Strain-level epidemiology of microbial communities and the human microbiome

1The importance of phenotypes and physiology characteristic of specific microbial strains has been recognized as early as the nineteenth century. Robert Koch’s postulates, for example, differentiate between disease-causing “pathogens” and benign but closely related microbial variants. While the surprising differences between otherwise similar microbial strains has thus been critical in infectious disease management and microbiology for centuries, it has only recently become accessible in the context of microbial communities and their ecology. It remains technically challenging to detect and differentiate among closely related microbial strains within communities, and we will discuss several high-throughput culture-independent and culture-based methods for doing so here. More importantly, though, the beginning of such work has shown strain variation in the human microbiome to be as important in the structure, function, immunology, and epidemiology of our “normal” microbial residents as it is in the definition of pathogenicity (Box 1).

2Particularly within communities that are by definition collections of heterogeneous cells, it has proven to be technically challenging to detect and differentiate among cells containing such closely related but highly variable genomes. Indeed, it is not yet clear how clonally most microbial lineages remain within typical in vivo communities. This suggests both basic questions about the generation and maintenance of closely related genome variants in any microbial community, and also pressing translational questions regarding the personalization and health consequences of strains in the human microbiome. Because of the extensive genetic and genomic (i.e., functional) differences between even closely related microbial strains, work to date has only rarely been powered to associate “commensal” microbial strains with their health consequences. Here, we thus review the ecology and effects known to date for microbial strain variants carried within the human microbiome, quantitative methods for their detection and epidemiology, and potential next steps including characterization of their surprisingly large pangenomic content of biochemical dark matter.

3Culture-based comparative genetics of isolates has been a mainstay of microbial characterization for decades, and along with culture-independent techniques, it is increasingly important in an era of high-throughput “culturomics” and creative isolation methods. Especially for human pathogens that are both of clinical interest and relatively easily culturable, hundreds or thousands of genomes have been used in some cases to compare strains and their transmission, associate SNV and structural variation to microbial or host phenotype, and define the genetic and evolutionary architectures of species and other clades. Metagenomic methods have the unique ability to extend these strain-specific investigations to almost any environment or microbe, while leveraging the insights already built up using isolate genomics. In particular, if a “strain” is considered to be a clonal genotype, it must correspond to a specific set of genes and resulting functionality. This functional perspective on strains has captured a wide range of operational architectures, since some processes are well-conserved across entire clades (e.g., butyrate production in Faecalibacterium prausnitzii). Others, conversely, are highly variable even within specific benign or pathogenic species—Escherichia coli in the gut being the most prominent example.

4The gut is the greatest reservoir of biomass in the human microbiome, the body’s largest immune exposure, the most well-studied contributor to microbiome-linked disease, and one of the most ecologically diverse human-associated microbial habitats. It is also the source of several of the most canonical examples of radically different microbial physiology among closely related strains, such as the benign E. coli variants carried in most guts as compared to acute pathogens such as enterohemorrhagic E. coli (EHEC) O157:H7, long-term risks such as colorectal cancer in association with colibactin production in pks + E. coli, or the probiotic E. coli Nissle 1917. Isolate cultures have identified other strain-specific characteristics associated with evolutionary advantages ranging from increased virulence, mobility, nutrient acquisition, antibiotic resistance, and defense.

5Strains abundant in the infant gut are only rarely abundant in maternal microbiomes and are often replaced within the first 1–2 years of life. Their similarity to maternal, familial, or generally environmental strains is also itself highly variable and species-specific, but even small structural variants may be crucial in immune programming during temporally specific developmental windows. Like developmental variants of human gene products, such as hemoglobin forms, this dynamism in early life has functional consequences: Bifidobacterium longum, for example, is selected for human milk oligosaccharide (HMO) utilization in breastfeeding infants, whereas closely related B. longum strains in the adult gut frequently possess the capacity to ferment carbohydrates, but not HMOs. Strains abundant in the infant gut are only rarely abundant in maternal microbiomes and are often replaced within the first 1–2 years of life, but even small structural variants may be crucial in immune programming during temporally specific developmental windows. Ultimately, microbial strain variants affect not only host and individual microbes’ physiology, but also the ecology and phylogenetics of the overall gut community: Helicobacter pylori is one of the best-known examples of resident microbial genetic variation paralleling that of human host populations, but this has recently been shown to be the case for multiple subsets of the gut microbiome, such as Prevotella copri or Eubacterium rectale. This leads to linkages between the evolution and diversification of gut microbial community strains and host migration, geography, and lifestyle.

