Search

RNA editing, the endogenous modification of nucleic acids, is known to be altered in genes with important neurological function in schizophrenia (SCZ). However, the global profile and molecular functions of disease-associated RNA editing remain unclear. Here, we analyzed RNA editing in postmortem brains of four SCZ cohorts and uncovered a significant and reproducible trend of hypoediting in patients of European descent. We report a set of SCZ-associated editing sites via WGCNA analysis, shared across cohorts. Using massively parallel reporter assays and bioinformatic analyses, we observed that differential 3′ untranslated region (3′UTR) editing sites affecting host gene expression were enriched for mitochondrial processes. Furthermore, we characterized the impact of two recoding sites in the mitofusin 1 ( MFN1 ) gene and showed their functional relevance to mitochondrial fusion and cellular apoptosis. Our study reveals a global reduction of editing in SCZ and a compelling link between editing and mitochondrial function in the disease.

Supplementary Figures S1-11 and Tables S1-2 are included in this file

RNA editing modifies single nucleotides of RNAs, regulating primary protein structure and protein abundance. In recent years, the diversity of proteins and complexity of gene regulation associated with RNA editing dysregulation has been increasingly appreciated in oncology. Large-scale shifts in editing have been observed in bulk tumors across various cancer types. However, RNA editing in single cells and individual cell types within tumors has not been explored. By profiling editing in single cells from lung adenocarcinoma biopsies, we found that the increased editing trend of bulk lung tumors was unique to cancer cells. Elevated editing levels were observed in cancer cells resistant to targeted therapy, and editing sites associated with drug response were enriched. Consistent with the regulation of antiviral pathways by RNA editing, higher editing levels in cancer cells were associated with reduced antitumor innate immune response, especially levels of natural killer cell infiltration. In addition, the level of RNA editing in cancer cells was positively associated with somatic point mutation burden. This observation motivated the definition of a new metric, RNA editing load, reflecting the amount of RNA mutations created by RNA editing. Importantly, in lung cancer, RNA editing load was a stronger predictor of patient survival than DNA mutations. This study provides the first single cell dissection of editing in cancer and highlights the significance of RNA editing load in cancer prognosis.Significance:RNA editing analysis in single lung adenocarcinoma cells uncovers RNA mutations that correlate with tumor mutation burden and cancer innate immunity and reveals the amount of RNA mutations that strongly predicts patient survival.See related commentary by Luo and Liang, p. 351

RNA editing modifies single nucleotides of RNAs, regulating primary protein structure and protein abundance. In recent years, the diversity of proteins and complexity of gene regulation associated with RNA editing dysregulation has been increasingly appreciated in oncology. Large-scale shifts in editing have been observed in bulk tumors across various cancer types. However, RNA editing in single cells and individual cell types within tumors has not been explored. By profiling editing in single cells from lung adenocarcinoma biopsies, we found that the increased editing trend of bulk lung tumors was unique to cancer cells. Elevated editing levels were observed in cancer cells resistant to targeted therapy, and editing sites associated with drug response were enriched. Consistent with the regulation of antiviral pathways by RNA editing, higher editing levels in cancer cells were associated with reduced antitumor innate immune response, especially levels of natural killer cell infiltration. In addition, the level of RNA editing in cancer cells was positively associated with somatic point mutation burden. This observation motivated the definition of a new metric, RNA editing load, reflecting the amount of RNA mutations created by RNA editing. Importantly, in lung cancer, RNA editing load was a stronger predictor of patient survival than DNA mutations. This study provides the first single cell dissection of editing in cancer and highlights the significance of RNA editing load in cancer prognosis. Significance: RNA editing analysis in single lung adenocarcinoma cells uncovers RNA mutations that correlate with tumor mutation burden and cancer innate immunity and reveals the amount of RNA mutations that strongly predicts patient survival. See related commentary by Luo and Liang, p. 351

To identify genetic pathways for addiction, we analyzed intravenous self-administration of cocaine or saline in a panel of 84 inbred and recombinant inbred mouse strains over 10 days. We integrated the behavior data with RNA-Seq data from the medial frontal cortex and nucleus accumbens from 41 strains. The self-administration of cocaine and saline showed distinct genetic bases. We maximized power to map loci for cocaine intake by using a linear mixed model to account for this longitudinal phenotype while correcting for population structure. A total of 15 unique significant loci were identified in the genome-wide association study (GWAS). A transcriptome-wide association study (TWAS) highlighted theTrpv2ion channel as a key locus for cocaine self-administration from the GWAS. In addition, 17 genes supplementary to the GWAS were identified includingArhgef26, Slc18b1andSlco5a1. We found numerous instances where alternate splice site selection or RNA editing altered transcript abundance. Our work emphasizes the importance ofTrpv2, a known cannabinoid receptor, for the response to cocaine as well as identifying further relevant loci.

The zebra finch (Taeniopygia guttata), a representative oscine songbird species, has been widely studied to investigate behavioral neuroscience, most notably the neurobiological basis of vocal learning, a rare trait shared in only a few animal groups including humans. In 2019, an updated zebra finch genome annotation (bTaeGut1_v1.p) was released from the Ensembl database and is substantially more comprehensive than the first version published in 2010. In this study, we utilized the publicly available RNA-seq data generated from Illumina-based short-reads and PacBio single-molecule real-time (SMRT) long-reads to assess the bird transcriptome. To analyze the high-throughput RNA-seq data, we adopted a hybrid bioinformatic approach combining short and long-read pipelines. From our analysis, we added 220 novel genes and 8,134 transcript variants to the Ensembl annotation, and predicted a new proteome based on the refined annotation. We further validated 18 different novel proteins by using mass-spectrometry data generated from zebra finch caudal telencephalon tissue. Our results provide additional resources for future studies of zebra finches utilizing this improved bird genome annotation and proteome.

PODXL, a protein that is dysregulated in multiple cancers, plays an important role in promoting cancer metastasis. In this study, we report that RNA editing promotes the inclusion of a PODXL alternative exon. The resulting edited PODXL long isoform is more prone to protease digestion and has the strongest effects on reducing cell migration and cisplatin chemoresistance among the three PODXL isoforms (short, unedited long, and edited long isoforms). Importantly, the editing level of the PODXL recoding site and the inclusion level of the PODXL alternative exon are strongly associated with overall patient survival in Kidney Renal Clear Cell Carcinoma (KIRC). Supported by significant enrichment of exonic RNA editing sites in alternatively spliced exons, we hypothesize that exonic RNA editing sites may enhance proteomic diversity through alternative splicing, in addition to amino acid changes, a previously under-appreciated aspect of RNA editing function.

Human stem-cell-derived models provide the promise of accelerating our understanding of brain disorders, but not knowing whether they possess the ability to mature beyond mid- to late-fetal stages potentially limits their utility. We leveraged a directed differentiation protocol to comprehensively assess maturation in vitro. Based on genome-wide analysis of the epigenetic clock and transcriptomics, as well as RNA editing, we observe that three-dimensional human cortical organoids reach postnatal stages between 250 and 300 days, a timeline paralleling in vivo development. We demonstrate the presence of several known developmental milestones, including switches in the histone deacetylase complex and NMDA receptor subunits, which we confirm at the protein and physiological levels. These results suggest that important components of an intrinsic in vivo developmental program persist in vitro. We further map neurodevelopmental and neurodegenerative disease risk genes onto in vitro gene expression trajectories to provide a resource and webtool (Gene Expression in Cortical Organoids, GECO) to guide disease modeling.

