Detecting of Genetic Diversity and Genetic Structure of Four Released Seeding Populations of Triangular Bream (Magalobrama terminalis) in the Middle and Lower Reaches of Qiantang River Based on SSR Marker

doi:10.11924/j.issn.1000-6850.casb2023-0395

Abstract

Abstract:

In order to detect genetic diversity and genetic structure of released populations of triangular bream (Magalobrama terminalis) in the middle and lower reaches of the Qiantang River, in this study, 192 samples from 4 Released seeding populations were amplified and sequenced based on 9 pairs of fluorescence-labeled SSR primers using capillary electrophoresis genotyping technique. The results showed that the average number of alleles (Na) was (20±2)-(24±2) per locus for each population, and the average number of effective alleles (Ne) was (10.76±1.51)-(13.66±1.33). The number of Na and Ne were the largest in Yuhang population, and the two indexes were similar in Jiande and Tonglu population. The average value of Shannon’s index (I) was (2.55±0.13)-(2.83±0.09) for four populations. The average value of observed heterozygosity (Ho) was (0.83±0.08)-(0.94±0.03) and that of expected heterozygosity (He) was (0.88±0.02)-(0.92±0.01), which showed heterozygote excess in most loci for four populations. The number of private alleles of Yuhang, Jiande, Tonglu and Fuyang populations was 54, 9, 15 and 17 respectively. The genetic diversity of the four released populations was relatively high, and some populations showed heterozygote excess and allelic loss. AMOVA analysis showed that 98.7% of genetic variation was between individuals within populations and 1.3% of genetic variation was among populations. The genetic differentiation index was 0.013-0.017, which revealed mild genetic differentiation (0<F_ST<0.05) among four released populations. Both structure analysis and PCoA analysis got the consistent result, indicating that all individuals could be divided into two genetic clusters. However, each cluster contained individuals from four populations, and there was no significant difference in the individual origin of each cluster. This study detected the germplasm genetic status of four released populations of M. terminalis, which would provide better guidance in the stock enhancement of M. terminalis in the future.

Key words: Magalobrama terminalis, released seeding, population, microsatellites, genetic diversity, genetic structure

ZHANG Minying, FANG Di’an, ZHOU Yanfeng, KUANG Zhen, REN Long, ZHENG Yuchen, XU Dongpo. Detecting of Genetic Diversity and Genetic Structure of Four Released Seeding Populations of Triangular Bream (Magalobrama terminalis) in the Middle and Lower Reaches of Qiantang River Based on SSR Marker[J]. Chinese Agricultural Science Bulletin, 2024, 40(13): 157-164.

Figures/Tables 8

References 36

