Using effective subnetworks to predict selected properties of gene networks

Background

Difficulties associated with implementing gene therapy are caused by the complexity of the underlying regulatory networks. The forms of interactions between the hundreds of genes, proteins, and metabolites in these networks are not known very accurately. An alternative approach is to limit consideration to genes on the network. Steady state measurements of these influence networks can be obtained from DNA microarray experiments. However, since they contain a large number of nodes, the computation of influence networks requires a prohibitively large set of microarray experiments. Furthermore, error estimates of the network make verifiable predictions impossible.

Methodology/Principal Findings

Here, we propose an alternative approach. Rather than attempting to derive an accurate model of the network, we ask what questions can be addressed using lower dimensional, highly simplified models. More importantly, is it possible to use such robust features in applications? We first identify a small group of genes that can be used to affect changes in other nodes of the network. The reduced effective empirical subnetwork (EES) can be computed using steady state measurements on a small number of genetically perturbed systems. We show that the EES can be used to make predictions on expression profiles of other mutants, and to compute how to implement pre-specified changes in the steady state of the underlying biological process. These assertions are verified in a synthetic influence network. We also use previously published experimental data to compute the EES associated with an oxygen deprivation network of E.coli, and use it to predict gene expression levels on a double mutant. The predictions are significantly different from the experimental results for less than of genes.

Conclusions/Significance

The constraints imposed by gene expression levels of mutants can be used to address a selected set of questions about a gene network.     

Uncovering genomic causes of co-morbidity in epilepsy: gene-driven phenotypic characterization of rare microdeletions

Background

Patients with epilepsy often suffer from other important conditions. The existence of such co-morbidities is frequently not recognized and their relationship with epilepsy usually remains unexplained.

Methodology/Principal Findings

We describe three patients with common, sporadic, non-syndromic epilepsies in whom large genomic microdeletions were found during a study of genetic susceptibility to epilepsy. We performed detailed gene-driven clinical investigations in each patient. Disruption of the function of genes in the deleted regions can explain co-morbidities in these patients.

Conclusions/Significance

Co-morbidities in patients with epilepsy can be part of a genomic abnormality even in the absence of (known) congenital malformations or intellectual disabilities. Gene-driven phenotype examination can also reveal clinically significant unsuspected condition.     

Reducing amplification artifacts in high multiplex amplicon sequencing by using molecular barcodes

Background

PCR amplicon sequencing has been widely used as a targeted approach for both DNA and RNA sequence analysis. High multiplex PCR has further enabled the enrichment of hundreds of amplicons in one simple reaction. At the same time, the performance of PCR amplicon sequencing can be negatively affected by issues such as high duplicate reads, polymerase artifacts and PCR amplification bias. Recently researchers have made some good progress in addressing these shortcomings by incorporating molecular barcodes into PCR primer design. So far, most work has been demonstrated using one to a few pairs of primers, which limits the size of the region one can analyze.

Results

We developed a simple protocol, which enables the use of molecular barcodes in high multiplex PCR with hundreds of amplicons. Using this protocol and reference materials, we demonstrated the applications in accurate variant calling at very low fraction over a large region and in targeted RNA quantification. We also evaluated the protocol’s utility in profiling FFPE samples.

Conclusions

We demonstrated the successful implementation of molecular barcodes in high multiplex PCR, with multiplex scale many times higher than earlier work. We showed that the new protocol combines the benefits of both high multiplex PCR and molecular barcodes, i.e. the analysis of a very large region, low DNA input requirement, very good reproducibility and the ability to detect as low as 1 % mutations with minimal false positives (FP).

Electronic supplementary material

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

Inferring gene-phenotype associations via global protein complex network propagation

Background

Phenotypically similar diseases have been found to be caused by functionally related genes, suggesting a modular organization of the genetic landscape of human diseases that mirrors the modularity observed in biological interaction networks. Protein complexes, as molecular machines that integrate multiple gene products to perform biological functions, express the underlying modular organization of protein-protein interaction networks. As such, protein complexes can be useful for interrogating the networks of phenome and interactome to elucidate gene-phenotype associations of diseases.

Methodology/Principal Findings

