韩国的Eleocharis ussuriensis Zinserl. 遗传变异与空间遗传结构：对生态监测与资源管理的启示

《Genes》：Genetic Variation and Spatial Genetic Structure of Eleocharis ussuriensis Zinserl. in South Korea: Implications for Ecological Monitoring and Resource Management Eun-Hye Kim, Kang-Rae Kim, Mi-Hwa Lee, Jaeduk Goh and Jeong-Nam Yu

　　

背景/目标：Eleocharis ussuriensis Zinserl. 是广泛分布于东北亚的多年生河流莎草，也是韩国淡水植被中的主要组成部分。然而，该物种在不同水文生境下的种内遗传结构仍不明确。本研究旨在基于全基因组数据开发新的简单重复序列(SSR)标记，并调查韩国E. ussuriensis 种群的遗传变异和群体结构。方法：研究人员从全基因组序列数据中开发了21个新的微卫星(SSR)标记，并将其应用于分析来自6个种群的120个个体。通过等位基因频率指标评估遗传多样性、分化和基因流，并利用地理坐标获得的空间信息进一步评估群体遗传结构。结果：共检测到201个等位基因，平均多态性信息含量(PIC)为0.759，表明标记具有较高信息性。种群的平均遗传多样性显示观察杂合度(HO=0.360)和期望杂合度(HE=0.281)均处于中等水平，多位点基因型比例(G/N)在种群间介于0.30至1.00之间。遗传分化显著(FST=0.373–0.669; Jost’s D=0.540–0.997)。Mantel检验显示遗传分化与地理距离显著相关(r=0.67, p<0.001)。基于等位基因频率和空间明确的方法均表明种群间存在遗传结构。结论：结果显示出种群间遗传结构的空间趋势，反映了等位基因在各区域的分布模式。这些发现为E. ussuriensis 的遗传变异提供了基础信息，并可能有助于更好地理解其生态动态。

1. 研究背景与目的

河流生态系统在气候变化驱动的快速环境变化中日益显现其关键生态功能。从景观生态学角度看，河流系统作为网状结构，通过水文驱动的扩散路径连接个体与基因。水文扰动（如洪水）可能通过影响扩散、定居和存活过程，导致河岸植被形成异质性遗传结构。河岸植物与河流动态紧密相关，受水质、水位波动、流速等环境变量直接影响。在频繁受干扰或环境不稳定的生境中，有性繁殖机会可能减少，有限基因型的重复克隆扩张可降低种群水平的遗传多样性。尽管针对河岸系统的生态功能与物种多样性已有大量研究，但对分布广泛、非受威胁植物的种内遗传变异关注相对较少，群体遗传研究多聚焦于濒危或稀有类群。

Eleocharis ussuriensis Zinserl. 是东北亚湿地（含韩国、中国、日本）的多年生河岸物种，属于莎草科(Cyperaceae)，风媒传粉，通常生长在静水边缘或水库沿岸。其通过根茎生长形成密集克隆群落，有助于沉积物稳定和河岸生态系统结构维持。其提取物因具有皮肤美白作用在韩国被注册用于化妆品。然而，尽管其可能具有较高的扩散潜力，但仍不清楚跨水文景观的有效基因流是否得以维持，以及这种潜在扩散是否会转化为不同水文生境种群间稳定的遗传连通性。

莎草科分类学上因形态相似、诊断特征有限而具有挑战性。在莎草科中，Eleocharis 属由于形态简单、关键性状变异较大，在物种内和物种间存在显著变异。针对这些分类学局限，需要可靠的分子工具。虽然已在相关物种（如E. dulcis）中开展基于微卫星(SSR)的分析，且有研究表明微卫星标记可在近缘种间转移，但其适用性有限且依赖背景。即使在同属内，跨物种扩增的成功率也常为中等水平（约60%），且多态性水平在种间差异较大，降低了标记的可靠性与信息性。迄今尚未开发出E. ussuriensis 的物种特异性微卫星标记，因此开发物种特异性标记对获取准确、稳健的遗传信息至关重要。

基于此，本研究旨在从全基因组序列数据中开发针对E. ussuriensis 的物种特异性微卫星(SSR)标记，以提高标记稳定性，克服有限基因组区域衍生标记的局限性。利用这些SSR标记，评估遗传多样性、种群分化与空间遗传结构，建立种群鉴别的遗传框架，为该物种的后期管理与利用提供遗传基础信息。

2. 关键技术方法

研究人员于2023年5月在韩国6个河岸种群中采集了120株E. ussuriensis 个体样本（每种群20株）。采样点选自独立流域系统与不同生境，包括4个流水环境（YC、HC、GS、US）和2个静水环境（JC、BS）种群。样本采集时个体间距>5 m，以避免采集克隆生长或有限花粉扩散产生的遗传相关个体。

从组装的全基因组数据中鉴定SSR位点，使用MISA工具，设置最小重复次数阈值4-10。经Primer3生成初始引物对，通过SnapGenev7.1.0验证结合特异性，在每种群6个个体上筛选引物，保留清晰、可重复扩增的位点。最终验证的21个微卫星标记被提交至GenBank。

