Disease-associated astrocyte epigenetic memory promotes CNS pathology

Disease-associated astrocyte subsets contribute to the pathology of neurologic diseases, including multiple sclerosis and experimental autoimmune encephalomyelitis1-8 (EAE), an experimental model for multiple sclerosis. However, little is known about the stability of these astrocyte subsets and their ability to integrate past stimulation events. Here we report the identification of an epigenetically controlled memory astrocyte subset that exhibits exacerbated pro-inflammatory responses upon rechallenge. Specifically, using a combination of single-cell RNA sequencing, assay for transposase-accessible chromatin with sequencing, chromatin immunoprecipitation with sequencing, focused interrogation of cells by nucleic acid detection and sequencing, and cell-specific in vivo CRISPR-Cas9-based genetic perturbation studies we established that astrocyte memory is controlled by the metabolic enzyme ATP-citrate lyase (ACLY), which produces acetyl coenzyme A (acetyl-CoA) that is used by histone acetyltransferase p300 to control chromatin accessibility. The number of ACLY+p300+ memory astrocytes is increased in acute and chronic EAE models, and their genetic inactivation ameliorated EAE. We also detected the pro-inflammatory memory phenotype in human astrocytes in vitro; single-cell RNA sequencing and immunohistochemistry studies detected increased numbers of ACLY+p300+ astrocytes in chronic multiple sclerosis lesions. In summary, these studies define an epigenetically controlled memory astrocyte subset that promotes CNS pathology in EAE and, potentially, multiple sclerosis. These findings may guide novel therapeutic approaches for multiple sclerosis and other neurologic diseases.

Disease-associated astrocyte epigenetic memory promotes CNS pathology

Astrocytes play important roles in the central nervous system (CNS) physiology and pathology. Indeed, astrocyte subsets defined by specific transcriptional activation states contribute to the pathology of neurologic diseases, including multiple sclerosis (MS) and its pre-clinical model experimental autoimmune encephalomyelitis (EAE) 1-8 . However, little is known about the stability of these disease-associated astrocyte subsets, their regulation, and whether they integrate past stimulation events to respond to subsequent challenges. Here, we describe the identification of an epigenetically controlled memory astrocyte subset which exhibits exacerbated pro-inflammatory responses upon re-challenge. Specifically, using a combination of single-cell RNA sequencing (scRNA-seq), assay for transposase-accessible chromatin with sequencing (ATAC-seq), chromatin immunoprecipitation with sequencing (ChIP-seq), focused interrogation of cells by nucleic acid detection and sequencing (FIND-seq), and cell-specific in vivo CRISPR/Cas9-based genetic perturbation studies we established that astrocyte memory is controlled by the metabolic enzyme ATP citrate lyase (ACLY), which produces acetyl coenzyme A (acetyl-CoA) used by the histone acetyltransferase p300 to control chromatin accessibility. ACLY + p300 + memory astrocytes are increased in acute and chronic EAE models; the genetic targeting of ACLY + p300 + astrocytes using CRISPR/Cas9 ameliorated EAE. We also detected responses consistent with a pro-inflammatory memory phenotype in human astrocytes in vitro ; scRNA-seq and immunohistochemistry studies detected increased ACLY + p300 + astrocytes in chronic MS lesions. In summary, these studies define an epigenetically controlled memory astrocyte subset that promotes CNS pathology in EAE and, potentially, MS. These findings may guide novel therapeutic approaches for MS and other neurologic diseases.

Microfluidics-free single-cell genomics with templated emulsification

Current single-cell RNA-sequencing approaches have limitations that stem from the microfluidic devices or fluid handling steps required for sample processing. We develop a method that does not require specialized microfluidic devices, expertise or hardware. Our approach is based on particle-templated emulsification, which allows single-cell encapsulation and barcoding of cDNA in uniform droplet emulsions with only a vortexer. Particle-templated instant partition sequencing (PIP-seq) accommodates a wide range of emulsification formats, including microwell plates and large-volume conical tubes, enabling thousands of samples or millions of cells to be processed in minutes. We demonstrate that PIP-seq produces high-purity transcriptomes in mouse-human mixing studies, is compatible with multiomics measurements and can accurately characterize cell types in human breast tissue compared to a commercial microfluidic platform. Single-cell transcriptional profiling of mixed phenotype acute leukemia using PIP-seq reveals the emergence of heterogeneity within chemotherapy-resistant cell subsets that were hidden by standard immunophenotyping. PIP-seq is a simple, flexible and scalable next-generation workflow that extends single-cell sequencing to new applications.

