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.

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.

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.

The impact of prenatal maternal stress due to potentially traumatic events on child temperament: A systematic review.

The objective of the current study was to complete a systematic review of the relationship between prenatal maternal stress due to potentially traumatic events (PTEs) and child temperament. Eligible studies through June 2020 were identified utilizing a search strategy in PubMed and PsycInfo. Included studies examined associations between prenatal maternal stress due to PTE and child temperament. Two independent coders extracted study characteristics and three coders assessed study quality. Of the 1969 identified studies, 20 met full inclusion criteria. Studies were classified on two dimensions: (1) disaster-related stress and (2) intimate partner violence during pregnancy. For disaster-related prenatal maternal stress, 75% (nine out of 12) of published reports found associations with increased child negative affectivity, 50% (five out of 10) also noted associations with lower effortful control/regulation, and 38% (three out of eight) found associations with lower positive affectivity. When considering prenatal intimate partner violence stress, 80% (four out of five) of published reports found associations with higher child negative affectivity, 67% (four out of six) found associations with lower effortful control/regulation, and 33% (one out of three) found associations with lower positive affectivity. Prenatal maternal stress due to PTEs may impact the offspring's temperament, especially negative affectivity. Mitigating the effects of maternal stress in pregnancy is needed in order to prevent adverse outcomes on the infant's socioemotional development.

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.

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 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.

Genome-wide interaction studies identify sex-specific risk alleles for nonsyndromic orofacial clefts.

Nonsyndromic cleft lip with or without cleft palate (NSCL/P) is the most common craniofacial birth defect in humans and is notable for its apparent sexual dimorphism where approximately twice as many males are affected as females. The sources of this disparity are largely unknown, but interactions between genetic and sex effects are likely contributors. We examined gene-by-sex (G × S) interactions in a worldwide sample of 2,142 NSCL/P cases and 1,700 controls recruited from 13 countries. First, we performed genome-wide joint tests of the genetic (G) and G × S effects genome-wide using logistic regression assuming an additive genetic model and adjusting for 18 principal components of ancestry. We further interrogated loci with suggestive results from the joint test ( p < 1.00 × 10 -5 ) by examining the G × S effects from the same model. Out of the 133 loci with suggestive results ( p < 1.00 × 10 -5 ) for the joint test, we observed one genome-wide significant G × S effect in the 10q21 locus (rs72804706; p = 6.69 × 10 -9 ; OR = 2.62 CI [1.89, 3.62]) and 16 suggestive G × S effects. At the intergenic 10q21 locus, the risk of NSCL/P is estimated to increase with additional copies of the minor allele for females, but the opposite effect for males. Our observation that the impact of genetic variants on NSCL/P risk differs for males and females may further our understanding of the genetic architecture of NSCL/P and the sex differences underlying clefts and other birth defects.

A Genome-wide Association Study of Nonsyndromic Cleft Palate Identifies an Etiologic Missense Variant in GRHL3.

Cleft palate (CP) is a common birth defect occurring in 1 in 2,500 live births. Approximately half of infants with CP have a syndromic form, exhibiting other physical and cognitive disabilities. The other half have nonsyndromic CP, and to date, few genes associated with risk for nonsyndromic CP have been characterized. To identify such risk factors, we performed a genome-wide association study of this disorder. We discovered a genome-wide significant association with a missense variant in GRHL3 (p.Thr454Met [c.1361C>T]; rs41268753; p = 4.08 × 10(-9)) and replicated the result in an independent sample of case and control subjects. In both the discovery and replication samples, rs41268753 conferred increased risk for CP (OR = 8.3, 95% CI 4.1-16.8; OR = 2.16, 95% CI 1.43-3.27, respectively). In luciferase transactivation assays, p.Thr454Met had about one-third of the activity of wild-type GRHL3, and in zebrafish embryos, perturbed periderm development. We conclude that this mutation is an etiologic variant for nonsyndromic CP and is one of few functional variants identified to date for nonsyndromic orofacial clefting. This finding advances our understanding of the genetic basis of craniofacial development and might ultimately lead to improvements in recurrence risk prediction, treatment, and prognosis.

Association of low-frequency genetic variants in regulatory regions with nonsyndromic orofacial clefts.

