Decoding the Human Face: Progress and Challenges in Understanding the Genetics of Craniofacial Morphology.

Variations in the form of the human face, which plays a role in our individual identities and societal interactions, have fascinated scientists and artists alike. Here, we review our current understanding of the genetics underlying variation in craniofacial morphology and disease-associated dysmorphology, synthesizing decades of progress on Mendelian syndromes in addition to more recent results from genome-wide association studies of human facial shape and disease risk. We also discuss the various approaches used to phenotype and quantify facial shape, which are of particular importance due to the complex, multipartite nature of the craniofacial form. We close by discussing how experimental studies have contributed and will further contribute to our understanding of human genetic variation and then proposing future directions and applications for the field.

Genome scan reveals several loci associated with torus palatinus.

OBJECTIVE: Torus Palatinus (TP) is a common trait with an unclear aetiology. Although prior studies suggest a hereditary component, the genetic factors that influence TP risk remain unknown. The purpose of this study is to identify genetic variants associated with TP. MATERIALS AND METHODS: We assessed the TP status of 829 individuals from various ancestral backgrounds using 3D palate scans. We then carried out a genome-wide association study (GWAS) to identify common variants associated with TP. We also performed gene-based tests across the exome to investigate the role of low-frequency coding variants. RESULTS: Our GWAS did not identify any genome-wide significant signals but identified suggestive associations including hits on chromosomes 2, 5 and 17 with p-values less than 5 × 10-6. Candidate genes at these suggestive loci have been implicated in normal-range craniofacial features, syndromes with facial and oral malformations, and bone density. We did not find evidence that low-frequency coding variants influence TP risk. In addition, we failed to replicate associations identified in prior genetic studies of TP. CONCLUSION: These findings suggest that multiple genes likely influence the development of TP. Independent replication will be required to confirm our suggestive associations.

3D facial phenotyping by biometric sibling matching used in contemporary genomic methodologies.

The analysis of contemporary genomic data typically operates on one-dimensional phenotypic measurements (e.g. standing height). Here we report on a data-driven, family-informed strategy to facial phenotyping that searches for biologically relevant traits and reduces multivariate 3D facial shape variability into amendable univariate measurements, while preserving its structurally complex nature. We performed a biometric identification of siblings in a sample of 424 children, defining 1,048 sib-shared facial traits. Subsequent quantification and analyses in an independent European cohort (n = 8,246) demonstrated significant heritability for a subset of traits (0.17-0.53) and highlighted 218 genome-wide significant loci (38 also study-wide) associated with facial variation shared by siblings. These loci showed preferential enrichment for active chromatin marks in cranial neural crest cells and embryonic craniofacial tissues and several regions harbor putative craniofacial genes, thereby enhancing our knowledge on the genetic architecture of normal-range facial variation.

Genome scans of facial features in East Africans and cross-population comparisons reveal novel associations.

Facial morphology is highly variable, both within and among human populations, and a sizable portion of this variation is attributable to genetics. Previous genome scans have revealed more than 100 genetic loci associated with different aspects of normal-range facial variation. Most of these loci have been detected in Europeans, with few studies focusing on other ancestral groups. Consequently, the degree to which facial traits share a common genetic basis across diverse sets of humans remains largely unknown. We therefore investigated the genetic basis of facial morphology in an East African cohort. We applied an open-ended data-driven phenotyping approach to a sample of 2,595 3D facial images collected on Tanzanian children. This approach segments the face into hierarchically arranged, multivariate features that capture the shape variation after adjusting for age, sex, height, weight, facial size and population stratification. Genome scans of these multivariate shape phenotypes revealed significant (p < 2.5 × 10-8) signals at 20 loci, which were enriched for active chromatin elements in human cranial neural crest cells and embryonic craniofacial tissue, consistent with an early developmental origin of the facial variation. Two of these associations were in highly conserved regions showing craniofacial-specific enhancer activity during embryological development (5q31.1 and 12q21.31). Six of the 20 loci surpassed a stricter threshold accounting for multiple phenotypes with study-wide significance (p < 6.25 × 10-10). Cross-population comparisons indicated 10 association signals were shared with Europeans (seven sharing the same associated SNP), and facilitated fine-mapping of causal variants at previously reported loci. Taken together, these results may point to both shared and population-specific components to the genetic architecture of facial variation.

