Spatial ecology of Haemophilus and Aggregatibacter in the human oral cavity.

Haemophilus and Aggregatibacter are two of the most common bacterial genera in the human oral cavity, encompassing both commensals and pathogens of substantial ecological and medical significance. In this study, we conducted a metapangenomic analysis of oral Haemophilus and Aggregatibacter species to uncover genomic diversity, phylogenetic relationships, and habitat specialization within the human oral cavity. Using three metrics-pangenomic gene content, phylogenomics, and average nucleotide identity (ANI)-we first identified distinct species and sub-species groups among these genera. Mapping of metagenomic reads then revealed clear patterns of habitat specialization, such as Aggregatibacter species predominantly in dental plaque, a distinctive Haemophilus parainfluenzae sub-species group on the tongue dorsum, and H. sp. HMT-036 predominantly in keratinized gingiva and buccal mucosa. In addition, we found that supragingival plaque samples contained predominantly only one out of the three taxa, H. parainfluenzae, Aggregatibacter aphrophilus, and A. sp. HMT-458, suggesting independent niches or a competitive relationship. Functional analyses revealed the presence of key metabolic genes, such as oxaloacetate decarboxylase, correlated with habitat specialization, suggesting metabolic versatility as a driving force. Additionally, heme synthesis distinguishes H. sp. HMT-036 from closely related Haemophilus haemolyticus, suggesting that the availability of micronutrients, particularly iron, was important in the evolutionary ecology of these species. Overall, our study exemplifies the power of metapangenomics to identify factors that may affect ecological interactions within microbial communities, including genomic diversity, habitat specialization, and metabolic versatility. IMPORTANCE: Understanding the microbial ecology of the mouth is essential for comprehending human physiology. This study employs metapangenomics to reveal that various Haemophilus and Aggregatibacter species exhibit distinct ecological preferences within the oral cavity of healthy individuals, thereby supporting the site-specialist hypothesis. Additionally, it was observed that the gene pool of different Haemophilus species correlates with their ecological niches. These findings shed light on the significance of key metabolic functions in shaping microbial distribution patterns and interspecies interactions in the oral ecosystem.

Site-specialization of human oral Gemella species.

Gemella species are core members of the human oral microbiome in healthy subjects and are regarded as commensals, although they can cause opportunistic infections. Our objective was to evaluate the site-specialization of Gemella species among various habitats within the mouth by combining pangenomics and metagenomics. With pangenomics, we identified genome relationships and categorized genes as core and accessory to each species. With metagenomics, we identified the primary oral habitat of individual genomes. Our results establish that the genomes of three species, G. haemolysans, G. sanguinis and G. morbillorum, are abundant and prevalent in human mouths at different oral sites: G. haemolysans on buccal mucosa and keratinized gingiva; G. sanguinis on tongue dorsum, throat, and tonsils; and G. morbillorum in dental plaque. The gene-level basis of site-specificity was investigated by identifying genes that were core to Gemella genomes at a specific oral site but absent from other Gemella genomes. The riboflavin biosynthesis pathway was present in G. haemolysans genomes associated with buccal mucosa but absent from the rest of the genomes. Overall, metapangenomics show that Gemella species have clear ecological preferences in the oral cavity of healthy humans and provides an approach to identifying gene-level drivers of site specificity.

Site Specialization of Human Oral Veillonella Species.

Veillonella species are abundant members of the human oral microbiome with multiple interspecies commensal relationships. Examining the distribution patterns of Veillonella species across the oral cavity is fundamental to understanding their oral ecology. In this study, we used a combination of pangenomic analysis and oral metagenomic information to clarify Veillonella taxonomy and to test the site specialist hypothesis for the Veillonella genus, which contends that most oral bacterial species are adapted to live at specific oral sites. Using isolate genome sequences combined with shotgun metagenomic sequence data, we showed that Veillonella species have clear, differential site specificity: Veillonella parvula showed strong preference for supra- and subgingival plaque, while closely related V. dispar, as well as more distantly related V. atypica, preferred the tongue dorsum, tonsils, throat, and hard palate. In addition, the provisionally named Veillonella sp. Human Microbial Taxon 780 showed strong site specificity for keratinized gingiva. Using comparative genomic analysis, we identified genes associated with thiamine biosynthesis and the reductive pentose phosphate cycle that may enable Veillonella species to occupy their respective habitats. IMPORTANCE Understanding the microbial ecology of the mouth is fundamental for understanding human physiology. In this study, metapangenomics demonstrated that different Veillonella species have clear ecological preferences in the oral cavity of healthy humans, validating the site specialist hypothesis. Furthermore, the gene pool of different Veillonella species was found to be reflective of their ecology, illuminating the potential role of vitamins and carbohydrates in determining Veillonella distribution patterns and interspecies interactions.

