💡这项突破是什么?
近日，北京大学生物医学前沿创新中心汤富酬研究员、文路副研究员与清华大学基础医学院纪家葵研究员合作，在Science杂志在线发表了题为“Single-cell long-read sequencing of chromatin accessibility unveils allele-specific chromatin states in preimplantation mouse embryos”的研究论文。该研究利用单细胞起始的scNanoATAC-seq2技术，构建了涵盖小鼠着床前发育全过程的高精度单细胞染色质可及性图谱，系统鉴定了不同发育阶段、不同谱系细胞的关键转录因子，揭示了调控合子基因组激活和谱系分化（上胚层、原始内胚层和滋养外胚层）的顺式调控网络。此外，还解析了雌性胚胎中父本X染色体的印记失活与重新激活，以及非经典印记基因的表观调控基础。这些发现为理解包括人类在内的哺乳动物早期胚胎发育的分子机制提供了新的线索。
1. 揭示了父本X染色体印记失活和重新激活的表观调控机制
哺乳动物雌性胚胎中两条X染色体中的一条会失活，以平衡雌雄之间的X连锁基因剂量。小鼠胚胎中父本X染色体在早期胚胎发育过程中经历印记失活，随后在多能性的上胚层细胞中被重新激活。该研究利用单细胞起始的scNanoATAC-seq2技术，在雌性小鼠胚胎中首次实现了对Xist和Tsix等关键调控基因的染色质可及性的等位基因特异性分析。Xist是X染色体失活的关键调控基因，其表达促进所在X染色体的失活；Tsix是Xist的反义链转录本，与Xist竞争结合Xite、Jpx、Six、Tsx和Linx等重要调控基因，其表达主要抑制所在X染色体的失活，即促进或维持该条X染色体的激活状态。
利用两种近交系小鼠交配产生的杂合胚胎进行分析，该研究发现，在4细胞胚胎阶段开始的父本X染色体印记失活过程中，Xist结构域特异性增强父本X染色体在该区域的染色质开放性，从而促进整条父本X染色体的印记失活。而此时Tsix结构域保持父本和母本X染色体在该区域的对称开放性，不参与父本X染色体的特异性印记失活（图2）。发育到早期桑椹胚阶段时，这一状态发生逆转，Xist结构域逐渐转变为父本和母本X染色体在该区域的对称开放性，而Tsix结构域逐渐特异性增强母本X染色体在该区域的染色质开放性，保持母本X染色体处于活跃状态，而仍然使整条父本X染色体维持印记失活状态（图2）。此后，在晚期囊胚阶段，在两个胚外谱系的细胞（滋养外胚层和原始内胚层）中保持这一由Tsix结构域主导的父本X染色体印记失活状态。而在多能性的上胚层细胞中Tsix结构域重新逆转回父本和母本X染色体在该区域的对称开放性，印记失活的父本X染色体被重新激活，但是其激活程度仍然略低于母本X染色体。
图2. Xist与Tsix结构域调控父本X染色体印记失活和重新激活
2. 鉴定了调控合子基因组激活和早期胚胎两次谱系命运决定的关键转录因子
对于结合基序已知的转录因子，通过分析其结合位点的染色质状态，ATAC-seq可以准确鉴定每个转录因子的调控活性，准确判断对应关键生物学事件的主导转录因子。通过对从受精卵到晚期囊胚的十个关键发育阶段的分析，包括受精卵、早期2细胞、晚期2细胞、4细胞、早期8细胞、晚期8细胞、16细胞（早期桑葚胚）、晚期桑椹胚期、早期囊胚和晚期囊胚，该研究准确鉴定了该过程中每个发育阶段和每个分化谱系的关键转录因子（图3）。scNanoATAC-seq2分析提示，这些转录因子的基因体和转录因子下游结合位点上的表观激活程度在发育过程中协同变化。这为后续通过基因敲除等手段对每个主导转录因子的功能进行精细分析奠定了基础，也为研究着床前胚胎的每个发育阶段的核心转录因子网络提供了表观组学依据。
图3.调控小鼠着床前胚胎各阶段发育的主导转录因子
3. 揭示了主要重复元件类别在早期胚胎发育过程中的表观调控特点
哺乳动物基因组中大约一半都是由重复元件组成的。一些重要类型的重复元件在基因组中有数千个序列几乎一模一样的拷贝，分散分布在整个基因组中。例如全长LINE1重复元件拷贝的长度在6 kb左右，同一亚家族内两个拷贝之间的序列相似度能够达到99%以上。二代短读段测序很难区分这些不同的拷贝，因而无法判断其染色质状态（图4）。该研究发现大部分全长LINE1在卵裂期染色质开放度下降，而后在囊胚期上升。然而，存在157个调控模式相反的全长LINE1拷贝，其染色质开放度在卵裂期上升，而后在囊胚期下降，可能涉及常染色质-异染色质的区室形成。
MERVL属于长末端重复序列（long terminal repeat, LTR）家族，全长约6kb，中间为内源逆转录病毒（endogenous retrovirus, ERV）蛋白编码序列，两端为辅助转座和顺式调控的MT2-mm元件。先前的研究表明，这类转座子在合子基因组激活过程中被广泛激活并转录，人为沉默MERVL会导致小鼠着床前发育异常（Sakashita et al., 2023）。由于小鼠基因组中570个完整MERVL元件的序列高度相似，短读段测序技术难以实现唯一比对。利用5 kb以上的长读段（涵盖MERVL侧翼的非重复序列），scNanoATAC-seq2能够实现MERVL来源DNA片段的唯一比对，从而精确描绘每个拷贝的MERVL的表观激活特征。在晚期2细胞阶段发生的合子基因组激活中，MERVL的表观激活发生于两端的MT2-mm重复元件，而内部的ERV编码序列所在的染色质保持相对关闭状态。此外，逐个拷贝分析的染色质可及性表明，这些序列相似的MERVL元件之间仍存在表观激活程度的差异。在scNanoATAC-seq2鉴定的活跃MERVL元件附近（30 kb以内），相应的靶基因呈现ZGA表达激活（如Zscan4c、Zscan4d和Sp110）；相反，表观遗传沉默的MERVL元件附近基因的表达则在ZGA过程中未激活。
图4. scNanoATAC-seq2鉴定序列高度相似的不同拷贝重复元件的染色质状态
4. 鉴定了一批新的小鼠早期胚胎非经典印记基因
scNanoATAC-seq2技术能够高效区分父、母本等位基因，适合研究发育中的基因印记。在C57×DBA和C57×CAST杂合胚胎中，该研究分别识别了47%和97%的染色质可及性的亲本来源。利用这些数据，该研究描绘了染色质可及性介导的非经典印记的动态变化（不依赖DNA甲基化）。结果显示，非经典印记强度在着床前发育过程中逐渐消失，并在晚期囊胚谱系分化后降至最低水平。该研究最晚在8细胞胚胎阶段观测到较强的等位基因特异性染色质可及性。由此，通过反交胚胎验证，研究人员在该发育阶段鉴定出325个小鼠非经典印记基因。其中266个基因是相较于Inoue等人的研究新鉴定出来的父源非经典印记基因。Inoue等人的研究是基于孤雌和孤雄小鼠桑葚胚的染色质可及性差异，而不是基于正常胚胎发育过程的等位基因特异性染色质可及性分析得出的（Inoue et al., 2017）（图5）。
图5. scNanoATAC-seq2鉴定小鼠早期胚胎非经典印记基因
5. 发现16细胞胚胎阶段已经出现胚内和胚外谱系分化的染色质可及性特征
该研究鉴定了早期囊胚中内细胞团和滋养外胚层谱系特异性的染色质开放区域特征，并且以此为基础解析其他发育阶段的胚内异质性。结果表明，小鼠胚胎发育至16细胞阶段时，首次出现了表观基因组层面的内细胞团和滋养外胚层命运分化特征（图6）。这提示在囊胚腔出现之前，囊胚中两种谱系的分化特征已然在16细胞胚胎（早期桑葚胚）阶段的染色质可及性层面出现。
图6. 小鼠胚胎发育16细胞期出现细胞谱系分化的染色质可及性特征
总之，该研究利用单细胞起始的scNanoATAC-seq2技术，构建了涵盖小鼠着床前发育全过程的高精度单细胞染色质可及性图谱，系统鉴定了不同发育阶段、不同谱系细胞的关键转录因子，揭示了调控合子基因组激活和谱系分化（上胚层、原始内胚层和滋养外胚层）的顺式调控网络。此外，还解析了雌性胚胎中父本X染色体的印记失活与重新激活，以及非经典印记基因的表观调控基础。这些发现为理解包括人类在内的哺乳动物早期胚胎发育的分子机制提供了新的线索。
北京大学生物医学前沿创新中心汤富酬研究员、文路副研究员与清华大学基础医学院纪家葵研究员为该论文的共同通讯作者。清华大学基础医学院博士生李孟瑶博士（已毕业）与北京大学生物医学前沿创新中心博士生蒋振寰为该论文的并列第一作者。
文章链接：
https://www.science.org/doi/10.1126/science.adp4319
来源：北京大学生物医学前沿创新中心
