skip to main content


The NSF Public Access Repository (NSF-PAR) system and access will be unavailable from 11:00 PM ET on Thursday, June 13 until 2:00 AM ET on Friday, June 14 due to maintenance. We apologize for the inconvenience.

Title: Whole-genome sequencing reveals de-novo mutations associated with nonsyndromic cleft lip/palate

The majority (85%) of nonsyndromic cleft lip with or without cleft palate (nsCL/P) cases occur sporadically, suggesting a role for de novo mutations (DNMs) in the etiology of nsCL/P. To identify high impact protein-altering DNMs that contribute to the risk of nsCL/P, we conducted whole-genome sequencing (WGS) analyses in 130 African case-parent trios (affected probands and unaffected parents). We identified 162 high confidence protein-altering DNMs some of which are based on available evidence, contribute to the risk of nsCL/P. These include novel protein-truncating DNMs in theACTL6A, ARHGAP10, MINK1, TMEM5andTTNgenes; as well as missense variants inACAN, DHRS3, DLX6, EPHB2, FKBP10, KMT2D, RECQL4, SEMA3C, SEMA4D, SHH, TP63,andTULP4. Many of these protein-altering DNMs were predicted to be pathogenic. Analysis using mouse transcriptomics data showed that some of these genes are expressed during the development of primary and secondary palate. Gene-set enrichment analysis of the protein-altering DNMs identified palatal development and neural crest migration among the few processes that were significantly enriched. These processes are directly involved in the etiopathogenesis of clefting. The analysis of the coding sequence in the WGS data provides more evidence of the opportunity for novel findings in the African genome.

