The phenotypic landscape of essential human genes.

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.

Optical Pooled Screens in Human Cells.

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.

Commensal Microbiota Promote Lung Cancer Development via γδ T Cells.

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.

Massively multiplexed nucleic acid detection with Cas13.

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 samples1-3 while simultaneously testing for many pathogens4-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 reagents7 self-organize in a microwell array8 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 health9-11.

A genome-wide optical pooled screen reveals regulators of cellular antiviral responses.

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.

DENT-seq for genome-wide strand-specific identification of DNA single-strand break sites with single-nucleotide resolution.

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.

Pooled genetic screens with image-based profiling.

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.

Random access DNA memory using Boolean search in an archival file storage system.

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 106 files, which thereby enables one in 106N selection capability using N optical channels. Our strategy thereby offers a scalable concept for random access of archival files in large-scale molecular datasets.

Nucleic acid cleavage with a hyperthermophilic Cas9 from an uncultured Ignavibacterium.

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.

Hydrogenotrophic methanogenesis in archaeal phylum Verstraetearchaeota reveals the shared ancestry of all methanogens.

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.

Massively parallel screening of synthetic microbial communities.

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.

Quantification of somatic mutation flow across individual cell division events by lineage sequencing.

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.

Targeting individual cells by barcode in pooled sequence libraries.

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.

Combinatorial drug discovery in nanoliter droplets.

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.

CloneSifter: enrichment of rare clones from heterogeneous cell populations.

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.

A Bayesian model of acquisition and clearance of bacterial colonization incorporating within-host variation.

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.

Single-particle trajectories reveal two-state diffusion-kinetics of hOGG1 proteins on DNA.

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.

Virtual microfluidics for digital quantification and single-cell sequencing.

We have developed hydrogel-based virtual microfluidics as a simple and robust alternative to complex engineered microfluidic systems for the compartmentalization of nucleic acid amplification reactions. We applied in-gel digital multiple displacement amplification (dMDA) to purified DNA templates, cultured bacterial cells and human microbiome samples in the virtual microfluidics system, and demonstrated whole-genome sequencing of single-cell MDA products with excellent coverage uniformity and markedly reduced chimerism compared with products of liquid MDA reactions.

Trends in Antibiotic Susceptibility in Staphylococcus aureus in Boston, Massachusetts, from 2000 to 2014.

The rate of infection by methicillin-resistant Staphylococcus aureus (MRSA) has declined over the past decade, but it is unclear whether this represents a decline in S. aureus infections overall. To evaluate the trends in the annual rates of infection by S. aureus subtypes and mean antibiotic resistance, we conducted a 15-year retrospective observational study at two tertiary care institutions in Boston, MA, of 31,753 adult inpatients with S. aureus isolated from clinical specimens. We inferred the gain and loss of methicillin resistance through genome sequencing of 180 isolates from 2016. The annual rates of infection by S. aureus declined from 2003 to 2014 by 4.2% (2.7% to 5.6%), attributable to an annual decline in MRSA of 10.9% (9.3% to 12.6%). Penicillin-susceptible S. aureus (PSSA) increased by 6.1% (4.2% to 8.1%) annually, and rates of methicillin-susceptible penicillin-resistant S. aureus (MSSA) did not change. Resistance in S. aureus decreased from 2000 to 2014 by 0.8 antibiotics (0.7 to 0.8). Within common MRSA clonal complexes, 3/14 MSSA and 2/21 PSSA isolates arose from the loss of resistance-conferring genes. Overall, in two tertiary care institutions in Boston, MA, a decline in S. aureus infections has been accompanied by a shift toward increased antibiotic susceptibility. The rise in PSSA makes penicillin an increasingly viable treatment option.

Molecular sled sequences are common in mammalian proteins.

Recent work revealed a new class of molecular machines called molecular sleds, which are small basic molecules that bind and slide along DNA with the ability to carry cargo along DNA. Here, we performed biochemical and single-molecule flow stretching assays to investigate the basis of sliding activity in molecular sleds. In particular, we identified the functional core of pVIc, the first molecular sled characterized; peptide functional groups that control sliding activity; and propose a model for the sliding activity of molecular sleds. We also observed widespread DNA binding and sliding activity among basic polypeptide sequences that implicate mammalian nuclear localization sequences and many cell penetrating peptides as molecular sleds. These basic protein motifs exhibit weak but physiologically relevant sequence-nonspecific DNA affinity. Our findings indicate that many mammalian proteins contain molecular sled sequences and suggest the possibility that substantial undiscovered sliding activity exists among nuclear mammalian proteins.