Next Article in Journal
Higher Circulating Levels of Neutrophils and Basophils Are Linked to Myopic Retinopathy
Next Article in Special Issue
The Epigenetics of Migraine
Previous Article in Journal
p53 Family in Resistance to Targeted Therapy of Melanoma
Previous Article in Special Issue
Gene Expression and DNA Methylation in Human Papillomavirus Positive and Negative Head and Neck Squamous Cell Carcinomas
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Association between Downstream Taste Signaling Genes, Oral Microbiome, and Severe Early Childhood Caries

by
Vivianne Cruz de Jesus
1,2,
Betty-Anne Mittermuller
1,2,3,
Pingzhao Hu
2,4,5,
Robert J. Schroth
1,2,3,6 and
Prashen Chelikani
1,2,4,7,*
1
Manitoba Chemosensory Biology Research Group, Department of Oral Biology, University of Manitoba, Winnipeg, MB R3E 0W2, Canada
2
Children’s Hospital Research Institute of Manitoba (CHRIM), Winnipeg, MB R3E 3P4, Canada
3
Department of Preventive Dental Science, University of Manitoba, Winnipeg, MB R3E 0W2, Canada
4
Department of Biochemistry and Medical Genetics, University of Manitoba, Winnipeg, MB R3E 3N4, Canada
5
Department of Biochemistry, Western University, London, ON N6A 5C1, Canada
6
Department of Pediatrics and Child Health, University of Manitoba, Winnipeg, MB R3A 1S1, Canada
7
Department of Physiology and Pathophysiology, University of Manitoba, Winnipeg, MB R3E 0J9, Canada
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2023, 24(1), 81; https://doi.org/10.3390/ijms24010081
Submission received: 22 November 2022 / Revised: 14 December 2022 / Accepted: 19 December 2022 / Published: 21 December 2022
(This article belongs to the Special Issue Genetic and Epigenetic Control of Disease Occurrence)

Abstract

:
Polymorphisms in taste receptor genes have been shown to play a role in early childhood caries (ECC), a multifactorial, biofilm-mediated disease. This study aimed to evaluate associations between severe-ECC (S-ECC), the oral microbiome, and variants in genes that encode components of the G protein-coupled receptor (GPCR) signaling cascade involved in taste sensation. A total of 176 children (88 caries-free; 88 with S-ECC) were recruited. Analyses of 16S and ITS1 rRNA microbial genes and seven (GNAQ, GNAS, GNAT3, GNAI2, RAC1, RALB, and PLCB2) human genes were pursued using next-generation sequencing. Regression analyses were performed to evaluate associations between genetic variants, S-ECC, and the supragingival plaque microbiome. Results suggest that PLCB2 rs2305645 (T), rs1869901 (G), and rs2305649 (G) alleles had a protective effect on S-ECC (rs2305645, odds ratio (OR) = 0.27 (95% confidence interval (CI): 0.14–0.51); rs1869901, OR = 0.34 (95% CI: 0.20–0.58); and rs2305649, OR = 0.43 (95% CI: 0.26–0.71)). Variants in GNAQ, GNAS, GNAT3, PLCB2, RALB, and RAC1 were associated with oral fungal and bacterial community composition. This study revealed that three loci at PLCB2 are significantly associated with S-ECC. Variants in multiple genes were associated with the composition of dental biofilm. These findings contribute to the current knowledge about the role of genetics in S-ECC.

1. Introduction

