Weighted gene coexpression network analysis strategies applied to mouse weight

Gene co-expression networks associated with carcass traits reveal new pathways for muscle and fat deposition in Nelore cattle

BACKGROUND

Positively correlated with carcass weight and animal growth, the ribeye area (REA) and the backfat thickness (BFT) are economic important carcass traits, which impact directly on producer's payment. The selection of these traits has not been satisfactory since they are expressed later in the animal's life and multigene regulated. So, next-generation technologies have been applied in this area to improve animal's selection and better understand the molecular mechanisms involved in the development of these traits. Correlation network analysis, performed by tools like WGCNA (Weighted Correlation Network Analysis), has been used to explore gene-gene interactions and gene-phenotype correlations. Thus, this study aimed to identify putative candidate genes and metabolic pathways that regulate REA and BFT by constructing a gene co-expression network using WGCNA and RNA sequencing data, to better understand genetic and molecular variations behind these complex traits in Nelore cattle.

RESULTS

The gene co-expression network analysis, using WGCNA, were built using RNA-sequencing data normalized by transcript per million (TPM) from 43 Nelore steers. Forty-six gene clusters were constructed, between them, three were positively correlated (p-value< 0.1) to the BFT (Green Yellow, Ivory, and Light Yellow modules) and, one cluster was negatively correlated (p-value< 0.1) with REA (Salmon module). The enrichment analysis performed by DAVID and WebGestalt (FDR 5%) identified eight Gene Ontology (GO) terms and three KEGG pathways in the Green Yellow module, mostly associated with immune response and inflammatory mechanisms. The enrichment of the Salmon module demonstrated 19 GO terms and 21 KEGG pathways, related to muscle energy metabolism, lipid metabolism, muscle degradation, and oxidative stress diseases. The Ivory and Light yellow modules have not shown significant results in the enrichment analysis.

CONCLUSION

With this study, we verified that inflammation and immune response pathways modulate the BFT trait. Energy and lipid metabolism pathways, highlighting fatty acid metabolism, were the central pathways associated with REA. Some genes, as RSAD2, EIF2AK2, ACAT1, and ACSL1 were considered as putative candidate related to these traits. Altogether these results allow us to a better comprehension of the molecular mechanisms that lead to muscle and fat deposition in bovine.

A Scalable Permutation Approach Reveals Replication and Preservation Patterns of Network Modules in Large Datasets

Network modules-topologically distinct groups of edges and nodes-that are preserved across datasets can reveal common features of organisms, tissues, cell types, and molecules. Many statistics to identify such modules have been developed, but testing their significance requires heuristics. Here, we demonstrate that current methods for assessing module preservation are systematically biased and produce skewed p values. We introduce NetRep, a rapid and computationally efficient method that uses a permutation approach to score module preservation without assuming data are normally distributed. NetRep produces unbiased p values and can distinguish between true and false positives during multiple hypothesis testing. We use NetRep to quantify preservation of gene coexpression modules across murine brain, liver, adipose, and muscle tissues. Complex patterns of multi-tissue preservation were revealed, including a liver-derived housekeeping module that displayed adipose- and muscle-specific association with body weight. Finally, we demonstrate the broader applicability of NetRep by quantifying preservation of bacterial networks in gut microbiota between men and women.

Identification of glioblastoma gene prognosis modules based on weighted gene co-expression network analysis

Background

Glioblastoma multiforme, the most prevalent and aggressive brain tumour, has a poor prognosis. The molecular mechanisms underlying gliomagenesis remain poorly understood. Therefore, molecular research, including various markers, is necessary to understand the occurrence and development of glioma.

Method

Weighted gene co-expression network analysis (WGCNA) was performed to construct a gene co-expression network in TCGA glioblastoma samples. Gene ontology (GO) and pathway-enrichment analysis were used to identify significance of gene modules. Cox proportional hazards regression model was used to predict outcome of glioblastoma patients.

Results

We performed weighted gene co-expression network analysis (WGCNA) and identified a gene module (yellow module) related to the survival time of TCGA glioblastoma samples. Then, 228 hub genes were calculated based on gene significance (GS) and module significance (MS). Four genes (OSMR + SOX21 + MED10 + PTPRN) were selected to construct a Cox proportional hazards regression model with high accuracy (AUC = 0.905). The prognostic value of the Cox proportional hazards regression model was also confirmed in GSE16011 dataset (GBM: n = 156).

Conclusion

We developed a promising mRNA signature for estimating overall survival in glioblastoma patients.

Electronic supplementary material

The online version of this article (10.1186/s12920-018-0407-1) contains supplementary material, which is available to authorized users.

Comparison of gene co-networks analysis provide a systems view of rice (Oryza sativa L.) response to Tilletia horrida infection

The biotrophic soil-borne fungus Tilletia horrida causes rice kernel smut, an important disease affecting the production of rice male sterile lines in most hybrid rice growing regions of the world. There are no successful ways of controlling this disease and there has been little study of mechanisms of resistance to T. horrida. Based on transcriptional data of different infection time points, we found 23, 782 and 23, 718 differentially expressed genes (fragments per kilobase of transcript sequence per million, FPKM >1) in Jiangcheng 3A (resistant to T. horrida) and 9311A (susceptible to T. horrida), respectively. In order to illuminate the differential responses of the two rice male sterile lines to T. horrida, we identified gene co-expression modules using the method of weighted gene co-expression network analysis (WGCNA) and compared the different biological functions of gene co-expression networks in key modules at different infection time points. The results indicated that gene co-expression networks in the two rice genotypes were different and that genes contained in some modules of the two groups may play important roles in resistance to T. horrida, such as DTH8 and OsHop/Sti1a. Furthermore, these results provide a global view of the responses of two different phenotypes to T. horrida, and assist our understanding of the regulation of expression changes after T. horrida infection.

Gene Expression Signatures Can Aid Diagnosis of Sexually Transmitted Infection-Induced Endometritis in Women

Sexually transmitted infection (STI) of the upper reproductive tract can result in inflammation and infertility. A biomarker of STI-induced upper tract inflammation would be significant as many women are asymptomatic and delayed treatment increases risk of sequelae. Blood mRNA from 111 women from three cohorts was profiled using microarray. Unsupervised analysis revealed a transcriptional profile that distinguished 9 cases of STI-induced endometritis from 18 with cervical STI or uninfected controls. Using a hybrid feature selection algorithm we identified 21 genes that yielded maximal classification accuracy within our training dataset. Predictive accuracy was evaluated using an independent testing dataset of 5 cases and 10 controls. Sensitivity was evaluated in a separate test set of 12 women with asymptomatic STI-induced endometritis in whom cervical burden was determined by PCR; and specificity in an additional test set of 15 uninfected women with pelvic pain due to unknown cause. Disease module preservation was assessed in 42 women with a clinical diagnosis of pelvic inflammatory disease (PID). We also tested the ability of the biomarker to discriminate STI-induced endometritis from other diseases. The biomarker was 86.7% (13/15) accurate in correctly distinguishing cases from controls in the testing dataset. Sensitivity was 83.3% (5/6) in women with high cervical Chlamydia trachomatis burden and asymptomatic endometritis, but 0% (0/6) in women with low burden. Specificity in patients with non-STI-induced pelvic pain was 86.7% (13/15). Disease modules were preserved in all 8 biomarker predicted cases. The 21-gene biomarker was highly discriminatory for systemic infections, lupus, and appendicitis, but wrongly predicted tuberculosis as STI-induced endometritis in 52.4%. A 21-gene biomarker can identify asymptomatic women with STI-induced endometritis that places them at risk for chronic disease development and discriminate STI-induced endometritis from non-STI pelvic pain and other diseases.

WISH-R- a fast and efficient tool for construction of epistatic networks for complex traits and diseases

Background

Genetic epistasis is an often-overlooked area in the study of the genomics of complex traits. Genome-wide association studies are a useful tool for revealing potential causal genetic variants, but in this context, epistasis is generally ignored. Data complexity and interpretation issues make it difficult to process and interpret epistasis. As the number of interaction grows exponentially with the number of variants, computational limitation is a bottleneck. Gene Network based strategies have been successful in integrating biological data and identifying relevant hub genes and pathways related to complex traits. In this study, epistatic interactions and network-based analysis are combined in the Weighted Interaction SNP hub (WISH) method and implemented in an efficient and easy to use R package.

