Search

High-throughput phenotypic screens using biochemical perturbations and high-content readouts are constrained by limitations of scale. To address this, we establish a method of pooling exogenous perturbations followed by computational deconvolution to reduce required sample size, labor and cost. We demonstrate the increased efficiency of compressed experimental designs compared to conventional approaches through benchmarking with a bioactive small-molecule library and a high-content imaging readout. We then apply compressed screening in two biological discovery campaigns. In the first, we use early-passage pancreatic cancer organoids to map transcriptional responses to a library of recombinant tumor microenvironment protein ligands, uncovering reproducible phenotypic shifts induced by specific ligands distinct from canonical reference signatures and correlated with clinical outcome. In the second, we identify the pleotropic modulatory effects of a chemical compound library with known mechanisms of action on primary human peripheral blood mononuclear cell immune responses. In sum, our approach empowers phenotypic screens with information-rich readouts to advance drug discovery efforts and basic biological inquiry., (© 2024. The Author(s).)

Bacterial vaginosis (BV), a common syndrome characterized by Lactobacillus-deficient vaginal microbiota, is associated with adverse health outcomes. BV often recurs after standard antibiotic therapy in part because antibiotics promote microbiota dominance by Lactobacillus iners instead of Lactobacillus crispatus, which has more beneficial health associations. Strategies to promote L. crispatus and inhibit L. iners are thus needed. We show that oleic acid (OA) and similar long-chain fatty acids simultaneously inhibit L. iners and enhance L. crispatus growth. These phenotypes require OA-inducible genes conserved in L. crispatus and related lactobacilli, including an oleate hydratase (ohyA) and putative fatty acid efflux pump (farE). FarE mediates OA resistance, while OhyA is robustly active in the vaginal microbiota and enhances bacterial fitness by biochemically sequestering OA in a derivative form only ohyA-harboring organisms can exploit. OA promotes L. crispatus dominance more effectively than antibiotics in an in vitro BV model, suggesting a metabolite-based treatment approach., Competing Interests: Declaration of interests M.Z., S.M.B., P.C.B., and D.S.K. are co-inventors on a patent related to this work. P.C.B. serves as a consultant and equity holder of companies in the microfluidics and life sciences industries, including 10x Genomics, GALT/Isolation Bio, Celsius Therapeutics, Next Gen Diagnostics, Cache DNA, Concerto Biosciences, Amber Bio, Stately, Ramona Optics, and Bifrost Biosystems. D.S.K. serves as an equity holder of Day Zero Diagnostics., (Copyright © 2024 The Authors. Published by Elsevier Inc. All rights reserved.)

Full-length RNA-sequencing methods using long-read technologies can capture complete transcript isoforms, but their throughput is limited. We introduce multiplexed arrays isoform sequencing (MAS-ISO-seq), a technique for programmably concatenating complementary DNAs (cDNAs) into molecules optimal for long-read sequencing, increasing the throughput >15-fold to nearly 40 million cDNA reads per run on the Sequel IIe sequencer. When applied to single-cell RNA sequencing of tumor-infiltrating T cells, MAS-ISO-seq demonstrated a 12- to 32-fold increase in the discovery of differentially spliced genes., (© 2023. The Author(s), under exclusive licence to Springer Nature America, Inc.)

Proton leakage from organelles is a common signal for noncanonical light chain 3B (LC3B) lipidation and inflammasome activation, processes induced upon stimulator of interferon genes (STING) activation. On the basis of structural analysis, we hypothesized that human STING is a proton channel. Indeed, we found that STING activation induced a pH increase in the Golgi and that STING reconstituted in liposomes enabled transmembrane proton transport. Compound 53 (C53), a STING agonist that binds the putative channel interface, blocked STING-induced proton flux in the Golgi and in liposomes. STING-induced LC3B lipidation and inflammasome activation were also inhibited by C53, suggesting that STING's channel activity is critical for these two processes. Thus, STING's interferon-induction function can be decoupled from its roles in LC3B lipidation and inflammasome activation.

The infection of mammalian cells by viruses and innate immune responses to infection are spatiotemporally organized processes. Cytosolic RNA sensors trigger nuclear translocation of the transcription factor interferon regulatory factor 3 (IRF3) and consequent induction of host immune responses to RNA viruses. Previous genetic screens for factors involved in viral sensing did not resolve changes in the subcellular localization of host or viral proteins. Here, we increased the throughput of our optical pooled screening technology by over fourfold. This allowed us to carry out a genome-wide CRISPR knockout screen using high-resolution multiparameter imaging of cellular responses to Sendai virus infection coupled with in situ cDNA sequencing by synthesis (SBS) to identify 80,408 single guide RNAs (sgRNAs) in 10,366,390 cells-over an order of magnitude more genomic perturbations than demonstrated previously using an in situ SBS readout. By ranking perturbations using human-designed and deep learning image feature scores, we identified regulators of IRF3 translocation, Sendai virus localization, and peroxisomal biogenesis. Among the hits, we found that ATP13A1, an ER-localized P5A-type ATPase, is essential for viral sensing and is required for targeting of mitochondrial antiviral signaling protein (MAVS) to mitochondrial membranes where MAVS must be localized for effective signaling through retinoic acid-inducible gene I (RIG-I). The ability to carry out genome-wide pooled screens with complex high-resolution image-based phenotyping dramatically expands the scope of functional genomics approaches.

