为抢救性保护该物种，中国科学院昆明植物研究所极小种群野生植物综合保护团队马永鹏研究组联合中国林业科学院资源昆虫研究所马宏课题组、中国环境科学研究院刘勇波研究组和云南省农业科学院花卉研究所王继华研究团队，对朱红大杜鹃以及近缘广布物种马缨杜鹃进行了一系列保护基因组学比较分析。该研究首次获得朱红大杜鹃染色体水平的高质量全基因组，基因组大小677 Mb，contig N50为34 Mb，scaffold N50 为53 Mb，组装质量大大超过了当前杜鹃花属所有完成测序的物种，并注释了38,280个蛋白编码基因。通过与其他物种进行系统发育研究表明，杜鹃花目与山茱萸目在大约于100 Ma前发生分化。

图2 朱红大杜鹃遗传多样性和有效种群大小变化。（a）17个木本植物的遗传多样性比较；（b）基于Stairway plot的朱红大杜鹃有效种群大小变化；（c）基于fastsimcoal2的朱红大杜鹃有效种群大小变化，T代表时间，深蓝色图形的大小及上面的数字代表有效种群大小。

为进一步解析朱红大杜鹃的濒危原因，我们对朱红大杜鹃进行了重测序，群体历史分析结果表明，历史上经历3次严重遗传瓶颈，与冰期-间冰期作用、“共和运动”等地质历史事件一致。最后一次遗传瓶颈之后，朱红大杜鹃的有效种群大小逐渐恢复。但种群调查大小仅数百株，这可能与人类活动如水库建设、采挖和修路等有关。同时，我们对同属广布种马缨杜鹃进行了重测序，并与朱红大杜鹃进行了比较分析，结果表明，朱红大杜鹃显著积累了更多数量的纯合有害突变；朱红大杜鹃的长片段纯合（Runs of homozygosity，ROH）程度显著高于马缨杜鹃，表明朱红大杜鹃近交严重。此外，与马缨杜鹃等其他16个具有染色体水平全基因组信息的木本植物相比，朱红大杜鹃的遗传多样性非常低。

图3. 遗传负荷和近交。桔红色柱子代表朱红大杜鹃，蓝绿色柱子代表马缨杜鹃。（a）有害突变数量比较。（b）杂合和纯合有害突变数量比较。（c）长纯合片段比较。

总之，人类活动导致的生境丧失、遗传多样性极低、地质历史事件导致的遗传瓶颈、近交和与热适应相关基因的有害突变，是朱红大杜鹃“极小种群”形成和维持的主要原因。最后，本文对朱红大杜鹃抢救性保护行动提出了建议。相关结果以“Chromosome-level genome assembly and population genetic analysis of a critically endangered rhododendron provide insights into its conservation”为题发表在植物学三大国际主流期刊、中科院1区期刊The Plant Journal上（链接地址：https://doi.org/10.1111/tpj.15399）。期刊编辑给予文章较高评价：“This is an outstanding manuscript and the more detailed and appropriate experimental design makes this a very valuable and scientifically sound study”。

该研究得到生态环境部生物多样性调查与评估、第二次青藏高原综合科学考察研究——青藏高原旗舰、极度濒危和球根宿根类植物的调查评估与综合保护（杜鹃专题）、国家自然科学基金、云南省中青年学术与技术带头人等项目的资助。

作者：马永鹏、刘德团（云南省极小种群野生植物综合保护重点实验室）