Previously reported placebo-response-associated variants do not predict patient outcomes in inflammatory disease Phase III trial placebo arms.

In clinical trials, a placebo response refers to improvement in disease symptoms arising from the psychological effect of receiving a treatment rather than the actual treatment under investigation. Previous research has reported genomic variation associated with the likelihood of observing a placebo response, but these studies have been limited in scope and have not been validated. Here, we analyzed whole-genome sequencing data from 784 patients undergoing placebo treatment in Phase III Asthma or Rheumatoid Arthritis trials to assess the impact of previously reported variation on patient outcomes in the placebo arms and to identify novel variants associated with the placebo response. Contrary to expectations based on previous reports, we did not observe any statistically significant associations between genomic variants and placebo treatment outcome. Our findings suggest that the biological origin of the placebo response is complex and likely to be variable between disease areas.

The neuronal sortilin-related receptor SORL1 is genetically associated with Alzheimer disease.

The recycling of the amyloid precursor protein (APP) from the cell surface via the endocytic pathways plays a key role in the generation of amyloid beta peptide (Abeta) in Alzheimer disease. We report here that inherited variants in the SORL1 neuronal sorting receptor are associated with late-onset Alzheimer disease. These variants, which occur in at least two different clusters of intronic sequences within the SORL1 gene (also known as LR11 or SORLA) may regulate tissue-specific expression of SORL1. We also show that SORL1 directs trafficking of APP into recycling pathways and that when SORL1 is underexpressed, APP is sorted into Abeta-generating compartments. These data suggest that inherited or acquired changes in SORL1 expression or function are mechanistically involved in causing Alzheimer disease.

Homologous chromosomes make contact at the sites of double-strand breaks in genes in somatic G0/G1-phase human cells.

Double-strand DNA breaks (DSBs) are continuously induced in cells by endogenously generated free radicals and exogenous genotoxic agents such as ionizing radiation. DSBs activate the kinase activity in sensor proteins such as ATM and DNA-PK, initiating a complex DNA damage response that coordinates various DNA repair pathways to restore genomic integrity. In this study, we report the unexpected finding that homologous chromosomes contact each other at the sites of DSBs induced by either radiation or the endonuclease I-PpoI in human somatic cells. Contact involves short segments of homologous chromosomes and is centered on a DSB in active genes but does not occur at I-PpoI sites in intergenic DNA. I-PpoI-induced contact between homologous genes is abrogated by the transcriptional inhibitors actinomycin D and α-amanitin and requires the kinase activity of ATM but not DNA-PK. Our findings provide documentation of a common transcription-related and ATM kinase-dependent mechanism that induces contact between allelic regions of homologous chromosomes at sites of DSBs in human somatic cells.

Replication of genome wide association identified candidate genes confirm the role of common and rare variants in PAX7 and VAX1 in the etiology of nonsyndromic CL(P).

Following recent genome wide association studies (GWAS), significant genetic associations have been identified for several genes with nonsyndromic cleft lip with or without cleft palate (CL(P)). To replicate two of these GWAS signals, we investigated the role of common and rare variants in the PAX7 and VAX1 genes. TaqMan genotyping was carried out for SNPs in VAX1 and PAX7 and transmission disequilibrium test (TDT) was performed to test for linkage and association in each population. Direct sequencing in and around the PAX7 and VAX1 genes in 1,326 individuals of European and Asian ancestry was done. The TDT analysis showed strong associations with markers in VAX1 (rs7078160, P = 2.7E-06 and rs475202, P = 0.0002) in a combined sample of Mongolian and Japanese CL(P) case-parent triads. Analyses using parent-of-origin effects showed significant excess transmission of the minor allele from both parents with the effect in the mothers (P = 6.5E-05, OR (transmission) = 1.91) more striking than in the fathers (P = 0.004, OR (transmission) = 1.67) for VAX1 marker rs7078160 in the combined Mongolian and Japanese samples when all cleft types were combined. The rs6659735 trinucleotide marker in PAX7 was significantly associated with all the US cleft groups combined (P = 0.007 in all clefts and P = 0.02 in CL(P)). Eight rare missense mutations found in PAX7 and two rare missense mutations in VAX1. Our study replicated previous GWAS findings for markers in VAX1 in the Asian population, and identified rare variants in PAX7 and VAX1 that may contribute to the etiology of CL(P). Determining the role of rare variants clearly warrants further investigation.

