Sub-species niche specialization in the oral microbiome is associated with nasopharyngeal carcinoma risk

We explored the relationship between NPC and the oral microbiome using V3-4 (341f-805r) 16s rRNA sequencing in a study of 499 NPC patients and 495 population-based age and sex frequency-matched controls from an endemic area of Southern China. We found a significant reduction in community richness in cases compared to controls. Differences in the overall microbial community structure between cases and controls could not be explained by other potential confounders; disease status explained 5 times more variation in the Unweighted UniFrac distance than the next most explanatory variable. In feature-based analyses, we identified a pair of co-excluding Granulicatella adiacens (G. adicans) amplicon sequence variants (ASVs) which were strongly associated with NPC status and differed by a single nucleotide. The G. adicans variant an individual carried was also associated with the overall microbial community based on beta diversity. Co-occurrence analysis suggested the two G. adicans ASVs sit at the center of two co-excluding clusters of closely related organisms. Our results suggest there are differences in the oral microbiome between NPC patients and healthy controls and these may be associated with both a loss of microbial diversity and niche specialization among closely related commensals.