A novel autophagy-related lncRNA prognostic risk model for breast cancer

Long non-coding RNAs (lncRNAs) are well known as crucial regulators to breast cancer development and are implicated in controlling autophagy. LncRNAs are also emerging as valuable prognostic factors for breast cancer patients. It is critical to identify autophagy-related lncRNAs with prognostic value in breast cancer. In this study, we identified autophagy-related lncRNAs in breast cancer by constructing a co-expression network of autophagy-related mRNAs-lncRNAs from The Cancer Genome Atlas (TCGA). We evaluated the prognostic value of these autophagy-related lncRNAs by univariate and multivariate Cox proportional hazards analyses and eventually obtained a prognostic risk model consisting of 11 autophagy-related lncRNAs (U62317.4, LINC01016, LINC02166, C6orf99, LINC00992, BAIAP2-DT, AC245297.3, AC090912.1, Z68871.1, LINC00578 and LINC01871). The risk model was further validated as a novel independent prognostic factor for breast cancer patients based on the calculated risk score by Kaplan-Meier analysis, univariate and multivariate Cox regression analyses and time-dependent receiver operating characteristic (ROC) curve analysis. Moreover, based on the risk model, the low-risk and high-risk groups displayed different autophagy and oncogenic statues by principal component analysis (PCA) and Gene Set Enrichment Analysis (GSEA) functional annotation. Taken together, these findings suggested that the risk model of the 11 autophagy-related lncRNAs has significant prognostic value for breast cancer and might be autophagy-related therapeutic targets in clinical practice.     

Identification of autophagy-related long non-coding RNA prognostic signature for breast cancer

Autophagy-related long non-coding RNAs (lncRNAs) disorders are related to the occurrence and development of breast cancer. The purpose of this study is to explore whether autophagy-related lncRNA can predict the prognosis of breast cancer patients. The autophagy-related lncRNAs prognostic signature was constructed by Least Absolute Shrinkage and Selection Operator (LASSO) Cox regression. We identified five autophagy-related lncRNAs (MAPT-AS1, LINC01871, AL122010.1, AC090912.1, AC061992.1) associated with prognostic value, and they were used to construct an autophagy-related lncRNA prognostic signature (ALPS) model. ALPS model offered an independent prognostic value (HR = 1.664, 1.381-2.006), where this risk score of the model was significantly related to the TNM stage, ER, PR and HER2 status in breast cancer patients. Nomogram could be utilized to predict survival for patients with breast cancer. Principal component analysis and Sankey Diagram results indicated that the distribution of five lncRNAs from the ALPS model tends to be low-risk. Gene set enrichment analysis showed that the high-risk group was enriched in autophagy and cancer-related pathways, and the low-risk group was enriched in regulatory immune-related pathways. These results indicated that the ALPS model composed of five autophagy-related lncRNAs could predict the prognosis of breast cancer patients.     

An eight-long noncoding RNA expression signature for colorectal cancer patients’ prognosis

Long noncoding RNAs (lncRNAs) have recently emerged as important biomarkers of cancer progression. Here, we proposed to develop a lncRNA-based signature with a prognostic value for colorectal cancer (CRC) overall survival (OS). Through mining microarray datasets, we analyzed the lncRNA expression profiles of 122 patients with CRC from Gene Expression Omnibus. Associations between lncRNA and CRC OS were firstly evaluated through univariate Cox regression analysis. A random survival forest method was applied for further screening of the lncRNA signature, which resulted in eight lncRNAs, including PEG3-AS1, LOC100505715, MINCR, DBH-AS1, LINC00664, FAM224A, LOC642852, and LINC00662. Combination of the eight lncRNAs weighted by their multivariate Cox regression coefficients formed a prognostic signature, through which, we could divide the 122 patients with CRC into two subgroups with significantly different OS. Good robustness of the lncRNA signature's prognostic value was verified through an independent data set consisting of 55 patients with CRC. In addition, gene set enrichment analysis indicated the potential association between high prognostic value and oxygen metabolism-related processes. This result should indicate that lncRNAs could be a useful signature for CRC prognosis.     

A novel autophagy-related lncRNA survival model for lung adenocarcinoma

