Dynamic ASXL1 Exon Skipping and Alternative Circular Splicing in Single Human Cells.

Circular RNAs comprise a poorly understood new class of noncoding RNA. In this study, we used a combination of targeted deletion, high-resolution splicing detection, and single-cell sequencing to deeply probe ASXL1 circular splicing. We found that efficient circular splicing required the canonical transcriptional start site and inverted AluSx elements. Sequencing-based interrogation of isoforms after ASXL1 overexpression identified promiscuous linear splicing between all exons, with the two most abundant non-canonical linear products skipping the exons that produced the circular isoforms. Single-cell sequencing revealed a strong preference for either the linear or circular ASXL1 isoforms in each cell, and found the predominant exon skipping product is frequently co-expressed with its reciprocal circular isoform. Finally, absolute quantification of ASXL1 isoforms confirmed our findings and suggests that standard methods overestimate circRNA abundance. Taken together, these data reveal a dynamic new view of circRNA genesis, providing additional framework for studying their roles in cellular biology.

Directed Chemical Evolution with an Outsized Genetic Code.

The first demonstration that macromolecules could be evolved in a test tube was reported twenty-five years ago. That breakthrough meant that billions of years of chance discovery and refinement could be compressed into a few weeks, and provided a powerful tool that now dominates all aspects of protein engineering. A challenge has been to extend this scientific advance into synthetic chemical space: to enable the directed evolution of abiotic molecules. The problem has been tackled in many ways. These include expanding the natural genetic code to include unnatural amino acids, engineering polyketide and polypeptide synthases to produce novel products, and tagging combinatorial chemistry libraries with DNA. Importantly, there is still no small-molecule analog of directed protein evolution, i.e. a substantiated approach for optimizing complex (≥ 10^9 diversity) populations of synthetic small molecules over successive generations. We present a key advance towards this goal: a tool for genetically-programmed synthesis of small-molecule libraries from large chemical alphabets. The approach accommodates alphabets that are one to two orders of magnitude larger than any in Nature, and facilitates evolution within the chemical spaces they create. This is critical for small molecules, which are built up from numerous and highly varied chemical fragments. We report a proof-of-concept chemical evolution experiment utilizing an outsized genetic code, and demonstrate that fitness traits can be passed from an initial small-molecule population through to the great-grandchildren of that population. The results establish the practical feasibility of engineering synthetic small molecules through accelerated evolution.

Evolutionary Conservation and Diversification of Puf RNA Binding Proteins and Their mRNA Targets.

Reprogramming of a gene's expression pattern by acquisition and loss of sequences recognized by specific regulatory RNA binding proteins may be a major mechanism in the evolution of biological regulatory programs. We identified that RNA targets of Puf3 orthologs have been conserved over 100-500 million years of evolution in five eukaryotic lineages. Focusing on Puf proteins and their targets across 80 fungi, we constructed a parsimonious model for their evolutionary history. This model entails extensive and coordinated changes in the Puf targets as well as changes in the number of Puf genes and alterations of RNA binding specificity including that: 1) Binding of Puf3 to more than 200 RNAs whose protein products are predominantly involved in the production and organization of mitochondrial complexes predates the origin of budding yeasts and filamentous fungi and was maintained for 500 million years, throughout the evolution of budding yeast. 2) In filamentous fungi, remarkably, more than 150 of the ancestral Puf3 targets were gained by Puf4, with one lineage maintaining both Puf3 and Puf4 as regulators and a sister lineage losing Puf3 as a regulator of these RNAs. The decrease in gene expression of these mRNAs upon deletion of Puf4 in filamentous fungi (N. crassa) in contrast to the increase upon Puf3 deletion in budding yeast (S. cerevisiae) suggests that the output of the RNA regulatory network is different with Puf4 in filamentous fungi than with Puf3 in budding yeast. 3) The coregulated Puf4 target set in filamentous fungi expanded to include mitochondrial genes involved in the tricarboxylic acid (TCA) cycle and other nuclear-encoded RNAs with mitochondrial function not bound by Puf3 in budding yeast, observations that provide additional evidence for substantial rewiring of post-transcriptional regulation. 4) Puf3 also expanded and diversified its targets in filamentous fungi, gaining interactions with the mRNAs encoding the mitochondrial electron transport chain (ETC) complex I as well as hundreds of other mRNAs with nonmitochondrial functions. The many concerted and conserved changes in the RNA targets of Puf proteins strongly support an extensive role of RNA binding proteins in coordinating gene expression, as originally proposed by Keene. Rewiring of Puf-coordinated mRNA targets and transcriptional control of the same genes occurred at different points in evolution, suggesting that there have been distinct adaptations via RNA binding proteins and transcription factors. The changes in Puf targets and in the Puf proteins indicate an integral involvement of RNA binding proteins and their RNA targets in the adaptation, reprogramming, and function of gene expression.

