基于线粒体16SrRNA探讨渤海海区矮拟帽贝的种群生态遗传

史博洋， 李媛媛， 王儒晓， 曲江勇, 2*， 郭承华*， 刘传林， 侯建海， 房辉， 刘楚童

(1.烟台大学生命科学学院，山东烟台264005； 2. 海南热带海洋学院，海南两栖爬行动物研究重点实验室，海南三亚572000)



基于线粒体16SrRNA探讨渤海海区矮拟帽贝的种群生态遗传

史博洋1， 李媛媛1， 王儒晓1， 曲江勇1, 2*， 郭承华1*， 刘传林1， 侯建海1， 房辉1， 刘楚童1

(1.烟台大学生命科学学院，山东烟台264005； 2. 海南热带海洋学院，海南两栖爬行动物研究重点实验室，海南三亚572000)

采用聚合酶链式反应和直接测序的方法研究了渤海海区15个矮拟帽贝Patelloidapygmaea种群遗传多样性与环境因子的关系。180个个体的16S rRNA部分序列(615 bp)共发现264个多态性位点，153个转换，114个颠换，界定32个单倍型。所有个体碱基含量依次为C(18.12%)、T(27.74%)、A(33.07%)、G(21.07%)；单倍型多样性为0.835 1±0.019 1，核苷酸多样性为0.449±0.214。矮拟帽贝大陆沿岸种群核苷酸多样性(0.191 8±0.091 7)高于岛屿种群(0.107 3±0.081 6)。种群历史动态结果显示，矮拟帽贝经历了种群近期扩张。与环境因子相关性研究发现，遗传距离、序列变异、核苷酸多样性与日照时数、日照百分率分别呈显著正相关(P<0.05)；形态相关分析结果表明，变异主成分(壳长、壳宽)与年均降水量呈显著负相关(P<0.05)；年均降水量、日照时数、日照百分率是影响矮拟帽贝形态变异、遗传多样性的主导因子，环境越稳定地区种群的遗传多样性越高。

矮拟帽贝； 16S rRNA； 生态遗传； 遗传多样性

环境因子在物种的进化过程中发挥着重要作用，物种的分布、行为、生理和遗传结构通常都会受到环境因子的影响，探究环境因子与生物多样性的相关性是生态遗传研究的重要内容之一(Scheiner，1993)。通常种群占据不同地理环境并呈现梯度性遗传变化，这种结果主要受自然选择与环境变迁影响(Jin & Liu，2008)，遗传结构分化和自然环境变迁之间的相互作用是目前生态地理学关注的焦点。同类物种生活在不同环境下(温度、降水、日照等)的地理种群遗传多样性存在地理变异，如荒漠沙蜥Phrynocephalusprzewalskii的遗传多样性与纬度、降水量显著正相关，大石鸡Alectorismagna、雉鸡甘肃亚种Phasianuscolchicusstrauchi的遗传多样性与温度呈显著负相关，以及同种生物因环境而产生的形态改变影响基因频率等(Michod，1981；Nevo，1981；曲江勇，刘迺发，2012)。环境对生物的影响(环境的演替、种群遗传结构的多变)伴随着表型与基因型的变化继而改变种群结构，环境促进演替使种群遗传分化趋于平衡稳定(Slatkin，1987)。

分子标记作为研究种群遗传结构和地理分化的重要工具，目前已广泛用于解决动物生态遗传等问题(Wayne & Robert，2004)。动物线粒体DNA具有分子量小、进化速度快、结构简单、表现为母系遗传等特点，非常适于动物群体遗传学和系统进化研究(Canapaetal.，2000)，其中16S rRNA基因进化速率适中，适合种群水平的遗传多样性及遗传结构研究(Hondaetal.，2002)。

矮拟帽贝Patelloidapygmaea隶属软体动物门Mollusca腹足纲Gastropoda原始腹足目Archaeogastropoda笠贝科Acmaeidae(Surhoneetal.，2010)，中国习见种类，沿岸、岛屿均有分布。矮拟帽贝遍布整个潮间带(海岸带以上到浅层潮下带)，栖息附着在大型水生植物、激浪区礁岩、珊瑚藻及其他贝壳等表面，主要摄食藻类与浮游植物(Lindberg，1990；Takenori & Takashi，1993；Kirkendale & Meyer，2004；Nakano & Ozawa，2007)。目前矮拟帽贝研究主要集中在系统学和分类学等领域，包括种水平鉴定、形态变异、遗传多样性等(Nakaietal.，2006；Gonzálezwevaretal.，2010)，种群系统地理研究较少。国内研究重点多为具有较大经济价值的贝类，分布广泛但无经济价值的潮间带贝类的系统地理格局受历史因素、现今环境和生态因素的影响，其遗传进化、系统地理分布格局未得到深入研究(Nietal.，2012)。

