Search

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.

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

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.

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

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

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

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

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

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

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.

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.

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

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.

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

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

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.

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

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.

Background: There are numerous anecdotal reports of ventricular arrhythmias during spaceflight; however, it is not known whether spaceflight or microgravity systematically increases the risk of cardiac dysrhythmias. Microvolt T wave alternans (MTWA) testing compares favorably with other noninvasive risk stratifiers and invasive electrophysiological testing in patients as a predictor of sudden cardiac death, ventricular tachycardia, and ventricular fibrillation. We hypothesized that simulated microgravity leads to an increase in MTWA. Methods: Twenty-four healthy male subjects underwent 9 to 16 days of head-down tilt bed rest (HDTB). MTWA was measured before and after the bed rest period during bicycle exercise stress. For the purposes of this study, we defined MTWA outcome to be positive if sustained MTWA was present with an onset heart rate ≤125 bpm. During various phases of HDTB, the following were also performed: daily 24-hour urine collections, serum electrolytes and catecholamines, and cardiovascular system identification (measure of autonomic function). Results: Before HDTB, 17% of the subjects were MTWA positive [95%CI: (0.6%, 37%)]; after HDTB, 42% of the subjects were MTWA positive [95%CI: (23%, 63%)] (P = 0.03). The subjects who were MTWA positive after HDTB compared with MTWA negative subjects had an increased versus decreased sympathetic responsiveness (P = 0.03) and serum norepinephrine levels (P = 0.05), and a trend toward higher potassium excretion (P = 0.06) after bed rest compared to baseline. Conclusions: HDTB leads to an increase in MTWA, providing the first evidence that simulated microgravity has a measurable effect on electrical repolarization processes. Possible contributing factors include loss in potassium and changes in sympathetic function.

We introduce a model order selection criterion called signal prediction error (SPE) for the identification of a linear regression model, which can be an adequate representation of a resting physiologic system. SPE is an estimate of the prediction error variance due only to model estimation error and not unobserved noise, which distinguishes it from the widely used final prediction error (FPE). We then present a theoretical analysis of SPE, which predicts that its ability to select correctly the model order is more dependent on the signal-to-noise ratio (SNR) and less dependent on the number of data samples available for analysis. We next propose a heuristic procedure based on SPE (called SPE/sup D/) to improve its robustness to SNR levels. We then demonstrate, through simulated physiologic data at high SNR levels, that SPE will be equivalent to consistent model order selection criteria for long data records but will become superior to FPE and other model order selection criteria as the size of the data record decreases. The simulated data results also show that SPE/sup D/ is indeed a significant improvement over SPE in terms of robustness to SNR. Finally, we demonstrate the applicability of SPE and SPE/sup D/ to actual cardio-respiratory-renal data.

Orthostatic intolerance (OI) is a major problem after spaceflight. Its etiology remains uncertain, but reports have pointed toward an individual susceptibility to OI. We hypothesized that individual predisposition plays an important role in post-bed rest OI. Twenty-four healthy male subjects were equilibrated on a constant diet, after which they underwent tilt-stand test (pre-TST). They then completed 14-16 days of head-down-tilt bed rest, and 14 of the subjects underwent repeat tilt-stand test (post-TST). During various phases, the following were performed: 24-h urine collections and hormonal measurements, plethysmography, and cardiovascular system identification (a noninvasive method to assess autonomic function and separately quantify parasympathetic and sympathetic responsiveness). Development of presyncope or syncope defined OI. During pre-TST, 11 subjects were intolerant and 13 were tolerant. At baseline, intolerant subjects had lower serum aldosterone ( P < 0.01), higher excretion of potassium ( P = 0.01), lower leg venous compliance ( P = 0.03), higher supine parasympathetic responsiveness ( P = 0.02), and lower standing sympathetic responsiveness ( P = 0.048). Of the 14 subjects who completed post-TST, 9 were intolerant and 5 were tolerant. Intolerant subjects had lower baseline serum cortisol ( P = 0.03) and a higher sodium level ( P = 0.02) compared with tolerant subjects. Thus several physiological characteristics were associated with increased susceptibility to OI. We propose a new model for OI, whereby individuals with greater leg venous compliance recruit compensatory mechanisms (activation of the renin-angiotensin-aldosterone system and sympathetic nervous system, and withdrawal of the parasympathetic nervous system) in the face of daily postural challenges, which places them at an advantage to face orthostatic stress. With head-down-tilt bed rest, the stimulus to recruit compensatory mechanisms disappears, and differences between the two subgroups attenuate.

Background: During postnatal maturation of mammalian heart the neonatal cardiomyocytes undergo dramatic changes including complete maturation, proliferation arrest, and terminal exit from the cell cycle (CC). However, transcriptome-wide analysis of CC programs has not been performed in perinatal stages among different cardiac chambers. In particular, the contribution of alternative RNA splicing to the chamber-specific CC activities is unexplored Design/Methods: To achieve comprehensive analysis of differential expression (DE) and alternative splicing (AS) of CC-related genes in left ventricle (LV) versus right ventricle (RV) during maturation deep RNA-seq was performed on male newborn mouse LV and RV at 3 time points of perinatal transition: P0, P3 and P7. Reads were mapped to mouse Transcriptome, and to mouse Genome. Transcriptome-Wide difference in inclusion of individual exons was performed using MATS. DE genes and AS variants were defined as those with fold change ≥2, at expression level ≥3 RBKM and a false discovery rate ≤0.05. Significant gene ontology (GO) terms were determined at P-value ≤0.05. Levels of expression were validated using qRT-PCR Results: Altogether, 2116 DE genes and 1162 AS events were observed. Among them, 109 CC-related genes were further analyzed. Distinct temporal patterns of DE and GO enrichment of CC genes in LV vs. RV during maturation were identified. Chamber -specific induction of genes involved in mitosis, karyokinesis, and cytokinesis was found at P7. RNA Splicing analysis of CC genes revealed 77 AS events. Skipping exon accounted for nearly half splicing events. Among 30 spliced exon variants, significant chamber-and temporal-specific inclusion were observed. Interestingly, the majority of AS variants exhibited opposing patterns of exon usage in RV vs. LV at p7 Conclusions: Our findings suggest novel molecular basis for chamber-specific programming of cellular proliferation and maturation in neonatal heart, including potential splicing regulation of dynamic exon enrichment of cell cycle related genes. Further functional studies to decipher putative splicing regulators of CC programming in LV vs. RV during maturation will likely lead to novel chamber-specific regenerative and therapeutic targets.

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.