Identification of messenger and long noncoding RNAs associated with gallbladder cancer via gene expression profile analysis

The present study aimed to investigate the long noncoding RNAs (lncRNAs) and messenger RNAs (mRNAs) involved in the progression of gallbladder cancer and explore the potential physiopathologic mechanisms of gallbladder cancer in terms of competing endogenous RNAs (ceRNAs). The original lncRNA and mRNA expression profile data (nine gallbladder cancer tissues samples and nine normal gallbladder samples) in GSE76633 was downloaded from the Gene Expression Omnibus database. Differentially expressed mRNAs and lncRNAs between gallbladder cancer tissue and normal control were selected and the pathways in which they are involved were analyzed using bioinformatics analyses. MicroRNAs (miRNAs) were also predicted based on the differentially expressed mRNAs. Finally, the co-expression relation between lncRNA and mRNA was analyzed and the ceRNA network was constructed by combining the lncRNA-miRNA, miRNA-mRNA, and lncRNA-mRNA pairs. Overall, 373 significantly differentially expressed mRNAs and 47 lncRNAs were identified between cancer and normal tissue samples. The upregulated genes were significantly enriched in the extracellular matrix (ECM)-receptor interaction pathway, while the downregulated genes were involved in the complement and coagulation cascades. Altogether, 128 co-expression relations between lncRNA and mRNA were obtained. In addition, 196 miRNA-mRNA regulatory relations and 145 miRNA-lncRNA relation pairs were predicted. Finally, the lncRNA-miRNA-gene ceRNA network was constructed by combining the three types of relation pairs, such as XLOC_011309-miR-548c-3p-SPOCK1 and XLOC_012588-miR-765-CEACAM6. mRNAs and lncRNAs may be involved in gallbladder cancer progression via ECM-receptor interaction pathways and the complement and coagulation cascades. Moreover, ceRNAs such as XLOC_011309-miR-548c-3p-SPOCK1 and XLOC_012588-miR-765-CEACAM6 can also be implicated in the pathogenesis of gallbladder cancer.     

Sl‐lncRNA15492 interacts with Sl‐miR482a and affects Solanum lycopersicum immunity against Phytophthora infestans

Long non‐coding RNAs (lncRNAs) are involved in the resistance of plants to infection by pathogens via interactions with microRNAs (miRNAs). Long non‐coding RNAs are cleaved by miRNAs to produce phased small interfering RNAs (phasiRNAs), which, as competing endogenous RNAs (ceRNAs), function as decoys for mature miRNAs, thus inhibiting their expression, and contain pre‐miRNA sequences to produce mature miRNAs. However, whether lncRNAs and miRNAs mediate other molecular mechanisms during plant resistance to pathogens is unknown. In this study, as a positive regulator, Sl‐lncRNA15492 from tomato (Solanum lycopersicum Zaofen No. 2) plants affected tomato resistance to Phytophthora infestans. Gain‐ and loss‐of‐function experiments and RNA ligase‐mediated 5′‐amplification of cDNA ends (RLM‐5′ RACE) also revealed that Sl‐miR482a was negatively involved in tomato resistance by targeting Sl‐NBS‐LRR genes and that silencing of Sl‐NBS‐LRR1 decreased tomato resistance. Sl‐lncRNA15492 inhibited the expression of mature Sl‐miR482a, whose precursor was located within the antisense sequence of Sl‐lncRNA15492. Further degradome analysis and additional RLM‐5′ RACE experiments verified that mature Sl‐miR482a could also cleave Sl‐lncRNA15492. These results provide a mechanism by which lncRNAs might inhibit precursor miRNA expression through antisense strands of lncRNAs, and demonstrate that Sl‐lncRNA15492 and Sl‐miR482a mutually inhibit the maintenance of Sl‐NBS‐LRR1 homeostasis during tomato resistance to P. infestans.     

Oestrogen induces epithelial‐mesenchymal transition in endometriosis via circ_0004712/miR‐148a‐3p sponge function

