Changing expression profiles of mRNA,lncRNA, circRNA,and miRNA in lung tissue reveal the pathophysiological of bronchopulmonary dysplasia (BPD) in mouse model

New perinatal care technologies have improved the survival rate of preterm neonates, but the prevalence of bronchopulmonary dysplasia (BPD), one of the most intractable problems in neonatal intensive care unit (NICU), remains unchanged. In present study, high-throughput sequencing (HTS) was performed to detect the expression profiles of long noncoding RNAs (lncRNAs), messenger RNAs (mRNAs), circular RNAs (circRNAs), and microRNAs (miRNAs) in hyperoxia-induced BPD mouse model. Significant differentially expressed RNAs were selected and clustered between the BPD group and the control group. The results revealed that expressions of 1778 lncRNAs, 1240 mRNAs, 97 circRNAs, and 201 miRNAs were significantly altered in the BPD group. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) were performed to predict the potential functions of differentially expressed RNAs. lncRNA-mRNA and circRNA-miRNA coexpression networks were constructed to detect their association with the pathogenesis of BPD. Our study provides a systematic perspective on the potential function of RNAs during BPD.     

小鼠正常黑素细胞和小鼠B16黑素瘤细胞的环状RNAs表达谱

黑素瘤是一种多发于皮肤的恶性肿瘤,因其侵袭性强,预后差等特点一直是科研人员关注的热点。环状RNAs(circRNAs)是一种新型内源性非编码RNA,广泛参与动物生长发育、细胞分化和信号转导等生理过程,但circRNAs在黑素瘤细胞内的分子机制尚未被充分解析。本研究以小鼠(C57BL/6J)正常黑素细胞及B16黑素瘤细胞为研究对象,采用二代测序技术分析两种细胞间circRNAs表达特性。测序结果显示,小鼠正常黑素细胞和黑素瘤细胞中共有851个circRNAs,其中195个差异表达circRNAs(DECs)。GO及KEGG数据库注释发现,DECs的来源基因主要参与细胞周期(cell cycle)、紧密结合(tight junction)、Rap1信号通路(Rap1 signaling pathway)、TGF-beta信号通路(TGF-beta signaling pathway)等与细胞增殖、迁移相关的信号通路;探究发现,黑素瘤细胞中显著性高表达的circE2F5(circ-3:14578602|14606309)通过上调E2F5的表达促进黑素瘤细胞增殖。circRNA靶基因预测发现,...     

Identification and characterization of circular RNA in lactating mammary glands from two breeds of sheep with different milk production profiles using RNA-Seq

CircRNA is a specific type of non-coding RNA that has been shown to have an important role in mammary gland (MG) activity, but no study of MG circRNA activity in sheep so far. In this study, the expression profile of circRNAs was investigated using RNA-Seq in MG parenchyma at peak lactation from Small-Tailed Han sheep and Gansu Alpine Merino sheep with phenotypic differences in milk yield and components. A total of 4, 906 circRNAs were found and 33 of these were differentially expressed between breeds. GO and KEGG results showed that the parental genes of differentially expressed circRNAs were mainly enriched in heterocyclic compound binding, kinase activity, adherens junction, the TGF-β signaling pathway, and the MAPK signaling pathway. This study provides an overview of circRNA expression in the ovine MG and the interaction between some key circRNAs and their target miRNAs. It improves our knowledge of the role of circRNA in sheep milk synthesis.     

结节性甲状腺肿circRNA-miRNA-mRNA调控网络的构建

基于生物信息分析筛选结节性甲状腺肿中差异表达的环状RNA(circRNA),并揭示circRNA-miRNA-mRNA调控网络在结节性甲状腺肿中的作用。从GEO数据库中检索结节性甲状腺肿组织基因芯片数据,利用R软件筛选出差异表达的circRNA。联合多个生物信息数据库预测差异表达circRNA下游的miRNA及mRNA, 并对靶mRNA进行GO及KEGG富集分析。利用STRING在线数据库及Cytoscape软件筛选核心基因。确定了2个circRNA,42个miRNA及546个mRNA。GO及KEGG富集分析表明靶mRNA主要涉及细胞生长及基因表达调控过程。基于Cytoscape软件筛选出了14个核心基因(SP1、IGF1R、RPS6KB1、SMAD2、SMAD3、SMAD4、VEGFA、CCND1、CDK2、HSPA4、HIF1A、CREB1,NR3C1和STAT5A)。最终基于2个circRNA、11个miRNA和14个核心mRNA构建了circRNA-miRNA-mRNA调控网络。结节性甲状腺肿组织中异常表达的circRNA及相关的circRNA-miRNA-mRNA调控网络可能成为结节性甲状腺肿诊断与治疗的新靶点。     