Results

The WISH R package (WISH-R) was developed to calculate epistatic interactions on a genome-wide level based on genomic data. It is easy to use and install, and works on regular genomic data. The package filters data based on linkage disequilibrium and calculates epistatic interaction coefficients between SNP pairs based on a parallelized efficient linear model and generalized linear model implementations. Normalized epistatic coefficients are analyzed in a network framework, alleviating multiple testing issues and integrating biological signal to identify modules and pathways related to complex traits. Functions for visualizing results and testing runtimes are also provided.

Conclusion

The WISH-R package is an efficient implementation for analyzing genome-wide epistasis for complex diseases and traits. It includes methods and strategies for analyzing epistasis from initial data filtering until final data interpretation. WISH offers a new way to analyze genomic data by combining epistasis and network based analysis in one method and provides options for visualizations. This alleviates many of the existing hurdles in the analysis of genomic interactions.

Electronic supplementary material

The online version of this article (10.1186/s12859-018-2291-2) contains supplementary material, which is available to authorized users.

Integrated pathogen load and dual transcriptome analysis of systemic host-pathogen interactions in severe malaria

The pathogenesis of infectious diseases depends on the interaction of host and pathogen. In Plasmodium falciparum malaria, host and parasite processes can be assessed by dual RNA sequencing of blood from infected patients. We performed dual transcriptome analyses on samples from 46 malaria-infected Gambian children to reveal mechanisms driving the systemic pathophysiology of severe malaria. Integrating these transcriptomic data with estimates of parasite load and detailed clinical information allowed consideration of potentially confounding effects due to differing leukocyte proportions in blood, parasite developmental stage, and whole-body pathogen load. We report hundreds of human and parasite genes differentially expressed between severe and uncomplicated malaria, with distinct profiles associated with coma, hyperlactatemia, and thrombocytopenia. High expression of neutrophil granule-related genes was consistently associated with all severe malaria phenotypes. We observed severity-associated variation in the expression of parasite genes, which determine cytoadhesion to vascular endothelium, rigidity of infected erythrocytes, and parasite growth rate. Up to 99% of human differential gene expression in severe malaria was driven by differences in parasite load, whereas parasite gene expression showed little association with parasite load. Coexpression analyses revealed interactions between human and P. falciparum, with prominent co-regulation of translation genes in severe malaria between host and parasite. Multivariate analyses suggested that increased expression of granulopoiesis and interferon-γ-related genes, together with inadequate suppression of type 1 interferon signaling, best explained severity of infection. These findings provide a framework for understanding the contributions of host and parasite to the pathogenesis of severe malaria and identifying new treatments.

Network-Based Predictors of Progression in Head and Neck Squamous Cell Carcinoma

The heterogeneity in head and neck squamous cell carcinoma (HNSCC) has made reliable stratification extremely challenging. Behavioral risk factors such as smoking and alcohol consumption contribute to this heterogeneity. To help elucidate potential mechanisms of progression in HNSCC, we focused on elucidating patterns of gene interactions associated with tumor progression. We performed de-novo gene co-expression network inference utilizing 229 patient samples from The Cancer Genome Atlas (TCGA) previously annotated by Bornstein et al. (2016). Differential network analysis allowed us to contrast progressor and non-progressor cohorts. Beyond standard differential expression (DE) analysis, this approach evaluates changes in gene expression variance (differential variability DV) and changes in covariance, which we denote as differential wiring (DW). The set of affected genes was overlaid onto the co-expression network, identifying 12 modules significantly enriched in DE, DV, and/or DW genes. Additionally, we identified modules correlated with behavioral measures such as alcohol consumption and smoking status. In the module enriched for differentially wired genes, we identified network hubs including IL10RA, DOK2, APBB1IP, UBASH3A, SASH3, CELF2, TRAF3IP3, GIMAP6, MYO1F, NCKAP1L, WAS, FERMT3, SLA, SELPLG, TNFRSF1B, WIPF1, AMICA1, PTPN22; the network centrality and progression specificity of these genes suggest a potential role in tumor evolution mechanisms related to inflammation and microenvironment. The identification of this network-based gene signature could be further developed to guide progression stratification, highlighting how network approaches may help improve clinical research end points and ultimately aid in clinical utility.

Effects of threshold on the topology of gene co-expression networks

Several developments regarding the analysis of gene co-expression profiles using complex network theory have been reported recently. Such approaches usually start with the construction of an unweighted gene co-expression network, therefore requiring the selection of a suitable threshold defining which pairs of vertices will be connected. We aimed at addressing such an important problem by suggesting and comparing five different approaches for threshold selection. Each of the methods considers a respective biologically-motivated criterion for electing a potentially suitable threshold. A set of 21 microarray experiments from different biological groups was used to investigate the effect of applying the five proposed criteria to several biological situations. For each experiment, we used the Pearson correlation coefficient to measure the relationship between each gene pair, and the resulting weight matrices were thresholded considering several values, generating respective adjacency matrices (co-expression networks). Each of the five proposed criteria was then applied in order to select the respective threshold value. The effects of these thresholding approaches on the topology of the resulting networks were compared by using several measurements, and we verified that, depending on the database, the impact on the topological properties can be large. However, a group of databases was verified to be similarly affected by most of the considered criteria. Based on such results, it can be suggested that when the generated networks present similar measurements, the thresholding method can be chosen with greater freedom. If the generated networks are markedly different, the thresholding method that better suits the interests of each specific research study represents a reasonable choice.

Transcription Factor Co-expression Networks of Adipose RNA-Seq Data Reveal Regulatory Mechanisms of Obesity

BACKGROUND

Transcription Factors (TFs) control actuation of genes in the genome and are key mediators of complex processes such as obesity. Master Regulators (MRs) are the genes at the top of a regulation hierarchy which regulate other genes.

OBJECTIVE

To elucidate clusters of highly co-expressed TFs (modules), involved pathways, highly inter-connected TFs (hub-TFs) and MRs leading to obesity and leanness, using porcine model for human obesity.

METHODS

We identified 817 expressed TFs in RNA-Sequencing dataset representing extreme degrees of obesity (DO; lean, obese). We built a single Weighted Transcription Factor Co-expression Network (WTFCN) and TF sub-networks (based on the DO). Hub-TFs and MRs (using iRegulon) were identi-fied in biologically relevant WTFCNs modules.

RESULTS

Single WTFCN detected the Red module significantly associated with DO (P < 0.03). This module was enriched for regulation processes in the immune system, e.g.: Immune system process (Padj = 2.50E-06) and metabolic lifestyle disorders, e.g. Circadian rhythm - mammal pathway (Padj = 2.33E-11). Detected MR, hub-TF SPI1 was involved in obesity, immunity and osteoporosis. Within the obese sub-network, the Red module suggested possible associations with immunity, e.g. TGF-beta signaling pathway (Padj = 1.73E-02) and osteoporosis, e.g. Osteoclast differentiation (Padj = 1.94E-02). Within the lean sub-network, the Magenta module displayed associations with type 2 diabetes, obesity and os-teoporosis e.g. Notch signaling pathway (Padj = 2.40E-03), osteoporosis e.g. hub-TF VDR (a prime candidate gene for osteoporosis).

CONCLUSION

Our results provide insights into the regulatory network of TFs and biologically relevant hub TFs in obesity.

Unsupervised, Statistically Based Systems Biology Approach for Unraveling the Genetics of Complex Traits: A Demonstration with Ethanol Metabolism

BACKGROUND

A statistical pipeline was developed and used for determining candidate genes and candidate gene coexpression networks involved in 2 alcohol (i.e., ethanol [EtOH]) metabolism phenotypes, namely alcohol clearance and acetate area under the curve in a recombinant inbred (RI) (HXB/BXH) rat panel. The approach was also used to provide an indication of how EtOH metabolism can impact the normal function of the identified networks.

METHODS

RNA was extracted from alcohol-naïve liver tissue of 30 strains of HXB/BXH RI rats. The reconstructed transcripts were quantitated, and data were used to construct gene coexpression modules and networks. A separate group of rats, comprising the same 30 strains, were injected with EtOH (2 g/kg) for measurement of blood EtOH and acetate levels. These data were used for quantitative trait loci (QTL) analysis of the rate of EtOH disappearance and circulating acetate levels. The analysis pipeline required calculation of the module eigengene values, the correction of these values with EtOH metabolism rates and acetate levels across the rat strains, and the determination of the eigengene QTLs. For a module to be considered a candidate for determining phenotype, the module eigengene values had to have significant correlation with the strain phenotypic values and the module eigengene QTLs had to overlap the phenotypic QTLs.