6One of the most crucial environmental factors related to this in the gut is diet, both acutely and over evolutionary time scales. However, the specifics of this relationship have been difficult to tease apart in human populations, due to the challenges of measuring diverse human diets, the confounding of long-term diet with other environmental factors, and the complexity of diet-microbial biochemical interactions. Indeed, diet represents only one aspect of gut microbial interaction with our biochemical environment, with several examples identified to date of strain-specific metabolism of drugs such as digoxin, metformin, acetaminophen, and potentially many others. With respect to diet itself, De Filippis et al., for example, found a greater abundance of P. copri among participants more closely adhering to a Mediterranean-style diet enriched with olive oil, fish, fruits, and vegetables. In contrast, Kovatcheva-Datchary et al. observed that even on the same barley-rich diet, Prevotella was only enriched among select participants, potentially in a strain-specific manner. De Filippis et al. later found similar heterogeneity among individuals on low-fat diets. Other examples include strains of short-chain fatty acid (SCFA)-producing bacteria with differential responses to fiber-enriched diets. Perhaps one of the most extreme examples of diet-linked strain specificity in the gut are among probiotic organisms such as Lactobacillus and Bifidobacterium, for which strains characteristic of fermented foods are highly distinct from those more typically resident in the human gut. The health consequences of probiotics can also be strain-specific dependent either on the strain context of the microbiome being entered, or on the strain of the probiotic organisms, e.g., the recently proposed ability of some bifidobacteria to facilitate cancer immunotherapy.

7While many studies have linked overall microbiome structure or microbial species enrichments to gastrointestinal (GI) or systemic disease, relatively few have identified strain-specific microbial variants associated with these diseases. The inflammatory bowel diseases (IBD) are among the best-studied chronic gastrointestinal conditions with respect to the microbiome, and in IBD, subspecies of E. coli and Ruminococcus gnavus have each been associated with disease severity. Hall et al. noted a particular subpopulation of R. gnavus strains more abundant in the IBD gut, enriched for adaptations to oxidative stress response, adhesion, and the utilization of iron and mucus. Bacteroides fragilis strains exhibit divergent behaviors leading to differential IgA induction in mouse models of IBD and have been associated with host immunomodulatory effects in monocolonization. While there are decades of work demonstrating the effects of such variants during animal monocolonization, understanding their effects in the human gut remains challenging, since the equivalent of a human genome-wide association study for most microbial community genetic variants (i.e., those not of very high penetrance) would be challenging, given the degree of multiple hypothesis testing necessary to account for the underlying microbial genetic variability.

8Studies of systemic disease outside of the gastrointestinal tract have also suggested functional roles for specific gut microbial strains. New-onset rheumatoid arthritis patients appear to be enriched for P. copri in the gut in some populations, for example, with evidence that this P. copri subset may be functionally or phylogenetically distinct. Obesity and type 2 diabetes (T2D) have shown relatively weak taxonomic or functional shifts in the gut microbiome overall, but again using mice to avoid challenges in human population structure, specific strains of Akkermansia muciniphila proved to be causal in alleviating these metabolic conditions. In human subjects, at least one study found SNPs specific to Bacteroides coprocola subpopulations within a T2D patient group. More broadly, strain-specific promotion of several SCFA producers, including Bifidobacterium spp., Eubacterium spp., and Lactobacillus spp., was selectively enriched by dietary fiber in a randomized clinical trial, improving T2D parameters.

