Autism-associated CHD8 controls reactive gliosis and neuroinflammation via remodeling chromatin in astrocytes.

Reactive changes of glial cells during neuroinflammation impact brain disorders and disease progression. Elucidating the mechanisms that control reactive gliosis may help us to understand brain pathophysiology and improve outcomes. Here, we report that adult ablation of autism spectrum disorder (ASD)-associated CHD8 in astrocytes attenuates reactive gliosis via remodeling chromatin accessibility, changing gene expression. Conditional Chd8 deletion in astrocytes, but not microglia, suppresses reactive gliosis by impeding astrocyte proliferation and morphological elaboration. Astrocyte Chd8 ablation alleviates lipopolysaccharide-induced neuroinflammation and septic-associated hypothermia in mice. Astrocytic CHD8 plays an important role in neuroinflammation by altering the chromatin landscape, regulating metabolic and lipid-associated pathways, and astrocyte-microglia crosstalk. Moreover, we show that reactive gliosis can be directly mitigated in vivo using an adeno-associated virus (AAV)-mediated Chd8 gene editing strategy. These findings uncover a role of ASD-associated CHD8 in the adult brain, which may warrant future exploration of targeting chromatin remodelers in reactive gliosis and neuroinflammation in injury and neurological diseases.

Low copy CRISPR-Cas13d mitigates collateral RNA cleavage.

While CRISPR-Cas13 systems excel in accurately targeting RNA, the potential for collateral RNA degradation poses a concern for therapeutic applications and limits broader adoption for transcriptome perturbations. We evaluate the extent to which collateral RNA cleavage occurs when Rfx Cas13d is delivered via plasmid transfection or lentiviral transduction and find that collateral activity only occurs with high levels of Rfx Cas13d expression. Using transcriptome-scale and combinatorial CRISPR pooled screens in cell lines with low-copy Rfx Cas13d, we find high on-target knockdown, without extensive collateral activity regardless of the expression level of the target gene. In contrast, transfection of Rfx Cas13d, which yields higher nuclease expression, results in collateral RNA degradation. Further, our analysis of a high-fidelity Cas13 variant uncovers a marked decrease in on-target efficiency, suggesting that its reduced collateral activity may be due to an overall diminished nuclease capability.

Genome-wide CRISPR screens in spheroid culture reveal that the tumor suppressor LKB1 inhibits growth via the PIKFYVE lipid kinase.

The tumor suppressor LKB1 is a serine/threonine protein kinase that is frequently mutated in human lung adenocarcinoma (LUAD). LKB1 regulates a complex signaling network that is known to control cell polarity and metabolism; however, the pathways that mediate the tumor-suppressive activity of LKB1 are incompletely defined. To identify mechanisms of LKB1-mediated growth suppression, we developed a spheroid-based cell culture assay to study LKB1-dependent growth. We then performed genome-wide CRISPR screens in spheroidal culture and found that LKB1 suppresses growth, in part, by activating the PIKFYVE lipid kinase. Finally, we used chemical inhibitors and a pH-sensitive reporter to determine that LKB1 impairs growth by promoting the internalization of wild-type EGFR in a PIKFYVE-dependent manner.

Enhancer-targeted CRISPR-Activation Rescues Haploinsufficient Autism Susceptibility Genes.

Autism Spectrum Disorder (ASD) is a highly heritable condition with diverse clinical presentations. Approximately 20% of ASD's genetic susceptibility is imparted by de novo mutations of major effect, most of which cause haploinsufficiency. We mapped enhancers of two high confidence autism genes - CHD8 and SCN2A and used CRISPR-based gene activation (CRISPR-A) in hPSC-derived excitatory neurons and cerebral forebrain organoids to correct the effects of haploinsufficiency, taking advantage of the presence of a wildtype allele of each gene and endogenous gene regulation. We found that CRISPR-A induced a sustained increase in CHD8 and SCN2A expression in treated neurons and organoids, with rescue of gene expression levels and mutation-associated phenotypes, including gene expression and physiology. These data support gene activation via targeting enhancers of haploinsufficient genes, as a therapeutic intervention in ASD and other neurodevelopmental disorders.

Non-linear transcriptional responses to gradual modulation of transcription factor dosage.