RESULTS

Of the 658 transcript coexpression modules generated from liver RNA sequencing data, a single module satisfied all criteria for being a candidate for determining the alcohol clearance trait. This module contained 2 alcohol dehydrogenase genes, including the gene whose product was previously shown to be responsible for the majority of alcohol elimination in the rat. This module was also the only module identified as a candidate for influencing circulating acetate levels. This module was also linked to the process of generation and utilization of retinoic acid as related to the autonomous immune response.

CONCLUSIONS

We propose that our analytical pipeline can successfully identify genetic regions and transcripts which predispose a particular phenotype and our analysis provides functional context for coexpression module components.

Genetic Drivers of Pancreatic Islet Function

The majority of gene loci that have been associated with type 2 diabetes play a role in pancreatic islet function. To evaluate the role of islet gene expression in the etiology of diabetes, we sensitized a genetically diverse mouse population with a Western diet high in fat (45% kcal) and sucrose (34%) and carried out genome-wide association mapping of diabetes-related phenotypes. We quantified mRNA abundance in the islets and identified 18,820 expression QTL. We applied mediation analysis to identify candidate causal driver genes at loci that affect the abundance of numerous transcripts. These include two genes previously associated with monogenic diabetes (PDX1 and HNF4A), as well as three genes with nominal association with diabetes-related traits in humans (FAM83E, IL6ST, and SAT2). We grouped transcripts into gene modules and mapped regulatory loci for modules enriched with transcripts specific for α-cells, and another specific for δ-cells. However, no single module enriched for β-cell-specific transcripts, suggesting heterogeneity of gene expression patterns within the β-cell population. A module enriched in transcripts associated with branched-chain amino acid metabolism was the most strongly correlated with physiological traits that reflect insulin resistance. Although the mice in this study were not overtly diabetic, the analysis of pancreatic islet gene expression under dietary-induced stress enabled us to identify correlated variation in groups of genes that are functionally linked to diabetes-associated physiological traits. Our analysis suggests an expected degree of concordance between diabetes-associated loci in the mouse and those found in human populations, and demonstrates how the mouse can provide evidence to support nominal associations found in human genome-wide association mapping.

Discovery of Blood Transcriptional Endotypes in Women with Pelvic Inflammatory Disease

Sexually transmitted infections with Chlamydia trachomatis and/or Neisseria gonorrhoeae and rates of pelvic inflammatory disease (PID) in women continue to rise, with reinfection being common because of poor adaptive immunity. Diagnosis remains imprecise, and pathogenesis data are derived primarily from monoinfection of mice with C. trachomatis or N. gonorrhoeae By comparing blood mRNA responses of women with C. trachomatis- and/or N. gonorrhoeae-induced PID and histologic endometritis with those from women with C. trachomatis and/or N. gonorrhoeae infection limited to their cervix and asymptomatic uninfected women determined via microarray, we discovered important pathogenic mechanisms in PID and response differences that provide a pathway to biomarker discovery. Women with N. gonorrhoeae- and/or C. trachomatis-induced PID exhibit overexpression of myeloid cell genes and suppression of protein synthesis, mitochondrial oxidative phosphorylation, and T cell-specific genes. Coinfected women exhibited the greatest activation of cell death pathways and suppression of responses essential for adaptive immunity. Women solely infected with C. trachomatis expressed elevated levels of type I and type II IFN genes, and enhanced type I IFN-induced chemokines in cervical secretions were associated with ascension of C. trachomatis to the endometrium. Blood microarrays reveal discrete pathobiological endotypes in women with PID that are driven by pathogen invasion of the upper genital tract.

Genetical genomics of growth in a chicken model

BACKGROUND

The genetics underlying body mass and growth are key to understanding a wide range of topics in biology, both evolutionary and developmental. Body mass and growth traits are affected by many genetic variants of small effect. This complicates genetic mapping of growth and body mass. Experimental intercrosses between individuals from divergent populations allows us to map naturally occurring genetic variants for selected traits, such as body mass by linkage mapping. By simultaneously measuring traits and intermediary molecular phenotypes, such as gene expression, one can use integrative genomics to search for potential causative genes.

RESULTS

In this study, we use linkage mapping approach to map growth traits (N = 471) and liver gene expression (N = 130) in an advanced intercross of wild Red Junglefowl and domestic White Leghorn layer chickens. We find 16 loci for growth traits, and 1463 loci for liver gene expression, as measured by microarrays. Of these, the genes TRAK1, OSBPL8, YEATS4, CEP55, and PIP4K2B are identified as strong candidates for growth loci in the chicken. We also show a high degree of sex-specific gene-regulation, with almost every gene expression locus exhibiting sex-interactions. Finally, several trans-regulatory hotspots were found, one of which coincides with a major growth locus.

CONCLUSIONS

These findings not only serve to identify several strong candidates affecting growth, but also show how sex-specificity and local gene-regulation affect growth regulation in the chicken.

Transcriptomic signatures of cellular and humoral immune responses in older adults after seasonal influenza vaccination identified by data-driven clustering

PBMC transcriptomes after influenza vaccination contain valuable information about factors affecting vaccine responses. However, distilling meaningful knowledge out of these complex datasets is often difficult and requires advanced data mining algorithms. We investigated the use of the data-driven Weighted Gene Correlation Network Analysis (WGCNA) gene clustering method to identify vaccine response-related genes in PBMC transcriptomic datasets collected from 138 healthy older adults (ages 50-74) before and after 2010-2011 seasonal trivalent influenza vaccination. WGCNA separated the 14,197 gene dataset into 15 gene clusters based on observed gene expression patterns across subjects. Eight clusters were strongly enriched for genes involved in specific immune cell types and processes, including B cells, T cells, monocytes, platelets, NK cells, cytotoxic T cells, and antiviral signaling. Examination of gene cluster membership identified signatures of cellular and humoral responses to seasonal influenza vaccination, as well as pre-existing cellular immunity. The results of this study illustrate the utility of this publically available analysis methodology and highlight genes previously associated with influenza vaccine responses (e.g., CAMK4, CD19), genes with functions not previously identified in vaccine responses (e.g., SPON2, MATK, CST7), and previously uncharacterized genes (e.g. CORO1C, C8orf83) likely related to influenza vaccine-induced immunity due to their expression patterns.

Application of Weighted Gene Co-expression Network Analysis for Data from Paired Design

Investigating how genes jointly affect complex human diseases is important, yet challenging. The network approach (e.g., weighted gene co-expression network analysis (WGCNA)) is a powerful tool. However, genomic data usually contain substantial batch effects, which could mask true genomic signals. Paired design is a powerful tool that can reduce batch effects. However, it is currently unclear how to appropriately apply WGCNA to genomic data from paired design. In this paper, we modified the current WGCNA pipeline to analyse high-throughput genomic data from paired design. We illustrated the modified WGCNA pipeline by analysing the miRNA dataset provided by Shiah et al. (2014), which contains forty oral squamous cell carcinoma (OSCC) specimens and their matched non-tumourous epithelial counterparts. OSCC is the sixth most common cancer worldwide. The modified WGCNA pipeline identified two sets of novel miRNAs associated with OSCC, in addition to the existing miRNAs reported by Shiah et al. (2014). Thus, this work will be of great interest to readers of various scientific disciplines, in particular, genetic and genomic scientists as well as medical scientists working on cancer.

Integrated omics dissection of proteome dynamics during cardiac remodeling

Transcript abundance and protein abundance show modest correlation in many biological models, but how this impacts disease signature discovery in omics experiments is rarely explored. Here we report an integrated omics approach, incorporating measurements of transcript abundance, protein abundance, and protein turnover to map the landscape of proteome remodeling in a mouse model of pathological cardiac hypertrophy. Analyzing the hypertrophy signatures that are reproducibly discovered from each omics data type across six genetic strains of mice, we find that the integration of transcript abundance, protein abundance, and protein turnover data leads to 75% gain in discovered disease gene candidates. Moreover, the inclusion of protein turnover measurements allows discovery of post-transcriptional regulations across diverse pathways, and implicates distinct disease proteins not found in steady-state transcript and protein abundance data. Our results suggest that multi-omics investigations of proteome dynamics provide important insights into disease pathogenesis in vivo.

