删除或更新信息,请邮件至freekaoyan#163.com(#换成@)

联合应用RNA-seq和ATAC-seq寻找FOXQ1转录因子的下游靶基因

本站小编 哈尔滨工业大学/2020-03-06

联合应用RNA-seq和ATAC-seq寻找FOXQ1转录因子的下游靶基因

伊梦杰1,2,张锦锦3,郭强2,唐慧2

(1.昆明理工大学 医学院, 昆明 650504;2.云南省第一人民医院,昆明理工大学附属医院,云南省消化内镜临床医学中心,云南省临床病毒学重点实验室,昆明市肿瘤分子与免疫防治重点实验室,昆明 650032;3.中国科学院昆明动物研究所,昆明 650223)



摘要:

结直肠癌(Colorectal cancer, CRC)是一种全球高发的恶性肿瘤,发病原因复杂且预后较差。近年来发现叉头框Q1(Forkhead box Q1, FOXQ1)基因作为一类核转录因子在结直肠癌中高表达,可控制下游基因转录活性。本实验拟探究CRC细胞中FOXQ1的转录调控功能并寻找其下游基因。方法:(1)构建低表达FOXQ1基因的稳定转染CRC细胞株;(2)应用RNA-seq检测FOXQ1敲低前后表达量显著差异的基因;(3)应用转座酶可接近性核染色质区域测序分析(Assay for Transposase-Accessible Chromatin using sequencing, ATAC-seq)检测FOXQ1敲低前后细胞染色质易接近性的变化;(4)进一步对FOXQ1敲低前后的RNA-seq和ATAC-seq数据进行一系列生物信息学分析,寻找CRC中FOXQ1转录调控的潜在下游基因。结果:应用RNA-seq筛选出了敲低FOXQ1后表达显著差异的基因EI24、TLR2、SMAD3,通过联合分析两细胞系的测序结果,发现FOXQ1基因敲低后,在DLD1和SW480两个细胞系中染色质易接近性均增强且表达量均上调的基因有61个,染色质易接近性均减弱且表达量均下调的基因有70个,且EI24、TLR2、SMAD3基因均位于重叠分析结果中,其中TLR2、SMAD3基因的染色质区域有明显变化,而EI24基因的染色质区域变化不明显。通过代谢通路分析找到了EI24、TLR2、SMAD3基因所富集的代谢通路。其中SMAD3、TLR2基因在炎症性肠病(Inflammatory bowel disease , IBD)通路中显著富集。EI24基因在p53信号通路(p53 signaling pathway)通路中显著富集。结论:基于染色质易接近性的变化和转录水平的研究发现:敲低FOXQ1基因对CRC细胞系中染色质的开放情况有较大的影响,且影响FOXQ1转录调控的下游基因的表达。找到了FOXQ1敲低后在SW480、DLD1中均发生变化的基因,为丰富FOXQ1转录因子的下游调控网络提供了研究基础。

关键词:  ATAC-seq  RNA-seq  结直肠癌  靶基因

DOI:10.12113/j.issn.1672-5565.201906003

分类号:R735.3

文献标识码:A

基金项目:国家自然科学基金项目(No.6,2);云南省卫生和计划生育委员会医学学科带头人培养基金(No.D-201642);云南省临床病毒学重点实验室基金(No.2018DG010),昆明市肿瘤分子与免疫防治重点实验室基金(No.2018-1-A-17334).



Combining ATAC-seq and RNA-seq to find downstream target genes of FOXQ1 transcription factor

YI Mengjie1,2, ZHANG Jinjin3, GUO Qiang2, TANG Hui2

(1. Medical School, Kunming University of Science and Technology, Kunming 650504, China; 2. Kunming Key Laboratory of CancerMolecular and Immunological Control, Yunnan Provincial Key Laboratory of Virology, Institute of Basic Medical Sciences, Affiliated Hospital of Kunming University of Science and Technology , First Peoples Hospital of Yunnan Province, Kunming 650032, China; 3. Kunming Institute of Zoology. CAS, Kunming 650223, China)

Abstract:

Colorectal cancer (CRC) is a kind of malignant tumor with high incidence in the world, which has complicated causes and poor prognosis. In recent years, Forkhead box Q1 (FOXQ1) gene has been found to be highly expressed in CRC as a nuclear transcription factor, which can control the transcriptional activity of downstream genes. This study aims aims to explore the transcriptional regulation function of FOXQ1in CRC cells and search for its downstream genes. Methods of the study include 1)constructing constructing stable transfected CRC cell lines with knockdown FOXQ1gene, 2) applying RNA-seq to detect genes with significant differences in expression level before and after FOXQ1, 3) assaying for the Transposase-Accessible Chromatin using sequencing (ATAC-seq) for the determination of Chromatin accessibility in cells before and after FOXQ1knockdown; and 4) conducting a series of bioinformatics analysis on the RNA-seq and ATAC-seq data before and after FOXQ1knockdown to search for potential downstream genes of FOXQ1in CRC. EI24, TLR2,and SMAD3genes were picked out using RNA-seq, which expressed significant difference. Through joint analysis of two cell lines sequencing results, we found that after knocking down FOXQ1,there there were 61 genes with enhanced chromatin accessibility and up-regulation in both cell lines of DLD1 and SW480, and 70 genes with reduced chromatin accessibility and down-regulated expression.EI24, TLR2,and SMAD3genes were located in the overlapping analysis. The chromatin region of TLR2and SMAD3genes showed significant changes, while that of EI24gene did not change significantly. Metabolic pathways enriched by EI24, TLR2,and SMAD3 geneswere found by metabolic pathway analysis. SMAD3and TLR2genes were significantly enriched in the Inflammatory Bowel Disease (IBD) pathway, and EI24gene was significantly enriched in the p53 signaling pathway. Based on the changes in chromatin accessibility and transcriptional levels, it was found that knocking down the FOXQ1gene had a great effect on the chromatin opening in CRC cell lines and it influenced the expression of downstream genes for FOXQ1transcriptional regulation. Genes that changed in SW480 and DLD1 after FOXQ1knockdown were found, which provides a research basis for enriching the downstream regulatory network of FOXQ1 transcription factors.