Alternative splicing is an RNA processing mechanism that affects most genes in human, contributing to disease mechanisms and phenotypic diversity. The regulation of splicing involves an intricate network of cis-regulatory elements and trans-acting factors. Due to their high sequence specificity, cis-regulation of splicing can be altered by genetic variants, significantly affecting splicing outcomes. Recently, multiple methods have been applied to understanding the regulatory effects of genetic variants on splicing. However, it is still challenging to go beyond apparent association to pinpoint functional variants. To fill in this gap, we utilized large-scale data sets of the Genotype-Tissue Expression (GTEx) project to study genetically modulated alternative splicing (GMAS) via identification of allele-specific splicing events. We demonstrate that GMAS events are shared across tissues and individuals more often than expected by chance, consistent with their genetically driven nature. Moreover, although the allelic bias of GMAS exons varies across samples, the degree of variation is similar across tissues versus individuals. Thus, genetic background drives the GMAS pattern to a similar degree as tissue-specific splicing mechanisms. Leveraging the genetically driven nature of GMAS, we developed a new method to predict functional splicing-altering variants, built upon a genotype-phenotype concordance model across samples. Complemented by experimental validations, this method predicted >1000 functional variants, many of which may alter RNA-protein interactions. Lastly, 72% of GMAS-associated SNPs were in linkage disequilibrium with GWAS-reported SNPs, and such association was enriched in tissues of relevance for specific traits/diseases. Our study enables a comprehensive view of genetically driven splicing variations in human tissues.

Single-cell RNA sequencing (scRNA-seq) data contain rich information at the gene, transcript, and nucleotide levels. Most analyses of scRNA-seq have focused on gene expression profiles, and it remains challenging to extract nucleotide variants and isoform-specific information. Here, we present scAllele, an integrative approach that detects single-nucleotide variants, insertions, deletions, and their allelic linkage with splicing patterns in scRNA-seq. We demonstrate that scAllele achieves better performance in identifying nucleotide variants than other commonly used tools. In addition, the read-specific variant calls by scAllele enables allele-specific splicing analysis, a unique feature not afforded by other methods. Applied to a lung cancer scRNA-seq dataset, scAllele identified variants with strong allelic linkage to alternative splicing, some of which are cancer specific and enriched in cancer-relevant pathways. scAllele represents a versatile tool to uncover multilayer information and previously unidentified biological insights from scRNA-seq data.

Using third-generation sequencers, long-read RNA-seq is increasingly applied in transcriptomic studies given its major advantage in characterizing full-length transcripts. A number of methods have been developed to analyze this new type of data for transcript isoforms and their abundance. Another application, which is significantly under-explored, is to identify and analyze single nucleotide variants (SNVs) in the RNA. Identification of SNVs, such as genetic mutations or RNA editing sites, is fundamental to many biomedical questions. In long-read RNA-seq, SNV analysis presents significant challenges, due to the well-known relatively high error rates of the third-generation sequencers. Here, we present the first study to detect and analyze RNA editing sites in long-read RNA-seq. Our new method, L-GIREMI, effectively handles sequencing errors and biases in the reads, and uses a model-based approach to score RNA editing sites. Applied to PacBio long-read RNA-seq data, L-GIREMI affords a high accuracy in RNA editing identification. In addition, the unique advantage of long reads allowed us to uncover novel insights about RNA editing occurrences in single molecules and double-stranded RNA (dsRNA) structures. L-GIREMI provides a valuable means to study RNA nucleotide variants in long-read RNA-seq.

Post-transcriptional regulation plays a key role in transcriptome reprogramming during cardiac pathogenesis. In previous studies, we have identified that cardiac enriched RNA-binding protein, RBFox1 plays key role in cardiac hypertrophy through mRNA alternative splicing regulation in nuclei. However, RBFox1 gene also generates a cytosolic isoform (RBFox1c), suggesting additional functions of post-transcriptional regulation in heart. In adult heart, RBFox1c mRNA constituted ~ 40% of total RBFox1 level but was significantly repressed in pressure-overloaded failing mouse heart. Using CRISPR-Cas9 technology, we have established an isoform specific RBFox1c-cKO mouse. At baseline inactivation of RBFox1c led to decreased cardiac function along with induction of cardiac fibrosis. RBFox1c-cKO mice also showed macrophages infiltration into myocardium post 7days MI. In contrast, restoration of RBFox1c expression in adult intact hearts significantly reduced cardiac fibrosis post stress. RNA-seq analyses in RBFox1c expressing cardiomyocytes showed that RBFox1c specifically suppressed the expression of pro-inflammatory genes. Secondly, CLIP-Seq analysis and targeted RNA-IP showed that RBFox1c could directly interact with inflammatory pathway mRNAs. These results suggested the inflammatory mRNAs are direct downstream targets regulated by RBFox1c. Using both in vitro cultured cardiomyocytes and intact mouse hearts, we demonstrated that expression of RBFox1c reduces pro-inflammatory mRNA expression at baseline and upon hypertrophy stimulation. Lastly, we characterized the interactome of RBFox1c through proteomic analysis and found RBFox1c specifically interacted with a component of the RNA NMD machinery-Upf1. RBFox1c interaction with Upf1 in cardiomyocytes was diminished upon hypertrophic stress. Furthermore, by inactivation of Upf1 via siRNA, we demonstrated that RBFox1c mediated repression of proinflammatory genes was Upf1 dependent.RBFox1 regulates cardiac transcriptome reprogramming in two post-transcriptional processes via distinct isoforms. While the RBFox1n regulates RNA splicing, the RBFox1c functions through targeted mRNA repression of proinflammatory genes by recruitment of Upf1 mediated RNA degradation.

Motivation Alternative polyadenylation (polyA) sites near the 3′ end of a pre-mRNA create multiple mRNA transcripts with different 3′ untranslated regions (3′ UTRs). The sequence elements of a 3′ UTR are essential for many biological activities such as mRNA stability, sub-cellular localization, protein translation, protein binding and translation efficiency. Moreover, numerous studies in the literature have reported the correlation between diseases and the shortening (or lengthening) of 3′ UTRs. As alternative polyA sites are common in mammalian genes, several machine learning tools have been published for predicting polyA sites from sequence data. These tools either consider limited sequence features or use relatively old algorithms for polyA site prediction. Moreover, none of the previous tools consider RNA secondary structures as a feature to predict polyA sites. Results In this paper, we propose a new deep learning model, called DeepPASTA, for predicting polyA sites from both sequence and RNA secondary structure data. The model is then extended to predict tissue-specific polyA sites. Moreover, the tool can predict the most dominant (i.e. frequently used) polyA site of a gene in a specific tissue and relative dominance when two polyA sites of the same gene are given. Our extensive experiments demonstrate that DeepPASTA signisficantly outperforms the existing tools for polyA site prediction and tissue-specific relative and absolute dominant polyA site prediction. Availability and implementation https://github.com/arefeen/DeepPASTA Supplementary information Supplementary data are available at Bioinformatics online.