Genome-Centric Dynamics Shape the Diversity of Oral Bacterial Populations.

Two major viewpoints have been put forward for how microbial populations change, differing in whether adaptation is driven principally by gene-centric or genome-centric processes. Longitudinal sampling at microbially relevant timescales, i.e., days to weeks, is critical for distinguishing these mechanisms. Because of its significance for both microbial ecology and human health and its accessibility and high level of curation, we used the oral microbiota to study bacterial intrapopulation genome dynamics. Metagenomes were generated by shotgun sequencing of total community DNA from the healthy tongues of 17 volunteers at four to seven time points obtained over intervals of days to weeks. We obtained 390 high-quality metagenome-assembled genomes (MAGs) defining population genomes from 55 genera. The vast majority of genes in each MAG were tightly linked over the 2-week sampling window, indicating that the majority of the population's genomes were temporally stable at the MAG level. MAG-defined populations were composed of up to 5 strains, as determined by single-nucleotide-variant frequencies. Although most were stable over time, individual strains carrying over 100 distinct genes that rose from low abundance to dominance in a population over a period of days were detected. These results indicate a genome-wide as opposed to a gene-level process of population change. We infer that genome-wide selection of ecotypes is the dominant mode of adaptation in the oral populations over short timescales. IMPORTANCE The oral microbiome represents a microbial community of critical relevance to human health. Recent studies have documented the diversity and dynamics of different bacteria to reveal a rich, stable ecosystem characterized by strain-level dynamics. However, bacterial populations and their genomes are neither monolithic nor static; their genomes are constantly evolving to lose, gain, or alter their functional potential. To better understand how microbial genomes change in complex communities, we used culture-independent approaches to reconstruct the genomes (MAGs) for bacterial populations that approximated different species, in 17 healthy donors' mouths over a 2-week window. Our results underscored the importance of strain-level dynamics, which agrees with and expands on the conclusions of previous research. Altogether, these observations reveal patterns of genomic dynamics among strains of oral bacteria occurring over a matter of days.

Corncob structures in dental plaque reveal microhabitat taxon specificity.

BACKGROUND: The human mouth is a natural laboratory for studying how bacterial communities differ across habitats. Different bacteria colonize different surfaces in the mouth-teeth, tongue dorsum, and keratinized and non-keratinized epithelia-despite the short physical distance between these habitats and their connection through saliva. We sought to determine whether more tightly defined microhabitats might have more tightly defined sets of resident bacteria. A microhabitat may be characterized, for example, as the space adjacent to a particular species of bacterium. Corncob structures of dental plaque, consisting of coccoid bacteria bound to filaments of Corynebacterium cells, present an opportunity to analyze the community structure of one such well-defined microhabitat within a complex natural biofilm. Here, we investigate by fluorescence in situ hybridization and spectral imaging the composition of the cocci decorating the filaments. RESULTS: The range of taxa observed in corncobs was limited to a small subset of the taxa present in dental plaque. Among four major groups of dental plaque streptococci, two were the major constituents of corncobs, including one that was the most abundant Streptococcus species in corncobs despite being relatively rare in dental plaque overall. Images showed both Streptococcus types in corncobs in all individual donors, suggesting that the taxa have different ecological roles or that mechanisms exist for stabilizing the persistence of functionally redundant taxa in the population. Direct taxon-taxon interactions were observed not only between the Streptococcus cells and the central corncob filament but also between Streptococcus cells and the limited subset of other plaque bacteria detected in the corncobs, indicating species ensembles involving these taxa as well. CONCLUSIONS: The spatial organization we observed in corncobs suggests that each of the microbial participants can interact with multiple, albeit limited, potential partners, a feature that may encourage the long-term stability of the community. Additionally, our results suggest the general principle that a precisely defined microhabitat will be inhabited by a small and well-defined set of microbial taxa. Thus, our results are important for understanding the structure and organizing principles of natural biofilms and lay the groundwork for future work to modulate and control biofilms for human health. Video Abstract.