We proposed a technique called RWPCN (Random Walker on Protein Complex Network) for predicting and prioritizing disease genes. The basis of RWPCN is a protein complex network constructed using existing human protein complexes and protein interaction network. To prioritize candidate disease genes for the query disease phenotypes, we compute the associations between the protein complexes and the query phenotypes in their respective protein complex and phenotype networks. We tested RWPCN on predicting gene-phenotype associations using leave-one-out cross-validation; our method was observed to outperform existing approaches. We also applied RWPCN to predict novel disease genes for two representative diseases, namely, Breast Cancer and Diabetes.

Conclusions/Significance

Guilt-by-association prediction and prioritization of disease genes can be enhanced by fully exploiting the underlying modular organizations of both the disease phenome and the protein interactome. Our RWPCN uses a novel protein complex network as a basis for interrogating the human phenome-interactome network. As the protein complex network can capture the underlying modularity in the biological interaction networks better than simple protein interaction networks, RWPCN was found to be able to detect and prioritize disease genes better than traditional approaches that used only protein-phenotype associations.     

Genomic Responses during Acute Human Anaphylaxis Are Characterized by Upregulation of Innate Inflammatory Gene Networks

Background

Systemic spread of immune activation and mediator release is required for the development of anaphylaxis in humans. We hypothesized that peripheral blood leukocyte (PBL) activation plays a key role.

Objective

To characterize PBL genomic responses during acute anaphylaxis.

Methods

PBL samples were collected at three timepoints from six patients presenting to the Emergency Department (ED) with acute anaphylaxis and six healthy controls. Gene expression patterns were profiled on microarrays, differentially expressed genes were identified, and network analysis was employed to explore underlying mechanisms.

Results

Patients presented with moderately severe anaphylaxis after oral aspirin (2), peanut (2), bee sting (1) and unknown cause (1). Two genes were differentially expressed in patients compared to controls at ED arrival, 67 genes at 1 hour post-arrival and 2,801 genes at 3 hours post-arrival. Network analysis demonstrated that three inflammatory modules were upregulated during anaphylaxis. Notably, these modules contained multiple hub genes, which are known to play a central role in the regulation of innate inflammatory responses. Bioinformatics analyses showed that the data were enriched for LPS-like and TNF activation signatures.

Conclusion

PBL genomic responses during human anaphylaxis are characterized by dynamic expression of innate inflammatory modules. Upregulation of these modules was observed in patients with different reaction triggers. Our findings indicate a role for innate immune pathways in the pathogenesis of human anaphylaxis, and the hub genes identified in this study represent logical candidates for follow-up studies.     

Enrichment Map: A Network-Based Method for Gene-Set Enrichment Visualization and Interpretation

Background

Gene-set enrichment analysis is a useful technique to help functionally characterize large gene lists, such as the results of gene expression experiments. This technique finds functionally coherent gene-sets, such as pathways, that are statistically over-represented in a given gene list. Ideally, the number of resulting sets is smaller than the number of genes in the list, thus simplifying interpretation. However, the increasing number and redundancy of gene-sets used by many current enrichment analysis software works against this ideal.

Principal Findings

To overcome gene-set redundancy and help in the interpretation of large gene lists, we developed “Enrichment Map”, a network-based visualization method for gene-set enrichment results. Gene-sets are organized in a network, where each set is a node and edges represent gene overlap between sets. Automated network layout groups related gene-sets into network clusters, enabling the user to quickly identify the major enriched functional themes and more easily interpret the enrichment results.

Conclusions

