图片设计源于“八卦”，在此代表人类来自精子的父源基因组与来自卵细胞的母源基因组在受精后的“阴阳融合”。整幅画面传递着人类受精卵在发育过程中，父源基因组（金色蛇）染色质快速打开并且其开放程度迅速超过母源基因组（红色蛇），并将这一逆转的开放程度不对称状态一直维持到4-细胞胚胎阶段这一含义。

2018年6月18日，北京大学北京未来基因诊断高精尖创新中心、生命科学学院生物动态光学成像中心汤富酬课题组携手北京大学第三医院乔杰课题组合作在国际知名学术期刊《Nature Cell Biology》在线发表了题为"Single-cell Multi-omics Sequencing of Human Early Embryos"的研究论文。利用汤富酬课题组发展的高精度的单细胞多组学测序技术“single-cell COOL-seq”(Cell Research, 2017)，首次在单细胞分辨率，解析了人类着床前胚胎发育过程中DNA甲基化组和染色质状态组的重编程过程，以及染色质状态与DNA甲基化之间的相互关系等关键生物学特征。

汤富酬教授课题组和乔杰教授课题组长期紧密合作，致力于人类着床前胚胎发育过程中基因表达的表观遗传学调控机理方面的研究。采用汤富酬课题组发展的单细胞RNA-seq转录组测序技术，于2013年绘制了完整的人类着床前胚胎的高精度单细胞转录组图谱ADDIN EN.CITEADDIN EN.CITE.DATA(Nature Structural & Molecular Biology, 2013)；之后利用汤富酬课题组发展的微量细胞DNA甲基化组高通量测序技术，于2014年在国际上首次实现了对人类早期胚胎发育过程中DNA甲基化组重编程的系统研究ADDIN EN.CITEADDIN EN.CITE.DATA(Nature, 2014)。2018年该团队利用单细胞DNA甲基化高通量测序技术，首次在单细胞分辨率深入解析了人类着床前胚胎发育的DNA甲基化组图谱(Nature Genetics, 2018)，揭示了人类着床前胚胎的去甲基化重编程，实际上是高度有序的大规模去甲基化和局部的DNA加甲基化动态平衡的结果，并揭示了父母源基因组差异甲基化等特征。

为了进一步在极限分辨率研究染色质状态在DNA甲基化重编程过程中的动态重构过程，该研究利用汤富酬课题组发展的国际领先的单细胞多组学测序技术（single-cell COOL-seq），在单细胞、单碱基分辨率系统地描绘了人类着床前胚胎发育过程中，各个发育阶段表观基因组多个层面的动态变化。对于存在高比例非整倍体的样品（例如人类着床前胚胎、人类癌症样品等），利用少量细胞的表观基因组测序方法，例如ChIP-seq、ATAC-seq、DNase-seq等，得到的结果有可能会被混杂在样品中的异常的非整体细胞所混淆。因此该研究利用汤富酬课题组发展的scCOOL-seq技术，可以在单细胞水平有效地区分整倍体和非整倍体细胞，在排除非整倍体细胞后、利用整倍体胚胎的单细胞数据，更为精准地反映人类着床前胚胎发育过程中表观基因组多层面的动态变化。该研究的主要发现有：

（1）受精后的19个小时以内，高度特化的精子和卵子都经历了大规模的染色质重构过程，来自精子的父源基因组染色质被迅速打开，而来自卵细胞的母源基因组染色质的开放程度降低（精子的平均染色质开放程度为5%，卵子为42%，受精卵为34%）。随后父母源基因组染色质的开放程度同时回落，直至合子基因激活的8-细胞胚胎阶段之后染色质开放程度再次增加，到桑椹胚阶段染色质开放程度达到最高点（2-细胞的平均染色质开放程度为32%, 4-细胞为30%, 8-细胞为27%，桑椹胚时期为40%；图1a）。利用该研究鉴定出的61,403个近端染色质开放区域进行非监督层次聚类分析发现，和小鼠胚胎相似，在人类早期胚胎发育过程中，近端染色质开放区域具有强烈的发育阶段特异性，更重要的是近端染色质开放区域在合子基因激活时期，即4-细胞到8-细胞时期，发生了最剧烈的染色质重构过程，着床前胚胎聚类分为两支，一支是4-细胞阶段之前的胚胎，另外一支是8-细胞阶段之后的胚胎（图1b）。

