Co-expression Network Analysis of Biomarkers for Adrenocortical Carcinoma

Molecular and Genetics Perspectives on Primary Adrenocortical Hyperfunction Disorders

Adrenocortical disorders encompass a broad spectrum of conditions ranging from benign hyperplasia to malignant tumors, significantly disrupting hormone balance and causing a variety of clinical manifestations. By leveraging next-generation sequencing and in silico analyses, recent studies have uncovered the genetic and molecular pathways implicated in these transitions. In this review, we explored the molecular and genetic alterations in adrenocortical disorders, with a particular focus on the transitions from normal adrenal function to hyperfunction. The insights gained are intended to enhance diagnostic and therapeutic strategies, offering up-to-date knowledge for managing these complex conditions effectively.

Gene network landscape of mouse splenocytes reveals integrin complex as the A151 ODN-responsive hub molecule in the immune transcriptome

Homeostatic restoration of an inflammatory response requires quenching of the immune system after pathogen threats vanish. A continued assault orchestrated by host defense results in tissue destruction or autoimmunity. A151 is the epitome of synthetic oligodeoxynucleotides (ODNs) that curb the immune response by a subset of white corpuscles through repetitive telomere-derived TTAGGG sequences. Currently, the genuine effect of A151 on the immune cell transcriptome remains unknown. Here, we leveraged an integrative approach where weighted gene co-expression network analysis (WGCNA), differential gene expression analysis, and gene set enrichment analysis (GSEA) of our in-house microarray datasets aided our understanding of how A151 ODN suppresses the immune response in mouse splenocytes. Our bioinformatics results, together with experimental validations, indicated that A151 ODN acts on components of integrin complexes, Itgam and Itga6, to interfere with immune cell adhesion and thereby suppresses the immune response in mice. Moreover, independent lines of evidence in this work converged on the observation that cell adhesion by integrin complexes serves as a focal point for cellular response to A151 ODN treatment in immune cells. Taken together, the outcome of this study sheds light on the molecular basis of immune suppression by a clinically useful DNA-based therapeutic agent.

Comprehensive Characterization of the Regulatory Landscape of Adrenocortical Carcinoma: Novel Transcription Factors and Targets Associated with Prognosis

We reconstructed a transcriptional regulatory network for adrenocortical carcinoma (ACC) using transcriptomic and clinical data from The Cancer Genome Atlas (TCGA)-ACC cohort. We investigated the association of transcriptional regulatory units (regulons) with overall survival, molecular phenotypes, and immune signatures. We annotated the ACC regulons with cancer hallmarks and assessed single sample regulon activities in the European Network for the Study of Adrenal Tumors (ENSAT) cohort. We found 369 regulons associated with overall survival and subdivided them into four clusters: RC1 and RC2, associated with good prognosis, and RC3 and RC4, associated with worse outcomes. The RC1 and RC3 regulons were highly correlated with the 'Steroid Phenotype,' while the RC2 and RC4 regulons were highly correlated with a molecular proliferation signature. We selected two regulons, NR5A1 (steroidogenic factor 1, SF-1) and CENPA (Centromeric Protein A), that were consistently associated with overall survival for further downstream analyses. The CENPA regulon was the primary regulator of MKI-67 (a marker of proliferation KI-67), while the NR5A1 regulon is a well-described transcription factor (TF) in ACC tumorigenesis. We also found that the ZBTB4 (Zinc finger and BTB domain-containing protein 4) regulon, which is negatively associated with CENPA in our transcriptional regulatory network, is also a druggable anti-tumorigenic TF. We anticipate that the ACC regulons may be used as a reference for further investigations concerning the complex molecular interactions in ACC tumors.

Identification of Survival-Related Genes in Acute Myeloid Leukemia (AML) Based on Cytogenetically Normal AML Samples Using Weighted Gene Coexpression Network Analysis

The prognosis of acute myeloid leukemia (AML) remains a challenge. In this study, we applied the weighted gene coexpression network analysis (WGCNA) to find survival-specific genes in AML based on 42 adult CN-AML samples from The Cancer Genome Atlas (TCGA) database. Eighteen hub genes (ABCA13, ANXA3, ARG1, BTNL8, C11orf42, CEACAM1, CEACAM3, CHI3L1, CRISP2, CYP4F3, GPR84, HP, LTF, MMP8, OLR1, PADI2, RGL4, and RILPL1) were found to be related to AML patient survival time. We then compared the hub gene expression levels between AML peripheral blood (PB) samples ( $n=162$ ) and control healthy whole blood samples ( $n=337$ ). Seventeen of the hub genes showed lower expression levels in AML PB samples. The gene expression analysis was also done among AML BM (bone marrow) samples of different stages: diagnosis ( $n=142$ ), posttreatment ( $n=42$ ), and recurrent ( $n=12$ ) stages. The results showed a significant increase of ANXA3, CEACM1, RGL4, RILPL1, and HP in posttreatment samples compared to diagnosis and/or recurrent samples. Transcription factor (TF) prediction of the hub genes suggested LTF as the top hit, overlapping 10 hub genes, while LTF itself is just one of the hub genes. Also, 3671 correlation links were shown between 128 mRNAs and 209 lncRNAs found in survival time-related modules. Generally, we identified candidate mRNA biomarkers based on CN-AML data which can be extensively used in AML prognosis. In addition, we mapped their potential regulatory mechanisms with correlated lncRNAs, providing new insights into potential targets for therapies in AML.

Root and shoot biology of Arabidopsis halleri dissected by WGCNA: an insight into the organ pivotal pathways and genes of an hyperaccumulator

Arabidopsis halleri is a hyperaccumulating pseudo-metallophyte and an emerging model to explore molecular basis of metal tolerance and hyperaccumulation. In this regard, understanding of interacting genes can be a crucial aspect as these interactions regulate several biological functions at molecular level in response to multiple signals. In this current study, we applied a weighted gene co-expression network analysis (WGCNA) on root and shoot RNA-seq data of A. halleri to predict the related scale-free organ specific co-expression networks, for the first time. A total of 19,653 genes of root and 18,081 genes of shoot were grouped into 14 modules and subjected to GO and KEGG enrichment analysis. "Photosynthesis" and "photosynthesis-antenna proteins" were identified as the most enriched and common pathway to both root and shoot. Whereas "glucosinolate biosynthesis," "autophagy," and "SNARE interactions in vesicular transport" were specific to root, and "circadian rhythm" was found to be enriched only in shoot. Later, hub and bottleneck genes were identified in each module by using cytoHubba plugin based on Cytoscape and scoring the relevance of each gene to the topology of the network. The modules with the most significant differential expression pattern across control and treatment (Cd-Zn treatment) were selected and their hub and bottleneck genes were screened to validate their possible involvement in heavy metal stress. Moreover, we combined the analysis of co-expression modules together with protein-protein interactions (PPIs), confirming some genes as potential candidates in plant heavy metal stress and as biomarkers. The results from this analysis shed the light on the pivotal functions to the hyperaccumulative trait of A. halleri, giving perspective to new paths for future research on this species.

