现代人类中声音与面部解剖相关基因的差异性DNA甲基化

Differential DNA methylation of vocal and facial anatomy genes in modern humans.

作者信息David Gokhman, Malka Nissim-Rafinia, Lily Agranat-Tamir, Genevieve Housman, Raquel García-Pérez, Esther Lizano, Olivia Cheronet, Swapan Mallick, Maria A Nieves-Colón, Heng Li, Songül Alpaslan-Roodenberg, Mario Novak, Hongcang Gu, Jason M Osinski, Manuel Ferrando-Bernal, Pere Gelabert, Iddi Lipende, Deus Mjungu, Ivanela Kondova, Ronald Bontrop, Ottmar Kullmer, Gerhard Weber, Tal Shahar, Mona Dvir-Ginzberg, Marina Faerman, Ellen E Quillen, Alexander Meissner, Yonatan Lahav, Leonid Kandel, Meir Liebergall, María E Prada, Julio M Vidal, Richard M Gronostajski, Anne C Stone, Benjamin Yakir, Carles Lalueza-Fox, Ron Pinhasi, David Reich, Tomas Marques-Bonet, Eran Meshorer, Liran Carmel
PMID32132541
发布时间2020-03-04
DOI10.1038/s41467-020-15020-6

实验完整度

高

基于69个骨骼甲基化图谱和6个黑猩猩图谱,包含DMR检测、基因功能富集、表达关联及原代软骨细胞验证等多层证据。

主要模型

现代人类骨骼样本(股骨、颅骨等) 尼安德特人和丹尼索瓦人骨骼样本 黑猩猩骨骼样本 人原代软骨细胞

重点核对

DMR检测阈值:甲基化差异>50%,至少覆盖50个CpG 样本过滤:使用59个额外人类图谱过滤群体内变异,保留组间完全无重叠的位点 谱系分配:利用黑猩猩甲基化数据通过简约法确定DMR起源分支 表达关联:在22种人类组织中分析DMR甲基化与基因表达的相关性 原代软骨细胞验证:直接测量部分DMR的甲基化水平

摘要

调控元件的变化被认为是表型分化的关键驱动因素。然而,识别导致人类特异性性状的调控元件变化一直极具挑战性。在本研究中,我们利用63个重建和实验测量的古代及现代人类DNA甲基化图谱,以及6只黑猩猩的图谱,来检测在现代人类与尼安德特人和丹尼索瓦人分化后可能出现的差异性甲基化区域。我们发现,与面部和声道解剖结构相关的基因经历了特别广泛的甲基化变化。具体而言,我们在一个面部和声音相关基因网络(SOX9、ACAN、COL2A1、NFIX和XYLT1)中识别出广泛的超甲基化。我们提出,这些抑制模式出现在与尼安德特人和丹尼索瓦人分化之后,并且它们可能在塑造现代人类面部和声道方面发挥了关键作用。

实验结论

提炼研究问题、关键发现与证据,快速把握文章的核心贡献。

研究问题
现代人类与古人类(尼安德特人和丹尼索瓦人)之间DNA甲基化差异如何贡献于人类特异性性状的演化?
核心机制
在现代人类谱系中,面部和声道相关基因(SOX9、ACAN、COL2A1、NFIX和XYLT1)启动子和增强子区域发生超甲基化,导致这些基因表达下调,从而可能影响面部和喉部骨骼形态。
主要证据
通过对69个骨骼甲基化图谱(包括现代人类、古人类和黑猩猩)的比较,鉴定出873个现代人类特有的DMR,这些DMR富集于面部和喉部相关基因;SOX9、ACAN、COL2A1、NFIX和XYLT1在現代人類中呈现超甲基化,且甲基化水平与表达负相关(如NFIX);原代软骨细胞中验证了部分DMR的超甲基化。
研究意义
研究结果表明,DNA甲基化变化可能在现代人类面部和声道解剖结构的演化中发挥了作用,为从表观遗传角度研究人类特异性性状的起源提供了新途径。

研究路径

按研究推进顺序梳理实验设计、验证步骤与关键观察。

1

生成骨骼DNA甲基化图谱

获取现代人类、古人类和黑猩猩的骨骼甲基化数据,以比较不同人种群体的甲基化差异。

重建了8个古人类和现代人类的甲基化图谱,产生了2个现代人类股骨WGBS图谱,收集了52个公开的现代人类骨骼甲基化芯片数据,并产生了6个黑猩猩的甲基化图谱(1个WGBS、1个RRBS和4个850K芯片)。