Endometriosis is a common, chronic gynaecologic disease affecting up to 10% of women in their reproductive age and leading to pain and infertility. Oestrogen (E₂)‐induced epithelial‐mesenchymal transition (EMT) process has been considered as a key factor of endometriosis development. Recently, the dysregulated circular RNAs (circRNAs) have been discovered in endometriosis tissues. However, the molecular mechanism of circRNAs on the E₂‐induced EMT process in endometriosis is still unknown. Here, we demonstrated that circ_0004712 up‐regulated by E₂ treatment in endometrial epithelial cells. Knock‐down the expression of circ_0004712 significantly suppressed E₂‐induced cell migration activity. Meanwhile, we identified miR‐148a‐3p as a potential target miRNA of circ_0004712. Inhibited the expression of miR‐148a‐3p could recovered the effect of circ_0004712 knock‐down in E₂‐treated endometrial epithelial. Furthermore, Western blot assay showed that E₂ treatment could increase the expression and activity of β‐catenin, snail and N‐cadherin and reduce the expression of E‐cadherin. The expression and activity of β‐catenin pathway were recovered by circ_0004712 knock‐down or miR‐148a‐3p overexpression. Altogether, the results demonstrate that circ_0004712/miR‐148a‐3p plays an important role in E₂‐induced EMT process in the development of endometriosis, and the molecular mechanism may be associated with the β‐catenin pathway. This work highlighted the importance of circRNAs in the development of endometriosis and provide a new biomarker for diagnosis and therapies.     

Involvement of a novel circularRNA,hsa_circ_0000520, attenuates tumorigenesis of cervical cancer cell through competitively binding with miR‐146b‐3p

The implication of circular RNAs (circRNAs) in the pathogenesis of human cervical cancer (CC) has been demonstrated by numerous of researches, nevertheless, the whole regulatory network of circRNAs in CC remains unclear. In the present study, two GSE data sets (GSE113696 and GSE102686) were enrolled to analysed different expressed circRNA. We found that hsa_circ_0000520(circ_0000520) was decreased in CC tissues and cell lines. Functional studies indicated circ_0000520 overexpression in vitro repressed CC cell proliferation, invasion and migration, while promoted CC cell apoptosis. Moreover, circ_0000520 overexpression in vivo repressed CC tumour growth. Mechanismly, circ_0000520 and PAX5 were revealed to directly bind to miR‐146b‐3p, and circ_0000520 could indirectly regulate PAX5 by sponging miR‐146b‐3p. In conclusion, c irc_0000520 repressed CC progression in vitro and in vivo by sponging miR‐146b‐3p to release PAX5.     

Long noncoding RNA MALAT1 enhances the docetaxel resistance of prostate cancer cells via miR‐145‐5p‐mediated regulation of AKAP12

Our present work was aimed to study on the regulatory role of MALAT1/miR‐145‐5p/AKAP12 axis on docetaxel (DTX) sensitivity of prostate cancer (PCa) cells. The microarray data (GSE33455) to identify differentially expressed lncRNAs and mRNAs in DTX‐resistant PCa cell lines (DU‐145‐DTX and PC‐3‐DTX) was retrieved from the Gene Expression Omnibus (GEO) database. QRT‐PCR analysis was performed to measure MALAT1 expression in DTX‐sensitive and DTX‐resistant tissues/cells. The human DTX‐resistant cell lines DU145‐PTX and PC3‐DTX were established as in vitro cell models, and the expression of MALAT1, miR‐145‐5p and AKAP12 was manipulated in DTX‐sensitive and DTX‐resistant cells. Cell viability was examined using MTT assay and colony formation methods. Cell apoptosis was assessed by TUNEL staining. Cell migration and invasion was determined by scratch test (wound healing) and Transwell assay, respectively. Dual‐luciferase assay was applied to analyse the target relationship between lncRNA MALAT1 and miR‐145‐5p, as well as between miR‐145‐5p and AKAP12. Tumour xenograft study was undertaken to confirm the correlation of MALAT1/miR‐145‐5p/AKAP12 axis and DTX sensitivity of PCa cells in vivo. In this study, we firstly notified that the MALAT1 expression levels were up‐regulated in clinical DTX‐resistant PCa samples. Overexpressed MALAT1 promoted cell proliferation, migration and invasion but decreased cell apoptosis rate of PCa cells in spite of DTX treatment. We identified miR‐145‐5p as a target of MALAT1. MiR‐145‐5p overexpression in PC3‐DTX led to inhibited cell proliferation, migration and invasion as well as reduced chemoresistance to DTX, which was attenuated by MALAT1. Moreover, we determined that AKAP12 was a target of miR‐145‐5p, which significantly induced chemoresistance of PCa cells to DTX. Besides, it was proved that MALAT1 promoted tumour cell proliferation and enhanced DTX‐chemoresistance in vivo. There was an lncRNA MALAT1/miR‐145‐5p/AKAP12 axis involved in DTX resistance of PCa cells and provided a new thought for PCa therapy.     

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.     

