Highly efficient generation of isogenic pluripotent stem cell models using prime editing.

The recent development of prime editing (PE) genome engineering technologies has the potential to significantly simplify the generation of human pluripotent stem cell (hPSC)-based disease models. PE is a multicomponent editing system that uses a Cas9-nickase fused to a reverse transcriptase (nCas9-RT) and an extended PE guide RNA (pegRNA). Once reverse transcribed, the pegRNA extension functions as a repair template to introduce precise designer mutations at the target site. Here, we systematically compared the editing efficiencies of PE to conventional gene editing methods in hPSCs. This analysis revealed that PE is overall more efficient and precise than homology-directed repair of site-specific nuclease-induced double-strand breaks. Specifically, PE is more effective in generating heterozygous editing events to create autosomal dominant disease-associated mutations. By stably integrating the nCas9-RT into hPSCs we achieved editing efficiencies equal to those reported for cancer cells, suggesting that the expression of the PE components, rather than cell-intrinsic features, limit PE in hPSCs. To improve the efficiency of PE in hPSCs, we optimized the delivery modalities for the PE components. Delivery of the nCas9-RT as mRNA combined with synthetically generated, chemically-modified pegRNAs and nicking guide RNAs improved editing efficiencies up to 13-fold compared with transfecting the PE components as plasmids or ribonucleoprotein particles. Finally, we demonstrated that this mRNA-based delivery approach can be used repeatedly to yield editing efficiencies exceeding 60% and to correct or introduce familial mutations causing Parkinson's disease in hPSCs.

From muscles to nerves, our body is formed of many kinds of cells which can each respond slightly differently to the same harmful genetic changes. Understanding the exact relationship between mutations and cell-type specific function is essential to better grasp how conditions such as Parkinson's disease or amyotrophic lateral sclerosis progress and can be treated. Stem cells could be an important tool in that effort, as they can be directed to mature into many cell types in the laboratory. Yet it remains difficult to precisely introduce disease-relevant mutations in these cells. To remove this obstacle, Li et al. focused on prime editing, a cutting-edge 'search and replace' approach which can introduce new genetic information into a specific DNA sequence. However, it was unclear whether this technique could be used to efficiently create stem cell models of human diseases. A first set of experiments showed that prime editing is superior to conventional approaches when generating mutated genes in stem cells. Li et al. then further improved the efficiency and precision of the method by tweaking how prime editing components are delivered into the cells. The refined approach could be harnessed to quickly generate large numbers of stem cells carrying mutations associated with Parkinson's disease; crucially, prime editing could then also be used to revert a mutated gene back to its healthy form. The improved prime editing approach developed by Li et al. removes a major hurdle for scientists hoping to use stem cells to study genetic diseases. This could potentially help to unlock progress in how we understand and ultimately treat these conditions.

JUM is a computational method for comprehensive annotation-free analysis of alternative pre-mRNA splicing patterns.

Alternative pre-mRNA splicing (AS) greatly diversifies metazoan transcriptomes and proteomes and is crucial for gene regulation. Current computational analysis methods of AS from Illumina RNA-sequencing data rely on preannotated libraries of known spliced transcripts, which hinders AS analysis with poorly annotated genomes and can further mask unknown AS patterns. To address this critical bioinformatics problem, we developed a method called the junction usage model (JUM) that uses a bottom-up approach to identify, analyze, and quantitate global AS profiles without any prior transcriptome annotations. JUM accurately reports global AS changes in terms of the five conventional AS patterns and an additional "composite" category composed of inseparable combinations of conventional patterns. JUM stringently classifies the difficult and disease-relevant pattern of intron retention (IR), reducing the false positive rate of IR detection commonly seen in other annotation-based methods to near-negligible rates. When analyzing AS in RNA samples derived from Drosophila heads, human tumors, and human cell lines bearing cancer-associated splicing factor mutations, JUM consistently identified approximately twice the number of novel AS events missed by other methods. Computational simulations showed JUM exhibits a 1.2 to 4.8 times higher true positive rate at a fixed cutoff of 5% false discovery rate. In summary, JUM provides a framework and improved method that removes the necessity for transcriptome annotations and enables the detection, analysis, and quantification of AS patterns in complex metazoan transcriptomes with superior accuracy.

