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

CRISPR/Cas9-mediated knock-in of BRCA1/2 mutations restores response to olaparib in pancreatic cancer cell lines

Pancreatic cancer is one of the most aggressive diseases with a very poor outcome. Olaparib, a PARP inhibitor, as maintenance therapy showed benefits in patients with metastatic pancreatic adenocarcinoma bearing germline BRCA1/2 mutations. However, germline BRCA mutation has been described in only 4-7% of patients with pancreatic adenocarcinoma. A CRISPR/Cas9-mediated system was used to knock-in the c.763G > T p.(Glu255*) and c.2133C > A p.(Cys711*) mutations in cell lines to obtain truncated BRCA1 and BRCA2 proteins, respectively. A CRISPR/Cas9 ribonucleoprotein complex was assembled for each mutation and transfected into two pancreatic cell lines (T3M4 and Capan-2) and into a breast cancer cell lines (MCF7) as control. BRCA protein levels were significantly decreased in all BRCA-depleted cells (P < 0.05), proving the transfection efficiency of our CRISPR/Cas9 systems. As expected, the calculated olaparib IC50 were significantly reduced for all cell lines harbored BRCA1 or BRCA2 mutations compared to wild-type BRCA1/2 cells (P < 0.01). Furthermore, we observed a higher induction of apoptosis after 72 h olaparib treatment in BRCA-depleted cells than in wild-type cells. This strategy might offer new insights into the management of patients with pancreatic cancer and open up new perspectives based on the in vivo use of CRISPR/Cas9 strategy.

Landscape of enhancer disruption and functional screen in melanoma cells

BACKGROUND

The high mutation rate throughout the entire melanoma genome presents a major challenge in stratifying true driver events from the background mutations. Numerous recurrent non-coding alterations, such as those in enhancers, can shape tumor evolution, thereby emphasizing the importance in systematically deciphering enhancer disruptions in melanoma.

RESULTS

Here, we leveraged 297 melanoma whole-genome sequencing samples to prioritize highly recurrent regions. By performing a genome-scale CRISPR interference (CRISPRi) screen on highly recurrent region-associated enhancers in melanoma cells, we identified 66 significant hits which could have tumor-suppressive roles. These functional enhancers show unique mutational patterns independent of classical significantly mutated genes in melanoma. Target gene analysis for the essential enhancers reveal many known and hidden mechanisms underlying melanoma growth. Utilizing extensive functional validation experiments, we demonstrate that a super enhancer element could modulate melanoma cell proliferation by targeting MEF2A, and another distal enhancer is able to sustain PTEN tumor-suppressive potential via long-range interactions.

CONCLUSIONS

Our study establishes a catalogue of crucial enhancers and their target genes in melanoma growth and progression, and illuminates the identification of novel mechanisms of dysregulation for melanoma driver genes and new therapeutic targeting strategies.

Epilepsy and sudden unexpected death in epilepsy in a mouse model of human SCN1B-linked developmental and epileptic encephalopathy

Voltage-gated sodium channel β1 subunits are essential proteins that regulate excitability. They modulate sodium and potassium currents, function as cell adhesion molecules and regulate gene transcription following regulated intramembrane proteolysis. Biallelic pathogenic variants in SCN1B, encoding β1, are linked to developmental and epileptic encephalopathy 52, with clinical features overlapping Dravet syndrome. A recessive variant, SCN1B-c.265C>T, predicting SCN1B-p.R89C, was homozygous in two children of a non-consanguineous family. One child was diagnosed with Dravet syndrome, while the other had a milder phenotype. We identified an unrelated biallelic SCN1B-c.265C>T patient with a clinically more severe phenotype than Dravet syndrome. We used CRISPR/Cas9 to knock-in SCN1B-p.R89C to the mouse Scn1b locus (Scn1bR89/C89). We then rederived the line on the C57BL/6J background to allow comparisons between Scn1bR89/R89 and Scn1bC89/C89 littermates with Scn1b+/+ and Scn1b-/- mice, which are congenic on C57BL/6J, to determine whether the SCN1B-c.265C>T variant results in loss-of-function. Scn1bC89/C89 mice have normal body weights and ∼20% premature mortality, compared with severely reduced body weight and 100% mortality in Scn1b-/- mice. β1-p.R89C polypeptides are expressed in brain at comparable levels to wild type. In heterologous cells, β1-p.R89C localizes to the plasma membrane and undergoes regulated intramembrane proteolysis similar to wild type. Heterologous expression of β1-p.R89C results in sodium channel α subunit subtype specific effects on sodium current. mRNA abundance of Scn2a, Scn3a, Scn5a and Scn1b was increased in Scn1bC89/C89 somatosensory cortex, with no changes in Scn1a. In contrast, Scn1b-/- mouse somatosensory cortex is haploinsufficient for Scn1a, suggesting an additive mechanism for the severity of the null model via disrupted regulation of another Dravet syndrome gene. Scn1bC89/C89 mice are more susceptible to hyperthermia-induced seizures at post-natal Day 15 compared with Scn1bR89/R89 littermates. EEG recordings detected epileptic discharges in young adult Scn1bC89/C89 mice that coincided with convulsive seizures and myoclonic jerks. We compared seizure frequency and duration in a subset of adult Scn1bC89/C89 mice that had been exposed to hyperthermia at post-natal Day 15 versus a subset that were not hyperthermia exposed. No differences in spontaneous seizures were detected between groups. For both groups, the spontaneous seizure pattern was diurnal, occurring with higher frequency during the dark cycle. This work suggests that the SCN1B-c.265C>T variant does not result in complete loss-of-function. Scn1bC89/C89 mice more accurately model SCN1B-linked variants with incomplete loss-of-function compared with Scn1b-/- mice, which model complete loss-of-function, and thus add to our understanding of disease mechanisms as well as our ability to develop new therapeutic strategies.

The W101C KCNJ5 Mutation Induces Slower Pacing by Constitutively Active GIRK Channels in hiPSC-Derived Cardiomyocytes

Mutations in the KCNJ5 gene, encoding one of the major subunits of cardiac G-protein-gated inwardly rectifying K+ (GIRK) channels, have been recently linked to inherited forms of sinus node dysfunction. Here, the pathogenic mechanism of the W101C KCNJ5 mutation underlying sinus bradycardia in a patient-derived cellular disease model of sinus node dysfunction (SND) was investigated. A human-induced pluripotent stem cell (hiPSCs) line of a mutation carrier was generated, and CRISPR/Cas9-based gene targeting was used to correct the familial mutation as a control line. Both cell lines were further differentiated into cardiomyocytes (hiPSC-CMs) that robustly expressed GIRK channels which underly the acetylcholine-regulated K+ current (IK,ACh). hiPSC-CMs with the W101C KCNJ5 mutation (hiPSCW101C-CM) had a constitutively active IK,ACh under baseline conditions; the application of carbachol was able to increase IK,ACh, further indicating that not all available cardiac GIRK channels were open at baseline. Additionally, hiPSCW101C-CM had a more negative maximal diastolic potential (MDP) and a slower pacing frequency confirming the bradycardic phenotype. Of note, the blockade of the constitutively active GIRK channel with XAF-1407 rescued the phenotype. These results provide further mechanistic insights and may pave the way for the treatment of SND patients with GIRK channel dysfunction.

Neuron Navigator 1 (Nav1) regulates the response to cocaine in mice

Genetic variation accounts for much of the risk for developing a substance use disorder, but the underlying genetic factors and their genetic effector mechanisms are mostly unknown. Inbred mouse strains exhibit substantial and heritable differences in the extent of voluntary cocaine self-administration. Computational genetic analysis of cocaine self-administration data obtained from twenty-one inbred strains identified Nav1, a member of the neuron navigator family that regulates dendrite formation and axonal guidance, as a candidate gene. To test this genetic hypothesis, we generated and characterized Nav1 knockout mice. Consistent with the genetic prediction, Nav1 knockout mice exhibited increased voluntary cocaine intake and had increased motivation for cocaine consumption. Immunohistochemistry, electrophysiology, and transcriptomic studies were performed as a starting point for investigating the mechanism for the Nav1 knockout effect. Nav1 knockout mice had a reduced inhibitory synapse density in their cortex, increased excitatory synaptic transmission in their cortex and hippocampus, and increased excitatory neurons in a deep cortical layer. Collectively, our results indicate that Nav1 regulates the response to cocaine, and we identified Nav1 knockout induced changes in the excitatory and inhibitory synaptic balance in the cortex and hippocampus that could contribute to this effect.