Most taste modalities are sensed through G protein-coupled receptors (GPCRs). They represent the largest receptor superfamily in humans and are composed of seven canonical transmembrane-spanning proteins [1]. GPCRs can elicit cellular and physiological responses by converting external stimuli into intracellular signals by recognizing and binding sensory input and ligands, including odors, tastes, bioamines, lipids, peptides, and proteins [1]. Chemosensory receptors such as odorant and taste receptors account for over half of the GPCR superfamily. Several studies have highlighted the diverse roles of these chemosensory receptors and have been extensively reviewed [2,3,4,5].
The signal transduction of GPCRs involves different components such as cell surface receptors, GTP-binding proteins (G proteins), second messengers, and effector molecules. GPCR-mediated calcium (Ca2+) signaling is essential for several physiological functions and dysregulations in this process are associated with many diseases, such as Alzheimer’s, Huntington’s disease, and Parkinson’s disease [6]. In a previous study, we showed that GPCRs and other proteins involved in detecting different taste signals, such as epithelial sodium channels, otopetrin 1, and carbonic anhydrase 6, are associated with severe early childhood caries (S-ECC), and the oral microbial composition of young children [7]. Here, we tested the hypothesis that variants in downstream components of the G protein-mediated GPCR signaling cascade may have similar effects.
Variants in genes that encode components of the GPCR signaling cascade may interfere with the host-microbial interactions mediated by taste receptors, as well as with taste preferences and dietary choices, leading to increased risk of S-ECC. Furthermore, Ca2+ is a well-known intracellular second messenger involved in GPCR signaling. In the context of oral health, Ca2+ signaling is important for tooth formation and components of the GPCR signaling cascade may be involved in the development of tooth and bone defects [8,9,10]. For instance, a study showed that increased Gαq activity in ameloblasts is linked to dental fluorosis, a tooth defect caused by increased ingestion of fluoride during enamel formation [11]. Activity of another G protein, Gαs, has been shown to modulate bone health and development [12]. Sweet taste receptors T1R2 and T1R3 have been associated with modulation of bone development, with reports of T1R2 and T1R3 knockout mice showing increased bone mass [13]. Additionally, defects in free fatty acid receptors FFAR1 and FFAR4 have been implicated in decreased bone formation and increased bone resorption in mice models [8]. As there are similarities in the process of bone and tooth formation, proteins involved in bone defects may also influence the formation of tooth defects [14]. Furthermore, tooth defects, along with dietary choices, which can be determined by taste preferences mediated by taste receptors, are known risk factors for early childhood caries (ECC) [15,16,17,18].
The canonical bitter taste signaling cascade mediated by heterotrimeric G-proteins suggests that the activation of taste receptors by a ligand leads to the activation of the intracellular heterotrimeric G-protein complex, Gαβγ. Then, the Gβγ dimer activates phospholipase C β2 (PLCβ2), which leads to Ca2+ release from the endoplasmic reticulum. Meanwhile, the Gαgustducin subunit activates phosphodiesterase to decrease cyclic AMP (cAMP) levels [2,4,19,20,21]. Studies have suggested that multiple Gα subunits, such as Gαi, Gαs, and Gαq, expressed in taste buds may also be involved in taste receptor signaling [22,23].
Upon activation of bitter taste receptors (T2Rs), the released second messenger Ca2+ binds to Ca2+-binding proteins, such as calmodulin, to continue the calcium-dependent signal transduction. Research has shown that calmodulin can interact with small GTPases such as Rac1 and RalB, suggesting that they may be involved in many calcium/calmodulin-mediated intracellular signaling pathways [24,25]. Recently, it was shown that Rac1 GTPase activity was regulated by quinine with the signal mediated by G-protein and T2R4 [26].
The objective of this study was to analyze the association between variants in downstream taste signaling genes, the oral microbiome, and S-ECC. While there are dozens of downstream proteins that might be potentially involved in human taste signaling, based on the previous published literature and to keep this objective feasible, the focus was only on seven downstream proteins. Thus, associations between variants in genes encoding cell signaling components, such as heterotrimeric G proteins (GNAT3, GNAS, GNAI2, and GNAQ), a downstream effector (PLCB2), and small GTPases (RAC1 and RALB) were assessed. The findings suggest that variants in the PLCB2 gene are significantly associated with S-ECC and that genetic variants of other signaling components are associated with the composition of the oral microbiota.

2. Results

2.1. Association between Genetic Variants and S-ECC

Eighty-eight caries-free children (45.34 ± 14.69 months old, 43 females, 45 males) and 88 children with S-ECC (44.88 ± 11.84 months old, 51 females, 37 males) were recruited. No significant differences in age and sex were identified between the groups (p > 0.05).
After correcting for multiple testing, the case–control allelic association analysis showed that three variants in PLCB2 were significantly associated with S-ECC (adj. p < 0.05, Table 1). Logistic regression analysis, adjusting for sex and age, showed that children with PLCB2 single nucleotide polymorphisms (SNPs) rs2305645 (T, OR = 0.27, 95% CI = 0.14–0.51), rs1869901 (G, OR = 0.34, 95% CI = 0.20–0.58), and rs2305649 (G, OR = 0.43, 95% CI = 0.26–0.71) had lower odds of developing S-ECC (i.e., they are protective SNPs; adj. p < 0.05, Figure 1).

2.2. Association between Bacterial, Fungal Microbiota and Host Genetic Variants

The results of the differential abundance analysis between the S-ECC and caries-free groups have been published previously [7]. Results from the linear regression model showed that eleven variants were positively or negatively correlated with the abundance of bacterial or fungal taxa. No significant associations were observed between the host variants and alpha (within samples) Shannon diversity (Adj. p > 0.05). None of the three PLCB2 SNPs associated with S-ECC significantly correlated with microbial species after adjustment for multiple testing. The significant associations are shown in Table 2.

3. Discussion

