Search

Rare coding variants are a current focus in studies of complex disease. Previously, at least 68 rare coding variants were reported from candidate gene sequencing studies in non-syndromic cleft lip and palate ( NSCL/P), a common birth defect. Advances in sequencing technology have now resulted in thousands of sequenced exomes, providing a large resource for comparative genetic studies. We collated rare coding variants reported to contribute to NSCL/P and compared them to variants identified from control exome databases to determine if some might be rare but benign variants. Seventy-one percentage of the variants described as etiologic for NSCL/P were not present in the exome data, suggesting that many likely contribute to disease. Our results strongly support a role for rare variants previously reported in the majority of NSCL/P candidate genes but diminish support for variants in others. However, because clefting is a complex trait it is not possible to be definitive about the role of any particular variant for its risk for NSCL/P. [ABSTRACT FROM AUTHOR]

The go-to cognitive enhancers of today are those that are widely available rather than optimal for the user, including drugs typically prescribed for treatment of ADHD and sleep disturbances such as narcolepsy. While highly effective in their intended therapeutic role, performance gains in healthy populations are modest at best and profoundly inconsistent across subgroups and individuals. We propose a method for in silico screening of possible novel cognitive enhancers followed by high-throughput in vivo and in vitro validation. The proposed method uses gene expression data to evaluate the the collection of activated or suppressed signaling pathways in tissues or neurons of the cognitively enhanced brain. An algorithm maps expression data onto signaling pathways and quantifies their individual activation strength. The collective pathways and their activation form, a biological fingerprint of cognitive enhancement (or any other condition of interest). Drugs can then be screened and ranked based on their ability to minimize, mimic, or exaggerate pathway activation or suppression within that cloud. Using this approach, one may predict the efficacy of many drugs that may enhance various aspects of cognition.

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 (G × S) interaction testing. There were 13 loci significant for G × S interactions, with the top finding in LTBP1 (RR = 3.37 [2.04-5.56], p = 1.93 × 10 -6 ). LTBP1 plays a role in regulating TGF-β 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 G × S 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., (© 2024. The Author(s), under exclusive licence to Springer-Verlag GmbH Germany, part of Springer Nature.)

Congenital heart defects (CHDs) are the most common structural birth defect and are present in 40-50% of children born with Down syndrome (DS). To characterize the genetic architecture of DS-associated CHD, we sequenced genomes of a multiethnic group of children with DS and a CHD (n=886: atrioventricular septal defects (AVSD), n=438; atrial septal defects (ASD), n=122; ventricular septal defects (VSD), n=170; other types of CHD, n=156) and DS with a structurally normal heart (DS+NH, n=572). We performed four GWAS for common variants (MAF>0.05) comparing DS with CHD, stratified by CHD-subtype, to DS+NH controls. Although no SNP achieved genome-wide significance, multiple loci in each analysis achieved suggestive significance (p<2×10 -6 ). Of these, the 1p35.1 locus (near RBBP4 ) was specifically associated with ASD risk and the 5q35.2 locus (near MSX2 ) was associated with any type of CHD. Each of the suggestive loci contained one or more plausible candidate genes expressed in the developing heart. While no SNP replicated (p<2×10 -6 ) in an independent cohort of DS+CHD (DS+CHD: n=229; DS+NH: n=197), most SNPs that were suggestive in our GWASs remained suggestive when meta-analyzed with the GWASs from the replication cohort. These results build on previous work to identify genetic modifiers of DS-associated CHD.

Orofacial cleft (OFC) is a common human congenital anomaly. Epithelial-specific RNA splicing regulators ESRP1 and ESRP2 regulate craniofacial morphogenesis and their disruption result in OFC in zebrafish, mouse and humans. Using esrp1/2 mutant zebrafish and murine Py2T cell line models, we functionally tested the pathogenicity of human ESRP1/2 gene variants. We found that many variants predicted by in silico methods to be pathogenic were functionally benign. Esrp1 also regulates the alternative splicing of Ctnnd1 and these genes are co-expressed in the embryonic and oral epithelium. In fact, over-expression of ctnnd1 is sufficient to rescue morphogenesis of epithelial-derived structures in esrp1/2 zebrafish mutants. Additionally, we identified 13 CTNND1 variants from genome sequencing of OFC cohorts, confirming CTNND1 as a key gene in human OFC. This work highlights the importance of functional assessment of human gene variants and demonstrates the critical requirement of Esrp-Ctnnd1 acting in the embryonic epithelium to regulate palatogenesis., (© 2024. The Author(s).)

