Large-scale mapping and mutagenesis of human transcriptional effector domains.

Human gene expression is regulated by more than 2,000 transcription factors and chromatin regulators1,2. Effector domains within these proteins can activate or repress transcription. However, for many of these regulators we do not know what type of effector domains they contain, their location in the protein, their activation and repression strengths, and the sequences that are necessary for their functions. Here, we systematically measure the effector activity of more than 100,000 protein fragments tiling across most chromatin regulators and transcription factors in human cells (2,047 proteins). By testing the effect they have when recruited at reporter genes, we annotate 374 activation domains and 715 repression domains, roughly 80% of which are new and have not been previously annotated3-5. Rational mutagenesis and deletion scans across all the effector domains reveal aromatic and/or leucine residues interspersed with acidic, proline, serine and/or glutamine residues are necessary for activation domain activity. Furthermore, most repression domain sequences contain sites for small ubiquitin-like modifier (SUMO)ylation, short interaction motifs for recruiting corepressors or are structured binding domains for recruiting other repressive proteins. We discover bifunctional domains that can both activate and repress, some of which dynamically split a cell population into high- and low-expression subpopulations. Our systematic annotation and characterization of effector domains provide a rich resource for understanding the function of human transcription factors and chromatin regulators, engineering compact tools for controlling gene expression and refining predictive models of effector domain function.

Decoupling of catalysis and transition state analog binding from mutations throughout a phosphatase revealed by high-throughput enzymology.

Using high-throughput microfluidic enzyme kinetics (HT-MEK), we measured over 9,000 inhibition curves detailing impacts of 1,004 single-site mutations throughout the alkaline phosphatase PafA on binding affinity for two transition state analogs (TSAs), vanadate and tungstate. As predicted by catalytic models invoking transition state complementary, mutations to active site and active-site-contacting residues had highly similar impacts on catalysis and TSA binding. Unexpectedly, most mutations to more distal residues that reduced catalysis had little or no impact on TSA binding and many even increased tungstate affinity. These disparate effects can be accounted for by a model in which distal mutations alter the enzyme's conformational landscape, increasing the occupancy of microstates that are catalytically less effective but better able to accommodate larger transition state analogs. In support of this ensemble model, glycine substitutions (rather than valine) were more likely to increase tungstate affinity (but not more likely to impact catalysis), presumably due to increased conformational flexibility that allows previously disfavored microstates to increase in occupancy. These results indicate that residues throughout an enzyme provide specificity for the transition state and discriminate against analogs that are larger only by tenths of an Ångström. Thus, engineering enzymes that rival the most powerful natural enzymes will likely require consideration of distal residues that shape the enzyme's conformational landscape and fine-tune active-site residues. Biologically, the evolution of extensive communication between the active site and remote residues to aid catalysis may have provided the foundation for allostery to make it a highly evolvable trait.

A bead-based method for high-throughput mapping of the sequence- and force-dependence of T cell activation.

Adaptive immunity relies on T lymphocytes that use αß T cell receptors (TCRs) to discriminate among peptides presented by major histocompatibility complex molecules (pMHCs). Identifying pMHCs capable of inducing robust T cell responses will not only enable a deeper understanding of the mechanisms governing immune responses but could also have broad applications in diagnosis and treatment. T cell recognition of sparse antigenic pMHCs in vivo relies on biomechanical forces. However, in vitro screening methods test potential pMHCs without force and often at high (nonphysiological) pMHC densities and thus fail to predict potent agonists in vivo. Here, we present a technology termed BATTLES (biomechanically assisted T cell triggering for large-scale exogenous-pMHC screening) that uses biomechanical force to initiate T cell triggering for peptides and cells in parallel. BATTLES displays candidate pMHCs on spectrally encoded beads composed of a thermo-responsive polymer capable of applying shear loads to T cells, facilitating exploration of the force- and sequence-dependent landscape of T cell responses. BATTLES can be used to explore basic T cell mechanobiology and T cell-based immunotherapies.

On the dependent recognition of some long zinc finger proteins.

