Genome-wide quantification of RNA flow across subcellular compartments reveals determinants of the mammalian transcript life cycle.

Dissecting the regulatory mechanisms controlling mammalian transcripts from production to degradation requires quantitative measurements of mRNA flow across the cell. We developed subcellular TimeLapse-seq to measure the rates at which RNAs are released from chromatin, exported from the nucleus, loaded onto polysomes, and degraded within the nucleus and cytoplasm in human and mouse cells. These rates varied substantially, yet transcripts from genes with related functions or targeted by the same transcription factors and RNA-binding proteins flowed across subcellular compartments with similar kinetics. Verifying these associations uncovered a link between DDX3X and nuclear export. For hundreds of RNA metabolism genes, most transcripts with retained introns were degraded by the nuclear exosome, while the remaining molecules were exported with stable cytoplasmic lifespans. Transcripts residing on chromatin for longer had extended poly(A) tails, whereas the reverse was observed for cytoplasmic mRNAs. Finally, machine learning identified molecular features that predicted the diverse life cycles of mRNAs.

Saccharomyces cerevisiae Mus81-Mms4 prevents accelerated senescence in telomerase-deficient cells.

Alternative lengthening of telomeres (ALT) in human cells is a conserved process that is often activated in telomerase-deficient human cancers. This process exploits components of the recombination machinery to extend telomere ends, thus allowing for increased proliferative potential. Human MUS81 (Mus81 in Saccharomyces cerevisiae) is the catalytic subunit of structure-selective endonucleases involved in recombination and has been implicated in the ALT mechanism. However, it is unclear whether MUS81 activity at the telomere is specific to ALT cells or if it is required for more general aspects of telomere stability. In this study, we use S. cerevisiae to evaluate the contribution of the conserved Mus81-Mms4 endonuclease in telomerase-deficient yeast cells that maintain their telomeres by mechanisms akin to human ALT. Similar to human cells, we find that yeast Mus81 readily localizes to telomeres and its activity is important for viability after initial loss of telomerase. Interestingly, our analysis reveals that yeast Mus81 is not required for the survival of cells undergoing recombination-mediated telomere lengthening, i.e. for ALT itself. Rather we infer from genetic analysis that Mus81-Mms4 facilitates telomere replication during times of telomere instability. Furthermore, combining mus81 mutants with mutants of a yeast telomere replication factor, Rrm3, reveals that the two proteins function in parallel to promote normal growth during times of telomere stress. Combined with previous reports, our data can be interpreted in a consistent model in which both yeast and human MUS81-dependent nucleases participate in the recovery of stalled replication forks within telomeric DNA. Furthermore, this process becomes crucial under conditions of additional replication stress, such as telomere replication in telomerase-deficient cells.

DNA polymerases δ and λ cooperate in repairing double-strand breaks by microhomology-mediated end-joining in Saccharomyces cerevisiae.

Maintenance of genome stability is carried out by a suite of DNA repair pathways that ensure the repair of damaged DNA and faithful replication of the genome. Of particular importance are the repair pathways, which respond to DNA double-strand breaks (DSBs), and how the efficiency of repair is influenced by sequence homology. In this study, we developed a genetic assay in diploid Saccharomyces cerevisiae cells to analyze DSBs requiring microhomologies for repair, known as microhomology-mediated end-joining (MMEJ). MMEJ repair efficiency increased concomitant with microhomology length and decreased upon introduction of mismatches. The central proteins in homologous recombination (HR), Rad52 and Rad51, suppressed MMEJ in this system, suggesting a competition between HR and MMEJ for the repair of a DSB. Importantly, we found that DNA polymerase delta (Pol δ) is critical for MMEJ, independent of microhomology length and base-pairing continuity. MMEJ recombinants showed evidence that Pol δ proofreading function is active during MMEJ-mediated DSB repair. Furthermore, mutations in Pol δ and DNA polymerase 4 (Pol λ), the DNA polymerase previously implicated in MMEJ, cause a synergistic decrease in MMEJ repair. Pol λ showed faster kinetics associating with MMEJ substrates following DSB induction than Pol δ. The association of Pol δ depended on RAD1, which encodes the flap endonuclease needed to cleave MMEJ intermediates before DNA synthesis. Moreover, Pol δ recruitment was diminished in cells lacking Pol λ. These data suggest cooperative involvement of both polymerases in MMEJ.

Landscape of target:guide homology effects on Cas9-mediated cleavage.