This work builds on the previous study assessing the role of variants in genes encoding taste receptors, ion channels, and the oral microbiome in S-ECC risk/protection [7]. It was previously identified that polymorphisms in taste-related genes are associated with S-ECC and that raised the question of whether variants in genes encoding downstream taste signaling components and other proteins involved in GPCR signaling could also play a role.
PLCβ2 is a critical component of the taste signaling cascade. A recent study looked at the expression of PLCβ2 in oral cells, using published data from single-cell RNA-seq studies, and reported that PLCβ2 was expressed in almost all endothelial, epithelial, fibroblast, and immune cell subsets evaluated [5]. PLCβ2 is responsible for the generation of the second messenger inositol-1,4,5—triphosphate (IP3), which activates IP3 receptors in the endoplasmic reticulum (ER), releasing Ca2+ into the cytoplasm [27]. Therefore, mutations in PLCB2 can have an important effect on Ca2+-mediated physiological functions. The findings from this study suggest that PLCB2 rs2305645 (T), rs1869901 (G), and rs2305649 (G) are protective SNPs against S-ECC. Interestingly, rs1869901 has been associated with autism spectrum disorder [28].
In this study, the association between mutations in downstream taste signaling components and the oral fungal and bacterial community composition was also evaluated. The results showed that variants in GNAQ, GNAS, GNAT3, RAC1, RALB, and PLCB2 were correlated with the relative abundances of bacterial taxa, while two RAC1 SNPs were correlated with the relative abundance of fungi. Although none of the three PLCB2 variants that were associated with S-ECC were correlated with oral microbial composition, another PLCB2 variant (rs72731486) was negatively associated with an unclassified Streptococcus species. An unclassified fungus from the order Malasseziales was previously associated with healthy dental plaque [7] and here it was shown to be negatively correlated with the RAC1 rs3729790 variant.
Though the canonical taste signaling pathway suggests that taste receptors activate Gαgustducin, the co-expression of multiple Gα subunits in taste buds have suggested that other Gα subunits may be involved in taste sensation [22]. Furthermore, recent studies suggest that taste receptors are involved in multiple functions, unrelated to taste sensation, in which other Gα subunits may be involved. For instance, a study showed that the T2Rs expressed in human airway smooth muscle, where their activation leads to relaxation and bronchodilation, couple to Gαi instead of Gαgustducin [29]. Therefore, it is possible that variants in the genes encoding Gα subunits may affect the taste receptor signal transduction and other aspects of oral health.
Variants in small G proteins may affect receptors that mediate microbial recognition. RalB is a small GTPase, a member of the Ras GTPase superfamily. It is activated by calcium/calmodulin interactions and is involved in cytoskeleton rearrangement and vesicle trafficking. RalB has been implicated in several human cancers such as oral squamous cell carcinoma [30]. Interestingly, RalB has been shown to play a role in innate immunity, helping to trigger innate immune pathways after activation of Toll-like receptors (TLRs) by viruses in human epithelial cells [31]. Furthermore, RalB activation after microbial stimuli in macrophages can induce autophagy, a key process for the clearance of intracellular pathogens [32]. This agrees with our finding that RALB variants may be associated with the composition of the oral microbiome. Rac1 (Ras-related C3 botulinum toxin substrate 1) is also a small G protein and it belongs to the family of Rho GTPases. It plays a role in various cellular functions such as actin cytoskeletal reorganization, cell cycle regulation, and movement [26,33]. Rac1 is activated by a variety of receptors, including GPCRs. Recent studies from our group showed that there is a link between Rac1 and T2Rs [26,34]. This link may involve the role of Rac1 in actin cytoskeletal reorganization and T2R-mediated internalization of Gram-positive bacteria, as well as the T2R4- and G protein-dependent inhibitory effect of Rac1 activity by the known T2R agonist quinine [26,34]. It is possible that variants in RAC1 could affect the T2R-mediated internalization of microbes, which would justify the correlation of RAC1 variants with the abundance of oral bacteria and fungi.
It is important to note that this study is not free of limitations. There was no information about some of the potentially confounding variables for the association observed, such as lifestyle and socioeconomic status. Despite having a small sample size, significant associations were detected. This is the first study to suggest an association between mutations in PLCB2 and decreased odds of S-ECC. Further work is needed to characterize this association and to investigate the correlations between variants in genes encoding G proteins and the composition of bacterial and fungal communities in dental plaque.

4. Materials and Methods

4.1. Study Design

A total of 176 children younger than 72 months of age with S-ECC or caries-free controls were recruited in Winnipeg, MB, Canada. The inclusion criteria were children younger than 72 months of age; child has S-ECC (American Academy of Pediatric Dentistry definition) [35], or is caries-free (dmft index = 0, i.e., no decayed, missing, or filled primary tooth surface; no incipient lesions). The exclusion criteria were age (>72 months); child with caries but who does not satisfy the case definition of S-ECC; current use of antibiotics. Further details about the study population and methods used for the evaluation of variants in candidate genes and association between host variants and the oral microbiome have been described previously [7].
Supragingival plaque and oral swab samples were collected for the study of the plaque microbiome (bacteriome and mycobiome) and the genetic variants in the seven genes involved in taste and calcium signaling, respectively (Table 3). The samples were stored in RNAprotect solution (Qiagen, Hilden, Germany) and kept at −80 °C until further analysis.
Written informed consent was obtained from all parents or legal caregivers and the study protocol was approved by the University of Manitoba’s Health Research Ethics Board (HREB HS23754-393 H2020:150) and the Misericordia Health Centre.