Inactive S. aureus Cas9 downregulates alpha-synuclein and reduces mtDNA damage and oxidative stress levels in human stem cell model of Parkinson's disease

Parkinson's disease (PD) is one of the most common neurodegenerative diseases, but no disease modifying therapies have been successful in clinical translation presenting a major unmet medical need. A promising target is alpha-synuclein or its aggregated form, which accumulates in the brain of PD patients as Lewy bodies. While it is not entirely clear which alpha-synuclein protein species is disease relevant, mere overexpression of alpha-synuclein in hereditary forms leads to neurodegeneration. To specifically address gene regulation of alpha-synuclein, we developed a CRISPR interference (CRISPRi) system based on the nuclease dead S. aureus Cas9 (SadCas9) fused with the transcriptional repressor domain Krueppel-associated box to controllably repress alpha-synuclein expression at the transcriptional level. We screened single guide (sg)RNAs across the SNCA promoter and identified several sgRNAs that mediate downregulation of alpha-synuclein at varying levels. CRISPRi downregulation of alpha-synuclein in iPSC-derived neuronal cultures from a patient with an SNCA genomic triplication showed functional recovery by reduction of oxidative stress and mitochondrial DNA damage. Our results are proof-of-concept in vitro for precision medicine by targeting the SNCA gene promoter. The SNCA CRISPRi approach presents a new model to understand safe levels of alpha-synuclein downregulation and a novel therapeutic strategy for PD and related alpha-synucleinopathies.

Prediction of CRISPR-Cas9 off-target activities with mismatches and indels based on hybrid neural network

The CRISPR/Cas9 system has significantly advanced the field of gene editing, yet its clinical application is constrained by the considerable challenge of off-target effects. Although numerous deep learning models for off-target prediction have been proposed, most struggle to effectively extract the nuanced features of guide RNA (gRNA) and DNA sequence pairs and to mitigate information loss during data transmission within the model. To address these limitations, we introduce a novel Hybrid Neural Network (HNN) model that employs a parallelized network structure to fully extract pertinent features from different positions and types of bases in the sequence to minimize information loss. Notably, this study marks the first application of word embedding techniques to extract information from sequence pairs that contain insertions and deletions (Indels). Comprehensive evaluation across diverse datasets indicates that our proposed model outperforms existing state-of-the-art prediction methods in off-target prediction. The datasets and source codes supporting this study can be found at https://github.com/Yang-k955/CRISPR-HW.

RNA-binding deficient TDP-43 drives cognitive decline in a mouse model of TDP-43 proteinopathy

TDP-43 proteinopathies including frontotemporal lobar degeneration (FTLD) and amyotrophic lateral sclerosis (ALS) are neurodegenerative disorders characterized by aggregation and mislocalization of the nucleic acid-binding protein TDP-43 and subsequent neuronal dysfunction. Here, we developed endogenous models of sporadic TDP-43 proteinopathy based on the principle that disease-associated TDP-43 acetylation at lysine 145 (K145) alters TDP-43 conformation, impairs RNA-binding capacity, and induces downstream mis-regulation of target genes. Expression of acetylation-mimic TDP-43K145Q resulted in stress-induced nuclear TDP-43 foci and loss of TDP-43 function in primary mouse and human-induced pluripotent stem cell (hiPSC)-derived cortical neurons. Mice harboring the TDP-43K145Q mutation recapitulated key hallmarks of FTLD, including progressive TDP-43 phosphorylation and insolubility, TDP-43 mis-localization, transcriptomic and splicing alterations, and cognitive dysfunction. Our study supports a model in which TDP-43 acetylation drives neuronal dysfunction and cognitive decline through aberrant splicing and transcription of critical genes that regulate synaptic plasticity and stress response signaling. The neurodegenerative cascade initiated by TDP-43 acetylation recapitulates many aspects of human FTLD and provides a new paradigm to further interrogate TDP-43 proteinopathies.

Reduced nephron endowment in the common Six2-TGC tg mouse line is due to Six3 misexpression by aberrant enhancer-promoter interactions in the transgene

Lifelong kidney function relies on the complement of nephrons generated during mammalian development from a mesenchymal nephron progenitor cell (NPC) population. Low nephron endowment confers increased susceptibility to chronic kidney disease. We asked whether reduced nephron numbers in the popular Six2TGC transgenic mouse line 1 was due to disruption of a regulatory gene at the integration site or to ectopic expression of a gene(s) contained within the transgene. Targeted locus amplification identified integration of the Six2TGC transgene within an intron of Cntnap5a on chr1. We generated Hi-C datasets from NPCs isolated from the Six2TGC tg/+ mice, the Cited1 CreERT2/+ control mice, and the Six2TGC tg/+ ; Tsc1 +/Flox,2 mice that exhibited restored nephron number compared with Six2TGC tg/+ mice, and mapped the precise integration of Six2TGC and Cited1 CreERT2 transgenes to chr1 and chr14, respectively. No changes in topology, accessibility, or expression were observed within the 50-megabase region centered on Cntnap5a in Six2TGC tg/+ mice compared with control mice. By contrast, we identified an aberrant regulatory interaction between a Six2 distal enhancer and the Six3 promoter contained within the transgene. Increasing the Six2TGC tg to Six2 locus ratio or removing one Six2 allele in Six2TGC tg/+ mice, caused severe renal hypoplasia. Furthermore, CRISPR disruption of Six3 within the transgene ( Six2TGC Δ Six3CT ) restored nephron endowment to wildtype levels and abolished the stoichiometric effect. Data from genetic and biochemical studies together suggest that in Six2TGC, SIX3 interferes with SIX2 function in NPC renewal through its C-terminal domain.

Significance

Using high-resolution chromatin conformation and accessibility datasets we mapped the integration site of two popular transgenes used in studies of nephron progenitor cells and kidney development. Aberrant enhancer-promoter interactions drive ectopic expression of Six3 in the Six2TGC tg line which was correlated with disruption of nephrogenesis. Disruption of Six3 within the transgene restored nephron numbers to control levels; further genetic and biochemical studies suggest that Six3 interferes with Six2 -mediated regulation of NPC renewal.

HDAC11 inhibition triggers bimodal thermogenic pathways to circumvent adipocyte catecholamine resistance

Stimulation of adipocyte β-adrenergic receptors (β-ARs) induces expression of uncoupling protein 1 (UCP1), promoting nonshivering thermogenesis. Association of β-ARs with a lysine-myristoylated form of A kinase-anchoring protein 12 (AKAP12, also known as gravin-α) is required for downstream signaling that culminates in UCP1 induction. Conversely, demyristoylation of gravin-α by histone deacetylase 11 (HDAC11) suppresses this pathway. Whether inhibition of HDAC11 in adipocytes is sufficient to drive UCP1 expression independently of β-ARs is not known. Here, we demonstrate that adipocyte-specific deletion of HDAC11 in mice leads to robust induction of UCP1 in adipose tissue (AT), resulting in increased body temperature. These effects are mimicked by treating mice in vivo or human AT ex vivo with an HDAC11-selective inhibitor, FT895. FT895 triggers biphasic, gravin-α myristoylation-dependent induction of UCP1 protein expression, with a noncanonical acute response that is posttranscriptional and independent of protein kinase A (PKA), and a delayed response requiring PKA activity and new Ucp1 mRNA synthesis. Remarkably, HDAC11 inhibition promotes UCP1 expression even in models of adipocyte catecholamine resistance where β-AR signaling is blocked. These findings define cell-autonomous, multimodal roles for HDAC11 as a suppressor of thermogenesis, and highlight the potential of inhibiting HDAC11 to therapeutically alter AT phenotype independently of β-AR stimulation.

ATF7IP2/MCAF2 directs H3K9 methylation and meiotic gene regulation in the male germline

