Gene networks and haloperidol-induced catalepsy

Coexpression of the discoidin domain receptor 1 gene with oligodendrocyte-related and schizophrenia risk genes in the developing and adult human brain

BACKGROUND

Discoidin domain receptor tyrosine kinase 1 (DDR1) is present in multiple types of epithelial cells and is highly expressed in the nervous system. Previous studies have revealed that DDR1 is involved in schizophrenia (SCZ). Although the expression of DDR1 in oligodendrocytes has been described, its role in brain myelination is not well understood. In this study, we aimed to explore the coexpression network of DDR1 in the human brain and to compare the list of DDR1 coexpressing genes with the list of genes containing single nucleotide polymorphisms (SNPs) that are associated with SCZ.

MATERIALS AND METHODS

We used a weighted gene coexpression network analysis (WGCNA) of a dataset from four brain areas (the dorsolateral prefrontal cortex, primary motor cortex, hippocampus, and striatum) and from four different intervals (I) of life (I-1 = 10-38 weeks postconception, I-2 ≥0 to < 6 years, I-3 ≥ 6 to < 40 years, and I-4 ≥ 40 years of age). We compared the list of genes that are associated with SCZ in the GWAS Catalog with the list of genes coexpressing with DDR1 in each interval.

RESULTS

Our study revealed that DDR1 was coexpressed with oligodendrocyte-related genes mainly in I-2 (adjP = 5.66e-24) and I-3 (adjP = 2.8e-114), which coincided with the coexpression of DDR1 with myelination-related genes (adjP = 9.04e-03 and 2.51e-08, respectively). DDR1 was also coexpressed with astrocyte-related genes in I-1 (adjP = 1.11e-71), I-2 (adjP = 2.12e-20) and I-4 (adjP = 9.93e-52) and with type 2 microglia-related genes in I-1 (adjP = 2.84e-08), I-2 (adjP = 5.68e-16) and I-4 (adjP = 3.66e-10). Moreover, we observed significant enrichment of SCZ susceptibility genes within the coexpression modules containing DDR1 in I-1 and I-4 (P = 1e-04 and 0.0037, respectively), during which the DDR1 module showed the highest association with the astrocytes.

CONCLUSIONS

Our study confirmed that DDR1 is coexpressed with oligodendrocyte- and myelin-related genes in the human brain but suggests that DDR1 may contribute mainly to SCZ risk through its role in other glial cell types, such as astrocytes.

Brain gene co-expression networks link complement signaling with convergent synaptic pathology in schizophrenia

The most significant common variant association for schizophrenia (SCZ) reflects increased expression of the complement component 4A (C4A). Yet, it remains unclear how C4A interacts with other SCZ risk genes or whether the complement system more broadly is implicated in SCZ pathogenesis. Here, we integrate several existing, large-scale genetic and transcriptomic datasets to interrogate the functional role of the complement system and C4A in the human brain. Unexpectedly, we find no significant genetic enrichment among known complement system genes for SCZ. Conversely, brain co-expression network analyses using C4A as a seed gene reveal that genes downregulated when C4A expression increases exhibit strong and specific genetic enrichment for SCZ risk. This convergent genomic signal reflects synaptic processes, is sexually dimorphic and most prominent in frontal cortical brain regions, and is accentuated by smoking. Overall, these results indicate that synaptic pathways-rather than the complement system-are the driving force conferring SCZ risk.

Alignment of the transcriptome with individual variation in animals selectively bred for High Drinking-In-the-Dark (HDID)

Among animals at risk for excessive ethanol consumption such as the HDID selected mouse lines, there is considerable individual variation in the amount of ethanol consumed and the associated blood ethanol concentrations (BECs). For the HDID lines, this variation occurs even though the residual genetic variation associated with the DID phenotype has been largely exhausted and thus is most likely associated with epigenetic factors. Here we focus on the question of whether the genes associated with individual variation in HDID-1 mice are different from those associated with selection (risk) (Iancu et al., 2013). Thirty-three HDID-1 mice were phenotyped for their BECs at the end of a standard DID trial, were sacrificed 3 weeks later, and RNA-Seq was used to analyze the striatal transcriptome. The data obtained illustrate that there is considerable overlap of the risk and variation gene sets, both focused on the fine-tuning of synaptic plasticity.

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.

