Small molecule regulators of microRNAs identified by high-throughput screen coupled with high-throughput sequencing.

MicroRNAs (miRNAs) regulate fundamental biological processes by silencing mRNA targets and are dysregulated in many diseases. Therefore, miRNA replacement or inhibition can be harnessed as potential therapeutics. However, existing strategies for miRNA modulation using oligonucleotides and gene therapies are challenging, especially for neurological diseases, and none have yet gained clinical approval. We explore a different approach by screening a biodiverse library of small molecule compounds for their ability to modulate hundreds of miRNAs in human induced pluripotent stem cell-derived neurons. We demonstrate the utility of the screen by identifying cardiac glycosides as potent inducers of miR-132, a key neuroprotective miRNA downregulated in Alzheimer's disease and other tauopathies. Coordinately, cardiac glycosides downregulate known miR-132 targets, including Tau, and protect rodent and human neurons against various toxic insults. More generally, our dataset of 1370 drug-like compounds and their effects on the miRNome provides a valuable resource for further miRNA-based drug discovery.

Small Molecule Regulators of microRNAs Identified by High-Throughput Screen Coupled with High-Throughput Sequencing.

MicroRNAs (miRNAs) regulate fundamental biological processes by silencing mRNA targets and are dysregulated in many diseases. Therefore, miRNA replacement or inhibition can be harnessed as potential therapeutics. However, existing strategies for miRNA modulation using oligonucleotides and gene therapies are challenging, especially for neurological diseases, and none have yet gained clinical approval. We explore a different approach by screening a biodiverse library of small molecule compounds for their ability to modulate hundreds of miRNAs in human induced pluripotent stem cell-derived neurons. We demonstrate the utility of the screen by identifying cardiac glycosides as potent inducers of miR-132, a key miRNA downregulated in Alzheimer's disease and other tauopathies. Coordinately, cardiac glycosides downregulate known miR-132 targets, including Tau, and protect rodent and human neurons against various toxic insults. More generally, our dataset of 1370 drug-like compounds and their effects on the miRNome provide a valuable resource for further miRNA-based drug discovery.

Decreasing miRNA sequencing bias using a single adapter and circularization approach.

The ability to accurately quantify all the microRNAs (miRNAs) in a sample is important for understanding miRNA biology and for development of new biomarkers and therapeutic targets. We develop a new method for preparing miRNA sequencing libraries, RealSeq®-AC, that involves ligating the miRNAs with a single adapter and circularizing the ligation products. When compared to other methods, RealSeq®-AC provides greatly reduced miRNA sequencing bias and allows the identification of the largest variety of miRNAs in biological samples. This reduced bias also allows robust quantification of miRNAs present in samples across a wide range of RNA input levels.

Coordinated tissue-specific regulation of adjacent alternative 3' splice sites in C. elegans.

Adjacent alternative 3' splice sites, those separated by ≤18 nucleotides, provide a unique problem in the study of alternative splicing regulation; there is overlap of the cis-elements that define the adjacent sites. Identification of the intron's 3' end depends upon sequence elements that define the branchpoint, polypyrimidine tract, and terminal AG dinucleotide. Starting with RNA-seq data from germline-enriched and somatic cell-enriched Caenorhabditis elegans samples, we identify hundreds of introns with adjacent alternative 3' splice sites. We identify 203 events that undergo tissue-specific alternative splicing. For these, the regulation is monodirectional, with somatic cells preferring to splice at the distal 3' splice site (furthest from the 5' end of the intron) and germline cells showing a distinct shift toward usage of the adjacent proximal 3' splice site (closer to the 5' end of the intron). Splicing patterns in somatic cells follow C. elegans consensus rules of 3' splice site definition; a short stretch of pyrimidines preceding an AG dinucleotide. Splicing in germline cells occurs at proximal 3' splice sites that lack a preceding polypyrimidine tract, and in three instances the germline-specific site lacks the AG dinucleotide. We provide evidence that use of germline-specific proximal 3' splice sites is conserved across Caenorhabditis species. We propose that there are differences between germline and somatic cells in the way that the basal splicing machinery functions to determine the intron terminus.

Co-option of the piRNA pathway for germline-specific alternative splicing of C. elegans TOR.