Identification of Four Novel Prognostic Biomarkers and Construction of Two Nomograms in Adrenocortical Carcinoma: A Multi-Omics Data Study via Bioinformatics and Machine Learning Methods

Background: Adrenocortical carcinoma (ACC) is an orphan tumor which has poor prognoses. Therefore, it is of urgent need for us to find candidate prognostic biomarkers and provide clinicians with an accurate method for survival prediction of ACC via bioinformatics and machine learning methods.

Methods: Eight different methods including differentially expressed gene (DEG) analysis, weighted correlation network analysis (WGCNA), protein-protein interaction (PPI) network construction, survival analysis, expression level comparison, receiver operating characteristic (ROC) analysis, and decision curve analysis (DCA) were used to identify potential prognostic biomarkers for ACC via seven independent datasets. Linear discriminant analysis (LDA), K-nearest neighbor (KNN), support vector machine (SVM), and time-dependent ROC were performed to further identify meaningful prognostic biomarkers (MPBs). Cox regression analyses were performed to screen factors for nomogram construction.

Results: We identified nine hub genes correlated to prognosis of patients with ACC. Furthermore, four MPBs (ASPM, BIRC5, CCNB2, and CDK1) with high accuracy of survival prediction were screened out, which were enriched in the cell cycle. We also found that mutations and copy number variants of these MPBs were associated with overall survival (OS) of ACC patients. Moreover, MPB expressions were associated with immune infiltration level. Two nomograms [OS-nomogram and disease-free survival (DFS)-nomogram] were established, which could provide clinicians with an accurate, quick, and visualized method for survival prediction.

Conclusion: Four novel MPBs were identified and two nomograms were constructed, which might constitute a breakthrough in treatment and prognosis prediction of patients with ACC.

eEF1A1 promotes colorectal cancer progression and predicts poor prognosis of patients

Colorectal cancer (CRC) is a major leading cause of cancer mortality worldwide in which dysregulated protein synthesis plays an etiologic role. The eukaryotic elongation factor 1 A1 (eEF1A1) exerts significant effects on protein synthesis by contributing to peptide chain extension. Whereas its role in CRC remains to be investigated. In this study, we found that the mRNA and protein levels of eEF1A1 were significantly upregulated in CRC cell lines and tissues. Elevated expression of eEF1A1 was correlated with shorter overall survival in 94 CRC patients. The inhibition of proliferation and cell cycle block were observed in CRC cells after eEF1A1 downregulation. Mechanistically, weighted gene correlation network analysis and further Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis suggested that mitogen-activated protein kinases (MAPKs) signaling pathways were significantly enriched in high-eEF1A1 expression group, and the levels of phosphorylated p38/JNK/ERK MAPK were dramatically decreased after eEF1A1 downregulation. Overexpression of eEF1A1 in CRC correlated with a poor prognosis. Collectively, this study determined the oncogenic role of eEF1A1 in CRC proliferation and tumorigenesis. eEF1A1 might be a promising therapeutic target and prognostic biomarker in CRC.

Pathological stage-associated non-coding RNA long intergenic non-protein coding RNA 1234 (LINC01234) participation in cell cycle regulation in adrenocortical carcinoma through bromodomain-containing protein 4 (BRD4) expression mediation via sponging microRNA (miR)-140-3p

Many researches indicated that long non-coding RNAs (lncRNAs) were involved in the malignant progression of tumors, including Adrenocortical Carcinoma (ACC). However, as for most lncRNAs, their biological behaviors and molecular mechanism remain unclear in ACC. In the present research, weighted gene co-expression network analysis (WGCNA) was used to identify pathologically relevant gene, including lncRNAs. By comparing their expressions in GSE61359 tumors and normal controls, long intergenic non-protein coding RNA 1234 (LINC01234) was selected to investigate the clinical significance, biological function, and mechanism in ACC. Data mining revealed that LINC01234 expression was significantly up-regulated in ACC patients, and a shorter survival time presents in patients with higher LINC01234 expression compared to that in patients with lower LINC01234 expression. Further, LINC01234 silencing resulted in cells growth arrest in vitro and in vivo. Mechanism studies suggested that LINC01234 silencing induced cell cycle arrest, and bromodomain-containing protein 4 (BRD4) overexpression could restore this phenomenon. Further research showed that LINC01234 could mediate BRD4 expression through competitively sequestering microRNA (miR)-140-3p, as evidenced by the positive correlation of LINC01234 with BRD4 and inverse correlation with miR-140-3p expression. Luciferase activity assay also verified the targeting relationship between LINC01234, BRD4 and miR-140-3p. And up-regulated LINC01234 in ACC cells significantly reversed the degradation of BRD4 by miR-140-3p. Collectively, we deduce that LINC01234 functions as a ceRNA to regulate BRD4 expression by sponging miR-140-3p in ACC progress. Our findings have the potential to provide a new target for the diagnosis and treatment of ACC.

Spindle and Kinetochore-Associated Complex is Associated With Poor Prognosis in Adrenocortical Carcinoma

INTRODUCTION

The spindle and kinetochore-associated (SKA) complex, composed of three subunits (SKA1, SKA2, and SKA3), stabilizes spindle microtubule attachment to the kinetochore (KT) in the middle stage of mitosis. High expression of this complex is associated with poor prognosis for several tumors. However, the potential role of SKA complex overexpression in rare malignant diseases, such as adrenocortical carcinoma (ACC), has not been well investigated.

MATERIALS AND METHODS

In this study, we used several databases to explore the relationship between SKA subunit expression and prognosis in ACC patients. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genome (KEGG) databases were used to analyze enriched pathways in ACC.

RESULTS

The results suggest that each of the three SKA subunits are overexpressed in ACC and that high expression is correlated with poor patient prognosis. Overexpression of the SKA complex is associated with the expression of organelle fission, nuclear division, and chromosome segregation pathways. Furthermore, differential expression of hub genes for proteins that interact physically or functionally with the SKA complex (CCNB2, UBE2C, BUB1B, TPX2, CCNA2, CDCA8, CCNB1, MELK, TOP2A, and KIF2C) revealed additional potential biomarkers for ACC.

CONCLUSIONS

Our findings provide additional understanding of the mechanisms of ACC and suggest an approach for biomarker discovery using publicly available resources.

Overexpression of SKA Complex Is Associated With Poor Prognosis in Gliomas

