De novo mutations mediate phenotypic switching in an opportunistic human lung pathogen.

Bacteria evolving within human hosts encounter selective tradeoffs that render mutations adaptive in one context and deleterious in another. Here, we report that the cystic fibrosis-associated pathogen Burkholderia dolosa overcomes in-human selective tradeoffs by acquiring successive point mutations that alternate phenotypes. We sequenced the whole genomes of 931 respiratory isolates from two recently infected patients and an epidemiologically-linked, chronically-infected patient. These isolates are contextualized using 112 historical genomes from the same outbreak strain. Within both newly infected patients, diverse parallel mutations that disrupt O-antigen expression quickly arose, comprising 29% and 63% of their B. dolosa communities by 3 years. The selection for loss of O-antigen starkly contrasts with our previous observation of parallel O-antigen-restoring mutations after many years of chronic infection in the historical outbreak. Experimental characterization revealed that O-antigen loss increases uptake in immune cells while decreasing competitiveness in the mouse lung. We propose that the balance of these pressures, and thus whether O-antigen expression is advantageous, depends on tissue localization and infection duration. These results suggest that mutation-driven alternation during infection may be more frequent than appreciated and is underestimated without dense temporal sampling.

High-throughput RNA isoform sequencing using programmed cDNA concatenation.

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.

Compressed phenotypic screens for complex multicellular models and high-content assays.

High-throughput phenotypic screens leveraging biochemical perturbations, high-content readouts, and complex multicellular models could advance therapeutic discovery yet remain constrained by limitations of scale. To address this, we establish a method for compressing screens by pooling perturbations followed by computational deconvolution. Conducting controlled benchmarks with a highly bioactive small molecule library and a high-content imaging readout, we demonstrate increased efficiency for compressed experimental designs compared to conventional approaches. To prove generalizability, we apply compressed screening to examine transcriptional responses of patient-derived pancreatic cancer organoids to a library of tumor-microenvironment (TME)-nominated recombinant protein ligands. Using single-cell RNA-seq as a readout, we uncover reproducible phenotypic shifts induced by ligands that correlate with clinical features in larger datasets and are distinct from reference signatures available in public databases. In sum, our approach enables phenotypic screens that interrogate complex multicellular models with rich phenotypic readouts to advance translatable drug discovery as well as basic biology.

The future of antibiotics begins with discovering new combinations.

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.

Plasma from patients with bacterial sepsis or severe COVID-19 induces suppressive myeloid cell production from hematopoietic progenitors in vitro.

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.

A validated lineage-derived somatic truth data set enables benchmarking in cancer genome analysis.

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.

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.

Single-Cell Transcriptional Profiling of Cells Derived From Regenerating Alveolar Ducts.

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.

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.

An immune-cell signature of bacterial sepsis.

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.

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.

Multiplexed enrichment and genomic profiling of peripheral blood cells reveal subset-specific immune signatures.

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.

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.

Sliding on DNA: From Peptides to Small Molecules.

Many DNA binding proteins utilize one-dimensional (1D) diffusion along DNA to accelerate their DNA target recognition. Although 1D diffusion of proteins along DNA has been studied for decades, a quantitative understanding is only beginning to emerge and few chemical tools are available to apply 1D diffusion as a design principle. Recently, we discovered that peptides can bind and slide along DNA-even transporting cargo along DNA. Such molecules are known as molecular sleds. Here, to advance our understanding of structure-function relationships governing sequence nonspecific DNA interaction of natural molecular sleds and to explore the potential for controlling sliding activity, we test the DNA binding and sliding activities of chemically modified peptides and analogs, and show that synthetic small molecules can slide on DNA. We found new ways to control molecular sled activity, novel small-molecule synthetic sleds, and molecular sled activity in N-methylpyrrole/N-methylimidazole polyamides that helps explain how these molecules locate rare target sites.

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.

Molecular sled is an eleven-amino acid vehicle facilitating biochemical interactions via sliding components along DNA.

