引言
桫椤科隶属于真蕨目，是现存最古老的树形蕨类谱系之一，广泛分布于热带、亚热带及南半球部分温带地区。该科植物不仅是森林生态系统稳定性的重要生物指示器，其独特的树形习性、相对完整的化石记录和进化连续性，也使其成为研究植物形态进化、环境适应机制及应对气候变化响应的理想模型。此外，其产生的次生代谢物具有抗氧化、抗炎、抗肿瘤等生物活性，显示出潜在的药用价值。遗传多样性是物种应对环境变化、维持长期进化适应潜力的基础，也是制定濒危物种保护策略的核心依据。近年来，现代分子生物学的进步极大地推动了基于组学（如基因组学、转录组学等）的桫椤科研究。本文旨在系统评述多组学技术在桫椤科遗传多样性研究中的应用与成就，重点解析测序技术和功能基因如何阐明其对生物和非生物胁迫（包括气候变化、生境破碎和种间竞争）的遗传多样性、隐性谱系分化及适应性进化机制，总结当前保护应用，并展望未来研究方向。
研究技术与方法学的演变
桫椤科遗传多样性研究的技术方法经历了四个阶段的演变。
早期研究阶段（2000年前）：主要依赖形态学标记、细胞学标记（如染色体核型分析）和生化标记（如同工酶/等位酶）。这些方法操作简便，但分辨率低，易受环境干扰，难以准确量化遗传变异。
传统DNA分子标记阶段（2000–2015年）：随机扩增多态性DNA（RAPD）、扩增片段长度多态性（AFLP）、简单序列重复区间（ISSR）、简单序列重复（SSR）和叶绿体DNA（cpDNA）序列标记等技术的应用，将研究推向分子水平。其中，SSR标记因其共显性、高多态性和高重复性，成为解析种群精细遗传结构的“金标准”。cpDNA序列（如atpB-rbcL, trnL-trnF）则因母系遗传和结构保守，是研究系统发育和谱系地理学的理想工具。
基因组简化测序阶段（2016–2020年）：二代测序（NGS）技术促进了如限制性位点相关DNA测序（RAD-seq）和基因分型测序（GBS）等基因组简化测序方法的应用。这些技术单次运行可获得成千上万个单核苷酸多态性（SNP）位点，提高了研究的分辨率和全基因组覆盖度，推动研究进入全基因组水平。
多组学与整合分析阶段（2021年至今）：全基因组测序技术的突破为解析核心遗传机制提供了根本支撑。例如，已完成桫椤（Alsophila spinulosa）染色体级别的基因组组装（6.23 Gb）。全转录组测序成为挖掘功能基因的核心手段。此外，生态位模型与遗传多样性的整合，在宏观分布与微观遗传学之间架起了桥梁。蛋白质组学、代谢组学和微生物组学等组学技术进一步拓展了研究边界。
分子系统发育与分类学修订
分子数据为澄清桫椤科内部长期存在的分类争议提供了关键支持。系统发育分析支持将桫椤科分为四个稳定的分支：白桫椤属（Sphaeropteris）、广义桫椤属（Alsophila s.l.）、狭义桫椤属（Alsophila s.s.）和黑桫椤属（Gymnosphaera）。其中，白桫椤属是其他三个分支的基部分支。狭义桫椤属和黑桫椤属形成一个高支持度的单系分支。研究支持将黑桫椤属恢复为独立的属。分子数据还有助于界定有争议的物种，例如基于cpDNA trnL-F序列分析，支持将 Alsophila austroyunnanensis 更名为 Gymnosphaera austroyunnanensis。对形态性状的系统发育意义进行了重新评估，发现鳞片形态（等边的vs.具边缘的）在主要类群界定中具有关键系统发育意义，而囊群盖形态的信号较弱。
遗传多样性的核心发现
质体基因组多样性：桫椤科叶绿体基因组具有典型的四分体结构，整体保守，但也存在物种特异性结构变异。例如，笔筒树（S. brunoniana）缺少trnV-UAC基因。筛选高核苷酸多样性的基因（如atpI, ycf2）为开发DNA条形码提供了支持。
核基因组中的隐性多样性与杂交：核基因组信息揭示了广泛存在的隐性多样性。例如，基于RAD-seq的研究发现，三个广布的黑桫椤属物种各自由两个高度分化的谱系组成，可能对应隐存种。杂交和基因渗入在桫椤科进化中扮演复杂角色，检测到种间及属间广泛的杂交和基因渗入信号，渗入基因富集在“胁迫响应”、“木质部发育”等适应性通路中。多倍体化（如异源四倍体 Gymnosphaera metteniana）是新物种形成的重要途径。
种群遗传结构与历史动态：地理隔离是驱动桫椤科种群遗传分化的关键因素。例如，琼州海峡导致桫椤（A. spinulosa）海南种群与大陆种群间存在高遗传分化。第四纪冰期-间冰期旋回对其分布格局和遗传多样性产生了深远影响，云南、四川峨眉山、海南吊罗山等地可能曾是其冰期避难所。局域尺度上，生活史性状（如有限的配子体扩散、广泛的克隆生长）显著影响种群遗传结构。
适应性遗传多样性：桫椤科的树形形态与独特的功能基因家族密切相关。全基因组水平鉴定出桫椤有186个bHLH基因，数量显著高于草本蕨类，这些基因在木质部特异性高表达，在茎发育和木质化中起核心调控作用。胁迫抗性和光适应基因的多样性为其应对复杂环境提供了保障。例如，ASR基因家族在干旱生境种群中的等位基因频率显著高于湿润生境种群。次生代谢相关基因（如苯丙烷代谢通路中的PAL, C4H, 4CL）的变异，赋予其重要的生态和药用价值，如独特的酚类化合物桫椤素（alsophilin）具有抗真菌活性。
生殖系统、微生物组与遗传多样性：生殖系统通过调控遗传变异的产生和传递，深刻影响遗传多样性。克隆繁殖广泛存在，会导致种群内遗传多样性显著降低。异交物种具有较高的种群内遗传多样性和较低的种群间遗传分化。微生物组与宿主遗传多样性密切相关，桫椤根、茎、叶的内生微生物群落具有显著的组织特异性，根系相关的伯克氏菌（Burkholderia）与耐旱性相关。
遗传多样性的地理格局与驱动机制
全球尺度上，物种多样性和遗传多样性的热点集中在美洲热带（亚马逊盆地、中美洲云雾林）和亚太热带-亚热带地区（东南亚群岛、中国南部、澳大利亚东北部）。非洲大陆因气候干旱化和生境破碎化，多样性显著匮乏。区域尺度上，中国南部及东南亚可分为华南、西南、岭南和稀有特有物种分布区四种典型区域，其遗传多样性水平和种群分化程度各异。
遗传多样性格局主要由温度、降水和生境连通性塑造。极端低温和高热波动抑制孢子萌发和幼苗存活，而热带-亚热带地区稳定温暖的条件是遗传多样性热点形成的基础。最暖季降水量通过影响水分可获得性来调控桫椤科生长和孢子传播效率。城市化扩张和农业开垦导致生境“碎片化”，阻碍孢子传播，形成“遗传岛屿”。遗传变异积累直接受交配系统和扩散能力调控，异交物种种群内遗传多样性高，种群间分化低；而自交或近交种群则限制遗传变异丰富度。孢子传播效率受风速、地形和适宜斑块连通性制约。第四纪冰期-间冰期旋回产生了决定性影响，冰期迫使桫椤科退缩至“避难所”，冰期后扩张的前缘种群受奠基者效应和遗传漂变影响，遗传多样性显著低于核心避难所种群。生物因素也参与调控，种间竞争、杂交以及内生菌群落间接影响种群生长和繁殖效率。此外，桫椤科巨大的基因组、低基因突变率以及多倍体化伴随的基因沉默等现象，也从遗传特性本身制约了多样性积累。
遗传多样性研究在桫椤科保护中的应用
指导迁地保护与精准种质保存：遗传多样性研究为迁地保护提供了核心科学依据。种质收集应优先考虑保留了更多适应性基因和独特遗传谱系的避难所种群。采样需满足最小有效种群大小，并通过SSR、AFLP等分子标记覆盖种群内主要遗传谱系。对于具有隐性物种分化的类群，每个隐性谱系应视为独立的保护单元。
辅助准确物种界定与科学划定保护单元：cpDNA片段、核SNP和SSR标记可准确区分近缘种、隐存种和自然杂交个体，为界定保护目标提供支持。结合遗传分化程度、适应性基因差异和基因流强度，可将种群划分为进化显著单元（ESU）和管理单元（MU）。
解析致危机制与科学评估灭绝风险：通过量化等位基因丰富度、杂合度等参数，可直接测量种群的适应进化潜力，较低的遗传多样性与较高的灭绝风险强相关。等位基因频率分布和杂合度偏差检验可检测种群是否经历过遗传瓶颈或长期近交。种间杂交和基因渗入对一些濒危物种构成威胁，需分子标记鉴定杂交个体并制定隔离保护措施。结合种群生态学研究，可实现“遗传因子-生态因子”协同的致危机制解析。整合遗传多样性参数与生境适宜性模型，可构建“遗传-生态”二维风险评估体系。
优化自然保护区网络与监测保护成效：将遗传多样性热点与物种分布模型预测的高适宜区叠加，识别“双高”区域，为保护区选址提供科学依据。对于未纳入保护地的边缘种群，需建立保护点或生态廊道。应建立基于遗传参数的长期监测体系，定期检测等位基因丰富度、基因流强度等指标以评估保护措施成效。结合转录组学，可监测胁迫响应和适应性基因的表达变化，为环境影响提供早期预警。
整合表型-基因型关联实现生境匹配保护：通过关联表型变异和遗传分化，可实现精准的生境调控。就地保护需明确关键生境因子的阈值（如林冠透光率30%–50%，土壤pH 5.5–6.5）。迁地保护选址应建立“生境相似性指数”，优先选择指数≥80%的区域，并结合光适应和胁迫抗性相关基因表达特征提高成活率。
支持可持续利用与优化人工繁殖技术：在严格保护野生种群的前提下，可基于遗传背景筛选功能性化合物含量高的种群，适度开发药用保健资源，建立“人工栽培-成分提取-收益反哺保护”模式。在人工繁殖中，应优先选择不同地理种群或谱系的个体进行杂交，以保障后代具有足够的变异。需优化孢子消毒、培养基配方等条件，模拟原生地的温湿度和光周期，以提高人工繁殖个体的遗传适应性和定植后的成活率。
研究展望
未来研究应聚焦四个方向：短期（≤3年）扩大组学研究覆盖度，填补白桫椤属等类群的遗传空白，开展跨大陆比较研究，构建全球尺度遗传多样性格局。中期（3–5年）加强多组学整合与功能验证，构建“基因组-转录组-表观组-代谢组-微生物组”整合分析框架，利用CRISPR–Cas9基因编辑等技术验证基因功能，阐明遗传变异影响表型适应性的分子通路。长期（>5年）深化保护遗传学与实践的结合，将遗传多样性热点、隐存谱系等信息纳入保护区网络优化规划，制定精准的迁地保护采样策略，建立“自然生境+迁地”双库保护标准化流程。持续创新技术方法与研究维度，引入单细胞测序、表观组学等，并加强与近缘蕨类科的比较组学研究。