The spindle and kinetochore-associated complex is composed of three members: SKA1, SKA2, and SKA3. It is necessary for stabilizing spindle microtubules attaching to kinetochore (KT) in the middle stage of mitosis. The SKA complex is associated with poor prognosis in several human cancers. However, the role of SKA complex in rare malignant diseases, such as gliomas, has not been fully investigated. We investigated several databases, including Oncomine, UALCAN, and cBioPortal to explore the expression profile and prognostic significance of SKA complex in patients with gliomas. Gene ontology and Kyoto Encyclopedia of Genes and Genome pathways were used to analyze the potential enriched pathways. The genes co-expressed with SKA complex were identified and used for developing a protein-protein interaction (PPI) network using the STRING database. We found a significant overexpression of the mRNA levels of SKA1, SKA2, and SKA3 in patients with glioma patients. Higher expression of SKA1 and SKA3, but not SKA2, was significantly correlated with shorter overall survival of patients with glioma. In glioma, SKA complex was found to be involved in nuclear division, chromosome segregation, and DNA replication. The results of PPI network identified 10 hub genes (CCNB2, UBE2C, BUB1B, TPX2, CCNA2, CCNB1, MELK, TOP2A, PBK, and KIF11), all of which were overexpressed and negatively associated with prognosis of patients with glioma. In conclusion, our study sheds new insights into the biological role and prognostic significance of SKA complex in glioma.

Endocrine and Neuroendocrine Tumors Special Issue—Checkpoint Inhibitors for Adrenocortical Carcinoma and Metastatic Pheochromocytoma and Paraganglioma: Do They Work?

Adrenocortical cancers and metastatic pheochromocytomas are the most common malignancies originating in the adrenal glands. Metastatic paragangliomas are extra-adrenal tumors that share similar genetic and molecular profiles with metastatic pheochromocytomas and, subsequently, these tumors are studied together. Adrenocortical cancers and metastatic pheochromocytomas and paragangliomas are orphan diseases with limited therapeutic options worldwide. As in any other cancers, adrenocortical cancers and metastatic pheochromocytomas and paragangliomas avoid the immune system. Hypoxia-pseudohypoxia, activation of the PD-1/PD-L1 pathway, and/or microsatellite instability suggest that immunotherapy with checkpoint inhibitors could be a therapeutic option for patients with these tumors. The results of clinical trials with checkpoint inhibitors for adrenocortical carcinoma or metastatic pheochromocytoma or paraganglioma demonstrate limited benefits; nevertheless, these results also suggest interesting mechanisms that might enhance clinical responses to checkpoint inhibitors. These mechanisms include the normalization of tumor vasculature, modification of the hormonal environment, and vaccination with specific tumor antigens. Combinations of checkpoint inhibitors with classical therapies, such as chemotherapy, tyrosine kinase inhibitors, radiopharmaceuticals, and/or novel therapies, such as vaccines, should be evaluated in clinical trials.

Weighted Gene Coexpression Network Analysis in Mouse Livers following Ischemia-Reperfusion and Extensive Hepatectomy

In mouse models, the recovery of liver volume is mainly mediated by the proliferation of hepatocytes after partial hepatectomy that is commonly accompanied with ischemia-reperfusion. The identification of differently expressed genes in liver following partial hepatectomy benefits the better understanding of the molecular mechanisms during liver regeneration (LR) with appliable clinical significance. Briefly, studying different gene expression patterns in liver tissues collected from the mice group that survived through extensive hepatectomy will be of huge critical importance in LR than those collected from the mice group that survived through appropriate hepatectomy. In this study, we performed the weighted gene coexpression network analysis (WGCNA) to address the central candidate genes and to construct the free-scale gene coexpression networks using the identified dynamic different expressive genes in liver specimens from the mice with 85% hepatectomy (20% for seven-day survial rate) and 50% hepatectomy (100% for seven-day survial rate under ischemia-reperfusion condition compared with the sham group control mice). The WGCNA combined with Kyoto Encyclopedia of Genes and Genomes (KEGG) and Gene Ontology (GO) enrichment analyses pinpointed out the apparent distinguished importance of three gene expression modules: the blue module for apoptotic process, the turquoise module for lipid metabolism, and the green module for fatty acid metabolic process in LR following extensive hepatectomy. WGCNA analysis and protein-protein interaction (PPI) network construction highlighted FAM175B, OGT, and PDE3B were the potential three hub genes in the previously mentioned three modules. This work may help to provide new clues to the future fundamental study and treatment strategy for LR following liver injury and hepatectomy.

Identification of a prognostic chemoresistance-related gene signature associated with immune microenvironment in breast cancer

Breast cancer is the most common form of cancer among women globally, and chemoresistance is a major challenge to disease treatment that is associated with a poor prognosis. This study was formulated to identify a reliable prognostic biosignature capable of predicting the survival of patients with chemoresistant breast cancer (CRBC) and evaluating the associated tumor immune microenvironment. Through a series of protein-protein interaction and weighted correlation network analyses, genes that were significantly associated with breast cancer chemoresistance were identified. Moreover, univariate Cox regression and lasso-penalized Cox regression analyses were employed to generate a prognostic model, and the prognostic utility of this model was then assessed using time-dependent receiver operating characteristic (ROC) and Kaplan-Meier survival curves. Finally, The CIBERSORT and ESTIMATE algorithms were additionally leveraged to assess relationships between the tumor immune microenvironment and patient prognostic signatures. Overall, a multigenic prognostic biosignature capable of predicting CRBC patient risk was successfully developed based on bioinformatics analysis and in vitro experiments. This biosignature was able to stratify CRBC patients into high- and low-risk subgroups. ROC curves also revealed that this biosignature achieved high diagnostic efficiency, and multivariate regression analyses indicated that this risk signature was an independent risk factor linked to CRBC patient outcomes. In addition, this signature was associated with the infiltration of the tumor microenvironment by multiple immune cell types. In conclusion, the chemoresistance-associated prognostic gene signature developed herein was able to effectively evaluate the prognosis of CRBC patients and to reflect the overall composition of the tumor immune microenvironment.

A Hypoxia Signature for Predicting Prognosis and Tumor Immune Microenvironment in Adrenocortical Carcinoma

Adrenocortical carcinoma (ACC) is a rare malignancy with dismal prognosis. Hypoxia is one of characteristics of cancer leading to tumor progression. For ACC, however, no reliable prognostic signature on the basis of hypoxia genes has been built. Our study aimed to develop a hypoxia-associated gene signature in ACC. Data of ACC patients were obtained from TCGA and GEO databases. The genes included in hypoxia risk signature were identified using the Cox regression analysis as well as LASSO regression analysis. GSEA was applied to discover the enriched gene sets. To detect a possible connection between the gene signature and immune cells, the CIBERSORT technique was applied. In ACC, the hypoxia signature including three genes (CCNA2, COL5A1, and EFNA3) was built to predict prognosis and reflect the immune microenvironment. Patients with high-risk scores tended to have a poor prognosis. According to the multivariate regression analysis, the hypoxia signature could be served as an independent indicator in ACC patients. GSEA demonstrated that gene sets linked to cancer proliferation and cell cycle were differentially enriched in high-risk classes. Additionally, we found that PDL1 and CTLA4 expression were significantly lower in the high-risk group than in the low-risk group, and resting NK cells displayed a significant increase in the high-risk group. In summary, the hypoxia risk signature created in our study might predict prognosis and evaluate the tumor immune microenvironment for ACC.

