Network-based differential gene expression analysis suggests cell cycle related genes regulated by E2F1 underlie the molecular difference between smoker and non-smoker lung adenocarcinoma

The MPRINT Hub Data, Model, Knowledge and Research Coordination Center: Bridging the gap in maternal-pediatric therapeutics research through data integration and pharmacometrics

Maternal and pediatric populations have historically been considered "therapeutic orphans" due to their limited inclusion in clinical trials. Physiologic changes during pregnancy and lactation and growth and maturation of children alter pharmacokinetics (PK) and pharmacodynamics (PD) of drugs. Precision therapy in these populations requires knowledge of these effects. Efforts to enhance maternal and pediatric participation in clinical studies have increased over the past few decades. However, studies supporting precision therapeutics in these populations are often small and, in isolation, may have limited impact. Integration of data from various studies, for example through physiologically based pharmacokinetic/pharmacodynamic (PBPK/PD) modeling or bioinformatics approaches, can augment the value of data from these studies, and help identify gaps in understanding. To catalyze research in maternal and pediatric precision therapeutics, the Obstetric and Pediatric Pharmacology and Therapeutics Branch of the Eunice Kennedy Shriver National Institute of Child Health and Human Development (NICHD) established the Maternal and Pediatric Precision in Therapeutics (MPRINT) Hub. Herein, we provide an overview of the status of maternal-pediatric therapeutics research and introduce the Indiana University-Ohio State University MPRINT Hub Data, Model, Knowledge and Research Coordination Center (DMKRCC), which aims to facilitate research in maternal and pediatric precision therapeutics through the integration and assessment of existing knowledge, supporting pharmacometrics and clinical trials design, development of new real-world evidence resources, educational initiatives, and building collaborations among public and private partners, including other NICHD-funded networks. By fostering use of existing data and resources, the DMKRCC will identify critical gaps in knowledge and support efforts to overcome these gaps to enhance maternal-pediatric precision therapeutics.

WDNE: an integrative graphical model for inferring differential networks from multi-platform gene expression data with missing values

The mechanisms controlling biological process, such as the development of disease or cell differentiation, can be investigated by examining changes in the networks of gene dependencies between states in the process. High-throughput experimental methods, like microarray and RNA sequencing, have been widely used to gather gene expression data, which paves the way to infer gene dependencies based on computational methods. However, most differential network analysis methods are designed to deal with fully observed data, but missing values, such as the dropout events in single-cell RNA-sequencing data, are frequent. New methods are needed to take account of these missing values. Moreover, since the changes of gene dependencies may be driven by certain perturbed genes, considering the changes in gene expression levels may promote the identification of gene network rewiring. In this study, a novel weighted differential network estimation (WDNE) model is proposed to handle multi-platform gene expression data with missing values and take account of changes in gene expression levels. Simulation studies demonstrate that WDNE outperforms state-of-the-art differential network estimation methods. When applied WDNE to infer differential gene networks associated with drug resistance in ovarian tumors, cell differentiation and breast tumor heterogeneity, the hub genes in the estimated differential gene networks can provide important insights into the underlying mechanisms. Furthermore, a Matlab toolbox, differential network analysis toolbox, was developed to implement the WDNE model and visualize the estimated differential networks.

Identification of Key Modules and Hub Genes of Annulus Fibrosus in Intervertebral Disc Degeneration