图1：人类着床前胚胎（整倍体胚胎）发育过程中染色质开放程度的动态变化（a）以及近端（基因启动子区域）染色质开放区域的非监督层次聚类分析（b）

（2）首次发现不同于小鼠受精卵，人类受精卵在发育过程中父源基因组染色质快速打开并且其开放程度迅速超过母源基因组，并将这一逆转的父母源基因组染色质开放程度不对称状态一直维持到4-细胞胚胎阶段。该研究首次实现了在人类着床前胚胎的单个细胞内，将父母源基因组信息精确分开并进行DNA甲基化和染色质开放程度的分析（图2a）。在此基础上，该研究发现，人类受精卵在受精后19个小时左右，其父源基因组的染色质已经比母源基因组更为开放，并且这一不对称状态一直维持到4-细胞胚胎阶段才结束（图2b）。这与小鼠中完全不同，小鼠胚胎在受精后很快每个胚胎细胞中其父源基因组染色质的开放程度变得跟母源基因组的一样，并且此后一直维持这种父母源基因组染色质开放程度的对称状态。人类受精卵中父源基因组更快速打开并且开放程度迅速超过母源基因组，这可能更有利于DNA去甲基化酶在父源基因组上的结合，从而加速父源基因组的去甲基化进程以及其他表观遗传学重编程进程。

（3）首次发现不同于小鼠胚胎，人类着床前胚胎具有更为开放、松散的染色质结构。该研究首次实现了对人类与小鼠着床前胚胎父母源基因组染色质开放程度的系统定量比较（图2c）。基于同一技术平台的绝对定量比较，该研究发现与小鼠胚胎相比，在对应的胚胎阶段，人类着床前胚胎具有更开放松散的染色质结构，而这一松散的染色质结构有可能更容易出现基因组不稳定性，进而导致人类着床前胚胎发育到囊胚阶段的比率远低于小鼠（人类胚胎的囊胚发育率只有40-60%，而小鼠胚胎的囊胚发育率高达95%以上）。小鼠胚胎的父源基因组和母源基因组在受精后被迅速打开，之后染色质开放程度回落，到2-细胞时期达到最低点；与此不同，人类胚胎母源基因组在受精后其染色质开放程度即开始下降，至8-细胞阶段降至最低点。而人类胚胎父源基因组在受精后先被迅速打开，之后染色质开放程度逐渐回落，到8-细胞阶段降至最低点。人类和小鼠染色质状态重构模式的差异刚好与这两个物种合子基因激活时间的差异吻合。

（4）首次发现受精后，人类雌性胚胎中父源X染色体迅速去甲基化并被激活，到2-细胞阶段后，父源X染色体的DNA甲基化水平已经远低于母源X染色体，并在此后一直维持这一逆转的不对称性。该研究可以精确地区分雌性和雄性胚胎，并在此基础上系统研究了人类雌性着床前胚胎中父母源X染色体的DNA甲基化和染色质开放程度的差异（图2d）。每个雌性细胞中有两条X染色体，其中一条来自精子（父源X染色体），另外一条来自卵细胞（母源X染色体）。scCOOL-seq技术可以基于杂合单核苷酸多态性（SNP）信息在每个单细胞中精确地区分父母源X染色体。该研究表明，与全基因组的DNA甲基化重编程和染色质状态重构特征一致，受精后，父源X染色体迅速去甲基化并被激活，到2-细胞阶段之后，父源X染色体的DNA甲基化水平已经显著低于母源X染色体的甲基化水平，并在此后一直维持这一不对称状态。与此同时，父源X染色体的染色质状态也迅速变得比母源X染色体更为开放，并一直维持到了4-细胞阶段。