Genome-wide association study of multiethnic nonsyndromic orofacial cleft families identifies novel loci specific to family and phenotypic subtypes.

Nonsyndromic orofacial clefts (OFCs) are among the most common craniofacial birth defects worldwide, and known to exhibit phenotypic and genetic heterogeneity. Cleft lip plus cleft palate (CLP) and cleft lip only (CL) are commonly combined together as one phenotype (CL/P), separately from cleft palate alone. In comparison, our study analyzes CL and CLP separately. A sample of 2218 CL and CLP cases, 4537 unaffected relatives of cases, and 2673 pure controls with no family history of OFC were selected from the Pittsburgh Orofacial Cleft (Pitt-OFC) multiethnic study.genome-wide association studies were run for seven specific phenotypes created based on the cleft type(s) observed within these families, as well as the combined CL/P phenotype. Five novel genome-wide significant associations, 3q29 (rs62284390), 5p13.2 (rs609659), 7q22.1 (rs6465810), 19p13.3 (rs628271), and 20q13.33 (rs2427238), and nine associations (p ≤ 1.0E-05) within previously confirmed OFC loci-PAX7, IRF6, FAM49A, DCAF4L2, 8q24.21, ARID3B, NTN1, TANC2 and the WNT9B:WNT3 gene cluster-were observed. We also found that single nucleotide polymorphisms within a subset of the associated loci, both previously known and novel, differ substantially in terms of their effects across cleft- or family-specific phenotypes, indicating not only etiologic differences between CL and CLP, but also genetic heterogeneity within each of the two OFC subtypes.

Using data-driven phenotyping to investigate the impact of sex on 3D human facial surface morphology.

The effects of sex on human facial morphology have been widely documented. Because sexual dimorphism is relevant to a variety of scientific and applied disciplines, it is imperative to have a complete and accurate account of how and where male and female faces differ. We apply a comprehensive facial phenotyping strategy to a large set of existing 3D facial surface images. We investigate facial sexual dimorphism in terms of size, shape, and shape variance. We also assess the ability to correctly assign sex based on shape, both for the whole face and for subregions. We applied a predefined data-driven segmentation to partition the 3D facial surfaces of 2446 adults into 63 hierarchically linked regions, ranging from global (whole face) to highly localized subparts. Each facial region was then analyzed with spatially dense geometric morphometrics. To describe the major modes of shape variation, principal components analysis was applied to the Procrustes aligned 3D points comprising each of the 63 facial regions. Both nonparametric and permutation-based statistics were then used to quantify the facial size and shape differences and visualizations were generated. Males were significantly larger than females for all 63 facial regions. Statistically significant sex differences in shape were also seen in all regions and the effects tended to be more pronounced for the upper lip and forehead, with more subtle changes emerging as the facial regions became more granular. Males also showed greater levels of shape variance, with the largest effect observed for the central forehead. Classification accuracy was highest for the full face (97%), while most facial regions showed an accuracy of 75% or greater. In summary, sex differences in both size and shape were present across every part of the face. By breaking the face into subparts, some shape differences emerged that were not apparent when analyzing the face as a whole. The increase in facial shape variance suggests possible evolutionary origins and may offer insights for understanding congenital facial malformations. Our classification results indicate that a high degree of accuracy is possible with only parts of the face, which may have implications for biometrics applications.

Genetic Analyses of Enamel Hypoplasia in Multiethnic Cohorts.