Gene Coexpression Networks Drive and Predict Reproductive Effects in Daphnia in Response to Environmental Disturbances

Increasing effects of anthropogenic stressors and those of natural origin on aquatic ecosystems have intensified the need for predictive and functional models of their effects. Here, we use gene expression patterns in combination with weighted gene coexpression networks and generalized additive models to predict effects on reproduction in the aquatic microcrustacean Daphnia. We developed models to predict effects on reproduction upon exposure to different cyanobacteria, different insecticides and binary mixtures of cyanobacteria and insecticides. Models developed specifically for groups of stressors (e.g., either cyanobacteria or insecticides) performed better than general models developed on all data. Furthermore, models developed using in silico generated mixture gene expression profiles from single stressor data were able to better predict effects on reproduction compared to models derived from the mixture exposures themselves. Our results highlight the potential of gene expression data to quantify effects of complex exposures at higher level organismal effects without prior mechanistic knowledge or complex exposure data.

MultiDCoX: Multi-factor analysis of differential co-expression

Background

Differential co-expression (DCX) signifies change in degree of co-expression of a set of genes among different biological conditions. It has been used to identify differential co-expression networks or interactomes. Many algorithms have been developed for single-factor differential co-expression analysis and applied in a variety of studies. However, in many studies, the samples are characterized by multiple factors such as genetic markers, clinical variables and treatments. No algorithm or methodology is available for multi-factor analysis of differential co-expression.

Results

We developed a novel formulation and a computationally efficient greedy search algorithm called MultiDCoX to perform multi-factor differential co-expression analysis. Simulated data analysis demonstrates that the algorithm can effectively elicit differentially co-expressed (DCX) gene sets and quantify the influence of each factor on co-expression. MultiDCoX analysis of a breast cancer dataset identified interesting biologically meaningful differentially co-expressed (DCX) gene sets along with genetic and clinical factors that influenced the respective differential co-expression.

Conclusions

MultiDCoX is a space and time efficient procedure to identify differentially co-expressed gene sets and successfully identify influence of individual factors on differential co-expression.

Electronic supplementary material

The online version of this article (10.1186/s12859-017-1963-7) contains supplementary material, which is available to authorized users.

Sugar metabolism reprogramming in a non-climacteric bud mutant of a climacteric plum fruit during development on the tree

We investigated sugar metabolism in leaves and fruits of two Japanese plum (Prunus salicina Lindl.) cultivars, the climacteric Santa Rosa and its bud sport mutant the non-climacteric Sweet Miriam, during development on the tree. We previously characterized differences between the two cultivars. Here, we identified key sugar metabolic pathways. Pearson coefficient correlations of metabolomics and transcriptomic data and weighted gene co-expression network analysis (WGCNA) of RNA sequencing (RNA-Seq) data allowed the identification of 11 key sugar metabolism-associated genes: sucrose synthase, sucrose phosphate synthase, cytosolic invertase, vacuolar invertase, invertase inhibitor, α-galactosidase, β-galactosidase, galactokinase, trehalase, galactinol synthase, and raffinose synthase. These pathways were further assessed and validated through the biochemical characterization of the gene products and with metabolite analysis. Our results demonstrated the reprogramming of sugar metabolism in both leaves and fruits in the non-climacteric plum, which displayed a shift towards increased sorbitol synthesis. Climacteric and non-climacteric fruits showed differences in their UDP-galactose metabolism towards the production of galactose and raffinose, respectively. The higher content of galactinol, myo-inositol, raffinose, and trehalose in the non-climacteric fruits could improve the ability of the fruits to cope with the oxidative processes associated with fruit ripening. Overall, our results support a relationship between sugar metabolism, ethylene, and ripening behavior.

Comparative assessment of differential network analysis methods

Differential network analysis (DiNA) denotes a recent class of network-based Bioinformatics algorithms which focus on the differences in network topologies between two states of a cell, such as healthy and disease, to identify key players in the discriminating biological processes. In contrast to conventional differential analysis, DiNA identifies changes in the interplay between molecules, rather than changes in single molecules. This ability is especially important in cases where effectors are changed, e.g. mutated, but their expression is not. A number of different DiNA approaches have been proposed, yet a comparative assessment of their performance in different settings is still lacking. In this paper, we evaluate 10 different DiNA algorithms regarding their ability to recover genetic key players from transcriptome data. We construct high-quality regulatory networks and enrich them with co-expression data from four different types of cancer. Next, we assess the results of applying DiNA algorithms on these data sets using a gold standard list (GSL). We find that local DiNA algorithms are generally superior to global algorithms, and that all DiNA algorithms outperform conventional differential expression analysis. We also assess the ability of DiNA methods to exploit additional knowledge in the underlying cellular networks. To this end, we enrich the cancer-type specific networks with known regulatory miRNAs and compare the algorithms performance in networks with and without miRNA. We find that including miRNAs consistently and considerably improves the performance of almost all tested algorithms. Our results underline the advantages of comprehensive cell models for the analysis of -omics data.

Estimation of the proteomic cancer co-expression sub networks by using association estimators

In this study, the association estimators, which have significant influences on the gene network inference methods and used for determining the molecular interactions, were examined within the co-expression network inference concept. By using the proteomic data from five different cancer types, the hub genes/proteins within the disease-associated gene-gene/protein-protein interaction sub networks were identified. Proteomic data from various cancer types is collected from The Cancer Proteome Atlas (TCPA). Correlation and mutual information (MI) based nine association estimators that are commonly used in the literature, were compared in this study. As the gold standard to measure the association estimators’ performance, a multi-layer data integration platform on gene-disease associations (DisGeNET) and the Molecular Signatures Database (MSigDB) was used. Fisher's exact test was used to evaluate the performance of the association estimators by comparing the created co-expression networks with the disease-associated pathways. It was observed that the MI based estimators provided more successful results than the Pearson and Spearman correlation approaches, which are used in the estimation of biological networks in the weighted correlation network analysis (WGCNA) package. In correlation-based methods, the best average success rate for five cancer types was 60%, while in MI-based methods the average success ratio was 71% for James-Stein Shrinkage (Shrink) and 64% for Schurmann-Grassberger (SG) association estimator, respectively. Moreover, the hub genes and the inferred sub networks are presented for the consideration of researchers and experimentalists.

Weighted gene co-expression network analysis of expression data of monozygotic twins identifies specific modules and hub genes related to BMI

Background

The therapeutic management of obesity is challenging, hence further elucidating the underlying mechanisms of obesity development and identifying new diagnostic biomarkers and therapeutic targets are urgent and necessary. Here, we performed differential gene expression analysis and weighted gene co-expression network analysis (WGCNA) to identify significant genes and specific modules related to BMI based on gene expression profile data of 7 discordant monozygotic twins.

Results

In the differential gene expression analysis, it appeared that 32 differentially expressed genes (DEGs) were with a trend of up-regulation in twins with higher BMI when compared to their siblings. Categories of positive regulation of nitric-oxide synthase biosynthetic process, positive regulation of NF-kappa B import into nucleus, and peroxidase activity were significantly enriched within GO database and NF-kappa B signaling pathway within KEGG database. DEGs of NAMPT, TLR9, PTGS2, HBD, and PCSK1N might be associated with obesity. In the WGCNA, among the total 20 distinct co-expression modules identified, coral1 module (68 genes) had the strongest positive correlation with BMI (r = 0.56, P = 0.04) and disease status (r = 0.56, P = 0.04). Categories of positive regulation of phospholipase activity, high-density lipoprotein particle clearance, chylomicron remnant clearance, reverse cholesterol transport, intermediate-density lipoprotein particle, chylomicron, low-density lipoprotein particle, very-low-density lipoprotein particle, voltage-gated potassium channel complex, cholesterol transporter activity, and neuropeptide hormone activity were significantly enriched within GO database for this module. And alcoholism and cell adhesion molecules pathways were significantly enriched within KEGG database. Several hub genes, such as GAL, ASB9, NPPB, TBX2, IL17C, APOE, ABCG4, and APOC2 were also identified. The module eigengene of saddlebrown module (212 genes) was also significantly correlated with BMI (r = 0.56, P = 0.04), and hub genes of KCNN1 and AQP10 were differentially expressed.

