Human genes influence the interaction between Streptococcus mutans and host caries susceptibility: a genome-wide association study in children with primary dentition

Abstract

Streptococcus mutans is a well-known cause of dental caries, due to its acidogenicity, aciduricity, and ability to synthesize exopolysaccharides in dental plaques. Intriguingly, not all children who carry S. mutans manifest caries, even with similar characteristics in oral hygiene, diet, and other environmental factors. This phenomenon suggests that host susceptibility potentially plays a role in the development of dental caries; however, the association between host genetics, S. mutans, and dental caries remains unclear. Therefore, this study examined the influence of host gene-by-S. mutans interaction on dental caries. Genome-wide association analyses were conducted in 709 US children (<13 years old), using the dbGap database acquired from the center for oral health research in appalachia (COHRA) and the Iowa Head Start programmes (GEIRS). A generalized estimating equation was used to examine the gene-by-S. mutans interaction effects on the outcomes (decayed and missing/filled primary teeth due to caries). Sequentially, the COHRA and GEIRS data were used to identify potential interactions and replicate the findings. Three loci at the genes interleukin 32 (IL32), galactokinase 2 (GALK2), and CUGBP, Elav-like family member 4 (CELF4) were linked to S. mutans carriage, and there was a severity of caries at a suggestive significance level among COHRA children (P < 9 × 10⁻⁵), and at a nominal significance level among GEIRS children (P = 0.047–0.001). The genetic risk score that combined the three loci also significantly interacted with S. mutans (P < 0.000 1). Functional analyses indicated that the identified genes are involved in the host immune response, galactose carbohydrate metabolism, and food-rewarding system, which could potentially be used to identify children at high risk for caries and to develop personalized caries prevention strategies.

References

[1] Dye BA, Li X, Thorton-Evans G. Oral health disparities as determined by selected healthy people 2020 oral health objectives for the United States, 2009–2010. NCHS. Data Brief., 2012, 104: 1-8.

[2] Dye BA, . Trends in oral health status: United States, 1988–1994 and 1999–2004. Vital. Health Stat., 2007, 11: 1-92.

[3] GBD 2016 Disease and Injury Incidence and Prevalence Collaborators. Global, regional, and national incidence, prevalence, and years lived with disability for 328 diseases and injuries for 195 countries, 1990-2016: a systematic analysis for the global burden of disease study 2016. Lancet, 2017, 390: 1211-1259.

[4] Tanzer JM. Dental caries is a transmissible infectious disease: the keyes and fitzgerald revolution. J. Dent. Res., 1995, 74: 1536-1542.

[5] Bowen WH. Dental caries—not just holes in teeth! a perspective. Mol. Oral. Microbiol., 2016, 31: 228-233.

[6] Takahashi N, Nyvad B. The role of bacteria in the caries process: ecological perspectives. J. Dent. Res., 2011, 90: 294-303.

[7] van Houte J. Role of micro-organisms in caries etiology. J. Dent. Res., 1994, 73: 672-681.

[8] Banas JA. Virulence properties of Streptococcus mutans. Front. Biosci., 2004, 9: 1267-1277.

[9] Bowen WH, Koo H. Biology of Streptococcus mutans-derived glucosyltransferases: role in extracellular matrix formation of cariogenic biofilms. Caries Res., 2011, 45: 69-86.

[10] Kanasi E, . Microbial risk markers for childhood caries in pediatricians’ offices. J. Dent. Res., 2010, 89: 378-383.

[11] Kowash MB, Curzon ME, Hart P. Association of salivary Streptococcus mutans with caries in young children: effect of dental health education on salivary levels. Eur. J. Paediatr. Dent., 2002, 3: 199-204.

[12] Burt BA, Loesche WJ, Eklund SA. Stability of selected plaque species and their relationship to caries in a child population over 2 years. Caries Res., 1985, 19: 193-200.

[13] Divaris K. Precision dentistry in early childhood: the central role of genomics. Dent. Clin. North Am., 2017, 61: 619-625.

