Cas13d-mediated isoform-specific RNA knockdown with a unified computational and experimental toolbox.

Alternative splicing is an essential mechanism for diversifying proteins, in which mature RNA isoforms produce proteins with potentially distinct functions. Two major challenges in characterizing the cellular function of isoforms are the lack of experimental methods to specifically and efficiently modulate isoform expression and computational tools for complex experimental design. To address these gaps, we developed and methodically tested a strategy which pairs the RNA-targeting CRISPR/Cas13d system with guide RNAs that span exon-exon junctions in the mature RNA. We performed a high-throughput essentiality screen, quantitative RT-PCR assays, and PacBio long read sequencing to affirm our ability to specifically target and robustly knockdown individual RNA isoforms. In parallel, we provide computational tools for experimental design and screen analysis. Considering all possible splice junctions annotated in GENCODE for multi-isoform genes and our gRNA efficacy predictions, we estimate that our junction-centric strategy can uniquely target up to 89% of human RNA isoforms, including 50,066 protein-coding and 11,415 lncRNA isoforms. Importantly, this specificity spans all splicing and transcriptional events, including exon skipping and inclusion, alternative 5' and 3' splice sites, and alternative starts and ends.

Proximity-dependent recruitment of Polycomb repressive complexes by the lncRNA Airn.

During mouse embryogenesis, expression of the long non-coding RNA (lncRNA) Airn leads to gene repression and recruitment of Polycomb repressive complexes (PRCs) to varying extents over a 15-Mb domain. The mechanisms remain unclear. Using high-resolution approaches, we show in mouse trophoblast stem cells that Airn expression induces long-range changes to chromatin architecture that coincide with PRC-directed modifications and center around CpG island promoters that contact the Airn locus even in the absence of Airn expression. Intensity of contact between the Airn lncRNA and chromatin correlated with underlying intensity of PRC recruitment and PRC-directed modifications. Deletion of CpG islands that contact the Airn locus altered long-distance repression and PRC activity in a manner that correlated with changes in chromatin architecture. Our data imply that the extent to which Airn expression recruits PRCs to chromatin is controlled by DNA regulatory elements that modulate proximity of the Airn lncRNA product to its target DNA.

Using RNA Sequencing and Spike-in RNAs to Measure Intracellular Abundance of lncRNAs and mRNAs.

Long noncoding RNAs (lncRNAs) play essential roles in normal physiology and in disease but their mechanisms of action can be challenging to identify. For mechanistic studies, it is often useful to know a lncRNA's intracellular abundance, i.e., approximately how many molecules of the lncRNA are present in a typical cell of a cell-type of interest. At least two approaches have been used to approximate lncRNA intracellular abundance: single-molecule sensitivity RNA fluorescence in situ hybridization (smFISH) and single-gene, calibrated reverse-transcription followed by quantitative PCR (RT-qPCR). However, like all experimental approaches, these methods have their limitations. smFISH, when analyzed using diffraction-limited microscopy, can underestimate intracellular abundance, especially for lncRNAs that accumulate in focused subcellular regions. Calibrated RT-qPCR may return inaccurate estimates of abundance because individual PCR amplicons spaced across the length of a transcript can vary in their efficiency of reverse transcription. Here, we describe a sequencing-based approach that is straightforward, orthogonal to smFISH and RT-qPCR, and can be used to approximate the intracellular abundance for most expressed long RNAs (lncRNAs and mRNAs) in a cell type of interest. Firstly, the average weight of total RNA per cell for the cell type of interest is estimated by replicate rounds of RNA purification from a known number of cells. Secondly, an rRNA-depletion RNA-Seq protocol is performed after adding spike-in control RNAs to a known quantity of total cellular RNA. Lastly, by comparing read counts per transcript to a standard curve derived from the spiked-in RNAs, the intracellular abundance for each transcript is estimated. The sequencing-based approach provides a powerful complement to existing methods, particularly in situations where it is desirable to quantify the abundance of multiple lncRNAs and/or mRNAs simultaneously.

