巴西西南海岸广盐性多毛类Laeonereis acuta（环节动物门：沙蚕科）微卫星标记的开发与表征

《Biochemical Genetics》：Development and Characterization of Microsatellite Markers for the Euryhaline Polychaete Laeonereis acuta (Annelida: Nereididae) in the Southwestern Brazilian Coast

　　

Laeonereis acuta是一种多毛类（polychaete）物种，通常在河口与海岸潟湖环境中具有较高丰度。由于该物种与污染栖息地密切相关，因此常被用于生态毒理学（ecotoxicology）研究。此外，其分布于空间上不连续且环境变异性较高的生境，使其成为研究局部适应（local adaptation）、遗传景观（genetic landscape）以及物种形成早期阶段（early stages of speciation）的适宜模型。本研究旨在为采自巴西西南部3个海岸潟湖的L. acuta开发引物并对微卫星标记进行表征。基于40个个体的基因分型（genotyping），共表征了10个位点（loci）。各位点等位基因（alleles）数范围为2至19。在5个位点检测到空等位基因（null alleles）证据，但当分别对各海岸潟湖进行分析时，其频率有所降低。当将全部个体视为单一种群时，5个位点显示为正且显著的F_IS值，7个位点偏离哈迪-温伯格平衡（Hardy-Weinberg equilibrium，HWE）。Maricá与Guarapina在若干位点表现出杂合子过剩（heterozygote excess），而Jaconé则显示出种群遗传隔离（population genetic isolation）的证据。这10个微卫星位点均具有多态性（polymorphic），适用于L. acuta的种群遗传分析，尽管这些模式未必能够代表其他地理区域。这些标记可通过阐明对污染物的生理响应是否与种群间遗传分化相关，从而促进生态毒理学研究。此外，它们还为探究不连续环境中的遗传结构与连通性提供了有价值的工具。

本文发表于《Biochemical Genetics》，主要围绕广盐性多毛类Laeonereis acuta微卫星标记的开发与群体遗传学应用潜力展开。研究背景在于，物种地理分布及其在异质环境中的存续，既受气候变化、盐度梯度等非生物因子影响，也受生活史特征、生态互作及资源需求等生物因子制约；而个体应对这些因素的能力又与其遗传组成密切相关。自然选择与遗传漂变可降低种群内遗传多样性并增加种群间分化，突变与基因流则有助于提高多样性，其中基因流还可促进种群基因库同质化。因此，在研究河口与海岸潟湖等空间不连续环境中的局部适应、种群连通性与谱系分化早期阶段时，亟需高多态性的分子标记。

L. acuta常见于河口和海岸潟湖，丰度高，并且由于与受污染栖息地关系密切，已被广泛应用于生态毒理学研究。论文指出，不同种群的生理与生态毒理学响应可能因局部适应及生活史性状差异而不同，而这些差异通常反映种群遗传变异。因此，明确其种群遗传结构和连通性，不仅有助于理解等位基因共享、遗传断裂与早期物种形成过程，也有助于判断污染响应是否与遗传分化相关。然而，此前L. acuta缺乏可直接用于精细尺度群体遗传研究的微卫星工具，这构成了本研究开展的直接动因。

研究人员首先针对L. acuta开发微卫星引物，并在巴西东南部里约热内卢州的3个海岸潟湖种群中进行验证与表征，包括Maricá、Guarapina与Jaconé。研究结果表明，共获得10个可用微卫星位点，且多数位点具有良好多态性，能够有效揭示细尺度种群遗传分化。总体分析显示，多个位点存在空等位基因信号、F_IS升高以及偏离哈迪-温伯格平衡的现象；但在分潟湖分析时，这些现象有所缓解，提示总体样本合并分析可能受种群结构影响。进一步比较发现，Maricá与Guarapina之间未检测到显著遗传分化，而Jaconé与前两者之间存在显著分化，并表现出更高近交水平，提示其种群更为封闭。由此，论文证明所开发的微卫星位点可用于L. acuta种群遗传结构、连通性及生态毒理学问题的研究。

研究的重要意义主要体现在三个方面。其一，为L. acuta这一河口—潟湖优势类群建立了可重复使用、成本相对较低的分子标记体系。其二，这些标记能够识别空间不连续环境中的当代遗传结构与细尺度连通性，从而补充线粒体标记在时间尺度上的局限。其三，该工具可拓展至生态毒理学场景，用于评估受遗传隔离的种群是否对环境胁迫因子表现出差异化生理响应。

在技术方法方面，研究人员从Maricá潟湖1个个体提取DNA，采用GS FLX 454（Roche）高通量测序（HTS，高通量测序）筛选微卫星序列，并利用QDD 3.1进行位点识别与引物设计。随后在巴西东南部3个海岸潟湖来源的40个个体样本中进行PCR扩增与基因分型，采用M13尾引物荧光标记策略降低分型成本，并在SeqStudio测序平台上完成片段分析。数据分析包括空等位基因与分型误差检测、连锁不平衡分析、等位基因数与稀释等位基因丰富度估计、杂合度、F_IS、哈迪-温伯格平衡（HWE）、成对F_ST计算，以及基于最大似然法的亲缘关系推断。

以下结合论文结果与讨论部分的小标题内容进行概括性解读。

Results and Discussion

首先，在全部40个完成基因分型的个体中，研究人员在10个微卫星位点共检测到74个等位基因，整体缺失数据比例为12.5%。除Lae_5存在口吃型分型误差（stuttering）证据外，其余位点未发现大片段等位基因脱扣（large-allele dropout）或明显评分错误。5个位点检测到潜在空等位基因，但在按各海岸潟湖分别分析时，其频率下降。这一结果表明，某些偏离可能与样本合并后的Wahlund效应有关，而不完全是位点本身质量问题。与此同时，位点间未检测到显著连锁不平衡（linkage disequilibrium），说明这些位点总体上可作为相对独立的遗传标记用于后续分析。

