Abstract
Cancers of the digestive system are malignant diseases. Our study focused on colon cancer, esophageal cancer (ESCC), rectal cancer, gastric cancer (GC), and rectosigmoid junction cancer to identify possible biomarkers for these diseases. The transcriptome data were downloaded from the TCGA database (The Cancer Genome Atlas Program), and a network was constructed using the WGCNA algorithm. Two significant modules were found, and coexpression networks were constructed. CytoHubba was used to identify hub genes of the two networks. GO analysis suggested that the network genes were involved in metabolic processes, biological regulation, and membrane and protein binding. KEGG analysis indicated that the significant pathways were the calcium signaling pathway, fatty acid biosynthesis, and pathways in cancer and insulin resistance. Some of the most significant hub genes were hsa-let-7b-3p, hsa-miR-378a-5p, hsa-miR-26a-5p, hsa-miR-382-5p, and hsa-miR-29b-2-5p and SECISBP2 L, NCOA1, HERC1, HIPK3, and MBNL1, respectively. These genes were predicted to be associated with the tumor prognostic reference for this patient population.
Introduction
Cancer is one of the most common and lethal disease types in humans, and the molecular mechanisms governing cancer progression have not been elucidated to date (; Yang et al., 2020; Zhang Z.M. et al., 2020). Tumors of the digestive system are a serious problem for human health, and esophageal tumors, colorectal tumors, gastric tumors and rectosigmoid-junction tumors are among the ten most common cancers worldwide. These cancers cause the most deaths, and the incidence rate is increasing year by year (, , ; ; ; ; ; Zeng et al., 2020a; Zhang L. et al., 2020). Therefore, the study of digestive system tumors has important practical significance.
Digestive system tumors are a common malignancy in the clinic. Colorectal cancer is the combination of colon and rectal cancer cases (; ). Colorectal cancer ranks third in men and second in women in the rankings of cancer incidence and is commonly diagnosed despite the increased death risk worldwide (Torre et al., 2015; ; ). Colorectal cancer (CRC) is derived from the accumulation of epigenetic and genetic changes that lead to cancer-related deaths worldwide. With significant advancements in the early diagnosis of colorectal cancer, the cancer mortality rate has increased in the United States (; ; ). In European countries, colorectal cancer also has a high incidence, and the annual rate is increasing worldwide (; ). Although the cancer-related death rate increases every year, treatments can be administered at the early stages of cancer (; Zhang Y. et al., 2019). Other digestive system tumors, such as esophageal cancer (ESCC), have been reported as the eighth most common cancer worldwide with a high cancer-related death rate (; ). According to previous studies, almost 400,000 patients die from this disease among 450,000 diagnosed patients, the cancer-related death rate is approximately equal to 87.8% worldwide, and the 5-year survival rate is as low as approximately 25% (, ). The most common cancer in Asia is gastric cancer (GC), which is the fifth ranking cancer leading to the third highest cancer-related death rate in the world. At stages I and IIA, the GC 5-year survival rate is 81.8–93.6%, but at stage IIIC, it is very low at only 17.9% (; ).
Recently, network analysis has become a popular method for large-scale data (; Zeng X. et al., 2019; ; ; Zeng et al., 2020b), and the WGCNA is a popular method for these coexpression network constructs (; Zeng W. et al., 2019; ). Five steps were included in this method: 1. a gene coexpression network was constructed by an adjacency matrix; 2. network modules were identified by the hierarchical cluster method; 3. modules with related phenotypes were analyzed; 4. module relationships were analyzed; and 5. key modules were found. The WGCNA method is widely used in biological networks and genomic data mining analysis, and also gives rise to network-based meta-analysis techniques ().
The theoretical basis of the ceRNA (competitive endogenous RNA) hypothesis is that lncRNAs interact directly to regulate the expression of targeted genes and indirectly combine with miRNA sequences via the common miRNA response elements (MREs) of lncRNAs (). Cytoplasmic lncRNAs mainly affect mRNA reliability and translation mechanisms by binding miRNAs in the ceRNA regulatory network (). Recently, analysis of lncRNA, miRNA, and mRNA networks has been identified to be involved in the progression of cancer (Zhang X. et al., 2017; Zeng et al., 2018), including endogenous cancer, hepatocellular carcinoma cancer and other malignant tumors, using the WGCNA method (Wang et al., 2010; Tay et al., 2014; Wei et al., 2014, 2018, 2019a,b; Thomson and Dinger, 2016; , ; ; ; ; ; ). To date, there are few studies on digestive system malignant tumors. Therefore, integrated analysis of the regulatory functions of digestive system malignant tumors for lncRNA-miRNA-mRNA interaction networks requires large sample data and methods (; Zeng et al., 2016; Zou et al., 2016; ; ; ).
Competitive endogenous RNAs (ceRNAs) represent a novel gene expression regulation model and have attracted much attention from the academic community in recent years (Zhao et al., 2020). Compared with the miRNA regulation network, the ceRNA regulation network is more sophisticated and complex because the ceRNA regulation network includes more RNA molecules, such as mRNAs, pseudogenes, long non-coding RNAs (lncRNAs) and miRNAs. The ceRNA network provides a new perspective for transcriptomic and biological research.
In our studies, we used RNA-seq data from TCGA for free-scale gene coexpression network construction of digestive system malignant tumors. This study was able to provide novel biomarkers for these malignant tumors. Therefore, potential biomarkers were identified by using bioinformatics methods for integrative analysis based on the large amount of RNA-seq. According to the results of the analysis of malignant tumors, NCOA1, HERC1, HIPK3, MBNL1, hsa-let-7b-3p, hsa-miR-378a-5p, and hsa-miR-26a-5p were predictive or prognostic factors for malignant tumors.
Data and Methods
Date Collection
For data analysis in this study, gdc-client 1.5.0 ww1 and TCGAbiolinks 2.16.3 tools2 () were used from the TCGA database. A total of 437 samples of mRNA data, 388 samples of miRNA data and 385 samples of clinical data for colon cancer were obtained from TCGA, including colon cancer samples and normal samples. For colon cancer, there were over 500 individuals, including 480 tumor and 41 normal samples. Eighty-seven samples of mRNA data, 97 samples of miRNA data, and 87 samples of clinical data for ESCC were obtained from TCGA, including esophageal cancer (ESCC) samples and normal samples. There were 90 samples of mRNA data, 79 samples of miRNA data and 87 samples of clinical data for rectal cancer obtained from TCGA, including ESCC samples and normal samples. A total of 373 samples of mRNA data, 452 samples of miRNA data and 406 samples of clinical data for gastric cancer (GC) were obtained from TCGA, including ESCC samples and normal samples. Sixty-six samples of mRNA data, 63 samples of miRNA data and 63 samples of clinical data for rectosigmoid junction cancer were obtained from TCGA, including ESCC samples and normal samples. The guidelines of this study are from the TCGA website3.
Computational Analysis of RNA-Seq Data
The level three data of these cancers from TCGA by Illumina HiSeq 2000 platform (Illumina Inc., San Diego, California, United States) including miRNA, lncRNA, and mRNA data were analyzed by Digital Gene Expression Data package (edgeR 3.30.3)4 and limma (limma 3.44.3)5 package in R. The genes were annotated gene symbols through the Ensembl database6 (Zerbino et al., 2018). According to these analyses, absolute log2(fold change) ≥ 2.0 and FDR ≤ 0.01 were set to screen the significant differentially expressed genes of mRNAs (DEGs), lncRNAs (DElncRNAs) and miRNAs (DEmiRNAs) between tumor tissues and normal tissues. Non-cancer-specific expression genes were filtered, and upregulated or downregulated genes were saved. The package ggplot2 3.3.27 was used in RStudio to generate volcano plots of upregulated or downregulated genes. Meanwhile, the VennDiagram R package8 was used for the Venn diagram of the datasets.
Construction of the Coexpression Network
The WGCNA package was used to analyze the data as described previously (). For the different conditions of tumor tissues, some samples may deviate from the actual situation. These tissue samples will affect the accuracy of the results. Therefore, these samples were removed before sample analysis. The function of goodSamplesGenes was used to verify if there were many missing values for each sample and to remove these samples from the total data. Then, the function hculst was used for cluster analysis of tissue samples. Appropriate soft thresholds need to be selected to build the network when the WGCNA method is performed. The soft threshold value should not only reach a scale-free fitting index of more than 0.9 but also requires greater mean connectivity of the network. Thus, 4 was selected for the soft threshold value. A good soft threshold is selected based on this rule. At the same time, the size of the gene module will be divided, and at least 30 genes will be identified for one module. If the correlation coefficient is greater than 0.75 for two modules, then these two modules will be merged into one module. To enhance the productivity of the modules, a cutoff (< 0.25) was selected to combine similar modules. After these modules were divided, the sample traits were associated with each module such that the most important associated modules for the traits were searched. Cytoscape 3.8.1 was used for the visualization of gene coexpression.
Gene Ontology and Pathway Enrichment Analysis
The DAVID online tool was applied for the GO (Gene Ontology) and KEGG analysis and functional annotation of DEmRNAs. The DAVID9 database website was used for annotation and GO enrichment. KEGG pathways and GO terms enriched adopt adjusted P-value ≤ 0.05. WebGestalt10 is used for enrichment analysis that supports three well-established and complementary methods for enrichment analysis. The genes of the network were analyzed by the WebGestalt website. The R package of GOplot (version 1.0.2) was used for the GO analysis of the DEmRNAs of the five cancers.
Hub Gene Identification and Validation
To identify the gene connectivity, Pearson’s correlation was used for the test. In general, hub genes existed in modules closely linked to traits (cor.geneTraitSignificance > 0.2), the modules with hub genes had high connectivity (cor.geneModuleMembership > 0.8), and the packages of cytoHubba and MCODE in Cytoscape 3.8.1 were used to search the hub genes.
Results
Identification and Validation of Differentially Expressed mRNAs, miRNAs, and lncRNAs
To identify the differentially expressed mRNAs, lncRNAs and miRNAs, the standard of log2(fold change) ≥ 1 and q-value < 0.05 was used, and the EdgeR package of R (; ; ) was adopted to calculate separately for the normal samples and cancer samples of five cancers. For colon cancer, 4982 DEmRNAs, 149 DEmiRNAs and 1636 DElncRNAs were identified. There were 1409 differentially expressed genes in lncRNAs, 78 DEmiRNAs and 3659 differentially expressed mRNAs genes found between healthy and cancer-treated samples for esophageal cancer. A total of 3337 DEmRNAs, 168 DEmiRNAs and 871 DElncRNAs were identified between the rectosigmoid junction cancer tissues and matched normal control tissues; in addition, 2438 DEmRNAs, 156 DEmiRNAs and 607 DElncRNAs were identified for rectal cancer, and 4325 DEmRNAs, 62 DEmiRNAs and 2120 DElncRNAs were identified for gastric cancer compared with normal control samples. The distributions of differentially expressed DEmRNAs, DEmiRNAs and DElncRNAs were identified through volcano plots (Figure 1). A number of differences were found, and 87 DElncRNAs, 6 DEmiRNAs and 602 DEmRNAs were shared among the five cancers as shown in Figure 2.
FIGURE 1
FIGURE 2
GO and KEGG Pathway Analysis of DEmRNAs
GO and KEGG methods were adopted to investigate the annotation of DEmRNAs for the five cancers (Figures 3, 4). Standards of log2(fold change) ≥ 3 and q-value < 0.01 were used for the analysis.
FIGURE 3
FIGURE 4
A total of 766 DEmRNAs were primarily enriched in biological processes (BP) such as organic anion transport for colon cancer (GO:0015711) and cell components (CC) such as the apical part of the cell (GO:0045177) and the apical plasma membrane (GO:0016324) as well as ligand activity (GO:0048018) and signaling receptor activator activity (GO:0030546) of molecular function (MF) receptor for colon cancer (Figure 3). In addition, KEGG enrichment analysis revealed that the significant pathways were “bile secretion,” “drug metabolism - cytochrome P450,” “chemical carcinogenesis,” “metabolism of xenobiotics by cytochrome P450” and “retinol metabolism,” which are associated with the progression of colon cancer (Figure 4).
To investigate the enrichment in esophageal cancer, 448 DEmRNAs were used and enriched in muscle contraction (GO:0006936), the muscle system process (GO:0003012), digestion (GO:0007586), the presynaptic membrane (GO:0042734), the intrinsic component of the presynaptic membrane (GO:0098889) of biological processes (BP) and cell components (CC) (Figure 3). In addition, “gastric acid secretion,” “neuroactive ligand-receptor interaction,” “protein digestion and absorption,” “IL-17 signaling pathway,” and “calcium signaling pathway” were the top 5 pathways of the KEGG analysis (Figure 4).
Rectosigmoid junction cancer was also analyzed by GO and KEGG enrichment; 623 DEmRNAs were screened out, and the top 5 GO enrichments were the regulation of membrane potential (GO:0042391), the transmembrane transporter complex (GO:1902495), the transporter complex (GO:1990351), muscle contraction (GO:0006936), and the ion channel complex (GO:0034702) of biological processes (BP) and cell components (CC) (Figure 3). The top 5 KEGG pathways were “neuroactive ligand-receptor interaction,” “nicotine addiction,” “vascular smooth muscle contraction,” “pancreatic secretion,” and “cAMP signaling pathway” (Figure 4).
Rectal cancer was also analyzed by GO and KEGG enrichment; 691 DEmRNAs were screened out, and the top 5 GO enrichments were the regulation of membrane potential (GO:0042391), the regulation of transsynaptic signaling (GO:0099177), the synaptic membrane (GO:0097060), the modulation of chemical synaptic transmission (GO:0050804), and the muscle system process (GO:0003012) of biological processes (BP) and cell components (CC) (Figure 3). The top 5 KEGG pathways were “neuroactive ligand-receptor interaction,” “circadian entrainment,” “cGMP-PKG signaling pathway,” “glutamatergic synapse” and “cAMP signaling pathway” (Figure 4).
Gastric cancer was also subjected to GO and KEGG enrichment analyses. A total of 637 DEmRNAs were screened out, and the top 5 GO enrichments were cornification (GO:0070268), digestion (GO:0007586), antimicrobial humoral response (GO:0019730) (; ), keratinocyte differentiation (GO:0030216), and skin development (GO:0043588) of biological processes (BP) (Figure 3). The top 5 KEGG pathways were “Protein digestion and absorption,” “Bile secretion,” “Fat digestion and absorption,” “Chemical carcinogenesis,” and “Metabolism of xenobiotics by cytochrome P450” (Figure 4).
Construction of the WGCNA Network
To explore further the key biological genes related to cancers, the WGCNA method was applied to select the DEmRNAs, DEmiRNAs and DElncRNAs. For further analysis, a soft threshold (β = 4) was adopted to guarantee a scale-free network with high scale independence and low mean connectivity (near 0) (Figure 5A). DEGs from the five cancers were divided into several modules by cluster analysis (Figure 5B). Approximately 43 modules were generated for the five cancers. The module trait relationship is shown in Figure 5C. The dark green module related to the colon, rectosigmoid junction and rectal tumors were the deepest (cor = 0.21, 0.16, and 0.52 P = 4E-11, 6e-07, and 4e-67). For the other two cancers, the esophageal and gastric cancers, the brown module was chosen according to the correlation (cor = 0.46 and 0.46 P = 4E-52, 3e-52). Dark green and brown modules were selected for further analysis (Figure 5C). The high correlation and high P-value indicated that these modules are suitable for hub gene identification in these cancers.
FIGURE 5
Identification of the Hub Genes in the Five Cancers
The dark green module was selected for the edges signifying the correlations in colon cancer, rectosigmoid junction cancer and rectal cancer with the filter criterion of a weight value greater than 0.02 by the WGCNA algorithm. A total of 629 edges and 61 nodes were obtained and input into Cytoscape 3.8.1 (Figure 6A). A network was generated by Cytoscape software, and the hub genes were obtained from the network by applying CytoHubba. The most significant hub gene network was discovered by CytoHubba, as shown in Figure 6B. The identified hub miRNAs included hsa-let-7b-3p, hsa-miR-378a-5p, hsa-miR-26a-5p, hsa-miR-382-5p, and hsa-miR-29b-2-5p. In addition, the brown module was selected for the edges signifying the correlations in esophageal cancer and gastric cancer with the filter criterion of a weight value greater than 0.12 by the WGCNA algorithm. A total of 1027 edges and 279 nodes were obtained and input into Cytoscape 3.8.1 (Figure 6C). The most significant hub genes were discovered by CytoHubba as shown in Figure 6D. The identified hub genes were SECISBP2 L, NCOA1, HERC1, HIPK3, and MBNL1.
FIGURE 6
We discover that most of the genes in dark green module that were significantly correlated with colon cancer, rectosigmoid junction cancer and rectal cancer were miRNAs. This dominance of miRNAs in dark green module suggests the possibility that lncRNAs play a significant role in these cancers through regulation of coding genes in key pathways. miRNA can act by guiding histone modifiers and chromatin modifiers to regulate transcription and play crucial roles in cell differentiation that ultimately determine cell fate.
Functional Annotation of the Module of Interest and the Network Genes
GO and KEGG analyses of the two modules of dark green and brown with 61 and 279 genes, respectively, in WebGestalt are shown in Figure 7. These were involved in molecular and functional biological processes and cellular components. For the dark green module, the genes were associated with metabolic process, membrane and protein binding, membrane among others (Figure 7A), and the genes in the brown module were enriched in biological regulation, membrane, and protein binding among others (Figure 7B). Analysis of KEGG pathways for the dark green module, the calcium signaling pathway, fatty acid biosynthesis and the chemokine signaling pathway were significant among all the pathways (Figure 7C), and pathways in cancer, insulin resistance and renal cell carcinoma were significant pathways in the brown module (Figure 7D).
FIGURE 7
Discussion
Cancer of the digestive system is a common cancer worldwide. Recently, RNAs have been reported to be related to cancers. The WGCNA method was applied to analyze the related modules with different tumors and to identify the important genes of these tumors.
ESCC remains a serious burden on the health system worldwide, and some research has suggested that RNAs are associated with ESCC. For example, 2,046 circRNAs were frequently altered in ESCC tissues (). Several induced circRNAs were identified in radioresistant ESCC cells compared to normal ESCC cells (). Eight lncRNAs could be practical and reliable prognostic tools for esophageal cancers (). Additionally, the special functions of RNAs have been clearly demonstrated. For the analysis in gastric cancer (GC), 9 miRNAs, 41 lncRNAs, and 10 mRNAs were suggested and significant for GC (). In addition, 15 lncRNAs were identified to be significantly related to the clinical features of colon cancer (). In rectal cancer, COL1A1 and MZB1, as the key genes of rectal cancer, can interact with other genes correlated with shorter survival for patients (Wu et al., 2019). However, these studies presented research on only one cancer for the digestion system. Therefore, there is still an urgent need to identify the hub genes in digestive system cancers.
In this study, data from the TCGA website were adopted, and DEmRNAs, 168 DEmiRNAs, and 871 DElncRNAs were identified for the different tumors with the RNA-sequencing method and integrated bioinformatics analysis. Based on the WGCNA algorithm, the modules associated with different tumors were generated under the appropriate thresholds. Therefore, a network related to the different cancers was constructed in this study. Based on the network analysis, two modules, dark green and brown, were selected for further research. CytoHubba from Cytoscape 3.8.1 was used for this analysis, and several hub genes were chosen for the prediction.
Twenty hub genes were associated with the staging of esophageal cancer and gastric cancer, such as SECISBP2 L, NCOA1, HERC1, HIPK3, and MBNL1. SECISBP2 L, which encodes the SECIS binding protein 2 (SBP2), is a selenoprotein that exists in yeast, fungi and higher plants (). However, the molecular mechanism of this protein has not been elucidated. NCOA1 is nuclear receptor coactivator 1 (also known as SRC-1) belonging to the SRC gene family (; Xu et al., 2009). The abnormal expression level of SRCs has triggered cancers, such as breast, endometrial and ovarian cancers. The SRC-3 and SRC-1 genes have high expression levels in breast cancer. SRC-1 regulates polyoma enhancer activator 3 (PEA3) and promotes breast-to-lung metastasis via transcription factors (; Walsh et al., 2012; ). HERC1 is a HECT and RLD domain containing E3 ubiquitin protein ligase family member 1 in humans. E3 ubiquitin ligases play a critical role in catalyzing ubiquitin transfer from E2 enzymes to the substrate in the ubiquitylation system. Similar to most E3 ubiquitin ligases, HERC1 has been shown to be a potential target for cancer therapy (; ). Homeodomain-interacting protein kinase 3 (HIPK3) is encoded by the HIPK3 gene in humans. In NSCLC tissues, low expression of HIPK3 was associated with poor survival rates. Therefore, HIPK3 is considered a valuable biomarker for the survival of NSCLC patients (). Muscleblind-like splicing regulator 1 (MBNL1), encoded by the MBNL1 gene, is an RNA splicing protein and essential for MLL-rearranged leukemia cell growth (). The isoforms of the MBNL1 protein affect cancer development and are targets for drug development (). The results of the network indicated that these genes are also related to esophageal cancer and gastric cancer, and further research is necessary to determine the molecular mechanisms of these genes.
We eventually obtained 20 miRNAs associated with the staging of esophageal cancer and gastric cancer (), including hsa-let-7b-3p, hsa-miR-378a-5p, hsa-miR-26a-5p, hsa-miR-382-5p, and hsa-miR-29b-2-5p. Previous studies indicate that gastric cancer stem-like cells are derived from MKN-45- by the expression change of hsa-let-7b-3p (). Moreover, microRNA-378a-5p increased the expression level in melanoma cells and is a novel positive regulator of melanoma progression (Tang et al., 2018; Tupone et al., 2020). The expression of hsa-miR-26a-5p was downregulated in breast cancer tissues and appeared to have a poor prognosis for patients (). However, overexpression of hsa-miR-382-5p increased oral squamous cell carcinoma cell invasion and migration (Sun et al., 2019). These studies suggested that miRNAs are also biomarkers for various cancers; however, the prediction of the interaction network indicated that these miRNAs are also related to colon cancer, rectosigmoid junction cancer and rectal cancer ().
Statements
Data availability statement
Publicly available datasets were analyzed in this study. This data can be found here: https://www.cancer.gov/about-nci/organization/ccg/research/structural-genomics/tcga.
Author contributions
LX and LZ designed the research. ZC and ZS performed the research. ZC, ZS, ZZ, and DZ analyzed the data. ZC wrote the manuscript. All authors read and approved the manuscript.
Funding
This work was supported by the National Natural Science Foundation of China (Nos. 62002244, 62002242, and 62001311).
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.
Footnotes
1.^https://github.com/NCI-GDC/gdc-client/
2.^http://www.bioconductor.org/packages/release/bioc/html/TCGAbiolinks.html
3.^http://cancergenome.nih.gov/publications/publicationguidelines
4.^http://www.bioconductor.org/packages/release/bioc/html/edgeR.html
5.^http://www.bioconductor.org/packages/release/bioc/html/limma.html
6.^https://uswest.ensembl.org/index.html
7.^https://cloud.r-project.org/package=ggplot2
8.^https://cran.r-project.org/web/packages/VennDiagram/
References
1
AbbasG.KrasnaM. (2017). Overview of esophageal cancer.Annals of Cardiothoracic Surgery6131–136. 10.21037/acs.2017.03.03
2
AmerA.NagahA.TianT.ZhangX. (2020). Mutation Mechanisms of Breast Cancer among the Female Population in China.Current Bioinformatics15253–259. 10.2174/1574893615666191220141548
3
AminM. B.GreeneF. L.EdgeS. B.ComptonC. C.GershenwaldJ. E.BrooklandR. K.et al (2017). The Eighth Edition AJCC Cancer Staging Manual: Continuing to build a bridge from a population-based to a more “personalized” approach to cancer staging.CA: A Cancer Journal for Clinicians6793–99. 10.3322/caac.21388
4
BaiY.LongJ.LiuZ.LinJ.HuangH.WangD.et al (2019). Comprehensive analysis of a ceRNA network reveals potential prognostic cytoplasmic lncRNAs involved in HCC progression.Journal of Cellular Physiology23418837–18848. 10.1002/jcp.28522
5
BrayF.FerlayJ.SoerjomataramI.SiegelR. L.TorreL. A.JemalA. (2018). Global cancer statistics 2018: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries.CA Cancer J Clin68394–424. 10.3322/caac.21492
6
CallahanK.PonceC.CrossC.SyF.PinheiroP. (2019). Low colorectal cancer survival in the Mountain West state of Nevada: A population-based analysis.PLoS One14:e0221337. 10.1371/journal.pone.0221337
7
CaoZ.PanX.YangY.HuangY.ShenH. B. (2018). The lncLocator: a subcellular localization predictor for long non-coding RNAs based on a stacked ensemble classifier.Bioinformatics342185–2194. 10.1093/bioinformatics/bty085
8
ChenY.PalB.VisvaderJ. E.SmythG. K. (2017). Differential methylation analysis of reduced representation bisulfite sequencing experiments using edgeR.F1000Res62055. 10.12688/f1000research.13196.2
9
ChengL. (2019). Computational and Biological Methods for Gene Therapy.Current Gene Therapy19210–210.
10
ChengL.HuY.SunJ.ZhouM.JiangQ. (2018). DincRNA: a comprehensive web-based bioinformatics toolkit for exploring disease associations and ncRNA function.Bioinformatics341953–1956. 10.1093/bioinformatics/bty002
11
ChengL.WangP.TianR.WangS.GuoQ.LuoM.et al (2019a). LncRNA2Target v2.0: a comprehensive database for target genes of lncRNAs in human and mouse.Nucleic Acids Res47D140–D144. 10.1093/nar/gky1051
12
ChengL.ZhaoH.WangP.ZhouW.LuoM.LiT.et al (2019b). Computational Methods for Identifying Similar Diseases.Molecular therapy. Nucleic acids18590–604.
13
ColapricoA.SilvaT. C.OlsenC.GarofanoL.CavaC.GaroliniD.et al (2016). TCGAbiolinks: an R/Bioconductor package for integrative analysis of TCGA data.Nucleic Acids Res44e71. 10.1093/nar/gkv1507
14
DojaM. N.KaurI.AhmadT. (2020). Current State of the Art for Survival Prediction in Cancer Using Data Mining Techniques.Current Bioinformatics15174–186. 10.2174/1574893614666190902152142
15
DongY.-M.BiJ.-H.HeQ.-E.SongK. (2020). ESDA: An Improved Approach to Accurately Identify Human snoRNAs for Precision Cancer Therapy.Current Bioinformatics1534–40. 10.2174/1574893614666190424162230
16
DonovanJ.CopelandP. R. (2009). Evolutionary history of selenocysteine incorporation from the perspective of SECIS binding proteins.BMC Evolutionary Biology9:229. 10.1186/1471-2148-9-229
17
FerlayJ.ShinH.-R.BrayF.FormanD.ParkinD. (2010). Estimates of worldwide burden of cancer in 2008: GLOBOCAN.International journal of cancer. Journal international du cancer1272893–2917. 10.1002/ijc.25516
18
FerlayJ.SoerjomataramI.DikshitR.EserS.MathersC.RebeloM.et al (2015). Cancer incidence and mortality worldwide: sources, methods and major patterns in GLOBOCAN 2012.Int J Cancer136E359–E386. 10.1002/ijc.29210
19
FitzmauriceC.AbateD.AbbasiN.AbbastabarH.Abd-AllahF.Abdel-RahmanO.et al (2019). Global, Regional, and National Cancer Incidence, Mortality, Years of Life Lost, Years Lived With Disability, and Disability-Adjusted Life-Years for 29 Cancer Groups, 1990 to 2017: A Systematic Analysis for the Global Burden of Disease Study.JAMA Oncol51749–1768. 10.1001/jamaoncol.2019.2996
20
FitzmauriceC.AkinyemijuT. F.Al LamiF. H.AlamT.Alizadeh-NavaeiR.AllenC.et al (2018). Global, Regional, and National Cancer Incidence, Mortality, Years of Life Lost, Years Lived With Disability, and Disability-Adjusted Life-Years for 29 Cancer Groups, 1990 to 2016: A Systematic Analysis for the Global Burden of Disease Study.JAMA Oncol41553–1568. 10.1001/jamaoncol.2018.2706
21
FitzmauriceC.AllenC.BarberR. M.BarregardL.BhuttaZ. A.BrennerH.et al (2017). Global, Regional, and National Cancer Incidence, Mortality, Years of Life Lost, Years Lived With Disability, and Disability-Adjusted Life-years for 32 Cancer Groups, 1990 to 2015: A Systematic Analysis for the Global Burden of Disease Study.JAMA Oncol3524–548. 10.1001/jamaoncol.2016.5688
22
GhoshA.YanH. (2020). Stability Analysis at Key Positions of EGFR Related to Non-small Cell Lung Cancer.Current Bioinformatics15260–267. 10.2174/1574893614666191212112026
23
HongZ.ZengX.WeiL.LiuX. J. B. (2020). Identifying enhancer–promoter interactions with neural network based on pre-trained DNA vectors and attention mechanism.Bioinformatics361037–1043.
24
HuangZ.-M.GeH.-F.YangC.-C.CaiY.ChenZ.TianW.-Z.et al (2019). MicroRNA-26a-5p inhibits breast cancer cell growth by suppressing RNF6 expression.The Kaohsiung Journal of Medical Sciences35467–473. 10.1002/kjm2.12085
25
IliopoulosA. C.BeisG.ApostolouP.PapasotiriouI. (2020). Complex Networks, Gene Expression and Cancer Complexity: A Brief Review of Methodology and Applications.Current Bioinformatics15629–655. 10.2174/1574893614666191017093504
26
ItskovichS. S.GurunathanA.ClarkJ.BurwinkelM.WunderlichM.BergerM. R.et al (2020). MBNL1 regulates essential alternative RNA splicing patterns in MLL-rearranged leukemia.Nature Communications112369. 10.1038/s41467-020-15733-8
27
JiJ.TangJ.XiaK.-J.JiangR. (2019). LncRNA in Tumorigenesis Microenvironment.Current Bioinformatics14640–641. 10.2174/157489361407190917161654
28
JiangL.DingY.TangJ.GuoF. (2018). MDA-SKF: Similarity Kernel Fusion for Accurately Discovering miRNA-Disease Association.Frontiers in Genetics9:1–13. 10.3389/fgene.2018.00618
29
JiangL.XiaoY.DingY.TangJ.GuoF. (2019). FKL-Spa-LapRLS: an accurate method for identifying human microRNA-disease association.BMC Genomics19:11–25.
30
JiangQ.HaoY.WangG.JuanL.ZhangT.TengM.et al (2010). Prioritization of disease microRNAs through a human phenome-microRNAome network.BMC Syst Biol4(Suppl. 1):S2. 10.1186/1752-0509-4-S1-S2
31
JiangQ.WangG.JinS.LiY.WangY. (2013). Predicting human microRNA-disease associations based on support vector machine.Int J Data Min Bioinform8282–293.
32
JiangQ.WangJ.WangY.MaR.WuX.LiY. (2014). TF2LncRNA: identifying common transcription factors for a list of lncRNA genes from ChIP-Seq data.Biomed Res Int2014317642. 10.1155/2014/317642
33
JiangQ.WangY.HaoY.JuanL.TengM.ZhangX.et al (2009). miR2Disease: a manually curated database for microRNA deregulation in human disease.Nucleic Acids Res37D98–D104. 10.1093/nar/gkn714
34
JinS.ZengX.FangJ.LinJ.ChanS. Y.ErzurumS. C.et al (2019). A network-based approach to uncover microRNA-mediated disease comorbidities and potential pathobiological implications.NPJ systems biology and applications51–11.
35
JinS.ZengX.XiaF.HuangW.LiuX. (2020). Application of deep learning methods in biological networks.Briefings in Bioinformatics2bbaa043. 10.1093/bib/bbaa043
36
LangfelderP.HorvathS. (2008). WGCNA: an R package for weighted correlation network analysis.BMC Bioinformatics9:559. 10.1186/1471-2105-9-559
37
LangfelderP.MischelP. S.HorvathS. (2013). When Is Hub Gene Selection Better than Standard Meta-Analysis?PLoS One8:e61505. 10.1371/journal.pone.0061505
38
LawC.AlhamdooshM.SuS.SmythG.RitchieM. (2016). RNA-seq analysis is easy as 1-2-3 with limma, Glimma and edgeR.F1000Research51408. 10.12688/f1000research.9005.1
39
LiJ.ChangM.GaoQ.SongX.GaoZ. (2020). Lung Cancer Classification and Gene Selection by Combining Affinity Propagation Clustering and Sparse Group Lasso.Current Bioinformatics15703–712. 10.2174/1574893614666191017103557
40
LiW.LiuJ.ZhaoH. (2020). Identification of a nomogram based on long non-coding RNA to improve prognosis prediction of esophageal squamous cell carcinoma.Aging121512–1526. 10.18632/aging.102697
41
LiZ.ZhangT.LeiH.WeiL.LiuY.ShiY.et al (2020). Research on Gastric Cancer’s Drug-resistant Gene Regulatory Network Model.Current Bioinformatics15225–234. 10.2174/1574893614666190722102557
42
LiS.JiangL.TangJ.GaoN.GuoF. (2020). Kernel Fusion Method for Detecting Cancer Subtypes via Selecting Relevant Expression Data.Front. Genet.11:979. 10.3389/fgene.2020.00979
43
LiangC.ChangluQ.HeZ.TongzeF.XueZ. (2019). gutMDisorder: a comprehensive database for dysbiosis of the gut microbiota in disorders and interventions.Nucleic Acids Research487603.
44
LinM.LiX.GuoH.JiF.YeL.MaX.et al (2019). Identification of Bone Metastasis-associated Genes of Gastric Cancer by Genome-wide Transcriptional Profiling.Current Bioinformatics1462–69. 10.2174/1574893612666171121154017
45
LiuX.HongZ.LiuJ.LinY.Rodríguez-PatónA.ZouQ.et al (2020). Computational methods for identifying the critical nodes in biological networks.Briefings in Bioinformatics21486–497.
46
LiuY.QianL.YangJ.HuangH.FengJ.LiX.et al (2018). The expression level and prognostic value of HIPK3 among non-small-cell lung cancer patients in China.OncoTargets and therapy117459–7469. 10.2147/OTT.S166878
47
LiuY.ZengX.HeZ.ZouQ. (2017). Inferring MicroRNA-Disease Associations by Random Walk on a Heterogeneous Network with Multiple Data Sources.IEEE/ACM Transactions on Computational Biology and Bioinformatics14905–915. 10.1109/tcbb.2016.2550432
48
LongJ.BaiY.YangX.LinJ.YangX.WangD.et al (2019). Construction and comprehensive analysis of a ceRNA network to reveal potential prognostic biomarkers for hepatocellular carcinoma.Cancer Cell International1990. 10.1186/s12935-019-0817-y
49
MaX.XiB.ZhangY.ZhuL.SuiX.TianG.et al (2020). A Machine Learning-based Diagnosis of Thyroid Cancer Using Thyroid Nodules Ultrasound Images.Current Bioinformatics15349–358. 10.2174/1574893614666191017091959
50
MazaE. (2016). In Papyro Comparison of TMM (edgeR), RLE (DESeq2), and MRN Normalization Methods for a Simple Two-Conditions-Without-Replicates RNA-Seq Experimental Design.Front Genet7:164. 10.3389/fgene.2016.00164
51
NadiaRamanaJ. (2020). The Human OncoBiome Database: A Database of Cancer Microbiome Datasets.Current Bioinformatics15472–477. 10.2174/1574893614666190902152727
52
OñateS. A.TsaiS. Y.TsaiM. J.O’MalleyB. W. (1995). Sequence and characterization of a coactivator for the steroid hormone receptor superfamily.Science2701354–1357. 10.1126/science.270.5240.1354
53
OzkanA.IsgorS. B.SengulG.IsgorY. G. (2019). Benchmarking Classification Models for Cell Viability on Novel Cancer Image Datasets.Current Bioinformatics14108–114. 10.2174/1574893614666181120093740
54
PanH.GuoC.PanJ.GuoD.SongS.ZhouY.et al (2019). Construction of a Competitive Endogenous RNA Network and Identification of Potential Regulatory Axis in Gastric Cancer.Frontiers in Oncology9:912. 10.3389/fonc.2019.00912
55
PawaN.ArulampalamT.NortonJ. D. (2011). Screening for colorectal cancer: established and emerging modalities.Nat Rev Gastroenterol Hepatol8711–722. 10.1038/nrgastro.2011.205
56
QianW.FengY.LiJ.PengW.GuQ.ZhangZ.et al (2019). Construction of ceRNA networks reveals differences between distal and proximal colon cancers.Oncol Rep413027–3040. 10.3892/or.2019.7083
57
QinL.XuY.XuY.MaG.LiaoL.WuY.et al (2015). NCOA1 promotes angiogenesis in breast tumors by simultaneously enhancing both HIF1α- and AP-1-mediated VEGFa transcription.Oncotarget623890–23904. 10.18632/oncotarget.4341
58
RahmanR.SchmaltzC.JacksonC. S.SimoesE. J.Jackson-ThompsonJ.IbdahJ. A. (2015). Increased risk for colorectal cancer under age 50 in racial and ethnic minorities living in the United States.Cancer Med41863–1870. 10.1002/cam4.560
59
RedmondA. M.BaneF. T.StaffordA. T.McIlroyM.DillonM. F.CrottyT. B.et al (2009). Coassociation of estrogen receptor and p160 proteins predicts resistance to endocrine treatment; SRC-1 is an independent predictor of breast cancer recurrence.Clin Cancer Res152098–2106. 10.1158/1078-0432.Ccr-08-1649
60
RustgiA.El-SeragH. (2014). Esophageal Carcinoma.The New England journal of medicine3712499–2509. 10.1056/NEJMra1314530
61
RustgiA.El-SeragH. B. (2015). Esophageal carcinoma.N Engl J Med3721472–1473. 10.1056/NEJMc1500692
62
SalehiZ.AkramiH.SisakhtnezhadS. (2018). The effect of placenta growth factor knockdown on hsa-miR-22-3p, hsa-let-7b-3p, hsa-miR-451b, and hsa-mir-4290 expressions in MKN-45-derived gastric cancer stem-like cells.Middle East Journal of Cancer9113–122.
63
SardinaD. S.AlaimoS.FerroA.PulvirentiA.GiugnoR. (2017). A novel computational method for inferring competing endogenous interactions.Brief Bioinform181071–1081. 10.1093/bib/bbw084
64
SchneiderT.Martinez-MartinezA.Cubillos-RojasM.BartronsR.VenturaF.RosaJ. L. (2018). The E3 ubiquitin ligase HERC1 controls the ERK signaling pathway targeting C-RAF for degradation.Oncotarget931531–31548. 10.18632/oncotarget.25847
65
ShenC.JiangL.DingY.TangJ.GuoF. (2019). LPI-KTASLP: Prediction of lncRNA-Protein Interaction by Semi-Supervised Link Learning with Multivariate Information.IEEE Access713486–13496.
66
SlatteryM. L.WolffR. K.HerrickJ. S.CaanB. J.PotterJ. D. (2007). IL6 genotypes and colon and rectal cancer.Cancer Causes & Control181095–1105. 10.1007/s10552-007-9049-x
67
SongJ.LuY.SunW.HanM.ZhangY.ZhangJ. (2019). Changing expression profiles of lncRNAs, circRNAs and mRNAs in esophageal squamous carcinoma.Oncology Letters185363–5373. 10.3892/ol.2019.10880
68
SuH.LinF.DengX.ShenL.FangY.FeiZ.et al (2016). Profiling and bioinformatics analyses reveal differential circular RNA expression in radioresistant esophageal cancer cells.Journal of translational medicine14225–225. 10.1186/s12967-016-0977-7
69
SunL.-P.XuK.CuiJ.YuanD.-Y.ZouB.LiJ.et al (2019). Cancer-associated fibroblast-derived exosomal miR-382-5p promotes the migration and invasion of oral squamous cell carcinoma.Oncology reports421319–1328. 10.3892/or.2019.7255
70
TangW.WanS.YangZ.TeschendorffA. E.ZouQ. (2018). Tumor origin detection with tissue-specific miRNA and DNA methylation markers.Bioinformatics34398–406. 10.1093/bioinformatics/btx622
71
TayY.RinnJ.PandolfiP. P. (2014). The multilayered complexity of ceRNA crosstalk and competition.Nature505344–352. 10.1038/nature12986
72
ThomsonD. W.DingerM. E. (2016). Endogenous microRNA sponges: evidence and controversy.Nature Reviews Genetics17272–283. 10.1038/nrg.2016.20
73
TorreL. A.BrayF.SiegelR. L.FerlayJ.Lortet-TieulentJ.JemalA. (2015). Global cancer statistics, 2012.CA Cancer J Clin6587–108. 10.3322/caac.21262
74
TuponeM. G.D’AguannoS.Di MartileM.ValentiniE.DesideriM.TrisciuoglioD.et al (2020). microRNA-378a-5p iS a novel positive regulator of melanoma progression.Oncogenesis922. 10.1038/s41389-020-0203-6
75
WalshC. A.QinL.TienJ. C.-Y.YoungL. S.XuJ. (2012). The function of steroid receptor coactivator-1 in normal tissues and cancer.International journal of biological sciences8470–485. 10.7150/ijbs.4125
76
WangJ.LiuX.WuH.NiP.GuZ.QiaoY.et al (2010). CREB up-regulates long non-coding RNA, HULC expression through interaction with microRNA-372 in liver cancer.Nucleic Acids Res385366–5383. 10.1093/nar/gkq285
77
WeiL.LiaoM.GaoY.JiR.HeZ.ZouQ. (2014). Improved and Promising Identification of Human MicroRNAs by Incorporating a High-Quality Negative Set.IEEE/ACM Transactions on Computational Biology and Bioinformatics11192–201. 10.1109/tcbb.2013.146
78
WeiL.SuR.WangB.LiX.ZouQ.GaoX. (2019a). Integration of deep feature representations and handcrafted features to improve the prediction of N-6-methyladenosine sites.Neurocomputing3243–9. 10.1016/j.neucom.2018.04.082
79
WeiL.XingP.ShiG.JiZ.ZouQ. (2019b). Fast Prediction of Protein Methylation Sites Using a Sequence-Based Feature Selection Technique.Ieee-Acm Transactions on Computational Biology and Bioinformatics161264–1273. 10.1109/tcbb.2017.2670558
80
WeiL.ZhouC.ChenH.SongJ.SuR. (2018). ACPred-FL: a sequence-based predictor using effective feature representation to improve the prediction of anti-cancer peptides.Bioinformatics344007–4016. 10.1093/bioinformatics/bty451
81
WuW.YangZ.LongF.LuoL.DengQ.WuJ.et al (2019). COL1A1 and MZB1 as the hub genes influenced the proliferation, invasion, migration and apoptosis of rectum adenocarcinoma cells by weighted correlation network analysis.Bioorganic Chemistry95103457. 10.1016/j.bioorg.2019.103457
82
XuJ.WuR. C.O’MalleyB. W. (2009). Normal and cancer-related functions of the p160 steroid receptor co-activator (SRC) family.Nat Rev Cancer9615–630. 10.1038/nrc2695
83
YangL.GaoH.WuK.ZhangH.LiC.TangL. (2020). Identification of Cancerlectins By Using Cascade Linear Discriminant Analysis and Optimal g-gap Tripeptide Composition.Current Bioinformatics15528–537. 10.2174/1574893614666190730103156
84
ZengW.WangF.MaY.LiangX.ChenP. (2019). Dysfunctional Mechanism of Liver Cancer Mediated by Transcription Factor and Non-coding RNA.Current Bioinformatics14100–107. 10.2174/1574893614666181119121916
85
ZengX.ZhuS.LiuX.ZhouY.NussinovR.ChengF. J. B. (2019). deepDR: a network-based deep learning approach to in silico drug repositioning.Bioinformatics355191–5198.
86
ZengX.LinY.HeY.LvL.MinX. (2020a). Deep collaborative filtering for prediction of disease genes.IEEE/ACM Transactions on Computational Biology and Bioinformatics171639–1647.
87
ZengX.ZhuS.LuW.LiuZ.HuangJ.ZhouY.et al (2020b). Target identification among known drugs by deep learning from heterogeneous networks.Chemical Science111775–1797.
88
ZengX.ZhangX.ZouQ. (2016). Integrative approaches for predicting microRNA function and prioritizing disease-related microRNA using biological interaction networks.Briefings in Bioinformatics17193–203. 10.1093/bib/bbv033
89
ZengX. X.LiuL.LuL. Y.ZouQ. (2018). Prediction of potential disease-associated microRNAs using structural perturbation method.Bioinformatics342425–2432. 10.1093/bioinformatics/bty112
90
ZerbinoD. R.AchuthanP.AkanniW.AmodeM. R.BarrellD.BhaiJ.et al (2018). Ensembl 2018.Nucleic Acids Research46D754–D761. 10.1093/nar/gkx1098
91
ZhangL.HeY.SongH.WangX.LuN.SunL.et al (2020). Elastic Net Regularized Softmax Regression Methods for Multi-subtype Classification in Cancer.Current Bioinformatics15212–224. 10.2174/1574893613666181112141724
92
ZhangZ. M.TanJ. X.WangF.DaoF. Y.ZhangZ. Y.LinH. (2020). Early Diagnosis of Hepatocellular Carcinoma Using Machine Learning Method.Front Bioeng Biotechnol8:254. 10.3389/fbioe.2020.00254
93
ZhangX.ZouQ.Rodriguez-PatonA.ZengX. (2017). Meta-path methods for prioritizing candidate disease miRNAs.IEEE/ACM Transactions on Computational Biology and Bioinformatics16283–291.
94
ZhangY.KouC.WangS.ZhangY. (2019). Genome-wide Differential-based Analysis of the Relationship between DNA Methylation and Gene Expression in Cancer.Current Bioinformatics14783–792. 10.2174/1574893614666190424160046
95
ZhaoT.HuY.PengJ.ChengL. (2020). DeepLGP: a novel deep learning method for prioritizing lncRNA target genes.Bioinformatics364466–4472. 10.1093/bioinformatics/btaa428
96
ZouQ.LiJ.SongL.ZengX.WangG. (2016). Similarity computation strategies in the microRNA-disease network: a survey.Briefings in Functional Genomics1555–64. 10.1093/bfgp/elv024
Summary
Keywords
cancer, long non-coding RNA, WGCNA, MicroRNA, hub genes
Citation
Chen Z, Shen Z, Zhang Z, Zhao D, Xu L and Zhang L (2021) RNA-Associated Co-expression Network Identifies Novel Biomarkers for Digestive System Cancer. Front. Genet. 12:659788. doi: 10.3389/fgene.2021.659788
Received
28 January 2021
Accepted
25 February 2021
Published
26 March 2021
Volume
12 - 2021
Edited by
Wei Jiang, Nanjing University of Aeronautics and Astronautics, China
Reviewed by
Yi Xiong, Shanghai Jiao Tong University, China; Lei Deng, Central South University, China
Updates
Copyright
© 2021 Chen, Shen, Zhang, Zhao, Xu and Zhang.
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: Lei Xu, csleixu@szpt.edu.cnLijun Zhang, c7zlj@szpt.edu.cn
†These authors have contributed equally to this work
This article was submitted to RNA, 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.