Recently, we showed the adenovirus proteinase interacts productively with its protein substrates in vitro and in vivo in nascent virus particles via one-dimensional diffusion along the viral DNA. The mechanism by which this occurs has heretofore been unknown. We show sliding of these proteins along DNA occurs on a new vehicle in molecular biology, a 'molecular sled' named pVIc. This 11-amino acid viral peptide binds to DNA independent of sequence. pVIc slides on DNA, exhibiting the fastest one-dimensional diffusion constant, 26±1.8 × 10(6) (bp)(2) s(-1). pVIc is a 'molecular sled,' because it can slide heterologous cargos along DNA, for example, a streptavidin tetramer. Similar peptides, for example, from the C terminus of ß-actin or NLSIII of the p53 protein, slide along DNA. Characteristics of the 'molecular sled' in its milieu (virion, nucleus) have implications for how proteins in the nucleus of cells interact and imply a new form of biochemistry, one-dimensional biochemistry.

A microfluidic platform enabling single-cell RNA-seq of multigenerational lineages.

We introduce a microfluidic platform that enables off-chip single-cell RNA-seq after multi-generational lineage tracking under controlled culture conditions. We use this platform to generate whole-transcriptome profiles of primary, activated murine CD8+ T-cell and lymphocytic leukemia cell line lineages. Here we report that both cell types have greater intra- than inter-lineage transcriptional similarity. For CD8+ T-cells, genes with functional annotation relating to lymphocyte differentiation and function--including Granzyme B--are enriched among the genes that demonstrate greater intra-lineage expression level similarity. Analysis of gene expression covariance with matched measurements of time since division reveals cell type-specific transcriptional signatures that correspond with cell cycle progression. We believe that the ability to directly measure the effects of lineage and cell cycle-dependent transcriptional profiles of single cells will be broadly useful to fields where heterogeneous populations of cells display distinct clonal trajectories, including immunology, cancer, and developmental biology.

Laser microdissection of the alveolar duct enables single-cell genomic analysis.

Complex tissues such as the lung are composed of structural hierarchies such as alveoli, alveolar ducts, and lobules. Some structural units, such as the alveolar duct, appear to participate in tissue repair as well as the development of bronchioalveolar carcinoma. Here, we demonstrate an approach to conduct laser microdissection of the lung alveolar duct for single-cell PCR analysis. Our approach involved three steps. (1) The initial preparation used mechanical sectioning of the lung tissue with sufficient thickness to encompass the structure of interest. In the case of the alveolar duct, the precision-cut lung slices were 200 µm thick; the slices were processed using near-physiologic conditions to preserve the state of viable cells. (2) The lung slices were examined by transmission light microscopy to target the alveolar duct. The air-filled lung was sufficiently accessible by light microscopy that counterstains or fluorescent labels were unnecessary to identify the alveolar duct. (3) The enzymatic and microfluidic isolation of single cells allowed for the harvest of as few as several thousand cells for PCR analysis. Microfluidics based arrays were used to measure the expression of selected marker genes in individual cells to characterize different cell populations. Preliminary work suggests the unique value of this approach to understand the intra- and intercellular interactions within the regenerating alveolar duct.

Regulation of a viral proteinase by a peptide and DNA in one-dimensional space: I. binding to DNA AND to hexon of the precursor to protein VI, pVI, of human adenovirus.

The precursor to adenovirus protein VI, pVI, is a multifunctional protein with different roles early and late in virus infection. Here, we focus on two roles late in infection, binding of pVI to DNA and to the major capsid protein hexon. pVI bound to DNA as a monomer independent of DNA sequence with an apparent equilibrium dissociation constant, K(d)((app)), of 46 nm. Bound to double-stranded DNA, one molecule of pVI occluded 8 bp. Upon the binding of pVI to DNA, three sodium ions were displaced from the DNA. A ΔG(0)(0) of -4.54 kcal/mol for the nonelectrostatic free energy of binding indicated that a substantial component of the binding free energy resulted from nonspecific interactions between pVI and DNA. The proteolytically processed, mature form of pVI, protein VI, also bound to DNA; its K(d)((app)) was much higher, 307 nm. The binding assays were performed in 1 mm MgCl(2) because in the absence of magnesium, the binding to pVI or protein VI to DNA was too tight to determine a K(d)((app)). Three molecules of pVI bound to one molecule of the hexon trimer with an equilibrium dissociation constant K(d)((app)) of 1.1 nm.