9One of the most complex conditions bridging the gut microbiome, gastrointestinal, and systemic health has proven to be cancer. Particularly in colorectal cancer (CRC), specific microbial strain functionality can be readily shown to be locally causal, such as DNA-damaging production of colibactin by pks + E. coli as introduced above or B. fragilis toxin. Other microbes such as CRC-specific lineages of Fusobacterium nucleatum have been identified more recently, with mechanisms such as Fap2-mediated binding to host Gal-GalNAc or immunomodulation via TIGIT mediating both their carcinogenicity and their differentiation from typical oral F. nucleatum strains. Other mechanisms of microbial influence on GI or systemic cancer remain less well-understood, with strong evidence of resident microbial effects on immunotherapy responsiveness, but as yet few strain-specific culprits. Likewise, limited studies have shown intratumoral bacteria within and outside of the colon to be capable of direct metabolism of chemotherapeutics such as gemcitabine, with potentially many more such microbe-chemical interactions waiting to be discovered.

10While the strain epidemiology of the gut microbiome is perhaps best developed, similar examples exist of the effects of “commensal” and pathogenic strains throughout the human body habitat. As with the gut, the most extreme examples are those of well-studied pathogens, such as resistant variants of Staphylococcus aureus in the skin and nasal microbiomes. More recently, combinations of culture-independent and high-throughput culture-based methods have exposed within-subject pathogen evolution over the course of months to years. In these cases, as with pks + E. coli, resistance functionality such as mecA can be attributed to just one or a few loci that are genetically variable among strains via mobile chromosomal or plasmid-encoded elements. More unexpectedly, however, recent findings have pointed to correspondingly strain-specific interactions with non-pathogenic commensals, such as coporphyrin III production by some Cutibacterium (formerly Propionibacterium) strains inducing Staphylococcus biofilm formation. Indeed, due to their biogeographical heterogeneity relative to the gut, exposed topographical surfaces such as the skin, nasopharynx, and lung are among the few body areas where detailed ecology and persistence of multiple competing strains within an individual has been directly observed, e.g., among S. epidermidis strains in psoriasis.

11Conversely, deep differentiation of strains within an individual is technically more challenging in the vaginal microbiome. Instead, this environment has revealed extensive subspecies heterogeneity between hosts within the dominant Lactobacillus and other species of the vagina, again raising issues regarding the exact definition of strains and species among different microbial clades. Specifically, analysis of the intraspecific diversity of vaginally dominant lactobacilli such as L. jensenii, L. iners, L. gasserii, and L. crispatus is complicated by the systematics of the clade, which has been under scrutiny for reorganization based on both isolate and culture-independent genomics. Nevertheless, vaginal Lactobacillus and other strains can be reasonably stable within individuals over time, with particularly large environmental changes such as pregnancy inducing shifts over the course of gestation. As in the gut, such genetic variation between strains can affect health, such as in the determinants of pathogenicity in E. coli causing urinary tract infections. In examples from even more acute infectious disease, strain-specific Lactobacillus bioactivity can itself contribute to risk of sexually transmitted infection acquisition such as HIV, both due to direct microbial biochemistry and its effect on host immunity.

12Finally, oral microbiology has historically provided some of the first and most striking examples of phenotypic heterogeneity between closely related microbial isolates, and this trend holds true in the era of culture-independent sequencing and whole-community studies as well. Indeed, some of the earliest large population-scale surveys of the microbiome found oral site tropism to be a strong driver of subspecies differentiation, with stable genetic differences among related microbial colonizers of different surfaces—including different teeth—within the same mouth. These potentially adaptive, highly niche-specific variants have begun to be explored at scale, remaining stable within individual up to hundreds of days within subjects, but revealing extensive long-term plasticity between members of clades such as the Neisseria. While there is extensive ongoing work regarding the role of overall oral microbial ecology in conditions from periodontitis to pancreatic cancer and heart disease, the ecological and genomic diversity of the oral microbiota has led to limited strain-specific associations to date. Several have been suggested for, e.g., Streptococcus variants in caries or F. nucleatum in association with oral cancer—suggesting intriguing links with its role in CRC. These include sufficient detail to implicate microbial processes such as polyamine biosynthesis, motility and chemotaxis, and immunostimulation (e.g., LPS and flagellar components), but without yet a clear picture of the many possible strains across which these functions may be distributed in the complex oral environment.