Many eukaryotic genes contain embedded antisense transcripts and repetitive sequences of unknown function. We report that male germline-specific expression of an antisense transcript contained in an intron of C. elegans Target of Rapamycin (TOR, let-363) is associated with (1) accumulation of endo-small interfering RNAs (siRNAs) against an embedded Helitron transposon and (2) activation of an alternative 3' splice site of TOR. The germline-specific Argonaute proteins PRG-1 and CSR-1, which participate in self/nonself RNA recognition, antagonistically regulate the generation of these endo-siRNAs, TOR mRNA levels, and 3' splice-site selection. Supply of exogenous double-stranded RNA against the region of sense/antisense overlap reverses changes in TOR expression and splicing and suppresses the progressive multigenerational sterility phenotype of prg-1 mutants. We propose that recognition of a "nonself" intronic transposon by endo-siRNAs/the piRNA system provides physiological regulation of expression and alternative splicing of a host gene that, in turn, contributes to the maintenance of germline function across generations.

Alternative Splicing Regulation of Cancer-Related Pathways in Caenorhabditis elegans: An In Vivo Model System with a Powerful Reverse Genetics Toolbox.

Alternative splicing allows for the generation of protein diversity and fine-tunes gene expression. Several model systems have been used for the in vivo study of alternative splicing. Here we review the use of the nematode Caenorhabditis elegans to study splicing regulation in vivo. Recent studies have shown that close to 25% of genes in the worm genome undergo alternative splicing. A big proportion of these events are functional, conserved, and under strict regulation either across development or other conditions. Several techniques like genome-wide RNAi screens and bichromatic reporters are available for the study of alternative splicing in worms. In this review, we focus, first, on the main studies that have been performed to dissect alternative splicing in this system and later on examples from genes that have human homologs that are implicated in cancer. The significant advancement towards understanding the regulation of alternative splicing and cancer that the C. elegans system has offered is discussed.

Epigenomic analysis detects widespread gene-body DNA hypomethylation in chronic lymphocytic leukemia.

We have extensively characterized the DNA methylomes of 139 patients with chronic lymphocytic leukemia (CLL) with mutated or unmutated IGHV and of several mature B-cell subpopulations through the use of whole-genome bisulfite sequencing and high-density microarrays. The two molecular subtypes of CLL have differing DNA methylomes that seem to represent epigenetic imprints from distinct normal B-cell subpopulations. DNA hypomethylation in the gene body, targeting mostly enhancer sites, was the most frequent difference between naive and memory B cells and between the two molecular subtypes of CLL and normal B cells. Although DNA methylation and gene expression were poorly correlated, we identified gene-body CpG dinucleotides whose methylation was positively or negatively associated with expression. We have also recognized a DNA methylation signature that distinguishes new clinico-biological subtypes of CLL. We propose an epigenomic scenario in which differential methylation in the gene body may have functional and clinical implications in leukemogenesis.

Co-regulation of alternative splicing by diverse splicing factors in Caenorhabditis elegans.

Regulation of alternative splicing is controlled by pre-mRNA sequences (cis-elements) and trans-acting protein factors that bind them. The combinatorial interactions of multiple protein factors with the cis-elements surrounding a given alternative splicing event lead to an integrated splicing decision. The mechanism of multifactorial splicing regulation is poorly understood. Using a splicing-sensitive DNA microarray, we assayed 352 Caenorhabditis elegans alternative cassette exons for changes in embryonic splicing patterns between wild-type and 12 different strains carrying mutations in a splicing factor. We identified many alternative splicing events that are regulated by multiple splicing factors. Many splicing factors have the ability to behave as splicing repressors for some alternative cassette exons and as splicing activators for others. Unexpectedly, we found that the ability of a given alternative splicing factor to behave as an enhancer or repressor of a specific splicing event can change during development. Our observations that splicing factors can change their effects on a substrate during development support a model in which combinatorial effects of multiple factors, both constitutive and developmentally regulated ones, contribute to the overall splicing decision.

HRP-2, the Caenorhabditis elegans homolog of mammalian heterogeneous nuclear ribonucleoproteins Q and R, is an alternative splicing factor that binds to UCUAUC splicing regulatory elements.

