Uncovering oral Neisseria tropism and persistence using metagenomic sequencing

Microbial epidemiology and population genomics have previously been carried out near-exclusively for organisms grown in vitro. Metagenomics helps to overcome this limitation, but it is still challenging to achieve strain-level characterization of microorganisms from culture-independent data with sufficient resolution for epidemiological modelling. Here, we have developed multiple complementary approaches that can be combined to profile and track individual microbial strains. To specifically profile highly recombinant neisseriae from oral metagenomes, we integrated four metagenomic analysis techniques: single nucleotide polymorphisms in the clade's core genome, DNA uptake sequence signatures, metagenomic multilocus sequence typing and strain-specific marker genes. We applied these tools to 520 oral metagenomes from the Human Microbiome Project, finding evidence of site tropism and temporal intra-subject strain retention. Although the opportunistic pathogen Neisseria meningitidis is enriched for colonization in the throat, N. flavescens and N. subflava populate the tongue dorsum, and N. sicca, N. mucosa and N. elongata the gingival plaque. The buccal mucosa appeared as an intermediate ecological niche between the plaque and the tongue. The resulting approaches to metagenomic strain profiling are generalizable and can be extended to other organisms and microbiomes across environments. Strain level analysis of metagenomic data defines the tissue tropism and persistence of human oral Neisseria species.