13It is not our goal here to summarize the many methods that have been used to differentiate among microbial strains in culture over decades of microbiology, so we will focus in this review mainly on culture-independent techniques, as well as some high-throughput culture-based methods appropriate for microbial communities (Fig. 1). In both of these categories, many strain definition methods rely on sequencing: assembly of culture-based isolates, or amplicon-based, shotgun metagenomic, or single-cell culture-independent approaches. Other molecular assays, particularly mass spectrometry (MS)-based proteomics, can be applied to strain-type either isolates or communities. This is also true for MS- or NMR-based metabolomics or metabolic flux measurements. Of course, microbial culture physiology and direct imaging has been used to differentiate among strains since the earliest microbiology, and in some cases, these time-tested methods can be applied to communities as well.

14The first breakthroughs in microbial strain identification from whole-community sequencing—like the first community-wide applications of sequencing generally—came from marker gene approaches relying on amplification of 16S rRNA gene variable regions (amplicon or “16S” sequencing, Table 1). In many cases, amplicon-based technologies struggle to differentiate closely related microbial strains, due both to technical (sequencing error, amplification noise, bioinformatics approximations) and biological (lack of nucleotide variants in the amplified regions) limitations. Once data generation platforms reached the fidelity necessary to preserve amplicon biological variation when present, however, several computational approaches emerged to classify such sequences in the most strain-specific manner possible. Oligotyping and Minimum Entropy Decomposition (MED) rely on semi-supervised and unsupervised classification, respectively, of variant positions within otherwise-identical 16S amplicons that show statistically unusual distributions across sample sets (and are thus unlikely due to technical factors). Other types of sub-operational taxonomic unit (OTU) clustering have subsequently extended this intuition to “exact” or “amplicon” sequence variants (ESVs or ASVs, respectively) using statistical error modeling (e.g., DADA2) or filtering before or after sequence identity clustering (e.g., Deblur or UNOISE2). Strain-resolved 16S amplicons have been used with methods like these to very specifically link, e.g., Porphyromonas asaccharolytica ATCC 25260 and Parvimonas micra ATCC 33270 to CRC, for example, or to assess the temporal stability of strains in the gut. With additional data generation efforts, they can also generally be extended to multiple - or non-16S amplicons, such as the VaST system for identifying a minimum group of target loci for amplification. While SNV diversity in sub-regions of the genome is typically highly correlated with that across the genome, the presence or absence of at least one reliably detected SNV within a single amplified 16S variable region can be so precise as to become highly clade- and protocol-specific.

15Notably, the earliest forms of full-length 16S rRNA gene sequencing avoided many of these issues by capturing biological variation across the entire locus with high fidelity, and this has recently become true again in higher throughput with the advancement of “long-read” technologies. Three main platforms can currently provide such long-reads: Pacific Biosciences, Oxford Nanopore, and linked-read analogs such as products from 10X Genomics and Loop Genomics. The extreme fidelity offered by Pacific Biosciences circular consensus sequencing (CCS) has been perhaps best-studied in this context, readily differentiating between single-nucleotide variants (SNVs, although sometimes not insertions or deletions) when they exist anywhere across the 16S rRNA gene locus between strains. Conversely, while Oxford Nanopore’s extremely cost-effective MinION can provide essentially full-length 16S rRNA gene reads, its error rates have restricted strain-specific applications to cases in which no other sequences highly homologous to microbes of interest are present in a community. Finally, several protocols now exist facilitating “simulated” long- or linked-reads on a variety of platforms, but those which have reached commercial viability are yet to be formally evaluated for amplicon profiling of microbial communities. Similarly, these technologies can sometimes be applied to entire microbial genomes isolated from single cells (e.g., via sorting or microfluidics) or from cross-linked genome copies. This abrogates the need for true metagenomic assembly or binning, as described below, although again with few quantitative studies of these emerging technologies in existence for whole-community profiling at the strain level.

16Overall, shotgun metagenomic approaches provide a richer profile of microbial communities’ genetic compositions, as they can in principle identify structural or SNVs anywhere within any microbe’s genome (Table 1). Two broad classes of analyses are currently able to identify microbial strains, the first based on the alignment of metagenomic nucleotides (typically unassembled) to a reference set of genes or genomes. This is generally efficient and sensitive, but of course only possible when sufficiently similar reference genomes (or prior metagenomic assemblies) exist to permit direct mapping of metagenomic reads. Notably, “sufficiently similar” references need not be particularly high-identity with respect to a target metagenome. Instead, they must simply permit sufficient genome-wide mapping to identify SNVs or structural variants unique to strains in the community, which can be successful at up to several tens of percent overall nucleotide divergence.