How to Differentiate Benign from Malignant Adrenocortical Tumors?

Simple Summary

Adrenocortical carcinoma is a rare cancer with a poor prognosis. Adrenal tumors are, however, commonly identified in clinical practice. Discrimination between benign and malignant adrenal tumors is of great importance to determine the appropriate treatment and follow-up strategy. This review summarizes the current diagnostic strategies and challenges to distinguish benign from malignant adrenal lesions. We will focus both on radiological and biochemical assessments, enabling diagnosis of the adrenal lesion preoperatively, and on histopathological and a wide variety of molecular assessments that can be done after surgical removal of the adrenal lesion. Furthermore, new non-invasive strategies such as liquid biopsies, in which blood samples are used to study circulating tumor cells, tumor DNA and microRNA, will be addressed in this review.

Abstract

Adrenocortical carcinoma (ACC) is a rare cancer with a poor prognosis. Adrenal incidentalomas are, however, commonly identified in clinical practice. Discrimination between benign and malignant adrenal tumors is of great importance considering the large differences in clinical behavior requiring different strategies. Diagnosis of ACC starts with a thorough physical examination, biochemical evaluation, and imaging. Computed tomography is the first-level imaging modality in adrenal tumors, with tumor size and Hounsfield units being important features for determining malignancy. New developments include the use of urine metabolomics, also enabling discrimination of ACC from adenomas preoperatively. Postoperatively, the Weiss score is used for diagnosis of ACC, consisting of nine histopathological criteria. Due to known limitations as interobserver variability and lack of accuracy in borderline cases, much effort has been put into new tools to diagnose ACC. Novel developments vary from immunohistochemical markers and pathological scores, to markers at the level of DNA, methylome, chromosome, or microRNA. Molecular studies have provided insights into the most promising and most frequent alterations in ACC. The use of liquid biopsies for diagnosis of ACC is studied, although in a small number of patients, requiring further investigation. In this review, current diagnostic modalities and challenges in ACC will be addressed.

Integrative Identification of Crucial Genes Associated With Plant Hormone-Mediated Bud Dormancy in Prunus mume

Prunus mume is an important ornamental woody plant with winter-flowering property, which is closely related to bud dormancy. Despite recent scientific headway in deciphering the mechanism of bud dormancy in P. mume, the overall picture of gene co-expression regulating P. mume bud dormancy is still unclear. Here a total of 23 modules were screened by weighted gene co-expression network analysis (WGCNA), of which 12 modules were significantly associated with heteroauxin, abscisic acid (ABA), and gibberellin (GA), including GA1, GA3, and GA4. The yellow module, which was positively correlated with the content of ABA and negatively correlated with the content of GA, was composed of 1,426 genes, among which 156 transcription factors (TFs) were annotated with transcriptional regulation function. An enrichment analysis revealed that these genes are related to the dormancy process and plant hormone signal transduction. Interestingly, the expression trends of PmABF2 and PmABF4 genes, the core members of ABA signal transduction, were positively correlated with P. mume bud dormancy. Additionally, the PmSVP gene had attracted lots of attention because of its co-expression, function enrichment, and expression level. PmABF2, PmABF4, and PmSVP were the genes with a high degree of expression in the co-expression network, which was upregulated by ABA treatment. Our results provide insights into the underlying molecular mechanism of plant hormone-regulated dormancy and screen the hub genes involved in bud dormancy in P. mume.

Identification and Validation of a Novel Prognosis Prediction Model in Adrenocortical Carcinoma by Integrative Bioinformatics Analysis, Statistics, and Machine Learning

Adrenocortical carcinoma (ACC) is a rare malignancy with poor prognosis. Thus, we aimed to establish a potential gene model for prognosis prediction of patients with ACC. First, weighted gene co-expression network (WGCNA) was constructed to screen two key modules (blue: P = 5e-05, R^2 = 0.65; red: P = 4e-06, R^2 = -0.71). Second, 93 survival-associated genes were identified. Third, 11 potential prognosis models were constructed, and two models were further selected. Survival analysis, receiver operating characteristic curve (ROC), Cox regression analysis, and calibrate curve were performed to identify the best model with great prognostic value. Model 2 was further identified as the best model [training set: P < 0.0001; the area under curve (AUC) value was higher than in any other models showed]. We further explored the prognostic values of genes in the best model by analyzing their mutations and copy number variations (CNVs) and found that MKI67 altered the most (12%). CNVs of the 14 genes could significantly affect the relative mRNA expression levels and were associated with survival of ACC patients. Three independent analyses indicated that all the 14 genes were significantly associated with the prognosis of patients with ACC. Six hub genes were further analyzed by constructing a PPI network and validated by AUC and concordance index (C-index) calculation. In summary, we constructed and validated a prognostic multi-gene model and found six prognostic biomarkers, which may be useful for predicting the prognosis of ACC patients.

Identification of collagen genes related to immune infiltration and epithelial-mesenchymal transition in glioma

Background

Gliomas account for the majority of fatal primary brain tumors, and there is much room for research in the underlying pathogenesis, the multistep progression of glioma, and how to improve survival. In our study, we aimed to identify potential biomarkers or therapeutic targets of glioma and study the mechanism underlying the tumor progression.

Methods

We downloaded the microarray datasets (GSE43378 and GSE7696) from the Gene Expression Omnibus (GEO) database. Then, we used weighted gene co-expression network analysis (WGCNA) to screen potential biomarkers or therapeutic targets related to the tumor progression. ESTIMATE (Estimation of STromal and Immune cells in MAlignant Tumors using Expression data) algorithm and TIMER (Tumor Immune Estimation Resource) database were used to analyze the correlation between the selected genes and the tumor microenvironment. Real-time reverse transcription polymerase chain reaction was used to measure the selected gene. Transwell and wound healing assays were used to measure the cell migration and invasion capacity. Western blotting was used to test the expression of epithelial-mesenchymal transition (EMT) related markers.

Results

We identified specific module genes that were positively correlated with the WHO grade but negatively correlated with OS of glioma. Importantly, we identified that 6 collagen genes (COL1A1, COL1A2, COL3A1, COL4A1, COL4A2, and COL5A2) could regulate the immunosuppressive microenvironment of glioma. Moreover, we found that these collagen genes were significantly involved in the EMT process of glioma. Finally, taking COL3A1 as a further research object, the results showed that knockdown of COL3A1 significantly inhibited the migration, invasion, and EMT process of SHG44 and A172 cells.

