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.

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.

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.

Lipin 1 regulates cellular lipid homeostasis through roles in glycerolipid synthesis (through phosphatidic acid phosphatase activity) and transcriptional coactivation. Lipin 1-deficient individuals exhibit episodic disease symptoms that are triggered by metabolic stress, such as stress caused by prolonged fasting. We sought to identify critical lipin 1 activities during fasting. We determined that lipin 1 deficiency induces widespread alternative mRNA splicing in liver during fasting, much of which is normalized by refeeding. The role of lipin 1 in mRNA splicing was largely independent of its enzymatic function. We identified interactions between lipin 1 and spliceosome proteins, as well as a requirement for lipin 1 to maintain homeostatic levels of spliceosome small nuclear RNAs and specific RNA splicing factors. In fasted Lpin1-/- liver, we identified a correspondence between alternative splicing of phospholipid biosynthetic enzymes and dysregulated phospholipid levels; splicing patterns and phospholipid levels were partly normalized by feeding. Thus, lipin 1 influences hepatic lipid metabolism through mRNA splicing, as well as through enzymatic and transcriptional activities, and fasting exacerbates the deleterious effects of lipin 1 deficiency on metabolic homeostasis.

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.

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.

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.

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.

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.

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.

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.

© 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.

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.

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.

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.

Ventricular chamber growth and development during perinatal circulatory transition is critical for functional adaptation of the heart. However, the chamber-specific programs of neonatal heart growth are poorly understood. We used integrated systems genomic and functional biology analyses of the perinatal chamber specific transcriptome and we identified Wnt11 as a prominent regulator of chamber-specific proliferation. Importantly, downregulation of Wnt11 expression was associated with cyanotic congenital heart defect (CHD) phenotypes and correlated with O2 saturation levels in hypoxemic infants with Tetralogy of Fallot (TOF). Perinatal hypoxia treatment in mice suppressed Wnt11 expression and induced myocyte proliferation more robustly in the right ventricle, modulating Rb1 protein activity. Wnt11 inactivation was sufficient to induce myocyte proliferation in perinatal mouse hearts and reduced Rb1 protein and phosphorylation in neonatal cardiomyocytes. Finally, downregulated Wnt11 in hypoxemic TOF infantile hearts was associated with Rb1 suppression and induction of proliferation markers. This study revealed a previously uncharacterized function of Wnt11-mediated signaling as an important player in programming the chamber-specific growth of the neonatal heart. This function influences the chamber-specific development and pathogenesis in response to hypoxia and cyanotic CHDs. Defining the underlying regulatory mechanism may yield chamber-specific therapies for infants born with CHDs.

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.

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.

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.

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

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.

Although major genetic networks controlling early liver specification and morphogenesis are known, the mechanisms responsible for postnatal hepatic maturation are poorly understood. Here we employ global analyses of the mouse liver transcriptome to demonstrate that postnatal remodelling of the liver is accompanied by large-scale transcriptional and post-transcriptional transitions that are cell-type-specific and temporally coordinated. Combining detailed expression analyses with gain- and loss-of-function studies, we identify epithelial splicing regulatory protein 2 (ESRP2) as a conserved regulatory factor that controls the neonatal-to-adult switch of ∼20% of splice isoforms in mouse and human hepatocytes. The normal shift in splicing coincides tightly with dramatic postnatal induction of ESRP2 in hepatocytes. We further demonstrate that forced expression of ESRP2 in immature mouse and human hepatocytes is sufficient to drive a reciprocal shift in splicing and causes various physiological abnormalities. These findings define a direct role for ESRP2 in the generation of conserved repertoires of adult splice isoforms that facilitate terminal differentiation and maturation of hepatocytes., Alternative RNA splicing is important during organismal development. Here, the authors perform RNA-Seq on mouse and human liver samples to provide a comprehensive view of splicing events during liver development and growth, and identify Espr2 as a main regulator of these splicing processes.

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.

Human saliva emerged as a research material in as early as the 17th century when investigators sought to understand the basis for salivary secretion (Garrett, 1975). Over the centuries, the focus of salivary research has evolved greatly and a wide range of topics has been examined (Garrett, 1975; Schipper et al., 2007). It is now known that the functions of saliva include at least lubrication, digestion of food, remineralization, prevention of demineralization, protection against microbial and viral infection, speech facilitation, and maintenance of oral and general health (Schipper et al., 2007; Malathi et al., 2014).