4.2. DNA Sequencing and Data Analysis

The paired-end targeted V4-16S rRNA and ITS1 rRNA amplicon sequencing (MiSeq PE250, Illumina Inc., San Diego, CA, USA) was performed by Génome Québec Innovation Centre (Montréal, Canada), using DNA extracted from supragingival plaque samples. The data were analyzed using the QIIME2 (v2018.11) pipeline [36] and various R packages (“qiime2R” v0.99.13, “phyloseq” v1.30.0, “DESeq2” v1.26.0), as previously described [7]. The negative binomial Wald test implemented in the “DESeq2” R package [37] was used for differential abundance analysis, controlling the false discovery rate (FDR) for multiple comparisons, and adjusting for sex and batch effect (library preparation, sequencing runs, DNA extraction batches). The Shannon diversity index was calculated using the R package “phyloseq” (version 1.30.0). It was used to detect associations between host genetics and microbial diversity.
The paired-end targeted sequencing of the seven genes (Table 3), using DNA extracted from oral swabs, was performed by Génome Québec Innovation Centre (Montréal, Canada). The Fluidigm Access Array technology (Fluidigm, South San Francisco, CA, USA) was used for library preparation and the NovaSeq6000 SP PE250 (Illumina Inc., San Diego, CA, USA) platform was used for sequencing. The primers used are listed in Table S1.
The genetic sequence data were analyzed using the GATK (v4.2.0, Broad Institute) best practices pipeline, Picardtools (v2.25.0), VCFtools (v0.1.15), and PLINK (v1.9) [38,39,40]. A total of 662 genetic variants (SNPs and INDELS) were called. Variants with >5% missing genotype call rate, mean sequencing depth < 10x, significantly different (p < 0.00001) missing data rate between cases and controls, Hardy–Weinberg equilibrium (HWE) p  <  0.05 in controls, and MAF < 0.01 were removed. Samples with a genotype failure rate > 0.5 were also removed. After quality control, 53 variants and 174 samples were included in downstream analyses. The filtered variants were then annotated using SnpEff 5.0e [41]. A case–control allelic association analysis was used to evaluate the association between host genetic variants and S-ECC using PLINK. The allelic chi-squared (χ2) test was used to identify significant differences in allele frequencies. Regression models were further used to analyze the associations between host genetic variants and S-ECC, and host genetic variants and the oral microbiome with an additive genetic model using PLINK, adjusting for age, sex, and microbiome sequencing batch [7].

Supplementary Materials

The supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/ijms24010081/s1.

Author Contributions

Conceptualization, V.C.d.J., R.J.S., and P.C.; methodology, V.C.d.J., P.H., R.J.S., and P.C.; formal analysis, V.C.d.J.; investigation, V.C.d.J., B.-A.M., and R.J.S.; resources, B.-A.M., R.J.S., and P.C.; data curation, V.C.d.J., and B.-A.M.; writing—original draft, V.C.d.J.; writing—review and editing, V.C.d.J., P.H., R.J.S., and P.C.; visualization, V.C.d.J.; supervision, R.J.S., and P.C.; project administration, B.-A.M.; funding acquisition, P.H., R.J.S., and P.C. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by an operating project grant from the Canadian Institute of Health Research (Grant No. PJT-159731).

Institutional Review Board Statement

The study was conducted according to the guidelines of the Declaration of Helsinki and approved by the University of Manitoba’s Health Research Ethics Board (HREB, HS23754-H2020:150).

Informed Consent Statement

Informed consent was obtained from all parents or legal caregivers of the subjects involved in the study.

Data Availability Statement

De-identified raw 16S and ITS1 rRNA gene sequences derived from the supragingival plaque samples are deposited at NCBI Sequence Read Archive (SRA) Repository (accession number PRJNA555320).

Acknowledgments

