中国东部广布植物芦苇（Phragmites australis）的遗传多样性与群体分化特征研究

《Frontiers in Plant Science》：Genetic diversity and differentiation characteristics of the widespread plant Phragmites australis in eastern China

　　

中文专业摘要翻译： 引言： 芦苇（Phragmites australis）是一种全球广布、生产力极高的挺水大型植物，在湿地生态系统中具有重要的生态功能。然而，中国东部地区芦苇的遗传多样性与群体分化特征尚未得到充分表征。 方法： 本研究采用基于RAD-seq（限制性位点相关DNA测序，Restriction-site Associated DNA Sequencing）的简化基因组测序技术，对中国东部10个群体采集的98个芦苇个体进行测序，以表征群体间的遗传变异、群体结构及地理分化特征。 结果： 共鉴定出730,413个高质量单核苷酸多态性位点（SNPs，Single Nucleotide Polymorphisms）。群体基因组分析显示，中国东部群体整体遗传多样性较高，且存在显著的群体结构，芦苇的广泛地理分布与多个遗传簇的形成相关。功能富集分析进一步表明，与植物–病原互作及糖代谢相关的通路显著富集，提示了与所采样群体间遗传变异相关的潜在生物学过程。在所采样群体中，龙岐湾（LQW）群体表现出最为独特的遗传特征，其次为白洋淀（BYD）和西枝湖（XZH）群体。相比之下，微山湖（WSH）和盘锦（PJ）群体的遗传分化程度较低，其基因组变异模式与参考基因组更为相似。值得注意的是，长江流域恰好与遗传结构的南北转变区域相吻合。 讨论： 总体而言，这些发现为中国东部芦苇的群体遗传结构、进化历史及地理分化提供了新见解，并为其进化动态与生态适应研究提供了有价值的基因组资源。

论文主体内容解读

一、研究背景与科学问题

芦苇（Phragmites australis）隶属于禾本科芦苇属，是一种多年生被子植物，也是全球分布最广的湿地植物之一。作为高生产力物种，芦苇在湿地生态系统中通过初级生产、养分循环、栖息地形成和污染物去除等过程发挥着重要生态作用。由于其广泛的逆境耐受性和高生物量产出，芦苇还被广泛应用于生态修复工程，包括构建人工湿地以去除水体中的重金属及其他污染物。在中国，芦苇还用于传统医药及多种工业应用，其生态功能、遗传多样性与进化特征近年来日益受到研究关注。

已有遗传学与谱系地理学研究揭示了芦苇存在显著的遗传变异与地理结构。近期一项基于全基因组RAD-seq的研究在全球分布范围内鉴定出六个主要遗传谱系，并证明多倍化、杂交和基因流共同塑造了其复杂的进化历史，其中存在一个独特的中国谱系及多个谱系间的广泛遗传混合。基于叶绿体的谱系地理分析鉴定出多个地理分布的单倍型，而基于SSR（简单序列重复标记，Simple Sequence Repeat）的研究也证实了中国芦苇群体间存在遗传变异。尽管中国多个重要湿地地区（包括黄河三角洲、渤海湾、闽江口和长江口）的芦苇群体已从生态和遗传角度得到研究，但对中国东部群体遗传变异的系统性全基因组分析仍然有限，尤其是地理格局及潜在环境因素对群体遗传分化的影响尚不清楚。

芦苇具有复杂的生活史策略，可能强烈影响其群体遗传结构。作为多年生克隆植物，芦苇可通过广泛的根状茎生长快速扩张，在湿地生境中长期持续存在并占据优势；同时，芦苇也能通过风媒传粉进行有性繁殖，产生大量轻质种子，可借助风、水和动物传播。有性繁殖与克隆繁殖的共存可能产生复杂的遗传模式，而水文连通性、生境破碎化和地理隔离等环境因素可能进一步调节群体间的基因流与遗传分化。因此，理解芦苇的遗传结构需同时考虑其繁殖策略与生态特征。

针对上述知识空白，本研究从中国东部十个代表性地区采集芦苇群体，采用RAD-seq简化基因组测序并结合参考基因组比对，表征基于SNP的遗传变异，旨在探究地理隔离如何与中国东部芦苇的遗传分化和群体结构相关联。研究人员假设地理隔离促进了群体间的遗传分化，并预测地理距离越远的群体遗传分化程度越高，符合基因流随地理距离增加而降低的预期。

二、主要技术方法

