Identification of lncRNA‐associated differential subnetworks in oesophageal squamous cell carcinoma by differential co‐expression analysis

Differential expression analysis has led to the identification of important biomarkers in oesophageal squamous cell carcinoma (ESCC). Despite enormous contributions, it has not harnessed the full potential of gene expression data, such as interactions among genes. Differential co‐expression analysis has emerged as an effective tool that complements differential expression analysis to provide better insight of dysregulated mechanisms and indicate key driver genes. Here, we analysed the differential co‐expression of lncRNAs and protein‐coding genes (PCGs) between normal oesophageal tissue and ESCC tissues, and constructed a lncRNA‐PCG differential co‐expression network (DCN). DCN was characterized as a scale‐free, small‐world network with modular organization. Focusing on lncRNAs, a total of 107 differential lncRNA‐PCG subnetworks were identified from the DCN by integrating both differential expression and differential co‐expression. These differential subnetworks provide a valuable source for revealing lncRNA functions and the associated dysfunctional regulatory networks in ESCC. Their consistent discrimination suggests that they may have important roles in ESCC and could serve as robust subnetwork biomarkers. In addition, two tumour suppressor genes (AL121899.1 and ELMO2), identified in the core modules, were validated by functional experiments. The proposed method can be easily used to investigate differential subnetworks of other molecules in other cancers.     

Long noncoding RNA HAND2-AS1 inhibits cancer cell proliferation,migration, and invasion in esophagus squamous cell carcinoma by regulating microRNA-21

Long noncoding RNA (lncRNA) HAND2-AS1 is a well-characterized tumor suppressor in several types of malignancies, while its role in esophagus squamous cell carcinoma (ESCC) is unknown. In this study, we found that lncRNA HAND2-AS1 was downregulated, while microRNA-21 ( miRNA-21) was upregulated in tumor tissues than in adjacent healthy tissues of ESCC patients. Expression levels of lncRNA HAND2-AS1 and miRNA-21 were significantly and inversely correlated in tumor tissues but not in healthy tissues. Plasma levels of lncRNA HAND2-AS1 were lower in ESCC patients than in healthy controls, and downregulation of plasma lncRNA HAND2-AS1 distinguished early stage ESCC patients from healthy controls. lncRNA HAND2-AS1 overexpression resulted in downregulation of miRNA-21 in cells of ESCC cell lines and inhibited cell proliferation, migration, and invasion. miRNA-21 overexpression failed to affect lncRNA HAND2-AS1 expression but significantly attenuated the inhibitory effect of lncRNA HAND2-AS1 overexpression on cancer cell proliferation, migration, and invasion. Therefore, lncRNA HAND2-AS1 may inhibit cancer cell proliferation, migration, and invasion in ESCC by regulating miRNA-21.     

LncRNA MEG3 inhibits the progression of prostate cancer by modulating miR‐9‐5p/QKI‐5 axis

This study was designed to detecting the influences of lncRNA MEG3 in prostate cancer. Aberrant lncRNAs expression profiles of prostate cancer were screened by microarray analysis. The qRT‐PCR and Western blot were employed to investigating the expression levels of lncRNA MEG3, miR‐9‐5p and QKI‐5. The luciferase reporter assay was utilized to testifying the interactions relationship among these molecules. Applying CCK‐8 assay, wound healing assay, transwell assay and flow cytometry in turn, the cell proliferation, migration and invasion abilities as well as apoptosis were measured respectively. LncRNA MEG3 was a down‐regulated lncRNA in prostate cancer tissues and cells and could inhibit the expression of miR‐9‐5p, whereas miR‐9‐5p down‐regulated QKI‐5 expression. Overexpressed MEG3 and QKI‐5 could decrease the abilities of proliferation, migration and invasion in prostate cancer cells effectively and increased the apoptosis rate. On the contrary, miR‐9‐5p mimics presented an opposite tendency in prostate cancer cells. Furthermore, MEG3 inhibited tumour growth and up‐regulated expression of QKI‐5 in vivo. LncRNA MEG3 was a down‐regulated lncRNA in prostate cancer and impacted the abilities of cell proliferation, migration and invasion, and cell apoptosis rate, this regulation relied on regulating miR‐9‐5p and its targeting gene QKI‐5.     