Integrative Analysis of Genetic, Genomic, and Phenotypic Data for Ethanol Behaviors: A Network-Based Pipeline for Identifying Mechanisms and Potential Drug Targets

Complex behavioral traits, such as alcohol abuse, are caused by an interplay of genetic and environmental factors, producing deleterious functional adaptations in the central nervous system. The long-term behavioral consequences of such changes are of substantial cost to both the individual and society. Substantial progress has been made in the last two decades in understanding elements of brain mechanisms underlying responses to ethanol in animal models and risk factors for alcohol use disorder (AUD) in humans. However, treatments for AUD remain largely ineffective and few medications for this disease state have been licensed. Genome-wide genetic polymorphism analysis (GWAS) in humans, behavioral genetic studies in animal models and brain gene expression studies produced by microarrays or RNA-seq have the potential to produce nonbiased and novel insight into the underlying neurobiology of AUD. However, the complexity of such information, both statistical and informational, has slowed progress toward identifying new targets for intervention in AUD. This chapter describes one approach for integrating behavioral, genetic, and genomic information across animal model and human studies. The goal of this approach is to identify networks of genes functioning in the brain that are most relevant to the underlying mechanisms of a complex disease such as AUD. We illustrate an example of how genomic studies in animal models can be used to produce robust gene networks that have functional implications, and to integrate such animal model genomic data with human genetic studies such as GWAS for AUD. We describe several useful analysis tools for such studies: ComBAT, WGCNA, and EW_dmGWAS. The end result of this analysis is a ranking of gene networks and identification of their cognate hub genes, which might provide eventual targets for future therapeutic development. Furthermore, this combined approach may also improve our understanding of basic mechanisms underlying gene x environmental interactions affecting brain functioning in health and disease.

NeuN+ neuronal nuclei in non-human primate prefrontal cortex and subcortical white matter after clozapine exposure

Increased neuronal densities in subcortical white matter have been reported for some cases with schizophrenia. The underlying cellular and molecular mechanisms remain unresolved. We exposed 26 young adult macaque monkeys for 6 months to either clozapine, haloperidol or placebo and measured by structural MRI frontal gray and white matter volumes before and after treatment, followed by observer-independent, flow-cytometry-based quantification of neuronal and non-neuronal nuclei and molecular fingerprinting of cell-type specific transcripts. After clozapine exposure, the proportion of nuclei expressing the neuronal marker NeuN increased by approximately 50% in subcortical white matter, in conjunction with a more subtle and non-significant increase in overlying gray matter. Numbers and proportions of nuclei expressing the oligodendrocyte lineage marker, OLIG2, and cell-type specific RNA expression patterns, were maintained after antipsychotic drug exposure. Frontal lobe gray and white matter volumes remained indistinguishable between antipsychotic-drug-exposed and control groups. Chronic clozapine exposure increases the proportion of NeuN+ nuclei in frontal subcortical white matter, without alterations in frontal lobe volumes or cell type-specific gene expression. Further exploration of neurochemical plasticity in non-human primate brain exposed to antipsychotic drugs is warranted.

Ethical and public policy challenges for pharmacogenomics

It is timely to consider the ethical and social questions raised by progress in pharmacogenomics, based on the current importance of pharmacogenomics for avoidance of predictable side effects of drugs, and for correct choice of medications in certain cancers. It has been proposed that the entire population be genotyped for drug-metabolizing enzyme polymorphisms, as a measure that would prevent many untoward and dangerous drug reactions. Pharmacologic treatment targeting based on genomics of disease can be expected to increase greatly in the coming years. Policy and ethical issues exist on consent for large-scale genomic pharmacogenomic data collection, public vs corporate ownership of genomic research results, testing efficacy and safety of drugs used for rare genomic indications, and accessibility of treatments based on costly research that is applicable to relatively few patients. In major psychiatric disorders and intellectual deficiency, rare and de novo deletion or duplication of chromosomal segments (copy number variation), in the aggregate, are common causes of increased risk. This implies that the policy problems of pharmacogenomics will be particularly important for the psychiatric disorders.

