Limb development genes underlie variation in human fingerprint patterns.

Fingerprints are of long-standing practical and cultural interest, but little is known about the mechanisms that underlie their variation. Using genome-wide scans in Han Chinese cohorts, we identified 18 loci associated with fingerprint type across the digits, including a genetic basis for the long-recognized "pattern-block" correlations among the middle three digits. In particular, we identified a variant near EVI1 that alters regulatory activity and established a role for EVI1 in dermatoglyph patterning in mice. Dynamic EVI1 expression during human development supports its role in shaping the limbs and digits, rather than influencing skin patterning directly. Trans-ethnic meta-analysis identified 43 fingerprint-associated loci, with nearby genes being strongly enriched for general limb development pathways. We also found that fingerprint patterns were genetically correlated with hand proportions. Taken together, these findings support the key role of limb development genes in influencing the outcome of fingerprint patterning.

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.

TWIST1 interacts with ß/δ-catenins during neural tube development and regulates fate transition in cranial neural crest cells.

Cell fate determination is a necessary and tightly regulated process for producing different cell types and structures during development. Cranial neural crest cells (CNCCs) are unique to vertebrate embryos and emerge from the neural plate borders into multiple cell lineages that differentiate into bone, cartilage, neurons and glial cells. We have previously reported that Irf6 genetically interacts with Twist1 during CNCC-derived tissue formation. Here, we have investigated the mechanistic role of Twist1 and Irf6 at early stages of craniofacial development. Our data indicate that TWIST1 is expressed in endocytic vesicles at the apical surface and interacts with ß/δ-catenins during neural tube closure, and Irf6 is involved in defining neural fold borders by restricting AP2α expression. Twist1 suppresses Irf6 and other epithelial genes in CNCCs during the epithelial-to-mesenchymal transition (EMT) process and cell migration. Conversely, a loss of Twist1 leads to a sustained expression of epithelial and cell adhesion markers in migratory CNCCs. Disruption of TWIST1 phosphorylation in vivo leads to epidermal blebbing, edema, neural tube defects and CNCC-derived structural abnormalities. Altogether, this study describes a previously uncharacterized function of mammalian Twist1 and Irf6 in the neural tube and CNCCs, and provides new target genes for Twist1 that are involved in cytoskeletal remodeling.

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.

Pleiotropy method reveals genetic overlap between orofacial clefts at multiple novel loci from GWAS of multi-ethnic trios.

Based on epidemiologic and embryologic patterns, nonsyndromic orofacial clefts- the most common craniofacial birth defects in humans- are commonly categorized into cleft lip with or without cleft palate (CL/P) and cleft palate alone (CP), which are traditionally considered to be etiologically distinct. However, some evidence of shared genetic risk in IRF6, GRHL3 and ARHGAP29 regions exists; only FOXE1 has been recognized as significantly associated with both CL/P and CP in genome-wide association studies (GWAS). We used a new statistical approach, PLACO (pleiotropic analysis under composite null), on a combined multi-ethnic GWAS of 2,771 CL/P and 611 CP case-parent trios. At the genome-wide significance threshold of 5 × 10-8, PLACO identified 1 locus in 1q32.2 (IRF6) that appears to increase risk for one OFC subgroup but decrease risk for the other. At a suggestive significance threshold of 10-6, we found 5 more loci with compelling candidate genes having opposite effects on CL/P and CP: 1p36.13 (PAX7), 3q29 (DLG1), 4p13 (LIMCH1), 4q21.1 (SHROOM3) and 17q22 (NOG). Additionally, we replicated the recognized shared locus 9q22.33 (FOXE1), and identified 2 loci in 19p13.12 (RAB8A) and 20q12 (MAFB) that appear to influence risk of both CL/P and CP in the same direction. We found locus-specific effects may vary by racial/ethnic group at these regions of genetic overlap, and failed to find evidence of sex-specific differences. We confirmed shared etiology of the two OFC subtypes comprising CL/P, and additionally found suggestive evidence of differences in their pathogenesis at 2 loci of genetic overlap. Our novel findings include 6 new loci of genetic overlap between CL/P and CP; 3 new loci between pairwise OFC subtypes; and 4 loci not previously implicated in OFCs. Our in-silico validation showed PLACO is robust to subtype-specific effects, and can achieve massive power gains over existing approaches for identifying genetic overlap between disease subtypes. In summary, we found suggestive evidence for new genetic regions and confirmed some recognized OFC genes either exerting shared risk or with opposite effects on risk to OFC subtypes.

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.

