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

Transcriptomic changes in the posterior pallium of male zebra finches associated with social niche conformance

Animals plastically adjust their physiological and behavioural phenotypes to conform to their social environment-social niche conformance. The degree of sexual competition is a critical part of the social environment to which animals adjust their phenotypes, but the underlying genetic mechanisms are poorly understood. We conducted a study to investigate how differences in sperm competition risk affect the gene expression profiles of the testes and two brain areas (posterior pallium and optic tectum) in breeding male zebra finches (Taeniopygia castanotis). In this pre-registered study, we investigated a large sample of 59 individual transcriptomes. We compared two experimental groups: males held in single breeding pairs (low sexual competition) versus those held in two pairs (elevated sexual competition) per breeding cage. Using weighted gene co-expression network analysis (WGCNA), we observed significant effects of the social treatment in all three tissues. However, only the treatment effects found in the pallium were confirmed by an additional randomisation test for statistical robustness. Likewise, the differential gene expression analysis revealed treatment effects only in the posterior pallium (ten genes) and optic tectum (six genes). No treatment effects were found in the testis at the single gene level. Thus, our experiments do not provide strong evidence for transcriptomic adjustment specific to manipulated sperm competition risk. However, we did observe transcriptomic adjustments to the manipulated social environment in the posterior pallium. These effects were polygenic rather than based on few individual genes with strong effects. Our findings are discussed in relation to an accompanying paper using the same animals, which reports behavioural results consistent with the results presented here.

Functional convergence of genomic and transcriptomic architecture underlies schooling behaviour in a live-bearing fish

The organization and coordination of fish schools provide a valuable model to investigate the genetic architecture of affiliative behaviours and dissect the mechanisms underlying social behaviours and personalities. Here we used replicate guppy selection lines that vary in schooling propensity and combine quantitative genetics with genomic and transcriptomic analyses to investigate the genetic basis of sociability phenotypes. We show that consistent with findings in collective motion patterns, experimental evolution of schooling propensity increased the sociability of female, but not male, guppies when swimming with unfamiliar conspecifics. This finding highlights a relevant link between coordinated motion and sociability for species forming fission-fusion societies in which both group size and the type of social interactions are dynamic across space and time. We further show that alignment and attraction, the two major traits forming the sociability personality axis in this species, showed heritability estimates at the upper end of the range previously described for social behaviours, with important variation across sexes. The results from both Pool-seq and RNA-seq data indicated that genes involved in neuron migration and synaptic function were instrumental in the evolution of sociability, highlighting a crucial role of glutamatergic synaptic function and calcium-dependent signalling processes in the evolution of schooling.

Metabolic changes in the brain and blood of rats following acoustic trauma, tinnitus and hyperacusis

It has been increasingly recognized that tinnitus is likely to be generated by complex network changes. Acoustic trauma that causes tinnitus induces significant changes in multiple metabolic pathways in the brain. However, it is not clear whether those metabolic changes in the brain could also be reflected in blood samples and whether metabolic changes could discriminate acoustic trauma, hyperacusis and tinnitus. We analyzed brain and serum metabolic changes in rats following acoustic trauma or a sham procedure using metabolomics. Hearing levels were recorded before and after acoustic trauma and behavioral measures to quantify tinnitus and hyperacusis were conducted at 4 weeks following acoustic trauma. Tissues from 11 different brain regions and serum samples were collected at about 3 months following acoustic trauma. Among the acoustic trauma animals, eight exhibited hyperacusis-like behavior and three exhibited tinnitus-like behavior. Using Gas chromatography-mass spectrometry and multivariate statistical analysis, significant metabolic changes were found in acoustic trauma animals in both the brain and serum samples with a number of metabolic pathways significantly perturbated. Furthermore, metabolic changes in the serum were able to differentiate sham from acoustic trauma animals, as well as sham from hyperacusis animals, with high accuracy. Our results suggest that serum metabolic profiling in combination with machine learning analysis may be a promising approach for identifying biomarkers for acoustic trauma, hyperacusis and potentially, tinnitus.

