Single Cell Gene Co-Expression Network Reveals FECH/CROT Signature as a Prognostic Marker

Systematic review and integrated analysis of prognostic gene signatures for prostate cancer patients

Prostate cancer (PC) is one of the most common cancers in men and becoming the second leading cause of cancer fatalities. At present, the lack of effective strategies for prognosis of PC patients is still a problem to be solved. Therefore, it is significant to identify potential gene signatures for PC patients' prognosis. Here, we summarized 71 different prognostic gene signatures for PC and concluded 3 strategies for signature construction after extensive investigation. In addition, 14 genes frequently appeared in 71 different gene signatures, which enriched in mitotic and cell cycle. This review provides extensive understanding and integrated analysis of current prognostic signatures of PC, which may help researchers to construct gene signatures of PC and guide future clinical treatment.

Systematic single cell RNA sequencing analysis reveals unique transcriptional regulatory networks of Atoh1-mediated hair cell conversion in adult mouse cochleae

Regeneration of mammalian cochlear hair cells (HCs) by modulating molecular pathways or transcription factors is a promising approach to hearing restoration; however, immaturity of the regenerated HCs in vivo remains a major challenge. Here, we analyzed a single cell RNA sequencing (scRNA-seq) dataset during Atoh1-induced supporting cell (SC) to hair cell (HC) conversion in adult mouse cochleae (Yamashita et al. (2018)) using multiple high-throughput sequencing analytical tools (WGCNA, SCENIC, ARACNE, and VIPER). Instead of focusing on differentially expressed genes, we established independent expression modules and confirmed the existence of multiple conversion stages. Gene regulatory network (GRN) analysis uncovered previously unidentified key regulators, including Nhlh1, Lhx3, Barhl1 and Nfia, that guide converted HC differentiation. Comparison of the late-stage converted HCs with the scRNA-seq data from neonatal mouse cochleae (Kolla et al. (2020)) revealed that they closely resemble postnatal day 1 wild-type OHCs, in contrast to other developmental stages. Using ARACNE and VIPER, we discovered multiple key regulators likely to promote conversion to a more mature OHC-like state, including Zbtb20, Nfia, Zmiz1, Gm14418, Bhlhe40, Six2, Fosb and Klf9. Our findings provide insights into the regulation of HC regeneration in adult mammalian cochleae in vivo and demonstrate an approach for analyzing GRNs in large scRNA-seq datasets.

Construction and validation of a fatty acid metabolism-related gene signature for predicting prognosis and therapeutic response in patients with prostate cancer

Background

Reprogramming of fatty acid metabolism is a newly-identified hallmark of malignancy. However, no studies have systematically investigated the fatty acid metabolism related-gene set in prostate cancer (PCa).

Methods

A cohort of 381 patients with gene expression and clinical data from The Cancer Genome Atlas was used as the training set, while another cohort of 90 patients with PCa from GEO (GSE70769) was used as the validation set. Differentially expressed fatty acid metabolism-related genes were subjected to least absolute shrinkage and selection operator (LASSO)-Cox regression to establish a fatty acid metabolism-related risk score. Associations between the risk score and clinical characteristics, immune cell infiltration, tumor mutation burden (TMB), tumor immune dysfunction and exclusion (TIDE) score, and response to chemotherapy were analyzed. Finally, the expression level of genes included in the model was validated using real-time PCR.

Results

A prognostic risk model based on five fatty acid metabolism related genes (ALDH1A1, CPT1B, CA2, CROT, and NUDT19) were constructed. Tumors with higher risk score were associated with larger tumor size, lymph node involvement, higher Gleason score, and poorer biochemical recurrence (BCR)-free survival. Furthermore, the high- and low-risk tumors exhibited distinct immune cell infiltration features and immune-related pathway activation. High-risk tumors were associated with favorable response to immunotherapy as indicated by high TMB and low TIDE score, but poor response to bicalutamide and docetaxel chemotherapy.

Conclusion

This study established a fatty acid metabolism-related gene signature which was predictive of BCR and response to chemotherapy and immunotherapy, providing a novel therapeutic biomarker for PCa.

Co-expression of fibrotic genes in inflammatory bowel disease; A localized event?