研究人员于2024年5月至9月从中国东部八个省份和一个直辖市采集了98个芦苇个体，涵盖十个自然群体，采样点包括西枝湖（XZH）、盘锦（PJ）、白洋淀（BYD）、独流减河口（DLJHK）、黄河口（HHK）、微山湖（WSH）、长江口（CJK）、闽江口（MJK）、龙岐湾（LQW）和防城港（FCG）。基因组DNA采用CTAB法提取，使用EcoRI限制性内切酶构建RAD-seq文库，并在Illumina NovaSeq X Plus平台进行双端测序。原始读段经Fastp v0.23.4质控后，使用BWA-MEM比对至芦苇染色体级参考基因组GCF_958298935.1，随后遵循GATK Best Practices流程，使用Sentieon HaplotypeCaller进行变异检测。变异位点经质量过滤和群体遗传过滤（平均测序深度≥4×且≤1000×、位点缺失率<30%、MAF≥0.05），并通过PLINK进行连锁不平衡（LD）剪枝，最终获得高质量SNP数据集。下游分析包括SnpEff功能注释与GO/KEGG富集分析、遗传多样性参数计算、基于Gephi和ArcGIS的地理遗传分化网络可视化、Mantel检验的隔离-by-距离（IBD）分析、FastTree系统发育树构建、ADMIXTURE群体结构分析、PCA主成分分析、亲缘关系分析以及PopLDdecay连锁不平衡衰减分析。

三、研究结果

3.1 简化基因组测序与变异鉴定

对98个芦苇样本进行简化基因组测序，共产生100.33 Gb测序读段，Q30率为91.57%，GC含量为46.52%。经序列聚类和变异检测，鉴定出55,054,643个标签、13,429,377个SNP和1,487,273个InDel（插入缺失变异，Insertion-Deletion）。经完整过滤流程和LD剪枝后，保留730,413个高质量SNP用于下游群体遗传分析。环形基因组图显示，24条染色体上的基因、SNP和InDel分布不均，多个基因组区域呈现较高的变异集中度。

3.2 多态性位点的功能富集分析

KEGG富集分析显示，植物–病原互作通路具有最高的Gene Ratio，其次为三个糖代谢相关通路，包括淀粉和蔗糖代谢、糖酵解/糖异生以及氨基糖和核苷酸糖代谢，这些通路均表现出较高的?log10(P值)，表明相关基因显著富集。GO富集分析中，分子功能类别以催化活性富集最强，其次为小分子结合和阴离子结合；细胞组分层面，与遗传变异相关的基因在细胞质组分中显著富集，超分子聚合物、超分子纤维和超分子复合物也显著富集；生物学过程方面，磷代谢相关过程（包括磷代谢过程、含磷酸盐化合物代谢过程和磷酸化）显著富集，碳代谢相关过程（包括碳水化合物代谢过程和有机物分解代谢过程）同样显著富集。

3.3 遗传多样性分析

群体遗传多样性参数显示，防城港（FCG）群体拥有最多的私有SNP（1409个），而微山湖（WSH）和盘锦（PJ）群体最少（分别为27个和66个）。龙岐湾（LQW）群体的观察杂合度（H_o）、期望杂合度（H_e）、多态信息含量（PIC）、Shannon信息指数（I）和最小等位基因频率（MAF）均为最低，分别为0.040、0.028、0.022、0.041和0.022，其次为白洋淀（BYD）和西枝湖（XZH）。所有十个群体的近交系数（F_IS）均为负值（范围?0.225至?0.015），表明观察杂合度普遍高于期望杂合度。长江口（CJK）群体的H_o最高（0.336），盘锦（PJ）和微山湖（WSH）群体的H_e（均为0.240）、PIC（均为0.196）和I（均为0.370）最高。核苷酸多样性（π）除LQW（0.0325）外，各群体介于0.131至0.256之间，PJ和WSH群体π值最高。群体分化分析显示，PJ与WSH之间的F_ST最低（0.052），BYD与LQW之间的F_ST最高（0.498），LQW与其他群体的F_ST均大于0.2。IBD分析显示，包含所有群体对时，遗传分化与地理距离呈弱正相关但不显著（Mantel r = 0.323，P = 0.122）；去除所有涉及LQW的群体对比较后，相关性强度显著提高，呈中等到强正相关且边缘显著（Mantel r = 0.458，P = 0.059）。

3.4 群体遗传关系与结构分析