Regulation of laryngeal squamous cell cancer progression by the lncRNA RP11‐159K7.2/miR‐206/DNMT3A axis

Long non‐coding RNAs (lncRNAs), which are longer than 200 nt, have been proved to play a role in promoting or inhibiting cancer progression. The following study investigated the role and underlying mechanisms of lncRNA RP11‐159K7.2 in laryngeal squamous cell carcinoma (LSCC) progression. Briefly, in situ hybridization (ISH) and real‐time quantitative PCR (RT‐qPCR) showed higher expression of RP11‐159K7.2 in LSCC tissues and cell lines. Patients with low expression level of RP11‐159K7.2 lived longer compared to those with high expression of RP11‐159K7.2 (χ² = 39.111, ***P < 0.001). Multivariate Cox regression analysis suggested that lncRNA RP11‐159K7.2 was an independent prognostic factor for LSCC patients (HR = 2.961, ***P < 0.001). Furthermore, to investigate the potential involvement of RP11‐159K7.2 in the development of LSCC, we knocked out the expression of endogenous RP11‐159K7.2 in TU‐212 cells and AMC‐HN‐8 cells via CRISPR/Cas9 double vector lentiviral system. RP11‐159K7.2 knockout decreased LSCC cell growth and invasion both in vitro and in vivo. Mechanically, we found that RP11‐159K7.2 could positively regulate the expression of DNMT3A by sponging miR‐206. In addition, a feedback loop was also discovered between DNMT3A and miR‐206. To sum up, these findings suggest that lncRNA RP11‐159K7.2 could be used as a potential biomarker for prognosis and treatment of LSCC.     

Comprehensive bioinformatics analysis identifies lncRNA HCG22 as a migration inhibitor in esophageal squamous cell carcinoma

Esophageal cancer is one of the most lethal malignancies worldwide, and esophageal squamous cell carcinoma (ESCC) is the dominant histological type. However, the long noncoding RNA (lncRNA) alterations in ESCC have not been elucidated to date. In this study, reliable databases from Gene Expression Omnibus (GEO), which analyzed lncRNA expression in ESCC tumor tissues and adjacent normal tissues were searched, and common differentially expressed lncRNAs and genes were analyzed. Next, cis- trans analysis was performed to predict the underlying relationships between altered lncRNAs and mRNAs, and the lncRNA-mRNA regulatory network was established. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses of altered lncRNA-related genes were performed. The promising lncRNA HCG22 was validated by quantitative polymerase chain reaction (qPCR), and clinicopathological data were collected to identify the relationship between lncRNA HCG22 expression level and clinical features. Finally, Transwell assays were performed to explore the biological functions of lncRNA HCG22 in ESCC cells. Two hundred forty-one lncRNAs and 835 mRNAs were observed to be remarkably altered between ESCC tumor tissues and adjacent normal tissues. The lncRNA-mRNA regulatory network showed the coexpression association between lncRNA HCG22 and SPINK7 and ADAMTS12. GO and KEGG analyses showed that HCG22 and ADAMTS12 had potential biological functions in the cell migration of ESCC. The downregulation of lncRNA HCG22 in ESCC tumor tissues was validated by qPCR, and the clinicopathological data showed a noticeable correlation between lncRNA HCG22 expression level and the ESCC differentiational degree and clinical TNM stage. Kaplan-Meier analysis showed that patients with ESCC having low lncRNA HCG22 expression in ESCC tissues had considerably shorter overall survival compared with patients with ESCC having high lncRNA HCG22 expression. Following Transwell assays confirmed the migratory role of lncRNA HCG22 in ESCC cells. In conclusion, lncRNA HCG22 was downregulated in ESCC tissues and can be a migration inhibitor of ESCC cells, and SPINK7 and ADAMTS12 are promising to be the regulatory targets of lncRNA HCG22.     

LncRNA ABHD11‐AS1 promotes the development of endometrial carcinoma by targeting cyclin D1