Correct splice site recognition is critical in pre-mRNA splicing. We find that almost all of a diverse panel of exonic splicing silencer (ESS) elements alter splice site choice when placed between competing sites, consistently inhibiting use of intron-proximal 5' and 3' splice sites. Supporting a general role for ESSs in splice site definition, we found that ESSs are both abundant and highly conserved between alternative splice site pairs and that mutation of ESSs located between natural alternative splice site pairs consistently shifted splicing toward the intron-proximal site. Some exonic splicing enhancers (ESEs) promoted use of intron-proximal 5' splice sites, and tethering of hnRNP A1 and SF2/ASF proteins between competing splice sites mimicked the effects of ESS and ESE elements, respectively. Further, we observed that specific subsets of ESSs had distinct effects on a multifunctional intron retention reporter and that one of these subsets is likely preferred for regulation of endogenous intron retention events. Together, our findings provide a comprehensive picture of the functions of ESSs in the control of diverse types of splicing decisions.

BACKGROUND Extracellular RNAs (exRNAs) in human body fluids are emerging as effective biomarkers for detection of diseases. Saliva, as the most accessible and noninvasive body fluid, has been shown to harbor exRNA biomarkers for several human diseases. However, the entire spectrum of exRNA from saliva has not been fully characterized. METHODS Using high-throughput RNA sequencing (RNA-Seq), we conducted an in-depth bioinformatic analysis of noncoding RNAs (ncRNAs) in human cell-free saliva (CFS) from healthy individuals, with a focus on microRNAs (miRNAs), piwi-interacting RNAs (piRNAs), and circular RNAs (circRNAs). RESULTS Our data demonstrated robust reproducibility of miRNA and piRNA profiles across individuals. Furthermore, individual variability of these salivary RNA species was highly similar to those in other body fluids or cellular samples, despite the direct exposure of saliva to environmental impacts. By comparative analysis of >90 RNA-Seq data sets of different origins, we observed that piRNAs were surprisingly abundant in CFS compared with other body fluid or intracellular samples, with expression levels in CFS comparable to those found in embryonic stem cells and skin cells. Conversely, miRNA expression profiles in CFS were highly similar to those in serum and cerebrospinal fluid. Using a customized bioinformatics method, we identified >400 circRNAs in CFS. These data represent the first global characterization and experimental validation of circRNAs in any type of extracellular body fluid. CONCLUSIONS Our study provides a comprehensive landscape of ncRNA species in human saliva that will facilitate further biomarker discoveries and lay a foundation for future studies related to ncRNAs in human saliva.

Tumor metastasis remains the major cause of cancer-related death, but its molecular basis is still not well understood. Here we uncovered a splicing-mediated pathway that is essential for breast cancer metastasis. We show that the RNA-binding protein heterogeneous nuclear ribonucleoprotein M (hnRNPM) promotes breast cancer metastasis by activating the switch of alternative splicing that occurs during epithelial–mesenchymal transition (EMT). Genome-wide deep sequencing analysis suggests that hnRNPM potentiates TGFβ signaling and identifies CD44 as a key downstream target of hnRNPM. hnRNPM ablation prevents TGFβ-induced EMT and inhibits breast cancer metastasis in mice, whereas enforced expression of the specific CD44 standard (CD44s) splice isoform overrides the loss of hnRNPM and permits EMT and metastasis. Mechanistically, we demonstrate that the ubiquitously expressed hnRNPM acts in a mesenchymal-specific manner to precisely control CD44 splice isoform switching during EMT. This restricted cell-type activity of hnRNPM is achieved by competition with ESRP1, an epithelial splicing regulator that binds to the same cis-regulatory RNA elements as hnRNPM and is repressed during EMT. Importantly, hnRNPM is associated with aggressive breast cancer and correlates with increased CD44s in patient specimens. These findings demonstrate a novel molecular mechanism through which tumor metastasis is endowed by the hnRNPM-mediated splicing program.