Science∣汤富酬研究组与合作者揭示人类结直肠癌发病新机制
2018年11月30日，北大-清华生命科学联合中心（CLS）、北京未来基因诊断高精尖创新中心（ICG）汤富酬研究组，与北医三院付卫研究组、乔杰研究组(北大-清华生命科学联合中心），在《Science》期刊发表了题为“Single-cell Multi-omics Sequencing and Analyses of Human Colorectal Cancer”的研究论文。该研究在国际上首次从单细胞分辨率、多组学水平深入解析了人类结直肠癌在发生和转移过程中，基因组拷贝数变异、DNA甲基化异常及基因表达改变的特点及相互关系。
结直肠癌是发病率和死亡率很高的一种癌症。基因组不稳定性、表观遗传学异常及基因表达紊乱等是结直肠癌典型的分子特征。肿瘤组织具有很强的异质性，以往在大量细胞水平的研究只能反映肿瘤组织的系综平均情况，而无法精准区分不同的肿瘤亚克隆以及肿瘤微环境中的其他细胞类型。一般的单细胞单组学测序技术只能分别研究单细胞基因组、单细胞DNA甲基化组、或者单细胞转录组，不能做到在一个单细胞中同时平行分析多个组学层面。汤富酬研究组对该课题组已发表的scTrio-seq（single-cell triple omics sequencing technique）技术进行了大幅度的优化，极大地提升了对单个细胞转录组的检测效率和DNA甲基化位点（CpG位点）的覆盖度，达到与利用完整单细胞进行高精度转录组分析、或者利用完整单细胞进行高精度DNA甲基化组分析相近的效果。在HEp-2细胞系中，单个完整细胞的高精度转录组测序平均可以检测到9,158个基因的表达，而scTrio-seq2的转录组部分平均可以检测到8,106（88%）个基因的表达；DNA甲基化组方面，scTrio-seq2的DNA甲基化组部分和单个完整细胞的高精度DNA甲基化组测序（scBS-seq）均可以检测到1600多万个CpG位点。该研究对患有结直肠癌并伴随转移的12例患者的原发位、淋巴结转移位、肝转移位等样品进行了多点取样，对1,800多个细胞进行了高精度的单细胞多组学测序（scTrio-seq2），共产生了7.6 Tb高质量的测序数据。scTrio-seq2的DNA甲基化组部分，平均每个细胞的测序量为4.1 Gb，平均覆盖到全基因组范围的870多万个CpG位点；scTrio-seq2的转录组部分，平均每个细胞测序量为0.9 Gb，平均可覆盖到3700多个基因。
图1. 课题研究思路及多点取样示意图
该研究的主要发现有：
（1）基于DNA甲基化组数据推断出的基因组拷贝数变异精准鉴定了肿瘤细胞的亚克隆。该研究精准鉴定出了结直肠癌患者单个癌细胞的染色体拷贝数变异谱，并利用其高精度的染色体内断点信息进行了肿瘤细胞的谱系追踪。对于其中5例患者，在每个患者肿瘤组织中鉴定出了多个（6-13个）肿瘤细胞亚克隆，发现原发位肿瘤的亚克隆结构通常比转移位肿瘤更复杂。
图2. CRC01患者癌细胞的单细胞染色体拷贝数变异谱
图3. CRC01患者癌细胞的单细胞染色体拷贝数变异的断点
图4. CRC01患者的肿瘤亚克隆结构示意图
（2）系统解析了肿瘤细胞的DNA甲基化异质性。在全基因组水平，检测到肿瘤细胞的DNA甲基化水平普遍低于癌旁的正常上皮细胞，而且DNA甲基化降低的程度存在很大的异质性。与癌旁的正常上皮细胞相比，单个肿瘤细胞全基因组DNA甲基化水平可降低7-36%。肿瘤细胞特异性DNA低甲基化基因组区域显著富集在LTR、LINE等基因组重复序列和异染色质区域，而肿瘤细胞特异性DNA高甲基化的基因组区域显著富集在CpG岛和常染色质区域。
图5. 不同结直肠癌患者及取样部位的单细胞DNA甲基化水平
（3）阐明了DNA甲基化谱和遗传谱系的相互关系。利用鉴定出的肿瘤细胞遗传谱系信息，发现肿瘤细胞的基因组DNA甲基化谱与遗传谱系高度一致，同一块肿瘤组织中同一谱系的肿瘤细胞DNA甲基化水平高度相似，而不同谱系的肿瘤细胞DNA甲基化水平可存在很大差异。所以肿瘤组织基因组DNA甲基化的强烈异质性主要来自同一个患者肿瘤内不同亚克隆之间的DNA甲基化差异，而不是同一个亚克隆内部不同细胞间的DNA甲基化差异。
图6. CRC01患者不同遗传谱系的癌细胞的单细胞DNA甲基化水平
（4）阐释了DNA甲基化谱和基因表达谱的相互关系。在单细胞水平，基因启动子区域的DNA甲基化与相应基因的表达呈现强烈的负相关关系，而基因区的DNA甲基化与相应基因的表达呈现正相关关系。对于测序单细胞数量较多的4例结直肠癌患者，其中1例患者肿瘤细胞的转录组亚群、遗传亚群（亚克隆）、和DNA甲基化亚群三者之间呈现出一一对应关系；而在另外3例患者中，转录组亚群、遗传亚群（亚克隆）、以及DNA甲基化亚群三者之间的关系较为复杂，没有简单的一一对应关系。
图7. 肿瘤细胞的遗传谱系、DNA甲基化、基因表达三者之间单细胞分辨率的对应关系
（5）首次在单细胞分辨率追踪了同一遗传谱系的肿瘤细胞在转移过程中DNA甲基化及基因表达的变化情况。对于同一个癌症患者同一遗传谱系的肿瘤细胞在从原发灶到转移灶的转移过程中其全基因组DNA甲基化水平相对稳定，而基因组局部区域的DNA甲基化水平可能发生较大变化。上皮-间充质转化相关的重要基因在转移前后的肿瘤细胞中其RNA表达水平和DNA甲基化水平均没有发生显著改变。
图8. 同一遗传谱系的肿瘤细胞在转移过程中的单细胞分辨率的DNA甲基化水平变化
（6）揭示了结直肠癌细胞DNA去甲基化的共同特点。在全部进行了单细胞基因组DNA甲基化测序的10个患者中，肿瘤细胞的DNA去甲基化程度与相应基因组区域在癌旁正常上皮细胞中的DNA甲基化水平呈现正相关关系，也就是说，在正常上皮细胞中DNA甲基化越高的基因组区域，其DNA去甲基化越强烈（DNA甲基化降低的绝对值和降低的相对比例都更大）。DNA去甲基化的程度与基因组重复序列L1以及癌旁正常组织中H3K9me3标记的异染色质区域的分布密度呈现强烈的正相关关系，而与开放的染色质区域以及H3K4me3标记的基因组区域呈现强烈的负相关关系。与正常的早期胚胎发育过程中的两次大规模基因组DNA去甲基化相反，相比于L2基因组重复序列，进化上更为年轻、更为活跃的L1基因组重复序列经历了更强烈的DNA去甲基化。
图9. CRC01患者癌细胞相比于癌旁细胞的DNA甲基化水平
