中国海洋湖沼学会主办。
文章信息
- 刘良国, 陈海康, 许倩, 杨品红, 谢春华, 李梦军. 2016.
- LIU Liang-Guo, CHEN Hai-Kang, XU Qian, YANG Pin-Hong, XIE Chun-Hua, LI Meng-Jun. 2016.
- 基于微卫星标记的洞庭青鲫与三个鲫品系群体遗传多样性分析
- GENETIC ANALYSIS OF Dongtingking crucian carp and other THREE strains of crucian carp BASED ON MICROSATELLITE MARKERS
- 海洋与湖沼, 47(3): 647-654
- Oceanologia et Limnologia Sinica, 47(3): 647-654.
- http://dx.doi.org/10.11693/hyhz20160300059
-
文章历史
- 收稿日期:2016-03-17
- 改回日期:2016-04-07
洞庭青鲫(Dongtingking crucian carp,Carassius auratus var. Dongtingking)是近年来本文作者在洞庭湖水系澧水北民湖水域发现的一个具有两性繁殖方式的天然二倍体鲫鱼新品系(杨品红等,2005,2007; 吴珊等,2006)。由于其在经济性状上具有生长快、产量高、抗逆性强、味道鲜美等特点,受到养殖户和广大消费者喜爱,发展前景广阔。在生物学特性上,与国内已报道的其它鲫品系性状相比,洞庭青鲫与彭泽鲫和普通野鲫尤为相似: 彭泽鲫原产于江西省彭泽县一带,最初的研究结果表明,彭泽鲫也是一种二倍体两性生殖鱼类(傅永进,1996),但随后报道的彭泽鲫均为染色体数目在150条左右的雌核发育鱼(Zhou et al,2002; 刘良国等,2004),在染色体倍性和生殖方式上,洞庭青鲫均与最初报道的彭泽鲫相似; 同时,我们采用核型分析方法对洞庭湖水系普通野鲫群体进行批量的染色体数目统计,发现在洞庭湖水系中,同时存在两种在形态上极为相似、染色体数为100和150 ±的不同倍性野鲫群体,而对洞庭青鲫染色体组型分析(杨品红等,2007)表明它与洞庭湖水系野生二倍体鲫倍性一致。洞庭青鲫作为一个新型的鲫鱼品系,它与彭泽鲫、二倍体和三倍体野鲫到底存在怎样的遗传关系,这将是本论文重点关注的问题。
微卫星(Simple Sequence Repeat,SSR,单一的序列重复)分子标记在真核生物基因组中广泛存在,因其具有多态性高、共显性遗传、可重复性强和高度特异性等优点,现已广泛应用于群体遗传多样性(周盼等,2011; 宋春妮等,2011)、基因连锁与遗传图谱构建(Miyao et al,1996)、品种鉴定、亲本分析、纯系检验等(Stockburge et al,2000; 刘静霞等,2003; Bessert et al,2003; 文萍等,2015)研究领域。为了了解洞庭青鲫与彭泽鲫、二倍体和三倍体野鲫的分子遗传多样性,揭示它们之间的遗传关系,本实验选用12个微卫星标记对洞庭青鲫与彭泽鲫、二倍体和三倍体野鲫的遗传多样性进行检测分析,实验结果对于鲫的多倍体起源演化研究,以及二倍体鲫的种质资源保护与选育种实践具有重要意义。
1 材料与方法 1.1 材料洞庭青鲫(F6)和彭泽鲫取自大湖股份湖南省水产工程技术研究中心实验基地,野生二倍体和三倍体鲫采自洞庭湖水系的沅水和澧水。实验鱼采回后,首先通过活体染色体制片的方法,进一步批量确定不同鲫品系的染色体倍性,然后,每个鲫品系群体各取样30尾,分别剪其背部肌肉用无水乙醇保存备用。
1.2 基因组DNA提取每尾鲫样品分别取1.0g肌肉用于基因组DNA提取,提取方法参照常规的酚-氯仿抽提程序进行。DNA提取后,用0.8%琼脂糖凝胶电泳进行质量检测,同时用紫外分光光度计估算其浓度,调整浓度至40ng/μL,4°C冰箱保存备用。
1.3 微卫星引物本实验选用在银鲫(周莉等,2001; 孙效文等,2009)和红鲫(李建中等,2005; 鲁双庆等,2005)中均能扩增出清晰的DNA谱带,且重复稳定性好的12对微卫星引物序列(表 1),由北京鼎国昌盛生物技术有限责任公司合成为蓝色(FAM)荧光标记引物。
引物名称 | 引物序列(5′→3′) | 退火温度(°C) | 等位基因大小(bp) |
MFW4 | F: TCCAAGTCAGTTTAATCACCG R: GGGAAGCGTTGACAACAAGC | 54 | 160-200 |
MFW16 | F: GTCCATTGTGTCAAGATAGAC R: TCTTCATTTCAGGCT GCAAAG | 52 | 350-366 |
MFW18 | F: GTCCCTGGTAGTGAGTGAGT R: GCGTTGACTTGTTTTATAGTAG | 56 | 272-286 |
MFW19 | F: GAATCCTCCATCATGCAAAC R: GCACAAACTCCACATTGTGCC | 58 | 186-224 |
MFW20 | F: CAGTGAGACGATTACCTTGG R: GTGAGCAGCCCACATTGAAG | 55 | 145-216 |
MFW23 | F: GTATAATTGGGAGTTTTAGGG R: CAGGTTTATCTCCCTTCTAG | 57 | 100-120 |
SCM4 | F: TTCACTAGAACAGCAGGAGT R: ACAATGGTAACAGGGACAC | 52 | 98-138 |
SCM8 | F: GCTGTTTGCTCAGGGAC R: TGGAGGAGGGTTAGGGT | 51 | 246-258 |
SCM10 | F: GAGCCGCTCTTACCTTGATG R: GTCTCATCCGCTCCAAACTG | 60 | 142-160 |
SCM12 | F: AAGCACCGGAGATCATTCTG R: GGGGCCTCTCTGTTATTCCT | 60 | 238-250 |
SCM13 | F: ACCGTTTTGAGAGCGAGTGT R: AACCAGGTGTGACACAATGC | 60 | 160-172 |
SCM17 | F:TGAATCAGTCACCAAGTCAGAG R: CGACAGACGCATGAGGTAGA | 59 | 108-146 |
F. 正向引物; R. 反向引物 |
PCR反应体系为25μL,其中40ng/μL DNA模板2.0μL,10mmol/L引物各0.5μL,10mmol/L dNTPs 0.5μL,10×PCR buffer 2.5μL,Taq DNA polymease 0.5μL(2U/μL),ddH2O 18.5μL。PCR扩增程序为: 95°C 5min,接着94°C 30s,退火30s(退火温度见表 1),72°C 30s,共35个循环,最后72°C 10min。将带有蓝色(FAM)荧光标记的PCR产物与红色分子量内标ROX 500混合,采用ABI3730X Genetic Analyzer对其进行毛细管电泳检测,得到扩增片段大小、峰图等原始实验数据; 同时,将扩增产物在4%的变性聚丙烯酰胺凝胶上电泳(50W恒功率2.5h),ABI PRISM 377测序仪扫描,得到电泳图。
1.5 数据统计与分析对毛细管电泳测序仪检测得到的原始扩增条带进行统计,将每一条带视为该位点的1个等位基因,利用Genepop V4.0和EXCEL MS TK-2007软件进行分析,计算并统计12个微卫星基因座位的等位基因数和等位基因频率、观测杂合度(observed heterozygosity,Ho)、期望杂合度(expected heterozygosity,He)和多态信息含量(polymorphism information content,PIC)。
通过Lynch(1990)的方法计算群体间的遗传相似性系数和遗传距离,根据遗传距离,再以MEGA 4.0软件包中的非加权配对算术平均法(unweighted pair group methods using arithmetic average,UPGMA)对四个鲫品系群体进行聚类分析。
2 结果与分析 2.1 微卫星PCR扩增结果及等位基因频率12对微卫星引物在4个鲫品系群体中均得到较好的扩增结果,图 1为引物MFW4和MFW20在4个群体中的部分变性聚丙烯酰胺凝胶电泳图谱。在12个基因座位中,每个座位检测到的等位基因数4—12个,平均等位基因数6.58,共检测出78个等位基因,其中洞庭青鲫群体有27个,二倍体野鲫群体49个,三倍体野鲫群体54个,彭泽鲫群体34个,78个等位基因中有14个等位基因为4个群体所共有。与二倍体和三倍体野鲫群体相比,洞庭青鲫、彭泽鲫等养殖群体的等位基因数明显偏少,基因座位MFW4的1,MFW19的2,MFW23的3,SCM4的5、7,SCM8的1,SCM10的3、5、7,SCM13的5,SCM17的1、2、8等一些等位基因在两个野生鲫群体中均有出现,而在洞庭青鲫和彭泽鲫养殖群体中缺失。12个基因座位在4个鲫品系群体中的等位基因频率见表 2。
座 位 | 等位基因 | 群 体 | ||||
DCC | 2nCC | 3nCC | PCC | |||
MFW4 | 1 | 0.200 | 0.100 | |||
2 | 0.500 | 0.100 | 0.500 | 0.500 | ||
3 | 0.200 | 0.300 | 0.300 | |||
4 | 0.500 | 0.200 | 0.200 | |||
5 | 0.100 | |||||
6 | 0.200 | |||||
7 | 0.100 | |||||
MFW16 | 1 | 0.600 | 0.200 | 0.300 | ||
2 | 0.200 | |||||
3 | 0.200 | 1.000 | 0.800 | 0.600 | ||
4 | 0.100 | |||||
MFW18 | 1 | 0.100 | ||||
2 | 0.400 | 0.100 | 0.100 | |||
3 | 0.100 | |||||
4 | 0.600 | 0.400 | 0.300 | 0.500 | ||
5 | 0.400 | 0.200 | 0.400 | 0.400 | ||
MFW19 | 1 | 0.200 | 0.100 | 0.300 | ||
2 | 0.200 | 0.300 | ||||
3 | 0.300 | |||||
4 | 0.200 | |||||
5 | 1.000 | 0.200 | 0.100 | 0.700 | ||
6 | 0.100 | |||||
7 | 0.200 | |||||
8 | 0.100 | |||||
MFW20 | 1 | 0.200 | 0.400 | 0.300 | 0.100 | |
2 | 0.200 | 0.100 | 0.300 | |||
3 | 0.200 | 0.100 | ||||
4 | 0.100 | |||||
5 | 0.200 | 0.100 | 0.300 | |||
6 | 0.200 | 0.300 | 0.200 | 0.300 | ||
7 | 0.100 | 0.300 | ||||
MFW23 | 1 | 0.500 | 0.100 | 0.200 | ||
2 | 0.500 | 0.800 | 0.600 | 0.500 | ||
3 | 0.200 | 0.200 | ||||
4 | 0.300 | |||||
5 | 0.100 | |||||
SCM8 | 1 | 0.200 | 0.200 | |||
2 | 0.200 | |||||
3 | 0.500 | 0.400 | 0.300 | 0.200 | ||
4 | 0.200 | |||||
5 | 0.500 | 0.500 | 0.800 | |||
SCM4 | 1 | 0.500 | 0.300 | 0.400 | 0.800 | |
2 | 0.100 | 0.100 | ||||
3 | 0.200 | |||||
4 | 0.100 | |||||
5 | 0.200 | 0.100 | ||||
6 | 0.100 | |||||
7 | 0.200 | 0.200 | ||||
8 | 0.400 | 0.100 | 0.100 | 0.100 | ||
SCM10 | 1 | 0.100 | ||||
2 | 0.500 | 0.200 | 0.100 | 0.200 | ||
3 | 0.400 | 0.300 | ||||
4 | 0.300 | 0.300 | ||||
5 | 0.300 | 0.100 | ||||
6 | 0.500 | 0.500 | ||||
7 | 0.100 | 0.100 | ||||
SCM12 | 1 | 0.100 | ||||
2 | 0.300 | 0.300 | 0.300 | |||
3 | 0.500 | 0.300 | 0.500 | 0.200 | ||
4 | 0.400 | |||||
5 | 0.500 | 0.100 | 0.500 | |||
SCM13 | 1 | 1.000 | 0.200 | 0.600 | ||
2 | 0.100 | |||||
3 | 0.200 | 0.300 | 0.400 | |||
4 | 0.400 | |||||
5 | 0.300 | 0.200 | ||||
6 | 0.300 | |||||
SCM17 | 1 | 0.200 | 0.100 | |||
2 | 0.200 | 0.100 | ||||
3 | 0.100 | |||||
4 | 0.100 | |||||
5 | 0.500 | 0.100 | 0.100 | 0.500 | ||
6 | 0.200 | |||||
7 | 0.200 | |||||
8 | 0.200 | 0.100 | ||||
9 | 0.100 | |||||
10 | 0.200 | |||||
11 | 0.500 | 0.200 | ||||
12 | 0.300 |
12个微卫星位点在同一鲫品系群体,以及同一个微卫星位点在不同鲫品系群体的遗传变异程度均存在较大差异(表 3)。洞庭青鲫、野生二倍体和三倍体鲫、彭泽鲫群体的平均观测杂合度分别为0.633、0.750、0.800、0.717,平均期望杂合度分别为0.502、0.713、0.757、0.602,平均多态信息含量分别为0.364、0.599、0.637、0.470,其中野生二倍体和三倍体鲫群体的平均观测杂合度、期望杂合度和多态信息含量均高于洞庭青鲫和彭泽鲫养殖群体,表现为高度多态,而洞庭青鲫和彭泽鲫养殖群体表现为中度多态,部分座位如洞庭青鲫的MFW19和SCM13两个位点PIC值为0,遗传多态性消失。
群体 | 参数 | 微卫星座位 | 平均值 | |||||||||||
MFW4 | MFW16 | MFW18 | MFW19 | MFW20 | MFW23 | SCM4 | SCM8 | SCM10 | SCM12 | SCM13 | SCM17 | |||
DCC | Ho | 1.000 | 0.400 | 0.400 | 0.000 | 0.800 | 1.000 | 1.000 | 1.000 | 1.000 | 1.000 | 0.000 | 1.000 | 0.750±0.056 |
He | 0.556 | 0.622 | 0.533 | 0.000 | 0.889 | 0.556 | 0.644 | 0.556 | 0.556 | 0.556 | 0.000 | 0.556 | 0.502±0.073 | |
PIC | 0.375 | 0.499 | 0.365 | 0.000 | 0.768 | 0.375 | 0.492 | 0.375 | 0.375 | 0.375 | 0.000 | 0.375 | 0.364±0.023 | |
2nCC | Ho | 1.000 | 0.000 | 0.800 | 0.800 | 1.000 | 0.400 | 1.000 | 0.400 | 1.000 | 0.600 | 0.600 | 1.000 | 0.717±0.058 |
He | 0.911 | 0.000 | 0.711 | 0.889 | 0.800 | 0.356 | 0.889 | 0.800 | 0.778 | 0.733 | 0.778 | 0.911 | 0.713±0.078 | |
PIC | 0.794 | 0.000 | 0.563 | 0.768 | 0.676 | 0.269 | 0.772 | 0.672 | 0.645 | 0.586 | 0.645 | 0.794 | 0.599±0.053 | |
3nCC | Ho | 0.800 | 0.000 | 0.800 | 0.600 | 1.000 | 0.400 | 0.200 | 0.200 | 0.800 | 0.800 | 1.000 | 1.000 | 0.633±0.062 |
He | 0.711 | 0.356 | 0.800 | 0.867 | 0.844 | 0.644 | 0.822 | 0.689 | 0.867 | 0.711 | 0.822 | 0.956 | 0.757±0.045 | |
PIC | 0.581 | 0.269 | 0.676 | 0.748 | 0.720 | 0.535 | 0.701 | 0.548 | 0.748 | 0.581 | 0.692 | 0.844 | 0.637±0.033 | |
PCC | Ho | 1.000 | 0.400 | 1.000 | 0.600 | 1.000 | 1.000 | 0.400 | 0.400 | 1.000 | 1.000 | 0.800 | 1.000 | 0.800±0.052 |
He | 0.689 | 0.600 | 0.644 | 0.467 | 0.800 | 0.689 | 0.378 | 0.356 | 0.689 | 0.689 | 0.533 | 0.689 | 0.602±0.040 | |
PIC | 0.548 | 0.466 | 0.492 | 0.332 | 0.666 | 0.548 | 0.314 | 0.269 | 0.548 | 0.548 | 0.365 | 0.548 | 0.470±0.027 |
根据等位基因频率计算得到不同鲫品系群体之间的遗传相似性系数和遗传距离(表 4),计算结果表明,4个鲫品系群体间的平均遗传相似系数在0.421— 0.787,平均遗传距离为0.213—0.579。其中洞庭青鲫与彭泽鲫群体之间、野生二倍体和三倍体鲫群体之间的遗传相似性系数较高,分别为0.787、0.621,遗传距离较小,分别为0.213、0.379; 而洞庭青鲫、彭泽鲫与野生二倍体鲫的遗传相似性系数较低,分别为0.421、0.482,遗传距离较大,分别为0.579、0.518。
群体 | DCC | 2nCC | 3nCC | PCC |
DCC | 0.421 | 0.543 | 0.787 | |
2nCC | 0.579 | 0.621 | 0.482 | |
3nCC | 0.457 | 0.379 | 0.591 | |
PCC | 0.213 | 0.518 | 0.409 |
根据群体间的遗传距离矩阵(表 4),采用MEGA 4.0软件包中的UPGMA聚类法分析4个鲫品系群体间的亲缘关系(图 2),图 2显示,4个鲫品系群体明显分为两个大支,其中,洞庭青鲫与彭泽鲫聚为一支,野生二倍体与三倍体鲫聚为另一支,表明洞庭青鲫与彭泽鲫养殖群体的亲缘关系较近,野生二倍体与三倍体鲫的亲缘关系较近。
3 讨论遗传多样性是指地球上生物所携带的各种遗传信息的总和,它包括生物遗传变异大小和遗传变异分布格局(即群体遗传结构)等内容。它是物种适应多变的环境条件、维持生存和进化的基础,物种的遗传多样性越丰富,其对环境的生存适应能力和进化潜力就越强(O’Connell et al,1997; 张文静等,2003; 朱晓东等,2007)。由于微卫星标记遵循孟德尔定律呈共显性遗传,在检测种群异质性方面有明显优势(Bártfai et al,2003),因而能更多的揭示群体的遗传变异水平。本研究利用12对微卫星标记对洞庭青鲫等4个鲫品系群体进行了比较分析,12个微卫星位点共检测到78个等位基因,其中有14个等位基因为4个鲫品系群体所共有,这些共有等位基因频率相对较高,推测是鲫进化过程中最原始、最保守的一部分,对于维持鲫物种的繁衍和稳定有着重要作用。其余等位基因,如在野生鲫群体中存在而在养殖鲫群体中丢失,或只在洞庭青鲫和彭泽鲫养殖鲫群体中出现而野生鲫群体中未见,或只在某一个鲫品系群体中出现等情况(表 2),可能是进化过程中由于基因的插入、缺失等突变机制所造成。
群体平均基因杂合度表示被检测位点上群体中杂合子的频率,它是衡量群体杂合程度的指标(周盼等,2011)。群体平均基因杂合度越低,反映该群体的遗传一致性越高,即群体的遗传变异越少,遗传多样性越低(Leberg,2002; 张燕萍等,2013)。张志伟等(2006)报道了江苏境内草鱼1个野生群体和2个养殖群体的平均观测杂合度为0.5286—0.7000,平均期望杂合度为0.5235—0.6939; 周盼等(2011)报道了三江水系4个野生草鱼群体的平均观测杂合度为0.7114— 0.8045,平均期望杂合度为0.6245—0.6739。本研究结果表明4个鲫品系群体的平均观测杂合度(Ho)介于0.633—0.800,平均期望杂合度(He)介于0.502— 0.757,与以上两个报道的结果相近,说明本实验选用的微卫星标记适合鲫品系群体的遗传多样性研究。从表 3可见,4个鲫品系群体的平均观测杂合度与期望杂合度分析结果较为一致,其中,三倍体野鲫的平均观测杂合度最高,为0.800,二倍体野鲫次之,为0.750,这说明洞庭湖水系野鲫群体的遗传多态性丰富,选育种潜力较大; 洞庭青鲫群体的平均观测杂合度最低,为0.633,说明洞庭青鲫养殖群体的遗传多样性水平较低,这与洞庭青鲫这一优良鲫鱼品系的多代的选育纯化有关,这一结果也与吴珊等(2006)对洞庭青鲫形态性状的分析结果一致: 即洞庭青鲫子代与亲代不仅没有性状分离现象,而且基因有进一步纯化的趋势。
多态信息含量(PIC)是群体内遗传变异的量度,可以作为微卫星位点变异程度指标(周盼等,2011)。依据Botstein等(1980)的划分标准: 当PIC>0.5时,该位点为高度多态位点; 当0.25<PIC<0.5时,为中度多态位点; 当PIC<0.25时,为低度多态位点。从本研究12个微卫星位点来看(表 3),高度多态位点主要集中在二倍体和三倍体野鲫群体,二倍体野鲫平均PIC为0.599,高度多态位点占总位点的83.3%,三倍体野鲫平均PIC为0.637,高度多态位点占总位点的91.7%,说明洞庭湖水系二倍体和三倍体野鲫群体的遗传多样性丰富; 而中、低度多态位点主要集中于洞庭青鲫和彭泽鲫养殖群体,洞庭青鲫平均PIC为0.364,中度多态位点占总位点的75.0%,MFW19和SCM13两个位点的PIC值为0,彭泽鲫平均PIC为0.470,中度和高度多态位点占总位点的比例各占50.0%,以上情况说明,与基因杂合度分析的结果一致,相比于野生鲫群体,洞庭青鲫和彭泽鲫养殖群体的遗传多样性降低,尤以洞庭青鲫遗传多样性下降为甚。
群体间的遗传距离是遗传变异的尺度,遗传距离的大小预示着亲缘关系的远近。本研究用微卫星实验数据分析得到的遗传距离和树状图与先前ISSR的分析结果(刘良国等,2015)基本一致,即洞庭青鲫与彭泽鲫的遗传距离最小,树状图上聚为一支,二者亲缘关系较近,主要原因可能是洞庭青鲫与彭泽鲫均起源于天然水域,因其都具有优良的生产性状而被多代选育,在相似的人工养殖环境下,选择压力相近所致; 二倍体和三倍体野鲫在树状图上聚为另一支,表明二者亲缘关系也较近,其原因可能要追溯到三倍体野鲫的起源上,洞庭湖水系二倍体和三倍体野鲫生活于同一水体,二者在形态上极为相似(刘良国等,2012),在极端环境(如温度的剧变等)条件下,二倍体鲫有可能产生二倍体卵子,这些二倍体卵子再与正常的单倍体精子杂交形成三倍体,因而在基因座位上,三倍体野鲫表现出与二倍体野鲫的相似性。
据沈俊宝等(1983)和王蕊芳等(1988)报道,由于环境急剧变化和天然杂交等因素,自然水体中的二倍体鲫可以演化成三倍体,而三倍体鲫具有雌核发育生殖特性,繁殖力和适应性强,种群优势明显,这可能是当前天然水域中二倍体鲫数量远少于三倍体鲫的重要原因。洞庭青鲫是近年来发现并选育出的生长速度快、营养价值高、遗传性状稳定的二倍体鲫品系,是自然界中十分珍贵的鲫鱼种质资源,具有很高的经济价值和科研价值。由于洞庭青鲫选育的时间较短,至今才经过6代,将来能否保持其二倍体稳定的遗传性状还未为可知。根据上述三倍体鲫的遗传起源推测,二倍体洞庭青鲫亦有可能在遭遇极端环境条件和杂交因素影响下,部分个体特化为行雌核发育方式的三倍体,因三倍体鲫的雌核发育生殖方式比二倍体鲫两性生殖方式具有明显的繁殖优势,在人工选育的条件下,新的洞庭青鲫三倍体种群将会越来越多,而二倍体洞庭青鲫种群将会逐渐消失。因此,加强对洞庭湖水系二倍体野鲫和洞庭青鲫种质资源的保护已经越来越迫切,笔者认为,可以从以下几方面来加强保护:(1)建立洞庭湖水系二倍体鲫种质资源保护区,如在洞庭青鲫的原产地澧水北民湖等水域建立洞庭青鲫省级及以上水产种质资源保护区,加强二倍体鲫种质资源管理;(2)设置洞庭湖水系禁渔期,加强渔业资源管理,杜绝电鱼、毒鱼等毁灭性的捕捞;(3)通过增殖放流进行洞庭青鲫和野生二倍体鲫种群的恢复;(4)加强对三倍体雌核发育鲫如银鲫、彭泽鲫的人工饲养管理,以防三倍体鲫逃逸到天然水体,导致二倍体鲫种群的减少或消失。