Automated Analysis and Classification of Histological Tissue Features by Multi-Dimensional Microscopic Molecular Profiling.

Characterization of the molecular attributes and spatial arrangements of cells and features within complex human tissues provides a critical basis for understanding processes involved in development and disease. Moreover, the ability to automate steps in the analysis and interpretation of histological images that currently require manual inspection by pathologists could revolutionize medical diagnostics. Toward this end, we developed a new imaging approach called multidimensional microscopic molecular profiling (MMMP) that can measure several independent molecular properties in situ at subcellular resolution for the same tissue specimen. MMMP involves repeated cycles of antibody or histochemical staining, imaging, and signal removal, which ultimately can generate information analogous to a multidimensional flow cytometry analysis on intact tissue sections. We performed a MMMP analysis on a tissue microarray containing a diverse set of 102 human tissues using a panel of 15 informative antibody and 5 histochemical stains plus DAPI. Large-scale unsupervised analysis of MMMP data, and visualization of the resulting classifications, identified molecular profiles that were associated with functional tissue features. We then directly annotated H&E images from this MMMP series such that canonical histological features of interest (e.g. blood vessels, epithelium, red blood cells) were individually labeled. By integrating image annotation data, we identified molecular signatures that were associated with specific histological annotations and we developed statistical models for automatically classifying these features. The classification accuracy for automated histology labeling was objectively evaluated using a cross-validation strategy, and significant accuracy (with a median per-pixel rate of 77% per feature from 15 annotated samples) for de novo feature prediction was obtained. These results suggest that high-dimensional profiling may advance the development of computer-based systems for automatically parsing relevant histological and cellular features from molecular imaging data of arbitrary human tissue samples, and can provide a framework and resource to spur the optimization of these technologies.

miR-142 regulates the tumorigenicity of human breast cancer stem cells through the canonical WNT signaling pathway.

MicroRNAs (miRNAs) are important regulators of stem and progenitor cell functions. We previously reported that miR-142 and miR-150 are upregulated in human breast cancer stem cells (BCSCs) as compared to the non-tumorigenic breast cancer cells. In this study, we report that miR-142 efficiently recruits the APC mRNA to an RNA-induced silencing complex, activates the canonical WNT signaling pathway in an APC-suppression dependent manner, and activates the expression of miR-150. Enforced expression of miR-142 or miR-150 in normal mouse mammary stem cells resulted in the regeneration of hyperproliferative mammary glands in vivo. Knockdown of endogenous miR-142 effectively suppressed organoid formation by BCSCs and slowed tumor growth initiated by human BCSCs in vivo. These results suggest that in some tumors, miR-142 regulates the properties of BCSCs at least in part by activating the WNT signaling pathway and miR-150 expression.

Transcriptome-wide mapping of pseudouridines: pseudouridine synthases modify specific mRNAs in S. cerevisiae.