Enamel hypoplasia causes reduction in the thickness of affected enamel and is one of the most common dental anomalies. This defect is caused by environmental and/or genetic factors that interfere with tooth formation, emphasizing the importance of investigating enamel hypoplasia on an epidemiological and genetic level. A genome-wide association of enamel hypoplasia was performed in multiple cohorts, overall comprising 7,159 individuals ranging in age from 7-82 years. Mixed-models were used to test for genetic association while simultaneously accounting for relatedness and genetic population structure. Meta-analysis was then performed. More than 5 million single-nucleotide polymorphisms were tested in individual cohorts. Analyses of the individual cohorts and meta-analysis identified association signals close to genome-wide significance (P < 510-8), and many suggestive association signals (510-8 < P < 510-6) near genes with plausible roles in tooth/enamel development. The strongest association signal (P = 1.5710-9) was observed near BMP2K in one of the individual cohorts. Additional suggestive signals were observed near genes with plausible roles in tooth development in the meta-analysis, such as SLC4A4 which can influence enamel hypoplasia. Additional human genetic studies are needed to replicate these results and functional studies in model systems are needed to validate our findings.

Correction: Genetic variants in pachyonychia congenita-associated keratins increase susceptibility to tooth decay.

[This corrects the article DOI: 10.1371/journal.pgen.1007168.].

Decreased DNA Methylation of RGMA is Associated with Intracranial Hypertension After Severe Traumatic Brain Injury: An Exploratory Epigenome-Wide Association Study.

BACKGROUND: Cerebral edema and intracranial hypertension are major contributors to unfavorable prognosis in traumatic brain injury (TBI). Local epigenetic changes, particularly in DNA methylation, may influence gene expression and thus host response/secondary injury after TBI. It remains unknown whether DNA methylation in the central nervous system is associated with cerebral edema severity or intracranial hypertension post TBI. We sought to identify epigenome-wide DNA methylation patterns associated with these forms of secondary injury after TBI. METHODS: We obtained genome-wide DNA methylation profiles of DNA extracted from ventricular cerebrospinal fluid samples at three different postinjury time points from a prospective cohort of patients with severe TBI (n = 89 patients, 254 samples). Cerebral edema and intracranial pressure (ICP) measures were clustered to generate composite end points of cerebral edema and ICP severity. We performed an unbiased epigenome-wide association study (EWAS) to test associations between DNA methylation at 419,895 cytosine-phosphate-guanine (CpG) sites and cerebral edema/ICP severity categories. Given inflated p values, we conducted permutation tests for top CpG sites to filter out potential false discoveries. RESULTS: Our data-driven hierarchical clustering across six cerebral edema and ICP measures identified two groups differing significantly in ICP based on the EWAS-identified CpG site cg22111818 in RGMA (Repulsive guidance molecule A, permutation p = 4.20 × 10-8). At 3-4 days post TBI, patients with severe intracranial hypertension had significantly lower levels of methylation at cg22111818. CONCLUSIONS: We report a novel potential relationship between intracranial hypertension after TBI and an acute, nonsustained reduction in DNA methylation at cg22111818 in the RGMA gene. To our knowledge, this is the largest EWAS in severe TBI. Our findings are further strengthened by previous findings that RGMA modulates axonal repair in other central nervous system disorders, but a role in intracranial hypertension or TBI has not been previously identified. Additional work is warranted to validate and extend these findings, including assessment of its possible role in risk stratification, identification of novel druggable targets, and ultimately our ability to personalize therapy in TBI.

Heritability Analysis in Twins Indicates a Genetic Basis for Velopharyngeal Morphology.