The human genome contains about 800 C2H2 zinc finger proteins (ZFPs), and most of them are composed of long arrays of zinc fingers. Standard ZFP recognition model asserts longer finger arrays should recognize longer DNA-binding sites. However, recent experimental efforts to identify in vivo ZFP binding sites contradict this assumption, with many exhibiting short motifs. Here we use ZFY, CTCF, ZIM3, and ZNF343 as examples to address three closely related questions: What are the reasons that impede current motif discovery methods? What are the functions of those seemingly unused fingers and how can we improve the motif discovery algorithms based on long ZFPs' biophysical properties? Using ZFY, we employed a variety of methods and find evidence for 'dependent recognition' where downstream fingers can recognize some previously undiscovered motifs only in the presence of an intact core site. For CTCF, high-throughput measurements revealed its upstream specificity profile depends on the strength of its core. Moreover, the binding strength of the upstream site modulates CTCF's sensitivity to different epigenetic modifications within the core, providing new insight into how the previously identified intellectual disability-causing and cancer-related mutant R567W disrupts upstream recognition and deregulates the epigenetic control by CTCF. Our results establish that, because of irregular motif structures, variable spacing and dependent recognition between sub-motifs, the specificities of long ZFPs are significantly underestimated, so we developed an algorithm, ModeMap, to infer the motifs and recognition models of ZIM3 and ZNF343, which facilitates high-confidence identification of specific binding sites, including repeats-derived elements. With revised concept, technique, and algorithm, we can discover the overlooked specificities and functions of those 'extra' fingers, and therefore decipher their broader roles in human biology and diseases.

Phosphatase specificity principles uncovered by MRBLE:Dephos and global substrate identification.

Phosphoprotein phosphatases (PPPs) regulate major signaling pathways, but the determinants of phosphatase specificity are poorly understood. This is because methods to investigate this at scale are lacking. Here, we develop a novel in vitro assay, MRBLE:Dephos, that allows multiplexing of dephosphorylation reactions to determine phosphatase preferences. Using MRBLE:Dephos, we establish amino acid preferences of the residues surrounding the dephosphorylation site for PP1 and PP2A-B55, which reveals common and unique preferences. To compare the MRBLE:Dephos results to cellular substrates, we focused on mitotic exit that requires extensive dephosphorylation by PP1 and PP2A-B55. We use specific inhibition of PP1 and PP2A-B55 in mitotic exit lysates coupled with phosphoproteomics to identify more than 2,000 regulated sites. Importantly, the sites dephosphorylated during mitotic exit reveal key signatures that are consistent with MRBLE:Dephos. Furthermore, integration of our phosphoproteomic data with mitotic interactomes of PP1 and PP2A-B55 provides insight into how binding of phosphatases to substrates shapes dephosphorylation. Collectively, we develop novel approaches to investigate protein phosphatases that provide insight into mitotic exit regulation.

Alleviating Cell Lysate-Induced Inhibition to Enable RT-PCR from Single Cells in Picoliter-Volume Double Emulsion Droplets.

Microfluidic droplet assays enable single-cell polymerase chain reaction (PCR) and sequencing analyses at unprecedented scales, with most methods encapsulating cells within nanoliter-sized single emulsion droplets (water-in-oil). Encapsulating cells within picoliter double emulsion (DE) (water-in-oil-in-water) allows sorting droplets with commercially available fluorescence-activated cell sorter (FACS) machines, making it possible to isolate single cells based on phenotypes of interest for downstream analyses. However, sorting DE droplets with standard cytometers requires small droplets that can pass FACS nozzles. This poses challenges for molecular biology, as prior reports suggest that reverse transcription (RT) and PCR amplification cannot proceed efficiently at volumes below 1 nL due to cell lysate-induced inhibition. To overcome this limitation, we used a plate-based RT-PCR assay designed to mimic reactions in picoliter droplets to systematically quantify and ameliorate the inhibition. We find that RT-PCR is blocked by lysate-induced cleavage of nucleic acid probes and primers, which can be efficiently alleviated through heat lysis. We further show that the magnitude of inhibition depends on the cell type, but that RT-PCR can proceed in low-picoscale reaction volumes for most mouse and human cell lines tested. Finally, we demonstrate one-step RT-PCR from single cells in 20 pL DE droplets with fluorescence quantifiable via FACS. These results open up new avenues for improving picoscale droplet RT-PCR reactions and expanding microfluidic droplet-based single-cell analysis technologies.

How duplicated transcription regulators can diversify to govern the expression of nonoverlapping sets of genes.

