Conserved function of lincRNAs in vertebrate embryonic development despite rapid sequence evolution.

Thousands of long intervening noncoding RNAs (lincRNAs) have been identified in mammals. To better understand the evolution and functions of these enigmatic RNAs, we used chromatin marks, poly(A)-site mapping and RNA-Seq data to identify more than 550 distinct lincRNAs in zebrafish. Although these shared many characteristics with mammalian lincRNAs, only 29 had detectable sequence similarity with putative mammalian orthologs, typically restricted to a single short region of high conservation. Other lincRNAs had conserved genomic locations without detectable sequence conservation. Antisense reagents targeting conserved regions of two zebrafish lincRNAs caused developmental defects. Reagents targeting splice sites caused the same defects and were rescued by adding either the mature lincRNA or its human or mouse ortholog. Our study provides a roadmap for identification and analysis of lincRNAs in model organisms and shows that lincRNAs play crucial biological roles during embryonic development with functionality conserved despite limited sequence conservation.

SPOCD1 is an essential executor of piRNA-directed de novo DNA methylation.

In mammals, the acquisition of the germline from the soma provides the germline with an essential challenge: the need to erase and reset genomic methylation1. In the male germline, RNA-directed DNA methylation silences young, active transposable elements2-4. The PIWI protein MIWI2 (PIWIL4) and its associated PIWI-interacting RNAs (piRNAs) instruct DNA methylation of transposable elements3,5. piRNAs are proposed to tether MIWI2 to nascent transposable element transcripts; however, the mechanism by which MIWI2 directs the de novo methylation of transposable elements is poorly understood, although central to the immortality of the germline. Here we define the interactome of MIWI2 in mouse fetal gonocytes undergoing de novo genome methylation and identify a previously unknown MIWI2-associated factor, SPOCD1, that is essential for the methylation and silencing of young transposable elements. The loss of Spocd1 in mice results in male-specific infertility but does not affect either piRNA biogenesis or the localization of MIWI2 to the nucleus. SPOCD1 is a nuclear protein whose expression is restricted to the period of de novo genome methylation. It co-purifies in vivo with DNMT3L and DNMT3A, components of the de novo methylation machinery, as well as with constituents of the NURD and BAF chromatin remodelling complexes. We propose a model whereby tethering of MIWI2 to a nascent transposable element transcript recruits repressive chromatin remodelling activities and the de novo methylation apparatus through SPOCD1. In summary, we have identified a previously unrecognized and essential executor of mammalian piRNA-directed DNA methylation.

lncRNAs in development and differentiation: from sequence motifs to functional characterization.

The number of long noncoding RNAs (lncRNAs) with characterized developmental and cellular functions continues to increase, but our understanding of the molecular mechanisms underlying lncRNA functions, and how they are dictated by RNA sequences, remains limited. Relatively short, conserved sequence motifs embedded in lncRNA transcripts are often important determinants of lncRNA localization, stability and interactions. Identifying such RNA motifs remains challenging due to the substantial length of lncRNA transcripts and the rapid evolutionary turnover of lncRNA sequences. Nevertheless, the recent discovery of specific RNA elements, together with their experimental interrogation, has enabled the first step in classifying heterogeneous lncRNAs into sub-groups with similar molecular mechanisms and functions. In this Review, we focus on lncRNAs with roles in development, cell differentiation and normal physiology in vertebrates, and we discuss the sequence elements defining their functions. We also summarize progress on the discovery of regulatory RNA sequence elements, as well as their molecular functions and interaction partners.

Genetic tools to dissect functions of long noncoding RNAs.

LncRNAs represent an abundant group of noncoding transcripts, some of which carry out important regulatory functions. To survey the biological and molecular roles of lncRNAs, reliable strategies for their genetic inactivation are required. Several lncRNA features make them challenging to target by genome editing. First, lncRNA loci often span large genomic distances. As such, full or partial deletion alleles are not always easy to generate and interpret as they might affect DNA regulatory elements. Second, in contrast to proteins, lncRNA transcripts are usually resistant to the minimally invasive approach of point substitutions. Third, lncRNA sequences exhibit rapid evolutionary turnover, impeding prediction and targeting of the specific functional sequence elements. Nonetheless, advances in genome editing and comparative genomics have expanded the repertoire of genetic strategies to dissect lncRNA functions in model organisms and cell lines. In this review, we discuss several approaches that have been used to generate lncRNA mutant alleles, focusing on vertebrate lncRNAs. We also briefly highlight comparative genomics approaches to identify conserved lncRNA sequence motifs, which represent attractive target sequences to abrogate lncRNA functions and to pinpoint functional contributions of these elements.

