Genetic manipulation of Patescibacteria provides mechanistic insights into microbial dark matter and the epibiotic lifestyle.

Patescibacteria, also known as the candidate phyla radiation (CPR), are a diverse group of bacteria that constitute a disproportionately large fraction of microbial dark matter. Its few cultivated members, belonging mostly to Saccharibacteria, grow as epibionts on host Actinobacteria. Due to a lack of suitable tools, the genetic basis of this lifestyle and other unique features of Patescibacteira remain unexplored. Here, we show that Saccharibacteria exhibit natural competence, and we exploit this property for their genetic manipulation. Imaging of fluorescent protein-labeled Saccharibacteria provides high spatiotemporal resolution of phenomena accompanying epibiotic growth, and a transposon-insertion sequencing (Tn-seq) genome-wide screen reveals the contribution of enigmatic Saccharibacterial genes to growth on their hosts. Finally, we leverage metagenomic data to provide cutting-edge protein structure-based bioinformatic resources that support the strain Southlakia epibionticum and its corresponding host, Actinomyces israelii, as a model system for unlocking the molecular underpinnings of the epibiotic lifestyle.

Episymbiotic Saccharibacteria TM7x modulates the susceptibility of its host bacteria to phage infection and promotes their coexistence.

Bacteriophages (phages) play critical roles in modulating microbial ecology. Within the human microbiome, the factors influencing the long-term coexistence of phages and bacteria remain poorly investigated. Saccharibacteria (formerly TM7) are ubiquitous members of the human oral microbiome. These ultrasmall bacteria form episymbiotic relationships with their host bacteria and impact their physiology. Here, we showed that during surface-associated growth, a human oral Saccharibacteria isolate (named TM7x) protects its host bacterium, a Schaalia odontolytica strain (named XH001) against lytic phage LC001 predation. RNA-Sequencing analysis identified in XH001 a gene cluster with predicted functions involved in the biogenesis of cell wall polysaccharides (CWP), whose expression is significantly down-regulated when forming a symbiosis with TM7x. Through genetic work, we experimentally demonstrated the impact of the expression of this CWP gene cluster on bacterial-phage interaction by affecting phage binding. In vitro coevolution experiments further showed that the heterogeneous populations of TM7x-associated and TM7x-free XH001, which display differential susceptibility to LC001 predation, promote bacteria and phage coexistence. Our study highlights the tripartite interaction between the bacterium, episymbiont, and phage. More importantly, we present a mechanism, i.e., episymbiont-mediated modulation of gene expression in host bacteria, which impacts their susceptibility to phage predation and contributes to the formation of "source-sink" dynamics between phage and bacteria in biofilm, promoting their long-term coexistence within the human microbiome.

Localized microbially induced inflammation influences distant healthy tissues in the human oral cavity.

Variation in human immune response to the same bacterial or viral pathogen is well established in the literature. Variation in immune response to microbial challenge has also been observed within the human oral cavity. Our recent study focused on characterizing observed variations in microbially induced gingival inflammation-resulting in three distinct clinical Inflammatory Responder Types (IRTs): High-IRT, Low-IRT, and Slow-IRT. Here, we applied a high-resolution temporal multiomic analysis during microbially induced inflammation in order to characterize the effects of localized oral inflammation on distant healthy tissues in young healthy adults. Our results highlight a nonlocalized subclinical effect with alterations in proinflammatory host mediators and an ecological shift toward dysbiosis within the subgingival microbiome in an IRT-dependent manner-despite maintained oral hygiene. Our results provide mechanistic insight into how healthy tissues within humans are influenced by distant localized inflammation and may ultimately become susceptible to disease.

Acquisition of the arginine deiminase system benefits epiparasitic Saccharibacteria and their host bacteria in a mammalian niche environment.