The duplication of transcription regulators can elicit major regulatory network rearrangements over evolutionary timescales. However, few examples of duplications resulting in gene network expansions are understood in molecular detail. Here we show that four Candida albicans transcription regulators that arose by successive duplications have differentiated from one another by acquiring different intrinsic DNA-binding specificities, different preferences for half-site spacing, and different associations with cofactors. The combination of these three mechanisms resulted in each of the four regulators controlling a distinct set of target genes, which likely contributed to the adaption of this fungus to its human host. Our results illustrate how successive duplications and diversification of an ancestral transcription regulator can underlie major changes in an organism's regulatory circuitry.

DeCoDe: degenerate codon design for complete protein-coding DNA libraries.

MOTIVATION: High-throughput protein screening is a critical technique for dissecting and designing protein function. Libraries for these assays can be created through a number of means, including targeted or random mutagenesis of a template protein sequence or direct DNA synthesis. However, mutagenic library construction methods often yield vastly more nonfunctional than functional variants and, despite advances in large-scale DNA synthesis, individual synthesis of each desired DNA template is often prohibitively expensive. Consequently, many protein-screening libraries rely on the use of degenerate codons (DCs), mixtures of DNA bases incorporated at specific positions during DNA synthesis, to generate highly diverse protein-variant pools from only a few low-cost synthesis reactions. However, selecting DCs for sets of sequences that covary at multiple positions dramatically increases the difficulty of designing a DC library and leads to the creation of many undesired variants that can quickly outstrip screening capacity. RESULTS: We introduce a novel algorithm for total DC library optimization, degenerate codon design (DeCoDe), based on integer linear programming. DeCoDe significantly outperforms state-of-the-art DC optimization algorithms and scales well to more than a hundred proteins sharing complex patterns of covariation (e.g. the lab-derived avGFP lineage). Moreover, DeCoDe is, to our knowledge, the first DC design algorithm with the capability to encode mixed-length protein libraries. We anticipate DeCoDe to be broadly useful for a variety of library generation problems, ranging from protein engineering attempts that leverage mutual information to the reconstruction of ancestral protein states. AVAILABILITY AND IMPLEMENTATION: github.com/OrensteinLab/DeCoDe. CONTACT: yaronore@bgu.ac.il. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Diversification of DNA binding specificities enabled SREBP transcription regulators to expand the repertoire of cellular functions that they govern in fungi.

The Sterol Regulatory Element Binding Proteins (SREBPs) are basic-helix-loop-helix transcription regulators that control the expression of sterol biosynthesis genes in higher eukaryotes and some fungi. Surprisingly, SREBPs do not regulate sterol biosynthesis in the ascomycete yeasts (Saccharomycotina) as this role was handed off to an unrelated transcription regulator in this clade. The SREBPs, nonetheless, expanded in fungi such as the ascomycete yeasts Candida spp., raising questions about their role and evolution in these organisms. Here we report that the fungal SREBPs diversified their DNA binding preferences concomitantly with an expansion in function. We establish that several branches of fungal SREBPs preferentially bind non-palindromic DNA sequences, in contrast to the palindromic DNA motifs recognized by most basic-helix-loop-helix proteins (including SREBPs) in higher eukaryotes. Reconstruction and biochemical characterization of the likely ancestor protein suggest that an intrinsic DNA binding promiscuity in the family was resolved by alternative mechanisms in different branches of fungal SREBPs. Furthermore, we show that two SREBPs in the human commensal yeast Candida albicans drive a transcriptional cascade that inhibits a morphological switch under anaerobic conditions. Preventing this morphological transition enhances C. albicans colonization of the mammalian intestine, the fungus' natural niche. Thus, our results illustrate how diversification in DNA binding preferences enabled the functional expansion of a family of eukaryotic transcription regulators.

Comprehensive, high-resolution binding energy landscapes reveal context dependencies of transcription factor binding.

Transcription factors (TFs) are primary regulators of gene expression in cells, where they bind specific genomic target sites to control transcription. Quantitative measurements of TF-DNA binding energies can improve the accuracy of predictions of TF occupancy and downstream gene expression in vivo and shed light on how transcriptional networks are rewired throughout evolution. Here, we present a sequencing-based TF binding assay and analysis pipeline (BET-seq, for Binding Energy Topography by sequencing) capable of providing quantitative estimates of binding energies for more than one million DNA sequences in parallel at high energetic resolution. Using this platform, we measured the binding energies associated with all possible combinations of 10 nucleotides flanking the known consensus DNA target interacting with two model yeast TFs, Pho4 and Cbf1. A large fraction of these flanking mutations change overall binding energies by an amount equal to or greater than consensus site mutations, suggesting that current definitions of TF binding sites may be too restrictive. By systematically comparing estimates of binding energies output by deep neural networks (NNs) and biophysical models trained on these data, we establish that dinucleotide (DN) specificities are sufficient to explain essentially all variance in observed binding behavior, with Cbf1 binding exhibiting significantly more nonadditivity than Pho4. NN-derived binding energies agree with orthogonal biochemical measurements and reveal that dynamically occupied sites in vivo are both energetically and mutationally distant from the highest affinity sites.