Genomic loci associated with common traits and diseases are typically non-coding and likely impact gene expression, sometimes coinciding with rare loss-of-function variants in the target gene. However, our understanding of how gradual changes in gene dosage affect molecular, cellular, and organismal traits is currently limited. To address this gap, we induced gradual changes in gene expression of four genes using CRISPR activation and inactivation. Downstream transcriptional consequences of dosage modulation of three master trans-regulators associated with blood cell traits (GFI1B, NFE2, and MYB) were examined using targeted single-cell multimodal sequencing. We showed that guide tiling around the TSS is the most effective way to modulate cis gene expression across a wide range of fold-changes, with further effects from chromatin accessibility and histone marks that differ between the inhibition and activation systems. Our single-cell data allowed us to precisely detect subtle to large gene expression changes in dozens of trans genes, revealing that many responses to dosage changes of these three TFs are non-linear, including non-monotonic behaviours, even when constraining the fold-changes of the master regulators to a copy number gain or loss. We found that the dosage properties are linked to gene constraint and that some of these non-linear responses are enriched for disease and GWAS genes. Overall, our study provides a straightforward and scalable method to precisely modulate gene expression and gain insights into its downstream consequences at high resolution.

Next-generation forward genetic screens: uniting high-throughput perturbations with single-cell analysis.

Programmable genome-engineering technologies, such as CRISPR (clustered regularly interspaced short palindromic repeats) nucleases and massively parallel CRISPR screens that capitalize on this programmability, have transformed biomedical science. These screens connect genes and noncoding genome elements to disease-relevant phenotypes, but until recently have been limited to individual phenotypes such as growth or fluorescent reporters of gene expression. By pairing massively parallel screens with high-dimensional profiling of single-cell types/states, we can now measure how individual genetic perturbations or combinations of perturbations impact the cellular transcriptome, proteome, and epigenome. We review technologies that pair CRISPR screens with single-cell multiomics and the unique opportunities afforded by extending pooled screens using deep multimodal phenotyping.

Prediction of on-target and off-target activity of CRISPR-Cas13d guide RNAs using deep learning.

Transcriptome engineering applications in living cells with RNA-targeting CRISPR effectors depend on accurate prediction of on-target activity and off-target avoidance. Here we design and test ~200,000 RfxCas13d guide RNAs targeting essential genes in human cells with systematically designed mismatches and insertions and deletions (indels). We find that mismatches and indels have a position- and context-dependent impact on Cas13d activity, and mismatches that result in G-U wobble pairings are better tolerated than other single-base mismatches. Using this large-scale dataset, we train a convolutional neural network that we term targeted inhibition of gene expression via gRNA design (TIGER) to predict efficacy from guide sequence and context. TIGER outperforms the existing models at predicting on-target and off-target activity on our dataset and published datasets. We show that TIGER scoring combined with specific mismatches yields the first general framework to modulate transcript expression, enabling the use of RNA-targeting CRISPRs to precisely control gene dosage.

Essential transcription factors for induced neuron differentiation.

Neurogenins are proneural transcription factors required to specify neuronal identity. Their overexpression in human pluripotent stem cells rapidly produces cortical-like neurons with spiking activity and, because of this, they have been widely adopted for human neuron disease models. However, we do not fully understand the key downstream regulatory effectors responsible for driving neural differentiation. Here, using inducible expression of NEUROG1 and NEUROG2, we identify transcription factors (TFs) required for directed neuronal differentiation by combining expression and chromatin accessibility analyses with a pooled in vitro CRISPR-Cas9 screen targeting all ~1900 TFs in the human genome. The loss of one of these essential TFs (ZBTB18) yields few MAP2-positive neurons. Differentiated ZBTB18-null cells have radically altered gene expression, leading to cytoskeletal defects and stunted neurites and spines. In addition to identifying key downstream TFs for neuronal differentiation, our work develops an integrative multi-omics and TFome-wide perturbation platform to rapidly characterize essential TFs for the differentiation of any human cell type.

LKB1 suppresses growth and promotes the internalization of EGFR through the PIKFYVE lipid kinase.