Alternative splicing is regulated by cis sequences in the pre-mRNA that serve as binding sites for trans-acting alternative splicing factors. In a previous study, we used bioinformatics and molecular biology to identify and confirm that the intronic hexamer sequence UCUAUC is a nematode alternative splicing regulatory element. In this study, we used RNA affinity chromatography to identify trans factors that bind to this sequence. HRP-2, the Caenorhabditis elegans homolog of human heterogeneous nuclear ribonucleoproteins Q and R, binds to UCUAUC in the context of unc-52 intronic regulatory sequences as well as to RNAs containing tandem repeats of this sequence. The three Us in the hexamer are the most important determinants of this binding specificity. We demonstrate, using RNA interference, that HRP-2 regulates the alternative splicing of two genes, unc-52 and lin-10, both of which have cassette exons flanked by an intronic UCUAUC motif. We propose that HRP-2 is a protein responsible for regulating alternative splicing through binding interactions with the UCUAUC sequence.

Global analysis of alternative splicing uncovers developmental regulation of nonsense-mediated decay in C. elegans.

Alternative splicing coupled to nonsense-mediated decay (AS-NMD) is a mechanism for post-transcriptional regulation of gene expression. We analyzed the global effects of mutations in seven genes of the C. elegans NMD pathway on AS isoform ratios. We find that mutations in two NMD factors, smg-6 and smg-7, have weaker global effects relative to mutations in other smg genes. We did an in-depth analysis of 12 pre-mRNA splicing factor genes that are subject to AS-NMD. For four of these, changes in the ratio of alternatively spliced isoforms during development are caused by developmentally regulated inhibition of NMD, and not by changes in alternative splicing. Using sucrose gradient analysis of mRNAs undergoing translation, we find several examples of NMD-dependent enrichment of premature termination codon (PTC) isoforms in the monosome fraction. In contrast, we present evidence of two genes for which the PTC-containing isoforms are found in polysomes and have a translational profile similar to non-PTC-containing transcripts from the same gene. We propose that NMD of certain alternatively spliced isoforms is regulated, and that some stabilized NMD targets may be translated.

Alternative splicing and the steady-state ratios of mRNA isoforms generated by it are under strong stabilizing selection in Caenorhabditis elegans.

Evolutionary studies indicate that a high proportion of alternative splicing (AS) events are species-specific; just 28% of minor-form alternatively spliced exons are conserved between mice and humans. We employed a splicing-sensitive microarray to study the evolution of allele-specific AS in nematodes. We compared splicing levels among five distinct Caenorhabditis elegans lines. Our results indicate that AS is less variable between natural isolates (NIs) from England, Hawaii, and Australia than when compared with mutation accumulation lines (6% vs. 21%, respectively, vary compared with N2). This suggests that strong stabilizing selection shapes the evolution of the ratios of isoforms generated by AS in C. elegans. When we analyzed some of the splicing changes between the NIs, we found examples of changes in both cis and trans that lead to alterations in gene-specific AS. This indicates that both these mechanisms for changing AS are employed along the path toward speciation in nematodes.

Alternative splicing regulation during C. elegans development: splicing factors as regulated targets.

Alternative splicing generates protein diversity and allows for post-transcriptional gene regulation. Estimates suggest that 10% of the genes in Caenorhabditis elegans undergo alternative splicing. We constructed a splicing-sensitive microarray to detect alternative splicing for 352 cassette exons and tested for changes in alternative splicing of these genes during development. We found that the microarray data predicted that 62/352 (approximately 18%) of the alternative splicing events studied show a strong change in the relative levels of the spliced isoforms (>4-fold) during development. Confirmation of the microarray data by RT-PCR was obtained for 70% of randomly selected genes tested. Among the genes with the most developmentally regulated alternatively splicing was the hnRNP F/H splicing factor homolog, W02D3.11 - now named hrpf-1. For the cassette exon of hrpf-1, the inclusion isoform comprises 65% of hrpf-1 steady state messages in embryos but only 0.1% in the first larval stage. This dramatic change in the alternative splicing of an alternative splicing factor suggests a complex cascade of splicing regulation during development. We analyzed splicing in embryos from a strain with a mutation in the splicing factor sym-2, another hnRNP F/H homolog. We found that approximately half of the genes with large alternative splicing changes between the embryo and L1 stages are regulated by sym-2 in embryos. An analysis of the role of nonsense-mediated decay in regulating steady-state alternative mRNA isoforms was performed. We found that 8% of the 352 events studied have alternative isoforms whose relative steady-state levels in embryos change more than 4-fold in a nonsense-mediated decay mutant, including hrpf-1. Strikingly, 53% of these alternative splicing events that are affected by NMD in our experiment are not obvious substrates for NMD based on the presence of premature termination codons. This suggests that the targeting of splicing factors by NMD may have downstream effects on alternative splicing regulation.

