GENETIC DIVERSITY OF THE BREEDING POPULATIONS OF GIANT FRESHWATER PRAWN MACROBRACHIUM ROSENBERGII
-
摘要: 为深入了解罗氏沼虾(Macrobrachium rosenbergii)不同种质资源的遗传背景, 研究采用微卫星标记和线粒体基因相结合, 对罗氏沼虾4个育种群体, 即选育3代的数丰核心群体(SF)、引进的正大群体(ZD)、正大和数丰杂交的群体(ZDS)、正大子代和数丰杂交的群体(ZD2S)的遗传多样性进行了研究。150个个体的微卫星分析结果显示, 7个微卫星位点均表现出高度多态性(PIC>0.5), 4个群体平均的等位基因数(Na)、期望杂合度(He)和多态信息含量(PIC)分别为19.43、0.8980和0.8867。SF群体的平均He (0.874)和PIC (0.854)最高, ZD2S的He (0.863)和PIC (0.834)其次, ZDS群体的He (0.798)和PIC (0.761)最低。4个群体间的遗传分化指数(FST)为0.04166—0.10438, 处于中低水平的遗传分化, 其中, ZD和ZDS的遗传分化最大, FST为0.10438。线粒体COⅠ和12S rRNA基因组合序列分析结果显示, 149个个体共识别27个单倍型, 平均单倍型多样性(Hd)和核苷酸多样性(π)分别为0.846和0.00313。在4个群体中, ZD2S群体的Hd值最高, 其次为SF群体, ZDS的最低; 对于π值, SF群体的最高, 其次为ZD2S群体, ZD群体的则最低, 该结果与微卫星的结果基本一致。但基于线粒体基因的群体间遗传分化较小, FST值为–0.02226—0.07310, 小于微卫星估计的结果。微卫星和线粒体基因一致表明, SF和ZD2S两个群体均保持着较高的遗传多样性, 具有进一步选育的潜力。Abstract: To learn the genetic background of the giant freshwater prawn (Macrobrachium rosenbergii), four different breeding populations, including Shufeng nucleus breeding population selected for three generations (SF), introduced Zhengda population (ZD), hybrid population between the introduced Zhengda population and Shufeng population (ZDS), and hybrid population between the progeny of the introduced Zhengda population and Shufeng population (ZD2S), were investigated by using both microsatellite loci and mitochondrial COⅠand 12S rRNA genes as molecular markers. The microsatellite data of 150 individuals revealed that all the seven loci had high diversity with PIC>0.5. The average number of alleles (Na), expected heterozygosity (He) and polymorphism information content (PIC) were 19.43, 0.8980 and 0.8867, respectively. Among seven loci of four populations, the average He (0.874) and PIC (0.854) of the SF were the highest, followed by He (0.863) and PIC (0.834) of the ZD2S, and the lowest He (0.798) and PIC (0.761) of the ZDS. The pairwised genetic differentiation indices (FST) among the four populations was 0.04166—0.10438, suggesting low to moderate genetic differentiations. The FST between ZD and ZDS was the highest with a value of 0.10438. Based on the concatenated mitochondrial COⅠ and 12S rRNA gene sequences, 149 individuals contained 27 haplotypes, with the average haplotype diversity (Hd) and nucleotide diversity (π) of 0.846 and 0.00313, respectively. Among four populations, the Hd of the ZD2S was the highest, followed by the SF, and ZDS was the lowest. For the π values, the SF was the highest, followed by ZD2S, and ZD was the lowest, which was generally consistent with the genetic diversity results analyzed based on microsatellite data. However, the genetic differentiation among four populations based on mitochondrial genes was low with FST values of –0.02226—0.07310, which is much lower than those based on microsatellite. These results indicated that both SF and ZD2S maintain high genetic diversity and have the potential for further selective breeding of Macrobrachium rosenbergii.
-
Keywords:
- Macrobrachium rosenbergii /
- Microsatellite loci /
- COⅠ /
- 12S rRNA /
- Genetic diversity /
- Breeding
-
罗氏沼虾(Macrobrachium rosenbergii)原产于东南亚的热带和亚热带地区, 是淡水虾类养殖中个体最大的品种[1]。中国大陆自1976年引进以来, 已产生了可观的经济效益并拥有良好的发展前景。但长期以来, 由于人工养殖中亲本管理不严谨(如近交等)导致种质衰退严重(如生长缓慢、病害频发等)[2], 尤其是近年来发生的长不大、性成熟提前的“铁虾”, 使罗氏沼虾产业遭受了巨大危害。因此, 进一步培育满足产业与市场需求的经济性状优良的罗氏沼虾新品种或新品系迫在眉睫。遗传多样性是物种进化的基础, 与进化潜力和适应力呈正相关, 遗传多样性越高的群体, 其生产性状如活力、生长速度、产卵量、环境适应力及抗病性能等提高的潜力越大[3]。因此, 遗传多样性水平不仅制约着生物的适应性进化, 也是选择育种的前提。
微卫星(Microsatellite)即短串联重复序列(Short tandem repeats, STRs)或简单重复序列(Simple sequence repeats, SSR), 是一种广泛应用于群体遗传多样性检测及亲缘关系鉴定的分子标记[4, 5]。线粒体DNA (mitochondrial DNA, mtDNA)是一种核外的遗传物质, 因具有结构简单、进化速度适中、易扩增等优点, 被认为是研究群体遗传多样性和遗传结构的有力工具[6, 7]。目前, 针对罗氏沼虾不同野生或养殖群体, 已有一些采用微卫星标记[3, 8—12]或线粒体基因[13—15]进行遗传多样性的研究报道, 但这些研究多集中在不同的养殖群体上, 亲本来源背景未知, 且部分研究采用的群体较少, 如钟丹丹等[3]仅研究了2个群体。另外, 在已有的线粒体基因研究中, 各群体覆盖的个体数极少, 如杨学明等[13] 3个群体仅17个个体, 姚茜等[14] 3个群体仅16个个体, 且研究采用的基因片段大多数较短, 约400—500 bp, 提供的信息量有限。目前罗氏沼虾中还未见将微卫星和线粒体基因相结合对群体遗传多样性进行研究的报道。微卫星和线粒体分别是细胞核和细胞质遗传, 将二者相结合进行研究, 能更全面地反映生物的遗传多样性, 应用较为广泛, 如东北林蛙(Rana dybowskii)[16]、中华蜜蜂(Apiscerana cerana)[17]、梨小食心虫(Grapholita molesta)[18]等物种的遗传多样性研究等。
本研究将微卫星和线粒体基因相结合, 对罗氏沼虾不同种质资源的育种群体进行遗传多样性研究, 期望全面了解国内已有种质资源的遗传背景, 为优良品种选育提供参考资料。
1. 材料与方法
1.1 实验材料
实验材料为取自罗氏沼虾育种基地(江苏数丰水产种业有限公司)培育及收集的不同来源的4个育种群体, 包括: 选育3代的数丰核心群体(SF)、引进的正大群体(ZD)、正大和数丰杂交的群体(ZDS)、正大子代和数丰杂交的群体(ZD2S)。数丰核心群体因群体数量较大选取60个个体, 其他每个群体各取30个个体, 分别剪取虾尾肌肉组织, 用95%的乙醇保存, 并冻存于–20℃冰箱, 供后续DNA提取。
1.2 基因组DNA提取、PCR扩增和测序
肌肉组织基因组DNA提取采用高盐抽提法, 用1.2%琼脂糖凝胶电泳检测提取DNA的完整性, 微量分光光度计检测浓度和纯度后, 用灭菌水稀释至50 ng/μL供后续PCR扩增实验。
采用7对荧光修饰标记(HEX或FAM)的微卫星引物对各群体进行PCR扩增。微卫星引物来自文献[3, 8]及本团队开发的引物。引物序列、目的片段大小及扩增条件见表 1。PCR反应体系为20 μL, 包括: 2×Taq Mix (TaKaRa) 10 μL、10 μmol/L的上下游引物各1 μL、灭菌水 7 μL和DNA 1 μL。PCR反应程序为: 95℃预变性3min, 然后94℃变性30s、54—62℃退火30s、72℃延伸45s, 重复35次循环, 最后72℃延伸10min。
表 1 用于本研究的微卫星引物相关信息Table 1. Primers used in the present study位点Locus 引物序列Primer sequence (5′—3′) 退火温度Tm 片段大小Size range (bp) 来源Resource 6CT55 FAM-CAGCTCTAACCTGATTGAAAGAC 55 220—264 钟丹丹等[3] GCAGGCATTATCGTTACTTCTCC 12AG57 HEX-GGAACAAATGAAGACTTGGATGC 57.6 232—294 钟丹丹等[3] CTCTCTCTCCCTCAAGGTGTTGG 27AG62 FAM-CAACTCTATGTTTCGGCATTTGG 62 224—276 孙成飞等[8] GGGGAATTTTACCGATGTTTCTG 28GT62 HEX-CCACCTACCGTACATTCCCAAAC 62 203—291 孙成飞等[8] CGGGGCGACTTTTAGTATCGAC 29AG62 FAM-CAAGGCTCGTGTCTCTTGTTTC 62 274—322 孙成飞等[8] GCTTGTACTTGTTCAGCTTTTGC 47AAG54 HEX-AAGAGGATTTGGAGCGATTGG 54 169—229 本研究 CTGAGTAAGATACGACGCCTTC 50AG54 FAM-GTAATGAACAGCACGAAAAGGAAG 54 274—322 本研究 TCTGCGTTATTTTGAGTTTGGTATG 线粒体COⅠ和12S rRNA基因的引物根据罗氏沼虾GenBank中的线粒体基因组全序列(登录号为AY659990和NC006880)进行设计。PCR反应体系为50 μL, 其中包括: 灭菌水31 μL、10×buffer 5 μL、2.5 mmol/L的dNTP(TaKaRa) 4 μL、25 mmol/L的MgC12(TaKaRa) 3 μL、10 μmol/L的上下游引物分别为 1 μL、5 U/μL的Taq酶(TaKaRa) 0.2 μL和DNA 4.8 μL。PCR反应采用Touchdown程序: 95℃预变性5min; 然后95℃变性30s、56—62℃退火30s、72℃延伸30s, 重复10次循环, 每次循环降低1℃; 然后再95℃变性30s、52℃退火30s、72℃延伸30s, 重复25次循环; 最后72℃延伸20min。
PCR扩增产物用1.5%琼脂糖凝胶进行电泳检测。微卫星分型使用ABI3730XL全自动DNA测序仪进行, 获得每个个体在不同微卫星位点的等位基因数据, 采用软件GeneMarker V2.2.0对扩增产物片段大小进行读取。线粒体基因采用与PCR扩增相同的引物进行双向测序。以上测序工作均由武汉天一辉远生物科技公司完成。
1.3 数据处理
对获得的微卫星位点分型数据进行校正。利用Genepop 4.0软件进行哈迪温伯格平衡(Hardy-Weinberg equilibrium, HWE)检验。采用Cervus 3.0软件计算等位基因数(Na)、观测杂合度(Ho)、期望杂合度(He)和多态信息量(PIC)。利用Arlequin 3.5软件来估计各种群间遗传分化指数(FST), 进行分子变异分析(AMOVA), 包括群体内以及群体间变异。
线粒体基因序列用Clustal X2.1软件进行比对, 并在Seaview中进行手工校正。采用MEGA7.0.26对序列的基本信息进行统计, 如序列长度、碱基组成、变异位点的数目和简约信息位点数目等。利用软件DnaSP 5.0和Arlequin进行单倍型数目(Number of haplotype, H)、单倍型多样性(Haplotype diversity, Hd)和核苷酸多样性(Nucleotidae diversity, π)等参数估计, 并采用最大似然 (ML)法对单倍型序列进行聚类分析。
2. 结果
2.1 基于微卫星位点的罗氏沼虾育种群体遗传多样性
经检验, 7个微卫星位点都极显著偏离哈迪温伯格平衡(P<0.01)。所有群体所有微卫星位点的等位基因数(Na)为14—24个, 平均每个位点为19.43个等位基因, 平均观测杂合度(Ho)、期望杂合度(He)和多态信息含量(PIC)分别为0.615、0.898和0.887。在7个位点中, 座位29AG62的He和PIC都最高, 分别为0.919和0.910; 座位27AG62的He和PIC值均最低, 分别为0.878和0.864。4个群体中, SF群体平均每个座位的He和PIC都最高, 分别为0.874和0.854; ZD2S的其次, He和PIC分别为0.863和0.834; 而ZDS群体的平均He和PIC都最低, 分别为0.798和0.761 (表 2)。
表 2 各微卫星位点在4个罗氏沼虾育种群体中的遗传多样性参数Table 2. Genetic diversity parameters of each microsatellite locus in the four breeding populations of Macrobrachium rosenbergii位点Locus 所有群体All stocks SF (n=60) ZD (n=30) Na Ho He PIC Na Ho He PIC Na Ho He PIC 12AG57 19 0.460 0.879 0.867 17 0.350 0.858 0.837 13 0.400 0.900 0.874 47AAG54 14 0.873 0.891 0.878 13 0.867 0.851 0.827 9 0.800 0.818 0.782 27AG62 18 0.386 0.878 0.864 14 0.549 0.857 0.836 10 0.200 0.817 0.779 29AG62 24 0.813 0.919 0.910 19 0.717 0.907 0.892 14 0.867 0.897 0.871 50AG54 20 0.604 0.908 0.898 16 0.508 0.872 0.851 11 0.867 0.847 0.814 28GT62 21 0.456 0.904 0.893 15 0.424 0.882 0.863 17 0.367 0.899 0.874 6CT55 20 0.711 0.907 0.897 16 0.683 0.890 0.873 10 0.567 0.772 0.737 平均值Mean 19.439 0.615 0.898 0.887 15.71 0.585 0.874 0.854 12.00 0.581 0.850 0.819 位点Locus ZDS (n=30) ZD2S (n=30) Na Ho He PIC Na Ho He PIC 12AG57 10 0.600 0.755 0.718 14 0.600 0.867 0.842 47AAG54 8 1.000 0.848 0.812 12 0.833 0.865 0.835 27AG62 5 0.182 0.555 0.499 10 0.448 0.816 0.781 29AG62 16 0.900 0.910 0.887 20 0.867 0.937 0.917 50AG54 11 0.567 0.880 0.852 12 0.567 0.875 0.846 28GT62 9 0.767 0.840 0.803 12 0.300 0.877 0.849 6CT55 8 0.900 0.801 0.756 9 0.724 0.806 0.770 平均值Mean 9.57 0.702 0.798 0.761 12.71 0.620 0.863 0.834 注: 表中各遗传多样性参数简写如下, Na. 等位基因数; Ho. 观察杂合度; He. 期望杂合度; PIC. 多态信息含量Note: the abbreviation of the genetic diversity parameters as the following. Na. number of alleles; Ho. observed heterozygosity; He. expected heterozygosity; PIC. polymorphism information content 各育种群体间及群体内分子变异分析(AMOVA)结果表明, 罗氏沼虾各育种群体的变异主要来自群体内的变异, 占总变异的93.59%, 而来自群体间的变异较小, 仅占总变异的6.31% (表格未列出)。各群体间的遗传分化指数(FST)见表 3, 群体间的遗传分化均达极显著水平。其中, ZD和ZDS群体间的遗传分化最大, FST值为0.10438, 而SF和ZD2S间的分化最小, FST值为0.04166。ZD和ZD2S间及ZD和SF群体间的FST值分别为0.06539和0.07281。
表 3 基于微卫星位点的罗氏沼虾各育种群体间遗传分化指数Table 3. Genetic differentiation index (FST) of the breeding populations of Macrobrachium rosenbergii based on SSRZDS ZD2S SF ZD ZDS 0.00000 ZD2S 0.06451** 0.00000 SF 0.05520** 0.04166** 0.00000 ZD 0.10438** 0.06539** 0.07281** 0.00000 注: **表示P≤0.01Note: ** represents P≤0.01 2.2 基于线粒体基因的罗氏沼虾育种群体遗传多样性
共获得4个群体117个个体的线粒体COⅠ基因及137个个体的12S rRNA基因序列。COⅠ基因比对长度为1247 bp, 12S rRNA基因比对长度为798 bp。COⅠ基因共检测到21个变异位点, 其中14个为简约信息位点; 117个个体共享9个单倍型, 单倍型多样性(Hd)为0.425, 核苷酸多样性(π)为0.00371(表 4)。在798 bp的12S rRNA基因中, 共检测到10个变异位点, 其中9个为简约信息位点; 137个个体共识别12个单倍型, Hd为0.660, π为0.0023(表 4)。
表 4 罗氏沼虾育种群体线粒体COⅠ、12S rRNA基因及两个基因组合序列的遗传多样性参数Table 4. Genetic diversity parameters of mitochondrial COⅠ, 12S rRNA and their concatenated gene sequences in each breeding population of Macrobrachium rosenbergii群体Population COⅠ 12S rRNA COⅠ+ 12S rRNA N H Hd π N H Hd π N H Hd π ZDS 19 3 0.205 0.00193 28 2 0.138 0.00086 29 5 0.594 0.001 ZD2S 26 5 0.603 0.00455 26 4 0.637 0.00197 30 11 0.876 0.0035 ZD 25 2 0.080 0.00006 28 3 0.500 0.00078 30 7 0.713 0.0004 SF 47 5 0.530 0.00512 55 10 0.598 0.00293 60 19 0.835 0.0042 总计Total 117 9 0.425 0.00371 137 12 0.660 0.00230 149 27 0.846 0.0031 注: 表中各参数简写如下, N. 样本数; H. 单倍型数; Hd. 单倍型多样性; π. 核苷酸多样性Note: the abbreviation of each parameter as the following. N. sample size; H. number of haplotype; Hd. haplotype diversity; π. nucleotidae diversity 将两个基因进行组合分析, 获得149个个体2045 bp的组合序列, 含31个变异位点, 其中27个为简约信息位点; 共识别27个单倍型, Hd值和π值分别为0.846和0.00313; 在4个群体中, ZD2S的Hd值最高(0.876), 其次为SF群体(0.835), ZDS的Hd值最低(0.594); 对于π值, SF群体的最高(0.00422), 其次为ZD2S(0.00349), ZD群体的最低(0.00044, 表 4)。
在两个基因组合序列获得的27个单倍型中, Hap6包含的个体数最多, 有46个, 出现频率为30.87%, 分布于ZD2S、SF和ZD三个群体中; 其次为单倍型Hap2, 包含29个个体, 出现频率为19.46%, 涵盖了所有4个群体。Hap1和Hap9分别包含15和13个个体, 其他23个单倍型包含的个体数在1—9个, 有17个单倍型只包含1个个体, 为特有单倍型。在各群体中单倍型数量不一, 其中ZDS群体29个个体定义了5个单倍型, ZD2S群体30个个体定义了11个单倍型, SF群体60个个体定义了19个单倍型, ZD群体30个个体定义了7个单倍型。各群体单倍型分布情况见表 5。
表 5 罗氏沼虾各育种群体COⅠ和12S rRNA基因组合序列单倍型在各群体中的分布Table 5. Distribution of combined COⅠ and 12S rRNA sequence haplotypes in each breeding population of Macrobrachium rosenbergii单倍型
HaplotypeZDS ZD2S SF ZD 总个体数
Total相对频率Relative frequency (%) Hap1 10 2 3 15 10.07 Hap2 16 5 2 6 29 19.46 Hap3 1 1 0.67 Hap4 1 1 0.67 Hap5 1 1 0.67 Hap6 9 22 15 46 30.87 Hap7 1 1 0.67 Hap8 1 1 0.67 Hap9 2 9 2 13 8.72 Hap10 1 1 0.67 Hap11 3 1 4 2.68 Hap12 1 1 0.67 Hap13 1 1 0.67 Hap14 1 1 0.67 Hap15 2 2 1.34 Hap16 1 1 0.67 Hap17 1 1 0.67 Hap18 1 1 0.67 Hap19 2 2 1.34 Hap20 1 2 3 6 4.03 Hap21 1 1 0.67 Hap22 6 6 4.03 Hap23 1 2 4 2 9 6.04 Hap24 1 1 0.67 Hap25 1 1 0.67 Hap26 1 1 0.67 Hap27 1 1 0.67 总个体数Total 29 30 60 30 149 100 基于最大似然法构建的27个单倍型聚类关系见图 1。所有单倍型聚为两个分支, Clade A和Clade B。Clade A包含16个单倍型, 主要是来自SF群体的个体, 同时有部分ZD2S的个体及个别ZD和ZDS的个体; ZD群体的单倍型主要聚类在Clade B, 该分支也包含较多的SF和ZD2S群体及部分ZDS群体的单倍型。所有单倍型间的遗传距离介于0.001—0.014, 两个分支Clade A和Clade B间的平均遗传距离为0.009。从整体上看, 单倍型间的遗传差异较小。
![]()
图 1 罗氏沼虾COⅠ和12S rRNA基因组合序列单倍型间的聚类分析节点处的数值代表支持率, 括弧中代表对应单倍型包含的种群Figure 1. Phylogenetic tree for the combined COⅠ and 12S rRNA sequence haplotypes of Macrobrachium rosenbergiiValues at the nodes correspond to the support values, and the populations included by the corresponding haplotype are listed in the parentheses基于COⅠ和12S rRNA基因组合序列的各群体两两间遗传分化指数(FST)结果表明, 群体间遗传分化较小, FST值为–0.02226—0.07310。ZD2S群体和ZDS群体间遗传分化最大, FST值为0.07310, 且分化显著; ZDS和ZD群体间及ZDS和SF群体间的FST值分别为0.03899和0.01419, 遗传分化较小, 且分化不显著; ZD和ZDS群体间、ZD和SF群体间及ZD2S和SF群体间的FST值为负值(表 6), 暗示群体间几乎不存在分化。
表 6 基于COⅠ和12S rRNA基因组合序列的罗氏沼虾育种群体间遗传分化指数Table 6. Genetic differentiation index (FST) of the breeding populations of Macrobrachium rosenbergii based on the combined COⅠ and 12S rRNA sequencesZDS ZD ZD2S SF ZDS 0 ZD 0.03899 0 ZD2S 0.07310* –0.02226 0 SF 0.01419 –0.01611 –0.00472 0 注: *表示P≤0.05Note: * represents P≤0.05 3. 讨论
3.1 基于微卫星位点的多态性及群体遗传分化
在本研究中哈迪温伯格平衡(HWE)检验结果显示, 7个微卫星座位都显著偏离哈迪温伯格平衡(P<0.05), 这与罗氏沼虾的其他微卫星研究结果类似[3]。这可能是因为研究对象为人工选育群体, 群体间存在非随机交配等导致偏离哈迪温伯格平衡[19]。另一方面, 这也可能与无效等位基因相关[20]。
多态信息含量(PIC) 能反映群体的遗传变异程度和位点多样性等[21]。根据Botstein等[22]的划分标准, PIC≥0.5为高度多态性; 0.25≤PIC<0.5为中度多态性; PIC< 0.25为低度多态性。在本研究结果中, 7个微卫星位点的PIC值除ZDS群体27AG62位点为0.499外, 其余各群体每个位点的PIC值均大于0.5, 为高度多态性。在4个群体中, SF群体平均每个座位的PIC最高(0.854), ZD2S的其次(0.834)。
杂合度是衡量群体遗传多样性的重要参数。本研究中的4个罗氏沼虾群体均表现出较高的杂合度, 其中SF群体平均每个位点的期望杂合度(He)最高, 为0.874, ZD2S群体的其次, 为0.863。从整体上, 所有群体所有位点的平均观测杂合度(Ho 0.615)小于平均期望杂合度(0.898), 表明这4个群体可能因选育丢失了部分杂合子。杂合子缺失将导致纯合子增加, 提高有害等位基因纯合的几率, 近交衰退几率升高, 从而降低物种对环境变化的适应性[23]。因此, 在缺乏外来优势野生型群体作为育种亲本时, 选择杂合子较多的群体进行育种可以降低近交衰退几率。
遗传分化指数(FST)是反映各群体间遗传分化的重要指标[24]。本研究基于7个微卫星位点分析的结果表明, 除SF与ZD2S群体的FST值稍低于0.05外, 其余两两群体间, FST值都在0.05—0.15, 为中度分化[24], 且分化均达极显著水平, 这与孙成飞等[8]的研究结果类似。选择遗传分化指数大的群体间进行杂交, 可发挥杂种优势, 减缓育种群体之间近亲繁殖带来的影响。
3.2 基于线粒体基因的遗传多样性及群体遗传分化
在本研究中, 线粒体COⅠ、12S rRNA及两个基因的组合序列分析结果一致表明, 罗氏沼虾SF群体的核苷酸多样性(π)和ZD2S群体单倍型多样性最高(Hd), 与微卫星位点分析的结果较为一致。4个育种群体COⅠ基因的单倍型多样性(Hd)为0.080—0.603, 核苷酸多样(π)为0.00006—0.00512, 12S rRNA基因的Hd为0.500—0.637, π为0.00078—0.00293, 两个基因组合序列的Hd为0.594—0.876, π为0.00044—0.00422, 均高于Thanh等[15]中罗氏沼虾越南2个野生种群及中国10个养殖群体16S rRNA基因的Hd (0—0.303)和π (0—0.0008), 但类似于其COⅠ基因的Hd (0.129—0.694)和π(0.0003—0.0073)。在本研究中, SF群体是通过收集国内包括罗氏沼虾“南太湖2号”在内的不同种质资源进行选育3代的群体, 亲本来源较为丰富。引进的ZD群体与SF群体杂交后也获得较高的遗传多样性, 这在本研究的微卫星和线粒体基因分析结果中得到进一步证明。另一方面, Thanh等[15]的研究采用的16S rRNA及COⅠ基因片段较短(分别仅395和498 bp), 提供的信息量有限。
从两个基因组合序列的单倍型分布来看, 27个单倍型中具有优势单倍型Hap6, 包含46个个体, 出现频率为30.87%, 分布于除ZDS群体以外的3个群体; 其次占优势的单倍型为Hap2, 包含29个个体, 出现频率为19.46%, 4个群体中均有出现; 其他各单倍型包含的个体数1—15不等。由此说明, 单倍型在各群体中分布不均, 优势单倍型与其他单倍型包含的个体数量相差极大。从各群体的单倍型多样性来看, SF群体包含的单倍型数目最多, 共有19个单倍型, ZD2S群体次之, 共11个。通常来说, 优势单倍型的个体适应能力强, 拥有更大的繁殖群体, 而低频率的单倍型则因环境适应能力差, 拥有的繁殖群体较小, 更容易对整个种群的遗传多样性造成不利影响[25]。另一方面, 优势单倍型体现了前期人工选择过程中的育种方向, 定向选育过程中保留了成活率高、生长性状优良的个体, 从而使优势单倍型包含的群体数量逐渐扩大。另外, 单倍型的聚类分析结果显示, SF和ZD2S群体拥有的单倍型, 在两个分支中均有较多的分布, 也体现出这两个群体具有丰富的多样性。
3.3 微卫星和线粒体基因的群体遗传多样性及遗传分化比较
本研究基于微卫星位点的等位基因数目(Na)、期望杂合度(He)和多态信息含量(PIC)等遗传多样性参数估计结果表明, 在4个育种群体中, SF群体的遗传多样性最高, 其次为ZD2S群体。该结果与基于线粒体COⅠ和12S rRNA基因组合序列的单倍型多样性(Hd)、核苷酸多样性(π)及单倍型聚类分析等结果较为一致, 均表明SF和ZD2S群体的遗传多样性较为丰富, 具有进一步选育的潜力。
从群体间的遗传分化来看, 微卫星位点分析结果显示各群体两两间均存在极显著的遗传分化, 且几乎均为中等程度分化, 而两个线粒体基因组合序列分析则显示群体间的遗传分化较小, 除ZD2S群体和ZDS群体的FST值(0.07310)大于0.05属于中度分化外, 其余均较小或为负值, 几乎不存在分化。这表明微卫星标记用于检测近缘种群间的遗传分化时比线粒体基因具有更高的灵敏度。
综合微卫星和线粒体基因分析结果, 在罗氏沼虾4个育种群体中, SF和ZD2S两个群体具有较高的遗传多样性, 可作为良好的选育材料。但已有的这4个群体间遗传分化不大, 若能引进遗传多样性更为丰富的、与这4个群体具有高度分化的种质资源, 将更有利于优良品种的选育。
-
![]()
图 1 罗氏沼虾COⅠ和12S rRNA基因组合序列单倍型间的聚类分析
节点处的数值代表支持率, 括弧中代表对应单倍型包含的种群
Figure 1. Phylogenetic tree for the combined COⅠ and 12S rRNA sequence haplotypes of Macrobrachium rosenbergii
Values at the nodes correspond to the support values, and the populations included by the corresponding haplotype are listed in the parentheses
表 1 用于本研究的微卫星引物相关信息
Table 1 Primers used in the present study
位点Locus 引物序列Primer sequence (5′—3′) 退火温度Tm 片段大小Size range (bp) 来源Resource 6CT55 FAM-CAGCTCTAACCTGATTGAAAGAC 55 220—264 钟丹丹等[3] GCAGGCATTATCGTTACTTCTCC 12AG57 HEX-GGAACAAATGAAGACTTGGATGC 57.6 232—294 钟丹丹等[3] CTCTCTCTCCCTCAAGGTGTTGG 27AG62 FAM-CAACTCTATGTTTCGGCATTTGG 62 224—276 孙成飞等[8] GGGGAATTTTACCGATGTTTCTG 28GT62 HEX-CCACCTACCGTACATTCCCAAAC 62 203—291 孙成飞等[8] CGGGGCGACTTTTAGTATCGAC 29AG62 FAM-CAAGGCTCGTGTCTCTTGTTTC 62 274—322 孙成飞等[8] GCTTGTACTTGTTCAGCTTTTGC 47AAG54 HEX-AAGAGGATTTGGAGCGATTGG 54 169—229 本研究 CTGAGTAAGATACGACGCCTTC 50AG54 FAM-GTAATGAACAGCACGAAAAGGAAG 54 274—322 本研究 TCTGCGTTATTTTGAGTTTGGTATG 
下载: 导出CSV
表 2 各微卫星位点在4个罗氏沼虾育种群体中的遗传多样性参数
Table 2 Genetic diversity parameters of each microsatellite locus in the four breeding populations of Macrobrachium rosenbergii
位点Locus 所有群体All stocks SF (n=60) ZD (n=30) Na Ho He PIC Na Ho He PIC Na Ho He PIC 12AG57 19 0.460 0.879 0.867 17 0.350 0.858 0.837 13 0.400 0.900 0.874 47AAG54 14 0.873 0.891 0.878 13 0.867 0.851 0.827 9 0.800 0.818 0.782 27AG62 18 0.386 0.878 0.864 14 0.549 0.857 0.836 10 0.200 0.817 0.779 29AG62 24 0.813 0.919 0.910 19 0.717 0.907 0.892 14 0.867 0.897 0.871 50AG54 20 0.604 0.908 0.898 16 0.508 0.872 0.851 11 0.867 0.847 0.814 28GT62 21 0.456 0.904 0.893 15 0.424 0.882 0.863 17 0.367 0.899 0.874 6CT55 20 0.711 0.907 0.897 16 0.683 0.890 0.873 10 0.567 0.772 0.737 平均值Mean 19.439 0.615 0.898 0.887 15.71 0.585 0.874 0.854 12.00 0.581 0.850 0.819
下载: 导出CSV
位点Locus ZDS (n=30) ZD2S (n=30) Na Ho He PIC Na Ho He PIC 12AG57 10 0.600 0.755 0.718 14 0.600 0.867 0.842 47AAG54 8 1.000 0.848 0.812 12 0.833 0.865 0.835 27AG62 5 0.182 0.555 0.499 10 0.448 0.816 0.781 29AG62 16 0.900 0.910 0.887 20 0.867 0.937 0.917 50AG54 11 0.567 0.880 0.852 12 0.567 0.875 0.846 28GT62 9 0.767 0.840 0.803 12 0.300 0.877 0.849 6CT55 8 0.900 0.801 0.756 9 0.724 0.806 0.770 平均值Mean 9.57 0.702 0.798 0.761 12.71 0.620 0.863 0.834 注: 表中各遗传多样性参数简写如下, Na. 等位基因数; Ho. 观察杂合度; He. 期望杂合度; PIC. 多态信息含量Note: the abbreviation of the genetic diversity parameters as the following. Na. number of alleles; Ho. observed heterozygosity; He. expected heterozygosity; PIC. polymorphism information content
下载: 导出CSV
表 3 基于微卫星位点的罗氏沼虾各育种群体间遗传分化指数
Table 3 Genetic differentiation index (FST) of the breeding populations of Macrobrachium rosenbergii based on SSR
ZDS ZD2S SF ZD ZDS 0.00000 ZD2S 0.06451** 0.00000 SF 0.05520** 0.04166** 0.00000 ZD 0.10438** 0.06539** 0.07281** 0.00000 注: **表示P≤0.01Note: ** represents P≤0.01
下载: 导出CSV
表 4 罗氏沼虾育种群体线粒体COⅠ、12S rRNA基因及两个基因组合序列的遗传多样性参数
Table 4 Genetic diversity parameters of mitochondrial COⅠ, 12S rRNA and their concatenated gene sequences in each breeding population of Macrobrachium rosenbergii
群体Population COⅠ 12S rRNA COⅠ+ 12S rRNA N H Hd π N H Hd π N H Hd π ZDS 19 3 0.205 0.00193 28 2 0.138 0.00086 29 5 0.594 0.001 ZD2S 26 5 0.603 0.00455 26 4 0.637 0.00197 30 11 0.876 0.0035 ZD 25 2 0.080 0.00006 28 3 0.500 0.00078 30 7 0.713 0.0004 SF 47 5 0.530 0.00512 55 10 0.598 0.00293 60 19 0.835 0.0042 总计Total 117 9 0.425 0.00371 137 12 0.660 0.00230 149 27 0.846 0.0031 注: 表中各参数简写如下, N. 样本数; H. 单倍型数; Hd. 单倍型多样性; π. 核苷酸多样性Note: the abbreviation of each parameter as the following. N. sample size; H. number of haplotype; Hd. haplotype diversity; π. nucleotidae diversity
下载: 导出CSV
表 5 罗氏沼虾各育种群体COⅠ和12S rRNA基因组合序列单倍型在各群体中的分布
Table 5 Distribution of combined COⅠ and 12S rRNA sequence haplotypes in each breeding population of Macrobrachium rosenbergii
单倍型
HaplotypeZDS ZD2S SF ZD 总个体数
Total相对频率Relative frequency (%) Hap1 10 2 3 15 10.07 Hap2 16 5 2 6 29 19.46 Hap3 1 1 0.67 Hap4 1 1 0.67 Hap5 1 1 0.67 Hap6 9 22 15 46 30.87 Hap7 1 1 0.67 Hap8 1 1 0.67 Hap9 2 9 2 13 8.72 Hap10 1 1 0.67 Hap11 3 1 4 2.68 Hap12 1 1 0.67 Hap13 1 1 0.67 Hap14 1 1 0.67 Hap15 2 2 1.34 Hap16 1 1 0.67 Hap17 1 1 0.67 Hap18 1 1 0.67 Hap19 2 2 1.34 Hap20 1 2 3 6 4.03 Hap21 1 1 0.67 Hap22 6 6 4.03 Hap23 1 2 4 2 9 6.04 Hap24 1 1 0.67 Hap25 1 1 0.67 Hap26 1 1 0.67 Hap27 1 1 0.67 总个体数Total 29 30 60 30 149 100
下载: 导出CSV
表 6 基于COⅠ和12S rRNA基因组合序列的罗氏沼虾育种群体间遗传分化指数
Table 6 Genetic differentiation index (FST) of the breeding populations of Macrobrachium rosenbergii based on the combined COⅠ and 12S rRNA sequences
ZDS ZD ZD2S SF ZDS 0 ZD 0.03899 0 ZD2S 0.07310* –0.02226 0 SF 0.01419 –0.01611 –0.00472 0 注: *表示P≤0.05Note: * represents P≤0.05
下载: 导出CSV
-
[1] 潘家模. 罗氏沼虾养殖新技术 [M]. 上海: 上海科学技术出版社, 1994: 1-6 Pan J M. New Cultured Technique of Macrobranchium rosenbergii [M]. Shanghai: Science and Technology Press, 1994: 1-6
[2] 陈雪峰, 杨国梁, 孔杰, 等. 人工养殖与选育对罗氏沼虾遗传多样性的影响 [J]. 水生生物学报, 2012, 36(5): 866-873. Chen X F, Yang G L, Kong J, et al. Effect of artificial culture and selective breeding on the genetic diversity of Macrobrachium rosenbergii [J]. Acta Hydrobiologica Sinica, 2012, 36(5): 866-873.
[3] 钟丹丹, 林勇, 宾石玉, 等. 两个罗氏沼虾种群的遗传多样性研究 [J]. 广东农业科学, 2015, 42(24): 140-145. doi: 10.3969/j.issn.1004-874X.2015.24.027 Zhong D D, Ln Y, Bin S Y, et al. Genetic diversity of two giant freshwater prawn Macrobrachium rosenbergii populations [J]. Guangdong Agricultural Sciences, 2015, 42(24): 140-145. doi: 10.3969/j.issn.1004-874X.2015.24.027
[4] 陈万光, 贾晓惠. 微卫星标记在水产动物遗传研究方面的应用 [J]. 安徽农业科学, 2009, 37(12): 5659-5660. doi: 10.3969/j.issn.0517-6611.2009.12.034 Chen W G, Jia X H. Application of microsatellite markers in genetic research of aquatic animals [J]. Journal of Anhui Agricultural Sciences, 2009, 37(12): 5659-5660. doi: 10.3969/j.issn.0517-6611.2009.12.034
[5] 张于光, 李迪强, 肖启明, 等. 微卫星技术及在动物遗传多样性研究中的应用 [J]. 湖南农业大学学报, 2001, 27(5): 410-414. Zhang Y G, Li D Q, Xiao Q M, et al. Microsatellites and their application to the genetic diversity in animal [J]. Journal of Hunan Agricultural University, 2001, 27(5): 410-414.
[6] 方李宏, 薛俊增, 董双林. 甲壳动物线粒体DNA的研究 [J]. 海洋湖沼通报, 2004(2): 59-65. doi: 10.3969/j.issn.1003-6482.2004.02.009 Fang L H, Xue J Z, Dong S L. A review on the studies of mitochondrial DNA (mtDNA) of crustacean [J]. Transactions of Oceanology and Limnology, 2004(2): 59-65. doi: 10.3969/j.issn.1003-6482.2004.02.009
[7] 申欣, 孙名安. 长臂虾科线粒体基因组特征分析及分子标记探讨 [J]. 水产科学, 2011, 30(6): 347-351. doi: 10.3969/j.issn.1003-1111.2011.06.008 Shen X, Sun M A. Analysis of mitochondrial genome characteristics and exploration of molecular markers in Palaemonidae [J]. Fisheries Science, 2011, 30(6): 347-351. doi: 10.3969/j.issn.1003-1111.2011.06.008
[8] 孙成飞, 叶星, 董浚健, 等. 罗氏沼虾6个养殖群体遗传多样性的微卫星分析 [J]. 南方水产科学, 2015, 11(2): 20-26. doi: 10.3969/j.issn.2095-0780.2015.02.003 Sun C F, Ye X, Dong J J, et al. Genetic diversity analysis of six cultured populations of Macrobrachium rosenbergii using microsatellite markers [J]. Southern Marine Science, 2015, 11(2): 20-26. doi: 10.3969/j.issn.2095-0780.2015.02.003
[9] 陈佳毅. 罗氏沼虾亲虾群体遗传结构分析及育苗参数比较 [D]. 扬州: 扬州大学, 2016: 12-13 Chen J Y. Genetic structrue analysis and larvae rearing parameters comparison of Macrobranchium rosenbergii [D]. Yangzhou: Yangzhou University, 2016: 12-13
[10] Chareontawee K, Poompuang S, Na-Nakorn U, et al. Genetic diversity of hatchery stocks of giant freshwater prawn (Macrobrachium rosenbergii) in Thailand [J]. Aquaculture, 2007, 271(1-4): 121-129. doi: 10.1016/j.aquaculture.2007.07.001
[11] Khan S R, Akter H, Sultana N, et al. Genetic diversity in three river populations of the giant freshwater prawn (Macrobrachium rosenbergii) in Bangladesh assessed by microsatellite DNA markers [J]. International Journal of Agriculture & Biology, 2014, 16(1): 195-200.
[12] Thanh H N, Liu Q, Zhao L, et al. Genetic diversity of the cultured giant freshwater prawn (Macrobrachium rosenbergii) in China based on microsatellite markers [J]. Biochemmical Systematics and Ecology, 2015(59): 144-154. doi: 10.1016/j.bse.2014.12.023
[13] 杨学明, 郭亚芬, 陈福艳, 等. 罗氏沼虾3个群体线粒体COI基因的序列差异和遗传标记研究 [J]. 遗传, 2006, 28(5): 540-544. doi: 10.3321/j.issn:0253-9772.2006.05.008 Yang X M, Guo Y F, Chen F Y, et al. Intraspecific DNA sequence polymorphism and genetic markers in the mitochondrial COI gene from three populations of Macrobranchium rosenbergii [J]. Hereditas, 2006, 28(5): 540-544. doi: 10.3321/j.issn:0253-9772.2006.05.008
[14] 姚茜, 杨频, 陈立侨, 等. 罗氏沼虾三群体线粒体D-Loop基因序列差异的初步研究 [J]. 水产学报, 2007, 31(S1): 18-22. Yao Q, Yang P, Chen L Q, et al. Interspecific DNA sequence polymorphism in the mitochondrial D-loop gene from three populations of Macrobranchium rosenbergii [J]. Journal of Fisheries of China, 2007, 31(S1): 18-22.
[15] Thanh H N, Liu Q, Zhao L, et al. Genetic diversity of cultured populations of giant freshwater prawn (Macrobrachium rosenbergii) in China using mtDNA COI and 16S rDNA markers [J]. Biochemical Systematics and Ecology, 2015, 62: 261-269. doi: 10.1016/j.bse.2015.09.011
[16] 吕广祺. 基于微卫星DNA和线粒体COI基因对辽宁地区东北林蛙遗传多样性的研究 [D]. 沈阳: 沈阳农业大学, 2017: 36 Lü G Q. Genetic diversity analysis of the Dybowski’s frog (Rana dybowskii) populations in Liaoning based on microsatellite and COI gene [D]. Shenyang: Shenyang Agricultural University, 2017: 36
[17] 李婉玫. 基于微卫星和线粒体DNA的陕西秦巴山区中华蜜蜂遗传多样性研究 [D]. 西安: 陕西师范大学, 2017: 63-64 Li W M. Genetic diversity analysis of the Apis cerana cerana Fabricius populations in Qinba, Shanxi based on microsatellite and mitochondrial DNA [D]. Xi’an: Shanxi Normal University, 2017: 63-64
[18] 郑燕. 基于微卫星和线粒体标记的梨小食心虫种群遗传多样性和遗传结构研究 [D]. 杨凌: 西北农林科技大学, 2014: 159 Zhen Y. Genetic diversity and genetic structure of oriental fruit moth, Grapholita molesta (L.) analyzed by microsatellite markers and DNA gene sequences [D]. Yangling: Northwest A & F University, 2014: 159
[19] 黄代新, 杨庆恩. 卡方检验和精确检验在HWE检验中应用 [J]. 法医学杂志, 2004(2): 116-119. doi: 10.3969/j.issn.1004-5619.2004.02.019 Huang D X, Yang Q E. Application of chi-square test and exact test in Hardy-Weinberg equilibrium testing [J]. Journal of Forensic Medicine, 2004(2): 116-119. doi: 10.3969/j.issn.1004-5619.2004.02.019
[20] 文亚峰, Kentaro U, 韩文军, 等. 微卫星标记中的无效等位基因 [J]. 生物多样性, 2013, 21(1): 117-126. doi: 10.3724/SP.J.1003.2013.10133 Wen Y F, Kentaro U, Han W J, et al. Null alleles in microsatellite markers [J]. Biodiversity Science, 2013, 21(1): 117-126. doi: 10.3724/SP.J.1003.2013.10133
[21] Shete S, Tiwari H, Elston R C. On estimating the heterozygosity and polymorphism information content value [J]. Theoretical Population Biology, 2000, 57(3): 265-271. doi: 10.1006/tpbi.2000.1452
[22] Botstein D, White R L, Skolnick M, et al. Construction of a genetic linkage map in man using restriction fragment length polymorphism [J]. American Journal of Human Genetics, 1980, 32(3): 314-331.
[23] 杨锐, 喻子牛. 山东沿海褶牡蛎与太平洋牡蛎等位基因酶的遗传变异 [J]. 水产学报, 2000, 24(2): 130-133. Yang R, Yu Z N. Allozyme variation within Crassostrea plicatula and Crassostrea gigas from Shandong coastal waters [J]. Journal of Fisheries of China, 2000, 24(2): 130-133.
[24] Wright S. Evolution and the Genetics of Populations [C]. Invariability Within and Among Natural Populations. Chicago: University of Chicago Press, 1978: 4
[25] 李大命, 张彤晴, 唐晟凯, 等. 基于线粒体 Cytb 基因和 D-loop 区序列的洪泽湖湖鲚遗传多样性分析 [J]. 江苏农业科学, 2018, 46(20): 36-39. Li D M, Zhang T Q, Tang S K, et al. Study on genetic diversity of Coilia nasus population in Lake Hongze based on mtDNA Cytb gene and D-loop sequences [J]. Jiangsu Agricultural Sciences, 2018, 46(20): 36-39.
-
期刊类型引用(3)
1. 刘欣，刘良，刘镇豪，秦强，张绍斌，张富斌. 巴河春季浮游动物群落结构及其与环境因子关系初探——以巴州区段为例. 天津农业科学. 2024(04): 66-73 . 
百度学术
2. 熊雪平，翟东东，刘红艳，陈元元，熊飞. 汉阳水网中5个连通湖泊的原生动物群落结构与物种多样性. 海洋与湖沼. 2024(06): 1466-1477 . 百度学术
3. 张紫薇，陈召莹，张甜娜，周石磊，崔建升，罗晓. 岗南水库沉积物好氧反硝化菌群落时空分布特征. 环境科学. 2022(01): 314-328 . 百度学术
其他类型引用(3)
计量
- 文章访问数: 1870
- HTML全文浏览量: 538
- PDF下载量: 86
- 被引次数: 6