Epigenetic influence of KAT6B and HDAC4 in the development of skeletal malocclusion.

INTRODUCTION: Genetic influences on the development of malocclusion include heritable effects on both masticatory muscles and jaw skeletal morphology. Beyond genetic variations, however, the characteristics of muscle and bone are also influenced by epigenetic mechanisms that produce differences in gene expression. We studied 2 enzymes known to change gene expressions through histone modifications, chromatin-modifying histone acetyltransferase KAT6B and deacetylase HDAC4, to determine their associations with musculoskeletal variations in jaw deformation malocclusions. METHODS: Samples of masseter muscle were obtained from subjects undergoing orthognathic surgery from 6 malocclusion classes based on skeletal sagittal and vertical dysplasia. The muscles were characterized for fiber type properties by immunohistochemistry, and their total RNA was isolated for gene expression studies by microarray analysis and quantitative real-time polymerase chain reaction. RESULTS: Gene expressions for fast isoforms of myosins and contractile regulatory proteins and for KAT6B and HDAC4 were severalfold greater in masseter muscles from a patient with a deepbite compared with one with an open bite, and genes related to exercise and activity did not differ substantially. In the total population, expressions of HDAC4 (P = 0.03) and KAT6B (P = 0.004) were significantly greater in subjects with sagittal Class III than in Class II malocclusion, whereas HDAC4 tended to correlate negatively with slow myosin type I and positively with fast myosin gene, especially type IIX. CONCLUSIONS: These data support other published reports of epigenetic regulation in the determination of skeletal muscle fiber phenotypes and bone growth. Further investigations are needed to elucidate how this regulatory model might apply to musculoskeletal development and malocclusion.

Genome-wide association scan of dental caries in the permanent dentition.

BACKGROUND: Over 90% of adults aged 20 years or older with permanent teeth have suffered from dental caries leading to pain, infection, or even tooth loss. Although caries prevalence has decreased over the past decade, there are still about 23% of dentate adults who have untreated carious lesions in the US. Dental caries is a complex disorder affected by both individual susceptibility and environmental factors. Approximately 35-55% of caries phenotypic variation in the permanent dentition is attributable to genes, though few specific caries genes have been identified. Therefore, we conducted the first genome-wide association study (GWAS) to identify genes affecting susceptibility to caries in adults. METHODS: Five independent cohorts were included in this study, totaling more than 7000 participants. For each participant, dental caries was assessed and genetic markers (single nucleotide polymorphisms, SNPs) were genotyped or imputed across the entire genome. Due to the heterogeneity among the five cohorts regarding age, genotyping platform, quality of dental caries assessment, and study design, we first conducted genome-wide association (GWA) analyses on each of the five independent cohorts separately. We then performed three meta-analyses to combine results for: (i) the comparatively younger, Appalachian cohorts (N = 1483) with well-assessed caries phenotype, (ii) the comparatively older, non-Appalachian cohorts (N = 5960) with inferior caries phenotypes, and (iii) all five cohorts (N = 7443). Top ranking genetic loci within and across meta-analyses were scrutinized for biologically plausible roles on caries. RESULTS: Different sets of genes were nominated across the three meta-analyses, especially between the younger and older age cohorts. In general, we identified several suggestive loci (P-value ≤ 10E-05) within or near genes with plausible biological roles for dental caries, including RPS6KA2 and PTK2B, involved in p38-depenedent MAPK signaling, and RHOU and FZD1, involved in the Wnt signaling cascade. Both of these pathways have been implicated in dental caries. ADMTS3 and ISL1 are involved in tooth development, and TLR2 is involved in immune response to oral pathogens. CONCLUSIONS: As the first GWAS for dental caries in adults, this study nominated several novel caries genes for future study, which may lead to better understanding of cariogenesis, and ultimately, to improved disease predictions, prevention, and/or treatment.

Heritable factors play a major role in determining host responses to Wuchereria bancrofti infection in an isolated South Pacific island population.