其次，在位点多态性方面，各位点等位基因数介于2至19之间，平均为7.4，显示出较高的遗传变异水平。其中二核苷酸重复位点的等位基因数较低，三核苷酸位点居中，而四核苷酸和五核苷酸位点均表现出10个及以上等位基因。论文将这一范围与其他多毛类及若干无脊椎动物中已报道的微卫星多态性水平进行对照，认为其总体处于合理范围内。这说明本研究筛选出的位点在L. acuta中具备足够分辨率，可满足种群遗传分析需求。

再次，在总体遗传平衡特征方面，除Lae_06和Lae_14外，其余位点的观察杂合度（observed heterozygosity）均低于期望杂合度（expected heterozygosity）。7个位点显著偏离HWE，5个位点呈现显著且为正的F_IS值，总体F_IS为0.178（p<0.0001）。研究人员指出，这种模式可能来源于非随机交配，例如近交，也可能反映样本中包含遗传结构化种群；同时，小样本量也可能影响HWE检验的统计特性。结合后续分群分析，论文更倾向于认为总体样本中存在明显种群结构，从而导致合并分析时的杂合子缺失与平衡偏离。

在种群间遗传分化方面，成对F_ST结果显示Maricá与Guarapina之间无显著遗传分化，而Jaconé与Maricá、Guarapina之间均存在显著分化。这说明前两者可能维持较高基因流，而Jaconé则相对孤立。进一步按潟湖分别分析后，每个潟湖中仅少数位点偏离HWE，且空等位基因主要集中于Jaconé与Guarapina的个别位点。共享等位基因数量总体多于私有等位基因，但Jaconé与Maricá在部分位点拥有较多私有等位基因，提示局部遗传分化已开始形成。该结果支持了海岸潟湖空间不连续性能够促进细尺度遗传结构形成的观点。

在各潟湖内部遗传组成方面，Maricá和Guarapina的多个位点表现出杂合子过剩，且总体F_IS分别为?0.105和?0.009，均不显著；而Jaconé则显示更高近交水平，总体F_IS为0.309（p<0.001）。论文将此与各潟湖的地貌水文特征联系起来：Maricá与Guarapina同属一个海岸潟湖系统，并由水道相连，可能促进基因流；Jaconé则依赖狭窄人工水道与另一受阻潟湖相连，水体交换受限、封闭性更高，因此其较高近交水平与地貌背景一致。也就是说，微卫星标记不仅揭示了种群间差异，也能够捕捉局部水文连通性对遗传结构的影响。

关于遗传多样性，3个潟湖检测到的总等位基因数分别为Maricá 46个、Guarapina 38个、Jaconé 47个。稀释等位基因丰富度在5个位点上以Jaconé更高，在4个位点上以Maricá更高。研究原先根据系统开放程度推测，较封闭潟湖的遗传多样性可能较低，而较开放的Guarapina可能因移民潜力较高而拥有更高多样性，但结果并未支持这一简单预期。由此可见，L. acuta的遗传多样性格局并不完全由潟湖开放度决定。

为解释Guarapina中较低遗传多样性，研究人员进一步开展了种群内亲缘关系分析。结果显示，Guarapina中半同胞（half-siblings）、全同胞（full-siblings）和亲子（parent-offspring）关系比例高于另两个潟湖，说明其采样个体之间存在更明显的亲缘聚集。这种亲缘结构可能压低局地遗传多样性估计。种群间亲缘比较则显示，大多数个体对属于无关关系，尤其是涉及Jaconé的跨种群比较中，无关配对占绝对多数，进一步支持Jaconé与其他潟湖之间连通性较低。

论文还将这些结果与生活史特征及既往分子研究进行对照。考虑到近缘种L. culveri具有在管中孵卵和直接发育（direct development）的特征，L. acuta预期也具有较低扩散潜力，并在小空间尺度上形成较高亲缘性。因此，即使潟湖与海水连通，也未必意味着L. acuta种群会发生显著移民输入。研究人员指出，本研究基于微卫星得到的当代细尺度遗传结构模式，与同一区域线粒体COI数据所反映的格局存在差异；这种差异很可能源于两类标记的时间分辨率不同：微卫星更敏感于近期和局地过程，而线粒体标记整合的是更长演化时间尺度上的信号。

讨论部分总体强调，所开发的10个微卫星位点虽然在总体样本分析中显示出一定程度的空等位基因、HWE偏离和F_IS升高，但分群分析表明这些问题在局地尺度上有所减弱，且位点总体多态性良好、无显著连锁不平衡，因此仍然适合作为L. acuta种群遗传学研究工具。研究同时指出，本研究观察到的遗传模式未必可以直接外推至其他地理区域，因此未来仍需在更广泛空间尺度上验证这些标记的适用性。总体而言，这套标记对于解析不连续环境中的遗传结构、近交、连通性及生态毒理响应的群体差异均具有应用价值。

研究结论部分可译为：本研究报道了多毛类L. acuta的10个微卫星位点的表征结果，并提供了用于扩增的引物序列。多数位点具有多态性，且被证明适用于检测细尺度种群遗传亚结构，凸显了其在生态学与进化研究中的应用价值。这些标记还可支持生态毒理学研究，使研究人员能够评估遗传隔离种群是否对环境胁迫因子表现出不同的生理响应。此外，将这些位点应用于巴西东南部以外的河口与海岸潟湖区域，有助于阐明生活于空间不连续环境中的种群之间遗传连通性的模式与水平。