We developed a novel technique, called pseudouridine site identification sequencing (PSI-seq), for the transcriptome-wide mapping of pseudouridylation sites with single-base resolution from cellular RNAs based on the induced termination of reverse transcription specifically at pseudouridines following CMCT treatment. PSI-seq analysis of RNA samples from S. cerevisiae correctly detected all of the 43 known pseudouridines in yeast 18S and 25S ribosomal RNA with high specificity. Moreover, application of PSI-seq to the yeast transcriptome revealed the presence of site-specific pseudouridylation within dozens of mRNAs, including RPL11a, TEF1, and other genes implicated in translation. To identify the mechanisms responsible for mRNA pseudouridylation, we genetically deleted candidate pseudouridine synthase (Pus) enzymes and reconstituted their activities in vitro. These experiments demonstrated that the Pus1 enzyme was necessary and sufficient for pseudouridylation of RPL11a mRNA, whereas Pus4 modified TEF1 mRNA, and Pus6 pseudouridylated KAR2 mRNA. Finally, we determined that modification of RPL11a at Ψ -68 was observed in RNA from the related yeast S. mikitae, and Ψ -239 in TEF1 mRNA was maintained in S. mikitae as well as S. pombe, indicating that these pseudouridylations are ancient, evolutionarily conserved RNA modifications. This work establishes that site-specific pseudouridylation of eukaryotic mRNAs is a genetically programmed RNA modification that naturally occurs in multiple yeast transcripts via distinct mechanisms, suggesting that mRNA pseudouridylation may provide an important novel regulatory function. The approach and strategies that we report here should be generally applicable to the discovery of pseudouridylation, or other RNA modifications, in diverse biological contexts.

Distinct stages of the translation elongation cycle revealed by sequencing ribosome-protected mRNA fragments.

During translation elongation, the ribosome ratchets along its mRNA template, incorporating each new amino acid and translocating from one codon to the next. The elongation cycle requires dramatic structural rearrangements of the ribosome. We show here that deep sequencing of ribosome-protected mRNA fragments reveals not only the position of each ribosome but also, unexpectedly, its particular stage of the elongation cycle. Sequencing reveals two distinct populations of ribosome footprints, 28-30 nucleotides and 20-22 nucleotides long, representing translating ribosomes in distinct states, differentially stabilized by specific elongation inhibitors. We find that the balance of small and large footprints varies by codon and is correlated with translation speed. The ability to visualize conformational changes in the ribosome during elongation, at single-codon resolution, provides a new way to study the detailed kinetics of translation and a new probe with which to identify the factors that affect each step in the elongation cycle.DOI: http://dx.doi.org/10.7554/eLife.01257.001.

Circular RNA is expressed across the eukaryotic tree of life.

An unexpectedly large fraction of genes in metazoans (human, mouse, zebrafish, worm, fruit fly) express high levels of circularized RNAs containing canonical exons. Here we report that circular RNA isoforms are found in diverse species whose most recent common ancestor existed more than one billion years ago: fungi (Schizosaccharomyces pombe and Saccharomyces cerevisiae), a plant (Arabidopsis thaliana), and protists (Plasmodium falciparum and Dictyostelium discoideum). For all species studied to date, including those in this report, only a small fraction of the theoretically possible circular RNA isoforms from a given gene are actually observed. Unlike metazoans, Arabidopsis, D. discoideum, P. falciparum, S. cerevisiae, and S. pombe have very short introns (∼ 100 nucleotides or shorter), yet they still produce circular RNAs. A minority of genes in S. pombe and P. falciparum have documented examples of canonical alternative splicing, making it unlikely that all circular RNAs are by-products of alternative splicing or 'piggyback' on signals used in alternative RNA processing. In S. pombe, the relative abundance of circular to linear transcript isoforms changed in a gene-specific pattern during nitrogen starvation. Circular RNA may be an ancient, conserved feature of eukaryotic gene expression programs.

Cell-type specific features of circular RNA expression.

Thousands of loci in the human and mouse genomes give rise to circular RNA transcripts; at many of these loci, the predominant RNA isoform is a circle. Using an improved computational approach for circular RNA identification, we found widespread circular RNA expression in Drosophila melanogaster and estimate that in humans, circular RNA may account for 1% as many molecules as poly(A) RNA. Analysis of data from the ENCODE consortium revealed that the repertoire of genes expressing circular RNA, the ratio of circular to linear transcripts for each gene, and even the pattern of splice isoforms of circular RNAs from each gene were cell-type specific. These results suggest that biogenesis of circular RNA is an integral, conserved, and regulated feature of the gene expression program.

Quantitative proteomic analysis reveals concurrent RNA-protein interactions and identifies new RNA-binding proteins in Saccharomyces cerevisiae.