BACKGROUND: It is increasingly recognized that host genetic factors may play an important role in determining the outcome of filarial infections. To test this hypothesis in bancroftian lymphatic filariasis, pedigree data were collected twice during an 18-year period from an isolated Polynesian population living on a Pacific island where lymphatic filariasis is endemic. METHODS: Using variance-component analysis, we examined the contribution of shared genetic and environmental effects on host clinical and immune responses to filarial infection, along with potential confounding determinants. RESULTS: Sex was found to have a negligible influence on heritability estimates, but shared-household effects accounted for up to 32% of host variability. After accounting for these shared-household effects, heritability estimates suggested that levels of microfilariae and circulating adult worm antigen, as well as host eosinophil and immunoglobulin G antibody responses to larval and adult worm antigens, were highly heritable (range of heritability estimates, 0.15-0.84). CONCLUSIONS: These data provide evidence of a key role for genetic factors in determining the host response to filarial infections in humans and emphasize the complexity of the relationships among the host, parasite, and environment.

Uterine leiomyomata and decreased height: a common HMGA2 predisposition allele.

Uterine leiomyomata (UL) are the most common female pelvic tumors and the primary indication for hysterectomy in the United States. We assessed genetic liability for UL by a known embryonic proliferation modulator, HMGA2, in 248 families ascertained through medical record-confirmed affected sister-pairs. Using a (TC)( n ) repeat in the 5' UTR and 17 SNPs spanning HMGA2, permutation-based association tests identified a significant increase in transmission of a single TC repeat allele (TC227) with UL (allele-specific P = 0.00005, multiple testing corrected min-P = 0.0049). The hypothesis that TC227 is a pathogenic variant is supported by a trend towards higher HMGA2 expression in TC227 allele-positive compared with non-TC227 UL tissue as well as by absence of culpable exonic sequence variants. HMGA2 has also been suggested recently by three genome-wide SNP studies to influence human height variation, and our examination of the affected sister-pair families revealed a significant association of TC227 with decreased height (allele-specific P = 0.00033, multiple testing corrected min-P = 0.016). Diminished stature and elevated risk of UL development have both been correlated with an earlier age of menarche, which may be the biological mechanism for TC227 effects as a tendency of women with TC227 to have an earlier onset of menarche was identified in our study population. These results indicate HMGA2 has a role in two growth-related phenotypes, UL predisposition and height, of which the former may affect future medical management decisions for many women.

Enabling genome-wide association testing with multiple diseases and no healthy controls.

MOTIVATION: While large-scale whole genome sequencing is feasible the high costs compel investigators to focus on disease subjects. As a result large sequencing datasets of samples with different diseases are often readily available, but not healthy controls to contrast them with. While it is possible to perform an association study using only diseases, the associations could be driven by a disease acting as a control and not the focal disease. METHODS: We developed a genotype-on-phenotype reverse regression with a Bayesian spike and slab prior to enable association testing in datasets with multiple diseases. This method, referred to as revreg, flagged associations (both common and rare) that were driven by diseases that were not of primary interest. RESULTS: Based on simulations, revreg had 80% power to detect an odds ratio of 1.74 for common variants (3500 samples total) and 3.73 for rare variants (14,000 samples total), with minimal type I error. For common variants, we tested this method on 3657 whole genome sequenced samples aimed at discovering variants associated with disease risk of Chronic Obstructive Pulmonary Disease using three other diseases as controls. We demonstrated detection of six highly significant associations likely due to Age-Related Macular Degeneration. In an exome dataset of 8836 samples aimed at characterizing rare variants associated with disease risk of Asthma, using five other diseases as controls, we detected and removed genic regions due to AMD (C3, CFH, CFHR5, CFI, and DNMT3A) and RA (KRTAP13-4).

The impact of race as a risk factor for symptom severity and age at diagnosis of uterine leiomyomata among affected sisters.

OBJECTIVE: The objective of the study was to identify risk factors for uterine leiomyomata (UL) in a racially diverse population of women with a family history of UL, and to evaluate their contribution to disease severity and age at diagnosis. STUDY DESIGN: We collected and analyzed epidemiologic data from 285 sister pairs diagnosed with UL. Risk factors for UL-related outcomes were compared among black (n = 73) and white (n = 212) sister pairs using univariate and multivariate regression models. RESULTS: Black women reported an average age at diagnosis of 5.3 years younger (SE, 1.1; P < .001) and were more likely to report severe disease (odds ratio, 5.22; 95% confidence interval, 1.99-13.7, P < .001) than white women of similar socioeconomic status. CONCLUSION: Self-reported race is a significant factor in the severity of UL among women with a family history of UL. Differences in disease presentation between races likely reflect underlying genetic heterogeneity. The affected sister-pair study design can address both epidemiological and genetic hypotheses about UL.