Indirect Correlative Light and Electron Microscopy (iCLEM): A Novel Pipeline for Multiscale Quantification of Structure From Molecules to Organs.

Correlative light and electron microscopy (CLEM) methods are powerful methods that combine molecular organization (from light microscopy) with ultrastructure (from electron microscopy). However, CLEM methods pose high cost/difficulty barriers to entry and have very low experimental throughput. Therefore, we have developed an indirect correlative light and electron microscopy (iCLEM) pipeline to sidestep the rate-limiting steps of CLEM (i.e., preparing and imaging the same samples on multiple microscopes) and correlate multiscale structural data gleaned from separate samples imaged using different modalities by exploiting biological structures identifiable by both light and electron microscopy as intrinsic fiducials. We demonstrate here an application of iCLEM, where we utilized gap junctions and mechanical junctions between muscle cells in the heart as intrinsic fiducials to correlate ultrastructural measurements from transmission electron microscopy (TEM), and focused ion beam scanning electron microscopy (FIB-SEM) with molecular organization from confocal microscopy and single molecule localization microscopy (SMLM). We further demonstrate how iCLEM can be integrated with computational modeling to discover structure-function relationships. Thus, we present iCLEM as a novel approach that complements existing CLEM methods and provides a generalizable framework that can be applied to any set of imaging modalities, provided suitable intrinsic fiducials can be identified.

Emergent activity, heterogeneity, and robustness in a calcium feedback model of the sinoatrial node.

The sinoatrial node (SAN) is the primary pacemaker of the heart. SAN activity emerges at an early point in life and maintains a steady rhythm for the lifetime of the organism. The ion channel composition and currents of SAN cells can be influenced by a variety of factors. Therefore, the emergent activity and long-term stability imply some form of dynamical feedback control of SAN activity. We adapt a recent feedback model-previously utilized to describe control of ion conductances in neurons-to a model of SAN cells and tissue. The model describes a minimal regulatory mechanism of ion channel conductances via feedback between intracellular calcium and an intrinsic target calcium level. By coupling a SAN cell to the calcium feedback model, we show that spontaneous electrical activity emerges from quiescence and is maintained at steady state. In a 2D SAN tissue model, spatial variability in intracellular calcium targets lead to significant, self-organized heterogeneous ion channel expression and calcium transients throughout the tissue. Furthermore, multiple pacemaking regions appear, which interact and lead to time-varying cycle length, demonstrating that variability in heart rate is an emergent property of the feedback model. Finally, we demonstrate that the SAN tissue is robust to the silencing of leading cells or ion channel knockouts. Thus, the calcium feedback model can reproduce and explain many fundamental emergent properties of activity in the SAN that have been observed experimentally based on a minimal description of intracellular calcium and ion channel regulatory networks.

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.

Diversity of cells and signals in the cardiovascular system.

This white paper is the outcome of the seventh UC Davis Cardiovascular Research Symposium on Systems Approach to Understanding Cardiovascular Disease and Arrhythmia. This biannual meeting aims to bring together leading experts in subfields of cardiovascular biomedicine to focus on topics of importance to the field. The theme of the 2022 Symposium was 'Cell Diversity in the Cardiovascular System, cell-autonomous and cell-cell signalling'. Experts in the field contributed their experimental and mathematical modelling perspectives and discussed emerging questions, controversies, and challenges in examining cell and signal diversity, co-ordination and interrelationships involved in cardiovascular function. This paper originates from the topics of formal presentations and informal discussions from the Symposium, which aimed to develop a holistic view of how the multiple cell types in the cardiovascular system integrate to influence cardiovascular function, disease progression and therapeutic strategies. The first section describes the major cell types (e.g. cardiomyocytes, vascular smooth muscle and endothelial cells, fibroblasts, neurons, immune cells, etc.) and the signals involved in cardiovascular function. The second section emphasizes the complexity at the subcellular, cellular and system levels in the context of cardiovascular development, ageing and disease. Finally, the third section surveys the technological innovations that allow the interrogation of this diversity and advancing our understanding of the integrated cardiovascular function and dysfunction.

Rare variant modifier analysis identifies variants in SEC24D associated with orofacial cleft subtypes.