Objective: To describe the genomic analysis and epidemiologic response related to a slow and prolonged methicillin-resistant Staphylococcus aureus (MRSA) outbreak., Design: Prospective observational study., Setting: Neonatal intensive care unit (NICU)., Methods: We conducted an epidemiologic investigation of a NICU MRSA outbreak involving serial baby and staff screening to identify opportunities for decolonization. Whole-genome sequencing was performed on MRSA isolates., Results: A NICU with excellent hand hygiene compliance and longstanding minimal healthcare-associated infections experienced an MRSA outbreak involving 15 babies and 6 healthcare personnel (HCP). In total, 12 cases occurred slowly over a 1-year period (mean, 30.7 days apart) followed by 3 additional cases 7 months later. Multiple progressive infection prevention interventions were implemented, including contact precautions and cohorting of MRSA-positive babies, hand hygiene observers, enhanced environmental cleaning, screening of babies and staff, and decolonization of carriers. Only decolonization of HCP found to be persistent carriers of MRSA was successful in stopping transmission and ending the outbreak. Genomic analyses identified bidirectional transmission between babies and HCP during the outbreak., Conclusions: In comparison to fast outbreaks, outbreaks that are "slow and sustained" may be more common to units with strong existing infection prevention practices such that a series of breaches have to align to result in a case. We identified a slow outbreak that persisted among staff and babies and was only stopped by identifying and decolonizing persistent MRSA carriage among staff. A repeated decolonization regimen was successful in allowing previously persistent carriers to safely continue work duties.

Detection of microbial cell-free DNA (cfDNA) circulating in the bloodstream has emerged as a promising new approach for diagnosing infection. Microbial diagnostics based on cfDNA require assays that can detect rare and highly fragmented pathogen nucleic acids. We now report WATSON (Whole-genome Assay using Tiled Surveillance Of Nucleic acids), a method to detect low amounts of pathogen cfDNA that couples pooled amplification of genomic targets tiled across the genome with pooled CRISPR/Cas13-based detection of these targets. We demonstrate that this strategy of tiling improves cfDNA detection compared to amplification and detection of a single targeted locus. WATSON can detect cfDNA from Mycobacterium tuberculosis in plasma of patients with active pulmonary tuberculosis, a disease that urgently needs accurate, minimally-invasive, field-deployable diagnostics. We thus demonstrate the potential for translating WATSON to a lateral flow platform. WATSON demonstrates the ability to capitalize on the strengths of targeting microbial cfDNA to address the need for point-of-care diagnostic tests for infectious diseases., (© 2023. The Author(s).)

Cells sense a wide variety of signals and respond by adopting complex transcriptional states. Most single-cell profiling is carried out today at cellular baseline, blind to cells' potential spectrum of functional responses. Exploring the space of cellular responses experimentally requires access to a large combinatorial perturbation space. Single-cell genomics coupled with multiplexing techniques provide a useful tool for characterizing cell states across several experimental conditions. However, current multiplexing strategies require programmatic handling of many samples in macroscale arrayed formats, precluding their application in large-scale combinatorial analysis. Here, we introduce StimDrop, a method that combines antibody-based cell barcoding with parallel droplet processing to automatically formulate cell population × stimulus combinations in a microfluidic device. We applied StimDrop to profile the effects of 512 sequential stimulation conditions on human dendritic cells. Our results demonstrate that priming with viral ligands potentiates hyperinflammatory responses to a second stimulus, and show transcriptional signatures consistent with this phenomenon in myeloid cells of patients with severe COVID-19.

Understanding the basis for cellular growth, proliferation, and function requires determining the roles of essential genes in diverse cellular processes, including visualizing their contributions to cellular organization and morphology. Here, we combined pooled CRISPR-Cas9-based functional screening of 5,072 fitness-conferring genes in human HeLa cells with microscopy-based imaging of DNA, the DNA damage response, actin, and microtubules. Analysis of >31 million individual cells identified measurable phenotypes for >90% of gene knockouts, implicating gene targets in specific cellular processes. Clustering of phenotypic similarities based on hundreds of quantitative parameters further revealed co-functional genes across diverse cellular activities, providing predictions for gene functions and associations. By conducting pooled live-cell screening of ∼450,000 cell division events for 239 genes, we additionally identified diverse genes with functional contributions to chromosome segregation. Our work establishes a resource detailing the consequences of disrupting core cellular processes that represents the functional landscape of essential human genes., Competing Interests: Declaration of interests P.C.B. is a consultant to and/or equity holder in companies in the life sciences industries including 10X Genomics, GALT, Celsius Therapeutics, Next Generation Diagnostics, Cache DNA, Concerto Biosciences, and Stately. P.C.B.’s laboratory receives research funding from Calico Life Sciences and Merck for work related to genetic screening. The Broad Institute and MIT have filed U.S. patent applications on work described here and may seek to license the technology., (Copyright © 2022 Elsevier Inc. All rights reserved.)