To investigate the expression, role and mechanism of action of long non‐coding RNA (lncRNA) ABHD11‐AS1 in endometrial carcinoma. The expression of lncRNA ABHD11‐AS1 was quantified by qRT‐PCR in human endometrial carcinoma (n = 89) and normal endometrial tissues (n = 27). LncRNA ABHD11‐AS1 was stably overexpressed or knocked‐down in endometrial carcinoma cell lines to examine the cellular phenotype and expression of related molecules. Compared to normal endometrial tissue, lncRNA ABHD11‐AS1 was significantly overexpressed in endometrial carcinoma. Overexpression of lncRNA ABHD11‐AS1 promoted the proliferation, G1‐S progression, invasion and migration of endometrial cancer cells; inhibited apoptosis; up‐regulated cyclin D1, CDK1, CDK2, CDK4, Bcl‐xl and VEGFA; and down‐regulated p16, while ABHD11‐AS1 down‐regulation has the opposite effect. RNA pull down demonstrated that lncRNA ABHD11‐AS1 binds directly to cyclin D1. Knockdown of cyclin D1 can reverse the effect of ABHD11‐AS1. Overexpression of lncRNA ABHD11‐AS1 increased the tumorigenicity and up‐regulated cyclin D1 in an in vivo model of endometrial cancer in nude mice. LncRNA ABHD11‐AS1 functions as an oncogene to promote cell proliferation and invasion in endometrial carcinoma by positively targeting cyclin D1.     

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.     

Functional analysis of RPS27 mutations and expression in melanoma

Next‐generation sequencing has enabled genetic and genomic characterization of melanoma to an unprecedent depth. However, the high mutational background plus the limited depth of coverage of whole‐genome sequencing performed on cutaneous melanoma samples make the identification of novel driver mutations difficult. We sought to explore the somatic mutation portfolio in exonic and gene regulatory regions in human melanoma samples, for which we performed targeted sequencing of tumors and matched germline DNA samples from 89 melanoma patients, identifying known and novel recurrent mutations. Two recurrent mutations found in the RPS27 promoter associated with decreased RPS27 mRNA levels in vitro. Data mining and IHC analyses revealed a bimodal pattern of RPS27 expression in melanoma, with RPS27‐low patients displaying worse prognosis. In vitro characterization of RPS27‐high and RPS27‐low melanoma cell lines, as well as loss‐of‐function experiments, demonstrated that high RPS27 status provides increased proliferative and invasive capacities, while low RPS27 confers survival advantage in low attachment and resistance to therapy. Additionally, we demonstrate that 10 other cancer types harbor bimodal RPS27 expression, and in those, similarly to melanoma, RPS27‐low expression associates with worse clinical outcomes. RPS27 promoter mutation could thus represent a mechanism of gene expression modulation in melanoma patients, which may have prognostic and predictive implications.     

Genome‐wide characterization and developmental expression profiling of long non‐coding RNAs in Sogatella furcifera

The genome‐wide characterization of long non‐coding RNA (lncRNA) in insects demonstrates their importance in fundamental biological processes. Essentially, an in‐depth understanding of the functional repertoire of lncRNA in insects is pivotal to insect resources utilization and sustainable pest control. Using a custom bioinformatics pipeline, we identified 1861 lncRNAs encoded by 1852 loci in the Sogatella furcifera genome. We profiled lncRNA expression in different developmental stages and observed that the expression of lncRNAs is more highly temporally restricted compared to protein‐coding genes. More up‐regulated Sogatella furcifera lncRNA expressed in the embryo, 4th and 5th instars, suggesting that increased lncRNA levels may play a role in these developmental stages. We compared the relationship between the expression of Sogatella furcifera lncRNA and its nearest protein gene and found that lncRNAs were more correlated to their downstream coding neighbors on the opposite strand. Our genome‐wide profiling of lncRNAs in Sogatella furcifera identifies exciting candidates for characterization of lncRNAs, and also provides information on lncRNA regulation during insect development.     

MiR‐130b promotes the progression of oesophageal squamous cell carcinoma by targeting SASH1