Saccharibacteria are a group of widespread and genetically diverse ultrasmall bacteria with highly reduced genomes that belong to the Candidate Phyla Radiation. Comparative genomic analyses suggest convergent evolution of key functions enabling the adaptation of environmental Saccharibacteria to mammalian microbiomes. Currently, our understanding of this environment-to-mammal niche transition within Saccharibacteria and their obligate episymbiotic association with host bacteria is limited. Here, we identified a complete arginine deiminase system (ADS), found in further genome streamlined mammal-associated Saccharibacteria but missing in their environmental counterparts, suggesting acquisition during environment-to-mammal niche transition. Using TM7x, the first cultured Saccharibacteria strain from the human oral microbiome and its host bacterium Actinomyces odontolyticus, we experimentally tested the function and impact of the ADS. We demonstrated that by catabolizing arginine and generating adenosine triphosphate, the ADS allows metabolically restrained TM7x to maintain higher viability and infectivity when disassociated from the host bacterium. Furthermore, the ADS protects TM7x and its host bacterium from acid stress, a condition frequently encountered within the human oral cavity due to bacterial metabolism of dietary carbohydrates. Intriguingly, with a restricted host range, TM7x forms obligate associations with Actinomyces spp. lacking the ADS but not those carrying the ADS, suggesting the acquired ADS may also contribute to partner selection for cooperative episymbiosis within a mammalian microbiome. These data present experimental characterization of a mutualistic interaction between TM7x and their host bacteria, and illustrate the benefits of acquiring a novel pathway in the transition of Saccharibacteria to mammalian microbiomes.

Human variation in gingival inflammation.

Oral commensal bacteria actively participate with gingival tissue to maintain healthy neutrophil surveillance and normal tissue and bone turnover processes. Disruption of this homeostatic host-bacteria relationship occurs during experimental gingivitis studies where it has been clearly established that increases in the bacterial burden increase gingival inflammation. Here, we show that experimental gingivitis resulted in three unique clinical inflammatory phenotypes (high, low, and slow) and reveal that interleukin-1ß, a reported major gingivitis-associated inflammatory mediator, was not associated with clinical gingival inflammation in the slow response group. In addition, significantly higher levels of Streptococcus spp. were also unique to this group. The low clinical response group was characterized by low concentrations of host mediators, despite similar bacterial accumulation and compositional characteristics as the high clinical response group. Neutrophil and bone activation modulators were down-regulated in all response groups, revealing novel tissue and bone protective responses during gingival inflammation. These alterations in chemokine and microbial composition responses during experimental gingivitis reveal a previously uncharacterized variation in the human host response to a disruption in gingival homeostasis. Understanding this human variation in gingival inflammation may facilitate the identification of periodontitis-susceptible individuals. Overall, this study underscores the variability in host responses in the human population arising from variations in host immune profiles (low responders) and microbial community maturation (slow responders) that may impact clinical outcomes in terms of destructive inflammation.

Exploitation of a Bacterium-Encoded Lytic Transglycosylase by a Human Oral Lytic Phage To Facilitate Infection.

Bacteriophages (phages) are an integral part of the human oral microbiome. Their roles in modulating bacterial physiology and shaping microbial communities have been discussed but remain understudied due to limited isolation and characterization of oral phage. Here, we report the isolation of LC001, a lytic phage targeting human oral Schaalia odontolytica (formerly known as Actinomyces odontolyticus) strain XH001. We showed that LC001 attached to and infected surface-grown, but not planktonic, XH001 cells, and it displayed remarkable host specificity at the strain level. Whole-genome sequencing of spontaneous LC001-resistant, surface-grown XH001 mutants revealed that the majority of the mutants carry nonsense or frameshift mutations in XH001 gene APY09_05145 (renamed ltg-1), which encodes a putative lytic transglycosylase (LT). The mutants are defective in LC001 binding, as revealed by direct visualization of the significantly reduced attachment of phage particles to the XH001 spontaneous mutants compared that to the wild type. Meanwhile, targeted deletion of ltg-1 produced a mutant that is defective in LC001 binding and resistant to LC001 infection even as surface-grown cells, while complementation of ltg-1 in the mutant background restored the LC001-sensitive phenotype. Intriguingly, similar expression levels of ltg-1 were observed in surface-grown and planktonic XH001, which displayed LC001-binding and nonbinding phenotypes, respectively. Furthermore, the overexpression of ltg-1 failed to confer an LC001-binding and -sensitive phenotype to planktonic XH001. Thus, our data suggested that rather than directly serving as a phage receptor, ltg-1-encoded LT may increase the accessibility of phage receptor, possibly via its enzymatic activity, by cleaving the peptidoglycan structure for better receptor exposure during peptidoglycan remodeling, a function that can be exploited by LC001 to facilitate infection. IMPORTANCE The evidence for the presence of a diverse and abundant phage population in the host-associated oral microbiome came largely from metagenomic analysis or the observation of virus-like particles within saliva/plaque samples, while the isolation of oral phage and investigation of their interaction with bacterial hosts are limited. Here, we report the isolation of LC001, the first lytic phage targeting oral Schaalia odontolytica. Our study suggested that LC001 may exploit the host bacterium-encoded lytic transglycosylase function to gain access to the receptor, thus facilitating its infection.