图2：人类着床前胚胎发育过程中父母源基因组DNA甲基化与染色质开放程度的不对称分布（a, b），人鼠物种间父母源基因组染色体开放程度的比较（c）以及单个雌性胚胎细胞中父母源DNA甲基化与染色质开放程度的差异（d）

（5）首次在单细胞分辨率，揭示了人类着床前胚胎发育过程中，各种基因组元件DNA去甲基化和染色质状态重构过程的异步性（图3a,b）。该研究发现，人类着床前胚胎发育过程中，基因启动子区域DNA甲基化异质性强烈的基因和染色质状态异质性强烈的基因是两类不同的基因。例如从4-细胞到桑椹胚阶段，启动子区域染色质状态有强烈异质性的基因其启动子DNA甲基化没有显著的异质性（图3c中I类基因），这些基因主要参与染色体组装、染色质的修饰和细胞周期。这暗示了染色质重构相关基因在这些胚胎时期的一部分细胞中活跃表达，而这种表达异质性主要是由这些基因启动子区域染色质状态的异质性驱动的。从受精卵到囊胚阶段，启动子区域DNA甲基化有强烈异质性的基因其启动子染色质状态没有显著的异质性（图3c中II类基因），这些基因主要参与免疫炎症反应和化学刺激的感知。这暗示了免疫炎症反应和化学刺激相关基因的DNA去甲基化在这些时期的胚胎的一部分细胞中还没有完成。

图3：各个元件DNA甲基化异质性（a）和染色质状态异质性（b），以及基因启动子区DNA甲基化异质性和染色质状态异质性的联系（c）

（6）通过转录抑制实验，首次发现，在人类着床前胚胎发育过程中，持续转录对于大量基因维持其启动子区域染色质持续处于开放状态具有重要作用。该研究发现，在人类和小鼠着床前胚胎发育过程中，合子基因激活伴随着启动子区域大片段（长度大于300bp）染色质开放区域比例的大幅度增加（图4）。使用RNA聚合酶抑制剂α-Amanitin处理受精卵、广谱抑制其转录活动，这些胚胎发育到8-细胞阶段时1700多个（35%）基因会失去启动子区域的染色质开放状态。这些依赖于持续转录的启动子大片段染色质开放区域对应的基因，很多都在胚胎发育过程中发挥重要作用，这说明持续转录对于维持大量胚胎发育重要基因启动子区域持续处于开放状态起了关键作用。

图4：持续转录对于染色质开放状态的维持是必须的

（7）对人类着床前胚胎进行了大规模的基因组调控区域预测，共得到37000多个候选增强子区域。利用scCOOL-seq数据，该研究结合染色质开放程度和DNA甲基化信息，进行了人类和小鼠的候选增强子预测。该研究发现对于人类和小鼠着床前胚胎中的染色质处于开放状态的候选增强子，只有700多个是在这两个物种中同时处于染色质开放状态，说明在着床前胚胎发育过程中处于活跃状态的增强子具有很强的物种特异性。与此同时，该研究对人类着床前胚胎中鉴定出的接近40万个染色质开放区域（NDR）进行了转录因子结合富集分析，发现相对于滋养外胚层细胞，多能性的内细胞团细胞中的远端（距离基因转录起始位点2kb以外）染色质开放区域更富集转录因子ZNF281的结合序列；而相对于内细胞团细胞，滋养外胚层细胞中的远端染色质开放区域更富集转录因子CDX2、TFAP2C（AP2γ）的结合序列，这为今后研究人类着床前胚胎发育过程中转录因子调控下游基因的表达提供了基础。