Multiomic Single Cell Sequencing Identifies Stemlike Nature of Mixed Phenotype Acute Leukemia and Provides Novel Risk Stratification

Mixed phenotype acute leukemia (MPAL) is a leukemia whose biologic drivers are poorly understood, therapeutic strategy remains unclear, and prognosis is poor. We performed multiomic single cell (SC) profiling of 14 newly diagnosed adult MPAL patients to characterize the immunophenotypic, genetic, and transcriptional landscapes of MPAL. We show that neither genetic profile nor transcriptome reliably correlate with specific MPAL immunophenotypes. However, progressive acquisition of mutations is associated with increased expression of immunophenotypic markers of immaturity. Using SC transcriptional profiling, we find that MPAL blasts express a stem cell-like transcriptional profile distinct from other acute leukemias and indicative of high differentiation potential. Further, patients with the highest differentiation potential demonstrated inferior survival in our dataset. A gene set score, MPAL95, derived from genes highly enriched in this cohort, is applicable to bulk RNA sequencing data and was predictive of survival in an independent patient cohort, suggesting utility for clinical risk stratification.

Droplet-based forward genetic screening of astrocyte-microglia cross-talk

Cell-cell interactions in the central nervous system play important roles in neurologic diseases. However, little is known about the specific molecular pathways involved, and methods for their systematic identification are limited. Here, we developed a forward genetic screening platform that combines CRISPR-Cas9 perturbations, cell coculture in picoliter droplets, and microfluidic-based fluorescence-activated droplet sorting to identify mechanisms of cell-cell communication. We used SPEAC-seq (systematic perturbation of encapsulated associated cells followed by sequencing), in combination with in vivo genetic perturbations, to identify microglia-produced amphiregulin as a suppressor of disease-promoting astrocyte responses in multiple sclerosis preclinical models and clinical samples. Thus, SPEAC-seq enables the high-throughput systematic identification of cell-cell communication mechanisms.

HIV silencing and cell survival signatures in infected T cell reservoirs

Rare CD4 T cells that contain HIV under antiretroviral therapy represent an important barrier to HIV cure1-3, but the infeasibility of isolating and characterizing these cells in their natural state has led to uncertainty about whether they possess distinctive attributes that HIV cure-directed therapies might exploit. Here we address this challenge using a microfluidic technology that isolates the transcriptomes of HIV-infected cells based solely on the detection of HIV DNA. HIV-DNA+ memory CD4 T cells in the blood from people receiving antiretroviral therapy showed inhibition of six transcriptomic pathways, including death receptor signalling, necroptosis signalling and antiproliferative Gα12/13 signalling. Moreover, two groups of genes identified by network co-expression analysis were significantly associated with HIV-DNA+ cells. These genes (n = 145) accounted for just 0.81% of the measured transcriptome and included negative regulators of HIV transcription that were higher in HIV-DNA+ cells, positive regulators of HIV transcription that were lower in HIV-DNA+ cells, and other genes involved in RNA processing, negative regulation of mRNA translation, and regulation of cell state and fate. These findings reveal that HIV-infected memory CD4 T cells under antiretroviral therapy are a distinctive population with host gene expression patterns that favour HIV silencing, cell survival and cell proliferation, with important implications for the development of HIV cure strategies.

Identification of astrocyte regulators by nucleic acid cytometry

Multiple sclerosis is a chronic inflammatory disease of the central nervous system1. Astrocytes are heterogeneous glial cells that are resident in the central nervous system and participate in the pathogenesis of multiple sclerosis and its model experimental autoimmune encephalomyelitis2,3. However, few unique surface markers are available for the isolation of astrocyte subsets, preventing their analysis and the identification of candidate therapeutic targets; these limitations are further amplified by the rarity of pathogenic astrocytes. Here, to address these challenges, we developed focused interrogation of cells by nucleic acid detection and sequencing (FIND-seq), a high-throughput microfluidic cytometry method that combines encapsulation of cells in droplets, PCR-based detection of target nucleic acids and droplet sorting to enable in-depth transcriptomic analyses of cells of interest at single-cell resolution. We applied FIND-seq to study the regulation of astrocytes characterized by the splicing-driven activation of the transcription factor XBP1, which promotes disease pathology in multiple sclerosis and experimental autoimmune encephalomyelitis4. Using FIND-seq in combination with conditional-knockout mice, in vivo CRISPR-Cas9-driven genetic perturbation studies and bulk and single-cell RNA sequencing analyses of samples from mouse experimental autoimmune encephalomyelitis and humans with multiple sclerosis, we identified a new role for the nuclear receptor NR3C2 and its corepressor NCOR2 in limiting XBP1-driven pathogenic astrocyte responses. In summary, we used FIND-seq to identify a therapeutically targetable mechanism that limits XBP1-driven pathogenic astrocyte responses. FIND-seq enables the investigation of previously inaccessible cells, including rare cell subsets defined by unique gene expression signatures or other nucleic acid markers.