Conclusions

In summary, our study demonstrated that collagen genes play an important role in regulating the immunosuppressive microenvironment and EMT process of glioma and could serve as potential therapeutic targets for glioma management.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12935-021-01982-0.

ESR1 as a recurrence-related gene in intrahepatic cholangiocarcinoma: a weighted gene coexpression network analysis

Background

Intrahepatic cholangiocarcinoma (iCCA) is the second most common malignant hepatic tumor and has a high postoperative recurrence rate and a poor prognosis. The key roles of most tumor recurrence-associated molecules in iCCA remain unclear. This study aimed to explore hub genes related to the postsurgical recurrence of iCCA.

Method

Differentially expressed genes (DEGs) between iCCA samples and normal liver samples were screened from The Cancer Genome Atlas (TCGA) database and used to construct a weighted gene coexpression network. Module-trait correlations were calculated to identify the key module related to recurrence in iCCA patients. Genes in the key module were subjected to functional enrichment analysis, and candidate hub genes were filtered through coexpression and protein–protein interaction (PPI) network analysis. Validation studies were conducted to detect the “real” hub gene. Furthermore, the biological functions and the underlying mechanism of the real hub gene in iCCA tumorigenesis and progression were determined via in vitro experiments.

Results

A total of 1019 DEGs were filtered and used to construct four coexpression modules. The red module, which showed the highest correlations with the recurrence status, family history, and day to death of patients, was identified as the key module. Gene Ontology (GO) enrichment and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses demonstrated that genes in the red module were enriched in genes and pathways related to tumorigenesis and tumor progression. We performed validation studies and identified estrogen receptor 1 (ESR1), which significantly impacted the prognosis of iCCA patients, as the real hub gene related to the recurrence of iCCA. The in vitro experiments demonstrated that ESR1 overexpression significantly suppressed cell proliferation, migration, and invasion, whereas ESR1 knockdown elicited opposite effects. Further investigation into the mechanism demonstrated that ESR1 acts as a tumor suppressor by inhibiting the JAK/STAT3 signaling pathway.

Conclusions

ESR1 was identified as the real hub gene related to the recurrence of iCCA that plays a critical tumor suppressor role in iCCA progression. ESR1 significantly impacts the prognosis of iCCA patients and markedly suppresses cholangiocarcinoma cell proliferation, migration and invasion by inhibiting JAK/STAT3 signaling pathway.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12935-021-01929-5.

Expression patterns and prognostic value of m6A RNA methylation regulators in adrenocortical carcinoma

Adrenocortical carcinoma (ACC) is considered a rare cancer with poor prognosis. We used public datasets from The Cancer Genome Atlas (TCGA) and Genotype-Tissue Expression (GTEx) databases to assess the relationships between N6-methyladenosine (m6A)-related genes and ACC.We used the Wilcoxon signed-rank test to compare m6A-related gene expression in ACC tissues with that in normal tissues. Then, ACC patients were grouped based on a cluster analysis of m6A-related gene expression. m6A-related genes that were significantly associated with survival were incorporated into a risk signature, and 2 groups were divided according to median risk score. Fisher exact tests were utilized to analyze differences in clinical variables between groups. We compared the overall survival (OS) rates of the groups by means of Kaplan-Meier curves and Cox regression analyses.We found that RBM15, ZC3H3, YTDHF1, YTDHF2, and ALBH5 were overexpressed in ACC and that KIAA1429, YTHDC1, HNRNPC, WTAP, METTL3, and FTO were down regulated in ACC. In addition, membership in cluster 2 or the high-risk group was associated with advanced clinical factors and poor prognosis. The univariable and multivariable Cox regression analyses showed that risk score can be considered an independent prognostic factor for ACC.We found that the expression of m6A-related genes could be used as an independent prognostic factor in ACC. However, the current study has some limitations, and further studies of m6A-related genes in ACC are needed.

The Role of Biomarkers in Adrenocortical Carcinoma: A Review of Current Evidence and Future Perspectives

Adrenocortical carcinoma (ACC) is a rare endocrine malignancy arising from the adrenal cortex often with unexpected biological behavior. It can occur at any age, with two peaks of incidence: in the first and between fifth and seventh decades of life. Although ACC are mostly hormonally active, precursors and metabolites, rather than end products of steroidogenesis are produced by dedifferentiated and immature malignant cells. Distinguishing the etiology of adrenal mass, between benign adenomas, which are quite frequent in general population, and malignant carcinomas with dismal prognosis is often unfeasible. Even after pathohistological analysis, diagnosis of adrenocortical carcinomas is not always straightforward and represents a great challenge for experienced and multidisciplinary expert teams. No single imaging method, hormonal work-up or immunohistochemical labelling can definitively prove the diagnosis of ACC. Over several decades' great efforts have been made in finding novel reliable and available diagnostic and prognostic factors including steroid metabolome profiling or target gene identification. Despite these achievements, the 5-year mortality rate still accounts for approximately 75% to 90%, ACC is frequently diagnosed in advanced stages and therapeutic options are unfortunately limited. Therefore, imperative is to identify new biological markers that can predict patient prognosis and provide new therapeutic options.

Weighted gene co-expression network analysis identifies modules and functionally enriched pathways in the lactation process

The exponential growth in knowledge has resulted in a better understanding of the lactation process in a wide variety of animals. However, the underlying genetic mechanisms are not yet clearly known. In order to identify the mechanisms involved in the lactation process, various mehods, including meta-analysis, weighted gene co-express network analysis (WGCNA), hub genes identification, gene ontology (GO), and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment at before peak (BP), peak (P), and after peak (AP) stages of the lactation processes have been employed. A total of 104, 85, and 26 differentially expressed genes were identified based on PB vs. P, BP vs. AP, and P vs. AP comparisons, respectively. GO and KEGG pathway enrichment analysis revealed that DEGs were significantly enriched in the "ubiquitin-dependent ERAD" and the "chaperone cofactor-dependent protein refolding" in BP vs. P and P vs. P, respectively. WGCNA identified five significant functional modules related to the lactation process. Moreover, GJA1, AP2A2, and NPAS3 were defined as hub genes in the identified modules, highlighting the importance of their regulatory impacts on the lactation process. The findings of this study provide new insights into the complex regulatory networks of the lactation process at three distinct stages, while suggesting several candidate genes that may be useful for future animal breeding programs. Furthermore, this study supports the notion that in combination with a meta-analysis, the WGCNA represents an opportunity to achieve a higher resolution analysis that can better predict the most important functional genes that might provide a more robust bio-signature for phenotypic traits, thus providing more suitable biomarker candidates for future studies.