Analysing the relationship between lncRNA and protein‐coding gene and the role of lncRNA as ceRNA in pulmonary fibrosis

Long non‐coding RNAs (lncRNAs) are involved in various pathophysiologic processes and human diseases. However, their dynamics and corresponding functions in pulmonary fibrosis remain poorly understood. In this study, portions of lncRNAs adjacent or homologous to protein‐coding genes were determined by searching the UCSC genome bioinformatics database. This was found to be potentially useful for exploring lncRNA functions in disease progression. Previous studies showed that competing endogenous RNA (ceRNA) hypothesis is another method to predict lncRNA function. However, little is known about the function of ceRNA in pulmonary fibrosis. In this study, we selected two differentially expressed lncRNAs MRAK088388 and MRAK081523 to explore their regulatory mechanisms. MRAK088388 and MRAK081523 were analysed as long‐intergenic non‐coding RNAs (lincRNAs), and identified as orthologues of mouse lncRNAs AK088388 and AK081523, respectively. qRT‐PCR and in situ hybridization (ISH) showed that they were significantly up‐regulated, and located in the cytoplasm of interstitial lung cells. We also showed that MRAK088388 and N4bp2 had the same miRNA response elements (MREs) for miR‐200, miR‐429, miR‐29, and miR‐30, whereas MRAK081523 and Plxna4 had the same MREs for miR‐218, miR‐141, miR‐98, and let‐7. Moreover, the expression levels of N4bp2 and Plxna4 significantly increased in fibrotic rats, and were highly correlated with those of MRAK088388 and MRAK081523, respectively. Among their shared miRNAs, miR‐29b‐3p and let‐7i‐5p decreased in the model group, and were negatively correlated with the expression of MRAK088388 and MRAK081523, respectively. MRAK088388 and MRAK081523 could regulate N4bp2 and Plxna4 expression by sponging miR‐29b‐3p and let‐7i‐5p, respectively, and possessed regulatory functions as ceRNAs. Thus, our study may provide insights into the functional interactions of lncRNA, miRNA and mRNA, and lead to new theories for the pathogenesis and treatment of pulmonary fibrosis.     

MiRNA‐145‐5p prevents differentiation of oligodendrocyte progenitor cells by regulating expression of myelin gene regulatory factor

The roles of specific microRNAs (miRNA) in oligodendrocyte (OL) differentiation have been studied in depth. However, miRNAs in OL precursors and oligodendrocyte progenitor cells (OPCs) have been less extensively investigated. MiR‐145‐5p is highly expressed in OPCs relative to differentiating OLs, suggesting this miRNA may serve a function specifically in OPCs. Knockdown of miR‐145‐5p in primary OPCs led to spontaneous differentiation, as evidenced by an increased proportion of MAG⁺ cells, increased cell ramification, and upregulation of multiple myelin genes including MYRF, TPPP, and MAG, and OL cell cycle exit marker Cdkn1c. Supporting this transition to a differentiating state, proliferation was reduced in miR‐145‐5p knockdown OPCs. Further, knockdown of miR‐145‐5p in differentiating OLs showed enhanced differentiation, with increased branching, myelin membrane production, and myelin gene expression. We identified several OL‐specific genes targeted by miR‐145‐5p that exhibited upregulation with miR‐145‐5p knockdown, including myelin gene regulatory factor (MYRF), that could be regulating the prodifferentiation phenotype in both miR‐145 knockdown OPCs and OLs. Indeed, spontaneous differentiation with knockdown of miR‐145‐5p was fully rescued by concurrent knockdown of MYRF. However, proliferation rate was only partially rescued with MYRF knockdown, and overexpression of miR‐145‐5p in OPCs increased proliferation rate without affecting expression of already lowly expressed differentiation genes. Taken together, these data suggest that in OPCs miR‐145‐5p both prevents differentiation at least in part by preventing expression of MYRF and promotes proliferation via as‐yet‐unidentified mechanisms. These findings clarify the need for differential regulation of miR‐145‐5p between OPCs and OLs and may have further implications in demyelinating diseases such as multiple sclerosis where miR‐145‐5p is dysregulated.     