2

识别差异性甲基化区域(DMR)

检测现代人类与古人类之间甲基化差异显著的区域,并确定其起源分支。

开发了针对古甲基化图谱的DMR检测算法,采用严格阈值(甲基化差异>50%,至少50个CpG),通过三步比较和黑猩猩外群进行谱系分配,并使用额外样本过滤群体内变异。

3

分析DMR相关基因的功能富集

探究现代人类特有DMR相关基因的生物学功能,特别是与解剖结构的关系。

将DMR映射到基因,使用Gene ORGANizer工具分析基因-器官关联,并进行GO富集分析、GWAS catalog重叠分析以及表达数据关联分析。

4

验证面部和喉部基因的富集

排除技术偏差和混杂因素,确认面部和喉部基因富集的可靠性。

进行了多种正交验证,包括:分析基因表达模式、重叠已知增强子、比较骨骼相关基因背景、置换检验、以及排除基因长度、细胞组成等因素的影响。

5

定量分析甲基化变化程度

评估面部和声音相关基因中甲基化变化的程度,并识别甲基化变化热点。

以100kb窗口扫描基因组,计算现代人类特有CpG的比例,比较不同基因类的变化程度,并识别变化最显著的基因(ACAN、SOX9、COL2A1、XYLT1、NFIX)。

6

研究甲基化与表达的关系

探究DMR甲基化水平是否与基因表达相关,特别是NFIX。

在22种人类组织中分析DMR甲基化与表达的相关性,利用DNMT3A诱导的MCF-7细胞甲基化数据验证NFIX表达受甲基化调控,并通过bisulfite-PCR在人类和黑猩猩颅骨样本中验证NFIX超甲基化。

7

原代软骨细胞验证

在相关细胞类型中验证DMR的甲基化模式。

从骨关节炎患者中分离原代软骨细胞,直接测量部分DMR的甲基化水平。

研究方法

按研究目的归类文中使用的方法,便于定位所需技术。

产品清单

实验环节名称品牌货号
QIAamp DNA Investigator试剂盒Qiagen56504
EpiTect亚硫酸氢盐试剂盒QIAGEN--
EZ DNA甲基化闪电试剂盒Zymo ResearchD5030
FastStart高保真PCR系统Sigma03553400001
凝胶/PCR DNA迷你试剂盒RBCYDF100
CloneJET PCR克隆试剂盒Thermo ScientificK1231
DH5α感受态细胞Invitrogen18265017
RBC质粒小提试剂盒RBCYPD100
Applied Biosystems 3730xl基因分析仪Applied Biosystems--
II型胶原酶WorthingtonLS004177
GenElute哺乳动物基因组DNA小提试剂盒SigmaG1N350
Illumina Infinium MethylationEPIC芯片Illumina--
Illumina HiSeq 2000测序仪Illumina--
Illumina NextSeq500测序仪Illumina--
Illumina HiSeqX10测序仪Illumina--

关键环节

汇总复现实验时建议重点确认的条件及原文阅读提示。

环节核对要点
骨骼样本甲基化图谱生成
样本来源、年龄、性别、骨骼类型、测序深度、甲基化数据产生技术
阅读提示:Methods中‘Present-day human bone DNA methylation maps’、‘Reconstructing ancient DNA methylation maps’及Supplementary Data 1
DMR检测
甲基化差异阈值(>50%)、最小CpG数(50)、平滑窗口(25 CpG)、FDR阈值(<0.05)、谱系分配方法
阅读提示:Methods中‘DMR detection’、‘DMR-detection algorithm’及Supplementary Methods
基因功能富集分析
基因定义(DMR重叠基因体或启动子上游5kb)、Gene ORGANizer参数、富集显著性阈值(FDR<0.05)
阅读提示:Methods中‘Gene ORGANizer analysis’及Supplementary Data 4
甲基化与表达关联分析
表达数据来源(22种组织)、相关性计算方法(Pearson)、FDR阈值(<0.05)
阅读提示:Methods中‘Computing correlation between methylation and expression’及Supplementary Data 2
原代软骨细胞验证
细胞来源(骨关节炎患者)、样本量(n=8)、甲基化检测方法(文中未明确说明)
阅读提示:Methods中‘Human primary chondrocyte validation’