The velopharyngeal mechanism is comprised of several muscular components that act in a coordinated manner to control airflow through the nose and mouth. Proper velopharyngeal function is essential for normal speech, swallowing, and breathing. The genetic basis of normal-range velopharyngeal morphology is poorly understood. The purpose of this study was to estimate the heritability of velopharyngeal dimensions.We measured five velopharyngeal variables (velar length, velar thickness, effective velar length, levator muscle length and pharyngeal depth) from MRIs of 155 monozygotic and 208 dizygotic twin pairs and then calculated heritability for these traits using a structural equation modeling approach.The heritability estimates were statistically significant (95% confidence intervals excluded zero) and ranged from 0.19 to 0.46. There was also evidence of significant genetic correlations between pairs of traits, pointing to the influence of common genetic effects.These results indicate that genetic factors influence variation in clinically relevant velopharyngeal structures.

Genome-Wide Association Analysis of Longitudinal Bone Mineral Content Data From the Iowa Bone Development Study.

The foundation for osteoporosis risk is, in part, established during childhood, adolescence, and young adulthood, all periods of development when bone mass is acquired rapidly. The relative quantity of bone mass accrued is influenced by both lifestyle and genetic factors, although the genetic component is not yet well understood. The purpose of this study was to use a genome-wide association (GWA) analysis to discover single nucleotide polymorphisms (SNPs) associated with: (1) the sex-specific hip bone mineral content at approximately the age of 19 when the amount of bone accrued is near its peak; and (2) the sex-specific rate of hip bone mineral content accrual during the adolescent growth spurt. The Iowa Bone Development Study, a longitudinal cohort study exploring bone health in children, adolescents, and young adults was the source of data. From this cohort, n = 364 (190 females, 174 males) participants were included in GWA analyses to address (1) and n = 258 participants (125 females and 133 males) were included in GWA analyses to address (2). Twenty SNPS were detected having p < 1.0 × 10-5. Of most biologic relevance were 2 suggestive SNPs (rs2051756 and rs2866908) detected in an intron of the DKK2 gene through the GWA analysis that explored peak bone mass in females.

Genetic variants in pachyonychia congenita-associated keratins increase susceptibility to tooth decay.

Pachyonychia congenita (PC) is a cutaneous disorder primarily characterized by nail dystrophy and painful palmoplantar keratoderma. PC is caused by mutations in KRT6A, KRT6B, KRT6C, KRT16, and KRT17, a set of keratin genes expressed in the nail bed, palmoplantar epidermis, oral mucosal epithelium, hair follicle and sweat gland. RNA-seq analysis revealed that all PC-associated keratins (except for Krt6c that does exist in the mouse genome) are expressed in the mouse enamel organ. We further demonstrated that these keratins are produced by ameloblasts and are incorporated into mature human enamel. Using genetic and intraoral examination data from 573 adults and 449 children, we identified several missense polymorphisms in KRT6A, KRT6B and KRT6C that lead to a higher risk for dental caries. Structural analysis of teeth from a PC patient carrying a p.Asn171Lys substitution in keratin-6a (K6a) revealed disruption of enamel rod sheaths resulting in altered rod shape and distribution. Finally, this PC-associated substitution as well as more frequent caries-associated SNPs, found in two of the KRT6 genes, that result in p.Ser143Asn substitution (rs28538343 in KRT6B and rs151117600 in KRT6C), alter the assembly of K6 filaments in ameloblast-like cells. These results identify a new set of keratins involved in tooth enamel formation, distinguish novel susceptibility loci for tooth decay and reveal additional clinical features of pachyonychia congenita.

Investigating the shared genetics of non-syndromic cleft lip/palate and facial morphology.