A growing body of evidence supports the existence of an extensive network of RNA-binding proteins (RBPs) whose combinatorial binding affects the post-transcriptional fate of every mRNA in the cell-yet we still do not have a complete understanding of which proteins bind to mRNA, which of these bind concurrently, and when and where in the cell they bind. We describe here a method to identify the proteins that bind to RNA concurrently with an RBP of interest, using quantitative mass spectrometry combined with RNase treatment of affinity-purified RNA-protein complexes. We applied this method to the known RBPs Pab1, Nab2, and Puf3. Our method significantly enriched for known RBPs and is a clear improvement upon previous approaches in yeast. Our data reveal that some reported protein-protein interactions may instead reflect simultaneous binding to shared RNA targets. We also discovered more than 100 candidate RBPs, and we independently confirmed that 77% (23/30) bind directly to RNA. The previously recognized functions of the confirmed novel RBPs were remarkably diverse, and we mapped the RNA-binding region of one of these proteins, the transcriptional coactivator Mbf1, to a region distinct from its DNA-binding domain. Our results also provided new insights into the roles of Nab2 and Puf3 in post-transcriptional regulation by identifying other RBPs that bind simultaneously to the same mRNAs. While existing methods can identify sets of RBPs that interact with common RNA targets, our approach can determine which of those interactions are concurrent-a crucial distinction for understanding post-transcriptional regulation.

Improved discovery of molecular interactions in genome-scale data with adaptive model-based normalization.

BACKGROUND: High throughput molecular-interaction studies using immunoprecipitations (IP) or affinity purifications are powerful and widely used in biology research. One of many important applications of this method is to identify the set of RNAs that interact with a particular RNA-binding protein (RBP). Here, the unique statistical challenge presented is to delineate a specific set of RNAs that are enriched in one sample relative to another, typically a specific IP compared to a non-specific control to model background. The choice of normalization procedure critically impacts the number of RNAs that will be identified as interacting with an RBP at a given significance threshold - yet existing normalization methods make assumptions that are often fundamentally inaccurate when applied to IP enrichment data. METHODS: In this paper, we present a new normalization methodology that is specifically designed for identifying enriched RNA or DNA sequences in an IP. The normalization (called adaptive or AD normalization) uses a basic model of the IP experiment and is not a variant of mean, quantile, or other methodology previously proposed. The approach is evaluated statistically and tested with simulated and empirical data. RESULTS AND CONCLUSIONS: The adaptive (AD) normalization method results in a greatly increased range in the number of enriched RNAs identified, fewer false positives, and overall better concordance with independent biological evidence, for the RBPs we analyzed, compared to median normalization. The approach is also applicable to the study of pairwise RNA, DNA and protein interactions such as the analysis of transcription factors via chromatin immunoprecipitation (ChIP) or any other experiments where samples from two conditions, one of which contains an enriched subset of the other, are studied.

The yeast Rab GTPase Ypt1 modulates unfolded protein response dynamics by regulating the stability of HAC1 RNA.

The unfolded protein response (UPR) is a conserved mechanism that mitigates accumulation of unfolded proteins in the ER. The yeast UPR is subject to intricate post-transcriptional regulation, involving recruitment of the RNA encoding the Hac1 transcription factor to the ER and its unconventional splicing. To investigate the mechanisms underlying regulation of the UPR, we screened the yeast proteome for proteins that specifically interact with HAC1 RNA. Protein microarray experiments revealed that HAC1 interacts specifically with small ras GTPases of the Ypt family. We characterized the interaction of HAC1 RNA with one of these proteins, the yeast Rab1 homolog Ypt1. We found that Ypt1 protein specifically associated in vivo with unspliced HAC1 RNA. This association was disrupted by conditions that impaired protein folding in the ER and induced the UPR. Also, the Ypt1-HAC1 interaction depended on IRE1 and ADA5, the two genes critical for UPR activation. Decreasing expression of the Ypt1 protein resulted in a reduced rate of HAC1 RNA decay, leading to significantly increased levels of both unspliced and spliced HAC1 RNA, and delayed attenuation of the UPR, when ER stress was relieved. Our findings establish that Ypt1 contributes to regulation of UPR signaling dynamics by promoting the decay of HAC1 RNA, suggesting a potential regulatory mechanism for linking vesicle trafficking to the UPR and ER homeostasis.

Circular RNAs are the predominant transcript isoform from hundreds of human genes in diverse cell types.