Key words:  ATAC-seq  RNA-seq  Colorectal cancer  Target genes


伊梦杰, 张锦锦, 郭强, 唐慧. 联合应用RNA-seq和ATAC-seq寻找FOXQ1转录因子的下游靶基因[J]. 生物信息学, 2019, 17(4): 227-236. DOI: 10.12113/j.issn.1672-5565.201906003.
YI Mengjie, ZHANG Jinjin, GUO Qiang, TANG Hui. Combining ATAC-seq and RNA-seq to find downstream target genes of FOXQ1 transcription factor[J]. Chinese Journal of Bioinformatics, 2019, 17(4): 227-236. DOI: 10.12113/j.issn.1672-5565.201906003.
基金项目 国家自然科学基金项目(No.81502556, 81860522);云南省卫生和计划生育委员会医学学科带头人培养基金(No.D-201642);云南省临床病毒学重点实验室基金(No.2018DG010),昆明市肿瘤分子与免疫防治重点实验室基金(No.2018-1-A-17334) 通信作者 唐慧, 女, 博士, 教授, 硕士生导师, 研究方向:肿瘤免疫.E-mail: htang1122@aliyun.com 作者简介 伊梦杰, 女, 硕士研究生, 研究方向:肿瘤免疫.E-mail:15762182860@163.com 文章历史 收稿日期: 2019-06-12 修回日期: 2019-07-02


Contents            Abstract            Full text            Figures/Tables            PDF


联合应用RNA-seq和ATAC-seq寻找FOXQ1转录因子的下游靶基因
伊梦杰1,2, 张锦锦3, 郭强2, 唐慧2     
1. 昆明理工大学 医学院,昆明 650504;
2. 云南省第一人民医院,昆明理工大学附属医院,云南省消化内镜临床医学中心, 云南省临床病毒学重点实验室,昆明市肿瘤分子与免疫防治重点实验室,昆明 650032;
3. 中国科学院昆明动物研究所,昆明 650223

收稿日期: 2019-06-12; 修回日期: 2019-07-02
基金项目: 国家自然科学基金项目(No.81502556, 81860522);云南省卫生和计划生育委员会医学学科带头人培养基金(No.D-201642);云南省临床病毒学重点实验室基金(No.2018DG010),昆明市肿瘤分子与免疫防治重点实验室基金(No.2018-1-A-17334)
作者简介:伊梦杰, 女, 硕士研究生, 研究方向:肿瘤免疫.E-mail:15762182860@163.com
通信作者: 唐慧, 女, 博士, 教授, 硕士生导师, 研究方向:肿瘤免疫.E-mail: htang1122@aliyun.com


摘要: 结直肠癌(Colorectal cancer, CRC)是一种全球高发的恶性肿瘤,发病原因复杂且预后较差。近年来发现叉头框Q1(Forkhead box Q1, FOXQ1)基因作为一类核转录因子在结直肠癌中高表达,可控制下游基因转录活性。本实验拟探究CRC细胞中FOXQ1的转录调控功能并寻找其下游基因。方法:(1)构建低表达FOXQ1基因的稳定转染CRC细胞株; (2)应用RNA-seq检测FOXQ1敲低前后表达量显著差异的基因; (3)应用转座酶可接近性核染色质区域测序分析(Assay for Transposase-Accessible Chromatin using sequencing, ATAC-seq)检测FOXQ1敲低前后细胞染色质易接近性的变化; (4)进一步对FOXQ1敲低前后的RNA-seq和ATAC-seq数据进行一系列生物信息学分析,寻找CRC中FOXQ1转录调控的潜在下游基因。结果:应用RNA-seq筛选出了敲低FOXQ1后表达显著差异的基因EI24、TLR2、SMAD3,通过联合分析两细胞系的测序结果,发现FOXQ1基因敲低后,在DLD1和SW480两个细胞系中染色质易接近性均增强且表达量均上调的基因有61个,染色质易接近性均减弱且表达量均下调的基因有70个,且EI24、TLR2、SMAD3基因均位于重叠分析结果中,其中TLR2、SMAD3基因的染色质区域有明显变化,而EI24基因的染色质区域变化不明显。通过代谢通路分析找到了EI24、TLR2、SMAD3基因所富集的代谢通路。其中SMAD3、TLR2基因在炎症性肠病(Inflammatory bowel disease, IBD)通路中显著富集。EI24基因在p53信号通路(p53 signaling pathway)通路中显著富集。结论:基于染色质易接近性的变化和转录水平的研究发现:敲低FOXQ1基因对CRC细胞系中染色质的开放情况有较大的影响,且影响FOXQ1转录调控的下游基因的表达。找到了FOXQ1敲低后在SW480、DLD1中均发生变化的基因,为丰富FOXQ1转录因子的下游调控网络提供了研究基础。
关键词: ATAC-seq    RNA-seq    结直肠癌    靶基因    
Combining ATAC-seq and RNA-seq to find downstream target genes of FOXQ1 transcription factor
YI Mengjie1,2, ZHANG Jinjin3, GUO Qiang2, TANG Hui2     
1. Medical School, Kunming University of Science and Technology, Kunming 650504, China;
2. Kunming Key Laboratory of CancerMolecular and Immunological Control, Yunnan Provincial Key Laboratory of Virology, Institute of Basic Medical Sciences, Affiliated Hospital of Kunming University of Science and Technology, First People's Hospital of Yunnan Province, Kunming 650032, China;
3. Kunming Institute of Zoology. CAS, Kunming 650223, China