Strategies for genetic inactivation of long noncoding RNAs in zebrafish.

The number of annotated long noncoding RNAs (lncRNAs) continues to grow; however, their functional characterization in model organisms has been hampered by the lack of reliable genetic inactivation strategies. While partial or full deletions of lncRNA loci disrupt lncRNA expression, they do not permit the formal association of a phenotype with the encoded transcript. Here, we examined several alternative strategies for generating lncRNA null alleles in zebrafish and found that they often resulted in unpredicted changes to lncRNA expression. Removal of the transcription start sites (TSSs) of lncRNA genes resulted in hypomorphic mutants, due to the usage of either constitutive or tissue-specific alternative TSSs. Deletions of short, highly conserved lncRNA regions can also lead to overexpression of truncated transcripts. In contrast, knock-in of a polyadenylation signal enabled complete inactivation of malat1, the most abundant vertebrate lncRNA. In summary, lncRNA null alleles require extensive in vivo validation, and we propose insertion of transcription termination sequences as the most reliable approach to generate lncRNA-deficient zebrafish.

Comparative analyses of super-enhancers reveal conserved elements in vertebrate genomes.

Super-enhancers (SEs) are key transcriptional drivers of cellular, developmental, and disease states in mammals, yet the conservational and regulatory features of these enhancer elements in nonmammalian vertebrates are unknown. To define SEs in zebrafish and enable sequence and functional comparisons to mouse and human SEs, we used genome-wide histone H3 lysine 27 acetylation (H3K27ac) occupancy as a primary SE delineator. Our study determined the set of SEs in pluripotent state cells and adult zebrafish tissues and revealed both similarities and differences between zebrafish and mammalian SEs. Although the total number of SEs was proportional to the genome size, the genomic distribution of zebrafish SEs differed from that of the mammalian SEs. Despite the evolutionary distance separating zebrafish and mammals and the low overall SE sequence conservation, ∼42% of zebrafish SEs were located in close proximity to orthologs that also were associated with SEs in mouse and human. Compared to their nonassociated counterparts, higher sequence conservation was revealed for those SEs that have maintained orthologous gene associations. Functional dissection of two of these SEs identified conserved sequence elements and tissue-specific expression patterns, while chromatin accessibility analyses predicted transcription factors governing the function of pluripotent state zebrafish SEs. Our zebrafish annotations and comparative studies show the extent of SE usage and their conservation across vertebrates, permitting future gene regulatory studies in several tissues.

Transposon-driven transcription is a conserved feature of vertebrate spermatogenesis and transcript evolution.

Spermatogenesis is associated with major and unique changes to chromosomes and chromatin. Here, we sought to understand the impact of these changes on spermatogenic transcriptomes. We show that long terminal repeats (LTRs) of specific mouse endogenous retroviruses (ERVs) drive the expression of many long non-coding transcripts (lncRNA). This process occurs post-mitotically predominantly in spermatocytes and round spermatids. We demonstrate that this transposon-driven lncRNA expression is a conserved feature of vertebrate spermatogenesis. We propose that transposon promoters are a mechanism by which the genome can explore novel transcriptional substrates, increasing evolutionary plasticity and allowing for the genesis of novel coding and non-coding genes. Accordingly, we show that a small fraction of these novel ERV-driven transcripts encode short open reading frames that produce detectable peptides. Finally, we find that distinct ERV elements from the same subfamilies act as differentially activated promoters in a tissue-specific context. In summary, we demonstrate that LTRs can act as tissue-specific promoters and contribute to post-mitotic spermatogenic transcriptome diversity.

Expanding the microRNA targeting code: functional sites with centered pairing.

