Search

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.

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.

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.

In the mammalian heart, the left ventricle (LV) rapidly becomes more dominant in size and function over the right ventricle (RV) after birth. The molecular regulators responsible for this chamber-specific differential growth are largely unknown. We found that cardiomyocytes in the neonatal mouse RV had lower proliferation, more apoptosis, and a smaller average size compared with the LV. This chamber-specific growth pattern was associated with a selective activation of p38 mitogen-activated protein kinase (MAPK) activity in the RV and simultaneous inactivation in the LV. Cardiomyocyte-specific deletion of both the Mapk14 and Mapk11 genes in mice resulted in loss of p38 MAPK expression and activity in the neonatal heart. Inactivation of p38 activity led to a marked increase in cardiomyocyte proliferation and hypertrophy but diminished cardiomyocyte apoptosis, specifically in the RV. Consequently, the p38-inactivated hearts showed RV-specific enlargement postnatally, progressing to pulmonary hypertension and right heart failure at the adult stage. Chamber-specific p38 activity was associated with differential expression of dual-specific phosphatases (DUSPs) in neonatal hearts, including DUSP26. Unbiased transcriptome analysis revealed that IRE1α/XBP1-mediated gene regulation contributed to p38 MAPK-dependent regulation of neonatal cardiomyocyte proliferation and binucleation. These findings establish an obligatory role of DUSP/p38/IRE1α signaling in cardiomyocytes for chamber-specific growth in the postnatal heart.

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 datasets 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 vs. 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.

microRNAs (miRNAs) are small non-coding RNAs that play critical roles in gene regulation. The presence of miRNAs in extracellular biofluids is increasingly recognized. However, most previous characterization of extracellular miRNAs focused on their overall expression levels. Alternative sequence isoforms and modifications of miRNAs were rarely considered in the extracellular space. Here, we developed a highly accurate bioinformatic method, called miNTA, to identify 3’ non-templated additions (NTAs) of miRNAs using small RNA-sequencing data. Using miNTA, we conducted an in-depth analysis of miRNA 3’ NTA profiles in 1047 extracellular RNA-sequencing data sets of 4 types of biofluids. This analysis identified hundreds of miRNAs with 3’ uridylation or adenylation, with the former being more prevalent. Among these miRNAs, up to 53% (22%) had an average 3’ uridylation (adenylation) level of at least 10% in a specific biofluid. Strikingly, we found that 3’ uridylation levels enabled segregation of different types of biofluids, more effectively than overall miRNA expression levels. This observation suggests that 3’ NTA levels possess fluid-specific information relatively robust to batch effects. In addition, we observed that extracellular miRNAs with 3’ uridylations are enriched in processes related to angiogenesis, apoptosis, and inflammatory response, and this type of modification may stabilize base-pairing between miRNAs and their target genes. Together, our study provides a comprehensive landscape of miRNA NTAs in human biofluids, which paves way for further biomarker discoveries. The insights generated in our work built a foundation for future functional, mechanistic, and translational discoveries.

Recent studies revealed global shifts in RNA editing, the modification of RNA sequences, across many cancers. Besides a few sites implicated in tumorigenesis or metastasis, most tumor-associated sites, predominantly in noncoding regions, have unknown function. Here, we characterize editing profiles between epithelial (E) and mesenchymal (M) phenotypes in seven cancer types, as epithelial-mesenchymal transition (EMT) is a key paradigm for metastasis. We observe distinct editing patterns between E and M tumors and EMT induction upon loss of ADAR enzymes in cultured cells. E-M differential sites are highly enriched in genes involved in immune and viral processes, some of which regulate mRNA abundance of their respective genes. We identify a novel mechanism in which ILF3 preferentially stabilizes edited transcripts. Among editing-dependent ILF3 targets is the transcript encoding PKR, a crucial player in immune response. Our study demonstrates the broad impact of RNA editing in cancer and relevance of editing to cancer-related immune pathways.

