中华皮肤科杂志 ›› 2019, Vol. 52 ›› Issue (10): 729-735.doi: 10.35541/cjd.20190001

• 论著 • 上一篇    下一篇

基于高通量RNA转录组测序技术分析特应性皮炎患者皮损与非皮损组织差异表达基因

陈丽洁    梁景耀    张锡宝    邵蕾    潘清丽    何素玲    刘玉梅    王建琴   

  1. 广州医科大学皮肤病研究所  广州市皮肤病防治所  510095
  • 收稿日期:2018-11-02 修回日期:2019-07-03 发布日期:2019-09-30
  • 通讯作者: 王建琴;刘玉梅 E-mail:jianqinwang@foxmail.com; 1109180611@ qq.com

Identification of differentially expressed genes in lesional versus nonlesional skin of patients with atopic dermatitis by using high-throughput transcriptome-wide RNA sequencing

Chen Lijie, Liang Jingyao, Zhang Xibao, Shao Lei, Pan Qingli, He Suling, Liu Yumei, Wang Jianqin   

  1. Institute of Dermatology, Guangzhou Medical University, Department of Dermatology, Guangzhou Institute of Dermatology, Guangzhou 510095, China
  • Received:2018-11-02 Revised:2019-07-03 Published:2019-09-30
  • Contact: Wang Jianqin; Liu Yumei E-mail:jianqinwang@foxmail.com; 1109180611@ qq.com

摘要: 【摘要】 目的 探讨中重度特应性皮炎(AD)患者皮损和非皮损组织转录组差异表达基因,并阐述差异表达基因在AD发病机制中的作用。方法 2016年7 - 10月于广州市皮肤病防治所/广州医科大学皮肤病研究所门诊收集5例汉族AD患者的皮损和非皮损组织,运用二代高通量RNA转录组测序(RNA-seq)技术进行测序,筛选差异表达基因,并对其进行GO功能注释及 KEGG通路分析。实时荧光定量PCR(qRT-PCR)验证候选基因在皮损和非皮损组织中的表达差异。结果 10个样品测序平均获得10.96 Gb数据,共检测到21 729条基因,其中19 268条为已知基因,2 545条为预测的新基因;共检测出23 153个新转录本,其中18 889个属于已知蛋白编码基因新的可变剪接亚型,2 545 个属于新的蛋白编码基因的转录本,其他1 719个属于长链非编码RNA。78条基因在皮损和非皮损中表达差异有统计学意义,其中67条在皮损中高表达,11条低表达,包括已知与AD炎症(CXCL1/2/8、IL6/IL1β、MMP1、SERPINB4、S100A2、GZMB、OASL、OSM)、屏障(KRT16、FABP5、CYP1A1)、角质形成细胞分化(IL-20)等相关的基因。GO分析显示,72条差异基因获得功能注释。KEGG 通路分析表明,差异表达基因共富集到132条信号通路中,其中13条通路显著富集,包括白细胞介素17通路、NOD样受体信号通路、Toll样受体信号通路等。qRT-PCR 结果显示,候选基因CXCL1、KRT6A、IL36A、SERPINB4和 PSAPL1的mRNA表达与转录组测序结果趋势一致。结论 本研究在转录组水平上筛选出AD患者皮损和非皮损组织差异表达基因及相关重要调节信号通路,发现汉族AD患者白细胞介素17通路最为富集,为深入探讨AD发病机制提供重要的基础。

关键词: 皮炎, 特应性; 转录组; 序列比对; 白细胞介素17; RNA-seq

Abstract: 【Abstract】 Objective To identify differentially expressed genes in the transcriptome of the lesional versus nonlesional skin tissues of patients with moderate and severe atopic dermatitis(AD), and to elucidate their roles in the pathogenesis of AD. Methods From July to October in 2016, lesional and nonlesional skin tissues were obtained from 5 outpatients of Han nationality with AD in Guangzhou Institute of Dermatology, Institute of Dermatology, Guangzhou Medical University. The next-generation high-throughput transcriptome-wide RNA sequencing (RNA-seq) was performed to identify differentially expressed genes, which were subjected to GO function annotation and KEGG pathway analysis. Real-time fluorescence-based quantitative PCR (qRT-PCR) was conducted to verify differences in candidate gene expression between lesional and nonlesional skin tissues. Results An average of 10.96 GBs sequence reads were acquired among 10 samples. A total of 21 729 genes were detected, including 19 268 known genes and 2 545 predicted novel genes. A total of 23 153 new transcripts were detected, of which 18 889 were new alternative splicing subtypes of known protein-coding genes, 2 545 were transcripts belonging to new protein-coding genes, and the remaining 1 719 belonged to long-stranded non-coding RNA. Totally, 78 differentially expressed genes were identified between the lesional and nonlesional skin tissues, including 69 upregulated and 11 downregulated genes in the lesional skin tissues. Among them, there were several genes known to be associated with AD inflammation (CXCL1/2/8, IL6/IL1β, MMP1, SERPINB4, S100A2, GZMB, OASL, OSM) and barrier (KRT16, FABP5, CYP1A1) and keratinocyte differentiation (IL-20). GO analysis revealed that functions of 72 differentially expressed genes could be annotated. KEGG pathway analysis showed that the differentially expressed genes were grouped into 132 signaling pathways, of which 13 were significantly enriched, including the interleukin (IL)-17 pathway, NOD-like receptor signaling pathway, Toll-like receptor signaling pathway, etc. qRT-PCR showed that the mRNA expression levels of candidate genes CXCL1, KRT6A, IL36A, SERPINB4 and PSAPL1 was consistent with the transcriptome sequencing results. Conclusions Differentially expressed genes and related important regulatory signaling pathways were identified between the lesional and nonlesional skin tissues of patients with AD at the transcriptional level, and the IL-17 pathway was found to be mostly enriched in AD lesions in patients of Han nationality. These findings provide an important basis for further study on the pathogenesis of AD.

Key words: Dermatitis, atopic, Transcriptome, Sequence alignment, Interleukin-17, RNA-seq

引用本文

陈丽洁 梁景耀 张锡宝 邵蕾 潘清丽 何素玲 刘玉梅 王建琴. 基于高通量RNA转录组测序技术分析特应性皮炎患者皮损与非皮损组织差异表达基因[J]. 中华皮肤科杂志, 2019,52(10):729-735. doi:10.35541/cjd.20190001

Chen Lijie, Liang Jingyao, Zhang Xibao, Shao Lei, Pan Qingli, He Suling, Liu Yumei, Wang Jianqin. Identification of differentially expressed genes in lesional versus nonlesional skin of patients with atopic dermatitis by using high-throughput transcriptome-wide RNA sequencing[J]. Chinese Journal of Dermatology, 2019, 52(10): 729-735.doi:10.35541/cjd.20190001