Site-tropism of streptococci in the oral microbiome.

A detailed understanding of where bacteria localize is necessary to advance microbial ecology and microbiome-based therapeutics. The site-specialist hypothesis predicts that most microbes in the human oral cavity have a primary habitat type within the mouth where they are most abundant. We asked whether this hypothesis accurately describes the distribution of the members of the genus Streptococcus, a clinically relevant taxon that dominates most oral sites. Prior analysis of 16S rRNA gene sequencing data indicated that some oral Streptococcus clades are site-specialists while others may be generalists. However, within complex microbial populations composed of numerous closely related species and strains, such as the oral streptococci, genome-scale analysis is necessary to provide the resolution to discriminate closely related taxa with distinct functional roles. Here, we assess whether individual species within this genus are specialists using publicly available genomic sequence data that provide species-level resolution. We chose a set of high-quality representative genomes for human oral Streptococcus species. Onto these genomes, we mapped shotgun metagenomic sequencing reads from supragingival plaque, tongue dorsum, and other sites in the oral cavity. We found that every abundant Streptococcus species in the healthy human oral cavity showed strong site-tropism and that even closely related species such as S. mitis, S. oralis, and S. infantis specialized in different sites. These findings indicate that closely related bacteria can have distinct habitat distributions in the absence of dispersal limitation and under similar environmental conditions and immune regimes. Substantial overlap between the core genes of these three species suggests that site-specialization is determined by subtle differences in genomic content.

Spatial scale in analysis of the dental plaque microbiome.

Ecologists have long recognized the importance of spatial scale in understanding structure-function relationships among communities of organisms within their environment. Here, we review historical and contemporary studies of dental plaque community structure in the context of three distinct scales: the micro (1-10 µm), meso (10-100 µm) and macroscale (100 µm to ≥1 cm). Within this framework, we analyze the compositional nature of dental plaque at the macroscale, the molecular interactions of microbes at the microscale, and the emergent properties of dental plaque biofilms at the mesoscale. Throughout our analysis of dental plaque across spatial scales, we draw attention to disease and health-associated structure-function relationships and include a discussion of host immune involvement in the mesoscale structure of periodontal disease-associated biofilms. We end with a discussion of two filamentous organisms, Fusobacterium nucleatum and Corynebacterium matruchotii, and their relevant contributions in structuring dental plaque biofilms.

No man's land: Species-specific formation of exclusion zones bordering Actinomyces graevenitzii microcolonies in nanoliter cultures.

To survive within complex environmental niches, including the human host, bacteria have evolved intricate interspecies communities driven by competition for limited nutrients, cooperation via complementary metabolic proficiencies, and establishment of homeostatic relationships with the host immune system. The study of such complex, interdependent relationships is often hampered by the challenges of culturing many bacterial strains in research settings and the limited set of tools available for studying the dynamic behavior of multiple bacterial species at the microscale. Here, we utilize a microfluidic-based co-culture system and time-lapse imaging to characterize dynamic interactions between Streptococcus species, Staphylococcus aureus, and Actinomyces species. Co-culture of Streptococcus cristatus or S. salivarius in nanoliter compartments with Actinomyces graevenitzii revealed localized exclusion of Streptococcus and Staphylococcus from media immediately surrounding A. graevenitzii microcolonies. This community structure did not occur with S. mitis or S. oralis strains or in co-cultures containing other Actinomycetaceae species such as S. odontolyticus or A. naeslundii. Moreover, fewer neutrophils were attracted to compartments containing both A. graevenitzii and Staphylococcus aureus than to an equal number of either species alone, suggesting a possible survival benefit together during immune responses.

Metapangenomics of the oral microbiome provides insights into habitat adaptation and cultivar diversity.

