Reference-based analysis of lung single-cell sequencing reveals a transitional profibrotic macrophage.

Tissue fibrosis is a major cause of mortality that results from the deposition of matrix proteins by an activated mesenchyme. Macrophages accumulate in fibrosis, but the role of specific subgroups in supporting fibrogenesis has not been investigated in vivo. Here, we used single-cell RNA sequencing (scRNA-seq) to characterize the heterogeneity of macrophages in bleomycin-induced lung fibrosis in mice. A novel computational framework for the annotation of scRNA-seq by reference to bulk transcriptomes (SingleR) enabled the subclustering of macrophages and revealed a disease-associated subgroup with a transitional gene expression profile intermediate between monocyte-derived and alveolar macrophages. These CX3CR1+SiglecF+ transitional macrophages localized to the fibrotic niche and had a profibrotic effect in vivo. Human orthologs of genes expressed by the transitional macrophages were upregulated in samples from patients with idiopathic pulmonary fibrosis. Thus, we have identified a pathological subgroup of transitional macrophages that are required for the fibrotic response to injury.

Flow-Cytometric Analysis and Purification of Airway Epithelial-Cell Subsets.

The human airway epithelium is essential in homeostasis, and epithelial dysfunction contributes to chronic airway disease. Development of flow-cytometric methods to characterize subsets of airway epithelial cells will enable further dissection of airway epithelial biology. Leveraging single-cell RNA-sequencing data in combination with known cell type-specific markers, we developed panels of antibodies to characterize and isolate the major airway epithelial subsets (basal, ciliated, and secretory cells) from human bronchial epithelial-cell cultures. We also identified molecularly distinct subpopulations of secretory cells and demonstrated cell subset-specific expression of low-abundance transcripts and microRNAs that are challenging to analyze with current single-cell RNA-sequencing methods. These new tools will be valuable for analyzing and separating airway epithelial subsets and interrogating airway epithelial biology.

Accurate Bulk Quantitation of Droplet Digital Polymerase Chain Reaction.

Droplet digital PCR provides superior accuracy for nucleic acid quantitation. The requirement of microfluidics to generate and analyze the emulsions, however, is a barrier to its adoption, particularly in low resource settings or clinical laboratories. Here, we report a novel method to prepare ddPCR droplets by vortexing and readout of the results by bulk analysis of recovered amplicons. We demonstrate the approach by accurately quantitating SARS-CoV-2 sequences using entirely bulk processing and no microfluidics. Our approach for quantitating reactions should extend to all digital assays that generate amplicons, including digital PCR and LAMP conducted in droplets, microchambers, or nanoliter wells. More broadly, our approach combines important attributes of ddPCR, including enhanced accuracy and robustness to inhibition, with the high-volume sample processing ability of quantitative PCR.

A comparative analysis of single cell and droplet-based FACS for improving production phenotypes: Riboflavin overproduction in Yarrowia lipolytica.

Evolutionary approaches to strain engineering inherently require the identification of suitable selection techniques for the product and phenotype of interest. In this work, we undertake a comparative analysis of two related but functionally distinct methods of high-throughput screening: traditional single cell fluorescence activated cell sorting (single cell FACS) and microdroplet-enabled FACS (droplet FACS) using water/oil/water (w/o/w) emulsions. To do so, we first engineer and evolve the non-conventional yeast Yarrowia lipolytica for high extracellular production of riboflavin (vitamin B2), an innately fluorescent product. Following mutagenesis and adaptive evolution, a direct parity-matched comparison of these two selection strategies was conducted. Both single cell FACS and droplet FACS led to significant increases in total riboflavin titer (32 and 54 fold relative to the parental PO1f strain, respectively). However, single cell FACS favored intracellular riboflavin accumulation (with only 70% of total riboflavin secreted) compared with droplet FACS that favored extracellular product accumulation (with 90% of total riboflavin secreted). We find that for the test case of riboflavin, the extent of secretion and total production were highly correlated. The resulting differences in production modes and levels clearly demonstrate the significant impact that selection approaches can exert on final evolutionary outcomes in strain engineering. Moreover, we note that these results provide a cautionary tale when intracellular read-outs of product concentration (including signals from biosensors) are used as surrogates for total production of potentially secreted products. In this regard, these results demonstrate that extracellular production is best assayed through an encapsulation technique when performing high throughput screening.