H3K9 tri-methylation (H3K9me3) plays emerging roles in gene regulation, beyond its accumulation on pericentric constitutive heterochromatin. It remains a mystery why and how H3K9me3 undergoes dynamic regulation in male meiosis. Here, we identify a novel, critical regulator of H3K9 methylation and spermatogenic heterochromatin organization: the germline-specific protein ATF7IP2 (MCAF2). We show that, in male meiosis, ATF7IP2 amasses on autosomal and X pericentric heterochromatin, spreads through the entirety of the sex chromosomes, and accumulates on thousands of autosomal promoters and retrotransposon loci. On the sex chromosomes, which undergo meiotic sex chromosome inactivation (MSCI), the DNA damage response pathway recruits ATF7IP2 to X pericentric heterochromatin, where it facilitates the recruitment of SETDB1, a histone methyltransferase that catalyzes H3K9me3. In the absence of ATF7IP2, male germ cells are arrested in meiotic prophase I. Analyses of ATF7IP2-deficient meiosis reveal the protein's essential roles in the maintenance of MSCI, suppression of retrotransposons, and global upregulation of autosomal genes. We propose that ATF7IP2 is a downstream effector of the DDR pathway in meiosis that coordinates the organization of heterochromatin and gene regulation through the spatial regulation of SETDB1-mediated H3K9me3 deposition.

A fusion framework of deep learning and machine learning for predicting sgRNA cleavage efficiency

CRISPR/Cas9 system is a powerful tool for genome editing. Numerous studies have shown that sgRNAs can strongly affect the efficiency of editing. However, it is still not clear what rules should be followed for designing sgRNA with high cleavage efficiency. At present, several machine learning or deep learning methods have been developed to predict the cleavage efficiency of sgRNAs, however, the prediction accuracy of these tools is still not satisfactory. Here we propose a fusion framework of deep learning and machine learning, which first deals with the primary sequence and secondary structure features of the sgRNAs using both convolutional neural network (CNN) and recurrent neural network (RNN), and then uses the features extracted by the deep neural network to train a conventional machine learning model with LGBM. As a result, the new approach overwhelmed previous methods. The Spearman's correlation coefficient between predicted and measured sgRNA cleavage efficiency of our model (0.917) is improved by over 5% compared with the most advanced method (0.865), and the mean square error reduces from 7.89 × 10-3 to 4.75 × 10-3. Finally, we developed an online tool, CRISep (http://www.cuilab.cn/CRISep), to evaluate the availability of sgRNAs based on our models.

Myb-like transcription factors have epistatic effects on circadian clock function but additive effects on plant growth

The functions of closely related Myb-like repressor and Myb-like activator proteins within the plant circadian oscillator have been well-studied as separate groups, but the genetic interactions between them are less clear. We hypothesized that these repressors and activators would interact additively to regulate both circadian and growth phenotypes. We used CRISPR-Cas9 to generate new mutant alleles and performed physiological and molecular characterization of plant mutants for five of these core Myb-like clock factors compared with a repressor mutant and an activator mutant. We first examined circadian clock function in plants likely null for both the repressor proteins, CIRCADIAN CLOCK ASSOCIATED 1 (CCA1) and LATE ELONGATED HYPOCOTYL (LHY), and the activator proteins, REVEILLE 4 (RVE4), REVEILLE (RVE6), and REVEILLE (RVE8). The rve468 triple mutant has a long period and flowers late, while cca1 lhy rve468 quintuple mutants, similarly to cca1 lhy mutants, have poor circadian rhythms and flower early. This suggests that CCA1 and LHY are epistatic to RVE4, RVE6, and RVE8 for circadian clock and flowering time function. We next examined hypocotyl elongation and rosette leaf size in these mutants. The cca1 lhy rve468 mutants have growth phenotypes intermediate between cca1 lhy and rve468 mutants, suggesting that CCA1, LHY, RVE4, RVE6, and RVE8 interact additively to regulate growth. Together, our data suggest that these five Myb-like factors interact differently in regulation of the circadian clock versus growth. More generally, the near-norm al seedling phenotypes observed in the largely arrhythmic quintuple mutant demonstrate that circadian-regulated output processes, like control of hypocotyl elongation, do not always depend upon rhythmic oscillator function.

Precise Correction of Lhcgr Mutation in Stem Leydig Cells by Prime Editing Rescues Hereditary Primary Hypogonadism in Mice

Hereditary primary hypogonadism (HPH), caused by gene mutation related to testosterone synthesis in Leydig cells, usually impairs male sexual development and spermatogenesis. Genetically corrected stem Leydig cells (SLCs) transplantation may provide a new approach for treating HPH. Here, a novel nonsense-point-mutation mouse model (LhcgrW495X ) is first generated based on a gene mutation relative to HPH patients. To verify the efficacy and feasibility of SLCs transplantation in treating HPH, wild-type SLCs are transplanted into LhcgrW495X mice, in which SLCs obviously rescue HPH phenotypes. Through comparing several editing strategies, optimized PE2 protein (PEmax) system is identified as an efficient and precise approach to correct the pathogenic point mutation in Lhcgr. Furthermore, delivering intein-split PEmax system via lentivirus successfully corrects the mutation in SLCs from LhcgrW495X mice ex vivo. Gene-corrected SLCs from LhcgrW495X mice exert ability to differentiate into functional Leydig cells in vitro. Notably, the transplantation of gene-corrected SLCs effectively regenerates Leydig cells, recovers testosterone production, restarts sexual development, rescues spermatogenesis, and produces fertile offspring in LhcgrW495X mice. Altogether, these results suggest that PE-based gene editing in SLCs ex vivo is a promising strategy for HPH therapy and is potentially leveraged to address more hereditary diseases in reproductive system.

Interpretable neural architecture search and transfer learning for understanding CRISPR/Cas9 off-target enzymatic reactions

Finely-tuned enzymatic pathways control cellular processes, and their dysregulation can lead to disease. Creating predictive and interpretable models for these pathways is challenging because of the complexity of the pathways and of the cellular and genomic contexts. Here we introduce Elektrum, a deep learning framework which addresses these challenges with data-driven and biophysically interpretable models for determining the kinetics of biochemical systems. First, it uses in vitro kinetic assays to rapidly hypothesize an ensemble of high-quality Kinetically Interpretable Neural Networks (KINNs) that predict reaction rates. It then employs a novel transfer learning step, where the KINNs are inserted as intermediary layers into deeper convolutional neural networks, fine-tuning the predictions for reaction-dependent in vivo outcomes. Elektrum makes effective use of the limited, but clean in vitro data and the complex, yet plentiful in vivo data that captures cellular context. We apply Elektrum to predict CRISPR-Cas9 off-target editing probabilities and demonstrate that Elektrum achieves state-of-the-art performance, regularizes neural network architectures, and maintains physical interpretability.

Tissue-specific knockout in Drosophila neuromuscular system reveals ESCRT's role in formation of synapse-derived extracellular vesicles

Tissue-specific gene knockout by CRISPR/Cas9 is a powerful approach for characterizing gene functions in animal development. However, this approach has been successfully applied in only a small number of Drosophila tissues. The Drosophila motor nervous system is an excellent model system for studying the biology of neuromuscular junction (NMJ). To expand tissue-specific CRISPR to the Drosophila motor system, here we present a CRISPR-mediated tissue-restricted mutagenesis (CRISPR-TRiM) toolkit for knocking out genes in motoneurons, muscles, and glial cells. We validated the efficacy of this toolkit by knocking out known genes in each tissue, demonstrated its orthogonal use with the Gal4/UAS binary expression system, and showed simultaneous knockout of multiple redundant genes. Using these tools, we discovered an essential role for SNARE pathways in NMJ maintenance. Furthermore, we demonstrate that the canonical ESCRT pathway suppresses NMJ bouton growth by downregulating the retrograde Gbb signaling. Lastly, we found that axon termini of motoneurons rely on ESCRT-mediated intra-axonal membrane trafficking to lease extracellular vesicles at the NMJ.

Benchmarking deep learning methods for predicting CRISPR/Cas9 sgRNA on- and off-target activities

In silico design of single guide RNA (sgRNA) plays a critical role in clustered regularly interspaced, short palindromic repeats/CRISPR-associated protein 9 (CRISPR/Cas9) system. Continuous efforts are aimed at improving sgRNA design with efficient on-target activity and reduced off-target mutations. In the last 5 years, an increasing number of deep learning-based methods have achieved breakthrough performance in predicting sgRNA on- and off-target activities. Nevertheless, it is worthwhile to systematically evaluate these methods for their predictive abilities. In this review, we conducted a systematic survey on the progress in prediction of on- and off-target editing. We investigated the performances of 10 mainstream deep learning-based on-target predictors using nine public datasets with different sample sizes. We found that in most scenarios, these methods showed superior predictive power on large- and medium-scale datasets than on small-scale datasets. In addition, we performed unbiased experiments to provide in-depth comparison of eight representative approaches for off-target prediction on 12 publicly available datasets with various imbalanced ratios of positive/negative samples. Most methods showed excellent performance on balanced datasets but have much room for improvement on moderate- and severe-imbalanced datasets. This study provides comprehensive perspectives on CRISPR/Cas9 sgRNA on- and off-target activity prediction and improvement for method development.