Dual-trait selection for ethanol consumption and withdrawal: genetic and transcriptional network effects

BACKGROUND

Data from C57BL/6J (B6) × DBA/2J (D2) F2 intercrosses (B6xD2 F2 ), standard and recombinant inbred strains, and heterogeneous stock mice indicate that a reciprocal (or inverse) genetic relationship exists between alcohol consumption and withdrawal severity. Furthermore, some genetic studies have detected reciprocal quantitative trait loci (QTLs) for these traits. We used a novel mouse model developed by simultaneous selection for both high alcohol consumption/low withdrawal and low alcohol consumption/high withdrawal and analyzed the gene expression and genome-wide genotypic differences.

METHODS

Randomly chosen third selected generation (S3 ) mice (N = 24/sex/line), bred from a B6xD2 F2 , were genotyped using the Mouse Universal Genotyping Array, which provided 2,760 informative markers. QTL analysis used a marker-by-marker strategy with the threshold for a significant log of the odds (LOD) set at 10. Gene expression in the ventral striatum was measured using the Illumina Mouse 8.2 array. Differential gene expression and the weighted gene co-expression network analysis (WGCNA) were implemented.

RESULTS

Significant QTLs for consumption/withdrawal were detected on chromosomes (Chr) 2, 4, 9, and 12. A suggestive QTL mapped to Chr 6. Some of the QTLs overlapped with known QTLs mapped for 1 of the traits individually. One thousand seven hundred and forty-five transcripts were detected as being differentially expressed between the lines; there was some overlap with known withdrawal genes (e.g., Mpdz) located within QTL regions. WGCNA revealed several modules of co-expressed genes showing significant effects in both differential expression and intramodular connectivity; a module richly annotated with kinase-related annotations was most affected.

CONCLUSIONS

Marked effects of selection on expression and network structure were detected. QTLs overlapping with differentially expressed genes on Chr 2 (distal) and 4 suggest that these are cis-eQTLs (Chr 2: Kif3b, Kcnq2; Chr 4: Mpdz, Snapc3). Other QTLs identified were on Chr 2 (proximal), 9, and 12. Network results point to involvement of kinase-related mechanisms and outline the need for further efforts such as interrogation of noncoding RNAs.

GeneWeaver: finding consilience in heterogeneous cross-species functional genomics data

A persistent challenge lies in the interpretation of consensus and discord from functional genomics experimentation. Harmonizing and analyzing this data will enable investigators to discover relations of many genes to many diseases, and from many phenotypes and experimental paradigms to many diseases through their genomic substrates. The GeneWeaver.org system provides a platform for cross-species integration and interrogation of heterogeneous curated and experimentally derived functional genomics data. GeneWeaver enables researchers to store, share, analyze, and compare results of their own genome-wide functional genomics experiments in an environment containing rich companion data obtained from major curated repositories, including the Mouse Genome Database and other model organism databases, along with derived data from highly specialized resources, publications, and user submissions. The data, largely consisting of gene sets and putative biological networks, are mapped onto one another through gene identifiers and homology across species. A versatile suite of interactive tools enables investigators to perform a variety of set analysis operations to find consilience among these often noisy experimental results. Fast algorithms enable real-time analysis of large queries. Specific applications include prioritizing candidate genes for quantitative trait loci, identifying biologically valid mouse models and phenotypic assays for human disease, finding the common biological substrates of related diseases, classifying experiments and the biological concepts they represent from empirical data, and applying patterns of genomic evidence to implicate novel genes in disease. These results illustrate an alternative to strict emphasis on replicability, whereby researchers classify experimental results to identify the conditions that lead to their similarity.

Cosplicing network analysis of mammalian brain RNA-Seq data utilizing WGCNA and Mantel correlations