An evolutionary metabolic engineering approach for enhancing lipogenesis in Yarrowia lipolytica.

Lipogenic organisms provide an ideal platform for biodiesel and oleochemical production. Through our previous rational metabolic engineering efforts, lipogenesis titers in Yarrowia lipolytica were significantly enhanced. However, the resulting strain still suffered from decreased biomass generation rates. Here, we employ a rapid evolutionary metabolic engineering approach linked with a floating cell enrichment process to improve lipogenesis rates, titers, and yields. Through this iterative process, we were able to ultimately improve yields from our prior strain by 55% to achieve production titers of 39.1g/L with upwards of 76% of the theoretical maximum yield of conversation. Isolated cells were saturated with up to 87% lipid content. An average specific productivity of 0.56g/L/h was achieved with a maximum instantaneous specific productivity of 0.89g/L/h during the lipid production phase in fermentation. Genomic sequencing of the evolved strains revealed a link between a decrease/loss of function mutation of succinate semialdehyde dehydrogenase, uga2, suggesting the importance of gamma-aminobutyric acid assimilation in lipogenesis. This linkage was validated through gene deletion experiments. This work presents an improved host strain that can serve as a platform for efficient oleochemical production.

Surveying the lipogenesis landscape in Yarrowia lipolytica through understanding the function of a Mga2p regulatory protein mutant.

Lipogenic organisms represent great starting points for metabolic engineering of oleochemical production. While previous engineering efforts were able to significantly improve lipid production in Yarrowia lipolytica, the lipogenesis landscape, especially with respect to regulatory elements, has not been fully explored. Through a comparative genomics and transcriptomics approach, we identified and validated a mutant mga2 protein that serves as a regulator of desaturase gene expression and potent lipogenesis factor. The resulting strain is enriched in unsaturated fatty acids. Comparing the underlying mechanism of this mutant to other previously engineered strains suggests that creating an imbalance between glycolysis and the TCA cycle can serve as a driving force for lipogenesis when combined with fatty acid catabolism overexpressions. Further comparative transcriptomics analysis revealed both distinct and convergent rewiring associated with these different genotypes. Finally, by combining metabolic engineering targets, it is possible to further engineer a strain containing the mutant mga2 gene to a lipid production titer of 25g/L.

Increasing expression level and copy number of a Yarrowia lipolytica plasmid through regulated centromere function.

Manipulating gene expression using different types of plasmids (centromeric, 2-micron, and autonomously replicating sequence) can expand expression levels and enable a quick characterization of combinations, both of which are highly desired for metabolic engineering applications. However, for the powerful industrial workhorse Yarrowia lipolytica, only a low-copy CEN plasmid is available. In this study, we engineered the CEN plasmid from Y. lipolytica by fusing different promoters upstream of the centromeric region to regulate its function and expand its range. By doing this, we successfully improved the copy number and gene expression at plasmid level by 80% and enabled a dynamic range of nearly 2.7-fold. This improvement was seen to be independent from both the promoter and gene used in the expression cassette. These results present a better starting point for more potent plasmids in Y. lipolytica.

G1/S restriction point coordinates phasic gene expression and cell differentiation.