Endogenous expression of inactive lysine deacetylases reveals deacetylation-dependent cellular mechanisms

Acetylation of lysine residues is an important and common post-translational regulatory mechanism occurring on thousands of non-histone proteins. Lysine deacetylases (KDACs or HDACs) are a family of enzymes responsible for removing acetylation. To identify the biological mechanisms regulated by individual KDACs, we created HT1080 cell lines containing chromosomal point mutations, which endogenously express either KDAC6 or KDAC8 having single inactivated catalytic domain. Engineered HT1080 cells expressing inactive KDA6 or KDAC8 domains remained viable and exhibited enhanced acetylation on known substrate proteins. RNA-seq analysis revealed that many changes in gene expression were observed when KDACs were inactivated, and that these gene sets differed significantly from knockdown and knockout cell lines. Using GO ontology, we identified several critical biological processes associated specifically with catalytic activity and others attributable to non-catalytic interactions. Treatment of wild-type cells with KDAC-specific inhibitors Tubastatin A and PCI-34051 resulted in gene expression changes distinct from those of the engineered cell lines, validating this approach as a tool for evaluating in-cell inhibitor specificity and identifying off-target effects of KDAC inhibitors. Probing the functions of specific KDAC domains using these cell lines is not equivalent to doing so using previously existing methods and provides novel insight into the catalytic functions of individual KDACs by investigating the molecular and cellular changes upon genetic inactivation.

Compound heterozygous IFT81 variations in a skeletal ciliopathy patient cause Bardet-Biedl syndrome-like ciliary defects

Owing to their crucial roles in development and homeostasis, defects in cilia cause ciliopathies with diverse clinical manifestations. The intraflagellar transport (IFT) machinery, containing the IFT-A and IFT-B complexes, mediates not only the intraciliary bidirectional trafficking but also import and export of ciliary proteins together with the kinesin-2 and dynein-2 motor complexes. The BBSome, containing eight subunits encoded by causative genes of Bardet-Biedl syndrome (BBS), connects the IFT machinery to ciliary membrane proteins to mediate their export from cilia. Although mutations in subunits of the IFT-A and dynein-2 complexes cause skeletal ciliopathies, mutations in some IFT-B subunits are also known to cause skeletal ciliopathies. We here show that compound heterozygous variations of an IFT-B subunit, IFT81, found in a patient with skeletal ciliopathy cause defects in its interactions with other IFT-B subunits, and in ciliogenesis and ciliary protein trafficking when one of the two variants was expressed in IFT81-knockout (KO) cells. Notably, we found that IFT81-KO cells expressing IFT81(Δ490-519), which lacks the binding site for the IFT25-IFT27 dimer, causes ciliary defects reminiscent of those found in BBS cells and those in IFT74-KO cells expressing a BBS variant of IFT74, which forms a heterodimer with IFT81. In addition, IFT81-KO cells expressing IFT81(Δ490-519) in combination with the other variant, IFT81 (L645*), which mimics the cellular conditions of the above skeletal ciliopathy patient, demonstrated essentially the same phenotype as those expressing only IFT81(Δ490-519). Thus, our data indicate that BBS-like defects can be caused by skeletal ciliopathy variants of IFT81.

Using CRISPR/Cas9 to identify genes required for mechanosensory neuron development and function

Tunicates are marine, non-vertebrate chordates that comprise the sister group to the vertebrates. Most tunicates have a biphasic lifecycle that alternates between a swimming larva and a sessile adult. Recent advances have shed light on the neural basis for the tunicate larva's ability to sense a proper substrate for settlement and initiate metamorphosis. Work in the highly tractable laboratory model tunicate Ciona robusta suggests that sensory neurons embedded in the anterior papillae transduce mechanosensory stimuli to trigger larval tail retraction and initiate the process of metamorphosis. Here, we take advantage of the low-cost and simplicity of Ciona by using tissue-specific CRISPR/Cas9-mediated mutagenesis to screen for genes potentially involved in mechanosensation and metamorphosis, in the context of an undergraduate 'capstone' research course. This small screen revealed at least one gene, Vamp1/2/3, which appears crucial for the ability of the papillae to trigger metamorphosis. We also provide step-by-step protocols and tutorials associated with this course, in the hope that it might be replicated in similar CRISPR-based laboratory courses wherever Ciona are available.

Genetic repair of a human induced pluripotent cell line from patient with Dutch-type cerebral amyloid angiopathy

Dutch-type cerebral amyloid angiopathy (D-CAA), also known as hereditary cerebral haemorrhage with amyloidosis-Dutch type (HCHWA-D), is an autosomal dominant disorder caused by a G to C transversion in codon 693 of the amyloid precursor protein (APP) that results in a Gln-to-Glu amino acid substitution. CRISPR-Cas9 editing was used for genetic correction of the mutation in a human induced pluripotent stem cell (hiPSC-) line established previously. The isogenic hiPSCs generated showed typical pluripotent stem cell morphology, expressed all markers of undifferentiated state, displayed a normal karyotype and had the capacity to differentiate into the three germ layers.

A streamlined guide RNA screening system for genome editing in Sorghum bicolor

BACKGROUND

Genome editing tools derived from clustered regularly interspaced short palindromic repeats (CRISPR) systems have been developed for generating targeted mutations in plants. Although these tools hold promise for rapid crop improvement, target-specific guide RNAs exhibit variable activity. To improve genome editing, a rapid and precise method for evaluating their efficiency is necessary.

RESULTS

Here we report an efficient system for screening single guide RNAs (sgRNAs) for genome editing in sorghum using a transient protoplast transfection assay. Protoplasts were isolated from leaves from sorghum plants cultivated under three different conditions. Cultivation for three days of continuous darkness following seven days with a 16-h light and 8-h dark photoperiod resulted in the highest yield of viable protoplasts and the highest protoplast transfection efficiency. We tested both plasmid-mediated and ribonucleoprotein-based delivery to protoplasts, via polyethylene glycol-mediated transfection, of CRISPR components targeting the sorghum genome. The frequencies of small insertions and deletions induced by a set of sgRNAs targeting four endogenous sorghum genes were analyzed via targeted deep sequencing. Our screening system induced indels in sorghum protoplasts at frequencies of up to 77.8% (plasmid) and 18.5% (RNP). The entire screening system was completed within 16 days.

CONCLUSIONS

The screening system optimized in this study for predicting sgRNA activity for genome editing in sorghum is efficient and straightforward. This system will reduce the time and effort needed for sorghum genome editing.

ATF3 induction prevents precocious activation of skeletal muscle stem cell by regulating H2B expression

Skeletal muscle stem cells (also called satellite cells, SCs) are important for maintaining muscle tissue homeostasis and damage-induced regeneration. However, it remains poorly understood how SCs enter cell cycle to become activated upon injury. Here we report that AP-1 family member ATF3 (Activating Transcription Factor 3) prevents SC premature activation. Atf3 is rapidly and transiently induced in SCs upon activation. Short-term deletion of Atf3 in SCs accelerates acute injury-induced regeneration, however, its long-term deletion exhausts the SC pool and thus impairs muscle regeneration. The Atf3 loss also provokes SC activation during voluntary exercise and enhances the activation during endurance exercise. Mechanistically, ATF3 directly activates the transcription of Histone 2B genes, whose reduction accelerates nucleosome displacement and gene transcription required for SC activation. Finally, the ATF3-dependent H2B expression also prevents genome instability and replicative senescence in SCs. Therefore, this study has revealed a previously unknown mechanism for preserving the SC population by actively suppressing precocious activation, in which ATF3 is a key regulator.

Mutagenesis of the melon Prv gene by CRISPR/Cas9 breaks papaya ringspot virus resistance and generates an autoimmune allele with constitutive defense responses