An Efficient Method for Electroporation of Small Interfering RNAs into ENCODE Project Tier 1 GM12878 and K562 Cell Lines.

The Encyclopedia of DNA Elements (ENCODE) Project aims to identify all functional sequence elements in the human genome sequence by use of high-throughput DNA/cDNA sequencing approaches. To aid the standardization, comparison, and integration of data sets produced from different technologies and platforms, the ENCODE Consortium selected several standard human cell lines to be used by the ENCODE Projects. The Tier 1 ENCODE cell lines include GM12878, K562, and H1 human embryonic stem cell lines. GM12878 is a lymphoblastoid cell line, transformed with the Epstein-Barr virus, that was selected by the International HapMap Project for whole genome and transcriptome sequencing by use of the Illumina platform. K562 is an immortalized myelogenous leukemia cell line. The GM12878 cell line is attractive for the ENCODE Projects, as it offers potential synergy with the International HapMap Project. Despite the vast amount of sequencing data available on the GM12878 cell line through the ENCODE Project, including transcriptome, chromatin immunoprecipitation-sequencing for histone marks, and transcription factors, no small interfering siRNA-mediated knockdown studies have been performed in the GM12878 cell line, as cationic lipid-mediated transfection methods are inefficient for lymphoid cell lines. Here, we present an efficient and reproducible method for transfection of a variety of siRNAs into the GM12878 and K562 cell lines, which subsequently results in targeted protein depletion.

Denaturation and electrophoresis of RNA with formaldehyde.

Electrophoretic size fractionation can be used to denature and separate large mRNA molecules (0.5-10 kb) on formaldehyde-containing agarose gels. Formaldehyde contains a carbonyl group that reacts to form Schiff bases with the imino or amino groups of guanine, adenine, and cytosine. These covalent adducts prevent normal base pairing and maintain the RNA in a denatured state. Because these adducts are unstable, formaldehyde must be present in the gel to maintain the RNA in the denatured state. This protocol describes the preparation of an agarose gel with formaldehyde and its setup in a horizontal electrophoresis apparatus. RNA samples are prepared and denatured in a solution of formamide and formaldehyde and, with 0.5- to 10-kb size markers, subjected to electrophoresis through the gel. Following electrophoresis, the gel is stained to visualize RNA markers or rRNA using one of several different types of stains.

5'-end labeling of RNA with [γ-32P]ATP and T4 polynucleotide kinase.

This protocol uses T4 polynucleotide kinase to catalyze the transfer of a radiolabeled, terminal (γ) phosphate of ATP to the 5'-hydroxyl terminus of a DNA or RNA molecule. The reaction is very efficient and hence is used as a general method for phosphorylating polynucleotides or oligonucleotides.

A segmental genomic duplication generates a functional intron.

An intron is an extended genomic feature whose function requires multiple constrained positions-donor and acceptor splice sites, a branch point, a polypyrimidine tract and suitable splicing enhancers-that may be distributed over hundreds or thousands of nucleotides. New introns are therefore unlikely to emerge by incremental accumulation of functional sub-elements. Here we demonstrate that a functional intron can be created de novo in a single step by a segmental genomic duplication. This experiment recapitulates in vivo the birth of an intron that arose in the ancestral jawed vertebrate lineage nearly half-a-billion years ago.

An SF1 affinity model to identify branch point sequences in human introns.