17Broadly speaking, four classes of reference-based community strain identification algorithms currently exist. The first identifies the one or more reference genotypes closest to those in a given community, with quantification based on some algorithm for ambiguity-resolved read mapping (e.g., PathoScope, Sigma). The second identifies the dominant, potentially novel genotype (strain) per species; these include StrainPhlAn, MetaMLST, MetaSNV, and others. These generally require deeper sequencing (up to 10× or more coverage of the strains to be targeted) and differ in their choice of which reference sequences to map against (e.g., complete genomes vs. universal core genes vs. species-specific marker genes) and the method and stringency of SNV identification. A third class of reference-based methods will further attempt to identify multiple strains per species within a metagenome, such as ConStrains or DESMAN, requiring even deeper coverage and more stringent noise removal to prevent false positives. Finally, fourth, methods that rely on structural rather than SNV variants are generally more sensitive (appropriate for community members as rare as ~1× or lower coverage) and include PanPhlan (which can be combined with gene-targeted functional profilers such as HUMAnN), MIDAS, and others.

18Alternatively, when sufficiently similar reference genomes are not available, metagenomic assembly can be used for highly novel strain discovery. There is an inherent tension in assembly-based metagenomic strain profiling, as most assemblers seek to identify a single consensus sequence for each contig and require > 1× coverage of an entire genome (or region) to do so. This is appropriate when a single strain dominates its nearby phylogenetic space within a community, in which case less-common strains can be found by mapping metagenomic reads back to, e.g., a binned assembly and identifying nucleotide or structural variants roughly as one would within complete genomes. However, in the presence of too many closely related strains within a community, such a consensus sequence is not achievable in the first place, and most assemblers will not be able to provide a contig appropriate for mapping. Even when possible, this process can be further complicated by the high ecological and technical variability of microbial community assemblies, resulting in diverse coverage and confidence (dependent on sequencing depth and population strain admixture) and benefitting from manual inspection of putative variants. Algorithms facilitating this process include Latent Strain Analysis (LSA), which can refine strain-level taxonomy using covariant clusters across multiple related (e.g., longitudinal) samples. Similarly, DESMAN uses statistical models not unlike those for ASV calling in amplicon data to identify variant genotypes well-supported across multiple samples’ co-assembly. In a very few cases to date, strain variants within microbial communities have been identified via analogous differences in metatranscriptomic gene expression quantification, such as strain-specific variation in Eggerthella lenta metabolism of the cardiac drug digoxin.

19Whether from reference sequences or assemblies, SNV versus structural approaches are often complementary and can provide unique information regarding the same underlying community: SNVs (when detectable) identify finer-grained phylogenetic and evolutionary differences, but can be difficult to interpret functionally, whereas structural variants (i.e., gain or loss of full genes or genomic regions) have a lower limit of detection within communities and can speak directly to the biochemical roles of the affected genes (when known, Fig. 2). Unsurprisingly, each approach can provide different strengths and weaknesses. Structural variation can be captured well by reference-based approaches, which are sensitive to unique gene (non-)detection. However, it is very difficult to identify rearrangements (rather than gains or losses) using such techniques, and these are better identified by assembly-based methods instead (when they can be reliably differentiated from, e.g., chimeric assembly errors). Conversely, SNV variation can be well-captured by either reference- or assembly-based approaches—the former more sensitively for organisms with representative isolates, the latter less sensitively but for novel organisms—and by either pangenome or whole-genome mapping approaches, depending where the most uniquely identifying polymorphisms occur. Finally, both structural variation and, to a lesser extent, nucleotide variation are particularly driven in microbial communities by mechanisms of genetic mobility, including all forms of lateral transfer, gene gain/loss, mobile elements, plasmids, and phage integration.

