意大利蜜蜂工蜂中肠的环状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在意大利蜜蜂中肠发育及胁迫响应过程中的功能奠定了基础。     

竞争性内源RNA的生物学功能及其调控

最新研究表明,RNA之间可以通过竞争结合共同的microRNA反应元件(microRNA response element, MRE)实现相互调节,这种调控模式构成竞争性内源RNA(Competing endogenous RNA, ceRNA)。已发现的ceRNA包括蛋白编码mRNA和非编码RNA,其中后者包括假基因转录物、长链非编码RNA(Long non-coding RNA, lncRNA)、环状RNA(Circular RNA, circRNA)等。文章主要从ceRNA分类的角度,阐述各类ceRNA构成的调控网络发挥的生物学功能在病理和生理相关过程中的作用,以及可能影响ceRNA调控有效性的因素。     

Circular RNA-associated ceRNA network involved in HIF-1 signalling in triple-negative breast cancer: circ_0047303 as a potential key regulator

The aggressive and highly metastatic nature of triple-negative breast cancer (TNBC) causes patients to suffer from the poor outcome. HIF-1 signalling pathway is a prominent pathway that contributes to angiogenesis and metastasis progression in tumours. On the contrary, the undeniable importance of circular RNAs (circRNAs) as multifunctional non-coding RNAs (ncRNAs) has been identified in breast cancer. These ncRNAs owing to their high stability and specificity have been becoming a hotspot in cancer researches. circRNAs act as competing endogenous RNAs (ceRNAs) and compete with mRNAs for shared miRNAs, thus modulate gene expression. Since the most dysregulated biological functions in TNBC are associated with cellular invasion, understanding the molecular pathogenesis of these processes is a crucial step towards the development of new treatment approaches. The purpose of this study is to undermine the circRNA-associated ceRNA network involved in HIF-1 signalling in TNBC using an integrative bioinformatics approach. In the next step, the novel circ_0047303-mediated ceRNA regulatory axes have been extracted and validated across TNBC samples. We show that circ_0047303 has the highest degree in the circRNA-associated ceRNA network and shows a significant up-expression in TNBC. Moreover, our results suggest that circ_0047303 could mediate the upregulation of key angiogenesis-related genes, including HIF-1, EIF4E2 and VEGFA in TNBC through sponging the tumour-suppressive miRNAs. The circ_0047303 could be a promising molecular biomarker and/or therapeutic target for TNBC.     

SRNAome and degradome sequencing analysis reveals specific regulation of sRNA in response to chilling injury in tomato fruit

Plant genomes encode diverse small RNA classes that function in distinct gene‐silencing pathways. To elucidate the intricate regulation of microRNAs (miRNAs) and endogenous small‐interfering RNAs (siRNAs) in response to chilling injury in tomato fruit, the deep sequencing and bioinformatic methods were combined to decipher the small RNAs landscape in the control and chilling‐injured groups. Except for the known miRNAs and ta‐siRNAs, 85 novel miRNAs and 5 ta‐siRNAs members belonging to 3 TAS families (TAS5, TAS9 and TAS10) were identified, 34 putative phased small RNAs and 740 cis/trans‐natural antisense small‐interfering RNAs (nat‐siRNAs) were also found in our results which enriched the tomato small RNAs repository. A large number of genes targeted by those miRNAs and siRNAs were predicted to be involved in the chilling injury responsive process and five of them were verified via degradome sequencing. Based on the above results, a regulatory model that comprehensively reveals the relationships between the small RNAs and their targets was set up. This work provides a foundation for further study of the regulation of miRNAs and siRNAs in the plant in response to chilling injury.     

Comprehensive analysis of circRNA expression profiles and circRNA-associated competing endogenous RNA networks in the development of mouse thymus

The thymus plays an irreplaceable role as a primary lymphoid organ. However, the complicate processes of its development and involution are incompletely understood. Accumulating evidence indicates that non-coding RNAs play key roles in the regulation of biological development. At present, the studies of the circRNA profiles and of circRNA-associated competing endogenous RNAs (ceRNAs) in the thymus are still scarce. Here, deep-RNA sequencing was used to study the biological mechanisms underlying the development process (from 2-week-old to 6-week-old) and the recession process (from 6-week-old to 3-month-old) of the mouse thymus. It was found that 196 circRNAs, 233 miRNAs and 3807 mRNAs were significantly dysregulated. The circRNA-associated ceRNA networks were constructed in the mouse thymus, which were mainly involved in early embryonic development and the proliferation and division of T cells. Taken together, these results elucidated the regulatory roles of ceRNAs in the development and involution processes of the mouse thymus.     

circRNA-associated ceRNA network construction reveals the circRNAs involved in the progression and prognosis of breast cancer

Recently, increasing evidences show that circular RNAs (circRNAs) are important regulators of various diseases, especially cancer. However, the regulatory role and the potential mechanism of action of circRNAs in breast cancer remain largely unknown. In this study, weighted gene co-expression network analysis was conducted with the differentially expressed miRNAs and mRNAs in breast cancer from The Cancer Genome Atlas database to identify the key modules associated with the carcinogenesis of breast cancer. In the significant turquoise and brown modules, 22 miRNAs and 1877 mRNAs were identified, respectively. Then, We compared and predicted the target genes and performed survival analysis to identify the miRNAs and mRNAs related to the prognosis of breast cancer. A circRNA-related competitive endogenous RNA network was identified by database co-screening, and deleted in liver cancer 1 (DLC1) was identified as a key gene. Finally, to assess how genes in key modules and key genes contribute to the development of breast cancer, relevant pathway information was obtained through DAVID and Gene Set Enrichment Analysis. These data demonstrated that three circRNAs (hsa-circ-0083373, hsa-circ-0083374, and hsa-circ-0083375) that regulate DLC1 expression via hsa-mir-511 and are involved in the pathogenesis and development of breast cancer.     

