Efficient site-specific integration of large genes in mammalian cells via continuously evolved recombinases and prime editing.

Methods for the targeted integration of genes in mammalian genomes suffer from low programmability, low efficiencies or low specificities. Here we show that phage-assisted continuous evolution enhances prime-editing-assisted site-specific integrase gene editing (PASSIGE), which couples the programmability of prime editing with the ability of recombinases to precisely integrate large DNA cargoes exceeding 10 kilobases. Evolved and engineered Bxb1 recombinase variants (evoBxb1 and eeBxb1) mediated up to 60% donor integration (3.2-fold that of wild-type Bxb1) in human cell lines with pre-installed recombinase landing sites. In single-transfection experiments at safe-harbour and therapeutically relevant sites, PASSIGE with eeBxb1 led to an average targeted-gene-integration efficiencies of 23% (4.2-fold that of wild-type Bxb1). Notably, integration efficiencies exceeded 30% at multiple sites in primary human fibroblasts. PASSIGE with evoBxb1 or eeBxb1 outperformed PASTE (for 'programmable addition via site-specific targeting elements', a method that uses prime editors fused to recombinases) on average by 9.1-fold and 16-fold, respectively. PASSIGE with continuously evolved recombinases is an unusually efficient method for the targeted integration of genes in mammalian cells.

Millennial drinks: acidity, fluoride content, and enamel softening.

This in vitro study aimed to evaluate the acidity and fluoride content of beverages commonly consumed by millennials and the enamel-softening effect of these drinks on tooth enamel. The study included 13 beverages in 4 categories: energy (sports) drink, flavored sparkling water, kombucha, and other (an unsweetened iced tea, a vegetable-fruit juice blend, and a soft drink). The acidity was measured with a pH/ion meter, and the fluoride concentration was measured with a combined fluoride electrode coupled to the meter (n = 10 measurements per beverage). The Vickers hardness number of extracted molars was measured before and after a 30-minute immersion in 4 representative beverages via 2 immersion protocols (n = 10 per beverage per protocol): (1) immersion in the beverage only and (2) immersion alternating between the beverage and artificial saliva every other minute. The pH and fluoride concentrations of the beverages ranged from 2.652 to 4.242 and from 0.0033 to 0.6045 ppm, respectively. One-way analysis of variance (ANOVA) revealed that all differences between beverages in pH values were statistically significant, as were the majority of differences in fluoride concentrations (P < 0.001). The beverages and the 2 immersion methods significantly affected enamel softening (2-way ANOVA, P = 0.0001 to 0.033). The representative energy drink (pH 2.990; 0.0102 ppm fluoride) caused the greatest enamel softening followed by the representative kombucha (pH 2.820; 0.2036 ppm fluoride). The representative flavored sparkling water (pH 4.066; 0.0098 ppm fluoride) caused significantly less enamel softening than the energy drink and kombucha. A root beer (pH 4.185; 0.6045 ppm fluoride) had the least enamel softening effect. All tested beverages were acidic and had a pH below 4.5; only some contained fluoride. Flavored sparkling water, likely due to its higher pH, caused less enamel softening than the tested energy drink and kombucha. The fluoride content of kombucha and root beer lower their enamel-softening effects. It is imperative that consumers be aware of the erosive potential of beverages they consume.

LKB1 drives stasis and C/EBP-mediated reprogramming to an alveolar type II fate in lung cancer.

LKB1 is among the most frequently altered tumor suppressors in lung adenocarcinoma. Inactivation of Lkb1 accelerates the growth and progression of oncogenic KRAS-driven lung tumors in mouse models. However, the molecular mechanisms by which LKB1 constrains lung tumorigenesis and whether the cancer state that stems from Lkb1 deficiency can be reverted remains unknown. To identify the processes governed by LKB1 in vivo, we generated an allele which enables Lkb1 inactivation at tumor initiation and subsequent Lkb1 restoration in established tumors. Restoration of Lkb1 in oncogenic KRAS-driven lung tumors suppressed proliferation and led to tumor stasis. Lkb1 restoration activated targets of C/EBP transcription factors and drove neoplastic cells from a progenitor-like state to a less proliferative alveolar type II cell-like state. We show that C/EBP transcription factors govern a subset of genes that are induced by LKB1 and depend upon NKX2-1. We also demonstrate that a defining factor of the alveolar type II lineage, C/EBPα, constrains oncogenic KRAS-driven lung tumor growth in vivo. Thus, this key tumor suppressor regulates lineage-specific transcription factors, thereby constraining lung tumor development through enforced differentiation.

