Gene set-based analysis of genome-wide association study (GWAS) data has recently emerged as a useful approach to examine the joint effects of multiple risk loci in complex human diseases or phenotypes. Dental caries is a common, chronic, and complex disease leading to a decrease in quality of life worldwide. In this study, we applied the approaches of gene set enrichment analysis to a major dental caries GWAS dataset, which consists of 537 cases and 605 controls. Using four complementary gene set analysis methods, we analyzed 1331 Gene Ontology (GO) terms collected from the Molecular Signatures Database (MSigDB). Setting false discovery rate (FDR) threshold as 0.05, we identified 13 significantly associated GO terms. Additionally, 17 terms were further included as marginally associated because they were top ranked by each method, although their FDR is higher than 0.05. In total, we identified 30 promising GO terms, including ‘Sphingoid metabolic process,’ ‘Ubiquitin protein ligase activity,’ ‘Regulation of cytokine secretion,’ and ‘Ceramide metabolic process.’ These GO terms encompass broad functions that potentially interact and contribute to the oral immune response related to caries development, which have not been reported in the standard single marker based analysis. Collectively, our gene set enrichment analysis provided complementary insights into the molecular mechanisms and polygenic interactions in dental caries, revealing promising association signals that could not be detected through single marker analysis of GWAS data.
Citation: Wang Q, Jia P, Cuenco KT, Zeng Z, Feingold E, Marazita ML, et al. (2013) Association Signals Unveiled by a Comprehensive Gene Set Enrichment Analysis of Dental Caries Genome-Wide Association Studies. PLoS ONE 8(8): e72653. doi:10.1371/journal.pone.0072653
Editor: Kai Wang, University of Southern California, United States of America
Received: April 30, 2013; Accepted: July 11, 2013; Published: August 14, 2013
Copyright: © 2013 Wang et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Funding: This work was supported by grants from National Institute of Dental and Craniofacial Research, National Institutes of Health (R03DE022093, R01DE014899, and U01DE018903). Dental caries GWAS dataset use was granted through the NIH GWAS Data Repository (dbGaP, accession number phs000095.v1.p1) Data Access Committee [request #7562-2] to ZZ. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Competing interests: The authors have declared that no competing interests exist.
Dental caries (also known as tooth decay or a cavity) is simply defined as a procedure that causes destruction and demineralization of hard tooth tissues such as enamel, dentin, and cementum. It is a highly pervasive chronic disease whose etiology is complex and multifactorial, with contributions from numerous factors, including microbial flora, salivary flow and composition, and fluoride exposure, among others. There has been increasing evidence of genetic components contributing to caries susceptibility , , , , . Benefiting from high-throughput genotyping technologies (up to a few million single nucleotide polymorphism (SNP) biomarkers on a single chip), genome-wide association studies (GWAS) have recently been employed to search for genetic susceptibility related to dental caries , , , among hundreds of other complex diseases and phenotypes . These dental caries GWA studies identified several loci and genes, such as ACTN2, LYZL2, and AJAP1 , . In these GWA studies, the statistical analyses of association signals are typically conducted for single markers, limiting the power to identify potential truly associated genes that may have been missed due to the multiple test adjustment necessary to control the false discovery rate (FDR). Recently, interrogating the joint effects of multiple risk loci or genes through the gene set-based analysis of GWAS data has become one popular complementary approach to single marker association tests . Gene set analysis of GWAS data has been successfully applied to many diseases or traits (see recent reviews , ), including schizophrenia , major depressive disorder , type II diabetes , , , Crohn’s disease , and several types of cancer , , , . However, to our knowledge, no such studies have been reported for gene set analysis of association data for human caries to date. In this work, we performed a comprehensive gene set analysis of GWAS data for dental caries, aiming to broaden our understanding of the role of interactions between genes for this worldwide disease.
Over the past several years, many gene set analysis methods have been proposed, which were extensively summarized in a recent review . These methods address two different null hypotheses on their tests of disease associations: 1) competitive null hypothesis (Q1), which tests whether the genes within a gene set show the same association magnitude compared to the genes outside the gene set; and 2) self-contained null hypothesis (Q2), which tests whether the genes within a gene set are associated with the disease phenotype. When the real causal genes are included in only a few gene sets, the two tests may have similar results. Nevertheless, the competitive tests may be less powerful when the causal genes are shared by multiple gene sets. Apart from the difference in the null hypotheses tested, each method has its own strengths and limitations, and no single proposed strategy outperforms all the others .
In this study, we employed four representative methods to conduct gene set enrichment analyses for dental caries, among which two perform competitive tests (Association List Go AnnoTatOR (ALIGATOR)  and GenGen ) and the other two perform self-contained tests (SNP ratio test (SRT)  and the mixed model ). The GWAS data was collected from a recent dental caries association study , and the Gene Ontology (GO) annotation database  was the source for candidate gene sets. Our study integrated the results from different approaches and reported 13 significantly associated and 17 marginally associated GO terms. To our knowledge, this is the first comprehensive gene set analysis for dental caries, or generally for dental health, to date. Our findings provide biological insights into the potential molecular mechanisms underlying dental caries, which helps to improve our understanding of dental caries beyond the single marker level.
Materials and Methods
We retrieved the dental caries GWAS data  from dbGaP (http://www.ncbi.nlm.nih.gov/gap) through approved access (dbGaP accession number: phs000095.v1.p1). A total of 4,020 individuals in this dataset have both genotype and phenotype data. We focused on the phenotype of ‘total primary tooth caries.’ In this dataset, the total primary tooth caries with white spots is described by the continuous variable ‘Prim_d1ft’ and the dichotomized variable ‘CAT1_PRIM_D1FT.’ By definition, individuals with disease are those with Prim_d1ft ≥ 1 (CAT1_PRIM_D1FT = 1) and controls are those with Prim_d1ft = 0 (CAT1_PRIM_D1FT = 0). Subjects who were between 3 and 12 years old at the time of dental exam were included. A total of 537 cases and 605 controls, among which there are 588 males and 554 females, formed our working dataset. The samples were genotyped on the Illumina platform Human610_Quadv1_B (Illumina, Inc.). Quality checks conducted in the original study as provided by dbGaP resulted in 589,735 SNPs for the following analyses.
Gene set annotation
The Molecular Signatures Database (MSigDB)  collects annotated gene sets from multiple sources. We downloaded the GO annotation  from MSigDB (version 3.0, C5) for gene set enrichment analysis. To avoid biological functions that are too narrowly or too broadly defined, only gene sets containing ≥ 5 and ≤ 250 genes were included in the following analyses. As a result, 1,331 GO terms passed the criteria, and the average number of genes per term was 44.
Logistic regression was performed for association test of each of the 589,735 SNPs with CAT1_PRIM_D1FT using the GWAS analysis tool PLINK . The variable “age at time of dental exam” was taken as a covariate in the regression. The overall genomic inflation factor was 1.031. We denoted the test statistic of each SNP as (i = 1,2,…,L, where L is total number of SNPs) and the p-value as (a higher indicates a lower ). A SNP was mapped to a gene if it is located in the gene region or within 20 kb upstream or downstream of the gene. We applied this criterion based on the previous studies , , . The SNP-gene mapping resulted in 20,756 protein coding genes based on the human reference assembly hg18.
Gene set enrichment studies for GWAS data have been proposed for several years. However, no single strategy outperforms all the others to date. To alleviate the potential biases in different statistical algorithms, we chose four representative methods to perform the gene set enrichment analysis in this study. These methods are GenGen , ALIGATOR , SRT , and the mixed model . The first two methods are used to test competitive null hypothesis (Q1), while the others are used to test self-contained null hypothesis (Q2) . We briefly describe the methods below. More details can be found in the original publications.
GenGen  is adapted from the Gene Set Enrichment Analysis (GSEA) method  that was originally designed to analyze gene expression data. The first step of this approach is to assign each gene a significance value (j = 1,2,…,N, where N is the total number of genes) with the most significant that can be mapped to this gene. Next, all the genes are ranked in descending order of , denoted by . Third, for a given gene set S consisting of genes, an enrichment score (ES) is computed using a weighted Kolmogorov-Smirnov-like running-sum statistic as follows:
where . Finally the significance of ES(S) is evaluated using a permutation test by shuffling the labels of cases and controls so that the linkage disequilibrium (LD) structures among SNPs are conserved.
The algorithm ALIGATOR  executes a SNP-based resampling procedure, which can effectively reduce the potential biases from gene size, SNP density, and LD structure. ALIGATOR defines a set of significantly associated SNPs through a predefined cutoff (e.g., p-value = 0.05). It maps these significant SNPs to genes, which are in turn denoted as significant genes, and counts the number of significant genes for each gene set. Then, the algorithm performs a SNP-based resampling, during which SNPs are selected and mapped to genes until the number of significant genes generated by the resampling process is the same as in the original case. Resample genes are mapped to gene sets in the same way as in actual cases, and the numbers of significant genes per gene set are recorded. In our analysis, we performed resampling 10,000 times. Finally, an empirical p-value is computed for each gene set by summing the number of resampling datasets that have a higher number of significant genes than the real case.
The SNP ratio test  similarly defines a set of SNPs that are significantly associated with the disease through a predefined p-value threshold . For a gene set S, the proportion of significant SNPs is first computed as , where M is the total number of SNPs mapped to all the genes in S. Then, the p-value of is estimated using a permutation by shuffling the case/control status among samples.
The mixed model  employs a hierarchical generalized linear model for gene set analysis. For each gene set, the mixed model includes the chi-square statistic (transformed from p-value) for each SNP as the outcome variable, random gene effects, and an intercept corresponding to the overall association with disease for all SNPs in the gene set. The statistical significance of the gene set is determined based on p-value for the intercept.
Both GenGen and SRT require permutation procedure by shuffling the case/control labels to determine the significance level. We generated permutation data with the same parameters and covariates 1,000 times for the use of these two algorithms. To correct multiple comparisons, the Benjamini-Hochberg method  was used to control the FDR.
We performed gene set enrichment analyses of dental caries GWAS data using four statistical methods (GenGen, ALIGATOR, SRT, and the mixed model) and GO annotation terms as the gene set pool. Setting FDR < 0.05 as the criterion to determine the statistical significance, the mixed model identified the largest number of GO terms that are statistically associated with dental caries, i.e., a total of 9 GO terms. The GenGen method claimed 4 significant GO terms, whereas no significant results could be found by either ALIGATOR or SRT (Tables 1 and 2). Interestingly, the GO terms identified by GenGen are all related to secretion or regulation of secretion: ‘Protein secretion,’ ‘Cytokine secretion,’ ‘Regulation of protein secretion,’ and ‘Regulation of cytokine secretion.’ The mixed model identified several GO terms that are related to neural development (‘Regulation of axonogenesis,’ ‘Regulation of neurogenesis,’ ‘Axonogenesis,’ and ‘Central nervous system development’) and three GO terms that are related to ligase activities (‘Ligase activity forming carbon nitrogen bonds,’ ‘Ubiquitin protein ligase activity,’ and ‘Small conjugating protein ligase activity’).
We further examined the genes that contributed to the association of these GO terms with dental caries. Genes that contained at least one SNP with its p-value < 0.05 calculated from the GWAS dataset were defined as “contributing genes.” Table 1 shows the contributing genes for the 13 associated GO terms. Some gene sets showing similar biological functions share many contributing genes. For example, ‘Regulation of axonogenesis’ and ‘Regulation of neurogenesis’ shared seven genes, including some interesting genes such as ROBO2 and SLIT2 (see Discussion section). Notably, the gene set ‘Cell matrix junction’ that was identified by the mixed model contains gene ACTN2, which was reported in the original GWAS dataset with suggestive evidence for association, but failed to meet the genome-wide significance (p-value < 10−7) . Our finding confirmed this result based on single SNP analysis of the original GWAS data at the gene set level. To further examine whether the association of this gene set with dental caries is driven by gene ACTN2, we excluded this gene and performed the same gene set analysis using the mixed model approach. Interestingly, the gene set ‘Cell matrix junction’ remained significant (FDR = 0.007) even without the gene ACTN2, indicating that there are additional informative genes in this gene set that contributed to the association.
Although ALIGATOR and SRT reported no significant GO terms under the criterion FDR < 0.05, several gene sets had reasonably low p-values before multiple testing correction and underwent further investigation. The high FDR values are likely due to the inherent characteristics of the algorithms used for these approaches, which is a phenomenon noticed in previous studies . To better explore the results of ALIGATOR and SRT, we adopted the strategy proposed in . Specifically, among the four methods we applied, the largest number of gene sets at FDR < 0.05 was 9, as reported by the mixed model approach. Therefore, we accordingly selected the top 9 gene sets ranked by their raw p-values and denoted them as candidate gene sets for each of the corresponding approaches (Table 2). Note that all the gene sets selected in this way have nominally significant p-values (within a range of 0 − 0.017). Among them, one gene set was identified by three methods, and four other gene sets were identified by two methods. Interestingly, ALIGATOR reported all five of the GO terms that can be identified by at least two strategies. Of especial note, the gene set ‘Sphingoid metabolic process’ was ranked as the most significant by the results from both ALIGATOR and SRT. The four other gene sets included ‘Ligase activity forming carbon nitrogen bonds,’ which was discovered by ALIGATOR, GenGen, and the mixed model, ‘Ubiquitin protein ligase activity’ by ALIGATOR and the mixed model, ‘Regulation of cytokine secretion’ by ALIGATOR and GenGen, and ‘Ceramide metabolic process’ by ALIGATOR and SRT. Note that GenGen and ALIGATOR are methods to investigate the competitive null hypothesis (Q1), and SRT and the mixed model are used for the self-contained null hypothesis (Q2). We saw from Table 2 that four gene sets were identified for both Q1 and Q2: ‘Sphingoid metabolic process,’ ‘Ligase activity forming carbon nitrogen bonds,’ ‘Ubiquitin protein ligase activity,’ and ‘Ceramide metabolic process.’ In total, we listed 30 top GO terms in Table 2.
In addition, we examined the set sizes (i.e., the number of genes) of the gene sets identified by each method. The sizes of the gene sets identified by the mixed model were greater than that of other methods. The median value of set sizes for the top 9 GO terms identified by the mixed model was 43, whereas the corresponding numbers were 18 for GenGen, 24 for ALIGATOR, and 23 for SRT, respectively. Meanwhile, the SNP density (represented by median number of SNPs per gene) in the GO terms discovered by four approaches are similar, i.e., 13, 12, 13, and 12 for the mixed model, GenGen, ALIGATOR, and SRT, respectively.
We further examined the association signals of the genes that resided in the 30 top GO terms reported by four different methods. A gene was considered nominally significant if it contained at least one SNP with its p-value < 0.05. Using this criterion, we found 383 nominally significant genes, among which 36 were involved in at least 4 GO terms (Table 3). The complete description of all the 383 significant genes was shown in the supplementary materials (Table S1). We used the Ingenuity Pathway Analysis (IPA, http://www.ingenuity.com, accessed in January, 2013) software to further investigate the phenotype annotations of these nominally significant genes. We searched the IPA using “dental” as the keyword in the category of ‘Functions and Diseases’ and obtained 122 related function annotation items. Eight of the 383 nominally significant genes were found in the dental related Ingenuity annotations: PBX3, PBX1, BCOR, GLI2, SHH, DIAPH1, SOX3, and RECQL4. They are mainly related to the Ingenuity functions ‘dental development’ and ‘dental disorder’ (Table S1).Of special note, association between BCOR and pit-and-fissure surface caries has been found in a recently published GWAS in the permanent dentition . However, it failed to be detected in primary caries through the genome-wide, single-marker analysis approach .
With many GWAS datasets having been released, gene set enrichment analysis was proposed as an important and complementary approach to the traditional single marker analysis of GWAS data. Compared to single marker analysis, gene set analysis focuses more on biological functions of gene products as well as their regulation in the cellular systems. Thus, this strategy has advantages in revealing potential molecular mechanisms underlying diseases or traits. In addition, both real and simulation studies indicated that gene set enrichment analysis could increase the power of detecting association signals , . In this study, we conducted a comprehensive gene set analysis for dental caries GWAS data . Applying four methods (GenGen, ALIGATOR, SRT, and the mixed model), we identified 30 GO terms that were significantly or marginally associated with dental caries (Table 2). Among them, five gene sets were identified by at least two enrichment methods (i.e. ‘Ligase activity forming carbon nitrogen bonds,’ ‘Regulation of cytokine secretion,’ ‘Ceramide metabolic process,’ ‘Sphingoid metabolic process,’ and ‘Ubiquitin protein ligase activity’). While definitive roles for the gene sets cannot be identified as sufficient to cause cariogenesis, the five GO terms are plausible factors for disease. These terms encompass broad functions that potentially interact and contribute to the oral immune response to caries-related organisms. The oral environment contains bacteria that may lead to a host inflammatory response eliciting cytokines , , . This inflammatory response involves the sphingolipds, of which sphingoids and ceramides are constituent components released during the response , . Anaerobic organisms present in the oral cavity thrive under hypoxic conditions, which have been observed to stimulate cytokine production regulated by ubiquitin protein ligases , , , . While no direct action from carbon nitrogen bond ligases is identified within the immune response pathways, it is possible that they function in a parallel maintenance mechanism for the immune-related pathways.
We further investigated the 36 identified genes associated with the 30 top GO terms (Table 3) for their potential overlap with caries development. For each gene, the GeneCards (http://www.genecards.org/) entry (summaries and function) and OMIM (http://omim.org/) entry were queried to summarize gene functions. GeneCards aliases were also searched for in OMIM. A query for gene name and each alias cross-listed with “caries”, “tooth”, and “dentin” was conducted in PubMed to further assess known genetic roles related to dental caries.
Based on our gene-based literature search, five genes from either ligase activity (gene: WWP2, RNF217), neuronal development (gene: ROBO2, SLIT2), or cytokine/protein secretion (gene: INS) gene sets listed in Table 2 might be potentially associated with dental traits. Only the cytokine/protein secretion term was identified by more than 2 gene set enrichment methods. WWP2 is a member of ligase activity pathways and functions as a ligase for and mediates degradation of PTEN, whose gene is expressed in mouse oral development , . RNF217 is located at 6q22.31, a genomic region reported to be associated with oral cleft . ROBO2 is a receptor for SLIT2 and possibly SLIT1. SLIT1 and SLIT2 appear to work cooperatively to establish anatomical midlines during neuronal development and establishment of olfactory organization . Gene SLIT1 is also expressed in the primary and secondary enamel knots during molar tooth cusp formation . INS may impact caries through insulin sensitivity  or more controversially through the activation of dentin-related genes , . Insulin receptor binding sites are present on rat incisors . None of these relationships are “smoking guns” for caries development, but the gene sets and the subset of tooth-related genes raise interesting possible mechanisms for caries. These contributing genes encompass multiple functions or biological processes related to tooth development or dental caries, suggesting that our gene set enrichment analysis was effective and the findings were insightful to the understanding of molecular mechanisms of disease at the system level.
Although the genetic research has been applied to dental caries for a long time (see a recent review ), interpretation of the results remains challenging. In our gene set enrichment analysis, few GO terms or genes we identified exhibit explicit roles for caries development. One possible reason is the complex characteristics of dental caries. While many caries risk factors have been reported, few of them have been rigorously replicated or confirmed . Thus, the predefined gene sets may be too general to play definitive functions in cariogenesis.
In this study, four popular gene set analysis methods, i.e., GenGen, ALIGATOR, SRT, and the mixed model, were applied to a real GWAS dataset. Although our primary interest is to unveil the genetic components of dental caries, these results also provided a comprehensive benchmark resource to compare these methods. We only observed limited consistency among the outputs of different algorithms. The inconsistency is not unexpected, mainly because different methods employ different intrinsic strategies and may test different null hypotheses (i.e., competitive vs. self-contained null hypothesis). In addition, different ways to preprocess GWAS data might influence the enrichment results. For example, one important step in performing gene set analysis of GWAS data is to map SNPs to genes and compute a gene-based statistical value. Typically, only a subset of SNPs within a gene plays roles in the disease, yet taking all the SNPs into account will likely reduce the test power. However, in practical applications, it is difficult to find the most relevant SNPs for gene set analysis. Many approaches, like GenGen, denote the most significant SNP as gene’s representative, which may exclude important additional SNPs if a gene has more than one association signal. Using ALIGATOR, all SNPs mapped to a gene are consulted, and a gene is defined as significant if it harbors at least one nominally significant SNP, requiring a predefined threshold that may be chosen arbitrarily. Therefore, the analysis results from ALIGATOR could be sensitive to the choice of threshold in different data sets . Similarly, in the SRT method, all SNPs mapped to a gene are considered, and this approach also requires a preselected threshold to define the associated SNPs. One advantage in SRT is its incorporated permutation test by randomly swapping case/control labels among samples to reduce the sensitivity driven by the choice of threshold. In contrast, the mixed model approach accounts for the p-values of all the SNPs mapped to a gene without requiring predefined thresholds. Thus, this method avoids potential arbitrary definitions and quantitatively leverages the information of all SNPs.
One limitation in this study is the FDR values attained using the four methods are quite different from each other. The top 9 GO terms identified by ALIGATOR had an FDR value of 1. The situation is better in SRT, but the top GO terms also hardly reach a noteworthy FDR significance level. The high FDR could be the result of several factors. One is the inherent drawbacks of the tools used. For example, two GO terms, ‘Ligase activity forming carbon nitrogen bonds’ and ‘Ubiquitin protein ligase activity,’ were ranked as the fourth and the fifth most significant gene sets, respectively, in both results by ALIGATOR and the mixed model. However, their FDR values differed substantially in the two results. Another possible reason for this high FDR might be attributed to the incomplete information in the current annotation databases, especially for some phenotypes without much molecular biology knowledge. In contrast to most common diseases such as cancer, the functional annotation for dental caries has been very limited so far. In fact, we also performed a gene set enrichment analysis using the canonical pathways from KEGG , a widely used pathway database. There were only a small number of KEGG pathways eligible for our analysis (181 pathways with ≥ 5 and ≤ 250 genes), and none were significant KEGG pathways at FDR < 0.05. The failed detection of promising pathways for dental caries reflected that most, if not all, genes in the current version of the KEGG database are not thoroughly annotated. Another limitation in pathway annotations is that we used an old version of the GO term set (MSigDB, version 3.0, C5). New versions of GO data were released during our data analysis, which now included more than 10,000 GO terms (06/26/2013 release). However, major efforts are needed to process the redundancy of genes in GO terms, as processed in version 3.0, C5, to avoid an over-adjustment through multiple testing correction. This work, as well as more robust pathway enrichment analysis in future, may help better define dental caries pathways.
In summary, we applied four representative gene set enrichment analysis methods to currently available dental caries GWAS data. Our work, to date, is the first gene set enrichment study for this worldwide disease. We reported 13 significantly associated and 17 marginally associated GO terms as likely involved in dental caries via their gene functions. The findings provided insights and interpretations into the underlying biological process for dental caries. Our study mainly focused on genetic signals in GWAS data. In future work, an integration of other genetic and genomic information (such as gene expression, linkage scan and protein-protein interaction network , , evidence from multiple species , and multi-dimensional functional module analysis ) may open new avenues to understand the etiology of dental caries.
List of genes that are nominally significant with dental caries in the 30 top GO terms. This table includes 383 nominally significant genes that appeared in the top 30 GO terms shown in Table 2. Genes that have at least one SNP with a p-value < 0.05 are regarded as nominally significant genes (without multiple testing correction).
We thank all the participants and investigators in the original GWAS that made this work possible.
Conceived and designed the experiments: Zhongming Zhao LW PJ MLM EF KTC. Performed the experiments: QW PJ LW. Analyzed the data: QW PJ LW KTC. Wrote the paper: QW PJ Zhongming Zhao LW EF KTC Zhen Zeng MLM.
- 1. Azevedo LF, Pecharki GD, Brancher JA, Cordeiro CA Jr, Medeiros KG, et al. (2010) Analysis of the association between lactotransferrin (LTF) gene polymorphism and dental caries. J Appl Oral Sci 18: 166–170. doi: 10.1590/s1678-77572010000200011
- 2. Wendell S, Wang X, Brown M, Cooper ME, DeSensi RS, et al. (2010) Taste genes associated with dental caries. J Dent Res 89: 1198–1202. doi: 10.1177/0022034510381502
- 3. Kang SW, Yoon I, Lee HW, Cho J (2011) Association between AMELX polymorphisms and dental caries in Koreans. Oral Dis 17: 399–406. doi: 10.1111/j.1601-0825.2010.01766.x
- 4. Valarini N, Maciel SM, Moura SK, Poli-Frederico RC (2012) Association of dental caries with HLA class II allele in brazilian adolescents. Caries Res 46: 530–535. doi: 10.1159/000341188
- 5. Tannure PN, Kuchler EC, Falagan-Lotsch P, Amorim LM, Raggio Luiz R, et al. (2012) MMP13 polymorphism decreases risk for dental caries. Caries Res 46: 401–407. doi: 10.1159/000339379
- 6. Zeng Z, Shaffer JR, Wang X, Feingold E, Weeks DE, et al. (2013) Genome-wide association studies of Pit-and-Fissure- and Smooth-surface Caries in permanent dentition. J Dent Res 92: 432–437. doi: 10.1177/0022034513481976
- 7. Shaffer JR, Wang X, Feingold E, Lee M, Begum F, et al. (2011) Genome-wide association scan for childhood caries implicates novel genes. J Dent Res 90: 1457–1462. doi: 10.1177/0022034511422910
- 8. Shaffer JR, Feingold E, Wang X, Lee M, Tcuenco K, et al. (2013) GWAS of dental caries patterns in the permanent dentition. J Dent Res 92: 38–44. doi: 10.1177/0022034512463579
- 9. Hindorff LA, Sethupathy P, Junkins HA, Ramos EM, Mehta JP, et al. (2009) Potential etiologic and functional implications of genome-wide association loci for human diseases and traits. Proc Natl Acad Sci USA 106: 9362–9367. doi: 10.1073/pnas.0903103106
- 10. Wang L, Jia P, Wolfinger RD, Chen X, Grayson BL, et al. (2011) An efficient hierarchical generalized linear mixed model for pathway analysis of genome-wide association studies. Bioinformatics 27: 686–692. doi: 10.1093/bioinformatics/btq728
- 11. Wang K, Li M, Hakonarson H (2010) Analysing biological pathways in genome-wide association studies. Nat Rev Genet 11: 843–854. doi: 10.1038/nrg2884
- 12. Wang L, Jia P, Wolfinger RD, Chen X, Zhao Z (2011) Gene set analysis of genome-wide association studies: methodological issues and perspectives. Genomics 98: 1–8. doi: 10.1016/j.ygeno.2011.04.006
- 13. Jia P, Wang L, Meltzer HY, Zhao Z (2010) Common variants conferring risk of schizophrenia: a pathway analysis of GWAS data. Schizophr Res 122: 38–42. doi: 10.1016/j.schres.2010.07.001
- 14. Kao CF, Jia P, Zhao Z, Kuo PH (2012) Enriched pathways for major depressive disorder identified from a genome-wide association study. Int J Neuropsychopharmacol 15: 1401–1411. doi: 10.1017/s1461145711001891
- 15. Elbers CC, van Eijk KR, Franke L, Mulder F, van der Schouw YT, et al. (2009) Using genome-wide pathway analysis to unravel the etiology of complex diseases. Genet Epidemiol 33: 419–431. doi: 10.1002/gepi.20395
- 16. Perry JRB, McCarthy MI, Hattersley AT, Zeggini E, Weedon MN, et al. (2009) Interrogating type 2 diabetes genome-wide association data using a biological pathway-based approach. Diabetes 58: 1463–1467. doi: 10.2337/db08-1378
- 17. Zhong H, Yang X, Kaplan LM, Molony C, Schadt EE (2010) Integrating pathway analysis and genetics of gene expression for genome-wide association studies. Am J Hum Genet 86: 581–591. doi: 10.1016/j.ajhg.2010.02.020
- 18. Chen X, Wang L, Hu B, Guo M, Barnard J, et al. (2010) Pathway-based analysis for genome-wide association studies using supervised principal components. Genet Epidemiol 34: 716–724. doi: 10.1002/gepi.20532
- 19. Li D, Duell EJ, Yu K, Risch HA, Olson SH, et al. (2012) Pathway analysis of genome-wide association study data highlights pancreatic development genes as susceptibility factors for pancreatic cancer. Carcinogenesis 33: 1384–1390. doi: 10.1093/carcin/bgs151
- 20. Menashe I, Figueroa JD, Garcia-Closas M, Chatterjee N, Malats N, et al. (2012) Large-scale pathway-based analysis of bladder cancer genome-wide association data from five studies of European background. PLoS One 7: e29396. doi: 10.1371/journal.pone.0029396
- 21. Chen LS, Hutter CM, Potter JD, Liu Y, Prentice RL, et al. (2010) Insights into colon cancer etiology via a regularized approach to gene set analysis of GWAS data. Am J Hum Genet 86: 860–871. doi: 10.1016/j.ajhg.2010.04.014
- 22. Menashe I, Maeder D, Garcia-Closas M, Figueroa JD, Bhattacharjee S, et al. (2010) Pathway analysis of breast cancer genome-wide association study highlights three pathways and one canonical signaling cascade. Cancer Res 70: 4453–4459. doi: 10.1158/0008-5472.can-09-4502
- 23. Holmans P, Green EK, Pahwa JS, Ferreira MA, Purcell SM, et al. (2009) Gene ontology analysis of GWA study data sets provides insights into the biology of bipolar disorder. Am J Hum Genet 85: 13–24. doi: 10.1016/j.ajhg.2009.05.011
- 24. Wang K, Li M, Bucan M (2007) Pathway-based approaches for analysis of genomewide association studies. Am J Hum Genet 81: 1278–1283. doi: 10.1086/522374
- 25. O'Dushlaine C, Kenny E, Heron EA, Segurado R, Gill M, et al. (2009) The SNP ratio test: pathway analysis of genome-wide association datasets. Bioinformatics 25: 2762–2763. doi: 10.1093/bioinformatics/btp448
- 26. Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, et al. (2000) Gene ontology: tool for the unification of biology. Nat Genet 25: 25–29. doi: 10.1038/75556
- 27. Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, et al. (2005) Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci USA 102: 15545–15550. doi: 10.1073/pnas.0506580102
- 28. Purcell S, Neale B, Todd-Brown K, Thomas L, Ferreira MA, et al. (2007) PLINK: a tool set for whole-genome association and population-based linkage analyses. Am J Hum Genet 81: 559–575. doi: 10.1086/519795
- 29. Fehringer G, Liu G, Briollais L, Brennan P, Amos CI, et al. (2012) Comparison of pathway analysis approaches using lung cancer GWAS data sets. PLoS One 7: e31816. doi: 10.1371/journal.pone.0031816
- 30. Jia P, Wang L, Meltzer HY, Zhao Z (2011) Pathway-based analysis of GWAS datasets: effective but caution required. Int J Neuropsychopharmacol 14: 567–572. doi: 10.1017/s1461145710001446
- 31. Benjamini Y, Hochberg Y (1995) Controlling the false discovery rate - a practical and powerful approach to multiple testing. J R Stat Soc Series B Stat Methodol 57: 289–300.
- 32. Chino T, Santer DM, Giordano D, Chen C, Li C, et al. (2009) Effects of oral commensal and pathogenic bacteria on human dendritic cells. Oral Microbiol Immunol 24: 96–103. doi: 10.1111/j.1399-302x.2008.00478.x
- 33. Peyyala R, Kirakodu SS, Novak KF, Ebersole JL (2012) Oral microbial biofilm stimulation of epithelial cell responses. Cytokine 58: 65–72. doi: 10.1016/j.cyto.2011.12.016
- 34. Wilson M, Seymour R, Henderson B (1998) Bacterial perturbation of cytokine networks. Infect Immun 66: 2401–2409.
- 35. Fischer CL, Drake DR, Dawson DV, Blanchette DR, Brogden KA, et al. (2012) Antibacterial activity of sphingoid bases and fatty acids against Gram-positive and Gram-negative bacteria. Antimicrob Agents Chemother 56: 1157–1161. doi: 10.1128/aac.05151-11
- 36. Nixon GF (2009) Sphingolipids in inflammation: pathological implications and potential therapeutic targets. Br J Pharmacol 158: 982–993. doi: 10.1111/j.1476-5381.2009.00281.x
- 37. Fuchs SY (2012) Ubiquitination-mediated regulation of interferon responses. Growth Factors 30: 141–148. doi: 10.3109/08977194.2012.669382
- 38. Grant MM, Kolamunne RT, Lock FE, Matthews JB, Chapple IL, et al. (2010) Oxygen tension modulates the cytokine response of oral epithelium to periodontal bacteria. J Clin Periodontol 37: 1039–1048. doi: 10.1111/j.1600-051x.2010.01622.x
- 39. Lawrence DW, Kornbluth J (2012) E3 ubiquitin ligase NKLAM is a macrophage phagosome protein and plays a role in bacterial killing. Cell Immunol 279: 46–52. doi: 10.1016/j.cellimm.2012.09.004
- 40. Cho KW, Cho SW, Lee JM, Lee MJ, Gang HS, et al. (2008) Expression of phosphorylated forms of ERK, MEK, PTEN and PI3K in mouse oral development. Gene Expr Patterns 8: 284–290. doi: 10.1016/j.gep.2007.12.001
- 41. Maddika S, Kavela S, Rani N, Palicharla VR, Pokorny JL, et al. (2011) WWP2 is an E3 ubiquitin ligase for PTEN. Nat Cell Biol 13: 728–733. doi: 10.1038/ncb2240
- 42. Vieira AR, McHenry TG, Daack-Hirsch S, Murray JC, Marazita ML (2008) Candidate gene/loci studies in cleft lip/palate and dental anomalies finds novel susceptibility genes for clefts. Genet Med 10: 668–674. doi: 10.1097/gim.0b013e3181833793
- 43. Nguyen-Ba-Charvet KT, Plump AS, Tessier-Lavigne M, Chedotal A (2002) Slit1 and slit2 proteins control the development of the lateral olfactory tract. J Neurosci 22: 5473–5480.
- 44. Loes S, Luukko K, Kvinnsland IH, Kettunen P (2001) Slit1 is specifically expressed in the primary and secondary enamel knots during molar tooth cusp formation. Mech Dev 107: 155–157. doi: 10.1016/s0925-4773(01)00454-3
- 45. Mozaffari MS, Abdelsayed R, Zakhary I, El-Salanty M, Liu JY, et al. (2011) Submandibular gland and caries susceptibility in the obese Zucker rat. J Oral Pathol Med 40: 194–200. doi: 10.1111/j.1600-0714.2010.00965.x
- 46. Pourdeihimi B, Keene JJ Jr, Gargiulo AW, Toto PD (1980) Insulin I125 distribution within oral tissues. J Periodontol 51: 199–205. doi: 10.1902/jop.19188.8.131.52
- 47. Torres PU, Prie D, Beck L, De Brauwere D, Leroy C, et al. (2009) Klotho gene, phosphocalcic metabolism, and survival in dialysis. J Ren Nutr 19: 50–56. doi: 10.1053/j.jrn.2008.10.018
- 48. Martineau-Doize B, McKee MD, Warshawsky H, Bergeron JJ (1986) In vivo demonstration by radioautography of binding sites for insulin in liver, kidney, and calcified tissues of the rat. Anat Rec 214: 130–140. doi: 10.1002/ar.1092140205
- 49. Werneck RI, Mira MT, Trevilatto PC (2010) A critical review: an overview of genetic influence on dental caries. Oral Dis 16: 613–623. doi: 10.1111/j.1601-0825.2010.01675.x
- 50. Kanehisa M, Goto S (2000) KEGG: kyoto encyclopedia of genes and genomes. Nucleic Acids Res 28: 27–30. doi: 10.1093/nar/28.1.27
- 51. Sun J, Jia P, Fanous AH, Webb BT, van den Oord EJ, et al. (2009) A multi-dimensional evidence-based candidate gene prioritization approach for complex diseases-schizophrenia as a case. Bioinformatics 25: 2595–2602. doi: 10.1093/bioinformatics/btp428
- 52. Jia P, Wang L, Fanous AH, Pato CN, Edwards TL, et al. (2012) Network-assisted investigation of combined causal signals from genome-wide association studies in schizophrenia. PLoS Comput Biol 8: e1002587. doi: 10.1371/journal.pcbi.1002587
- 53. Zhao Z, Guo AY, van den Oord EJ, Aliev F, Jia P, et al. (2012) Multi-species data integration and gene ranking enrich significant results in an alcoholism genome-wide association study. BMC Genomics 13 Suppl 8S16. doi: 10.1186/1471-2164-13-s8-s16
- 54. Zhang S, Liu CC, Li W, Shen H, Laird PW, et al. (2012) Discovery of multi-dimensional modules by integrative analysis of cancer genomic data. Nucleic Acids Res 40: 9379–9391. doi: 10.1093/nar/gks725