The essential biological cofactor heme is synthesized in cells in the Fe(II) form. Oxidized Fe(III) heme is specifically required for processing primary transcripts of microRNAs (pri-miRNAs) by the RNA-binding protein DGCR8, a core component of the Microprocessor complex. It is unknown how readily available Fe(III) heme is in the largely reducing environment in human cells and how changes in cellular Fe(III) heme availability alter microRNA (miRNA) expression. Here we address the first question by characterizing DGCR8 mutants with various degrees of deficiency in heme-binding. We observed a strikingly simple correlation between Fe(III) heme affinityin vitroand the Microprocessor activity in HeLa cells, with the heme affinity threshold for activation estimated to be between 0.6-5 pM under typical cell culture conditions. The threshold is strongly influenced by cellular heme synthesis and uptake. We suggest that the threshold reflects a labile Fe(III) heme pool in cells. Based on our understanding of DGCR8 mutants, we reanalyzed recently reported miRNA sequencing data and conclude that heme is generally required for processing canonical pri-miRNAs, that heme modulates the specificity of Microprocessor, and that cellular heme level and differential DGCR8 heme occupancy alter the expression of distinct groups of miRNAs in a hierarchical fashion. Overall, our study provides the first glimpse of a labile Fe(III) heme pool important for a fundamental physiological function and reveal principles governing how Fe(III) heme modulates miRNA maturation at a genomic scale. We also discuss potential states and biological significance of the labile Fe(III) heme pool.

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.

RNA editing introduces single nucleotide changes to RNA, thus potentially diversifying gene expression. Recent studies have reported significant changes in RNA editing profiles in disease and development. The functional consequences of these widespread alterations remain elusive because of the unknown function of most RNA editing sites. Here, we carry out a comprehensive analysis of A-to-I editomes in human populations. Surprisingly, we observe highly similar editing profiles across populations despite striking differences in the expression levels of ADAR genes. Striving to explain this discrepancy, we uncover a functional mechanism of A-to-I editing in regulating mRNA abundance. We show that A-to-I editing stabilizes RNA secondary structures and reduces the accessibility of AGO2-miRNA to target sites in mRNAs. The editing-dependent stabilization of mRNAs in turn alters the observed editing levels in the stable RNA repertoire. Our study provides valuable insights into the functional impact of RNA editing in human cells., RNA A-to-I editing introduces single nucleotide changes to RNA, but its role in cells remains unclear. Here, the authors analyse A-to-I editomes in human samples and find that A-to-I editing stabilizes RNA secondary structures and reduces the accessibility of AGO2-miRNA to target sites in mRNAs.

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.

© 2019 Elsevier Inc. The Extracellular RNA Communication Consortium (ERCC) was launched to accelerate progress in the new field of extracellular RNA (exRNA) biology and to establish whether exRNAs and their carriers, including extracellular vesicles (EVs), can mediate intercellular communication and be utilized for clinical applications. Phase 1 of the ERCC focused on exRNA/EV biogenesis and function, discovery of exRNA biomarkers, development of exRNA/EV-based therapeutics, and construction of a robust set of reference exRNA profiles for a variety of biofluids. Here, we present progress by ERCC investigators in these areas, and we discuss collaborative projects directed at development of robust methods for EV/exRNA isolation and analysis and tools for sharing and computational analysis of exRNA profiling data. The Extracellular RNA Communication Consortium (ERCC) presents progress toward understanding the biology of extracellular RNAs and their use as biomarkers and therapeutics.

Allele-specific protein-RNA binding is an essential aspect that may reveal functional genetic variants (GVs) mediating post-transcriptional regulation. Recently, genome-wide detection of in vivo binding of RNA-binding proteins is greatly facilitated by the enhanced crosslinking and immunoprecipitation (eCLIP) method. We developed a new computational approach, called BEAPR, to identify allele-specific binding (ASB) events in eCLIP-Seq data. BEAPR takes into account crosslinking-induced sequence propensity and variations between replicated experiments. Using simulated and actual data, we show that BEAPR largely outperforms often-used count analysis methods. Importantly, BEAPR overcomes the inherent overdispersion problem of these methods. Complemented by experimental validations, we demonstrate that the application of BEAPR to ENCODE eCLIP-Seq data of 154 proteins helps to predict functional GVs that alter splicing or mRNA abundance. Moreover, many GVs with ASB patterns have known disease relevance. Overall, BEAPR is an effective method that helps to address the outstanding challenge of functional interpretation of GVs., Differential binding of RNA-binding proteins mediated by genetic variants (GVs) can influence posttranscriptional regulation. Here, the authors develop BEAPR, a computational approach to identify allele-specific binding events in eCLIP-Seq data.

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.