Spatial structure in biology, spanning molecular, organellular, cellular, tissue, and organismal scales, is encoded through a combination of genetic and epigenetic factors in individual cells. Microscopy remains the most direct approach to exploring the intricate spatial complexity defining biological systems and the structured dynamic responses of these systems to perturbations. Genetic screens with deep single-cell profiling via image features or gene expression programs have the capacity to show how biological systems work in detail by cataloging many cellular phenotypes with one experimental assay. Microscopy-based cellular profiling provides information complementary to next-generation sequencing (NGS) profiling and has only recently become compatible with large-scale genetic screens. Optical screening now offers the scale needed for systematic characterization and is poised for further scale-up. We discuss how these methodologies, together with emerging technologies for genetic perturbation and microscopy-based multiplexed molecular phenotyping, are powering new approaches to reveal genotype-phenotype relationships., (© 2022 The Authors. Published under the terms of the CC BY 4.0 license.)

Design of nucleic acid-based viral diagnostics typically follows heuristic rules and, to contend with viral variation, focuses on a genome's conserved regions. A design process could, instead, directly optimize diagnostic effectiveness using a learned model of sensitivity for targets and their variants. Toward that goal, we screen 19,209 diagnostic-target pairs, concentrated on CRISPR-based diagnostics, and train a deep neural network to accurately predict diagnostic readout. We join this model with combinatorial optimization to maximize sensitivity over the full spectrum of a virus's genomic variation. We introduce Activity-informed Design with All-inclusive Patrolling of Targets (ADAPT), a system for automated design, and use it to design diagnostics for 1,933 vertebrate-infecting viral species within 2 hours for most species and within 24 hours for all but three. We experimentally show that ADAPT's designs are sensitive and specific to the lineage level and permit lower limits of detection, across a virus's variation, than the outputs of standard design techniques. Our strategy could facilitate a proactive resource of assays for detecting pathogens., (© 2022. The Author(s).)

The coronavirus disease 2019 (COVID-19) pandemic has demonstrated a clear need for high-throughput, multiplexed and sensitive assays for detecting severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) and other respiratory viruses and their emerging variants. Here, we present a cost-effective virus and variant detection platform, called microfluidic Combinatorial Arrayed Reactions for Multiplexed Evaluation of Nucleic acids (mCARMEN), which combines CRISPR-based diagnostics and microfluidics with a streamlined workflow for clinical use. We developed the mCARMEN respiratory virus panel to test for up to 21 viruses, including SARS-CoV-2, other coronaviruses and both influenza strains, and demonstrated its diagnostic-grade performance on 525 patient specimens in an academic setting and 166 specimens in a clinical setting. We further developed an mCARMEN panel to enable the identification of 6 SARS-CoV-2 variant lineages, including Delta and Omicron, and evaluated it on 2,088 patient specimens with near-perfect concordance to sequencing-based variant classification. Lastly, we implemented a combined Cas13 and Cas12 approach that enables quantitative measurement of SARS-CoV-2 and influenza A viral copies in samples. The mCARMEN platform enables high-throughput surveillance of multiple viruses and variants simultaneously, enabling rapid detection of SARS-CoV-2 variants., (© 2022. The Author(s).)

Rapid and accurate diagnosis of infections is fundamental to individual patient care and public health management. Nucleic acid detection methods are critical to this effort, but are limited either in the breadth of pathogens targeted or by the expertise and infrastructure required. We present here a high-throughput system that enables rapid identification of bacterial pathogens, bCARMEN, which utilizes: (1) modular CRISPR-Cas13-based nucleic acid detection with enhanced sensitivity and specificity; and (2) a droplet microfluidic system that enables thousands of simultaneous, spatially multiplexed detection reactions at nanoliter volumes; and (3) a novel preamplification strategy that further enhances sensitivity and specificity. We demonstrate bCARMEN is capable of detecting and discriminating 52 clinically relevant bacterial species and several key antibiotic resistance genes. We further develop a simple proof of principle workflow using stabilized reagents and cell phone camera optical readout, opening up the possibility of a rapid point-of-care multiplexed bacterial pathogen identification and antibiotic susceptibility testing., (© The Author(s) 2022. Published by Oxford University Press on behalf of the National Academy of Sciences.)

Discovery of the genetic components underpinning fundamental and disease-related processes is being rapidly accelerated by combining efficient, programmable genetic engineering with phenotypic readouts of high spatial, temporal and/or molecular resolution. Microscopy is a fundamental tool for studying cell biology, but its lack of high-throughput sequence readouts hinders integration in large-scale genetic screens. Optical pooled screens using in situ sequencing provide massively scalable integration of barcoded lentiviral libraries (e.g., CRISPR perturbation libraries) with high-content imaging assays, including dynamic processes in live cells. The protocol uses standard lentiviral vectors and molecular biology, providing single-cell resolution of phenotype and engineered genotype, scalability to millions of cells and accurate sequence reads sufficient to distinguish >10 6 perturbations. In situ amplification takes ~2 d, while sequencing can be performed in ~1.5 h per cycle. The image analysis pipeline provided enables fully parallel automated sequencing analysis using a cloud or cluster computing environment., (© 2022. The Author(s), under exclusive licence to Springer Nature Limited.)