Orofacial clefts (OFCs) are common congenital birth defects with various etiologies, including genetic variants. Online Mendelian Inheritance in Man (OMIM) annotated several hundred genes involving OFCs. Furthermore, several hundreds of de novo variants (DNVs) have been identified from individuals with OFCs. Some DNVs are related to known OFC genes or pathways, but there are still many DNVs whose relevance to OFC development is unknown. To explore novel gene functions and their cellular expression profiles, we focused on DNVs in genes that were not listed in OMIM. We collected 960 DNVs in 853 genes from published studies and curated these genes, based on the DNVs' deleteriousness, into 230 and 23 genes related to cleft lip with or without cleft palate (CL/P) and cleft palate only (CPO), respectively. For comparison, we curated 178 CL/P and 277 CPO genes from OMIM. In CL/P, the pathways enriched in DNV and OMIM genes were significantly overlapped (p = 0.002). Single-cell RNA sequencing (scRNA-seq) analysis of mouse lip development revealed that both gene sets had abundant expression in the ectoderm (DNV genes: adjusted p = 0.032, OMIM genes: adjusted p < 0.0002), while only DNV genes were enriched in the endothelium (adjusted p = 0.032). Although we did not achieve significant findings using CPO gene sets, which was mainly due to the limited number of DNV genes, scRNA-seq analysis implicated various expression patterns among DNV and OMIM genes. Our results suggest that combinatory pathway and scRNA-seq data analyses are helpful for contextualizing genes in OFC development., Competing Interests: Declaration of interests The authors declare no competing interests., (Copyright © 2024 The Author(s). Published by Elsevier Inc. All rights reserved.)

Orofacial cleft (OFC) is a common human congenital anomaly. Epithelial-specific RNA splicing regulators ESRP1 and ESRP2 regulate craniofacial morphogenesis and their disruption result in OFC in zebrafish, mouse and humans. Using esrp1/2 mutant zebrafish and murine Py2T cell line models, we functionally tested the pathogenicity of human ESRP1/2 gene variants. We found that many variants predicted by in silico methods to be pathogenic were functionally benign. Esrp1 also regulates the alternative splicing of Ctnnd1 and these genes are co-expressed in the embryonic and oral epithelium. In fact, over-expression of ctnnd1 is sufficient to rescue morphogenesis of epithelial-derived structures in esrp1/2 zebrafish mutants. Additionally, we identified 13 CTNND1 variants from genome sequencing of OFC cohorts, confirming CTNND1 as a key gene in human OFC. This work highlights the importance of functional assessment of human gene variants and demonstrates the critical requirement of Esrp-Ctnnd1 acting in the embryonic epithelium to regulate palatogenesis.

Orofacial clefts (OFCs) are common, affecting 1:1000 live births. OFCs occur across a phenotypic spectrum - including cleft lip (CL), cleft lip and palate (CLP), or cleft palate (CP) - and can be further subdivided based on laterality, severity, or specific structures affected. Herein we review what is known about the genetic architecture underlying each of these subtypes, considering both shared and subtype-specific risks. While there are more known genetic similarities between CL and CLP than CP, recent research supports both shared and subtype-specific genetic risk factors within and between phenotypic classifications of OFCs. Larger sample sizes and deeper phenotyping data will be of increasing importance for the discovery of novel genetic risk factors for OFCs and various subtypes going forward., Competing Interests: Declaration of interests The authors have no conflicts of interest to disclose., (Copyright © 2024 Elsevier Ltd. All rights reserved.)

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., (© 2023. The Author(s), under exclusive licence to Springer-Verlag GmbH Germany, part of Springer Nature.)

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., (© 2023 Wiley Periodicals LLC.)

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., Competing Interests: Conflict of Interest The authors declare no conflicts of interests., (Copyright © 2023 The Authors. Published by Elsevier Inc. All rights reserved.)

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., Competing Interests: M.P.E. is a member of the Human Genetics and Genomics Advances Editorial Board., (© 2023 The Authors.)

Craniofacial disorders arise in early pregnancy and are one of the most common congenital defects. To fully understand how craniofacial disorders arise, it is essential to characterize gene expression during the patterning of the craniofacial region. To address this, we performed bulk and single-cell RNA-seq on human craniofacial tissue from 4-8 weeks post conception. Comparisons to dozens of other human tissues revealed 239 genes most strongly expressed during craniofacial development. Craniofacial-biased developmental enhancers were enriched +/- 400 kb surrounding these craniofacial-biased genes. Gene co-expression analysis revealed that regulatory hubs are enriched for known disease causing genes and are resistant to mutation in the normal healthy population. Combining transcriptomic and epigenomic data we identified 539 genes likely to contribute to craniofacial disorders. While most have not been previously implicated in craniofacial disorders, we demonstrate this set of genes has increased levels of de novo mutations in orofacial clefting patients warranting further study., (© 2023. The Author(s).)