We thank the participants, the parents and caregivers, the recruitment team, and the Misericordia Health Centre.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Lagerström, M.C.; Schiöth, H.B. Structural Diversity of G Protein-Coupled Receptors and Significance for Drug Discovery. Nat. Rev. Drug Discov. 2008, 7, 339–357. [Google Scholar] [CrossRef] [PubMed]
  2. Jaggupilli, A.; Howard, R.; Upadhyaya, J.D.; Bhullar, R.P.; Chelikani, P. Bitter Taste Receptors: Novel Insights into the Biochemistry and Pharmacology. Int. J. Biochem. Cell Biol. 2016, 77, 184–196. [Google Scholar] [CrossRef] [PubMed]
  3. Liu, L.; Li, S.; Zheng, J.; Bu, T.; He, G.; Wu, J. Safety Considerations on Food Protein-Derived Bioactive Peptides. Trends Food Sci. Technol. 2020, 96, 199–207. [Google Scholar] [CrossRef]
  4. Shaik, F.A.; Singh, N.; Arakawa, M.; Duan, K.; Bhullar, R.P.; Chelikani, P. Bitter Taste Receptors: Extraoral Roles in Pathophysiology. Int. J. Biochem. Cell Biol. 2016, 77, 197–204. [Google Scholar] [CrossRef]
  5. Xi, R.; Zheng, X.; Tizzano, M. Role of Taste Receptors in Innate Immunity and Oral Health. J. Dent. Res. 2022, 101, 002203452210779. [Google Scholar] [CrossRef]
  6. Pchitskaya, E.; Popugaeva, E.; Bezprozvanny, I. Calcium Signaling and Molecular Mechanisms Underlying Neurodegenerative Diseases. Cell Calcium 2018, 70, 87–94. [Google Scholar] [CrossRef]
  7. de Jesus, V.C.; Mittermuller, B.; Hu, P.; Schroth, R.J.; Chelikani, P. Genetic Variants in Taste Genes Play a Role in Oral Microbial Composition and Severe Early Childhood Caries. iScience 2022, 25, 105489. [Google Scholar] [CrossRef]
  8. Luo, J.; Sun, P.; Siwko, S.; Liu, M.; Xiao, J. The Role of GPCRs in Bone Diseases and Dysfunctions. Bone Res. 2019, 7, 19. [Google Scholar] [CrossRef] [Green Version]
  9. Parry, D.A.; Smith, C.E.L.; El-Sayed, W.; Poulter, J.A.; Shore, R.C.; Logan, C.V.; Mogi, C.; Sato, K.; Okajima, F.; Harada, A.; et al. Mutations in the PH-Sensing G-Protein-Coupled Receptor GPR68 Cause Amelogenesis Imperfecta. Am. J. Hum. Genet. 2016, 99, 984–990. [Google Scholar] [CrossRef] [Green Version]
  10. Seymen, F.; Zhang, H.; Kasimoglu, Y.; Koruyucu, M.; Simmer, J.P.; Hu, J.C.-C.; Kim, J.-W. Novel Mutations in GPR68 and SLC24A4 Cause Hypomaturation Amelogenesis Imperfecta. J. Pers. Med. 2022, 12, 13. [Google Scholar] [CrossRef]
  11. Zhang, Y.; Kim, J.-Y.; Horst, O.; Nakano, Y.; Zhu, L.; Radlanski, R.J.; Ho, S.; Besten, P.K.D. Fluorosed Mouse Ameloblasts Have Increased SATB1 Retention and Gαq Activity. PLoS ONE 2014, 9, e103994. [Google Scholar] [CrossRef] [PubMed]
  12. Zhang, L.; Sugamori, K.S.; Claridge, C.; dela Cruz, A.; Grynpas, M.D.; Mitchell, J. Overexpression of GαS in Murine Osteoblasts In Vivo Leads to Increased Bone Mass and Decreased Bone Quality. J. Bone Miner. Res. 2017, 32, 2171–2181. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Simon, B.R.; Learman, B.S.; Parlee, S.D.; Scheller, E.L.; Mori, H.; Cawthorn, W.P.; Ning, X.; Krishnan, V.; Ma, Y.L.; Tyrberg, B.; et al. Sweet Taste Receptor Deficient Mice Have Decreased Adiposity and Increased Bone Mass. PLoS ONE 2014, 9, e86454. [Google Scholar] [CrossRef]
  14. Yu, T.; Klein, O.D. Molecular and Cellular Mechanisms of Tooth Development, Homeostasis and Repair. Development 2020, 147, dev184754. [Google Scholar] [CrossRef] [PubMed]
  15. Cruvinel, V.R.N.; Gravina, D.B.L.; Azevedo, T.D.P.L.; de Rezende, C.S.; Bezerra, A.C.B.; Toledo, O.A. de Prevalence of Enamel Defects and Associated Risk Factors in Both Dentitions in Preterm and Full Term Born Children. J. Appl. Oral Sci. 2012, 20, 310–317. [Google Scholar] [CrossRef] [Green Version]
  16. Oliveira, A.F.B.; Chaves, A.M.B.; Rosenblatt, A. The Influence of Enamel Defects on the Development of Early Childhood Caries in a Population with Low Socioeconomic Status: A Longitudinal Study. Caries Res. 2006, 40, 296–302. [Google Scholar] [CrossRef] [PubMed]
  17. Pierce, A.; Singh, S.; Lee, J.; Grant, C.; de Jesus, V.C.; Schroth, R.J. The Burden of Early Childhood Caries in Canadian Children and Associated Risk Factors. Front. Public Health 2019, 7, 328. [Google Scholar] [CrossRef] [Green Version]
  18. Wendell, S.; Wang, X.; Brown, M.; Cooper, M.E.; DeSensi, R.S.; Weyant, R.J.; Crout, R.; McNeil, D.W.; Marazita, M.L. Taste Genes Associated with Dental Caries. J. Dent. Res. 2010, 89, 1198–1202. [Google Scholar] [CrossRef]
  19. Caicedo, A.; Pereira, E.; Margolskee, R.F.; Roper, S.D. Role of the G-Protein Subunit Alpha-Gustducin in Taste Cell Responses to Bitter Stimuli. J. Neurosci. 2003, 23, 9947–9952. [Google Scholar] [CrossRef] [Green Version]
  20. Singh, N.; Bhagirath, A.Y.; Medapati, M.R.; de Jesus, V.C.; Chelikani, P. Chemistry and Biological Mechanisms of Peptides that Modulate Taste. In Food Proteins and Peptides: Emerging Biofunctions, Food and Biomaterial Applications; Royal Society of Chemistry: London, UK, 2021; p. 31. [Google Scholar]
  21. Lee, A.; Owyang, C. Sugars, Sweet Taste Receptors, and Brain Responses. Nutrients 2017, 9, 653. [Google Scholar] [CrossRef]
  22. von Molitor, E.; Riedel, K.; Krohn, M.; Hafner, M.; Rudolf, R.; Cesetti, T. Sweet Taste Is Complex: Signaling Cascades and Circuits Involved in Sweet Sensation. Front. Hum. Neurosci. 2021, 15, 667709. [Google Scholar] [CrossRef] [PubMed]
  23. Wong, G.T.; Gannon, K.S.; Margolskee, R.F. Transduction of Bitter and Sweet Taste by Gustducin. Nature 1996, 381, 796–800. [Google Scholar] [CrossRef] [PubMed]
  24. Clough, R.R.; Sidhu, R.S.; Bhullar, R.P. Calmodulin Binds RalA and RalB and Is Required for the Thrombin-Induced Activation of Ral in Human Platelets. J. Biol. Chem. 2002, 277, 28972–28980. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Elsaraj, S.M.; Bhullar, R.P. Regulation of Platelet Rac1 and Cdc42 Activation through Interaction with Calmodulin. Biochim. Et Biophys. Acta (BBA)-Mol. Cell Res. 2008, 1783, 770–778. [Google Scholar] [CrossRef] [Green Version]
  26. Sidhu, C.; Jaggupilli, A.; Chelikani, P.; Bhullar, R.P. Regulation of Rac1 GTPase Activity by Quinine through G-Protein and Bitter Taste Receptor T2R4. Mol. Cell. Biochem. 2017, 426, 129–136. [Google Scholar] [CrossRef]
  27. Behrens, M.; Meyerhof, W. Signaling in the Chemosensory Systems. Cell. Mol. Life Sci. 2006, 63, 1501–1509. [Google Scholar] [CrossRef]
  28. Skafidas, E.; Testa, R.; Zantomio, D.; Chana, G.; Everall, I.P.; Pantelis, C. Predicting the Diagnosis of Autism Spectrum Disorder Using Gene Pathway Analysis. Mol. Psychiatry 2014, 19, 504–510. [Google Scholar] [CrossRef] [Green Version]
  29. Kim, D.; Woo, J.A.; Geffken, E.; An, S.S.; Liggett, S.B. Coupling of Airway Smooth Muscle Bitter Taste Receptors to Intracellular Signaling and Relaxation Is via Gαi1,2,3. Am. J. Respir. Cell Mol. Biol. 2017, 56, 762–771. [Google Scholar] [CrossRef]
  30. Zinatizadeh, M.R.; Zarandi, P.K.; Keshavarz-Fathi, M.; Yousefi, M.H.; Rezaei, N. The Role of Ral Signaling and Post Translational Modifications (PTMs) of Ras in Cancer. Genome Instab. Dis. 2022, 3, 22–32. [Google Scholar] [CrossRef]
  31. Chien, Y.; Kim, S.; Bumeister, R.; Loo, Y.M.; Kwon, S.W.; Johnson, C.L.; Balakireva, M.G.; Romeo, Y.; Kopelovich, L.; Gale, M.; et al. RalB GTPase-Mediated Activation of the IκB Family Kinase TBK1 Couples Innate Immune Signaling to Tumor Cell Survival. Cell 2006, 127, 157–170. [Google Scholar] [CrossRef]
  32. Shi, C.S.; Shenderov, K.; Huang, N.N.; Kabat, J.; Abu-Asab, M.; Fitzgerald, K.A.; Sher, A.; Kehrl, J.H. Activation of Autophagy by Inflammatory Signals Limits IL-1β Production by Targeting Ubiquitinated Inflammasomes for Destruction. Nat. Immunol. 2012, 13, 255–263. [Google Scholar] [CrossRef] [PubMed]
  33. Jaffe, A.B.; Hall, A. RHO GTPASES: Biochemistry and Biology. Annu. Rev. Cell Dev. Biol. 2005, 21, 247–269. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Medapati, M.R.; Bhagirath, A.Y.; Singh, N.; Schroth, R.J.; Bhullar, R.P.; Duan, K.; Chelikani, P. Bitter Taste Receptor T2R14 Modulates Gram-Positive Bacterial Internalization and Survival in Gingival Epithelial Cells. IJMS 2021, 22, 9920. [Google Scholar] [CrossRef]
  35. AAPD. Policy on Early Childhood Caries (ECC): Classifications, Consequences, and Preventive Strategies. Ref. Man. Pediatr. Dent. Chic. Ill. Am. Acad. Pediatr. Dent. 2017, 79–81. [Google Scholar]
  36. Bolyen, E.; Rideout, J.R.; Dillon, M.R.; Bokulich, N.A.; Abnet, C.C.; Al-Ghalith, G.A.; Alexander, H.; Alm, E.J.; Arumugam, M.; Asnicar, F.; et al. Reproducible, Interactive, Scalable and Extensible Microbiome Data Science Using QIIME 2. Nat. Biotechnol. 2019, 37, 852–857. [Google Scholar] [CrossRef]
  37. Love, M.I.; Huber, W.; Anders, S. Moderated Estimation of Fold Change and Dispersion for RNA-Seq Data with DESeq2. Genome Biol. 2014, 15, 550. [Google Scholar] [CrossRef] [Green Version]
  38. Danecek, P.; Auton, A.; Abecasis, G.; Albers, C.A.; Banks, E.; DePristo, M.A.; Handsaker, R.E.; Lunter, G.; Marth, G.T.; Sherry, S.T.; et al. The Variant Call Format and VCFtools. Bioinformatics 2011, 27, 2156–2158. [Google Scholar] [CrossRef] [PubMed]
  39. Purcell, S.; Neale, B.; Todd-Brown, K.; Thomas, L.; Ferreira, M.A.R.; Bender, D.; Maller, J.; Sklar, P.; de Bakker, P.I.W.; Daly, M.J.; et al. PLINK: A Tool Set for Whole-Genome Association and Population-Based Linkage Analyses. Am. J. Hum. Genet. 2007, 81, 559–575. [Google Scholar] [CrossRef] [Green Version]
  40. Van der Auwera, G.; O’Connor, B.D. Genomics in the Cloud: Using Docker, GATK, and WDL in Terra, 1st ed.; O’Reilly Media: Sebastopol, CA, USA, 2020. [Google Scholar]
  41. Cingolani, P.; Patel, V.M.; Coon, M.; Nguyen, T.; Land, S.J.; Ruden, D.M.; Lu, X. Using Drosophila Melanogaster as a Model for Genotoxic Chemical Mutational Studies with a New Program, SnpSift. Front. Genet. 2012, 3, 35. [Google Scholar] [CrossRef]