本研究以矮拟帽贝线粒体16S rRNA部分序列，探究矮拟帽贝的生态遗传，揭示种群遗传多样性随环境因子的变异规律，为中国矮拟帽贝种群遗传研究提供理论依据。

1 材料与方法

1.1样品采集、DNA提取及序列获取

矮拟帽贝采样点见图1，大陆沿岸种群包括大连、庄河、营口、锦州、葫芦岛、秦皇岛、莱州、蓬莱、烟台；岛屿种群包括长岛、砣矶岛、大钦岛、南隍城岛、北隍城岛。采集后测量个体的壳长、壳宽、壳高、壳周长、长斜边、短斜边、壁厚等；返回实验室称量壳干重、计算头端角度与尾端角度(表1)，无水乙醇保存备用。

图1 矮拟帽贝地理种群采集点Fig. 1 The sampling sites of Patelloida pygmaea in the Bohai Sea

DNA提取：DNeasy©组织&血液试剂盒(GTIN: 04053228006077，Lot：151046088，GIAGEN，德国)。PCR扩增：选用16S rRNA通用引物(华大基因合成)，16Sar-L：5′-CGCCTGTTTATCAAAAACAT-3′；16SbrH：5′-CCGGTCTGAACTCAGATCACGT-3′(Folmeretal.，1994)，PCR扩增反应体系(25 μL)：12.5 μL Taq MIX(0.625 U Taq酶、0.2 mmol·L-1dNTP、1.5 mmol·L-1MgCl2；Lot：A3701A，TaKaRa，中国大连)；0.15 μL (100 μM)上下游引物、1 μL DNA模板(80～100 ng)。PCR反应条件设置：94 ℃预变性2 min；94 ℃变性30 s，49 ℃退火40 s，72 ℃延伸50 s，35个循环；72 ℃延伸5 min。PCR结束后取4 μL于1%琼脂糖凝胶电泳检测(1×TBE、5 V/cm恒压、Goldenview染色)(琼脂糖，Lot：142075，BIOWEST，西班牙；Goldenview，Lot：OHO9RL，济南力戈科技有限公司)。PCR产物送北京六合华大公司进行ABI PRIMS 3730双向测序。

1.2环境因子

选取经度、纬度、年均温度、温度变异系数、年均降水量、降水量变异系数、年均日照时数、日照率作为环境因子(表2)，从国家气象数据共享服务网获取10个地点近30年气象数据，并计算气温、降水量变异系数(砣矶岛、大钦岛、小钦岛、北隍城岛、南隍城岛等地区气象站因缺少相关数据，未被国家气象数据共享服务网收录)。

1.3数据分析

双向测序结果用Conting进行拼接，通过CLUSTAL_X 1.83进行序列比对，单倍型用DnaSP v5统计(Librado & Rozas，2009)。通过ARLEQUIN 3.5计算群体多样性指数：单倍型多样性(h)、核苷酸多样度(π)、序列平均核苷酸差异数(k)、扩张系数(Tau)、有效种群大小(θ)、种群分化系数(Fst)；AMOVA分析(Excoffier & Lischer，2010)；通过MEGA7计算遗传距离(D)与序列变异(S)(Kumaretal.，2016)；样品所测壳长、壳宽、壳高、壳周长、短斜边、长斜边、头端角度、尾端角度、干重、壁厚等形态数据用平均数±标准差表示，形态变异采用主成分分析与相关性分析，以上统计分析均在SPSS 19.0(Arbuckle，2010)完成。

2 结果

2.1序列变异与遗传多样性

180个矮拟帽贝样本的mtDNA 16S rRNA部分序列中，共发现了264个多态位点，碱基组成为C(18.12%)、T(27.74%)、A(33.07%)、G(21.07%)，转换/颠换比为1.34(缺失位点少于5%)，界定32种单倍型。总体单倍型多样性为0.835 1，核苷酸多样性为0.449，核苷酸两两差异度为118.533。除庄河、蓬莱外，其他地理种群单倍型多样性较高(表3)；大陆沿岸各种群的核苷酸多样性、遗传距离与序列变异整体高于岛屿种群，总体核苷酸多样性差异较大(π大陆沿岸=0.197 8、π岛屿=0.107 3)，扩张系数与种群遗传多样性表明种群存在近期扩张(总体CI 95%：Tau=426.109、θ0=0、θ1=inf)；大陆沿岸种群与岛屿种群Fst=0.099(P<0.001)；AMOVA分析可知变异主要来自种群内(66.93%，P<0.001)(表4)。