The tumor suppressor LKB1 is a serine/threonine protein kinase that is frequently mutated in human lung adenocarcinoma (LUAD). LKB1 regulates a complex signaling network that is known to control cell polarity and metabolism; however, the pathways that mediate the tumor suppressive activity of LKB1 are incompletely defined. To identify mechanisms of LKB1- mediated growth suppression we developed a spheroid-based cell culture assay to study LKB1- dependent growth. Using this assay, along with genome-wide CRISPR screens and validation with orthogonal methods, we discovered that LKB1 suppresses growth, in part, by activating the PIKFYVE lipid kinase, which promotes the internalization of wild-type EGFR. Our findings reveal a new mechanism of regulation of EGFR, which may have implications for the treatment of LKB1 -mutant LUAD.

Heterozygous deletion of the autism-associated gene CHD8 impairs synaptic function through widespread changes in gene expression and chromatin compaction.

Whole-exome sequencing of autism spectrum disorder (ASD) probands and unaffected family members has identified many genes harboring de novo variants suspected to play a causal role in the disorder. Of these, chromodomain helicase DNA-binding protein 8 (CHD8) is the most recurrently mutated. Despite the prevalence of CHD8 mutations, we have little insight into how CHD8 loss affects genome organization or the functional consequences of these molecular alterations in neurons. Here, we engineered two isogenic human embryonic stem cell lines with CHD8 loss-of-function mutations and characterized differences in differentiated human cortical neurons. We identified hundreds of genes with altered expression, including many involved in neural development and excitatory synaptic transmission. Field recordings and single-cell electrophysiology revealed a 3-fold decrease in firing rates and synaptic activity in CHD8+/- neurons, as well as a similar firing-rate deficit in primary cortical neurons from Chd8+/- mice. These alterations in neuron and synapse function can be reversed by CHD8 overexpression. Moreover, CHD8+/- neurons displayed a large increase in open chromatin across the genome, where the greatest change in compaction was near autism susceptibility candidate 2 (AUTS2), which encodes a transcriptional regulator implicated in ASD. Genes with changes in chromatin accessibility and expression in CHD8+/- neurons have significant overlap with genes mutated in probands for ASD, intellectual disability, and schizophrenia but not with genes mutated in healthy controls or other disease cohorts. Overall, this study characterizes key molecular alterations in genome structure and expression in CHD8+/- neurons and links these changes to impaired neuronal and synaptic function.

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

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

Discovery of target genes and pathways at GWAS loci by pooled single-cell CRISPR screens.

Most variants associated with complex traits and diseases identified by genome-wide association studies (GWAS) map to noncoding regions of the genome with unknown effects. Using ancestrally diverse, biobank-scale GWAS data, massively parallel CRISPR screens, and single-cell transcriptomic and proteomic sequencing, we discovered 124 cis-target genes of 91 noncoding blood trait GWAS loci. Using precise variant insertion through base editing, we connected specific variants with gene expression changes. We also identified trans-effect networks of noncoding loci when cis target genes encoded transcription factors or microRNAs. Networks were themselves enriched for GWAS variants and demonstrated polygenic contributions to complex traits. This platform enables massively parallel characterization of the target genes and mechanisms of human noncoding variants in both cis and trans.

Genome-wide CRISPR/Cas9 screens reveal shared and cell-specific mechanisms of resistance to SHP2 inhibition.

SHP2 (PTPN11) acts upstream of SOS1/2 to enable RAS activation. Allosteric SHP2 inhibitors (SHP2i) in the clinic prevent SHP2 activation, block proliferation of RTK- or cycling RAS mutant-driven cancers, and overcome "adaptive resistance." To identify SHP2i resistance mechanisms, we performed genome-wide CRISPR/Cas9 knockout screens on two SHP2i-sensitive cell lines, recovering genes expected to cause resistance (NF1, PTEN, CDKN1B, LZTR1, and RASA2) and novel targets (INPPL1, MAP4K5, epigenetic modifiers). We screened 14 additional lines with a focused CRISPR library targeting common "hits" from the genome-wide screens. LZTR1 deletion conferred resistance in 12/14 lines, followed by MAP4K5 (8/14), SPRED2/STK40 (6/14), and INPPL1 (5/14). INPPL1, MAP4K5, or LZTR1 deletion reactivated ERK signaling. INPPL1-mediated sensitization to SHP2i required its NPXY motif but not lipid phosphatase activity. MAP4K5 acted upstream of MEK through a kinase-dependent target(s); LZTR1 had cell-dependent effects on RIT and RAS stability. INPPL1, MAP4K5, or LZTR1 deletion also conferred SHP2i resistance in vivo. Defining the SHP2i resistance landscape could suggest effective combination approaches.