Interspecies interactions shape the structure and function of microbial communities. In particular, positive, growth-promoting interactions can substantially affect the diversity and productivity of natural and engineered communities. However, the prevalence of positive interactions and the conditions in which they occur are not well understood. To address this knowledge gap, we used kChip, an ultrahigh-throughput coculture platform, to measure 180,408 interactions among 20 soil bacteria across 40 carbon environments. We find that positive interactions, often described to be rare, occur commonly and primarily as parasitisms between strains that differ in their carbon consumption profiles. Notably, nongrowing strains are almost always promoted by strongly growing strains (85%), suggesting a simple positive interaction–mediated approach for cultivation, microbiome engineering, and microbial consortium design.

DNA is an ultrahigh-density storage medium that could meet exponentially growing worldwide demand for archival data storage if DNA synthesis costs declined sufficiently and if random access of files within exabyte-to-yottabyte-scale DNA data pools were feasible. Here, we demonstrate a path to overcome the second barrier by encapsulating data-encoding DNA file sequences within impervious silica capsules that are surface labelled with single-stranded DNA barcodes. Barcodes are chosen to represent file metadata, enabling selection of sets of files with Boolean logic directly, without use of amplification. We demonstrate random access of image files from a prototypical 2-kilobyte image database using fluorescence sorting with selection sensitivity of one in 10 6 files, which thereby enables one in 10 6N selection capability using N optical channels. Our strategy thereby offers a scalable concept for random access of archival files in large-scale molecular datasets., (© 2021. The Author(s), under exclusive licence to Springer Nature Limited.)

Antibiotic resistance is a worldwide and growing clinical problem. With limited drug development in the antibacterial space, combination therapy has emerged as a promising strategy to combat multidrug-resistant bacteria. Antibacterial combinations can improve antibiotic efficacy and suppress antibacterial resistance through independent, synergistic, or even antagonistic activities. Combination therapies are famously used to treat viral and mycobacterial infections and cancer. However, antibacterial combinations are only now emerging as a common treatment strategy for other bacterial infections owing to challenges in their discovery, development, regulatory approval, and commercial/clinical deployment. Here, we focus on discovery-where the sheer scale of combinatorial chemical spaces represents a significant challenge-and discuss how combination therapy can impact the treatment of bacterial infections. Despite these challenges, recent advancements, including new in silico methods, theoretical frameworks, and microfluidic platforms, are poised to identify the new and efficacious antibacterial combinations needed to revitalize the antibacterial drug pipeline., (© 2021 The Authors. Annals of the New York Academy of Sciences published by Wiley Periodicals LLC on behalf of New York Academy of Sciences.)

Bacterial sepsis and severe COVID-19 share similar clinical manifestations and are both associated with dysregulation of the myeloid cell compartment. We previously reported an expanded CD14 + monocyte state, MS1, in patients with bacterial sepsis and validated expansion of this cell subpopulation in publicly available transcriptomics data. Here, using published datasets, we show that the gene expression program associated with MS1 correlated with sepsis severity and was up-regulated in monocytes from patients with severe COVID-19. To examine the ontogeny and function of MS1 cells, we developed a cellular model for inducing CD14 + MS1 monocytes from healthy bone marrow hematopoietic stem and progenitor cells (HSPCs). We found that plasma from patients with bacterial sepsis or COVID-19 induced myelopoiesis in HSPCs in vitro and expression of the MS1 gene program in monocytes and neutrophils that differentiated from these HSPCs. Furthermore, we found that plasma concentrations of IL-6, and to a lesser extent IL-10, correlated with increased myeloid cell output from HSPCs in vitro and enhanced expression of the MS1 gene program. We validated the requirement for these two cytokines to induce the MS1 gene program through CRISPR-Cas9 editing of their receptors in HSPCs. Using this cellular model system, we demonstrated that induced MS1 cells were broadly immunosuppressive and showed decreased responsiveness to stimulation with a synthetic RNA analog. Our in vitro study suggests a potential role for systemic cytokines in inducing myelopoiesis during severe bacterial or SARS-CoV-2 infection., (Copyright © 2021 The Authors, some rights reserved; exclusive licensee American Association for the Advancement of Science. No claim to original U.S. Government Works. Distributed under a Creative Commons Attribution License 4.0 (CC BY).)

DNA single-strand breaks (SSBs), or "nicks," are the most common form of DNA damage. Oxidative stress, endogenous enzyme activities, and other processes cause tens of thousands of nicks per cell per day. Accumulation of nicks, caused by high rates of occurrence or defects in repair enzymes, has been implicated in multiple diseases. However, improved methods for nick analysis are needed to characterize the mechanisms of these processes and learn how the location and number of nicks affect cells, disease progression, and health outcomes. In addition to natural processes, including DNA repair, leading genome editing technologies rely on nuclease activity, including nick generation, at specific target sites. There is currently a pressing need for methods to study off-target nicking activity genome-wide to evaluate the side effects of emerging genome editing tools on cells and organisms. Here, we developed a new method, DENT-seq, for efficient strand-specific profiling of nicks in complex DNA samples with single-nucleotide resolution and low false-positive rates. DENT-seq produces a single deep sequence data set enriched for reads near nick sites and establishes a readily detectable mutational signal that allows for determination of the nick site and strand with single-base resolution at penetrance as low as one strand per thousand. We apply DENT-seq to profile the off-target activity of the Nb.BsmI nicking endonuclease and an engineered spCas9 nickase. DENT-seq will be useful in exploring the activity of engineered nucleases in genome editing and other biotechnological applications as well as spontaneous and therapeutic-associated strand breaks., (© 2021 Elacqua et al.; Published by Cold Spring Harbor Laboratory Press.)