[1]	陈宜瑜. 中国动物志:硬骨鱼纲鲤形目(中卷)[M]. 北京: 科学出版社,1998:200-208.
[2]	陈马康, 童合一, 俞泰济, 等. 钱塘江鱼类资源[M]. 上海: 上海科学技术文献出版社,1990.
[3]	童勇义. 钱塘江下游渔业资源开发利用现状及制约因素浅析[J]. 中国水产, 2001(4):18-19.
[4]	冯晓宇. 浙江省三角鲂研究进展及示范推广情况[J]. 科学养鱼, 2009, 25(10):40-41.
[5]	戴杨鑫, 冯晓宇, 储忝江, 等. 铜山水库人工繁殖三角鲂种群的生长情况[J]. 浙江农业科学, 2020, 61(2):365-367.
[6]	张敏莹, 王邢艳, 张婉平, 等. 基于微卫星亲权鉴定对钱塘江中下游三角鲂增殖放流效果的研究[J]. 农业生物技术学报, 2020, 28(7):1260-1268.
[7]	ARAKI H, SCHMID C. Is hatchery stocking a help or harm? Evidence, limitations and future directions in ecological and genetic surveys[J]. Aquaculture, 2010, 308(1):S2-S11.
[8]	李思发, 朱泽闻, 邹曙明, 等. 鲂属团头鲂、三角鲂及广东鲂种间遗传关系及种内遗传差异[J]. 动物学报, 2002(3):339-345.
[9]	童富淡, 刘红云. 三角鲂IGF-I基因的分子克隆及其序列分析[J]. 中国农业科学, 2004(10):1564-1570.
[10]	刘凯, 谢楠, 郭炜, 等. 三角鲂MHCⅠα基因全长cDNA克隆与生物信息学分析[J]. 浙江农业学报, 2022, 34(6):1162-1174. doi: 10.3969/j.issn.1004-1524.2022.06.07
[11]	苏晓磊, 郑国栋, 蒋霞云, 等. 三角鲂×翘嘴鲌、团头鲂×翘嘴鲌两种杂交后代微卫星遗传结构分析[J]. 水生生物学报, 2019, 43(2):264-271.
[12]	郭洪洪, 郑国栋, 吴成宾, 等. 三角鲂(♀)×翘嘴鲌(♂)杂交F_1的生长性能及形态差异[J]. 水产学报, 2018, 42(10):1572-1581.
[13]	王邢艳, 徐东坡, 张婉平, 等. 基于COⅠ基因研究钱塘江三角鲂亲本、放流和自然捕捞群体的遗传多样性与遗传分化[J]. 大连海洋大学学报, 2021, 36(4):603-611.
[14]	张敏莹, 方弟安, 周彦锋, 等. 基于微卫星标记的钱塘江中下游三角鲂亲本群体和自然捕捞群体遗传多样性及遗传结构[J]. 大连海洋大学学报, 2022, 37(5):775-783.
[15]	刘凯, 冯晓宇, 马恒甲, 等. 钱塘江三角鲂线粒体基因组测序及其结构特征分析[J]. 浙江农业学报, 2020, 32(9):1591-1608. doi: 10.3969/j.issn.1004-1524.2020.09.08
[16]	刘凯, 谢楠, 冯晓宇, 等. 三角鲂肝脏EST序列微卫星信息分析[J]. 江西农业学报, 2016, 28(6):87-91.
[17]	聂竹兰. 三角鲂转录组分析与不同地理种群遗传多样性研究[D]. 武汉: 华中农业大学, 2014.
[18]	郭爱环, 练青平, 罗伟, 等. 钱塘江三角鲂的生长特性及种群资源现状评价[J]. 渔业研究, 2022, 44(5):459-466. doi: 10.14012/j.cnki.fjsc.2022.05.004
[19]	孙效文, 张晓锋, 赵莹莹, 等. 水产生物微卫星标记技术研究进展及其应用[J]. 中国水产科学, 2008(4):689-703.
[20]	PEAKALL R, SMOUSE P E. GenAlEx 6: Genetic analysis in excel population genetic software for teaching and research[J]. Molecular ecology notes, 2006, 6(1):288-295. doi: 10.1111/men.2006.6.issue-1 URL
[21]	EXCOFFIE R L, LISCHE R H E L. Arlequin suite ver3.5: A new series of programs to perform population genetics analyses under Linux and Windows[J]. Molecular ecology resources, 2010, 10(3):564-567. doi: 10.1111/men.2010.10.issue-3 URL
[22]	SMOUSE P E, PEAKALL R. Spatial autocorrelation analysis of individual multiallele and multilocus genetic structure[J]. Heredity, 1999, 82(5):561-573. doi: 10.1038/sj.hdy.6885180
[23]	EARL D A, VONHOLDT B M. STRUCTURE HARVESTER: A website and program for visualizing STRUCTURE output and implementing the Evanno method[J]. Conservation genetics resources, 2012, 4(2):359-361. doi: 10.1007/s12686-011-9548-7 URL
[24]	EVANNO G, REGNAUT S, GOUDET J. Detecting the number of clusters of individuals using the software structure: A simulation study[J]. Molecular ecology, 2005, 14(8):2611-2620. doi: 10.1111/j.1365-294X.2005.02553.x pmid: 15969739
[25]	JAKOBSSON M, ROSENBERG N A. Clumpp: A cluster matching and permutation program for dealing with label switching and multimodality in analysis of population structure[J]. Bioinformatics, 2007, 23(14):1801-1806. doi: 10.1093/bioinformatics/btm233 pmid: 17485429
[26]	SWIFT M L. GraphPad prism, data analysis, and scientific graphing[J]. Journal of chemical information and computer sciences, 1997, 37(2):411-412. doi: 10.1021/ci960402j URL
[27]	BALLOUX F, LUGON-MOULIN N. The estimation of population differentiation with microsatellite markers[J]. Molecular ecology, 2002, 11(2):155-165. doi: 10.1046/j.0962-1083.2001.01436.x pmid: 11856418
[28]	李莉, 孙振兴, 扬树德, 等. 用微卫星标记分析皱纹盘鲍群体的遗传变异[J]. 遗传, 2006, 28(12):1549-1554.
[29]	DEWOODY J A, AVISE J C. Microsatellite variation in marine, freshwater and anadromous fishes compared with other animals[J]. Journal of fish biology, 2000, 56(3):461-473. doi: 10.1111/jfb.2000.56.issue-3 URL
[30]	WEIR B S, COCKERHAM C C. Estimation F-statistics for the analysis of population structure[J]. Evolution, 1984,38:1358-1370.
[31]	李盛文, 贾智英, 柏盈盈, 等. 松浦红镜鲤保种群体的遗传结构[J]. 中国水产科学, 2014, 21(1):67-74.
[32]	吉磊, 区又君, 李加. 卵形鲳鲹3个养殖群体的微卫星多态性分析[J]. 热带海洋学报, 2011, 30(3):62-68.
[33]	LUIKART G, CORNUET J M. Estimating the effective number of breeders from heterozygote excess in progeny[J]. Genetics, 1999, 151(3):1211-1216. doi: 10.1093/genetics/151.3.1211 pmid: 10049936
[34]	SPENCER C C, NEIGEL J E, LEBERG P L. Experimental evaluation of the usefulness of microsatellite DNA for detecting demographic bottlenecks[J]. Molecular ecology, 2000, 9(10):1517-1528. doi: 10.1046/j.1365-294x.2000.01031.x pmid: 11050547
[35]	唐首杰, 杨洁, 王成辉, 等. 驯养、选育条件下尼罗罗非鱼群体的选择压力分析[J]. 中国水产科学, 2016, 23(4):900-913.
[36]	ARAKI H, COOPER B, BLOUIN M S. Genetic effects of captive breeding cause a rapid, cumulative fitness decline in the wild[J]. Science, 2007,318:100-103.