Efficient combinatorial targeting of RNA transcripts in single cells with Cas13 RNA Perturb-seq.

Pooled CRISPR screens coupled with single-cell RNA-sequencing have enabled systematic interrogation of gene function and regulatory networks. Here, we introduce Cas13 RNA Perturb-seq (CaRPool-seq), which leverages the RNA-targeting CRISPR-Cas13d system and enables efficient combinatorial perturbations alongside multimodal single-cell profiling. CaRPool-seq encodes multiple perturbations on a cleavable CRISPR array that is associated with a detectable barcode sequence, allowing for the simultaneous targeting of multiple genes. We compared CaRPool-seq to existing Cas9-based methods, highlighting its unique strength to efficiently profile combinatorially perturbed cells. Finally, we apply CaRPool-seq to perform multiplexed combinatorial perturbations of myeloid differentiation regulators in an acute myeloid leukemia (AML) model system and identify extensive interactions between different chromatin regulators that can enhance or suppress AML differentiation phenotypes.

Expanding the Scope of Bacterial CRISPR Activation with PAM-Flexible dCas9 Variants.

CRISPR-Cas transcriptional tools have been widely applied for programmable regulation of complex biological networks. In comparison to eukaryotic systems, bacterial CRISPR activation (CRISPRa) has stringent target site requirements for effective gene activation. While genes may not always have an NGG protospacer adjacent motif (PAM) at the appropriate position, PAM-flexible dCas9 variants can expand the range of targetable sites. Here we systematically evaluate a panel of PAM-flexible dCas9 variants for their ability to activate bacterial genes. We observe that dxCas9-NG provides a high dynamic range of gene activation for sites with NGN PAMs while dSpRY permits modest activity across almost any PAM. Similar trends were observed for heterologous and endogenous promoters. For all variants tested, improved PAM-flexibility comes with the trade-off that CRISPRi-mediated gene repression becomes less effective. Weaker CRISPR interference (CRISPRi) gene repression can be partially rescued by expressing multiple sgRNAs to target many sites in the gene of interest. Our work provides a framework to choose the most effective dCas9 variant for a given set of gene targets, which will further expand the utility of CRISPRa/i gene regulation in bacterial systems.

Cancer genes disfavoring T cell immunity identified via integrated systems approach.

Adoptive T cell therapies (ACT) have been curative for a limited number of cancer patients. The sensitization of cancer cells to T cell killing may expand the benefit of these therapies for more patients. To this end, we use a three-step approach to identify cancer genes that disfavor T cell immunity. First, we profile gene transcripts upregulated by cancer under selection pressure from T cell killing. Second, we identify potential tumor gene targets and pathways that disfavor T cell killing using signaling pathway activation libraries and genome-wide loss-of-function CRISPR-Cas9 screens. Finally, we implement pharmacological perturbation screens to validate these targets and identify BIRC2, ITGAV, DNPEP, BCL2, and ERRα as potential ACT-drug combination candidates. Here, we establish that BIRC2 limits antigen presentation and T cell recognition of tumor cells by suppressing IRF1 activity and provide evidence that BIRC2 inhibition in combination with ACT is an effective strategy to increase efficacy.

Recurrent somatic mutations as predictors of immunotherapy response.

Immune checkpoint blockade (ICB) has transformed the treatment of metastatic cancer but is hindered by variable response rates. A key unmet need is the identification of biomarkers that predict treatment response. To address this, we analyzed six whole exome sequencing cohorts with matched disease outcomes to identify genes and pathways predictive of ICB response. To increase detection power, we focus on genes and pathways that are significantly mutated following correction for epigenetic, replication timing, and sequence-based covariates. Using this technique, we identify several genes (BCLAF1, KRAS, BRAF, and TP53) and pathways (MAPK signaling, p53 associated, and immunomodulatory) as predictors of ICB response and develop the Cancer Immunotherapy Response CLassifiEr (CIRCLE). Compared to tumor mutational burden alone, CIRCLE led to superior prediction of ICB response with a 10.5% increase in sensitivity and a 11% increase in specificity. We envision that CIRCLE and more broadly the analysis of recurrently mutated cancer genes will pave the way for better prognostic tools for cancer immunotherapy.

Invited review: Effect of subacute ruminal acidosis on gut health of dairy cows.