该研究利用国际领先的单细胞多组学高通量测序技术，深度解析了人类着床前胚胎发育过程中，单个整倍体细胞中DNA甲基化组重编程伴随的染色质状态的重构过程，并系统地研究了人类胚胎父母源基因组在单个细胞中DNA甲基化与染色质开放程度的不对称分布。首次定量比较了人类和小鼠胚胎中父母源基因组的染色质开放程度，发现了人类胚胎不同于小鼠胚胎的物种特异性的关键表观遗传学特征，例如：与小鼠胚胎相比，人类着床前胚胎的染色质状态更为松散；不同于小鼠胚胎，人类着床前胚胎中父源基因组染色质状态在卵裂阶段早期比母源基因组更为开放。这为今后人们继续研究人、鼠两个物种早期胚胎中表观遗传学的异同提供了理论基础，明确了利用小鼠作为模式生物研究哺乳动物早期胚胎发育的优势和局限性。同时为研究非整倍体等发育异常的胚胎、以及相关的反复流产等临床医学问题提供了新思路和新研究手段。

北京大学北京未来基因诊断高精尖创新中心博士生李琳、博士生高云，四川大学研究员郭帆，以及北京大学第三医院博士生任一昕为该论文的并列第一作者；汤富酬教授和乔杰教授为该论文的共同通讯作者。该项目得到了国家自然科学基金、国家重大科学研究计划、北京市科学技术委员会、北京未来基因诊断高精尖创新中心（ICG）、生命科学联合中心（CLS）等的资助。

原文链接：https://www.nature.com/articles/s41556-018-0123-2

参考文献：

ADDIN EN.REFLISTGuo, F., Li, L., Li, J., Wu, X., Hu, B., Zhu, P., Wen, L., and Tang, F. (2017). Single-cell multi-omics sequencing of mouse early embryos and embryonic stem cells. Cell research.

Yan, L., Yang, M., Guo, H., Yang, L., Wu, J., Li, R., Liu, P., Lian, Y., Zheng, X., Yan, J., et al.(2013). Single-cell RNA-Seq profiling of human preimplantation embryos and embryonic stem cells. Nat Struct Mol Biol20, 1131-1139.

Guo, H., Zhu, P., Yan, L., Li, R., Hu, B., Lian, Y., Yan, J., Ren, X., Lin, S., Li, J., et al.(2014). The DNA methylation landscape of human early embryos. Nature511, 606-610.

Zhu, P., Guo, H., Ren, Y., Hou, Y., Dong, J., Li, R., Lian, Y., Fan, X., Hu, B., Gao, Y., et al.(2018). Single-cell DNA methylome sequencing of human preimplantation embryos. Nature genetics.

2018年6月4日，北京大学北京未来基因诊断高精尖创新中心、生命科学学院生物动态光学成像中心汤富酬课题组携手中国科学院生物物理研究所王晓群课题组、北京大学第三医院乔杰课题组联合在学术期刊《Cell Research》上在线发表了题为“Spatial transcriptomic survey of human embryonic cerebral cortex by single-cell RNA-seq analysis”的文章。该研究利用单细胞转录组测序技术，绘制了人类胚胎中晚期（22-23周龄）大脑全皮层及脑干共22个脑区的高精度单细胞转录组图谱，系统解析了全皮层细胞类型多样性，揭示了不同脑区细胞组成的差异性以及同种类型神经元在不同脑区的基因差异表达特点，阐明了胚胎期不同脑区发育成熟的异步性特征。这是继2018年3月14日该合作团队在《Nature》发表对人类大脑前额叶发育调控工作之后在空间维度上对人类大脑发育的进一步系统分析。通过与其他课题组已发表的人类成年大脑神经元单细胞转录组数据进行系统的比较分析，确定了胚胎期大脑神经元细胞基因表达的关键特征。此外，该研究还系统分析了已知的自闭症风险基因在各个脑区兴奋性神经元中的表达特点和区域特异性，推测了胚胎发育时期对自闭症发生最密切相关的脑区以及潜在的风险基因，为研究神经系统疾病的发生机理提供了参考。