Characteristics of the competition among RNAs for the binding of shared miRNAs

Competing endogenous RNAs (ceRNAs) are RNAs that share common miRNA binding sites and compete with each other for the miRNA association at these sites. The observation of this phenomenon in the cells altered the view of the miRNA target RNAs from molecules that are passively controlled by miRNAs to molecules that also modulate the miRNAs activity. In this review, we build a general profile of ceRNAS characteristics in order to facilitate ceRNAs identification by researchers. The information summarized here contains an actualized list of previously reported ceRNAs and classes of RNAs that can participate in this type of interaction, the expression behavior and characteristics of ceRNAs and miRNAs in the context of competition, the influence of the shared MREs/miRNAs numbers and the miRNA binding strength on the competition, reports on competition between RNAs in different subcellular localizations and the concept that ceRNAs may form a huge regulatory network in the cell.     

Long noncoding RNA MEG3 play an important role in osteosarcoma development through sponging microRNAs

More and more evidence indicate long noncoding RNAs (lncRNAs) as competing endogenous RNAs (ceRNAs) to indirectly regulate messenger RNAs (mRNAs) by acting as microRNA (miRNA) sponges, which represents a novel layer of gene regulation that plays a critical role in the development of cancers. However, functional roles and regulatory mechanisms of lncRNA-mediated ceRNAs network in osteosarcoma are still largely unknown. Here, we comprehensively compared the expression profiles of mRNAs, lncRNAs, and miRNAs between osteosarcoma and normal samples from the Gene Expression Omnibus (GEO) to elaborate related latent mechanisms. Two lncRNAs, ie, LINC01560 and MEG3, were identified to be aberrantly expressed. Importantly, MEG3 was considered as a promising diagnostic biomarker and therapeutic target for patients with osteosarcoma according to the Kaplan-Meier analysis of another independent osteosarcoma data set from the Cancer Genome Atlas (P = 0.05). Eventually, we successfully established a dysregulated lncRNA-related ceRNA network, including one osteosarcoma-specific lncRNA, three miRNAs and four mRNAs. In conclusion, this study should be beneficial for improving our understanding of the lncRNA-mediated ceRNA regulatory mechanisms in the pathogenesis of osteosarcoma and providing it with novel candidate diagnostic and therapeutic biomarkers.     

Genome-wide identification of a competing endogenous RNA network in cholangiocarcinoma

Cholangiocarcinoma (CCA) is the second widespread liver tumor with relatively poor survival. Increasing evidence in recent studies showed long noncoding RNAs (lncRNAs) exert a crucial impact on the development and progression of CCA based on the mechanism of competing endogenous RNAs (ceRNAs). However, functional roles and regulatory mechanisms of lncRNA-regulated ceRNA in CCA, are only partially understood. The expression profile of messenger RNAs (mRNAs), lncRNAs, and microRNAs (miRNAs) downloaded from The Cancer Genome Atlas were comprehensively investigated. Differential expression of these three types of RNA between CCA and corresponding precancerous tissues were screened out for further analysis. On the basis of interactive information generated from miRDB, miRTarBase, TargetScan, and miRcode public databases, we then constructed an mRNA-miRNA-lncRNA regulatory network. Kyoto Encyclopedia of Genes and Genomes and Gene Ontology analyses were conducted to identify the biological function of the ceRNA network involved in CCA. As a result, 2883 mRNAs, 136 miRNAs, and 993 lncRNAs were screened out as differentially expressed RNAs in CCA. In addition, a ceRNA network in CCA was constructed, composing of 50 up and 27 downregulated lncRNAs, 14 up and 7 downregulated miRNAs, 29 up and 25 downregulated mRNAs. Finally, gene set enrichment and pathway analysis indicated our CCA-specific ceRNA network was related with cancer-related pathway and molecular function. In conclusion, our research identified a novel lncRNA-related ceRNA network in CCA, which might act as a potential therapeutic target for patients with CCA.     

Identification and characterization of tumorigenic circular RNAs in cervical cancer

Circular RNAs (circRNAs) are a distinctive family of ncRNAs, and they function as key regulators in the initiation, development and progression of various diseases. However, the regulatory roles of circRNAs in the tumorigenesis of cervical cancer (CC) have not been fully understood. In this study, we identified a set of circRNAs in CC and paired normal tissues, using RNA sequencing data, and found that the cancer and normal tissues could be told apart by those differentially expressed (DE) circRNAs, indicating that circRNA expression profiles in CC were significantly different from those in the normal tissues. Meanwhile, the upregulated genes in CC were enriched in inflammation-related pathways, and the correlation analysis between the DE circRNAs and genes revealed that the abundance of DE circRNAs was positively related to the expression of their host genes. However, the expression of TGFBR2 and KDM4C were found to exhibit a negative correlation with their corresponding circRNAs. Furthermore, we also predicted the interactions between circRNAs and proteins, and constructed a competing endogenous RNA (ceRNA) network. Specifically, hsa_circ_0001495 was predicted to interact with ESRP2, and acted as a sponge by competing for miRNAs with TBL1XR1. Functionally, hsa_circ_0001495 was predicted to regulate epithelial cell proliferation and NOTCH signaling via ESRP2 and TBL1XR1, respectively. We also evaluated the prognostic values of downstream target genes of selected circRNAs, using clinical records of CC patients. In summary, the present study provided some regulatory circRNAs involved in CC tumorigenesis based on bioinformatics approaches, which brought strong evidences for researchers to further explore their biological and clinical values.