位点	荧光标记	引物序列		重复序列	退火温度/℃
位点	荧光标记	F	R	重复序列	退火温度/℃
SJ4	FAM	TAAAGGGAAATTCTGGT	AGTATAAGTTGAGTGGGTG	(ATCT)₂₅	50
SJ9	HEX	GCCTGACAGTCTTCTGC	GCTATCCGATTATCATTTAC	(AC)₁₃	59
SJ14	FAM	AAAAAAAGAACGAGAGGG	TGATGATTTGGAGGAAGT	(GA)₁₅	50
SJ18	HEX	CCATCCAACACGCCGACTC	GCGGGAAATCACTTCATAACTCAA	(AC)₁₅	61.5
SJ21	TAM	AGTATAAGTTGAGTGGGTG	TAAAGGGAAATTCTGGT	(ATCT)₂₅	50
SJ22	FAM	ATGTATTGGGTTGAGGTT	GAGCTATGGACTCCGTTAT	(TG)₁₄	53
SJ23	FAM	CTTACAGACTCCGACAGG	ATCCACGACTTCCAGAAC	(AC)₁₂	57
SJ28	TAM	TGGAGTTAGTGTCCGCTTGT	AGGATACGGGTGAGTTCG	(TG)₁₃	56
SJ31	TAM	GAACCTGAATGAGGACATAAAGGT	TAACGCACAGAATTTCACATAACC	(TG)₈	53