Enrichment Map is a significant advance in the interpretation of enrichment analysis. Any research project that generates a list of genes can take advantage of this visualization framework. Enrichment Map is implemented as a freely available and user friendly plug-in for the Cytoscape network visualization software (http://baderlab.org/Software/EnrichmentMap/).     

Expression Profiling of a Genetic Animal Model of Depression Reveals Novel Molecular Pathways Underlying Depressive-Like Behaviours

Background

The Flinders model is a validated genetic rat model of depression that exhibits a number of behavioural, neurochemical and pharmacological features consistent with those observed in human depression.

Principal Findings

In this study we have used genome-wide microarray expression profiling of the hippocampus and prefrontal/frontal cortex of Flinders Depression Sensitive (FSL) and control Flinders Depression Resistant (FRL) lines to understand molecular basis for the differences between the two lines. We profiled two independent cohorts of Flinders animals derived from the same colony six months apart, each cohort statistically powered to allow independent as well as combined analysis. Using this approach, we were able to validate using real-time-PCR a core set of gene expression differences that showed statistical significance in each of the temporally distinct cohorts, representing consistently maintained features of the model. Small but statistically significant increases were confirmed for cholinergic (chrm2, chrna7) and serotonergic receptors (Htr1a, Htr2a) in FSL rats consistent with known neurochemical changes in the model. Much larger gene changes were validated in a number of novel genes as exemplified by TMEM176A, which showed 35-fold enrichment in the cortex and 30-fold enrichment in hippocampus of FRL animals relative to FSL.

Conclusions

These data provide significant insights into the molecular differences underlying the Flinders model, and have potential relevance to broader depression research.     

Regulation of gene expression with thyroid hormone in rats with myocardial infarction

Introduction

The expression of hundreds of genes is altered in response to left ventricular (LV) remodeling following large transmural myocardial infarction (MI). Thyroid hormone (TH) improves LV remodeling and cardiac performance after MI. However, the molecular basis is unknown.

Methods

MI was produced by ligation of the left anterior descending coronary artery in female SD rats. Rats were divided into the following groups: (1) Sham MI, (2) MI, and (3) MI+T4 treatment (T4 pellet 3.3 mg, 60 days release, implanted subcutaneously immediately following MI). Four weeks after surgery, total RNA was isolated from LV non-infarcted areas for microarray analysis using the Illumina RatRef-12 Expression BeadChip Platform.

Results

Signals were detected in 13,188 genes (out of 22,523), of which the expression of 154 genes were decreased and the expression of 200 genes were increased in MI rats compared with Sham MI rats (false discovery rate (FDR) <0.05). Compared to MI rats, T4 treatment decreased expression of 27 genes and increased expression of 28 genes. In particular, 6 genes down-regulated by MI and 12 genes up-regulated by MI were reversed by T4. Most of the 55 genes altered by T4 treatment are in the category of molecular function under binding (24) and biological processes which includes immune system process (9), multi-organism process (5) and biological regulation (19) nonexclusively.

Conclusions

These results suggest that altered expression of genes for molecular function and biological process may be involved in the beneficial effects of thyroid hormone treatment following MI in rats.     

Categorizer: a tool to categorize genes into user-defined biological groups based on semantic similarity

Background

Communalities between large sets of genes obtained from high-throughput experiments are often identified by searching for enrichments of genes with the same Gene Ontology (GO) annotations. The GO analysis tools used for these enrichment analyses assume that GO terms are independent and the semantic distances between all parent–child terms are identical, which is not true in a biological sense. In addition these tools output lists of often redundant or too specific GO terms, which are difficult to interpret in the context of the biological question investigated by the user. Therefore, there is a demand for a robust and reliable method for gene categorization and enrichment analysis.

Results

We have developed Categorizer, a tool that classifies genes into user-defined groups (categories) and calculates p-values for the enrichment of the categories. Categorizer identifies the biologically best-fit category for each gene by taking advantage of a specialized semantic similarity measure for GO terms. We demonstrate that Categorizer provides improved categorization and enrichment results of genetic modifiers of Huntington’s disease compared to a classical GO Slim-based approach or categorizations using other semantic similarity measures.

Conclusion

Categorizer enables more accurate categorizations of genes than currently available methods. This new tool will help experimental and computational biologists analyzing genomic and proteomic data according to their specific needs in a more reliable manner.     

A network-based gene expression signature informs prognosis and treatment for colorectal cancer patients

Background

Several studies have reported gene expression signatures that predict recurrence risk in stage II and III colorectal cancer (CRC) patients with minimal gene membership overlap and undefined biological relevance. The goal of this study was to investigate biological themes underlying these signatures, to infer genes of potential mechanistic importance to the CRC recurrence phenotype and to test whether accurate prognostic models can be developed using mechanistically important genes.

Methods and Findings

We investigated eight published CRC gene expression signatures and found no functional convergence in Gene Ontology enrichment analysis. Using a random walk-based approach, we integrated these signatures and publicly available somatic mutation data on a protein-protein interaction network and inferred 487 genes that were plausible candidate molecular underpinnings for the CRC recurrence phenotype. We named the list of 487 genes a NEM signature because it integrated information from Network, Expression, and Mutation. The signature showed significant enrichment in four biological processes closely related to cancer pathophysiology and provided good coverage of known oncogenes, tumor suppressors, and CRC-related signaling pathways. A NEM signature-based Survival Support Vector Machine prognostic model was trained using a microarray gene expression dataset and tested on an independent dataset. The model-based scores showed a 75.7% concordance with the real survival data and separated patients into two groups with significantly different relapse-free survival (p = 0.002). Similar results were obtained with reversed training and testing datasets (p = 0.007). Furthermore, adjuvant chemotherapy was significantly associated with prolonged survival of the high-risk patients (p = 0.006), but not beneficial to the low-risk patients (p = 0.491).

Conclusions

The NEM signature not only reflects CRC biology but also informs patient prognosis and treatment response. Thus, the network-based data integration method provides a convergence between biological relevance and clinical usefulness in gene signature development.     

Long-Term Effects of Temporal Lobe Epilepsy on Local Neural Networks: A Graph Theoretical Analysis of Corticography Recordings

Purpose

Pharmaco-resistant temporal lobe epilepsy (TLE) is often treated with surgical intervention at some point. As epilepsy surgery is considered a last resort by most physicians, a long history of epileptic seizures prior to surgery is not uncommon. Little is known about the effects of ongoing TLE on neural functioning. A better understanding of these effects might influence the moment of surgical intervention. Functional connectivity (interaction between spatially distributed brain areas) and network structure (integration and segregation of information processing) are thought to be essential for optimal brain functioning. We report on the impact of TLE duration on temporal lobe functional connectivity and network characteristics.

Methods

Functional connectivity of the temporal lobe at the time of surgery was assessed by means of interictal electrocorticography (ECoG) recordings of 27 TLE patients by using the phase lag index (PLI). Graphs (abstract network representations) were reconstructed from the PLI matrix and characterized by the clustering coefficient C (local clustering), the path length L (overall network interconnectedness), and the “small world index” S (network configuration).

Results

Functional connectivity (average PLI), clustering coefficients, and the small world index were negatively correlated with TLE duration in the broad frequency band (0.5–48 Hz).

Discussion

Temporal lobe functional connectivity is lower in patients with longer TLE history, and longer TLE duration is correlated with more random network configuration. Our findings suggest that the neural networks of TLE patients become more pathological over time, possibly due to temporal lobe changes associated with long-standing lesional epilepsy.     

The Role of Resting State Networks in Focal Neocortical Seizures

Objective

The role of resting state functional networks in epilepsy is incompletely understood. While some pathologic diagnoses have been shown to have maintained but altered resting state connectivity, others have implicated resting state connectivity in disease progression. However little is known about how these resting state networks influence the behavior of a focal neocortical seizure.

Methods

Using data taken from invasively monitored patients with intractable focal neocortical epilepsy, we evaluated network connectivity (as determined by oscillatory covariance of the slow cortical potential (<0.5 Hz)) as it relates to neocortical seizure foci both in the interictal and ictal states.

Results

Similar to what has been shown in the past for sleep and anesthesia, electophysiologic resting state networks that are defined by this slow cortical potential covariance maintain their topographic correlation structure throughout an ictal event. Moreover, in the context of focal epilepsy in which the seizure has a specific site of onset, seizure propagation is not chaotic or random. Rather, the seizure (reflected by an elevation of high frequency power) preferentially propagates along the network that contains the seizure onset zone.

Significance

Taken together, these findings further undergird the fundamental role of resting state networks, provide novel insights into the network-influenced behavior of seizures, and potentially identify additional targets for surgical disconnection including informing the location for the completion of multiple subpial transections (MSPTs).