BACKGROUND: The increasing availability of microbial genomes and environmental shotgun metagenomes provides unprecedented access to the genomic differences within related bacteria. The human oral microbiome with its diverse habitats and abundant, relatively well-characterized microbial inhabitants presents an opportunity to investigate bacterial population structures at an ecosystem scale. RESULTS: Here, we employ a metapangenomic approach that combines public genomes with Human Microbiome Project (HMP) metagenomes to study the diversity of microbial residents of three oral habitats: tongue dorsum, buccal mucosa, and supragingival plaque. For two exemplar taxa, Haemophilus parainfluenzae and the genus Rothia, metapangenomes reveal distinct genomic groups based on shared genome content. H. parainfluenzae genomes separate into three distinct subgroups with differential abundance between oral habitats. Functional enrichment analyses identify an operon encoding oxaloacetate decarboxylase as diagnostic for the tongue-abundant subgroup. For the genus Rothia, grouping by shared genome content recapitulates species-level taxonomy and habitat preferences. However, while most R. mucilaginosa are restricted to the tongue as expected, two genomes represent a cryptic population of R. mucilaginosa in many buccal mucosa samples. For both H. parainfluenzae and the genus Rothia, we identify not only limitations in the ability of cultivated organisms to represent populations in their native environment, but also specifically which cultivar gene sequences are absent or ubiquitous. CONCLUSIONS: Our findings provide insights into population structure and biogeography in the mouth and form specific hypotheses about habitat adaptation. These results illustrate the power of combining metagenomes and pangenomes to investigate the ecology and evolution of bacteria across analytical scales.

Oral Microbiome Geography: Micron-Scale Habitat and Niche.

The mouth presents a multiplicity of local environments in communication with one another via saliva. The spatial organization of microbes within the mouth is shaped by opposing forces in dynamic equilibrium-salivary flow and adhesion, shedding and colonization-and by interactions among and between microbes and the host. Here we review recent evidence confirming that oral microbes are specialized for individual habitats within the mouth and that microbial habitats and niches are defined by micron-scale gradients in combination with short- and long-range interactions. Micron-scale structure illuminates the roles of individual taxa and provides insight into their community ecology and potential pathogenicity.

Spatial Ecology of the Human Tongue Dorsum Microbiome.

A fundamental question in microbial ecology is how microbes are spatially organized with respect to each other and their host. A test bed for examining this question is the tongue dorsum, which harbors a complex and important microbial community. Here, we use multiplexed fluorescence spectral imaging to investigate the organization of the tongue microbiome at micron to hundred-micron scales. We design oligonucleotide probes for taxa both abundant and prevalent, as determined by sequence analysis. Imaging reveals a highly structured spatial organization of microbial consortia, ranging in linear dimension from tens to hundreds of microns. The consortia appear to develop from a core of epithelial cells, with taxa clustering in domains suggestive of clonal expansion. Quantitative proximity analysis provides the basis for a model of tongue dorsum microbiome organization and dynamics. Our work illustrates how high-resolution analysis of micron-scale organization provides insights into physiological functions and microbiome-host interactions.

Semi-blind sparse affine spectral unmixing of autofluorescence-contaminated micrographs.

MOTIVATION: Spectral unmixing methods attempt to determine the concentrations of different fluorophores present at each pixel location in an image by analyzing a set of measured emission spectra. Unmixing algorithms have shown great promise for applications where samples contain many fluorescent labels; however, existing methods perform poorly when confronted with autofluorescence-contaminated images. RESULTS: We propose an unmixing algorithm designed to separate fluorophores with overlapping emission spectra from contamination by autofluorescence and background fluorescence. First, we formally define a generalization of the linear mixing model, called the affine mixture model (AMM), that specifically accounts for background fluorescence. Second, we use the AMM to derive an affine nonnegative matrix factorization method for estimating fluorophore endmember spectra from reference images. Lastly, we propose a semi-blind sparse affine spectral unmixing (SSASU) algorithm that uses knowledge of the estimated endmembers to learn the autofluorescence and background fluorescence spectra on a per-image basis. When unmixing real-world spectral images contaminated by autofluorescence, SSASU greatly improved proportion indeterminacy as compared to existing methods for a given relative reconstruction error. AVAILABILITY AND IMPLEMENTATION: The source code used for this paper was written in Julia and is available with the test data at https://github.com/brossetti/ssasu.

Biogeography of the Oral Microbiome: The Site-Specialist Hypothesis.