Barcoded viral tracing of single-cell interactions in central nervous system inflammation

Cell-cell interactions control the physiology and pathology of the central nervous system (CNS). To study astrocyte cell interactions in vivo, we developed rabies barcode interaction detection followed by sequencing (RABID-seq), which combines barcoded viral tracing and single-cell RNA sequencing (scRNA-seq). Using RABID-seq, we identified axon guidance molecules as candidate mediators of microglia-astrocyte interactions that promote CNS pathology in experimental autoimmune encephalomyelitis (EAE) and, potentially, multiple sclerosis (MS). In vivo cell-specific genetic perturbation EAE studies, in vitro systems, and the analysis of MS scRNA-seq datasets and CNS tissue established that Sema4D and Ephrin-B3 expressed in microglia control astrocyte responses via PlexinB2 and EphB3, respectively. Furthermore, a CNS-penetrant EphB3 inhibitor suppressed astrocyte and microglia proinflammatory responses and ameliorated EAE. In summary, RABID-seq identified microglia-astrocyte interactions and candidate therapeutic targets.

Targeted Single-Cell RNA and DNA Sequencing With Fluorescence-Activated Droplet Merger

Analyzing every cell in a diverse sample provides insight into population-level heterogeneity, but abundant cell types dominate the analysis and rarer populations are scarcely represented in the data. To focus on specific cell types, the current paradigm is to physically isolate subsets of interest prior to analysis; however, it remains difficult to isolate and then single-cell sequence such populations because of compounding losses. Here, we describe an alternative approach that selectively merges cells with reagents to achieve enzymatic reactions without having to physically isolate cells. We apply this technique to perform single-cell transcriptome and genome sequencing of specific cell subsets. Our method for analyzing heterogeneous populations obviates the need for pre- or post-enrichment and simplifies single-cell workflows, making it useful for other applications in single-cell biology, combinatorial chemical synthesis, and drug screening.

MAFG-driven astrocytes promote CNS inflammation

Multiple sclerosis is a chronic inflammatory disease of the CNS1. Astrocytes contribute to the pathogenesis of multiple sclerosis2, but little is known about the heterogeneity of astrocytes and its regulation. Here we report the analysis of astrocytes in multiple sclerosis and its preclinical model experimental autoimmune encephalomyelitis (EAE) by single-cell RNA sequencing in combination with cell-specific Ribotag RNA profiling, assay for transposase-accessible chromatin with sequencing (ATAC-seq), chromatin immunoprecipitation with sequencing (ChIP-seq), genome-wide analysis of DNA methylation and in vivo CRISPR-Cas9-based genetic perturbations. We identified astrocytes in EAE and multiple sclerosis that were characterized by decreased expression of NRF2 and increased expression of MAFG, which cooperates with MAT2α to promote DNA methylation and represses antioxidant and anti-inflammatory transcriptional programs. Granulocyte-macrophage colony-stimulating factor (GM-CSF) signalling in astrocytes drives the expression of MAFG and MAT2α and pro-inflammatory transcriptional modules, contributing to CNS pathology in EAE and, potentially, multiple sclerosis. Our results identify candidate therapeutic targets in multiple sclerosis.

Microfluidic bead encapsulation above 20 kHz with triggered drop formation

Microsphere beads are functionalized with oligonucleotides, antibodies, and other moieties to enable specific detection of analytes. Droplet microfluidics leverages this for single-molecule or -cell analysis by pairing beads and targets in water-in-oil droplets. Pairing is achieved with devices operating in the dripping regime, limiting throughput. Here, we describe a pairing method that uses beads to trigger the breakup of a jet into monodispersed droplets. We use the method to pair 105 Human T cells with polyacrylamide beads ten times faster than methods operating in the dripping regime. Our method improves the throughput of bead-based droplet workflows, enabling analysis of large populations and the detection of rare events.