Background: Intervertebral disc degeneration impairs the quality of patients lives. Even though there has been development of many therapeutic strategies, most of them remain unsatisfactory due to the limited understanding of the mechanisms that underlie the intervertebral disc degeneration. Questions/purposes: This study is meant to identify the key modules and hub genes related to the annulus fibrosus in intervertebral disc degeneration (IDD) through: (1) constructing a weighted gene co-expression network; (2) identifying key modules and hub genes; (3) verifying the relationships of key modules and hub genes with IDD; and (4) confirming the expression pattern of hub genes in clinical samples. Methods: The Gene Expression Omnibus provided 24 sets of annulus fibrosus microarray data. Differentially expressed genes between the annulus fibrosus of degenerative and non-degenerative intervertebral disc samples have gone through the Gene Ontology (GO) and pathway analysis. The construction of a gene network and classification of genes into different modules were conducted through performing Weighted Gene Co-expression Network Analysis. The identification of modules and hub genes that were most related to intervertebral disc degeneration was proceeded. In order to verify the relationships of the module and hub genes with intervertebral disc degeneration, Ingenuity Pathway Analysis was operated. Clinical samples were adopted to help verify the hub gene expression profile. Results: One thousand one hundred ninety differentially expressed genes were identified. Terms and pathways associated with intervertebral disc degeneration were presented by GO and pathway analysis. The construction of a Weighted Gene Coexpression Network was completed and clustering differentially expressed genes into four modules was also achieved. The module with the lowest P-value and the highest absolute correlation coefficient was selected and its relationship with intervertebral disc degeneration was confirmed by Ingenuity Pathway Analysis. The identification of hub genes and the confirmation of their expression profile were also realized. Conclusions: This study generated a comprehensive overview of the gene networks underlying annulus fibrosus in intervertebral disc degeneration. Clinical Relevance: Modules and hub genes identified in this study are highly associated with intervertebral disc degeneration, and may serve as potential therapeutic targets for intervertebral disc degeneration.

Assessing the Role of Long Noncoding RNA in Nucleus Accumbens in Subjects With Alcohol Dependence

BACKGROUND

Long noncoding RNA (lncRNA) have been implicated in the etiology of alcohol use. Since lncRNA provide another layer of complexity to the transcriptome, assessing their expression in the brain is the first critical step toward understanding lncRNA functions in alcohol use and addiction. Thus, we sought to profile lncRNA expression in the nucleus accumbens (NAc) in a large postmortem alcohol brain sample.

METHODS

LncRNA and protein-coding gene (PCG) expressions in the NAc from 41 subjects with alcohol dependence (AD) and 41 controls were assessed via a regression model. Weighted gene coexpression network analysis was used to identify lncRNA and PCG networks (i.e., modules) significantly correlated with AD. Within the significant modules, key network genes (i.e., hubs) were also identified. The lncRNA and PCG hubs were correlated via Pearson correlations to elucidate the potential biological functions of lncRNA. The lncRNA and PCG hubs were further integrated with GWAS data to identify expression quantitative trait loci (eQTL).

RESULTS

At Bonferroni adj. p-value ≤ 0.05, we identified 19 lncRNA and 5 PCG significant modules, which were enriched for neuronal and immune-related processes. In these modules, we further identified 86 and 315 PCG and lncRNA hubs, respectively. At false discovery rate (FDR) of 10%, the correlation analyses between the lncRNA and PCG hubs revealed 3,125 positive and 1,860 negative correlations. Integration of hubs with genotype data identified 243 eQTLs affecting the expression of 39 and 204 PCG and lncRNA hubs, respectively.

CONCLUSIONS

Our study identified lncRNA and gene networks significantly associated with AD in the NAc, coordinated lncRNA and mRNA coexpression changes, highlighting potentially regulatory functions for the lncRNA, and our genetic (cis-eQTL) analysis provides novel insights into the etiological mechanisms of AD.

Identifying drug candidates for hepatocellular carcinoma based on differentially expressed genes