Across species and tissues and especially in the mammalian brain, production of gene isoforms is widespread. While gene expression coordination has been previously described as a scale-free coexpression network, the properties of transcriptome-wide isoform production coordination have been less studied. Here we evaluate the system-level properties of cosplicing in mouse, macaque, and human brain gene expression data using a novel network inference procedure. Genes are represented as vectors/lists of exon counts and distance measures sensitive to exon inclusion rates quantifies differences across samples. For all gene pairs, distance matrices are correlated across samples, resulting in cosplicing or cotranscriptional network matrices. We show that networks including cosplicing information are scale-free and distinct from coexpression. In the networks capturing cosplicing we find a set of novel hubs with unique characteristics distinguishing them from coexpression hubs: heavy representation in neurobiological functional pathways, strong overlap with markers of neurons and neuroglia, long coding lengths, and high number of both exons and annotated transcripts. Further, the cosplicing hubs are enriched in genes associated with autism spectrum disorders. Cosplicing hub homologs across eukaryotes show dramatically increasing intronic lengths but stable coding region lengths. Shared transcription factor binding sites increase coexpression but not cosplicing; the reverse is true for splicing-factor binding sites. Genes with protein-protein interactions have strong coexpression and cosplicing. Additional factors affecting the networks include shared microRNA binding sites, spatial colocalization within the striatum, and sharing a chromosomal folding domain. Cosplicing network patterns remain relatively stable across species.

Splicing landscape of the eight collaborative cross founder strains

Background

The Collaborative Cross (CC) is a large panel of genetically diverse recombinant inbred mouse strains specifically designed to provide a systems genetics resource for the study of complex traits. In part, the utility of the CC stems from the extensive genome-wide annotations of founder strain sequence and structural variation. Still missing, however, are transcriptome-specific annotations of the CC founder strains that could further enhance the utility of this resource.

Results

We provide a comprehensive survey of the splicing landscape of the 8 CC founder strains by leveraging the high level of alternative splicing within the brain. Using deep transcriptome sequencing, we found that a majority of the splicing landscape is conserved among the 8 strains, with ~65% of junctions being shared by at least 2 strains. We, however, found a large number of potential strain-specific splicing events as well, with an average of ~3000 and ~500 with ≥3 and ≥10 sequence read coverage, respectively, within each strain. To better understand strain-specific splicing within the CC founder strains, we defined criteria for and identified high-confidence strain-specific splicing events. These splicing events were defined as exon-exon junctions 1) found within only one strain, 2) with a read coverage ≥10, and 3) defined by a canonical splice site. With these criteria, a total of 1509 high-confidence strain-specific splicing events were identified, with the majority found within two of the wild-derived strains, CAST and PWK. Strikingly, the overwhelming majority, 94%, of these strain-specific splicing events are not yet annotated. Strain-specific splicing was also located within genomic regions recently reported to be over- and under-represented within CC populations.

Conclusions

Phenotypic characterization of CC populations is increasing; thus these results will not only aid in further elucidating the transcriptomic architecture of the individual CC founder strains, but they will also help in guiding the utilization of the CC populations in the study of complex traits. This report is also the first to establish guidelines in defining and identifying strain-specific splicing across different mouse strains.

Electronic supplementary material

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

Multi-omic network signatures of disease

