Despite the growing recognition of a host genetic effect on shaping gut microbiota composition, the genetic determinants of oral microbiota remain largely unexplored, especially in the context of oral diseases. Here, we performed a microbiome genome-wide association study in 2 independent cohorts of patients with oral squamous cell carcinoma (OSCC, n = 144 and 67) and an additional group of noncancer individuals (n = 104). Besides oral bacterial dysbiosis and signatures observed in OSCC, associations of 3 loci with the abundance of genus-level taxa and 4 loci with β diversity measures were detected (q < 0.05) at the discovery stage. The most significant hit (rs10906082 with the genus Lachnoanaerobaculum, P = 3.55 × 10–9 at discovery stage) was replicated in a second OSCC cohort. Moreover, the other 2 taxonomical associations, rs10973953 with the genus Kingella (P = 1.38 × 10–9) and rs4721629 with the genus Parvimonas (P = 3.53 × 10–8), were suggestive in the meta-analysis combining 2 OSCC cohorts. Further pathway analysis revealed that these loci were enriched for genes in regulation of oncogenic and angiogenic responses, implicating a genetic anchor to the oral microbiome in estimation of casual relationships with OSCC. Our findings delineate the role of host genotypes in influencing the structure of oral microbial communities.
AkiyamaMOkadaYKanaiMTakahashiAMomozawaYIkedaMIwataNIkegawaSHirataMMatsudaK, et al. 2017. Genome-wide association study identifies 112 new loci for body mass index in the Japanese population. Nat Genet. 49(10):1458–1467.
2.
BlekhmanRGoodrichJKHuangKSunQBukowskiRBellJTSpectorTDKeinanALeyREGeversD, et al. 2015. Host genetic variation impacts microbiome composition across human body sites. Genome Biol. 16(1):191.
3.
BonderMJKurilshikovATigchelaarEFMujagicZImhannFVilaAVDeelenPVatanenTSchirmerMSmeekensSP, et al. 2016. The effect of host genetics on the gut microbiome. Nat Genet. 48(11):1407–1412.
CorbyPMBretzWAHartTCSchorkNJWesselJLyons-WeilerJPasterBJ. 2007. Heritability of oral microbial species in caries-active and caries-free twins. Twin Res Hum Genet. 10(6):821–828.
6.
DickinsonREDallolABiecheIKrexDMortonDMaherERLatifF. 2004. Epigenetic inactivation of SLIT3 and SLIT1 genes in human cancers. Br J Cancer. 91(12):2071–2078.
7.
FurquimCPSoaresGMRibeiroLLAzcarate-PerilMAButzNRoachJMossKBonfimCTorres-PereiraCCTelesFR. 2017. The salivary microbiome and oral cancer risk: a pilot study in Fanconi anemia. J Dent Res. 96(3):292–299.
8.
GomezAEspinozaJLHarkinsDMLeongPSafferyRBockmannMTorralbaMKuelbsCKodukulaRInmanJ, et al. 2017. Host genetic control of the oral microbiome in health and disease. Cell Host Microbe. 22(3):269–278.e3.
9.
GoodrichJKDavenportERBeaumontMJacksonMAKnightROberCSpectorTDBellJTClarkAGLeyRE. 2016. Genetic determinants of the gut microbiome in UK twins. Cell Host Microbe. 19(5):731–743.
10.
Guerrero-PrestonRGodoy-VitorinoFJedlickaARodriguez-HilarioAGonzalezHBondyJLawsonFFolawiyoOMichailidiCDziedzicA, et al. 2016. 16S rRNA amplicon sequencing identifies microbiota associated with oral cancer, human papilloma virus infection and surgical treatment. Oncotarget. 7(32):51320–51334.
11.
HuWLiXChengRKeJLiuYMaMCaoYLiuD. 2019. NAV2 facilitates invasion of cutaneous melanoma cells by targeting SNAI2 through the GSK-3β/β-catenin pathway. Arch Dermatol Res. 311(5):399–410.
12.
HughesDABacigalupeRWangJRuhlemannMCTitoRYFalonyGJoossensMVieira-SilvaSHenckaertsLRymenansL, et al. 2020. Genome-wide associations of human gut microbiome variation and implications for causal inference analyses. Nat Microbiol. 5(9):1079–1087.
13.
IgartuaCDavenportERGiladYNicolaeDLPintoJOberC. 2017. Host genetic variation in mucosal immunity pathways influences the upper airway microbiome. Microbiome. 5(1):16.
14.
IshibashiTHayakawaHSekiguchiM. 2003. A novel mechanism for preventing mutations caused by oxidation of guanine nucleotides. EMBO Rep. 4(5):479–483.
15.
IshidaSKatoKTanakaMOdamakiTKuboRMitsuyamaEXiaoJZYamaguchiRUematsuSImotoS, et al. 2020. Genome-wide association studies and heritability analysis reveal the involvement of host genetics in the Japanese gut microbiota. Commun Biol. 3(1):686.
16.
IshiguroHShimokawaTTsunodaTTanakaTFujiiYNakamuraYFurukawaY. 2002. Isolation of HELAD1, a novel human helicase gene up-regulated in colorectal carcinomas. Oncogene. 21(41):6387–6394.
17.
JamiesonSMTsaiPKondratyevMKBudhaniPLiuASenzerNNChioreanEGJalalSINemunaitisJJKeeD, et al. 2018. Evofosfamide for the treatment of human papillomavirus–negative head and neck squamous cell carcinoma. JCI Insight. 3(16):e122204.
18.
JemalASiegelRWardEHaoYXuJMurrayTThunMJ. 2008. Cancer statistics, 2008. CA Cancer J Clin. 58(2):71–96.
19.
KessenbrockKPlaksVWerbZ. 2010. Matrix metalloproteinases: regulators of the tumor microenvironment. Cell. 141(1):52–67.
20.
KnightsDSilverbergMSWeersmaRKGeversDDijkstraGHuangHTylerADvan SommerenSImhannFStempakJM, et al. 2014. Complex host genetics influence the microbiome in inflammatory bowel disease. Genome Med. 6(12):107.
21.
KoldeRFranzosaEARahnavardGHallABVlamakisHStevensCDalyMJXavierRJHuttenhowerC. 2018. Host genetic variation and its microbiome interactions within the human microbiome project. Genome Med. 10(1):6.
22.
KurilshikovAMedina-GomezCBacigalupeRRadjabzadehDWangJDemirkanALe RoyCIRaygoza GarayJAFinnicumCTLiuX, et al. 2021. Large-scale association analyses identify host factors influencing human gut microbiome composition. Nat Genet. 53(2):156–165.
23.
La RosaGRMGattusoGPedullaERapisardaENicolosiDSalmeriM. 2020. Association of oral dysbiosis with oral cancer development. Oncol Lett. 19(4):3045–3058.
24.
LeeWHChenHMYangSFLiangCPengCYLinFMTsaiLLWuBCHsinCHChuangCY, et al. 2017. Bacterial alterations in salivary microbiota and their association in oral cancer. Sci Rep. 7(1):16540.
25.
LiYXieJXuXWangJAoFWanYZhuY. 2013. MicroRNA-548 down-regulates host antiviral response via direct targeting of IFN-λ1. Protein Cell. 4(2):130–141.
26.
LimYTotsikaMMorrisonMPunyadeeraC. 2017. Oral microbiome: a new biomarker reservoir for oral and oropharyngeal cancers. Theranostics. 7(17):4313–4321.
27.
LiuCXLiYObermoeller-McCormickLMSchwartzALBuG. 2001. The putative tumor suppressor LRP1B, a novel member of the low density lipoprotein (LDL) receptor family, exhibits both overlapping and distinct properties with the LDL receptor-related protein. J Biol Chem. 276(31):28889–28896.
28.
NarayanGGoparajuCArias-PulidoHKaufmannAMSchneiderADurstMMansukhaniMPothuriBMurtyVV. 2006. Promoter hypermethylation-mediated inactivation of multiple slit-robo pathway genes in cervical cancer progression. Mol Cancer. 5:16.
29.
NonesKWaddellNSongSPatchAMMillerDJohnsAWuJKassahnKSWoodDBaileyP, et al. 2014. Genome-wide DNA methylation patterns in pancreatic ductal adenocarcinoma reveal epigenetic deregulation of SLIT-ROBO, ITGA2 and MET signaling. Int J Cancer. 135(5):1110–1118.
30.
Norskov-LauritsenN. 2014. Classification, identification, and clinical significance of Haemophilus and Aggregatibacter species with host specificity for humans. Clin Microbiol Rev. 27(2):214–240.
31.
PapapostolouAKroffkeBTatakisDNNagarajaHNKumarPS. 2011. Contribution of host genotype to the composition of health-associated supragingival and subgingival microbiomes. J Clin Periodontol. 38(6):517–524.
32.
PiriyapongsaJJordanIK. 2007. A family of human microRNA genes from miniature inverted-repeat transposable elements. PLoS One. 2(2):e203.
33.
RothschildDWeissbrodOBarkanEKurilshikovAKoremTZeeviDCosteaPIGodnevaAKalkaINBarN, et al. 2018. Environment dominates over host genetics in shaping human gut microbiota. Nature. 555(7695):210–215.
34.
ScepanovicPHodelFMondotSPartulaVByrdAHammerCAlanioCBergstedtJPatinETouvierM, et al. 2019. A comprehensive assessment of demographic, environmental, and host genetic associations with gut microbiome diversity in healthy individuals. Microbiome. 7(1):130.
TannoTFujiwaraASakaguchiKTanakaKTakenakaSTsuyamaS. 2007. Slit3 regulates cell motility through Rac/Cdc42 activation in lipopolysaccharide-stimulated macrophages. FEBS Lett. 581(5):1022–1026.
37.
TurpinWEspin-GarciaOXuWSilverbergMSKevansDSmithMIGuttmanDSGriffithsAPanaccioneROtleyA, et al. 2016. Association of host genome with intestinal microbial composition in a large healthy cohort. Nat Genet. 48(11):1413–1417.
38.
WangJThingholmLBSkiecevicieneJRauschPKummenMHovJRDegenhardtFHeinsenFARuhlemannMCSzymczakS, et al. 2016. Genome-wide association analysis identifies variation in vitamin D receptor and other host factors influencing the gut microbiota. Nat Genet. 48(11):1396–1406.
39.
YangSFHuangHDFanWLJongYJChenMKHuangCNChuangCYKuoYLChungWHSuSC. 2018. Compositional and functional variations of oral microbiota associated with the mutational changes in oral cancer. Oral Oncol. 77:1–8.
40.
ZehnerMMarschallALBosESchloetelJGKreerCFehrenschildDLimmerAOssendorpFLangTKosterAJ, et al. 2015. The translocon protein Sec61 mediates antigen transport from endosomes in the cytosol for cross-presentation to CD8(+) T cells. Immunity. 42(5):850–863.
Supplementary Material
Please find the following supplemental material available below.
For Open Access articles published under a Creative Commons License, all supplemental material carries the same license as the article it is associated with.
For non-Open Access articles published, all supplemental material carries a non-exclusive license, and permission requests for re-use of supplemental material or any part of supplemental material shall be sent directly to the copyright owner as specified in the copyright notice associated with the article.