Abstract: Colorectal cancer (CRC) is a kind of malignant tumor with high incidence in the world, which has complicated causes and poor prognosis. In recent years, Forkhead box Q1 (FOXQ1) gene has been found to be highly expressed in CRC as a nuclear transcription factor, which can control the transcriptional activity of downstream genes. This study aims aims to explore the transcriptional regulation function of FOXQ1 in CRC cells and search for its downstream genes. Methods of the study include 1)constructing constructing stable transfected CRC cell lines with knockdown FOXQ1 gene, 2) applying RNA-seq to detect genes with significant differences in expression level before and after FOXQ1, 3) assaying for the Transposase-Accessible Chromatin using sequencing (ATAC-seq) for the determination of Chromatin accessibility in cells before and after FOXQ1 knockdown; and 4) conducting a series of bioinformatics analysis on the RNA-seq and ATAC-seq data before and after FOXQ1 knockdown to search for potential downstream genes of FOXQ1 in CRC. EI24, TLR2, and SMAD3 genes were picked out using RNA-seq, which expressed significant difference. Through joint analysis of two cell lines sequencing results, we found that after knocking down FOXQ1, there there were 61 genes with enhanced chromatin accessibility and up-regulation in both cell lines of DLD1 and SW480, and 70 genes with reduced chromatin accessibility and down-regulated expression.EI24, TLR2, and SMAD3 genes were located in the overlapping analysis. The chromatin region of TLR2 and SMAD3 genes showed significant changes, while that of EI24 gene did not change significantly. Metabolic pathways enriched by EI24, TLR2, and SMAD3 genes were found by metabolic pathway analysis. SMAD3 and TLR2 genes were significantly enriched in the Inflammatory Bowel Disease (IBD) pathway, and EI24 gene was significantly enriched in the p53 signaling pathway. Based on the changes in chromatin accessibility and transcriptional levels, it was found that knocking down the FOXQ1 gene had a great effect on the chromatin opening in CRC cell lines and it influenced the expression of downstream genes for FOXQ1 transcriptional regulation. Genes that changed in SW480 and DLD1 after FOXQ1 knockdown were found, which provides a research basis for enriching the downstream regulatory network of FOXQ1 transcription factors.
Key Words: ATAC-seq    RNA-seq    Colorectal cancer    Target genes    
结直肠癌(Colorectal cancer, CRC)是一种由多致病因素导致且预后很差的恶性消化道肿瘤,是世界第三大常见肿瘤,其死亡率在所有癌症中居于第四位(位于肺癌、肝癌、胃癌之后)[1]。叉头框Q1(Forkhead box Q1, FOXQ1)是叉头框(Forkhead box, FOX)基因家族的成员之一,基因定位于6p23-25,编码含403个氨基酸的FOXQl蛋白,作为一类核转录因子,可以稳定结合到靶基因启动子区域的GC盒等核心元件,控制下游基因转录活性从而发挥生物效应[2]。2010年Kaneda等[3]研究发现FOXQ1在CRC中异常高表达,近年来大量研究亦证实FOXQ1在卵巢癌、乳腺癌、膀胱移形细胞癌、胃癌、肝癌、非小细胞肺癌和神经脑胶质瘤等多种肿瘤中异常表达[4-6],与多种肿瘤的发生、发展密切相关,在许多肿瘤中具有明确的促肿瘤生长转移的功能。因此,鉴定新的FOXQ1下游靶基因,丰富由FOXQ1参与并介导的信号通路信息,有望为肿瘤的靶向治疗提供新的靶点。

ATAC-seq是一种检测染色质易接近性的测序方法,可以检测样本间染色质易接近性的变化情况。ATAC-seq只需要很少的细胞(50 000个)就能检测基因组中所有活跃的调控序列[7],DNA探针(作为转座子发挥作用)通过酶促反应(转座酶Tn5)被整合到基因组的开放区域,然后通过测序来鉴定这些区域[8]。本实验联合应用RNA-seq与ATAC-seq寻找由FOXQ1敲低引起的染色质易接近性改变所导致的表达改变基因,即FOXQ1转录因子调控的潜在下游基因。

1 材料和方法1.1 细胞、仪器和试剂本实验所用细胞系均购于中科院上海细胞库,DMEM高糖培养液、RPMI 1640培养液购于美国Corning公司,胎牛血清购于美国Gibco公司。Puromycin购于北京索莱宝科技有限公司,FOXQ1一抗购于abcam公司,HRP-Rb-anti-goat二抗购于Cell Signaling Technology公司,β-actin一抗、HRP-goat-anti-mouse二抗均购于Proteintech公司。RNAzol RT RNA Isolation Reagent购于美国MRC公司,引物由宝生物工程有限公司合成。核酸定量仪为Thermo公司产品,LightCycler 480实时荧光定量仪为Roche公司产品,WB垂直基础电泳仪为伯乐公司产品。