Microbial communities are complex and dynamic, composed of hundreds of taxa interacting across multiple spatial scales. Advances in sequencing and imaging technology have led to great strides in understanding both the composition and the spatial organization of these complex communities. In the human mouth, sequencing results indicate that distinct sites host microbial communities that not only are distinguishable but to a meaningful degree are composed of entirely different microbes. Imaging suggests that the spatial organization of these communities is also distinct. Together, the literature supports the idea that most oral microbes are site specialists. A clear understanding of microbiota structure at different sites in the mouth enables mechanistic studies, informs the generation of hypotheses, and strengthens the position of oral microbiology as a model system for microbial ecology in general.

Systematic evasion of the restriction-modification barrier in bacteria.

Bacteria that are recalcitrant to genetic manipulation using modern in vitro techniques are termed genetically intractable. Genetic intractability is a fundamental barrier to progress that hinders basic, synthetic, and translational microbiology research and development beyond a few model organisms. The most common underlying causes of genetic intractability are restriction-modification (RM) systems, ubiquitous defense mechanisms against xenogeneic DNA that hinder the use of genetic approaches in the vast majority of bacteria and exhibit strain-level variation. Here, we describe a systematic approach to overcome RM systems. Our approach was inspired by a simple hypothesis: if a synthetic piece of DNA lacks the highly specific target recognition motifs for a host's RM systems, then it is invisible to these systems and will not be degraded during artificial transformation. Accordingly, in this process, we determine the genome and methylome of an individual bacterial strain and use this information to define the bacterium's RM target motifs. We then synonymously eliminate RM targets from the nucleotide sequence of a genetic tool in silico, synthesize an RM-silent "SyngenicDNA" tool, and propagate the tool as minicircle plasmids, termed SyMPL (SyngenicDNA Minicircle Plasmid) tools, before transformation. In a proof-of-principle of our approach, we demonstrate a profound improvement (five orders of magnitude) in the transformation of a clinically relevant USA300 strain of Staphylococcus aureus This stealth-by-engineering SyngenicDNA approach is effective, flexible, and we expect in future applications could enable microbial genetics free of the restraints of restriction-modification barriers.

Lévy-like movement patterns of metastatic cancer cells revealed in microfabricated systems and implicated in vivo.

Metastatic cancer cells differ from their non-metastatic counterparts not only in terms of molecular composition and genetics, but also by the very strategy they employ for locomotion. Here, we analyzed large-scale statistics for cells migrating on linear microtracks to show that metastatic cancer cells follow a qualitatively different movement strategy than their non-invasive counterparts. The trajectories of metastatic cells display clusters of small steps that are interspersed with long "flights". Such movements are characterized by heavy-tailed, truncated power law distributions of persistence times and are consistent with the Lévy walks that are also often employed by animal predators searching for scarce prey or food sources. In contrast, non-metastatic cancerous cells perform simple diffusive movements. These findings are supported by preliminary experiments with cancer cells migrating away from primary tumors in vivo. The use of chemical inhibitors targeting actin-binding proteins allows for "reprogramming" the Lévy walks into either diffusive or ballistic movements.

Preservation of three-dimensional spatial structure in the gut microbiome.

Preservation of three-dimensional structure in the gut is necessary in order to analyze the spatial organization of the gut microbiota and gut luminal contents. In this study, we evaluated preparation methods for mouse gut with the goal of preserving micron-scale spatial structure while performing fluorescence imaging assays. Our evaluation of embedding methods showed that commonly used media such as Tissue-Tek Optimal Cutting Temperature (OCT) compound, paraffin, and polyester waxes resulted in redistribution of luminal contents. By contrast, a hydrophilic methacrylate resin, Technovit H8100, preserved three-dimensional organization. Our mouse intestinal preparation protocol optimized using the Technovit H8100 embedding method was compatible with microbial fluorescence in situ hybridization (FISH) and other labeling techniques, including immunostaining and staining with both wheat germ agglutinin (WGA) and 4', 6-diamidino-2-phenylindole (DAPI). Mucus could be visualized whether the sample was fixed with paraformaldehyde (PFA) or with Carnoy's fixative. The protocol optimized in this study enabled simultaneous visualization of micron-scale spatial patterns formed by microbial cells in the mouse intestines along with biogeographical landmarks such as host-derived mucus and food particles.