There is increasing evidence that genetic risk variants for non-syndromic cleft lip/palate (nsCL/P) are also associated with normal-range variation in facial morphology. However, previous analyses are mostly limited to candidate SNPs and findings have not been consistently replicated. Here, we used polygenic risk scores (PRS) to test for genetic overlap between nsCL/P and seven biologically relevant facial phenotypes. Where evidence was found of genetic overlap, we used bidirectional Mendelian randomization (MR) to test the hypothesis that genetic liability to nsCL/P is causally related to implicated facial phenotypes. Across 5,804 individuals of European ancestry from two studies, we found strong evidence, using PRS, of genetic overlap between nsCL/P and philtrum width; a 1 S.D. increase in nsCL/P PRS was associated with a 0.10 mm decrease in philtrum width (95% C.I. 0.054, 0.146; P = 2x10-5). Follow-up MR analyses supported a causal relationship; genetic variants for nsCL/P homogeneously cause decreased philtrum width. In addition to the primary analysis, we also identified two novel risk loci for philtrum width at 5q22.2 and 7p15.2 in our Genome-wide Association Study (GWAS) of 6,136 individuals. Our results support a liability threshold model of inheritance for nsCL/P, related to abnormalities in development of the philtrum.

The Influence of Sex and Ancestry on Three-Dimensional Palate Shape.

ABSTRACT: Modern human palate shape has been reported to vary by sex and ancestry, but limitations in the methods used to quantify shape and in population coverage have led to inconsistent findings. In the present study, the authors aim to characterize the effects of sex and ancestry on normal-range three-dimensional palate shape through landmark-based morphometrics.Three-dimensional digital dental casts were obtained and landmarked from 794 adults of European (n = 429), African (n = 295), and East Asian (n = 70) ancestry. Principal component analysis was conducted to identify patterns of shape variation present in our cohort, and canonical variates analysis was performed to test for shape differences between sexes and ancestries.Principal component analysis showed that 3 principal components, explaining 76.52% of variance, linked higher palatal vault with either a relative reduction in anteroposterior or mediolateral dimensions. Canonical variates analysis showed that males had wider and shorter palates with more posteriorly located maximum vault depth than females. Individuals of African ancestry, having higher vaults with more posteriorly located maximal depths, also had wider and shorter palates, whereas individuals of European ancestry had narrower and longer palates with more anteriorly located maximum vault depths. Individuals of East Asian ancestry showed the shallowest vaults.It was found that both sex and ancestry influence palate shape, suggesting a possible genetic component underlying this variation. Additionally, our findings indicate that vault height tends to co-vary with anteroposterior or mediolateral dimensions. Further investigation of these morphological patterns may shed light on possible links to common congenital anomalies such as orofacial clefting.

Parents of Children With Nonsyndromic Orofacial Clefting Show Altered Palate Shape.

OBJECTIVE: The unaffected relatives of individuals with nonsyndromic orofacial clefts have been shown to exhibit subtle craniofacial differences compared with the general population. Here, we investigate whether these morphological differences extend to the shape of the palate. DESIGN: We conducted a geometric morphometric analysis to compare palate shape in the clinically unaffected parents of children with nonsyndromic cleft lip with or without cleft palate and adult controls of European, Asian, and African ancestry. We conducted pairwise group comparisons using canonical variates analysis, and then confirmed and characterized findings of shape differences using Euclidean distance matrix analysis. RESULTS: Significant differences in palate shape were detected in unaffected mothers (but not fathers) compared to demographically matched controls. The differences in shape were ancestry-specific; mothers of Asian-derived and African-derived ancestry showed wider and shorter palates with higher posterior palatal vaults, while mothers of European-derived ancestry showed narrower palates with higher anterior palatal vaults. CONCLUSIONS: Our findings suggest that altered palate shape is a subclinical phenotypic feature, which may be indicative of elevated orofacial cleft risk. The risk phenotype varied by sex and ancestry, suggesting possible etiologic heterogeneity among demographic groups. Understanding the genetic basis of these informative palate shape traits may reveal new genes and pathways relevant to nonsyndromic orofacial clefting.

Replication of GWAS significant loci in a sub-Saharan African Cohort with early childhood caries: a pilot study.