Transcriptome of Epibiont Saccharibacteria Nanosynbacter lyticus Strain TM7x During the Establishment of Symbiosis.

Saccharibacteria Nanosynbacter lyticus strain TM7x is a member of the broadly distributed candidate phylum radiation. These bacteria have ultrasmall cell sizes, have reduced genomes, and live as epibionts on the surfaces of other bacteria. The mechanisms by which they establish and maintain this relationship are not yet fully understood. The transcriptomes of the epibiont TM7x and its host bacteria Schaalia odontolytica strain XH001 were captured across the establishment of symbiosis during both the initial interaction and stable symbiosis. The results showed a dynamic interaction with large shifts in gene expression for both species between the initial encounter and stable symbiosis, notably in transporter genes. During stable symbiosis, the host XH001 showed higher gene expression for peptidoglycan biosynthesis, mannosylation, cell cycle and stress-related genes, whereas it showed lower expression of chromosomal partitioning genes. This was consistent with the elongated cell shape seen in XH001 infected with TM7x and our discovery that infection resulted in thickened cell walls. Within TM7x, increased pili, type IV effector genes, and arginine catabolism/biosynthesis gene expression during stable symbiosis implied a key role for these functions in the interaction. Consistent with its survival and persistence in the human microbiome as an obligate epibiont with reduced de novo biosynthetic capacities, TM7x also showed higher levels of energy production and peptidoglycan biosynthesis, but lower expression of stress-related genes, during stable symbiosis. These results imply that TM7x and its host bacteria keep a delicate balance in order to sustain an episymbiotic lifestyle. IMPORTANCE Nanosynbacter lyticus type strain TM7x is the first cultivated member of the Saccharibacteria and the candidate phyla radiation (CPR). It was discovered to be ultrasmall in cell size with a highly reduced genome that establishes an obligate epibiotic relationship with its host bacterium. The CPR is a large, monophyletic radiation of bacteria with reduced genomes that includes Saccharibacteria. The vast majority of the CPR have yet to be cultivated, and our insights into these unique organisms to date have been derived from only a few Saccharibacteria species. Being obligate parasites, it is unknown how these ultrasmall Saccharibacteria, which are missing many de novo biosynthetic pathways, are maintained at a high prevalence within the human microbiome as well as in the environment.

Klebsiella and Providencia emerge as lone survivors following long-term starvation of oral microbiota.

It is well-understood that many bacteria have evolved to survive catastrophic events using a variety of mechanisms, which include expression of stress-response genes, quiescence, necrotrophy, and metabolic advantages obtained through mutation. However, the dynamics of individuals leveraging these abilities to gain a competitive advantage in an ecologically complex setting remain unstudied. In this study, we observed the saliva microbiome throughout the ecological perturbation of long-term starvation, allowing only the species best equipped to access and use the limited resources to survive. During the first several days, the community underwent a death phase that resulted in a ∼50-100-fold reduction in the number of viable cells. Interestingly, after this death phase, only three species, Klebsiella pneumoniae, Klebsiella oxytoca, and Providencia alcalifaciens, all members of the family Enterobacteriaceae, appeared to be transcriptionally active and recoverable. Klebsiella are significant human pathogens, frequently resistant to multiple antibiotics, and recently, ectopic colonization of the gut by oral Klebsiella was documented to induce dysbiosis and inflammation. MetaOmics analyses provided several leads for further investigation regarding the ecological success of the Enterobacteriaceae. The isolates accumulated single nucleotide polymorphisms in known growth advantage in stationary phase alleles and produced natural products closely resembling antimicrobial cyclic depsipeptides. The results presented in this study suggest that pathogenic Enterobacteriaceae persist much longer than their more benign neighbors in the salivary microbiome when faced with starvation. This is particularly significant, given that hospital surfaces contaminated with oral fluids, especially sinks and drains, are well-established sources of outbreaks of drug-resistant Enterobacteriaceae.

Rapid evolution of decreased host susceptibility drives a stable relationship between ultrasmall parasite TM7x and its bacterial host.