To better understand dynamic disease processes, integrated multi-omic methods are needed, yet comparing different types of omic data remains difficult. Integrative solutions benefit experimenters by eliminating potential biases that come with single omic analysis. We have developed the methods needed to explore whether a relationship exists between co-expression network models built from transcriptomic and proteomic data types, and whether this relationship can be used to improve the disease signature discovery process. A naïve, correlation based method is utilized for comparison. Using publicly available infectious disease time series data, we analyzed the related co-expression structure of the transcriptome and proteome in response to SARS-CoV infection in mice. Transcript and peptide expression data was filtered using quality scores and subset by taking the intersection on mapped Entrez IDs. Using this data set, independent co-expression networks were built. The networks were integrated by constructing a bipartite module graph based on module member overlap, module summary correlation, and correlation to phenotypes of interest. Compared to the module level results, the naïve approach is hindered by a lack of correlation across data types, less significant enrichment results, and little functional overlap across data types. Our module graph approach avoids these problems, resulting in an integrated omic signature of disease progression, which allows prioritization across data types for down-stream experiment planning. Integrated modules exhibited related functional enrichments and could suggest novel interactions in response to infection. These disease and platform-independent methods can be used to realize the full potential of multi-omic network signatures. The data (experiment SM001) are publically available through the NIAID Systems Virology (https://www.systemsvirology.org) and PNNL (http://omics.pnl.gov) web portals. Phenotype data is found in the supplementary information. The ProCoNA package is available as part of Bioconductor 2.13.

Coexpression and cosplicing network approaches for the study of mammalian brain transcriptomes

Next-generation sequencing experiments have demonstrated great potential for transcriptome profiling. While transcriptome sequencing greatly increases the level of biological detail, system-level analysis of these high-dimensional datasets is becoming essential. We illustrate gene network approaches to the analysis of transcriptional data, with particular focus on the advantage of RNA-Seq technology compared to microarray platforms. We introduce a novel methodology for constructing cosplicing networks, based on distance measures combined with matrix correlations. We find that the cosplicing network is distinct and complementary to the coexpression network, although it shares the scale-free properties. In the cosplicing network, we find a set of novel hubs that have unique characteristics distinguishing them from coexpression hubs: they are heavily represented in neurobiological functional pathways and have strong overlap with markers of neurons and neuroglia, long-coding lengths, and high number of both exons and annotated transcripts. We also find that gene networks are plastic in the face of genetic and environmental pressures.

Correlated mRNAs and miRNAs from co-expression and regulatory networks affect porcine muscle and finally meat properties

Background

Physiological processes aiding the conversion of muscle to meat involve many genes associated with muscle structure and metabolic processes. MicroRNAs regulate networks of genes to orchestrate cellular functions, in turn regulating phenotypes.

Results

We applied weighted gene co-expression network analysis to identify co-expression modules that correlated to meat quality phenotypes and were highly enriched for genes involved in glucose metabolism, response to wounding, mitochondrial ribosome, mitochondrion, and extracellular matrix. Negative correlation of miRNA with mRNA and target prediction were used to select transcripts out of the modules of trait-associated mRNAs to further identify those genes that are correlated with post mortem traits.

Conclusions

Porcine muscle co-expression transcript networks that correlated to post mortem traits were identified. The integration of miRNA and mRNA expression analyses, as well as network analysis, enabled us to interpret the differentially-regulated genes from a systems perspective. Linking co-expression networks of transcripts and hierarchically organized pairs of miRNAs and mRNAs to meat properties yields new insight into several biological pathways underlying phenotype differences. These pathways may also be diagnostic for many myopathies, which are accompanied by deficient nutrient and oxygen supply of muscle fibers.

Opportunities for bioinformatics in the classification of behavior and psychiatric disorders

A bioinformatics approach to behavioral neuroscience provides both unique opportunities and challenges for research on behavior. A major challenge has been to describe, define, and discriminate among abstract behavioral processes, in large part by distinguishing among the biological mechanisms of unique but not entirely discrete, entities of behavior. Understanding the complexity of neurobiology and behavior requires integration of data across diverse biological systems, types of data, and levels of scale. With the perspective and application of bioinformatics, we can uncover the relationships among these systems and take steps forward in realizing the common and distinct bases of psychiatric disease.

Selection for drinking in the dark alters brain gene coexpression networks

BACKGROUND

Heterogeneous stock (HS/NPT) mice have been used to create lines selectively bred in replicate for elevated drinking in the dark (DID). Both selected lines routinely reach a blood ethanol (EtOH) concentration (BEC) of 1.00 mg/ml or greater at the end of the 4-hour period of access in Day 2. The mechanisms through which genetic differences influence DID are currently unclear. Therefore, the current study examines the transcriptome, the first stage at which genetic variability affects neurobiology. Rather than focusing solely on differential expression (DE), we also examine changes in the ways that gene transcripts collectively interact with each other, as revealed by changes in coexpression patterns.

METHODS

Naïve mice (N = 48/group) were genotyped using the Mouse Universal Genotyping Array, which provided 3,683 informative markers. Quantitative trait locus (QTL) analysis used a marker-by-marker strategy with the threshold for a significant logarithm of odds (LOD) set at 10.6. Gene expression in the ventral striatum was measured using the Illumina Mouse 8.2 array. Differential gene expression and the weighted gene coexpression network analysis (WGCNA) were implemented largely as described elsewhere.

RESULTS

Significant QTLs for elevated BECs after DID were detected on chromosomes 4, 14, and 16; the latter 2 were associated with gene-poor regions. None of the QTLs overlapped with known QTLs for EtOH preference drinking. Ninety-four transcripts were detected as being differentially expressed in both selected lines versus HS controls; there was no overlap with known preference genes. The WGCNA revealed 2 modules as showing significant effects of both selections on intramodular connectivity. A number of genes known to be associated with EtOH phenotypes (e.g., Gabrg1, Glra2, Grik1, Npy2r, and Nts) showed significant changes in connectivity.

CONCLUSIONS

We found marked and consistent effects of selection on coexpression patterns; DE changes were more modest and less concordant. The QTLs and differentially expressed genes detected here are distinct from the preference phenotype. This is consistent with behavioral data and suggests that the DID and preference phenotypes are markedly different genetically.

Differential network analysis reveals genetic effects on catalepsy modules

We performed short-term bi-directional selective breeding for haloperidol-induced catalepsy, starting from three mouse populations of increasingly complex genetic structure: an F₂ intercross, a heterogeneous stock (HS) formed by crossing four inbred strains (HS4) and a heterogeneous stock (HS-CC) formed from the inbred strain founders of the Collaborative Cross (CC). All three selections were successful, with large differences in haloperidol response emerging within three generations. Using a custom differential network analysis procedure, we found that gene coexpression patterns changed significantly; importantly, a number of these changes were concordant across genetic backgrounds. In contrast, absolute gene-expression changes were modest and not concordant across genetic backgrounds, in spite of the large and similar phenotypic differences. By inferring strain contributions from the parental lines, we are able to identify significant differences in allelic content between the selected lines concurrent with large changes in transcript connectivity. Importantly, this observation implies that genetic polymorphisms can affect transcript and module connectivity without large changes in absolute expression levels. We conclude that, in this case, selective breeding acts at the subnetwork level, with the same modules but not the same transcripts affected across the three selections.

Genes, behavior and next-generation RNA sequencing

Advances in next-generation sequencing suggest that RNA-Seq is poised to supplant microarray-based approaches for transcriptome analysis. This article briefly reviews the use of microarrays in the brain-behavior context and then illustrates why RNA-Seq is a superior strategy. Compared with microarrays, RNA-Seq has a greater dynamic range, detects both coding and noncoding RNAs, is superior for gene network construction, detects alternative spliced transcripts, detects allele specific expression and can be used to extract genotype information, e.g. nonsynonymous coding single nucleotide polymorphisms. Examples of where RNA-Seq has been used to assess brain gene expression are provided. Despite the advantages of RNA-Seq, some disadvantages remain. These include the high cost of RNA-Seq and the computational complexities associated with data analysis. RNA-Seq embraces the complexity of the transcriptome and provides a mechanism to understand the underlying regulatory code; the potential to inform the brain-behavior relationship is substantial.

Utilizing RNA-Seq data for de novo coexpression network inference

MOTIVATION

RNA-Seq experiments have shown great potential for transcriptome profiling. While sequencing increases the level of biological detail, integrative data analysis is also important. One avenue is the construction of coexpression networks. Because the capacity of RNA-Seq data for network construction has not been previously evaluated, we constructed a coexpression network using striatal samples, derived its network properties and compared it with microarray-based networks.

RESULTS

The RNA-Seq coexpression network displayed scale-free, hierarchical network structure. We detected transcripts groups (modules) with correlated profiles; modules overlap distinct ontology categories. Neuroanatomical data from the Allen Brain Atlas reveal several modules with spatial colocalization. The network was compared with microarray-derived networks; correlations from RNA-Seq data were higher, likely because greater sensitivity and dynamic range. Higher correlations result in higher network connectivity, heterogeneity and centrality. For transcripts present across platforms, network structure appeared largely preserved. From this study, we present the first RNA-Seq data de novo network inference.