BACKGROUND: Early childhood caries (ECC) is a rapidly progressing form of dental infection and a significant public health problem, especially among socially and economically disadvantaged populations. This study aimed to assess the risk factors for ECC among a cohort of Sub-Saharan African children and to determine the role of genetics in the etiology of ECC. METHODS: A sample of 691 children (338 with ECC, 353 without ECC, age < 6 years) was recruited from schools in Lagos, Nigeria. Socio-demographic, dental services utilization and infant dietary data were obtained with interviewer-administered questionnaire. Oral examination was conducted using the WHO oral health diagnostic criteria. Saliva samples were collected from the children for genetic analysis. Single nucleotide polymorphisms were selected from previous study for genotyping. Genetic association analyses to investigate the role of genetics in the etiology of ECC was done. Bivariate comparisons and Multivariate logistic regression analyses were conducted to assess associations between ECC and predictor variables, p < 0.05. RESULTS: Of the 338 children with ECC, 64 (18.9%) had Severe-Early Childhood Caries (S-ECC). Children aged 48-59 months comprised the highest proportion of subjects with ECC (165; 48.8%) and S-ECC (24; 37.5%) while female subjects had higher dt (3.13 ± 2.56) and dmft values 3.27 ± 2.64. ECC was significantly more prevalent among children who were breastfed at night ≥ 12 months (OR 3.30; CI 0.39, 4.75), those with no previous dental visit (OR 1.71; CI 0.24, 2.77), those who used sweetened pacifiers (OR 1.85; CI 0.91, 3.79) and those who daily consumed sugar-sweetened drinks/snacks (OR 1.35; CI 0.09, 18.51). A suggestive increased risk for ECC (OR 1.26, p = 0. 0.0397) was observed for the genetic variant rs11239282 on chromosome 10. We also observed a suggestive reduced risk for ECC (OR 0.80, p = 0.03) for the rs131777 on chromosome 22. None of the genetic variants were significant after correction for multiple testing (Bonferroni p value p = 0.004). CONCLUSIONS: Prolonged night-time breastfeeding, poor utilization of dental services and daily consumption of sugar were risk factors for ECC. Larger sample size is needed to confirm the results of the genetic analysis and to conduct genome wide studies in order to discover new risk loci for ECC.

Genome-Wide Association Study (GWAS) of dental caries in diverse populations.

BACKGROUND: Dental caries is one of the most common chronic diseases and is influenced by a complex interplay of genetic and environmental factors. Most previous genetic studies of caries have focused on identifying genes that contribute to dental caries in specific ethnic groups, usually of European descent. METHODS: The aim of this study is to conduct a genome-wide association study (GWAS) to identify associations affecting susceptibility to caries in a large multiethnic population from Argentina, the Philippines, Guatemala, Hungary, and the USA, originally recruited for studies of orofacial clefts (POFC, N = 3686). Ages of the participants ranged from 2 to 12 years for analysis of the primary dentition, and 18-60 years for analysis of the permanent dentition. For each participant, dental caries was assessed by counts of decayed and filled teeth (dft/DFT) and genetic variants (single nucleotide polymorphisms, SNPs) were genotyped or imputed across the entire genome. Caries was analyzed separately for the primary and permanent dentitions, with age, gender, and presence/absence of any type of OFC treated as covariates. Efficient Mixed-Model Association eXpedited (EMMAX) was used to test genetic association, while simultaneously accounting for relatedness and stratification. RESULTS: We identified several suggestive loci (5 × 10-8 < P < 5 × 10-6) within or near genes with plausible biological roles for dental caries, including a cluster of taste receptor genes (TAS2R38, TAS2R3, TAS2R4, TASR25) on chromosome 7 for the permanent dentition analysis, and DLX3 and DLX4 on chromosome 17 for the primary dentition analysis. Genome-wide significant results were seen with SNPs in the primary dentition only; however, none of the identified genes near these variants have known roles in cariogenesis. CONCLUSION: The results of this study warrant further investigation and may lead to a better understanding of cariogenesis in diverse populations, and help to improve dental caries prediction, prevention, and/or treatment in future.