表1 15个矮拟帽贝种群形态数据的平均数与标准差Table 1 Mean values and standard deviation of the morphological data of 15 Patelloida pygmaea populations

表2 10个矮拟帽贝地理种群近30年环境因子Table 2 Environmental factors of 10 Patelloida pygmaea populations in recent 30 years

表3 15个矮拟帽贝地理种群遗传多样性Table 3 Genetic diversity of 15 Patelloida pygmaea populations

表4 矮拟帽贝分子方差分析Table 4 AMOVA analysis of Patelloida pygmaea

表5 矮拟帽贝遗传多样性、形态与环境因子相关性分析Table 5 Relationship correlation of genetic diversity and morphology of Patelloida pygmaea with environmental factors

注Note:*P<0.05.

2.2遗传多样性与环境因子间的相关性

所有环境因子经正态分布检验，结果显示均为正态分布，通过Person相关性检验(表5)，发现矮拟帽贝种群遗传距离、序列变异、核苷酸多样性与日照时数、日照百分率呈显著正相关(P<0.05)；其他数据间无显著相关性(P>0.05)。主成分分析显示，壳长、壳宽为变异关键因子(贡献值：74.6%、18.4%)，壳长、壳宽与环境因子相关性分析结果表明：壳长、壳宽与年均降水量呈显著负相关(r=-0.703/-0.728，P=0.023/0.017)。

3 讨论

对矮拟帽贝15个种群的线粒体16S rRNA基因序列研究表明，各地理种群较高的单倍型多样性(0.486～1.000)和较低的核苷酸多样性(0.001 6～0.185 5)，符合Grant和Bowen(1998)提出的海洋类生物较高的单倍型多样性和较低的核苷酸多样性模式，该类型种群可能经历过历史扩张事件(Bowenetal.，2001)，Tau值与θ值也表明该物种存在近期历史扩张；大陆沿岸种群与岛屿种群遗传差异显著，大陆沿岸种群遗传多样性整体高于岛屿种群，且分化显著，暗示存在显著的陆岛遗传模式，大陆沿岸种群更易扩张与进化，推测大陆沿岸种群更为古老，岛屿种群基因频率受陆地基因流影响显著(Finston & Peck，1995；曲若竹等，2004)。

气候因子是影响生物进化的重要因子，其中温度和降水量则是影响陆生动物进化的2个最重要生态因子，对外温动物与内温动物均具有显著影响，且气候对不同地区的同类物种的影响程度也不同。渤海不同地理种群矮拟帽贝遗传结构研究结果表明，16S rRNA序列的遗传距离、序列变异、核苷酸多样性与日照时数、日照百分率呈显著正相关，说明日照是影响矮拟帽贝遗传多样性的重要因素，日照时数越大，遗传多样性越高。温度是影响贝类遗传变异的关键环境因子(Kimetal.，2009)，海表温度作为海水热力状况的物理量，受太阳辐射影响显著(Joos & Spahni，2008)，推测矮拟帽贝作为潮间带生物，栖息地以岩石海床为主，日照改变海面温度继而影响种群遗传多样性，且环境越不稳定，遗传多样性越低，与海洋无脊椎动物生态研究结果基本一致(Monarietal.，2007；Yuetal.，2009)。

形态变异主成分(壳长、壳宽)与环境因子相关性分析表明，降水量是影响形态变异的主要因素，与遗传存在差异的主要原因是自然选择导致物种基因遗传与形态变异的模式不同，即选择压力对不同的类型作用不同(Correetal.，2002)。形态量度和其他环境因子的相关性趋势与遗传分析一致，说明表型通常是受基因控制和影响，自然选择对表型作用实际上是作用于基因或基因型，即环境因素可以通过影响个体表型进而影响基因型频率(Huangetal.，2007)。稳定环境中，较小的自然选择压力使不同的基因型被最大限度地保留，因此遗传多样性较高；在气候不稳定区域，种群保留适应多变环境的基因型而在选择压力下淘汰其他基因型，遗传多样性会降低。

曲江勇, 刘迺发. 2012. 环境因子对雉鸡甘肃亚种种群遗传多样性的影响[J]. 四川动物, 31(4): 518-523.