Introduction

Extracellular matrix turnover, a ubiquitous dynamic biological process, can be diverted to fibrosis. The latter can affect the intestine as a serious complication of Inflammatory Bowel Diseases (IBD) and is resistant to current pharmacological interventions. It embosses the need for out-of-the-box approaches to identify and target molecular mechanisms of fibrosis.

Methods and results

In this study, a novel mRNA sequencing dataset of 22 pairs of intestinal biopsies from the terminal ileum (TI) and the sigmoid of 7 patients with Crohn's disease, 6 with ulcerative colitis and 9 control individuals (CI) served as a validation cohort of a core fibrotic transcriptomic signature (FIBSig), This signature, which was identified in publicly available data (839 samples from patients and healthy individuals) of 5 fibrotic disorders affecting different organs (GI tract, lung, skin, liver, kidney), encompasses 241 genes and the functional pathways which derive from their interactome. These genes were used in further bioinformatics co-expression analyses to elucidate the site-specific molecular background of intestinal fibrosis highlighting their involvement, particularly in the terminal ileum. We also confirmed different transcriptomic profiles of the sigmoid and terminal ileum in our validation cohort. Combining the results of these analyses we highlight 21 core hub genes within a larger single co-expression module, highly enriched in the terminal ileum of CD patients. Further pathway analysis revealed known and novel inflammation-regulated, fibrogenic pathways operating in the TI, such as IL-13 signaling and pyroptosis, respectively.

Discussion

These findings provide a rationale for the increased incidence of fibrosis at the terminal ileum of CD patients and highlight operating pathways in intestinal fibrosis for future evaluation with mechanistic and translational studies.

scGENA: A Single-Cell Gene Coexpression Network Analysis Framework for Clustering Cell Types and Revealing Biological Mechanisms

Single-cell RNA-sequencing (scRNA-seq) is a recent high-throughput technique that can measure gene expression, reveal cell heterogeneity, rare and complex cell populations, and discover cell types and their relationships. The analysis of scRNA-seq data is challenging because of transcripts sparsity, replication noise, and outlier cell populations. A gene coexpression network (GCN) analysis effectively deciphers phenotypic differences in specific states by describing gene–gene pairwise relationships. The underlying gene modules with different coexpression patterns partially bridge the gap between genotype and phenotype. This study presents a new framework called scGENA (single-cell gene coexpression network analysis) for GCN analysis based on scRNA-seq data. Although there are several methods for scRNA-seq data analysis, we aim to build an integrative pipeline for several purposes that cover primary data preprocessing, including data exploration, quality control, normalization, imputation, and dimensionality reduction of clustering as downstream of GCN analysis. To demonstrate this integrated workflow, an scRNA-seq dataset of the human diabetic pancreas with 1600 cells and 39,851 genes was implemented to perform all these processes in practice. As a result, scGENA is demonstrated to uncover interesting gene modules behind complex diseases, which reveal biological mechanisms. scGENA provides a state-of-the-art method for gene coexpression analysis for scRNA-seq data.

Metabolic-related gene signature model forecasts biochemical relapse in primary prostate cancer

Metabolism plays an important role in the pathogenesis of prostate cancer (PCa). Hence, we explored candidate metabolic-related genes attributed to biochemical relapse (BCR) of PCa. Gene expression profile and clinical parameters were downloaded from GSE70769 as a "training set". Using univariate Cox and LASSO-COX regression models, risk scores (RSs) were constructed. Kaplan-Meier (K-M) survival and time-dependent receiver operating characteristic (t-ROC) curves were employed. Univariate and multivariate Cox models were utilized to validate prognostic factors for biochemical relapse-free survival (BCRFS). Nomogram was plotted to facilitate clinical application. The dataset obtained from GSE70768 served as "validation set". RSs were constructed by using 7 metabolic-related genes. RSs could significantly predict 1, 3, 5-year BCRFS (AUCs for training set: 0.810-0.836; AUC for validation set: 0.673-0.827). Nomograms could effectively predicted BCRFS (training set: C-index=0.831; validation set: C-index=0.737). RSs model is an independent prognostic factor for BCR, holding greater predictive value than traditional clinicopathological parameters. Clinical Relevance- We built the prognostic nomogram based on metabolic-related gene signatures and clinicopathological features. The nomogram might further optimize biochemical relapse risk stratification for prostate cancer patients with crucial accuracy.

