BioMed Research International
Volume 2017 (2017), Article ID 3859582, 7 pages
https://doi.org/10.1155/2017/3859582

Research Article

Cancer-Related Triplets of mRNA-lncRNA-miRNA Revealed by Integrative Network in Uterine Corpus Endometrial Carcinoma

Chenglin Liu,¹ Yu-Hang Zhang,² Qinfang Deng,³ Yixue Li,1,4,5 Tao Huang,² Songwen Zhou,³ and Yu-Dong Cai⁶

¹School of Life Sciences and Biotechnology, Shanghai Jiaotong University, Shanghai 200240, China
²Institute of Health Sciences, Shanghai Institutes for Biological Sciences, Chinese Academy of Sciences, Shanghai 200031, China
³Department of Medical Oncology, Shanghai Pulmonary Hospital, Cancer Institute, Tongji University Medical School, Shanghai 200433, China
⁴Key Lab of Computational Biology, CAS-MPG Partner Institute for Computational Biology, Shanghai Institutes for Biological Sciences, Chinese Academy of Sciences, Shanghai 200031, China
⁵Collaborative Innovation Center for Genetics and Development, Fudan University, Shanghai 200433, China
⁶School of Life Sciences, Shanghai University, Shanghai 200444, China

Correspondence should be addressed to Tao Huang; tohuangtao@126.com, Songwen Zhou; zhou_songwen@126.com, and Yu-Dong Cai; cai_yud@126.com

Received 2 September 2016; Revised 28 September 2016; Accepted 22 November 2016; Published 8 February 2017

Academic Editor: Kazuhisa Nishizawa

Copyright © 2017 Chenglin Liu et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

The regulation of transcriptome expression level is a complex process involving multiple-level interactions among molecules such as protein coding RNA (mRNA), long noncoding RNA (lncRNA), and microRNA (miRNA), which are essential for the transcriptome stability and maintenance and regulation of body homeostasis. The availability of multilevel expression data enables a comprehensive view of the regulatory network. In this study, we analyzed the coding and noncoding gene expression profiles of 301 patients with uterine corpus endometrial carcinoma (UCEC). A new method was proposed to construct a genome-wide integrative network based on variance inflation factor (VIF) regression method. The cross-regulation relations of mRNA, lncRNA, and miRNA were then selected based on clique-searching algorithm from the network, when any two molecules of the three were shown as interacting according to the integrative network. Such relation, which we call the mRNA-lncRNA-miRNA triplet, demonstrated the complexity in transcriptome regulation process. Finally, six UCEC-related triplets were selected in which the mRNA participates in endometrial carcinoma pathway, such as CDH1 and TP53. The multi-type RNAs are proved to be cross-regulated as to each of the six triplets according to literature. All the triplets demonstrated the association with the initiation and progression of UCEC. Our method provides a comprehensive strategy for the investigation of transcriptome regulation mechanism.

1. Introduction

Uterine corpus endometrial carcinoma (UCEC) develops from the cells of the inner lining of the uterus, which is one of the most common female genital cancer threatening the health of women all over the world [1, 2]. Only counting 2012, approximately 320,000 women have been diagnosed and about 76,000 people have died of UCEC, according to incomplete statistics [3]. Most commonly, UCEC occurs in postmenopausal women, due to the unstable level of estrogen after menopause [4]. Smoking, high blood pressure, and being overweight also indirectly relate to uterus diseases via various regulation mechanisms [5–7]. In addition, genetic disorders also contribute to the development of UCEC and associate it with other diseases such as Lynch syndrome and colon cancer [8, 9]. A potential inherited tendency shows in UCEC with an increased risk in women with a family history of endometrial cancer [10]. Clinical diagnosis is according to the symptoms such as postmenopausal vaginal bleeding, enlarged uterus, low abdominal pain, and pelvic cramping [11–13]. The understanding of regulatory network could help to investigate its mechanism and benefit the diagnosis and treatment of UCEC.

The multi-type-molecular regulatory network, especially the interaction network between coding RNAs (mRNAs) and noncoding RNAs, has gained many interests in recent years. Previous reports have revealed that the microRNAs (miRNAs) which are small size noncoding RNAs of about 22 nucleotides and long noncoding RNAs (lncRNAs) which contain more than 200 nucleotides cross-regulate their expression levels and comodulate the expression of mRNAs. On the other hand, mRNAs also affect the expression of noncoding RNAs in specific ways [14, 15]. For example, the long intergenic noncoding RNA lincRNA-p21 has been reported to be downregulated by miRNA let-7. The binding of lincRNA-p21 to JUNB and CTNNB1 mRNAs results in the repression of JunB and -catenin translation [16]. Another experiment has shown that the depletion of lncRNA highly upregulated in liver cancer (HULC) results in significant deregulation of several genes involved in liver cancer. This lncRNA is upregulated by CREB mRNA which is underregulated by miR372 [17]. Such interaction, which we call the mRNA-lncRNA-miRNA triplet, is essential for the maintenance and regulation of body homeostasis. The aberrance of any of its molecules may influence the stability of multilevel expression and affect the tumorigenesis accordingly.