Profiling and integrated analysis of differentially expressed circRNAs as novel biomarkers for breast cancer

Breast cancer (BC) is a globally common cancer with the highest and increasing morbidity and mortality among females. Novel biomarkers are warranted to be discovered for the early detection, treatment, and prognosis of BC. In this study, we investigated the profiles of differentially expressed (DE) circular RNAs (circRNAs) by competing endogenous RNAs (ceRNA) microarray to construct a genome-wide circRNA profile. Then, we performed Gene Ontology (GO) analysis and Kyoto Encyclopedia of Gene and Genome (KEGG) pathway analysis of the host genes (HGs) of circRNAs. A total of 4,370 DE circRNAs were detected and GO and KEGG analysis showed that they were significantly associated with cell cycle, DNA replication, BC, and familial BC. We validated the differential circRNAs and relevant HGs through quantitative real-time polymerase chain reaction and constructed a putative circRNA–microRNA–messenger RNA regulatory network. Eight circRNAs, including hsa_circ_0069094, hsa_circ_0062558, hsa_circ_0074026, hsa_circ_0079876, hsa_circ_0017536, hsa_circ_0023302, hsa_circ_0017650, and hsa_circ_0017545, were validated significantly DE in BC tissue and associated with TNM staging, lymph node infiltration, and Ki67. Hsa_circ_0069094, hsa_circ_0079876, hsa_circ_0017650, and hsa_circ_0017526 were upregulated in plasma. This study revealed the general expression characteristics of specific DE circRNAs in BC and hsa_circ_0069094, hsa_circ_0079876, hsa_circ_0017650, and hsa_circ_0017526 might be promising candidate targets.     

A novel strategy of identifying circRNA biomarkers in cardiovascular disease by meta-analysis

Circular RNAs (circRNAs) are stable and abundantly expressed in vivo but are abnormally expressed in several diseases. This study aimed to identify circRNAs acting as potential biomarkers for cardiovascular disease (CVD). Research were retrieved from the articles published by September 2018 in eight databases to compare circRNA expression profiles between CVD and non-CVD in human and animal models. Meta-analysis under a random effects model was conducted. Subgroup analysis of tissue, species, and disease-specific circRNAs was examined. Sensitivity analysis was performed to explain the uncertainty among all studies. Diagnostic accuracy of circRNAs in CVD was analyzed to testify the discriminative ability. Bioinformatics analysis including Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis was conducted. Among 6,284 differentially expressed circRNAs from 32 original studies, only 322 circRNAs were reported in three or more studies. The meta-analysis identified 63 significantly dysregulated circRNAs, 44 upregulated and 19 downregulated. Among the tissue-specific or disease-specific circRNAs identified in the subgroup analysis, two circRNAs (circCDKN2BAS and circMACF1) showed the potential to be circulating biomarkers for CVD. Sensitivity analysis demonstrated 69% of circRNAs were in conformity with the overall analysis. The pooled diagnostic odds ratio was 2.94 (95% confidence interval [CI], 2.35–3.58), and the overall area under the curve value was 0.86 (95% CI, 0.83–0.89). GO and KEGG enrichment analyses indicated that the target genes of circRNAs participate in cardiogenesis-related processes and pathways. This study demonstrates circRNAs have a high diagnostic value as potential biomarkers for CVD, and two candidate circRNAs, circCDKN2BAS and circMACF1, are potential circulating biomarkers for CVD diagnosis and treatment.     

Identification of downregulated circRNAs from tissue and plasma of patients with gastric cancer and construction of a circRNA-miRNA-mRNA network