Identification of co-expression hub genes for ferroptosis in kidney renal clear cell carcinoma based on weighted gene co-expression network analysis and The Cancer Genome Atlas clinical data

Renal clear cell carcinoma (KIRC) is one of the most common tumors worldwide and has a high mortality rate. Ferroptosis is a major mechanism of tumor occurrence and development, as well as important for prognosis and treatment of KIRC. Here, we conducted bioinformatics analysis to identify KIRC hub genes that target ferroptosis. By Weighted gene co-expression network analysis (WGCNA), 11 co-expression-related genes were screened out. According to Kaplan Meier's survival analysis of the data from the gene expression profile interactive analysis database, it was identified that the expression levels of two genes, PROM2 and PLIN2, are respectively related to prognosis. In conclusion, our findings indicate that PROM2 and PLIN2 may be effective new targets for the treatment and prognosis of KIRC.

Guidelines for bioinformatics of single-cell sequencing data analysis in Alzheimer's disease: review, recommendation, implementation and application

Alzheimer's disease (AD) is the most common form of dementia, characterized by progressive cognitive impairment and neurodegeneration. Extensive clinical and genomic studies have revealed biomarkers, risk factors, pathways, and targets of AD in the past decade. However, the exact molecular basis of AD development and progression remains elusive. The emerging single-cell sequencing technology can potentially provide cell-level insights into the disease. Here we systematically review the state-of-the-art bioinformatics approaches to analyze single-cell sequencing data and their applications to AD in 14 major directions, including 1) quality control and normalization, 2) dimension reduction and feature extraction, 3) cell clustering analysis, 4) cell type inference and annotation, 5) differential expression, 6) trajectory inference, 7) copy number variation analysis, 8) integration of single-cell multi-omics, 9) epigenomic analysis, 10) gene network inference, 11) prioritization of cell subpopulations, 12) integrative analysis of human and mouse sc-RNA-seq data, 13) spatial transcriptomics, and 14) comparison of single cell AD mouse model studies and single cell human AD studies. We also address challenges in using human postmortem and mouse tissues and outline future developments in single cell sequencing data analysis. Importantly, we have implemented our recommended workflow for each major analytic direction and applied them to a large single nucleus RNA-sequencing (snRNA-seq) dataset in AD. Key analytic results are reported while the scripts and the data are shared with the research community through  GitHub. In summary, this comprehensive review provides insights into various approaches to analyze single cell sequencing data and offers specific guidelines for study design and a variety of analytic directions. The review and the accompanied software tools will serve as a valuable resource for studying cellular and molecular mechanisms of AD, other diseases, or biological systems at the single cell level.

The miR-33a-5p/CROT axis mediates ovarian cancer cell behaviors and chemoresistance via the regulation of the TGF-β signal pathway

Due to the lack of symptoms and detection biomarkers at the early stage, most patients with ovarian cancer (OC) are diagnosed at an advanced stage and often face chemoresistance and relapse. Hence, defining detection biomarkers and mechanisms of chemoresistance is imperative. A previous report of a cDNA microarray analysis shows a potential association of carnitine O-octanoyltransferase (CROT) with taxane resistance but the biological function of CROT in OC remains unknown. The current study explored the function and regulatory mechanism of CROT on cellular behavior and paclitaxel (PTX)-resistance in OC. We found that CROT was downregulated in OC tissues and PTX-resistant cells. Furthermore, CROT expression was negatively correlated with the prognosis of OC patients. Overexpression of CROT inhibited the OC cell proliferation, migration, invasion, and colony formation, arrested the cell cycle at the G2/M phase, and promoted cell apoptosis. In addition, miR-33a-5p bound directly to the 3'UTR of CROT to negatively regulate the expression of CROT and promoted OC cell growth. Finally, overexpression of CROT decreased the phosphorylation of Smad2, whereas knockdown of CROT increased the nuclear translocation of Smad2 and Smad4, two transducer proteins of TGF-β signaling, indicating that CROT is a tumor suppressor which mediates OC cell behaviors through the TGF-β signaling pathway. Thus, targeting the miR-33a-5p/CROT axis may have clinical potential for the treatment of patients with OC.