基于过滤后SNP数据构建的近似最大似然（aML）系统发育树显示，来自北方群体（包括BYD、XZH、HHK、PJ、WSH和DLJHK）的个体大体聚为一支，而南方群体（包括CJK、FCG和MJK）形成独立分支，LQW独立于南北主要分支之外，呈现明显的南北遗传格局。三维主成分分析（3D-PCA）中，PC1、PC2和PC3分别解释11.4%、7.2%和6.4%的遗传方差，南北群体沿主成分明显分离，LQW形成独立聚类，与系统发育树结果一致。成对亲缘关系分析显示，大多数个体对属于无关或远亲类别，FCG群体内部存在大量近亲信号，而BYD、PJ和WSH群体内部亲缘关系较弱。ADMIXTURE分析中，交叉验证误差在K=9时达到最小值，表明十个群体中存在九个遗传簇。在K=8–10时，个体呈现多个遗传簇的复杂组合，PJ和WSH遗传组成较为相似，XZH包含与PJ和WSH相似的遗传组分，南方群体（特别是CJK、FCG和MJK）的遗传组成与北方群体明显不同。

3.5 连锁不平衡分析

连锁不平衡（LD）衰减分析显示，所有群体的LD水平（r2）随SNP位点间物理距离增加而逐渐下降，但衰减速率存在显著差异。BYD、XZH、FCG和HHK群体LD衰减较慢，在较大物理距离上仍维持较高r2值；PJ和WSH群体LD衰减最快，r2值在最初几百千碱基内迅速下降；CJK、DLJHK、MJK和LQW群体呈现中等衰减模式。半衰减距离在群体间差异较大，PJ最短（167.898 kb），而BYD、XZH、FCG和HHK最长（1000 kb）。

四、讨论与结论

讨论部分指出，LQW群体位于深圳，其独特的遗传模式和较低的遗传多样性可能与该地区密集的人类活动相关，历史引种或人为介导的传播可能促成了这一格局，但需进一步采样和群体历史分析加以验证。BYD位于雄安新区，是华北最大的淡水湿地，其独特的遗传特征可能反映地理隔离、局部群体历史或环境条件的影响。XZH的分化提示局部群体过程可能促进了群体分化。PJ群体采自辽河口湿地，其相对较高的遗传多样性与连续生境、环境异质性和群体连通性相关。WSH群体来自亚洲最大的草甸型湖泊湿地微山湖，其遗传特征为湖泊型芦苇群体的进化历史提供了新见解。

功能富集分析中，植物–病原互作通路的富集提示病原相关防御机制在芦苇群体中可能发挥作用，不同谱系间病原易感性差异表明宿主遗传背景可能影响植物–病原互作。糖代谢相关通路的富集与芦苇蔗糖合成、运输和利用密切关联于碳分配和根状茎发育，支持其多年生生长策略，同时盐度等环境条件下可溶性糖、淀粉积累和糖酵解活性的变化凸显了碳水化合物代谢在渗透调节和植物生理响应中的重要性。磷代谢相关基因的富集提示养分相关过程的变异可能促成群体间遗传分化。

遗传多样性空间格局方面，PJ群体较高的遗传多样性可能归因于其大面积连续生境、环境异质性和群体连通性。所有群体观察杂合度普遍高于期望杂合度导致的负F_IS值，可能反映异交和基因流的作用，但芦苇复杂的多倍体基因组背景、二倍体参考基因组和二倍体基因型模型可能引入杂合度和F_IS估计的不确定性，等位基因剂量未估计且旁系同源区域未专门过滤，这些限制需在解释负F_IS值时加以考虑。长江流域与遗传结构的南北转变相吻合，但应被解释为与遗传结构相关的地理过渡区，而非已证实的基因流物理屏障。

ADMIXTURE分析揭示的九个遗传簇反映了芦苇作为广布种的复杂群体遗传结构，其广泛地理分布使群体暴露于多样化环境条件，促进了空间遗传结构的形成与维持。这一模式与玉米（最优K=15）和鹰嘴豆（最优K=13）等其他广布物种的研究一致，支持广布物种因长期适应异质环境和地理隔离群体而可能蕴含复杂遗传结构的假说。

研究结论指出，本研究对中国东部芦苇群体产生了100.33 Gb高质量测序数据，鉴定出13,429,377个SNP和1,487,273个InDel。功能富集分析将植物–病原互作和碳水化合物代谢确定为与群体间遗传变异和分化相关的潜在生物学过程。群体遗传分析揭示了地理群体间的显著分化，LQW、BYD和XZH群体呈现独特遗传结构，PJ和WSH共享最高的H_e、PIC和Shannon指数值，CJK的H_o和MAF值最高。沿纬度梯度观察到明显的南北遗传结构转变，长江流域与该转变区域吻合。十个群体中检测到九个遗传簇，反映了芦苇作为广布种的遗传多样性和复杂进化历史。该研究为芦苇群体的空间遗传结构和遗传分化提供了新见解，并为未来进化历史、群体连通性及湿地群体保护管理研究提供了基因组资源。