A fully adjusted two-stage procedure for rank-normalization in genetic association studies.

When testing genotype-phenotype associations using linear regression, departure of the trait distribution from normality can impact both Type I error rate control and statistical power, with worse consequences for rarer variants. Because genotypes are expected to have small effects (if any) investigators now routinely use a two-stage method, in which they first regress the trait on covariates, obtain residuals, rank-normalize them, and then use the rank-normalized residuals in association analysis with the genotypes. Potential confounding signals are assumed to be removed at the first stage, so in practice, no further adjustment is done in the second stage. Here, we show that this widely used approach can lead to tests with undesirable statistical properties, due to both combination of a mis-specified mean-variance relationship and remaining covariate associations between the rank-normalized residuals and genotypes. We demonstrate these properties theoretically, and also in applications to genome-wide and whole-genome sequencing association studies. We further propose and evaluate an alternative fully adjusted two-stage approach that adjusts for covariates both when residuals are obtained and in the subsequent association test. This method can reduce excess Type I errors and improve statistical power.

A systematic genetic analysis and visualization of phenotypic heterogeneity among orofacial cleft GWAS signals.

Phenotypic heterogeneity is a hallmark of complex traits, and genetic studies of such traits may focus on them as a single diagnostic entity or by analyzing specific components. For example, in orofacial clefting (OFC), three subtypes-cleft lip (CL), cleft lip and palate (CLP), and cleft palate (CP) have been studied separately and in combination. To further dissect the genetic architecture of OFCs and how a given associated locus may be contributing to distinct subtypes of a trait we developed a framework for quantifying and interpreting evidence of subtype-specific or shared genetic effects in complex traits. We applied this technique to create a "cleft map" of the association of 30 genetic loci with three OFC subtypes. In addition to new associations, we found loci with subtype-specific effects (e.g., GRHL3 [CP], WNT5A [CLP]), as well as loci associated with two or all three subtypes. We cross-referenced these results with mouse craniofacial gene expression datasets, which identified additional promising candidate genes. However, we found no strong correlation between OFC subtypes and expression patterns. In aggregate, the cleft map revealed that neither subtype-specific nor shared genetic effects operate in isolation in OFC architecture. Our approach can be easily applied to any complex trait with distinct phenotypic subgroups.

Consortium-based genome-wide meta-analysis for childhood dental caries traits.

Prior studies suggest dental caries traits in children and adolescents are partially heritable, but there has been no large-scale consortium genome-wide association study (GWAS) to date. We therefore performed GWAS for caries in participants aged 2.5-18.0 years from nine contributing centres. Phenotype definitions were created for the presence or absence of treated or untreated caries, stratified by primary and permanent dentition. All studies tested for association between caries and genotype dosage and the results were combined using fixed-effects meta-analysis. Analysis included up to 19 003 individuals (7530 affected) for primary teeth and 13 353 individuals (5875 affected) for permanent teeth. Evidence for association with caries status was observed at rs1594318-C for primary teeth [intronic within ALLC, odds ratio (OR) 0.85, effect allele frequency (EAF) 0.60, P 4.13e-8] and rs7738851-A (intronic within NEDD9, OR 1.28, EAF 0.85, P 1.63e-8) for permanent teeth. Consortium-wide estimated heritability of caries was low [h2 of 1% (95% CI: 0%: 7%) and 6% (95% CI 0%: 13%) for primary and permanent dentitions, respectively] compared with corresponding within-study estimates [h2 of 28% (95% CI: 9%: 48%) and 17% (95% CI: 2%: 31%)] or previously published estimates. This study was designed to identify common genetic variants with modest effects which are consistent across different populations. We found few single variants associated with caries status under these assumptions. Phenotypic heterogeneity between cohorts and limited statistical power will have contributed; these findings could also reflect complexity not captured by our study design, such as genetic effects which are conditional on environmental exposure.