Silent transcription intervals and translational bursting lead to diverse phenotypic switching

For complex process of gene expression, we use theoretical analysis and stochastic simulations to study the phenotypic diversity induced by silent transcription intervals and translational bursting.

Epigenetic loss of heterogeneity from low to high grade localized prostate tumours

Identifying precise molecular subtypes attributable to specific stages of localized prostate cancer has proven difficult due to high levels of heterogeneity. Bulk assays represent a population-average, which mask the heterogeneity that exists at the single-cell level. In this work, we sequence the accessible chromatin regions of 14,424 single-cells from 18 flash-frozen prostate tumours. We observe shared chromatin features among low-grade prostate cancer cells are lost in high-grade tumours. Despite this loss, high-grade tumours exhibit an enrichment for FOXA1, HOXB13 and CDX2 transcription factor binding sites, indicating a shared trans-regulatory programme. We identify two unique genes encoding neuronal adhesion molecules that are highly accessible in high-grade prostate tumours. We show NRXN1 and NLGN1 expression in epithelial, endothelial, immune and neuronal cells in prostate cancer using cyclic immunofluorescence. Our results provide a deeper understanding of the active gene regulatory networks in primary prostate tumours, critical for molecular stratification of the disease.

Identification of hub genes in colorectal cancer based on weighted gene co-expression network analysis and clinical data from The Cancer Genome Atlas

Colorectal cancer (CRC) is one of the most common tumors worldwide and is associated with high mortality. Here we performed bioinformatics analysis, which we validated using immunohistochemistry in order to search for hub genes that might serve as biomarkers or therapeutic targets in CRC. Based on data from The Cancer Genome Atlas (TCGA), we identified 4832 genes differentially expressed between CRC and normal samples (1562 up-regulated and 3270 down-regulated in CRC). Gene ontology (GO) analysis showed that up-regulated genes were enriched mainly in organelle fission, cell cycle regulation, and DNA replication; down-regulated genes were enriched primarily in the regulation of ion transmembrane transport and ion homeostasis. Weighted gene co-expression network analysis (WGCNA) identified eight gene modules that were associated with clinical characteristics of CRC patients, including brown and blue modules that were associated with cancer onset. Analysis of the latter two hub modules revealed the following six hub genes: adhesion G protein-coupled receptor B3 (BAI3, also known as ADGRB3), cyclin F (CCNF), cytoskeleton-associated protein 2 like (CKAP2L), diaphanous-related formin 3 (DIAPH3), oxysterol binding protein-like 3 (OSBPL3), and RERG-like protein (RERGL). Expression levels of these hub genes were associated with prognosis, based on Kaplan–Meier survival analysis of data from the Gene Expression Profiling Interactive Analysis database. Immunohistochemistry of CRC tumor tissues confirmed that OSBPL3 is up-regulated in CRC. Our findings suggest that CCNF, DIAPH3, OSBPL3, and RERGL may be useful as therapeutic targets against CRC. BAI3 and CKAP2L may be novel biomarkers of the disease.

Identification of a three-gene-based prognostic model in multiple myeloma using bioinformatics analysis

Background

Multiple myeloma (MM), the second most hematological malignancy, has high incidence and remains incurable till now. The pathogenesis of MM is poorly understood. This study aimed to identify novel prognostic model for MM on gene expression profiles.

Methods

Gene expression datas of MM (GSE6477, GSE136337) were downloaded from Gene Expression Omnibus (GEO) database. The differentially expressed genes (DEGs) in GSE6477 between case samples and normal control samples were screened by the limma package. Meanwhile, enrichment analysis was conducted, and a protein-protein interaction (PPI) network of these DEGs was established by STRING and cytoscape software. Co-expression modules of genes were built by Weighted Correlation Network Analysis (WGCNA). Key genes were identified both from hub genes and the DEGs. Univariate and multivariate Cox congression were performed to screen independent prognostic genes to construct a predictive model. The predictive power of the model was evaluated by Kaplan–Meier curve and time-dependent receiver operating characteristic (ROC) curves. Finally, univariate and multivariate Cox regression analyse were used to investigate whether the prognostic model could be independent of other clinical parameters.