The connection between circular RNAs (circRNAs) and gastric cancer has been reported widely in recent years. However, previous studies have focused mainly on circRNAs from gastric cancer tissue. The objectives of the present study were to detect dysregulated circRNAs from both tissue and plasma of patients with gastric cancer and to explore their potential roles in the pathogenesis of gastric cancer. Expression profiles of circRNAs were obtained from the Gene Expression Omnibus (GEO) and analyzed using the GEO2R tool to identify differential expressed circRNAs. The significance threshold was set as |log2 (fold change)| > 2 and adjusted P < .05. The microRNA (miRNA) binding sites of the differentially expressed circRNAs were predicted using the Circular RNA Interactome web tool. TargetScan and the miRNet database were used to predict the miRNA target genes. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathway analyses were performed using Database for Annotation Visualization and Integrated Discovery. Hub genes were identified and a network was constructed with Cytoscape. The overall survival rates for the selected miRNAs and messenger RNAs were evaluated by Kaplan-Meier Plotter. A total of three downregulated circRNAs (hsa_circ_0001190, hsa_circ_0036287, and hsa_circ_0048607) were identified in this study. Six miRNAs and eight hub genes met the significance criteria and were selected for further analysis. A circRNA-miRNA-hub gene network was constructed based on three circRNAs, six miRNAs, and eight hub genes. Evaluation of overall survival rates for the hub genes showed that low expression levels of GADD45A, PPP1CB, PJA2, and KLF2 were associated with poor overall survival. This study identified potential novel plasma circRNA biomarkers and provides insights into the underlying mechanisms of gastric cancer pathogenesis.     

Characteristics of circular RNA expression of pulmonary macrophages in mice with sepsis‐induced acute lung injury

Circular RNAs (circRNAs) make up a large class of non‐coding RNAs and play important roles in the pathology of a variety of diseases. However, their roles in pulmonary macrophage polarization after sepsisinduced lung injury is unknown. In this study, mice were divided into two groups: Sham control group and cecal ligation and puncture (CLP)‐induced ALI group. Macrophages were isolated from lung homogenates 24 hours after SCLP/CLP. We started with RNA‐seq of circRNA changes in macrophages and validated by RT‐PCR in the following experiments. A total of 4318 circRNAs were detected in the two groups. Of these, 11 and 126 circRNAs were found to be significantly upregulated and downregulated, respectively, compared to the control (p≤0.05, Fold Change ≥2). Differentially expressed circRNAs with a high foldchange (fold‐change >4, P<0.05) were selected for validation by qRT‐PCR, 10 of which were verified. Furthermore, the most differentially expressed circRNAs within all the comparisons were annotated in detail with circRNA/miRNA interaction information using miRNA target prediction software. The network of circRNA‐miRNA‐mRNA was illustrated by cytoscape software. Gene ontology analyses indicated the upregulated circRNAs were involved in the multiple biological functions such as regulation of mitochondrion distribution and Notch binding, while the down‐regulated circRNAs mainly involved in the biological process as histone H3K27 methylation. KEGG pathway analysis revealed TGF‐beta signaling pathway was related to the upregulated circRNAs. The present study provides a novel insight into the roles of circRNAs in pulmonary macrophage differentiation and polarization post septic lung injury.     

Identification of circRNA and mRNA expression profiles and functional networks of vascular tissue in lipopolysaccharide‐induced sepsis

Sepsis is the most common cause of death in intensive care units. This study investigated the circular RNA (circRNA) and mRNA expression profiles and functional networks of the aortic tissue in sepsis. We established a lipopolysaccharide (LPS)‐induced rat sepsis model. High‐throughput sequencing was performed on the aorta tissue to identify differentially expressed (DE) circRNAs and mRNAs, which were validated by real‐time quantitative polymerase chain reaction (RT‐qPCR). Bioinformatic analysis was carried out and coding and non‐coding co‐expression (CNC) and competing endogenous RNA (ceRNA) regulatory networks were constructed to investigate the mechanisms. In total, 373 up‐regulated and 428 down‐regulated circRNAs and 2063 up‐regulated and 2903 down‐regulated mRNAs were identified. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses of mRNAs showed that the down‐regulated genes were mainly enriched in the process of energy generation. CNC and ceRNA regulatory networks were constructed with seven DE circRNAs. The results of functional enrichment analysis of CNC target genes revealed the important role of circRNAs in inflammatory response. The ceRNA network also highlighted the significant enrichment in calcium signalling pathway. Significant alterations in circRNAs and mRNAs were observed in the aortic tissue of septic rats. In addition, CNC and ceRNA networks were established.     