LKB1 inactivation modulates chromatin accessibility to drive metastatic progression.

Metastasis is the leading cause of cancer-related deaths and enables cancer cells to compromise organ function by expanding in secondary sites. Since primary tumours and metastases often share the same constellation of driver mutations, the mechanisms that drive their distinct phenotypes are unclear. Here we show that inactivation of the frequently mutated tumour suppressor gene LKB1 (encoding liver kinase B1) has evolving effects throughout the progression of lung cancer, which leads to the differential epigenetic re-programming of early-stage primary tumours compared with late-stage metastases. By integrating genome-scale CRISPR-Cas9 screening with bulk and single-cell multi-omic analyses, we unexpectedly identify LKB1 as a master regulator of chromatin accessibility in lung adenocarcinoma primary tumours. Using an in vivo model of metastatic progression, we further show that loss of LKB1 activates the early endoderm transcription factor SOX17 in metastases and a metastatic-like sub-population of cancer cells within primary tumours. The expression of SOX17 is necessary and sufficient to drive a second wave of epigenetic changes in LKB1-deficient cells that enhances metastatic ability. Overall, our study demonstrates how the downstream effects of an individual driver mutation can change throughout cancer development, with implications for stage-specific therapeutic resistance mechanisms and the gene regulatory underpinnings of metastatic evolution.

High-throughput single-cell chromatin accessibility CRISPR screens enable unbiased identification of regulatory networks in cancer.

Chromatin accessibility profiling can identify putative regulatory regions genome wide; however, pooled single-cell methods for assessing the effects of regulatory perturbations on accessibility are limited. Here, we report a modified droplet-based single-cell ATAC-seq protocol for perturbing and evaluating dynamic single-cell epigenetic states. This method (Spear-ATAC) enables simultaneous read-out of chromatin accessibility profiles and integrated sgRNA spacer sequences from thousands of individual cells at once. Spear-ATAC profiling of 104,592 cells representing 414 sgRNA knock-down populations reveals the temporal dynamics of epigenetic responses to regulatory perturbations in cancer cells and the associations between transcription factor binding profiles.

ArchR is a scalable software package for integrative single-cell chromatin accessibility analysis.