Results

GSE6477, including 101 case and 15 normal control, were screened as the datasets. A total of 178 DEGs were identified, including 59 up-regulated and 119 down-regulated genes. In WGCNA analysis, module black and module purple were the most relevant modules with cancer traits, and 92 hub genes in these two modules were selected for further analysis. Next, 47 genes were chosen both from the DEGs and hub genes as key genes. Three genes (LYVE1, RNASE1, and RNASE2) were finally screened by univariate and multivariate Cox regression analyses and used to construct a risk model. In addition, the three-gene prognostic model revealed independent and accurate prognostic capacity in relation to other clinical parameters for MM patients.

Conclusion

In summary, we identified and constructed a three-gene-based prognostic model that could be used to predict overall survival of MM patients.

Inferring Cell Subtypes and LncRNA Function by a Cell-Specific CeRNA Network in Breast Cancer

Single-cell RNA sequencing is a powerful tool to explore the heterogeneity of breast cancer. The identification of the cell subtype that responds to estrogen has profound significance in breast cancer research and treatment. The transcriptional regulation of estrogen is an intricate network involving crosstalk between protein-coding and non-coding RNAs, which is still largely unknown, particularly at the single cell level. Therefore, we proposed a novel strategy to specify cell subtypes based on a cell-specific ceRNA network (CCN). The CCN was constructed by integrating a cell-specific RNA-RNA co-expression network (RCN) with an existing ceRNA network. The cell-specific RCN was built based on single cell expression profiles with predefined reference cells. Heterogeneous cell subtypes were inferred by enriching RNAs in CCN to the estrogen response hallmark. Edge biomarkers were identified in the early estrogen response subtype. Topological analysis revealed that NEAT1 was a hub lncRNA for the early response subtype, and its ceRNAs could predict patient survival. Another hub lncRNA, DLEU2, could potentially be involved in GPCR signaling, based on CCN. The CCN method that we proposed here facilitates the inference of cell subtypes from a network perspective and explores the function of hub lncRNAs, which are promising targets for RNA-based therapeutics.

Single-cell network biology for resolving cellular heterogeneity in human diseases

Understanding cellular heterogeneity is the holy grail of biology and medicine. Cells harboring identical genomes show a wide variety of behaviors in multicellular organisms. Genetic circuits underlying cell-type identities will facilitate the understanding of the regulatory programs for differentiation and maintenance of distinct cellular states. Such a cell-type-specific gene network can be inferred from coregulatory patterns across individual cells. Conventional methods of transcriptome profiling using tissue samples provide only average signals of diverse cell types. Therefore, reconstructing gene regulatory networks for a particular cell type is not feasible with tissue-based transcriptome data. Recently, single-cell omics technology has emerged and enabled the capture of the transcriptomic landscape of every individual cell. Although single-cell gene expression studies have already opened up new avenues, network biology using single-cell transcriptome data will further accelerate our understanding of cellular heterogeneity. In this review, we provide an overview of single-cell network biology and summarize recent progress in method development for network inference from single-cell RNA sequencing (scRNA-seq) data. Then, we describe how cell-type-specific gene networks can be utilized to study regulatory programs specific to disease-associated cell types and cellular states. Moreover, with scRNA data, modeling personal or patient-specific gene networks is feasible. Therefore, we also introduce potential applications of single-cell network biology for precision medicine. We envision a rapid paradigm shift toward single-cell network analysis for systems biology in the near future.

Gene regulatory networks reconstructed from single-cell RNA sequencing datasets are allowing researchers to better understand the molecular circuits and cell states that contribute to complex human disease. Junha Cha and Insuk Lee from Yonsei University in Seoul, South Korea, review the concept of ‘single-cell network biology’, which involves using computational algorithms on genetic expression data from thousands of cells to infer functional interactions in various biological contexts. This systems biology approach to analyzing the profiles of messenger RNA in single cells is helping researchers discover new signaling pathways that could serve as disease biomarkers or therapeutic targets. In the future, patient-specific models of personal gene networks could explain why certain genetic variants affect disease risk. This research could also eventually lead to new types of individualized medical treatments.