1.2 方法1.2.1 实验分组根据课题组前期工作并结合文献报道,选取在CRC细胞系中FOXQ1表达量较高的SW480和DLD1细胞系进行FOXQ1基因的敲低实验,实验分为两组共4株细胞。FOXQ1基因敲低组:SW480-shFOXQ1和DLD1-shFOXQ1;对照组:SW480-shControl和DLD1-shControl。

1.2.2 FOXQ1基因低表达稳定细胞系的构建1.2.2.1 慢病毒表达质粒的构建、筛选与扩增根据FOXQ1序列,得到siRNA靶序列,设计合成3对shRNA干扰序列(见表 1),命名为shRNA-A、shRNA-B、shRNA-C,并设计一条对照序列shRNA-control;应用Addgene的pSPAX2慢病毒包装系统,PLKO.1-puro-shFOXQ1质粒包装构建获得lenti-shFOXQ1慢病毒,每组重组质粒经测序鉴定无误后,用去内毒素大提试剂盒提取质粒,测浓度;酶切鉴定无误后,在HEK293T细胞中扩增,之后用LB培养基筛选、扩增[9]

表1(Table 1)
表 1 shRNA引物序列Table 1 shRNA primer sequences The name of the primer Base sequence (5'to3')

shRNA-A F CCGGCTCCATCAAACGTGCCTTACTCGAGTAAGGCACGTTTGATGGAGTTTTTG

shRNA-A R AATTCAAAAACTCCATCAAACGTGCCTTACTCGAGTAAGGCACGTTTGATGGAG

shRNA-B F CCGGGCTATTGACCGATGCTTCACTCGAGTGAAGCATCGGTCAATAGCTTTTTG

shRNA-B R AATTCAAAAAGCTATTGACCGATGCTTCATTCTCGAGTGAAGCATCGGTCAATAGC

shRNA-C F CCGGTCTTGAAAGCAAGTGTGATTCTCGAGAATCACACTTGCTTTCAAGTTTTTG

shRNA-C R AATTCAAAAACTTGAAAGCAAGTGTGATTCTCGAGAATCACACTTGCTTTCAAGAA



表 1 shRNA引物序列Table 1 shRNA primer sequences


1.2.2.2 慢病毒的感染及阳性细胞的筛选将待感染的SW480和DLD1细胞铺在24孔板中培养,生长密度达70%~80%时将病毒浓缩液加至细胞中,培养12~15 h后更换为完全培养基。待细胞长满时,将24孔板中的细胞传代至六孔板中进行初步扩大培养,并加入含有1 μg/ml嘌呤霉素的培养基筛选阳性细胞。筛选获得的阳性细胞分别应用qRT-PCR和WB在mRNA和蛋白水平验证FOXQ1基因的敲低效率[10]

1.2.3 RNA-seq1.2.3.1 收集细胞将稳定转染的细胞传代培养稳定生长后,分别收集生长状态良好的DLD1、SW480基因敲低组与对照组细胞进行后续RNA-seq实验。

1.2.3.2 RNA的提取和检测提取基因敲低组与对照组细胞总RNA,用Nanodrop 2000检测所提总RNA的浓度和纯度后用琼脂糖凝胶电泳检测RNA的完整性,最后利用Agilent 2100测定RIN值。

1.2.3.3 RNA的富集和测序使用结合有poly-T寡核苷酸的磁珠从总RNA中分离出含有poly-A的mRNA,加入片段化缓冲剂将其打成片段,将片段化的mRNA逆转录成cDNA并纯化;经末端修复、poly-A添加、测序接头连接及AMPure XP beads筛选后得出大小合适的片段,进行PCR扩增,建立测序文库并进行文库质控[11]

1.2.3.4 上机测试。对构建合格的测序文库进行双末端(Paired-end)测序,本测序实验组和对照组均设置了三次生物学重复,测序工作由上海嘉因生物科技有限公司完成。

1.2.3.5 RNA-seq的数据分析RNA-seq数据下机后首先进行空载去除、接头去除等数据预处理后产出原始数据,紧接着利用FastQC软件进行数据质量控制,利用生物信息学软件STAR、HTSeq和DESeq2对测序结果进行参考序列比对、表达量统计、差异基因筛选等分析[12-14],并通过R软件绘制基因聚类分析图、火山图等。分别建立DLD1和SW480细胞FOXQ1敲低前后的转录谱。

1.2.4 ATAC-seq1.2.4.1 收集细胞分别收集DLD1、SW480基因敲低组和对照组细胞,计数50 000个细胞,离心去上清后依次用预冷的PBS、lysis buffer悬浮细胞500g,4 ℃离心去除上清液,立即进行转座反应[15]

1.2.4.2 转座反应与纯化确保细胞始终置于冰上,配置转座反应体系悬浮细胞,37 ℃孵育30 min,立即用Qiagen MinElute PCR Purification Kit纯化DNA,之后用10μl elution buffer洗脱。

1.2.4.3 PCR扩增配置PCR反应体系循环扩增,之后用Qiagen MinElute PCR Purification Kit纯化DNA[8]

1.2.4.4 上机测试库检合格后,把不同文库按照有效浓度及目标下机数据量的需求合并后进行Illumina HiSeq测序。测序基于边合成边测序(Sequencing by Synthesis)的原理进行,在序的流动池中加入四种荧光标记的dNTP、DNA聚合酶以及接头引物进行扩增,在每一个测序簇延伸互补链时,每加入一个被荧光标记的dNTP就能释放出相对应的荧光,测序仪通过捕获荧光信号,并通过计算机软件将光信号转化为测序峰,从而获得待测片段的序列信息。本测序实验由上海嘉因生物科技有限公司完成。