As one of the most common structural birth defects, orofacial clefts (OFCs) have been studied for decades, and recent studies have demonstrated that there are genetic differences between the different phenotypic presentations of OFCs. However, the contribution of rare genetic variation genome-wide to different subtypes of OFCs has been understudied, with most studies focusing on common genetic variation or rare variation within targeted regions of the genome. Therefore, we used whole-genome sequencing data from the Gabriella Miller Kids First Pediatric Research Program to conduct a gene-based burden analysis to test for genetic modifiers of cleft lip (CL) vs cleft lip and palate (CLP). We found that there was a significantly increased burden of rare variants in SEC24D in CL cases compared to CLP cases (p = 6.86 [Formula: see text] 10-7). Of the 15 variants within SEC24D, 53.3% were synonymous, but overlapped a known craniofacial enhancer. We then tested whether these variants could alter predicted transcription factor binding sites (TFBS), and found that the rare alleles destroyed binding sites for 9 transcription factors (TFs), including Pax1 (p = 0.0009), and created binding sites for 23 TFs, including Pax6 (p = 6.12 [Formula: see text] 10-5) and Pax9 (p = 0.0001), which are known to be involved in normal craniofacial development, suggesting a potential mechanism by which these synonymous variants could have a functional impact. Overall, this study indicates that rare genetic variation may contribute to the phenotypic heterogeneity of OFCs and suggests that regulatory variation may also contribute and warrant further investigation in future studies of genetic variants controlling risk to OFC.

Genome-wide Enrichment of De Novo Coding Mutations in Orofacial Cleft Trios.

Although de novo mutations (DNMs) are known to increase an individual's risk of congenital defects, DNMs have not been fully explored regarding orofacial clefts (OFCs), one of the most common human birth defects. Therefore, whole-genome sequencing of 756 child-parent trios of European, Colombian, and Taiwanese ancestry was performed to determine the contributions of coding DNMs to an individual's OFC risk. Overall, we identified a significant excess of loss-of-function DNMs in genes highly expressed in craniofacial tissues, as well as genes associated with known autosomal dominant OFC syndromes. This analysis also revealed roles for zinc-finger homeobox domain and SOX2-interacting genes in OFC etiology.

Modeling incomplete penetrance in long QT syndrome type 3 through ion channel heterogeneity: an in silico population study.

Many cardiac diseases are characterized by an increased late sodium current, including heart failure, hypertrophic cardiomyopathy, and inherited long QT syndrome type 3 (LQT3). The late sodium current in LQT3 is caused by a gain-of-function mutation in the voltage-gated sodium channel Nav1.5. Despite a well-defined genetic cause of LQT3, treatment remains inconsistent because of incomplete penetrance of the mutation and variability of antiarrhythmic efficacy. Here, we investigate the relationship between LQT3-associated mutation incomplete penetrance and variability in ion channel expression, simulating a population of 1,000 individuals with the O'Hara-Rudy model of the human ventricular myocyte. We first simulate healthy electrical activity (i.e., in the absence of a mutation) and then incorporate heterozygous expression for three LQT3-associated mutations (Y1795C, I1768V, and ΔKPQ), to directly compare the effects of each mutation on individuals across a diverse population. For all mutations, we find that susceptibility, defined by either the presence of an early afterdepolarization (EAD) or prolonged action potential duration (APD), primarily depends on the balance between the conductance of IKr and INa, for which individuals with a higher IKr-to-INa ratio are less susceptible. Furthermore, we find distinct differences across the population, observing individuals susceptible to zero, one, two, or all three mutations. Individuals tend to be less susceptible with an appropriate balance of repolarizing currents, typically via increased IKs or IK1. Interestingly, the more critical repolarizing current is mutation specific. We conclude that the balance between key currents plays a significant role in mutant-specific presentation of the disease phenotype in LQT3.NEW & NOTEWORTHY An in silico population approach investigates the relationship between variability in ion channel expression and gain-of-function mutations in the voltage-gated sodium channel associated with the congenital disorder long QT syndrome type 3 (LQT3). We find that ion channel variability can contribute to incomplete penetrance of the mutation, with mutant-specific differences in ion channel conductances leading to susceptibility to proarrhythmic action potential duration prolongation or early afterdepolarizations.

Long-term changes in heart rate and electrical remodeling contribute to alternans formation in heart failure: a patient-specific in silico study.