采用M13-tailed法进行PCR扩增，产物经GeneScan? 500 ROX和Hi-Di? 甲酰胺混合，变性后在ABI 3730xl DNA分析仪上分析，用GeneMarker? v2.6.7进行等位基因分型。通过MICROCHECKER v2.2.3检测可能的基因分型错误。用PopGene v1.32计算遗传多样性指标，Cervus v3.0.7计算PIC。用R包UpSetR v1.4.0可视化等位基因频率数据。种群分化（F-统计量）用hierfstat v0.5-11和FSTAT v2.9.4估算，用BOTTLENECK v1.2.02检测近期瓶颈。遗传分化与基因流通过AMOVA、Jost’s D、PCoA、Mantel检验评估。用STRUCTURE v2.3.4和TESS3r分析遗传结构，通过交叉熵准则确定最优K值。

3. 研究结果

3.1. 微卫星位点变异

21个微卫星位点共检测到201个等位基因，每个位点等位基因数4-18，平均9.6。PIC值范围0.465-0.882，平均0.759，显示高度多态性。位点特异性F_IS值范围-0.935-0.683，F_ST值范围0.339-0.887，平均0.639，表明种群间遗传分化强。

3.2. 遗传多样性与瓶颈

6个种群中观测等位基因数范围33-70，其中148个为种群特异性等位基因。多态位点比例(P_0.95)范围33.3-76.2。有效等位基因数/位点(A_e/L)范围1.4-2.2。观察杂合度(H_O)范围0.155-0.645，期望杂合度(H_E)范围0.163-0.432。多位点基因型比例(G/N)范围0.30-1.00。平均H_O=0.362，平均H_E=0.284，平均A_e/L=1.7，平均G/N=0.80。

近交系数(F_IS)范围-0.840-0.076，多数种群呈现杂合子过量。HC种群在所有三种突变模型下均表现强烈瓶颈信号。GS、US、BS种群在IAM和TPM下检测到显著杂合子过量。YC、JC种群未检测到显著瓶颈信号。

3.3. 遗传分化与基因流

AMOVA显示62%的遗传变异存在于种群间。配对F_ST值范围0.373-0.669，Jost’s D值范围0.540-0.997。Mantel检验显示遗传分化与地理距离显著相关(r=0.67, p<0.001)。

PCoA前两轴解释总遗传变异的50.2%，前三维累计解释62.1%。YC、HC、JC聚为一组，US、BS在左侧，GS明显分离。方向性迁移分析表明YC、HC、JC间基因流较强（YC与HC间迁移系数=1.00），US、BS间存在中度相互迁移，GS种群基因流较低。

3.4. 空间与遗传结构

STRUCTURE和TESS3r分析均支持最优K=6。个体主要被分配到对应采样种群的遗传簇，显示明确的种群水平遗传结构。

4. 讨论与结论

讨论 所开发的多个微卫星位点呈现高多态性，但部分位点F_IS值极端，偏离Hardy–Weinberg平衡，这可能与莎草科植物有限的传粉和种子扩散、克隆繁殖及区域特异性种群分化有关。E. ussuriensis 的遗传多样性（H_O=0.362, H_E=0.284）低于相近物种E. parvula，与莎草科其他物种如Carex scabrifolia、C. kobomugi相当或略低。直接比较遗传多样性时需谨慎，因标记系统与分析方法的差异。

不同种群遗传多样性差异可能与生境环境相关。四个流水种群中，HC、GS、US遗传多样性低于YC，多位点基因型比例(0.30-0.40)低于总体均值，且检测到瓶颈信号，表明这些种群可能经历近期群体波动或干扰。这些种群位于易受夏季季风洪水影响的河岸环境，水文干扰与人为河道改造可能影响种群动态并促进克隆繁殖。YC表现出高的多位点基因型比例与均衡杂合度，而HC呈现明显的杂合子过量与瓶颈信号，显示其近期群体不稳定。JC种群位于较老的水库，受管理缓冲，未检测到明确瓶颈信号，但遗传多样性相对较低，反映有限的种群扩张。BS种群遗传多样性较高，位于受保护的饮用水水库，稳定生境与独特水文特征可能有助于维持和积累遗传变异。

AMOVA结果显示62%的遗传变异存在于种群间，配对F_ST(0.373-0.669)和Jost’s D(0.540-0.997)表明分化显著。显著的距离隔离(IBD)模式显示地理距离推动遗传分化，但分化程度超出单纯距离影响，反映有限的等位基因共享和种群特异性等位基因存在。方向性迁移分析显示YC、HC、JC间基因流强，US与BS间存在中度相互迁移，GS则呈现基因流较低，表明基因流空间异质性强。

遗传分化与有限基因流在等位基因频率和空间聚类分析中一致体现，每个种群形成明显遗传簇。E. ussuriensis 的斑块分布与克隆生长型可能促进此模式，有限的有效扩散与局部繁殖可增强种群间遗传分化。显著IBD模式支持地理距离与局地过程共同塑造其空间遗传结构。

结论 研究人员开发了21个新的基因组衍生的SSR标记，用于评估韩国6个E. ussuriensis 种群的遗传变异、分化和空间结构。这些标记呈现足够多态性和可重复性，为评估该物种群体水平遗传模式提供了可靠工具。

遗传与克隆多样性在种群间存在差异，显示生境背景与局部种群历史影响遗传组成。种群间存在遗传分化，并伴随显著的距离隔离(IBD)模式。尽管位点间无直接水文连通性，基因流与地理距离及位点特异性群体过程相关。这些模式反映在空间遗传结构中，显示种群形成明显群体的趋势。这些发现表明E. ussuriensis 区域种群呈现异质性遗传特征。

总体而言，本研究为E. ussuriensis 的遗传变异提供了初步见解，有助于更好地理解其生态动态。开发的SSR标记可促进未来比较研究，为基于证据的保护与资源管理策略提供信息。