Elements at the 5' end of Xist harbor SPEN-independent transcriptional antiterminator activity.

The Xist lncRNA requires Repeat A, a conserved RNA element located in its 5' end, to induce gene silencing during X-chromosome inactivation. Intriguingly, Repeat A is also required for production of Xist. While silencing by Repeat A requires the protein SPEN, how Repeat A promotes Xist production remains unclear. We report that in mouse embryonic stem cells, expression of a transgene comprising the first two kilobases of Xist (Xist-2kb) causes transcriptional readthrough of downstream polyadenylation sequences. Readthrough required Repeat A and the ∼750 nucleotides downstream, did not require SPEN, and was attenuated by splicing. Despite associating with SPEN and chromatin, Xist-2kb did not robustly silence transcription, whereas a 5.5-kb Xist transgene robustly silenced transcription and read through its polyadenylation sequence. Longer, spliced Xist transgenes also induced robust silencing yet terminated efficiently. Thus, in contexts examined here, Xist requires sequence elements beyond its first two kilobases to robustly silence transcription, and the 5' end of Xist harbors SPEN-independent transcriptional antiterminator activity that can repress proximal cleavage and polyadenylation. In endogenous contexts, this antiterminator activity may help produce full-length Xist RNA while rendering the Xist locus resistant to silencing by the same repressive complexes that the lncRNA recruits to other genes.

lncRNA-Induced Spread of Polycomb Controlled by Genome Architecture, RNA Abundance, and CpG Island DNA.

Long noncoding RNAs (lncRNAs) cause Polycomb repressive complexes (PRCs) to spread over broad regions of the mammalian genome. We report that in mouse trophoblast stem cells, the Airn and Kcnq1ot1 lncRNAs induce PRC-dependent chromatin modifications over multi-megabase domains. Throughout the Airn-targeted domain, the extent of PRC-dependent modification correlated with intra-nuclear distance to the Airn locus, preexisting genome architecture, and the abundance of Airn itself. Specific CpG islands (CGIs) displayed characteristics indicating that they nucleate the spread of PRCs upon exposure to Airn. Chromatin environments surrounding Xist, Airn, and Kcnq1ot1 suggest common mechanisms of PRC engagement and spreading. Our data indicate that lncRNA potency can be tightly linked to lncRNA abundance and that within lncRNA-targeted domains, PRCs are recruited to CGIs via lncRNA-independent mechanisms. We propose that CGIs that autonomously recruit PRCs interact with lncRNAs and their associated proteins through three-dimensional space to nucleate the spread of PRCs in lncRNA-targeted domains.

A piggyBac-based toolkit for inducible genome editing in mammalian cells.

We describe the development and application of a novel series of vectors that facilitate CRISPR-Cas9-mediated genome editing in mammalian cells, which we call CRISPR-Bac. CRISPR-Bac leverages the piggyBac transposon to randomly insert CRISPR-Cas9 components into mammalian genomes. In CRISPR-Bac, a single piggyBac cargo vector containing a doxycycline-inducible Cas9 or catalytically dead Cas9 (dCas9) variant and a gene conferring resistance to Hygromycin B is cotransfected with a plasmid expressing the piggyBac transposase. A second cargo vector, expressing a single-guide RNA (sgRNA) of interest, the reverse-tetracycline TransActivator (rtTA), and a gene conferring resistance to G418, is also cotransfected. Subsequent selection on Hygromycin B and G418 generates polyclonal cell populations that stably express Cas9, rtTA, and the sgRNA(s) of interest. We show that CRISPR-Bac can be used to knock down proteins of interest, to create targeted genetic deletions with high efficiency, and to activate or repress transcription of protein-coding genes and an imprinted long noncoding RNA. The ratio of sgRNA-to-Cas9-to-transposase can be adjusted in transfections to alter the average number of cargo insertions into the genome. sgRNAs targeting multiple genes can be inserted in a single transfection. CRISPR-Bac is a versatile platform for genome editing that simplifies the generation of mammalian cells that stably express the CRISPR-Cas9 machinery.