Conclusion

We identified significant genes and specific modules potentially related to BMI based on the gene expression profile data of monozygotic twins. The findings may help further elucidate the underlying mechanisms of obesity development and provide novel insights to research potential gene biomarkers and signaling pathways for obesity treatment. Further analysis and validation of the findings reported here are important and necessary when more sample size is acquired.

Electronic supplementary material

The online version of this article (10.1186/s12864-017-4257-6) contains supplementary material, which is available to authorized users.

Gene Network Construction from Microarray Data Identifies a Key Network Module and Several Candidate Hub Genes in Age-Associated Spatial Learning Impairment

As humans age many suffer from a decrease in normal brain functions including spatial learning impairments. This study aimed to better understand the molecular mechanisms in age-associated spatial learning impairment (ASLI). We used a mathematical modeling approach implemented in Weighted Gene Co-expression Network Analysis (WGCNA) to create and compare gene network models of young (learning unimpaired) and aged (predominantly learning impaired) brains from a set of exploratory datasets in rats in the context of ASLI. The major goal was to overcome some of the limitations previously observed in the traditional meta- and pathway analysis using these data, and identify novel ASLI related genes and their networks based on co-expression relationship of genes. This analysis identified a set of network modules in the young, each of which is highly enriched with genes functioning in broad but distinct GO functional categories or biological pathways. Interestingly, the analysis pointed to a single module that was highly enriched with genes functioning in “learning and memory” related functions and pathways. Subsequent differential network analysis of this “learning and memory” module in the aged (predominantly learning impaired) rats compared to the young learning unimpaired rats allowed us to identify a set of novel ASLI candidate hub genes. Some of these genes show significant repeatability in networks generated from independent young and aged validation datasets. These hub genes are highly co-expressed with other genes in the network, which not only show differential expression but also differential co-expression and differential connectivity across age and learning impairment. The known function of these hub genes indicate that they play key roles in critical pathways, including kinase and phosphatase signaling, in functions related to various ion channels, and in maintaining neuronal integrity relating to synaptic plasticity and memory formation. Taken together, they provide a new insight and generate new hypotheses into the molecular mechanisms responsible for age associated learning impairment, including spatial learning.

Identification of novel therapeutic target genes and pathway in pancreatic cancer by integrative analysis

BACKGROUND

Gene alterations are crucial to the molecular pathogenesis of pancreatic cancer. The present study was designed to identify the potential candidate genes in the pancreatic carcinogenesis.

METHODS

Gene Expression Omnibus database (GEO) datasets of pancreatic cancer tissue were retrieval and the differentially expressed genes (DEGs) from individual microarray data were merged. Gene Ontology (GO) enrichment analysis, Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis, protein-protein interaction (PPI) networks, and gene coexpression analysis were performed.

RESULTS

Three GEO datasets, including 74 pancreatic cancer samples and 55 controls samples were selected. A total of 2325 DEGs were identified, including 1383 upregulated and 942 downregulated genes. The GO terms for molecular functions, biological processes, and cellular component were protein binding, small molecule metabolic process, and integral to membrane, respectively. The most significant pathway in KEGG analysis was metabolic pathways. PPI network analysis indicated that the significant hub genes including cytochrome P450, family 2, subfamily E, polypeptide 1 (CYP2E1), mitogen-activated protein kinase 3 (MAPK3), and phospholipase C, gamma 1 (PLCG1). Gene coexpression network analysis identified 4 major modules, and the potassium channel tetramerization domain containing 10 (KCTD10), kin of IRRE like (KIRREL), dipeptidyl-peptidase 10 (DPP10), and unc-80 homolog (UNC80) were the hub gene of each modules, respectively.

CONCLUSION

Our integrative analysis provides a comprehensive view of gene expression patterns associated with the pancreatic carcinogenesis.

Differentiated transcriptional signatures in the maize landraces of Chiapas, Mexico

Background

Landrace farmers are the keepers of crops locally adapted to the environments where they are cultivated. Patterns of diversity across the genome can provide signals of past evolution in the face of abiotic and biotic change. Understanding this rich genetic resource is imperative especially since diversity can provide agricultural security as climate continues to shift.

Results

Here we employ RNA sequencing (RNA-seq) to understand the role that conditions that vary across a landscape may have played in shaping genetic diversity in the maize landraces of Chiapas, Mexico. We collected landraces from three distinct elevational zones and planted them in a midland common garden. Early season leaf tissue was collected for RNA-seq and we performed weighted gene co-expression network analysis (WGCNA). We then used association analysis between landrace co-expression module expression values and environmental parameters of landrace origin to elucidate genes and gene networks potentially shaped by environmental factors along our study gradient. Elevation of landrace origin affected the transcriptome profiles. Two co-expression modules were highly correlated with temperature parameters of landrace origin and queries into their ‘hub’ genes suggested that temperature may have led to differentiation among landraces in hormone biosynthesis/signaling and abiotic and biotic stress responses. We identified several ‘hub’ transcription factors and kinases as candidates for the regulation of these responses.

Conclusions

These findings indicate that natural selection may influence the transcriptomes of crop landraces along an elevational gradient in a major diversity center, and provide a foundation for exploring the genetic basis of local adaptation. While we cannot rule out the role of neutral evolutionary forces in the patterns we have identified, combining whole transcriptome sequencing technologies, established bioinformatics techniques, and common garden experimentation can powerfully elucidate structure of adaptive diversity across a varied landscape. Ultimately, gaining such understanding can facilitate the conservation and strategic utilization of crop genetic diversity in a time of climate change.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-4005-y) contains supplementary material, which is available to authorized users.

Transcriptomic profiling in muscle and adipose tissue identifies genes related to growth and lipid deposition

Growth performance and meat quality are important traits for the pig industry and consumers. Adipose tissue is the main site at which fat storage and fatty acid synthesis occur. Therefore, we combined high-throughput transcriptomic sequencing in adipose and muscle tissues with the quantification of corresponding phenotypic features using seven Chinese indigenous pig breeds and one Western commercial breed (Yorkshire). We obtained data on 101 phenotypic traits, from which principal component analysis distinguished two groups: one associated with the Chinese breeds and one with Yorkshire. The numbers of differentially expressed genes between all Chinese breeds and Yorkshire were shown to be 673 and 1056 in adipose and muscle tissues, respectively. Functional enrichment analysis revealed that these genes are associated with biological functions and canonical pathways related to oxidoreductase activity, immune response, and metabolic process. Weighted gene coexpression network analysis found more coexpression modules significantly correlated with the measured phenotypic traits in adipose than in muscle, indicating that adipose regulates meat and carcass quality. Using the combination of differential expression, QTL information, gene significance, and module hub genes, we identified a large number of candidate genes potentially related to economically important traits in pig, which should help us improve meat production and quality.

Temporal network analysis identifies early physiological and transcriptomic indicators of mild drought in Brassica rapa

The dynamics of local climates make development of agricultural strategies challenging. Yield improvement has progressed slowly, especially in drought-prone regions where annual crop production suffers from episodic aridity. Underlying drought responses are circadian and diel control of gene expression that regulate daily variations in metabolic and physiological pathways. To identify transcriptomic changes that occur in the crop Brassica rapa during initial perception of drought, we applied a co-expression network approach to associate rhythmic gene expression changes with physiological responses. Coupled analysis of transcriptome and physiological parameters over a two-day time course in control and drought-stressed plants provided temporal resolution necessary for correlation of network modules with dynamic changes in stomatal conductance, photosynthetic rate, and photosystem II efficiency. This approach enabled the identification of drought-responsive genes based on their differential rhythmic expression profiles in well-watered versus droughted networks and provided new insights into the dynamic physiological changes that occur during drought.

Transcriptional Changes in the Mouse Retina after Ocular Blast Injury: A Role for the Immune System