Long non-coding RNA (lncRNA) is an important regulatory factor in the development of lung adenocarcinoma, which is related to the control of autophagy. LncRNA can also be used as a biomarker of prognosis in patients with lung adenocarcinoma. Therefore, it is important to determine the prognostic value of autophagy-related lncRNA in lung adenocarcinoma. In this study, autophagy-related mRNAs-lncRNAs were screened from lung adenocarcinoma and a co-expression network of autophagy-related mRNAs-lncRNAs was constructed by using The Cancer Genome Atlas (TCGA). The univariate and multivariate Cox proportional hazard analyses were used to evaluate the prognostic value of the autophagy-related lncRNAs and finally obtained a survival model composed of 11 autophagy-related lncRNAs. Through Kaplan-Meier analysis, univariate and multivariate Cox regression analysis and time-dependent receiver operating characteristic (ROC) curve analysis, it was further verified that the survival model was a new independent prognostic factor for patients with lung adenocarcinoma. In addition, based on the survival model, gene set enrichment analysis (GSEA) was used to illustrate the function of genes in low-risk and high-risk groups. These 11 lncRNAs were GAS6-AS1, AC106047.1, AC010980.2, AL034397.3, NKILA, AL606489.1, HLA-DQB1-AS1, LINC01116, LINC01806, FAM83A-AS1 and AC090559.1. The hazard ratio (HR) of the risk score was 1.256 (1.196-1.320) (P < .001) in univariate Cox regression analysis and 1.215 (1.149-1.286) (P < .001) in multivariate Cox regression analysis. And the AUC value of the risk score was 0.809. The 11 autophagy-related lncRNA survival models had important predictive value for the prognosis of lung adenocarcinoma and may become clinical autophagy-related therapeutic targets.     

Development and validation of a five-lncRNA signature with prognostic value in colon cancer

Dysregulation of long noncoding RNAs (lncRNAs) has been found in a large number of human cancers, including colon cancer. Therefore, the implementation of potential lncRNAs biomarkers with prognostic prediction value are very much essential. GSE39582 data set was downloaded from database of Gene Expression Omnibus. Re-annotation analysis of lncRNA expression profiles was performed by NetAffx annotation files. Univariate and multivariate Cox proportional analyses helped select prognostic lncRNAs. Algorithm of random survival forest-variable hunting (RSF-VH) together with stepwise multivariate Cox proportional analysis were performed to establish lncRNA signature. The log-rank test was carried out to analyze and compare the Kaplan-Meier survival curves of patients’ overall survival (OS). Receiver operating characteristic (ROC) analysis was used for comparing the survival prediction regarding its specificity and sensitivity based on lncRNA risk score, followed by calculating the values of area under the curve (AUC). The single-sample GSEA (ssGSEA) analysis was used to describe biological functions associated with this signature. Finally, to determine the robustness of this model, we used the validation sets including GSE17536 and The Cancer Genome Atlas data set. After re-annotation analysis of lncRNAs, a total of 14 lncRNA probes were obtained by univariate and multivariate Cox proportional analysis. Then, the RSF-VH algorithm and stepwise multivariate Cox analysis helped to build a five-lncRNA prognostic signature for colon cancer. The patients in group with high risk showed an obviously shorter survival time compared with patients in group with low risk with AUC of 0.75. In addition, the five-lncRNA signature can be used to independently predict the survival of patients with colon cancer. The ssGSEA analysis revealed that pathways such as extracellular matrix-receptor interaction was activated with an increase in risk score. These findings determined the strong power of prognostic prediction value of this five-lncRNA signature for colon cancer.     

Dysregulated ceRNA network modulated by copy number variation-driven lncRNAs in breast cancer: A comprehensive analysis