位点	荧光标记	引物序列		重复序列	退火温度/℃
位点	荧光标记	F	R	重复序列	退火温度/℃
SJ4	FAM	TAAAGGGAAATTCTGGT	AGTATAAGTTGAGTGGGTG	(ATCT)₂₅	50
SJ9	HEX	GCCTGACAGTCTTCTGC	GCTATCCGATTATCATTTAC	(AC)₁₃	59
SJ14	FAM	AAAAAAAGAACGAGAGGG	TGATGATTTGGAGGAAGT	(GA)₁₅	50
SJ18	HEX	CCATCCAACACGCCGACTC	GCGGGAAATCACTTCATAACTCAA	(AC)₁₅	61.5
SJ21	TAM	AGTATAAGTTGAGTGGGTG	TAAAGGGAAATTCTGGT	(ATCT)₂₅	50
SJ22	FAM	ATGTATTGGGTTGAGGTT	GAGCTATGGACTCCGTTAT	(TG)₁₄	53
SJ23	FAM	CTTACAGACTCCGACAGG	ATCCACGACTTCCAGAAC	(AC)₁₂	57
SJ28	TAM	TGGAGTTAGTGTCCGCTTGT	AGGATACGGGTGAGTTCG	(TG)₁₃	56
SJ31	TAM	GAACCTGAATGAGGACATAAAGGT	TAACGCACAGAATTTCACATAACC	(TG)₈	53

群体	位点	等位基因数Na	有效等位基因数Ne	Shannon’s信息指数I	观测杂合度Ho	期望杂合度He	固定指数F	私有等位基因数
YH	SJ4	26	14.05	2.94	0.79	0.93	0.15	54
	SJ9	23	14.58	2.88	0.81	0.93	0.13
	SJ14	15	7.88	2.32	1.00	0.87	-0.15
	SJ18	30	19.78	3.17	1.00	0.95	-0.05
	SJ21	25	9.29	2.67	0.77	0.89	0.14
	SJ22	20	9.54	2.58	0.92	0.90	-0.02
	SJ23	25	15.11	2.93	1.00	0.93	-0.07
	SJ28	27	17.45	3.06	0.98	0.94	-0.04
	SJ31	28	15.26	2.98	1.00	0.93	-0.07
	均值	24±2	13.66±1.33	2.83±0.09	0.92±0.03	0.92±0.01	0.00±0.04
JD	SJ4	18	11.16	2.64	0.75	0.91	0.18	9
	SJ9	15	7.27	2.24	0.75	0.86	0.13
	SJ14	15	7.33	2.25	1.00	0.86	-0.16
	SJ18	29	18.14	3.11	1.00	0.94	-0.06
	SJ21	18	7.58	2.42	0.75	0.87	0.14
	SJ22	18	6.02	2.21	0.98	0.83	-0.17
	SJ23	24	14.77	2.91	0.94	0.93	-0.01
	SJ28	15	7.81	2.37	1.00	0.87	-0.15
	SJ31	31	19.38	3.16	1.00	0.95	-0.05
	均值	20±2	11.05±1.71	2.59±0.13	0.91±0.04	0.89±0.01	0.02±0.05
TL	SJ4	14	6.97	2.18	0.85	0.86	0.00	15
	SJ9	14	10.15	2.46	0.92	0.90	-0.02
	SJ14	14	7.25	2.21	1.00	0.86	-0.16
	SJ18	32	17.66	3.18	0.98	0.94	-0.04
	SJ21	15	4.81	2.03	0.75	0.79	0.05
	SJ22	19	10.17	2.56	0.98	0.90	-0.09
	SJ23	23	12.84	2.78	0.98	0.92	-0.06
	SJ28	17	9.33	2.44	0.98	0.89	-0.10
	SJ31	31	17.66	3.10	1.00	0.94	-0.06
	均值	20±2	10.76±1.51	2.55±0.13	0.94±0.03	0.89±0.02	0.05±0.02