The prognosis for patients with advanced hepatocellular carcinoma (HCC) is extremely poor, mainly due to rapid progression and a paucity of effective drugs. Genome-wide analysis allows for potential drugs to be explored based on differentially expressed genes (DEGs). However, drug candidates and DEGs in HCC are largely unknown. In this study, we investigated DEGs and prognostication using The Cancer Genome Atlas (TCGA), the International Cancer Genome Consortium (ICGC), the Gene Expression Omnibus (GEO), and immunohistochemical staining. Protein-protein interaction networks between DEGs were also analyzed to clarify 12 hub genes and query online databases for potential HCC therapeutic drugs. We found that 885 of 3219 DEGs from a TCGA dataset were associated with prognosis. We clarified 12 hub genes that were overexpressed in tumor samples and significantly associated with poor overall survival (OS) in HCC patients. These findings were validated using GEO and ICGC cohorts. Moreover, promising drug candidates targeted against HCC were predicted using online databases. Collectively, the upregulation of 12 hub genes was associated with poor prognosis for patients with HCC, and focusing on their expression may advance efforts towards targeted HCC therapies.

Gene expression network analysis of lymph node involvement in colon cancer identifies AHSA2, CDK10, and CWC22 as possible prognostic markers

Colon cancer has been well studied using a variety of molecular techniques, including whole genome sequencing. However, genetic markers that could be used to predict lymph node (LN) involvement, which is the most important prognostic factor for colon cancer, have not been identified. In the present study, we compared LN(+) and LN(−) colon cancer patients using differential gene expression and network analysis. Colon cancer gene expression data were obtained from the Cancer Genome Atlas and divided into two groups, LN(+) and LN(−). Gene expression networks were constructed using LASSO (Least Absolute Shrinkage and Selection Operator) regression. We identified hub genes, such as APBB1, AHSA2, ZNF767, and JAK2, that were highly differentially expressed. Survival analysis using selected hub genes, such as AHSA2, CDK10, and CWC22, showed that their expression levels were significantly associated with the survival rate of colon cancer patients, which indicates their possible use as prognostic markers. In addition, protein-protein interaction network, GO enrichment, and KEGG pathway analysis were performed with selected hub genes from each group to investigate the regulatory relationships between hub genes and LN involvement in colon cancer; these analyses revealed differences between the LN(−) and LN(+) groups. Our network analysis may help narrow down the search for novel candidate genes for the treatment of colon cancer, in addition to improving our understanding of the biological processes underlying LN involvement. All R implementation codes are available at journal website as Supplementary Materials.

GeneSurrounder: network-based identification of disease genes in expression data

Background

A key challenge of identifying disease–associated genes is analyzing transcriptomic data in the context of regulatory networks that control cellular processes in order to capture multi-gene interactions and yield mechanistically interpretable results. One existing category of analysis techniques identifies groups of related genes using interaction networks, but these gene sets often comprise tens or hundreds of genes, making experimental follow-up challenging. A more recent category of methods identifies precise gene targets while incorporating systems-level information, but these techniques do not determine whether a gene is a driving source of changes in its network, an important characteristic when looking for potential drug targets.

Results

We introduce GeneSurrounder, an analysis method that integrates expression data and network information in a novel procedure to detect genes that are sources of dysregulation on the network. The key idea of our method is to score genes based on the evidence that they influence the dysregulation of their neighbors on the network in a manner that impacts cell function. Applying GeneSurrounder to real expression data, we show that our method is able to identify biologically relevant genes, integrate pathway and expression data, and yield more reproducible results across multiple studies of the same phenotype than competing methods.

Conclusions

Together these findings suggest that GeneSurrounder provides a new avenue for identifying individual genes that can be targeted therapeutically. The key innovation of GeneSurrounder is the combination of pathway network information with gene expression data to determine the degree to which a gene is a source of dysregulation on the network. By prioritizing genes in this way, our method provides insights into disease mechanisms and suggests diagnostic and therapeutic targets. Our method can be used to help biologists select among tens or hundreds of genes for further validation. The implementation in R is available at github.com/sahildshah1/gene-surrounder.

Electronic supplementary material

The online version of this article (10.1186/s12859-019-2829-y) contains supplementary material, which is available to authorized users.

A network-based signature to predict the survival of non-smoking lung adenocarcinoma

Background

A substantial increase in the number of non-smoking lung adenocarcinoma (LAC) patients has been drawing extensive attention in the past decade. However, effective biomarkers, which could guide the precise treatment, are still limited for identifying high-risk patients. Here, we provide a network-based signature to predict the survival of non-smoking LAC.

