1. 领域背景与文献
文献英文标题：Ancient genomics study reveals low HLA diversity in eastern hunter-gatherers；
发表期刊：BMC Genomics；影响因子：未公开；研究领域：古基因组学、人类群体遗传学、免疫遗传学
古人类群体遗传学领域中，东欧狩猎采集者（EHG）是公元前14000至5000年活跃于欧洲的重要群体，其基因组主要由来自西伯利亚的古北亚人（ANE）祖先成分，以及少量源自Villabruna集群的旧石器时代成分构成，与西欧狩猎采集者（WHG）共享约14000年前的共同祖先。EHG与WHG、早期农民（EF）存在地理接触但未发生大规模基因混合，公元前5000年后，草原牧民的扩张导致EHG作为独立文化群体消失，但其基因成分通过草原牧民与当地农民的融合，持续保留在现代欧洲人群中（如Y染色体单倍群R1b1a1）。当前领域研究热点聚焦于古人类免疫相关基因的演化机制，尤其是人类白细胞抗原（HLA）基因，该基因家族是适应性免疫的核心介质，其多样性与病原体压力密切相关，然而此前针对EHG的人类白细胞抗原研究仅能基于标签SNP进行间接推断，存在标记不全、种群特异性不足的局限性，无法准确揭示EHG的人类白细胞抗原多样性特征及演化驱动因素。未解决的核心问题包括EHG的人类白细胞抗原多样性水平如何、与现代人群的差异由何种因素导致、病原体暴露是否参与塑造其人类白细胞抗原等位基因库等。在此背景下，本研究以俄罗斯西部Sakhtysh遗址的EHG样本为研究对象，通过古基因组测序、病原体筛查及人类白细胞抗原靶向富集分析，首次获得EHG的直接人类白细胞抗原基因型数据，旨在明确EHG的遗传连续性、病原体暴露谱及人类白细胞抗原多样性演化特征，为古人类免疫基因演化研究填补关键空白。
2. 文献综述解析
作者对领域现有研究的分类维度涵盖EHG群体遗传学研究、古人类人类白细胞抗原基因型研究、病原体压力与人类白细胞抗原演化关联研究三个方向，系统梳理了各方向的研究进展与局限性，进而明确本研究的创新定位。
EHG群体遗传学研究方面，已明确EHG的地理分布范围（从波罗的海、黑海至乌拉尔山脉）、生存策略（渔猎采集）、基因组构成及与其他古人类群体的互动模式，证实EHG对现代欧洲人群的基因贡献，但针对单一遗址长期时间尺度的遗传连续性研究仍较为有限；古人类人类白细胞抗原基因型研究方面，已报道新石器时代农民的人类白细胞抗原基因型数据，发现其与现代欧洲人群存在显著差异，但针对狩猎采集者的研究仅能依赖标签SNP进行等位基因频率推断，该方法受限于标记覆盖度不足及种群特异性差异，无法获得准确的人类白细胞抗原基因型；病原体压力与人类白细胞抗原演化关联研究方面，已建立人类白细胞抗原多样性与病原体多样性正相关的理论框架，但针对古狩猎采集者的病原体谱数据匮乏，无法直接验证病原体压力对人类白细胞抗原演化的驱动作用。
针对现有研究中缺乏EHG直接人类白细胞抗原基因型数据的核心空白，本研究采用靶向富集测序技术，首次获得27名EHG个体的人类白细胞抗原基因型（覆盖I类和II类共6个位点），解决了标签SNP推断的局限性；同时结合群体遗传连续性分析、病原体筛查数据，系统揭示EHG低人类白细胞抗原多样性的特征及可能驱动因素（小种群规模、低病原体多样性），为古人类免疫基因演化的研究提供了直接的实验证据，弥补了狩猎采集者群体在古人类免疫遗传学研究中的缺失环节。
3. 研究思路总结与详细解析
本研究的整体框架以Sakhtysh遗址EHG群体为研究对象，核心目标是明确该群体的遗传连续性、病原体暴露谱及人类白细胞抗原多样性演化特征；核心科学问题包括EHG群体在2000年文化变迁中是否保持遗传连续性、其人类白细胞抗原多样性低于现代人群的驱动因素是什么、病原体暴露与人类白细胞抗原等位基因频率是否存在关联；技术路线遵循“样本采集与测序→群体遗传分析→病原体筛查→人类白细胞抗原分型与多样性分析→结论推导”的闭环逻辑，通过多组学数据整合实现研究目标。
3.1 样本采集与基因组测序
实验目的：获取Sakhtysh遗址古人类样本的高质量基因组数据，为后续群体遗传、病原体及人类白细胞抗原分析提供基础。
方法细节：从Sakhtysh遗址的43名个体中采集84件骨骼样本（包括岩骨和牙齿），严格遵循古DNA实验的洁净室操作规范进行DNA提取，采用部分尿嘧啶-DNA糖基化酶处理的方法制备文库，对所有样本进行全基因组鸟枪测序（Illumina HiSeq 6000平台，2×100 bp读长）；同时对31名保留岩骨的个体进行人类白细胞抗原区域的靶向富集测序，以获得高覆盖度的人类白细胞抗原序列数据。
结果解读：43个样本中，32个样本的1240K SNP panel分型位点数量≥20000，符合群体遗传分析的质量标准；结合已发表的4个Sakhtysh样本数据，最终纳入36个样本进行后续分析；通过检测DNA末端C→T替换比例、线粒体DNA污染水平、男性样本X染色体杂合性等指标，确认所有样本的古DNA真实性，并通过性染色体与常染色体的测序深度比例确定个体的分子性别。
产品关联：文献未提及具体实验产品，领域常规使用古DNA提取试剂盒、Illumina测序文库制备试剂盒、人类白细胞抗原靶向富集探针组等。
3.2 群体遗传连续性与亲缘关系分析
实验目的：验证Sakhtysh群体在Lyalovo至Volosovo文化阶段（跨度约2000年）的遗传连续性，分析个体间的亲缘关系模式。
方法细节：对36个样本进行主成分分析（PCA），将古人类样本投影到现代西欧亚人群的PCA空间中；通过f3统计分析Sakhtysh个体与其他古人类群体的遗传漂移共享水平，采用f4统计验证不同阶段个体的遗传相似性；使用无监督Admixture分析解析样本的祖先成分构成；采用READv2软件分析个体间的亲缘关系，限定检测二级以内的亲属关系。
结果解读：PCA结果显示，所有Sakhtysh个体均聚类在EHG群体的分布范围内，与俄罗斯其他遗址的EHG样本（如Sidelkino、Popovo）遗传距离较近；f3/f4统计结果表明，不同文化阶段的Sakhtysh个体与参考个体的遗传漂移水平无显著差异，Admixture分析显示所有个体的祖先成分（ANE和WHG）保持稳定，证明该群体在2000年的文化变迁中维持了遗传连续性；亲缘关系分析未发现一级亲属对，仅检测到3对二级亲属对，且均属于Late Volosovo阶段，提示该遗址的埋葬群体并非紧密的亲缘群体。
3.3 病原体筛查与分析
实验目的：检测Sakhtysh个体的病原体暴露情况，为人类白细胞抗原演化的病原体压力驱动假说提供直接证据。
方法细节：使用MEGAN Alignment Tool（MALT）对全基因组测序数据进行宏基因组病原体筛查，设置90%的序列相似性阈值，比对细菌和病毒参考基因组数据库；对阳性样本将测序数据映射到对应病原体的参考基因组，验证古DNA的特征（如末端脱氨模式）以确认感染的真实性。
结果解读：在2个样本中检测到猪红斑丹毒丝菌（Erysipelothrix rhusiopathiae）的DNA，序列显示典型的古DNA脱氨特征，推测个体曾发生败血症感染，该病原体为人畜共患病原体，可能通过渔猎活动中的皮肤伤口感染；4个样本检测到人细小病毒B19，5个样本检测到乙型肝炎病毒（HBV），其中1个样本同时感染两种病毒，HBV基因组与“西欧亚新石器至青铜时代”（WENBA）谱系最接近，支持该病毒谱系的传播早于农民扩张的假说；此外还检测到9种与牙菌斑、牙周炎相关的口腔细菌。
产品关联：文献未提及具体实验产品，领域常规使用MEGAN、BWA等宏基因组分析软件。
3.4 人类白细胞抗原靶向富集、分型与多样性分析
实验目的：获得Sakhtysh个体的准确人类白细胞抗原基因型，分析其多样性特征、与其他群体的差异及潜在的选择压力。
方法细节：对31个样本的人类白细胞抗原靶向富集测序数据，采用TARGT pipeline进行人类白细胞抗原分型（覆盖I类HLA-A、HLA-B、HLA-C及II类HLA-DPB1、HLA-DQB1、HLA-DRB1位点），同时使用OptiType软件对I类人类白细胞抗原分型结果进行验证，仅保留两种方法一致的分型结果；计算Shannon多样性指数评估人类白细胞抗原多样性，通过主坐标分析（基于FST值）比较Sakhtysh与其他群体的遗传分化，采用Ewens-Watterson检验评估人类白细胞抗原位点的选择压力。
结果解读：27个样本获得有效人类白细胞抗原基因型，其中13个样本完成6个位点的全部分型；Shannon多样性指数分析显示，Sakhtysh群体的人类白细胞抗原多样性显著低于现代欧洲人群，甚至低于新石器时代农民群体；HLA-B*27:05等位基因在Sakhtysh群体中的频率约为50%（n=27，文献未明确提供该数据的统计学显著性），该等位基因在现代人群中相对罕见；主坐标分析显示，Sakhtysh与新石器时代农民的遗传分化程度（FST值）小于与现代人群的分化；Ewens-Watterson检验显示，HLA-B位点存在定向选择的趋势，HLA-C、HLA-DRB1、HLA-DQB1位点存在平衡选择的趋势，但所有检验结果均未达到统计学显著性（文献未明确提供P值）。
产品关联：文献未提及具体实验产品，领域常规使用人类白细胞抗原分型软件（TARGT、OptiType）及群体遗传分析软件（R语言hierfstat包）。
4. Biomarker研究及发现成果解析
本研究中涉及的Biomarker为人类白细胞抗原等位基因（包括I类和II类共6个位点的等位基因），作为反映古人类免疫适应特征的分子标志物，其筛选与验证逻辑为“靶向富集测序获得高覆盖度人类白细胞抗原序列→双软件联合分型确认基因型→与不同时期群体比较分析多样性与频率特征→检验选择压力明确演化驱动因素”，形成完整的研究链条。
明确Biomarker类型为人类白细胞抗原等位基因，包括I类的HLA-A、HLA-B、HLA-C及II类的HLA-DPB1、HLA-DQB1、HLA-DRB1位点的等位基因；筛选与验证逻辑为，首先对Sakhtysh遗址的古人类样本进行人类白细胞抗原区域靶向富集测序，获得高覆盖度的人类白细胞抗原序列，然后采用TARGT pipeline和OptiType两种软件进行分型，仅保留两种方法一致的结果以确保基因型准确性，随后将Sakhtysh群体的人类白细胞抗原基因型与新石器时代农民、现代欧洲人群进行比较，分析多样性特征与等位基因频率差异，最后通过Ewens-Watterson检验评估选择压力，明确人类白细胞抗原等位基因的演化驱动因素。
Biomarker的来源为Sakhtysh遗址古人类的岩骨样本，通过靶向富集测序获得人类白细胞抗原序列数据；验证方法采用两种独立的人类白细胞抗原分型软件（TARGT和OptiType）联合确认，确保分型结果的准确性；特异性方面，分型结果针对人类白细胞抗原的6个核心位点，覆盖I类和II类免疫基因，能够准确反映群体的免疫遗传特征；敏感性方面，31个进行靶向富集的样本中27个获得有效基因型，分型成功率约为87%（n=31）。
本研究的核心成果包括三个方面，一是发现Sakhtysh EHG群体的人类白细胞抗原多样性显著低于现代欧洲人群，甚至低于新石器时代农民群体，推测该特征可能由小种群规模和低病原体多样性共同驱动（文献未明确提供统计学P值）；二是发现HLA-B27:05等位基因在Sakhtysh群体中的频率约为50%（n=27），该等位基因在现代人群中相对罕见，推测其高频率可能与古代病毒感染压力相关，而新石器时代后选择压力的转变（从病毒到细菌感染）导致其频率下降；三是发现DRB108:01-DQB1*04:02单倍型在Sakhtysh群体中的频率约为40%（n=27），该单倍型可能源自EHG与WHG的共同祖先，为古人类群体的共同起源提供了免疫遗传学证据。此外，本研究首次获得EHG的直接人类白细胞抗原基因型数据，填补了狩猎采集者群体在古人类免疫遗传学研究中的空白，为理解人类免疫基因的演化历程提供了关键的古数据支撑。