1.2.4.5 ATAC-seq的数据分析首先对ATAC-seq下机数据进行预处理,如FastQC质量控制,去除duplicate sequences,去除blacklist reads[16],原始序列比对,比对后去除重复序列和细胞器序列等,再利用生物信息学软件BWA、Macs2、Homer、Deep tools等对ATAC-seq测序结果进行参考序列比对分析、差异结合位点检测、特征峰在全基因组上的分布注释等,建立DLD1和SW480细胞FOXQ1敲低前后的染色质易接近性变化图谱。

1.2.5 RNA-seq和ATAC-seq数据的关联分析分别将DLD1和SW480两个细胞系中RNA-seq的差异基因和ATAC-seq的差异表达峰进行关联分析,得到两组关联分析结果,以明确两组细胞中染色质易接近性变化区域对下游基因的调控功能,并找出每种细胞实验组与对照组相比差异染色质易接近性区域可能调控的下游基因。

1.2.6 两个细胞系关联分析数据的重叠分析从DLD1和SW480的关联分析结果中找出两细胞共有的重叠部分,即获得FOXQ1敲低后引起的染色质易接近性改变所导致的表达改变基因。

1.2.7 代谢通路分析将重叠分析结果中所获基因基于代谢通路数据库(Kyoto Encyclopedia of Genes and Genomes, KEGG)进行代谢通路注释,得到注释基因参与的所有代谢通路名称,采用Fisher检验计算代谢通路的显著性水平(P < 0.05),从而筛选出注释基因富集的显著性代谢通路[17]

2 结果分析2.1 在mRNA和蛋白水平验证DLD1和SW480细胞中FOXQ1的敲低效率qRT-PCR和WB实验结果(见图 1)表明成功构建了FOXQ1基因敲低组DLD1-sh-FOXQ1、SW480-sh-FOXQ1以及对照组DLD1-sh-Control和SW480-sh-Control。FOXQ1基因敲低组在mRNA和蛋白水平与对照组相比FOXQ1基因的表达量均显著降低(***P < 0.001)。

图 1(Figure 1)
图 1 稳定低表达FOXQ1的DLD1和SW480细胞验证(***P < 0.001)Figure 1 Verification of the FOXQ1 expression quantity in DLD1 and SW480 cells which knock down FOXQ1 gene steadily (***P < 0.001)


2.2 RNA-seqRNA-seq结果表明,DLD1细胞FOXQ1基因敲低组与对照组相比表达显著上调的基因有215个,表达显著下调的基因有131个;SW480细胞FOXQ1敲低组与对照组相比表达显著上调的基因有171个,表达显著下调的基因有358个(表达差异基因的筛选阈值为FDR < 0.05,log2FC>1)。对这些表达显著差异基因进行分析,发现在FOXQ1基因敲低后,与侵袭、自噬相关的EI24基因在表达下调基因中位居前列,与先天免疫应答和炎症反应相关的TLR2基因以及与迁移相关的SMAD3基因在表达上调的基因中位居前列。推测这些基因的表达改变与CRC的发生、侵袭、发展等密切相关,将进一步结合ATAC-seq分析这几种基因在FOXQ1敲低后,染色质开放区域的变化情况(见图 2)。

图 2(Figure 2)
图 2 RNA-seq基因聚类分析和火山图Figure 2 RNA-seq cluster analysis and volcano map


2.3 ATAC-seqATAC-seq结果表明DLD1-shFOXQ1组有39 146个特征峰,DLD1-shControl组有49 381个特征峰,其中,基因敲低组与对照组相比染色质易接近性增强的差异表达峰有2 385个,染色质易接近性减弱的差异表达峰有6 205个。SW480-shFOXQ1组有38 962个特征峰,SW480-shControl组有42 244个特征峰,其中,基因敲低组与对照组相比染色质易接近性增强的差异表达峰有4 563个,染色质易接近性减弱的差异表达峰有3 733个(见图 3)。

图 3(Figure 3)
图 3 ATAC-seq分析图Figure 3 ATAC-seq analysis注:(a)左侧最高峰代表无核小体片段,对应开放的染色质区域,右侧是单核小体峰,同时具有非核小体片段和单核小体片段说明数据质量良好;(b)回帖序列主要富集在转录起始区(TSS);(c)ATAC-seq序列主要富集在峰中心附近,说明数据良好;(d)将注释结果在全基因组功能性区域上的分布进行统计分析,根据其分布了解结合位点可能具有的功能。


2.4 RNA-seq与ATAC-seq的关联分析本实验将RNA-seq结果中的差异基因和ATAC-seq结果中差异染色质开放区域进行关联分析,预测由FOXQ1基因敲低所导致的染色质开放区域改变引起的转录因子的结合能力的改变,及最终导致下游基因的表达上调或下调。结果显示,DLD1染色质易接近性减弱区域对下游基因表达具有抑制作用(见图 4(a));染色质易接近性增强区域对下游基因表达具有促进作用(见图 4(b))。在SW480细胞中也观察到相同的结果(见图 4(c)、图 4(d))[4]

图 4(Figure 4)
图 4 染色质功能预测Figure 4 Prediction of chromatin function注:红线代表上调基因,紫线代表下调基因,黑色虚线代表没有表达改变的基因作为背景。