To study target sequence specificity, selectivity, and reaction kinetics of Streptococcus pyogenes Cas9 activity, we challenged libraries of random variant targets with purified Cas9::guide RNA complexes in vitro. Cleavage kinetics were nonlinear, with a burst of initial activity followed by slower sustained cleavage. Consistent with other recent analyses of Cas9 sequence specificity, we observe considerable (albeit incomplete) impairment of cleavage for targets mutated in the PAM sequence or in 'seed' sequences matching the proximal 8 bp of the guide. A second target region requiring close homology was located at the other end of the guide::target duplex (positions 13-18 relative to the PAM). Sequences flanking the guide+PAM region had measurable (albeit modest) effects on cleavage. In addition, the first-base Guanine constraint commonly imposed by gRNA expression systems has little effect on overall cleavage efficiency. Taken together, these studies provide an in vitro understanding of the complexities of Cas9-gRNA interaction and cleavage beyond the general paradigm of site determination based on the 'seed' sequence and PAM.

Target-dependent nickase activities of the CRISPR-Cas nucleases Cpf1 and Cas9.

Clustered regularly interspaced short palindromic repeats (CRISPR) machineries are prokaryotic immune systems that have been adapted as versatile gene editing and manipulation tools. We found that CRISPR nucleases from two families, Cpf1 (also known as Cas12a) and Cas9, exhibit differential guide RNA (gRNA) sequence requirements for cleavage of the two strands of target DNA in vitro. As a consequence of the differential gRNA requirements, both Cas9 and Cpf1 enzymes can exhibit potent nickase activities on an extensive class of mismatched double-stranded DNA (dsDNA) targets. These properties allow the production of efficient nickases for a chosen dsDNA target sequence, without modification of the nuclease protein, using gRNAs with a variety of patterns of mismatch to the intended DNA target. In parallel to the nicking activities observed with purified Cas9 in vitro, we observed sequence-dependent nicking for both perfectly matched and partially mismatched target sequences in a Saccharomyces cerevisiae system. Our findings have implications for CRISPR spacer acquisition, off-target potential of CRISPR gene editing/manipulation, and tool development using homology-directed nicking.

Cooperation between non-essential DNA polymerases contributes to genome stability in Saccharomyces cerevisiae.

DNA polymerases influence genome stability through their involvement in DNA replication, response to DNA damage, and DNA repair processes. Saccharomyces cerevisiae possess four non-essential DNA polymerases, Pol λ, Pol η, Pol ζ, and Rev1, which have varying roles in genome stability. In order to assess the contribution of the non-essential DNA polymerases in genome stability, we analyzed the pol4Δ rev1Δ rev3Δ rad30Δ quadruple mutant in microhomology mediated repair, due to recent studies linking some of these DNA polymerases to this repair pathway. Our results suggest that the length and quality of microhomology influence both the overall efficiency of repair and the involvement of DNA polymerases. Furthermore, the non-essential DNA polymerases demonstrate overlapping and redundant functions when repairing double-strand breaks using short microhomologies containing mismatches. Then, we examined genome-wide mutation accumulation in the pol4Δ rev1Δ rev3Δ rad30Δ quadruple mutant compared to wild type cells. We found a significant decrease in the overall rate of mutation accumulation in the quadruple mutant cells compared to wildtype, but an increase in frameshift mutations and a shift towards transversion base-substitution with a preference for G:C to T:A or C:G. Thus, the non-essential DNA polymerases have an impact on the nature of the mutational spectrum. The sequence and functional homology shared between human and S. cerevisiae non-essential DNA polymerases suggest these DNA polymerases may have a similar role in human cells.

A Small RNA Isolation and Sequencing Protocol and Its Application to Assay CRISPR RNA Biogenesis in Bacteria.

Next generation high-throughput sequencing has enabled sensitive and unambiguous analysis of RNA populations in cells. Here, we describe a method for isolation and strand-specific sequencing of small RNA pools from bacteria that can be multiplexed to accommodate multiple biological samples in a single experiment. Small RNAs are isolated by polyacrylamide gel electrophoresis and treated with T4 polynucleotide kinase. This allows for 3' adapter ligation to CRISPR RNAs, which don't have pre-existing 3'-OH ends. Pre-adenylated adapters are then ligated using T4 RNA ligase 1 in the absence of ATP and with a high concentration of polyethylene glycol (PEG). The 3' capture step enables precise determination of the 3' ends of diverse RNA molecules. Additionally, a random hexamer in the ligated adapter helps control for potential downstream amplification bias. Following reverse-transcription, the cDNA product is circularized and libraries are prepared by PCR. We show that the amplified library need not be visible by gel electrophoresis for efficient sequencing of the desired product. Using this method, we routinely prepare RNA sequencing libraries from minute amounts of purified small RNA. This protocol is tailored to assay for CRISPR RNA biogenesis in bacteria through sequencing of mature CRISPR RNAs, but can be used to sequence diverse classes of small RNAs. We also provide a fully worked example of our data processing pipeline, with instructions for running the provided scripts.

High-Throughput Characterization of Cascade type I-E CRISPR Guide Efficacy Reveals Unexpected PAM Diversity and Target Sequence Preferences.