Breast cancer is a malignancy harmful to physical and mental health in women, with quite high mortality. Copy number variations (CNVs) are vital factors affecting the progression of breast cancer. Detecting CNVs in breast cancer to predict the prognosis of patients has become a promising approach to accurate treatment in recent years. The differential analysis was performed on CNVs of long noncoding RNAs (lncRNAs) as well as the expression of lncRNAs, microRNAs (miRNAs) and mRNAs in normal tissue and breast tumor tissue based on The Cancer Genome Atlas (TCGA) database. The CNV-driven lncRNAs were identified by the Kruskal–Wallis test. Meanwhile, a competitive endogenous RNA (ceRNA) network regulated by CNV-driven lncRNA was constructed. As the Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses revealed, the mRNAs in the dysregulated ceRNA network were mainly enriched in the biological functions and signaling pathways, including the Focal Adhesion-PI3K-Akt–mTOR-signaling pathway, the neuronal system, metapathway biotransformation Phase I and II and blood circulation, etc. The relationship between the CNVs of five lncRNAs and their gene expression in the ceRNA network was analyzed via a chi-square test, which confirmed that except for LINC00243, the expression of four lncRNAs was notably correlated with the CNVs. The survival analysis revealed that only the copy number gain of LINC00536 was evidently related to the poor prognosis of patients. The CIBERSORT algorithm showed that five lncRNAs were correlated with the abundance of immune cell infiltration and immune checkpoints. In a word, by analyzing CNV-driven lncRNAs and the ceRNA network regulated by these lncRNAs, this study explored the mechanism of breast cancer and provided novel insights into new biomarkers.     

Biomarker expression analysis in different age groups revealed age was a risk factor for breast cancer

The relationship between age and breast cancer is ambiguous. Here, we analyzed the differential expression pattern of long noncoding RNAs (lncRNAs) and messenger RNAs (mRNAs) in different age groups to provide an effective association between age and breast cancer risk at the molecular level. We integrated the microarray information from the Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) data sets. The patients were divided into young ( < 50 years) and old ( ≥ 50 years) age groups and evaluated by differential gene expression, weighted gene correlation network analysis (WGCNA), functional enrichment analyses, and coexpression analysis. To determine their potential clinical significance, univariate Cox regression analysis and survival assessment were conducted. We identified two lncRNAs (AL139280.1 and AP000851.1) and three mRNAs (MT1M, HBB, and TFPI2) as the risk markers, and Gene set enrichment analysis (GSEA) focusing on a single gene revealed that "pyrimidine metabolism," "cell cycle," and "P53 signaling pathway" were coenriched. These data demonstrated that age may be a risk factor for breast carcinogenesis and prognosis and provide an in-depth molecular characterization based on the expression patterns of lncRNAs and mRNAs.     

Genome-wide screening and cohorts validation identifying novel lncRNAs as prognostic biomarkers for clear cell renal cell carcinoma

Clear cell renal cell carcinoma (ccRCC) is the main subtype of renal cell carcinoma with varied prognosis. We aimed to identify and assess the possible prognostic long noncoding RNA (lncRNA) biomarkers. LncRNAs expression data and corresponding clinical information of 619 ccRCC patients were downloaded from The Cancer Genome Atlas (TCGA) and International Cancer Genome Consortium (ICGC) databases. Differentially expressed genes analysis, univariate Cox regression, the least absolute shrinkage and selection operator Cox regression model were utilized to identify hub lncRNAs. Multivariate Cox regression was used to establish the risk model. Statistical analysis was performed using R 3.5.3. The expression value of five lncRNAs and the risk-score levels were significantly associated with a survival prognosis of ccRCC patients (all P < .001). In the TCGA validation cohort, the area under the curve (AUC) for the integrated nomogram was 0.905 and 0.91 for 3-, 5-year prediction separately. The AUC reached up to 0.757 in an independent ICGC cohort. Besides, the calibration plots also illustrated well curve-fitting between observation values and predictive values. Weighted gene co-expression network analysis and subsequent pathway analysis revealed that the PI3K-Akt-mTOR and hypoxia-inducible factor signaling crosstalk might function as the most essential mechanisms related to the five-lncRNAs signature. Our study suggested that lncRNA AC009654.1, AC092490.2, LINC00524, LINC01234, and LINC01885 were significantly associated with ccRCC prognosis. The prognostic model based on this five lncRNA may predict the overall survival of ccRCC.     

Knockdown of long non‐coding RNA LINC00176 suppresses ovarian cancer progression by BCL3‐mediated down‐regulation of ceruloplasmin