人类胚胎期大脑全皮层示意图

该研究发现人类大脑皮层的区域差异性首先体现在各种类型细胞的组成比例不同，特别是神经元和胶质细胞组成比例的显著不同上。在额叶区域以兴奋性神经元为主而在颞叶区域以抑制性神经元为主，皮层下侧的脑区则主要富集神经胶质细胞，这为不同脑区神经信号网络的差异提供了细胞基础。兴奋性神经元与起源自原始神经节的抑制性神经元在发育调控上有着截然不同的特征，前者是一个连续成熟的异步分化发生过程，而后者则在发育上比较同步，同一发育时期大脑皮层中不同抑制性神经元的成熟度比较一致，异质性很小。通过对兴奋性神经元细胞发育过程的深入分析，鉴定了调控大脑皮层神经元成熟的关键候选转录因子，为体外神经元定向诱导分化提供了参考。

为了全面解析人类胚胎期神经元细胞与成体神经元细胞的基因表达差异，该研究系统分析比较了多个已经发表的单细胞转录组数据集，发现大脑皮层中的兴奋性神经元在胚胎期与成年期的基因表达差异较小，而抑制性神经元在胚胎期与成年期之间差异很大，说明抑制性神经元在胚胎发育晚期及出生后发生了更为剧烈的基因表达变化。即使胚胎-成体差异程度存在不同，这两种神经元在所研究的胚胎阶段都还不具备成体大脑皮层神经元的分层分布特征。当把所分析的所有类型的神经元细胞做聚类分析时，发现胚胎-成体阶段间的基因表达差异远大于同一发育阶段不同类型神经元细胞之间的基因表达差异，并且胚胎期兴奋性神经元与抑制性神经元之间的基因表达差异大于成年期这二者之间的差异。

兴奋性神经元分化成熟过程的调控、大脑区域成熟差异、以及胚胎期与与成体神经元的比较分析

以往的不同脑区基因表达差异的相关研究，通常是基于大量细胞平均分析的结果来研究脑区差异表达基因，然而这种方式无法确定这样的脑区特异性分子标记具体是哪种细胞表达的，以及这样的区域差异是否仅仅由于不同脑区相同细胞类型的组成比例不同导致的。而该研究在单细胞分辨率上鉴定出的不同神经元细胞类型中比较了不同脑区的差异表达基因，发现了未成熟兴奋性神经元、成熟兴奋性神经元、以及抑制性神经元在不同脑区的基因表达差异，这些区域特异性表达的基因中很多与已知的神经系统疾病或者大脑区域化功能发育基因密切相关，进一步揭示了神经系统功能形成的脑区特异性调控机制。原有的研究发现自闭症的发生与妊娠中期的投射神经元发育关系密切，该研究在此基础上进一步提示已知的自闭症风险基因与较为成熟的兴奋性神经元特异性表达基因的相关性更高，而且这些自闭症风险基因与大脑更靠近背侧的皮层区域的发育更为相关，提示这些背侧皮层区域的神经元发育异常可能更容易引发自闭症。

区域细胞组成差异，区域特异基因验证及自闭症风险基因互作网络

人类大脑皮层包含了数十亿个神经元。近年来利用单细胞测序手段研究神经系统个别脑区多样性的工作已经广泛开展，但是对人类胚胎期的大脑全皮层的各个主要区域进行系统比较分析的研究此前尚未见报导。大脑皮层不同区域的神经元细胞存在明显的功能差异，比如听觉神经主要分布在颞叶区域，枕叶区域是视觉信号传输的主要区域。这些区域功能差异必然与脑区的结构和细胞类型组成密切相关，因而研究发育过程中不同脑区的基因表达差异对于认识大脑区域特异性功能的形成具有重要参考价值。该项研究通过高精度的单细胞转录组测序，结合系统的生物信息学分析和深度数据挖掘，对全皮层的细胞类型和基因表达差异进行了深入、精准的探索和分析，为解答“全皮层区域功能形成的调控”这一关键问题的后续研究提供了高精度的单细胞图谱，也为后续神经系统发育研究提供了参考。

