留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

西江鲫的遗传多样性与群体结构

匡天旭 帅方敏 陈蔚涛 李新辉

引用本文:
Citation:

西江鲫的遗传多样性与群体结构

    作者简介: 匡天旭(1991 — ),男,硕士研究生,从事鱼类分子生物学研究。E-mail: 1462019455@qq.com;
    通讯作者: 李新辉, lxhui01@aliyun.com
  • 中图分类号: S 932.4

Genetic diversity and population structure of Carassius auratus in Xijiang River

    Corresponding author: Xinhui LI, lxhui01@aliyun.com
  • CLC number: S 932.4

  • 摘要: 为了解珠江鲫(Carassius auratus)遗传种质资源现状,文章以西江为例,利用线粒体细胞色素c氧化酶亚基Ⅰ (COⅠ)基因研究了西江8个地理群体的遗传多样性和群体结构。结果表明,西江鲫群体的遗传多样性处于较高水平。鲫群体存在显著遗传分化(FST=0.164,P<0.001),地理距离和自身生活习性可能是导致鲫群体产生遗传分化的重要因素。此外,中性检验和错配分布分析表明西江鲫没有经历近期的群体扩张。
  • 图 1  鲫样本采集地点

    Figure 1.  Sampling site of C.auratus

    图 2  基于COⅠ基因构建的单倍型NJ系统发育树

    Figure 2.  Neighbor-Joining tree based on COⅠ gene

    图 3  基于COⅠ基因构建的鲫单倍型网络图

    Figure 3.  Median-joining haplotype network of C.auratusbased on COⅠ gene

    图 4  地理距离与两两群体间遗传分化的关系

    Figure 4.  Relationship between geographic distance and pairwise differentiation

    图 5  西江鲫群体核苷酸错配分布分析

    Figure 5.  Mismatch distribution analysis of C.auratus in Xijiang River

    表 1  基于COⅠ基因序列的鲫群体遗传多样性参数

    Table 1.  Genetic diversity ind in C.auratus populations based on COⅠ sequences

    群体
    population
    经纬度
    latitude and longitude
    样本数
    numbr of samples
    单倍型数
    number of haplotypes
    单倍型多样性指数 (Hd)
    haplotype diversity
    核苷酸多样性指数 (π)
    nucleotide diversity
    鲁布革 Lubuge 104°31′45″E, 24°44′46″N 17 5 0.809±0.057 0.002 8±0.000 7
    八渡 Badu 105°48′26″E, 24°42′8″N 17 6 0.721±0.087 0.002 7±0.000 9
    天峨 Tian'e 108°52′22″E, 23°48′43″N 7 3 0.714±0.127 0.001 4±0.000 3
    大化 Dahua 107°59′16″E, 23°44′5″N 9 6 0.917±0.073 0.006 5±0.001 3
    合山 Heshan 110°04′19″E, 23°24′16″N 19 8 0.836±0.057 0.008 7±0.000 9
    桂平 Guiping 110°53′6″E, 23°21′46″N 21 8 0.833±0.055 0.005 8±0.001 1
    藤县 Tengxian 112°27′33″E, 23°4′54″N 18 9 0.895±0.048 0.006 4±0.001 0
    肇庆 Zhaoqing 110°04′20″E, 23°24′15″N 23 12 0.901±0.041 0.008 1±0.000 9
    合计 total 131 25 0.900±0.017 0.006 7±0.000 4
    下载: 导出CSV

    表 2  基于COⅠ基因对珠江干流鲫的AMOVA分析

    Table 2.  AMOVA analysis of C.auratus populations based on COⅠ gene

    变异来源
    source of variation
    自由度
    df
    平方和
    SS
    变异组分
    variance components
    变异百分数
    percentage of variation
    FST P
    群体间 among populations 7 52.692 0.354 46 Va 16.40 0.163 97 0.000 00
    群体内 within populations 123 222.301 1.807 32 Vb 83.60
    总计 total 130 274.992 2.161 78
    下载: 导出CSV

    表 3  基于COⅠ序列的野生鲫群体间的遗传分化系数 (FST)

    Table 3.  FST analysis among populations of C.auratusbased on COⅠ sequences

    群体
    population
    鲁布革
    Lubuge
    八渡
    Badu
    天峨
    Tian'e
    大化
    Dahua
    合山
    Heshan
    桂平
    Guiping
    藤县
    Tengxian
    鲁布革 Lubuge
    八渡 Badu 0.010 27
    天峨 Tian'e – 0.020 33 0.098 00*
    大化 Dahua 0.124 63* 0.176 13* 0.099 99
    合山 Heshan 0.123 91* 0.165 03* 0.118 59 0.061 15
    桂平 Guiping 0.033 80 0.094 36* 0.043 29 0.081 67* 0.081 05*
    藤县 Tengxian 0.279 76 * 0.319 31* 0.339 04* 0.227 87* 0.130 90* 0.144 19*
    肇庆 Zhaoqing 0.279 36 * 0.313 49* 0.307 59* 0.190 04* 0.148 04* 0.150 29* – 0.019 77
     注:*. P<0.05
    下载: 导出CSV

    表 4  基于COⅠ基因对野生鲫群体的中性检验

    Table 4.  Neutrality test of C.auratus populations based on COⅠ gene

    群体
    population
    Tajima's D Fu's Fs
    D P Fs P
    鲁布革 Lubuge – 0.461 0.344 0.136 0.565
    八渡 Badu – 1.550 0.054 – 1.005 0.227
    天峨 Tian'e 0.206 0.598 – 0.237 0.241
    大化 Dahua – 0.326 0.412 – 0.438 0.338
    合山 Heshan 0.075 0.576 0.757 0.666
    桂平 Guiping – 0.645 0.278 – 0.096 0.505
    藤县 Tengxian – 0.838 0.205 – 1.042 0.318
    肇庆 Zhaoqing – 0.785 0.244 – 1.869 0.226
    总体 total – 0.540 0.339 – 0.474 0.386
    下载: 导出CSV
  • [1] 帅方敏, 李新辉, 刘乾甫, 等. 珠江水系鱼类群落多样性空间分布格局[J]. 生态学报, 2017, 37(9): 3182-3192
    [2] 李捷, 李新辉, 贾晓平, 等. 西江鱼类群落多样性及其演变[J]. 中国水产科学, 2010, 17(2): 298-311
    [3] WARD R D. Genetics in fisheries management[J]. Hydrobiologia, 2000, 420: 191-201
    [4] STENSETH N C, ROUYER T. Ecology-destabilized fish stocks[J]. Nature, 2008, 452(7189): 825-826
    [5] BARASA J E, ABILA R, GROBLER J P, et al. High genetic diversity and population differentiation in Clarias gariepinus of Yala Swamp: evidence from mitochondrial DNA sequences[J]. J Fish Biol, 2016, 89(6): 2557-2570
    [6] WAN Y, ZHOU C H, OUYANG S, et al. Genetic diversity analysis of mitochondrial DNA control region in artificially propagated Chinese sucker Myxocyprinus asiaticus[J]. Mitochondr DNA, 2015, 26(4): 514-519
    [7] 桂建芳. 鱼类性别和生殖的遗传基础及其人工控制[M]. 北京: 科学出版社, 2007: 107-124.
    [8] RUMISHA C, LEERMAKERS M, MDEGELA R H, et al. Bioaccumulation and public health implications of trace metals in edible tissues of the crustaceans Scylla serrata and Penaeus monodon from the Tanzanian coast[J]. Environ Monit Assess, 2017, 189(10): 529
    [9] LIU Z J, CORDES J F. DNA marker technologies and their applications in aquaculture genetics[J]. Aquaculture, 2004, 238(1/2/3/4): 1-37
    [10] AVISE J C. Molecular markers, natural history, and evolution[M]. 2nd ed. London: Chapman and Hall, 2004: 684.
    [11] HABIB K A, JEONG D, MYOUNG J G, et al. Population genetic structure and demographic history of the fat greenling Hexagrammos otakii[J]. Gene Genom, 2011, 33(4): 413-423
    [12] SHEN Y J, GUAN L H, WANG D Q, et al. DNA barcoding and evaluation of genetic diversity in Cyprinidae fish in the midstream of the Yangtze River[J]. Ecol Evol, 2016, 6(9): 2702-2713
    [13] 刘连为, 许强华, 陈新军. 基于线粒体COⅠ和Cyt b基因序列的北太平洋柔鱼种群遗传结构研究[J]. 水产学报, 2012, 36(11): 1675-1684
    [14] THOMPSON J D, GIBSON T J, PLEWNIAK F, et al. The CLUSTAL_X Windows interface: flexible strategies for multiple sequence alignment aided by quality analysis tools[J]. Nucleic Acids Res, 1997, 25(24): 4876-4882
    [15] TAMURA K, PETERSON D, PETERSON N, et al. MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods[J]. Mol Biol Evol, 2011, 28(10): 2731-2739
    [16] FELSENSTEIN J. Confidence-limits on phylogenies: an approach using the bootstrap[J]. Evolution, 1985, 39(4): 783-791
    [17] BANDELT H J, FORSTER P, RÖHL A. Median-joining networks for inferring intraspecific phylogenies[J]. Mol Biol Evol, 1999, 16(1): 37-48
    [18] LIBRADO P, ROZAS J. DnaSP v5: a software for comprehensive analysis of DNA polymorphism data[J]. Bioinformatics, 2009, 25(11): 1451-1452
    [19] EXCOFFIER L, LAVAL G, SCHNEIDER S. Arlequin Ver. 3.01: an integrated software package for population genetics data analysis[J]. Evol Bioinform, 2005, 1: 47-50
    [20] WANG W, MA C Y, CHEN W, et al. Population genetic diversity of Chinese sea bass (Lateolabrax maculatus) from southeast coastal regions of China based on mitochondrial COⅠ gene sequences[J]. Biochem Syst Ecol, 2017, 71: 114-120
    [21] XU D D, LOU B, SHI L H, et al. Genetic diversity and population structure of Nibea albiflora in the China sea revealed by mitichondrial COⅠ sequences[J]. Biochem Syst Ecol, 2012, 45: 158-165
    [22] LI S F, YANG Q L, XU J W, et al. Genetic diversity and variation of mitochondrial DNA in native and introduced bighead carp[J]. Trans Am Fish Soc, 2010, 139(4): 937-946
    [23] GROOMBRIDGE J, JONES C, BRUFORD M, et al. 'Ghost' alleles of the Mauritius kestrel[J]. Nature, 2000, 403(6770): 616-616
    [24] GRANT W S, BOWEN B W. Shallow population histories in deep evolutionary lineages of marine fishes: insights from sardines and anchovies and lessons for conservation[J]. J Hered, 1998, 89(5): 415-426
    [25] 袁振兴, 代应贵, 刘伟, 等. 都柳江鲤、鲫和草鱼种群mtDNA控制区及遗传多样性分析[J]. 基因组学与应用生物学, 2017, 36(5): 1926-1934
    [26] 邓朝阳. 鲫4个群体线粒体Cyt b序列和D-loop区的比较分析[J]. 现代农业科技, 2015(9): 277-279, 283
    [27] MA C Y, CHENG Q Q, ZHANG Q Y, et al. Genetic variation of Coilia ectenes (Clupeiformes: Engraulidae) revealed by the complete cytochrome b sequences of mitochondrial DNA[J]. J Exp Mar Bio Ecol, 2010, 385(1/2): 14-19
    [28] HEWITT G M. Some genetic consequences of ice ages, and their role in divergence and speciation[J]. Biol J Linnean Soc, 1996, 58(3): 247-276
    [29] HEWITT G. The genetic legacy of the Quaternary ice ages[J]. Nature, 2000, 405(6789): 907-913
    [30] JU L, WANG H, JIANG D. Simulation of the last glacial maximum climate over East Asia with a regional climate model nested in a general circulation model[J]. Palaeogeogr Palaeoclimatol Palaeoecol, 2007, 248(3/4): 376-390
    [31] WEAVER L M, GAN S, QUIRINO B, et al. A comparison of the expression patterns of several senescence-associated genes in response to stress and hormone treatment[J]. Plant Mol Biol, 1998, 37(3): 455-469
    [32] CHEN W T, ZHONG Z X, DAI W, et al. Phylogeographic structure, cryptic speciation and demographic history of the sharpbelly (Hemiculter leucisculus), a freshwater habitat generalist from southern China[J]. BMC Evol Biol, 2017, 17(1): 216
  • [1] 李莉好喻达辉黄桂菊杜博符云童馨郭奕惠叶卫 . 吉富罗非鱼不同选育群体的遗传多样性. 南方水产科学, 2007, 3(5): 40-48.
    [2] 马国强高天翔孙典荣 . 基于16S rRNA 和Cyt b 基因序列探讨2 种梅童鱼的遗传分化. 南方水产科学, 2010, 6(2): 13-20. doi: 10.3969/j.issn.1673-2227.2010.02.003
    [3] 孙奉玉宋忠魁赵鹏聂振平苏琼王芳宇 . 广西沿海及其邻近海区拟穴青蟹群体遗传多样性的RAPD分析. 南方水产科学, 2012, 8(2): 30-35. doi: 10.3969/j.issn.2095-0780.2012.02.005
    [4] 郑德育郭易佳杨天燕高天翔郑瑶袁冬皓斯舒谨 . 基于线粒体ND2基因序列的少鳞遗传多样性研究. 南方水产科学, 2019, 15(5): 84-91. doi: 10.12131/20190042
    [5] 王春晓高风英卢迈新刘志刚朱华平叶星 . 2个尼罗罗非鱼群体GHSR基因5侧翼序列的多态性及其遗传多样性分析. 南方水产科学, 2015, 11(1): 18-25. doi: 10.3969/j.issn.2095-0780.2015.01.003
    [6] 杨喜书章群余帆洋吕金磊底晓丹邵伟军黄镇宇卢丽锋 . 华南6水系与澜沧江-湄公河攀鲈线粒体ND2基因的遗传多样性分析. 南方水产科学, 2017, 13(3): 43-50. doi: 10.3969/j.issn.2095-0780.2017.03.006
    [7] 王小玉喻达辉黄桂菊郭奕惠杜博龚世园王爱民 . 合浦珠母贝3个家系的AFLP标记分离与遗传多样性研究. 南方水产科学, 2007, 3(5): 54-60.
    [8] 李敏张鹏李玉芳陈森张魁孔啸兰陈作志 . 南海扁舵鲣种群遗传结构和遗传多样性评价. 南方水产科学, 2015, 11(5): 82-89. doi: 10.3969/j.issn.2095-0780.2015.05.010
    [9] 莫艳秀王晓清莫永亮 . 长吻鮠遗传多样性的RAPD分析. 南方水产科学, 2010, 6(6): 77-80. doi: 10.3969/j.issn.1673-2227.2010.06.014
    [10] 孙立元郭华阳朱彩艳马振华江世贵张殿昌 . 卵形鲳鲹育种群体遗传多样性分析. 南方水产科学, 2014, 10(2): 67-71. doi: doi:10.3969/j.issn.2095-0780.2014.02.010
    [11] 范嗣刚王婧璇黄桂菊刘宝锁郭奕惠喻达辉1 . 合浦珠母贝选育家系的遗传多样性分析. 南方水产科学, 2016, 12(5): 90-96. doi: 10.3969/j.issn.2095-0780.2016.05.011
    [12] 孙成飞叶星董浚键田园园梁健辉 . 罗氏沼虾6个养殖群体遗传多样性的微卫星分析. 南方水产科学, 2015, 11(2): 20-26. doi: 10.3969/j.issn.2095-0780.2015.02.003
    [13] 朱彩艳叶卫夏军红符云周发林江世贵 . 广东1个鲮原种群体的种质特征及遗传多样性分析. 南方水产科学, 2005, 1(4): 1-5.
    [14] 孙成飞谢汶峰胡婕董浚键田园园吴灶和叶星 . 大口黑鲈3个养殖群体的遗传多样性分析. 南方水产科学, 2019, 15(2): 64-71. doi: 10.12131/20180203
    [15] 吴文婷张磊宋志民赵翠汤晓荣刘涛陈伟洲 . 5个紫菜属物种丝状体的微卫星初步遗传分析. 南方水产科学, 2012, 8(4): 29-36. doi: 10.3969/j.issn.2095-0780.2012.04.005
    [16] 符云钟金香颉晓勇叶卫林碧海陈辉崇张汉华 . 罗非鱼3个养殖群体的遗传多样性及特异性AFLP标记研究. 南方水产科学, 2008, 4(6): 50-55.
    [17] 原居林朱俊杰王高学 . 秦岭细鳞鲑黑河种群和湑水河种群的遗传多样性分析. 南方水产科学, 2009, 5(3): 63-66. doi: 10.3969/j.issn.1673-2227.2009.03.011
    [18] 荣朝振祖国掌胡建华孙守旗孙棠丽 . 泥鳅线粒体DNA控制区结构分析及遗传多样性研究. 南方水产科学, 2011, 7(5): 55-62. doi: 10.3969/j.issn.2095-0780.2011.05.009
    [19] 颉晓勇苏天凤陈文张志颜远义江世贵 . 凡纳滨对虾6个养殖群体遗传多样性的比较分析. 南方水产科学, 2008, 4(6): 42-49.
    [20] 赵彦花区又君温久福李加儿周慧 . 基于微卫星标记的黄唇鱼遗传多样性研究. 南方水产科学, 2019, 15(4): 127-132. doi: 10.12131/20180261
  • 加载中