Ovarian cancer is a common malignancy among women with some clinically approved diagnostic coding gene biomarkers. However, long non‐coding RNAs (lncRNAs) have been indicated to play an important role in controlling tumorigenesis of ovarian cancer. Hereby, the aim of the study was to uncover the function of lncRNA LINC00176 in the development and progression of ovarian cancer by regulating ceruloplasmin (CP). Bioinformatics prediction in combination with RT‐qPCR analysis for the expression pattern of LINC00176 revealed that LINC00176 was highly expressed in ovarian cancer tissues as well as in ovarian cancer cell lines, respectively. LINC00176 was predominantly localized in the nucleus. Delivery of si‐LINC00176, oe‐LINC00176, si‐BCL3 and si‐CP plasmids was conducted to explore the effects of LINC00176 on ovarian cancer. Promoted proliferation, migration and invasion along with reduced apoptosis were observed in cells treated with oe‐LINC00176, while si‐BCL3 and si‐CP were able to block the promoting effects. Investigations with regard to the correlation between LINC00176 and promoter region of CP turned out to be positive via B‐cell CLL/lymphoma 3 (BCL3) by means of dual‐luciferase reporter gene assay, ChIP and RIP assays. Furthermore, oncogenic properties of the LINC00176/BCL3/CP axis were also demonstrated by tumour formation in vivo generated upon injecting cells in nude mice. Our results demonstrate that restored LINC00176 initiates tumorigenesis in ovarian cancer by increasing CP expression via recruiting BCL3, the mechanism of which represented a potential and promising therapeutic target for the disease.     

Identification of LINC01615 as potential metastasis-related long noncoding RNA in hepatocellular carcinoma

Hepatocellular carcinoma is one of the most prevalent and fatal cancers. Studying the long noncoding RNA (lncRNA) alterations in hepatocellular carcinoma may lead to new therapeutic strategies. We checked whether there were correlations between The Cancer Genome Atlas expression profiles of the differentially expressed lncRNAs and their DNA methylation status or the copy number variations for hepatocellular carcinoma. We obtained 41 lncRNAs that were differentially expressed between tumor and normal samples, and their DNA methylation status was negatively correlated with the expression levels. We identified five lncRNAs that were recurrently amplified or deleted in tumor samples, but none of them were associated with the messenger RNA (mRNA) expression levels. To obtain the biological function of these lncRNAs, the coexpressed mRNAs in the hepatocellular carcinoma were figured out. A total of 10 lncRNAs were highly correlated with at least one gene. Six out of the ten lncRNAs were already known to be related with cancer previously. LINC01615 had 72 coexpressed genes, and we carried out the gene ontology (GO) term enrichment for these protein-coding genes. The results suggested that these lncRNAs were associated with extracellular matrix organization. To summarize, we identified 41 potentially cancer-related lncRNAs. In particular, we proposed that LINC01615 potentially affected the extracellular matrix and had further impacts on the metastasis of hepatocellular carcinoma.     

The joint role of methylation and immune-related lncRNAs in ovarian cancer: Defining molecular subtypes and developing prognostic signature

IntroductionComplex outcome of ovarian cancer (OC) stems from the tumor immune microenvironment (TIME) influenced by genetic and epigenetic factors. This study aimed to comprehensively explored the subclasses of OC through lncRNAs related to both N6-methyladenosine (m6A)/N1-methyladenosine (m1A)/N7-methylguanosine (m7G)/5-methylcytosine (m5C) in terms of epigenetic variability and immune molecules and develop a new set of risk predictive systems.Material and methodsThe lncRNA data of OC were collected from TCGA. Spearman correlation analysis on lncRNA data of OC with immune-related gene expression and with m6A/m5C/m1A/m7G were respectively conducted. The m6A/m5C/m1A/m7G-related m6A/m5C/m1A/m7G related immune lncRNA subtypes were identified on the basis of the prognostic lncRNAs. Heterogeneity among subtypes was evaluated by tumor mutation analysis, tumor microenvironment (TME) component analysis, response to immune checkpoint blocked (ICB) and chemotherapeutic drugs. A risk predictive system was developed based on the results of Cox regression analysis and random survival forest analysis of the differences between each specific cluster and other clusters.ResultsThree m6A/m5C/m1A/m7G-related immune lncRNA subtypes of OC showing distinct differences in prognosis, mutation pattern, TIME components, immunotherapy and chemotherapy response were identified. A set of risk predictive system consisting of 10 lncRNA for OC was developed, according to which the risk score of samples in each OC dataset was calculated and risk type was defined.ConclusionsThis study classified three m6A/m5C/m1A/m7G-related immune lncRNA subtypes with distinct heterogeneous mutation patterns, TME components, ICB therapy and immune response, and provided a set of risk predictive system consisted of 10 lncRNA for OC.     