The majority of plant disease resistance (R) genes encode nucleotide binding-leucine-rich repeat (NLR) proteins. In melon, two closely linked NLR genes, Fom-1 and Prv, were mapped and identified as candidate genes that control resistance to Fusarium oxysporum f.sp. melonis races 0 and 2, and to papaya ringspot virus (PRSV), respectively. In this study, we validated the function of Prv and showed that it is essential for providing resistance against PRSV infection. We generated CRISPR/Cas9 [clustered regularly interspaced palindromic repeats (CRISPR)/CRISPR-associated protein 9] mutants using Agrobacterium-mediated transformation of a PRSV-resistant melon genotype, and the T1 progeny proved susceptible to PRSV, showing strong disease symptoms and viral spread upon infection. Three alleles having 144, 154, and ~3 kb deletions, respectively, were obtained, all of which caused loss of resistance. Interestingly, one of the Prv mutant alleles, prvΔ154, encoding a truncated product, caused an extreme dwarf phenotype, accompanied by leaf lesions, high salicylic acid levels, and defense gene expression. The autoimmune phenotype observed at 25 °C proved to be temperature dependent, being suppressed at 32 °C. This is a first report on the successful application of CRISPR/Cas9 to confirm R gene function in melon. Such validation opens up new opportunities for molecular breeding of disease resistance in this important vegetable crop.

Converging evidence from exome sequencing and common variants implicates target genes for osteoporosis

In this study, we leveraged the combined evidence of rare coding variants and common alleles to identify therapeutic targets for osteoporosis. We undertook a large-scale multiancestry exome-wide association study for estimated bone mineral density, which showed that the burden of rare coding alleles in 19 genes was associated with estimated bone mineral density (P < 3.6 × 10-7). These genes were highly enriched for a set of known causal genes for osteoporosis (65-fold; P = 2.5 × 10-5). Exome-wide significant genes had 96-fold increased odds of being the top ranked effector gene at a given GWAS locus (P = 1.8 × 10-10). By integrating proteomics Mendelian randomization evidence, we prioritized CD109 (cluster of differentiation 109) as a gene for which heterozygous loss of function is associated with higher bone density. CRISPR-Cas9 editing of CD109 in SaOS-2 osteoblast-like cell lines showed that partial CD109 knockdown led to increased mineralization. This study demonstrates that the convergence of common and rare variants, proteomics and CRISPR can highlight new bone biology to guide therapeutic development.

Sperm chromatin structure and reproductive fitness are altered by substitution of a single amino acid in mouse protamine 1

Conventional dogma presumes that protamine-mediated DNA compaction in sperm is achieved by electrostatic interactions between DNA and the arginine-rich core of protamines. Phylogenetic analysis reveals several non-arginine residues conserved within, but not across species. The significance of these residues and their post-translational modifications are poorly understood. Here, we investigated the role of K49, a rodent-specific lysine residue in protamine 1 (P1) that is acetylated early in spermiogenesis and retained in sperm. In sperm, alanine substitution (P1(K49A)) decreases sperm motility and male fertility-defects that are not rescued by arginine substitution (P1(K49R)). In zygotes, P1(K49A) leads to premature male pronuclear decompaction, altered DNA replication, and embryonic arrest. In vitro, P1(K49A) decreases protamine-DNA binding and alters DNA compaction and decompaction kinetics. Hence, a single amino acid substitution outside the P1 arginine core is sufficient to profoundly alter protein function and developmental outcomes, suggesting that protamine non-arginine residues are essential for reproductive fitness.

Pre-mRNA splicing order is predetermined and maintains splicing fidelity across multi-intronic transcripts

Combinatorially, intron excision within a given nascent transcript could proceed down any of thousands of paths, each of which would expose different dynamic landscapes of cis-elements and contribute to alternative splicing. In this study, we found that post-transcriptional multi-intron splicing order in human cells is largely predetermined, with most genes spliced in one or a few predominant orders. Strikingly, these orders were conserved across cell types and stages of motor neuron differentiation. Introns flanking alternatively spliced exons were frequently excised last, after their neighboring introns. Perturbations to the spliceosomal U2 snRNA altered the preferred splicing order of many genes, and these alterations were associated with the retention of other introns in the same transcript. In one gene, early removal of specific introns was sufficient to induce delayed excision of three proximal introns, and this delay was caused by two distinct cis-regulatory mechanisms. Together, our results demonstrate that multi-intron splicing order in human cells is predetermined, is influenced by a component of the spliceosome and ensures splicing fidelity across long pre-mRNAs.

Optimized minimal genome-wide human sgRNA library

Genome-wide clustered regularly interspaced short palindromic repeats (CRISPR)-based knockout screening is revolting the genetic analysis of a cellular or molecular phenotype in question but is challenged by the large size of single-guide RNA (sgRNA) library. Here we designed a minimal genome-wide human sgRNA library, H-mLib, which is composed of 21,159 sgRNA pairs assembled based on a dedicated selection strategy from all potential SpCas9/sgRNAs in the human genome. These sgRNA pairs were cloned into a dual-gRNA vector each targeting one gene, resulting in a compact library size nearly identical to the number of human protein-coding genes. The performance of the H-mLib was benchmarked to other CRISPR libraries in a proliferation screening conducted in K562 cells. We also identified groups of core essential genes and cell-type specific essential genes by comparing the screening results from the K562 and Jurkat cells. Together, the H-mLib exemplified high specificity and sensitivity in identifying essential genes while containing minimal library complexity, emphasizing its advantages and applications in CRISPR screening with limited cell numbers.

Dynein-2-driven intraciliary retrograde trafficking indirectly requires multiple interactions of IFT54 in the IFT-B complex with the dynein-2 complex

Within cilia, the dynein-2 complex needs to be transported as an anterograde cargo to achieve its role as a motor to drive retrograde trafficking of the intraflagellar transport (IFT) machinery containing IFT-A and IFT-B complexes. We previously showed that interactions of WDR60 and the DYNC2H1-DYNC2LI1 dimer of dynein-2 with multiple IFT-B subunits, including IFT54, are required for the trafficking of dynein-2 as an IFT cargo. However, specific deletion of the IFT54-binding site from WDR60 demonstrated only a minor effect on dynein-2 trafficking and function. We here show that the C-terminal coiled-coil region of IFT54, which participates in its interaction with the DYNC2H1-DYNC2LI1 dimer of dynein-2 and with IFT20 of the IFT-B complex, is essential for IFT-B function, and suggest that the IFT54 middle linker region between the N-terminal WDR60-binding region and the C-terminal coiled-coil is required for ciliary retrograde trafficking, probably by mediating the effective binding of IFT-B to the dynein-2 complex, and thereby ensuring dynein-2 loading onto the anterograde IFT trains. The results presented here agree with the notion predicted from the previous structural models that the dynein-2 loading onto the anterograde IFT train relies on intricate, multivalent interactions between the dynein-2 and IFT-B complexes.

Identification of a psychiatric risk gene NISCH at 3p21.1 GWAS locus mediating dendritic spine morphogenesis and cognitive function

BACKGROUND

Schizophrenia and bipolar disorder (BD) are believed to share clinical symptoms, genetic risk, etiological factors, and pathogenic mechanisms. We previously reported that single nucleotide polymorphisms spanning chromosome 3p21.1 showed significant associations with both schizophrenia and BD, and a risk SNP rs2251219 was in linkage disequilibrium with a human specific Alu polymorphism rs71052682, which showed enhancer effects on transcriptional activities using luciferase reporter assays in U251 and U87MG cells.

METHODS

CRISPR/Cas9-directed genome editing, real-time quantitative PCR, and public Hi-C data were utilized to investigate the correlation between the Alu polymorphism rs71052682 and NISCH. Primary neuronal culture, immunofluorescence staining, co-immunoprecipitation, lentiviral vector production, intracranial stereotaxic injection, behavioral assessment, and drug treatment were used to examine the physiological impacts of Nischarin (encoded by NISCH).

RESULTS

Deleting the Alu sequence in U251 and U87MG cells reduced mRNA expression of NISCH, the gene locates 180 kb from rs71052682, and Hi-C data in brain tissues confirmed the extensive chromatin contacts. These data suggested that the genetic risk of schizophrenia and BD predicted elevated NISCH expression, which was also consistent with the observed higher NISCH mRNA levels in the brain tissues from psychiatric patients compared with controls. We then found that overexpression of NISCH resulted in a significantly decreased density of mushroom dendritic spines with a simultaneously increased density of thin dendritic spines in primary cultured neurons. Intriguingly, elevated expression of this gene in mice also led to impaired spatial working memory in the Y-maze. Given that Nischarin is the target of anti-hypertensive agents clonidine and tizanidine, which have shown therapeutic effects in patients with schizophrenia and patients with BD in preliminary clinical trials, we demonstrated that treatment with those antihypertensive drugs could reduce NISCH mRNA expression and rescue the impaired working memory in mice.