Overexpression of PKMYT1 Facilitates Tumor Development and Is Correlated with Poor Prognosis in Clear Cell Renal Cell Carcinoma

BACKGROUND Protein kinase membrane-associated tyrosine/threonine (PKMYT1) has been found in many tumors, but its association with clear cell renal cell carcinoma (ccRCC) remains unclear. MATERIAL AND METHODS PKMYT1 expression in ccRCC was examined in the Cancer Genome Atlas (TCGA), Gene Expression Omnibus (GEO), and Tumor Immune Estimation Resource databases. The correlation between PKMYT1 expression and clinicopathological parameters was explored via the chi-square test. Receiver operating characteristic curves were used to estimate the diagnostic performance of PKMYT1. Kaplan-Meier curves, a Cox model, nomogram, time-dependent receiver operating characteristic curves, and decision curve analysis (DCA) were used to evaluate the prognostic value and clinical utility of PKMYT1. Genes coexpressed with PKMYT1 in ccRCC were identified based on TCGA, the gene expression profiling interactive, and cBioPortal. Gene Set Enrichment Analysis revealed biological pathways associated with PKMYT1 in ccRCC. RESULTS Weighted gene coexpression network analysis identified PKMYT1 as one of the genes most significantly correlated with progression of histological grade. PKMYT1 was significantly upregulated in ccRCC compared with normal tissue (P<0.001), with a trend toward differentiating between individuals with ccRCC and those who were healthy (area under the curve=0.942). High PKMYT1 expression was correlated with unsatisfactory survival (hazard ratio=1.67, P=0.001), indicating that it is a risk factor for ccRCC. A nomogram incorporating PKMYT1 level was created and showed a clinical net benefit. PKMYT1 was strongly positively correlated with the anti-silencing function of 1B histone chaperone (ASF1B) gene in ccRCC. CONCLUSIONS PKMYT1 is upregulated in ccRCC and its presence indicates poor prognosis, making it a potential therapeutic target for ccRCC.

MiR-581/SMAD7 Axis Contributes to Colorectal Cancer Metastasis: A Bioinformatic and Experimental Validation-Based Study

Metastasis is a well-known poor prognostic factor and primary cause of mortality in patients with colorectal cancer (CRC). Recently, with the progress of high through-put sequencing, aberrantly expressed non-coding RNAs (ncRNAs) were found to participate in the initiation and development of cancer. However, the mechanisms of ncRNA-mediated regulation of metastasis in CRC remain largely unknown. In this study, we systematically analyzed the expression network of microRNAs (miRNAs) and genes in CRC metastasis using bioinformatics, and discovered that the miR-581/SMAD7 axis could be a potential factor that drives CRC metastasis. A dual luciferase report assay and protein analysis confirmed the binding relationship between miR-581 and SMAD7. Further functional assays revealed that miR-581 inhibition could suppress cell proliferation and induce apoptosis in SW480 cells. Up-regulation or down-regulation of miR-581 could both affect cell invasion capacity and modulate epithelial to mesenchymal transition (EMT) via a SMAD7/TGFβ signaling pathway. In conclusion, our findings elucidated that miR-581/SMAD7 could be essential for CRC metastasis, and may serve as a potential therapeutic target for CRC patients.

CCNB2, TOP2A, and ASPM Reflect the Prognosis of Hepatocellular Carcinoma, as Determined by Weighted Gene Coexpression Network Analysis

Background

Hepatocellular carcinoma (HCC) is characterized by increased mortality and poor prognosis. We aimed to identify potential prognostic markers by weighted gene coexpression network analysis (WGCNA), to assist clinical outcome prediction and improve treatment decisions for HCC patients.

Methods

Prognosis-related gene modules were first established by WGCNA. Venn diagrams obtained intersection genes of module genes and differentially expressed genes. The Kaplan-Meier overall survival curves and disease-free survival curves of intersection genes were further analyzed on the Gene Expression Profiling Interactive Analysis website. Chi-square tests were performed to explore the associations between prognostic gene expressions and clinicopathological features.

Results

CCNB2, TOP2A, and ASPM were identified as both prognosis-related genes and differentially expressed genes. TOP2A (HR: 1.7, P = 0.003) and ASPM (HR: 1.8, P < 0.001) exhibited a significant difference between the high- and low-expression groups in the overall survival analysis, while CCNB2 (HR: 1.4, P = 0.052) was not statistically significant. CCNB2 (HR: 1.5, P = 0.006), TOP2A (HR: 1.7, P < 0.001), and ASPM (HR: 1.6, P = 0.003) were all statistically significant in the disease-free survival analysis. All three genes were significantly associated with race and fetoprotein values (P < 0.05). CCNB2 expression was associated with tumor stage (P = 0.01), and ASPM expression was associated with new tumor events (P = 0.03).

Conclusion

Overexpression of CCNB2, TOP2A, and ASPM are associated with poor prognosis, and these genes could serve as potential prognostic markers and therapeutic targets for HCC.

Identification of the key genes and pathways involved in the tumorigenesis and prognosis of kidney renal clear cell carcinoma

Kidney renal clear cell carcinoma (KIRC) is the most common renal cell carcinoma (RCC). However, patients with KIRC usually have poor prognosis due to limited biomarkers for early detection and prognosis prediction. In this study, we analysed key genes and pathways involved in KIRC from an array dataset including 26 tumour and 26 adjacent normal tissue samples. Weighted gene co-expression network analysis (WGCNA) was performed with the WGCNA package, and 20 modules were characterized as having the highest correlation with KIRC. The upregulated genes in the tumour samples are involved in the innate immune response, whereas the downregulated genes contribute to the cellular catabolism of glucose, amino acids and fatty acids. Furthermore, the key genes were evaluated through a protein-protein interaction (PPI) network combined with a co-expression network. The comparatively lower expression of AGXT, PTGER3 and SLC12A3 in tumours correlates with worse prognosis in KIRC patients, while higher expression of ALOX5 predicts reduced survival. Our integrated analysis illustrated the hub genes involved in KIRC tumorigenesis, shedding light on the development of prognostic markers. Further understanding of the function of the identified KIRC hub genes could provide deep insights into the molecular mechanisms of KIRC.

Weighted gene co-expression network analysis identified MYL9 and CNN1 are associated with recurrence in colorectal cancer