Measuring associations between the microbiota and repeated measures of continuous clinical variables using a lasso-penalized generalized linear mixed model.

BACKGROUND: Human microbiome studies in clinical settings generally focus on distinguishing the microbiota in health from that in disease at a specific point in time. However, microbiome samples may be associated with disease severity or continuous clinical health indicators that are often assessed at multiple time points. While the temporal data from clinical and microbiome samples may be informative, analysis of this type of data can be problematic for standard statistical methods. RESULTS: To identify associations between microbiota and continuous clinical variables measured repeatedly in two studies of the respiratory tract, we adapted a statistical method, the lasso-penalized generalized linear mixed model (LassoGLMM). LassoGLMM can screen for associated clinical variables, incorporate repeated measures of individuals, and address the large number of species found in the microbiome. As is common in microbiome studies, when the number of variables is an order of magnitude larger than the number of samples LassoGLMM can be imperfect in its variable selection. We overcome this limitation by adding a pre-screening step to reduce the number of variables evaluated in the model. We assessed the use of this adapted two-stage LassoGLMM for its ability to determine which microbes are associated with continuous repeated clinical measures.We found associations (retaining a non-zero coefficient in the LassoGLMM) between 10 laboratory measurements and 43 bacterial genera in the oral microbiota, and between 2 cytokines and 3 bacterial genera in the lung. We compared our associations with those identified by the Wilcoxon test after dichotomizing our outcomes and identified a non-significant trend towards differential abundance between high and low outcomes. Our two-step LassoGLMM explained more of the variance seen in the outcome of interest than other variants of the LassoGLMM method. CONCLUSIONS: We demonstrated a method that can account for the large number of genera detected in microbiome studies and repeated measures of clinical or longitudinal studies, allowing for the detection of strong associations between microbes and clinical measures. By incorporating the design strengths of repeated measurements and a prescreening step to aid variable selection, our two-step LassoGLMM will be a useful analytic method for investigating relationships between microbes and repeatedly measured continuous outcomes.

Analysis of protein-altering variants in telomerase genes and their association with MUC5B common variant status in patients with idiopathic pulmonary fibrosis: a candidate gene sequencing study.

BACKGROUND: Idiopathic pulmonary fibrosis (IPF) risk has a strong genetic component. Studies have implicated variations at several loci, including TERT, surfactant genes, and a single nucleotide polymorphism at chr11p15 (rs35705950) in the intergenic region between TOLLIP and MUC5B. Patients with IPF who have risk alleles at rs35705950 have longer survival from the time of IPF diagnosis than do patients homozygous for the non-risk allele, whereas patients with shorter telomeres have shorter survival times. We aimed to assess whether rare protein-altering variants in genes regulating telomere length are enriched in patients with IPF homozygous for the non-risk alleles at rs35705950. METHODS: Between Nov 1, 2014, and Nov 1, 2016, we assessed blood samples from patients aged 40 years or older and of European ancestry with sporadic IPF from three international phase 3 clinical trials (INSPIRE, CAPACITY, ASCEND), one phase 2 study (RIFF), and US-based observational studies (Vanderbilt Clinical Interstitial Lung Disease Registry and the UCSF Interstitial Lung Disease Clinic registry cohorts) at the Broad Institute (Cambridge, MA, USA) and Human Longevity (San Diego, CA, USA). We also assessed blood samples from non-IPF controls in several clinical trials. We did whole-genome sequencing to assess telomere length and identify rare protein-altering variants, stratified by rs35705950 genotype. We also assessed rare functional variation in TERT exons and compared telomere length and disease progression across genotypes. FINDINGS: We assessed samples from 1510 patients with IPF and 1874 non-IPF controls. 30 (3%) of 1046 patients with an rs35705950 risk allele had a rare protein-altering variant in TERT compared with 34 (7%) of 464 non-risk allele carriers (odds ratio 0·40 [95% CI 0·24-0·66], p=0·00039). Subsequent analyses identified enrichment of rare protein-altering variants in PARN and RTEL1, and rare variation in TERC in patients with IPF compared with controls. We expanded our study population to provide a more accurate estimation of rare variant frequency in these four loci, and to calculate telomere length. The proportion of patients with at least one rare variant in TERT, PARN, TERC, or RTEL1 was higher in patients with IPF than in controls (149 [9%] of 1739 patients vs 205 [2%] of 8645 controls, p=2·44 × 10-8). Patients with IPF who had a variant in any of the four identified telomerase component genes had telomeres that were 3·69-16·10% shorter than patients without a variant in any of the four genes and had an earlier mean age of disease onset than patients without one or more variants (65·1 years [SD 7·8] vs 67·1 years [7·9], p=0·004). In the placebo arms of clinical trials, shorter telomeres were significantly associated with faster disease progression (1·7% predicted forced vital capacity per kb per year, p=0·002). Pirfenidone had treatment benefit regardless of telomere length (p=4·24 × 10-8 for telomere length lower than the median, p=0·0044 for telomere length greater than the median). INTERPRETATION: Rare protein-altering variants in TERT, PARN, TERC, and RTEL1 are enriched in patients with IPF compared with controls, and, in the case of TERT, particularly in individuals without a risk allele at the rs35705950 locus. This suggests that multiple genetic factors contribute to sporadic IPF, which might implicate distinct mechanisms of pathogenesis and disease progression. FUNDING: Genentech, National Institutes of Health, Francis Family Foundation, Pulmonary Fibrosis Foundation, Nina Ireland Program for Lung Health, US Department of Veterans Affairs.