MiR‐130b and SAM and SH3 domain containing 1 (SASH1) play an important role in many types of human cancers. The aim of our research was to study their interactions in the process of the proliferation and aggressiveness of oesophageal squamous cell carcinoma (ESCC) cells. Microarray analysis was done to screen the differentially expressed genes in the ESCC tissues. miR‐130b and SASH1 mRNA levels in the ESCC tissues and cells were detected by qRT‐PCR. Dual luciferase reporter system was used to verify the target relationship between miR‐130b and SASH1. The effects of miR‐130b on SASH1 expression were explored by western blot in KYSE30 and TE1 cell lines. CCK‐8 assay, flow cytometry, Transwell, and wound healing assays were conducted to explore the effects of miR‐130b and SASH1 in vitro. In addition, in vivo experiments were conducted to study the roles of miR‐130b and SASH1. miR‐130b was highly expressed, while SASH1 was the opposite in both the ESCC tissues and cells. The expression of SASH1 was inhibited by the direct binding of miR‐130b. The inhibition of miR‐130b reduced the proliferation and aggressiveness of ESCC cells, while it also induced apoptosis and cell cycle arrest in the ESCC cells by suppressing SASH1. The in vivo assay suggested that the overexpression of miR‐130b promoted the growth of ESCC tumours. MiR‐130b was up‐regulated in the ESCC tumour tissues and cells, acting as a tumour promoter. A stimulating effect was demonstrated on ESCC cell growth and aggressiveness by suppressing SASH1, which is an anti‐oncogene.     

LNRLMI: Linear neighbour representation for predicting lncRNA‐miRNA interactions

LncRNA and miRNA are key molecules in mechanism of competing endogenous RNAs(ceRNA), and their interactions have been discovered with important roles in gene regulation. As supplementary to the identification of lncRNA‐miRNA interactions from CLIP‐seq experiments, in silico prediction can select the most potential candidates for experimental validation. Although developing computational tool for predicting lncRNA‐miRNA interaction is of great importance for deciphering the ceRNA mechanism, little effort has been made towards this direction. In this paper, we propose an approach based on linear neighbour representation to predict lncRNA‐miRNA interactions (LNRLMI). Specifically, we first constructed a bipartite network by combining the known interaction network and similarities based on expression profiles of lncRNAs and miRNAs. Based on such a data integration, linear neighbour representation method was introduced to construct a prediction model. To evaluate the prediction performance of the proposed model, k‐fold cross validations were implemented. As a result, LNRLMI yielded the average AUCs of 0.8475 ± 0.0032, 0.8960 ± 0.0015 and 0.9069 ± 0.0014 on 2‐fold, 5‐fold and 10‐fold cross validation, respectively. A series of comparison experiments with other methods were also conducted, and the results showed that our method was feasible and effective to predict lncRNA‐miRNA interactions via a combination of different types of useful side information. It is anticipated that LNRLMI could be a useful tool for predicting non‐coding RNA regulation network that lncRNA and miRNA are involved in.     

Long non‐coding RNA LINC01535 promotes cervical cancer progression via targeting the miR‐214/EZH2 feedback loop