RETRACTED: A 5' fragment of Xist can sequester RNA produced from adjacent genes on chromatin.

Xist requires Repeat-A, a protein-binding module in its first two kilobases (2kb), to repress transcription. We report that when expressed as a standalone transcript in mouse embryonic stem cells (ESCs), the first 2kb of Xist (Xist-2kb) does not induce transcriptional silencing. Instead, Xist-2kb sequesters RNA produced from adjacent genes on chromatin. Sequestration does not spread beyond adjacent genes, requires the same sequence elements in Repeat-A that full-length Xist requires to repress transcription and can be induced by lncRNAs with similar sequence composition to Xist-2kb. We did not detect sequestration by full-length Xist, but we did detect it by mutant forms of Xist with attenuated transcriptional silencing capability. Xist-2kb associated with SPEN, a Repeat-A binding protein required for Xist-induced transcriptional silencing, but SPEN was not necessary for sequestration. Thus, when expressed in mouse ESCs, a 5' fragment of Xist that contains Repeat-A sequesters RNA from adjacent genes on chromatin and associates with the silencing factor SPEN, but it does not induce transcriptional silencing. Instead, Xist-induced transcriptional silencing requires synergy between Repeat-A and additional sequence elements in Xist. We propose that sequestration is mechanistically related to the Repeat-A dependent stabilization and tethering of Xist near actively transcribed regions of chromatin.

Functional classification of long non-coding RNAs by k-mer content.

The functions of most long non-coding RNAs (lncRNAs) are unknown. In contrast to proteins, lncRNAs with similar functions often lack linear sequence homology; thus, the identification of function in one lncRNA rarely informs the identification of function in others. We developed a sequence comparison method to deconstruct linear sequence relationships in lncRNAs and evaluate similarity based on the abundance of short motifs called k-mers. We found that lncRNAs of related function often had similar k-mer profiles despite lacking linear homology, and that k-mer profiles correlated with protein binding to lncRNAs and with their subcellular localization. Using a novel assay to quantify Xist-like regulatory potential, we directly demonstrated that evolutionarily unrelated lncRNAs can encode similar function through different spatial arrangements of related sequence motifs. K-mer-based classification is a powerful approach to detect recurrent relationships between sequence and function in lncRNAs.

A survey of imprinted gene expression in mouse trophoblast stem cells.

Several hundred mammalian genes are expressed preferentially from one parental allele as the result of a process called genomic imprinting. Genomic imprinting is prevalent in extra-embryonic tissue, where it plays an essential role during development. Here, we profiled imprinted gene expression via RNA-Seq in a panel of six mouse trophoblast stem lines, which are ex vivo derivatives of a progenitor population that gives rise to the placental tissue of the mouse. We found evidence of imprinted expression for 48 genes, 31 of which had been described previously as imprinted and 17 of which we suggest as candidate imprinted genes. An equal number of maternally and paternally biased genes were detected. On average, candidate imprinted genes were more lowly expressed and had weaker parent-of-origin biases than known imprinted genes. Several known and candidate imprinted genes showed variability in parent-of-origin expression bias between the six trophoblast stem cell lines. Sixteen of the 48 known and candidate imprinted genes were previously or newly annotated noncoding RNAs and six encoded for a total of 60 annotated microRNAs. Pyrosequencing across our panel of trophoblast stem cell lines returned levels of imprinted expression that were concordant with RNA-Seq measurements for all eight genes examined. Our results solidify trophoblast stem cells as a cell culture-based experimental model to study genomic imprinting, and provide a quantitative foundation upon which to delineate mechanisms by which the process is maintained in the mouse.