Transcriptomic analyses of postmortem brains have begun to elucidate molecular abnormalities in autism spectrum disorder (ASD). However, a crucial pathway involved in synaptic development, RNA editing, has not yet been studied on a genome-wide scale. Here we profiled global patterns of adenosine-to-inosine (A-to-I) editing in a large cohort of postmortem brains of people with ASD. We observed a global bias for hypoediting in ASD brains, which was shared across brain regions and involved many synaptic genes. We show that the Fragile X proteins FMRP and FXR1P interact with RNA-editing enzymes (ADAR proteins) and modulate A-to-I editing. Furthermore, we observed convergent patterns of RNA-editing alterations in ASD and Fragile X syndrome, establishing this as a molecular link between these related diseases. Our findings, which are corroborated across multiple data sets, including dup15q (genomic duplication of 15q11.2-13.1) cases associated with intellectual disability, highlight RNA-editing dysregulation in ASD and reveal new mechanisms underlying this disorder.

BackgroundRNA editing generates modifications to the RNA sequences, thereby increasing protein diversity and shaping various layers of gene regulation. Recent studies have revealed global shifts in editing levels across many cancer types, as well as a few specific mechanisms implicating individual sites in tumorigenesis or metastasis. However, most tumor-associated sites, predominantly in noncoding regions, have unknown functional relevance.ResultsHere, we carry out integrative analysis of RNA editing profiles between epithelial and mesenchymal tumors, since epithelial-mesenchymal transition is a key paradigm for metastasis. We identify distinct editing patterns between epithelial and mesenchymal tumors in seven cancer types using TCGA data, an observation further supported by single-cell RNA sequencing data and ADAR perturbation experiments in cell culture. Through computational analyses and experimental validations, we show that differential editing sites between epithelial and mesenchymal phenotypes function by regulating mRNA abundance of their respective genes. Our analysis of RNA-binding proteins reveals ILF3 as a potential regulator of this process, supported by experimental validations. Consistent with the known roles of ILF3 in immune response, epithelial-mesenchymal differential editing sites are enriched in genes involved in immune and viral processes. The strongest target of editing-dependent ILF3 regulation is the transcript encoding PKR, a crucial player in immune and viral response.ConclusionsOur study reports widespread differences in RNA editing between epithelial and mesenchymal tumors and a novel mechanism of editing-dependent regulation of mRNA abundance. It reveals the broad impact of RNA editing in cancer and its relevance to cancer-related immune pathways.

Motivation RNA-sequencing (RNA-seq) enables global identification of RNA-editing sites in biological systems and disease. A salient step in many studies is to identify editing sites that statistically associate with treatment (e.g. case versus control) or covary with biological factors, such as age. However, RNA-seq has technical features that incumbent tests (e.g. t-test and linear regression) do not consider, which can lead to false positives and false negatives. Results In this study, we demonstrate the limitations of currently used tests and introduce the method, RNA-editing tests (REDITs), a suite of tests that employ beta-binomial models to identify differential RNA editing. The tests in REDITs have higher sensitivity than other tests, while also maintaining the type I error (false positive) rate at the nominal level. Applied to the GTEx dataset, we unveil RNA-editing changes associated with age and gender, and differential recoding profiles between brain regions. Availability and implementation REDITs are implemented as functions in R and freely available for download at https://github.com/gxiaolab/REDITs. The repository also provides a code example for leveraging parallelization using multiple cores.

Alcohol (ethanol, EtOH) consumption during pregnancy can result in fetal alcohol spectrum disorders (FASDs), which are characterized by prenatal and postnatal growth restriction and craniofacial dysmorphology. Recently, cell-derived extracellular vesicles, including exosomes and microvesicles containing several species of RNAs (exRNAs), have emerged as a mechanism of cell-to-cell communication. However, EtOH’s effects on the biogenesis and function of non-coding exRNAs during fetal development have not been explored. Therefore, we studied the effects of maternal EtOH exposure on the composition of exosomal RNAs in the amniotic fluid (AF) using rat fetal alcohol exposure (FAE) model. Through RNA-Seq analysis we identified and verified AF exosomal miRNAs with differential expression levels specifically associated with maternal EtOH exposure. Uptake of purified FAE AF exosomes by rBMSCs resulted in significant alteration of molecular markers associated with osteogenic differentiation of rBMSCs. We also determined putative functional roles for AF exosomal miRNAs (miR-199a-3p, miR-214-3p and let-7g) that are dysregulated by FAE in osteogenic differentiation of rBMSCs. Our results demonstrate that FAE alters AF exosomal miRNAs and that exosomal transfer of dysregulated miRNAs has significant molecular effects on stem cell regulation and differentiation. Our results further suggest the usefulness of assessing molecular alterations in AF exRNAs to study the mechanisms of FAE teratogenesis that should be further investigated by using an in vivo model.

BACKGROUND It was recently discovered that abundant and stable extracellular RNA (exRNA) species exist in bodily fluids. Saliva is an emerging biofluid for biomarker development for noninvasive detection and screening of local and systemic diseases. Use of RNA-Sequencing (RNA-Seq) to profile exRNA is rapidly growing; however, no single preparation and analysis protocol can be used for all biofluids. Specifically, RNA-Seq of saliva is particularly challenging owing to high abundance of bacterial contents and low abundance of salivary exRNA. Given the laborious procedures needed for RNA-Seq library construction, sequencing, data storage, and data analysis, saliva-specific and optimized protocols are essential. METHODS We compared different RNA isolation methods and library construction kits for long and small RNA sequencing. The role of ribosomal RNA (rRNA) depletion also was evaluated. RESULTS The miRNeasy Micro Kit (Qiagen) showed the highest total RNA yield (70.8 ng/mL cell-free saliva) and best small RNA recovery, and the NEBNext library preparation kits resulted in the highest number of detected human genes [5649–6813 at 1 reads per kilobase RNA per million mapped (RPKM)] and small RNAs [482–696 microRNAs (miRNAs) and 190–214 other small RNAs]. The proportion of human RNA-Seq reads was much higher in rRNA-depleted saliva samples (41%) than in samples without rRNA depletion (14%). In addition, the transfer RNA (tRNA)-derived RNA fragments (tRFs), a novel class of small RNAs, were highly abundant in human saliva, specifically tRF-4 (4%) and tRF-5 (15.25%). CONCLUSIONS Our results may help in selection of the best adapted methods of RNA isolation and small and long RNA library constructions for salivary exRNA studies.

Motivation The length of the 3′ untranslated region (3′ UTR) of an mRNA is essential for many biological activities such as mRNA stability, sub-cellular localization, protein translation, protein binding and translation efficiency. Moreover, correlation between diseases and the shortening (or lengthening) of 3′ UTRs has been reported in the literature. This length is largely determined by the polyadenylation cleavage site in the mRNA. As alternative polyadenylation (APA) sites are common in mammalian genes, several tools have been published recently for detecting APA sites from RNA-Seq data or performing shortening/lengthening analysis. These tools consider either up to only two APA sites in a gene or only APA sites that occur in the last exon of a gene, although a gene may generally have more than two APA sites and an APA site may sometimes occur before the last exon. Furthermore, the tools are unable to integrate the analysis of shortening/lengthening events with APA site detection. Results We propose a new tool, called TAPAS, for detecting novel APA sites from RNA-Seq data. It can deal with more than two APA sites in a gene as well as APA sites that occur before the last exon. The tool is based on an existing method for finding change points in time series data, but some filtration techniques are also adopted to remove change points that are likely false APA sites. It is then extended to identify APA sites that are expressed differently between two biological samples and genes that contain 3′ UTRs with shortening/lengthening events. Our extensive experiments on simulated and real RNA-Seq data demonstrate that TAPAS outperforms the existing tools for APA site detection or shortening/lengthening analysis significantly. Availability and implementation https://github.com/arefeen/TAPAS Supplementary information Supplementary data are available at Bioinformatics online.