LncRNA‐mRNA expression profiles and functional networks in osteoclast differentiation

Human osteoclasts are differentiated from CD14⁺ monocytes and are responsible for bone resorption. Long non‐coding RNAs (lncRNAs) have been proved to be significantly involved in multiple biologic processes, especially in cell differentiation. However, the effect of lncRNAs in osteoclast differentiation is less appreciated. In our study, RNA sequencing (RNA‐seq) was used to identify the expression profiles of lncRNAs and mRNAs in osteoclast differentiation. The results demonstrated that expressions of 1117 lncRNAs and 296 mRNAs were significantly altered after osteoclast differentiation. qRT‐PCR assays were performed to confirm the expression profiles, and the results were almost consistent with the RNA‐seq data. GO and KEGG analyses were used to predict the functions of these differentially expressed mRNA and lncRNAs. The Path‐net analysis demonstrated that MAPK pathway, PI3K‐AKT pathway and NF‐kappa B pathway played important roles in osteoclast differentiation. Co‐expression networks and competing endogenous RNA networks indicated that ENSG00000257764.2‐miR‐106a‐5p‐TIMP2 may play a central role in osteoclast differentiation. Our study provides a foundation to further understand the role and underlying mechanism of lncRNAs in osteoclast differentiation, in which many of them could be potential targets for bone metabolic disease.     

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.     

Hormone replacement therapy enhances IGF‐1 signaling in skeletal muscle by diminishing miR‐182 and miR‐223 expressions: a study on postmenopausal monozygotic twin pairs

MiRNAs are fine‐tuning modifiers of skeletal muscle regulation, but knowledge of their hormonal control is lacking. We used a co‐twin case–control study design, that is, monozygotic postmenopausal twin pairs discordant for estrogen‐based hormone replacement therapy (HRT) to explore estrogen‐dependent skeletal muscle regulation via miRNAs. MiRNA profiles were determined from vastus lateralis muscle of nine healthy 54–62‐years‐old monozygotic female twin pairs discordant for HRT (median 7 years). MCF‐7 cells, human myoblast cultures and mouse muscle experiments were used to confirm estrogen's causal role on the expression of specific miRNAs, their target mRNAs and proteins and finally the activation of related signaling pathway. Of the 230 miRNAs expressed at detectable levels in muscle samples, qPCR confirmed significantly lower miR‐182, miR‐223 and miR‐142‐3p expressions in HRT using than in their nonusing co‐twins. Insulin/IGF‐1 signaling emerged one common pathway targeted by these miRNAs. IGF‐1R and FOXO3A mRNA and protein were more abundantly expressed in muscle samples of HRT users than nonusers. In vitro assays confirmed effective targeting of miR‐182 and miR‐223 on IGF‐1R and FOXO3A mRNA as well as a dose‐dependent miR‐182 and miR‐223 down‐regulations concomitantly with up‐regulation of FOXO3A and IGF‐1R expression. Novel finding is the postmenopausal HRT‐reduced miRs‐182, miR‐223 and miR‐142‐3p expression in female skeletal muscle. The observed miRNA‐mediated enhancement of the target genes' IGF‐1R and FOXO3A expression as well as the activation of insulin/IGF‐1 pathway signaling via phosphorylation of AKT and mTOR is an important mechanism for positive estrogen impact on skeletal muscle of postmenopausal women.     