Figure 1. Forest plot for the results of logistic regression. Effect sizes are expressed on a log OR scale. Adj. p, p corrected for multiple testing by Bonferroni adjustment. Adjusted p < 0.05 were considered statistically significant.
Figure 1. Forest plot for the results of logistic regression. Effect sizes are expressed on a log OR scale. Adj. p, p corrected for multiple testing by Bonferroni adjustment. Adjusted p < 0.05 were considered statistically significant.
Ijms 24 00081 g001
Table 1. Allelic association of genetic variants with severe early childhood caries (S-ECC).
Table 1. Allelic association of genetic variants with severe early childhood caries (S-ECC).
GeneVariantsLocation (GRCh38)Effect AlleleFrequency of Effect Alleleχ2pORAdj. p *Type of Variant
Cases (S-ECC)Controls (Caries-Free)
PLCB2rs2305645chr15: 40303364T0.090.2719.051.28 × 10−50.270.0007Intron
rs1869901chr15: 40303426G0.170.3818.431.77 × 10−50.340.0009Intron
rs2305649chr15: 40297629G0.160.3413.960.00020.380.0099Intron
Note: * Adjusted p (corrected for multiple testing by Bonferroni adjustment test) less than 0.05 were considered statistically significant. Only significant associations are shown. χ2, basic allelic test Chi-square. The basic allelic test compares frequencies of alleles in cases versus controls. Chr, chromosome. OR, odds ratio.
Table 2. Association between host genetic variants and oral bacterial and fungal taxa.
Table 2. Association between host genetic variants and oral bacterial and fungal taxa.
GeneVariantLocationAlleleBETAAdj. pTaxa
BACTERIA (n = 174)
GNAQc.736-12T > Cchr9:77728679G−1.2410.0029Capnocytophaga sputigena
GNAQrs1478186975chr9:77728678G−2.1920.0311Actinomyces gerencseriae
RAC1rs1051504128chr7:6387120A6.9130.0095Unclassified bacteria
RAC1rs836478chr7:6392059C0.5340.0443Lachnospiraceae [G-3] bacterium HMT 100
GNASrs3730173chr20:58909879T−3.140.0208Lachnoanaerobaculum saburreum
RALBc.356A > Gchr2:120289612G−3.140.0208Lachnoanaerobaculum saburreum
RALBrs11545293chr2:120278757A0.72320.0311Bergeyella sp. HMT 907
PLCB2rs72731486chr15:40299054A−0.58030.0325Genus Streptococcus
GNAT3rs6975345chr7:80494683C2.7550.0351Prevotella salivae
FUNGI (n = 155)
RAC1rs836478chr7:6392059C2.5190.0324Genus Alternaria
RAC1rs3729790chr7:6387323A−2.4940.0329Order Malasseziales
Beta, regression coefficient. Adj. p < 0.05 are corrected for multiple testing by Bonferroni adjustment test.
Table 3. List of sequenced genes.
Table 3. List of sequenced genes.
Gene SymbolGene NameRefSeq IDLocationFunction
GNASG protein
subunit alpha s
NM_001077489Chr20: 58,891,364–58,911,192It encodes the guanine nucleotide-binding protein Gαs, which is involved in activation of adenylyl cyclase (AC), the enzyme that synthesizes cyclic adenosine monophosphate (cAMP) from adenosine triphosphate (ATP), and a variety of cellular responses. With relevance to dental and oral health, mutations in this gene are linked to bone defects.
GNAI2G protein
subunit alpha i2
NM_002070Chr3: 50,236,204–50,259,362It encodes the Gαi2 subunit, which is involved in hormonal regulation of AC. There is evidence of interaction between Gαi2 and bitter taste receptors (T2Rs).
GNAQG protein
subunit alpha q
NM_002072Chr9: 77,716,097–78,031,811It encodes the Gαq subunit, which couples GPCRs and PLCβ.
GNAT3G protein
subunit alpha
transducin 3
NM_001102386Chr7: 80,458,635–80,512,064The Gαgustducin encoded by this gene binds to taste receptors and is involved in the canonical taste signaling pathway.
PLCB2Phospholipase
C beta 2
NM_004573Chr15: 40,287,909–40,307,935It encodes PLCβ2, which catalyzes the hydrolysis of PIP2 to IP3, which elicits Ca2+ release from internal stores. PLCβ2 is involved in the canonical taste signaling pathway.
RAC1Rac family
small GTPase 1
NM_006908Chr7: 6,374,527–6,403,967It encodes a GTPase belonging to the RAS superfamily of small G proteins. Members of this superfamily have been shown to regulate a broad number of cellular events such as cytoskeletal dynamics and have a possible link with T2Rs.
RALBRAS-like proto-oncogene BNM_002881Chr2: 120,252,852–120,294,710It encodes a GTP-binding protein that is a member of the small GTPase superfamily and Ras family of proteins. It is involved in innate immunity and tumor growth.
Note: The starting and ending nucleotide positions are from the December 2013 (GRCh38/hg38) human assembly (http://genome.ucsc.edu (accessed on 14 March 2022)).
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

de Jesus, V.C.; Mittermuller, B.-A.; Hu, P.; Schroth, R.J.; Chelikani, P. Association between Downstream Taste Signaling Genes, Oral Microbiome, and Severe Early Childhood Caries. Int. J. Mol. Sci. 2023, 24, 81. https://doi.org/10.3390/ijms24010081

AMA Style

de Jesus VC, Mittermuller B-A, Hu P, Schroth RJ, Chelikani P. Association between Downstream Taste Signaling Genes, Oral Microbiome, and Severe Early Childhood Caries. International Journal of Molecular Sciences. 2023; 24(1):81. https://doi.org/10.3390/ijms24010081

Chicago/Turabian Style

de Jesus, Vivianne Cruz, Betty-Anne Mittermuller, Pingzhao Hu, Robert J. Schroth, and Prashen Chelikani. 2023. "Association between Downstream Taste Signaling Genes, Oral Microbiome, and Severe Early Childhood Caries" International Journal of Molecular Sciences 24, no. 1: 81. https://doi.org/10.3390/ijms24010081

APA Style

de Jesus, V. C., Mittermuller, B. -A., Hu, P., Schroth, R. J., & Chelikani, P. (2023). Association between Downstream Taste Signaling Genes, Oral Microbiome, and Severe Early Childhood Caries. International Journal of Molecular Sciences, 24(1), 81. https://doi.org/10.3390/ijms24010081

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop