《Ecosphere》：Glacial refugia, postglacial dynamics, and hybrid zones of Pinaceae in Eurasia captured from sedimentary ancient DNA
研究人员应用杂交捕获(hybridization capture)方法靶向云杉属(Picea)、冷杉属(Abies)、落叶松属(Larix)及松属(Pinus)完整叶绿体基因组(chloroplast genome, cpDNA)，对源自欧洲、亚洲及阿拉斯加19个沉积岩芯共238份样品（涵盖过去40,000年）提取的沉积古DNA(sedimentary ancient DNA, sedaDNA)进行分析。方法学改进提升了测序深度与基因组覆盖度，使低丰度类群可被检出并建立半定量针叶树时间序列。基于单核苷酸多态性(single-nucleotide polymorphism, SNP)的物种鉴定、混合分析(admixture analysis)及单倍型网络(haplotype network)重建推断种群动态。研究发现此前未被记录的末次冰期极大(Last Glacial Maximum, LGM)松科避难种群分布于西泰梅尔(Taymyr)、南雅库特(southern Yakutia)及远东西伯利亚。冰后重定殖(recolonization)始于晚冰期(Late Glacial)先锋种落叶松属(Larix)，随后为早—中全新世逐步扩大的云杉属(Picea)与松属(Pinus)。青藏高原东部晚冰期至全新世以粗枝云杉(Picea asperata)和紫果云杉(Picea purpurea)为优势种；欧洲冰后动态显示高纬度（科拉半岛Kola Peninsula）早为欧洲赤松(Pinus sylvestris)主导，阿尔卑斯山及迪纳拉山脉(Dinaric Alps)为瑞士五针松(Pinus cembra)与矮山松(Pinus mugo)，晚全新世后被欧洲云杉(Picea abies)（捷克、俄西南、可能南瑞典）取代，而德国北部仍以北欧赤松为主。研究人员还鉴定出早全新世科拉半岛潜在杂交云杉（可能为芬云杉Picea fennica），及欧—亚云杉(Picea abies–Picea obovata)、西伯利亚冷杉(Abies sibirica)–欧洲冷杉(Abies alba)、欧洲落叶松(Larix decidua)–西伯利亚落叶松(Larix sibirica)杂交带，西伯利亚落叶松种间混合及西伯利亚红松(Pinus sibirica)–偃松(Pinus pumila)渐渗。综上，松科sedaDNA杂交捕获为欧亚古植被动态提供高分辨率视角，可支撑未来古基因组研究与森林恢复管理。
《欧亚松科(Pinaceae)冰期避难所、冰后动态及杂交带的沉积古DNA(sedaDNA)杂交捕获研究》论文解读
研究背景与立项依据
北方针叶林（Boreal forests）约占全球森林面积27%，以松科(Pinaceae)常绿或落叶针叶树种（云杉属Picea、冷杉属Abies、落叶松属Larix、松属Pinus）为主，提供气候调节、冻土稳定及生物栖息地等关键生态系统服务。现代针叶林分布格局主要受第四纪冰川旋回——特别是末次冰期极大(Last Glacial Maximum, LGM, 23–19 cal ka BP)——而非仅现代气候塑造。欧洲松科类群传统认为主要幸存于南欧避难所（比利牛斯、阿尔卑斯、喀尔巴阡、南乌拉尔），西伯利亚则存留于乌拉尔、阿尔泰—萨彦、贝加尔及东北西伯利亚避难所，北美存留于东部白令及南部避难所；但高纬度欧亚是否存在未记载的LGM避难所、物种水平(species-level)冰后重定殖(postglacial recolonization)次序、以及避难所内不同松科物种接触形成杂交带(hybrid zone)并促进冰后扩张的机制均不明确。传统孢粉(pollen)与植物大化石(plant macrofossil)受产量/传播偏差及低丰度类群漏检限制，此前sedaDNA靶向捕获仅限单属（如Larix），缺乏基于全叶绿体基因组(chloroplast genome, cpDNA)SNP的物种鉴定与群体遗传分析。为此研究人员开展跨欧亚—阿拉斯加19个沉积档案的Pinaceae四大属cpDNA杂交捕获sedaDNA研究。本文发表于《Ecosphere》。
主要关键技术方法
研究人员采集欧洲、亚洲及美国阿拉斯加15个湖泊、2个古湖、1个洞穴与1个滑塌沉积共238份样品（跨度过去~40 ka），提取sedaDNA构建单链DNA文库；设计云杉属（新组装西伯利亚云杉Picea obovata cp基因组，GenBank: PV631374）、西伯利亚落叶松Larix gmelinii、偃松Pinus pumila、西伯利亚冷杉Abies sibirica四类探针进行单轮杂交捕获(hybridization capture)富集目标cpDNA，Illumina测序后经去重复、Kraken2分类、竞争比对至四属串联参考基因组，对Pinaceae reads做属级定量；进一步比对至远缘外群cp参考基因组，通过FreeBayes调用SNP，以属/种特异性SNP进行物种判定，并对各属SNP数据做ADMIXTURE混合分析及TCS单倍型网络重建，结合古代DNA损伤(mapDamage)验证真实性。
研究结果
Early migration of pine and spruce in northwestern Russia（西北俄罗斯松树与云杉的早期迁移）
科拉半岛(Kola Peninsula) Lake Imandra记录显示晚冰期—Younger Dryas微弱Pinaceae信号；早全新世(10.5 cal ka BP起)欧洲赤松(Pinus sylvestris)显著扩增并于5.7 ka BP达峰，后渐降；云杉低丰度出现(10.5–5.1 ka BP)，SNP鉴定以西伯利亚云杉(Picea obovata)变体为主，但ADMIXTURE识别出区别于欧/西伯利亚云杉的独特祖源，推测为二者杂交种芬云杉(Picea fennica)；该祖先成分亦见于芬兰Lake Nautaj?rvi中晚全新世、泰梅尔Lake Lama晚冰期及迪纳拉Zminje jezero，TCS网络显示Imandra中—晚全新世样本更近西伯利亚/南雅库特而非俄西南种群，表明东/中亚云杉西迁与北上的欧洲云杉在俄平原杂交形成Picea fennica并向科拉低密扩散。
Eastern spruce populations expanded to southwestern Russia in the Late Glacial（晚冰期东欧云杉向东欧西南部扩展）
俄特维尔州(Tver Oblast) Lake Ezzhinitsa晚冰期(13.5 cal ka BP)具高产云杉序列，SNP呈Picea abies–Picea obovata混合且ADMIXTURE/K均值归属西伯利亚云杉祖源（同Lake Lama/Lake Khamra），TCS亲缘关系近北/东亚群而非中欧，说明晚冰期西伯利亚或东源杂交云杉沿西伯利亚→北欧→西南俄路线分布；全新世样本祖源转为中心欧Picea abies型，反映早全新世中欧种群北迁定殖。
Spruce presence near the Laurentide ice margin in southern Sweden during the Late Glacial（晚冰期瑞典南部劳伦泰德冰盖边缘的云杉存在）
瑞典Attek?p与H?sseldala Port湿地晚—末冰期(14.8–10.7 cal ka BP)检出低丰度Picea reads及典型aDNA损伤特征，虽不足SNP定种但佐证晚冰期瑞典南部近冰缘已有Picea abies局部存活或极早迁入，与先前sedaDNA metabarcoding结论相符。
Spruce and fir presence in Czechia in the Late Glacial and Early Holocene（波西米亚晚冰期与早全新世云杉与冷杉出现）
捷克北波西米亚Otto M?rtsch洞穴晚冰期(12.9 cal ka BP, 280 reads) SNP鉴定为Picea abies，早于孢粉记录(~9 ka BP)；早全新世(10–9.4 ka BP)高通量显示Picea abies、瑞士五针松(Pinus cembra)、欧洲冷杉(Abies alba)共存，后者较孢粉早约3000年；南波西米亚?erné jezero晚冰期—早全新世同样提前检出Picea与Abies低信号。
Mid-to-Late Holocene expansion of Pinaceae across central and northern Europe and potential evidence of population contact（中—晚全新世松科在欧洲中—北部的扩张及种群接触证据）
?erné jezero自6.3 ka BP起Picea abies（3489–39,466 reads）与Abies alba大幅扩增，Pinus sylvestris次之，Larix decidua少量；瑞士Sulsseewli早全新世Pinus cembra占优，~4 ka BP后Picea abies二次扩增并超松属；德北Lake Gro? Glienicker全全新世Pinus sylvestris恒为优势；芬南Lake Nautaj?rvi与俄Lake Imandra中晚全新世Picea超越Pinus；TCS显示捷克、瑞士、德、芬、俄西欧云杉样本聚为中—南欧支，区别于俄西南晚冰期东源支，说明全新世西欧云杉北向协调扩张并与东源种群接触，Imandra渐被西欧谱系取代。
Potential low-abundance Abies sibirica populations may have reached the Kola Peninsula during the Middle Holocene expansion（中全新世扩张中西伯利亚冷杉低丰度种群可能达科拉半岛）
Lake Imandra与Lake Nautaj?rvi早—中全新世检出低丰度Abies，SNP以西伯利亚冷杉(Abies sibirica)变体为主，暗示早全新世Abies sibirica西迁遇北上欧洲冷杉(Abies alba)产生杂交并向北扩至科拉。
Holocene Pinus dynamics in Eurasia（欧亚松属全新世动态）
德北恒为Pinus sylvestris主导；芬兰早—中全新世Pinus→Picea演替；捷克早全新世Pinus cembra高丰度，瑞士Sulsseewli早—中全新世Pinus cembra为主后降；西伯利亚Lake Lama晚冰期—LGM检出低丰Pinus sylvestris（为已知最北LGM松分布），Lake Khamra中全新世后出现西伯利亚红松(Pinus sibirica)/Pinus cembra subsp. sibirica及与偃松(Pinus pumila)混合信号，Lake Bolshoe Toko东南西伯利亚以Pinus pumila为主偶见Pinus sibirica侵入，反映两分类群接触带杂交。
Holocene Pinaceae dynamics in the Swiss Alps（瑞士阿尔卑斯山松科全新世动态）
Sulsseewli早全新世Pinus cembra高峰，少量Picea abies与Abies alba；~4 ka BP起Picea abies大幅扩增超松属，伴随人类活动影响，晚全新世零星Larix decidua与Abies。
Widespread low-abundance larch populations and evidence of hybridization across Europe during the Holocene（欧洲广泛低丰度落叶松种群及全新世杂交证据）
科拉、德北、捷克、迪纳拉、芬南、瑞士均检出早—中—晚全新世低丰Larix（传统孢粉/大化石未见），SNP呈欧洲落叶松(Larix decidua)–西伯利亚落叶松(Larix sibirica)高度混合，ADMIXTURE可分欧/西伯利亚祖源，欧区归Larix decidua，西伯利亚Taymyr/Lake Khamra为Larix sibirica，暗示早全新世Larix sibirica西迁与Larix decidua接触杂交，欧区残存低密度种群。
Early presence and lineage shift of Picea abies, with sparse Larix decidua in Late Glacial pine-dominated woodlands of the Dinaric Alps（迪纳拉阿尔卑斯晚冰期欧洲云杉早现、谱系转换及早低丰落叶松）
黑山Zminje jezero Younger Dryas(12.1–11.8 ka BP)松属(Pinus mugo SNP鉴定)占优伴低Picea abies，11.1 ka BP起Picea abies扩增超松属；晚冰期样(11.8 ka BP)ADMIXTURE示东源杂交祖源（类Picea fennica），早全新世样转西欧Picea abies祖源，TCS证实晚冰期东源支→早全新世西欧支转换；另首报该地晚冰期低丰Larix decidua。
Glacial persistence and Holocene genetic shift of Pinaceae in northern Siberia（北西伯利亚松科冰期持存与全新世遗传转换）
泰梅尔西Lake Lama LGM(23–19.2 ka BP)持续检出Picea(最丰)>Larix>Pinus>Abies，为四大属最北LGM避难所；LGM—晚冰期Picea呈Picea abies–Picea obovata混合且ADMIXTURE偏西欧祖源（暗示西欧云杉或杂交种曾达中西伯利亚），Younger Dryas后转Picea obovata主导；Larix LGM—晚冰期为Larix gmelinii–Larix sibirica混合，11.5 ka BP起因Larix sibirica占优，早全新世Picea obovata扩增，中全新世Larix衰退、Picea续扩，反映YD气候事件对遗传谱系的筛选效应。
Southwestern Yakutia as a glacial refugium for Pinaceae and Holocene expansion patterns（西南雅库特作为松科冰期避难所及全新世扩张模式）
Lake Khamra LGM—晚冰期以Larix gmelinii为主伴低Picea/Pinus为避难种群，Younger Dryas小幅衰退，中—晚全新世Picea obovata、Larix sibirica（晚冰期后为Larix sibirica型）、Pinus sibirica及Abies sibirica扩增，说明西南雅库特为LGM多属避难所且全新世接纳西伯利亚种混交带东扩。
Pinaceae dynamics in the southeastern Siberian Mountains during the last 34 ka（东南西伯利亚山地过去34 ka松科动态）
Lake Bolshoe Toko自MIS 3起Larix(主要为Larix gmelinii与Larix cajanderi混合/交替优势)>Picea>Pinus，LGM持续持存，Picea主要为鱼鳞云杉Picea jezoensis(MIS 3最老层为Picea obovata型)，Pinus以偃松Pinus pumila为主伴Pinus sibirica渗入，反映东南西伯利亚为长期多属持存区及Larix gmelinii–Larix cajanderi杂交带。
Larix dynamics in the hybridization zone of Larix gmelinii and Larix cajanderi in the East Siberian Mountains over the last 40 ka（东西伯利亚Larix gmelinii–Larix cajanderi杂交带40 ka落叶松动态）
Oyumyakon高地Lake Ulu MIS 3–LGM Larix持存，早为Larix gmelinii主配Larix cajanderi变体，LGM前后转Larix cajanderi主导型，说明该区既是LGM避难所也是两落叶松杂交带，幸存类群偏东种（可能经适应性渐渗增强抗寒性）。
Northeastern Tibetan Plateau(Nianbaoyuze Mountains)（青藏高原东北部——年保玉则山区）
Lake Ximen LGM极低Picea reads鉴定为粗枝云杉Picea asperata，早全新世(9.4 ka BP)Picea asperata与紫果云杉Picea purpurea共扩后衰退，伴极低Larix/Pinus/Abies；晚全新世信号弱符合气候变冷干。
Southeastern Tibetan Plateau(Hengduan Mountains)（青藏高原东南部——横断山区）
Lake Naleng晚冰期极低Pinus/Picea/偶Abies，中全新世Picea asperata为主伴Picea purpurea，晚全新世衰退，与区域植被演变一致。
Picea forest dynamics in south-central Alaska: Early Holocene colonization and Late Holocene high-density expansion（南中阿拉斯加云杉林动态：早全新世定殖与晚全新世高密度扩张）
阿拉斯加基奈半岛Hidden Lake无晚冰期—早全新世前Pinaceae，8.7 cal ka BP起白云杉Picea glauca大量扩增并保持中全新世高位，3.5 ka BP再升且单倍型多样性增，晚全新世偶见低Larix/Pinus/Abies，符合孢粉记录并细化定殖时序。
讨论与结论总结
研究人员通过改良松科四属完整cpDNA杂交捕获sedaDNA流程获得高覆盖度数据，重建了欧亚—阿拉斯加过去40 ka Pinaceae属/种水平时空动态。方法相较传统孢粉/大化石能敏感检出低丰度及此前未记载种群，虽有暖期DNA保存偏差、远距花粉cpDNA捕获及近缘种偏向性等局限，结合古代损伤验证与阴性对照可有效甄别。
新识别西泰梅尔(Taymyr)、南雅库特(southern Yakutia)及远东西伯利亚为LGM非冰川区多属避难所，其中Lake Lama同时存Picea abies–Picea obovata及Larix gmelinii–Larix sibirica杂交种群，Lake Khamra以Larix gmelinii为主伴Picea/Pinus，Lake Bolshoe Toko含Picea jezoensis/Pinus pumila/Pinus sibirica，Lake Ulu为Larix gmelinii–Larix cajanderi杂交避难带，杂交可能增强LGM适应性与冰后扩张力。
冰后动态显示：欧洲高纬/高山先由先锋Pinus sylvestris或Pinus cembra定植，升温湿后Picea abies（捷克/俄/瑞/芬）或Picea obovata（西伯利亚）替代，德北仍Pinus sylvestris主导；欧区广泛检出自Larix decidua低丰种群并具Larix decidua–Larix sibirica混合信号；西伯利亚Larix为早演替种，温暖期Picea obovata/Picea jezoensis扩伸；青藏高原东部全新世以Picea asperata+Picea purpurea为主；阿拉斯加早全新世Picea glauca定殖晚全新世密林扩伸。
冰川旋回致隔离种群接触促杂交：欧俄平原早全新世西伯利亚/中亚云杉西迁遇北上Picea abies形成Picea fennica（杂交带或产物向科拉扩散）；Lake Lama LGM为Picea abies–Picea obovata杂交带后全新世西移出该带；欧区Larix示Larix decidua–Larix sibirica历史渐渗；西伯利亚Lake Lama LGM为Larix gmelinii–Larix sibirica杂交带，Lake Khamra晚冰期为Larix gmelinii近带，早全新世Larix sibirica西扩使两区全新世均入杂交带；Lake Ulu存Larix gmelinii–Larix cajanderi古老杂交与LGM偏向Larix cajanderi型保存；Abies sibirica早全新世西迁遇Abies alba北扩可能产杂交。
综上，Pinaceae sedaDNA杂交捕获为古植被动态、冰期避难所及杂交热点提供高分辨率证据，可支撑古基因组学与森林恢复管理研究。