The pathogenesis of filarial lymphedema: is it the worm or is it the host?

Our understanding of the pathogenesis of filarial lymphedema, although evolving, is still limited. Recurrent bacterial infections play a major role in the progression of lymphedema to elephantiasis, but the host and parasite factors that trigger disease development are not known. Field studies in Haiti show that lymphedema and host responses to parasite antigens cluster in families, consistent with the hypothesis that host genes influence lymphedema susceptibility. The recent recognition that filarial parasites harbor the endosymbiotic bacteria, Wolbachia, also raises questions about the potential contribution of the inflammatory response to Wolbachia antigens to lymphedema development. In this review, we discuss potential risk factors for lymphedema and try to integrate these in a model of pathogenesis.

Assessment of families for excess risk of lymphedema of the leg in a lymphatic filariasis-endemic area.

The influence of host genes on the distribution of lymphedema due to lymphatic filariasis is unknown. To assess this, pedigree and disease information were collected from lymphedema patients in a lymphatic filariasis-endemic area. These patients were female, with an average age of approximately 40 years, who were enrolled between June 1995 and July 1999 in a lymphedema treatment clinic, and from the rural Haitian community served by the clinic. Interviews were conducted between September 1998 and December 1999. Families with multiple lymphedema cases were of similar size, with an average of 15 members, as those families with only a single lymphedema case. We determined whether families observed to have multiple lymphedema cases had a higher prevalence of lymphedema than expected when stratified population estimates and family size were considered. Lymphedema of the leg was excessive in 15 of 43 families with multiple lymphedema cases. The number of families demonstrating excess disease was significantly different than was expected based on population estimates of lymphedema prevalence (P = 0.026). Families with multiple cases of lymphedema were not significantly larger in family size than families with a single lymphedema cases. Twelve of the 15 families had a male with lymphedema, which influenced the interpretation of the results. The significance of these results is discussed.

A family study of lymphedema of the leg in a lymphatic filariasis-endemic area.

The risk of filarial lymphedema may not be equivalent for all members of filaria-exposed populations. While evidence for a genetic factor that influences acquisition of infection has been growing, very little work has addressed whether there is a genetic basis to the development of disease due to lymphatic filariasis. We designed a family study of lymphedema in a rural community in Haiti to assess disease aggregation. Three hundred sixty-eight female patients sixteen years of age or older were enrolled at a lymphedema treatment clinic between June 1995 and December 1999. After applying additional eligibility criteria, 172 probands were enrolled into the family study for detailed pedigree collection between September 1998 and December 1999. Fifty-three lymphedema cases were identified among the probands' parents, full-siblings, children, half-siblings, and mating partners of the parents. Twelve of the 53 cases were among males. The proportion of cases occurring in a biologic parent of the proband was higher than in unrelated individuals married into the proband's family (P = 0.0010). This is the first large family study based on pedigrees to assess the familial aggregation of lymphedema due to filariasis. This family study will be useful to investigate the role of genes and environment in the development of filarial-related lymphedema.