[14] Shaffer JR, . Effects of enamel matrix genes on dental caries are moderated by fluoride exposures. Hum. Genet., 2015, 134: 159-167.

[15] Ballantine JL, . Exploring the genomic basis of early childhood caries: a pilot study. Int. J. Paediatr. Dent., 2018, 28: 217-225.

[16] Horikoshi M, . Genome-wide associations for birth weight and correlations with adult disease. Nature, 2016, 538: 248-252.

[17] Morrison J, . Genome-wide association study of dental caries in the hispanic communities health study/study of latinos (HCHS/SOL). Hum. Mol. Genet., 2016, 25: 807-816.

[18] Shaffer JR, . Genome-wide association scan for childhood caries implicates novel genes. J. Dent. Res., 2011, 90: 1457-1462.

[19] Zeng Z, . Genome-wide association studies of pit-and-fissure- and smooth-surface caries in permanent dentition. J. Dent. Res., 2013, 92: 432-437.

[20] Justice AE, . Genome-wide meta-analysis of 241,258 adults accounting for smoking behaviour identifies novel loci for obesity traits. Nat. Commun., 2017, 8

[21] Khawar MB, Abbasi MH, Sheikh N. IL-32: a novel pluripotent inflammatory interleukin, towards gastric inflammation, gastric cancer, and chronic rhino sinusitis. Mediators Inflamm., 2016, 2016: 8413768.

[22] Cregger R, Langworthy K, Salako N, Streckfus C. Relationship between salivarycytokines, and caries experience in children with different body mass indices. J. Dent. Oral. Health, 2017, 5: 2.

[23] Gornowicz A, . Pro-inflammatory cytokines in saliva of adolescents with dental caries disease. Ann. Agric. Environ. Med., 2012, 19: 711-716.

[24] Kim DH, . Intracellular interleukin-32γ mediates antiviral activity of cytokines against hepatitis B virus. Nat. Commun., 2018, 9

[25] Gorvel L, Korenfeld D, Tung T, Klechevsky E. Dendritic cell-derived IL-32α: a novel inhibitory cytokine of NK cell function. J. Immunol., 2017, 199: 1290-1300.

[26] Cho JS, . Toll-like receptor 4-mediated expression of interleukin-32 via the c-Jun N-terminal kinase/protein kinase B/cyclic adenosine monophosphate response element binding protein pathway in chronic rhinosinusitis with nasal polyps. Int. Forum Allergy Rhinol., 2016, 6: 1020-1028.

[27] Bai X, Dinarello CA, Chan ED. The role of interleukin-32 against tuberculosis. Cytokine, 2015, 76: 585-587.

[28] Dede FO, Balli U, Doğan SB, Güven B. Interleukin-32 levels in gingival crevicular fluid and saliva of patients with chronic periodontitis after periodontal treatment. J. Periodontal Res., 2017, 52: 397-407.

[29] Pastuszak I, O’Donnell J, Elbein AD. Identification of the GalNAc kinase amino acid sequence. J. Biol. Chem., 1996, 271: 23653-23656.

[30] Agnew A, Timson D. Mechanistic studies on human N-acetylgalactosamine kinase. J. Enzyme Inhib. Med. Chem., 2010, 25: 370-376.

[31] Riedl B, Schmid W. A stereoselective and flexible synthesis to access both enantiomers of N-acetylgalactosamine and peracetylated N-acetylidosamine. Beilstein J. Org. Chem., 2018, 14: 856-860.

[32] Blixt O, . Analysis of Tn antigenicity with a panel of new IgM and IgG1 monoclonal antibodies raised against leukemic cells. Glycobiology, 2012, 22: 529-542.

[33] Coelho AI, Berry GT, Rubio-Gozalbo ME. Galactose metabolism and health. Curr. Opin. Clin. Nutr. Metab. Care, 2015, 18: 422-427.

[34] Zeng L, Das S, Burne RA. Utilization of lactose and galactose by Streptococcus mutans: transport, toxicity, and carbon catabolite repression. J. Bacteriol., 2010, 192: 2434-2444.