Around one-quarter of bacterial diversity comprises a single radiation with reduced genomes, known collectively as the Candidate Phyla Radiation. Recently, we coisolated TM7x, an ultrasmall strain of the Candidate Phyla Radiation phylum Saccharibacteria, with its bacterial host Actinomyces odontolyticus strain XH001 from human oral cavity and stably maintained as a coculture. Our current work demonstrates that within the coculture, TM7x cells establish a long-term parasitic association with host cells by infecting only a subset of the population, which stay viable yet exhibit severely inhibited cell division. In contrast, exposure of a naïve A. odontolyticus isolate, XH001n, to TM7x cells leads to high numbers of TM7x cells binding to each host cell, massive host cell death, and a host population crash. However, further passaging reveals that XH001n becomes less susceptible to TM7x over time and enters a long-term stable relationship similar to that of XH001. We show that this reduced susceptibility is driven by rapid host evolution that, in contrast to many forms of phage resistance, offers only partial protection. The result is a stalemate where infected hosts cannot shed their parasites; nevertheless, parasite load is sufficiently low that the host population persists. Finally, we show that TM7x can infect and form stable long-term relationships with other species in a single clade of Actinomyces, displaying a narrow host range. This system serves as a model to understand how parasitic bacteria with reduced genomes such as those of the Candidate Phyla Radiation have persisted with their hosts and ultimately expanded in their diversity.

Recent advances in genomic DNA sequencing of microbial species from single cells.

The vast majority of microbial species remain uncultivated and, until recently, about half of all known bacterial phyla were identified only from their 16S ribosomal RNA gene sequence. With the advent of single-cell sequencing, genomes of uncultivated species are rapidly filling in unsequenced branches of the microbial phylogenetic tree. The wealth of new insights gained from these previously inaccessible groups is providing a deeper understanding of their basic biology, taxonomy and evolution, as well as their diverse roles in environmental ecosystems and human health.

Identification of PGN_1123 as the Gene Encoding Lipid A Deacylase, an Enzyme Required for Toll-Like Receptor 4 Evasion, in Porphyromonas gingivalis.

Removal of one acyl chain from bacterial lipid A by deacylase activity is a mechanism used by many pathogenic bacteria to evade the host's Toll-like receptor 4 (TLR4)-mediated innate immune response. In Porphyromonas gingivalis, a periodontal pathogen, lipid A deacylase activity converts a majority of the initially synthesized penta-acylated lipid A, a TLR4 agonist, to tetra-acylated structures, which effectively evade TLR4 sensing by being either inert or antagonistic at TLR4. In this paper, we report successful identification of the gene that encodes the P. gingivalis lipid A deacylase enzyme. This gene, PGN_1123 in P. gingivalis 33277, is highly conserved within P. gingivalis, and putative orthologs are phylogenetically restricted to the Bacteroidetes phylum. Lipid A of ΔPGN_1123 mutants is penta-acylated and devoid of tetra-acylated structures, and the mutant strain provokes a strong TLR4-mediated proinflammatory response, in contrast to the negligible response elicited by wild-type P. gingivalis Heterologous expression of PGN_1123 in Bacteroides thetaiotaomicron promoted lipid A deacylation, confirming that PGN_1123 encodes the lipid A deacylase enzyme.IMPORTANCE Periodontitis, commonly referred to as gum disease, is a chronic inflammatory condition that affects a large proportion of the population. Porphyromonas gingivalis is a bacterium closely associated with periodontitis, although how and if it is a cause for the disease are not known. It has a formidable capacity to dampen the host's innate immune response, enabling its persistence in diseased sites and triggering microbial dysbiosis in animal models of infection. P. gingivalis is particularly adept at evading the host's TLR4-mediated innate immune response by modifying the structure of lipid A, the TLR4 ligand. In this paper, we report identification of the gene encoding lipid A deacylase, a key enzyme that modifies lipid A to TLR4-evasive structures.

The Distinct Immune-Stimulatory Capacities of Porphyromonas gingivalis Strains 381 and ATCC 33277 Are Determined by the fimB Allele and Gingipain Activity.

