Integrating many co-splicing networks to reconstruct splicing regulatory modules

(Differential) Co-Expression Analysis of Gene Expression: A Survey of Best Practices

Analysis of gene expression data is widely used in transcriptomic studies to understand functions of molecules inside a cell and interactions among molecules. Differential co-expression analysis studies diseases and phenotypic variations by finding modules of genes whose co-expression patterns vary across conditions. We review the best practices in gene expression data analysis in terms of analysis of (differential) co-expression, co-expression network, differential networking, and differential connectivity considering both microarray and RNA-seq data along with comparisons. We highlight hurdles in RNA-seq data analysis using methods developed for microarrays. We include discussion of necessary tools for gene expression analysis throughout the paper. In addition, we shed light on scRNA-seq data analysis by including preprocessing and scRNA-seq in co-expression analysis along with useful tools specific to scRNA-seq. To get insights, biological interpretation and functional profiling is included. Finally, we provide guidelines for the analyst, along with research issues and challenges which should be addressed.

Transcriptional Networks of Microglia in Alzheimer's Disease and Insights into Pathogenesis

Microglia, the main immune cells of the central nervous system, are increasingly implicated in Alzheimer's disease (AD). Manifold transcriptomic studies in the brain have not only highlighted microglia's role in AD pathogenesis, but also mapped crucial pathological processes and identified new therapeutic targets. An important component of many of these transcriptomic studies is the investigation of gene expression networks in AD brain, which has provided important new insights into how coordinated gene regulatory programs in microglia (and other cell types) underlie AD pathogenesis. Given the rapid technological advancements in transcriptional profiling, spanning from microarrays to single-cell RNA sequencing (scRNA-seq), tools used for mapping gene expression networks have evolved to keep pace with the unique features of each transcriptomic platform. In this article, we review the trajectory of transcriptomic network analyses in AD from brain to microglia, highlighting the corresponding methodological developments. Lastly, we discuss examples of how transcriptional network analysis provides new insights into AD mechanisms and pathogenesis.

Bioinformatics challenges and perspectives when studying the effect of epigenetic modifications on alternative splicing

It is widely known that epigenetic modifications are important in regulating transcription, but several have also been reported in alternative splicing. The regulation of pre-mRNA splicing is important to explain proteomic diversity and the misregulation of splicing has been implicated in many diseases. Here, we give a brief overview of the role of epigenetics in alternative splicing and disease. We then discuss the bioinformatics methods that can be used to model interactions between epigenetic marks and regulators of splicing. These models can be used to identify alternative splicing and epigenetic changes across different phenotypes.

This article is part of a discussion meeting issue ‘Frontiers in epigenetic chemical biology’.

Gene co-expression analysis for functional classification and gene-disease predictions

Gene co-expression networks can be used to associate genes of unknown function with biological processes, to prioritize candidate disease genes or to discern transcriptional regulatory programmes. With recent advances in transcriptomics and next-generation sequencing, co-expression networks constructed from RNA sequencing data also enable the inference of functions and disease associations for non-coding genes and splice variants. Although gene co-expression networks typically do not provide information about causality, emerging methods for differential co-expression analysis are enabling the identification of regulatory genes underlying various phenotypes. Here, we introduce and guide researchers through a (differential) co-expression analysis. We provide an overview of methods and tools used to create and analyse co-expression networks constructed from gene expression data, and we explain how these can be used to identify genes with a regulatory role in disease. Furthermore, we discuss the integration of other data types with co-expression networks and offer future perspectives of co-expression analysis.

Co-expression networks reveal the tissue-specific regulation of transcription and splicing

Gene co-expression networks capture biologically important patterns in gene expression data, enabling functional analyses of genes, discovery of biomarkers, and interpretation of genetic variants. Most network analyses to date have been limited to assessing correlation between total gene expression levels in a single tissue or small sets of tissues. Here, we built networks that additionally capture the regulation of relative isoform abundance and splicing, along with tissue-specific connections unique to each of a diverse set of tissues. We used the Genotype-Tissue Expression (GTEx) project v6 RNA sequencing data across 50 tissues and 449 individuals. First, we developed a framework called Transcriptome-Wide Networks (TWNs) for combining total expression and relative isoform levels into a single sparse network, capturing the interplay between the regulation of splicing and transcription. We built TWNs for 16 tissues and found that hubs in these networks were strongly enriched for splicing and RNA binding genes, demonstrating their utility in unraveling regulation of splicing in the human transcriptome. Next, we used a Bayesian biclustering model that identifies network edges unique to a single tissue to reconstruct Tissue-Specific Networks (TSNs) for 26 distinct tissues and 10 groups of related tissues. Finally, we found genetic variants associated with pairs of adjacent nodes in our networks, supporting the estimated network structures and identifying 20 genetic variants with distant regulatory impact on transcription and splicing. Our networks provide an improved understanding of the complex relationships of the human transcriptome across tissues.

The determinants of alternative RNA splicing in human cells

Alternative splicing represents an important level of the regulation of gene function in eukaryotic organisms. It plays a critical role in virtually every biological process within an organism, including regulation of cell division and cell death, differentiation of tissues in the embryo and the adult organism, as well as in cellular response to diverse environmental factors. In turn, studies of the last decade have shown that alternative splicing itself is controlled by different mechanisms. Unfortunately, there is no clear understanding of how these diverse mechanisms, or determinants, regulate and constrain the set of alternative RNA species produced from any particular gene in every cell of the human body. Here, we provide a consolidated overview of alternative splicing determinants including RNA-protein interactions, epigenetic regulation via chromatin remodeling, coupling of transcription-to-alternative splicing, effect of secondary structures in pre-RNA, and function of the RNA quality control systems. We also extensively and critically discuss some mechanistic insights on coordinated inclusion/exclusion of exons during the formation of mature RNA molecules. We conclude that the final structure of RNA is pre-determined by a complex interplay between cis- and trans-acting factors. Altogether, currently available empirical data significantly expand our understanding of the functioning of the alternative splicing machinery of cells in normal and pathological conditions. On the other hand, there are still many blind spots that require further deep investigations.