Motivation: There is widespread interest in identifying genetic variants that exhibit parent-of-origin effects (POEs) wherein the effect of an allele on phenotype expression depends on its parental origin. POEs can arise from different phenomena including genomic imprinting and have been documented for many complex traits. Traditional tests for POEs require family data to determine parental origins of transmitted alleles. As most genome-wide association studies (GWAS) sample unrelated individuals (where allelic parental origin is unknown), the study of POEs in such datasets requires sophisticated statistical methods that exploit genetic patterns we anticipate observing when POEs exist. We propose a method to improve discovery of POE variants in large-scale GWAS samples that leverages potential pleiotropy among multiple correlated traits often collected in such studies. Our method compares the phenotypic covariance matrix of heterozygotes to homozygotes based on a Robust Omnibus Test. We refer to our method as the Parent of Origin Inference using Robust Omnibus Test (POIROT) of multiple quantitative traits., Results: Through simulation studies, we compared POIROT to a competing univariate variance-based method which considers separate analysis of each phenotype. We observed POIROT to be well-calibrated with improved power to detect POEs compared to univariate methods. POIROT is robust to non-normality of phenotypes and can adjust for population stratification and other confounders. Finally, we applied POIROT to GWAS data from the UK Biobank using BMI and two cholesterol phenotypes. We identified 338 genome-wide significant loci for follow-up investigation., Availability and Implementation: The code for this method is available at https://github.com/staylorhead/POIROT-POE., (© The Author(s) 2023. Published by Oxford University Press.)

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.

Mediation models are a set of statistical techniques that investigate the mechanisms that produce an observed relationship between an exposure variable and an outcome variable in order to deduce the extent to which the relationship is influenced by intermediate mediator variables. For a case-control study, the most common mediation analysis strategy employs a counterfactual framework that permits estimation of indirect and direct effects on the odds ratio scale for dichotomous outcomes, assuming either binary or continuous mediators. While this framework has become an important tool for mediation analysis, we demonstrate that we can embed this approach in a unified likelihood framework for mediation analysis in case-control studies that leverages more features of the data (in particular, the relationship between exposure and mediator) to improve efficiency of indirect effect estimates. One important feature of our likelihood approach is that it naturally incorporates cases within the exposure-mediator model to improve efficiency. Our approach does not require knowledge of disease prevalence and can model confounders and exposure-mediator interactions, and is straightforward to implement in standard statistical software. We illustrate our approach using both simulated data and real data from a case-control genetic study of lung cancer., (© 2022 John Wiley & Sons Ltd.)

Introduction: Orofacial clefts (OFCs) are common craniofacial birth defects with heterogeneous phenotype and etiology. Geneticists have applied nearly every available method and technology for further understanding of the genetic architectures of OFCs., Objective: This review describes the evidence for a genetic etiology in OFCs, statistical genetic approaches employed to identify genetic causes, and how the results have shaped our current understanding of the genetic architectures of syndromic and nonsyndromic OFCs., Conclusion: There has been rapid progress toward elucidating the genetic architectures of OFCs due to the availability of large collections of DNA samples from cases, controls, and families with OFCs and the consistent adoption of new methodologies and novel statistical approaches as they are developed. Genetic studies have identified rare and common variants influencing risk of OFCs in both Mendelian and complex forms of OFCs, blurring the distinction traditional categories used in genetic studies and clinical medicine., (© 2021 Wiley Periodicals LLC.)

Objective: This study assessed the feasibility of unpaid social media advertising to recruit participants affected with an orofacial cleft (OFC) for a genetic study., Design: This is a retrospective analysis of recruitment based on enrollment and participation in a genetic study. Participants completed a series of enrollment surveys, provided saliva samples, and completed postparticipation feedback surveys., Participants: Participants were eligible if they or a minor in their care were affected by an OFC, the affected participant was not adopted, and the mother of the affected individual had not taken antiseizure medication during pregnancy., Main Outcome Measures: Success of recruitment was evaluated from the number of enrolled participants and sample return rate., Results: In the first 12 months of recruitment, 313 individuals completed initial screening surveys; of these, 306 participants were eligible. A total of 263 individuals completed all online surveys and were sent DNA sample kits. One hundred sixty-two subject DNA samples were returned within 12 months of sending, for a return rate of 62%. Approximately two-thirds (66.3%) of all returned samples were sent back within the first 6 weeks after receiving DNA kits., Conclusions: Unpaid social media advertising enabled the recruitment of a large cohort of participants in a short time (12 months). The resulting study population was limited in racial and ethnic diversity, suggesting that other recruitment strategies will be needed for studies seeking specific demographic or socioeconomic groups. Nonetheless, social media recruitment was efficient and effective for recruiting participants for a genetic study in comparison to traditional clinic-based modes of recruitment.

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., (© 2022 Wiley Periodicals LLC.)