Pluripotent embryonic stem cells have a unique cell cycle structure with a suppressed G1/S restriction point and little differential expression across the cell cycle phases. Here, we evaluate the link between G1/S restriction point activation, phasic gene expression, and cellular differentiation. Expression analysis reveals a gain in phasic gene expression across lineages between embryonic days E7.5 and E9.5. Genetic manipulation of the G1/S restriction point regulators miR-302 and P27 respectively accelerates or delays the onset of phasic gene expression in mouse embryos. Loss of miR-302-mediated p21 or p27 suppression expedites embryonic stem cell differentiation, while a constitutive Cyclin E mutant blocks it. Together, these findings uncover a causal relationship between emergence of the G1/S restriction point with a gain in phasic gene expression and cellular differentiation.

Droplet-microfluidics-assisted sequencing of HIV proviruses and their integration sites in cells from people on antiretroviral therapy.

The human immunodeficiency virus (HIV) integrates its genome into that of infected cells and may enter an inactive state of reversible latency that cannot be targeted using antiretroviral therapy. Sequencing such a provirus and the adjacent host junctions in individual cells may elucidate the mechanisms of the persistence of infected cells, but this is difficult owing to the 150-million-fold higher amount of background human DNA. Here we show that full-length proviruses connected to their contiguous HIV-host DNA junctions can be assembled via a high-throughput microfluidic assay where droplet-based whole-genome amplification of HIV DNA in its native context is followed by a polymerase chain reaction (PCR) to tag droplets containing proviruses for sequencing. We assayed infected cells from people with HIV receiving suppressive antiretroviral therapy, resulting in the detection and sequencing of paired proviral genomes and integration sites, 90% of which were not recovered by commonly used nested-PCR methods. The sequencing of individual proviral genomes with their integration sites could improve the genetic analysis of persistent HIV-infected cell reservoirs.

Tuning gene expression in Yarrowia lipolytica by a hybrid promoter approach.

The development of strong and tunable promoter elements is necessary to enable metabolic and pathway engineering applications for any host organism. Here, we have expanded and generalized a hybrid promoter approach to produce libraries of high-expressing, tunable promoters in the nonconventional yeast Yarrowia lipolytica. These synthetic promoters are comprised of two modular components: the enhancer element and the core promoter element. By exploiting this basic promoter architecture, we have overcome native expression limitations and provided a strategy for both increasing the native promoter capacity and producing libraries for tunable gene expression in a cellular system with ill-defined genetic tools. In doing so, this work has created the strongest promoters ever reported for Y. lipolytica. Furthermore, we have characterized these promoters at the single-cell level through the use of a developed fluorescence-based assay as well as at the transcriptional and whole-cell levels. The resulting promoter libraries exhibited a range of more than 400-fold in terms of mRNA levels, and the strongest promoters in this set had 8-fold-higher fluorescence levels than those of typically used endogenous promoters. These results suggest that promoters in Y. lipolytica are enhancer limited and that this limitation can be partially or fully alleviated through the addition of tandem copies of upstream activation sequences (UASs). Finally, this work illustrates that tandem copies of UAS regions can serve as synthetic transcriptional amplifiers that may be generically used to increase the expression levels of promoters.

Accurate bulk quantitation of droplet digital PCR.

Droplet digital PCR provides superior accuracy in nucleic acid quantitation. The requirement of microfluidics to generate and analyze the emulsions, however, is a barrier to its adoption, particularly in low resource or clinical settings. Here, we report a novel method to prepare ddPCR droplets by vortexing and readout the results by bulk analysis of recovered amplicons. We demonstrate the approach by accurately quantitating SARS-CoV-2 sequences using entirely bulk processing and no microfluidics. Our approach for quantitating reactions should extend to all digital assays that generate amplicons, including digital PCR and LAMP conducted in droplets, microchambers, or nanoliter wells. More broadly, our approach combines important attributes of ddPCR, including enhanced accuracy and robustness to inhibition, with the high-volume sample processing ability of quantitative PCR.

Digital droplet PCR accurately quantifies SARS-CoV-2 viral load from crude lysate without nucleic acid purification.