Concentric electrodes improve microfluidic droplet sorting

Microfluidic droplet sorting allows selection of subpopulations of cells, nucleic acids, and biomolecules with soluble assays. Dielectrophoresis is widely used for sorting because it generates strong forces on droplets, actuates rapidly, and is easy to integrate into microfluidic chips. However, existing device designs apply a short force, limiting the deflection of droplets, and therefore the speed and reliability of sorting. We describe a concentric design that applies a long force, allowing large deflections and increased reliability. We demonstrate the utility of this design by sorting polydisperse emulsions, which are typically difficult to sort with high purity.

Single-Cell RT-PCR in Microfluidic Droplets with Integrated Chemical Lysis

Droplet microfluidics can identify and sort cells using digital reverse transcription polymerase chain reaction (RT-PCR) signals from individual cells. However, current methods require multiple microfabricated devices for enzymatic cell lysis and PCR reagent addition, making the process complex and prone to failure. Here, we describe a new approach that integrates all components into a single device. The method enables controlled exposure of isolated single cells to a high pH buffer, which lyses cells and inactivates reaction inhibitors but can be instantly neutralized with RT-PCR buffer. Using our chemical lysis approach, we distinguish individual cells' gene expression with data quality equivalent to more complex two-step workflows. Our system accepts cells and produces droplets ready for amplification, making single-cell droplet RT-PCR faster and more reliable.

Finding a helix in a haystack: nucleic acid cytometry with droplet microfluidics

Nucleic acids encode the information of life, programming cellular functions and dictating many biological outcomes. Differentiating between cells based on their nucleic acid programs is, thus, a powerful way to unravel the genetic bases of many phenotypes. This is especially important considering that most cells exist in heterogeneous populations, requiring them to be isolated before they can be studied. Existing flow cytometry techniques, however, are unable to reliably recover specific cells based on nucleic acid content. Nucleic acid cytometry is a new field built on droplet microfluidics that allows robust identification, sorting, and sequencing of cells based on specific nucleic acid biomarkers. This review highlights applications that immediately benefit from the approach, biological questions that can be addressed for the first time with it, and considerations for building successful workflows.

Efficient extraction of oil from droplet microfluidic emulsions

Droplet microfluidic techniques can perform large numbers of single molecule and cell reactions but often require controlled, periodic flow to merge, split, and sort droplets. Here, we describe a simple method to convert aperiodic flows into periodic ones. Using an oil extraction module, we efficiently remove oil from emulsions to readjust the droplet volume fraction, velocity, and packing, producing periodic flows. The extractor acts as a universal adaptor to connect microfluidic modules that do not operate under identical flow conditions, such as droplet generators, incubators, and merger devices.

Measurement of copy number variation in single cancer cells using rapid-emulsification digital droplet MDA

Uniform amplification of low input DNA is important for applications across biology, including single-cell genomics, forensic science, and microbial and viral sequencing. However, the requisite biochemical amplification methods are prone to bias, skewing sequence proportions and obscuring signals relating to copy number. Digital droplet multiple displacement amplification enables uniform amplification, but requires expert knowledge of microfluidics to generate monodisperse emulsions. In addition, existing microfluidic methods are tedious and labor intensive for preparing many samples. Here, we introduce rapid emulsification multiple displacement amplification, a method to generate monodisperse droplets with a hand-held syringe and hierarchical droplet splitter. While conventional microfluidic devices require >10 minutes to emulsify a sample, our system takes tens of seconds and yields data of equivalent quality. We demonstrate the approach by using it to accurately measure copy number variation in single cancer cells.

Perchlorate Reductase Is Distinguished by Active Site Aromatic Gate Residues