图(5)表(4)
计量
  • 文章访问数:  1908
  • HTML全文浏览量:  636
  • PDF下载量:  113
  • 被引次数: 0
出版历程
  • 收稿日期:  2018-05-17
  • 录用日期:  2018-05-30
  • 刊出日期:  2018-10-01

西江鲫的遗传多样性与群体结构

    作者简介:匡天旭(1991 — ),男,硕士研究生,从事鱼类分子生物学研究。E-mail: 1462019455@qq.com
    通讯作者: 李新辉, lxhui01@aliyun.com
  • 1. 中国水产科学研究院珠江水产研究所,广东 广州 510380
  • 2. 上海海洋大学海洋科学学院,上海 201306

摘要: 为了解珠江鲫(Carassius auratus)遗传种质资源现状,文章以西江为例,利用线粒体细胞色素c氧化酶亚基Ⅰ (COⅠ)基因研究了西江8个地理群体的遗传多样性和群体结构。结果表明,西江鲫群体的遗传多样性处于较高水平。鲫群体存在显著遗传分化(FST=0.164,P<0.001),地理距离和自身生活习性可能是导致鲫群体产生遗传分化的重要因素。此外,中性检验和错配分布分析表明西江鲫没有经历近期的群体扩张。

English Abstract

  • 珠江由东江、西江和北江组成,生境异质性高,是目前生物多样性研究的热点地区之一[1]。西江全长2 214 km,流域面积为3.53×105 km2,年径流量达2 240×108 m3,孕育了丰富的鱼类资源[2]。然而,由于受到水利设施、环境污染、过度捕捞等因素的影响,西江鱼类资源退化严重[3-6],因此了解西江鱼类的遗传多样性现状和群体结构对今后鱼类资源利用与保护尤为必要。

    鲫(Carassius auratus)属于鲤形目、鲤科、鲫属,为常见定居型鱼类,广泛分布于欧洲和亚洲诸多水系。由于其肉质细嫩鲜美,营养价值高,深受消费者喜爱,中国年消费量达200×104 t[7]。由于鲫具有较强的环境适应能力且易于被遗传改良,常被用于水产养殖[7]。鲫是珠江水系传统渔业的主要捕捞对象之一,但近几十年来其产量持续走低[1,8]。为了合理利用珠江鲫的种质资源,弄清该物种的遗传多样性现状和群体结构是关键。本研究利用线粒体细胞色素c氧化酶亚基Ⅰ (COⅠ)基因[9]作为分子标记[10],对西江水系鲫的8个地理群体共131尾样本开展遗传多样性和群体结构分析,旨在了解西江水系在高强度梯级开发、环境变化和人为干扰等条件下鲫的遗传多样性现状和群体结构,为鲫的资源利用和保护提供科学依据[11-12]

    • 2016年在西江干流共设8个采样点,各采样点间平均间隔约200 km (图1表1)。使用流刺网在8个采样站位共采集鲫样本131尾,剪取胸鳍或尾鳍保存于无水乙醇中,带回实验室后转移至 – 40 ℃冰箱内储存。

      图  1  鲫样本采集地点

      Figure 1.  Sampling site of C.auratus

    • 提取总DNA时,选取20 mg鳍条组织并将其铰成细小粒状,选用Axygen总基因组试剂盒进行提取,提取方法及步骤按照该试剂盒使用说明进行。使用通用引物(FishF1和FishR1)[13]对每尾样本进行PCR扩增。PCR反应体系为40 μL,10×buffer 4 μL,dNTPs (10 mmol·L–1) 2 μL,正反向引物各1 μL,Taq酶(5 U) 0.2 μL,基因组DNA 2 μL做反应模板,加超纯水补至40 μL。PCR反应条件为95 ℃预变性3 min;94 ℃变性30 s,54 ℃退火30 s,72 ℃延伸1 min,循环35次;最后72 ℃延伸10 min。扩增产物经1 %琼脂糖凝胶电泳检测后送至广州艾基生物技术有限公司进行双向测序,测序引物与扩增引物相同。

    • 使用Seqman软件对正反向测序结果进行拼接,运用Clustal X软件[14]对序列进行比对、编辑与剪切并辅以手工调整。使用MEGA 7.0 软件[15]计算遗传距离,模型选择Kimura双参数模型(K2P)。以鲤(Cyprinus carpio)作外类群,基于单倍型序列以邻接法(neighbor-joining,NJ)建树[16],使用1 000次重复抽样(bootstrap)检测其置信度,所用模型为K2P距离模型。使用Network 4.6软件[17]绘制单倍型网络图,并利用Dnasp 5.0软件 [18] 统计所有个体和每个群体的单倍型数,计算单倍型多样性(Hd)和核苷酸多样性(π)。使用Arlequin 3.5软件 [19],根据pairwise difference模型,估算两两群体间的遗传分化指数(FST),并基于Mantel检验分析地理距离是否与遗传分化存在联系。两两群体之间的地理距离大致通过采样点的经纬度坐标估算获得。本研究还对鲫群体进行分子方差分析(AMOVA),将干流的8个群体划分为一个组群,验证8个群体之间是否存在显著的遗传分化。最后,使用中性检验Fu's Fs和Tajima's D及核苷酸错配分布分析(mismatch distribution analysis)推断研究群体是否经历了近期的群体扩张。FST、AMOVA、中性检验和错配分布分析均采用1 000次重复。

    • 本研究成功获得了131尾鲫的COⅠ基因片段。在序列两端截齐后,序列长度为636 bp,其中包含了36个变异位点,定义了25个单倍型。所有序列均未发现插入和缺失。邻接树的支持率较低,并未发现西江鲫形成了显著的谱系分枝(图2)。单倍型网络图表明,15个单倍型(H3~H9、H11、H12、H14~H18、H21)被多个地理群体共享。出现频率最高的单倍型是H4,占所有检测个体的26.0%,并且地理分布最广,为8个地理群体共享;其次出现频率较高的是H15,占总个体的10.7% (14/131)。有10个单倍型为单个地理群体所特有,其中,单倍型H1、H2仅在上游南盘江八渡江段分布;单倍型H10仅在中游红水河大化江段分布;单倍型H13等5个样本仅在中游红水河合山江段分布;单倍型H19、H20仅在下游西江藤县江段分布;单倍型H22、H23、H24和H25仅在下游西江肇庆江段分布。各单倍型网络关系分解见图3

      图  2  基于COⅠ基因构建的单倍型NJ系统发育树

      Figure 2.  Neighbor-Joining tree based on COⅠ gene

      图  3  基于COⅠ基因构建的鲫单倍型网络图

      Figure 3.  Median-joining haplotype network of C.auratusbased on COⅠ gene

      鲫总的单倍型多样性为0.900±0.017,核苷酸多样性为0.006 7±0.000 4,平均核苷酸差异数为4.231。各个地理群体的单倍型多样性介于0.714~0.917,核苷酸多样性介于0.001 4~0.008 7。此外,本研究发现上游站位的单倍型多样性和核苷酸多样性均低于中下游站位(表1)。

      群体
      population
      经纬度
      latitude and longitude
      样本数
      numbr of samples
      单倍型数
      number of haplotypes
      单倍型多样性指数 (Hd)
      haplotype diversity
      核苷酸多样性指数 (π)
      nucleotide diversity
      鲁布革 Lubuge 104°31′45″E, 24°44′46″N 17 5 0.809±0.057 0.002 8±0.000 7
      八渡 Badu 105°48′26″E, 24°42′8″N 17 6 0.721±0.087 0.002 7±0.000 9
      天峨 Tian'e 108°52′22″E, 23°48′43″N 7 3 0.714±0.127 0.001 4±0.000 3
      大化 Dahua 107°59′16″E, 23°44′5″N 9 6 0.917±0.073 0.006 5±0.001 3
      合山 Heshan 110°04′19″E, 23°24′16″N 19 8 0.836±0.057 0.008 7±0.000 9
      桂平 Guiping 110°53′6″E, 23°21′46″N 21 8 0.833±0.055 0.005 8±0.001 1
      藤县 Tengxian 112°27′33″E, 23°4′54″N 18 9 0.895±0.048 0.006 4±0.001 0
      肇庆 Zhaoqing 110°04′20″E, 23°24′15″N 23 12 0.901±0.041 0.008 1±0.000 9
      合计 total 131 25 0.900±0.017 0.006 7±0.000 4

      表 1  基于COⅠ基因序列的鲫群体遗传多样性参数

      Table 1.  Genetic diversity ind in C.auratus populations based on COⅠ sequences

      AMOVA分析发现鲫群体间存在显著的遗传分化(FST=0.164, P=0.000),占总变异的16.40%(表2)。在估算两两群体之间的遗传分化系数时发现,鲁布革与大化、合山、藤县、肇庆群体间均存在显著的遗传分化(表3)。此外,八渡群体除与鲁布革群体未检测到显著遗传分化外,与其他江段群体都具有显著的遗传分化;藤县和肇庆2个群体与其他群体之间都有显著分化。但藤县和肇庆之间的FST为负值,说明两地理群体之间基因交流频繁。此外,藤县与距离较远的群体(鲁布革、八渡、天峨)之间的遗传分化值要明显大于距离较近的群体(大化、合山、桂平)。Mantel检验发现群体之间的地理距离与遗传分化呈现出显著的正相关关系(Mantel R=0.661,P=0.004,图4)。然而,部分距离较近的群体之间也存在显著的遗传分化,例如八渡与天峨之间,以及藤县与桂平之间。

      变异来源
      source of variation
      自由度
      df
      平方和
      SS
      变异组分
      variance components
      变异百分数
      percentage of variation
      FST P
      群体间 among populations 7 52.692 0.354 46 Va 16.40 0.163 97 0.000 00
      群体内 within populations 123 222.301 1.807 32 Vb 83.60
      总计 total 130 274.992 2.161 78

      表 2  基于COⅠ基因对珠江干流鲫的AMOVA分析

      Table 2.  AMOVA analysis of C.auratus populations based on COⅠ gene

      群体
      population
      鲁布革
      Lubuge
      八渡
      Badu
      天峨
      Tian'e
      大化
      Dahua
      合山
      Heshan
      桂平
      Guiping
      藤县
      Tengxian
      鲁布革 Lubuge
      八渡 Badu 0.010 27
      天峨 Tian'e – 0.020 33 0.098 00*
      大化 Dahua 0.124 63* 0.176 13* 0.099 99
      合山 Heshan 0.123 91* 0.165 03* 0.118 59 0.061 15
      桂平 Guiping 0.033 80 0.094 36* 0.043 29 0.081 67* 0.081 05*
      藤县 Tengxian 0.279 76 * 0.319 31* 0.339 04* 0.227 87* 0.130 90* 0.144 19*
      肇庆 Zhaoqing 0.279 36 * 0.313 49* 0.307 59* 0.190 04* 0.148 04* 0.150 29* – 0.019 77
       注:*. P<0.05

      表 3  基于COⅠ序列的野生鲫群体间的遗传分化系数 (FST)

      Table 3.  FST analysis among populations of C.auratusbased on COⅠ sequences

      图  4  地理距离与两两群体间遗传分化的关系

      Figure 4.  Relationship between geographic distance and pairwise differentiation

      Tajima's D和Fu's Fs在总的群体和单个地理群体中均为不显著负值(表4),表明西江鲫群体可能并未经历近期的群体扩张。另外,错配分布分析未显示单峰分布(图5),并不支持群体历史上有过扩张。

      群体
      population
      Tajima's D Fu's Fs
      D P Fs P
      鲁布革 Lubuge – 0.461 0.344 0.136 0.565
      八渡 Badu – 1.550 0.054 – 1.005 0.227
      天峨 Tian'e 0.206 0.598 – 0.237 0.241
      大化 Dahua – 0.326 0.412 – 0.438 0.338
      合山 Heshan 0.075 0.576 0.757 0.666
      桂平 Guiping – 0.645 0.278 – 0.096 0.505
      藤县 Tengxian – 0.838 0.205 – 1.042 0.318
      肇庆 Zhaoqing – 0.785 0.244 – 1.869 0.226
      总体 total – 0.540 0.339 – 0.474 0.386

      表 4  基于COⅠ基因对野生鲫群体的中性检验

      Table 4.  Neutrality test of C.auratus populations based on COⅠ gene

      图  5  西江鲫群体核苷酸错配分布分析

      Figure 5.  Mismatch distribution analysis of C.auratus in Xijiang River

    • 单倍型多样性(Hd)和核苷酸多样性(π)是衡量群体遗传多样性高低的重要指标[8]。诸多内外因素均会导致物种遗传多样性水平发生变化[20-23]。西江鲫总群体的Hd为0.900±0.017,π为0.006 7±0.000 4。依据Grant和Bowen[24]的标准,西江鲫的遗传多样性属于较高水平。袁振兴等[25]对线粒体DNA控制区 (D-loop)进行分析,发现贵州都柳江鲫的遗传多样性较高(Hd=0.892 6,π=0.006 32);邓朝阳[26]通过线粒体D-loop序列分析鲫的遗传多样性发现,分布于长江水系的野生鲫的遗传多样性同样较高(Hd =0.933,π=0.005 28)。另外一个有趣的发现是,鲁布革、八渡、天峨等群体的遗传多样性明显低于较下游的群体,大量的水利开发可能是导致这一结果的重要原因。由于大量的水电站建设,特别是中游红水河为十级梯级开发,已经将西江原河流截为多段,导致上游群体与中下游群体相对隔离,上游的鱼受泄洪及船闸排水的影响可迁移至下游江段,而中下游的鱼则难以迁移至较上游的区域,导致遗传多样性的差异分布[27]

      AMOVA分析和遗传分化系数(FST)估算表明西江鲫群体之间存在显著的遗传分化,而中下游群体具有的一些特有单倍型也从侧面反映出群体之间显著的遗传分化。FST估算发现距离较近的群体的遗传分化值往往小于距离远的群体,比如鲁布革群体与较近群体(八渡、天峨群体)之间的遗传分化值较小且不显著,而与相对较远的群体(鲁布革、天峨群体)之间的遗传分化值较高,达0.3。在Mantel检验中也发现群体之间遗传分化随着地理距离增大而显著增加。因此,地理距离是导致鲫遗传分化的一个重要因素。此外,本研究还发现部分距离较近的群体之间也存在显著的遗传分化,如藤县群体与桂平群体。自身的生活习性可能是导致距离较近群体间遗传分化的主要原因。鲫营定居性生活史,其迁徙能力有限,习惯生活于特定环境;而水利梯级开发将原河流截断,导致部分群体处于相对隔离的环境,限制了西江鲫群体之间的迁移扩散。

      大量研究表明,周期性的冰川循环(glacial cycles)对许多现存物种的遗传多样性和分布都有重要的影响,但是不同物种或不同栖息环境的物种往往对冰期循环有着不同的响应[28-29]。本研究利用中性检验和错配分布分析了西江鲫的群体历史动态,发现近期西江鲫群体保持相对稳定的规模,并未经历群体扩张。这可能是由于西江地处热带和亚热带区域,在更新世冰期期间(glacial period)气候温和[30-31],鱼类能够维持其相对稳定的生态位,进而能保持相对稳定的群体数量。有研究表明,珠江流域的䱗(Hemiculter leucisculus)也经历了类似的群体历史动态[32]

      本研究探讨了西江鲫种群的遗传多样性现状和群体结构,结果显示鲫的遗传多样性处于相对较高的水平,地理距离和自身特定的生活史可能促进了鲫的遗传分化,西江鲫群体数量保持了相对稳定的状态。本研究结果为鲫的资源利用与保护以及今后珠江水系其他鱼类遗传多样性和群体结构的研究提供理论参考依据。由于本研究的分子标记单一且采样范围有限,因此后续研究有必要选取更多的分子标记(微卫星标记、简化基因组)并加大采样范围来进一步阐释鲫群体的遗传多样性现状和群体结构。

参考文献 (32)

目录

    /

    返回文章
    返回