Splicing factor 1 (SF1) binds to the branch point sequence (BPS) of mammalian introns and is believed to be important for the splicing of some, but not all, introns. To help identify BPSs, particularly those that depend on SF1, we generated a BPS profile model in which SF1 binding affinity data, validated by branch point mapping, were iteratively incorporated into computational models. We searched a data set of 117,499 human introns for best matches to the SF1 Affinity Model above a threshold, and counted the number of matches at each intronic position. After subtracting a background value, we found that 87.9% of remaining high-scoring matches identified were located in a region upstream of 3'-splice sites where BPSs are typically found. Since U2AF65 recognizes the polypyrimidine tract (PPT) and forms a cooperative RNA complex with SF1, we combined the SF1 model with a PPT model computed from high affinity binding sequences for U2AF65. The combined model, together with binding site location constraints, accurately identified introns bound by SF1 that are candidates for SF1-dependent splicing.

Global analysis of alternative splicing regulation by insulin and wingless signaling in Drosophila cells.

BACKGROUND: Despite the prevalence and biological relevance of both signaling pathways and alternative pre-mRNA splicing, our knowledge of how intracellular signaling impacts on alternative splicing regulation remains fragmentary. We report a genome-wide analysis using splicing-sensitive microarrays of changes in alternative splicing induced by activation of two distinct signaling pathways, insulin and wingless, in Drosophila cells in culture. RESULTS: Alternative splicing changes induced by insulin affect more than 150 genes and more than 50 genes are regulated by wingless activation. About 40% of the genes showing changes in alternative splicing also show regulation of mRNA levels, suggesting distinct but also significantly overlapping programs of transcriptional and post-transcriptional regulation. Distinct functional sets of genes are regulated by each pathway and, remarkably, a significant overlap is observed between functional categories of genes regulated transcriptionally and at the level of alternative splicing. Functions related to carbohydrate metabolism and cellular signaling are enriched among genes regulated by insulin and wingless, respectively. Computational searches identify pathway-specific sequence motifs enriched near regulated 5' splice sites. CONCLUSIONS: Taken together, our data indicate that signaling cascades trigger pathway-specific and biologically coherent regulatory programs of alternative splicing regulation. They also reveal that alternative splicing can provide a novel molecular mechanism for crosstalk between different signaling pathways.

Polypyrimidine tract binding protein controls the transition from exon definition to an intron defined spliceosome.

The polypyrimidine tract binding protein (PTB) binds pre-mRNAs to alter splice-site choice. We characterized a series of spliceosomal complexes that assemble on a pre-mRNA under conditions of either PTB-mediated splicing repression or its absence. In the absence of repression, exon definition complexes that were assembled downstream of the regulated exon could progress to pre-spliceosomal A complexes and functional spliceosomes. Under PTB-mediated repression, assembly was arrested at an A-like complex that was unable to transition to spliceosomal complexes. Trans-splicing experiments indicated that, even when the U1 and U2 small nuclear ribonucleoprotein particles (snRNPs) are properly bound to the upstream and downstream exons, the presence of PTB prevents the interaction of the two exon complexes. Proteomic analyses of these complexes provide a new description of exon definition complexes, and indicate that splicing regulators can act on the transition between the exon definition complex and an intron-defined spliceosome.

Identification and analysis of U5 snRNA variants in Drosophila.

Distinct isoforms of spliceosomal RNAs may be involved in regulating pre-messenger RNA splicing in eukaryotic cells. During a large-scale effort to identify small noncoding RNAs in Drosophila, we isolated a U5 snRNA-like molecule containing a 5' segment identical to that of the canonical (major) U5 snRNA but with a variant Sm binding site and a distinct 3' hairpin sequence. Based on this finding, another six similar U5 snRNA-like sequences were identified within the Drosophila genome by sequence similarity to the invariant loop in the 5' half of U5. Interestingly, although all of these variants are expressed in vivo, each shows a distinct temporal expression profile during Drosophila development, and one is expressed primarily in fly heads. The presence of these U5 snRNA variants within RNP particles suggests their role in splicing and implies a possible connection to regulation of developmental and tissue-specific gene expression.