Colorectal cancer (CRC) is one of the most common carcinomas and the fourth leading cause of cancer-related death worldwide. One of the obstacles in the successful treatment of CRC is a high rate of recurrence. We aimed to construct weighted gene co-expression network analysis (WGCNA) to identify key modules and hub genes in association with recurrence in CRC patients. We firstly used the microarray data, GSE41258, to construct a co-expression network and identify gene modules. Furthermore, protein and protein interaction (PPI) network was also performed to screen hub genes. To validate the hub genes, an independent dataset GSE17536 was used for survival analyses. Additionally, another two databases were also performed to investigate the survival rates and expression levels of hub genes. Gene set enrichment analyses (GSEA) combined with gene ontology (GO) were performed to further explore function and mechanisms. In our study, the midnightblue module was identified to be significant, 15 hub genes were screened, four of which were identified as hub nodes in the PPI network. In the test dataset, we found higher expression of MYL9 and CNN1 were significantly associated with shorter survival time of CRC patients. GO analyses showed that MYL9 and CNN1 were enriched in “muscle system process” and “cytoskeletal protein binding”. GSEA found the two hub genes were enriched in “pathways in cancer” and “calcium signaling pathway”. In conclusion, our study demonstrated that MYL9 and CNN1 were hub genes associated with the recurrence of CRC, which may contribute to the improvement of recurrence-free survival time of CRC patients.

In silico identification of EP400 and TIA1 as critical transcription factors involved in human hepatocellular carcinoma relapse

Hepatocellular carcinoma (HCC) is the second leading cause of cancer-associated mortality worldwide. Transcription factors (TFs) are crucial proteins that regulate gene expression during cancer progression; however, the roles of TFs in HCC relapse remain unclear. To identify the TFs that drive HCC relapse, the present study constructed co-expression network and identified the Tan module the most relevant to HCC relapse. Numerous hub TFs (highly connected) were subsequently obtained from the Tan module according to the intra-module connectivity and the protein-protein interaction network connectivity. Next, E1A-binding protein p400 (EP400) and TIA1 cytotoxic granule associated RNA binding protein (TIA1) were identified as hub TFs differentially connected between the relapsed and non-relapsed subnetworks. In addition, zinc finger protein 143 (ZNF143) and Yin Yang 1 (YY1) were also identified by using the plugin iRegulon in Cytoscape as master upstream regulatory elements, which could potentially regulate expression of the genes and TFs of the Tan module, respectively. The Kaplan-Meier (KM) curves obtained from KMplot and Gene Expression Profiling Interactive Analysis tools confirmed that the high expression of EP400 and TIA1 were significantly associated with shorter relapse-free survival and disease-free survival of patients with HCC. Furthermore, the KM curves from the UALCAN database demonstrated that high EP400 expression significantly reduced the overall survival of patients with HCC. EP400 and TIA1 may therefore serve as potential prognostic and therapeutic biomarkers.

ACAT1 and Metabolism-Related Pathways Are Essential for the Progression of Clear Cell Renal Cell Carcinoma (ccRCC), as Determined by Co-expression Network Analysis

Kidney cancer ranks as one of the top 10 causes of cancer death; this cancer is difficult to detect, difficult to treat, and poorly understood. The most common subtype of kidney cancer is clear cell renal cell carcinoma (ccRCC) and its progression is influenced by complex gene interactions. Few clinical studies have investigated the molecular markers associated with the progression of ccRCC. In this study, we collected microarray profiles of 72 ccRCCs and matched normal samples to identify differentially expressed genes (DEGs). Then a weighted gene co-expression network analysis (WGCNA) was conducted to identify co-expressed gene modules. By relating all co-expressed modules to clinical features, we found that the brown module and Fuhrman grade had the highest correlation (r = -0.8, p = 1e-09). Thus, the brown module was regarded as a clinically significant module and subsequently analyzed. Functional annotation showed that the brown module focused on metabolism-related biological processes and pathways, such as fatty acid oxidation and amino acid metabolism. We then performed a protein-protein interaction (PPI) network to identify the hub nodes in the brown module. It is worth noting that only one candidate, acetyl-CoA acetyltransferase (ACAT1), was considered to be the final target most relevant to the Fuhrman grade of ccRCC, by applying the intersection of hub genes in the co-expressed network and the PPI network. ACAT1 was subsequently validated using another two external microarray datasets and the TCGA dataset. Intriguingly, validation results indicated that ACAT1 was negatively correlated with four grades of ccRCC, which was also consistent with our results from qRT-PCR analysis and immunohistochemistry staining of clinical samples. Overexpression of ACAT1 inhibited the proliferation and migration of human ccRCC cells in vitro. In addition, the Kaplan-Meier survival curve showed that patients with a lower expression of ACAT1 showed a significantly lower overall survival rate and disease-free survival rate, indicating that ACAT1 could act as a prognostic and recurrence/progression biomarker of ccRCC. In summary, we found and confirmed that ACAT1 might help to identify the progression of ccRCC, which might have important clinical implications for enhancing risk stratification, therapeutic decision, and prognosis prediction in ccRCC patients.

Identification of a blood-based 12-gene signature that predicts the severity of coronary artery stenosis: An integrative approach based on gene network construction, Support Vector Machine algorithm, and multi-cohort validation

BACKGROUND AND AIMS

We aimed to identify a blood-based gene expression score (GES) to predict the severity of coronary artery stenosis in patients with known or suspected coronary artery disease (CAD) by integrative use of gene network construction, Support Vector Machine (SVM) algorithm, and multi-cohort validation.

METHODS

In the discovery phase, a public blood-based microarray dataset of 110 patients with known CAD was analyzed by weighted gene coexpression network analysis and protein-protein interaction network analysis to identify candidate hub genes. In the training set with 151 CAD patients, bioinformatically identified hub genes were experimentally verified by real-time polymerase chain reaction, and statistically filtered with the SVM algorithm to develop a GES. Internal and external validation of GES was performed in patients with suspected CAD from two validation cohorts (n = 209 and 206).

RESULTS

The discovery phase screened 15 network-centric hub genes significantly correlated with the Duke CAD Severity Index. In the training cohort, 12 of 15 hub genes were filtered to construct a blood-based GES12, which showed good discrimination for higher modified Gensini scores (AUC: 0.798 and 0.812), higher Sullivan Extent scores (AUC: 0.776 and 0.778), and the presence of obstructive CAD (AUC: 0.834 and 0.792) in two validation cohorts. A nomogram comprising GES12, smoking status, hypertension status, low density lipoprotein cholesterol level, and body mass index further improved performance, with respect to discrimination, risk classification, and clinical utility, for prediction of coronary stenosis severity.

CONCLUSIONS

GES12 is useful in predicting the severity of coronary artery stenosis in patients with known or suspected CAD.

Screening and identification of key biomarkers in adrenocortical carcinoma based on bioinformatics analysis