[35] Wagnon JL, . CELF4 regulates translation and local abundance of a vast set of mRNAs, including genes associated with regulation of synaptic function. PLoS Genet., 2012, 8

[36] Yang Y, . Complex seizure disorder caused by Brunol4 deficiency in mice. PLoS Genet., 2007, 3

[37] Leigh Sarah-Jane, Morris Margaret J.. The role of reward circuitry and food addiction in the obesity epidemic: An update. Biological Psychology, 2018, 131: 31-42.

[38] Alonso-Alonso M, . Food reward system: current perspectives and future research needs. Nutr. Rev., 2015, 73: 296-307.

[39] Xiao J, . Association between oral candida and bacteriome in children with severe ECC. J. Dent. Res., 2018, 97: 1468-1476.

[40] Xiao J, . Candida albicans and early childhood caries: a systematic review and meta-analysis. Caries Res., 2018, 52: 102-112.

[41] Gross EL, . Beyond Streptococcus mutans: dental caries onset linked to multiple species by 16S rRNA community analysis. PLoS ONE, 2012, 7

[42] Johansson I, Witkowska E, Kaveh B, Holgerson PL, Tanner AC. The microbiome in populations with a low and high prevalence of caries. J. Dent. Res., 2016, 95: 80-86.

[43] Ma C, . Comparison of oral microbial profiles between children with severe early childhood caries and caries-free children using the human oral microbe identification microarray. PLoS ONE, 2015, 10

[44] Richards VP, . Microbiomes of site-specific dental plaques from children with different caries status. Infect. Immun., 2017, 85: e00106-e00117.

[45] Tanner AC, . Microbiota of severe early childhood caries before and after therapy. J. Dent. Res., 2011, 90: 1298-1305.

[46] Yang F, . Saliva microbiomes distinguish caries-active from healthy human populations. ISME J., 2012, 6: 1-10.

[47] Stanley BO, . Genetic association of MPPED2 and ACTN2 with dental caries. J. Dent. Res., 2014, 93: 626-632.

[48] Chang CC, . Second-generation PLINK: rising to the challenge of larger and richer datasets. Gigascience, 2015, 4

[49] Purcell S, . PLINK: a tool set for whole-genome association and population-based linkage analyses. Am. J. Hum. Genet., 2007, 81: 559-575.

[50] Chen MH, Yang Q. GWAF: an R package for genome-wide association analyses with family data. Bioinformatics, 2010, 26: 580-581.

[51] Zhao W, . Interaction between social/psychosocial factors and genetic variants on body mass index: a gene-environment interaction analysis in a longitudinal setting. Int. J. Environ. Res. Public Health, 2017, 14: 1153.

[52] Nettleton JA, . Gene x dietary pattern interactions in obesity: analysis of up to 68 317 adults of European ancestry. Hum. Mol. Genet., 2015, 24: 4728-4738.

[53] Storey JD, Tibshirani R. Statistical significance for genomewide studies. Proc. Natl Acad. Sci. USA, 2003, 100: 9440-9445.

[54] Ward LD, Kellis M. HaploReg: a resource for exploring chromatin states, conservation, and regulatory motif alterations within sets of genetically linked variants. Nucleic Acids Res, 2012, 40: D930-D934.

[55] Kuleshov MV, . Enrichr: a comprehensive gene set enrichment analysis web server 2016 update. Nucleic Acids Res., 2016, 44: W90-W97.

[56] Chen J, Bardes EE, Aronow BJ, Jegga AG. ToppGene suite for gene list enrichment analysis and candidate gene prioritization. Nucleic Acids Res., 2009, 37: W305-W311.

[57] Little J, . Strengthening the reporting of genetic association studies (STREGA): an extension of the STROBE statement. Eur. J. Clin. Invest., 2009, 39: 247-266.

Funding

U.S. Department of Health & Human Services | NIH | National Institute of Dental and Craniofacial Research (NIDCR)(K23DE027412)