Materials and methods

Gene expression profiles were downloaded from The Cancer Genome Atlas and Gene Expression Omnibus. Significant gene co-expression networks and hub genes were identified by Weighted Gene Co-expression Network Analysis. Potential mechanisms and pathways of co-expression networks were analyzed by Gene Ontology. The predictive signature was constructed by penalized Cox regression analysis and tested in two independent datasets.

Results

Two distinct co-expression modules were significantly correlated with the non-smoking status across 4 Gene Expression Omnibus datasets. Gene Ontology revealed that nuclear division and cell cycle pathways were main mechanisms of the blue module and that genes in the turquoise module were involved in lymphocyte activation and cell adhesion pathways. Seventeen genes were selected from hub genes at an optimal lambda value and built the prognostic signature. The prognostic signature distinguished the survival of non-smoking LAC (training: hazard ratio [HR]=3.696, 95% CI: 2.025–6.748, P<0.001; testing: HR=2.9, 95% CI: 1.322–6.789, P=0.006; HR=2.78, 95% CI: 1.658–6.654, P=0.022) and had moderate predictive abilities in the training and validation datasets.

Conclusion

The prognostic signature is a promising predictor of non-smoking LAC patients, which might benefit clinical practice and precision therapeutic management.

Meta-analysis of mRNA expression profiles to identify differentially expressed genes in lung adenocarcinoma tissue from smokers and non-smokers

Compared to other types of lung cancer, lung adenocarcinoma patients with a history of smoking have a poor prognosis during the treatment of lung cancer. How lung adenocarcinoma-related genes are differentially expressed between smoker and non-smoker patients has yet to be fully elucidated. We performed a meta-analysis of four publicly available microarray datasets related to lung adenocarcinoma tissue in patients with a history of smoking using R statistical software. The top 50 differentially expressed genes (DEGs) in smoking vs. non‑smoking patients are shown using heat maps. Additionally, we conducted KEGG and GO analyses. In addition, we performed a PPI network analysis for 8 genes that were selected during a previous analysis. We identified a total of 2,932 DEGs (1,806 upregulated, 1,126 downregulated) and five genes (CDC45, CDC20, ANAPC7, CDC6, ESPL1) that may link lung adenocarcinoma to smoking history. Our study may provide new insights into the complex mechanisms of lung adenocarcinoma in smoking patients, and our novel gene expression signatures will be useful for future clinical studies.

Screening of potential therapy targets for prostate cancer using integrated analysis of two gene expression profiles

The aim of the present study was to analyze potential therapy targets for prostate cancer using integrated analysis of two gene expression profiles. First, gene expression profiles GSE38241 and GSE3933 were downloaded from the Gene Expression Omnibus database. Differentially expressed genes (DEGs) between prostate cancer and normal control samples were identified using the Linear Models for Microarray Data package. Pathway enrichment analysis of DEGs was performed using Gene Ontology and the Kyoto Encyclopedia of Genes and Genomes. Furthermore, protein-protein interaction (PPI) networks of DEGs were constructed, on the basis of the Search Tool for the Retrieval of Interacting Genes/Proteins database. The Molecular Complex Detection was utilized to perform module analysis of the PPI networks. In addition, transcriptional regulatory networks were constructed on the basis of the associations between transcription factors (TFs) and target genes. A total of 529 DEGs were identified, including 129 upregulated genes that were primarily associated with to the cell cycle. Additionally, 400 downregulated genes were identified, which were principally enriched in the pathways associated with vascular smooth muscle contraction and focal adhesion. Cell Division Cycle Associated 8, Cell Division Cycle 45, Ubiquitin Conjugating Enzyme E2 C and Thymidine Kinase 1 were identified as hub genes in the upregulated sub-network. Furthermore, the upregulated TF E2F, and the downregulated TF Early Growth Response 1, were identified to be critical in the transcriptional regulatory networks. The identified DEGs and TFs may have critical roles in the progression of prostate cancer, and may be used as target molecules for treating prostate cancer.