Identification of a glycolysis-related lncRNA prognostic signature for clear cell renal cell carcinoma

Background: The present study investigated the independent prognostic value of glycolysis-related long noncoding (lnc)RNAs in clear cell renal cell carcinoma (ccRCC).Methods: A coexpression analysis of glycolysis-related mRNAs–long noncoding RNAs (lncRNAs) in ccRCC from The Cancer Genome Atlas (TCGA) was carried out. Clinical samples were randomly divided into training and validation sets. Univariate Cox regression and least absolute shrinkage and selection operator (LASSO) regression analyses were performed to establish a glycolysis risk model with prognostic value for ccRCC, which was validated in the training and validation sets and in the whole cohort by Kaplan–Meier, univariate and multivariate Cox regression, and receiver operating characteristic (ROC) curve analyses. Principal component analysis (PCA) and functional annotation by gene set enrichment analysis (GSEA) were performed to evaluate the risk model.Results: We identified 297 glycolysis-associated lncRNAs in ccRCC; of these, 7 were found to have prognostic value in ccRCC patients by Kaplan–Meier, univariate and multivariate Cox regression, and ROC curve analyses. The results of the GSEA suggested a close association between the 7-lncRNA signature and glycolysis-related biological processes and pathways.Conclusion: The seven identified glycolysis-related lncRNAs constitute an lncRNA signature with prognostic value for ccRCC and provide potential therapeutic targets for the treatment of ccRCC patients.     

Five key lncRNAs considered as prognostic targets for predicting pancreatic ductal adenocarcinoma

Pancreatic ductal adenocarcinoma (PDAC) has a poor prognosis, and the 5‐year survival rate was only 7.7%. To improve prognosis, a screening biomarker for early diagnosis of pancreatic cancer is in urgent need. Long non‐coding RNA (lncRNA) expression profiles as potential cancer prognostic biomarkers play critical roles in development of tumorigenesis and metastasis of cancer. However, lncRNA signatures in predicting the survival of a patient with PDAC remain unknown. In the current study, we try to identify potential lncRNA biomarkers and their prognostic values in PDAC. LncRNAs expression profiles and corresponding clinical information for 182 cases with PDAC were acquired from The Cancer Genome Atlas (TCGA). A total of 14 470 lncRNA were identified in the cohort, and 175 PDAC patients had clinical variables. We obtained 108 differential expressed lncRNA via R packages. Univariate and multivariate Cox proportional hazards regression, lasso regression was performed to screen the potential prognostic lncRNA. Five lncRNAs have been recognized to significantly correlate with OS. We established a linear prognostic model of five lncRNA (C9orf139, MIR600HG, RP5‐965G21.4, RP11‐436K8.1, and CTC‐327F10.4) and divided patients into high‐ and low‐risk group according to the prognostic index. The five lncRNAs played independent prognostic biomarkers of OS of PDAC patients and the AUC of the ROC curve for the five lncRNAs signatures prediction 5‐year survival was 0.742. In addition, targeted genes of MIR600HG, C9orf139, and CTC‐327F10.4 were explored and functional enrichment was also conducted. These results suggested that this five‐lncRNAs signature could act as potential prognostic biomarkers in the prediction of PDAC patient's survival.     

Downregulation of long noncoding RNA LINC00683 associated with unfavorable prognosis in prostate cancer based on TCGA

