哮喘是一种慢性呼吸道疾病,主要表现为反复发作的咳嗽、喘息、胸闷和呼吸困难等症状,病理特征包括气道重塑和慢性炎症等。由于儿童的肺功能发育和免疫功能均尚不健全,因此儿童哮喘的发病率较高,并呈逐年上升的趋势
[1-3]。过敏性哮喘是儿童哮喘的最主要类型,占儿童哮喘的80%以上
[4]。尘螨、蟑螂和宠物毛发是儿童过敏性哮喘最常见的室内诱因
[5]。儿童过敏性哮喘的发病机制十分复杂,目前认为主要与遗传、环境、免疫等因素有关
[6]。因此,积极探索儿童过敏性哮喘的发病机制,寻找有效的治疗方法和预防策略,具有重要意义。
由于高通量RNA测序的广泛应用和生物信息学的发展,在真核生物中越来越多的环状RNA(circular RNA, circRNA)已被鉴定出来
[7]。circRNA是包含外显子和/或内含子序列的稳定闭合RNA结构
[8]。研究发现,circRNA能够作为内源竞争RNA(competitive endogenous RNA, ceRNA),竞争性结合微小RNA(microRNA, miRNA),间接调控下游靶基因的功能,从而参与多种生物学过程
[9-10]。近年来已有研究表明,ceRNA在哮喘的气道重塑、气道炎症等方面发挥了重要作用。Lin等
[11]发现circHIPK3能够作为miR-326的海绵,通过靶向STIM1调节气道平滑肌细胞的增殖、迁移及凋亡,从而参与哮喘的气道重塑。Huang等
[12]研究发现,has_circ_0005519可能通过调节CD4
+T细胞中的hsa-let-7a-5p从而促进白细胞介素(interleukin, IL)-13和IL-6的表达,进而参与哮喘气道炎症。而关于ceRNA机制在儿童过敏性哮喘中的研究仍处于初级阶段,有待进一步探索。
基于目前的研究基础和现状,本研究通过对基因表达数据库和生物信息学方法的探索,挖掘儿童过敏性哮喘患者与正常儿童血液中显著差异表达的circRNA、miRNA和信使RNA(messenger RNA, mRNA),构建一个完整的circRNA-miRNA-mRNA调控网络,为深入研究circRNA在儿童过敏性哮喘中的作用及其潜在机制并优化治疗方案提供依据。
1 资料与方法
1.1 数据来源
通过GEO数据库(
https://www.ncbi.nlm.nih.gov/geo/),以“asthma”“children”为关键词检索儿童哮喘相关的基因芯片数据集,对相关的数据集进一步筛选,筛选标准:(1)样本为过敏性哮喘患儿组和健康儿童对照组;(2)年龄6~12岁,根据筛选标准下载获得哮喘患儿的3组相关数据集:circRNA数据集(GSE143192)、miRNA数据集(GSE120172)和mRNA数据集(GSE40732)。
1.2 筛选准则
经处理确定数据可用于分析后,使用R软件中的“limma”数据包完成差异分析,鉴定过敏性哮喘患儿和健康儿童间的差异表达基因。在GSE143192中adj.p.value<0.05以及|log2FC|≥1.5设置为阈值,鉴定DEcircRNA。在GSE120172和GSE40732中,将adj.p.value<0.05以及|log2FC|≥0.25设置为阈值,分别鉴定DEmiRNA和DEmRNA。随后,使用“ggplot2”和“ComplexHeatmap”数据包将结果可视化为火山图和聚类图。
1.3 构建circRNA-miRNA-mRNA调控网络
在circRNA、miRNA和mRNA数据集中提取出DEcircRNA、DEmiRNA和DEmRNA。将DEcircRNA通过ENCORI和circBank预测出相应的靶标miRNA,结合DEmiRNA进一步筛选目标miRNA。随后,通过miRDB、miRTarBase和TargetScan预测出相应的mRNA,再与DEmRNA通过venn图分析以识别目标mRNA。根据ceRNA理论,通过Cytoscape 3.9.0软件构建circRNA-miRNA-mRNA调控网络。
1.4 靶基因功能和通路富集分析
通过使用DAVID(
https://david.ncifcrf.gov/)对靶基因进行基因本体论(Gene Ontology, GO)的功能注释富集分析和京都基因和基因组百科全书(Kyoto Encyclopedia of Genes and Genomes, KEGG)通路富集分析,以此明确基因的功能,并筛选参与的相关生物学通路。
1.5 免疫细胞浸润及其与关键mRNA相关性的分析
利用CIBERSORT对样本中的22种免疫细胞浸润情况进行计算,筛选出P<0.05的样本进行后续分析。使用R软件及其“corrplot”“vioplot”数据包,进行Spearman相关性分析及绘制相关热图,并比较和可视化儿童过敏性哮喘患者与正常样本之间的22种免疫细胞浸润率。通过Pearson相关系数计算,分析关键mRNA与免疫细胞浸润之间的相关性,并借助R软件中的“ggpubr”数据包进行可视化。
1.6 研究对象
本研究选取就诊于复旦大学附属儿科医院的哮喘患儿20例,纳入标准参照《儿童支气管哮喘诊断与防治指南(2016年版)》
[13],诊断标准为复旦大学附属儿科医院呼吸科门诊确诊为哮喘并排除其他呼吸道疾病(如肺炎)及其他系统疾病的患儿。对照组选取健康体检的儿童20例,排除存在过敏病史及家族史的儿童。收集哮喘组和对照组儿童的年龄、性别和体重等相关资料,并按照性别与年龄进行匹配。抽取每位受试者1~2 mL血液标本,进行外周血单个核细胞(peripheral blood mononuclear cells, PBMC)的分离。本研究经复旦大学附属儿科医院伦理委员会批准(2020-81),研究对象监护人均签署知情同意书。
1.7 PBMC分离
将血液用PBS等体积稀释,缓慢加入等体积的Ficoll-Hypaque中,400 g缓升缓降离心30 min,吸取第2层单个核细胞悬浮液,用PBS洗2次,所获得的细胞即为PBMC。
1.8 RNA提取、逆转录和qPCR
采用Trizol法提取RNA,分别用Takara逆转录试剂盒和qPCR试剂盒进行RNA逆转录和实时定量PCR,引物序列见
表1。以β-actin为内参,数据分析采用2
-△△CT法。
1.9 统计学分析
数据采用R软件进行统计分析。采用CIBERSORT软件计算免疫细胞浸润百分率,Pearson相关系数分析关键mRNA与免疫细胞浸润的相关性。数据采用GraphPad Prism 8软件进行处理。正态分布的计量资料以均数±标准差()表示,组间比较采用两样本t检验。P<0.05代表差异有统计学意义。
2 结果
2.1 研究对象基本情况
共纳入哮喘患儿20例,其中男13例,女7例,年龄为(8.1±2.6)岁,体重指数(body mass index, BMI)为(17.5±2.9)kg/m2。对照组20例,其中男12例,女8例,年龄为(7.2±2.7)岁,BMI为(17.7±3.4)kg/m2。两组在性别、年龄和体重指数方面比较差异均无统计学意义(P>0.05)。
2.2 筛选儿童过敏性哮喘患者和对照者的DEcircRNA、DEmiRNA和DEmRNA
从GEO中下载儿童过敏性哮喘患者的circRNA、miRNA和mRNA表达谱数据。随后获取差异表达的circRNA、miRNA和mRNA。与对照样本相比,GSE143192共鉴定出130个DEcircRNA,包括81个上调和49个下调的circRNA;GSE120172共鉴定出40个DEmiRNA,包括16个上调和24个下调miRNA;GSE40732共鉴定出802个DEmRNA,包括602个上调和200个下调基因mRNA。为直观反映差异基因的表达情况,对差异基因进行了分层聚类分析,绘制了热图,图中颜色深浅代表基因表达的高低,红色代表上调,蓝色代表下调。绘制的DEM火山图中,红点表示上调的基因,蓝点表示下调的基因,灰点表示无显著差异。具体见
图1。
2.3 构建ceRNA调控网络
利用数据库预测DEcircRNA、DEmiRNA和DEmRNA之间的相互作用,构建circRNA-miRNA-mRNA的调控网络。利用ENCORI和circBank预测出287个与DEcircRNA相互作用的miRNA,并将预测得到的miRNA与DEmiRNA通过veen图分析,取交集获得14个miRNA,采用TargetScan和 miRDB数据库预测到4 090个DEmiRNA的靶基因,将预测得到的mRNA与DEmRNA取交集获得184个差异基因(
图2A~B)。将上述得到的DEcircRNA和DEmiRNA相互作用对以及DEmiRNA和DEmRNA相互作用对导入Cytoscape3.9.0软件中,最终获得由12个circRNA、7个miRNA和75个mRNA组成的94个基因ceRNA调控网络,选取其中的hsa_circ_0047447、MTARC1、RICTOR和GK5进行RT-qPCR验证,发现hsa_circ_0047447和MTARC1在哮喘组PBMC中的表达量高于对照组,而RICTOR和GK5在哮喘组PBMC中的表达量低于对照组(
图3)。
2.4 GO和KEGG富集分析
为了探索ceRNA网络中差异表达mRNA的潜在生物学功能,采用DAVID进行GO和KEGG分析。GO分析发现,生物学过程主要参与调节生长发育、调节细胞成分、与宿主相互作用、调节细胞极性、参与病毒的生命周期等;细胞组分主要包括:细胞迁移过程中形成的前沿区域、由细胞膜包围的细胞表面突出部分、轴突生长锥、细胞分裂等;分子功能主要为:结合胰岛素受体底物、结合外源蛋白、模式识别受体活性、转录辅抑制因子活性、核苷酸跨膜转运蛋白活性等(
图4A)。KEGG结果显示,Nod样受体信号通路、mTOR信号通路、胰岛素信号通路、肌醇磷酸代谢、内吞作用、hedgehog信号通路等(
图4B)。
2.5 分析免疫细胞浸润及其与关键mRNA的相关性
利用CIBERSORT对免疫细胞浸润水平进行分析,发现哮喘组的CD8
+ T细胞和静息NK细胞高于对照组,而对照组中的静息CD4
+记忆T细胞和活化肥大细胞高于哮喘组(
图5A)。免疫细胞之间的相关分析表明,在上调的免疫细胞中,CD8
+ T细胞与静息肥大细胞、活化CD4
+记忆T细胞呈较强的正相关性,而与静息CD4
+记忆T细胞、活化肥大细胞、浆细胞呈较强的负相关性。静息NK细胞与单核细胞、活化的树突状细胞、活化的CD4
+记忆T细胞呈较强的正相关性,与嗜酸性粒细胞、活化的NK细胞、静息树突状细胞呈较强的负相关性。下调免疫细胞中,静息CD4
+记忆T细胞与记忆B细胞、调节性T细胞、初始T细胞、M2巨噬细胞、CD8
+ T细胞呈较强的负相关性,活化肥大细胞与嗜酸性粒细胞、中性粒细胞呈较强的正相关性,与单核细胞、初始T细胞、CD8
+ T细胞呈较强的负相关性(
图5B)。通过对ceRNA调控网络中的75个mRNA与22种免疫细胞浸润的相关性分析,根据在上调免疫细胞CD8
+ T细胞、静息NK细胞中显著高表达,而在下调免疫细胞静息CD4
+记忆T细胞、活化肥大细胞中显著低表达,共筛选出KIAA0930、PPARGC1B、RAB5B、SNTB2、TRIM5和USP28这6个关键mRNA,进一步地,通过RT-qPCR检测发现KIAA0930、PPARGC1B和TRIM5在哮喘患儿PBMC中高表达,SNTB2低表达,而RAB5B和USP28在2组之间差异无统计学意义(
图6)。
3 讨论
哮喘是一种高异质性的常见慢性炎症性疾病
[14]。已有研究证实,儿童哮喘与过敏有关,并且受到遗传和环境因素的影响,这些因素决定了哮喘的易感性和严重程度
[15-16]。过敏性哮喘是儿童哮喘的最主要类型。近年来本课题组及其他学者的研究表明,circS100A11、LncRNA
PTPRE-
AS1等非编码RNA在儿童过敏性哮喘中发挥重要作用,可能成为治疗的靶点
[17-19]。因此本研究通过构建circRNA-miRNA-mRNA调控网络,进一步探讨非编码RNA在儿童过敏性哮喘中的作用,具有重要的意义。
本研究通过生物信息学分析,构建了哮喘患儿的ceRNA调控网络,并揭示了其在免疫浸润中的特征。我们从数据集中筛选出儿童过敏性哮喘和健康儿童之间的130个DEcircRNA、40个DEmiRNA和802个DEmRNA。基于RNA之间的相互作用,构建了包含12个circRNA、7个miRNA和75个mRNA的ceRNA网络。随后,通过RT-qPCR验证发现hsa_circ_0047447、MTARC1、RICTOR与预测结果一致,在哮喘患儿PBMC中的表达量高于正常儿童,而GK5与预测结果相反,在哮喘患儿PBMC中的表达量低于对照组,后续将通过扩大样本量进行进一步验证。我们对差异表达的mRNA进行了功能和通路富集分析,探索其潜在的生物学功能。GO结果表明,差异基因主要富集于细胞极性、细胞分裂和模式识别受体活性等。KEGG分析发现,ceRNA网络中的差异基因主要参与了mTOR信号通路、Nod样受体信号通路以及肌醇磷酸代谢等通路。在circRNA相关研究中,ceRNA机制被认为是circRNA发挥功能的最重要途径之一。Su等
[20]研究发现,哮喘患者血浆中高表达的circVPS33A,能够通过吸附miR-192-5p调节下游HMGB1的表达,从而发挥作用。Wang等
[21]发现哮喘儿童肺泡灌洗液中高表达的circZNF652吸附miR-452-5p调节下游JAK2/STAT6信号通路,从而促进杯状细胞的分化。circERBB2通过调控miR-98-5p/IGF1R轴参与哮喘的气道重塑
[22]。在小鼠哮喘模型中发现,mmu_circ_0001359可吸附miR-183-5p,促进FoxO1信号介导的M2型巨噬细胞活化;利用体内、外实验发现人工制备包含mmu_circ_0001359的外泌体可以降低炎症水平并减缓气道重塑
[23]。该研究提示我们,可以利用外泌体包裹circRNA,以ceRNA调控网络为靶点开发儿童过敏性哮喘的新型疗法。
既往研究发现,ceRNA网络在免疫浸润、免疫逃逸、免疫治疗等调节免疫反应中具有重要作用
[24-26]。在儿童过敏性哮喘中,免疫系统也扮演着重要的角色,因此,探讨ceRNA网络与免疫细胞之间的潜在关系具有重要的意义。免疫浸润分析显示,儿童过敏性哮喘患者CD8
+ T细胞和静息NK细胞高于对照组,而CD4
+记忆T细胞和活化肥大细胞低于对照组。ceRNA调控网络在T细胞中发挥重要作用,在哮喘患者CD4
+T细胞中circ_0005519表达上调,并可能与hsa-let-7a-5p结合,从而减轻对CD4
+ T细胞中IL-13/IL-6的抑制作用
[12]。circUBAP2介导的ceRNA系统通过调控CD8
+ T细胞的浸润调节胰腺癌
[27]。NK细胞是机体重要的免疫细胞,能够参与超敏反应和自身免疫性疾病的发生
[28]。Carpio-Pedroza等
[29]研究发现,哮喘患儿外周血NK细胞较健康对照组升高,并且分泌较多的IL-4和较少的γ干扰素,表明NK细胞能促进哮喘的炎症反应。肥大细胞是过敏性哮喘的主要效应细胞之一,在哮喘的炎症中具有重要作用。在严重哮喘患儿支气管活检样本中发现,黏膜下肥大细胞的数量增多
[30]。而本研究显示,哮喘患儿活化的肥大细胞浸润率较对照组低,这可能与哮喘患儿的不同发病阶段有关。
值得注意的是,通过对ceRNA网络中的75个mRNA与免疫细胞进行相关性分析,筛选出KIAA0930、PPARGC1B、RAB5B、SNTB2、TRIM5和USP28这6个关键mRNA。PPARGC1B是一种核激素受体,参与细胞激活、分化、增殖和凋亡,PPARGC1B基因中的遗传变异可能会影响PPAR的功能,进而可能导致患者的气道炎症和重塑
[31]。已有研究揭示,在大多数癌症中,USP28与NK T细胞浸润呈显著正相关,在非小细胞肺癌中,USP28在CD8
+ T细胞中高表达
[32]。Jimenez-Moyano等
[33]发现了非人类TRIM5变异在细胞免疫中的功能,TRIM5可以将先天病毒感应和CD8
+ T细胞激活结合起来,以削弱物种对逆转录病毒的感染能力。上述研究表明,TRIM5和USP28在一些疾病的细胞免疫中扮演了重要角色。本研究发现KIAA0930、PPARGC1B和TRIM5在哮喘患儿PBMC中显著高表达,SNTB2显著低表达,而RAB5B和USP28在两组之间没有显著性差异,后续我们将通过扩大样本量进一步验证,为儿童过敏性哮喘的治疗提供更加有力的证据。
综上所述,本研究通过生物信息学分析,筛选出儿童过敏性哮喘患者中差异表达的circRNA、miRNA和mRNA,构建了circRNA-miRNA-mRNA免疫相关的调控网络,并分析了ceRNA网络中差异表达基因的生物学功能和调控通路,为进一步探究儿童过敏性哮喘的ceRNA机制和免疫治疗提供更多的可能性。然而,本研究得出的结论需要通过体内、体外实验进一步验证,实验完善后有望为治疗儿童过敏性哮喘提供更可靠的依据。
上海市卫生健康委员会科研项目(20214Y0440)
国家自然科学基金青年项目(82100033)
国家自然科学基金面上项目(82470029)