Integrative analysis of circRNAs acting as ceRNAs involved in ethylene pathway in tomato

Circular RNAs (circRNAs) are a large class of non‐coding endogenous RNAs that could act as competing endogenous RNAs (ceRNAs) to terminate the mRNA targets' suppression of miRNAs. To elucidate the intricate regulatory roles of circRNAs in the ethylene pathway in tomato fruit, deep sequencing and bioinformatics methods were performed. After strict screening, a total of 318 circRNAs were identified. Among these circRNAs, 282 were significantly differentially expressed among wild‐type and sense‐/antisense‐LeERF1 transgenic tomato fruits. Besides, 1254 target genes were identified and a large amount of them were found to be involved in ethylene pathway. In addition, a sophisticated regulatory model consisting of circRNAs, target genes and ethylene was set up. Importantly, 61 circRNAs were found to be potential ceRNAs to combine with miRNAs and some of the miRNAs had been revealed to participate in the ethylene signaling pathway. This research further raised the possibility that the ethylene pathway in tomato fruit may be under the regulation of various circRNAs and provided a new perspective of the roles of circRNAs.     

意大利蜜蜂工蜂中肠的环状RNA及其调控网络分析

【目的】环状RNA(circRNA)在可变剪接、转录调控和来源基因的表达调控等方面具有重要功能。本研究旨在分析意大利蜜蜂Apis mellifera ligustica工蜂中肠circRNA的数量、种类、结构特征和作用,并通过构建和分析circRNA的调控网络探索circRNA的调控功能。【方法】在实验室条件下人工饲养意大利蜜蜂工蜂,利用circRNA-seq技术对意大利蜜蜂7和10日龄成年工蜂中肠样品进行深度测序。利用find_circ软件从质控后的数据中预测circRNA。通过BLAST比对GO和KEGG数据库,对circRNA的来源基因进行功能和代谢通路注释。利用TargetFinder软件预测circRNA靶向结合的miRNA及miRNA靶向结合的mRNA,通过Cytoscape v.3.2.1软件对circRNA-miRNA和circRNA-miRNA-mRNA调控网络进行构建及可视化。通过设计背靠背引物和线性扩增引物RT-PCR对预测出的circRNA进行验证。【结果】意大利蜜蜂工蜂中肠样品的测序平均得到136 463 071条clean reads,去除rRNA后各样品的anchor reads均在136 779 122条及以上。共预测出10 833个circRNA,长度主要介于15~1 000 nt;上述circRNA的类型丰富,其中已注释的外显子circRNA数量最多,分布在西方蜜蜂1号染色体的circRNA数量最多,其次为8号染色体。CircRNA的来源基因可注释到包括结合、细胞进程和细胞在内的45个GO条目,以及包括内吞作用、内质网蛋白加工及核糖体在内的121条KEGG代谢通路,表明circRNA在意大利蜜蜂工蜂中肠的生长、发育、新陈代谢和细胞生命活动等生物学过程中发挥重要作用。进一步构建circRNA-miRNA和circRNA-miRNA-mRNA调控网络,分析结果显示部分circRNA可能作为竞争性内源RNA吸附结合microRNA,从而调控基因的表达水平。最后,对随机选择的3个circRNA的RT-PCR结果验证了其真实存在。【结论】本研究对意大利蜜蜂工蜂中肠中的circRNA进行预测、分析及鉴定。研究结果提供了中肠circRNA的数量、种类、结构特征、作用和调控网络的信息,揭示了circRNA可能通过作用于来源基因和作为竞争性内源RNA在意大利蜜蜂工蜂中肠的生长发育和免疫防御中发挥作用,为深入研究circRNA在意大利蜜蜂中肠发育及胁迫响应过程中的功能奠定了基础。