more » « less
Award ID(s):
Author(s) / Creator(s):
; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; more » ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; « less
Publisher / Repository:
Nature Publishing Group
Date Published:
Journal Name:
Scientific Reports
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. 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 nearFOXP1(mm60;p‐value = .0002). Additionally, we observed a suggestive association between CL/P and a forebrain enhancer nearFOXE1(hs1717;p‐value = .001). These findings suggest that low‐frequency variants in craniofacial enhancer regions contribute to the complex etiology of nonsyndromic OFCs.

    more » « less
  2. Summary

    The chromodomain family member chromodomain 1 (CHD1) has been shown to have numerous critical molecular functions including transcriptional regulation, splicing, and DNA repair. Complete loss of function of this gene is not compatible with life. On the other hand, missense and copy number variants ofCHD1can result in intellectual disabilities and craniofacial malformations in human patients including cleft palate and Pilarowski–Bjornsson Syndrome. We have used the aquatic developmental model organismXenopus laevis, to determine a specific role for Chd1 in such cranioafcial disorders. Protein and gene knockdown techniques inXenopus, including antisense oligos and mosaic Crispr/Cas9‐mediated mutagenesis, recapitulated the craniofacial defects observed in humans. Further analysis indicated that embryos deficient in Chd1 had defects in cranial neural crest development and jaw cartilage morphology. Additionally, flow cytometry and immunohistochemistry revealed that decreased Chd1 resulted in increased in apoptosis in the developing head. Together, these experiments demonstrate that Chd1 is critical for fundamental processes and cell survival in craniofacial development. We also presented evidence that Chd1 is regulated by retinoic acid signaling during craniofacial development. Expression levels ofchd1mRNA, specifically in the head, were increased by RAR agonist exposure and decreased upon antagonist treatment. Subphenotypic levels of an RAR antagonist and Chd1 morpholinos synergized to result in orofacial defects. Further, RAR DNA binding sequences (RAREs) were detected inchd1regulatory regions by bioinformatic analysis. In summary, by combining human genetics and experiments in an aquatic model we now have a better understanding of the role of CHD1 in craniofacial disorders.

    more » « less
  3. Abstract

    Genome‐wide association studies (GWAS) of alcohol dependence (AD) have reliably identified variation within alcohol metabolizing genes (eg,ADH1B) but have inconsistently located other signals, which may be partially attributable to symptom heterogeneity underlying the disorder. We conducted GWAS of DSM‐IV AD (primary analysis), DSM‐IV AD criterion count (secondary analysis), and individual dependence criteria (tertiary analysis) among 7418 (1121 families) European American (EA) individuals from the Collaborative Study on the Genetics of Alcoholism (COGA). Trans‐ancestral meta‐analyses combined these results with data from 3175 (585 families) African‐American (AA) individuals from COGA. In the EA GWAS, three loci were genome‐wide significant: rs1229984 inADH1Bfor AD criterion count (P= 4.16E−11) andDesire to cut drinking(P= 1.21E−11); rs188227250 (chromosome 8,Drinking more than intended,P= 6.72E−09); rs1912461 (chromosome 15,Time spent drinking,P= 1.77E−08). In the trans‐ancestral meta‐analysis, rs1229984 was associated with multiple phenotypes and two additional loci were genome‐wide significant: rs61826952 (chromosome 1, DSM‐IV AD,P= 8.42E−11); rs7597960 (chromosome 2,Time spent drinking,P= 1.22E−08). Associations with rs1229984 and rs18822750 were replicated in independent datasets. Polygenic risk scores derived from the EA GWAS of AD predicted AD in two EA datasets (P < .01; 0.61%‐1.82% of variance). Identified novel variants (ie, rs1912461, rs61826952) were associated with differential central evoked theta power (loss − gain;P= .0037) and reward‐related ventral striatum reactivity (P= .008), respectively. This study suggests that studying individual criteria may unveil new insights into the genetic etiology of AD liability.

    more » « less
  4. Introduction:Heavy metal pollutants can have long lasting negative impacts on ecosystem health and can shape the evolution of species. The persistent and ubiquitous nature of heavy metal pollution provides an opportunity to characterize the genetic mechanisms that contribute to metal resistance in natural populations.

    Methods:We examined variation in resistance to copper, a common heavy metal contaminant, using wild collections of the model organismDrosophila melanogaster. Flies were collected from multiple sites that varied in copper contamination risk. We characterized phenotypic variation in copper resistance within and among populations using bulked segregant analysis to identify regions of the genome that contribute to copper resistance.

    Results and Discussion:Copper resistance varied among wild populations with a clear correspondence between resistance level and historical exposure to copper. We identified 288 SNPs distributed across the genome associated with copper resistance. Many SNPs had population-specific effects, but some had consistent effects on copper resistance in all populations. Significant SNPs map to several novel candidate genes involved in refolding disrupted proteins, energy production, and mitochondrial function. We also identified one SNP with consistent effects on copper resistance in all populations nearCG11825, a gene involved in copper homeostasis and copper resistance. We compared the genetic signatures of copper resistance in the wild-derived populations to genetic control of copper resistance in theDrosophilaSynthetic Population Resource (DSPR) and theDrosophilaGenetic Reference Panel (DGRP), two copper-naïve laboratory populations. In addition toCG11825, which was identified as a candidate gene in the wild-derived populations and previously in the DSPR, there was modest overlap of copper-associated SNPs between the wild-derived populations and laboratory populations. Thirty-one SNPs associated with copper resistance in wild-derived populations fell within regions of the genome that were associated with copper resistance in the DSPR in a prior study. Collectively, our results demonstrate that the genetic control of copper resistance is highly polygenic, and that several loci can be clearly linked to genes involved in heavy metal toxicity response. The mixture of parallel and population-specific SNPs points to a complex interplay between genetic background and the selection regime that modifies the effects of genetic variation on copper resistance.

    more » « less
  5. Abstract Background

    Cleft palate is one of the most prevalent birth defects. Mice are useful for studying palate development because of their morphological and genetic similarities to humans. In mice, palate development occurs between embryonic days (E)11.5 to 15.5. Single cell transcriptional profiles of palate cell populations have been a valuable resource for the craniofacial research community, but we lack a single cell transcriptional profile for anterior palate at E13.5, at the transition from proliferation to shelf elevation.


    A detailed single cell RNA sequencing analysis reveals heterogeneity in expression profiles of the cell populations of the E13.5 anterior palate. Hybridization chain reaction RNA fluorescent in situ hybridization (HCR RNA FISH) reveals epithelial populations segregate into layers. Mesenchymal populations spatially segregate into four domains. One of these mesenchymal populations expresses ligands and receptors distinct from the rest of the mesenchyme, suggesting that these cells have a unique function. RNA velocity analysis shows two terminal cell states that contribute to either the proximal or distal palatal regions emerge from a single progenitor pool.


    This single cell resolution expression data and detailed analysis from E13.5 anterior palate provides a powerful resource for mechanistic insight into secondary palate morphogenesis for the craniofacial research community.

    more » « less