Non-syndromic cleft lip with or without cleft palate (NSCL/P) is a common birth defect, affecting approximately 1 in 700 births. NSCL/P has complex etiology including several known genes and environmental factors; however, known genetic risk variants only account for a small fraction of the heritability of NSCL/P. It is commonly suggested that gene-by-environment (G×E) interactions may help explain some of the "missing" heritability of NSCL/P. We conducted a genome-wide G×E interaction study in cases and controls of European ancestry with three common maternal exposures during pregnancy: alcohol, smoking, and vitamin use using a two-stage design. After selecting 127 loci with suggestive 2df tests for gene and G x E effects, 40 loci showed significant G x E effects after correcting for multiple tests. Notable interactions included SNPs of 6q22 near VGLL2 with alcohol and 6p22.3 near PRL with smoking. These interactions could provide new insights into the etiology of CL/P and new opportunities to modify risk through behavioral changes., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest. The reviewer ZJ declared a past collaboration with several of the authors EL and MM to the handling editor., (Copyright © 2022 Carlson, Shaffer, Deleyiannis, Hecht, Wehby, Christensen, Feingold, Weinberg, Marazita and Leslie.)

Background: In a large pedigree with an unusual phenotype of spastic paraplegia or dystonia and autosomal dominant inheritance, linkage analysis previously mapped the disease to chromosome 2q24-2q31., Objective: The aim of this study is to identify the genetic cause and molecular basis of an unusual autosomal dominant spastic paraplegia and dystonia., Methods: Whole exome sequencing following linkage analysis was used to identify the genetic cause in a large family. Cosegregation analysis was also performed. An additional 384 individuals with spastic paraplegia or dystonia were screened for pathogenic sequence variants in the adenosine triphosphate (ATP) synthase membrane subunit C locus 3 gene (ATP5MC3). The identified variant was submitted to the "GeneMatcher" program for recruitment of additional subjects. Mitochondrial functions were analyzed in patient-derived fibroblast cell lines. Transgenic Drosophila carrying mutants were studied for movement behavior and mitochondrial function., Results: Exome analysis revealed a variant (c.318C > G; p.Asn106Lys) (NM&#95;001689.4) in ATP5MC3 in a large family with autosomal dominant spastic paraplegia and dystonia that cosegregated with affected individuals. No variants were identified in an additional 384 individuals with spastic paraplegia or dystonia. GeneMatcher identified an individual with the same genetic change, acquired de novo, who manifested upper-limb dystonia. Patient fibroblast studies showed impaired complex V activity, ATP generation, and oxygen consumption. Drosophila carrying orthologous mutations also exhibited impaired mitochondrial function and displayed reduced mobility., Conclusion: A unique form of familial spastic paraplegia and dystonia is associated with a heterozygous ATP5MC3 variant that also reduces mitochondrial complex V activity., (© 2021 International Parkinson and Movement Disorder Society.)

Orofacial clefts (OFCs) are common (1 in 700 births) congenital malformations that include a cleft lip (CL) and cleft lip and palate (CLP). These OFC subtypes are also heterogeneous themselves, with the CL occurring on the left, right, or both sides of the upper lip. Unilateral CL and CLP have a 2:1 bias towards left-sided clefts, suggesting a nonrandom process. Here, we performed a study of left- and right-sided clefts within the CL and CLP subtypes to better understand the genetic factors controlling cleft laterality. We conducted genome-wide modifier analyses by comparing cases that had right unilateral CL (RCL; N = 130), left unilateral CL (LCL; N = 216), right unilateral CLP (RCLP; N = 416), or left unilateral CLP (LCLP; N = 638), and identified a candidate region on 4q28, 400 kb downstream from FAT4, that approached genome-wide significance for LCL versus RCL (p = 8.4 × 10 -8 ). Consistent with its potential involvement as a genetic modifier of CL, we found that Fat4 exhibits a specific domain of expression in the mesenchyme of the medial nasal processes that form the median upper lip. Overall, these results suggest that the epidemiological similarities in left- to right-sided clefts in CL and CLP are not reflected in the genetic association results., (© 2021 Wiley Periodicals LLC.)

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., Competing Interests: The authors have declared that no competing interests exist.

Non-syndromic cleft lip with or without cleft palate (nsCL/P) is a common congenital facial malformation with a multifactorial etiology. Genome-wide association studies (GWASs) have identified multiple genetic risk loci. However, functional interpretation of these loci is hampered by the underrepresentation in public resources of systematic functional maps representative of human embryonic facial development. To generate novel insights into the etiology of nsCL/P, we leveraged published GWAS data on nsCL/P as well as available chromatin modification and expression data on mid-facial development. Our analyses identified five novel risk loci, prioritized candidate target genes within associated regions, and highlighted distinct pathways. Furthermore, the results suggest the presence of distinct regulatory effects of nsCL/P risk variants throughout mid-facial development and shed light on its regulatory architecture. Our integrated data provide a platform to advance hypothesis-driven molecular investigations of nsCL/P and other human facial defects., Competing Interests: The authors declare no competing interests., (© 2021 The Authors.)

Purpose: To understand the consequences of the 3q29 deletion on medical, neurodevelopmental, psychiatric, brain structural, and neurological sequalae by systematic evaluation of affected individuals. To develop evidence-based recommendations using these data for effective clinical care., Methods: Thirty-two individuals with the 3q29 deletion were evaluated using a defined phenotyping protocol and standardized data collection instruments., Results: Medical manifestations were varied and reported across nearly every organ system. The most severe manifestations were congenital heart defects (25%) and the most common were gastrointestinal symptoms (81%). Physical examination revealed a high proportion of musculoskeletal findings (81%). Neurodevelopmental phenotypes represent a significant burden and include intellectual disability (34%), autism spectrum disorder (38%), executive function deficits (46%), and graphomotor weakness (78%). Psychiatric illness manifests across the lifespan with psychosis prodrome (15%), psychosis (20%), anxiety disorders (40%), and attention deficit-hyperactivity disorder (ADHD) (63%). Neuroimaging revealed structural anomalies of the posterior fossa, but on neurological exam study subjects displayed only mild or moderate motor vulnerabilities., Conclusion: By direct evaluation of 3q29 deletion study subjects, we document common features of the syndrome, including a high burden of neurodevelopmental and neuropsychiatric phenotypes. Evidence-based recommendations for evaluation, referral, and management are provided to help guide clinicians in the care of 3q29 deletion patients.

Two large studies of case-parent trios ascertained through a proband with a non-syndromic orofacial cleft (OFC, which includes cleft lip and palate, cleft lip alone, or cleft palate alone) were used to test for possible gene-environment (G × E) interaction between genome-wide markers (both observed and imputed) and self-reported maternal exposure to smoking, alcohol consumption, and multivitamin supplementation during pregnancy. The parent studies were as follows: GENEVA, which included 1,939 case-parent trios recruited largely through treatment centers in Europe, the United States, and Asia, and 1,443 case-parent trios from the Pittsburgh Orofacial Cleft Study (POFC) also ascertained through a proband with an OFC including three major racial/ethnic groups (European, Asian, and Latin American). Exposure rates to these environmental risk factors (maternal smoking, alcohol consumption, and multivitamin supplementation) varied across studies and among racial/ethnic groups, creating substantial differences in power to detect G × E interaction, but the trio design should minimize spurious results due to population stratification. The GENEVA and POFC studies were analyzed separately, and a meta-analysis was conducted across both studies to test for G × E interaction using the 2 df test of gene and G × E interaction and the 1 df test for G × E interaction alone. The 2 df test confirmed effects for several recognized risk genes, suggesting modest G × E effects. This analysis did reveal suggestive evidence for G × Vitamin interaction for CASP9 on 1p36 located about 3 Mb from PAX7 , a recognized risk gene. Several regions gave suggestive evidence of G × E interaction in the 1 df test. For example, for G × Smoking interaction, the 1 df test suggested markers in MUSK on 9q31.3 from meta-analysis. Markers near SLCO3A1 also showed suggestive evidence in the 1 df test for G × Alcohol interaction, and rs41117 near RETREG1 (a.k.a. FAM134B ) also gave suggestive significance in the meta-analysis of the 1 df test for G × Vitamin interaction. While it remains quite difficult to obtain definitive evidence for G × E interaction in genome-wide studies, perhaps due to small effect sizes of individual genes combined with low exposure rates, this analysis of two large case-parent trio studies argues for considering possible G × E interaction in any comprehensive study of complex and heterogeneous disorders such as OFC., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2021 Zhang, Venkataraghavan, Hetmanski, Leslie, Marazita, Feingold, Weinberg, Ruczinski, Taub, Scott, Ray and Beaty.)