2.5 SW480与DLD1关联分析结果的重叠分析将两个细胞系的关联分析结果基因进行重叠分析。重叠部分为在两种细胞中染色质易接近性变化相同且调控的表达量变化趋势相同的基因。其中,DLD1细胞系中染色质易接近性减弱区域调控的下调基因有1 436个,SW480细胞系中染色质易接近性减弱区域调控的下调基因有407个(见图 5),两细胞系染色质易接近性减弱区域共同调控的表达下调基因有70个。DLD1细胞系中染色质易接近性增强区域调控的上调基因有792个,SW480细胞系中染色质易接近性增强区域调控的上调基因有531个(见图 5),两细胞系染色质易接近性增强区域共同调控的表达上调基因有61个。

图 5(Figure 5)
图 5 SW480与DLD1关联分析结果的重叠分析Figure 5 Overlapping analysis of correlation analysis results in SW480 and DLD1


在这两个细胞的交集基因中发现EI24、TLR2、SMAD3基因也在其中,所以观察这些基因在染色质易接近性方面发生的变化(见图 6~图 8)。其中,TLR2、SMAD3基因在ATAC-seq数据中基因敲低组与对照组相比染色质区域有较为明显的变化,而EI24基因的染色质区域变化很弱。

图 6(Figure 6)
图 6 ATAC-seq EI24基因比对Figure 6 ATAC-seq EI24 genome comparison


图 7(Figure 7)
图 7 ATAC-seq TLR2基因比对Figure 7 ATAC-seq TLR2 genome comparison


图 8(Figure 8)
图 8 ATAC-seq SMAD3基因比对Figure 8 ATAC-seq SMAD3 genome comparison


2.6 代谢通路分析将两个细胞系的重叠基因进行代谢通路分析,结果表明染色质易接近性增强区域调控的表达上调基因显著性富集的代谢通路包括炎症性肠病(Inflammatory bowel disease, IBD)、甘油磷脂代谢通路(Glycerophospholipid metabolism)、细胞周期代谢通路(Cell cycle)等。其中SMAD3、TLR2基因显著性富集在IBD代谢通路(见图 9(a))。染色质易接近性减弱区域调控的表达下调基因显著性富集的代谢通路包括p53信号通路(p53 signaling pathway)、TRP通道炎症调控通路(Inflammatory mediator regulation of TRP channels)、半胱氨酸和蛋氨酸代谢通路(Cysteine and methionine metabolism)等。EI24基因显著性富集在p53信号通路(图 9(b))。

图 9(Figure 9)
图 9 代谢通路分析Figure 9 Metabolic pathway analysis注:红色代表具有显著性(P < 0.05)蓝色代表不具有显著性(P > 0.05), 横坐标代表代谢通路名称,纵坐标代表富集程度.


有研究表明,患有IBD的患者与健康人群相比患结直肠癌的风险更高[18-19],当肠道系统与其微生物群之间的关系(包括屏障功能、免疫信号和代谢物)受到干扰后引起的慢性炎症是其发病的主要潜在原因[20]。其次,在CRC中P53信号通路也是一种重要的信号通路,p53抑癌基因突变是导致CRC发生的最主要原因之一,同时也是结直肠癌侵袭和转移的原因之一,还有研究认为p53突变在腺瘤-癌转移过程中也发挥重要作用[21]。进一步研究IBD、P53信号通路在CRC发生发展中的作用将具有重要意义。

3 讨论一般来说,染色质有“关闭”“开放”两种状态,处于“关闭”状态的染色质,在异染色质蛋白以及修饰酶的作用下,被包装成致密、紧凑的结构,阻遏转录因子等蛋白的结合,此时染色质处于沉默失去生物功能的阶段;而处于“开放”状态的染色质,具有不太紧致的结构,可招募转录因子等蛋白的结合,进而调控基因的表达水平。ATAC-seq作为一种绘制全基因组染色质可及性图谱的方法[8],利用超活性Tn5转座酶检测染色质的可接近性,是本实验的重要研究手段之一[7]。通过对ATAC-seq和RNA-seq数据的分析,初步确定了FOXQ1基因敲低后发生差异表达的基因,丰富了FOXQ1转录因子的下游调控网络。


参考文献
[1] INéS M, CRISTINA S D D, PRADILLA D A, et al. Colorectal carcinoma:A general overview and future perspectives in colorectal cancer[J]. International Journal of Molecular Sciences, 2017, 18(1): 197-236. DOI:10.3390/ijms18010197 (0)


[2] BIELLER A, PASCHE B, FRANK S, et al. Isolation and characterization of the human forkhead gene FOXQ1[J]. DNA & Cell Biology, 2001, 20(9): 555-561. DOI:10.1089/104454901317094963 (0)


[3] HIROYASU K, TOKUZO A, KAORU T, et al. FOXQ1 is overexpressed in colorectal cancer and enhances tumorigenicity and tumor growth[J]. Cancer Research, 2010, 70(5): 2053-2063. DOI:10.1158/0008-5472.CAN-09-2161 (0)


[4] WANG S, SUN H, MA J, et al. Target analysis by integration of transcriptome and ChIP-seq data with BETA[J]. Nature Protocols, 2013, 8(12): 2502-2515. DOI:10.1038/nprot.2013.150 (0)


[5] GAO M, SHIH I M, WANG T L. The role of forkhead box Q1 transcription factor in ovarian epithelial carcinomas[J]. International Journal of Molecular Sciences, 2012, 13(11): 13881-13893. DOI:10.3390/ijms131113881 (0)