Alternative transcription start (ATS) and alternative polyadenylation (APA) create alternative RNA isoforms and modulate many aspects of RNA expression and protein production. However, ATS and APA remain difficult to detect in RNA sequencing (RNA-seq). Here, we developed mountainClimber, a de novo cumulative-sum-based approach to identify ATS and APA as change points. Unlike many existing methods, mountainClimber runs on a single sample and identifies multiple ATS or APA sites anywhere in the transcript. We analyzed 2,342 GTEx samples (36 tissues, 215 individuals) and found that tissue type is the predominant driver of transcript end variations. 75% and 65% of genes exhibited differential APA and ATS across tissues, respectively. In particular, testis displayed longer 5' untranslated regions (UTRs) and shorter 3' UTRs, often in genes related to testis-specific biology. Overall, we report the largest study of transcript ends across human tissues to our knowledge. mountainClimber is available at github.com/gxiaolab/mountainClimber.

Adenosine-to-inosine (A-to-I) editing, mediated by the ADAR enzymes, diversifies the transcriptome by altering RNA sequences. Recent studies reported global changes in RNA editing in disease and development. Such widespread editing variations necessitate an improved understanding of the regulatory mechanisms of RNA editing. Here, we study the roles of >200 RNA-binding proteins (RBPs) in mediating RNA editing in two human cell lines. Using RNA-sequencing and global protein-RNA binding data, we identify a number of RBPs as key regulators of A-to-I editing. These RBPs, such as TDP-43, DROSHA, NF45/90 and Ro60, mediate editing through various mechanisms including regulation of ADAR1 expression, interaction with ADAR1, and binding to Alu elements. We highlight that editing regulation by Ro60 is consistent with the global up-regulation of RNA editing in systemic lupus erythematosus. Additionally, most key editing regulators act in a cell type-specific manner. Together, our work provides insights for the regulatory mechanisms of RNA editing., Giovanni Quinones-Valdez et al. examined the role of over 200 RNA-binding proteins in mediating A-to-I RNA editing. They identified several RNA-binding proteins that regulate ADAR1 expression, interaction, or binding with Alu elements in a cell type-specific manner.

Identification of functional genetic variants and elucidation of their regulatory mechanisms represent significant challenges of the post-genomic era. A poorly understood topic is the involvement of genetic variants in mediating post-transcriptional RNA processing, including alternative splicing. Thus far, little is known about the genomic, evolutionary, and regulatory features of genetically modulated alternative splicing (GMAS). Here, we systematically identified intronic tag variants for genetic modulation of alternative splicing using RNA-seq data specific to cellular compartments. Combined with our previous method that identifies exonic tags for GMAS, this study yielded 622 GMAS exons. We observed that GMAS events are highly cell type independent, indicating that splicing-altering genetic variants could have widespread function across cell types. Interestingly, GMAS genes, exons, and single-nucleotide variants (SNVs) all demonstrated positive selection or accelerated evolution in primates. We predicted that GMAS SNVs often alter binding of splicing factors, with SRSF1 affecting the most GMAS events and demonstrating global allelic binding bias. However, in contrast to their GMAS targets, the predicted splicing factors are more conserved than expected, suggesting that cis-regulatory variation is the major driving force of splicing evolution. Moreover, GMAS-related splicing factors had stronger consensus motifs than expected, consistent with their susceptibility to SNV disruption. Intriguingly, GMAS SNVs in general do not alter the strongest consensus position of the splicing factor motif, except the more than 100 GMAS SNVs in linkage disequilibrium with polymorphisms reported by genome-wide association studies. Our study reports many GMAS events and enables a better understanding of the evolutionary and regulatory features of this phenomenon.

Male vertebrate social displays vary from physically simple to complex, with the latter involving exquisite motor command of the body and appendages. Studies of these displays have, in turn, provided substantial insight into neuromotor mechanisms. The neotropical golden-collared manakin (Manacus vitellinus) has been used previously as a model to investigate intricate motor skills because adult males of this species perform an acrobatic and androgen-dependent courtship display. To support this behavior, these birds express elevated levels of androgen receptors (AR) in their skeletal muscles. Here we use RNA sequencing to explore how testosterone (T) modulates the muscular transcriptome to support male manakin courtship displays. In addition, we explore how androgens influence gene expression in the muscles of the zebra finch (Taenopygia guttata), a model passerine bird with a limited courtship display and minimal muscle AR. We identify androgen-dependent, muscle-specific gene regulation in both species. In addition, we identify manakin-specific effects that are linked to muscle use during the manakin display, including androgenic regulation of genes associated with muscle fiber contractility, cellular homeostasis, and energetic efficiency. Overall, our results point to numerous genes and gene networks impacted by androgens in male birds, including some that underlie optimal muscle function necessary for performing acrobatic display routines. Manakins are excellent models to explore gene regulation promoting athletic ability.

Background: RBFox1 is known to be an RNA splicing regulator with enriched expression in cardiac muscle. Loss of RBFox1 expression is a molecular hallmark associated with heart failure in mouse and human. Genetic manipulation of RBFox1 reveals a major function of RBFox1 in the pathogenesis of cardiac hypertrophy, fibrosis and dysfunction under pathological stresses. However, much of our current knowledge about RBFox1 focuses on the nuclear RBFox1 with a major impact on global alternative splicing changes. Yet, RBFox1 gene also generates a cytosolic isoform-RBFox1c through alternative splicing, and the specific function of RBFox1c has not been characterized. Goal and Methods: This study investigated the functional impact and the underlying mechanism of the RBFox1c in cardiac pathological remodeling and targeted gene regulation. Results: RBFox1c expression is significantly repressed in stressed cardiomyocytes in vivo and in vitro. RNA-Seq combined with IPA analyses from cardiomyocytes revealed RBFox1c but not the nucleus RBFox1 specifically suppressed pro-inflammatory genes. In the cardiac specific RBFox1 knockout mice, enhanced cardiac fibrosis is observed following I/R injury associated with elevated expression of pro-inflammatory genes. In contrast, cardiac specific expression of RBFox1c reduced cardiac fibrosis and inflammatory gene expression following pressure-overload and myocardial infarct injury associated with improved ejection fraction. Motif enrichment analysis identified significant enrichment of the RBFox1 binding motif in the 3’UTR of the RBFox1c regulated genes. We performed CLIP analysis followed by RT-PCR and observed RBFox1c interacted with inflammatory gene 3’UTR. Lastly, we explored the interactome of RBFox1c and found RBFox1c specifically interacted with a component of RNA decay machinery-Upf1. Indeed, while expression of RBFox1c repressed inflammatory gene expression, inactivation of Upf1 abolished this effect in cardiomyocytes. Conclusion: RBFox1 regulates cardiac transcriptome reprogramming at two post-transcriptional steps. The RBFox1 nuclei isoform regulates RNA splicing reprogramming, and the RBFox1c represses proinflammatory genes via recruitment of Upf1 mediated RNA degradation.