Spatial organization of a model 15-member human gut microbiota established in gnotobiotic mice.

Knowledge of the spatial organization of the gut microbiota is important for understanding the physical and molecular interactions among its members. These interactions are thought to influence microbial succession, community stability, syntrophic relationships, and resiliency in the face of perturbations. The complexity and dynamism of the gut microbiota pose considerable challenges for quantitative analysis of its spatial organization. Here, we illustrate an approach for addressing this challenge, using (i) a model, defined 15-member consortium of phylogenetically diverse, sequenced human gut bacterial strains introduced into adult gnotobiotic mice fed a polysaccharide-rich diet, and (ii) in situ hybridization and spectral imaging analysis methods that allow simultaneous detection of multiple bacterial strains at multiple spatial scales. Differences in the binding affinities of strains for substrates such as mucus or food particles, combined with more rapid replication in a preferred microhabitat, could, in principle, lead to localized clonally expanded aggregates composed of one or a few taxa. However, our results reveal a colonic community that is mixed at micrometer scales, with distinct spatial distributions of some taxa relative to one another, notably at the border between the mucosa and the lumen. Our data suggest that lumen and mucosa in the proximal colon should be conceptualized not as stratified compartments but as components of an incompletely mixed bioreactor. Employing the experimental approaches described should allow direct tests of whether and how specified host and microbial factors influence the nature and functional contributions of "microscale" mixing to the dynamic operations of the microbiota in health and disease.

Multiplexed Spectral Imaging of 120 Different Fluorescent Labels.

The number of fluorescent labels that can unambiguously be distinguished in a single image when acquired through band pass filters is severely limited by the spectral overlap of available fluorophores. The recent development of spectral microscopy and the application of linear unmixing algorithms to spectrally recorded image data have allowed simultaneous imaging of fluorophores with highly overlapping spectra. However, the number of distinguishable fluorophores is still limited by the unavoidable decrease in signal to noise ratio when fluorescence signals are fractionated over multiple wavelength bins. Here we present a spectral image analysis algorithm to greatly expand the number of distinguishable objects labeled with binary combinations of fluorophores. Our algorithm utilizes a priori knowledge about labeled specimens and imposes a binary label constraint on the unmixing solution. We have applied our labeling and analysis strategy to identify microbes labeled by fluorescence in situ hybridization and here demonstrate the ability to distinguish 120 differently labeled microbes in a single image.

Individuality, Stability, and Variability of the Plaque Microbiome.

Dental plaque is a bacterial biofilm composed of a characteristic set of organisms. Relatively little information from cultivation-independent, high-throughput analyses has been published on the temporal dynamics of the dental plaque microbiome. We used Minimum Entropy Decomposition, an information theory-based approach similar to oligotyping that provides single-nucleotide resolution, to analyze a previously published time series data set and investigate the dynamics of the plaque microbiome at various analytic and taxonomic levels. At both the genus and 97% Operational Taxonomic Unit (OTU) levels of resolution, the range of variation within each individual overlapped that of other individuals in the data set. When analyzed at the oligotype level, however, the overlap largely disappeared, showing that single-nucleotide resolution enables differentiation of individuals from one another without ambiguity. The overwhelming majority of the plaque community in all samples was made up of bacteria from a moderate number of plaque-typical genera, indicating that the overall community framework is shared among individuals. Each of these genera fluctuated in abundance around a stable mean that varied between individuals, with some genera having higher inter-individual variability than others. Thus, at the genus level, differences between individuals lay not in the identity of the major genera but in consistently differing proportions of these genera from mouth to mouth. However, at the oligotype level, we detected oligotype "fingerprints," a highly individual-specific set of persistently abundant oligotypes fluctuating around a stable mean over time. For example, within the genus Corynebacterium, more than a dozen oligotypes were detectable in each individual, of which a different subset reached high abundance in any given person. This pattern suggests that each mouth contains a subtly different community of organisms. We also compared the Chinese plaque community characterized here to previously characterized Western plaque communities, as represented by analyses of data emerging from the Human Microbiome Project, and found no major differences between Chinese and Western supragingival plaque. In conclusion, we found the plaque microbiome to be highly individualized at the oligotype level and characterized by stability of community membership, with variability in the relative abundance of community members between individuals and over time.