CONCLUSIONS

We identify a psychiatric risk gene NISCH at 3p21.1 GWAS locus influencing dendritic spine morphogenesis and cognitive function, and Nischarin may have potentials for future therapeutic development.

Recyclable CRISPR/Cas9 mediated gene disruption and deletions in Histoplasma

Targeted gene disruption is challenging in the dimorphic fungal pathogen Histoplasma due to the low frequency of homologous recombination. Transformed DNA is either integrated ectopically into the genome or maintained extra chromosomally by de novo addition of telomeric sequences. Based on a system developed in Blastomyces, we adapted a CRISPR/Cas9 system to facilitate targeted gene disruption in Histoplasma with high efficiency. We express a codon-optimized version of Cas9 as well as guide RNAs from a single ectopic vector carrying a selectable marker. Once the desired mutation is verified, one can screen for isolates that have lost the Cas9 vector by simply removing the selective pressure. Multiple mutations can then be generated in the same strain by retransforming the Cas9 vector carrying different guides. We used this system to disrupt a number of target genes including RYP2 and SRE1 where loss-of-function mutations could be monitored visually by colony morphology or color, respectively. Interestingly, expression of two guide RNAs targeting the 5' and 3' ends of a gene allowed isolation of deletion mutants where the sequence between the guide RNAs was removed from the genome. Whole-genome sequencing showed that the frequency of off-target mutations associated with the Cas9 nuclease was negligible. Finally, we increased the frequency of gene disruption by using an endogenous Histoplasma regulatory sequence to drive guide RNA expression. These tools transform our ability to generate targeted mutations in Histoplasma.

Golgi satellites are essential for polysialylation of NCAM and expression of LTP at distal synapses

The complex cytoarchitecture of neurons poses significant challenges for the maturation of synaptic membrane proteins. It is currently unclear whether locally secreted synaptic proteins bypass the Golgi or whether they traffic through Golgi satellites (GSs). Here, we create a transgenic GS reporter mouse line and show that GSs are widely distributed along dendrites and are capable of mature glycosylation, in particular sialylation. We find that polysialylation of locally secreted NCAM takes place at GSs. Accordingly, in mice lacking a component of trans-Golgi network-to-plasma membrane trafficking, we find fewer GSs and significantly reduced PSA-NCAM levels in distal dendrites of CA1 neurons that receive input from the temporoammonic pathway. Induction of long-term potentiation at those, but not more proximal, synapses is severely impaired. We conclude that GSs serve the need for local mature glycosylation of synaptic membrane proteins in distal dendrites and thereby contribute to rapid changes in synaptic strength.

Barley stripe mosaic virus-mediated somatic and heritable gene editing in barley (Hordeum vulgare L.)

Genome editing strategies in barley (Hordeum vulgare L.) typically rely on Agrobacterium-mediated genetic transformation for the delivery of required genetic reagents involving tissue culture techniques. These approaches are genotype-dependent, time-consuming, and labor-intensive, which hampers rapid genome editing in barley. More recently, plant RNA viruses have been engineered to transiently express short guide RNAs facilitating CRISPR/Cas9-based targeted genome editing in plants that constitutively express Cas9. Here, we explored virus-induced genome editing (VIGE) based on barley stripe mosaic virus (BSMV) in Cas9-transgenic barley. Somatic and heritable editing in the ALBOSTRIANS gene (CMF7) resulting in albino/variegated chloroplast-defective barley mutants is shown. In addition, somatic editing in meiosis-related candidate genes in barley encoding ASY1 (an axis-localized HORMA domain protein), MUS81 (a DNA structure-selective endonuclease), and ZYP1 (a transverse filament protein of the synaptonemal complex) was achieved. Hence, the presented VIGE approach using BSMV enables rapid somatic and also heritable targeted gene editing in barley.

Secondary structures that regulate mRNA translation provide insights for ASO-mediated modulation of cardiac hypertrophy

Translation of upstream open reading frames (uORFs) typically abrogates translation of main (m)ORFs. The molecular mechanism of uORF regulation in cells is not well understood. Here, we identified a double-stranded RNA (dsRNA) structure residing within the GATA4 uORF that augments uORF translation and inhibits mORF translation. Antisense oligonucleotides (ASOs) that disrupt this dsRNA structure promote mORF translation, while ASOs that base-pair immediately downstream (i.e., forming a bimolecular double-stranded region) of either the uORF or mORF start codon enhance uORF or mORF translation, respectively. Human cardiomyocytes and mice treated with a uORF-enhancing ASO showed reduced cardiac GATA4 protein levels and increased resistance to cardiomyocyte hypertrophy. We further show the general utility of uORF-dsRNA- or mORF- targeting ASO to regulate mORF translation for other mRNAs. Our work demonstrates a regulatory paradigm that controls translational efficiency and a useful strategy to alter protein expression and cellular phenotypes by targeting or generating dsRNA downstream of a uORF or mORF start codon.

Bullet points for discoveries

dsRNA within GATA4 uORF activates uORF translation and inhibits mORF translation. ASOs that target the dsRNA can either inhibit or enhance GATA4 mORF translation. ASOs can be used to impede hypertrophy in human cardiomyocytes and mouse hearts.uORF- and mORF-targeting ASOs can be used to control translation of multiple mRNAs.

A phloem-localized Arabidopsis metacaspase (AtMC3) improves drought tolerance

Increasing drought phenomena pose a serious threat to agricultural productivity. Although plants have multiple ways to respond to the complexity of drought stress, the underlying mechanisms of stress sensing and signaling remain unclear. The role of the vasculature, in particular the phloem, in facilitating inter-organ communication is critical and poorly understood. Combining genetic, proteomic and physiological approaches, we investigated the role of AtMC3, a phloem-specific member of the metacaspase family, in osmotic stress responses in Arabidopsis thaliana. Analyses of the proteome in plants with altered AtMC3 levels revealed differential abundance of proteins related to osmotic stress pointing into a role of the protein in water-stress-related responses. Overexpression of AtMC3 conferred drought tolerance by enhancing the differentiation of specific vascular tissues and maintaining higher levels of vascular-mediated transportation, while plants lacking the protein showed an impaired response to drought and inability to respond effectively to the hormone abscisic acid. Overall, our data highlight the importance of AtMC3 and vascular plasticity in fine-tuning early drought responses at the whole plant level without affecting growth or yield.

High-throughput sgRNA testing reveals rules for Cas9 specificity and DNA repair in tomato cells

CRISPR/Cas9 technology has the potential to significantly enhance plant breeding. To determine the specificity and the mutagenic spectrum of SpCas9 in tomato, we designed 89 g(uide) RNAs targeting genes of the tomato MYB transcription factor family with varying predicted specificities. Plasmids encoding sgRNAs and Cas9 were introduced into tomato protoplasts, and target sites as well as 224 predicted off-target sites were screened for the occurrence of mutations using amplicon sequencing. Algorithms for the prediction of efficacy of the sgRNAs had little predictive power in this system. The analysis of mutations suggested predictable identity of single base insertions. Off-target mutations were found for 13 out of 89 sgRNAs and only occurred at positions with one or two mismatches (at 14 and 3 sites, respectively). We found that PAM-proximal mismatches do not preclude low frequency off-target mutations. Off-target mutations were not found at all 138 positions that had three or four mismatches. We compared off-target mutation frequencies obtained with plasmid encoding sgRNAs and Cas9 with those induced by ribonucleoprotein (RNP) transfections. The use of RNPs led to a significant decrease in relative off-target frequencies at 6 out of 17, no significant difference at 9, and an increase at 2 sites. Additionally, we show that off-target sequences with insertions or deletions relative to the sgRNA may be mutated, and should be considered during sgRNA design. Altogether, our data help sgRNA design by providing insight into the Cas9-induced double-strand break repair outcomes and the occurrence of off-target mutations.

Detect-seq, a chemical labeling and biotin pull-down approach for the unbiased and genome-wide off-target evaluation of programmable cytosine base editors