Ocular blast injury is a major medical concern for soldiers and explosion victims due to poor visual outcomes. To define the changes in gene expression following a blast injury to the eye, we examined retinal ribonucleic acid (RNA) expression in 54 mouse strains 5 days after a single 50-psi overpressure air wave blast injury. We observe that almost 40% of genes are differentially expressed with a false discovery rate (FDR) of <0.001, even though the nominal changes in RNA expression are rather small. Moreover, we find through machine learning approaches that genetic networks related to the innate and acquired immune system are activated. Accompanied by lymphocyte invasion into the inner retina, blast injury also results in progressive loss of visual function and retinal ganglion cells (RGCs). Collectively, these data demonstrate how systems genetics can be used to put meaning to the transcriptome changes following ocular blast injury that eventually lead to blindness.

Co-Expression Network and Pathway Analyses Reveal Important Modules of miRNAs Regulating Milk Yield and Component Traits

Co-expression network analyses provide insights into the molecular interactions underlying complex traits and diseases. In this study, co-expression network analysis was performed to detect expression patterns (modules or clusters) of microRNAs (miRNAs) during lactation, and to identify miRNA regulatory mechanisms for milk yield and component traits (fat, protein, somatic cell count (SCC), lactose, and milk urea nitrogen (MUN)) via miRNA target gene enrichment analysis. miRNA expression (713 miRNAs), and milk yield and components (Fat%, Protein%, lactose, SCC, MUN) data of nine cows at each of six different time points (day 30 (D30), D70, D130, D170, D230 and D290) of an entire lactation curve were used. Four modules or clusters (GREEN, BLUE, RED and TURQUOISE) of miRNAs were identified as important for milk yield and component traits. The GREEN and BLUE modules were significantly correlated (|r| > 0.5) with milk yield and lactose, respectively. The RED and TURQUOISE modules were significantly correlated (|r| > 0.5) with both SCC and lactose. In the GREEN module, three abundantly expressed miRNAs (miR-148a, miR-186 and miR-200a) were most significantly correlated to milk yield, and are probably the most important miRNAs for this trait. DDR1 and DDHX1 are hub genes for miRNA regulatory networks controlling milk yield, while HHEX is an important transcription regulator for these networks. miR-18a, miR-221/222 cluster, and transcription factors HOXA7, and NOTCH 3 and 4, are important for the regulation of lactose. miR-142, miR-146a, and miR-EIA17-14144 (a novel miRNA), and transcription factors in the SMAD family and MYB, are important for the regulation of SCC. Important signaling pathways enriched for target genes of miRNAs of significant modules, included protein kinase A and PTEN signaling for milk yield, eNOS and Noth signaling for lactose, and TGF β, HIPPO, Wnt/β-catenin and cell cycle signaling for SCC. Relevant enriched gene ontology (GO)-terms related to milk and mammary gland traits included cell differentiation, G-protein coupled receptor activity, and intracellular signaling transduction. Overall, this study uncovered regulatory networks in which miRNAs interacted with each other to regulate lactation traits.

Epistatic Networks Jointly Influence Phenotypes Related to Metabolic Disease and Gene Expression in Diversity Outbred Mice

In this study, Tyler et al. analyzed the complex genetic architecture of metabolic disease-related traits using the Diversity Outbred mouse population

Genetic studies of multidimensional phenotypes can potentially link genetic variation, gene expression, and physiological data to create multi-scale models of complex traits. The challenge of reducing these data to specific hypotheses has become increasingly acute with the advent of genome-scale data resources. Multi-parent populations derived from model organisms provide a resource for developing methods to understand this complexity. In this study, we simultaneously modeled body composition, serum biomarkers, and liver transcript abundances from 474 Diversity Outbred mice. This population contained both sexes and two dietary cohorts. Transcript data were reduced to functional gene modules with weighted gene coexpression network analysis (WGCNA), which were used as summary phenotypes representing enriched biological processes. These module phenotypes were jointly analyzed with body composition and serum biomarkers in a combined analysis of pleiotropy and epistasis (CAPE), which inferred networks of epistatic interactions between quantitative trait loci that affect one or more traits. This network frequently mapped interactions between alleles of different ancestries, providing evidence of both genetic synergy and redundancy between haplotypes. Furthermore, a number of loci interacted with sex and diet to yield sex-specific genetic effects and alleles that potentially protect individuals from the effects of a high-fat diet. Although the epistatic interactions explained small amounts of trait variance, the combination of directional interactions, allelic specificity, and high genomic resolution provided context to generate hypotheses for the roles of specific genes in complex traits. Our approach moves beyond the cataloging of single loci to infer genetic networks that map genetic etiology by simultaneously modeling all phenotypes.

Network Analysis of Lung Transcriptomics Reveals a Distinct B-Cell Signature in Emphysema

RATIONALE

Chronic obstructive pulmonary disease (COPD) is characterized by chronic airflow limitation caused by a combination of airways disease (bronchiolitis) and parenchymal destruction (emphysema), whose relative proportion varies from patient to patient.

OBJECTIVES

To explore and contrast the molecular pathogenesis of emphysema and bronchiolitis in COPD.

METHODS

We used network analysis of lung transcriptomics (Affymetrix arrays) in 70 former smokers with COPD to compare differential expression and gene coexpression in bronchiolitis and emphysema.

MEASUREMENTS AND MAIN RESULTS

We observed that in emphysema (but not in bronchiolitis) (1) up-regulated genes were enriched in ontologies related to B-cell homing and activation; (2) the immune coexpression network had a central core of B cell-related genes; (3) B-cell recruitment and immunoglobulin transcription genes (CXCL13, CCL19, and POU2AF1) correlated with emphysema severity; (4) there were lymphoid follicles (CD20(+)IgM(+)) with active B cells (phosphorylated nuclear factor-κB p65(+)), proliferation markers (Ki-67(+)), and class-switched B cells (IgG(+)); and (5) both TNFRSF17 mRNA and B cell-activating factor protein were up-regulated. These findings were by and large reproduced in a group of patients with incipient emphysema and when patients with emphysema were matched for the severity of airflow limitation of those with bronchiolitis.

CONCLUSIONS

Our study identifies enrichment in B cell-related genes in patients with COPD with emphysema that is absent in bronchiolitis. These observations contribute to a better understanding of COPD pathobiology and may open new therapeutic opportunities for patients with COPD.

Preservation affinity in consensus modules among stages of HIV-1 progression

Background

Analysis of gene expression data provides valuable insights into disease mechanism. Investigating relationship among co-expression modules of different stages is a meaningful tool to understand the way in which a disease progresses. Identifying topological preservation of modular structure also contributes to that understanding.

Methods

HIV-1 disease provides a well-documented progression pattern through three stages of infection: acute, chronic and non-progressor. In this article, we have developed a novel framework to describe the relationship among the consensus (or shared) co-expression modules for each pair of HIV-1 infection stages. The consensus modules are identified to assess the preservation of network properties. We have investigated the preservation patterns of co-expression networks during HIV-1 disease progression through an eigengene-based approach.

Results

We discovered that the expression patterns of consensus modules have a strong preservation during the transitions of three infection stages. In particular, it is noticed that between acute and non-progressor stages the preservation is slightly more than the other pair of stages. Moreover, we have constructed eigengene networks for the identified consensus modules and observed the preservation structure among them. Some consensus modules are marked as preserved in two pairs of stages and are analyzed further to form a higher order meta-network consisting of a group of preserved modules. Additionally, we observed that module membership (MM) values of genes within a module are consistent with the preservation characteristics. The MM values of genes within a pair of preserved modules show strong correlation patterns across two infection stages.

Conclusions

We have performed an extensive analysis to discover preservation pattern of co-expression network constructed from microarray gene expression data of three different HIV-1 progression stages. The preservation pattern is investigated through identification of consensus modules in each pair of infection stages. It is observed that the preservation of the expression pattern of consensus modules remains more prominent during the transition of infection from acute stage to non-progressor stage. Additionally, we observed that the module membership values of genes are coherent with preserved modules across the HIV-1 progression stages.

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-017-1590-3) contains supplementary material, which is available to authorized users.

Differential Aging Analysis in Human Cerebral Cortex Identifies Variants in TMEM106B and GRN that Regulate Aging Phenotypes