Identification of Key Modules and Hub Genes of Keloids with Weighted Gene Coexpression Network Analysis

BACKGROUND

Keloid scarring impairs patients' quality of life, and although many therapeutic strategies have been developed, most remain unsatisfactory because of limited understanding of the mechanisms underlying keloid development.

METHODS

A microarray gene expression data set from keloid tissue was acquired from the Gene Expression Omnibus. Differentially expressed genes in fibroblasts and keratinocytes underwent functional annotation and pathway analysis. Weighted gene coexpression network analysis was applied to identify the gene targets of keloid scars within differentially expressed genes. Modules and hub genes for keloids were identified. Enrichment analysis was undertaken to verify the modules' and hub genes' relationship with keloids.

RESULTS

Enrichment analysis and pathway analysis showed gene ontology terms and pathways related to keloids. Each cell type generated three modules in weighted gene coexpression network analysis, with one module most related to keloids. Enrichment analysis showed that the modules concerned are enriched with terms related to keloids. Three hub genes were selected for fibroblasts and keratinocytes, and their relationship to keloids was verified. Immunohistochemical staining verified expression change of some hub genes.

CONCLUSIONS

This is the first study to describe the gene networks underlying keloids. Modules and hub genes generated in the present study are highly related to keloids and may identify novel therapeutic targets for treatment of keloids.

CLINICAL QUESTION/LEVEL OF EVIDENCE

Therapeutic, V.

Construction of protein interaction network involved in lung adenocarcinomas using a novel algorithm

Studies that only assess differentially-expressed (DE) genes do not contain the information required to investigate the mechanisms of diseases. A complete knowledge of all the direct and indirect interactions between proteins may act as a significant benchmark in the process of forming a comprehensive description of cellular mechanisms and functions. The results of protein interaction network studies are often inconsistent and are based on various methods. In the present study, a combined network was constructed using selected gene pairs, following the conversion and combination of the scores of gene pairs that were obtained across multiple approaches by a novel algorithm. Samples from patients with and without lung adenocarcinoma were compared, and the RankProd package was used to identify DE genes. The empirical Bayesian (EB) meta-analysis approach, the search tool for the retrieval of interacting genes/proteins database (STRING), the weighted gene coexpression network analysis (WGCNA) package and the differentially-coexpressed genes and links package (DCGL) were used for network construction. A combined network was also constructed with a novel rank-based algorithm using a combined score. The topological features of the 5 networks were analyzed and compared. A total of 941 DE genes were screened. The topological analysis indicated that the gene interaction network constructed using the WGCNA method was more likely to produce a small-world property, which has a small average shortest path length and a large clustering coefficient, whereas the combined network was confirmed to be a scale-free network. Gene pairs that were identified using the novel combined method were mostly enriched in the cell cycle and p53 signaling pathway. The present study provided a novel perspective to the network-based analysis. Each method has advantages and disadvantages. Compared with single methods, the combined algorithm used in the present study may provide a novel method to analyze gene interactions, with increased credibility.

Comparative study of computational methods for reconstructing genetic networks of cancer-related pathways

Network reconstruction is an important yet challenging task in systems biology. While many methods have been recently proposed for reconstructing biological networks from diverse data types, properties of estimated networks and differences between reconstruction methods are not well understood. In this paper, we conduct a comprehensive empirical evaluation of seven existing network reconstruction methods, by comparing the estimated networks with different sparsity levels for both normal and tumor samples. The results suggest substantial heterogeneity in networks reconstructed using different reconstruction methods. Our findings also provide evidence for significant differences between networks of normal and tumor samples, even after accounting for the considerable variability in structures of networks estimated using different reconstruction methods. These differences can offer new insight into changes in mechanisms of genetic interaction associated with cancer initiation and progression.