Programmable cytosine base editors show promising approaches for correcting pathogenic mutations; yet, their off-target effects have been of great concern. Detect-seq (dU-detection enabled by C-to-T transition during sequencing) is an unbiased, sensitive method for the off-target evaluation of programmable cytosine base editors. It profiles the editome by tracing the editing intermediate dU, which is introduced inside living cells and edited by programmable cytosine base editors. The genomic DNA is extracted, preprocessed and labeled by successive chemical and enzymatic reactions, followed by biotin pull-down to enrich the dU-containing loci for sequencing. Here, we describe a detailed protocol for performing the Detect-seq experiment, and a customized, open-source, bioinformatic pipeline for analyzing the characteristic Detect-seq data is also provided. Unlike those previous whole-genome sequencing-based methods, Detect-seq uses an enrichment strategy and hence is endowed with great sensitivity, a higher signal-to-noise ratio and no requirement for high sequencing depth. Furthermore, Detect-seq is widely applicable for both mitotic and postmitotic biological systems. The entire protocol typically takes 5 d from the genomic DNA extraction to sequencing and ~1 week for data analysis.

CRISPR/Cas9 editing reveals IRF8 regulated gene signatures restraining plasmablast differentiation

The transcription factor Interferon regulatory factor 8 (IRF8) is involved in maintaining B cell identity. However, how IRF8 regulates T cell independent B cell responses are not fully characterized. Here, an in vivo CRISPR/Cas9 system was optimized to generate Irf8-deficient murine B cells and used to determine the role of IRF8 in B cells responding to LPS stimulation. Irf8-deficient B cells more readily formed CD138⁺ plasmablasts in response to LPS with the principal dysregulation occurring at the activated B cell stage. Transcriptional profiling revealed an upregulation of plasma cell associated genes prematurely in activated B cells and a failure to repress the gene expression programs of IRF1 and IRF7 in Irf8-deficient cells. These data expand on the known roles of IRF8 in regulating B cell identity by preventing premature plasma cell formation and highlight how IRF8 helps evolve TLR responses away from the initial activation towards those driving humoral immunity.

A redundant transcription factor network steers spatiotemporal Arabidopsis triterpene synthesis

Plant specialized metabolites modulate developmental and ecological functions and comprise many therapeutic and other high-value compounds. However, the mechanisms determining their cell-specific expression remain unknown. Here we describe the transcriptional regulatory network that underlies cell-specific biosynthesis of triterpenes in Arabidopsis thaliana root tips. Expression of thalianol and marneral biosynthesis pathway genes depends on the phytohormone jasmonate and is limited to outer tissues. We show that this is promoted by the activity of redundant bHLH-type transcription factors from two distinct clades and coactivated by homeodomain factors. Conversely, the DOF-type transcription factor DAG1 and other regulators prevent expression of the triterpene pathway genes in inner tissues. We thus show how precise expression of triterpene biosynthesis genes is determined by a robust network of transactivators, coactivators and counteracting repressors.

Using traditional machine learning and deep learning methods for on- and off-target prediction in CRISPR/Cas9: a review

CRISPR/Cas9 (Clustered Regularly Interspaced Short Palindromic Repeats and CRISPR-associated protein 9) is a popular and effective two-component technology used for targeted genetic manipulation. It is currently the most versatile and accurate method of gene and genome editing, which benefits from a large variety of practical applications. For example, in biomedicine, it has been used in research related to cancer, virus infections, pathogen detection, and genetic diseases. Current CRISPR/Cas9 research is based on data-driven models for on- and off-target prediction as a cleavage may occur at non-target sequence locations. Nowadays, conventional machine learning and deep learning methods are applied on a regular basis to accurately predict on-target knockout efficacy and off-target profile of given single-guide RNAs (sgRNAs). In this paper, we present an overview and a comparative analysis of traditional machine learning and deep learning models used in CRISPR/Cas9. We highlight the key research challenges and directions associated with target activity prediction. We discuss recent advances in the sgRNA-DNA sequence encoding used in state-of-the-art on- and off-target prediction models. Furthermore, we present the most popular deep learning neural network architectures used in CRISPR/Cas9 prediction models. Finally, we summarize the existing challenges and discuss possible future investigations in the field of on- and off-target prediction. Our paper provides valuable support for academic and industrial researchers interested in the application of machine learning methods in the field of CRISPR/Cas9 genome editing.

Epigenome Editing Durability Varies Widely Across Cardiovascular Disease Target Genes

Background

Hepatic knockdown of the proprotein convertase subtilisin/kexin type 9 ( PCSK9 ) gene or the angiopoietin-like 3 ( ANGPTL3 ) gene has been demonstrated to reduce blood low-density lipoprotein cholesterol (LDL-C) levels, and hepatic knockdown of the angiotensinogen ( AGT ) gene has been demonstrated to reduce blood pressure. Genome editing can productively target each of these three genes in hepatocytes in the liver, offering the possibility of durable "one-and-done" therapies for hypercholesterolemia and hypertension. However, concerns around making permanent gene sequence changes via DNA strand breaks might hinder acceptance of these therapies. Epigenome editing offers an alternative approach to gene inactivation, via silencing of gene expression by methylation of the promoter region, but the long-term durability of epigenome editing remains to be established.

Methods

We assessed the ability of epigenome editing to durably reduce the expression of the human PCSK9, ANGPTL3 , and AGT genes in HuH-7 hepatoma cells. Using the CRISPRoff epigenome editor, we identified guide RNAs that produced efficient gene knockdown immediately after transfection. We assessed the durability of gene expression and methylation changes through serial cell passages.

Results

Cells treated with CRISPRoff and PCSK9 guide RNAs were maintained for up to 124 cell doublings and demonstrated durable knockdown of gene expression and increased CpG dinucleotide methylation in the promoter, exon 1, and intron 1 regions. In contrast, cells treated with CRISPRoff and ANGPTL3 guide RNAs experienced only transient knockdown of gene expression. Cells treated with CRISPRoff and AGT guide RNAs also experienced transient knockdown of gene expression; although initially there was increased CpG methylation throughout the early part of the gene, this methylation was geographically heterogeneous-transient in the promoter, and stable in intron 1.

Conclusions

This work demonstrates precise and durable gene regulation via methylation, supporting a new therapeutic approach for protection against cardiovascular disease via knockdown of genes such as PCSK9 . However, the durability of knockdown with methylation changes is not generalizable across target genes, likely limiting the therapeutic potential of epigenome editing compared to other modalities.

Ring Finger 149-Related Is an FGF/MAPK-Independent Regulator of Pharyngeal Muscle Fate Specification

During embryonic development, cell-fate specification gives rise to dedicated lineages that underlie tissue formation. In olfactores, which comprise tunicates and vertebrates, the cardiopharyngeal field is formed by multipotent progenitors of both cardiac and branchiomeric muscles. The ascidian Ciona is a powerful model to study cardiopharyngeal fate specification with cellular resolution, as only two bilateral pairs of multipotent cardiopharyngeal progenitors give rise to the heart and to the pharyngeal muscles (also known as atrial siphon muscles, ASM). These progenitors are multilineage primed, in as much as they express a combination of early ASM- and heart-specific transcripts that become restricted to their corresponding precursors, following oriented and asymmetric divisions. Here, we identify the primed gene ring finger 149 related (Rnf149-r), which later becomes restricted to the heart progenitors, but appears to regulate pharyngeal muscle fate specification in the cardiopharyngeal lineage. CRISPR/Cas9-mediated loss of Rnf149-r function impairs atrial siphon muscle morphogenesis, and downregulates Tbx1/10 and Ebf, two key determinants of pharyngeal muscle fate, while upregulating heart-specific gene expression. These phenotypes are reminiscent of the loss of FGF/MAPK signaling in the cardiopharyngeal lineage, and an integrated analysis of lineage-specific bulk RNA-seq profiling of loss-of-function perturbations has identified a significant overlap between candidate FGF/MAPK and Rnf149-r target genes. However, functional interaction assays suggest that Rnf149-r does not directly modulate the activity of the FGF/MAPK/Ets1/2 pathway. Instead, we propose that Rnf149-r acts both in parallel to the FGF/MAPK signaling on shared targets, as well as on FGF/MAPK-independent targets through (a) separate pathway(s).

Deep sampling of gRNA in the human genome and deep-learning-informed prediction of gRNA activities