[6] ZHANG H, MENG F, LIU G, et al. Forkhead transcription factor foxq1 promotes epithelial-mesenchymal transition and breast cancer metastasis[J]. Cancer Research, 2011, 71(4): 1292-1301. DOI:10.1158/0008-5472.CAN-10-2825 (0)


[7] POTT S, LIEB J D. Single-cell ATAC-seq: Strength in numbers[J]. Genome Biology, 2015, 16(1): 172. DOI:10.1186/s13059-015-0737-7 (0)


[8] BUENROSTRO J D, WU B, CHANG H Y, et al. ATAC-seq: A method for assaying chromatin accessibility genome-wide[J]. Current Protocols in Molecular Biology, 2015, 109(1): 21.9.1-21.9.9. DOI:10.1002/0471142727.mb2129s109 (0)


[9] 白璇, 唐慧, 郎丰超, 等. 慢病毒表达载体的构建及沉默FOXQ1基因在大肠癌细胞系DLD-1种的表达[J]. 世界华人消化杂志, 2014, 22(19): 2752-2757.
BAI Xuan, TANG Hui, LANG Fengchao, et al. Construction of lentiviral expression vectors to silence expression of FOXQ1 gene in colorectal cancer cell line DLD-1[J]. Shijie Huaren Xiaohua Zazhi, 2014, 22(19): 2752-2757. DOI:10.11569/wcjd.v22.i19.2752 (0)


[10] 刘捷, 施露露, 刘春燕, 等. 慢病毒载体沉默ADAMTS6人非小细胞肺癌稳转株的构建[J]. 现代生物医学进展, 2019, 19(8): 1401-1406.
LIU Jie, SHI Lulu, LIU Chunyan, et al. Construction of ADAMTS6 stable knockdown non-small cell lung cancer cell line using lentiviral vector[J]. Progress in Modern Biomedicine, 2019, 19(8): 1401-1406. DOI:10.13241/j.cnki.pmb.2019.08.001 (0)


[11] 李露双, 董文慧, 丁兴萃, 等. 麻竹笋转录组测序及苦涩味物质合成基因差异表达分析[J]. 林业科学研究, 2018, 31(4): 38-46.
LI Lushuang, DONG Wenhui, DING Xingcui, et al. Transcriptome sequencing and differential expression analysis of bitter and astringent substances biosynthesis related gene in dendrocalamus latiflorus[J]. Forest Reseach, 2018, 31(4): 38-46. DOI:10.13275/j.cnki.lykxyj.2018.04.006 (0)


[12] CHEN G, SHI T, SHI L. Characterizing and annotating the genome using RNA-seq data[J]. Science China-life Sciences, 2017, 60(2): 116-125. DOI:10.1007/s11427-015-0349-4 (0)


[13] DOBIN A, DAVIS C A, SCHLESINGER F, et al. STAR: Ultrafast universal RNA-seq aligner[J]. Bioinformatics, 2013, 29(1): 15-21. DOI:10.1093/bioinformatics/bts635 (0)


[14] ANDERS S, PYL P T, HUBER W. HTSeq--a Python framework to work with high-throughput sequencing data[J]. Bioinformatics, 2015, 31(2): 166-169. DOI:10.1093/bioinformatics/btu638 (0)


[15] LOVE M I, HUBER W, ANDERS S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2[J]. Genome Biology, 2014, 15(12): 550. DOI:10.1186/s13059-014-0550-8 (0)


[16] CORCES M R, TREVINO A E, HAMILTON E G, et al. An improved ATAC-seq protocol reduces background and enables interrogation of frozen tissues[J]. Nature Methods, 2017, 14(10): 959. DOI:10.1038/nmeth.4396 (0)


[17] GUO Y, YE F, SHENG Q, et al. Three-stage quality control strategies for DNA re-sequencing data[J]. Briefings in bioinformatics, 2014, 15(6): 879-889. DOI:10.1093/bib/bbt069 (0)


[18] DRAGHICI S, KHATRI P, TARCA A L, et al. A systems biology approach for pathway level analysis[J]. Genome Research, 2007, 17(10): 1537. DOI:10.1101/gr.6202607 (0)


[19] PULUSU S S R, LAWRANCE I C. Dysplasia and colorectal cancer surveillance in inflammatory bowel disease[J]. Expert Review Gastroenterology Hepatology, 2017, 11(8): 711-722. DOI:10.1080/17474124.2017.1327347 (0)


[20] YANG Y, JOBIN C. Novel insights into microbiome in colitis and colorectal cancer[J]. Current Opinion in Gastroenterology, 2017, 33(6): 422-427. DOI:10.1097/MOG.0000000000000399 (0)


[21] LI X L, ZHOU J, CHEN Z R, et al. P53 mutations in colorectal cancer-molecular pathogenesis and pharmacological reactivation[J]. World Journal of Gastroenterology, 2015, 21(1): 84-93. DOI:10.3748/wjg.v21.i1.84 (0)



