海洋生态与环境科学重点实验室知识库

探究生物不同地理群体的遗传结构，不仅可以加深对生物多样性的认识，也是合理有效管理生物资源的基础。目前对于日本鳗鲡（Anguilla japonica）群体遗传结构的研究大多基于中性遗传标记，本研究采用与基因相关的微卫星标记和基于全基因组扫描筛选的离散SNP标记开展研究，从新的角度解析日本鳗鲡的群体遗传结构。除了群体遗传结构的解析，适应性进化也是近几十年来进化生态学研究的热点。本地适应性是指本地群体的某些特征，使其相较于外来群体可更好地适应本地环境。一直以来都普遍认为，适应性进化是一个缓慢的过程，需要经过很长时间的积累，因此它可能不会对当前的生态过程造成影响。然而近年来发现，适应性进化可以发生在几个世代的时间尺度内。鳗鲡属鱼类随机交配的特征使它们成为研究微时间尺度上本地适应性问题理想的材料。欧洲鳗鲡（Anguilla anguilla）、美洲鳗鲡（Anguilla rostrata）的研究均已证明，尽管鳗鲡是随机交配物种，依然可以检测到其对本地环境适应的遗传学信号。对日本鳗鲡的本地适应性进行研究，将加深对于单一世代时间尺度内适应性进化的认识，也有助于深入理解全球气候变化对海洋鱼类的影响。从二十世纪八十年代开始，日本鳗鲡的种群衰退不断加剧。玻璃鳗（Glass eel）、黄鳗（Yellow eel），以及银鳗（Silver eel）的野生群体资源量均出现了严重衰退。在2016年公布的IUCN濒危物种红色名录中，日本鳗鲡已经被列为“濒危”。目前对于日本鳗鲡种群资源的保护已经迫在眉睫。解析日本鳗鲡的群体遗传结构将为其种群资源的管理和保护提供科学依据，而阐明日本鳗鲡对本地环境的适应机制有助于深入理解日本鳗鲡种群衰退的机制。
本研究的主要研究内容和结果如下：
1.      开发与基因相关的微卫星标记
我们开发了日本鳗鲡种群中与基因相关的微卫星标记。结果，总共有24个位点具有多态性，且它们在所有的地理群体都呈现出多态性。在这24个位点中，有18个是根据欧洲鳗鲡和美洲鳗鲡中受自然选择作用的基因开发，另外6个来自于随机挑选的转录本序列。每个位点的平均等位基因数为4.846到29.231。每个群体的平均观测杂合度为0.675到0.716，平均期望杂合度为0.746到0.769。经过Bonferroni校正 （α= 0.05, K = 13），群体和位点的312个组合中有29个显著偏离哈温平衡。而在连锁不平衡检验中，经过Bonferroni校正（α= 0.05, K = 276）所有的位点间均不存在连锁不平衡。在无效等位基因的检验中，有8个位点在大于3个群体（未在所有群体）中表现出无效等位基因的特征。
2.      开发SNP标记
选取两个群体各12尾玻璃鳗幼体进行酶切位点相关DNA测序（Restriction site-associated DNA sequencing，RAD-Seq）。将过滤后的序列根据个体特异的标签序列进行归类。使用BWA version 0.5.9将序列比对到日本鳗鲡的基因组草图上，使用SAMTOOLS检测到73,557个SNP。使用ARLEQUIN 3.5筛选可能受到自然选择作用的位点，共检测到250个离散SNP标记，在这其中有85个SNP标记可以成功进行功能注释。使用KASPar（KBiosciences Competitive Allele-Specific PCR genotyping system）技术对来自12个地理群体的373个样品在96个离散SNP位点（85个可以进行功能注释的离散SNP位点+随机挑选的另外11个未能进行功能注释的离散SNP位点）上进行基因分型。
3.      基于与基因相关的微卫星标记探究日本鳗鲡的群体遗传结构
基于所有群体计算的F_ST值为-0.001，且不显著（95% CI: -0.002 < F_ST < 0.000），这表明不同地理群体间并不存在显著的遗传分化。所有的两两群体间F_ST值经过FDR校正（α = 0.05, K = 78）后均不显著。STRUCTURE分析的最优K值为5。然而，在K = 5时，所有的群体呈现出类似的遗传组成，这也暗示不同地理群体间并不存在遗传结构。Mantel test的结果表明，13个地理群体间并不存在距离隔离模式（R² = 4.634e-03, P = 0.6233）。所有的结果均支持日本鳗鲡是随机交配物种。
4.      基于SNP标记探究日本鳗鲡的本地适应性
基于SNP标记，在12个地理群体的两两比较中，出现了部分显著的F_ST值，表明这12个地理群体间存在一定程度的适应性分化。基于SNP位点的STRUCTURE结果表明，中国台湾群体（MA、MY）和日本群体（MJ）的遗传组成类似，而中国大陆各群体之间的遗传组成类似。该结果也得到了DAPC结果的验证。
基于BAYENV2和SAMBADA两种软件筛选到受自然选择的位点相同，即AVPY01361442_1313和AVPY01251147_815。AVPY01361442_1313位点与SMC6基因相关，而AVPY01251147_815位点与taar13c基因相关。SMC6基因可能参与了玻璃鳗幼体对剧烈环境扰动的适应，通过DNA修复维持自身染色体的稳定。而taar13c基因很可能介导了日本鳗鲡对尸胺等二元胺的行为反应；然而，若想彻底阐明taar13c基因参与何种行为反应，尚需要开展玻璃鳗幼体对气味的行为学研究。
本研究验证了日本鳗鲡随机交配的假说，并检测到了其发生本地适应性的信号。基于离散SNP标记，发现了岛源群体（采自岛屿河口的群体）和陆源群体（采自大陆河口的群体）适应性分化的信号。研究结果可以为日本鳗鲡种群资源的保护和管理提供基础资料。