Pharmacological Evaluation of Drugs in Animal Models of Tinnitus

Despite the pressing need for effective drug treatments for tinnitus, currently, there is no single drug that is approved by the FDA for this purpose. Instead, a wide range of unproven over-the-counter tinnitus remedies are available on the market with little or no benefit for tinnitus but with potential harm and adverse effects. Animal models of tinnitus have played a critical role in exploring the pathophysiology of tinnitus, identifying therapeutic targets and evaluating novel and existing drugs for tinnitus treatment. This review summarises and compares the studies on pharmacological evaluation of tinnitus treatment in different animal models based on the pharmacological properties of the drug and provides insights into future directions for tinnitus drug discovery.

Early neurogenomic response associated with variation in guppy female mate preference

Understanding the evolution of mate choice requires dissecting the mechanisms of female preference, particularly how these differ among social contexts and preference phenotypes. Here we study the female neurogenomic response after only 10 minutes of mate exposure in both a sensory component (optic tectum) and a decision-making component (telencephalon) of the brain. By comparing the transcriptional response between females with and without preferences for colorful males, we identified unique neurogenomic elements associated with the female preference phenotype that are not present in females without preference. Network analysis revealed different properties for this response at the sensory-processing and the decision-making levels, and showed that this response is highly centralized in the telencephalon. Furthermore, we identified an additional set of genes that vary in expression across social contexts, beyond mate evaluation. We show that transcription factors among those loci are predicted to regulate the transcriptional response of the genes we found to be associated with female preference.

Progress With Nonhuman Animal Models of Addiction

Nonhuman animals have been major contributors to the science of the genetics of addiction. Given the explosion of interest in genetics, it is fair to ask, are we making reasonable progress toward our goals with animal models? I will argue that our goals are changing and that overall progress has been steady and seems likely to continue apace. Genetics tools have developed almost incredibly rapidly, enabling both more reductionist and more synthetic or integrative approaches. I believe that these approaches to making progress have been unbalanced in biomedical science, favoring reductionism, particularly in animal genetics. I argue that substantial, novel progress is also likely to come in the other direction, toward synthesis and abstraction. Another area in which future progress with genetic animal models seems poised to contribute more is the reconciliation of human and animal phenotypes, or consilience. The inherent power of the genetic animal models could be more profitably exploited. In the end, animal research has continued to provide novel insights about how genes influence individual differences in addiction risk and consequences. The rules of the genetics game are changing so fast that it is hard to remember how comparatively little we knew even a generation ago. Rather than worry about whether we have been wasting time and resources asking the questions we have been, we should look to the future and see if we can come up with some new ones. The valuable findings from the past will endure, and the sidetracks will be forgotten.

Brain Metabolic Changes in Rats following Acoustic Trauma

Acoustic trauma is the most common cause of hearing loss and tinnitus in humans. However, the impact of acoustic trauma on system biology is not fully understood. It has been increasingly recognized that tinnitus caused by acoustic trauma is unlikely to be generated by a single pathological source, but rather a complex network of changes involving not only the auditory system but also systems related to memory, emotion and stress. One obvious and significant gap in tinnitus research is a lack of biomarkers that reflect the consequences of this interactive "tinnitus-causing" network. In this study, we made the first attempt to analyse brain metabolic changes in rats following acoustic trauma using metabolomics, as a pilot study prior to directly linking metabolic changes to tinnitus. Metabolites in 12 different brain regions collected from either sham or acoustic trauma animals were profiled using a gas chromatography mass spectrometry (GC/MS)-based metabolomics platform. After deconvolution of mass spectra and identification of the molecules, the metabolomic data were processed using multivariate statistical analysis. Principal component analysis showed that metabolic patterns varied among different brain regions; however, brain regions with similar functions had a similar metabolite composition. Acoustic trauma did not change the metabolite clusters in these regions. When analyzed within each brain region using the orthogonal projection to latent structures discriminant analysis sub-model, 17 molecules showed distinct separation between control and acoustic trauma groups in the auditory cortex, inferior colliculus, superior colliculus, vestibular nucleus complex (VNC), and cerebellum. Further metabolic pathway impact analysis and the enrichment overview with network analysis suggested the primary involvement of amino acid metabolism, including the alanine, aspartate and glutamate metabolic pathways, the arginine and proline metabolic pathways and the purine metabolic pathway. Our results provide the first metabolomics evidence that acoustic trauma can induce changes in multiple metabolic pathways. This pilot study also suggests that the metabolomic approach has the potential to identify acoustic trauma-specific metabolic shifts in future studies where metabolic changes are correlated with the animal's tinnitus status.