The Porphyromonas gingivalis strain ATCC 33277 (33277) and 381 genomes are nearly identical. However, strain 33277 displays a significantly diminished capacity to stimulate host cell Toll-like receptor 2 (TLR2)-dependent signaling and interleukin-1ß (IL-1ß) production relative to 381, suggesting that there are strain-specific differences in one or more bacterial immune-modulatory factors. Genomic sequencing identified a single nucleotide polymorphism in the 33277 fimB allele (A→T), creating a premature stop codon in the 33277 fimB open reading frame relative to the 381 fimB allele. Gene exchange experiments established that the 33277 fimB allele reduces the immune-stimulatory capacity of this strain. Transcriptome comparisons revealed that multiple genes related to carboxy-terminal domain (CTD) family proteins, including the gingipains, were upregulated in 33277 relative to 381. A gingipain substrate degradation assay demonstrated that cell surface gingipain activity is higher in 33277, and an isogenic mutant strain deficient for the gingipains exhibited an increased ability to induce TLR2 signaling and IL-1ß production. Furthermore, 33277 and 381 mutant strains lacking CTD cell surface proteins were more immune-stimulatory than the parental wild-type strains, consistent with an immune-suppressive role for the gingipains. Our data show that the combination of an intact fimB allele and limited cell surface gingipain activity in P. gingivalis 381 renders this strain more immune-stimulatory. Conversely, a defective fimB allele and high-level cell surface gingipain activity reduce the capacity of P. gingivalis 33277 to stimulate host cell innate immune responses. In summary, genomic and transcriptomic comparisons identified key virulence characteristics that confer divergent host cell innate immune responses to these highly related P. gingivalis strains.

Discovery of a Novel Periodontal Disease-Associated Bacterium.

One of the world's most common infectious disease, periodontitis (PD), derives from largely uncharacterized communities of oral bacteria growing as biofilms (a.k.a. plaque) on teeth and gum surfaces in periodontal pockets. Bacteria associated with periodontal disease trigger inflammatory responses in immune cells, which in later stages of the disease cause loss of both soft and hard tissue structures supporting teeth. Thus far, only a handful of bacteria have been characterized as infectious agents of PD. Although deep sequencing technologies, such as whole community shotgun sequencing have the potential to capture a detailed picture of highly complex bacterial communities in any given environment, we still lack major reference genomes for the oral microbiome associated with PD and other diseases. In recent work, by using a combination of supervised machine learning and genome assembly, we identified a genome from a novel member of the Bacteroidetes phylum in periodontal samples. Here, by applying a comparative metagenomics read-classification approach, including 272 metagenomes from various human body sites, and our previously assembled draft genome of the uncultivated Candidatus Bacteroides periocalifornicus (CBP) bacterium, we show CBP's ubiquitous distribution in dental plaque, as well as its strong association with the well-known pathogenic "red complex" that resides in deep periodontal pockets.

A Linear Plasmid-Like Prophage of Actinomyces odontolyticus Promotes Biofilm Assembly.

The human oral cavity is home to a large number of bacteria and bacteriophages (phages). However, the biology of oral phages as members of the human microbiome is not well understood. Recently, we isolated Actinomyces odontolyticus subsp. actinosynbacter strain XH001 from the human oral cavity, and genomic analysis revealed the presence of an intact prophage named xhp1. Here, we demonstrated that xhp1 is a linear plasmid-like prophage, which is a newly identified phage of A. odontolyticus The prophage xhp1 genome is a 35-kb linear double-stranded DNA with 10-bp single-stranded, 3' cohesive ends. xhp1 exists extrachromosomally, with an estimated copy number of 5. Annotation of xhp1 revealed 54 open reading frames, while phylogenetic analysis suggests that it has limited similarity with other phages. xhp1 phage particles can be induced by mitomycin C and belong to the Siphoviridae family, according to transmission electron microscopic examination. The released xhp1 particles can reinfect the xhp1-cured XH001 strain and result in tiny blurry plaques. Moreover, xhp1 promotes XH001 biofilm formation through spontaneous induction and the release of host extracellular DNA (eDNA). In conclusion, we identified a linear plasmid-like prophage of A. odontolyticus, which enhances bacterial host biofilm assembly and could be beneficial to the host for its persistence in the oral cavity.IMPORTANCE The biology of phages as members of the human oral microbiome is understudied. Here, we report the characterization of xhp1, a novel linear plasmid-like prophage identified from a human oral isolate, Actinomyces odontolyticus subsp. actinosynbacter strain XH001. xhp1 can be induced and reinfect xhp1-cured XH001. The spontaneous induction of xhp1 leads to the lysis of a subpopulation of bacterial hosts and the release of eDNA that promotes biofilm assembly, thus potentially contributing to the persistence of A. odontolyticus within the oral cavity.