Construction of a 6-gene prognostic signature to assess prognosis of patients with pancreatic cancer

Pancreatic cancer (PaCa) is one of the most fatal cancers in the world. Although great efforts have made to explore the mechanisms of PaCa oncogenesis, the prognosis of PaCa patients is still unsatisfactory. Thus, it is imperative to further understand the potential carcinogenesis of PaCa and reliable prognostic models.The gene expression profile and clinical information of GSE21501 were downloaded from the Gene Expression Omnibus (GEO) database. Weighted gene co-expression network analysis (WGCNA) was applied to explore the potent genes associated with the overall survival (OS) events of PaCa patients. Cox regression model was applied to selecting prognostic genes and establish prognostic model. The prognostic values of six-gene signature were validated in TCGA-PAAD cohort.According to the WGCNA analysis, a total of 19 modules were identified and 115 hub genes in the mostly associated module were reserved for next analysis. According to the univariate and multivariate Cox regression analysis, we established a six-gene signature (FTSJ3, STAT1, STX2, CDX2, RASSF4, MACF1) which could effectively evaluate the overall survival (OS) of PaCa patients. In validated patients' cohorts, the six-gene signature exhibited excellent prognostic value in TCGA-PAAD cohort as well.We developed a six-gene signature to exactly predict OS of PaCa patients and provide a novel personalized strategy for evaluating prognosis. The findings may be contributed to medical customization and therapeutic decision in clinical practice.

Elucidation of Biological Networks across Complex Diseases Using Single-Cell Omics

Single-cell multimodal omics (scMulti-omics) technologies have made it possible to trace cellular lineages during differentiation and to identify new cell types in heterogeneous cell populations. The derived information is especially promising for computing cell-type-specific biological networks encoded in complex diseases and improving our understanding of the underlying gene regulatory mechanisms. The integration of these networks could, therefore, give rise to a heterogeneous regulatory landscape (HRL) in support of disease diagnosis and drug therapeutics. In this review, we provide an overview of this field and pay particular attention to how diverse biological networks can be inferred in a specific cell type based on integrative methods. Then, we discuss how HRL can advance our understanding of regulatory mechanisms underlying complex diseases and aid in the prediction of prognosis and therapeutic responses. Finally, we outline challenges and future trends that will be central to bringing the field of HRL in complex diseases forward.

Identification of four hub genes as promising biomarkers to evaluate the prognosis of ovarian cancer in silico

BACKGROUND

Ovarian cancer (OvCa) is one of the most fatal cancers among females in the world. With growing numbers of individuals diagnosed with OvCa ending in deaths, it is urgent to further explore the potential mechanisms of OvCa oncogenesis and development and related biomarkers.

METHODS

The gene expression profiles of GSE49997 were downloaded from the Gene Expression Omnibus (GEO) database. Weighted gene co-expression network analysis (WGCNA) was applied to explore the most potent gene modules associated with the overall survival (OS) and progression-free survival (PFS) events of OvCa patients, and the prognostic values of these genes were exhibited and validated based on data from training and validation sets. Next, protein-protein interaction (PPI) networks were built by GeneMANIA. Besides, enrichment analysis was conducted using DAVID website.

RESULTS

According to the WGCNA analysis, a total of eight modules were identified and four hub genes (MM > 0.90) in the blue module were reserved for next analysis. Kaplan-Meier analysis exhibited that these four hub genes were significantly associated with worse OS and PFS in the patient cohort from GSE49997. Moreover, we validated the short-term (4-years) and long-term prognostic values based on the GSE9891 data, respectively. Last, PPI networks analysis, Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis revealed several potential mechanisms of four hub genes and their co-operators participating in OvCa progression.

CONCLUSION

Four hub genes (COL6A3, CRISPLD2, FBN1 and SERPINF1) were identified to be associated with the prognosis in OvCa, which might be used as monitoring biomarkers to evaluate survival time of OvCa patients.