Double Emulsion Picoreactors for High-Throughput Single-Cell Encapsulation and Phenotyping via FACS.

In the past five years, droplet microfluidic techniques have unlocked new opportunities for the high-throughput genome-wide analysis of single cells, transforming our understanding of cellular diversity and function. However, the field lacks an accessible method to screen and sort droplets based on cellular phenotype upstream of genetic analysis, particularly for large and complex cells. To meet this need, we developed Dropception, a robust, easy-to-use workflow for precise single-cell encapsulation into picoliter-scale double emulsion droplets compatible with high-throughput screening via fluorescence-activated cell sorting (FACS). We demonstrate the capabilities of this method by encapsulating five standardized mammalian cell lines of varying sizes and morphologies as well as a heterogeneous cell mixture of a whole dissociated flatworm (5-25 µm in diameter) within highly monodisperse double emulsions (35 µm in diameter). We optimize for preferential encapsulation of single cells with extremely low multiple-cell loading events (<2% of cell-containing droplets), thereby allowing direct linkage of cellular phenotype to genotype. Across all cell lines, cell loading efficiency approaches the theoretical limit with no observable bias by cell size. FACS measurements reveal the ability to discriminate empty droplets from those containing cells with good agreement to single-cell occupancies quantified via microscopy, establishing robust droplet screening at single-cell resolution. High-throughput FACS screening of cellular picoreactors has the potential to shift the landscape of single-cell droplet microfluidics by expanding the repertoire of current nucleic acid droplet assays to include functional phenotyping.

Discovering epistatic feature interactions from neural network models of regulatory DNA sequences.

Motivation: Transcription factors bind regulatory DNA sequences in a combinatorial manner to modulate gene expression. Deep neural networks (DNNs) can learn the cis-regulatory grammars encoded in regulatory DNA sequences associated with transcription factor binding and chromatin accessibility. Several feature attribution methods have been developed for estimating the predictive importance of individual features (nucleotides or motifs) in any input DNA sequence to its associated output prediction from a DNN model. However, these methods do not reveal higher-order feature interactions encoded by the models. Results: We present a new method called Deep Feature Interaction Maps (DFIM) to efficiently estimate interactions between all pairs of features in any input DNA sequence. DFIM accurately identifies ground truth motif interactions embedded in simulated regulatory DNA sequences. DFIM identifies synergistic interactions between GATA1 and TAL1 motifs from in vivo TF binding models. DFIM reveals epistatic interactions involving nucleotides flanking the core motif of the Cbf1 TF in yeast from in vitro TF binding models. We also apply DFIM to regulatory sequence models of in vivo chromatin accessibility to reveal interactions between regulatory genetic variants and proximal motifs of target TFs as validated by TF binding quantitative trait loci. Our approach makes significant strides in improving the interpretability of deep learning models for genomics. Availability and implementation: Code is available at: https://github.com/kundajelab/dfim. Supplementary information: Supplementary data are available at Bioinformatics online.

Identification and characterization of a previously undescribed family of sequence-specific DNA-binding domains.

