Genome-Wide Association Study Reveals Multiple Loci Associated with Primary Tooth Development during Infancy
Tooth development is a highly heritable process which relates to other growth and developmental processes, and which interacts with the development of the entire craniofacial complex. Abnormalities of tooth development are common, with tooth agenesis being the most common developmental anomaly in humans. We performed a genome-wide association study of time to first tooth eruption and number of teeth at one year in 4,564 individuals from the 1966 Northern Finland Birth Cohort (NFBC1966) and 1,518 individuals from the Avon Longitudinal Study of Parents and Children (ALSPAC). We identified 5 loci at P<5×10−8, and 5 with suggestive association (P<5×10−6). The loci included several genes with links to tooth and other organ development (KCNJ2, EDA, HOXB2, RAD51L1, IGF2BP1, HMGA2, MSRB3). Genes at four of the identified loci are implicated in the development of cancer. A variant within the HOXB gene cluster associated with occlusion defects requiring orthodontic treatment by age 31 years.
Published in the journal:
. PLoS Genet 6(2): e32767. doi:10.1371/journal.pgen.1000856
Category:
Research Article
doi:
https://doi.org/10.1371/journal.pgen.1000856
Summary
Tooth development is a highly heritable process which relates to other growth and developmental processes, and which interacts with the development of the entire craniofacial complex. Abnormalities of tooth development are common, with tooth agenesis being the most common developmental anomaly in humans. We performed a genome-wide association study of time to first tooth eruption and number of teeth at one year in 4,564 individuals from the 1966 Northern Finland Birth Cohort (NFBC1966) and 1,518 individuals from the Avon Longitudinal Study of Parents and Children (ALSPAC). We identified 5 loci at P<5×10−8, and 5 with suggestive association (P<5×10−6). The loci included several genes with links to tooth and other organ development (KCNJ2, EDA, HOXB2, RAD51L1, IGF2BP1, HMGA2, MSRB3). Genes at four of the identified loci are implicated in the development of cancer. A variant within the HOXB gene cluster associated with occlusion defects requiring orthodontic treatment by age 31 years.
Introduction
Heritability of primary tooth emergence is estimated to be over 70% [1]. Abnormalities in tooth development are common with tooth agenesis alone affecting up to 10% of the population, ranking it as the most common developmental anomaly in humans [2]. Such abnormalities contribute to a variety of challenging and expensive orthodontic, prosthetic and surgical treatments and account for approximately 6% of all dental health care attendances [3]. Many genes implicated in primary dentition have regulatory functions important to several developmental processes in the embryo [4], and the developing tooth is a useful model for the study of organogenesis [5]. However, despite substantial research into tooth development in mice and human malformation syndromes [5], the genetic determinants of the normal variation in human tooth development have not been established.
To identify genetic loci regulating primary dentition we performed a general population based genome-wide association (GWA) study of tooth development in infancy among individuals from the 1966 Northern Finland Birth Cohort (NFBC1966) and the Avon Longitudinal Study of Parents and Children (ALSPAC). Specifically, we tested for associations with time to first tooth eruption and number of teeth by one year of age. These phenotypes are relevant to later tooth development because teeth largely acquire their final form at a very early age [6]. The availability of longitudinal birth cohort data allowed us to investigate life-course associations with dental occlusion defects.
Results
We tested 300,766 SNPs common to both studies (each used the Illumina platform). The analyses were adjusted for sex, gestational age and population structure (Materials and Methods). Results for the two cohorts were combined using fixed effects inverse variance meta-analysis. Five genetic loci were identified at genome-wide significance (P<5×10−8). Table 1 shows the top-ranking SNPs at each locus (see also Figure 1 and Figure 2, Figures S1, S2, S3). For all SNPs the allele associated with a delay in tooth eruption was associated with fewer teeth at the end of infancy. Table S1 shows details of the functions of genes linked to the identified loci.
The strongest association with both phenotypes was for SNP rs8079702, located 15 kb downstream of KCNJ2 (inward rectifier potassium channel 2) (P = 3.77×10−22 for time of first tooth, P = 1.24×10−14 for number of teeth; Table 1). There are no SNPs in KCNJ2 in our data, but rs8079702 had highest correlation with SNP rs4328485 which was the closest available SNP to KCNJ2 (r2 = 0.17; 1 kb away). KCNJ2 has been implicated in Pierre Robin sequence [7] and Andersen-Tawil syndrome [8], which show abnormalities in tooth development (missing teeth, delays in eruption) and are characterized by craniofacial anomalies such as narrowing of the jaw and cleft palate [8]. The second strongest association was for SNP rs5936487, located within the EDA (ectodermal dysplasia protein) gene (P = 6.18×10−11 for time of first tooth, P = 3.36×10−10 for number of teeth). EDA was fundamental in forming the first teeth in organisms [9], and mutations cause hypohidrotic ectodermal dysplasia (HED) and non-syndromic disorders of tooth agenesis [10].
The three remaining loci at genome-wide significance (P<5×10−8) have SNPs located within the genes RAD51L1 (RAD51-like1), IGF2BP1 (insulin-like growth factor 2 mRNA binding protein 1) and MSRB3 (methionine sulfoxide reductase B3). RAD51L1 is involved in DNA repair and a variant in the gene has been found to confer susceptibility to breast cancer [11]. It is responsible for protein kinase activity, and the injection of activators of protein kinase C (PKC) in rats causes delays in tooth eruption [12]. IGF2BP1 regulates the growth factor IGF2, and knockouts of the gene in mice suggest a role in organ development [13], while its expression is associated with ovarian cancer [14]. A microarray study in the developing mouse molar tooth found MSRB3 to be in the top 100 most expressed genes of 34,000 examined [15].
Each of the associated SNPs explain a small fraction of the residual phenotypic variation in time to first tooth (0.2%–1.6%, NFBC1966; 0.4%–1.5%, ALSPAC) and number of teeth by one year (0.2%–1.2%, NFBC1966; 0.5%–1.6%, ALSPAC), after controlling for sex and gestational age. Selecting the SNP with the most extreme signal for either phenotype to represent each locus (“top SNPs”), and analysing them together, the additive effects of these five top SNPs explain 2.9% of the variance of both tooth eruption time and number of teeth in the NFBC1966, and 4.2% and 4.0% of the variance in tooth eruption and number of teeth in ALSPAC. Without a suitable external replication cohort these estimates were derived in the two discovery cohorts and therefore may overestimate the true values due to the “winner's curse”. GWA studies have thus far explained only a small proportion of heritability [16], and our estimates are comparable with the variance explained in human height by a GWA study [17]. In order to identify variants with lower effect sizes or rarer variants larger sample sizes would be required.
We also summarized the predictive power of the five top SNPs by defining a ‘delayed tooth eruption’ measure as the number of alleles across the SNPs that delay tooth eruption. Figure 3 shows the number of delayed tooth eruption alleles against the mean of both time to first tooth eruption and number of teeth by one year in NFBC1966. Individuals with 8 or more delayed eruption alleles (10% of NFBC1966) have an average of 1.5 fewer teeth at 12 months, and later tooth eruption by 1.1 months, compared to individuals with 3 or fewer such alleles (11% of NFBC1966). Figure S4 shows the same plot for time to first tooth in ALSPAC.
In addition to the five loci attaining genome wide significance, there were 5 loci with SNPs that had P-values between 5×10−6 and 5×10−8 (Table 1). We investigated the biological functions of nearby genes to see if any of these loci were related to tooth development. These signals included SNP rs6504340, which is located between the developmental regulatory genes HOXB1 (homeobox B1) and HOXB2 (homeobox B2). Although previous studies have indicated that tooth development is independent of a Hox patterning program [18], Homeobox genes have recently been shown to be expressed in the dental mesenchyme in the pharyngeal teeth of bony fishes [9]. SNP rs6504340 lies 500 kb upstream of rs9674544 in IGF2BP1, but the two SNPs show almost no linkage disequilibrium with each other (r2 = 0.002 in NFBC1966 and r2 = 0.006 in ALSPAC). Furthermore, a test for association of rs6504340 conditional on rs9674544 was significant (P = 6.3×10−5 in NFBC1966 and P = 0.01 in ALSPAC; Materials and Methods), indicating that these represent two independent signals (Figure 1). We also identified three SNPs at 2q35, the most significant of which had r2 = 0.48 with a variant associated with breast cancer [11],[19], and SNP rs12424086 located close to the HMGA2 gene and 6 kb away from rs1042725, the SNP identified by a GWA study for adult and childhood height [20].
Given the influence of tooth development on dental occlusion, we hypothesized that genetic determinants of early tooth eruption may associate with dental occlusion later in life. We tested for associations between the SNP with the most extreme signal for either phenotype at each of the 10 identified loci and defects in occlusion requiring orthodontic treatment by the age of 31 years in the NFBC1966 (data not available in ALSPAC). A total of 611 individuals (13.5%) reported a defect in occlusion that had required orthodontic treatment. Of the 10 SNPs tested, SNP rs6504340 (HOXB gene cluster) gave a significant association, where each G allele (associated with delayed tooth eruption and lower number of teeth in infancy, Table 1), increased the odds of having an occlusal defect requiring orthodontic treatment by 35%, after adjusting for sex (odds ratio (OR) = 1.35, 95% CI = 1.16–1.57; P = 1.13×10−4; further adjustment for gestational age did not change the result). A smaller number of teeth at 1 year also predicted higher risk of orthodontic treatment (OR = 1.05, 95% CI = 1.01–1.09; P = 0.009). However, when number of teeth or time to first tooth were included in the model with dental occlusion as outcome, the associations with the G allele remained (P = 0.001, P = 1.71×10−4), suggesting an independent association between rs6504340 and dental occlusion.
Discussion
Teeth and several other organs have common growth and developmental pathways during early life [21]. The genes at the loci identified in our study have roles in organogenesis, growth and developmental processes, and cancer. Mutations in three of the genes lead to altered organogenesis and development; KCNJ2 (teeth, jaws, palates, ears, fingers, toes), EDA (teeth, hair, sweat glands, salivary glands) and IGF2BP1 (intestines) [8],[13],[22]. Of the loci at suggestive levels of significance, the HOXB gene cluster is an established regulator of development, and the HMGA2 gene has previously been associated with adult height [20]. Normal development and cancer both involve shifts between cell proliferation and differentiation [23] and genes regulating organ-specific growth are known to be involved in oncogenesis [24]. A previous study identified a common genetic link between an abnormal tooth development and cancer [25]. From our identified loci, IGF2BP1 and RAD51L1 have been implicated in cancer [11],[14] as have HOXB2, 2q35, and HMGA2 [19],[26],[27].
We provide the first detailed insight into the genetic architecture of primary dentition and our findings could have implications for the study of other developmental and organogenic processes. Exploiting the availability of longitudinal cohort data [28] we found an association between a variant within the HOXB gene cluster and the requirement for orthodontic treatment due to defective occlusion by the age of 31 years. Further GWA studies of developmental processes during infancy may establish whether the genetic determinants of infant development can contribute to the study of chronic diseases, such as cancer, that occur later in life.
Materials and Methods
Study population and phenotype description
The data was derived from two genome-wide scans of the geographically defined prospective birth cohorts; the NFBC1966 and ALSPAC. The NFBC1966 followed pregnancies in the two northernmost provinces of Finland with expected delivery dates in 1966. ALSPAC recruited mothers during pregnancy with expected dates of delivery between April 1991 and December 1992 from Bristol and the surrounding area in the South West of England. A total of 4,564 samples were available from the NFBC1966 and 1,518 from ALSPAC. In both cohorts, two separate measures of primary tooth development were collected: i) date of first tooth eruption (in months), and ii) number of teeth (measured at 12 months in NFBC1966 and 15 months in ALSPAC). In the NFBC1966 date of first tooth eruption and number of teeth was gathered by public health professionals during children's monthly visits to child welfare centers (parents carried a booklet where they had recorded the developmental milestones reached). In ALSPAC, parents reported the date of first tooth eruption and number of teeth at 15 months on a questionnaire. In order to ascertain the accuracy of the parental responses, a subsample were examined and validated by a dentist. Information on date of first tooth eruption was available for 4,523 individuals in the NFBC1966 (99% of available GWA samples) and 1396 (92%) in ALSPAC and for number of teeth, 4,326 (95%) in the NFBC1966 and 1,426 (94%) in ALSPAC. All aspects of the study were reviewed and approved by the Ethics Committee of the University of Oulu and the ALSPAC Law and Ethics Committee and by the respective local research committees. Participants (in NFBC1966) and parents (in ALSPAC) gave written informed consent.
Genotyping
The Illumina HumanCNV370-Duo DNA Analysis BeadChip was used for genotyping the NFBC1966, and Illumina HumanHap317K BeadChip for ALSPAC. The genotyping and quality control procedures have been described elsewhere [29],[30]. SNPs were excluded from the analysis if the call rate in the final sample was <95%, if there was a lack of Hardy-Weinberg Equilibrium (HWE) (P<10−4 in NFBC1966, P<5×10−7 in ALSPAC), or if the MAF was <1%. After quality control, 329,091 SNPs in NFBC1966 and 310,611 in ALSPAC were available. We report here the results from the 300,766 genotyped SNPs common to both studies.
Statistical analyses
Age of first tooth eruption in the NFBC1966 was recorded in months, such that the first tooth could have erupted at any time between the end of previous month and the end of the recorded month. In ALSPAC it was recorded to the nearest month and 3 individuals were recorded as having no teeth after 15 months. To account for the censoring in the two cohorts the outcome was analyzed using parametric survival analysis in the R software package 2.7.1.The Gaussian distribution gave a good fit to the data in both cohorts and was used to model the underlying event time. Number of teeth in the NFBC1966 was recorded at 12 months. In ALSPAC, measurements were taken at around 15 months but there was variability in the exact time of measurement, therefore the ALSPAC analysis was adjusted for age of measurement. Teeth typically erupt in pairs from the upper and lower jaw (75% of children had an even number of teeth in the NFBC1966), making the Poisson distribution inappropriate for modeling the number of teeth. Therefore ordinal logistic regression was used as implemented by the polr function in the R package. Analyses of the X chromosome treated males as homozygous females. The allele frequencies of the identified SNPs on the X chromosome did not differ significantly between the sexes. GWA analyses were adjusted for sex, gestational age and population stratification using principal components (PC). Each analysis was corrected for population stratification separately by including those of the top 10 PCs that were associated with the phenotype at P<0.05 [31]. For number of teeth, PCs 3, 6 and 9 were included in ALSPAC and none in the NFBC1966. For time to first tooth eruption no PCs were included in ALSPAC and PC 2 was included in the NFBC1966. After correction by PCs, the estimated variance inflation factors [32] for date of first tooth eruption were 1.039 and 1.047 in ALSPAC and NFBC1966 respectively, and 1.011 and 1.039 for number of teeth. Genomic control [32] was then used to correct the residual population stratification. The variance inflation factors from the meta-analyses were 1.012 for number of teeth and 1.015 date of first tooth eruption.
Results from the two studies were combined using fixed effects inverse variance meta-analysis [33]. Analyses were performed using the statistical package R and metaMapper (a meta-analysis software developed in-house). Conditional analyses were calculated using the likelihood ratio test comparing ordinal regression models, one including rs9674544 and the other including rs9674544 and rs6504340. Variance explained by each SNP was computed as 1 minus the ratio of variance of residuals of the model with age, gestational age and SNP to variance of residuals of the model with just age and gestational age. To correct for overfitting, each individual's phenotype was estimated from a model that did not include that individual. The total variance explained by the five loci reaching genome-wide significance was calculated similarly using the most associated SNPs for each phenotype at each locus. Additional tests for association with orthodontic treatment used the SNPs most associated with number of teeth at the 10 loci at P<5×10−6. Table 1 reports the top GWA signals at each of the ten loci (i.e. the SNP with the strongest association with either time to first tooth eruption or number of teeth at age 1 year).
URLs
Jackson Laboratory website, http://www.jax.org; NCBI, http://www.ncbi.nlm.nih.gov; R project, www.r-project.org; Uniprot, http://www.uniprot.org.
Supporting Information
Zdroje
1. HughesTE
BockmannMR
SeowK
GotjamanosT
GullyN
2007 Strong genetic control of emergence of human primary incisors. J Dent Res 86 1160 1165
2. NanciA
2008 Ten Cate's oral histology. Development, structure, and function, 7th ed. St Louis Mosby Elsevier, pg 96
3. AndersonR
RichmondS
ThomasDW
1999 Patient presentation at medical practices with dental problems: an analysis of the 1996 General Practice Morbidity Database for Wales. Br Dent J 186 297 300
4. ThesleffI
2000 Genetic basis of tooth development and dental defects. Acta Odontol Scand 58 191 194
5. TuckerA
SharpeP
2004 The cutting-edge of mammalian development; how the embryo makes teeth. Nat Rev Genet 5 499 508
6. KoussoulakouDS
MargaritisLH
KoussoulakosSL
2009 A curriculum vitae of teeth: evolution, generation, regeneration. Int J Biol Sci 5 226 243
7. BenkoS
FantesJA
AmielJ
KleinjanDJ
ThomasS
2009 Highly conserved non-coding elements on either side of SOX9 associated with Pierre Robin sequence. Nat Genet 41 359 364
8. YoonG
OberoiS
Tristani-FirouziM
EtheridgeSP
QuitaniaL
2006 Andersen-Tawil syndrome: prospective cohort analysis and expansion of the phenotype. Am J Med Genet A 140 312 321
9. FraserGJ
HulseyCD
BloomquistRF
UyesugiK
ManleyNR
2009 An ancient gene network is co-opted for teeth on old and new jaws. PLoS Biol 7 e31 doi:10.1371/journal.pbio.1000031
10. MuesGI
GriggsR
HartungAJ
WhelanG
BestLG
2009 From ectodermal dysplasia to selective tooth agenesis. Am J Med Genet A 149A 2037 2041
11. ThomasG
JacobsKB
KraftP
YeagerM
WacholderS
2009 A multistage genome-wide association study in breast cancer identifies two new risk alleles at 1p11.2 and 14q24.1 (RAD51L1). Nat Genet 41 579 584
12. WiseGE
YaoS
LiuD
2006 Injections of osteoprotegerin and PMA delay tooth eruption. Clin Anat 19 19 24
13. HansenTV
HammerNA
NielsenJ
MadsenM
DalbaeckC
2004 Dwarfism and impaired gut development in insulin-like growth factor II mRNA-binding protein 1-deficient mice. Mol Cell Biol 24 4448 4464
14. GuL
ShigemasaK
OhamaK
2004 Increased expression of IGF II mRNA-binding protein 1 mRNA is associated with an advanced clinical stage and poor prognosis in patients with ovarian cancer. Int J Oncol 24 671 678
15. PembertonTJ
LiF-Y
OkaS
Mendoza-FandinoGA
HsuY-H
2007 Identification of novel genes expressed during mouse tooth development by microarray gene expression analysis. Dev Dyn 236 2245 2257
16. ManolioTA
CollinsFS
CoxNJ
GoldsteinDB
HindorffLA
2009 Finding the missing heritability of complex diseases. Nature 461 747 753
17. WeedonMN
LangoH
LindgrenCM
WallaceC
EvansDM
2008 Genome-wide association analysis identifies 20 loci that influence adult height. Nat Genet 40 575 583
18. JamesCT
OhazamaA
TuckerAS
SharpePT
2002 Tooth development is independent of a Hox patterning programme. Dev Dyn 225 332 335
19. StaceySN
ManolescuA
SulemP
RafnarT
GudmundssonJ
2007 Common variants on chromosomes 2q35 and 16q2 confer susceptibility to estrogen receptor-positive breast cancer. Nat Genet 39 865 869
20. WeedonMN
LettreG
FreathyRM
LindgrenCM
VoightBF
2007 A common variant of HMGA2 is associated with adult and childhood height in the general population. Nat Genet 39 1245 1250
21. ThesleffI
VaahtokariA
PartanenAM
1995 Common molecular mechanisms regulating the development of teeth and other organs. Int J Dev Biol 39 35 50
22. MikkolaML
TheleffI
2003 Ectodysplasin signalling in development. Cytokine Growth Factor Rev 14 211 224
23. NunesFD
de AlmeidaFC
TucciR
de SousaSC
2003 Homeobox genes: a molecular link between development and cancer. Pesqui Odontol Bras 17 94 98
24. HallikasO
PalinK
SinjashinaN
RautiainenR
PartanenJ
2006 Genome-wide prediction of mammalian enhancers based on analysis of transcription-factor binding affinity. Cell 124 47 59
25. LammiL
ArteS
SomerM
JarvinenH
LahermoP
2004 Mutations in AXIN2 cause familial tooth agenesis and predispose to colorectal cancer. Am J Hum Genet 74 1043 1050
26. CollinsY
TanDF
PejovicT
MorG
QianF
2005 Identification of differentially expressed genes in clinically distinct groups of serous ovarian carcinomas using cDNA microarray. Int J Mol Med 14 43 53
27. GrierDG
ThompsonA
KwasniewskaA
McGonigleGJ
HallidayHL
2005 The pathophysiology of HOX genes and their role in cancer. J Pathol 205 154 171
28. ManolioTA
2009 Cohort studies and the genetics of complex disease. Nat Genet 41 5 6
29. SabattiC
ServiceSK
HartikainenAL
PoutaA
RipattiS
2009 Genome-wide association analysis of metabolic traits in a birth cohort from a founder population. Nat Genet 41 35 46
30. TimpsonNJ
TobiasJH
RichardsJB
SoranzoN
DuncanEL
2009 Common variants in the region around Osterix are associated with bone mineral density and growth in childhood. Hum Mol Genet 18 1510 1517
31. NovembreJ
StephensM
2008 Interpreting principal component analyses of spatial population genetic variation. Nat Genet 40 646 649
32. DevlinB
RoederK
1999 Genomic control for association studies. Biometrics 55 997 1004
33. SuttonAJ
AbramsKR
JonesDR
Sheldon TA SongF
2000 Methods for meta-analysis in medical research. Chichester, UK John Wiley & Sons, pg 58
Štítky
Genetika Reprodukční medicínaČlánek vyšel v časopise
PLOS Genetics
2010 Číslo 2
- Primární hyperoxalurie – aktuální možnosti diagnostiky a léčby
- Srdeční frekvence embrya může být faktorem užitečným v předpovídání výsledku IVF
- Akutní intermitentní porfyrie
- Vztah užívání alkoholu a mužské fertility
- Šanci na úspěšný průběh těhotenství snižují nevhodné hladiny progesteronu vznikající při umělém oplodnění
Nejčtenější v tomto čísle
- Genome-Wide Association Study in Asian Populations Identifies Variants in and Associated with Systemic Lupus Erythematosus
- Nuclear Pore Proteins Nup153 and Megator Define Transcriptionally Active Regions in the Genome
- The Genetic Interpretation of Area under the ROC Curve in Genomic Profiling
- Nucleoporins and Transcription: New Connections, New Questions