Most metazoan microRNA (miRNA) target sites have perfect pairing to the seed region, located near the miRNA 5' end. Although pairing to the 3' region sometimes supplements seed matches or compensates for mismatches, pairing to the central region has been known to function only at rare sites that impart Argonaute-catalyzed mRNA cleavage. Here, we present "centered sites," a class of miRNA target sites that lack both perfect seed pairing and 3'-compensatory pairing and instead have 11-12 contiguous Watson-Crick pairs to the center of the miRNA. Although centered sites can impart mRNA cleavage in vitro (in elevated Mg(2+)), in cells they repress protein output without consequential Argonaute-catalyzed cleavage. Our study also identified extensively paired sites that are cleavage substrates in cultured cells and human brain. This expanded repertoire of cleavage targets and the identification of the centered site type help explain why central regions of many miRNAs are evolutionarily conserved.

Coherent but overlapping expression of microRNAs and their targets during vertebrate development.

MicroRNAs (miRNAs) are small noncoding RNAs that direct post-transcriptional repression of protein-coding genes. In vertebrates, each highly conserved miRNA typically regulates hundreds of target mRNAs. However, the precise relationship between expression of the miRNAs and that of their targets has remained unclear, in part because of the scarcity of quantitative expression data at cellular resolution. Here we report quantitative analyses of mRNA levels in miRNA-expressing cells of the zebrafish embryo, capturing entire miRNA expression domains, purified to cellular resolution using fluorescent-activated cell sorting (FACS). Focus was on regulation by miR-206 and miR-133 in the developing somites and miR-124 in the developing central nervous system. Comparison of wild-type embryos and those lacking miRNAs revealed predicted targets that responded to the miRNAs and distinguished miRNA-mediated mRNA destabilization from other regulatory effects. For all three miRNAs examined, expression of the miRNAs and that of their predicted targets usually overlapped. A few targets were expressed at higher levels in miRNA-expressing cells than in the rest of the embryo, demonstrating that miRNA-mediated repression can act in opposition to other regulatory processes. However, for most targets expression was lower in miRNA-expressing cells than in the rest of the embryo, indicating that miRNAs usually operate in concert with the other regulatory machinery of the cell.

Extensive alternative polyadenylation during zebrafish development.

The post-transcriptional fate of messenger RNAs (mRNAs) is largely dictated by their 3' untranslated regions (3' UTRs), which are defined by cleavage and polyadenylation (CPA) of pre-mRNAs. We used poly(A)-position profiling by sequencing (3P-seq) to map poly(A) sites at eight developmental stages and tissues in the zebrafish. Analysis of over 60 million 3P-seq reads substantially increased and improved existing 3' UTR annotations, resulting in confidently identified 3' UTRs for >79% of the annotated protein-coding genes in zebrafish. mRNAs from most zebrafish genes undergo alternative CPA, with those from more than a thousand genes using different dominant 3' UTRs at different stages. These included one of the poly(A) polymerase genes, for which alternative CPA reinforces its repression in the ovary. 3' UTRs tend to be shortest in the ovaries and longest in the brain. Isoforms with some of the shortest 3' UTRs are highly expressed in the ovary, yet absent in the maternally contributed RNAs of the embryo, perhaps because their 3' UTRs are too short to accommodate a uridine-rich motif required for stability of the maternal mRNA. At 2 h post-fertilization, thousands of unique poly(A) sites appear at locations lacking a typical polyadenylation signal, which suggests a wave of widespread cytoplasmic polyadenylation of mRNA degradation intermediates. Our insights into the identities, formation, and evolution of zebrafish 3' UTRs provide a resource for studying gene regulation during vertebrate development.

The lincRNA JUNI regulates the stress-dependent induction of c-Jun, cellular migration and survival through the modulation of the DUSP14-JNK axis.