20Other molecular technologies for microbial strain typing in communities are often limited to microbes that can be cultured or otherwise isolated, although advances in (semi-)automated anaerobic culture and nanoculture have made this feasible in high throughput as well. Particularly in clinical microbiology, near-strain variant typing via mass spectrometry peptide fingerprinting is commonplace for pathogen isolates, due to its rapid turnaround time and low cost per individual sample relative to sequencing. The technology has some of the same caveats as ASV identification from sequence amplicons introduced above, however: amino acid variants must exist between the strains of interest in the profiled proteins, at a level detectable above experimental noise, and must be classifiable to a taxon of origin in a reference database or by clustering. While in principle the same types of strain-level protein variants could be detected using MALDI-TOF MS technologies in culture-independent community extracts, such applications remain extremely challenging, and instead, community proteomics are currently more commonly analyzed in a gene- or taxon-centric way.

21Conversely, microbial imaging—arguably the first method for differentiating strains—has made the high-throughput leap to whole communities in several culture-independent forms that are, under appropriate circumstances, able to provide strain-level identification. In some cases, this can mean literally direct microscopy of microfluidically separated (or nanocultured) cells, using automated cell isolation and image analysis. More molecular techniques include spectral or combinatorial fluorescent in situ hybridization (Combinatorial Labeling and Spectral Imaging or CLASI-FISH), which can currently identify over a dozen microbes within a community while maintaining spatial structure. Along with related techniques such as multilabel FISH (MiL-FISH), this relies on the presence of sufficient genetic variants at the FISH-probed loci (often 16S rRNA gene regions) to be differentially bound by spectrally distinct probes, but can in some cases be extended to living bacteria. This is also true for other microbial probe imaging methods such as flow cytometry or light sheet microscopy, which can retain viable cells, but require probes or genetically manipulated microbes with loci capable of distinguishing between closely related strains.

22While many of these methods are in part or whole culture-independent, it is difficult to understate the importance of the “culturomics” renaissance in separating and characterizing microbial strain isolates from communities including the human microbiome. While pathogen epidemiology has long relied on comparative genomics among up to tens of thousands of isolates, it has only recently become efficient to carry out large-scale isolation of commensal organisms from human populations or individuals. Doing so, however, opens up the ability to identify strain-level differences among isolates of the same species among individuals, within an individual microbiome at different spatial locations, or over time. Once isolated, of course, such microbial strains can be characterized by any number of standard methods, including differences among growth curves or media, chemical (e.g., antimicrobial) resistance, metabolic flux profiling, or amplicon or shotgun sequencing. Alternatively, whole-community culture via chemostat bioreactors provides an intermediate environment in which strains that are rare in situ can sweep to dominance, or be perturbed in a controlled manner, to amplify differential phenotypes or sequences that may otherwise remain below the limit of detection. Finally, culture-based and culture-independent strain identification techniques blur in the areas of single-cell microbial isolation and microcolony growth from communities. Microfluidic techniques in this vein include gel microdroplets (GMDs) for single-cell amplification or phenotyping, as well as microfluidic streak plates (MSPs) that combine the specificity of single cells with the biomass of streaked colonies (if desired).

23Particularly when considering culture-based and ex vivo/in vitro/model system assays, the combination of culture-independent high-throughput epidemiology with subsequent strain isolation or manipulation opens up a world of possibilities for characterizing novel health-relevant strains in the human microbiome. This review has taken an essentially “top down” perspective, akin to forward genetics, in which strain-specific features of interest (SNVs, gene cassettes, metabolism, etc.) are identified by various means from human population studies. Such an approach leads naturally to the subsequent biochemical characterization of these variants, either via isolation from primary samples or by in silico retrieval of homologous sequences or related strains from databases or repositories (e.g., ATCC, BEI, DSMZ). Primary samples can be characterized as an entire community via gnotobiotics or continuous culture, or individual isolate strains grown, characterized, or (when possible) genetically manipulated. Such approaches dovetail nicely with “bottom up” approaches (analogous to reverse genetics) that identify and characterize health-relevant strains by directly beginning with isolates and assessing their phenotypes in gnotobiotic mono- or combinatorial colonization or, when possible, human feeding or microbiota transplant clinical trials.