Adrenocortical carcinoma (ACC) is a rare malignancy with a poor prognosis. The presently available understanding of the pathogenesis of ACC is incomplete and the treatment options for patients with ACC are limited. Gene marker identification is required for accurate and timely diagnosis of the disease. In order to identify novel candidate genes associated with the occurrence and progression of ACC, the microarray datasets, GSE12368 and GSE19750, were obtained from Gene Expression Omnibus. Differentially expressed genes (DEGs) were identified, and functional enrichment analysis was performed. A protein-protein interaction network (PPI) was constructed to identify significantly altered modules, and module analysis was performed using Search Tool for the Retrieval of Interacting Genes and Cytoscape. A total of 228 DEGs were screened, consisting of 29 up and 199 downregulated genes. The enriched functions and pathways of the DEGs primarily included 'cell division', 'regulation of transcription involved in G1/S transition of mitotic cell cycle', 'G1/S transition of mitotic cell cycle', 'p53 signaling pathway' and 'oocyte meiosis'. A total of 14 hub genes were identified, and biological process analysis revealed that these genes were significantly enriched in cell division and mitotic cell cycle. Furthermore, survival analysis revealed that AURKA, TYMS, GINS1, RACGAP1, RRM2, EZH2, ZWINT, CDK1, CCNB1, NCAPG and TPX2 may be involved in the tumorigenesis, progression or prognosis of ACC. In conclusion, the 14 hub genes identified in the present study may aid researchers in elucidating the molecular mechanisms associated with the tumorigenesis and progression of ACC, and may be powerful and promising candidate biomarkers for the diagnosis and treatment of ACC.

Statistical genomics in rare cancer

Rare cancers make of more than 20% of cancer cases. Due to the rare nature, less research has been conducted on rare cancers resulting in worse outcomes for patients with rare cancers compared to common cancers. The ability to study rare cancers is impaired by the ability to collect a large enough set of patients to complete an adequately powered genomic study. In this manuscript we outline analytical approaches and public genomic datasets that have been used in genomic studies of rare cancers. These statistical analysis approaches and study designs include: gene set / pathway analyses, pedigree and consortium studies, meta-analysis or horizontal integration, and integration of multiple types of genomic information or vertical integration. We also discuss some of the publicly available resources that can be leveraged in rare cancer genomic studies.

Identification of Hub Genes Associated With Progression and Prognosis in Patients With Bladder Cancer

Given that most bladder cancers (BCs) are diagnosed in advanced stages with poor prognosis, this study aims to find novel biomarkers associated with the progression and prognosis in patients with BC. 1,779 differentially expressed genes (DEGs) between BC samples and normal bladder tissues were identified in total. Then, 24 DEGs were regarded as candidate hub genes by constructing a protein–protein interaction (PPI) network and a random forest model. Among them, six genes (BUB1B, CCNB1, CDK1, ISG15, KIF15, and RAD54L) were eventually identified by using five analysis methods (one-way Analysis of Variance analysis, spearman correlation analysis, distance correlation analysis, receiver operating characteristic curve, and expression values comparison), which were correlated with the progression and prognosis of BC. Moreover, the validation of hub genes was conducted based on GSE13507, Oncomine, and CBioPortal. Results of univariate Cox regression analysis showed that the expression levels of all the hub genes were influence features of overall survival (OS) and cancer specific survival (CSS) based on GSE13507, and we further established a six-gene signature based on the expression levels of the six genes and their Cox regression coefficients. This signature showed good potential for clinical application suggested by survival analysis (OS: Hazard Ratio = 0.484, 95%CI: 0.298–0.786; P = 0.0034; CSS: Hazard Ratio = 0.244, 95%CI: 0.121–0.493, P < 0.0001) and decision curve analysis. In conclusion, our study indicates that six hub genes have great predictive value for the prognosis and progression of BC and may contribute to the exploration of further basic and clinical research of BC.

Identification of biomarkers of chromophobe renal cell carcinoma by weighted gene co-expression network analysis

Background

Chromophobe renal cell carcinoma (ChRCC) is the second common subtype of non-clear cell renal cell carcinoma (nccRCC), which accounting for 4–5% of renal cell carcinoma (RCC). However, there is no effective bio-marker to predict clinical outcomes of this malignant disease. Bioinformatic methods may provide a feasible potential to solve this problem.

Methods

In this study, differentially expressed genes (DEGs) of ChRCC samples on The Cancer Genome Atlas database were filtered out to construct co-expression modules by weighted gene co-expression network analysis and the key module were identified by calculating module-trait correlations. Functional analysis was performed on the key module and candidate hub genes were screened out by co-expression and MCODE analysis. Afterwards, real hub genes were filter out in an independent dataset GSE15641 and validated by survival analysis.

Results

Overall 2215 DEGs were screened out to construct eight co-expression modules. Brown module was identified as the key module for the highest correlations with pathologic stage, neoplasm status and survival status. 29 candidate hub genes were identified. GO and KEGG analysis demonstrated most candidate genes were enriched in mitotic cell cycle. Three real hub genes (SKA1, ERCC6L, GTSE-1) were selected out after mapping candidate genes to GSE15641 and two of them (SKA1, ERCC6L) were significantly related to overall survivals of ChRCC patients.

Conclusions

In summary, our findings identified molecular markers correlated with progression and prognosis of ChRCC, which might provide new implications for improving risk evaluation, therapeutic intervention, and prognosis prediction in ChRCC patients.

Electronic supplementary material

The online version of this article (10.1186/s12935-018-0703-z) contains supplementary material, which is available to authorized users.

Identification of Crucial Genes and Pathways in Human Arrhythmogenic Right Ventricular Cardiomyopathy by Coexpression Analysis

As one common disease causing young people to die suddenly due to cardiac arrest, arrhythmogenic right ventricular cardiomyopathy (ARVC) is a disorder of heart muscle whose progression covers one complicated gene interaction network that influence the diagnosis and prognosis of it. In our research, differentially expressed genes (DEGs) were screened, and we established a weighted gene coexpression network analysis (WGCNA) and gene set net correlations analysis (GSNCA) for identifying crucial genes as well as pathways related to ARVC pathogenic mechanism (n = 12). In the research, the results demonstrated that there were 619 DEGs in total between non-failing donor myocardial samples and ARVC tissues (FDR < 0.05). WGCNA analysis identified the two gene modules (brown and turquoise) as being most significantly associated with ARVC state. Then the ARVC-related four key biological pathways (cytokine–cytokine receptor interaction, chemokine signaling pathway, neuroactive ligand receptor interaction, and JAK-STAT signaling pathway) and four hub genes (CXCL2, TNFRSF11B, LIFR, and C5AR1) in ARVC samples were further identified by GSNCA method. Finally, we used t-test and receiver operating characteristic (ROC) curves for validating hub genes, results showed significant differences in t-test and their AUC areas all greater than 0.8. Together, these results revealed that the new four hub genes as well as key pathways that might be involved into ARVC diagnosis. Even though further experimental validation is required for the implication by association, our findings demonstrate that the computational methods based on systems biology might complement the traditional gene-wide approaches, as such, might offer a new insight in therapeutic intervention within rare diseases of people like ARVC.