CircNr1h4 regulates the pathological process of renal injury in salt‐sensitive hypertensive mice by targeting miR‐155‐5p

Circular RNAs are a class of widespread and diverse endogenous RNAs that may regulate gene expression in various diseases, but their regulation and function in hypertensive renal injury remain unclear. In this study, we generated ribosomal‐depleted RNA sequencing data from normal mouse kidneys and from injured mouse kidneys induced by deoxycorticosterone acetate‐salt hypertension and identified at least 4900 circRNA candidates. A total of 124 of these circRNAs were differentially expressed between the normal and injured kidneys. Furthermore, we characterized one abundant circRNA, termed circNr1h4, which is derived from the Nr1h4 gene and significantly down‐regulated in the injured kidneys. RNA sequencing data and qPCR analysis also showed many microRNAs and mRNAs, including miR‐155‐5p and fatty acid reductase 1 (Far1), were differentially expressed between the normal and injured kidney and related to circNr1h4. In vitro, the silencing of circNr1h4 or overexpression of miR‐155‐5p significantly decreased Far1 levels and increased reactive oxygen species. Mechanistic investigations indicated that circNr1h4 acts as a competing endogenous RNA for miR‐155‐5p, leading to regulation of its target gene Far1. Our study provides novel insight into the molecular mechanisms underlying kidney injury in hypertension, which will be required to develop therapeutic strategies of targeting circRNAs for hypertensive kidney injury.     

Differential circular RNAs expression in ovary during oviposition in honey bees

Circular RNAs (circRNAs) are non-coding RNAs newly identified and play important roles in RNA regulation. The mechanism and function of circRNAs have been reported in some species. However, little is known regarding circRNAs in honey bees. In this study, we analyzed circRNAs through bioinformatics, and predicted 12,211 circRNAs in the ovary of honey bee queens. 1340, 175 and 100 circRNAs were differentially expressed in comparisons of egg-laying queens vs virgin queens, egg-laying inhibited queens vs egg-laying queens and egg-laying recovery queens vs egg-laying inhibited queens. Further, functional annotation of differentially expressed circRNAs revealed several pathways that are closely related to ovary activation and oviposition, including insulin secretion and calcium signaling pathways. Moreover, the potential interactions among circRNAs, miRNAs, lncRNAs and mRNAs were investigated. Ame_circ_0005197 and ame_circ_0016640 were observed to sponge several reproductive related miRNAs. These findings demonstrate that circRNAs have potential effects in ovary activation and oviposition of honey bees.     

Circular RNA circ_001842 plays an oncogenic role in renal cell carcinoma by disrupting microRNA‐502‐5p‐mediated inhibition of SLC39A14

Renal cell carcinoma (RCC) is a common urologic malignancy, and up to 30% of RCC patients present with locally advanced or metastatic disease at the time of initial diagnosis. Increasing evidence suggests that circular RNAs (circRNAs) serve as genomic regulatory molecules in various human cancers. Our initial in silico microarray‐based analysis identified that circRNA circ_001842 was highly expressed in RCC. Such up‐regulation of circ_001842 in RCC was experimentally validated in tissues and cell lines using RT‐qPCR. Thereafter, we attempted to identify the role of circ_001842 in the pathogenesis of RCC. Through a series of gain‐ and loss‐of function assays, cell biological functions were examined using colony formation assay, Transwell assay, annexin V‐FITC/PI‐labelled flow cytometry and scratch test. A high expression of circ_001842 in tissues was observed as associated with poor prognosis of RCC patients. circ_001842 was found to elevate SLC39A14 expression by binding to miR‐502‐5p, consequently resulting in augmented RCC cell proliferation, migration and invasion, as well as EMT in vitro and tumour growth in vivo. These observations imply the involvement of circ_001842 in RCC pathogenesis through a miR‐502‐5p‐dependent SLC39A14 mechanism, suggesting circ_001842 is a potential target for RCC treatment.