Prostate cancer (PCa) is the third most common reason of cancer-related deaths in men. Accumulating evidence has shown that dysregulation of long noncoding RNAs (lncRNAs) is closely related to cancer initiation and development. Although large numbers of lncRNAs have been discovered, knowledge regarding their function and physiological/pathological significance remains limited. In this study, we aimed to reveal functional lncRNAs and identify prognosis-related RNAs in PCa by analyzing data from The Cancer Genome Atlas (TCGA). To achieve this, an lncRNA-mRNA coexpression network was constructed by weighted correlation network analysis. Additionally, a subnetwork was extracted from this weighted correlation network, and seven lncRNAs were identified as core nodes. Further Kaplan-Meier survival analysis showed that three lncRNAs (LINC00683, LINC00857, and FENDRR) were significantly downregulated in PCa samples, and there was a strong positive correlation with patient survival. Importantly, LINC00683 has not been fully reported as related with PCa. Additionally, gene set enrichment analysis indicated that LINC00683 might be involved in cancer-related pathways such as the Wnt pathway. Based on the findings of this study, lncRNA LINC00683 is likely to provide a new diagnostic biomarker and therapeutic target for future PCa treatments.     

Development of a novel six-long noncoding RNA signature predicting survival of patients with bladder urothelial carcinoma

Bladder urothelial carcinoma is a malignant tumor with a high incidence in the uropoietic system. Considerable studies have shown that long noncoding RNA (lncRNA) plays an important role in the development and progression of bladder urothelial carcinoma. In this study, the lncRNA expression and clinical data of 377 bladder urothelial carcinoma patients were obtained from The Cancer Genome Atlas database and differentially expressed lncRNAs in cancer and normal groups were evaluated. Univariate COX and multivariate COX regression analyses of prognosis were performed on differentially expressed lncRNAs in the training data sets, six prognosis-related lncRNAs (LINC02195, LINC01484, LINC01468, SMC2-AS1, AC011298.1, and PTPRD-AS1) were assessed, and a six-lncRNA signature was constructed. The predictive capability of this six-lncRNA signature was validated in the testing data sets and entire data sets. The prognostic ability of the six-lncRNA signature was independent of other clinical elements after multivariate COX regression and stratified analyses of with other clinical elements. We performed functional enrichment analysis with the six prognosis-related lncRNAs. Results of functional enrichment revealed that these prognosis-related lncRNAs might promote the development and metastasis of bladder urothelial carcinoma. In summary, the six-lncRNA signature that we developed could effectively predict the prognosis of bladder urothelial carcinoma patients. This six-lncRNA signature might be a novel independent prognostic marker of bladder urothelial carcinoma. Moreover, it also provides novel insights into the mechanism of bladder urothelial carcinoma.     

Prognostic power of a lipid metabolism gene panel for diffuse gliomas

Lipid metabolism reprogramming plays important role in cell growth, proliferation, angiogenesis and invasion in cancers. However, the diverse lipid metabolism programmes and prognostic value during glioma progression remain unclear. Here, the lipid metabolism‐related genes were profiled using RNA sequencing data from The Cancer Genome Atlas (TCGA) and Chinese Glioma Genome Atlas (CGGA) database. Gene ontology (GO) and gene set enrichment analysis (GSEA) found that glioblastoma (GBM) mainly exhibited enrichment of glycosphingolipid metabolic progress, whereas lower grade gliomas (LGGs) showed enrichment of phosphatidylinositol metabolic progress. According to the differential genes of lipid metabolism between LGG and GBM, we developed a nine‐gene set using Cox proportional hazards model with elastic net penalty, and the CGGA cohort was used for validation data set. Survival analysis revealed that the obtained gene set could differentiate the outcome of low‐ and high‐risk patients in both cohorts. Meanwhile, multivariate Cox regression analysis indicated that this signature was a significantly independent prognostic factor in diffuse gliomas. Gene ontology and GSEA showed that high‐risk cases were associated with phenotypes of cell division and immune response. Collectively, our findings provided a new sight on lipid metabolism in diffuse gliomas.     

Identification of lncRNA–miRNA–mRNA regulatory network associated with epithelial ovarian cancer cisplatin-resistant