CDK phosphorylation inhibits the DNA-binding and ATP-hydrolysis activities of the Drosophila origin recognition complex.

Faithful propagation of eukaryotic chromosomes usually requires that no DNA segment be replicated more than once during one cell cycle. Cyclin-dependent kinases (Cdks) are critical for the re-replication controls that inhibit the activities of components of the pre-replication complexes (pre-RCs) following origin activation. The origin recognition complex (ORC) initiates the assembly of pre-RCs at origins of replication and Cdk phosphorylation of ORC is important for the prevention of re-initiation. Here we show that Drosophila melanogaster ORC (DmORC) is phosphorylated in vivo and is a substrate for Cdks in vitro. Cdk phosphorylation of DmORC subunits DmOrc1p and DmOrc2p inhibits the intrinsic ATPase activity of DmORC without affecting ATP binding to DmOrc1p. Moreover, Cdk phosphorylation inhibits the ATP-dependent DNA-binding activity of DmORC in vitro, thus identifying a novel determinant for DmORC-DNA interaction. DmORC is a substrate for both Cdk2 x cyclin E and Cdk1 x cyclin B in vitro. Such phosphorylation of DmORC by Cdk2 x cyclin E, but not by Cdk1 x cyclin B, requires an "RXL" motif in DmOrc1p. We also identify casein kinase 2 (CK2) as a kinase activity in embryonic extracts targeting DmORC for modification. CK2 phosphorylation does not affect ATP hydrolysis by DmORC but modulates the ATP-dependent DNA-binding activity of DmORC. These results suggest molecular mechanisms by which Cdks may inhibit ORC function as part of re-replication control and show that DmORC activity may be modulated in response to phosphorylation by multiple kinases.

Regulatory role of dADAR in ROS metabolism in Drosophila CNS.

Pre-mRNA adenosine deaminase (ADAR) is involved in many physiological processes by either directly converting adenosine to inosine in certain pre-mRNAs or indirectly regulating expression of certain genes. Mutations of Drosophila ADAR (dADAR) results in neuronal dysfunction and hypersensitivity to oxygen deprivation. Recently, we found that the mutant flies were very resistant to paraquat, a compound that generates free radicals. In order to further characterize the neuronal role of dADAR and understand the basis for the resistance to the oxidative stress, we investigated the effect of dADAR on the expression of genes encoding scavengers of cellular reactive oxygen species (ROS) in both dADAR mutant and overexpression flies. Our data show that the expression of the genes encoding known ROS scavengers [superoxide dismutase (SOD) and catalase] is not regulated by dADAR. However, the transcripts of genes encoding two potential ROS scavengers (dhd and Cyp4g1) were robustly increased in dADAR mutant flies, and conversely both were significantly decreased in dADAR overexpressing flies. Using dhd [encoding a Drosophila homolog of the mammalian protein thioredoxin (Trx)] transgenic flies, we confirmed that the resistance of dADAR mutant flies to paraquat resulted, at least partially, from the up-regulation of dhd gene in dADAR mutant flies. Our data not only confirm the importance of ADAR in maintenance of neuronal function but also reveal its regulatory role in the expression of genes encoding ROS scavengers.

Genome-wide analysis reveals an unexpected function for the Drosophila splicing factor U2AF50 in the nuclear export of intronless mRNAs.

The protein factor U2AF is an essential component required for pre-mRNA splicing. Mutations identified in the S. pombe large U2AF subunit were used to engineer transgenic Drosophila carrying temperature-sensitive U2AF large subunit alleles. Mutant recombinant U2AF heterodimers showed reduced polypyrimidine tract RNA binding at elevated temperatures. Genome-wide RNA profiling comparing wild-type and mutant strains identified more than 400 genes differentially expressed in the dU2AF50 mutant flies grown at the restrictive temperature. Surprisingly, almost 40% of the downregulated genes lack introns. Microarray analyses revealed that nuclear export of a large number of intronless mRNAs is impaired in Drosophila-cultured cells RNAi knocked down for dU2AF50. Immunopurification of nuclear RNP complexes showed that dU2AF50 associates with intronless mRNAs. These results reveal an unexpected role for the splicing factor dU2AF50 in the nuclear export of intronless mRNAs.