北京大学北京未来基因诊断高精尖创新中心、生命科学学院博士后范小英、博士生董骥、中科院生物物理所博士生钟穗娟，以及北医三院博士生魏源为该研究的并列第一作者，汤富酬教授、乔杰教授及王晓群教授为共同通讯作者。该项工作得到了国家自然科学基金委、科技部、北京未来基因诊断高精尖创新中心、生命科学联合中心等的支持。

文章链接：https://www.nature.com/articles/s41422-018-0053-3

消化系统是人体中最重要的器官系统之一，而消化道作为消化系统最重要的组成部分、其在食物消化、营养吸收、废物排泄、抵抗微生物入侵等多个方面起着极其重要的作用。消化道主要由食道、胃、小肠、大肠组成。揭示这四种器官在人类胚胎发育过程中的基因表达调控网络以及“细胞类型图谱”是发育生物学的重大问题。

2018年5月25日，北京大学北京未来基因诊断高精尖创新中心、生命科学学院生物动态光学成像中心汤富酬课题组、葛颢课题组，携手北京大学第三医院乔杰课题组联合在国际知名学术期刊《Nature Cell Biology》上在线发表了题为“Tracing the temporal-spatial transcriptome landscapes of the human fetal digestive tract using single-cell RNA-sequencing”的研究论文。在国际上首次在单细胞分辨率和全转录组水平，全面、系统、深入地阐明了食道、胃、小肠和大肠这四种器官在人类胚胎发育过程中的基因表达图谱及其信号调控机制，揭示了这四种器官不同细胞类型之间的精准发育路径和基因表达特征，并进一步详细解析了大肠从胎儿到成人的发育、成熟路径和关键生物学特征。

本项研究首先分离并获取了人类6到25周胚胎中食道、胃、小肠和大肠以及成年人的大肠总计5,277个单细胞，并对其进行了高精度的单细胞转录组测序，获得了多达5T的基因表达谱数据。针对这些数据该团队结合t-SNE，KNN，PCA，Monocle2和StemID等多种生物信息学分析方法进行深度数据挖掘，全面系统地对这四种器官的发育进程在基因表达调控层面进行了深入的解析。主要发现有：

（1）通过关键基因表达特征定义了这四种器官共40种不同的重要细胞类型，其中在人类胚胎发育阶段鉴定出30种不同的细胞类型。通过对胚胎发育时期的这30种细胞类型的关键生物学特征，特别是其独特发育路径，交互信号通路网络，细胞周期分布，营养代谢状态和营养吸收特点，以及层级转录因子调控网络等特征进行了细致翔实的分析后发现，其中具有干、祖细胞特征的细胞类型有19种，并且与其它11种已经分化成熟的细胞类型相比，其在胚胎发育的较早期就已经出现并且贯穿于整个胚胎发育时期，其细胞周期也呈现出活跃分裂增殖的特点。另一方面，整体的基因表达聚类分析显示出胚胎早期消化道的四种不同器官的细胞虽然已经具有明确的空间分布差异和不同的形态特征，但是仍然具有高度相似的基因表达特点，而胚胎晚期这四种不同器官的细胞基因表达差异非常明显，显示这些细胞类型功能成熟后基因表达的高度器官特异性和细胞类型特异性。同时小肠同其它三种消化道器官相比，营养吸收代谢的关键基因在胚胎晚期被大规模同步激活，说明处于胚胎晚期的小肠可能已经具备了基本的营养吸收功能。这也说明随着胚胎发育的进行，四种器官在功能上的差异性也逐渐展现出来。

人类胚胎发育过程中食道、胃、小肠、大肠的细胞类型图谱及其关键生物学特征