Orofacial clefts (OFCs) are among the most prevalent craniofacial birth defects worldwide and create a significant public health burden. The majority of OFCs are non-syndromic and vary in prevalence by ethnicity. Africans have the lowest prevalence of OFCs (~ 1/2,500), Asians have the highest prevalence (~1/500), Europeans and Latin Americans lie somewhere in the middle (~1/800 and 1/900, respectively). Thus, ethnicity appears to be a major determinant of the risk of developing OFC. The Pittsburgh Orofacial Clefts Multiethnic study was designed to explore this ethnic variance, comprising a large number of families and individuals (~12,000 individuals) from multiple populations worldwide: US and Europe, Asians, mixed Native American/Caucasians, and Africans. In this current study, we analyzed 2,915 OFC cases, 6,044 unaffected individuals related to the OFC cases, and 2,685 controls with no personal or family history of OFC. Participants were grouped by their ancestry into African, Asian, European, and Central and South American subsets, and genome-wide association run on the combined sample as well as the four ancestry-based groups. We observed 22 associations to cleft lip with or without cleft palate at 18 distinct loci with p -values < 1e-06, including 10 with genome-wide significance (<5e-08), in the combined sample and within ancestry groups. Three loci - 2p12 (rs62164740, p = 6.27e-07), 10q22.2 (rs150952246, p = 3.14e-07), and 10q24.32 (rs118107597, p = 8.21e-07) are novel. Nine were in or near known OFC loci - PAX7, IRF6, FAM49A, DCAF4L2 , 8q24.21, NTN1, WNT3-WNT9B, TANC2 , and RHPN2 . The majority of the associations were observed only in the combined sample, European, and Central and South American groups. We investigated whether the observed differences in association strength were (a) purely due to sample sizes, (b) due to systematic allele frequency difference at the population level, or (c) due to the fact certain OFC-causing variants confer different amounts of risk depending on ancestral origin, by comparing effect sizes to observed allele frequencies of the effect allele in our ancestry-based groups. While some of the associations differ due to systematic differences in allele frequencies between groups, others show variation in effect size despite similar frequencies across ancestry groups., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2021 Mukhopadhyay, Feingold, Moreno-Uribe, Wehby, Valencia-Ramirez, Muñeton, Padilla, Deleyiannis, Christensen, Poletta, Orioli, Hecht, Buxó, Butali, Adeyemo, Vieira, Shaffer, Murray, Weinberg, Leslie and Marazita.)