Interplay between Drosophila Bloom's syndrome helicase and Ku autoantigen during nonhomologous end joining repair of P element-induced DNA breaks.

P transposable elements in Drosophila are mobilized via a cut-and-paste mechanism. The broken DNA ends generated during transposition can be repaired via the homology-directed synthesis-dependent strand annealing or by nonhomologous end joining (NHEJ). Genetic studies have demonstrated an interaction between the gene (mus309, for mutagen-sensitive) encoding the Drosophila Bloom's syndrome helicase homolog (DmBLM) and the Ku70 gene, which is involved in NHEJ. We have used RNA interference (RNAi) to knock down expression of DmBLM and one or both of the Drosophila Ku subunits, DmKu70 or DmKu80. Our results show that upon reduction of DmKu, an increase in small deletions (1-49 bp) and large deletions (>/=50 bp) flanking the site of P element-induced breaks is observed, and a reduction in large deletions at these sites is found upon reduction of DmBLM. Moreover, double RNAi of DmKu and DmBLM results in an increase in small deletions characteristic of the DmKu RNAi and also partially suppresses the reduction in repair efficiency observed with DmKu RNAi. These results suggest that there are DNA double-strand break recognition and/or processing events involving DmKu and DmBLM that, when eliminated by RNAi, lead to deletions. Finally, these results raise the possibility that, unlike the situation in mammals, where BLM appears to function exclusively in the homologous repair pathway, in Drosophila, DmBLM may be directly involved in, or at least influence the double-strand break recognition that leads to the NHEJ repair pathway.

Drosophila melanogaster MNK/Chk2 and p53 regulate multiple DNA repair and apoptotic pathways following DNA damage.

We have used genetic and microarray analysis to determine how ionizing radiation (IR) induces p53-dependent transcription and apoptosis in Drosophila melanogaster. IR induces MNK/Chk2-dependent phosphorylation of p53 without changing p53 protein levels, indicating that p53 activity can be regulated without an Mdm2-like activity. In a genome-wide analysis of IR-induced transcription in wild-type and mutant embryos, all IR-induced increases in transcript levels required both p53 and the Drosophila Chk2 homolog MNK. Proapoptotic targets of p53 include hid, reaper, sickle, and the tumor necrosis factor family member EIGER: Overexpression of Eiger is sufficient to induce apoptosis, but mutations in Eiger do not block IR-induced apoptosis. Animals heterozygous for deletions that span the reaper, sickle, and hid genes exhibited reduced IR-dependent apoptosis, indicating that this gene complex is haploinsufficient for induction of apoptosis. Among the genes in this region, hid plays a central, dosage-sensitive role in IR-induced apoptosis. p53 and MNK/Chk2 also regulate DNA repair genes, including two components of the nonhomologous end-joining repair pathway, Ku70 and Ku80. Our results indicate that MNK/Chk2-dependent modification of Drosophila p53 activates a global transcriptional response to DNA damage that induces error-prone DNA repair as well as intrinsic and extrinsic apoptosis pathways.

Widespread predicted nonsense-mediated mRNA decay of alternatively-spliced transcripts of human normal and disease genes.

We have recently shown that a third of reliably-inferred alternative mRNA isoforms are candidates for nonsense-mediated mRNA decay (NMD), an mRNA surveillance system (Lewis et al., 2003; PROC: Natl Acad. Sci. USA, 100, 189-192). Rather than being translated to yield protein, these transcripts are expected to be degraded and may be subject to regulated unproductive splicing and translation (RUST). Our initial experimental studies are consistent with these predictions and suggest an unappreciated role for NMD in several human diseases.