Most human pre-mRNAs are spliced into linear molecules that retain the exon order defined by the genomic sequence. By deep sequencing of RNA from a variety of normal and malignant human cells, we found RNA transcripts from many human genes in which the exons were arranged in a non-canonical order. Statistical estimates and biochemical assays provided strong evidence that a substantial fraction of the spliced transcripts from hundreds of genes are circular RNAs. Our results suggest that a non-canonical mode of RNA splicing, resulting in a circular RNA isoform, is a general feature of the gene expression program in human cells.

Widespread mRNA association with cytoskeletal motor proteins and identification and dynamics of myosin-associated mRNAs in S. cerevisiae.

Programmed mRNA localization to specific subcellular compartments for localized translation is a fundamental mechanism of post-transcriptional regulation that affects many, and possibly all, mRNAs in eukaryotes. We describe here a systematic approach to identify the RNA cargoes associated with the cytoskeletal motor proteins of Saccharomyces cerevisiae in combination with live-cell 3D super-localization microscopy of endogenously tagged mRNAs. Our analysis identified widespread association of mRNAs with cytoskeletal motor proteins, including association of Myo3 with mRNAs encoding key regulators of actin branching and endocytosis such as WASP and WIP. Using conventional fluorescence microscopy and expression of MS2-tagged mRNAs from endogenous loci, we observed a strong bias for actin patch nucleator mRNAs to localize to the cell cortex and the actin patch in a Myo3- and F-actin dependent manner. Use of a double-helix point spread function (DH-PSF) microscope allowed super-localization measurements of single mRNPs at a spatial precision of 25 nm in x and y and 50 nm in z in live cells with 50 ms exposure times, allowing quantitative profiling of mRNP dynamics. The actin patch mRNA exhibited distinct and characteristic diffusion coefficients when compared to a control mRNA. In addition, disruption of F-actin significantly expanded the 3D confinement radius of an actin patch nucleator mRNA, providing a quantitative assessment of the contribution of the actin cytoskeleton to mRNP dynamic localization. Our results provide evidence for specific association of mRNAs with cytoskeletal motor proteins in yeast, suggest that different mRNPs have distinct and characteristic dynamics, and lend insight into the mechanism of actin patch nucleator mRNA localization to actin patches.

ESRRA-C11orf20 is a recurrent gene fusion in serous ovarian carcinoma.

Every year, ovarian cancer kills approximately 14,000 women in the United States and more than 140,000 women worldwide. Most of these deaths are caused by tumors of the serous histological type, which is rarely diagnosed before it has disseminated. By deep paired-end sequencing of mRNA from serous ovarian cancers, followed by deep sequencing of the corresponding genomic region, we identified a recurrent fusion transcript. The fusion transcript joins the 5' exons of ESRRA, encoding a ligand-independent member of the nuclear-hormone receptor superfamily, to the 3' exons of C11orf20, a conserved but uncharacterized gene located immediately upstream of ESRRA in the reference genome. To estimate the prevalence of the fusion, we tested 67 cases of serous ovarian cancer by RT-PCR and sequencing and confirmed its presence in 10 of these. Targeted resequencing of the corresponding genomic region from two fusion-positive tumor samples identified a nearly clonal chromosomal rearrangement positioning ESRRA upstream of C11orf20 in one tumor, and evidence of local copy number variation in the ESRRA locus in the second tumor. We hypothesize that the recurrent novel fusion transcript may play a role in pathogenesis of a substantial fraction of serous ovarian cancers and could provide a molecular marker for detection of the cancer. Gene fusions involving adjacent or nearby genes can readily escape detection but may play important roles in the development and progression of cancer.

Defining the specificity of cotranslationally acting chaperones by systematic analysis of mRNAs associated with ribosome-nascent chain complexes.

Polypeptides exiting the ribosome must fold and assemble in the crowded environment of the cell. Chaperones and other protein homeostasis factors interact with newly translated polypeptides to facilitate their folding and correct localization. Despite the extensive efforts, little is known about the specificity of the chaperones and other factors that bind nascent polypeptides. To address this question we present an approach that systematically identifies cotranslational chaperone substrates through the mRNAs associated with ribosome-nascent chain-chaperone complexes. We here focused on two Saccharomyces cerevisiae chaperones: the Signal Recognition Particle (SRP), which acts cotranslationally to target proteins to the ER, and the Nascent chain Associated Complex (NAC), whose function has been elusive. Our results provide new insights into SRP selectivity and reveal that NAC is a general cotranslational chaperone. We found surprising differential substrate specificity for the three subunits of NAC, which appear to recognize distinct features within nascent chains. Our results also revealed a partial overlap between the sets of nascent polypeptides that interact with NAC and SRP, respectively, and showed that NAC modulates SRP specificity and fidelity in vivo. These findings give us new insight into the dynamic interplay of chaperones acting on nascent chains. The strategy we used should be generally applicable to mapping the specificity, interplay, and dynamics of the cotranslational protein homeostasis network.