Cancer cells employ adaptive mechanisms to survive various stressors, including genotoxic drugs. Understanding the factors promoting survival is crucial for developing effective treatments. In this study, we unveil a previously unexplored long non-coding RNA, JUNI (JUN-DT, LINC01135), which is upregulated by genotoxic drugs through the activation of stress-activated MAPKs, JNK, and p38 and consequently exerts positive control over the expression of its adjacent gene product c-Jun, a well-known oncoprotein, which transduces signals to multiple transcriptional outputs. JUNI regulates cellular migration and has a crucial role in conferring cellular resistance to chemotherapeutic drugs or UV radiation. Depletion of JUNI markedly increases the sensitivity of cultured cells and spheroids to chemotherapeutic agents. We identified 57 proteins interacting with JUNI. The activity of one of them the MAPK phosphatase and inhibitor, DUSP14, is counteracted by JUNI, thereby, facilitating efficient JNK phosphorylation and c-Jun induction when cells are exposed to UV radiation. The antagonistic interplay with DUSP14 contributes not only to c-Jun induction but also augments the survival of UV-exposed cells. In summary, we introduce JUNI as a novel stress-inducible regulator of c-Jun, positioning it as a potential target for enhancing the sensitivity of cancer cells to chemotherapy.

Novel non-stimulants rescue hyperactive phenotype in an adgrl3.1 mutant zebrafish model of ADHD.

ADHD is a highly prevalent neurodevelopmental disorder. The first-line therapeutic for ADHD, methylphenidate, can cause serious side effects including weight loss, insomnia, and hypertension. Therefore, the development of non-stimulant-based therapeutics has been prioritized. However, many of these also cause other effects, most notably somnolence. Here, we have used a uniquely powerful genetic model and unbiased drug screen to identify novel ADHD non-stimulant therapeutics. We first found that adgrl3.1 null (adgrl3.1-/-) zebrafish larvae showed a robust hyperactive phenotype. Although the hyperactivity was rescued by three ADHD non-stimulant therapeutics, all interfered significantly with sleep. Second, we used wild-type zebrafish larvae to characterize a simple behavioral phenotype generated by atomoxetine and screened the 1200 compound Prestwick Chemical Library® for a matching behavioral profile resulting in 67 hits. These hits were re-assayed in the adgrl3.1-/-. Using the previously identified non-stimulants as a positive control, we identified four compounds that matched the effect of atomoxetine: aceclofenac, amlodipine, doxazosin, and moxonidine. We additionally demonstrated cognitive effects of moxonidine in mice using a T-maze spontaneous alternation task. Moxonidine, has high affinity for imidazoline 1 receptors. We, therefore, assayed a pure imidazoline 1 agonist, LNP599, which generated an effect closely matching other non-stimulant ADHD therapeutics suggesting a role for this receptor system in ADHD. In summary, we introduce a genetic model of ADHD in zebrafish and identify five putative therapeutics. The findings offer a novel tool for understanding the neural circuits of ADHD, suggest a novel mechanism for its etiology, and identify novel therapeutics.

adgrl3.1-deficient zebrafish show noradrenaline-mediated externalizing behaviors, and altered expression of externalizing disorder-candidate genes, suggesting functional targets for treatment.

Externalizing disorders (ED) are a cause of concern for public health, and their high heritability makes genetic risk factors a priority for research. Adhesion G-Protein-Coupled Receptor L3 (ADGRL3) is strongly linked to several EDs, and loss-of-function models have shown the impacts of this gene on several core ED-related behaviors. For example, adgrl3.1-/- zebrafish show high levels of hyperactivity. However, our understanding of the mechanisms by which this gene influences behavior is incomplete. Here we characterized, for the first time, externalizing behavioral phenotypes of adgrl3.1-/- zebrafish and found them to be highly impulsive, show risk-taking in a novel environment, have attentional deficits, and show high levels of hyperactivity. All of these phenotypes were rescued by atomoxetine, demonstrating noradrenergic mediation of the externalizing effects of adgrl3.1. Transcriptomic analyses of the brains of adgrl3.1-/- vs. wild-type fish revealed several differentially expressed genes and enriched gene clusters that were independent of noradrenergic manipulation. This suggests new putative functional pathways underlying ED-related behaviors, and potential targets for the treatment of ED.

Demarcation of Topologically Associating Domains Is Uncoupled from Enriched CTCF Binding in Developing Zebrafish.