Existing cancer benchmark data sets for human sequencing data use germline variants, synthetic methods, or expensive validations, none of which are satisfactory for providing a large collection of true somatic variation across a whole genome. Here we propose a data set, Lineage derived Somatic Truth (LinST), of short somatic mutations in the HT115 colon cancer cell-line, that are validated using a known cell lineage that includes thousands of mutations and a high confidence region covering 2.7 gigabases per sample.

Background: Many biological processes, such as cancer metastasis, organismal development, and acquisition of resistance to cytotoxic therapy, rely on the emergence of rare sub-clones from a larger population. Understanding how the genetic and epigenetic features of diverse clones affect clonal fitness provides insight into molecular mechanisms underlying selective processes. While large-scale barcoding with NGS readout has facilitated cellular fitness assessment at the population level, this approach does not support characterization of clones prior to selection. Single-cell genomics methods provide high biological resolution, but are challenging to scale across large populations to probe rare clones and are destructive, limiting further functional analysis of important clones., Results: Here, we develop CloneSifter, a methodology for tracking and enriching rare clones throughout their response to selection. CloneSifter utilizes a CRISPR sgRNA-barcode library that facilitates the isolation of viable cells from specific clones within the barcoded population using a sequence-specific retrieval reporter. We demonstrate that CloneSifter can measure clonal fitness of cancer cell models in vitro and retrieve targeted clones at abundance as low as 1 in 1883 in a heterogeneous cell population., Conclusions: CloneSifter provides a means to track and access specific and rare clones of interest across dynamic changes in population structure to comprehensively explore the basis of these changes.

The great majority of globally circulating pathogens go undetected, undermining patient care and hindering outbreak preparedness and response. To enable routine surveillance and comprehensive diagnostic applications, there is a need for detection technologies that can scale to test many samples 1-3 while simultaneously testing for many pathogens 4-6 . Here, we develop Combinatorial Arrayed Reactions for Multiplexed Evaluation of Nucleic acids (CARMEN), a platform for scalable, multiplexed pathogen detection. In the CARMEN platform, nanolitre droplets containing CRISPR-based nucleic acid detection reagents 7 self-organize in a microwell array 8 to pair with droplets of amplified samples, testing each sample against each CRISPR RNA (crRNA) in replicate. The combination of CARMEN and Cas13 detection (CARMEN-Cas13) enables robust testing of more than 4,500 crRNA-target pairs on a single array. Using CARMEN-Cas13, we developed a multiplexed assay that simultaneously differentiates all 169 human-associated viruses with at least 10 published genome sequences and rapidly incorporated an additional crRNA to detect the causative agent of the 2020 COVID-19 pandemic. CARMEN-Cas13 further enables comprehensive subtyping of influenza A strains and multiplexed identification of dozens of HIV drug-resistance mutations. The intrinsic multiplexing and throughput capabilities of CARMEN make it practical to scale, as miniaturization decreases reagent cost per test by more than 300-fold. Scalable, highly multiplexed CRISPR-based nucleic acid detection shifts diagnostic and surveillance efforts from targeted testing of high-priority samples to comprehensive testing of large sample sets, greatly benefiting patients and public health 9-11 .

Lung regeneration occurs in a variety of adult mammals after surgical removal of one lung (pneumonectomy). Previous studies of murine post-pneumonectomy lung growth have identified regenerative "hotspots" in subpleural alveolar ducts; however, the cell-types participating in this process remain unclear. To identify the single cells participating in post-pneumonectomy lung growth, we used laser microdissection, enzymatic digestion and microfluidic isolation. Single-cell transcriptional analysis of the murine alveolar duct cells was performed using the C1 integrated fluidic circuit (Fluidigm) and a custom PCR panel designed for lung growth and repair genes. The multi-dimensional data set was analyzed using visualization software based on the tSNE algorithm. The analysis identified 6 cell clusters; 1 cell cluster was present only after pneumonectomy. This post-pneumonectomy cluster was significantly less transcriptionally active than 3 other clusters and may represent a transitional cell population. A provisional cluster identity for 4 of the 6 cell clusters was obtained by embedding bulk transcriptional data into the tSNE analysis. The transcriptional pattern of the 6 clusters was further analyzed for genes associated with lung repair, matrix production, and angiogenesis. The data demonstrated that multiple cell-types (clusters) transcribed genes linked to these basic functions. We conclude that the coordinated gene expression across multiple cell clusters is likely a response to a shared regenerative microenvironment within the subpleural alveolar ducts., (Copyright © 2020 Ysasi, Bennett, Wagner, Valenzuela, Servais, Tsuda, Pyne, Li, Grimsby, Pokharel, Livak, Ackermann, Blainey and Mentzer.)