Distinct ribonucleoprotein reservoirs for microRNA and siRNA populations in C. elegans.

MicroRNAs (miRNAs) are regulatory molecules that share both biosynthetic derivation (cleavage from short hairpin precursor RNAs) and functional roles (downregulation of specific mRNAs through targeted degradation and/or translational inhibition). A distinct family of small RNAs, termed siRNAs, have some common characteristics but exhibit distinct modes of biosynthesis and function. In this study, we report procedures for purification of a predominant species of miRNA-containing ribonucleoprotein complexes from Caenorhabditis elegans and demonstrate that this population is distinct from the predominant pool of siRNA-containing ribonucleoprotein complexes. An observed miRNP-associated RNA population consisting predominantly (>95%) of miRNAs supported the unique identity of miRNPs as biological effectors within the cell, provided clean material for analysis of changes in miRNA spectra during development, and provided strong evidence of miRNA character for a number of novel small RNAs. Likewise, the RNA spectrum derived from partial siRNP purification was useful in defining functional characteristics of this more diverse population of small RNAs.

Intronic alternative splicing regulators identified by comparative genomics in nematodes.

Many alternative splicing events are regulated by pentameric and hexameric intronic sequences that serve as binding sites for splicing regulatory factors. We hypothesized that intronic elements that regulate alternative splicing are under selective pressure for evolutionary conservation. Using a Wobble Aware Bulk Aligner genomic alignment of Caenorhabditis elegans and Caenorhabditis briggsae, we identified 147 alternatively spliced cassette exons that exhibit short regions of high nucleotide conservation in the introns flanking the alternative exon. In vivo experiments on the alternatively spliced let-2 gene confirm that these conserved regions can be important for alternative splicing regulation. Conserved intronic element sequences were collected into a dataset and the occurrence of each pentamer and hexamer motif was counted. We compared the frequency of pentamers and hexamers in the conserved intronic elements to a dataset of all C. elegans intron sequences in order to identify short intronic motifs that are more likely to be associated with alternative splicing. High-scoring motifs were examined for upstream or downstream preferences in introns surrounding alternative exons. Many of the high-scoring nematode pentamer and hexamer motifs correspond to known mammalian splicing regulatory sequences, such as (T)GCATG, indicating that the mechanism of alternative splicing regulation is well conserved in metazoans. A comparison of the analysis of the conserved intronic elements, and analysis of the entire introns flanking these same exons, reveals that focusing on intronic conservation can increase the sensitivity of detecting putative splicing regulatory motifs. This approach also identified novel sequences whose role in splicing is under investigation and has allowed us to take a step forward in defining a catalog of splicing regulatory elements for an organism. In vivo experiments confirm that one novel high-scoring sequence from our analysis, (T)CTATC, is important for alternative splicing regulation of the unc-52 gene.

A computational screen for mammalian pseudouridylation guide H/ACA RNAs.

The box H/ACA RNA gene family is one of the largest non-protein-coding gene families in eukaryotes and archaea. Recently, we developed snoGPS, a computational screening program for H/ACA snoRNAs, and applied it to Saccharomyces cerevisiae. We report here results of extending our method to screen for H/ACA RNAs in multiple large genomes of related species, and apply it to the human, mouse, and rat genomes. Because of the 250-fold larger search space compared to S. cerevisiae, significant enhancements to our algorithms were required. Complementing extensive cloning experiments performed by others, our findings include the detection and experimental verification of seven new mammalian H/ACA RNAs and the prediction of 23 new H/ACA RNA pseudouridine guide assignments. These assignments include four for H/ACA RNAs previously classified as orphan H/ACA RNAs with no known targets. We also determined systematic syntenic conservation among human and mouse H/ACA RNAs. With this work, 82 of 97 ribosomal RNA pseudouridines and 18 of 32 spliceosomal RNA pseudouridines in mammals have been linked to H/ACA guide RNAs.