曲若竹, 侯林, 吕红丽, 等. 2004. 群体遗传结构中的基因流[J]. 遗传, 26(3): 377-382.

Arbuckle JL. 2010. IBM SPSS Amos 19 user’s guide[M]. Chicago: Amos Development Corporation.

Bowen BW, Bass AL, Rocha LA,etal. 2001. Phylogeography of the trumpetfishes(Aulostomus): ring species complex on a global scale[J]. Evolution, 55(5): 1029-1039.

Canapa A, Barucca M, Caputo V,etal. 2000. A molecular analysis of the systematics of three Antarctic bivalves[J]. Bolletino di Zoologia, 67(sup1): 127-132.

Corre VL, Roux F, Reboud X. 2002. DNA polymorphism at the FRIGIDA gene inArabidopsisthaliana: extensive nonsynonymous variation is consistent with local selection for flowering time[J]. Molecular Biology & Evolution, 19(8): 1261-1271.

Excoffier L, Lischer HEL. 2010. Arlequin suite ver 3.5: a new series of programs to perform population genetics analyses under Linux and Windows[J]. Molecular Ecology Resources, 10(3): 564-567.

Finston TL, Peck SB. 1995. Population structure and gene flow in Stomion: a species swarm of flightless beetles of the Galápagos Islands[J]. Heredity An International Journal of Genetics, 75(4): 390-397.

Folmer O, Black M, Hoeh W,etal. 1994. DNA primers for amplification of mitochondrial cytochrome c oxidase subunit I from diverse metazoan invertebrates[J]. Molecular Marine Biology & Biotechnology, 3(5): 294-299.

Gonzálezwevar CA, Nakano T, Caete JI,etal. 2010. Molecular phylogeny and historical biogeography ofNacella(Patellogastropoda: Nacellidae) in the Southern Ocean[J]. Molecular Phylogenetics & Evolution, 56(1): 115-124.

Grant W, Bowen BW. 1998. Shallow population histories in deep evolutionary lineages of marine fishes: insights from sardines and anchovies and lessons for conservation[J]. Journal of Heredity, 89(5): 415-426.

Honda M, Yasukawa Y, Ren H,etal. 2002. Phylogenetic relationships of the Asian box turtles of the genusCuorasensu lato (Reptilia: Bataguridae) inferred from mitochondrial DNA sequences[J]. Zoologicalence, 19(11): 1305-1312.

Huang Z, Liu N, Luo S,etal. 2007. Ecological genetics of rusty-necklaced partridge (Alectorismagna): environmental factors and population genetic variability correlations[J]. Korean Journal of Genetics, 29(2): 115-120.

Jin YT, Liu NF. 2008. Ecological genetics ofPhrynocephalusvlangaliion the North Tibetan (Qinghai) Plateau: correlation between environmental factors and population genetic variability[J]. Biochemical Genetics, 46(9-10): 598-604.

Joos F, Spahni R. 2008. Rates of change in natural and anthropogenic radiative forcing over the past 20,000 years[J]. Proceedings of the National Academy of Sciences of the United States of America, 105(5): 1425-1430.

Kim M, Ahn IY, Cheon J,etal. 2009. Molecular cloning and thermal stress-induced expression of a pi-class glutathione S-transferase (GST) in the Antarctic bivalveLaternulaelliptica[J]. Comparative Biochemistry & Physiology Part A Molecular & Integrative Physiology, 152(2): 207-213.

Kirkendale LA, Meyer CP. 2004. Phylogeography of thePatelloidaprofunda group (Gastropoda: Lottiidae): diversification in a dispersal-driven marine system[J]. Molecular Ecology, 13(9): 2749-2762.

Kumar S, Stecher G, Tamura K. 2016. MEGA7: molecular evolutionary genetics analysis version 7.0 for bigger datasets.[J]. Molecular Biology & Evolution, 33(7): 1870.

Librado P, Rozas J. 2009. DnaSP v5[J]. Bioinformatics, 25(11): 1451-1452.

Lindberg DR. 1990. Systematics ofPotamacmaeafluviatilis(Blanford): a brackish water patellogastropod (Patelloidinae: Lottiidae)[J]. Journal of Molluscan Studies, 56(2): 309-316.

Michod RE. 1981. Ecological genetics[J]. Science, 212(4491): 154-155.

Monari M, Matozzo V, Foschi J,etal. 2007. Effects of high temperatures on functional responses of haemocytes in the clamChameleagallina[J]. Fish & Shellfish Immunology, 22(1-2): 98-114.