Dysregulation of the immune response to bacterial infection can lead to sepsis, a condition with high mortality. Multiple whole-blood gene-expression studies have defined sepsis-associated molecular signatures, but have not resolved changes in transcriptional states of specific cell types. Here, we used single-cell RNA-sequencing to profile the blood of people with sepsis (n = 29) across three clinical cohorts with corresponding controls (n = 36). We profiled total peripheral blood mononuclear cells (PBMCs, 106,545 cells) and dendritic cells (19,806 cells) across all subjects and, on the basis of clustering of their gene-expression profiles, defined 16 immune-cell states. We identified a unique CD14 + monocyte state that is expanded in people with sepsis and validated its power in distinguishing these individuals from controls using public transcriptomic data from subjects with different disease etiologies and from multiple geographic locations (18 cohorts, n = 1,467 subjects). We identified a panel of surface markers for isolation and quantification of the monocyte state and characterized its epigenomic and functional phenotypes, and propose a model for its induction from human bone marrow. This study demonstrates the utility of single-cell genomics in discovering disease-associated cytologic signatures and provides insight into the cellular basis of immune dysregulation in bacterial sepsis.

Clustered regularly interspaced short palindromic repeats (CRISPR)-associated 9 (Cas9) systems have been effectively harnessed to engineer the genomes of organisms from across the tree of life. Nearly all currently characterized Cas9 proteins are derived from mesophilic bacteria, and canonical Cas9 systems are challenged by applications requiring enhanced stability or elevated temperatures. We discovered IgnaviCas9, a Cas9 protein from a hyperthermophilic Ignavibacterium identified through mini-metagenomic sequencing of samples from a hot spring. IgnaviCas9 is active at temperatures up to 100 °C in vitro, which enables DNA cleavage beyond the 44 °C limit of Streptococcus pyogenes Cas9 (SpyCas9) and the 70 °C limit of both Geobacillus stearothermophilus Cas9 (GeoCas9) and Geobacillus thermodenitrificans T12 Cas9 (ThermoCas9). As a potential application of this enzyme, we demonstrate that IgnaviCas9 can be used in bacterial RNA-seq library preparation to remove unwanted cDNA from 16s ribosomal rRNA without increasing the number of steps, thus underscoring the benefits provided by its exceptional thermostability in improving molecular biology and genomic workflows. IgnaviCas9 is an exciting addition to the CRISPR-Cas9 toolbox and expands its temperature range., Competing Interests: Competing interest statement: A patent concerning the technology disclosed in this publication has been filed., (Copyright © 2019 the Author(s). Published by PNAS.)

Profiling multiple omic layers in a single cell enables the discovery and analysis of biological phenomena that are not apparent from analysis of mono-omic data. While methods for multi-omic profiling have been reported, their adoption has been limited due to high cost and complex workflows. Here, we present a simple method for joint profiling of gene expression and chromatin accessibility in tens to hundreds of single cells. We assess the quality of resulting single cell ATAC- and RNA-seq data across three cell types, examine the link between accessibility and expression at the CD3G and FTH1 loci in human primary T cells and monocytes, and compare the accuracy of clustering solutions for mono-omic and combined data. The new method allows biological laboratories to perform simultaneous profiling of gene expression and chromatin accessibility using standard reagents and instrumentation. This technique, in conjunction with other advances in multi-omic profiling, will enable highly-resolved cell state classification and more specific mechanistic hypothesis generation than is possible with mono-omic analysis., Competing Interests: CONFLICTS OF INTEREST P.C.B. is an extramural faculty member of MIT’s Koch Institute for Integrative Cancer Research and a consultant to and equity holder in two companies in the microfluidics and single-cell industries, 10X Genomics and General Automation Lab Technologies. The Broad Institute and MIT may seek to commercialize aspects of this work.

Genetic screens are critical for the systematic identification of genes underlying cellular phenotypes. Pooling gene perturbations greatly improves scalability but is not compatible with imaging of complex and dynamic cellular phenotypes. Here, we introduce a pooled approach for optical genetic screens in mammalian cells. We use targeted in situ sequencing to demultiplex a library of genetic perturbations following image-based phenotyping. We screened a set of 952 genes across millions of cells for involvement in nuclear factor κB (NF-κB) signaling by imaging the translocation of RelA (p65) to the nucleus. Screening at a single time point across 3 cell lines recovered 15 known pathway components, while repeating the screen with live-cell imaging revealed a role for Mediator complex subunits in regulating the duration of p65 nuclear retention. These results establish a highly multiplexed approach to image-based screens of spatially and temporally defined phenotypes with pooled libraries., (Copyright © 2019 Elsevier Inc. All rights reserved.)

Synapses contain hundreds of distinct proteins whose heterogeneous expression levels are determinants of synaptic plasticity and signal transmission relevant to a range of diseases. Here, we use diffusible nucleic acid imaging probes to profile neuronal synapses using multiplexed confocal and super-resolution microscopy. Confocal imaging is performed using high-affinity locked nucleic acid imaging probes that stably yet reversibly bind to oligonucleotides conjugated to antibodies and peptides. Super-resolution PAINT imaging of the same targets is performed using low-affinity DNA imaging probes to resolve nanometer-scale synaptic protein organization across nine distinct protein targets. Our approach enables the quantitative analysis of thousands of synapses in neuronal culture to identify putative synaptic sub-types and co-localization patterns from one dozen proteins. Application to characterize synaptic reorganization following neuronal activity blockade reveals coordinated upregulation of the post-synaptic proteins PSD-95, SHANK3 and Homer-1b/c, as well as increased correlation between synaptic markers in the active and synaptic vesicle zones.