（2）通过胚胎发育的时间轴进行纵向分析，该研究发现了一系列在胚胎发育早期（6W-11W）、中期（12W-18W）和晚期（19W-25W）阶段特异性表达的标志基因，这为多能性干细胞体外分化形成各种有功能的消化道细胞类型提供了重要的候选标志基因以及基因表达谱的金标准。进一步的Gene Ontology分析发现在胚胎发育早期的消化道是以干、祖细胞的快速增殖和保持未分化状态为主要特征，而胚胎发育中期的消化道则初步建立了消化功能和免疫功能，提示在多能性干细胞体外定向分化过程中要密切关注免疫细胞对于消化道细胞的分化成熟可能发挥的重要调控作用。胚胎发育晚期的消化道则主要以组织细胞的结构特化和消化吸收功能特化和成熟为主要特征。另一方面，通过这四个器官的空间轴进行横向比较分析，发现在胚胎发育早期，四个消化道器官之间差异很小（不同器官之间差异表达的基因只有259个），来自不同器官的细胞其基本基因表达特征仍然非常相似，而到了胚胎发育晚期，四个消化道器官之间差异巨大（不同器官之间差异表达的基因增加到1193个），这也反映出消化道不同器官在发育过程中的功能特化和细胞类型差异不断积累的特点。

（3）在消化道四种器官的发育过程中，通过解析单细胞水平上的基因表达调控网络发现，对于发育起核心调控作用的HOX家族的基因，不同的HOX基因在不同的消化道器官中可能起着关键的区域化调节作用。在Hedgehog信号通路中，IHH配体可能起到了关键作用，并且发现间充质与上皮细胞之间存在旁分泌作用，同时在小肠和大肠还发现Wnt，FGF，TGF-β和BMP等信号通路协同调控着它们的发育过程。

人类胚胎发育过程中小肠的信号通路等调控模式示意图

（4）针对大肠的发育进程，该研究进一步获取了成年人大肠高精度的单细胞转录组数据，首次系统地研究了大肠从胎儿到成人的发育过程，成年人大肠中10种主要的细胞类型都被精确地鉴定出来。该研究系统地揭示了包括大肠干、祖细胞，肠上皮细胞，杯状细胞等功能成熟的细胞类型的核心基因表达模式和关键标志基因等重要生物学特征。在胎儿发育的后期，大肠中的细胞类型以及其基因表达模式与成年人阶段已经很接近，其中一些特定细胞类型的关键基因的表达也都趋于一致，说明在胚胎晚期大肠的主要细胞类型可能已经具备了初步的食物消化和营养吸收功能。

成年人大肠中的细胞类型（左图）及大肠从胚胎时期到成年人阶段的发育路径（右图）

该项研究首次系统、全面、深入地同时将消化系统最重要的食道，胃，小肠，大肠四个主要器官在单细胞转录组水平上进行了时空交叉的深度解析。这四种器官，尤其是小肠，一直被作为研究细胞类型之间的相互作用以及成体干细胞的干性维持和分化的经典模型。本研究的发现无疑为消化道生物学研究领域提供了全面翔实的发育细胞图谱数据，具有重大参考价值，同时也为体外多能性干细胞定向分化为消化道相关细胞类型建立了基因表达图谱金标准。另一方面，对于这四种消化道器官相关的癌症，特别是胃癌和结直肠癌的致病机理研究，此项研究结果也在基因表达层面上提供了重要参考。

北京大学北京未来基因诊断高精尖创新中心的博士后高帅、博士生汪睿，博士生李静云、以及北医三院闫丽盈博士为该文章并列第一作者，汤富酬教授、乔杰教授、及葛颢教授为该文章的共同通讯作者。此项工作得到了国家自然科学基金委、科技部、未来基因诊断高精尖创新中心（ICG）、以及生命科学联合中心（CLS）的支持。

原文链接：https://www.nature.com/articles/s41556-018-0105-4