Genome-wide scans have shown that common risk alleles for orofacial clefts (OFC) tend to be located in noncoding regulatory elements and cumulatively explain only part of the heritability of OFCs. Low-frequency variants may account for some of the "missing" heritability. Therefore, we scanned low-frequency variants located within putative craniofacial enhancers to identify novel OFC risk variants and implicate new regulatory elements in OFC pathogenesis. Analyses were performed in a multiethnic sample of 1,995 cases of cleft lip with or without cleft palate (CL/P), 221 cases with cleft palate (CP) only, and 1,576 unaffected controls. One hundred and nineteen putative craniofacial enhancers identified from ChIP-Seq studies in craniofacial tissues or cell lines contained multiple low-frequency (0.01-1%) variants, which we genotyped in participants using a custom Illumina panel. Two complementary statistical approaches, sequence kernel association test and combined multivariate and collapsing, were used to test association of the aggregated low-frequency variants across each enhancer region with CL/P and CP. We discovered a significant association between CP and a branchial arch enhancer near FOXP1 (mm60; p-value = .0002). Additionally, we observed a suggestive association between CL/P and a forebrain enhancer near FOXE1 (hs1717; p-value = .001). These findings suggest that low-frequency variants in craniofacial enhancer regions contribute to the complex etiology of nonsyndromic OFCs.

Identification of 16q21 as a modifier of nonsyndromic orofacial cleft phenotypes.

Orofacial clefts (OFCs) are common, complex birth defects with extremely heterogeneous phenotypic presentations. Two common subtypes-cleft lip alone (CL) and CL plus cleft palate (CLP)-are typically grouped into a single phenotype for genetic analysis (i.e., CL with or without cleft palate, CL/P). However, mounting evidence suggests there may be unique underlying pathophysiology and/or genetic modifiers influencing expression of these two phenotypes. To this end, we performed a genome-wide scan for genetic modifiers by directly comparing 450 CL cases with 1,692 CLP cases from 18 recruitment sites across 13 countries from North America, Central or South America, Asia, Europe, and Africa. We identified a region on 16q21 that is strongly associated with different cleft type (P = 5.611 × 10-8 ). We also identified significant evidence of gene-gene interactions between this modifier locus and two recognized CL/P risk loci: 8q21 and 9q22 (FOXE1) (P = 0.012 and 0.023, respectively). Single nucleotide polymorphism (SNPs) in the 16q21 modifier locus demonstrated significant association with CL over CLP. The marker alleles on 16q21 that increased risk for CL were found at highest frequencies among individuals with a family history of CL (P = 0.003). Our results demonstrate the existence of modifiers for which type of OFC develops and suggest plausible elements responsible for phenotypic heterogeneity, further elucidating the complex genetic architecture of OFCs.

A multi-ethnic genome-wide association study identifies novel loci for non-syndromic cleft lip with or without cleft palate on 2p24.2, 17q23 and 19q13.

Orofacial clefts (OFCs), which include non-syndromic cleft lip with or without cleft palate (CL/P), are among the most common birth defects in humans, affecting approximately 1 in 700 newborns. CL/P is phenotypically heterogeneous and has a complex etiology caused by genetic and environmental factors. Previous genome-wide association studies (GWASs) have identified at least 15 risk loci for CL/P. As these loci do not account for all of the genetic variance of CL/P, we hypothesized the existence of additional risk loci. We conducted a multiethnic GWAS in 6480 participants (823 unrelated cases, 1700 unrelated controls and 1319 case-parent trios) with European, Asian, African and Central and South American ancestry. Our GWAS revealed novel associations on 2p24 near FAM49A, a gene of unknown function (P = 4.22 × 10-8), and 19q13 near RHPN2, a gene involved in organizing the actin cytoskeleton (P = 4.17 × 10-8). Other regions reaching genome-wide significance were 1p36 (PAX7), 1p22 (ARHGAP29), 1q32 (IRF6), 8q24 and 17p13 (NTN1), all reported in previous GWASs. Stratification by ancestry group revealed a novel association with a region on 17q23 (P = 2.92 × 10-8) among individuals with European ancestry. This region included several promising candidates including TANC2, an oncogene required for development, and DCAF7, a scaffolding protein required for craniofacial development. In the Central and South American ancestry group, significant associations with loci previously identified in Asian or European ancestry groups reflected their admixed ancestry. In summary, we have identified novel CL/P risk loci and suggest new genes involved in craniofacial development, confirming the highly heterogeneous etiology of OFCs.

Genome-wide meta-analyses of nonsyndromic orofacial clefts identify novel associations between FOXE1 and all orofacial clefts, and TP63 and cleft lip with or without cleft palate.