Recently, the availability of large scaled multilevel expression data provides an opportunity to obtain the comprehensive map of the multi-type-molecular regulatory network. The Cancer Genome Atlas (TCGA) database [18–20], especially the TCGA long noncoding RNAs website, provides the whole-genome profiling of 301 UCEC patients including the expression levels of mRNA, lncRNAs, and miRNAs. Such multidimensional resources allow us to investigate the mRNA-lncRNA-miRNA interactions, understand the transcriptional characteristic of UCEC, and dig deeper into the essential genetic alterations, transcriptional regulations, and posttranscriptional mechanisms throughout its initiation and progression [19].

Here, a new method is built to systematically investigate the mRNA-lncRNA-miRNA interactions in UCEC based on the patient expression profiles downloaded from TCGA long noncoding RNA website. An integrative network of mRNAs, lncRNAs, and miRNAs is constructed using an accurate and extremely efficient algorithm, the variance inflation factor (VIF) regression method. Many mRNA-lncRNA-miRNA triplets, which depict the cross-regulation relations among mRNA, lncRNA, and miRNA, are detected by searching all cliques (that is complete subgraphs with all vertices adjacent to each other) consisting of these three elements. The clique searching problem is a fundamental topic in computer science, which is very important in clustering analysis based on density and grid of data elements [21], and many solutions have been proposed to improve the searching performance. At last, the detected triplets are screened for their biological functions, and three of them are determined as UCEC-related triplets according to KEGG database and published literature. All in all, the proposed algorithm can find out disease-associated transcriptional RNA (miRNAs, lncRNAs, and mRNAs) interactions and may contribute to reveal the potential posttranscriptional regulatory mechanisms of UCEC.

2. Material and Methods

2.1. Datasets