Long non‐coding RNAs (lncRNAs) have shown critical roles in multiple cancers via competitively binding common microRNAs. miR‐214 has been proved to play tumour suppressive roles in various cancers, including cervical cancer. In this study, we identified that lncRNA LINC01535 physically binds miR‐214, relieves the repressive roles of miR‐214 on its target EZH2, and therefore up‐regulates EZH2 protein expression. Intriguingly, we also found that EZH2 directly represses the expression of miR‐214. Thus, miR‐214 and EZH2 form double negative regulatory loop. Through up‐regulating EZH2, LINC01535 further represses miR‐214 expression. Functional experiments showed that enhanced expression of LINC01535 promotes cervical cancer cell growth, migration and invasion in vitro and cervical cancer xenograft growth in vivo. Reciprocally, LINC01535 knockdown suppresses cervical cancer cell growth, migration and invasion. Activation of the miR‐214/EZH2 regulatory loop by overexpression of miR‐214 or silencing of EZH2 reverses the roles of LINC01535 in promoting cervical canc`er cell growth, migration and invasion in vitro and cervical cancer xenograft growth in vivo. Clinically, LINC01535 is significantly up‐regulated in cervical cancer tissues and correlated with advanced clinical stage and poor prognosis. Moreover, the expression of LINC01535 is reversely associated with the expression of miR‐214 and positively associated with the expression of EZH2 in cervical cancer tissues. In conclusion, this study reveals that LINC01535 promotes cervical cancer progression via repressing the miR‐214/EZH2 regulatory loop.     

Long non‐coding RNA SNHG14 induces trastuzumab resistance of breast cancer via regulating PABPC1 expression through H3K27 acetylation

Currently, resistance to trastuzumab, a human epidermal growth factor receptor 2 (HER2) inhibitor, has become one major obstacle for improving the clinical outcome of patients with advanced HER2⁺ breast cancer. While cell behaviour can be modulated by long non‐coding RNAs (lncRNAs), the contributions of lncRNAs in progression and trastuzumab resistance of breast cancer are largely unknown. To this end, the involvement and regulatory functions of lncRNA SNHG14 in human breast cancer were investigated. RT‐qPCR assay showed that SNHG14 was up‐regulated in breast cancer tissues and associated with trastuzumab response. Gain‐ and loss‐of‐function experiments revealed that overexpression of SNHG14 promotes cell proliferation, invasion and trastuzumab resistance, whereas knockdown of SNHG14 showed an opposite effect. PABPC1 gene was identified as a downstream target of SNHG14, and PABPC1 mediates the SNHG14‐induced oncogenic effects. More importantly, ChIP assays demonstrated that lncRNA SNHG14 may induce PABPC1 expression through modulating H3K27 acetylation in the promoter of PABPC1 gene, thus resulting in the activation of Nrf2 signalling pathway. These data suggest that lncRNA SNHG14 promotes breast cancer tumorigenesis and trastuzumab resistance through regulating PABPC1 expression through H3K27 acetylation. Therefore, SNHG14 may serve as a novel diagnostic and therapeutic target for breast cancer patients.     

Atractylenolide II reverses the influence of lncRNA XIST/miR‐30a‐5p/ROR1 axis on chemo‐resistance of colorectal cancer cells

This investigation was conducted to elucidate whether atractylenolide II could reverse the role of lncRNA XIST/miR‐30a‐5p/ROR1 axis in modulating chemosensitivity of colorectal cancer cells. We totally collected 294 pairs of colorectal cancer tissues and adjacent normal tissues and also purchased colorectal cancer cell lines and human embryonic kidney cell line. 5‐fluorouracil, cisplatin, mitomycin and adriamycin were designated as the chemotherapies for colorectal cell lines, and atractylenolides were arranged as the Chinese drug. The expressions of XIST, miR‐30a‐5p and ROR1 were quantified with aid of qRT‐PCR or Western blot, and luciferase reporter gene assay was implemented to determine the relationships among XIST, miR‐30a‐5p and ROR1. Our results demonstrated that XIST and ROR1 expressions were dramatically up‐regulated, yet miR‐30a‐5p expression was down‐regulated within colorectal cancer tissues (P < 0.05). The overexpressed XIST and ROR1, as well as under‐expressed miR‐30a‐5p, were inclined to promote viability and proliferation of colorectal cells under the influence of chemo drugs (P < 0.05). In addition, XIST could directly target miR‐30a‐5p, and ROR1 acted as the targeted molecule of miR‐30a‐5p. Interestingly, atractylenolides not only switched the expressions of XIST, miR‐30a‐5p and ROR1 within colorectal cancer cells but also significantly intensified the chemosensitivity of colorectal cancer cells (P < 0.05). Finally, atractylenolide II was discovered to slow down the viability and proliferation of colorectal cancer cells (P < 0.05). In conclusion, the XIST/miR‐30a‐5p/ROR1 axis could be deemed as pivotal markers underlying colorectal cancer, and administration of atractylenolide II might improve the chemotherapeutic efficacy for colorectal cancer.