Nonsyndromic orofacial clefts (OFCs) are a heterogeneous group of common craniofacial birth defects with complex etiologies that include genetic and environmental risk factors. OFCs are commonly categorized as cleft lip with or without cleft palate (CL/P) and cleft palate alone (CP), which have historically been analyzed as distinct entities. Genes for both CL/P and CP have been identified via multiple genome-wide linkage and association studies (GWAS); however, altogether, known variants account for a minority of the estimated heritability in risk to these craniofacial birth defects. We performed genome-wide meta-analyses of CL/P, CP, and all OFCs across two large, multiethnic studies. We then performed population-specific meta-analyses in sub-samples of Asian and European ancestry. In addition to observing associations with known variants, we identified a novel genome-wide significant association between SNPs located in an intronic TP63 enhancer and CL/P (p = 1.16 × 10-8). Several novel loci with compelling candidate genes approached genome-wide significance on 4q21.1 (SHROOM3), 12q13.13 (KRT18), and 8p21 (NRG1). In the analysis of all OFCs combined, SNPs near FOXE1 reached genome-wide significance (p = 1.33 × 10-9). Our results support the highly heterogeneous nature of OFCs and illustrate the utility of meta-analysis for discovering new genetic risk factors.

Association studies of low-frequency coding variants in nonsyndromic cleft lip with or without cleft palate.

Nonsyndromic cleft lip with or without cleft palate (NSCL/P) is a group of common human birth defects with complex etiology. Although genome-wide association studies have successfully identified a number of risk loci, these loci only account for about 20% of the heritability of orofacial clefts. The "missing" heritability may be found in rare variants, copy number variants, or interactions. In this study, we investigated the role of low-frequency variants genotyped in 1995 cases and 1626 controls on the Illumina HumanCore + Exome chip. We performed two statistical tests, Sequence Kernel Association Test (SKAT) and Combined Multivariate and Collapsing (CMC) method using two minor allele frequency cutoffs (1% and 5%). We found that a burden of low-frequency coding variants in N4BP2, CDSN, PRTG, and AHRR were associated with increased risk of NSCL/P. Low-frequency variants in other genes were associated with decreased risk of NSCL/P. These results demonstrate that low-frequency variants contribute to the genetic etiology of NSCL/P.

Maternal Mental Health in Latinx/Hispanic Women after COVID-19 Pandemic.

The recent COVID-19 global emergency may have ripple effects on mental health of many people worldwide. This is especially true for populations like birthing and postpartum women where many changes to daily routines, access to medical care, work-related routines and socialization were experienced. This brief report presents data from an ongoing cohort study aiming to describe maternal mental health during the pandemic T12 (March 2020 to April 2021) with post-pandemic T2 (May 2022 to May 2023) of mothers followed in Puerto Rico. 47 out of 100 mothers have been recalled and assessed with psychosocial interviews (COPE-IUS) and assessments of anxiety (GAD-7) and depression (PHQ-9). Paired t-test revealed mean scores of depressions (PHQ-9) were significantly higher for T2 with a mean of 6.35 and a range of 4.4+/- than for T1 where mean was 5.15 (+/- 2.9), t=-1.954, df=45, p < .05. Similarly, anxiety scores (GAD-7) were significantly higher in T2 6.67 (4.2) than for T1 5.35 (3.7), t=-1.8, df=45, p < .05. Also, COPE-IUS Post-pandemic psychosocial interview results evidence that 80% of mothers do not feel the COVID-19 pandemic is a significant stressor at T2 and are able to maintain routine activities with no social distancing measures. However, reports of loneliness, sadness, worry, and fear continue to be present. Our findings point to the need to further identify other contributing factors to the deterioration of maternal mental health during the perinatal/peripartum period (pregnancy, birth, and postpartum) in Puerto Rico. Possibly the effects of repeated adversity that has been present in the island (multiple environmental stressors, history of traumatic experiences, and constant hardships) may all have cumulative impact over maternal mental health during the perinatal/peripartum period.

Genome-wide study of gene-by-sex interactions identifies risks for cleft palate.

Structural birth defects affect 3-4% of all live births and, depending on the type, tend to manifest in a sex-biased manner. Orofacial clefts (OFCs) are the most common craniofacial structural birth defects and are often divided into cleft lip with or without cleft palate (CL/P) and cleft palate only (CP). Previous studies have found sex-specific risks for CL/P, but these risks have yet to be evaluated in CP. CL/P is more common in males and CP is more frequently observed in females, so we hypothesized there would also be sex-specific differences for CP. Using a trio-based cohort, we performed sex-stratified genome-wide association studies (GWAS) based on proband sex followed by a genome-wide gene-by-sex (GxS) interaction testing. There were 13 loci significant for GxS interactions, with the top finding in LTBP1 (RR=3.37 [2.04 - 5.56], p=1.93x10 -6 ). LTBP1 plays a role in regulating TGF-B bioavailability, and knockdown in both mice and zebrafish lead to craniofacial anomalies. Further, there is evidence for differential expression of LTBP1 between males and females in both mice and humans. Therefore, we tested the association between the imputed genetically regulated gene expression of genes with significant GxS interactions and the CP phenotype. We found significant association for LTBP1 in cell cultured fibroblasts in female probands (p=0.0013) but not in males. Taken altogether, we show there are sex-specific risks for CP that are otherwise undetectable in a combined sex cohort, and LTBP1 is a candidate risk gene, particularly in females.

