研究了红水河野生鲢Hypophthalmichthys molitrix上游、中游和下游3个群体24个个体的线粒体控制区(D-loop)和细胞色素b (Cyt b)基因序列,通过PCR扩增与测序,分别获得了长度为820 bp和1117 bp的D-loop和Cyt b基因片段。分析表明,D-loop序列共定义了20个单倍型,存在74个多态位点,发生转换19次、颠换3次、插入/缺失52次。Cyt b序列共定义了9个单倍型,存在48个多态性位点,发生转换13次、颠换1次、插入/缺失34次。两序列分析结果表明,群体间和群体内遗传距离较小,遗传变异的分子变异等级不显著。上、中、下游均有共享单倍型,3个群体间遗传分化不大。 Sequences of mitochondrial D-loop and Cyt b gene of 24 individuals of Hypophthalmichthys molitrix which were sampled from three populations in Hongshui River were determined. Using PCR amplification and sequencing, 820 bp of D-loop partial sequences and 1117 bp of Cyt b partial se-quences were obtained. In D-loop sequences, a total of 74 polymorphic sites were defined; 20 dis-tinct haplotypes, 19 transitions, 3 transversions and 52 insert/lacuna loci were found. In Cyt b se-quences, 48 polymorphic sites and 9 haplotypes were defined; 13 transitions, 1 transversion and 34 insert/lacuna loci were found. Two-sequence analysis showed that, within groups and between groups, the genetic distance is small, and the genetic variation of molecular variation level is not significant. Upper, middle and lower reaches are shared haplotypes, and the genetic differentiation among three populations is little.
吴伟军,彭敏,王大鹏,施军,李育森,韩耀全,雷建军,何安尤*
广西水产科学研究院/广西水产遗传育种与健康养殖重点实验室,广西 南宁
收稿日期:2015年12月5日;录用日期:2015年12月18日;发布日期:2015年12月24日
研究了红水河野生鲢Hypophthalmichthys molitrix上游、中游和下游3个群体24个个体的线粒体控制区(D-loop)和细胞色素b (Cyt b)基因序列,通过PCR扩增与测序,分别获得了长度为820 bp和1117 bp的D-loop和Cyt b基因片段。分析表明,D-loop序列共定义了20个单倍型,存在74个多态位点,发生转换19次、颠换3次、插入/缺失52次。Cyt b序列共定义了9个单倍型,存在48个多态性位点,发生转换13次、颠换1次、插入/缺失34次。两序列分析结果表明,群体间和群体内遗传距离较小,遗传变异的分子变异等级不显著。上、中、下游均有共享单倍型,3个群体间遗传分化不大。
关键词 :鲢,红水河,D-Loop,Cyt b,遗传变异
鲢(Hypophthalmichthys molitrix)属鲤形目,鲤科,鲢亚科,鲢属。又名白鲢、跳鲢、鲢子。体侧扁、稍高,呈纺锤形,头大吻短,口宽、眼小,鳞细,体呈银白色,腹部有肉棱。鲢滤食浮游植物,生长速度快、产量高,在广西一般做为生态养殖品种。鲢属中上层鱼。春夏秋三季,绝大多数时间在水域的中上层游动觅食,冬季则潜至深水越冬。秋除浮游生物外,还大量地吃腐屑类饵料。其食欲与水温成正比,夏季水位越低,其摄食量越大:冬季则少吃少动。最适宜的水温为23℃~32℃。鲢性情活泼,喜欢跳跃,有逆流而上的习性。
线粒体D-环区(displacement loop region, D-loop),是线粒体DNA中的一段非编码区,位于tRNA Pro和tRNA Phe基因之间,是线粒体DNA序列和长度变异最大的区域,也是线粒体基因中进化最快的部分,适用于种群水平差异的检测,也可用于种间分析[
采集3个鲢群体共27条,分别为上游岩滩水库群体(shang)、中游百龙滩水库群体(zhong)和下游来宾市区江段群体(xia),三个群体采样区域有大坝分隔,大坝无过鱼设施,可认为群体间无遗传交流。每尾个体取0.1 g背部肌肉新鲜样品,于无水乙醇中保存。
通过液氮研磨、裂解液裂解、有机溶剂抽提,使进入水相的蛋白质与核酸分层,在RNA酶的作用下,降解RNA,得纯度较高的DNA样品并进行PCR测序。D-loop引物序列为:DL1 (5′-ACCCCTGGCTCCCAA AGC-3′),DH2 (5′-ATCTTAGCATCTTCAGTG-3′) [
使用SeqMan II软件拼接获得序列并人工校对,用CLUSTAL-X程序比对,DAMBE统计单倍型。用ARLEQUIN统计软件计算获得单倍型多样性(Haplotype diversity, H)和核苷酸多样性(Nucleotide diversity, π),群体间的分化指数(F-statistics, Fst) [
由表1可见,在27个820 bp的D-loop基因序列中共检测出20个单倍型。检测到4个单突变位点,70个简约信息位点数;在74个突变位点中:转换位点19个,颠换位点3个,插入/缺失位点52个,转换高于颠换。基因片段中,A + T的含量(67.3%)明显高于C + G的含量(32.7%)。29个1117 bp的Cyt b基因序列中共检测出9个单倍型,检测到1个单突变位点,47个简约信息位点;在48个突变位点中:转换位点13个,颠换位点1个,插入/缺失位点34个,A + T的含量(58.3%)高于C + G的含量(41.6%)。对比分析可知:Cyt b单突变位点、多态位点百分比和简约信息位点百分比分别为4.30%、4.21%和0.09%,均小于D-loop,后者分别为9.02%,8.54%和0.48%。
由表2可见,D-loop序列中,红水河3个鲢群体内的遗传距离分别为0.0322、0.0246和0.0150。上游群体与中游群体间遗传距离为0.0284,上游群体与下游群体间遗传距离为0.0246,中游群体和下游群体间遗传距离为0.0199。Cytb基因序列中,群体内的遗传距离分别为0.0164、0.0111、0.0036。上游群体与中游群体间遗传距离为0.0137,上游群体与下游群体间遗传距离为0.0107,中游群体和下游群体间遗传距离为0.0074。两序列的群体内遗传距离均为上游高于中游,中游又高于下游。
3个群体的D-loop和Cyt b序列的单倍型多样性H、平均核苷酸变异数K、和核苷酸多样性π见表3,下游群体的D-loop和Cyt b序列单倍型多样性高于上游和中游,可能与其采集到的样品较少,检测到的单倍型数也较少有关,可以作为例证的是,其平均核苷酸变异数K和核苷酸多样性π低于上游和中游。对比D-loop和Cyt b序列,D-loop的平均核苷酸变异数和核苷酸多样性在上、中下游均高于Cyt b序列。
基因序列 | 单倍型 | 比对碱基位点数 | 多态位点数(比例/%) | 简约信息位点数(比例/%) | 单突变位点数(比例/%) | 转换、颠换数 | A + T/% | C + G/% |
---|---|---|---|---|---|---|---|---|
D-loop | 20 | 820 | 74 (9.02%) | 70 (8.54%) | 4 (0.48%) | 19/3 | 67.3 | 32.7 |
Cyt b | 9 | 1117 | 48 (4.30%) | 47 (4.21%) | 1 (0.09%) | 13/1 | 58.3 | 41.6 |
表1. 鲢线粒体D-loop及Cyt b碱基组成及变异分析
遗传距离 | 上游(Shang) | 中游(Zhong) | 下游(Xia) |
---|---|---|---|
上游 (Shang) | 0.0322 | ||
0.0164 | |||
中游 (Zhong) | 0.0284 | 0.0246 | |
0.0137 | 0.0111 | ||
下游 (Xia) | 0.0246 | 0.0199 | 0.0150 |
0.0107 | 0.0074 | 0.0036 |
表2. 红水河鲢D-loop和Cyt b基因序列群体内和群体间遗传距离
注:上面的代表D-loop,下面的代表Cyt b (下同)。
群体 | 样本数 | 多态位点 | 单倍型数 | 单倍型多样性H | 平均核苷酸变异数K | 核苷酸多样性π |
---|---|---|---|---|---|---|
上游 (Shang) | 14 | 72 | 13 | 0.9890 | 24.9670 | 0.0309 |
48 | 7 | 0.8791 | 17.7252 | 0.0159 | ||
中游 (Zhong) | 10 | 66 | 8 | 0.9556 | 19.2889 | 0.0238 |
46 | 5 | 0.8444 | 12.1333 | 0.0109 | ||
下游 (Xia) | 3 | 18 | 3 | 1.0000 | 12.0000 | 0.0148 |
6 | 3 | 1.0000 | 4.0000 | 0.0036 |
表3. 3个鲢群体相关遗传多样性指标
表4和表5分别列出群体间的Fst、Nm值和群体内群体间的遗传变异的分子变异等级分析(AMOVA),从分析结果可见,3个群体间D-loop和Cyt b序列遗传差异均不显著(Fst = −0.01379或−0.03160,p > 0.05)。
以鲢D-loop和Cyt b基因序列为外群,用NJ法构建分子系统树,其拓扑结构见图1。可见,3个群体的个体混杂在一起,群体间未表现出聚类关系。
单倍型多样性(H)、平均核苷酸差异数(K)和核苷酸多样性(π)是反映物种遗传变异水平的3个重要参数,3个群体遗传多样性差别不大,群体内遗传距离表现为上游高于中游,中游高于下游。遗传变异的分子变异等级分析结果为不显著。24个样本的D-loop序列分属20个单倍型,而Cyt b序列分属9个单倍型。上、中、下游均有共享单倍型,NJ系统树也显示3个群体并未聚集为独立的分支。由于线粒体DNA为母系遗传,一般不发生重组,表明3个群体间遗传分化不大。
上游(Shang) | 中游(Zhong) | 下游(Xia) | |
---|---|---|---|
上游(Shang) | 0.0322 | ||
0.0164 | |||
中游(Zhong) | 0.0284 | 0.0246 | |
0.0137 | 0.0111 | ||
下游(Xia) | 0.0246 | 0.0199 | 0.0150 |
0.0107 | 0.0074 | 0.0036 |
表4. 3个鲢群体间的Fst和Nm
变异来源 | 自由度d.f | 平方和 | 方差组分 | 变异百分率 | F统计量 | P值 |
---|---|---|---|---|---|---|
群体间 | 2 | 22.140 | −0.16859 Va | −1.38 | −0.01379 | 0.46823 |
2 | 11.001 | −0.22186 Va | −3.16 | −0.03160 | 0.50635 | |
群体内 | 24 | 297.452 | 12.39385 Vb | 101.38 | ||
24 | 173.814 | 7.24226 Vb | 103.16 | |||
总体 | 26 | 319.593 | 12.22526 | |||
26 | 184.815 | 7.02040 |
表5. 3个鲢群体间的遗传变异的分子变异等级分析(AMOVA)
图1. 红水河3个鲢群体D-loop (左)和Cyt b (右)基因序列的NJ分子系统树
红水河综合利用规划实施后,红水河流域形成了9个首尾衔接的梯级水库,流水生境绝大部分已经萎缩,天然河道和库区流水江段大为缩短,多数流水性鱼类,特别是产漂流性卵的鱼类,在库区已不具备形成自然种群条件 [
对红水河野生鲢上游、中游和下游3个群体的24个个体的D-loop和Cyt b基因序列分析表明,D-loop序列的单倍型数、多态位点数均高于Cyt b序列。红水河鲢群体间和群体内遗传距离较小,遗传变异的分子变异等级不显著。上、中、下游均有共享单倍型,3个群体间遗传分化不大。可能是原本连通的群体分隔后尚未产生显著遗传分化,也可能是来自近缘亲本的增殖放流群体已经取代了野生群体。
社会公益研究项目(CXIF-2014-006);自治区直属公益性科研院所基本科研业务专项资金项目(编号:GXIF-2014-006)。
吴伟军,彭 敏,王大鹏,施 军,李育森,韩耀全,雷建军,雷建军,何安尤. 红水河鲢线粒体D-Loop和Cyt b基因序列分析 Comparison of Mitochondrial D-Loop and Cyt b Sequences of Hypophthalmichthys molitrix[J]. 水产研究, 2015, 02(04): 67-73. http://dx.doi.org/10.12677/OJFR.2015.24008