24As introduced above, the precise definition of “strain” is somewhat fluid throughout biology, let alone in microbiology or microbial community biology. While it has most often referred to a single colony isolate culture in the past, the introduction of technologies and tools for precisely resolved genetic variant identification within microbial communities has led to increased broadening of the term. It is now used with some frequency to mean a subspecies or intraspecific clade with relatively low genetic diversity, defined by core or pangenomic identity, nucleotide identity within an amplicon such as the 16S rRNA gene, or the other genotyping or phenotypic similarities described above. As has increasingly been discussed in the literature for microbial systematics overall, this suggests the need for a more quantitative definition of strains or subspecies clades, particularly within naturally variant microbial communities. In the absence of a single consensus definition, it is extremely useful for individual studies to define their use of “strain” up front when describing culture-based or (especially) culture-independent microbial community research.

25Regardless of their precise definition, several emerging technologies offer exciting new approaches for identifying, isolating, and characterizing health-relevant strains in the human microbiome. Historically, microbial genetic variants not associated with an overt, acute phenotype have gone largely undetected, until the relatively recent availability of whole-community profiling techniques by which they can be efficiently captured. Truly single-cell approaches reliant on individual microbial separation have been so far difficult to apply to human epidemiology, with methods for eukaryotic cells not transferring well at scale to the heterogeneity of microbial cell wall biochemistry and methods from environmental community profiling difficult to apply to matrices as diverse as human stool or skin. In addition to bioengineering for cell separation and lysis, advances in low-input, low-noise DNA isolation, amplification, and sequencing will help to address this challenge, as will nanoculture approaches that inherently amplify genomes in vivo. Such methods for capturing strains from the human microbiome go hand-in-hand with additional technologies for characterizing them at scale, including cheaper experimental systems such as gut-on-chip or organoid variants that sit in between single isolate culture and rich gnotobiotic models. Ultimately, understanding human microbiome biology will require not just the detection of specific microbial genetic variants in communities, but their introduction and manipulation, including the theoretical ability to genetically perturb any microbial strain either after or even before isolation from its host community.

26Even in the absence of such technology, extensive work remains to be done to characterize the microbial strain diversity in the human microbiome that has already been uncovered. Of the tens of millions of gene families identified within the human microbiome, some ~ 75% are not biochemically characterized by anything more than (in some cases remote) homology to reference sequences, and ~ 25% are not closely homologous to any isolate open reading frames. This astounding pool of biochemical dark matter may be unsurprising to microbial bioprospectors, who have mined primarily environmental communities for novel enzymatic and antimicrobial function for decades. As such, it represents a remarkable potential for new bioactive discovery in human health as well, since human-associated microbes could easily be enriched for protein and metabolite products that modulate host responses. In many of the examples described above, successful associations of SNV or structural variants in the microbiome with human phenotypes or environmental factors have led to genes of unknown function. Strain-level epidemiology in the human microbiome can thus help to prioritize the daunting task of identifying and characterizing the “most interesting” novel microbial variants and products of greatest relevance to health.

27Finally, the ways in which better techniques for strain characterization in the microbiome can benefit human health are themselves diverse. Cheap, rapid, and reproducible methods to quantify microbiome SNVs and genetic variants across human populations will allow the identification of precise microbial risk factors, much as did the standardization of human genetics platforms for genome-wide association studies (GWAS). Also analogously to GWAS, microbial strains can thus provide prognostic or diagnostic biomarkers for disease risk or diagnosis, or hints as to their underlying molecular mechanisms. This has been the case for decades in for comparative genetics microbial isolates, and as the number and depth of metagenomes continues to increase, it will undoubtedly become practical in microbial communities as well. Conversely, features of strains found to be bioactive can be used to develop novel interventions for health maintenance or therapy. These can range from better targeting of existing fecal microbiota transplant (FMT) technologies based on donor or recipient strain content, to the rational design of synthetic FMTs, treatment response prediction for FMTs or prebiotics, or the eventual administration of genetically modified organisms or communities. Recent work in strain-level epidemiology of microbial communities and the human microbiome is thus one of many important, ongoing efforts to realize the microbiome’s substantial translational potential.