Genome-wide association study of periodontal health measured by probing depth in adults ages 18-49 years.

The etiology of chronic periodontitis clearly includes a heritable component. Our purpose was to perform a small exploratory genome-wide association study in adults ages 18-49 years to nominate genes associated with periodontal disease-related phenotypes for future consideration. Full-mouth periodontal pocket depth probing was performed on participants (N = 673), with affected status defined as two or more sextants with probing depths of 5.5 mm or greater. Two variations of this phenotype that differed in how missing teeth were treated were used in analysis. More than 1.2 million genetic markers across the genome were genotyped or imputed and tested for genetic association. We identified ten suggestive loci (p-value ≤ 1E-5), including genes/loci that have been previously implicated in chronic periodontitis: LAMA2, HAS2, CDH2, ESR1, and the genomic region on chromosome 14q21-22 between SOS2 and NIN. Moreover, we nominated novel loci not previously implicated in chronic periodontitis or related pathways, including the regions 3p22 near OSBPL10 (a lipid receptor implicated in hyperlipidemia), 4p15 near HSP90AB2P (a heat shock pseudogene), 11p15 near GVINP1 (a GTPase pseudogene), 14q31 near SEL1L (an intracellular transporter), and 18q12 in FHOD3 (an actin cytoskeleton regulator). Replication of these results in additional samples is needed. This is one of the first research efforts to identify genetic polymorphisms associated with chronic periodontitis-related phenotypes by the genome-wide association study approach. Though small, efforts such this are needed in order to nominate novel genes and generate new hypotheses for exploration and testing in future studies.

Homologous chromosomes move and rapidly initiate contact at the sites of double-strand breaks in genes in G0-phase human cells.

We recently reported that homologous chromosomes make contact at the sites of double-strand breaks (DSBs) induced by ionizing radiation (IR) and the restriction endonuclease I-PpoI in G0/G1-phase somatic human cells. The contact involves short segments of homologous chromosomes and is centered on a DSB that occurs in a gene; contact does not occur at a DSB in intergenic DNA. Contact between homologous chromosomes is abrogated by inhibition of transcription and requires the kinase activity of ATM, but not DNA-PK. Here, we report additional insights into the mechanism underlying this novel phenomenon. We identify four patterns of homologous chromosome contact, and show that contact between homologous arms, but not centrosomes, is induced by IR. Significantly, we demonstrate that contact is induced by IR in non-proliferating, G0-phase human cells derived from tissue explants. Finally, we show that contact between homologous chromosomes is detectable as early as 5 min after IR. These results point to the existence of a mechanism that rapidly localizes homologous chromosome arms at sites of DSBs in genes in G0-phase human cells.

Multi-dimensional prioritization of dental caries candidate genes and its enriched dense network modules.

A number of genetic studies have suggested numerous susceptibility genes for dental caries over the past decade with few definite conclusions. The rapid accumulation of relevant information, along with the complex architecture of the disease, provides a challenging but also unique opportunity to review and integrate the heterogeneous data for follow-up validation and exploration. In this study, we collected and curated candidate genes from four major categories: association studies, linkage scans, gene expression analyses, and literature mining. Candidate genes were prioritized according to the magnitude of evidence related to dental caries. We then searched for dense modules enriched with the prioritized candidate genes through their protein-protein interactions (PPIs). We identified 23 modules comprising of 53 genes. Functional analyses of these 53 genes revealed three major clusters: cytokine network relevant genes, matrix metalloproteinases (MMPs) family, and transforming growth factor-beta (TGF-ß) family, all of which have been previously implicated to play important roles in tooth development and carious lesions. Through our extensive data collection and an integrative application of gene prioritization and PPI network analyses, we built a dental caries-specific sub-network for the first time. Our study provided insights into the molecular mechanisms underlying dental caries. The framework we proposed in this work can be applied to other complex diseases.

Association signals unveiled by a comprehensive gene set enrichment analysis of dental caries genome-wide association studies.

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.