Subacute ruminal acidosis (SARA) is assumed to be a common disease in high-yielding dairy cows. Despite this, the epidemiological evidence is limited by the lack of survey data. The prevalence of SARA has mainly been determined by measuring the pH of ruminal fluid collected using rumenocentesis. This may not be sufficiently accurate, because the symptoms of SARA are not solely due to ruminal pH depression, and ruminal pH varies among sites in the rumen, throughout a 24-h period, and among days. The impact of SARA has mainly been studied by conducting SARA challenges in cows, sheep, and goats based on a combination of feed restriction and high-grain feeding. The methodologies of these challenges vary considerably among studies. Variations include differences in the duration and amount of grain feeding, type of grain, amount and duration of feed restriction, number of experimental cows, and sensitivity of cows to SARA challenges. Grain-based SARA challenges affect gut health. These effects include depressing the pH in, and increasing the toxin content of, digesta. They also include altering the taxonomic composition of microbiota, reducing the functionality of the epithelia throughout the gastrointestinal tract (GIT), and a moderate inflammatory response. The effects on the epithelia include a reduction in its barrier function. Effects on microbiota include reductions in their richness and diversity, which may reduce their functionality and reflect dysbiosis. Changes in the taxonomic composition of gut microbiota throughout the GIT are evident at the phylum level, but less evident and more variable at the genus level. Effects at the phylum level include an increase in the Firmicutes to Bacteroidetes ratio. More studies on the effects of a SARA challenge on the functionality of gut microbiota are needed. The inflammatory response resulting from grain-based SARA challenges is innate and moderate and mainly consists of an acute phase response. This response is likely a combination of systemic inflammation and inflammation of the epithelia of the GIT. The systemic inflammation is assumed to be caused by translocation of immunogenic compounds, including bacterial endotoxins and bioamines, through the epithelia into the interior circulation. This translocation is increased by the increase in concentrations of toxins in digesta and a reduction of the barrier function of epithelia. Severe SARA can cause rumenitis, but moderate SARA may activate an immune response in the epithelia of the GIT. Cows grazing highly fermentable pastures with high sugar contents can also have a low ruminal pH indicative of SARA. This is not accompanied by an inflammatory response but may affect milk production and gut microbiota. Grain-based SARA affects several aspects of gut health, but SARA resulting from grazing high-digestible pastures and insufficient coarse fiber less so. We need to determine which method for inducing SARA is the most representative of on-farm conditions.

Effects of calcareous marine algae on milk production, feed intake, energy balance, mineral status, and inflammatory markers in transition dairy cows.

The objective of this experiment was to compare the effects of calcareous marine algae (CMA; Acid Buf, Celtic Sea Minerals) with a limestone-based control on feed intake, milk production, energy balance, serum mineral metabolites, and inflammatory markers in transition dairy cows. Twenty-two multiparous and 10 primiparous cows were assigned to 2 treatments from 25 d before expected parturition until 42 d postpartum. Cows were assigned to treatment according to a randomized complete block design based on parity, pre-experimental body condition score, previous 305-d milk yield, and either fat + protein yield (for multiparous cows) or predicted transmitting ability for milk yield and fat + protein yield (for primiparous cows). Cows were fed a negative dietary cation-anion difference [-50 mEq/kg] total mixed ration (TMR) based on corn silage, grass silage, and straw during the prepartum period and a 50:50 forage:concentrate TMR based on grass silage, corn silage, and concentrate during the postpartum period. The 2 dietary treatments consisted of a control (CON), which contained limestone as the primary calcium source, and CMA, in which limestone was replaced by CMA at 0.42% and 0.47% of dry matter for the pre- and postpartum periods, respectively. The dietary treatments were fed as 2 different concentrate pellets added to the TMR. Cows fed the CMA diet had higher dry matter intake in both the prepartum (+1.08 kg) and postpartum (+0.94 kg) periods compared with cows fed the CON diet. Fat yield (+0.11 kg), fat concentration (+0.43%), and 4% fat-corrected milk (+1.56 kg) were higher in cows fed CMA than in cows fed CON. The concentration of plasma serum amyloid A was reduced and that of serum P was increased on the CMA treatment compared with the CON treatment. These findings demonstrate the benefits of supplementing CMA to dairy cows during the transition period compared with a CON treatment containing limestone as the primary Ca source.