Individuals with chronic heart failure (CHF) have an increased risk of ventricular arrhythmias, which has been linked to pathological cellular remodeling and may also be mediated by changes in heart rate. Heart rate typically fluctuates on a timescale ranging from seconds to hours, termed heart rate variability (HRV). This variability is reduced in CHF, and this HRV reduction is associated with a greater risk for arrhythmias. Furthermore, variations in heart rate influence the formation of proarrhythmic alternans, a beat-to-beat alternation in the action potential duration (APD), or intracellular calcium (Ca). In this study, we investigate how long-term changes in heart rate and electrical remodeling associated with CHF influence alternans formation. We measure key statistical properties of the RR-interval sequences from ECGs of individuals with normal sinus rhythm (NSR) and CHF. Patient-specific RR-interval sequences and synthetic sequences (randomly generated to mimicking these statistical properties) are used as the pacing protocol for a discrete time-coupled map model that governs APD and intracellular Ca handling of a single cardiac myocyte, modified to account for pathological electrical remodeling in CHF. Patient-specific simulations show that beat-to-beat differences in APD vary temporally in both populations, with alternans formation more prevalent in CHF. Parameter studies using synthetic sequences demonstrate that increasing the autocorrelation time or mean RR-interval reduces APD alternations, whereas increasing the RR-interval standard deviation leads to higher alternans magnitudes. Importantly, we find that although both the CHF-associated changes in heart rate and electrical remodeling influence alternans formation, variations in heart rate may be more influential.NEW & NOTEWORTHY Using patient-specific data, we show that both the changes in heart rate and electrical remodeling associated with chronic heart failure influence the formation of proarrhythmic alternans in the heart.

Rare variants found in clinical gene panels illuminate the genetic and allelic architecture of orofacial clefting.

PURPOSE: Orofacial clefts (OFCs) are common birth defects including cleft lip, cleft lip and palate, and cleft palate. OFCs have heterogeneous etiologies, complicating clinical diagnostics because it is not always apparent if the cause is Mendelian, environmental, or multifactorial. Sequencing is not currently performed for isolated or sporadic OFCs; therefore, we estimated the diagnostic yield for 418 genes in 841 cases and 294 controls. METHODS: We evaluated 418 genes using genome sequencing and curated variants to assess their pathogenicity using American College of Medical Genetics criteria. RESULTS: 9.04% of cases and 1.02% of controls had "likely pathogenic" variants (P < .0001), which was almost exclusively driven by heterozygous variants in autosomal genes. Cleft palate (17.6%) and cleft lip and palate (9.09%) cases had the highest yield, whereas cleft lip cases had a 2.80% yield. Out of 39 genes with likely pathogenic variants, 9 genes, including CTNND1 and IRF6, accounted for more than half of the yield (4.64% of cases). Most variants (61.8%) were "variants of uncertain significance", occurring more frequently in cases (P = .004), but no individual gene showed a significant excess of variants of uncertain significance. CONCLUSION: These results underscore the etiological heterogeneity of OFCs and suggest sequencing could reduce the diagnostic gap in OFCs.

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.

Ephaptic Coupling as a Resolution to the Paradox of Action Potential Wave Speed and Discordant Alternans Spatial Scales in the Heart.

Previous computer simulations have suggested that existing models of action potential wave propagation in the heart are not consistent with observed wave propagation behavior. Specifically, computer models cannot simultaneously reproduce the rapid wave speeds and small spatial scales of discordant alternans patterns measured experimentally in the same simulation. The discrepancy is important, because discordant alternans can be a key precursor to the development of abnormal and dangerous rapid rhythms in the heart. In this Letter, we show that this paradox can be resolved by allowing so-called ephaptic coupling to play a primary role in wave front propagation in place of conventional gap-junction coupling. With this modification, physiological wave speeds and small discordant alternans spatial scales both occur with gap-junction resistance values that are more in line with those observed in experiments. Our theory thus also provides support to the hypothesis that ephaptic coupling plays an important role in normal wave propagation.

Rare variants found in multiplex families with orofacial clefts: Does expanding the phenotype make a difference?

Exome sequencing (ES) is now a relatively straightforward process to identify causal variants in Mendelian disorders. However, the same is not true for ES in families where the inheritance patterns are less clear, and a complex etiology is suspected. Orofacial clefts (OFCs) are highly heritable birth defects with both Mendelian and complex etiologies. The phenotypic spectrum of OFCs may include overt clefts and several subclinical phenotypes, such as discontinuities in the orbicularis oris muscle (OOM) in the upper lip, velopharyngeal insufficiency (VPI), microform clefts or bifid uvulas. We hypothesize that expanding the OFC phenotype to include these phenotypes can clarify inheritance patterns in multiplex families, making them appear more Mendelian. We performed exome sequencing to find rare, likely causal genetic variants in 31 multiplex OFC families, which included families with multiple individuals with OFCs and individuals with subclinical phenotypes. We identified likely causal variants in COL11A2, IRF6, SHROOM3, SMC3, TBX3, and TP63 in six families. Although we did not find clear evidence supporting the subclinical phenotype hypothesis, our findings support a role for rare variants in the etiology of OFCs.

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.