Human age-associated traits, such as cognitive decline, can be highly variable across the population, with some individuals exhibiting traits that are not expected at a given chronological age. Here we present differential aging (Δ-aging), an unbiased method that quantifies individual variability in age-associated phenotypes within a tissue of interest, and apply this approach to the analysis of existing transcriptome-wide cerebral cortex gene expression data from several cohorts totaling 1,904 autopsied human brain samples. We subsequently performed a genome-wide association study and identified the TMEM106B and GRN gene loci, previously associated with frontotemporal dementia, as determinants of Δ-aging in the cerebral cortex with genome-wide significance. TMEM106B risk variants are associated with inflammation, neuronal loss, and cognitive deficits, even in the absence of known brain disease, and their impact is highly selective for the frontal cerebral cortex of older individuals (>65 years). The methodological framework we describe can be broadly applied to the analysis of quantitative traits associated with aging or with other parameters.

Detection of statistically significant network changes in complex biological networks

Background

Biological networks contribute effectively to unveil the complex structure of molecular interactions and to discover driver genes especially in cancer context. It can happen that due to gene mutations, as for example when cancer progresses, the gene expression network undergoes some amount of localized re-wiring. The ability to detect statistical relevant changes in the interaction patterns induced by the progression of the disease can lead to the discovery of novel relevant signatures. Several procedures have been recently proposed to detect sub-network differences in pairwise labeled weighted networks.

Methods

In this paper, we propose an improvement over the state-of-the-art based on the Generalized Hamming Distance adopted for evaluating the topological difference between two networks and estimating its statistical significance. The proposed procedure exploits a more effective model selection criteria to generate p-values for statistical significance and is more efficient in terms of computational time and prediction accuracy than literature methods. Moreover, the structure of the proposed algorithm allows for a faster parallelized implementation.

Results

In the case of dense random geometric networks the proposed approach is 10-15x faster and achieves 5-10% higher AUC, Precision/Recall, and Kappa value than the state-of-the-art. We also report the application of the method to dissect the difference between the regulatory networks of IDH-mutant versus IDH-wild-type glioma cancer. In such a case our method is able to identify some recently reported master regulators as well as novel important candidates.

Conclusions

We show that our network differencing procedure can effectively and efficiently detect statistical significant network re-wirings in different conditions. When applied to detect the main differences between the networks of IDH-mutant and IDH-wild-type glioma tumors, it correctly selects sub-networks centered on important key regulators of these two different subtypes. In addition, its application highlights several novel candidates that cannot be detected by standard single network-based approaches.

Electronic supplementary material

The online version of this article (doi:10.1186/s12918-017-0412-6) contains supplementary material, which is available to authorized users.

Analysis of long noncoding RNA and mRNA using RNA sequencing during the differentiation of intramuscular preadipocytes in chicken

Long noncoding RNAs (lncRNAs) regulate metabolic tissue development and function, including adipogenesis. However, little is known about the function and profile of lncRNAs in intramuscular preadipocyte differentiation in chicken. Here, we identified lncRNAs in chicken intramuscular preadipocytes at different differentiation stages using RNA sequencing. A total of 1,311,382,604 clean reads and 25,435 lncRNAs were obtained from 12 samples. In total, 7,433 differentially expressed genes (4,698 lncRNAs and 2,735 mRNAs) were identified by pairwise comparison. These 7,433 differentially expressed genes were grouped into 11 clusters based on their expression patterns by K-means clustering. Using Weighted Gene Coexpression Network Analysis, we identified four stage-specific modules positively related to I0, I2, I4, and I6 stages and two stage-specific modules negatively related to I0 and I2 stages, respectively. Many well-known and novel pathways associated with intramuscular preadipocyte differentiation were identified. We also identified hub genes in each stage-specific module and visualized them in Cytoscape. Our analysis revealed many highly-connected genes, including XLOC_058593, BMP3, MYOD1, and LAMP3. This study provides a valuable resource for chicken lncRNA study and improves our understanding of the biology of preadipocyte differentiation in chicken.

A Systems Approach Implicates a Brain Mitochondrial Oxidative Homeostasis Co-expression Network in Genetic Vulnerability to Alcohol Withdrawal

Genetic factors significantly affect vulnerability to alcohol dependence (alcoholism). We previously identified quantitative trait loci on distal mouse chromosome 1 with large effects on predisposition to alcohol physiological dependence and associated withdrawal following both chronic and acute alcohol exposure in mice (Alcdp1 and Alcw1, respectively). We fine-mapped these loci to a 1.1–1.7 Mb interval syntenic with human 1q23.2-23.3. Alcw1/Alcdp1 interval genes show remarkable genetic variation among mice derived from the C57BL/6J and DBA/2J strains, the two most widely studied genetic animal models for alcohol-related traits. Here, we report the creation of a novel recombinant Alcw1/Alcdp1 congenic model (R2) in which the Alcw1/Alcdp1 interval from a donor C57BL/6J strain is introgressed onto a uniform, inbred DBA/2J genetic background. As expected, R2 mice demonstrate significantly less severe alcohol withdrawal compared to wild-type littermates. Additionally, comparing R2 and background strain animals, as well as reciprocal congenic (R8) and appropriate background strain animals, we assessed Alcw1/Alcdp1 dependent brain gene expression using microarray and quantitative PCR analyses. To our knowledge this includes the first Weighted Gene Co-expression Network Analysis using reciprocal congenic models. Importantly, this allows detection of co-expression patterns limited to one or common to both genetic backgrounds with high or low predisposition to alcohol withdrawal severity. The gene expression patterns (modules) in common contain genes related to oxidative phosphorylation, building upon human and animal model studies that implicate involvement of oxidative phosphorylation in alcohol use disorders (AUDs). Finally, we demonstrate that administration of N-acetylcysteine, an FDA-approved antioxidant, significantly reduces symptoms of alcohol withdrawal (convulsions) in mice, thus validating a phenotypic role for this network. Taken together, these studies support the importance of mitochondrial oxidative homeostasis in alcohol withdrawal and identify this network as a valuable therapeutic target in human AUDs.

WGCNA Application to Proteomic and Metabolomic Data Analysis

Progresses in mass spectrometric instrumentation and bioinformatics identification algorithms made over the past decades allow quantitative measurements of relative or absolute protein/metabolite amounts in cells in a high-throughput manner, which has significantly expedited the exploration into functions and dynamics of complex biological systems. However, interpretation of high-throughput data is often restricted by the limited availability of suitable computational methods and enough statistical power. While many computational methodologies have been developed in the past decades to address the issue, it becomes clear that network-focused rather than individual gene/protein-focused strategies would be more appropriate to obtain a complete picture of cellular responses. Recently, an R analytical package named as weighted gene coexpression network analysis (WGCNA) was developed and applied to high-throughput microarray or RNA-seq datasets since it provides a systems-level insights, high sensitivity to low abundance, or small fold changes genes without any information loss. The approach was also recently applied to proteomic and metabolomic data analysis. However, due to the fact that low coverage of the current proteomic and metabolomic analytical technologies, causing the format of datasets are often incomplete, the method needs to be modified so that it can be properly utilized for meaningful biologically interpretation. In this chapter, we provide a detailed introduction of the modified protocol and its tutorials for applying the WGCNA approach in analyzing proteomic and metabolomic datasets.

The bronchial epithelial cell bacterial microbiome and host response in patients infected with human immunodeficiency virus

BACKGROUND

Chronic Obstructive Pulmonary Disease (COPD) is an important comorbidity in patients living with human immunodeficiency virus (HIV). Previous bacterial microbiome studies have shown increased abundance of specific bacterium, like Tropheryma whipplei, and no overall community differences. However, the host response to the lung microbiome is unknown in patients infected with HIV.

METHODS

Two bronchial brush samples were obtained from 21 HIV-infected patients. One brush was used for bacterial microbiome analysis using the Illumina MiSeq^TM platform, while the other was used to evaluate gene expression patterns of the host using the Affymetrix Human Gene ST 2.0 array. Weighted gene co-expression network analysis was used to determine the relationship between the bacterial microbiome and host gene expression response.

RESULTS

The Shannon Diversity was inversely related to only one gene expression module (p = 0.02); whereas evenness correlated with five different modules (p ≤ 0.05). After FDR correction only the Firmicutes phylum was significantly correlated with any modules (FDR < 0.05). These modules were enriched for cilia, transcription regulation, and immune response. Specific operational taxonomic units (OTUs), such as OTU4 (Pasteurellaceae), were able to distinguish HIV patients with and without COPD and severe emphysema.

CONCLUSION

These data support the hypothesis that the bacterial microbiome in HIV lungs is associated with specific host immune responses. Whether or not these responses are also seen in non-HIV infected individuals needs to be addressed in future studies.