In eukaryotes, nascent RNA transcripts undergo an intricate series of RNA processing steps to achieve mRNA maturation. RNA editing and alternative splicing are two major RNA processing steps that can introduce significant modifications to the final gene products. By tackling these processes in isolation, recent studies have enabled substantial progress in understanding their global RNA targets and regulatory pathways. However, the interplay between individual steps of RNA processing, an essential aspect of gene regulation, remains poorly understood. By sequencing the RNA of different subcellular fractions, we examined the timing of adenosine-to-inosine (A-to-I) RNA editing and its impact on alternative splicing. We observed that >95% A-to-I RNA editing events occurred in the chromatin-associated RNA prior to polyadenylation. We report about 500 editing sites in the 3′ acceptor sequences that can alter splicing of the associated exons. These exons are highly conserved during evolution and reside in genes with important cellular function. Furthermore, we identified a second class of exons whose splicing is likely modulated by RNA secondary structures that are recognized by the RNA editing machinery. The genome-wide analyses, supported by experimental validations, revealed remarkable interplay between RNA editing and splicing and expanded the repertoire of functional RNA editing sites.

The epithelial-mesenchymal transition (EMT) is a fundamental developmental process that is abnormally activated in cancer metastasis. Dynamic changes in alternative splicing occur during EMT. ESRP1 and hnRNPM are splicing regulators that promote an epithelial splicing program and a mesenchymal splicing program, respectively. The functional relationships between these splicing factors in the genome-scale remain elusive. Comparing alternative splicing targets of hnRNPM and ESRP1 revealed that they co-regulate a set of cassette exon events, with the majority showing discordant splicing regulation. hnRNPM discordantly regulated splicing events show a positive correlation with splicing during EMT while concordant splicing events do not, highlighting the antagonistic role of hnRNPM and ESRP1 during EMT. Motif enrichment analysis near co-regulated exons identifies guanine-uridine rich motifs downstream of hnRNPM-repressed and ESRP1-enhanced exons, supporting a model of competitive binding to these cis-elements to antagonize alternative splicing. The set of co-regulated exons are enriched in genes associated with cell-migration and cytoskeletal reorganization, which are pathways associated with EMT. Splicing levels of co-regulated exons are associated with breast cancer patient survival and correlate with gene sets involved in EMT and breast cancer subtypes. These data identify complex modes of interaction between hnRNPM and ESRP1 in regulation of splicing in disease-relevant contexts.

The epithelial–mesenchymal transition (EMT) is a fundamental developmental process that is abnormally activated in cancer metastasis. Dynamic changes in alternative splicing occur during EMT. ESRP1 and hnRNPM are splicing regulators that promote an epithelial splicing program and a mesenchymal splicing program, respectively. The functional relationships between these splicing factors in the genome scale remain elusive. Comparing alternative splicing targets of hnRNPM and ESRP1 revealed that they coregulate a set of cassette exon events, with the majority showing discordant splicing regulation. Discordant splicing events regulated by hnRNPM show a positive correlation with splicing during EMT; however, concordant events do not, indicating the role of hnRNPM in regulating alternative splicing during EMT is more complex than previously understood. Motif enrichment analysis near hnRNPM–ESRP1 coregulated exons identifies guanine–uridine rich motifs downstream from hnRNPM-repressed and ESRP1-enhanced exons, supporting a general model of competitive binding to these cis-elements to antagonize alternative splicing. The set of coregulated exons are enriched in genes associated with cell migration and cytoskeletal reorganization, which are pathways associated with EMT. Splicing levels of coregulated exons are associated with breast cancer patient survival and correlate with gene sets involved in EMT and breast cancer subtyping. This study identifies complex modes of interaction between hnRNPM and ESRP1 in regulation of splicing in disease-relevant contexts.