Microbial communities have numerous potential applications in biotechnology, agriculture, and medicine. Nevertheless, the limited accuracy with which we can predict interspecies interactions and environmental dependencies hinders efforts to rationally engineer beneficial consortia. Empirical screening is a complementary approach wherein synthetic communities are combinatorially constructed and assayed in high throughput. However, assembling many combinations of microbes is logistically complex and difficult to achieve on a timescale commensurate with microbial growth. Here, we introduce the kChip, a droplets-based platform that performs rapid, massively parallel, bottom-up construction and screening of synthetic microbial communities. We first show that the kChip enables phenotypic characterization of microbes across environmental conditions. Next, in a screen of ∼100,000 multispecies communities comprising up to 19 soil isolates, we identified sets that promote the growth of the model plant symbiont Herbaspirillum frisingense in a manner robust to carbon source variation and the presence of additional species. Broadly, kChip screening can identify multispecies consortia possessing any optically assayable function, including facilitation of biocontrol agents, suppression of pathogens, degradation of recalcitrant substrates, and robustness of these functions to perturbation, with many applications across basic and applied microbial ecology., Competing Interests: Conflict of interest statement: P.C.B. is an extramural faculty member of MIT’s Koch Institute for Integrative Cancer Research and a consultant to and equity holder in two companies in the microfluidics industry, 10X Genomics and General Automation Lab Technologies. The Broad Institute and MIT may seek to commercialize aspects of this work, and related applications for intellectual property have been filed.

Bacterial populations that colonize a host can play important roles in host health, including serving as a reservoir that transmits to other hosts and from which invasive strains emerge, thus emphasizing the importance of understanding rates of acquisition and clearance of colonizing populations. Studies of colonization dynamics have been based on assessment of whether serial samples represent a single population or distinct colonization events. With the use of whole genome sequencing to determine genetic distance between isolates, a common solution to estimate acquisition and clearance rates has been to assume a fixed genetic distance threshold below which isolates are considered to represent the same strain. However, this approach is often inadequate to account for the diversity of the underlying within-host evolving population, the time intervals between consecutive measurements, and the uncertainty in the estimated acquisition and clearance rates. Here, we present a fully Bayesian model that provides probabilities of whether two strains should be considered the same, allowing us to determine bacterial clearance and acquisition from genomes sampled over time. Our method explicitly models the within-host variation using population genetic simulation, and the inference is done using a combination of Approximate Bayesian Computation (ABC) and Markov Chain Monte Carlo (MCMC). We validate the method with multiple carefully conducted simulations and demonstrate its use in practice by analyzing a collection of methicillin resistant Staphylococcus aureus (MRSA) isolates from a large recently completed longitudinal clinical study. An R-code implementation of the method is freely available at: https://github.com/mjarvenpaa/bacterial-colonization-model., Competing Interests: The authors have declared that no relevant competing interests exist.

Methanogenic archaea are major contributors to the global carbon cycle and were long thought to belong exclusively to the euryarchaeal phylum. Discovery of the methanogenesis gene cluster methyl-coenzyme M reductase (Mcr) in the Bathyarchaeota, and thereafter the Verstraetearchaeota, led to a paradigm shift, pushing back the evolutionary origin of methanogenesis to predate that of the Euryarchaeota. The methylotrophic methanogenesis found in the non-Euryarchaota distinguished itself from the predominantly hydrogenotrophic methanogens found in euryarchaeal orders as the former do not couple methanogenesis to carbon fixation through the reductive acetyl-CoA [Wood-Ljungdahl pathway (WLP)], which was interpreted as evidence for independent evolution of the two methanogenesis pathways. Here, we report the discovery of a complete and divergent hydrogenotrophic methanogenesis pathway in a thermophilic order of the Verstraetearchaeota, which we have named Candidatus Methanohydrogenales, as well as the presence of the WLP in the crenarchaeal order Desulfurococcales. Our findings support the ancient origin of hydrogenotrophic methanogenesis, suggest that methylotrophic methanogenesis might be a later adaptation of specific orders, and provide insight into how the transition from hydrogenotrophic to methylotrophic methanogenesis might have occurred., Competing Interests: Conflict of interest statement: S.R.Q. is a shareholder of Fluidigm Corporation. T.W. and M.S. are coauthors on a 2017 software assessment paper., (Copyright © 2019 the Author(s). Published by PNAS.)

Lung cancer is closely associated with chronic inflammation, but the causes of inflammation and the specific immune mediators have not been fully elucidated. The lung is a mucosal tissue colonized by a diverse bacterial community, and pulmonary infections commonly present in lung cancer patients are linked to clinical outcomes. Here, we provide evidence that local microbiota provoke inflammation associated with lung adenocarcinoma by activating lung-resident γδ T cells. Germ-free or antibiotic-treated mice were significantly protected from lung cancer development induced by Kras mutation and p53 loss. Mechanistically, commensal bacteria stimulated Myd88-dependent IL-1β and IL-23 production from myeloid cells, inducing proliferation and activation of Vγ6 + Vδ1 + γδ T cells that produced IL-17 and other effector molecules to promote inflammation and tumor cell proliferation. Our findings clearly link local microbiota-immune crosstalk to lung tumor development and thereby define key cellular and molecular mediators that may serve as effective targets in lung cancer intervention., (Copyright © 2018 Elsevier Inc. All rights reserved.)