Sequence-specific DNA-binding proteins are among the most important classes of gene regulatory proteins, controlling changes in transcription that underlie many aspects of biology. In this work, we identify a transcriptional regulator from the human fungal pathogen Candida albicans that binds DNA specifically but has no detectable homology with any previously described DNA- or RNA-binding protein. This protein, named White-Opaque Regulator 3 (Wor3), regulates white-opaque switching, the ability of C. albicans to switch between two heritable cell types. We demonstrate that ectopic overexpression of WOR3 results in mass conversion of white cells to opaque cells and that deletion of WOR3 affects the stability of opaque cells at physiological temperatures. Genome-wide chromatin immunoprecipitation of Wor3 and gene expression profiling of a wor3 deletion mutant strain indicate that Wor3 is highly integrated into the previously described circuit regulating white-opaque switching and that it controls a subset of the opaque transcriptional program. We show by biochemical, genetic, and microfluidic experiments that Wor3 binds directly to DNA in a sequence-specific manner, and we identify the set of cis-regulatory sequences recognized by Wor3. Bioinformatic analyses indicate that the Wor3 family arose more recently in evolutionary time than most previously described DNA-binding domains; it is restricted to a small number of fungi that include the major fungal pathogens of humans. These observations show that new families of sequence-specific DNA-binding proteins may be restricted to small clades and suggest that current annotations--which rely on deep conservation--underestimate the fraction of genes coding for transcriptional regulators.

Microfluidic affinity and ChIP-seq analyses converge on a conserved FOXP2-binding motif in chimp and human, which enables the detection of evolutionarily novel targets.

The transcription factor forkhead box P2 (FOXP2) is believed to be important in the evolution of human speech. A mutation in its DNA-binding domain causes severe speech impairment. Humans have acquired two coding changes relative to the conserved mammalian sequence. Despite intense interest in FOXP2, it has remained an open question whether the human protein's DNA-binding specificity and chromatin localization are conserved. Previous in vitro and ChIP-chip studies have provided conflicting consensus sequences for the FOXP2-binding site. Using MITOMI 2.0 microfluidic affinity assays, we describe the binding site of FOXP2 and its affinity profile in base-specific detail for all substitutions of the strongest binding site. We find that human and chimp FOXP2 have similar binding sites that are distinct from previously suggested consensus binding sites. Additionally, through analysis of FOXP2 ChIP-seq data from cultured neurons, we find strong overrepresentation of a motif that matches our in vitro results and identifies a set of genes with FOXP2 binding sites. The FOXP2-binding sites tend to be conserved, yet we identified 38 instances of evolutionarily novel sites in humans. Combined, these data present a comprehensive portrait of FOXP2's-binding properties and imply that although its sequence specificity has been conserved, some of its genomic binding sites are newly evolved.

Basic leucine zipper transcription factor Hac1 binds DNA in two distinct modes as revealed by microfluidic analyses.

A quantitative understanding of how transcription factors interact with genomic target sites is crucial for reconstructing transcriptional networks in vivo. Here, we use Hac1, a well-characterized basic leucine zipper (bZIP) transcription factor involved in the unfolded protein response (UPR) as a model to investigate interactions between bZIP transcription factors and their target sites. During the UPR, the accumulation of unfolded proteins leads to unconventional splicing and subsequent translation of HAC1 mRNA, followed by transcription of UPR target genes. Initial candidate-based approaches identified a canonical cis-acting unfolded protein response element (UPRE-1) within target gene promoters; however, subsequent studies identified a large set of Hac1 target genes lacking this UPRE-1 and containing a different motif (UPRE-2). Using a combination of unbiased and directed microfluidic DNA binding assays, we established that Hac1 binds in two distinct modes: (i) to short (6-7 bp) UPRE-2-like motifs and (ii) to significantly longer (11-13 bp) extended UPRE-1-like motifs. Using a genetic screen, we demonstrate that a region of extended homology N-terminal to the basic DNA binding domain is required for this dual site recognition. These results establish Hac1 as the first bZIP transcription factor known to adopt more than one binding mode and unify previously conflicting and discrepant observations of Hac1 function into a cohesive model of UPR target gene activation. Our results also suggest that even structurally simple transcription factors can recognize multiple divergent target sites of very different lengths, potentially enriching their downstream target repertoire.

Structure of the transcriptional network controlling white-opaque switching in Candida albicans.

The human fungal pathogen Candida albicans can switch between two phenotypic cell types, termed 'white' and 'opaque'. Both cell types are heritable for many generations, and the switch between the two types occurs epigenetically, that is, without a change in the primary DNA sequence of the genome. Previous work identified six key transcriptional regulators important for white-opaque switching: Wor1, Wor2, Wor3, Czf1, Efg1, and Ahr1. In this work, we describe the structure of the transcriptional network that specifies the white and opaque cell types and governs the ability to switch between them. In particular, we use a combination of genome-wide chromatin immunoprecipitation, gene expression profiling, and microfluidics-based DNA binding experiments to determine the direct and indirect regulatory interactions that form the switch network. The six regulators are arranged together in a complex, interlocking network with many seemingly redundant and overlapping connections. We propose that the structure (or topology) of this network is responsible for the epigenetic maintenance of the white and opaque states, the switching between them, and the specialized properties of each state.

