Original Research ARTICLE
The Identification and Analysis of mRNA–lncRNA–miRNA Cliques From the Integrative Network of Ovarian Cancer
- 1Department of Tumor Biological Treatment, The Third Affiliated Hospital of Soochow University, Changzhou, China
- 2Jiangsu Engineering Research Center for Tumor Immunotherapy, Changzhou, China
- 3Institute of Cell Therapy, Soochow University, Changzhou, China
- 4Shanghai Institute of Nutrition and Health, Shanghai Institutes for Biological Sciences (CAS), Shanghai, China
Ovarian cancer is one of the leading causes of cancer mortality in women. Since little clinical symptoms were shown in the early period of ovarian cancer, most patients were found in phases III–IV or with abdominal metastasis when diagnosed. The lack of effective early diagnosis biomarkers makes ovarian cancer difficult to screen. However, in essence, the fundamental problem is we know very little about the regulatory mechanisms during tumorigenesis of ovarian cancer. There are emerging regulatory factors, such as long noncoding RNAs (lncRNAs) and microRNAs (miRNAs), which have played important roles in cancers. Therefore, we analyzed the RNA-seq profiles of 407 ovarian cancer patients. An integrative network of 20,424 coding RNAs (mRNAs), 10,412 lncRNAs, and 742 miRNAs were construed with variance inflation factor (VIF) regression method. The mRNA–lncRNA–miRNA cliques were identified from the network and analyzed. Such promising cliques showed significant correlations with survival and stage of ovarian cancer and characterized the complex sponge regulatory mechanism, suggesting their contributions to tumorigenicity. Our results provided novel insights of the regulatory mechanisms among mRNAs, lncRNAs, and miRNAs and highlighted several promising regulators for ovarian cancer detection and treatment.
Ovarian cancer (OC) ranks the fifth leading cause of cancer mortality among gynecologic malignancies in women (Torre et al., 2018). According to the SEER registry, 22,240 new cases and 14,070 deaths of ovarian cancer have been estimated in the United States in 2018. Ovarian cancer, encompassing various tumors of ovarian origin, is a typical example of heterogeneous disease and classified into three categories based on the affected cells: epithelial cells, germ cells, and stromal cells, of which epithelial ovarian cancer (EOC) is the most common cause of death within gynecologic cancer (Vetter and Hays, 2018). Due to the insidious nature of ovarian cancer, manifesting with little to no clinical symptoms in the early period, 80% of patients were found in stage III or IV, which explains the low 5-year survival rate and poses significant therapeutic challenges (Torre et al., 2018). Although great progress has been made in screening methods such as transvaginal ultrasonography as well as serum biomarker CA125 and HE4 for early detection of ovarian cancer, diagnostic and treatment limitations still exist because of low sensitivity, high expenses, and/or inconvenience. Thus, it is imperative to identify new risk genes and their regulatory network in order to elucidate ovarian carcinogenetic mechanisms through genetic research.
A decade-long interest in the cross talk between coding RNAs (mRNAs) and noncoding RNAs has unveiled multilayer regulatory circuitry of cells including chromatin remodeling, protein stability, transcription, and mRNA turnover; the dysregulation of which could also lead to the occurrence of cancer. Increasing studies have demonstrated that microRNAs (miRNAs), which are small noncoding RNAs (∼22 nts), and long noncoding RNAs (lncRNAs), which are usually longer than 200 nucleotides, mutually regulated their expression levels and worked jointly to control the expression of mRNAs (Braconi et al., 2011; Fan et al., 2013; Martens-Uzunova et al., 2014; Liu et al., 2017a). Conversely, mRNAs also impose their influences on the expression of noncoding RNAs in a variety of ways (Helwak et al., 2013; Ribeiro et al., 2014; Sundaram and Veera Bramhachari, 2017). For example, a muscle-specific lncRNA linc-MD1 elevated the expression of MAML1 and MEF2C mRNAs by preventing miR-133 and miR-135 from binding to their targets (Cesana et al., 2011; Legnini et al., 2014). Linc-ROR has been reported to function as a key competing endogenous RNA to partially sequester miR-145-5p, which was activated by Nanog, Oct4, and Sox2, forming a feedback loop with core transcription factors and miRNAs to regulate embryonic stem cell maintenance and differentiation (Wang et al., 2013). The widely studied lncRNA HOTAIR acted as a scaffold to link several RNA-binding proteins including Ataxin-1 and Snurportin-1, whose stability was decreased by miR-34a in the human prostate cancer cells (Chiyomaru et al., 2013; Yoon et al., 2013). Such interactions are called mRNA–lncRNA–miRNA triplets; the aberration of which could destroy the gene expression patterns and promote tumorigenesis (Gao et al., 2015; Imig et al., 2015; Wang et al., 2015a).
With the development of Next Generation Sequencing (NGS), all kinds of RNAs, such as lncRNAs, miRNAs, mRNAs, Piwi-interacting RNAs (piRNAs), and circular RNAs (circRNAs), can be easily quantified. The complex regulatory mechanism among them can be investigated with advanced network analysis. For ovarian cancer, The Cancer Genome Atlas (TCGA) (The Cancer Genome Atlas Research Network et al., 2013) is a valuable dataset with the mRNA, lncRNA, and miRNA expression profiles of 407 patients. To construct the mRNA–lncRNA–miRNA regulatory network and the cancer-related mRNA, lncRNA, and miRNA triplets, we adopted variance inﬂation factor (VIF) regression method (Liu et al., 2017a). A comprehensive genome-wide mRNA, lncRNA, and miRNA regulatory network of ovarian cancer was built, and the substructures of mRNA–lncRNA–miRNA cliques were discovered on the network. These triplets were further filtered by associations with survival and ovarian cancer stage. At last, nine triplets were analyzed in detail. The mechanisms revealed by this work may help understand the tumorigenesis and may develop new treatment for ovarian cancer.
Ovarian Cancer RNA-Sequencing Datasets
The RNA-sequencing data of ovarian cancer are downloaded from Akrami et al. (2013) at http://larssonlab.org/tcga-lncrnas/datasets.php. It included the expression levels of 20,462 protein coding mRNAs, 10,419 lncRNAs, and 742 miRNAs in 407 ovarian cancer patients. The summary of this ovarian cancer dataset is shown in Table 1. As described by Akrami et al. (2013), the expression levels of these transcripts were quantified as RPKM (reads per kilobase per million mapped reads) by normalizing the mRNA length and library size after the read counts were generated using TopHat (Kim et al., 2013) and HTSeq-count (Anders et al., 2015). To enable the log2 transformation, the zero values of the expression data were replaced with the minimum nonzero RPKM values. The TCGA sample IDs and their clinical information were listed in Supplementary Material 1. Within the 407 samples, four samples did not have stage information, 20 samples were in stage II, 323 samples were in stage III, and 60 samples were in stage IV. The sample distribution was consistent with the incidences of ovarian cancer: most of them were in stages III and IV (Torre et al., 2018).
Table 1 Number of transcripts in expression data, constructed network, triplets of the 407 ovarian cancer patients.
The mRNA–lncRNA–miRNA Network Construction Based on Variance Inflation Factor Regression
Previous studies have shown that miRNAs and lncRNAs can interact with each other and regulate the expression of mRNAs jointly (Braconi et al., 2011; Fan et al., 2013; Martens-Uzunova et al., 2014; Liu et al., 2017a). To get a genome-wide view of such complex regulations among miRNAs, lncRNAs, and mRNAs, we adopted a method called variance inflation factor (VIF) regression (Lin et al., 2011) to construct the mRNA–lncRNA–miRNA network of ovarian cancer.
Since there were 31,623 transcripts in total, including 20,462 protein coding mRNAs, 10,419 lncRNAs, and 742 miRNAs, for each regression model that corresponded to a transcript, there will be 31,622 variables that may regulate it. It will be a very large regression model. The computational complexity is extremely high. Fortunately, the VIF algorithm is a highly efficient method and can select the possible regulators of the target RNA in a short time.
The optimization goal of the VIF algorithm is to estimate β that minimizes l0 norm (Lin et al., 2011), the penalized sum of squared errors between actual expression levels of the target transcript and predicted expression levels,
where y are actual expression levels of the target transcript in 407 ovarian cancer patients, x are the expression levels of the 31,622 candidate regulators, β is the vector of coefficient parameters, λ0 is the penalty, p is the number of candidate regulators, i.e., 31,622 in this study, is the number of coefficients that are not zero.
If we use the conventional method, there will be 2p possible combinations of candidate regulators to search for the best β. It is an innumerous space to search. The VIF algorithm takes a greed strategy. First, it evaluates the marginal correlation of each candidate regulator with the target transcript in a pre-sampled small dataset. Then, it searches the optimal candidate regulator subset with t-statistic correction by adding/removing one candidate regulator each time. The algorithm is not only efficient but also accurate. Its accuracy is comparable with the highly recognized regression method, LASSO. We used the R package VIF (Lin et al., 2011) to apply the VIF regression algorithm.
After the regression model was constructed, we used the adjust coefficient of determination R2 to evaluate the goodness-of-fit. It measured how well the predicted value approximates the actual values with adjustment of how complex the regression model was, i.e. how many regulators were included in the model. We only kept the regression models with the adjust R2 greater than 0.8. The target transcript and all the regulators in the optimal subset of each regression model constructed a regulatory association network with edge directions that the regulators regulated the target. Then, all the association networks constructed from selected regression models are combined to a single network.
The Identification and Analysis of mRNA–lncRNA–miRNA Cliques From the Network
As we mentioned before, the mRNA–lncRNA–miRNA triplet functions together and forms a clique on network. A mRNA–lncRNA–miRNA clique includes three fully connected nodes: one mRNA, one lncRNA, and one miRNA. Such structure has significant biological meanings and worth to be identified and analyzed. We used R package igraph (Gabor Csardi, 2006) to find the cliques of mRNA–lncRNA–miRNA triplets on the integrative network constructed with VIF regression. In a clique of mRNA–lncRNA–miRNA, there were regulations between lncRNA and mRNA, miRNA and mRNA, and lncRNA and miRNA. Since the transcriptome data were not time course, it was difficult to decide the direction of regulation accurately. We treated the VIF integrative network as an association network rather than a directional network.
The Survival Analysis of the RNAs Within the mRNA–lncRNA–miRNA Cliques
There were 406 ovarian cancer patients with both the overall survival time (months) and the survival status. We performed survival analysis using Cox proportional hazards regression model (Andersen and Gill, 1982) on these samples. For each RNA, the patients were divided into two groups: the patients with expression levels smaller than the median and the patients with expression levels greater than or equal to the median. The survival curves of these two groups of patients were plotted as Kaplan–Meier plot. The significance of survival difference between these two patient groups was evaluated with log rank test p value calculated by univariate survival analysis. For each triplet, the log rank p value was calculated using multivariate survival analysis of the mRNA, lncRNA, and miRNA. If the log rank test p value was smaller than 0.05, their survivals were considered as significantly different. The R package survival (https://CRAN.R-project.org/package=survival) was used to perform the survival analysis.
Results and Discussion
Construction of the Ovarian Cancer Integrative Network Based on VIF Regression
Previous studies have built various computational models to predict potential disease-related lncRNA/miRNA (Chen and Yan, 2013; Chen and Huang, 2017; Chen et al., 2017; Chen et al., 2018a; Chen et al., 2018b; Chen et al., 2018c) that showed high efficiency, accuracy, and stability. In this study, we firstly introduced cliques incorporating mRNA, lncRNA, and miRNA and constructed their genome-wide integrative network based on VIF regression (Figure 1). First, for each target transcript, a regression model was estimated using all the other 31,622 transcripts. The candidate regulators could be mRNA, lncRNA, or miRNA. Then, we filtered the 31,623 VIF regression models with cutoff of adjusted R2 greater than 0.8. Finally, we combined the filtered highly possible regression models to obtain a complex integrative network that was given in Supplementary Material 2. It had 381,085 regulations between 21,670 transcripts that included 16,667 coding mRNAs, 4,796 lncRNAs, and 207 miRNAs as shown in Table 1. It can be seen that about 81.5% measured coding mRNAs, 46.0% measured lncRNAs, and 27.9% measured miRNAs played roles in the integrative network. The lncRNAs may have played more functional regulation roles than miRNAs.
Identification of mRNA–lncRNA–miRNA Cliques Based on Network Structure Analysis
We analyzed the graph structure of ovarian cancer integrative network using R package igraph (Gabor Csardi, 2006) and identified 7,311 mRNA–lncRNA–miRNA cliques. Each clique included one mRNA, one lncRNA, and one miRNA. What is more, there were three interactions between mRNA and lncRNA, mRNA and miRNA, and miRNA and lncRNA. These 7,311 mRNA–lncRNA–miRNA triplets were given in Supplementary Material 3. To explore the biological significances of these triplets, we did survival analysis of each RNA and calculated the log rank test p value. If a RNA’s log rank test p value was smaller than 0.05, it was listed as 1, otherwise as 0 in Supplementary Material 3. What is more, we performed the analysis of variance (ANOVA) among different cancer stages and calculated the ANOVA p values. If a RNA’s ANOVA p value was smaller than 0.05, it was listed as 1, otherwise as 0 in Supplementary Material 3. For univariate survival analysis of RNAs, there were nine mRNAs and 15 miRNAs that showed significant associations with survival but no survival related lncRNAs. For multivariate survival analysis of triplets, there were 786 triplets with log rank p values smaller than 0.05. For stage analysis, there were 637 mRNAs, 464 lncRNA, and 316 miRNAs that showed significant associations with stage.
As shown in Table 1, within the 7,311 mRNA–lncRNA–miRNA triplets, there were 633 coding mRNAs, 1,184 lncRNAs, and 169 miRNAs. It can be seen that about 3.80% coding mRNAs, 24.7% lncRNAs, and 81.6% miRNAs on the integrative network function through the clique structures. The miRNAs relied more on the cliques than lncRNAs. The lncRNAs may have other regulatory mechanisms, while most miRNAs were involved in sponge mechanism (Furió-Tarí et al., 2016). For mRNAs, there are many other mechanisms, such transcription factors, alternative splicing, and A-to-I RNA-editing. Furthermore, they may be regulated by lncRNAs or miRNAs alone without the complete clique structure.
By categorizing these mRNA–lncRNA–miRNA triplets based on survival and stage association, we plotted nine triplets that showed significant association with survival and cancer stage in Figure 2. The Kaplan–Meier plots of the four RNAs associated with survival, MIR600HG, MIR4519, POLD3, and MTA1, were given in Figure 3. The high expression of MIR600HG and the low expression of MIR4519, POLD3, and MTA1 were associated with high risk.
Figure 2 Key mRNA–lncRNA–miRNA triplets on the integrative network of ovarian cancer. The red, green, and yellow nodes represented mRNA, lncRNA, and miRNA, respectively. Each mRNA–lncRNA–miRNA triplet formed a clique on the integrative network. The rectangle, triangle, and hexagon represented survival significant, cancer-stage significant, and both significant, respectively. In a clique of mRNA–lncRNA–miRNA, there were regulations between lncRNA and mRNA, miRNA and mRNA, and lncRNA and miRNA. The numbers on the edges were adjusted R2. (A) The mRNA–lncRNA–miRNA triplets of NFYA, RP11-401P9.4.1 and MIR600HG. (B) The mRNA–lncRNA–miRNA triplets of RBM5, RP11-401P9.4.1 and MIR600HG. (C) The mRNA–lncRNA–miRNA triplets of ALS2, RP11-401P9.4.1 and MIR600HG. (D) The mRNA–lncRNA–miRNA triplets of LAP3, RP11-401P9.4.1 and MIR600HG. (E) The mRNA–lncRNA–miRNA triplets of RAD52, RP11-401P9.4.1 and MIR600HG. (F) The mRNA–lncRNA–miRNA triplets of LUC7L, RP11-169K16.9.1 and MIR4519. (G) The mRNA–lncRNA–miRNA triplets of POLD3, RP11-736I10.1.1 and MIR130B. (H) The mRNA–lncRNA–miRNA triplets of POLD3, CTC-471C19.2.1 and MIR130B. (I) The mRNA–lncRNA–miRNA triplets of MTA1, RP11-492L8.2.1 and MIR218-1.
Figure 3 Kaplan–Meier plots of the four RNAs associated with survival. The Kaplan–Meier plots of MIR600HG, MIR4519, POLD3, and MTA1 were shown in A, B, C, and D, respectively. (A) The high expression of MIR600HG was associated with high risk. (B) The low expression of MIR4519 was associated with high risk. (C) The low expression of POLD3 was associated with high risk. (D) The low expression of MTA1 was associated with high risk.
Functional Analysis of Key mRNA–lncRNA–miRNA Triplets on the Ovarian Cancer Integrative Network
First, we focus on the mRNA–lncRNA–miRNA interaction and biological functions in Figures 2A–E, as they share the same lncRNA and miRNA. RP11-401P9.4.1, also known as lnc-NKD1-1:1, has been reported to be up regulated in multiple cancers including hepatocellular carcinoma (HCC) and clear cell renal cell carcinoma (ccRCC) (Ellinger et al., 2015; Zhang et al., 2015). Evolutionary clues in lncRNAs, especially structure, function, and expression not canonical sequences, have been validated among lncRNA homologs in vertebrate (Li et al., 2013; Xu et al., 2016; Nitsche and Stadler, 2017; Tu et al., 2018; Fico et al., 2019), underscoring the importance of conservative function in lncRNAs. In silico functional analysis has shown that lnc-NKD1-1:1 regulated the genes involved in mitosis and cell cycle, implying its important roles in the progression of cancer. Hence, it is reasonable to infer that lnc-NKD1-1:1 performs similar functions in ovarian cancer. MIR600HG has been proved as a prognostic biomarker in predicting survival of pancreatic ductal adenocarcinoma (PDAC) patients (Song et al., 2018), which is also consistent with our results (Supplementary Material 3). In addition, its function could be inferred by its significantly co-expressed gene scaffold adapter GRB2-associated binding protein 2 (GAB2), whose overexpression has been suggested to promote tumor angiogenesis and growth in ovarian cancer cells (Duckworth et al., 2016). Therefore, lnc-NKD1-1:1 may have its specific way to interact with MIR600HG for their similar contributions to promote the malignant process of cancer.
Next, we consider the five mRNAs in the triplets. NFYA, the regulatory subunit of the CCAAT binding transcription factor NFY, has been observed with elevated levels in EOC cells compared with ovarian surface epithelial cells (Garipov et al., 2013), knockdown of which suppressed growth and induced apoptosis of human EOC cells, emphasizing its requirement for ovarian cancer cell proliferation. As a component of pre-spliceosomal complexes, the RNA binding motif protein 5 (RBM5) has been demonstrated to regulate the alternative splicing process of several mRNAs (Nguyen et al., 2011). The upregulation of RBM5 has been found in ovarian and breast cancers (Farina et al., 2011), while decreased levels of RBM5 upon RAS activation were detected in various solid cancer types and exhibited its association with metastasis (Kim et al., 2010; Farina et al., 2011), suggesting its crucial roles in tumorigenesis. ALS2 encodes a ubiquitous protein Alsin, which is mainly expressed in the central nervous system (Gautam et al., 2016). It is important for maintaining cellular integrity, as mutations in ALS2 can lead to many motor neuron diseases (Yang et al., 2001). Although there is little literature about the connections between ALS2 and ovarian cancer, our results may provide new insights into the function of this gene, and further mechanisms need to be investigated. As a cell surface aminopeptidase, leucine aminopeptidase 3 (LAP3) is responsible for catalyzing the hydrolysis of leucine residues. Accumulating evidence has revealed that overexpression of LAP3 correlated with prognosis of several cancers, such as esophageal squamous cell carcinoma (ESCC), glioma, and HCC (He et al., 2015). Meanwhile, a recent study on ovarian cancer cells has indicated that suppression of LAP3 inhibited cell migration and invasion through its target genes fascin and MMP-2/9 (Wang et al., 2015b), underlying its important roles in regulating cancer cell metastasis. Previous studies have presented the involvement of radiation repair protein 52 (RAD52) in the cell cycle control as well as DNA repair, and depletion of RAD52 could cause synthetic lethality in BRCA1 mutant breast cancer cells (Hromas et al., 2017). Additionally, several groups have conferred RAD52 as novel risk loci in ovarian cancer patients (Stafford et al., 2017; Zhao et al., 2017), underlining its significant effects on controlling cancer cell activities. In conclusion, all elements in these five triplets contribute cohesively to the development of ovarian cancer through regulating cell cycle and growth.
Then, we explore the triplet mRNA LUC7L–lncRNA RP11-169K16.9.1–miRNA MIR4519 shown in Figure 2F. Apart from functioning in myogenesis, putative RNA-binding protein Luc7-like 1 (LUC7L) has recently been regarded as a predictor of breast cancer survival through transcriptional network construction and showed its correlation with metastasis (Crawford et al., 2008). MIR4519 contributes to tumorigenicity by interacting with its direct target N-ethylmaleimide-sensitive factor (NSF) and associated receptor SNARE. The latter have been demonstrated to be involved in cell migration and invasion by regulating vesicle trafficking and membrane fusion (Meng and Wang, 2015; Sun et al., 2016; Naydenov et al., 2018; Yoon and Munson, 2018). Furthermore, SNARE complex have been reported to affect cell apoptosis and proliferation, highlighting their great impacts on facilitating ovarian cancer cell viability (Sun et al., 2016). RP11-169K16.9.1, also called lnc-UQCRHL-1:1, is an antisense lncRNA derived from UQCRHL gene, whose genetic variants have been identified in breast cancer samples with specific DNA amplification region (Parris et al., 2018). Also, prior research has considered UQCRHL as a prognostic factor in HCC for its pivotal roles in mitochondrial respiration (Park et al., 2017). Thus, it can be inferred that lnc-UQCRHL-1:1 may exert key influences on cellular respiration, aberration of which is frequently found in tumor progression (Vyas et al., 2016). This triplet concentrates on the migration process that triggers ovarian cancer metastasis.
Subsequently, we analyze two triplets in Figures 2G, H, as they share the same mRNA and miRNA. Acting as a component of the DNA polymerase delta complex, POLD3 has been served as a prognostic biomarker in ovarian cancer (Willis et al., 2016), which is in accordance with our results (Supplementary Material 3) and indicates its possible roles in clinical application. Increasing evidence has shown that MIR130B participated in multidrug resistance in ovarian cancer by targeting NRP1 and CSF1, which influenced cell motility and adhesion (Yang et al., 2012; Chen et al., 2016). Furthermore, abnormal methylation levels of MIR130B were found in drug-resistant cell lines as well as ovarian cancer tissues (Yang et al., 2012). Although RP11-736I10.1.1 is a rare lncRNA with little literature report, its function could be speculated by its neighboring gene PLAC1, which is an X-linked trophoblast gene and has recently been identified as a serum biomarker for breast cancer (Baldwin et al., 2008; Yuan et al., 2018), suggesting its potential relationship with tumor initiation and progression. CTC-471C19.2.1, also named lnc-UBE2QL1-1:1, is originated from UBE2QL1, which has proved to interconnect with other genes enriched in cell cycle and metabolic pathways in HCC and renal cancer (Wake et al., 2013; Ye et al., 2016). Considering the aforementioned interspecies conservation in lncRNA functionality (Nitsche and Stadler, 2017; Fico et al., 2019), our results expand new links between this lncRNA and ovarian cancer. In summary, above coding and noncoding RNAs may serve as candidate genes that need more attention considering drug resistance in ovarian cancer.
The last triplet is mRNA MTA1–lncRNA RP11-492L8.2.1–miRNA MIR218-1 (Figure 2I). MTA1 is a metastasis-associated gene and has been observed to stimulate proliferation of EOC cells by virtue of enhancing DNA repair (Yang et al., 2014). Previous studies have revealed that MIR218-1 affected tumor angiogenesis and further has specific influences on cell migration and invasion in several cancers including cervical squamous cell carcinoma and prostate cancer (Yamamoto et al., 2013; Guan et al., 2016), stressing its great effects on determining cancer cell fate. RP11-492L8.2.1, also annotated as lnc-DMXL1-7, may interact with DMXL1, whose upregulation was found in leukemic cells (Koldehoff et al., 2008). Given these interactions, these three RNAs cluster together for their similar efforts to tumor initiation and progression mediated by regulating angiogenesis, proliferation, and migration of malignant cells.
Identification of Four Cliques Based on Network Structure Analysis
We have analyzed the mRNA–lncRNA–miRNA triplets, which were actually three cliques. If we increase the size of cliques to four, what will happen? To investigate the effects of clique size, we extracted the four cliques from the integrative network. There were 11,117 four cliques with at least one mRNA, lncRNA, and miRNA. The number of four cliques was greater than the number of triplets, which was 7,311. It suggested that some four cliques shared the same triplets. The 11,117 four cliques containing mRNA–lncRNA–miRNA triplet were given in Supplementary Material 4. To explore the biological significances of these four cliques, we also did survival analysis and ANOVA for each RNA. If a RNA’s log rank test p value was smaller than 0.05, it was listed as 1, otherwise as 0 in Supplementary Material 4. If an RNA’s ANOVA p value among different cancer stages was smaller than 0.05, it was listed as 1, otherwise as 0 in Supplementary Material 4.
There were 15 four cliques that were both associated with survival and cancer stages (Figure 4). The survival-associated RNAs in these four cliques were MIR600HG, POLD3, and CLPTM1L. MIR600HG and POLD3 have already been discussed. The Kaplan–Meier plots of CLPTM1L were given in Figure 5. The high expression of CLPTM1L was associated with high risk.
Figure 4 Four-cliques that were both associated with survival and cancer stage. The red, green, and yellow nodes represented mRNA, lncRNA, and miRNA, respectively. The rectangle, triangle, and hexagon represented survival significant, cancer-stage significant, and both significant, respectively. (A) Four-cliques of CTD-2142D14.1.1, MIR124-3, CLPTM1L and NSUN2. (B) Four-cliques of RP11-318I4.1.1, MIR124-3, CLPTM1L and NSUN2. (C) Four-cliques of CTC-471C19.2.1, MIR130B, POLD3 and FAM168A. (D) Four-cliques of RP11-736I10.1.1, MIR130B, POLD3 and FAM168A. (E) Four-cliques of RP11-401P9.4.1, MIR600HG, KRIT1 and RAD52. (F) Four-cliques of RP11-401P9.4.1, MIR600HG, RBM5 and KRIT1. (G) Four-cliques of RP11-401P9.4.1, MIR600HG, NFYA and KRIT1. (H) Four-cliques of RP11-401P9.4.1, MIR600HG, KRIT1 and LAP3. (I) Four-cliques of RP11-401P9.4.1, MIR600HG, RBM5 and RAD52. (J) Four-cliques of RP11-401P9.4.1, MIR600HG, NFYA and RAD52. (K) Four-cliques of RP11-401P9.4.1, MIR600HG, LAP3 and RAD52. (L) Four-cliques of RP11-401P9.4.1, MIR600HG, ALS2 and RAD52. (M) Four-cliques of RP11-401P9.4.1, MIR600HG, RBM5 and ALS2. (N) Four-cliques of RP11-401P9.4.1, MIR600HG, RBM5 and DPM1. (O) Four-cliques of RP11-401P9.4.1, MIR600HG, NFYA and ALS2.
Figure 5 Kaplan–Meier plot of CLPTM1L. The high expression of CLPTM1L was associated with high risk.
We then move on to the functional analysis of a few new genes in four cliquesas other genes have been discussed above (Figures 4I, J, K, L, M and O). Figures 4A, B shares the same mRNAs and miRNA. Cleft lip and palate transmembrane protein 1-like (CLPTM1L) was frequently found to be up regulated in cisplatin-resistant human ovarian cancer cells (Yamamoto et al., 2001; James et al., 2012) and also associated with cisplatin-induced apoptosis. Recent study has uncovered the role of CLPTM1L in modulating cell survival signaling and identified it as a lung cancer risk gene (James et al., 2012; James et al., 2014), which was also validated by our results (Figure 5). The tRNA methytransferase NSUN2 has been shown to stimulate cell growth in a cell cycle-dependent manner due to its dynamic expression levels throughout the cell cycle (Xing et al., 2015), suggesting its great roles in promoting tumorigenesis. Accumulating studies have revealed the function of MIR124-3 in inhibiting invasion and migration of malignant cells both in ovarian cancer and hepatocellular cancer (Zhang et al., 2013; Cai et al., 2017; Yuan et al., 2017). Although the roles of two lncRNAs CTD-2142D14.1.1 and RP11-318I4.1.1 have not been elucidated yet, it can be deduced by their significant co-expression with CLPTM1L and NSUN2 from TCGA ovarian cancer data, indicating their cooperation in contributing to the progression of ovarian cancer. As for Figures 4C and D, FAM168A, also known as tongue cancer resistance-related protein 1(TCRP1), has previously been found to mediate resistance to chemotherapy in several cancers including ovarian and lung cancers (Gu et al., 2017; Liu et al., 2017b), exerting similar function to the triplet in Figure 2H, which explained their clusters and emphasized their importance on drug resistance in ovarian cancer. Besides, FAM168A has been identified as a prognostic marker in renal cancer, endometrial cancer, and urothelial cancer (Gu et al., 2017). The new gene in Figures 4E, F, G, H is Krev interaction trapped protein 1 (KRIT1), which is responsible for cerebrovascular disease (Goitre et al., 2010). Consistent with other partners in this four cliques, KRIT1 has been shown to be involved in cell growth, DNA damage, angiogenesis, and cell death (Goitre et al., 2010; Orso et al., 2013), indicating that it may function as a potential mechanism of susceptibility to ovarian tumorigenesis. Figure 4N contains a new gene DPM1 that encodes dolichol-phosphate mannosyltransferase. Previous work has examined increased expression levels of DPM1 in primary EOC of various histotypes (Ramakrishna et al., 2010). Moreover, DPM1 has been considered as a prognostic marker in liver cancer. Altogether, these new genes in four cliques are capable of undertaking ovarian cancer risk candidate genes on account of their contributions to the development of ovarian cancer.
Identification of the Common Clique Between Ovarian Cancer and Uterine Corpus Endometrial Carcinoma
In our previous paper (Liu et al., 2017a), we have identified the cliques in uterine corpus endometrial carcinoma (UCEC). Both uterine corpus endometrial carcinoma and ovarian cancer were common female cancers. We compared the cliques in ovarian cancer with the cliques in UCEC and found one overlapped clique: mRNA KRTAP24-1–lncRNA LL22NC03-121E8.3.1–miRNA MIR2116. KRTAP24-1, also called KAP24.1, is a cuticular hair keratin-associated protein (KAP) and belongs to human KAP family because of its structure and location (Rogers et al., 2007). More than as the intermediate filament forming proteins of epithelial cells, keratins have recently been regarded as regulators of protein synthesis, cell motility, and membrane traffic and signaling (Karantza, 2011). Furthermore, keratins have been extensively involved in cancer cell invasion and metastasis as well as acted as diagnostic markers in epithelial cancers (Moll et al., 2008), suggesting its vital roles in the development of cancer. MIR2116 is a newly found miRNA expressed in EOC (Wyman et al., 2009) and associated with bacterial infection (Zhao et al., 2018; Modak et al., 2019). Our results expand its function in the pathogenesis of gynecological cancers. Although LL22NC03-121E8.3.1 is a novel lncRNA with little reports, we can speculate its function by its interacted protein fragile X mental retardation 1 (FMR1), which is expressed at high levels and considered as a prognostic marker in gynecological cancers including ovarian and endometrial cancers (Hagerman et al., 2017), indicating its potential roles in facilitating tumorigenesis. More recently, several lines of evidence have supported that the estrogen receptor (ER) was highly expressed in ovarian cancers and regulated cell growth, which attracted endocrine therapy for ovarian cancer treatment (Lee et al., 2014; Langdon et al., 2017). Meanwhile, the main histopathologic type of endometrial cancer (type I) is strongly related to unopposed estrogen (Lee et al., 2014), and hormones represent a promising target especially in advanced or recurrent endometrial cancer. Taken these into account, hormones become common features in the progression of ovarian and endometrial cancers. As this triplet is both found in OC and UCEC, it is reasonable to infer that these three RNAs have a strong relationship with hormone production through regulating their expression and stability, which needs more investigations. In addition, our results provide new hormone-related RNAs about gynecological oncogenesis.
The regulation roles of lncRNAs and miRNAs are largely unknown. Recent studies have shown that lncRNAs and miRNA can cross talk and create a competition for binding between miRNA, lncRNA, and regulatory target, which is called sponge effect. To investigate the functional relationship between miRNA, lncRNA, and mRNA on a genome-wide scale is essential for understanding the regulatory roles of lncRNA and miRNA. Therefore, we analyzed the RNA-seq profiles of 407 ovarian cancer patients and constructed an integrative network of 20,424 coding mRNAs, 10,412 lncRNAs, and 742 miRNAs using VIF regression. Then, all the mRNA–lncRNA–miRNA cliques were detected using R package igraph. Specifically, we analyzed the triplets that showed significant correlations with survival and stage of ovarian cancer. Our results not only provide novel insights of the regulatory mechanisms among mRNAs, lncRNAs, and miRNAs, but also shield light on the tumorigenesis mechanisms of ovarian cancer.
YZ and JJ conceived of the study and participated in its design and coordination. XZ and WH conducted the database search and analysis. TH and BX conducted the network construction. YZ conducted the biological analysis and wrote the first draft. The draft was improved through discussion and editing by all the authors, who read and approved the final manuscript.
This study was supported by the National Key Technology Support Program of China (2015BAI12B00), National Key R&D Program of China (2018YFC0910403), National Natural Science Foundation of China (31701111, 31701151), Shanghai Municipal Science and Technology Major Project (2017SHZDZX01), China Postdoctoral Science Foundation (2018M642305), Natural Science Foundation of Jiangsu Province (BK20170295), Youth Innovation Promotion Association of Chinese Academy of Sciences (2016245), and Basic Research Project of Changzhou (CJ20179024).
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
We would like to thank Xiangyin Kong and Landian Hu for their continuous enthusiasm and support. We would like to thank also Yudong Cai and Changping Wu for editing different versions of this manuscript.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2019.00751/full#supplementary-material
Supplementary Material 1 | The TCGA sample IDs and their clinical information.
Supplementary Material 2 | Constructed integrative network.
Supplementary Material 3 | mRNA–lncRNA–miRNA triplets.
Supplementary Material 4 | Four-cliques containing mRNA–lncRNA–miRNA triplets.
Akrami, R., Jacobsen, A., Hoell, J., Schultz, N., Sander, C., Larsson, E. (2013). Comprehensive analysis of long non-coding RNAs in ovarian cancer reveals global patterns and targeted DNA amplification. PLoS One 8 (11), e80306. doi: 10.1371/journal.pone.0080306
Baldwin, E. L., May, L. F., Justice, A. N., Martin, C. L., Ledbetter, D. H. (2008). Mechanisms and consequences of small supernumerary marker chromosomes: from Barbara McClintock to modern genetic-counseling issues. Am. J. Hum. Genet. 82 (2), 398–410. doi: 10.1016/j.ajhg.2007.10.013
Braconi, C., Kogure, T., Valeri, N., Huang, N., Nuovo, G., Costinean, S., et al. (2011). microRNA-29 can regulate expression of the long non-coding RNA gene MEG3 in hepatocellular cancer. Oncogene 30 (47), 4750–4756. doi: 10.1038/onc.2011.193
Cai, Q. Q., Dong, Y. W., Wang, R., Qi, B., Guo, J. X., Pan, J., et al. (2017). MiR-124 inhibits the migration and invasion of human hepatocellular carcinoma cells by suppressing integrin alphaV expression. Sci. Rep. 7, 40733. doi: 10.1038/srep40733
Cesana, M., Cacchiarelli, D., Legnini, I., Santini, T., Sthandier, O., Chinappi, M., et al. (2011). A long noncoding RNA controls muscle differentiation by functioning as a competing endogenous RNA. Cell 147 (2), 358–369. doi: 10.1016/j.cell.2011.09.028
Chen, C., Hu, Y., Li, L. (2016). NRP1 is targeted by miR-130a and miR-130b, and is associated with multidrug resistance in epithelial ovarian cancer based on integrated gene network analysis. Mol. Med. Rep. 13 (1), 188–196. doi: 10.3892/mmr.2015.4556
Chen, X., Wang, L., Qu, J., Guan, N. N., Li, J. Q. (2018a). Predicting miRNA-disease association based on inductive matrix completion. Bioinformatics 34 (24), 4256–4265. doi: 10.1093/bioinformatics/bty503
Chen, X., Xie, D., Wang, L., Zhao, Q., You, Z. H., Liu, H. (2018b). BNPMDA: bipartite network projection for miRNA-disease association prediction. Bioinformatics 34 (18), 3178–3186. doi: 10.1093/bioinformatics/bty333
Chen, X., Yan, C. C., Zhang, X., You, Z. H. (2017). Long non-coding RNAs and complex diseases: from experimental results to computational models. Brief Bioinform. 18 (4), 558–576. doi: 10.1093/bib/bbw060
Chen, X., Yin, J., Qu, J., Huang, L. (2018c). MDHGI: matrix decomposition and heterogeneous graph inference for miRNA-disease association prediction. PLoS Comput. Biol. 14 (8), e1006418. doi: 10.1371/journal.pcbi.1006418
Chiyomaru, T., Yamamura, S., Fukuhara, S., Yoshino, H., Kinoshita, T., Majid, S., et al. (2013). Genistein inhibits prostate cancer cell growth by targeting miR-34a and oncogenic HOTAIR. PLoS One 8 (8), e70372. doi: 10.1371/journal.pone.0070372
Crawford, N. P., Walker, R. C., Lukes, L., Officewala, J. S., Williams, R. W., Hunter, K. W. (2008). The diasporin pathway: a tumor progression-related transcriptional network that predicts breast cancer survival. Clin. Exp. Metastasis 25 (4), 357–369. doi: 10.1007/s10585-008-9146-6
Duckworth, C., Zhang, L., Carroll, S. L., Ethier, S. P., Cheung, H. W. (2016). Overexpression of GAB2 in ovarian cancer cells promotes tumor growth and angiogenesis by upregulating chemokine expression. Oncogene 35 (31), 4036–4047. doi: 10.1038/onc.2015.472
Ellinger, J., Alam, J., Rothenburg, J., Deng, M., Schmidt, D., Syring, I., et al. (2015). The long non-coding RNA lnc-ZNF180-2 is a prognostic biomarker in patients with clear cell renal cell carcinoma. Am. J. Cancer Res. 5 (9), 2799.
Fan, M., Li, X., Jiang, W., Huang, Y., Li, J., Wang, Z. (2013). A long non-coding RNA, PTCSC3, as a tumor suppressor and a target of miRNAs in thyroid cancer cells. Exp. Ther. Med. 5 (4), 1143–1146. doi: 10.3892/etm.2013.933
Farina, B., Fattorusso, R., Pellecchia, M. (2011). Targeting zinc finger domains with small molecules: solution structure and binding studies of the RanBP2-type zinc finger of RBM5. Chem. Biochem. 12 (18), 2837–2845. doi: 10.1002/cbic.201100582
Fico, A., Fiorenzano, A., Pascale, E., Patriarca, E. J., Minchiotti, G. (2019). Long non-coding RNA in stem cell pluripotency and lineage commitment: functions and evolutionary conservation. Cell. Mol. Life Sci. 76 (8), 1459–1471. doi: 10.1007/s00018-018-3000-z
Furió-Tarí, P., Tarazona, S., Gabaldón, T., Enright, A. J., Conesa, A. (2016). Spongescan: a web for detecting microRNA binding elements in lncRNA sequences. Nucleic Acids Res. 44, W176–W180. doi: 10.1093/nar/gkw443
Gao, Y., Meng, H., Liu, S., Hu, J., Zhang, Y., Jiao, T., et al. (2015). LncRNA-HOST2 regulates cell biological behaviors in epithelial ovarian cancer through a mechanism involving microRNA let-7b. Hum. Mol. Genet. 24 (3), 841–852. doi: 10.1093/hmg/ddu502
Garipov, A., Li, H., Bitler, B. G., Thapa, R. J., Balachandran, S., Zhang, R. (2013). NF-YA underlies EZH2 upregulation and is essential for proliferation of human epithelial ovarian cancer cells. Mol. Cancer. Res. 11 (4), 360–369. doi: 10.1158/1541-7786.MCR-12-0661
Gautam, M., Jara, J. H., Sekerkova, G., Yasvoina, M. V., Martina, M., Ozdinler, P. H. (2016). Absence of alsin function leads to corticospinal motor neuron vulnerability via novel disease mechanisms. Hum. Mol. Genet. 25 (6), 1074–1087. doi: 10.1093/hmg/ddv631
Goitre, L., Balzac, F., Degani, S., Degan, P., Marchi, S., Pinton, P., et al. (2010). KRIT1 regulates the homeostasis of intracellular reactive oxygen species. PLoS One 5 (7), e11786. doi: 10.1371/journal.pone.0011786
Gu, Y., Zhang, Z., Yin, J., Ye, J., Song, Y., Liu, H., et al. (2017). Epigenetic silencing of miR-493 increases the resistance to cisplatin in lung cancer by targeting tongue cancer resistance-related protein 1(TCRP1). J. Exp. Clin. Cancer Res. 36 (1), 114. doi: 10.1186/s13046-017-0582-5
Guan, B., Wu, K., Zeng, J., Xu, S., Mu, L., Yang, G., et al. (2016). Tumor-suppressive microRNA-218 inhibits tumor angiogenesis via targeting the mTOR component RICTOR in prostate cancer. Oncotarget 8 (5), 8162. doi: 10.18632/oncotarget.14131
He, X., Huang, Q., Qiu, X., Liu, X., Sun, G., Guo, J., et al. (2015). LAP3 promotes glioma progression by regulating proliferation, migration and invasion of glioma cells. Int. J. Biol. Macromol. 72 (72C), 1081–1089. doi: 10.1016/j.ijbiomac.2014.10.021
Hromas, R., Kim, H. S., Sidhu, G., Williamson, E., Jaiswal, A., Totterdale, T. A., et al. (2017). The endonuclease EEPD1 mediates synthetic lethality in RAD52-depleted BRCA1 mutant breast cancer cells. Breast Cancer Res. 19 (1), 122. doi: 10.1186/s13058-017-0912-8
Imig, J., Brunschweiger, A., Brummer, A., Guennewig, B., Mittal, N., Kishore, S., et al. (2015). miR-CLIP capture of a miRNA targetome uncovers a lincRNA H19-miR-106a interaction. Nat. Chem. Biol. 11 (2), 107–114. doi: 10.1038/nchembio.1713
James, M. A., Vikis, H. G., Tate, E., Rymaszewski, A. L., You, M. (2014). CRR9/CLPTM1L regulates cell survival signaling and is required for Ras transformation and lung tumorigenesis. Cancer Res. 74 (4), 1116–1127. doi: 10.1158/0008-5472.CAN-13-1617
James, M. A., Wen, W., Wang, Y., Byers, L. A., Heymach, J. V., Coombes, K. R., et al. (2012). Functional characterization of CLPTM1L as a lung cancer risk candidate gene in the 5p15.33 locus. PLoS One 7 (6), e36116. doi: 10.1371/journal.pone.0036116
Kim, D., Pertea, G., Trapnell, C., Pimentel, H., Kelley, R., Salzberg, S. L. (2013). TopHat2: accurate alignment of transcriptomes in the presence of insertions, deletions and gene fusions. Genome Biol. 14 (4), R36. doi: 10.1186/gb-2013-14-4-r36
Kim, Y. S., Hwan, J. D., Bae, S., Bae, D. H., Shick, W. A. (2010). Identification of differentially expressed genes using an annealing control primer system in stage III serous ovarian carcinoma. BMC Cancer 10, 576. doi: 10.1186/1471-2407-10-576
Koldehoff, M., Zakrzewski, J. L., Klein-Hitpass, L., Beelen, D. W., Elmaagacli, A. H. (2008). Gene profiling of growth factor independence 1B gene (Gfi-1B) in leukemic cells. Int. J. Hematol. 87 (1), 39–47. doi: 10.1007/s12185-007-0013-z
Lee, W. L., Yen, M. S., Chao, K. C., Yuan, C. C., Ng, H. T., Chao, H. T., et al. (2014). Hormone therapy for patients with advanced or recurrent endometrial cancer. J. Chin. Med. Assoc. 77 (5), 221–226. doi: 10.1016/j.jcma.2014.02.007
Legnini, I., Morlando, M., Mangiavacchi, A., Fatica, A., Bozzoni, I. (2014). A feedforward regulatory loop between HuR and the long noncoding RNA linc-MD1 controls early phases of myogenesis. Mol. Cell. 53 (3), 506–514. doi: 10.1186/1471-2458-13-507
Li, L., Liu, B., Wapinski, O., Tsai, M. C., Qu, K., Zhang, J., et al. (2013). Targeted disruption of hotair leads to homeotic transformation and gene derepression. Cell. Rep. 5 (1), 3–12. doi: 10.1016/j.celrep.2013.09.003
Liu, C., Zhang, Y. H., Deng, Q., Li, Y., Huang, T., Zhou, S., et al. (2017a). Cancer-related triplets of mRNA-lncRNA-miRNA revealed by integrative network in uterine corpus endometrial carcinoma. Biomed. Res. Int. 2017, 3859582. doi: 10.1155/2017/3859582
Liu, X., Feng, M., Zheng, G., Gu, Y., Wang, C., He, Z. (2017b). TCRP1 expression is associated with platinum sensitivity in human lung and ovarian cancer cells. Oncol. Lett. 13 (3), 1398–1405. doi: 10.3892/ol.2016.5534
Martens-Uzunova, E. S., Bottcher, R., Croce, C. M., Jenster, G., Visakorpi, T., Calin, G. A. (2014). Long noncoding RNA in prostate, bladder, and kidney cancer. Eur. Urol. 65 (6), 1140–1151. doi: 10.1016/j.eururo.2013.12.003
Meng, J., Wang, J. (2015). Role of SNARE proteins in tumourigenesis and their potential as targets for novel anti-cancer therapeutics. Biochim. Biophys. Acta 1856 (1), 1–12. doi: 10.1016/j.bbcan.2015.04.002
Modak, J. M., Roy-O’Reilly, M., Zhu, L., Staff, I., McCullough, L. D. (2019). Differential microribonucleic acid expression in cardioembolic stroke. J. Stroke Cerebrovasc. Dis. 28 (1), 121–124. doi: 10.1016/j.jstrokecerebrovasdis.2018.09.018
Naydenov, N. G., Joshi, S., Feygin, A., Saini, S., Litovchick, L., Ivanov, A. I. (2018). A membrane fusion protein, Ykt6, regulates epithelial cell migration via microRNA-mediated suppression of Junctional Adhesion Molecule A. Cell. Cycle 17 (14), 1812–1831. doi: 10.1080/15384101.2018.1496755
Nguyen, C. D., Mansfield, R. E., Leung, W., Vaz, P. M., Loughlin, F. E., Grant, R. P., et al. (2011). Characterization of a family of RanBP2-type Zinc fingers that can recognize single-stranded RNA. J. Mol. Biol. 407 (2), 273. doi: 10.1016/j.jmb.2010.12.041
Orso, F., Balzac, F., Marino, M., Lembo, A., Retta, S. F., Taverna, D. (2013). miR-21 coordinates tumor growth and modulates KRIT1 levels. Biochem. Biophys. Res. Commun. 438 (1), 90–96. doi: 10.1016/j.bbrc.2013.07.031
Park, E. R., Kim, S. B., Lee, J. S., Kim, Y. H., Lee, D. H., Cho, E. H., et al. (2017). The mitochondrial hinge protein, UQCRH, is a novel prognostic factor for hepatocellular carcinoma. Cancer Med. 6 (4), 749–760. doi: 10.1002/cam4.1042
Parris, T. Z., Rönnerman, E. W., Engqvist, H., Biermann, J., Truvé, K., Nemes, S., et al. (2018). Genome-wide multi-omics profiling of the 8p11-p12 amplicon in breast carcinoma. Oncotarget 9 (35), 24140–24154. doi: 10.18632/oncotarget.25329
Ramakrishna, M., Williams, L. H., Boyle, S. E., Bearfoot, J. L., Sridhar, A., Speed, T. P., et al. (2010). Identification of candidate growth promoting genes in ovarian cancer through integrated copy number and expression analysis. PLoS One 5 (4), e9983. doi: 10.1371/journal.pone.0009983.t001
Ribeiro, A. O., Schoof, C. R., Izzotti, A., Pereira, L. V., Vasques, L. R. (2014). MicroRNAs: modulators of cell identity, and their applications in tissue engineering. Microrna 3 (1), 45–53. doi: 10.2174/2211536603666140522003539
Rogers, M. A., Winter, H., Langbein, L., Wollschlager, A., Praetzel-Wunder, S., Jave-Suarez, L. F., et al. (2007). Characterization of human KAP24.1, a cuticular hair keratin-associated protein with unusual amino-acid composition and repeat structure. J. Invest. Dermatol. 127 (5), 1197–1204. doi: 10.1038/sj.jid.5700702
Song, J., Xu, Q., Zhang, H., Yin, X., Zhu, C., Zhao, K., et al. (2018). Five key lncRNAs considered as prognostic targets for predicting pancreatic ductal adenocarcinoma. J. Cell. Biochem. 119 (6), 4559–4569. doi: 10.1002/jcb.26598
Stafford, J. L., Dyson, G., Levin, N. K., Chaudhry, S., Rosati, R., Kalpage, H., et al. (2017). Reanalysis of BRCA1/2 negative high risk ovarian cancer patients reveals novel germline risk loci and insights into missing heritability. PLoS One 12 (6), e0178450. doi: 10.1371/journal.pone.0178450
Sundaram, G. M., Veera Bramhachari, P. (2017). Molecular interplay of pro-inflammatory transcription factors and non-coding RNAs in esophageal squamous cell carcinoma. Tumour Biol. 39 (6), 1–12. doi: 10.1177/1010428317705760
The Cancer Genome Atlas Research Network, Weinstein, J. N., Collisson, E. A., Mills, G. B., Shaw, K. R., Ozenberger, B. A., et al. (2013). The cancer genome atlas pan-cancer analysis project. Nat. Genet. 45 (10), 1113–1120. doi: 10.1038/ng.2764
Tu, J., Tian, G., Cheung, H. H., Wei, W., Lee, T. L. (2018). Gas5 is an essential lncRNA regulator for self-renewal and pluripotency of mouse embryonic stem cells and induced pluripotent stem cells. Stem. Cell Res. Ther. 9 (1), 71. doi: 10.1186/s13287-018-0813-5
Vetter, M. H., Hays, J. L. (2018). Use of targeted therapeutics in epithelial ovarian cancer: a review of current literature and future directions. Clin. Ther. 40 (3), 361–371. doi: 10.1016/j.clinthera.2018.01.012
Wake, N. C., Ricketts, C. J., Morris, M. R., Prigmore, E., Gribble, S. M., Skytte, A. B., et al. (2013). UBE2QL1 is disrupted by a constitutional translocation associated with renal tumor predisposition and is a novel candidate renal tumor suppressor gene. Hum. Mutat. 34 (12), 1650–1661. doi: 10.1002/humu.22433
Wang, X., Li, M., Wang, Z., Han, S., Tang, X., Ge, Y., et al. (2015a). Silencing of long noncoding RNA MALAT1 by miR-101 and miR-217 inhibits proliferation, migration, and invasion of esophageal squamous cell carcinoma cells. J. Biol. Chem. 290 (7), 3925–3935. doi: 10.1074/jbc.M114.596866
Wang, X., Shi, L., Deng, Y., Qu, M., Mao, S., Xu, L., et al. (2015b). Inhibition of leucine aminopeptidase 3 suppresses invasion of ovarian cancer cells through down-regulation of fascin and MMP-2/9. Eur. J. Pharmacol. 768, 116–122. doi: 10.1016/j.ejphar.2015.10.039
Wang, Y., Xu, Z., Jiang, J., Xu, C., Kang, J., Xiao, L., et al. (2013). Endogenous miRNA sponge lincRNA-RoR regulates Oct4, Nanog, and Sox2 in human embryonic stem cell self-renewal. Dev. Cell. 25 (1), 69–80. doi: 10.1016/j.devcel.2013.03.002
Willis, S., Villalobos, V. M., Gevaert, O., Abramovitz, M., Williams, C., Sikic, B. I., et al. (2016). Single gene prognostic biomarkers in ovarian cancer: a meta-analysis. PLoS One 11 (2), e0149183. doi: 10.1371/journal.pone.0149183
Wyman, S. K., Parkin, R. K., Mitchell, P. S., Fritz, B. R., O’Briant, K., Godwin, A. K., et al. (2009). Repertoire of microRNAs in epithelial ovarian cancer as determined by next generation sequencing of small RNA cDNA libraries. PLoS One 4 (4), e5311. doi: 10.1371/journal.pone.0005311
Xing, J., Yi, J., Cai, X., Tang, H., Liu, Z., Zhang, X., et al. (2015). NSun2 promotes cell growth via elevating cyclin-dependent Kinase 1 translation. Mol. Cell. Biol. 35 (23), 4043–4052. doi: 10.1128/MCB.00742-15
Xu, C., Zhang, Y., Wang, Q., Xu, Z., Jiang, J., Gao, Y., et al. (2016). Long non-coding RNA GAS5 controls human embryonic stem cell self-renewal by maintaining NODAL signalling. Nat. Commun. 7, 13287. doi: 10.1038/ncomms13287
Yamamoto, K., Okamoto, A., Isonishi, S., Ochiai, K., Ohtake, Y. (2001). A novel gene, CRR9, which was up-regulated in CDDP-resistant ovarian tumor cell line, was associated with apoptosis. Biochem. Biophys. Res. Commun. 280 (4), 1148–1154. doi: 10.1006/bbrc.2001.4250
Yamamoto, N., Kinoshita, T., Nohata, N., Itesako, T., Yoshino, H., Enokida, H., et al. (2013). Tumor suppressive microRNA-218 inhibits cancer cell migration and invasion by targeting focal adhesion pathways in cervical squamous cell carcinoma. Int. J. Oncol. 42 (5), 1523–1532. doi: 10.3892/ijo.2013.1851
Yang, C., Cai, J., Wang, Q., Tang, H., Cao, J., Wu, L., et al. (2012). Epigenetic silencing of miR-130b in ovarian cancer promotes the development of multidrug resistance by targeting colony-stimulating factor 1. Gynecol. Oncol. 124 (2), 325–334. doi: 10.1016/j.ygyno.2011.10.013
Yang, Q. Y., Li, J. H., Wang, Q. Y., Wu, Y., Qin, J. L., Cheng, J. J., et al. (2014). MTA1 promotes cell proliferation via DNA damage repair in epithelial ovarian cancer. Genet. Mol. Res. 13 (4), 10269–10278. doi: 10.4238/2014.December.4.21
Yang, Y., Hentati, A., Deng, H. X., Dabbagh, O., Sasaki, T., Hirano, M., et al. (2001). The gene encoding alsin, a protein with three guanine-nucleotide exchange factor domains, is mutated in a form of recessive amyotrophic lateral sclerosis. Nat. Genet. 29 (2), 160. doi: 10.1038/ng1001-160
Ye, C., Tao, R., Cao, Q., Zhu, D., Wang, Y., Wang, J., et al. (2016). Whole-genome DNA methylation and hydroxymethylation profiling for HBV-related hepatocellular carcinoma. Int. J. Oncol. 49 (2), 589–602. doi: 10.3892/ijo.2016.3535
Yoon, J. H., Abdelmohsen, K., Kim, J., Yang, X., Martindale, J. L., Tominaga-Yamanaka, K., et al. (2013). Scaffold function of long non-coding RNA HOTAIR in protein ubiquitination. Nat. Commun. 4, 2939. doi: 10.1038/ncomms3939
Yuan, L., Li, S., Zhou, Q., Wang, D., Zou, D., Shu, J., et al. (2017). MiR-124 inhibits invasion and induces apoptosis of ovarian cancer cells by targeting programmed cell death 6. Oncol. Lett. 14 (6), 7311–7317. doi: 10.3892/ol.2017.7157
Zhang, H., Zhu, C., Yi, Z., Ming, L., Wu, L., Yang, X., et al. (2015). Long non-coding RNA expression profiles of hepatitis C virus-related dysplasia and hepatocellular carcinoma. Oncotarget 6 (41), 43770. doi: 10.18632/oncotarget.6087
Zhao, Q., Yang, J., Li, L., Cao, D., Yu, M., Shen, K., et al. (2017). Germline and somatic mutations in homologous recombination genes among Chinese ovarian cancer patients detected using next-generation sequencing. J. Gynecol. Oncol. 28 (4), e39. doi: 10.3802/jgo.2017.28.e39
Keywords: variance inflation factor regression, mRNA–lncRNA–miRNA cliques, regulatory network construction, ovarian cancer, functional regulator
Citation: Zhou Y, Zheng X, Xu B, Hu W, Huang T and Jiang J (2019) The Identification and Analysis of mRNA–lncRNA–miRNA Cliques From the Integrative Network of Ovarian Cancer. Front. Genet. 10:751. doi: 10.3389/fgene.2019.00751
Received: 30 January 2019; Accepted: 17 July 2019;
Published: 21 August 2019.
Edited by:Marco Pellegrini, Institute of Computer Science and Telematics (IIT), Italy
Reviewed by:Suman Ghosal, National Institutes of Health (NIH), United States
Qi Zhao, Shenyang Aerospace University, China
Jia Qu, China University of Mining and Technology, China, in collaboration with reviewer QZ
Copyright © 2019 Zhou, Zheng, Xu, Hu, Huang and Jiang. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.