Precision-guided antimicrobial peptide as a targeted modulator of human microbial ecology.

One major challenge to studying human microbiome and its associated diseases is the lack of effective tools to achieve targeted modulation of individual species and study its ecological function within multispecies communities. Here, we show that C16G2, a specifically targeted antimicrobial peptide, was able to selectively kill cariogenic pathogen Streptococcus mutans with high efficacy within a human saliva-derived in vitro oral multispecies community. Importantly, a significant shift in the overall microbial structure of the C16G2-treated community was revealed after a 24-h recovery period: several bacterial species with metabolic dependency or physical interactions with S. mutans suffered drastic reduction in their abundance, whereas S. mutans' natural competitors, including health-associated Streptococci, became dominant. This study demonstrates the use of targeted antimicrobials to modulate the microbiome structure allowing insights into the key community role of specific bacterial species and also indicates the therapeutic potential of C16G2 to achieve a healthy oral microbiome.

Cultivation of a human-associated TM7 phylotype reveals a reduced genome and epibiotic parasitic lifestyle.

The candidate phylum TM7 is globally distributed and often associated with human inflammatory mucosal diseases. Despite its prevalence, the TM7 phylum remains recalcitrant to cultivation, making it one of the most enigmatic phyla known. In this study, we cultivated a TM7 phylotype (TM7x) from the human oral cavity. This extremely small coccus (200-300 nm) has a distinctive lifestyle not previously observed in human-associated microbes. It is an obligate epibiont of an Actinomyces odontolyticus strain (XH001) yet also has a parasitic phase, thereby killing its host. This first completed genome (705 kb) for a human-associated TM7 phylotype revealed a complete lack of amino acid biosynthetic capacity. Comparative genomics analyses with uncultivated environmental TM7 assemblies show remarkable conserved gene synteny and only minimal gene loss/gain that may have occurred as TM7x adapted to conditions within the human host. Transcriptomic and metabolomic profiles provided the first indications, to our knowledge, that there is signaling interaction between TM7x and XH001. Furthermore, the induction of TNF-α production in macrophages by XH001 was repressed in the presence of TM7x, suggesting its potential immune suppression ability. Overall, our data provide intriguing insights into the uncultivability, pathogenicity, and unique lifestyle of this previously uncharacterized oral TM7 phylotype.

hybridSPAdes: an algorithm for hybrid assembly of short and long reads.

MOTIVATION: Recent advances in single molecule real-time (SMRT) and nanopore sequencing technologies have enabled high-quality assemblies from long and inaccurate reads. However, these approaches require high coverage by long reads and remain expensive. On the other hand, the inexpensive short reads technologies produce accurate but fragmented assemblies. Thus, a hybrid approach that assembles long reads (with low coverage) and short reads has a potential to generate high-quality assemblies at reduced cost. RESULTS: We describe hybridSPAdes algorithm for assembling short and long reads and benchmark it on a variety of bacterial assembly projects. Our results demonstrate that hybridSPAdes generates accurate assemblies (even in projects with relatively low coverage by long reads) thus reducing the overall cost of genome sequencing. We further present the first complete assembly of a genome from single cells using SMRT reads. AVAILABILITY AND IMPLEMENTATION: hybridSPAdes is implemented in C++ as a part of SPAdes genome assembler and is publicly available at http://bioinf.spbau.ru/en/spades CONTACT: d.antipov@spbu.ru SUPPLEMENTARY INFORMATION: supplementary data are available at Bioinformatics online.

Human Oral Buccal Microbiomes Are Associated with Farmworker Status and Azinphos-Methyl Agricultural Pesticide Exposure.