Design automation of microfluidic single and double emulsion droplets with machine learning.

Droplet microfluidics enables kHz screening of picoliter samples at a fraction of the cost of other high-throughput approaches. However, generating stable droplets with desired characteristics typically requires labor-intensive empirical optimization of device designs and flow conditions that limit adoption to specialist labs. Here, we compile a comprehensive droplet dataset and use it to train machine learning models capable of accurately predicting device geometries and flow conditions required to generate stable aqueous-in-oil and oil-in-aqueous single and double emulsions from 15 to 250 µm at rates up to 12000 Hz for different fluids commonly used in life sciences. Blind predictions by our models for as-yet-unseen fluids, geometries, and device materials yield accurate results, establishing their generalizability. Finally, we generate an easy-to-use design automation tool that yield droplets within 3 µm (<8%) of the desired diameter, facilitating tailored droplet-based platforms and accelerating their utility in life sciences.

Individual dimers of the mitotic kinesin motor Eg5 step processively and support substantial loads in vitro.

Eg5, a member of the kinesin superfamily of microtubule-based motors, is essential for bipolar spindle assembly and maintenance during mitosis, yet little is known about the mechanisms by which it accomplishes these tasks. Here, we used an automated optical trapping apparatus in conjunction with a novel motility assay that employed chemically modified surfaces to probe the mechanochemistry of Eg5. Individual dimers, formed by a recombinant human construct Eg5-513-5His, stepped processively along microtubules in 8-nm increments, with short run lengths averaging approximately eight steps. By varying the applied load (with a force clamp) and the ATP concentration, we found that the velocity of Eg5 was slower and less sensitive to external load than that of conventional kinesin, possibly reflecting the distinct demands of spindle assembly as compared with vesicle transport. The Eg5-513-5His velocity data were described by a minimal, three-state model where a force-dependent transition follows nucleotide binding.

Versatility and stability optimization of flow-focusing droplet generators via quality metric-driven design automation.

Droplet generation is a fundamental component of droplet microfluidics, compartmentalizing biological or chemical systems within a water-in-oil emulsion. As adoption of droplet microfluidics expands beyond expert labs or integrated devices, quality metrics are needed to contextualize the performance capabilities, improving the reproducibility and efficiency of operation. Here, we present two quality metrics for droplet generation: performance versatility, the operating range of a single device, and stability, the distance of a single operating point from a regime change. Both metrics were characterized in silico and validated experimentally using machine learning and rapid prototyping. These metrics were integrated into a design automation workflow, DAFD 2.0, which provides users with droplet generators of a desired performance that are versatile or flow stable. Versatile droplet generators with stable operating points accelerate the development of sophisticated devices by facilitating integration of other microfluidic components and improving the accuracy of design automation tools.

Double emulsions as a high-throughput enrichment and isolation platform for slower-growing microbes.

Our understanding of in situ microbial physiology is primarily based on physiological characterization of fast-growing and readily-isolatable microbes. Microbial enrichments to obtain novel isolates with slower growth rates or physiologies adapted to low nutrient environments are plagued by intrinsic biases for fastest-growing species when using standard laboratory isolation protocols. New cultivation tools to minimize these biases and enrich for less well-studied taxa are needed. In this study, we developed a high-throughput bacterial enrichment platform based on single cell encapsulation and growth within double emulsions (GrowMiDE). We showed that GrowMiDE can cultivate many different microorganisms and enrich for underrepresented taxa that are never observed in traditional batch enrichments. For example, preventing dominance of the enrichment by fast-growing microbes due to nutrient privatization within the double emulsion droplets allowed cultivation of slower-growing Negativicutes and Methanobacteria from stool samples in rich media enrichment cultures. In competition experiments between growth rate and growth yield specialist strains, GrowMiDE enrichments prevented competition for shared nutrient pools and enriched for slower-growing but more efficient strains. Finally, we demonstrated the compatibility of GrowMiDE with commercial fluorescence-activated cell sorting (FACS) to obtain isolates from GrowMiDE enrichments. Together, GrowMiDE + DE-FACS is a promising new high-throughput enrichment platform that can be easily applied to diverse microbial enrichments or screens.