Nonsyndromic orofacial clefts (OFCs) are a common birth defect and are phenotypically heterogenous in the structure affected by the cleft - cleft lip (CL) and cleft lip and palate (CLP) - as well as other features, such as the severity of the cleft. Here, we focus on bilateral and unilateral clefts as one dimension of OFC severity, because the genetic architecture of these subtypes is not well understood. We tested for subtype-specific genetic associations in 44 bilateral CL (BCL) cases, 434 unilateral CL (UCL) cases, 530 bilateral CLP cases (BCLP), 1123 unilateral CLP (UCLP) cases, and unrelated controls (N = 1626), using a mixed-model approach. While no novel loci were found, the genetic architecture of UCL was distinct compared to BCL, with 44.03% of suggestive loci having different effects between the two subtypes. To further understand the subtype-specific genetic risk factors, we performed a genome-wide scan for modifiers and found a significant modifier locus on 20p11 (p=7.53×10 -9 ), 300kb downstream of PAX1 , that associated with higher odds of BCL vs. UCL, and replicated in an independent cohort (p=0.0018) with no effect in BCLP (p>0.05). We further found that this locus was associated with normal human nasal shape. Taken together, these results suggest bilateral and unilateral clefts may have different genetic architectures. Moreover, our results suggest BCL, the rarest form of OFC, may be genetically distinct from the other OFC subtypes. This expands our understanding of modifiers for OFC subtypes and further elucidates the genetic mechanisms behind the phenotypic heterogeneity in OFCs., Competing Interests: Declaration of Interests: None of the authors report any conflicts of interest.

The FaceBase Consortium was established by the National Institute of Dental and Craniofacial Research in 2009 as a 'big data' resource for the craniofacial research community. Over the past decade, researchers have deposited hundreds of annotated and curated datasets on both normal and disordered craniofacial development in FaceBase, all freely available to the research community on the FaceBase Hub website. The Hub has developed numerous visualization and analysis tools designed to promote integration of multidisciplinary data while remaining dedicated to the FAIR principles of data management (findability, accessibility, interoperability and reusability) and providing a faceted search infrastructure for locating desired data efficiently. Summaries of the datasets generated by the FaceBase projects from 2014 to 2019 are provided here. FaceBase 3 now welcomes contributions of data on craniofacial and dental development in humans, model organisms and cell lines. Collectively, the FaceBase Consortium, along with other NIH-supported data resources, provide a continuously growing, dynamic and current resource for the scientific community while improving data reproducibility and fulfilling data sharing requirements., Competing Interests: Competing interestsThe authors declare no competing or financial interests., (© 2020. Published by The Company of Biologists Ltd.)

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., (Copyright © 2020 American Society of Human Genetics. Published by Elsevier Inc. All rights reserved.)

Understanding the etiology of congenital disorders requires interdisciplinary research and close collaborations between clinicians, geneticists and developmental biologists. The pace of gene discovery has quickened due to advances in sequencing technology, resulting in a wealth of publicly available sequence data but also a gap between gene discovery and crucial mechanistic insights provided by studies in model systems. In this Spotlight, I highlight the opportunities for developmental biologists to engage with human geneticists and genetic resources to advance the study of congenital disorders., Competing Interests: Competing interestsThe authors declare no competing or financial interests., (© 2020. Published by The Company of Biologists Ltd.)