Interactions between Clustered Regularly Interspaced Short Palindromic Repeat (CRISPR) RNAs and CRISPR-associated (Cas) proteins form an RNA-guided adaptive immune system in prokaryotes. The adaptive immune system utilizes segments of the genetic material of invasive foreign elements in the CRISPR locus. The loci are transcribed and processed to produce small CRISPR RNAs (crRNAs), with degradation of invading genetic material directed by a combination of complementarity between RNA and DNA and in some cases recognition of adjacent motifs called PAMs (Protospacer Adjacent Motifs). Here we describe a general, high-throughput procedure to test the efficacy of thousands of targets, applying this to the Escherichia coli type I-E Cascade (CRISPR-associated complex for antiviral defense) system. These studies were followed with reciprocal experiments in which the consequence of CRISPR activity was survival in the presence of a lytic phage. From the combined analysis of the Cascade system, we found that (i) type I-E Cascade PAM recognition is more expansive than previously reported, with at least 22 distinct PAMs, with many of the noncanonical PAMs having CRISPR-interference abilities similar to the canonical PAMs; (ii) PAM positioning appears precise, with no evidence for tolerance to PAM slippage in interference; and (iii) while increased guanine-cytosine (GC) content in the spacer is associated with higher CRISPR-interference efficiency, high GC content (>62.5%) decreases CRISPR-interference efficiency. Our findings provide a comprehensive functional profile of Cascade type I-E interference requirements and a method to assay spacer efficacy that can be applied to other CRISPR-Cas systems.

On the Origin of Reverse Transcriptase-Using CRISPR-Cas Systems and Their Hyperdiverse, Enigmatic Spacer Repertoires.

Cas1 integrase is the key enzyme of the clustered regularly interspaced short palindromic repeat (CRISPR)-Cas adaptation module that mediates acquisition of spacers derived from foreign DNA by CRISPR arrays. In diverse bacteria, the cas1 gene is fused (or adjacent) to a gene encoding a reverse transcriptase (RT) related to group II intron RTs. An RT-Cas1 fusion protein has been recently shown to enable acquisition of CRISPR spacers from RNA. Phylogenetic analysis of the CRISPR-associated RTs demonstrates monophyly of the RT-Cas1 fusion, and coevolution of the RT and Cas1 domains. Nearly all such RTs are present within type III CRISPR-Cas loci, but their phylogeny does not parallel the CRISPR-Cas type classification, indicating that RT-Cas1 is an autonomous functional module that is disseminated by horizontal gene transfer and can function with diverse type III systems. To compare the sequence pools sampled by RT-Cas1-associated and RT-lacking CRISPR-Cas systems, we obtained samples of a commercially grown cyanobacterium-Arthrospira platensis Sequencing of the CRISPR arrays uncovered a highly diverse population of spacers. Spacer diversity was particularly striking for the RT-Cas1-containing type III-B system, where no saturation was evident even with millions of sequences analyzed. In contrast, analysis of the RT-lacking type III-D system yielded a highly diverse pool but reached a point where fewer novel spacers were recovered as sequencing depth was increased. Matches could be identified for a small fraction of the non-RT-Cas1-associated spacers, and for only a single RT-Cas1-associated spacer. Thus, the principal source(s) of the spacers, particularly the hypervariable spacer repertoire of the RT-associated arrays, remains unknown.IMPORTANCE While the majority of CRISPR-Cas immune systems adapt to foreign genetic elements by capturing segments of invasive DNA, some systems carry reverse transcriptases (RTs) that enable adaptation to RNA molecules. From analysis of available bacterial sequence data, we find evidence that RT-based RNA adaptation machinery has been able to join with CRISPR-Cas immune systems in many, diverse bacterial species. To investigate whether the abilities to adapt to DNA and RNA molecules are utilized for defense against distinct classes of invaders in nature, we sequenced CRISPR arrays from samples of commercial-scale open-air cultures of Arthrospira platensis, a cyanobacterium that contains both RT-lacking and RT-containing CRISPR-Cas systems. We uncovered a diverse pool of naturally occurring immune memories, with the RT-lacking locus acquiring a number of segments matching known viral or bacterial genes, while the RT-containing locus has acquired spacers from a distinct sequence pool for which the source remains enigmatic.

Rad59 regulates association of Rad52 with DNA double-strand breaks.

Homologous recombination among repetitive sequences is an important mode of DNA repair in eukaryotes following acute radiation exposure. We have developed an assay in Saccharomyces cerevisiae that models how multiple DNA double-strand breaks form chromosomal translocations by a nonconservative homologous recombination mechanism, single-strand annealing, and identified the Rad52 paralog, Rad59, as an important factor. We show through genetic and molecular analyses that Rad59 possesses distinct Rad52-dependent and -independent functions, and that Rad59 plays a critical role in the localization of Rad52 to double-strand breaks. Our analysis further suggests that Rad52 and Rad59 act in multiple, sequential processes that determine genome structure following acute exposure to DNA damaging agents.