Life science studies involving clustered regularly interspaced short palindromic repeat (CRISPR) editing generally apply the best-performing guide RNA (gRNA) for a gene of interest. Computational models are combined with massive experimental quantification on synthetic gRNA-target libraries to accurately predict gRNA activity and mutational patterns. However, the measurements are inconsistent between studies due to differences in the designs of the gRNA-target pair constructs, and there has not yet been an integrated investigation that concurrently focuses on multiple facets of gRNA capacity. In this study, we analyzed the DNA double-strand break (DSB)-induced repair outcomes and measured SpCas9/gRNA activities at both matched and mismatched locations using 926,476 gRNAs covering 19,111 protein-coding genes and 20,268 non-coding genes. We developed machine learning models to forecast the on-target cleavage efficiency (AIdit_ON), off-target cleavage specificity (AIdit_OFF), and mutational profiles (AIdit_DSB) of SpCas9/gRNA from a uniformly collected and processed dataset by deep sampling and massively quantifying gRNA capabilities in K562 cells. Each of these models exhibited superlative performance in predicting SpCas9/gRNA activities on independent datasets when benchmarked with previous models. A previous unknown parameter was also empirically determined regarding the "sweet spot" in the size of datasets used to establish an effective model to predict gRNA capabilities at a manageable experimental scale. In addition, we observed cell type-specific mutational profiles and were able to link nucleotidylexotransferase as the key factor driving these outcomes. These massive datasets and deep learning algorithms have been implemented into the user-friendly web service http://crispr-aidit.com to evaluate and rank gRNAs for life science studies.

Dynamics of Chromatin Accessibility During Hematopoietic Stem Cell Differentiation Into Progressively Lineage-Committed Progeny

Epigenetic mechanisms regulate the multilineage differentiation capacity of hematopoietic stem cells (HSCs) into a variety of blood and immune cells. Mapping the chromatin dynamics of functionally defined cell populations will shed mechanistic insight into 2 major, unanswered questions in stem cell biology: how does epigenetic identity contribute to a cell type's lineage potential, and how do cascades of chromatin remodeling dictate ensuing fate decisions? Our recent work revealed evidence of multilineage gene priming in HSCs, where open cis-regulatory elements (CREs) exclusively shared between HSCs and unipotent lineage cells were enriched for DNA binding motifs of known lineage-specific transcription factors. Oligopotent progenitor populations operating between the HSCs and unipotent cells play essential roles in effecting hematopoietic homeostasis. To test the hypothesis that selective HSC-primed lineage-specific CREs remain accessible throughout differentiation, we used ATAC-seq to map the temporal dynamics of chromatin remodeling during progenitor differentiation. We observed epigenetic-driven clustering of oligopotent and unipotent progenitors into distinct erythromyeloid and lymphoid branches, with multipotent HSCs and MPPs associating with the erythromyeloid lineage. We mapped the dynamics of lineage-primed CREs throughout hematopoiesis and identified both unique and shared CREs as potential lineage reinforcement mechanisms at fate branch points. Additionally, quantification of genome-wide peak count and size revealed overall greater chromatin accessibility in HSCs, allowing us to identify HSC-unique peaks as putative regulators of self-renewal and multilineage potential. Finally, CRISPRi-mediated targeting of ATACseq-identified putative CREs in HSCs allowed us to demonstrate the functional role of selective CREs in lineage-specific gene expression. These findings provide insight into the regulation of stem cell multipotency and lineage commitment throughout hematopoiesis and serve as a resource to test functional drivers of hematopoietic lineage fate.

Target-seq: single workflow for detection of genome integration site, DNA translocation and off-target events

Designed donor DNA delivery through viral or nonviral systems to target loci in the host genome is a critical step for gene therapy. Adeno-associated virus and lentivirus are leading vehicles for in vivo and ex vivo delivery of therapeutic genes due to their high delivery and editing efficiency. Nonviral editing tools, such as CRISPR/Cas9, are getting more attention for gene modification. However, there are safety concerns; for example, tumorigenesis due to off-target effects and DNA rearrangement. Analysis tools to detect and characterize on-target and off-target genome modification post editing in the host genome are pivotal for evaluating the success and safety of gene therapy. We developed Target-seq combined with different analysis tools to detect the genome integration site, DNA translocation and off-target events.

Cytoskeletal Keratins Are Overexpressed in a Zebrafish Model of Idiopathic Scoliosis

Idiopathic scoliosis (IS) is a three-dimensional rotation of the spine >10 degrees with an unknown etiology. Our laboratory established a late-onset IS model in zebrafish (Danio rerio) containing a deletion in kif7. A total of 25% of kif7^co63/co63 zebrafish develop spinal curvatures and are otherwise developmentally normal, although the molecular mechanisms underlying the scoliosis are unknown. To define transcripts associated with scoliosis in this model, we performed bulk mRNA sequencing on 6 weeks past fertilization (wpf) kif7^co63/co63 zebrafish with and without scoliosis. Additionally, we sequenced kif7^co63/co63, kif7^co63/+, and AB zebrafish (n = 3 per genotype). Sequencing reads were aligned to the GRCz11 genome and FPKM values were calculated. Differences between groups were calculated for each transcript by the t-test. Principal component analysis showed that transcriptomes clustered by sample age and genotype. kif7 mRNA was mildly reduced in both homozygous and heterozygous zebrafish compared to AB. Sonic hedgehog target genes were upregulated in kif7^co63/co63 zebrafish over AB, but no difference was detected between scoliotic and non-scoliotic mutants. The top upregulated genes in scoliotic zebrafish were cytoskeletal keratins. Pankeratin staining of 6 wpf scoliotic and non-scoliotic kif7^co63/co63 zebrafish showed increased keratin levels within the zebrafish musculature and intervertebral disc (IVD). Keratins are major components of the embryonic notochord, and aberrant keratin expression has been associated with intervertebral disc degeneration (IVDD) in both zebrafish and humans. The role of increased keratin accumulation as a molecular mechanism associated with the onset of scoliosis warrants further study.

Genome-Editing of FtsZ1 for Alteration of Starch Granule Size in Potato Tubers

Genome-editing has enabled rapid improvement for staple food crops, such as potato, a key beneficiary of the technology. In potato, starch contained within tubers represents the primary product for use in food and non-food industries. Starch granules are produced in the plastids of tubers with plastid size correlated with the size of starch grana. The division of plastids is controlled by proteins, including the tubulin-like GTPase FtsZ1. The altered expression of FtsZ1 has been shown to disrupt plastid division, leading to the production of "macro-plastid"-containing plants. These macro-chloroplast plants are characterized by cells containing fewer and enlarged plastids. In this work, we utilize CRISPR/Cas9 to generate FtsZ1 edited potato lines to demonstrate that genome-editing can be used to increase the size of starch granules in tubers. Altered plastid morphology was comparable to the overexpression of FtsZ1 in previous work in potato and other crops. Several lines were generated with up to a 1.98-fold increase in starch granule size that was otherwise phenotypically indistinguishable from wild-type plants. Further, starch paste from one of the most promising lines showed a 2.07-fold increase in final viscosity. The advantages of enlarged starch granules and the potential of CRISPR/Cas9-based technologies for food crop improvement are further discussed.

POLR1A variants underlie phenotypic heterogeneity in craniofacial, neural, and cardiac anomalies

Heterozygous pathogenic variants in POLR1A, which encodes the largest subunit of RNA Polymerase I, were previously identified as the cause of acrofacial dysostosis, Cincinnati-type. The predominant phenotypes observed in the cohort of 3 individuals were craniofacial anomalies reminiscent of Treacher Collins syndrome. We subsequently identified 17 additional individuals with 12 unique heterozygous variants in POLR1A and observed numerous additional phenotypes including neurodevelopmental abnormalities and structural cardiac defects, in combination with highly prevalent craniofacial anomalies and variable limb defects. To understand the pathogenesis of this pleiotropy, we modeled an allelic series of POLR1A variants in vitro and in vivo. In vitro assessments demonstrate variable effects of individual pathogenic variants on ribosomal RNA synthesis and nucleolar morphology, which supports the possibility of variant-specific phenotypic effects in affected individuals. To further explore variant-specific effects in vivo, we used CRISPR-Cas9 gene editing to recapitulate two human variants in mice. Additionally, spatiotemporal requirements for Polr1a in developmental lineages contributing to congenital anomalies in affected individuals were examined via conditional mutagenesis in neural crest cells (face and heart), the second heart field (cardiac outflow tract and right ventricle), and forebrain precursors in mice. Consistent with its ubiquitous role in the essential function of ribosome biogenesis, we observed that loss of Polr1a in any of these lineages causes cell-autonomous apoptosis resulting in embryonic malformations. Altogether, our work greatly expands the phenotype of human POLR1A-related disorders and demonstrates variant-specific effects that provide insights into the underlying pathogenesis of ribosomopathies.