CoSpliceNet: a framework for co-splicing network inference from transcriptomics data

Background

Alternative splicing has been proposed to increase transcript diversity and protein plasticity in eukaryotic organisms, but the extent to which this is the case is currently unclear, especially with regard to the diversification of molecular function. Eukaryotic splicing involves complex interactions of splicing factors and their targets. Inference of co-splicing networks capturing these types of interactions is important for understanding this crucial, highly regulated post-transcriptional process at the systems level.

Results

First, several transcript and protein attributes, including coding potential of transcripts and differences in functional domains of proteins, were compared between splice variants and protein isoforms to assess transcript and protein diversity in a biological system. Alternative splicing was shown to increase transcript and function-related protein diversity in developing Arabidopsis embryos. Second, CoSpliceNet, which integrates co-expression and motif discovery at splicing regulatory regions to infer co-splicing networks, was developed. CoSpliceNet was applied to temporal RNA sequencing data to identify candidate regulators of splicing events and predict RNA-binding motifs, some of which are supported by prior experimental evidence. Analysis of inferred splicing factor targets revealed an unexpected role for the unfolded protein response in embryo development.

Conclusions

The methods presented here can be used in any biological system to assess transcript diversity and protein plasticity and to predict candidate regulators, their targets, and RNA-binding motifs for splicing factors. CoSpliceNet is freely available at http://delasa.github.io/co-spliceNet/.

Electronic supplementary material

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

An Integrative Genomic Study Implicates the Postsynaptic Density in the Pathogenesis of Bipolar Disorder

Genome-wide association studies (GWAS) have identified several common variants associated with bipolar disorder (BD), but the biological meaning of these findings remains unclear. Integrative genomics-the integration of GWAS signals with gene expression data-may illuminate genes and gene networks that have key roles in the pathogenesis of BD. We applied weighted gene co-expression network analysis (WGCNA), which exploits patterns of co-expression among genes, to brain transcriptome data obtained by sequencing of poly-A RNA derived from postmortem dorsolateral prefrontal cortex from people with BD, along with age- and sex-matched controls. WGCNA identified 33 gene modules. Many of the modules corresponded closely to those previously reported in human cortex. Three modules were associated with BD, enriched for genes differentially expressed in BD, and also enriched for signals in prior GWAS of BD. Functional analysis of genes within these modules revealed significant enrichment of several functionally related sets of genes, especially those involved in the postsynaptic density (PSD). These results provide convergent support for the hypothesis that dysregulation of genes involved in the PSD is a key factor in the pathogenesis of BD. If replicated in larger samples, these findings could point toward new therapeutic targets for BD.

Advanced Applications of RNA Sequencing and Challenges

Next-generation sequencing technologies have revolutionarily advanced sequence-based research with the advantages of high-throughput, high-sensitivity, and high-speed. RNA-seq is now being used widely for uncovering multiple facets of transcriptome to facilitate the biological applications. However, the large-scale data analyses associated with RNA-seq harbors challenges. In this study, we present a detailed overview of the applications of this technology and the challenges that need to be addressed, including data preprocessing, differential gene expression analysis, alternative splicing analysis, variants detection and allele-specific expression, pathway analysis, co-expression network analysis, and applications combining various experimental procedures beyond the achievements that have been made. Specifically, we discuss essential principles of computational methods that are required to meet the key challenges of the RNA-seq data analyses, development of various bioinformatics tools, challenges associated with the RNA-seq applications, and examples that represent the advances made so far in the characterization of the transcriptome.

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.