相关话题/基因 临床 基金 免疫 细胞

  • 领限时大额优惠券,享本站正版考研考试资料!
    大额优惠券
    优惠券领取后72小时内有效,10万种最新考研考试考证类电子打印资料任你选。涵盖全国500余所院校考研专业课、200多种职业资格考试、1100多种经典教材,产品类型包含电子书、题库、全套资料以及视频,无论您是考研复习、考证刷题,还是考前冲刺等,不同类型的产品可满足您学习上的不同需求。 ...
    本站小编 Free壹佰分学习网 2022-09-19
  • 合成生物学基因设计软件:iGEM设计综述
    合成生物学基因设计软件:iGEM设计综述伍克煜1,刘峰江1,许浩1,张浩天1,王贝贝1,2(1.电子科技大学生命科学与技术学院,成都611731;2.电子科技大学信息生物学研究中心,成都611731)摘要:随着基因回路规模的扩大,和应用范围的拓展,传统的合成基因回路的设计思路面临着新的挑战。新合成基 ...
    本站小编 哈尔滨工业大学 2020-12-05
  • 蒺藜苜蓿多聚半乳糖醛酸酶基因家族的全基因组分析
    蒺藜苜蓿多聚半乳糖醛酸酶基因家族的全基因组分析胡龄予,伊霖晟,刘冬雪,王阳,杨欢欢,张鑫茹,杨爽,束永俊(哈尔滨师范大学生命科学与技术学院,哈尔滨150025)摘要:多聚半乳糖醛酸酶(Polygalacturonases,PGs)是一种果胶水解酶,参与果实成熟、器官脱落、花粉成熟等多个植物发育过程。 ...
    本站小编 哈尔滨工业大学 2020-12-05
  • 生物信息学方法筛选胶质母细胞瘤的核心基因
    生物信息学方法筛选胶质母细胞瘤的核心基因柯帅,王文波,廖红展,彭志柱,邱县生,唐慧敏,李清华,夏学巍(桂林医学院附属医院神经外科,广西桂林541000)摘要:多形性胶质母细胞瘤(GBM)是成人最常见的恶性神经上皮肿瘤,关于其诊断和治疗的靶点研究一直是困扰研究者的难题。采用生物信息学的方法对GBM的基 ...
    本站小编 哈尔滨工业大学 2020-12-05
  • 医学遗传学2.0:导致人类慢病的主因可能首先是人体共生微生物基因异常,其次才是人类基因异常
    医学遗传学2.0:导致人类慢病的主因可能首先是人体共生微生物基因异常,其次才是人类基因异常张成岗(军事科学院军事医学研究院辐射医学研究所,全军军事认知与心理卫生研究中心,北京100850)摘要:当前慢病高发的现实对“健康中国2030”战略目标的实现提出了巨大挑战。虽然众多医疗机构和政府管理部门付出巨 ...
    本站小编 哈尔滨工业大学 2020-12-05
  • Hsa-miR-210-5p靶基因预测及其相关信号通路的生物信息学分析
    Hsa-miR-210-5p靶基因预测及其相关信号通路的生物信息学分析蔡丹平,龙鼎新(南华大学船山学院,湖南衡阳421001)摘要:为深入研究miR-210-5p的调控机制及生物学功能提供理论机制,应用生物信息学方法分析miR-210-5p序列,预测其靶基因,用Veney2.1.0绘制韦恩图得到靶基 ...
    本站小编 哈尔滨工业大学 2020-12-05
  • 大鼠颅脑损伤后差异表达基因及miRNA研究
    大鼠颅脑损伤后差异表达基因及miRNA研究王家昕,王晓霞,李洁,苗泽远,倪爽,王子钰,苏立宁(河北北方学院基础医学院,河北张家口075000)摘要:研究背景创伤性脑损伤(Traumaticbraininjury,TBI)是致死率和致残率极高的外科疾患,我国在对于TBI的判断、治疗等方面还处于薄弱阶段 ...
    本站小编 哈尔滨工业大学 2020-12-05
  • IgA 肾病患者与膜性肾病患者外周血单核细胞mRNA分析
    IgA肾病患者与膜性肾病患者外周血单核细胞mRNA分析梁爽,凡奎,张燕,谢杨眉(四川省三台县人民医院肾内科,四川绵阳621100)摘要:为了寻找诊断、鉴别IgA肾病(IgAN)和膜性肾病(MN)的血液特异性标记物,利用公共数据库中的IgAN和MN患者的外周血单核细胞(PBMCs)的转录组表达谱数据集 ...
    本站小编 哈尔滨工业大学 2020-12-05
  • MutPrimerDesign:用于人类基因编码区域突变位点的引物设计程序
    MutPrimerDesign:用于人类基因编码区域突变位点的引物设计程序曹英豪,彭公信(中国医学科学院基础医学研究所&北京协和医学院基础医学院,北京100730)摘要:位于基因编码区的DNA突变与基因的功能密切相关。在已知人类基因编码区的突变位点时,如何在基因组上设计引物验证该突变是一个重要的问题 ...
    本站小编 哈尔滨工业大学 2020-12-05
  • 肺鳞状细胞癌发生的早期标志物及肿瘤预测模型
    肺鳞状细胞癌发生的早期标志物及肿瘤预测模型尚文慧,王晓曦,李晓琴,高斌(北京工业大学生命科学与生物工程学院,北京100124)摘要:选取癌症基因组图谱数据库的肺鳞状细胞癌(LungSquamousCellCarcinoma,LUSC)样本作为数据集,在全基因组的水平上研究肺鳞状细胞癌病人从正常到发病 ...
    本站小编 哈尔滨工业大学 2020-12-05
  • 毛竹SWEET基因家族的全基因组鉴定与分析
    毛竹SWEET基因家族的全基因组鉴定与分析李新然1,2,张智俊1,2,喻珮瑶1,2,陈家璐1,2,马瑞芳1,2,刘笑雨1,2(1.浙江农林大学省部共建亚热带森林培育国家重点实验室,杭州311300;2.浙江农林大学林业与生物技术学院,杭州311300)摘要:糖外排转运蛋白(Sugarswillev ...
    本站小编 哈尔滨工业大学 2020-12-05