Autism spectrum disorder (ASD) is a genetically complex, clinically heterogeneous neurodevelopmental disease. Recently, our understanding of the molecular abnormalities in ASD has been expanded through transcriptomic analyses of postmortem brains. However, a crucial molecular pathway involved in synaptic development, RNA editing, has not yet been studied on a genome-wide scale. Here, we profiled the global patterns of adenosine-to-inosine (A-to-I) editing in a large cohort of post-mortem ASD brains. Strikingly, we observed a global bias of hypo-editing in ASD brains, common to different brain regions and involving many genes with known neurobiological functions. Through genome-wide protein-RNA binding analyses and detailed molecular assays, we show that the Fragile X proteins, FMRP and FXR1P, interact with ADAR proteins and modulate A-to-I editing. Furthermore, we observed convergent patterns of RNA editing alterations in ASD and Fragile X syndrome, thus establishing RNA editing as a molecular link underlying these two highly related diseases. Our findings support a role for RNA editing dysregulation in ASD and highlight novel mechanisms for RNA editing regulation.

© International & American Associations for Dental Research 2018. Coevolution of the human host and its associated microbiota has led to sophisticated interactions to maintain a delicate homeostasis. Emerging evidence suggests that in addition to small molecules, peptides, and proteins, small regulatory noncoding RNAs (sRNAs) might play an important role in cross-domain interactions. In this study, we revealed the presence of diverse host transfer RNA–derived small RNAs (tsRNAs) among human salivary sRNAs. We selected 2 tsRNAs (tsRNA-000794 and tsRNA-020498) for further study based on their high sequence similarity to specific tRNAs from a group of Gram-negative oral bacteria, including Fusobacterium nucleatum, a key oral commensal and opportunistic pathogen. We showed that the presence of F. nucleatum triggers exosome-mediated release of tsRNA-000794 and tsRNA-020498 by human normal oral keratinocyte cells. Furthermore, both tsRNA candidates exerted a growth inhibition effect on F. nucleatum, likely through interference with bacterial protein biosynthesis, but did not affect the growth of Streptococcus mitis, a health-associated oral Gram-positive bacterium whose genome does not carry sequences bearing high similarity to either tsRNA. Our data provide the first line of evidence for the modulatory role of host-derived tsRNAs in the microbial-host interaction.

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.

Human speech is one of the few examples of vocal learning among mammals yet ~half of avian species exhibit this ability. Its neurogenetic basis is largely unknown beyond a shared requirement for FoxP2 in both humans and zebra finches. We manipulated FoxP2 isoforms in Area X, a song-specific region of the avian striatopallidum analogous to human anterior striatum, during a critical period for song development. We delineate, for the first time, unique contributions of each isoform to vocal learning. Weighted gene coexpression network analysis of RNA-seq data revealed gene modules correlated to singing, learning, or vocal variability. Coexpression related to singing was found in juvenile and adult Area X whereas coexpression correlated to learning was unique to juveniles. The confluence of learning and singing coexpression in juvenile Area X may underscore molecular processes that drive vocal learning in young zebra finches and, by analogy, humans., eLife digest Songbirds, much like in humans, have a critical period in youth when they are best at learning vocal communication skills. In birds, this is when they learn a song they will use later in life as a courtship song. In humans, this is when language skills are most easily learned. After this critical period ends, it is much harder for people to learn languages, and for certain bird species to learn their song. When birds sing every morning, the activity of a gene called FoxP2 drops, which causes a coordinated change in the activity of thousands of other genes. It is suspected that FoxP2 – and the changes it causes – could be a part of the molecular basis for vocal learning. FoxP2 is also known to play a role in speech in humans, and both birds and humans have a long and a short version of this gene. Previous research has shown that when the long version of the gene was altered so its activity would no longer decrease when birds were singing, the birds failed to learn their song. Moreover, humans with a mutation in the long version have problems with their speech. However, until now, it was not known if modifications to the short version had the same effect. Burkett et al. investigated whether there was a noticeable pattern in the effects of FoxP2 before and after the critical period in a songbird. The analysis found that during the critical period, a set of genes changed together as young birds learned to sing. This particular pattern disappeared as the birds aged and the critical period ended. Burkett et al. confirmed that when birds had the long version of FoxP2 altered, they were less able to learn. However, changing the short version of FoxP2 had little effect on learning but led to changes in the birds’ song. The genetic pathways identified in the experiments are known to be present in many different species, including humans. Related pathways have also been found to play a role in non-vocal learning in organisms as distantly related as rats and snails. This suggests that they could be acting as a blueprint for learning new skills. Few treatments for language impairments have been developed so far due to poor understanding of the molecular basis for vocal communication. The findings of this study could help to create new treatments for speech problems in people, such as children with autism or people with mutated versions of FoxP2.

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.