Dermatoglyphic patterns on the fingers often differ in syndromes and other conditions with a developmental component, compared to the general population. Previous literature on the relationship between orofacial clefts-the most common craniofacial birth defect in humans-and dermatoglyphics is inconsistent, with some studies reporting altered pattern frequencies and/or increased asymmetry and others failing to find differences. To investigate dermatoglyphics in orofacial clefting, we obtained dermatoglyphic patterns in a large multiethnic cohort of orofacial cleft cases (N = 367), their unaffected family members (N = 836), and controls (N = 299). We categorized fingerprint pattern types from males and females who participated at five sites of the Pittsburgh Orofacial Cleft study (Hungary, United States of America (Pennsylvania, Texas), Spain, and Argentina). We also calculated a pattern dissimilarity score for each individual as a measure of left-right asymmetry. We tested for group differences in the number of arches, ulnar and radial loops, and whorls on each individual's hands, and in the pattern dissimilarity scores using ANOVA. After taking sex and site differences into account, we did not find any significant pattern count differences between cleft and non-cleft individuals. Notably, we did observe increased pattern dissimilarity in individuals with clefts, compared to both their unaffected relatives and controls. Increased dermatoglyphic pattern dissimilarity in individuals with nonsyndromic orofacial clefts may reflect a generalized developmental instability., Competing Interests: The authors have declared that no competing interests exist.

Orofacial clefts (OFCs) are among the most prevalent craniofacial birth defects worldwide and create a significant public health burden. The majority of OFCs are non-syndromic, and the genetic etiology of non-syndromic OFCs is only partially determined. Here, we analyze whole genome sequence (WGS) data for association with risk of OFCs in European and Colombian families selected from a multicenter family-based OFC study. This is the first large-scale WGS study of OFC in parent-offspring trios, and a part of the Gabriella Miller Kids First Pediatric Research Program created for the study of childhood cancers and structural birth defects. WGS provides deeper and more specific genetic data than using imputation on present-day single nucleotide polymorphic (SNP) marker panels. Genotypes of case-parent trios at single nucleotide variants (SNV) and short insertions and deletions (indels) spanning the entire genome were called from their sequences using human GRCh38 genome assembly, and analyzed for association using the transmission disequilibrium test. Among genome-wide significant associations, we identified a new locus on chromosome 21 in Colombian families, not previously observed in other larger OFC samples of Latin American ancestry. This locus is situated within a region known to be expressed during craniofacial development. Based on deeper investigation of this locus, we concluded that it contributed risk for OFCs exclusively in the Colombians. This study reinforces the ancestry differences seen in the genetic etiology of OFCs, and underscores the need for larger samples when studying for OFCs and other birth defects in populations with diverse ancestry.

Standard methods for case-control association studies of rare variation often treat disease outcome as a dichotomous phenotype. However, both theoretical and experimental studies have demonstrated that subjects with a family history of disease can be enriched for risk variation relative to subjects without such history. Assuming family history information is available, this observation motivates the idea of replacing the standard dichotomous outcome variable used in case-control studies with a more informative ordinal outcome variable that distinguishes controls (0), sporadic cases (1), and cases with a family history (2), with the expectation that we should observe increasing number of risk variants with increasing category of the ordinal variable. To leverage this expectation, we propose a novel rare-variant association test that incorporates family history information based on our previous GAMuT framework for rare-variant association testing of multivariate phenotypes. We use simulated data to show that, when family history information is available, our new method outperforms standard rare-variant association methods, like burden and SKAT tests, that ignore family history. We further illustrate our method using a rare-variant study of cleft lip and palate., (Copyright © 2020 by the Genetics Society of America.)

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., (© 2019 Wiley Periodicals, Inc.)

Mutations in IRF6, TFAP2A and GRHL3 cause orofacial clefting syndromes in humans. However, Tfap2a and Grhl3 are also required for neurulation in mice. Here, we found that homeostasis of Irf6 is also required for development of the neural tube and associated structures. Over-expression of Irf6 caused exencephaly, a rostral neural tube defect, through suppression of Tfap2a and Grhl3 expression. Conversely, loss of Irf6 function caused a curly tail and coincided with a reduction of Tfap2a and Grhl3 expression in tail tissues. To test whether Irf6 function in neurulation was conserved, we sequenced samples obtained from human cases of spina bifida and anencephaly. We found two likely disease-causing variants in two samples from patients with spina bifida. Overall, these data suggest that the Tfap2a-Irf6-Grhl3 genetic pathway is shared by two embryologically distinct morphogenetic events that previously were considered independent during mammalian development. In addition, these data suggest new candidates to delineate the genetic architecture of neural tube defects and new therapeutic targets to prevent this common birth defect., (Published by Oxford University Press 2019.)

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., (© 2018 Wiley Periodicals, Inc.)