Pre- and peri-natal hurricane exposure alters DNA methylation patterns in children.

Hurricane Maria was the worst recorded natural disaster to affect Puerto Rico. Increased stress in pregnant women during and in the aftermath of the hurricane may have induced epigenetic changes in their infants, which could affect gene expression. Stage of gestation at the time of the event was associated with significant differences in DNA methylation in the infants, especially those who were at around 20-25 weeks of gestation when the hurricane struck. Significant differences in DNA methylation were also associated with maternal mental status assessed after the hurricane, and with property damage. Hurricane Maria could have long lasting consequences to children who were exposed to this disaster during pregnancy.

Trio-based GWAS identifies novel associations and subtype-specific risk factors for cleft palate.

Cleft palate (CP) is one of the most common craniofacial birth defects; however, there are relatively few established genetic risk factors associated with its occurrence despite high heritability. Historically, CP has been studied as a single phenotype, although it manifests across a spectrum of defects involving the hard and/or soft palate. We performed a genome-wide association study using transmission disequilibrium tests of 435 case-parent trios to evaluate broad risks for any cleft palate (ACP) (n = 435), and subtype-specific risks for any cleft soft palate (CSP), (n = 259) and any cleft hard palate (CHP) (n = 125). We identified a single genome-wide significant locus at 9q33.3 (lead SNP rs7035976, p = 4.24 × 10-8) associated with CHP. One gene at this locus, angiopoietin-like 2 (ANGPTL2), plays a role in osteoblast differentiation. It is expressed both in craniofacial tissue of human embryos and developing mouse palatal shelves. We found 19 additional loci reaching suggestive significance (p < 5 × 10-6), of which only one overlapped between groups (chromosome 17q24.2, ACP and CSP). Odds ratios for the 20 loci were most similar across all 3 groups for SNPs associated with the ACP group, but more distinct when comparing SNPs associated with either subtype. We also found nominal evidence of replication (p < 0.05) for 22 SNPs previously associated with orofacial clefts. Our study to evaluate CP risks in the context of its subtypes and we provide newly reported associations affecting the broad risk for CP as well as evidence of subtype-specific risks.

Trio-based GWAS identifies novel associations and subtype-specific risk factors for cleft palate.

Orofacial clefts (OFCs) are the most common craniofacial birth defects and are often categorized into two etiologically distinct groups: cleft lip with or without cleft palate (CL/P) and isolated cleft palate (CP). CP is highly heritable, but there are still relatively few established genetic risk factors associated with its occurrence compared to CL/P. Historically, CP has been studied as a single phenotype despite manifesting across a spectrum of defects involving the hard and/or soft palate. We performed GWAS using transmission disequilibrium tests using 435 case-parent trios to evaluate broad risks for any cleft palate (ACP, n=435), as well as subtype-specific risks for any cleft soft palate (CSP, n=259) and any cleft hard palate (CHP, n=125). We identified a single genome-wide significant locus at 9q33.3 (lead SNP rs7035976, p=4.24×10 -8 ) associated with CHP. One gene at this locus, angiopoietin-like 2 ( ANGPTL2 ), plays a role in osteoblast differentiation. It is expressed in craniofacial tissue of human embryos, as well as in the developing mouse palatal shelves. We found 19 additional loci reaching suggestive significance (p<5×10 -6 ), of which only one overlapped between groups (chromosome 17q24.2, ACP and CSP). Odds ratios (ORs) for each of the 20 loci were most similar across all three groups for SNPs associated with the ACP group, but more distinct when comparing SNPs associated with either the CSP or CHP groups. We also found nominal evidence of replication (p<0.05) for 22 SNPs previously associated with cleft palate (including CL/P). Interestingly, most SNPs associated with CL/P cases were found to convey the opposite effect in those replicated in our dataset for CP only. Ours is the first study to evaluate CP risks in the context of its subtypes and we provide newly reported associations affecting the broad risk for CP as well as evidence of subtype-specific risks.