CCCTC-binding factor (CTCF) is a conserved architectural protein that plays crucial roles in gene regulation and three-dimensional (3D) chromatin organization. To better understand mechanisms and evolution of vertebrate genome organization, we analyzed genome occupancy of CTCF in zebrafish utilizing an endogenously epitope-tagged CTCF knock-in allele. Zebrafish CTCF shares similar facets with its mammalian counterparts, including binding to enhancers, active promoters and repeat elements, and bipartite sequence motifs of its binding sites. However, we found that in vivo CTCF binding is not enriched at boundaries of topologically associating domains (TADs) in developing zebrafish, whereas TAD demarcation by chromatin marks did not differ from mammals. Our data suggest that general mechanisms underlying 3D chromatin organization, and in particular the involvement of CTCF in this process, differ between distant vertebrate species.

Author Correction: In-cell identification and measurement of RNA-protein interactions.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

In-cell identification and measurement of RNA-protein interactions.

Regulatory RNAs exert their cellular functions through RNA-binding proteins (RBPs). Identifying RNA-protein interactions is therefore key for a molecular understanding of regulatory RNAs. To date, RNA-bound proteins have been identified primarily through RNA purification followed by mass spectrometry. Here, we develop incPRINT (in cell protein-RNA interaction), a high-throughput method to identify in-cell RNA-protein interactions revealed by quantifiable luminescence. Applying incPRINT to long noncoding RNAs (lncRNAs), we identify RBPs specifically interacting with the lncRNA Firre and three functionally distinct regions of the lncRNA Xist. incPRINT confirms previously known lncRNA-protein interactions and identifies additional interactions that had evaded detection with other approaches. Importantly, the majority of the incPRINT-defined interactions are specific to individual functional regions of the large Xist transcript. Thus, we present an RNA-centric method that enables reliable identification of RNA-region-specific RBPs and is applicable to any RNA of interest.

MicroRNA degradation by a conserved target RNA regulates animal behavior.

microRNAs (miRNAs) repress target transcripts through partial complementarity. By contrast, highly complementary miRNA-binding sites within viral and artificially engineered transcripts induce miRNA degradation in vitro and in cell lines. Here, we show that a genome-encoded transcript harboring a near-perfect and deeply conserved miRNA-binding site for miR-29 controls zebrafish and mouse behavior. This transcript originated in basal vertebrates as a long noncoding RNA (lncRNA) and evolved to the protein-coding gene NREP in mammals, where the miR-29-binding site is located within the 3' UTR. We show that the near-perfect miRNA site selectively triggers miR-29b destabilization through 3' trimming and restricts its spatial expression in the cerebellum. Genetic disruption of the miR-29 site within mouse Nrep results in ectopic expression of cerebellar miR-29b and impaired coordination and motor learning. Thus, we demonstrate an endogenous target-RNA-directed miRNA degradation event and its requirement for animal behavior.

Evaluation of off-target and on-target scoring algorithms and integration into the guide RNA selection tool CRISPOR.

BACKGROUND: The success of the CRISPR/Cas9 genome editing technique depends on the choice of the guide RNA sequence, which is facilitated by various websites. Despite the importance and popularity of these algorithms, it is unclear to which extent their predictions are in agreement with actual measurements. RESULTS: We conduct the first independent evaluation of CRISPR/Cas9 predictions. To this end, we collect data from eight SpCas9 off-target studies and compare them with the sites predicted by popular algorithms. We identify problems in one implementation but found that sequence-based off-target predictions are very reliable, identifying most off-targets with mutation rates superior to 0.1 %, while the number of false positives can be largely reduced with a cutoff on the off-target score. We also evaluate on-target efficiency prediction algorithms against available datasets. The correlation between the predictions and the guide activity varied considerably, especially for zebrafish. Together with novel data from our labs, we find that the optimal on-target efficiency prediction model strongly depends on whether the guide RNA is expressed from a U6 promoter or transcribed in vitro. We further demonstrate that the best predictions can significantly reduce the time spent on guide screening. CONCLUSIONS: To make these guidelines easily accessible to anyone planning a CRISPR genome editing experiment, we built a new website ( http://crispor.org ) that predicts off-targets and helps select and clone efficient guide sequences for more than 120 genomes using different Cas9 proteins and the eight efficiency scoring systems evaluated here.