FY	SJ4	18	5.20	2.12	0.48	0.81	0.41	17
	SJ9	14	8.83	2.32	0.67	0.89	0.25
	SJ14	17	9.29	2.45	1.00	0.89	-0.12
	SJ18	35	23.88	3.34	1.00	0.96	-0.04
	SJ21	21	4.71	2.15	0.46	0.79	0.42
	SJ22	15	6.81	2.18	0.96	0.85	-0.12
	SJ23	24	9.80	2.67	1.00	0.90	-0.11
	SJ28	18	11.05	2.59	0.96	0.91	-0.05
	SJ31	34	23.04	3.31	0.94	0.96	0.02
	均值	22±3	11.40±2.38	2.57±0.16	0.83±0.08	0.88±0.02	0.07±0.07

群体	位点	等位基因数Na	有效等位基因数Ne	Shannon’s信息指数I	观测杂合度Ho	期望杂合度He	固定指数F	私有等位基因数
YH	SJ4	26	14.05	2.94	0.79	0.93	0.15	54
	SJ9	23	14.58	2.88	0.81	0.93	0.13
	SJ14	15	7.88	2.32	1.00	0.87	-0.15
	SJ18	30	19.78	3.17	1.00	0.95	-0.05
	SJ21	25	9.29	2.67	0.77	0.89	0.14
	SJ22	20	9.54	2.58	0.92	0.90	-0.02
	SJ23	25	15.11	2.93	1.00	0.93	-0.07
	SJ28	27	17.45	3.06	0.98	0.94	-0.04
	SJ31	28	15.26	2.98	1.00	0.93	-0.07
	均值	24±2	13.66±1.33	2.83±0.09	0.92±0.03	0.92±0.01	0.00±0.04
JD	SJ4	18	11.16	2.64	0.75	0.91	0.18	9
	SJ9	15	7.27	2.24	0.75	0.86	0.13
	SJ14	15	7.33	2.25	1.00	0.86	-0.16
	SJ18	29	18.14	3.11	1.00	0.94	-0.06
	SJ21	18	7.58	2.42	0.75	0.87	0.14
	SJ22	18	6.02	2.21	0.98	0.83	-0.17
	SJ23	24	14.77	2.91	0.94	0.93	-0.01
	SJ28	15	7.81	2.37	1.00	0.87	-0.15
	SJ31	31	19.38	3.16	1.00	0.95	-0.05
	均值	20±2	11.05±1.71	2.59±0.13	0.91±0.04	0.89±0.01	0.02±0.05
TL	SJ4	14	6.97	2.18	0.85	0.86	0.00	15
	SJ9	14	10.15	2.46	0.92	0.90	-0.02
	SJ14	14	7.25	2.21	1.00	0.86	-0.16
	SJ18	32	17.66	3.18	0.98	0.94	-0.04
	SJ21	15	4.81	2.03	0.75	0.79	0.05
	SJ22	19	10.17	2.56	0.98	0.90	-0.09
	SJ23	23	12.84	2.78	0.98	0.92	-0.06
	SJ28	17	9.33	2.44	0.98	0.89	-0.10
	SJ31	31	17.66	3.10	1.00	0.94	-0.06
	均值	20±2	10.76±1.51	2.55±0.13	0.94±0.03	0.89±0.02	0.05±0.02
FY	SJ4	18	5.20	2.12	0.48	0.81	0.41	17
	SJ9	14	8.83	2.32	0.67	0.89	0.25
	SJ14	17	9.29	2.45	1.00	0.89	-0.12
	SJ18	35	23.88	3.34	1.00	0.96	-0.04
	SJ21	21	4.71	2.15	0.46	0.79	0.42
	SJ22	15	6.81	2.18	0.96	0.85	-0.12
	SJ23	24	9.80	2.67	1.00	0.90	-0.11
	SJ28	18	11.05	2.59	0.96	0.91	-0.05
	SJ31	34	23.04	3.31	0.94	0.96	0.02
	均值	22±3	11.40±2.38	2.57±0.16	0.83±0.08	0.88±0.02	0.07±0.07

变异来源	自由度	离均差平方和	均方差	变异的标准误	占总变异百分率/%
群体间	3	21.026	7.009	0.057	1.3
群体内	188	800.646	4.259	4.259	98.7
总计	191	821.672		4.316	100