Perchlorate is an important ion on both Earth and Mars. Perchlorate reductase (PcrAB), a specialized member of the dimethylsulfoxide reductase superfamily, catalyzes the first step of microbial perchlorate respiration, but little is known about the biochemistry, specificity, structure, and mechanism of PcrAB. Here we characterize the biophysics and phylogeny of this enzyme and report the 1.86-Å resolution PcrAB complex crystal structure. Biochemical analysis revealed a relatively high perchlorate affinity (Km = 6 μm) and a characteristic substrate inhibition compared with the highly similar respiratory nitrate reductase NarGHI, which has a relatively much lower affinity for perchlorate (Km = 1.1 mm) and no substrate inhibition. Structural analysis of oxidized and reduced PcrAB with and without the substrate analog SeO3 (2-) bound to the active site identified key residues in the positively charged and funnel-shaped substrate access tunnel that gated substrate entrance and product release while trapping transiently produced chlorate. The structures suggest gating was associated with shifts of a Phe residue between open and closed conformations plus an Asp residue carboxylate shift between monodentate and bidentate coordination to the active site molybdenum atom. Taken together, structural and mutational analyses of gate residues suggest key roles of these gate residues for substrate entrance and product release. Our combined results provide the first detailed structural insight into the mechanism of biological perchlorate reduction, a critical component of the chlorine redox cycle on Earth.

Genetic dissection of chlorate respiration in Pseudomonas stutzeri PDA reveals syntrophic (per)chlorate reduction

Genes important for growth of Pseudomonas stutzeri PDA on chlorate were identified using a randomly DNA bar-coded transposon mutant library. During chlorate reduction, mutations in genes encoding the chlorate reductase clrABC, predicted molybdopterin cofactor chaperon clrD, molybdopterin biosynthesis and two genes of unknown function (clrE, clrF) had fitness defects in pooled mutant assays (Bar-seq). Markerless in-frame deletions confirmed that clrA, clrB and clrC were essential for chlorate reduction, while clrD, clrE and clrF had less severe growth defects. Interestingly, the key detoxification gene cld was essential for chlorate reduction in isogenic pure culture experiments, but showed only minor fitness defects in Bar-seq experiments. We hypothesized this was enabled through chlorite dismutation by the community, as most strains in the Bar-seq library contained an intact cld. In support of this, Δcld grew with wild-type PDA or ΔclrA, and purified Cld also restored growth to the Δcld mutant. Expanding on this, wild-type PDA and a Δcld mutant of the perchlorate reducer Azospira suillum PS grew on perchlorate in co-culture, but not individually. These results demonstrate that co-occurrence of cld and a chloroxyanion reductase within a single organism is not necessary and raises the possibility of syntrophic (per)chlorate respiration in the environment.

Chlorate reduction in Shewanella algae ACDC is a recently acquired metabolism characterized by gene loss, suboptimal regulation and oxidative stress

Previous work on respiratory chlorate reduction has biochemically identified the terminal reductase ClrABC and the chlorite detoxifying enzyme Cld. In Shewanella algae ACDC, genes encoding these enzymes reside on composite transposons whose core we refer to as the chlorate reduction composite transposon interior (CRI). To better understand this metabolism in ACDC, we used RNA-seq and proteomics to predict carbon and electron flow during chlorate reduction and posit that formate is an important electron carrier with lactate as the electron donor, but that NADH predominates on acetate. Chlorate-specific transcription of electron transport chain components or the CRI was not observed, but clr and cld transcription was attenuated by oxygen. The major chlorate-specific response related to oxidative stress and was indicative of reactive chlorine species production. A genetic system based on rpsL-streptomycin counter selection was developed to further dissect the metabolism, but ACDC readily lost the CRI via homologous recombination of the composite transposon's flanking insertion sequences. An engineered strain containing a single chromosomal CRI did not grow on chlorate, but overexpression of cld and its neighbouring cytochrome c restored growth. We postulate that the recently acquired CRI underwent copy-number expansion to circumvent insufficient expression of key genes in the pathway.

Toward a mechanistic understanding of anaerobic nitrate-dependent iron oxidation: balancing electron uptake and detoxification

The anaerobic oxidation of Fe(II) by subsurface microorganisms is an important part of biogeochemical cycling in the environment, but the biochemical mechanisms used to couple iron oxidation to nitrate respiration are not well understood. Based on our own work and the evidence available in the literature, we propose a mechanistic model for anaerobic nitrate-dependent iron oxidation. We suggest that anaerobic iron-oxidizing microorganisms likely exist along a continuum including: (1) bacteria that inadvertently oxidize Fe(II) by abiotic or biotic reactions with enzymes or chemical intermediates in their metabolic pathways (e.g., denitrification) and suffer from toxicity or energetic penalty, (2) Fe(II) tolerant bacteria that gain little or no growth benefit from iron oxidation but can manage the toxic reactions, and (3) bacteria that efficiently accept electrons from Fe(II) to gain a growth advantage while preventing or mitigating the toxic reactions. Predictions of the proposed model are highlighted and experimental approaches are discussed.