The advent of single-cell chromatin accessibility profiling has accelerated the ability to map gene regulatory landscapes but has outpaced the development of scalable software to rapidly extract biological meaning from these data. Here we present a software suite for single-cell analysis of regulatory chromatin in R (ArchR; https://www.archrproject.com/ ) that enables fast and comprehensive analysis of single-cell chromatin accessibility data. ArchR provides an intuitive, user-focused interface for complex single-cell analyses, including doublet removal, single-cell clustering and cell type identification, unified peak set generation, cellular trajectory identification, DNA element-to-gene linkage, transcription factor footprinting, mRNA expression level prediction from chromatin accessibility and multi-omic integration with single-cell RNA sequencing (scRNA-seq). Enabling the analysis of over 1.2 million single cells within 8 h on a standard Unix laptop, ArchR is a comprehensive software suite for end-to-end analysis of single-cell chromatin accessibility that will accelerate the understanding of gene regulation at the resolution of individual cells.

Finding needles in a haystack: dissecting tumor heterogeneity with single-cell transcriptomic and chromatin accessibility profiling.

Tumor evolution often results in a wealth of heterogeneous cancer cell types within a single tumor - heterogeneity that can include epigenetic and gene expression changes that are impossible to identify from histological features alone. The invasion of cancer cells into nearby healthy tissue, accompanied by the infiltration of responding immune cells, results in an even more complex architecture of tumor and non-tumor cells. However, bulk genomics-based methods can only assay the aggregate transcriptomic and epigenetic profiles across all of this rich cellular diversity. Such bulk averaging hides small subpopulations of tumor cells with unique phenotypes that might result in therapeutic resistance or metastatic progression. The advent of single-cell-based genomics assays for measuring transcription and chromatin accessibility - particularly scRNA-seq and scATAC-seq - has enabled the dissection of cell-types within tumors at a scale and resolution capable of unraveling the epigenetic and gene expression programs of rare and unique cellular subpopulations. This Review focuses on recent advances in scRNA-seq and scATAC-seq technologies and their application to cancer biology in the context of furthering our understanding of tumor heterogeneity.

Associations between sensory processing and electrophysiological and neurochemical measures in children with ASD: an EEG-MRS study.

BACKGROUND: Autism spectrum disorder (ASD) is associated with hyper- and/or hypo-sensitivity to sensory input. Spontaneous alpha power, which plays an important role in shaping responsivity to sensory information, is reduced across the lifespan in individuals with ASD. Furthermore, an excitatory/inhibitory imbalance has also been linked to sensory dysfunction in ASD and has been hypothesized to underlie atypical patterns of spontaneous brain activity. The present study examined whether resting-state alpha power differed in children with ASD as compared to TD children, and investigated the relationships between alpha levels, concentrations of excitatory and inhibitory neurotransmitters, and atypical sensory processing in ASD. METHODS: Participants included thirty-one children and adolescents with ASD and thirty-one age- and IQ-matched typically developing (TD) participants. Resting-state electroencephalography (EEG) was used to obtain measures of alpha power. A subset of participants (ASD = 16; TD = 16) also completed a magnetic resonance spectroscopy (MRS) protocol in order to measure concentrations of excitatory (glutamate + glutamine; Glx) and inhibitory (GABA) neurotransmitters. RESULTS: Children with ASD evidenced significantly decreased resting alpha power compared to their TD peers. MRS estimates of GABA and Glx did not differ between groups with the exception of Glx in the temporal-parietal junction. Inter-individual differences in alpha power within the ASD group were not associated with region-specific concentrations of GABA or Glx, nor were they associated with sensory processing differences. However, atypically decreased Glx was associated with increased sensory impairment in children with ASD. CONCLUSIONS: Although we replicated prior reports of decreased alpha power in ASD, atypically reduced alpha was not related to neurochemical differences or sensory symptoms in ASD. Instead, reduced Glx in the temporal-parietal cortex was associated with greater hyper-sensitivity in ASD. Together, these findings may provide insight into the neural underpinnings of sensory processing differences present in ASD.

Zmat3 Is a Key Splicing Regulator in the p53 Tumor Suppression Program.

Although TP53 is the most commonly mutated gene in human cancers, the p53-dependent transcriptional programs mediating tumor suppression remain incompletely understood. Here, to uncover critical components downstream of p53 in tumor suppression, we perform unbiased RNAi and CRISPR-Cas9-based genetic screens in vivo. These screens converge upon the p53-inducible gene Zmat3, encoding an RNA-binding protein, and we demonstrate that ZMAT3 is an important tumor suppressor downstream of p53 in mouse KrasG12D-driven lung and liver cancers and human carcinomas. Integrative analysis of the ZMAT3 RNA-binding landscape and transcriptomic profiling reveals that ZMAT3 directly modulates exon inclusion in transcripts encoding proteins of diverse functions, including the p53 inhibitors MDM4 and MDM2, splicing regulators, and components of varied cellular processes. Interestingly, these exons are enriched in NMD signals, and, accordingly, ZMAT3 broadly affects target transcript stability. Collectively, these studies reveal ZMAT3 as a novel RNA-splicing and homeostasis regulator and a key component of p53-mediated tumor suppression.

Helicase-inactivating BRIP1 mutation yields Fanconi anemia with microcephaly and other congenital abnormalities.

Fanconi anemia is a genetically and phenotypically heterogeneous disorder characterized by congenital anomalies, bone marrow failure, cancer, and sensitivity of chromosomes to DNA cross-linking agents. One of the 22 genes responsible for Fanconi anemia is BRIP1, in which biallelic truncating mutations lead to Fanconi anemia group J and monoallelic truncating mutations predispose to certain cancers. However, of the more than 1000 reported missense mutations in BRIP1, very few have been functionally characterized. We evaluated the functional consequence of BRIP1 p.R848H (c.2543G > A), which was homozygous in two cousins with low birth weight, microcephaly, upper limb abnormalities, and imperforate anus and for whom chromosome breakage analysis of patient cells revealed increased mitomycin C sensitivity. BRIP1 p.R848H alters a highly conserved residue in the catalytic DNA helicase domain. We show that BRIP1 p.R848H leads to a defect in helicase activity. Heterozygosity at this missense has been reported in multiple cancer patients but, in the absence of functional studies, classified as of unknown significance. Our results support that this mutation is pathogenic for Fanconi anemia in homozygotes and for increased cancer susceptibility in heterozygous carriers.

Refuges and ecological traps: Extreme drought threatens persistence of an endangered fish in intermittent streams.

Recent droughts raise global concern over potential biodiversity loss and mitigating impacts to vulnerable species has become a management priority. However, drought impacts on populations are difficult to predict, in part, because habitat refuges can buffer organisms from harsh environmental conditions. In a global change context, more extreme droughts may turn previously suitable habitats into ecological traps, where vulnerable species can no longer persist. Here, we explore the impacts of California's recent record-breaking drought on endangered juvenile Coho salmon. We estimated the variability of cumulative salmon survival using mark-recapture of nearly 20,000 tagged fish in intermittent stream pools during a 7-year period encompassing drought and non-drought conditions. We then determined the relative importance of physical habitat, streamflow, precipitation, landscape, and biological characteristics that may limit survival during drought. Our most striking result was an increase in the number of pools with reduced or zero survival during drought years and a coincident increase in spatial variability in survival among study reaches. In nearly half of the stream pools, salmon survival during drought was similar to mean survival of pools assessed during non-drought years, indicating some pools had remarkable resistance (ability to withstand disturbance) to extreme drought. Lower survival was most attributable to longer duration of disconnection between upstream and downstream habitats, a consequence of increasing drought severity. Our results not only suggest that many pools sustain juvenile salmon in non-drought years transition into ecological traps during drought but also highlight that some pools serve as refuges even under extreme drought conditions. Projected increases in drought severity that lead to longer droughts and greater habitat fragmentation could transform an increasing proportion of suitable habitats into ecological traps. Predicting future impacts of drought on Coho salmon and other sensitive species will require identification and protection of drought refuges and management strategies that prevent further habitat fragmentation.

CRISPR screens in cancer spheroids identify 3D growth-specific vulnerabilities.

Cancer genomics studies have identified thousands of putative cancer driver genes1. Development of high-throughput and accurate models to define the functions of these genes is a major challenge. Here we devised a scalable cancer-spheroid model and performed genome-wide CRISPR screens in 2D monolayers and 3D lung-cancer spheroids. CRISPR phenotypes in 3D more accurately recapitulated those of in vivo tumours, and genes with differential sensitivities between 2D and 3D conditions were highly enriched for genes that are mutated in lung cancers. These analyses also revealed drivers that are essential for cancer growth in 3D and in vivo, but not in 2D. Notably, we found that carboxypeptidase D is responsible for removal of a C-terminal RKRR motif2 from the α-chain of the insulin-like growth factor 1 receptor that is critical for receptor activity. Carboxypeptidase D expression correlates with patient outcomes in patients with lung cancer, and loss of carboxypeptidase D reduced tumour growth. Our results reveal key differences between 2D and 3D cancer models, and establish a generalizable strategy for performing CRISPR screens in spheroids to reveal cancer vulnerabilities.

Loss-of-function mutations in UDP-Glucose 6-Dehydrogenase cause recessive developmental epileptic encephalopathy.

Developmental epileptic encephalopathies are devastating disorders characterized by intractable epileptic seizures and developmental delay. Here, we report an allelic series of germline recessive mutations in UGDH in 36 cases from 25 families presenting with epileptic encephalopathy with developmental delay and hypotonia. UGDH encodes an oxidoreductase that converts UDP-glucose to UDP-glucuronic acid, a key component of specific proteoglycans and glycolipids. Consistent with being loss-of-function alleles, we show using patients' primary fibroblasts and biochemical assays, that these mutations either impair UGDH stability, oligomerization, or enzymatic activity. In vitro, patient-derived cerebral organoids are smaller with a reduced number of proliferating neuronal progenitors while mutant ugdh zebrafish do not phenocopy the human disease. Our study defines UGDH as a key player for the production of extracellular matrix components that are essential for human brain development. Based on the incidence of variants observed, UGDH mutations are likely to be a frequent cause of recessive epileptic encephalopathy.

Massively parallel single-cell chromatin landscapes of human immune cell development and intratumoral T cell exhaustion.

Understanding complex tissues requires single-cell deconstruction of gene regulation with precision and scale. Here, we assess the performance of a massively parallel droplet-based method for mapping transposase-accessible chromatin in single cells using sequencing (scATAC-seq). We apply scATAC-seq to obtain chromatin profiles of more than 200,000 single cells in human blood and basal cell carcinoma. In blood, application of scATAC-seq enables marker-free identification of cell type-specific cis- and trans-regulatory elements, mapping of disease-associated enhancer activity and reconstruction of trajectories of cellular differentiation. In basal cell carcinoma, application of scATAC-seq reveals regulatory networks in malignant, stromal and immune cells in the tumor microenvironment. Analysis of scATAC-seq profiles from serial tumor biopsies before and after programmed cell death protein 1 blockade identifies chromatin regulators of therapy-responsive T cell subsets and reveals a shared regulatory program that governs intratumoral CD8+ T cell exhaustion and CD4+ T follicular helper cell development. We anticipate that scATAC-seq will enable the unbiased discovery of gene regulatory factors across diverse biological systems.

Rhizobacterial Community Assembly Patterns Vary Between Crop Species.

Currently our limited understanding of crop rhizosphere community assembly hinders attempts to manipulate it beneficially. Variation in root communities has been attributed to plant host effects, soil type, and plant condition, but it is hard to disentangle the relative importance of soil and host without experimental manipulation. To examine the effects of soil origin and host plant on root associated bacterial communities we experimentally manipulated four crop species in split-plot mesocosms and surveyed variation in bacterial diversity by Illumina amplicon sequencing. Overall, plant species had a greater impact than soil type on community composition. While plant species associated with different Operational Taxonomic Units (OTUs) in different soils, plants tended to recruit bacteria from similar, higher order, taxonomic groups in different soils. However, the effect of soil on root-associated communities varied between crop species: Onion had a relatively invariant bacterial community while other species (maize and pea) had a more variable community structure. Dynamic communities could result from environment specific recruitment, differential bacterial colonization or reflect broader symbiont host range; while invariant community assembly implies tighter evolutionary or ecological interactions between plants and root-associated bacteria. Irrespective of mechanism, it appears both communities and community assembly rules vary between crop species.

Inherited thrombocytopenia associated with mutation of UDP-galactose-4-epimerase (GALE).

Severe thrombocytopenia, characterized by dysplastic megakaryocytes and intracranial bleeding, was diagnosed in six individuals from a consanguineous kindred. Three of the individuals were successfully treated by bone marrow transplant. Whole-exome sequencing and homozygosity mapping of multiple family members, coupled with whole-genome sequencing to reveal shared non-coding variants, revealed one potentially functional variant segregating with thrombocytopenia under a recessive model: GALE p.R51W (c.C151T, NM_001127621). The mutation is extremely rare (allele frequency = 2.5 × 10-05), and the likelihood of the observed co-segregation occurring by chance is 1.2 × 10-06. GALE encodes UDP-galactose-4-epimerase, an enzyme of galactose metabolism and glycosylation responsible for two reversible reactions: interconversion of UDP-galactose with UDP-glucose and interconversion of UDP-N-acetylgalactosamine with UDP-N-acetylglucosamine. The mutation alters an amino acid residue that is conserved from yeast to humans. The variant protein has both significantly lower enzymatic activity for both interconversion reactions and highly significant thermal instability. Proper glycosylation is critical to normal hematopoiesis, in particular to megakaryocyte and platelet development, as reflected in the presence of thrombocytopenia in the context of congenital disorders of glycosylation. Mutations in GALE have not previously been associated with thrombocytopenia. Our results suggest that GALE p.R51W is inadequate for normal glycosylation and thereby may impair megakaryocyte and platelet development. If other mutations in GALE are shown to have similar consequences, this gene may be proven to play a critical role in hematopoiesis.

De novo mutation in RING1 with epigenetic effects on neurodevelopment.

RING1 is an E3-ubiquitin ligase that is involved in epigenetic control of transcription during development. It is a component of the polycomb repressive complex 1, and its role in that complex is to ubiquitylate histone H2A. In a 13-year-old girl with syndromic neurodevelopmental disabilities, we identified a de novo mutation, RING1 p.R95Q, which alters a conserved arginine residue in the catalytic RING domain. In vitro assays demonstrated that the mutant RING1 retains capacity to catalyze ubiquitin chain formation, but is defective in its ability to ubiquitylate histone H2A in nucleosomes. Consistent with this in vitro effect, cells of the patient showed decreased monoubiquitylation of histone H2A. We modeled the mutant RING1 in Caenorhabditis elegans by editing the comparable amino acid change into spat-3, the suggested RING1 ortholog. Animals with either the missense mutation or complete knockout of spat-3 were defective in monoubiquitylation of histone H2A and had defects in neuronal migration and axon guidance. Relevant to our patient, animals heterozygous for either the missense or knockout allele also showed neuronal defects. Our results support three conclusions: mutation of RING1 is the likely cause of a human neurodevelopmental syndrome, mutation of RING1 can disrupt histone H2A ubiquitylation without disrupting RING1 catalytic activity, and the comparable mutation in C. elegans spat-3 both recapitulates the effects on histone H2A ubiquitylation and leads to neurodevelopmental abnormalities. This role for RING1 adds to our understanding of the importance of aberrant epigenetic effects as causes of human neurodevelopmental disorders.

Integrin-ß4 identifies cancer stem cell-enriched populations of partially mesenchymal carcinoma cells.

Neoplastic cells within individual carcinomas often exhibit considerable phenotypic heterogeneity in their epithelial versus mesenchymal-like cell states. Because carcinoma cells with mesenchymal features are often more resistant to therapy and may serve as a source of relapse, we sought to determine whether such cells could be further stratified into functionally distinct subtypes. Indeed, we find that a basal epithelial marker, integrin-ß4 (ITGB4), can be used to enable stratification of mesenchymal-like triple-negative breast cancer (TNBC) cells that differ from one another in their relative tumorigenic abilities. Notably, we demonstrate that ITGB4+ cancer stem cell (CSC)-enriched mesenchymal cells reside in an intermediate epithelial/mesenchymal phenotypic state. Among patients with TNBC who received chemotherapy, elevated ITGB4 expression was associated with a worse 5-year probability of relapse-free survival. Mechanistically, we find that the ZEB1 (zinc finger E-box binding homeobox 1) transcription factor activity in highly mesenchymal SUM159 TNBC cells can repress expression of the epithelial transcription factor TAp63α (tumor protein 63 isoform 1), a protein that promotes ITGB4 expression. In addition, we demonstrate that ZEB1 and ITGB4 are important in modulating the histopathological phenotypes of tumors derived from mesenchymal TNBC cells. Hence, mesenchymal carcinoma cell populations are internally heterogeneous, and ITGB4 is a mechanistically driven prognostic biomarker that can be used to identify the more aggressive subtypes of mesenchymal carcinoma cells in TNBC. The ability to rapidly isolate and mechanistically interrogate the CSC-enriched, partially mesenchymal carcinoma cells should further enable identification of novel therapeutic opportunities to improve the prognosis for high-risk patients with TNBC.

Supplemental Learning in the Laboratory: An Innovative Approach for Evaluating Knowledge and Method Transfer.

The Multi-Rule Quality Control System (MRQCS) is a tool currently employed by the Centers for Disease Control and Prevention (CDC) to evaluate and compare laboratory performance. We have applied the MRQCS to a comparison of instructor and computer-led pre-laboratory lectures for a supplemental learning experiment. Students in general chemistry and analytical chemistry from both two- and four-year institutions performed two laboratory experiments as part of their normal laboratory curriculum. The first laboratory experiment was a foundational learning experiment in which all the students were introduced to Beer-Lambert's Law and spectrophotometric light absorbance measurements. The foundational learning experiment was instructor-led only, and participant performance was evaluated against a mean characterized value. The second laboratory experiment was a supplemental learning experiment in which students were asked to build upon the methodology they learned in the foundational learning experiment and apply it to a different analyte. The instruction type was varied randomly into two delivery modes, participants receiving either instructor-led or computer-led pre-laboratory instruction. The MRQCS was applied and determined that no statistical difference was found to exist in the QC (quality control) passing rates between the participants in the instructor-led instruction and the participants in the computer-led instruction. These findings demonstrate the successful application of the MRQCS to evaluate knowledge and technology transfer.