The expression data of UCEC are obtained from TCGA long noncoding RNAs website (http://larssonlab.org/tcga-lncrnas/datasets.php), including the profiles of 20,462 protein coding genes, 10,419 lncRNA genes, and 742 miRNAs from 301 UCEC patients, as shown in Table 1. Specifically, the miRNA expression data are selected from the profiles of noncoding genes by their gene symbols. The expression levels are given as reads per kilobase per million (RPKM) values. The zero values of the expression data are set as the minimum nonzero RPKM of their corresponding sample for the allowance of log transformation. UCEC-related pathway information is adopted from KEGG (Kyoto Encyclopedia of Genes and Genomes, http://www.kegg.jp/) database with entry ID “hsa05213”. The pathway involves 52 genes including tumor protein coding gene TP53 and cadherin protein coding gene CDH1.

Table 1: Number of genes in raw data and constructed network of 301 UCEC patients.

mRNA
lncRNA
miRNA
Total

Raw data
20,462
10,419
742
31,623

Network
14,229
4,601
268
19,098

Triplet list
736
1,799
227
2,762

2.2. Integrative Network Construction

The interaction network is built by firstly determine the key factors (mRNA, lncRNA, and miRNA) affecting the expression level of each RNA molecule, respectively, and integrating them into a complete network after that. Hence, each RNA molecule is regarded as the dependent variable in one linear regression model, while all others are treated as the independent variables. In summary, 20,462 + 10,419 + 742 = 31,623 regression models are built where each one is based on 31,622 RNA expression features, and the integrative network is constructed after that.

Due to the large dataset in each regression model with far more features than observations (31,622 versus 301), an efficient regression and feature selection method, the variance inflation factor (VIF) regression algorithm [22] is utilized to select the optimal regulator set that is most related to each target RNA. The algorithm is designed to find the optimal that can minimize the penalized sum of squared errors,

where are observations and are predictors, , . Instead of searching over all subsets for the best , this algorithm evaluates the marginal correlations of each candidate predictors with the target factor using a small presampled set of data and searches the optimal subset by including -statistic correction procedure when adding or removing one variable at a time. The method has shown great efficiency but is also accurate compared to other methods such as LASSO and has comparable accuracy even compared to the most accurate but slowest regression method FoBa. The construction of the VIF regression models is based on program from http://cran.r-project.org/web/packages/VIF/.

Next, the goodness-of-fit for linear regression models is assessed by the adjust coefficient of determination (denoted as adjust ). The statistic measures how well the regression line approximates the real data points and can compare the regression model containing different number of regulators. In this paper, the regression model is retained only if it surpasses the adjust cutoff of 0.8. Regulation relations between the target RNA and its regulators are obtained from the retained regression models. These relations are further integrated as a comprehensive map for the mRNA-lncRNA-miRNA interaction. Note that the constructed network is an undirected graph. The edges are constructed if the two factors are connected by arcs of any direction. It is because the regression model can only identify the regulators of the target based on their gene expression associations but cannot determine if the regulators induce the perturbation of the target or vice versa without prior biological knowledge.

2.3. The mRNA-lncRNA-miRNA Triplet Detection

The detection of mRNA-lncRNA-miRNA triplets from the integrative network is a typical clique problem in computer science. Clique problem tries to search all complete subgraphs with all vertices connected to each other. Here, the size of subgraph is set as three, and the vertices of each subgraph are restricted to contain all of the three RNA types. The subgraphs, called mRNA-lncRNA-miRNA triplets, describe the relations of mRNA, lncRNA, and miRNA with each two of them coregulated according to the VIF regression model. The detection procedure is fulfilled by the “cliques” function in R package igraph [23].

Next, the UCEC-related triplets are further screened out if its mRNA participates in the hsa05213 pathway (endometrial cancer, Homo sapiens) according to KEGG database. These triplets are further analyzed for their interactions and biological functions as to UCEC according to literature.

3. Results and Discussion

3.1. Structure of the Integrative Network

The whole-genome integrative network of mRNA, lncRNA, and miRNA is constructed based on their cross-regulation relations using VIF regression. Totally, 19,098 factors are included in the integrative network, composed of 14229 coding mRNAs, 4,601 lncRNA, and 268 miRNA. On the network, each RNA is regulated by an average of 30 factors. The protein coding gene-gene interactions dominate the integrative network, as the expression levels of most coding genes are largely affected by only the coding mRNAs. Noncoding RNAs tend to have more interactions with noncoding RNAs instead of coding RNAs, which implies the extensive cross-talk of noncoding RNAs in their regulation of transcriptome and posttranscriptome. The details of the integrative network can be referred in Supplementary Material S1 available online at https://doi.org/10.1155/2017/3859582.

3.2. Candidate mRNA-lncRNA-miRNA Triplets of UCEC

By restricting the vertices types of clique problem to have all three RNA types, 14,416 mRNA-lncRNA-miRNA triplets are detected from the integrative network. These triplets involve 736 coding mRNAs, 1,799 lncRNAs, and 227 miRNAs, and provide a comprehensive map for the mRNA-lncRNA-miRNA interaction. The relatively small number of coding mRNAs compared to the noncoding RNAs indicates that many coding genes are coregulated by multiple lncRNAs and miRNAs. Extensive cross-talks exist in the regulatory process of noncoding RNAs, which also explain the complexity of transcriptome regulation process. The list of the detected triplets can be found in Supplementary Material S2.

Next, the triplet is considered as UCEC-related if its mRNA participates in hsa05213 endometrial cancer pathway. Note that the mRNA-lncRNA-miRNA cross-interaction is a very special interaction case that the genes in the selected triplet have little chance to be enriched in the pathway. However, studies have shown that the mutations in a pathway are mutual exclusive, and only one functional gene mutation is enough to perturb the pathway [24–26]. Hence, any triplet having overlapped genes with the pathway may contribute to the progression of cancer. Here, six triplets related to hsa05213 are detected and are retained for further analysis, as shown in Figure 1. The mRNA, lncRNA, and miRNA are labeled as red, green, and yellow, and the cross-regulation relations are shown as an undirected 3-vertex graph. Four of the six triplets, as shown in the first graph in Figure 1, involve the same mRNA and miRNA, but different lncRNAs, that is, mRNA CDH1-lncRNA (RP4-591L5.1, CTA.929C8.5.1, U47924.27.1, and AP006285.7.1)-miRNA miR128-1. The other two triplets are mRNA CDH1-lncRNA AP006285.7.1-miRNA miR126 and mRNA TP53-lncRNA CTD-2008N3.1.1-miRNA miR203, respectively.

Figure 1: UCEC-related triples screened out from the integrative network. The mRNA, lncRNA, and miRNA are labeled as red, green, and yellow. The cross-regulation relations are described as an undirected 3-vertex graph. The four triplets including the same mRNA and miRNA but different lncRNAs are presented in the top figure. The other two triplets are presented in the following.

3.3. Interaction and Biological Function of UCEC-Related Genes and Triplets

First, we focus on the mRNA-miRNA interaction and biological functions of the first set of triplets in Figure 1, which involves CDH1 and miR128-1, as mentioned above. CDH1 encodes a classical cadherin from cadherin superfamily, which is a calcium-dependent cell adhesion regulatory protein [27, 28]. CDH1 contributes the cell adhesion, mobility, and proliferation in specific microenvironment, especially in tumor [29]. As for UCEC, CDH1 contributes the initiation and invasion of endometrial cancer through its specific role in epithelial-mesenchymal transition (EMT) [30–32]. Additionally, miR128-1 has been proved to interact with the expression product of CDH1 cadherin and participate in the regulation of EMT in prostate cancer stem during the tumorigenesis [33, 34]. Apart from that, miR128-1 also participates in the regulation of progression and EMT in glioblastoma [35]. In fact, miR128-1 interacts with CDH1 coding protein cadherin via a specific upstream protein Bmi1 which is the direct target of miR128-1 [34, 35].

Next, we consider the four lncRNAs in the triplets. The first lncRNA RP4-591L5.1 is a crucial lncRNA which binds a specific miRNA miR218. MiR218 contributes the cellular chemosensitivity, migration, and invasion, which may further influence the cadherin regulation and associate with the function of miR128-1 [36, 37]. Another lncRNA, CTA.929C8.5.1, also called lnc-CRYBA4-7:1, has been predicted to be interacted with miR4268. MiR4268 is a rare miRNA with a special 3D structure. It has been proved to participate in the maintenance of stemness and may activate the initiation process of tumor in specific environment [38]. Additionally, as the stemness of cancer cells is associated with tumor migration and has specific relationship with the process of EMT [39, 40], the interaction of CTA.929C8.5.1 with miR4268 may affect the stemness of tumor cells and further have a specific influence on EMT, which explains its potential relationship with miR218 and gene CDH1 in the triplet. The third lncRNA U47924.27.1, also named lnc-PTPN6-1:1, is a unique lncRNA that is the target of several functional miRNAs, such as miR139, and may participate in the initiation of several tumors especially in hematopoietic malignancy [41]. Therefore, it is reasonable that, in UCEC, such lncRNA may play a similar way to interact with the miRNA and lncRNAs mentioned above and contribute to the tumor initiation and progression. Apart from that, U47924.27.1 is associated with PTEN, while CDH1 coding protein cadherin is also associated with PTEN cascades. Hence, this lncRNA may also have interaction with CDH1 [42]. The last lncRNA AP006285.7.1, also annotated as lnc-KRTAP5-4-1:1, has been proved to be associated with miR513a through sequence analysis. MiR513a may play its specific role in various cancer types and mainly regulate the proliferation of cells and contributes to the modeling of inflammation environment [43, 44]. Such regulatory functions may participate in EMT process and further promote the migration of tumor cells, which further explains the interaction of AP006285.7.1 with other factors in the triplet. In summary, all elements in the triplet contribute cohesively to the initiation and progression of UCEC and may exert influence on EMT process.

As to the triplet CDH1-lncRNA AP006285.7.1-miRNA miR126, we only investigate the functions of miR126 since the other two have been mentioned above. MiR126 has been proved as a predictive and diagnosis marker of esophageal cancer [45]. It is also associated with cell adhesion and migration and may contribute to the cadherin regulation in a similar way with other miRNAs (miR99a, miR200, etc.) [46, 47]. Therefore, miR126, lncRNA AP006285.7.1, and CDH1 can be clustered together because of their specific function and contribution to UCEC. This triplet focuses more on the cell adhesion instead of EMT progression and concentrates on the progression and migration process of the tumorigenesis of UCEC.

The next triplet is mRNA TP53-lncRNA CTD-2008N3.1.1-miRNA miR203. TP53 is the most famous tumor suppressor gene which generally contributes to every common type of tumor including endometrial cancer [48, 49]. As a multifunctional gene, TP53 also interacts with several crucial miRNAs (miR181, miR34a, miR520g, etc.) especially in various tumor tissues [50–53]. Consistent with our screen triplet, the interaction of TP53 and miR203 has been proved by several publications [54, 55]. Such interaction is quite crucial for certain kind of tumor especially for colon cancer [55]. As to lncRNA CTD-2008N3.1.1, which is also called lnc-CTD-2012 M11.2.1-1:1, it has been reported to associate with several miRNAs using computational prediction, which may have its specific way to interact with TP53 and miR203 in UCEC [56–58]. Additionally, CTD-2008N3.1.1 interacts with miR331 which participates in the tumorigenesis of various tumor types [56, 59–61]. Furthermore, CTD-2008N3.1.1 is a specific lncRNA originating from CTD sequence [62]. Since TP53 has been reported to be associated with several CTD structures in different tumor types, such screened lncRNA may also interact with TP53 and have its specific function in the process of UCEC initiation and progression [63, 64].

4. Conclusion

In summary, all of the selected triplets have been partially or fully confirmed to be associated with tumorigenesis especially in UCEC. Moreover, some of our preliminarily screened genes which are not included in these six triplets are also proved to be UCEC-associated and may have their specific function in the process of tumorigenesis. For example, miR204 is a specific miRNA in non-small cell lung cancer, which can specifically regulate the metastasis of tumor cells [65]. Such miRNA may have similar function in UCEC. Another miRNA, miR320, has been reported as a functional regulatory miRNA in stage I endometrioid endometrial carcinoma [66]. All in all, based on the expression profile of UCEC, our proposed method can cluster miRNAs, lncRNAs, and coding genes into functional interacted groups. Such an algorithm can also be applied to other cancer types and benefit the deeper understanding of the transcriptome regulatory mechanisms, and the cross-talk of multilevel RNAs such as miRNAs and lncRNAs. Additionally, the transcriptional level regulation network prediction helps to reveal the posttranscriptional regulation in tumors and other severe diseases.

Disclosure

Chenglin Liu, Yu-Hang Zhang, and Qinfang Deng are co-first authors.

Competing Interests

The authors declare that there is no conflict of interests regarding the publication of this manuscript.

Acknowledgments

This study was supported by the National Natural Science Foundation of China (Grant no. 81301994), the Natural Science Foundation of Shanghai Science and Technique Committee (Grant no. 13ZR1434700), Shanghai Sailing Program, and the Youth Innovation Promotion Association of Chinese Academy of Sciences (CAS) (Grant no. 2016245).

References

S. Lockwood, Contemporary Issues in Women's Cancers, Jones and Bartlett Publishers, Sudbury, Mass, USA, 2009.

A. Keen and E. Lennan, Women's Cancers, John Wiley & Sons, West Sussex, UK, 2011.

L. Ross, “2012 International Cancer Education Conference Proceedings,” Journal of Cancer Education, vol. 30, supplement 1, pp. 1–96, 2015.

Y. Chen, Q. Huang, Q. Chen et al., “The inflammation and estrogen metabolism impacts of polychlorinated biphenyls on endometrial cancer cells,” Toxicology in Vitro, vol. 29, no. 2, pp. 308–313, 2015.

M. Al-Zoughool, L. Dossus, R. Kaaks et al., “Risk of endometrial cancer in relationship to cigarette smoking: results from the EPIC study,” International Journal of Cancer, vol. 121, no. 12, pp. 2741–2747, 2007.

M. L. McCarroll, S. Armbruster, R. J. Pohle-Krauza et al., “Feasibility of a lifestyle intervention for overweight/obese endometrial and breast cancer survivors using an interactive mobile application,” Gynecologic Oncology, vol. 137, no. 3, pp. 508–515, 2015.

Y. Zhang, H. Liu, S. Yang, J. Zhang, L. Qian, and X. Chen, “Overweight, obesity and endometrial cancer risk: results from a systematic review and meta-analysis,” International Journal of Biological Markers, vol. 29, no. 1, pp. e21–e29, 2014.

J. T. Rabban, S. M. Calkins, A. N. Karnezis et al., “Association of tumor morphology with mismatch-repair protein status in older endometrial cancer patients: implications for universal versus selective screening strategies for lynch syndrome,” American Journal of Surgical Pathology, vol. 38, no. 6, pp. 793–800, 2014.

J. Y. Lee, H. J. Kim, E. H. Lee, H. W. Lee, J. Kim, and M. K. Kim, “One case of endometrial cancer occurrence: over 10 years after colon cancer in Lynch family,” Obstetrics & Gynecology Science, vol. 56, no. 6, pp. 408–411, 2013.

A. Dutt, H. B. Salvesen, H. Greulich, W. R. Sellers, R. Beroukhim, and M. Meyerson, “Somatic mutations are present in all members of the AKT family in endometrial carcinoma,” British Journal of Cancer, vol. 101, no. 7, pp. 1218–1219, 2009.

M. M. El Behery, M. H. Huda, E. M. Kamal, and A. E. Shehatta, “Diagnostic accuracy of uterine fluid lactate dehydrogenase isoenzyme activity profile and vaginal ultrasound in detecting endometrial cancer in women with postmenopausal bleeding,” Archives of Gynecology and Obstetrics, vol. 281, no. 4, pp. 717–721, 2010.

C. Neppe, R. Land, and A. Obermair, “Wrigley forceps to deliver a bulky uterus following a total laparoscopic hysterectomy for endometrial cancer,” Australian and New Zealand Journal of Obstetrics and Gynaecology, vol. 45, no. 5, pp. 444–445, 2005.

K. A. Donovan, A. R. Boyington, P. L. Judson, and J. F. Wyman, “Bladder and bowel symptoms in cervical and endometrial cancer survivors,” Psycho-Oncology, vol. 23, no. 6, pp. 672–678, 2014.

E. S. Martens-Uzunova, R. Böttcher, C. M. Croce, G. Jenster, T. Visakorpi, and G. A. Calin, “Long noncoding RNA in prostate, bladder, and kidney cancer,” European Urology, vol. 65, no. 6, pp. 1140–1151, 2014.

A. O. Ribeiro, C. R. G. Schoof, A. Izzotti, L. V. Pereira, and L. R. Vasques, “MicroRNAs: modulators of cell identity, and their applications in tissue engineering,” MicroRNA, vol. 3, no. 1, pp. 45–53, 2014.

J.-H. Yoon, K. Abdelmohsen, and M. Gorospe, “Functional interactions among microRNAs and long noncoding RNAs,” Seminars in Cell and Developmental Biology, vol. 34, pp. 9–14, 2014.

J.-H. Yoon, K. Abdelmohsen, S. Srikantan et al., “LincRNA-p21 suppresses target mRNA translation,” Molecular Cell, vol. 47, no. 4, pp. 648–655, 2012.

J. N. Weinstein, E. A. Collisson, G. B. Mills et al., “The Cancer Genome Atlas Pan-Cancer analysis project,” Nature Genetics, vol. 45, no. 10, pp. 1113–1120, 2013.

Cancer Genome Atlas Research Network, “Integrated genomic analyses of ovarian carcinoma,” Nature, vol. 474, no. 7353, pp. 609–615, 2011.

J. N. Weinstein, E. A. Collisson, G. B. Mills et al., “The cancer genome atlas pan-cancer analysis project,” Nature Genetics, vol. 45, no. 10, pp. 1113–1120, 2013.

S. Luccioli, E. Ben-Jacob, A. Barzilai, P. Bonifazi, and A. Torcini, “Clique of functional hubs orchestrates population bursts in developmentally regulated neural networks,” PLoS Computational Biology, vol. 10, no. 9, Article ID e1003823, 2014.

D. Lin, D. P. Foster, and L. H. Ungar, “VIF regression: a fast regression algorithm for large data,” Journal of the American Statistical Association, vol. 106, no. 493, pp. 232–247, 2011.

G. Csardi and T. Nepusz, “The igraph software package for complex network research,” InterJournal Complex Systems, Article ID 1695, 2006.

J. Zhao, S. Zhang, L.-Y. Wu, and X.-S. Zhang, “Efficient methods for identifying mutated driver pathways in cancer,” Bioinformatics, vol. 28, no. 22, pp. 2940–2947, 2012.

B. Vogelstein and K. W. Kinzler, “Cancer genes and the pathways they control,” Nature Medicine, vol. 10, no. 8, pp. 789–799, 2004.

C.-H. Yeang, F. McCormick, and A. J. Levine, “Combinatorial patterns of somatic gene mutations in cancer,” The FASEB Journal, vol. 22, no. 8, pp. 2605–2622, 2008.

M. Silies and C. Klämbt, “APC/CFzr/Cdh1-dependent regulation of cell adhesion controls glial migration in the Drosophila PNS,” Nature Neuroscience, vol. 13, no. 11, pp. 1357–1364, 2010.

Q.-W. Deng, B.-S. He, Y.-Q. Pan et al., “Roles of E-cadherin (CDH1) genetic variations in cancer risk: a meta-analysis,” Asian Pacific Journal of Cancer Prevention, vol. 15, no. 8, pp. 3705–3713, 2014.

Z. R. Qian, T. Sano, K. Yoshimoto et al., “Tumor-specific downregulation and methylation of the CDH13 (H-cadherin) and CDH1 (E-cadherin) genes correlate with aggressiveness of human pituitary adenomas,” Modern Pathology, vol. 20, no. 12, pp. 1269–1277, 2007.

E. Płuciennik, M. Nowakowska, K. Pospiech et al., “The role of WWOX tumor suppressor gene in the regulation of EMT process via regulation of CDH1-ZEB1-VIM expression in endometrial cancer,” International Journal of Oncology, vol. 46, no. 6, pp. 2639–2648, 2015.

S.-M. Hsiao, M.-W. Chen, C.-A. Chen et al., “The H3K9 methyltransferase G9a represses E-cadherin and is associated with myometrial invasion in endometrial cancer,” Annals of Surgical Oncology, vol. 22, pp. 1556–1565, 2015.

K. Banno, M. Yanokura, M. Iida, K. Masuda, and D. Aoki, “Carcinogenic mechanisms of endometrial cancer: involvement of genetics and epigenetics,” Journal of Obstetrics and Gynaecology Research, vol. 40, no. 8, pp. 1957–1967, 2014.

Y. Zhang, T. Chao, R. Li et al., “MicroRNA-128 inhibits glioma cells proliferation by targeting transcription factor E2F3a,” Journal of Molecular Medicine, vol. 87, no. 1, pp. 43–51, 2009.

R. Nanta, D. Kumar, D. Meeker et al., “NVP-LDE-225 (Erismodegib) inhibits epithelial-mesenchymal transition and human prostate cancer stem cell growth in NOD/SCID IL2Rγ null mice by regulating Bmi-1 and microRNA-128,” Oncogenesis, vol. 2, article no. e42, 2013.

J. Fu, M. Rodova, R. Nanta et al., “NPV-LDE-225 (Erismodegib) inhibits epithelial mesenchymal transition and self-renewal of glioblastoma initiating cells by regulating miR-21, miR-128, and miR-200,” Neuro-Oncology, vol. 15, no. 6, pp. 691–706, 2013.

R. Dong, H. Qiu, G. Du, Y. Wang, J. Yu, and C. Mao, “Restoration of microRNA-218 increases cellular chemosensitivity to cervical cancer by inhibiting cell-cycle progression,” Molecular Medicine Reports, vol. 10, no. 6, pp. 3289–3295, 2014.

B. Peng, D. Li, M. Qin et al., “MicroRNA218 inhibits glioma migration and invasion via inhibiting glioma-associated oncogene homolog 1 expression at N terminus,” Tumor Biology, vol. 35, no. 4, pp. 3831–3837, 2014.

L. A. Goff, J. Davila, M. R. Swerdel et al., “Ago2 immunoprecipitation identifies predicted MicroRNAs in human embryonic stem cells and neural precursors,” PLoS ONE, vol. 4, no. 9, Article ID e7192, 2009.

T. Saito and K. Mimori, “Cancer stemness and circulating tumor cells,” Nihon rinsho. Japanese journal of clinical medicine, vol. 73, no. 5, pp. 806–810, 2015.

A. W. Fender, J. M. Nutter, T. L. Fitzgerald, F. E. Bertrand, and G. Sigounas, “Notch-1 promotes stemness and epithelial to mesenchymal transition in colorectal cancer,” Journal of Cellular Biochemistry, vol. 116, no. 11, pp. 2517–2527, 2015.

M. F. Alemdehy, J. R. Haanstra, H. W. J. de Looper et al., “Interstrand cross-link induced miR139-3p and miR199a-3p have opposite roles in hematopoietic cell expansion and leukemic transformation,” Blood, vol. 125, no. 25, pp. 3937–3948, 2015.

Y. Chen, Y. Sun, L. Chen et al., “MiRNA-200c increases the sensitivity of breast cancer cells to doxorubicin through the suppression of E-cadherin-mediated PTEN/Akt signaling,” Molecular Medicine Reports, vol. 7, no. 5, pp. 1579–1584, 2013.

X. Chen, G. Zhao, F. Wang et al., “Upregulation of miR-513b inhibits cell proliferation, migration, and promotes apoptosis by targeting high mobility group-box 3 protein in gastric cancer,” Tumor Biology, vol. 35, no. 11, pp. 11081–11089, 2014.

A.-Y. Gong, R. Zhou, G. Hu et al., “MicroRNA-513 regulates B7-H1 translation and is involved in IFN-γ-induced B7-H1 expression in cholangiocytes,” Journal of Immunology, vol. 182, no. 3, pp. 1325–1333, 2009.

N. S. Sakai, E. Samia-Aly, M. Barbera, and R. C. Fitzgerald, “A review of the current understanding and clinical utility of miRNAs in esophageal cancer,” Seminars in Cancer Biology, vol. 23, no. 6, pp. 512–521, 2013.

D. Li, X. Li, W. Cao, Y. Qi, and X. Yang, “Antagonism of microRNA-99a promotes cell invasion and down-regulates E-cadherin expression in pancreatic cancer cells by regulating mammalian target of rapamycin,” Acta Histochemica, vol. 116, no. 5, pp. 723–729, 2014.

L. Romero-Pérez, M. Á. López-García, J. Díaz-Martín et al., “ZEB1 overexpression associated with E-cadherin and microRNA-200 downregulation is characteristic of undifferentiated endometrial carcinoma,” Modern Pathology, vol. 26, no. 11, pp. 1514–1524, 2013.

B. Leroy, M. Anderson, and T. Soussi, “TP53 mutations in human cancer: database reassessment and prospects for the next decade,” Human Mutation, vol. 35, no. 6, pp. 672–688, 2014.

M. Chmelarova, S. Kos, E. Dvorakova et al., “Importance of promoter methylation of GATA4 and TP53 genes in endometrioid carcinoma of endometrium,” Clinical Chemistry and Laboratory Medicine, vol. 52, no. 8, pp. 1229–1234, 2014.

F. Ganci, A. Sacconi, N. B. Ben-Moshe et al., “Expression of TP53 mutation-associated microRNAs predicts clinical outcome in head and neck squamous cell carcinoma patients,” Annals of Oncology, vol. 24, no. 12, pp. 3082–3088, 2013.

Y. K. Cheah, R. W. Cheng, S. K. Yeap, C. H. Khoo, and H. S. See, “Analysis of TP53 gene expression and p53 level of human hypopharyngeal FaDu (HTB-43) head and neck cancer cell line after microRNA-181a inhibition,” Genetics and Molecular Research, vol. 13, no. 1, pp. 1679–1683, 2014.

A. Dufour, G. Palermo, E. Zellmeier et al., “Inactivation of TP53 correlates with disease progression and low miR-34a expression in previously treated chronic lymphocytic leukemia patients,” Blood, vol. 121, no. 18, pp. 3650–3657, 2013.

Y. Zhang, L. Geng, G. Talmon, and J. Wang, “MicroRNA-520g confers drug resistance by regulating p21 expression in colorectal cancer,” Journal of Biological Chemistry, vol. 290, no. 10, pp. 6215–6225, 2015.

D. J. McKenna, S. S. McDade, D. Patel, and D. J. McCance, “MicroRNA 203 expression in keratinocytes is dependent on regulation of p53 levels by E6,” Journal of Virology, vol. 84, no. 20, pp. 10644–10652, 2010.

J. A. Li, Y. Chen, J. Zhao, F. Kong, and Y. Zhang, “miR-203 reverses chemoresistance in p53-mutated colon cancer cells through downregulation of Akt2 expression,” Cancer Letters, vol. 304, no. 1, pp. 52–59, 2011.

N. Léveillé, C. A. Melo, K. Rooijers et al., “Genome-wide profiling of p53-regulated enhancer RNAs uncovers a subset of enhancers controlled by a lncRNA,” Nature Communications, vol. 6, article no. 6520, 2015.

Y. Sánchez, V. Segura, O. Marín-Béjar et al., “Genome-wide analysis of the human p53 transcriptional network unveils a lncRNA tumour suppressor signature,” Nature communications, vol. 5, article 5812, p. 5812, 2014.

A. Zhang, M. Xu, and Y.-Y. Mo, “Role of the lncRNA-p53 regulatory network in cancer,” Journal of Molecular Cell Biology, vol. 6, no. 3, pp. 181–191, 2014.

X.-H. Liu, M. Sun, F.-Q. Nie et al., “Lnc RNA HOTAIR functions as a competing endogenous RNA to regulate HER2 expression by sponging miR-331-3p in gastric cancer,” Molecular Cancer, vol. 13, no. 1, article 92, 2014.

M. R. Epis, K. M. Giles, P. A. Candy, R. J. Webster, and P. J. Leedman, “MiR-331-3p regulates expression of neuropilin-2 in glioblastoma,” Journal of Neuro-Oncology, vol. 116, no. 1, pp. 67–75, 2014.

M. R. Epis, A. Barker, K. M. Giles, D. J. Beveridge, and P. J. Leedman, “The RNA-binding protein HuR opposes the repression of ERBB-2 gene expression by microRNA miR-331-3p in prostate cancer cells,” The Journal of Biological Chemistry, vol. 286, no. 48, pp. 41442–41454, 2011.

J.-H. Yoon, K. Abdelmohsen, and M. Gorospe, “Posttranscriptional gene regulation by long noncoding RNA,” Journal of Molecular Biology, vol. 425, no. 19, pp. 3723–3730, 2013.

T. Terakawa, H. Kenzaki, and S. Takada, “P53 searches on DNA by rotation-uncoupled sliding at C-terminal tails and restricted hopping of core domains,” Journal of the American Chemical Society, vol. 134, no. 35, pp. 14555–14562, 2012.

S. K. Mungamuri, S. Wang, J. J. Manfredi, W. Gu, and S. A. Aaronson, “Ash2L enables P53-dependent apoptosis by favoring stable transcription pre-initiation complex formation on its pro-apoptotic target promoters,” Oncogene, vol. 34, no. 19, pp. 2461–2470, 2015.

L. Shi, B. Zhang, X. Sun et al., “MiR-204 inhibits human NSCLC metastasis through suppression of NUAK1,” British Journal of Cancer, vol. 111, no. 12, pp. 2316–2327, 2014.

H. Xiong, Q. Li, S. Liu et al., “Integrated microRNA and mRNA transcriptome sequencing reveals the potential roles of miRNAs in stage I endometrioid endometrial carcinoma,” PLOS ONE, vol. 9, no. 10, Article ID e110163, 2014.