Abstract
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.
Introduction
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 inflation 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.
Methods
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
| mRNA | lncRNA | miRNA | Total | |
|---|---|---|---|---|
| Expression | 20,462 | 10,419 | 742 | 31,623 |
| Network | 16,667 | 4,796 | 207 | 21,670 |
| Triplets | 633 | 1,184 | 169 | 1,986 |
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.
Figure 1
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
Figure 3
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
Figure 5
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.
Conclusion
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.
Funding
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).
Statements
Author contributions
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.
Acknowledgments
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.
Conflict of interest
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.
Supplementary material
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 1The TCGA sample IDs and their clinical information.
Supplementary Material 2Constructed integrative network.
Supplementary Material 3mRNA–lncRNA–miRNA triplets.
Supplementary Material 4Four-cliques containing mRNA–lncRNA–miRNA triplets.
References
1
AkramiR.JacobsenA.HoellJ.SchultzN.SanderC.LarssonE. (2013). Comprehensive analysis of long non-coding RNAs in ovarian cancer reveals global patterns and targeted DNA amplification. PLoS One8 (11), e80306. doi: 10.1371/journal.pone.0080306
2
AndersS.PylP. T.HuberW. (2015). HTSeq–a Python framework to work with high-throughput sequencing data. Bioinformatics31 (2), 166–169. doi: 10.1093/bioinformatics/btu638
3
AndersenP. K.GillR. D. (1982). Cox’s regression model for counting processes: a large sample study. Ann. Statist.10 (4), 1100–1120. doi: 10.1214/aos/1176345976
4
BaldwinE. L.MayL. F.JusticeA. N.MartinC. L.LedbetterD. 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
5
BraconiC.KogureT.ValeriN.HuangN.NuovoG.CostineanS.et al. (2011). microRNA-29 can regulate expression of the long non-coding RNA gene MEG3 in hepatocellular cancer. Oncogene30 (47), 4750–4756. doi: 10.1038/onc.2011.193
6
CaiQ. Q.DongY. W.WangR.QiB.GuoJ. X.PanJ.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
7
CesanaM.CacchiarelliD.LegniniI.SantiniT.SthandierO.ChinappiM.et al. (2011). A long noncoding RNA controls muscle differentiation by functioning as a competing endogenous RNA. Cell147 (2), 358–369. doi: 10.1016/j.cell.2011.09.028
8
ChenC.HuY.LiL. (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
9
ChenX.HuangL. (2017). LRSSLMDA: laplacian regularized sparse subspace learning for miRNA-disease association prediction. PLoS Comput. Biol.13 (12), e1005912. doi: 10.1371/journal.pcbi.1005912
10
ChenX.WangL.QuJ.GuanN. N.LiJ. Q. (2018a). Predicting miRNA-disease association based on inductive matrix completion. Bioinformatics34 (24), 4256–4265. doi: 10.1093/bioinformatics/bty503
11
ChenX.XieD.WangL.ZhaoQ.YouZ. H.LiuH. (2018b). BNPMDA: bipartite network projection for miRNA-disease association prediction. Bioinformatics34 (18), 3178–3186. doi: 10.1093/bioinformatics/bty333
12
ChenX.YanC. C.ZhangX.YouZ. 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
13
ChenX.YanG. Y. (2013). Novel human lncRNA-disease association inference based on lncRNA expression profiles. Bioinformatics29 (20), 2617–2624. doi: 10.1093/bioinformatics/btt426
14
ChenX.YinJ.QuJ.HuangL. (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
15
ChiyomaruT.YamamuraS.FukuharaS.YoshinoH.KinoshitaT.MajidS.et al. (2013). Genistein inhibits prostate cancer cell growth by targeting miR-34a and oncogenic HOTAIR. PLoS One8 (8), e70372. doi: 10.1371/journal.pone.0070372
16
CrawfordN. P.WalkerR. C.LukesL.OfficewalaJ. S.WilliamsR. W.HunterK. W. (2008). The diasporin pathway: a tumor progression-related transcriptional network that predicts breast cancer survival. Clin. Exp. Metastasis25 (4), 357–369. doi: 10.1007/s10585-008-9146-6
17
DuckworthC.ZhangL.CarrollS. L.EthierS. P.CheungH. W. (2016). Overexpression of GAB2 in ovarian cancer cells promotes tumor growth and angiogenesis by upregulating chemokine expression. Oncogene35 (31), 4036–4047. doi: 10.1038/onc.2015.472
18
EllingerJ.AlamJ.RothenburgJ.DengM.SchmidtD.SyringI.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.
19
FanM.LiX.JiangW.HuangY.LiJ.WangZ. (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
20
FarinaB.FattorussoR.PellecchiaM. (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
21
FicoA.FiorenzanoA.PascaleE.PatriarcaE. J.MinchiottiG. (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
22
Furió-TaríP.TarazonaS.GabaldónT.EnrightA. J.ConesaA. (2016). Spongescan: a web for detecting microRNA binding elements in lncRNA sequences. Nucleic Acids Res.44, W176–W180. doi: 10.1093/nar/gkw443
23
Gabor CsardiT. N. (2006). The igraph software package for complex network research. IntJ. Complex Syst., 1695.
24
GaoY.MengH.LiuS.HuJ.ZhangY.JiaoT.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
25
GaripovA.LiH.BitlerB. G.ThapaR. J.BalachandranS.ZhangR. (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
26
GautamM.JaraJ. H.SekerkovaG.YasvoinaM. V.MartinaM.OzdinlerP. 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
27
GoitreL.BalzacF.DeganiS.DeganP.MarchiS.PintonP.et al. (2010). KRIT1 regulates the homeostasis of intracellular reactive oxygen species. PLoS One5 (7), e11786. doi: 10.1371/journal.pone.0011786
28
GuY.ZhangZ.YinJ.YeJ.SongY.LiuH.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
29
GuanB.WuK.ZengJ.XuS.MuL.YangG.et al. (2016). Tumor-suppressive microRNA-218 inhibits tumor angiogenesis via targeting the mTOR component RICTOR in prostate cancer. Oncotarget8 (5), 8162. doi: 10.18632/oncotarget.14131
30
HagermanR. J.Berry-KravisE.HazlettH. C.BaileyD. B.Jr.MoineH.KooyR. F.et al. (2017). Fragile X syndrome. Nat. Rev. Dis. Primers3, 17065. doi: 10.1038/nrdp.2017.65
31
HeX.HuangQ.QiuX.LiuX.SunG.GuoJ.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
32
HelwakA.KudlaG.DudnakovaT.TollerveyD. (2013). Mapping the human miRNA interactome by CLASH reveals frequent noncanonical binding. Cell153 (3), 654–665. doi: 10.1016/j.cell.2013.03.043
33
HromasR.KimH. S.SidhuG.WilliamsonE.JaiswalA.TotterdaleT. 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
34
ImigJ.BrunschweigerA.BrummerA.GuennewigB.MittalN.KishoreS.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
35
JamesM. A.VikisH. G.TateE.RymaszewskiA. L.YouM. (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
36
JamesM. A.WenW.WangY.ByersL. A.HeymachJ. V.CoombesK. R.et al. (2012). Functional characterization of CLPTM1L as a lung cancer risk candidate gene in the 5p15.33 locus. PLoS One7 (6), e36116. doi: 10.1371/journal.pone.0036116
37
KarantzaV. (2011). Keratins in health and cancer: more than mere epithelial cell markers. Oncogene30 (2), 127–138. doi: 10.1038/onc.2010.456
38
KimD.PerteaG.TrapnellC.PimentelH.KelleyR.SalzbergS. 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
39
KimY. S.HwanJ. D.BaeS.BaeD. H.ShickW. A. (2010). Identification of differentially expressed genes using an annealing control primer system in stage III serous ovarian carcinoma. BMC Cancer10, 576. doi: 10.1186/1471-2407-10-576
40
KoldehoffM.ZakrzewskiJ. L.Klein-HitpassL.BeelenD. W.ElmaagacliA. 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
41
LangdonS. P.GourleyC.GabraH.StanleyB. (2017). Endocrine therapy in epithelial ovarian cancer. Expert Rev. Anticancer Ther.17 (2), 109–117. doi: 10.1080/14737140.2017.1272414
42
LeeW. L.YenM. S.ChaoK. C.YuanC. C.NgH. T.ChaoH. 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
43
LegniniI.MorlandoM.MangiavacchiA.FaticaA.BozzoniI. (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
44
LiL.LiuB.WapinskiO.TsaiM. C.QuK.ZhangJ.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
45
LinD.FosterD. P.UngarL. H. (2011). VIF Regression: a fast regression algorithm for large data. J. Am. Stat. Assoc.106 (493), 232–247. doi: 10.1198/jasa.2011.tm10113
46
LiuC.ZhangY. H.DengQ.LiY.HuangT.ZhouS.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
47
LiuX.FengM.ZhengG.GuY.WangC.HeZ. (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
48
Martens-UzunovaE. S.BottcherR.CroceC. M.JensterG.VisakorpiT.CalinG. 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
49
MengJ.WangJ. (2015). Role of SNARE proteins in tumourigenesis and their potential as targets for novel anti-cancer therapeutics. Biochim. Biophys. Acta1856 (1), 1–12. doi: 10.1016/j.bbcan.2015.04.002
50
ModakJ. M.Roy-O’ReillyM.ZhuL.StaffI.McCulloughL. 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
51
MollR.DivoM.LangbeinL. (2008). The human keratins: biology and pathology. Histochem. Cell. Biol.129 (6), 705–733. doi: 10.1007/s00418-008-0435-6
52
NaydenovN. G.JoshiS.FeyginA.SainiS.LitovchickL.IvanovA. I. (2018). A membrane fusion protein, Ykt6, regulates epithelial cell migration via microRNA-mediated suppression of Junctional Adhesion Molecule A. Cell. Cycle17 (14), 1812–1831. doi: 10.1080/15384101.2018.1496755
53
NguyenC. D.MansfieldR. E.LeungW.VazP. M.LoughlinF. E.GrantR. 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
54
NitscheA.StadlerP. F. (2017). Evolutionary clues in lncRNAs. Wiley Interdiscip. Rev. RNA8 (1), 1–13. doi: 10.1002/wrna.1376
55
OrsoF.BalzacF.MarinoM.LemboA.RettaS. F.TavernaD. (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
56
ParkE. R.KimS. B.LeeJ. S.KimY. H.LeeD. H.ChoE. 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
57
ParrisT. Z.RönnermanE. W.EngqvistH.BiermannJ.TruvéK.NemesS.et al. (2018). Genome-wide multi-omics profiling of the 8p11-p12 amplicon in breast carcinoma. Oncotarget9 (35), 24140–24154. doi: 10.18632/oncotarget.25329
58
RamakrishnaM.WilliamsL. H.BoyleS. E.BearfootJ. L.SridharA.SpeedT. P.et al. (2010). Identification of candidate growth promoting genes in ovarian cancer through integrated copy number and expression analysis. PLoS One5 (4), e9983. doi: 10.1371/journal.pone.0009983.t001
59
RibeiroA. O.SchoofC. R.IzzottiA.PereiraL. V.VasquesL. R. (2014). MicroRNAs: modulators of cell identity, and their applications in tissue engineering. Microrna3 (1), 45–53. doi: 10.2174/2211536603666140522003539
60
RogersM. A.WinterH.LangbeinL.WollschlagerA.Praetzel-WunderS.Jave-SuarezL. 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
61
SongJ.XuQ.ZhangH.YinX.ZhuC.ZhaoK.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
62
StaffordJ. L.DysonG.LevinN. K.ChaudhryS.RosatiR.KalpageH.et al. (2017). Reanalysis of BRCA1/2 negative high risk ovarian cancer patients reveals novel germline risk loci and insights into missing heritability. PLoS One12 (6), e0178450. doi: 10.1371/journal.pone.0178450
63
SunQ.HuangX.ZhangQ.QuJ.ShenY.WangX.et al. (2016). SNAP23 promotes the malignant process of ovarian cancer. J. Ovarian Res.9 (1), 80. doi: 10.1186/s13048-016-0289-9
64
SundaramG. M.Veera BramhachariP. (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
65
The Cancer Genome Atlas Research NetworkWeinsteinJ. N.CollissonE. A.MillsG. B.ShawK. R.OzenbergerB. A.et al. (2013). The cancer genome atlas pan-cancer analysis project. Nat. Genet.45 (10), 1113–1120. doi: 10.1038/ng.2764
66
TorreL. A.TrabertB.DeSantisC. E.MillerK. D.SamimiG.RunowiczC. D.et al. (2018). Ovarian cancer statistics. CA. Cancer J. Clin.68 (4), 284–296 doi: 10.3322/caac.21456
67
TuJ.TianG.CheungH. H.WeiW.LeeT. 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
68
VetterM. H.HaysJ. 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
69
VyasS.ZaganjorE.HaigisM. C. (2016). Mitochondria and cancer. Biomed. Res. Int.166 (3), 555. doi: 10.1016/j.cell.2016.07.002
70
WakeN. C.RickettsC. J.MorrisM. R.PrigmoreE.GribbleS. M.SkytteA. 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
71
WangX.LiM.WangZ.HanS.TangX.GeY.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
72
WangX.ShiL.DengY.QuM.MaoS.XuL.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
73
WangY.XuZ.JiangJ.XuC.KangJ.XiaoL.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
74
WillisS.VillalobosV. M.GevaertO.AbramovitzM.WilliamsC.SikicB. I.et al. (2016). Single gene prognostic biomarkers in ovarian cancer: a meta-analysis. PLoS One11 (2), e0149183. doi: 10.1371/journal.pone.0149183
75
WymanS. K.ParkinR. K.MitchellP. S.FritzB. R.O’BriantK.GodwinA. K.et al. (2009). Repertoire of microRNAs in epithelial ovarian cancer as determined by next generation sequencing of small RNA cDNA libraries. PLoS One4 (4), e5311. doi: 10.1371/journal.pone.0005311
76
XingJ.YiJ.CaiX.TangH.LiuZ.ZhangX.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
77
XuC.ZhangY.WangQ.XuZ.JiangJ.GaoY.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
78
YamamotoK.OkamotoA.IsonishiS.OchiaiK.OhtakeY. (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
79
YamamotoN.KinoshitaT.NohataN.ItesakoT.YoshinoH.EnokidaH.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
80
YangC.CaiJ.WangQ.TangH.CaoJ.WuL.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
81
YangQ. Y.LiJ. H.WangQ. Y.WuY.QinJ. L.ChengJ. 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
82
YangY.HentatiA.DengH. X.DabbaghO.SasakiT.HiranoM.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
83
YeC.TaoR.CaoQ.ZhuD.WangY.WangJ.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
84
YoonJ. H.AbdelmohsenK.KimJ.YangX.MartindaleJ. L.Tominaga-YamanakaK.et al. (2013). Scaffold function of long non-coding RNA HOTAIR in protein ubiquitination. Nat. Commun.4, 2939. doi: 10.1038/ncomms3939
85
YoonT. Y.MunsonM. (2018). SNARE complex assembly and disassembly. Curr. Biol.28 (8), R397–R401. doi: 10.1016/j.cub.2018.01.005
86
YuanH.ChenV.BoisvertM.IsaacsC.GlazerR. I. (2018). PLAC1 as a serum biomarker for breast cancer. PLoS One13 (2), e0192106. doi: 10.1371/journal.pone.0192106
87
YuanL.LiS.ZhouQ.WangD.ZouD.ShuJ.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
88
ZhangH.ZhuC.YiZ.MingL.WuL.YangX.et al. (2015). Long non-coding RNA expression profiles of hepatitis C virus-related dysplasia and hepatocellular carcinoma. Oncotarget6 (41), 43770. doi: 10.18632/oncotarget.6087
89
ZhangH. W.WangQ. Y.ZhaoQ.DiW. (2013). MiR-124 inhibits the migration and invasion of ovarian cancer cells by targeting SphK1. J. Ovarian Res.6 (1), 1–9. doi: 10.1186/1757-2215-6-84
90
ZhaoQ.YangJ.LiL.CaoD.YuM.ShenK.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
91
ZhaoS.ChenQ.KangX.KongB.WangZ. (2018). Aberrantly expressed genes and miRNAs in slow transit constipation based on RNA-Seq analysis. Biomed. Res. Int.2018, 2617432. doi: 10.1155/2018/2617432
Summary
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
Volume
10 - 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
Updates
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.
*Correspondence: Jingting Jiang, jiangjingting@suda.edu.cn; Tao Huang, tohuangtao@126.com
This article was submitted to Bioinformatics and Computational Biology, a section of the journal Frontiers in Genetics
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.