In a longitudinal agricultural community cohort sampling of 65 adult farmworkers and 52 adult nonfarmworkers, we investigated agricultural pesticide exposure-associated changes in the oral buccal microbiota. We found a seasonally persistent association between the detected blood concentration of the insecticide azinphos-methyl and the taxonomic composition of the buccal swab oral microbiome. Blood and buccal samples were collected concurrently from individual subjects in two seasons, spring/summer 2005 and winter 2006. Mass spectrometry quantified blood concentrations of the organophosphate insecticide azinphos-methyl. Buccal oral microbiome samples were 16S rRNA gene DNA sequenced, assigned to the bacterial taxonomy, and analyzed after "centered-log-ratio" transformation to handle the compositional nature of the proportional abundances of bacteria per sample. Nonparametric analysis of the transformed microbiome data for individuals with and without azinphos-methyl blood detection showed significant perturbations in seven common bacterial taxa (>0.5% of sample mean read depth), including significant reductions in members of the common oral bacterial genus Streptococcus Diversity in centered-log-ratio composition between individuals' microbiomes was also investigated using principal-component analysis (PCA) to reveal two primary PCA clusters of microbiome types. The spring/summer "exposed" microbiome cluster with significantly less bacterial diversity was enriched for farmworkers and contained 27 of the 30 individuals who also had azinphos-methyl agricultural pesticide exposure detected in the blood. IMPORTANCE: In this study, we show in human subjects that organophosphate pesticide exposure is associated with large-scale significant alterations of the oral buccal microbiota composition, with extinctions of whole taxa suggested in some individuals. The persistence of this association from the spring/summer to the winter also suggests that long-lasting effects on the commensal microbiota have occurred. The important health-related outcomes of these agricultural community individuals' pesticide-associated microbiome perturbations are not understood at this time. Future investigations should index medical and dental records for common and chronic diseases that may be interactively caused by this association between pesticide exposure and microbiome alteration.

Metagenome and Metatranscriptome Analyses Using Protein Family Profiles.

Analyses of metagenome data (MG) and metatranscriptome data (MT) are often challenged by a paucity of complete reference genome sequences and the uneven/low sequencing depth of the constituent organisms in the microbial community, which respectively limit the power of reference-based alignment and de novo sequence assembly. These limitations make accurate protein family classification and abundance estimation challenging, which in turn hamper downstream analyses such as abundance profiling of metabolic pathways, identification of differentially encoded/expressed genes, and de novo reconstruction of complete gene and protein sequences from the protein family of interest. The profile hidden Markov model (HMM) framework enables the construction of very useful probabilistic models for protein families that allow for accurate modeling of position specific matches, insertions, and deletions. We present a novel homology detection algorithm that integrates banded Viterbi algorithm for profile HMM parsing with an iterative simultaneous alignment and assembly computational framework. The algorithm searches a given profile HMM of a protein family against a database of fragmentary MG/MT sequencing data and simultaneously assembles complete or near-complete gene and protein sequences of the protein family. The resulting program, HMM-GRASPx, demonstrates superior performance in aligning and assembling homologs when benchmarked on both simulated marine MG and real human saliva MG datasets. On real supragingival plaque and stool MG datasets that were generated from healthy individuals, HMM-GRASPx accurately estimates the abundances of the antimicrobial resistance (AMR) gene families and enables accurate characterization of the resistome profiles of these microbial communities. For real human oral microbiome MT datasets, using the HMM-GRASPx estimated transcript abundances significantly improves detection of differentially expressed (DE) genes. Finally, HMM-GRASPx was used to reconstruct comprehensive sets of complete or near-complete protein and nucleotide sequences for the query protein families. HMM-GRASPx is freely available online from http://sourceforge.net/projects/hmm-graspx.

Genome of the pathogen Porphyromonas gingivalis recovered from a biofilm in a hospital sink using a high-throughput single-cell genomics platform.

Although biofilms have been shown to be reservoirs of pathogens, our knowledge of the microbial diversity in biofilms within critical areas, such as health care facilities, is limited. Available methods for pathogen identification and strain typing have some inherent restrictions. In particular, culturing will yield only a fraction of the species present, PCR of virulence or marker genes is mainly focused on a handful of known species, and shotgun metagenomics is limited in the ability to detect strain variations. In this study, we present a single-cell genome sequencing approach to address these limitations and demonstrate it by specifically targeting bacterial cells within a complex biofilm from a hospital bathroom sink drain. A newly developed, automated platform was used to generate genomic DNA by the multiple displacement amplification (MDA) technique from hundreds of single cells in parallel. MDA reactions were screened and classified by 16S rRNA gene PCR sequence, which revealed a broad range of bacteria covering 25 different genera representing environmental species, human commensals, and opportunistic human pathogens. Here we focus on the recovery of a nearly complete genome representing a novel strain of the periodontal pathogen Porphyromonas gingivalis (P. gingivalis JCVI SC001) using the single-cell assembly tool SPAdes. Single-cell genomics is becoming an accepted method to capture novel genomes, primarily in the marine and soil environments. Here we show for the first time that it also enables comparative genomic analysis of strain variation in a pathogen captured from complex biofilm samples in a healthcare facility.