RNA editing generates post-transcriptional sequence changes that can be deduced from RNA-seq data, but detection typically requires matched genomic sequence or multiple related expression data sets. We developed the GIREMI tool (genome-independent identification of RNA editing by mutual information; https://www.ibp.ucla.edu/research/xiao/GIREMI.html) to predict adenosine-to-inosine editing accurately and sensitively from a single RNA-seq data set of modest sequencing depth. Using GIREMI on existing data, we observed tissue-specific and evolutionary patterns in editing sites in the human population.

Motivation Analysis of RNA sequencing (RNA-Seq) data in human saliva is challenging. Lack of standardization and unification of the bioinformatic procedures undermines saliva‘s diagnostic potential. Thus, it motivated us to perform this study. Results We applied principal pipelines for bioinformatic analysis of small RNA-Seq data of saliva of 98 healthy Korean volunteers including either direct or indirect mapping of the reads to the human genome using Bowtie1. Analysis of alignments to exogenous genomes by another pipeline revealed that almost all of the reads map to bacterial genomes. Thus, salivary exRNA has fundamental properties that warrant the design of unique additional steps while performing the bioinformatic analysis. Our pipelines can serve as potential guidelines for processing of RNA-Seq data of human saliva. Availability and implementation Processing and analysis results of the experimental data generated by the exceRpt (v4.6.3) small RNA-seq pipeline (github.gersteinlab.org/exceRpt) are available from exRNA atlas (exrna-atlas.org). Alignment to exogenous genomes and their quantification results were used in this paper for the analyses of small RNAs of exogenous origin.

Extensive pre-mRNA back-splicing generates numerous circular RNAs (circRNAs) in human transcriptome. However, the biological functions of these circRNAs remain largely unclear. Here we report that N6-methyladenosine (m6A), the most abundant base modification of RNA, promotes efficient initiation of protein translation from circRNAs in human cells. We discover that consensus m6A motifs are enriched in circRNAs and a single m6A site is sufficient to drive translation initiation. This m6A-driven translation requires initiation factor eIF4G2 and m6A reader YTHDF3, and is enhanced by methyltransferase METTL3/14, inhibited by demethylase FTO, and upregulated upon heat shock. Further analyses through polysome profiling, computational prediction and mass spectrometry reveal that m6A-driven translation of circRNAs is widespread, with hundreds of endogenous circRNAs having translation potential. Our study expands the coding landscape of human transcriptome, and suggests a role of circRNA-derived proteins in cellular responses to environmental stress.

To determine the role for mutations of MECP2 in Rett Syndrome, we generated isogenic lines of human iPSCs (hiPSCs), neural progenitor cells (NPCs), and neurons from patient fibroblasts with and without MECP2 expression in an attempt to recapitulate disease phenotypes in vitro. Molecular profiling uncovered neuronal specific gene expression changes including induction of a Senescence Associated Secretory Phenotype (SASP) program. Patient derived Neurons made without MECP2 show signs of stress, including induction of p53, and senescence. The induction of p53 appeared to affect dendritic branching in Rett neurons, as p53 inhibition restored dendritic complexity. These disease-in-a-dish data suggest that loss of MECP2 can lead to dendritic defects due to an increase in aspects of neuronal aging.

Mutations that lead to splicing defects can have severe consequences on gene function and cause disease. Here, we explore how human genetic variation affects exon recognition by developing a multiplexed functional assay of splicing using Sort-seq (MFASS). We assayed 27,733 variants in the Exome Aggregation Consortium (ExAC) within or adjacent to 2,198 human exons in the MFASS minigene reporter and found that 3.8% (1,050) of variants, most of which are extremely rare, led to large-effect splice-disrupting variants (SDVs). Importantly, we find that 83% of SDVs are located outside of canonical splice sites, are distributed evenly across distinct exonic and intronic regions, and are difficult to predict a priori. Our results indicate extant, rare genetic variants can have large functional effects on splicing at appreciable rates, even outside the context of disease, and MFASS enables their empirical assessment at scale.

Understanding protein-DNA interactions, specifically transcription factor (TF) and transcription factor binding site (TFBS) bindings, is crucial in deciphering gene regulation. The recent associated TF-TFBS pattern discovery combines one-sided motif discovery on both the TF and the TFBS sides. Using sequences only, it identifies the short protein-DNA binding cores available only in high-resolution 3D structures. The discovered patterns lead to promising subtype and disease analysis applications. While the related studies use either association rule mining or existing TFBS annotations, none has proposed any formal unified (both-sided) model to prioritize the top verifiable associated patterns. We propose the unified scores and develop an effective pipeline for associated TF-TFBS pattern discovery. Our stringent instance-level evaluations show that the patterns with the top unified scores match with the binding cores in 3D structures considerably better than the previous works, where up to 90 percent of the top 20 scored patterns are verified. We also introduce extended verification from literature surveys, where the high unified scores correspond to even higher verification percentage. The top scored patterns are confirmed to match the known WRKY binding cores with no available 3D structures and agree well with the top binding affinities of in vivo experiments.

RNA-sequencing (RNA-Seq) technologies hold enormous promise for novel discoveries in genomics and transcriptomics. In the past year, a surge of reports has analyzed RNA-Seq data to gain a global view of the RNA editome. Opposing results have been presented, giving rise to extensive debate surrounding one of the first such studies in which a daunting list of all 12 types of RNA–DNA differences (RDDs) were identified. Although a consensus is forming that some of the initial “paradigm-shifting” results of this study may be questionable, recent reports on this topic differed in terms of the number and relative abundance of each type of RDD. Many outstanding issues exist, most importantly, the choice of bioinformatic approaches. Here we discuss the critical data analysis and experimental design issues of such studies to enable improved systematic investigation of the largely unexplored frontier of single-nucleotide variants in RNA.

Background— Cardiac maturation during perinatal transition of heart is critical for functional adaptation to hemodynamic load and nutrient environment. Perturbation in this process has major implications in congenital heart defects. Transcriptome programming during perinatal stages is an important information but incomplete in current literature, particularly, the expression profiles of the long noncoding RNAs (lncRNAs) are not fully elucidated. Methods and Results— From comprehensive analysis of transcriptomes derived from neonatal mouse heart left and right ventricles, a total of 45 167 unique transcripts were identified, including 21 916 known and 2033 novel lncRNAs. Among these lncRNAs, 196 exhibited significant dynamic regulation along maturation process. By implementing parallel weighted gene co-expression network analysis of mRNA and lncRNA data sets, several lncRNA modules coordinately expressed in a developmental manner similar to protein coding genes, while few lncRNAs revealed chamber-specific patterns. Out of 2262 lncRNAs located within 50 kb of protein coding genes, 5% significantly correlate with the expression of their neighboring genes. The impact of Ppp1r1b-lncRNA on the corresponding partner gene Tcap was validated in cultured myoblasts. This concordant regulation was also conserved in human infantile hearts. Furthermore, the Ppp1r1b-lncRNA/Tcap expression ratio was identified as a molecular signature that differentiated congenital heart defect phenotypes. Conclusions— The study provides the first high-resolution landscape on neonatal cardiac lncRNAs and reveals their potential interaction with mRNA transcriptome during cardiac maturation. Ppp1r1b-lncRNA was identified as a regulator of Tcap expression, with dynamic interaction in postnatal cardiac development and congenital heart defects.

Previous studies had shown that the integration of genome wide expression profiles, in metabolic tissues, with genetic and phenotypic variance, provided valuable insight into the underlying molecular mechanisms. We used RNA-Seq to characterize hypothalamic transcriptome in 99 inbred strains of mice from the Hybrid Mouse Diversity Panel (HMDP), a reference resource population for cardiovascular and metabolic traits. We report numerous novel transcripts supported by proteomic analyses, as well as novel non coding RNAs. High resolution genetic mapping of transcript levels in HMDP, reveals both local and trans expression Quantitative Trait Loci (eQTLs) demonstrating 2 trans eQTL 'hotspots' associated with expression of hundreds of genes. We also report thousands of alternative splicing events regulated by genetic variants. Finally, comparison with about 150 metabolic and cardiovascular traits revealed many highly significant associations. Our data provide a rich resource for understanding the many physiologic functions mediated by the hypothalamus and their genetic regulation. DOI: http://dx.doi.org/10.7554/eLife.15614.001, eLife digest Metabolism is a term that describes all the chemical reactions that are involved in keeping a living organism alive. Diseases related to metabolism – such as obesity, heart disease and diabetes – are a major health problem in the Western world. The causes of these diseases are complex and include both environmental factors, such as diet and exercise, and genetics. Indeed, many genetic variants that contribute to obesity have been uncovered in both humans and mice. However, it is only dimly understood how these genetic variants affect the underlying networks of interacting genes that cause metabolic disorders. Measuring gene activity or expression, and tracing how genetic instructions are carried from DNA into RNA and proteins, can reliably identify groups of genes that correlate with metabolic traits in specific organs. This strategy was successfully used in previous studies to reveal new information about abnormalities linked to obesity in specific tissues such as the liver and fat tissues. It was also shown that this approach might suggest new molecules that could be targeted to treat metabolic disorders. A brain region called the hypothalamus is key to the control of metabolism, including feeding behavior and obesity. Hasin-Brumshtein et al. set out to explore gene expression in the hypothalamus of 99 different strains of mice, in the hope that the data will help identify new connections between gene expression and metabolism. This approach showed that thousands of new and known genes are expressed in the mouse hypothalamus, some of which coded for proteins, and some of which did not. Hasin-Brumshtein et al. uncovered two genetic variants that controlled the expression of hundreds of other genes. Further analysis then revealed thousands of genetic variants that regulated the expression of, and type of RNA (so-called "spliceforms") produced from neighboring genes. Also, the expression of many individual genes showed significant similarities with about 150 metabolic measurements that had been evaluated previously in the mice. This new dataset is a unique resource that can be coupled with different approaches to test existing ideas and develop new ones about the role of particular genes or genetic mechanisms in obesity. Future studies will likely focus on new genes that show strong associations with attributes that are relevant to metabolic disorders, such as insulin levels, weight and fat mass. DOI: http://dx.doi.org/10.7554/eLife.15614.002

Background: Left ventricle (LV) and right ventricle (RV) in mouse heart undergo dramatically different chamber-specific remodeling after birth, leading to rapid increase in LV vs. RV chamber size. However, the underlying regulatory mechanism mediating chamber specific remodeling process remains enigmatic. Results and Methods: In neonatal mouse heart, p38 MAP kinase activity is dynamically activated in a chamber specific manner. p38 activity is specifically elevated in RV comparing to LV at E18.5, postnatal day 3 (P3) and P7 stages whereas p38 activity is lower in both ventricles at P0 and P1. In mouse heart with cardiomyocyte specific-knockout of p38α and β (p38ab-cdKO), total p38 activity was diminished in both chambers. The p38ab-cdKO mice had significant neonatal lethality associated with RV specific chamber enlargement and significant increase in both RV wall thickness (RVW) and inner diameter of RV (RVID) as early as P3. Interestingly, p38 inactivation suppressed myocyte apoptotic activity specifically in RV while increased RV myocyte proliferation and hypertrophy during neonatal period. Unexpectedly, RNA-seq results implicated Xbp1 mediated transcriptional regulation significantly contributing to p38 dependent transcriptome reprogramming in RV. Indeed, IRE1α expression in neonatal cardiomyocyte is sufficient to induce proliferation in vitro. Furthermore, knockdown of Xbp1 blunted p38 inhibition-induced myocyte proliferation, suggesting that IRE1a/Xbp1 mediate p38 signaling in neonatal myocyte proliferation. Conclusion: Chamber-specific remodeling in neonatal heart involves temporally regulated and RV specific p38 MAP kinase activity. RV specific myocyte proliferation and hypertrophy concurrent with RV specific programmed myocyte death is orchestrated by two innate stress-response pathways, p38 and Xbp1.

Nutrition plays a significant role in the increasing prevalence of metabolic and brain disorders. Here we employ systems nutrigenomics to scrutinize the genomic bases of nutrient–host interaction underlying disease predisposition or therapeutic potential. We conducted transcriptome and epigenome sequencing of hypothalamus (metabolic control) and hippocampus (cognitive processing) from a rodent model of fructose consumption, and identified significant reprogramming of DNA methylation, transcript abundance, alternative splicing, and gene networks governing cell metabolism, cell communication, inflammation, and neuronal signaling. These signals converged with genetic causal risks of metabolic, neurological, and psychiatric disorders revealed in humans. Gene network modeling uncovered the extracellular matrix genes Bgn and Fmod as main orchestrators of the effects of fructose, as validated using two knockout mouse models. We further demonstrate that an omega-3 fatty acid, DHA, reverses the genomic and network perturbations elicited by fructose, providing molecular support for nutritional interventions to counteract diet-induced metabolic and brain disorders. Our integrative approach complementing rodent and human studies supports the applicability of nutrigenomics principles to predict disease susceptibility and to guide personalized medicine., Highlights • Fructose promotes transcriptomic and epigenomic reprogramming to perturb brain networks linking metabolism and brain function. • The extracellular matrix genes Bgn and Fmod emerge as key regulators of gene networks responsive to fructose. • The omega-3 fatty acid DHA reverses fructose-induced genomic and network reprogramming. Meng et al. report fructose as a powerful inducer of genomic and epigenomic variability with the capacity to reorganize gene networks critical for central metabolic regulation and neuronal processes in the brain; conversely, an omega-3 fatty acid, DHA, has the potential to normalize the genomic impact of fructose. Our findings help explain the pathogenic actions of fructose on prevalent metabolic and brain disorders and provide proof-of-concept for nutritional remedies supported by nutrigenomics evidence. Our integrative approach complementing rodent and human studies supports the applicability of nutrigenomics principles to predict disease susceptibility and to guide personalized medicine.

In mammals, small RNAs are important players in post-transcriptional gene regulation. While their roles in mRNA destabilization and translational repression are well appreciated, their involvement in endonucleolytic cleavage of target RNAs is poorly understood. Very few microRNAs are known to guide RNA cleavage. Endogenous small interfering RNAs are expected to induce target cleavage, but their target genes remain largely unknown. We report a systematic study of small RNA-mediated endonucleolytic cleavage in mouse through integrative analysis of small RNA and degradome sequencing data without imposing any bias toward known small RNAs. Hundreds of small cleavage-inducing RNAs and their cognate target genes were identified, significantly expanding the repertoire of known small RNA-guided cleavage events. Strikingly, both small RNAs and their target sites demonstrated significant overlap with retrotransposons, providing evidence for the long-standing speculation that retrotransposable elements in mRNAs are leveraged as signals for gene targeting. Furthermore, our analysis showed that the RNA cleavage pathway is also present in human cells but affecting a different repertoire of retrotransposons. These results show that small RNA-guided cleavage is more widespread than previously appreciated. Their impact on retrotransposons in non-coding regions shed light on important aspects of mammalian gene regulation.

Epigenetic reprogramming is a critical process of pathological gene induction during cardiac hypertrophy and remodeling, but the underlying regulatory mechanisms remain to be elucidated. Here we identified a heart-enriched long noncoding (lnc)RNA, named cardiac-hypertrophy-associated epigenetic regulator (Chaer), which is necessary for the development of cardiac hypertrophy. Mechanistically, Chaer directly interacts with the catalytic subunit of polycomb repressor complex 2 (PRC2). This interaction, which is mediated by a 66-mer motif in Chaer, interferes with PRC2 targeting to genomic loci, thereby inhibiting histone H3 lysine 27 methylation at the promoter regions of genes involved in cardiac hypertrophy. The interaction between Chaer and PRC2 is transiently induced after hormone or stress stimulation in a process involving mammalian target of rapamycin complex 1, and this interaction is a prerequisite for epigenetic reprogramming and induction of genes involved in hypertrophy. Inhibition of Chaer expression in the heart before, but not after, the onset of pressure overload substantially attenuates cardiac hypertrophy and dysfunction. Our study reveals that stress-induced pathological gene activation in the heart requires a previously uncharacterized lncRNA-dependent epigenetic checkpoint.

Motivation: Accurate identification of genetic variants such as single-nucleotide polymorphisms (SNPs) or RNA editing sites from RNA-Seq reads is important, yet challenging, because it necessitates a very low false-positive rate in read mapping. Although many read aligners are available, no single aligner was specifically developed or tested as an effective tool for SNP and RNA editing prediction. Results: We present RASER, an accurate read aligner with novel mapping schemes and index tree structure that aims to reduce false-positive mappings due to existence of highly similar regions. We demonstrate that RASER shows the best mapping accuracy compared with other popular algorithms and highest sensitivity in identifying multiply mapped reads. As a result, RASER displays superb efficacy in unbiased mapping of the alternative alleles of SNPs and in identification of RNA editing sites. Availability and implementation: RASER is written in C++ and freely available for download at https://github.com/jaegyoonahn/RASER. Contact: gxxiao@ucla.edu Supplementary information: Supplementary data are available at Bioinformatics online.

Rationale: Accurate and comprehensive de novo transcriptome profiling in heart is a central issue to better understand cardiac physiology and diseases. Although significant progress has been made in genome-wide profiling for quantitative changes in cardiac gene expression, current knowledge offers limited insights to the total complexity in cardiac transcriptome at individual exon level. Objective: To develop more robust bioinformatic approaches to analyze high-throughput RNA sequencing (RNA-Seq) data, with the focus on the investigation of transcriptome complexity at individual exon and transcript levels. Methods and Results: In addition to overall gene expression analysis, the methods developed in this study were used to analyze RNA-Seq data with respect to individual transcript isoforms, novel spliced exons, novel alternative terminal exons, novel transcript clusters (ie, novel genes), and long noncoding RNA genes. We applied these approaches to RNA-Seq data obtained from mouse hearts after pressure-overload–induced by transaortic constriction. Based on experimental validations, analyses of the features of the identified exons/transcripts, and expression analyses including previously published RNA-Seq data, we demonstrate that the methods are highly effective in detecting and quantifying individual exons and transcripts. Novel insights inferred from the examined aspects of the cardiac transcriptome open ways to further experimental investigations. Conclusions: Our work provided a comprehensive set of methods to analyze mouse cardiac transcriptome complexity at individual exon and transcript levels. Applications of the methods may infer important new insights to gene regulation in normal and disease hearts in terms of exon utilization and potential involvement of novel components of cardiac transcriptome.

Alternative splicing (AS) has emerged as a key mechanism that accounts for gene expression diversity in metazoan organisms. Splicing is tightly regulated by a repertoire of RNA and protein factors and RNA sequence elements that function in a cooperative manner. Systems-level experimental and computational approaches have been instrumental in establishing comprehensive profiles of transcript variants generated by AS. In addition, systems biology approaches are starting to define how combinatorial splicing regulation shapes the complex splicing phenotypes observed in different tissue types and developmental stages and under different conditions. Here, we review recent progress in these areas. Copyright © 2010 John Wiley & Sons, Inc. For further resources related to this article, please visit the WIREs website

Alternative splicing is a powerful mechanism affording extensive proteomic and regulatory diversity from a limited repertoire of genes. However, the extent to which alternative splicing has contributed to the evolution of primate species-specific characteristics has not been assessed previously. Using comparative genomics and quantitative microarray profiling, we performed the first global analysis of alternative splicing differences between humans and chimpanzees. Surprisingly, 6%–8% of profiled orthologous exons display pronounced splicing level differences in the corresponding tissues from the two species. Little overlap is observed between the genes associated with alternative splicing differences and the genes that display steady-state transcript level differences, indicating that these layers of regulation have evolved rapidly to affect distinct subsets of genes in humans and chimpanzees. The alternative splicing differences we detected are predicted to affect diverse functions including gene expression, signal transduction, cell death, immune defense, and susceptibility to diseases. Differences in expression at the protein level of the major splice variant of GlutathioneS-transferase omega-2 (GSTO2), which functions in the protection against oxidative stress and is associated with human aging-related diseases, suggests that this enzyme is less active in human cells compared with chimpanzee cells. The results of this study thus support an important role for alternative splicing in establishing differences between humans and chimpanzees.

Background: Fetal to neonatal transition of heart is an elaborate process, during which, neonatal cardiomyocytes undergo functional maturation and terminal exit from the cell cycle. However, transcriptome programming in neonatal cardiac chambers during perinatal stages is understudied. In particular, the changes in long non-coding RNAs (lncRNAs) in neonatal heart have not been explored. Objective: To achieve transcriptome-wide analysis of lncRNAs in neonatal left ventricle (LV) and right ventricle (RV) during maturation stages using deep RNA-Sequencing Methods: Deep RNA-sequencing was performed on male newborn mouse (C57 BL) LV and RV at 3 time points of perinatal circulatory transition: P0, P3 and P7. Reads were mapped to mouse genome (mm10). The lncRNAs annotated in NONCODE database were identified. Differentially expressed lncRNAs were defined as those with coefficient of variation ≥0.2, at a false discovery rate ≤0.05, and expressed at ≥3 RPKM in at least one sample. Correlated lncRNAs/ gene pairs were identified using Pearson’s (r2≥0.8, P≤0.05). A subset of LncRNAs/gene expression was validated using qRT-PCR. Results: Out of the 70, 983 observed unique lncRNAs, approximately 7000 were identified exhibiting significant variation during maturation windows with highly spatial-temporal dependent expression patterns, including approximately 5000 known and 2000 novel lncRNAs. Notably, 20% of these lncRNAs were located within 50 KB of a protein coding gene. Out of a total of 2400 lncRNAs/gene pairs, 10 % exhibited significantly concordant (lncRNA/gene) expression patterns. These correlated genes were significantly enriched in metabolism, cell cycle and contractility functional ontology. Interestingly, some of these lncRNAs exhibited concordance with their neighboring gene in human tissues with congenital heart defects, suggesting conserved, potentially significant, regulatory function. Conclusions: Transcriptome programming during neonatal heart maturation involves global changes in lncRNAs. Their expression concordance with neighboring protein coding genes implicates potential important regulatory role of lncRNAs in neonatal heart chamber specification and congenital diseases.

Background: Chamber-specific and temporally orchestrated perinatal cardiac growth and maturation is critical to functional adaptation of heart to the dramatic changes in hemodynamic load and nutrient environment, and perturbation of this process may have major implications in congenital heart defects (CHDs). However, the responsible molecular regulatory mechanisms are still poorly understood. Objective: To determine the factors that control chamber specificity of neonatal heart growth during perinatal maturation Methods: Total transcriptome was analyzed from deep RNA-seq in male newborn mouse (C57BL) left ventricle (LV) and right ventricle (RV) at postnatal day 0, 3 and 7. Functional validation was performed in Neonatal Rat Ventricular Myocytes (NRVMs) and intact neonatal mice. Results: Chamber specific- regulation of cell cycle genes was correlated with chamber specific- expression of the non-canonical Wnt11. Exogenous expression of Wnt11 induced maturation markers and suppressed cell cycle genes in myocytes. In contrast, targeted Wnt11 inhibition induced myocyte proliferation associated with upregulated cell cycle genes both in cultured NRVMs and in the intact heart. Notably, Wnt inhibition also significantly enhanced β-catenin phosphorylation without impacting JNK, PKCα, or p38- activities. Remarkably, both Wnt mediated signals and cell cycle markers were altered in heart samples from infants cyanotic CHDs. Furthermore, hypoxia exposure suppressed Wnt11 and induced cell cycle genes in cultured cardiomyocytes. Conclusions: Differential expression of Wnt11 and downstream cell cycle genes are novel molecular basis for chamber specific maturation. Wnt11 mediated-signaling has potential impact on cyanotic CHDs. Differential role of Wnt signaling in chamber specific post-natal heart maturation can lead to potential chamber-specific diagnoses and therapies for congenital heart diseases.

Pre-mRNA splicing is an essential RNA processing step in eukaryotes. Alternative splicing generates distinct spliced isoforms of the same gene, thereby dramatically increasing transcriptome diversity. Since most human genes undergo alternative splicing, this process contributes to a wide spectrum of biological functions in healthy and disease states. Splicing is closely regulated by various cis-regulatory elements and trans-factors. With the advent of high-throughput experimental technologies and bioinformatic algorithms, we now have powerful means to study alternative splicing globally and uncover its functional impact and regulatory mechanisms. As more RNA sequencing (RNA-Seq) data from normal and disease conditions are becoming available, many studies are underway to dissect global misregulation of splicing in diseases and develop novel splicing-targeted therapeutics. In this chapter, we first discuss the experimental and bioinformatic approaches for identification of alternative splicing, followed by a comprehensive review on the state-of-the-art methodologies to study splicing regulation. In addition, we discuss the current challenges and open questions in the RNA splicing field including gene expression kinetics, co-transcriptional splicing, and therapeutic approaches targeting splicing.