To construct a long noncoding RNA (lncRNA)–microRNA (miRNA)–messenger RNA (mRNA) regulatory network related to epithelial ovarian cancer (EOC) cisplatin-resistant, differentially expressed genes (DEGs), differentially expressed lncRNAs (DELs), and differentially expressed miRNAs (DEMs) between MDAH and TOV-112D cells lines were identified. Gene Ontology (GO) analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis were conducted to analyze the biological functions of DEGs. Downstream mRNAs or upstream lncRNAs for miRNAs were analyzed at miRTarBase 7.0 or DIANA-LncBase V2, respectively. A total of 485 significant DEGs, 85 DELs, and 5 DEMs were identified. Protein–protein interaction (PPI) network of DEGs contrains 81 nodes and 141 edges was constructed, and 25 hub genes related to EOC cisplatin-resistant were identified. Subsequently, a lncRNA–miRNA–mRNA regulatory network contains 4 lncRNAs, 4 miRNAs, and 35 mRNAs was established. Taken together, our study provided evidence concerning the alteration genes involved in EOC cisplatin-resistant, which will help to unravel the mechanisms underlying drug resistant.     

Identification of two molecular subtypes of dysregulated immune lncRNAs in ovarian cancer

Long non-coding RNA (lncRNA) has increasingly been identified as a key regulator in pathologies such as cancer. Multiple platforms were used for comprehensive analysis of ovarian cancer to identify molecular subgroups. However, lncRNA and its role in mapping the ovarian cancer subpopulation are still largely unknown. RNA-sequencing and clinical characteristics of ovarian cancer were acquired from The Cancer Genome Atlas database (TCGA). A total of 52 lncRNAs were identified as aberrant immune lncRNAs specific to ovarian cancer. We redefined two different molecular subtypes, C1(188) and C2(184 samples), in “iClusterPlus” R package, among which C2 grouped ovarian cancer samples have higher survival probability and longer median survival time (P <0.05) with activated IFN-gamma response, Wound Healing and Cytotoxic lymphocytes signal; 456 differentially expressed genes were acquired in C1 and C2 subtypes using limma (3.40.6) package, among which 419 were up-regulated and 37 were down-regulated, in TCGA dataset. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes (KEGG) functional enrichment analysis revealed that these genes were actively involved in ECM-receptor interaction, PI3K-Akt signaling pathway interaction KEGG pathway. Compared with the existing immune subtype, the Cluster2 sample showed a substantial increase in the proportion of the existing C2 immune subtype, accounting for 81.37%, which was associated with good prognosis. Our C1 subtype contains only 56.49% of the existing immune C1 and C4, which also explains the poor prognosis of C1. Furthermore, 52 immune-related lncRNAs were used to divide the TCGA-endometrial cancer and cervical cancer samples into two categories, and C2 had a good prognosis. The differentially expressed genes were highly correlated with immune-cell-related pathways. Based on lncRNA, two molecular subtypes of ovarian cancer were identified and had significant prognostic differences and immunological characteristics.     

LINC01101 and LINC00277 expression levels as novel factors in HPV‐induced cervical neoplasia

Recently long non‐coding RNAs were identified as new factors involved in gene expression regulation. To gain insight into expression pattern of these factors related to E7 HPV18 oncogene, this study uses HeLa cell culture transfected with E7‐siRNA. Gene expression profile was investigated using microarray analysis. After analysing the microarray results, we identified 15,387 RNA species differentially expressed in E7‐siRNA‐transfected cells compared with controls (fold change >2). The expression profiles of lncRNA species highlighted 731 lncRNAs and 203 lincRNAs. We selected two lincRNAs (LINC01101 and LINC00277) and we evaluated the expression profile in HPV‐induced neoplasia. Both lincRNAs investigated display a significantly reduced pattern of expression in cervical lesions and cancer, associated with clinical parameters. A connection between HPV presence and lincRNAs was noted. hrHPV‐positive samples exhibit significantly reduced LINC01101 and LINC00277 expression level (P < 0.05). These results provide new insights into involvement of lncRNA in HPV‐induced cervical cancer, enriching our understanding of their potential role in this pathology.