Extracting Stage-Specific and Dynamic Modules Through Analyzing Multiple Networks Associated with Cancer Progression

Determining the dynamics of pathways associated with cancer progression is critical for understanding the etiology of diseases. Advances in biological technology have facilitated the simultaneous genomic profiling of multiple patients at different clinical stages, thus generating the dynamic genomic data for cancers. Such data provide enable investigation of the dynamics of related pathways. However, methods for integrative analysis of dynamic genomic data are inadequate. In this study, we develop a novel nonnegative matrix factorization algorithm for dynamic modules ( NMF-DM), which simultaneously analyzes multiple networks for the identification of stage-specific and dynamic modules. NMF-DM applies the temporal smoothness framework by balancing the networks at the current stage and the previous stage. Experimental results indicate that the NMF-DM algorithm is more accurate than the state-of-the-art methods in artificial dynamic networks. In breast cancer networks, NMF-DM reveals the dynamic modules that are important for cancer stage transitions. Furthermore, the stage-specific and dynamic modules have distinct topological and biochemical properties. Finally, we demonstrate that the stage-specific modules significantly improve the accuracy of cancer stage prediction. The proposed algorithm provides an effective way to explore the time-dependent cancer genomic data.

Normalized lmQCM: An Algorithm for Detecting Weak Quasi-Cliques in Weighted Graph with Applications in Gene Co-Expression Module Discovery in Cancers

In this paper, we present a new approach for mining weighted networks to identify densely connected modules such as quasi-cliques. Quasi-cliques are densely connected subnetworks in a network. Detecting quasi-cliques is an important topic in data mining, with applications such as social network study and biomedicine. Our approach has two major improvements upon previous work. The first is the use of local maximum edges to initialize the search in order to avoid excessive overlaps among the modules, thereby greatly reducing the computing time. The second is the inclusion of a weight normalization procedure to enable discovery of "subtle" modules with more balanced sizes. We carried out careful tests on multiple parameters and settings using two large cancer datasets. This approach allowed us to identify a large number of gene modules enriched in both biological functions and chromosomal bands in cancer data, suggesting potential roles of copy number variations (CNVs) involved in the cancer development. We then tested the genes in selected modules with enriched chromosomal bands using The Cancer Genome Atlas data, and the results strongly support our hypothesis that the coexpression in these modules are associated with CNVs. While gene coexpression network analyses have been widely adopted in disease studies, most of them focus on the functional relationships of coexpressed genes. The relationship between coexpression gene modules and CNVs are much less investigated despite the potential advantage that we can infer from such relationship without genotyping data. Our new approach thus provides a means to carry out deep mining of the gene coexpression network to obtain both functional and genetic information from the expression data.

Reversal of dendritic phenotypes in 16p11.2 microduplication mouse model neurons by pharmacological targeting of a network hub

The architecture of dendritic arbors contributes to neuronal connectivity in the brain. Conversely, abnormalities in dendrites have been reported in multiple mental disorders and are thought to contribute to pathogenesis. Rare copy number variations (CNVs) are genetic alterations that are associated with a wide range of mental disorders and are highly penetrant. The 16p11.2 microduplication is one of the CNVs most strongly associated with schizophrenia and autism, spanning multiple genes possibly involved in synaptic neurotransmission. However, disease-relevant cellular phenotypes of 16p11.2 microduplication and the driver gene(s) remain to be identified. We found increased dendritic arborization in isolated cortical pyramidal neurons from a mouse model of 16p11.2 duplication (dp/+). Network analysis identified MAPK3, which encodes ERK1 MAP kinase, as the most topologically important hub in protein-protein interaction networks within the 16p11.2 region and broader gene networks of schizophrenia-associated CNVs. Pharmacological targeting of ERK reversed dendritic alterations associated with dp/+ neurons, outlining a strategy for the analysis and reversal of cellular phenotypes in CNV-related psychiatric disorders.

The Association Between Low Back Pain and Composition of IgG Glycome

Low back pain (LBP) is a common debilitating condition which aetiology and pathogenesis are poorly understood. We carried out a first so far analysis of associations between LBP and plasma IgG N-glycome in a sample of 4511 twins from TwinsUK database assessed for LBP, lumbar disc degeneration (LDD) as its possible cause, and IgG-glycan levels. Using weighted correlation network analysis, we established a correlation between LBP and glycan modules featured by glycans that either promote or block antibody-dependent cell-mediated cytotoxicity (ADCC). The levels of four glycan traits representing two of those modules were statistically significantly different in monozygotic twins discordant for LBP. Also, the trend to higher prevalence of systemic inflammatory disorders was shown for twins with low level of fucosylated glycans and high level of non-fucosylated glycans. Core fucosylation of IgG is a “safety switch” reducing ADCC, thus our results suggest the involvement of ADCC and associated inflammation in pathogenesis of LBP. No correlation between LDD scores and glycans was found assuming that the inflammation may not be a part of LDD. These data provide a new insight into understanding the complex pathophysiology of LBP and suggest glycan levels as a possible biomarker for inflammation-related subtypes of LBP.

Transcriptomic plasticity in brown adipose tissue contributes to an enhanced capacity for nonshivering thermogenesis in deer mice

For small mammals living at high altitude, aerobic heat generation (thermogenesis) is essential for survival during prolonged periods of cold, but is severely impaired under conditions of hypobaric hypoxia. Recent studies in deer mice (Peromyscus maniculatus) reveal adaptive enhancement of thermogenesis in high- compared to low-altitude populations under hypoxic cold stress, an enhancement that is attributable to modifications in the aerobic metabolism of muscles used in shivering. However, because small mammals rely heavily on nonshivering mechanisms for cold acclimatization, we tested for evidence of adaptive divergence in nonshivering thermogenesis (NST) under hypoxia. To do so, we measured NST and characterized transcriptional profiles of brown adipose tissue (BAT) in high- and low-altitude deer mice that were (i) wild-caught and acclimatized to their native altitude, and (ii) born and reared under common garden conditions at low elevation. We found that NST performance under hypoxia is enhanced in wild-caught, high-altitude deer mice, a difference that is associated with increased expression of coregulated genes that influence several physiological traits. These traits include vascularization and O2 supply to BAT, brown adipocyte proliferation and the uncoupling of oxidative phosphorylation from ATP synthesis in the generation of heat. Our results suggest that acclimatization to hypoxic cold stress is facilitated by enhancement of nonshivering heat production, which is driven by regulatory plasticity in a suite of genes that influence intersecting physiological pathways.

Differential gene network analysis for the identification of asthma-associated therapeutic targets in allergen-specific T-helper memory responses

BACKGROUND

Asthma is strongly associated with allergic sensitization, but the mechanisms that determine why only a subset of atopics develop asthma are not well understood. The aim of this study was to test the hypothesis that variations in allergen-driven CD4 T cell responses are associated with susceptibility to expression of asthma symptoms.

METHODS

The study population consisted of house dust mite (HDM) sensitized atopics with current asthma (n = 22), HDM-sensitized atopics without current asthma (n = 26), and HDM-nonsensitized controls (n = 24). Peripheral blood mononuclear cells from these groups were cultured in the presence or absence of HDM extract for 24 h. CD4 T cells were then isolated by immunomagnetic separation, and gene expression patterns were profiled on microarrays.

RESULTS

Differential network analysis of HDM-induced CD4 T cell responses in sensitized atopics with or without asthma unveiled a cohort of asthma-associated genes that escaped detection by more conventional data analysis techniques. These asthma-associated genes were enriched for targets of STAT6 signaling, and they were nested within a larger coexpression module comprising 406 genes. Upstream regulator analysis suggested that this module was driven primarily by IL-2, IL-4, and TNF signaling; reconstruction of the wiring diagram of the module revealed a series of hub genes involved in inflammation (IL-1B, NFkB, STAT1, STAT3), apoptosis (BCL2, MYC), and regulatory T cells (IL-2Ra, FoxP3). Finally, we identified several negative regulators of asthmatic CD4 T cell responses to allergens (e.g. IL-10, type I interferons, microRNAs, drugs, metabolites), and these represent logical candidates for therapeutic intervention.

CONCLUSION

Differential network analysis of allergen-induced CD4 T cell responses can unmask covert disease-associated genes and pin point novel therapeutic targets.