Stunned silence: gene expression programs in human cells infected with monkeypox or vaccinia virus.

Poxviruses use an arsenal of molecular weapons to evade detection and disarm host immune responses. We used DNA microarrays to investigate the gene expression responses to infection by monkeypox virus (MPV), an emerging human pathogen, and Vaccinia virus (VAC), a widely used model and vaccine organism, in primary human macrophages, primary human fibroblasts and HeLa cells. Even as the overwhelmingly infected cells approached their demise, with extensive cytopathic changes, their gene expression programs appeared almost oblivious to poxvirus infection. Although killed (gamma-irradiated) MPV potently induced a transcriptional program characteristic of the interferon response, no such response was observed during infection with either live MPV or VAC. Moreover, while the gene expression response of infected cells to stimulation with ionomycin plus phorbol 12-myristate 13-acetate (PMA), or poly (I-C) was largely unimpaired by infection with MPV, a cluster of pro-inflammatory genes were a notable exception. Poly(I-C) induction of genes involved in alerting the innate immune system to the infectious threat, including TNF-alpha, IL-1 alpha and beta, CCL5 and IL-6, were suppressed by infection with live MPV. Thus, MPV selectively inhibits expression of genes with critical roles in cell-signaling pathways that activate innate immune responses, as part of its strategy for stealthy infection.

Identification of RNA recognition elements in the Saccharomyces cerevisiae transcriptome.

Post-transcriptional regulation of gene expression, including mRNA localization, translation and decay, is ubiquitous yet still largely unexplored. How is the post-transcriptional regulatory program of each mRNA encoded in its sequence? Hundreds of specific RNA-binding proteins (RBPs) appear to play roles in mediating the post-transcriptional regulatory program, akin to the roles of specific DNA-binding proteins in transcription. As a step toward decoding the regulatory programs encoded in each mRNA, we focused on specific mRNA-protein interactions. We computationally analyzed the sequences of Saccharomyces cerevisiae mRNAs bound in vivo by 29 specific RBPs, identifying eight novel candidate motifs and confirming or extending six earlier reported recognition elements. Biochemical selections for RNA sequences selectively recognized by 12 yeast RBPs yielded novel motifs bound by Pin4, Nsr1, Hrb1, Gbp2, Sgn1 and Mrn1, and recovered the known recognition elements for Puf3, She2, Vts1 and Whi3. Most of the RNA elements we uncovered were associated with coherent mRNA expression changes and were significantly conserved in related yeasts, supporting their functional importance and suggesting that the corresponding RNA-protein interactions are evolutionarily conserved.

Three-dimensional tracking of single mRNA particles in Saccharomyces cerevisiae using a double-helix point spread function.

Optical imaging of single biomolecules and complexes in living cells provides a useful window into cellular processes. However, the three-dimensional dynamics of most important biomolecules in living cells remains essentially uncharacterized. The precise subcellular localization of mRNA-protein complexes plays a critical role in the spatial and temporal control of gene expression, and a full understanding of the control of gene expression requires precise characterization of mRNA transport dynamics beyond the optical diffraction limit. In this paper, we describe three-dimensional tracking of single mRNA particles with 25-nm precision in the x and y dimensions and 50-nm precision in the z dimension in live budding yeast cells using a microscope with a double-helix point spread function. Two statistical methods to detect intermittently confined and directed transport were used to quantify the three-dimensional trajectories of mRNA for the first time, using ARG3 mRNA as a model. Measurements and analysis show that the dynamics of ARG3 mRNA molecules are mostly diffusive, although periods of non-Brownian confinement and directed transport are observed. The quantitative methods detailed in this paper can be broadly applied to the study of mRNA localization and the dynamics of diverse other biomolecules in a wide variety of cell types.