Next generation sequencing technology and genomewide data analysis: Perspectives for retinal research

The advent of high throughput next generation sequencing (NGS) has accelerated the pace of discovery of disease-associated genetic variants and genomewide profiling of expressed sequences and epigenetic marks, thereby permitting systems-based analyses of ocular development and disease. Rapid evolution of NGS and associated methodologies presents significant challenges in acquisition, management, and analysis of large data sets and for extracting biologically or clinically relevant information. Here we illustrate the basic design of commonly used NGS-based methods, specifically whole exome sequencing, transcriptome, and epigenome profiling, and provide recommendations for data analyses. We briefly discuss systems biology approaches for integrating multiple data sets to elucidate gene regulatory or disease networks. While we provide examples from the retina, the NGS guidelines reviewed here are applicable to other tissues/cell types as well.

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.

Integromics network meta-analysis on cardiac aging offers robust multi-layer modular signatures and reveals micronome synergism

Background

The avalanche of integromics and panomics approaches shifted the deciphering of aging mechanisms from single molecular entities to communities of them. In this orientation, we explore the cardiac aging mechanisms – risk factor for multiple cardiovascular diseases - by capturing the micronome synergism and detecting longevity signatures in the form of communities (modules).

For this, we developed a meta-analysis scheme that integrates transcriptome expression data from multiple cardiac-specific independent studies in mouse and human along with proteome and micronome interaction data in the form of multiple independent weighted networks. Modularization of each weighted network produced modules, which in turn were further analyzed so as to define consensus modules across datasets that change substantially during lifespan. Also, we established a metric that determines - from the modular perspective - the synergism of microRNA-microRNA interactions as defined by significantly functionally associated targets.

Results

The meta-analysis provided 40 consensus integromics modules across mouse datasets and revealed microRNA relations with substantial collective action during aging. Three modules were reproducible, based on homology, when mapped against human-derived modules. The respective homologs mainly represent NADH dehydrogenases, ATP synthases, cytochrome oxidases, Ras GTPases and ribosomal proteins. Among various observations, we corroborate to the involvement of miR-34a (included in consensus modules) as proposed recently; yet we report that has no synergistic effect. Moving forward, we determined its age-related neighborhood in which HCN3, a known heart pacemaker channel, was included. Also, miR-125a-5p/-351, miR-200c/-429, miR-106b/-17, miR-363/-92b, miR-181b/-181d, miR-19a/-19b, let-7d/-7f, miR-18a/-18b, miR-128/-27b and miR-106a/-291a-3p pairs exhibited significant synergy and their association to aging and/or cardiovascular diseases is supported in many cases by a disease database and previous studies. On the contrary, we suggest that miR-22 has not substantial impact on heart longevity as proposed recently.

Conclusions

We revised several proteins and microRNAs recently implicated in cardiac aging and proposed for the first time modules as signatures. The integromics meta-analysis approach can serve as an efficient subvening signature tool for more-oriented better-designed experiments. It can also promote the combinational multi-target microRNA therapy of age-related cardiovascular diseases along the continuum from prevention to detection, diagnosis, treatment and outcome.

Electronic supplementary material

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

Integrative analysis of many RNA-seq datasets to study alternative splicing

Alternative splicing is an important gene regulatory mechanism that dramatically increases the complexity of the proteome. However, how alternative splicing is regulated and how transcription and splicing are coordinated are still poorly understood, and functions of transcript isoforms have been studied only in a few limited cases. Nowadays, RNA-seq technology provides an exceptional opportunity to study alternative splicing on genome-wide scales and in an unbiased manner. With the rapid accumulation of data in public repositories, new challenges arise from the urgent need to effectively integrate many different RNA-seq datasets for study alterative splicing. This paper discusses a set of advanced computational methods that can integrate and analyze many RNA-seq datasets to systematically identify splicing modules, unravel the coupling of transcription and splicing, and predict the functions of splicing isoforms on a genome-wide scale.

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.

Introduction to sequencing the brain transcriptome

High-throughput next-generation sequencing is now entering its second decade. However, it was not until 2008 that the first report of sequencing the brain transcriptome appeared (Mortazavi, Williams, Mccue, Schaeffer, & Wold, 2008). These authors compared short-read RNA-Seq data for mouse whole brain with microarray results for the same sample and noted both the advantages and disadvantages of the RNA-Seq approach. While RNA-Seq provided exon level resolution, the majority of the reads were provided by a small proportion of highly expressed genes and the data analysis was exceedingly complex. Over the past 6 years, there have been substantial improvements in both RNA-Seq technology and data analysis. This volume contains 11 chapters that detail various aspects of sequencing the brain transcriptome. Some of the chapters are very methods driven, while others focus on the use of RNA-Seq to study such diverse areas as development, schizophrenia, and drug abuse. This chapter briefly reviews the transition from microarrays to RNA-Seq as the preferred method for analyzing the brain transcriptome. 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, and can be used to extract genotype information, e.g., nonsynonymous coding single nucleotide polymorphisms. RNA-Seq embraces the complexity of the brain transcriptome and provides a mechanism to understand the underlying regulatory code; the potential to inform the brain-behavior-disease relationships is substantial.