Objective: We examined genomic variation potentially associated with the cortisol stress response in children having a painful medical procedure., Design: Children 4-10 years old having a peripheral intravenous line inserted provided saliva samples for evaluation of the cortisol response as a biological measure of distress: two on the day of the procedure and two at home on a nonstressful day for comparison values. Children and biological parents also provided samples for genotyping of variants with known or suspected association with the cortisol stress response. Analysis included child-only association and family-based transmission disequilibrium tests (TDTs)., Results: Genotype and phenotype data on the cortisol stress response were available from 326 children for child-only association analyses and 376 complete family trios for TDTs. Children were 50% female, an average of 7.5 years old, and mostly (83%) White/non-Hispanic. We identified four single-nucleotide polymorphisms (SNPs) potentially associated with the cortisol stress response: rs1176744 ( HTR3B), rs10062367 ( CRHBP), rs634479 ( OPRM1), and rs8030107 ( NTRK3). Family-based analysis identified a two-SNP haplotype in HTR1B suggestive for association with the cortisol response (rs6296, rs11568817). Allelic TDTs identified rs7897947 ( NFKB2) as potentially related to cortisol response., Conclusions: Findings provide preliminary evidence for genes potentially important in cortisol response to an acute stressor in children in the serotonin, dopamine, and brain-derived neurotrophic factor pathways, the hypothalamic-pituitary-adrenal axis, and the inflammatory response. Combined with analyses of related phenotypes and clinical data, these results could help identify patients at increased risk of adverse responses to painful medical procedures who might benefit from tailored interventions.

Motivation: De novo copy number deletions have been implicated in many diseases, but there is no formal method to date that identifies de novo deletions in parent-offspring trios from capture-based sequencing platforms., Results: We developed Minimum Distance for Targeted Sequencing (MDTS) to fill this void. MDTS has similar sensitivity (recall), but a much lower false positive rate compared to less specific CNV callers, resulting in a much higher positive predictive value (precision). MDTS also exhibited much better scalability., Availability and Implementation: MDTS is freely available as open source software from the Bioconductor repository., Supplementary Information: Supplementary data are available at Bioinformatics online., (© The Author(s) 2018. Published by Oxford University Press.)

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., (© 2018 Wiley Periodicals, Inc.)

Velopharyngeal dysfunction (VPD) occurs when the muscular soft palate (velum) and lateral pharyngeal walls are physically unable to separate the oral and nasal cavities during speech production leading to hypernasality and abnormal speech reduction. Because VPD is often associated with overt or submucous cleft palate, it could be present as a subclinical phenotype in families with a history of orofacial clefting. A key assumption to this model is that the overt and subclinical manifestations of the orofacial cleft phenotype exist on a continuum and therefore share common etiological factors. We performed a genome-wide association study in 976 unaffected relatives of isolated CP probands, 54 of whom had VPD. Five loci were significantly (p < 5 × 10 -8 ) associated with VPD: 3q29, 9p21.1, 12q21.31, 16p12.3 and 16p13.3. An additional 15 loci showing suggestive evidence of association with VPD were observed. Several genes known to be involved in orofacial clefting and craniofacial development are located in these regions, such as TFRC, PCYT1A, BNC2 and FREM1. Although further research is necessary, this could be an indication for a potential shared genetic architecture between VPD and cleft palate, and supporting the hypothesis that VPD is a subclinical phenotype of orofacial clefting.

We present association results from a large genome-wide association study of tooth agenesis (TA) as well as selective TA, including 1,944 subjects with congenitally missing teeth, excluding third molars, and 338,554 controls, all of European ancestry. We also tested the association of previously identified risk variants, for timing of tooth eruption and orofacial clefts, with TA. We report associations between TA and 9 novel risk variants. Five of these variants associate with selective TA, including a variant conferring risk of orofacial clefts. These results contribute to a deeper understanding of the genetic architecture of tooth development and disease. The few variants previously associated with TA were uncovered through candidate gene studies guided by mouse knockouts. Knowing the etiology and clinical features of TA is important for planning oral rehabilitation that often involves an interdisciplinary approach.

Genome-wide association scans of complex multipartite traits like the human face typically use preselected phenotypic measures. Here we report a data-driven approach to phenotyping facial shape at multiple levels of organization, allowing for an open-ended description of facial variation while preserving statistical power. In a sample of 2,329 persons of European ancestry, we identified 38 loci, 15 of which replicated in an independent European sample (n = 1,719). Four loci were completely new. For the others, additional support (n = 9) or pleiotropic effects (n = 2) were found in the literature, but the results reported here were further refined. All 15 replicated loci highlighted distinctive patterns of global-to-local genetic effects on facial shape and showed enrichment for active chromatin elements in human cranial neural crest cells, suggesting an early developmental origin of the facial variation captured. These results have implications for studies of facial genetics and other complex morphological traits.