The COVID-19 pandemic caused by the SARS-CoV-2 virus motivates diverse diagnostic approaches due to the novel causative pathogen, incompletely understood clinical sequelae, and limited availability of testing resources. Given the variability in viral load across and within patients, absolute viral load quantification directly from crude lysate is important for diagnosis and surveillance. Here, we investigate the use of digital droplet PCR (ddPCR) for SARS-CoV-2 viral load measurement directly from crude lysate without nucleic acid purification. We demonstrate ddPCR accurately quantifies SARS-CoV-2 standards from purified RNA and multiple sample matrices, including commonly utilized universal transport medium (UTM). In addition, we find ddPCR functions robustly at low input viral copy numbers on nasopharyngeal swab specimens stored in UTM without upfront RNA extraction. We also show ddPCR, but not qPCR, from crude lysate shows high concordance with viral load measurements from purified RNA. Our data suggest ddPCR offers advantages to qPCR for SARS-CoV-2 detection with higher sensitivity and robustness when using crude lysate rather than purified RNA as input. More broadly, digital droplet assays provide a potential method for nucleic acid measurement and infectious disease diagnosis with limited sample processing, underscoring the utility of such techniques in laboratory medicine.

Mapping enzyme catalysis with metabolic biosensing.

Enzymes are represented across a vast space of protein sequences and structural forms and have activities that far exceed the best chemical catalysts; however, engineering them to have novel or enhanced activity is limited by technologies for sensing product formation. Here, we describe a general and scalable approach for characterizing enzyme activity that uses the metabolism of the host cell as a biosensor by which to infer product formation. Since different products consume different molecules in their synthesis, they perturb host metabolism in unique ways that can be measured by mass spectrometry. This provides a general way by which to sense product formation, to discover unexpected products and map the effects of mutagenesis.

High Throughput Yeast Strain Phenotyping with Droplet-Based RNA Sequencing.

The powerful tools available to edit yeast genomes have made this microbe a valuable platform for engineering. While it is now possible to construct libraries of millions of genetically distinct strains, screening for a desired phenotype remains a significant obstacle. With existing screening techniques, there is a tradeoff between information output and throughput, with high-throughput screening typically being performed on one product of interest. Therefore, we present an approach to accelerate strain screening by adapting single cell RNA sequencing to isogenic picoliter colonies of genetically engineered yeast strains. To address the unique challenges of performing RNA sequencing on yeast cells, we culture isogenic yeast colonies within hydrogels and spheroplast prior to performing RNA sequencing. The RNA sequencing data can be used to infer yeast phenotypes and sort out engineered pathways. The scalability of our method addresses a critical obstruction in microbial engineering.

Digital droplet PCR accurately quantifies SARS-CoV-2 viral load from crude lysate without nucleic acid purification.

The COVID-19 pandemic caused by the SARS-CoV-2 virus motivates diverse diagnostic approaches due to the novel causative pathogen, incompletely understood clinical sequelae, and limited availability of testing resources. Given the variability in viral load across and within patients, absolute viral load quantification directly from crude lysate is important for diagnosis and surveillance. Here, we investigate the use of digital droplet PCR (ddPCR) for SARS-CoV-2 viral load measurement directly from crude lysate without nucleic acid purification. We demonstrate ddPCR accurately quantifies SARS-CoV-2 standards from purified RNA and multiple sample matrices, including commonly utilized universal transport medium (UTM). In addition, we find ddPCR functions robustly at low input viral copy numbers on nasopharyngeal swab specimens stored in UTM without upfront RNA extraction. We also show ddPCR, but not qPCR, from crude lysate shows high concordance with viral load measurements from purified RNA. Our data suggest ddPCR offers advantages to qPCR for SARS-CoV-2 detection with higher sensitivity and robustness when using crude lysate rather than purified RNA as input. More broadly, digital droplet assays provide a potential method for nucleic acid measurement and infectious disease diagnosis with limited sample processing, underscoring the utility of such techniques in laboratory medicine.

Linked optical and gene expression profiling of single cells at high-throughput.