Specialized immune cell subsets are involved in autoimmune disease, cancer immunity, and infectious disease through a diverse range of functions mediated by overlapping pathways and signals. However, subset-specific responses may not be detectable in analyses of whole blood samples, and no efficient approach for profiling cell subsets at high throughput from small samples is available. We present a low-input microfluidic system for sorting immune cells into subsets and profiling their gene expression. We validate the system's technical performance against standard subset isolation and library construction protocols and demonstrate the importance of subset-specific profiling through in vitro stimulation experiments. We show the ability of this integrated platform to identify subset-specific disease signatures by profiling four immune cell subsets in blood from patients with systemic lupus erythematosus (SLE) and matched control subjects. The platform has the potential to make multiplexed subset-specific analysis routine in many research laboratories and clinical settings.

Transcriptional profiling of thousands of single cells in parallel by RNA-seq is now routine. However, due to reliance on pooled library preparation, targeting analysis to particular cells of interest is difficult. Here, we present a multiplexed PCR method for targeted sequencing of select cells from pooled single-cell sequence libraries. We demonstrated this molecular enrichment method on multiple cell types within pooled single-cell RNA-seq libraries produced from primary human blood cells. We show how molecular enrichment can be combined with FACS to efficiently target ultra-rare cell types, such as the recently identified AXL+SIGLEC6+ dendritic cell (AS DC) subset, in order to reduce the required sequencing effort to profile single cells by 100-fold. Our results demonstrate that DNA barcodes identifying cells within pooled sequencing libraries can be used as targets to enrich for specific molecules of interest, for example reads from a set of target cells.

Mutation data reveal the dynamic equilibrium between DNA damage and repair processes in cells and are indispensable to the understanding of age-related diseases, tumor evolution, and the acquisition of drug resistance. However, available genome-wide methods have a limited ability to resolve rare somatic variants and the relationships between these variants. Here, we present lineage sequencing, a new genome sequencing approach that enables somatic event reconstruction by providing quality somatic mutation call sets with resolution as high as the single-cell level in subject lineages. Lineage sequencing entails sampling single cells from a population and sequencing subclonal sample sets derived from these cells such that knowledge of relationships among the cells can be used to jointly call variants across the sample set. This approach integrates data from multiple sequence libraries to support each variant and precisely assigns mutations to lineage segments. We applied lineage sequencing to a human colon cancer cell line with a DNA polymerase epsilon ( POLE ) proofreading deficiency (HT115) and a human retinal epithelial cell line immortalized by constitutive telomerase expression (RPE1). Cells were cultured under continuous observation to link observed single-cell phenotypes with single-cell mutation data. The high sensitivity, specificity, and resolution of the data provide a unique opportunity for quantitative analysis of variation in mutation rate, spectrum, and correlations among variants. Our data show that mutations arrive with nonuniform probability across sublineages and that DNA lesion dynamics may cause strong correlations between certain mutations., (© 2018 Brody et al.; Published by Cold Spring Harbor Laboratory Press.)

Combinatorial drug treatment strategies perturb biological networks synergistically to achieve therapeutic effects and represent major opportunities to develop advanced treatments across a variety of human disease areas. However, the discovery of new combinatorial treatments is challenged by the sheer scale of combinatorial chemical space. Here, we report a high-throughput system for nanoliter-scale phenotypic screening that formulates a chemical library in nanoliter droplet emulsions and automates the construction of chemical combinations en masse using parallel droplet processing. We applied this system to predict synergy between more than 4,000 investigational and approved drugs and a panel of 10 antibiotics against Escherichia coli , a model gram-negative pathogen. We found a range of drugs not previously indicated for infectious disease that synergize with antibiotics. Our validated hits include drugs that synergize with the antibiotics vancomycin, erythromycin, and novobiocin, which are used against gram-positive bacteria but are not effective by themselves to resolve gram-negative infections., Competing Interests: Conflict of interest statement: Broad Institute and MIT may seek to commercialize aspects of this work; related applications for intellectual property have been filed.

We reanalyze trajectories of hOGG1 repair proteins diffusing on DNA. A previous analysis of these trajectories with the popular mean-squared-displacement approach revealed only simple diffusion. Here, a new optimal estimator of diffusion coefficients reveals two-state kinetics of the protein. A simple, solvable model, in which the protein randomly switches between a loosely bound, highly mobile state and a tightly bound, less mobile state is the simplest possible dynamic model consistent with the data. It yields accurate estimates of hOGG1's (i) diffusivity in each state, uncorrupted by experimental errors arising from shot noise, motion blur and thermal fluctuations of the DNA; (ii) rates of switching between states and (iii) rate of detachment from the DNA. The protein spends roughly equal time in each state. It detaches only from the loosely bound state, with a rate that depends on pH and the salt concentration in solution, while its rates for switching between states are insensitive to both. The diffusivity in the loosely bound state depends primarily on pH and is three to ten times higher than in the tightly bound state. We propose and discuss some new experiments that take full advantage of the new tools of analysis presented here.