Nakai S, Miura O, Maki M,etal. 2006. Morphological and habitat divergence in the intertidal limpetPatelloidapygmaea[J]. Marine Biology, 149(3): 515-523.

Nakano T, Ozawa T. 2007. Worldwide phylogeography of limpets of the order Patellogastropoda: molecular, morphological and palaeontological evidence[J]. Journal of Molluscan Studies, 73(1): 79-99.

Nevo E. 1981. Genetic variation and climatic selection in the lizardAgamastellioin Israel and Sinai[J]. Theoretical and Applied Genetics, 60(6): 369-380.

Ni G, Li Q, Kong L,etal. 2012. Phylogeography of the bivalveTegillarcagranosain coastal China: implications for management and conservation[J]. Marine Ecology Progress, 452: 119-130.

Scheiner SM. 1993. Genetics and evolution of phenotypic plasticity[J]. Annual Review of Ecology & Systematics, 24(1): 35-68.

Slatkin M. 1987. Gene flow and the geographic structure of natural populations[J]. Science, 236(4803): 787.

Surhone LM, Tennoe MT, Henssonow SF. 2010.Patelloidapygmaea[M]. Saarbrücken: Betascript Publishing.

Takenori S, Takashi O. 1993. Anatomy and systematic position ofYayoiacmea, a new genus for Japanese tiny limpet “Collisella” oyamai Habe, 1955 (Gastropoda: Lottiidae)[J]. The Japanese Journal of Malacology, 52: 193-209.

Wayne, Robert K. 2004. Molecular markers, natural history and evolution[M]. Stamford: Sinauer Associates: 25-86.

Yu JH, Song JH, Choi MC,etal. 2009. Effects of water temperature change on immune function in surf clams,Mactraveneriformis(Bivalvia: Mactridae)[J]. Journal of Invertebrate Pathology, 102(1): 30-35.

PopulationEcologicalGeneticsofPatelloidapygmaeaBasedonMitochondrial16SrRNAFragmentSequencesintheBohaiSea

SHI Boyang1, LI Yuanyuan1, WANG Ruxiao1, QU Jiangyong1, 2*, GUO Chenghua1*, LIU Chuanlin1, HOU Jianhai1, FANG Hui1, LIU Chutong1

(1. College of Life Sciences, Yantai University, Yantai, Shandong Province 264005, China; 2. Hainan Key Laboratory for Herpetological Research, Hainan Tropical Ocean University, Sanya, Hainan Province 572000, China)

The relationship between genetic diversity and environmental factors of 15Patelloidapygmaeapopulations in the Bohai Sea was studied by PCR and DNA sequencing. A total of 264 polymorphic loci were found in the 16S rRNA partial sequence of the 180 individuals, and 153 transitions as well as 114 transversions were defined by 32 haplotypes. The base contents of C, T, A and G were 18.12%, 27.74%, 33.07%, and 21.07%, respectively. The haplotype diversity was 0.835 1±0.019 1, and the nucleotide diversity was 0.449±0.214. Moreover, the nucleotide diversity ofP.pygmaeapopulation in the coast (0.191 8±0.091 7) was higher than that of island (0.107 3±0.081 6), and the results of population demographic history indicated that the populations ofP.pygmaeamay have experienced a recent demographic expansion. The correlation analysis between environmental factors and genetic factors showed that the genetic distance, sequence variation and nucleotide diversity ofP.pygmaeawere significantly positive correlated with sunshine duration and sunshine percentage (P<0.05). The result of morphological correlation analysis showed that the main components (length and width) were significantly negative correlated with the annual average precipitation (P<0.05). Average annual precipitation, sunshine duration and sunshine percentage were the dominant factors ofP.pygmaeamorphological variation, and the population genetic diversity will be higher in the the area with high environmental stablility.

Patelloidapygmaea; 16S rRNA; ecological genetical; genetic diversity

2017-01-12接受日期：2017-03-28

国家自然科学基金地区基金项目(31460562)； 烟台大学博士科研启动经费(SM15B01)； 2016年国家级大学生创新创业训练计划项目(201611066009)； 烟台大学研究生科技创新基金项目； 烟台大学开放实验室项目

史博洋, 硕士研究生, 主要从事海洋生物学研究, E-mail:sby_hope@163.com

*通信作者Corresponding author, E-mail:qjy@ytu.edu.cn; gch@ytu.edu.cn

10.11984/j.issn.1000-7083.20170013

Q145； Q346

： A

： 1000-7083(2017)04-0386-06