Single-cell RNA sequencing has emerged as a powerful tool for characterizing cells, but not all phenotypes of interest can be observed through changes in gene expression. Linking sequencing with optical analysis has provided insight into the molecular basis of cellular function, but current approaches have limited throughput. Here, we present a high-throughput platform for linked optical and gene expression profiling of single cells. We demonstrate accurate fluorescence and gene expression measurements on thousands of cells in a single experiment. We use the platform to characterize DNA and RNA changes through the cell cycle and correlate antibody fluorescence with gene expression. The platform's ability to isolate rare cell subsets and perform multiple measurements, including fluorescence and sequencing-based analysis, holds potential for scalable multi-modal single-cell analysis.

High throughput gene expression profiling of yeast colonies with microgel-culture Drop-seq.

Yeast can be engineered into "living foundries" for non-natural chemical production by reprogramming them via a "design-build-test" cycle. While methods for "design" and "build" are relatively scalable and efficient, "test" remains a bottleneck, limiting the effectiveness of the procedure. Here we describe isogenic colony sequencing (ICO-seq), a massively-parallel strategy to assess the gene expression, and thus engineered pathway efficacy, of large numbers of genetically distinct yeast colonies. We use the approach to characterize opaque-white switching in 658 C. albicans colonies. By profiling the transcriptomes of 1642 engineered S. cerevisiae strains, we assess gene expression heterogeneity in a protein mutagenesis library. Our approach will accelerate synthetic biology by allowing facile and cost-effective transcriptional profiling of large numbers of genetically distinct yeast strains.

Combined aptamer and transcriptome sequencing of single cells.

The transcriptome and proteome encode distinct information that is important for characterizing heterogeneous biological systems. We demonstrate a method to simultaneously characterize the transcriptomes and proteomes of single cells at high throughput using aptamer probes and droplet-based single cell sequencing. With our method, we differentiate distinct cell types based on aptamer surface binding and gene expression patterns. Aptamers provide advantages over antibodies for single cell protein characterization, including rapid, in vitro, and high-purity generation via SELEX, and the ability to amplify and detect them with PCR and sequencing.

RNA-aptamers-in-droplets (RAPID) high-throughput screening for secretory phenotypes.

Synthetic biology and metabolic engineering seek to re-engineer microbes into "living foundries" for the production of high value chemicals. Through a "design-build-test" cycle paradigm, massive libraries of genetically engineered microbes can be constructed and tested for metabolite overproduction and secretion. However, library generation capacity outpaces the rate of high-throughput testing and screening. Well plate assays are flexible but with limited throughput, whereas droplet microfluidic techniques are ultrahigh-throughput but require a custom assay for each target. Here we present RNA-aptamers-in-droplets (RAPID), a method that greatly expands the generality of ultrahigh-throughput microfluidic screening. Using aptamers, we transduce extracellular product titer into fluorescence, allowing ultrahigh-throughput screening of millions of variants. We demonstrate the RAPID approach by enhancing production of tyrosine and secretion of a recombinant protein in Saccharomyces cerevisiae by up to 28- and 3-fold, respectively. Aptamers-in-droplets affords a general approach for evolving microbes to synthesize and secrete value-added chemicals.Screening libraries of genetically engineered microbes for secreted products is limited by the available assay throughput. Here the authors combine aptamer-based fluorescent detection with droplet microfluidics to achieve high throughput screening of yeast strains engineered for enhanced tyrosine or streptavidin production.

Draft Genome Sequence of the Oleaginous Yeast Yarrowia lipolytica PO1f, a Commonly Used Metabolic Engineering Host.

The draft genome sequence of the oleaginous yeast Yarrowia lipolytica stain PO1f, a commonly used metabolic engineering host, is presented here. The approximately 20.3-Mb genome sequence of PO1f will greatly facilitate research efforts in metabolic engineering of Yarrowia lipolytica for value-added chemical production.