Exposure of endothelial cells (ECs) to agents such as oxidized glycerophospholipids (oxGPs) and cytokines, known to accumulate in atherosclerotic lesions, perturbs the expression of hundreds of genes in ECs involved in inflammatory and other biological processes. We hypothesized that microRNAs (miRNAs) are involved in regulating the inflammatory response in human aortic endothelial cells (HAECs) in response to oxGPs and interleukin 1β (IL-1β). Using next-generation sequencing and RT-quantitative PCR, we characterized the profile of expressed miRNAs in HAECs pre- and postexposure to oxGPs. Using this data, we identified miR-21-3p and miR-27a-5p to be induced 3- to 4-fold in response to oxGP and IL-1β treatment compared with control treatment. Transient overexpression of miR-21-3p and miR-27a-5p resulted in the downregulation of 1,253 genes with 922 genes overlapping between the two miRNAs. Gene Ontology functional enrichment analysis predicted that the two miRNAs were involved in the regulation of nuclear factor κB (NF-κB) signaling. Overexpression of these two miRNAs leads to changes in p65 nuclear translocation. Using 3′ untranslated region luciferase assay, we identified 20 genes within the NF-κB signaling cascade as putative targets of miRs-21-3p and -27a-5p, implicating these two miRNAs as modulators of NF-κB signaling in ECs.

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.

Background— Although metabolic reprogramming is critical in the pathogenesis of heart failure, studies to date have focused principally on fatty acid and glucose metabolism. Contribution of amino acid metabolic regulation in the disease remains understudied. Methods and Results— Transcriptomic and metabolomic analyses were performed in mouse failing heart induced by pressure overload. Suppression of branched-chain amino acid (BCAA) catabolic gene expression along with concomitant tissue accumulation of branched-chain α-keto acids was identified as a significant signature of metabolic reprogramming in mouse failing hearts and validated to be shared in human cardiomyopathy hearts. Molecular and genetic evidence identified the transcription factor Krüppel-like factor 15 as a key upstream regulator of the BCAA catabolic regulation in the heart. Studies using a genetic mouse model revealed that BCAA catabolic defect promoted heart failure associated with induced oxidative stress and metabolic disturbance in response to mechanical overload. Mechanistically, elevated branched-chain α-keto acids directly suppressed respiration and induced superoxide production in isolated mitochondria. Finally, pharmacological enhancement of branched-chain α-keto acid dehydrogenase activity significantly blunted cardiac dysfunction after pressure overload. Conclusions— BCAA catabolic defect is a metabolic hallmark of failing heart resulting from Krüppel-like factor 15–mediated transcriptional reprogramming. BCAA catabolic defect imposes a previously unappreciated significant contribution to heart failure.

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.

Summary Most human genes produce multiple splicing isoforms with distinct functions. To systematically understand splicing regulation, we conducted an unbiased screen and identified >100 intronic splicing enhancers (ISEs) that were clustered by sequence similarity into six groups. All ISEs functioned in another cell type and heterologous introns, and their distribution and conservation patterns in different pre-mRNA regions are similar to exonic splicing silencers. Consistently all ISEs inhibited use of splice sites from exonic locations. The putative trans-factors of each ISE group were identified and validated. Five distinct ISE motifs were recognized by hnRNP H and F whose C-terminal domains were sufficient to render context-dependent activities of ISEs. The sixth group was controlled by factors that either activate or suppress splicing. This work provided a comprehensive picture of general ISE activities and provided new models of how a single element can function oppositely depending on its locations and binding factors.

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

Pre-mRNA splicing is regulated through the combinatorial activity of RNA motifs, including splice sites and splicing regulatory elements. Here we show that the activity of the G-run (polyguanine sequence) class of splicing enhancer elements is approx4-fold higher when adjacent to intermediate strength 5' splice sites (ss) than when adjacent to weak 5' ss, and approx1.3-fold higher relative to strong 5' ss. We observed this dependence on 5' ss strength in both splicing reporters and in global microarray and mRNA-Seq analyses of splicing changes following RNA interference against heterogeneous nuclear ribonucleoprotein (hnRNP) H, which cross-linked to G-runs adjacent to many regulated exons. An exon's responsiveness to changes in hnRNP H levels therefore depends in a complex way on G-run abundance and 5' ss strength. This pattern of activity enables G-runs and hnRNP H to buffer the effects of 5' ss mutations, augmenting both the frequency of 5' ss polymorphism and the evolution of new splicing patterns. Certain other splicing factors may function similarly., American Heart Association, Human Frontier Science Program (Strasbourg, France), National Institutes of Health (U.S.), National Science Foundation (U.S.) (equipment grant DBI-0821391)

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.