DNA methylation effects on Van der Woude Syndrome phenotypic variability.

OBJECTIVE: Van der Woude Syndrome (VWS) classically presents with combinations of lip pits (LP) and orofacial clefts, with marked phenotypic discordance even amongst individuals carrying the same mutation. Such discordance suggests a possible role for epigenetic factors as phenotypic modifiers. Both IRF6 , causal for 70% of VWS cases, and TP63 interact in a regulatory loop to coordinate epithelial proliferation and differentiation for palatogenesis. We hypothesize that differential DNA methylation (DNAm) in CpG sites within regulatory regions of IRF6 and TP63 are associated with VWS phenotypic discordance. METHODS: We measured DNAm levels of CpG sites located in the promoter regions of IRF6 and TP63 and in an IRF6 enhancer element (MCS9.7) in 83 individuals with VWS grouped within 5 phenotypes for primary analysis: 1=CL+/-P+LP, 2=CL+/-P, 3=CP+LP, 4=CP, 5=LP and 2 phenotypes for secondary analysis: 1=any cleft and LP, 2= any cleft without LP. DNA samples were bisulfite converted and pyrosequenced with target-specific primers. Methylation levels were compared amongst phenotypes. RESULTS: CpG sites in the IRF6 promoter showed statistically significant differences in methylation among phenotypic groups in both analyses (P<0.05). Individuals with any form of cleft (Groups 1-4) had significantly higher methylation levels than individuals with lip pits only (Group 5). In the secondary analysis, individuals in Group 1 (cleft+LP) had significantly higher methylation than Group 2 (cleft only). CONCLUSION: Results indicated that hypermethylation of the IRF6 promoter is associated with more severe phenotypes (any cleft +/- lip pits); thus, possibly impacting an already genetically weakened IRF6 protein and leading to a more severe phenotype.

Damaging Mutations in AFDN Contribute to Risk of Nonsyndromic Cleft Lip With or Without Cleft Palate.

Novel or rare damaging mutations have been implicated in the developmental pathogenesis of nonsyndromic cleft lip with or without cleft palate (nsCL ± P). Thus, we investigated the human genome for high-impact mutations that could explain the risk of nsCL ± P in our cohorts.We conducted next-generation sequencing (NGS) analysis of 130 nsCL ± P case-parent African trios to identify pathogenic variants that contribute to the risk of clefting. We replicated this analysis using whole-exome sequence data from a Brazilian nsCL ± P cohort. Computational analyses were then used to predict the mechanism by which these variants could result in increased risks for nsCL ± P.We discovered damaging mutations within the AFDN gene, a cell adhesion molecule (CAMs) that was previously shown to contribute to cleft palate in mice. These mutations include p.Met1164Ile, p.Thr453Asn, p.Pro1638Ala, p.Arg669Gln, p.Ala1717Val, and p.Arg1596His. We also discovered a novel splicing p.Leu1588Leu mutation in this protein. Computational analysis suggests that these amino acid changes affect the interactions with other cleft-associated genes including nectins (PVRL1, PVRL2, PVRL3, and PVRL4) CDH1, CTNNA1, and CTNND1.This is the first report on the contribution of AFDN to the risk for nsCL ± P in humans. AFDN encodes AFADIN, an important CAM that forms calcium-independent complexes with nectins 1 and 4 (encoded by the genes PVRL1 and PVRL4). This discovery shows the power of NGS analysis of multiethnic cleft samples in combination with a computational approach in the understanding of the pathogenesis of nsCL ± P.

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 the ACTL6A, ARHGAP10, MINK1, TMEM5 and TTN genes; as well as missense variants in ACAN, DHRS3, DLX6, EPHB2, FKBP10, KMT2D, RECQL4, SEMA3C, SEMA4D, SHH, TP63, and TULP4. 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.

Variant analyses of candidate genes in orofacial clefts in multi-ethnic populations.

OBJECTIVES: Cleft lip with/without cleft palate and cleft palate only is congenital birth defects where the upper lip and/or palate fail to fuse properly during embryonic facial development. Affecting ~1.2/1000 live births worldwide, these orofacial clefts impose significant social and financial burdens on affected individuals and their families. Orofacial clefts have a complex etiology resulting from genetic variants combined with environmental covariates. Recent genome-wide association studies and whole-exome sequencing for orofacial clefts identified significant genetic associations and variants in several genes. Of these, we investigated the role of common/rare variants in SHH, RORA, MRPL53, ACVR1, and GDF11. MATERIALS AND METHODS: We sequenced these five genes in 1255 multi-ethnic cleft lip with/without palate and cleft palate only samples in order to find variants that may provide potential explanations for the missing heritability of orofacial clefts. Rare and novel variants were further analyzed using in silico predictive tools. RESULTS: Ninteen total variants of interest were found, with variant types including stop-gain, missense, synonymous, intronic, and splice-site variants. Of these, 3 novel missense variants were found, one in SHH, one in RORA, and one in GDF11. CONCLUSION: This study provides evidence that variants in SHH, RORA, MRPL53, ACVR1, and GDF11 may contribute to risk of orofacial clefts in various populations.

Genome-Wide Scan for Parent-of-Origin Effects in a sub-Saharan African Cohort With Nonsyndromic Cleft Lip and/or Cleft Palate (CL/P).

OBJECTIVE: Nonsyndromic cleft lip and/or cleft palate (NSCL/P) have multifactorial etiology where genetic factors, gene-environment interactions, stochastic factors, gene-gene interactions, and parent-of-origin effects (POEs) play cardinal roles. POEs arise when the parental origin of alleles differentially impacts the phenotype of the offspring. The aim of this study was to identify POEs that can increase risk for NSCL/P in humans using a genome-wide dataset. METHODS: The samples (174 case-parent trios from Ghana, Ethiopia, and Nigeria) included in this study were from the African only genome wide association studies (GWAS) that was published in 2019. Genotyping of individual DNA using over 2 million multiethnic and African ancestry-specific single-nucleotide polymorphisms from the Illumina Multi-Ethnic Genotyping Array v2 15070954 A2 (genome build GRCh37/hg19) was done at the Center for Inherited Diseases Research. After quality control checks, PLINK was employed to carry out POE analysis employing the pooled subphenotypes of NSCL/P. RESULTS: We observed possible hints of POEs at a cluster of genes at a 1 mega base pair window at the major histocompatibility complex class 1 locus on chromosome 6, as well as at other loci encompassing candidate genes such as ASB18, ANKEF1, AGAP1, GABRD, HHAT, CCT7, DNMT3A, EPHA7, FOXO3, lncRNAs, microRNA, antisense RNAs, ZNRD1, ZFAT, and ZBTB16. CONCLUSION: Findings from our study suggest that some loci may increase the risk for NSCL/P through POEs. Additional studies are required to confirm these suggestive loci in NSCL/P etiology.

Non-random distribution of deleterious mutations in the DNA and protein-binding domains of IRF6 are associated with Van Der Woude syndrome.

BACKGROUND: The development of the face occurs during the early days of intrauterine life by the formation of facial processes from the first Pharyngeal arch. Derangement in these well-organized fusion events results in Orofacial clefts (OFC). Van der Woude syndrome (VWS) is one of the most common causes of syndromic cleft lip and/or palate accounting for 2% of all cases. Mutations in the IRF6 gene account for 70% of cases with the majority of these mutations located in the DNA-binding (exon 3, 4) or protein-binding domains (exon 7-9). The current study was designed to update the list of IRF6 variants reported for VWS by compiling all the published mutations from 2013 to date as well as including the previously unreported VWS cases from Africa and Puerto Rico. METHODS: We used PubMed with the search terms; "Van der Woude syndrome," "Popliteal pterygium syndrome," "IRF6," and "Orofacial cleft" to identify eligible studies. We compiled the CADD score for all the mutations to determine the percentage of deleterious variants. RESULTS: Twenty-one new mutations were identified from nine papers. The majority of these mutations were in exon 4. Mutations in exon 3 and 4 had CADD scores between 20 and 30 and mutations in exon 7-9 had CADD scores between 30 and 40. The presence of higher CADD scores in the protein-binding domain (exon 7-9) further confirms the crucial role played by this domain in the function of IRF6. In the new cases, we identified five IRF6 mutations, three novel missense mutations (p.Phe36Tyr, p.Lys109Thr, and p.Gln438Leu), and two previously reported nonsense mutations (p.Ser424*and p.Arg250*). CONCLUSION: Mutations in the protein and DNA-binding domains of IRF6 ranked among the top 0.1% and 1% most deleterious genetic mutations, respectively. Overall, these findings expand the range of VWS mutations and are important for diagnostic and counseling purposes.

Identification of Isthmin 1 as a Novel Clefting and Craniofacial Patterning Gene in Humans.

Orofacial clefts are one of the most common birth defects, affecting 1-2 per 1000 births, and have a complex etiology. High-resolution array-based comparative genomic hybridization has increased the ability to detect copy number variants (CNVs) that can be causative for complex diseases such as cleft lip and/or palate. Utilizing this technique on 97 nonsyndromic cleft lip and palate cases and 43 cases with cleft palate only, we identified a heterozygous deletion of Isthmin 1 in one affected case, as well as a deletion in a second case that removes putative 3' regulatory information. Isthmin 1 is a strong candidate for clefting, as it is expressed in orofacial structures derived from the first branchial arch and is also in the same "synexpression group" as fibroblast growth factor 8 and sprouty RTK signaling antagonist 1a and 2, all of which have been associated with clefting. CNVs affecting Isthmin 1 are exceedingly rare in control populations, and Isthmin 1 scores as a likely haploinsufficiency locus. Confirming its role in craniofacial development, knockdown or clustered randomly interspaced short palindromic repeats/Cas9-generated mutation of isthmin 1 in Xenopus laevis resulted in mild to severe craniofacial dysmorphologies, with several individuals presenting with median clefts. Moreover, knockdown of isthmin 1 produced decreased expression of LIM homeobox 8, itself a gene associated with clefting, in regions of the face that pattern the maxilla. Our study demonstrates a successful pipeline from CNV identification of a candidate gene to functional validation in a vertebrate model system, and reveals Isthmin 1 as both a new human clefting locus as well as a key craniofacial patterning gene.

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.

Exome sequencing provides additional evidence for the involvement of ARHGAP29 in Mendelian orofacial clefting and extends the phenotypic spectrum to isolated cleft palate.

BACKGROUND: Recent advances in genomics methodologies, in particular the availability of next-generation sequencing approaches have made it possible to identify risk loci throughout the genome, in particular the exome. In the current study, we present findings from an exome study conducted in five affected individuals of a multiplex family with cleft palate only. METHODS: The GEnome MINIng (GEMINI) pipeline was used to functionally annotate the single nucleotide polymorphisms, insertions and deletions. Filtering methods were applied to identify variants that are clinically relevant and present in affected individuals at minor allele frequencies (≤1%) in the 1000 Genomes Project single nucleotide polymorphism database, Exome Aggregation Consortium, and Exome Variant Server databases. The bioinformatics tool Systems Tool for Craniofacial Expression-Based Gene Discovery was used to prioritize cleft candidates in our list of variants, and Sanger sequencing was used to validate the presence of identified variants in affected and unaffected relatives. RESULTS: Our analyses approach narrowed the candidates down to the novel missense variant in ARHGAP29 (GenBank: NM_004815.3, NP_004806.3;c.1654T>C [p.Ser552Pro]. A functional assay in zebrafish embryos showed that the encoded protein lacks the activity possessed by its wild-type counterpart, and migration assays revealed that keratinocytes transfected with wild-type ARHGAP29 migrated faster than counterparts transfected with the p.Ser552Pro ARHGAP29 variant or empty vector (control). CONCLUSION: These findings reveal ARHGAP29 to be a regulatory protein essential for proper development of the face, identifies an amino acid that is key for this, and provides a potential new diagnostic tool.Birth Defects Research 109:27-37, 2017. © 2016 Wiley Periodicals, Inc.

Replication of genome wide association identified candidate genes confirm the role of common and rare variants in PAX7 and VAX1 in the etiology of nonsyndromic CL(P).

Following recent genome wide association studies (GWAS), significant genetic associations have been identified for several genes with nonsyndromic cleft lip with or without cleft palate (CL(P)). To replicate two of these GWAS signals, we investigated the role of common and rare variants in the PAX7 and VAX1 genes. TaqMan genotyping was carried out for SNPs in VAX1 and PAX7 and transmission disequilibrium test (TDT) was performed to test for linkage and association in each population. Direct sequencing in and around the PAX7 and VAX1 genes in 1,326 individuals of European and Asian ancestry was done. The TDT analysis showed strong associations with markers in VAX1 (rs7078160, P = 2.7E-06 and rs475202, P = 0.0002) in a combined sample of Mongolian and Japanese CL(P) case-parent triads. Analyses using parent-of-origin effects showed significant excess transmission of the minor allele from both parents with the effect in the mothers (P = 6.5E-05, OR (transmission) = 1.91) more striking than in the fathers (P = 0.004, OR (transmission) = 1.67) for VAX1 marker rs7078160 in the combined Mongolian and Japanese samples when all cleft types were combined. The rs6659735 trinucleotide marker in PAX7 was significantly associated with all the US cleft groups combined (P = 0.007 in all clefts and P = 0.02 in CL(P)). Eight rare missense mutations found in PAX7 and two rare missense mutations in VAX1. Our study replicated previous GWAS findings for markers in VAX1 in the Asian population, and identified rare variants in PAX7 and VAX1 that may contribute to the etiology of CL(P). Determining the role of rare variants clearly warrants further investigation.

Analysis of MLL2 gene in the first Brazilian family with Kabuki syndrome.

Most patients with Kabuki syndrome (KS) are the only person in their family with the condition. However, familial cases of KS have been described showing evidence that this syndrome can be inherited as a dominant trait with variable expressivity. We report on two related individuals with facial findings characteristic of KS. The proposita had arched eyebrows, long and upward slanting palpebral fissures, cleft lip and palate, retromicrognathia, brachydactyly of hands and feet, stubby fingers, nail hypoplasia, and prominent finger pads. Her mother had eyebrows with dispersed lateral half, long and upward slanting palpebral fissures, retrognathia, abnormal and posteriorly rotated ears, prominent finger pads, brachydactyly of feet, learning difficulties, and psychomotor development delay. DNA sequencing revealed a novel missense mutation in the MLL2 gene in both the proposita and her mother. The mutation (p.R5432Q) was found in the exon 51, within the SET domain of the gene, which confers methyltransferase activity on the protein. Therefore, the epigenetic and transcriptional regulatory properties of this protein may be altered and this suggests that the mutation is the cause of phenotype observed in both the patient and her mother. The clinical signs and the molecular evidence in this family further support the notion that KS is an autosomal dominant condition with variable expressivity. To our knowledge this is the first report of a Brazilian family with recurrence of this syndrome.

Deficiency of the cytoskeletal protein SPECC1L leads to oblique facial clefting.

Genetic mutations responsible for oblique facial clefts (ObFC), a unique class of facial malformations, are largely unknown. We show that loss-of-function mutations in SPECC1L are pathogenic for this human developmental disorder and that SPECC1L is a critical organizer of vertebrate facial morphogenesis. During murine embryogenesis, Specc1l is expressed in cell populations of the developing facial primordial, which proliferate and fuse to form the face. In zebrafish, knockdown of a SPECC1L homolog produces a faceless phenotype with loss of jaw and facial structures, and knockdown in Drosophila phenocopies mutants in the integrin signaling pathway that exhibit cell-migration and -adhesion defects. Furthermore, in mammalian cells, SPECC1L colocalizes with both tubulin and actin, and its deficiency results in defective actin-cytoskeleton reorganization, as well as abnormal cell adhesion and migration. Collectively, these data demonstrate that SPECC1L functions in actin-cytoskeleton reorganization and is required for proper facial morphogenesis.

A case of 3q29 microdeletion syndrome involving oral cleft inherited from a nonaffected mosaic parent: molecular analysis and ethical implications.

OBJECTIVE: The objective of this study was to use array comparative genomic hybridization to detect causal microdeletions in samples of subjects with cleft lip and palate. SUBJECTS: We analyzed DNA samples from a male patient and his parents seen during surgical screening for an Operation Smile medical mission in the Philippines. METHOD: We used Affymetrix® Genome-Wide Human SNP Array 6.0 followed by sequencing and quantitative polymerase chain reaction using SYBR Green I dye. RESULTS: We report the second case of 3q29 microdeletion syndrome including cleft lip with or without cleft palate and the first case of this microdeletion syndrome inherited from a phenotypically normal mosaic parent. CONCLUSIONS: Our findings confirm the usefulness of a comparative genomic hybridization to detect causal microdeletions and indicate that parental somatic mosaicism should be considered in healthy parents for genetic counseling of the families. We discuss important ethical implications of sharing health impact results from research studies with the participant families.

Identification of a microdeletion at the 7q33-q35 disrupting the CNTNAP2 gene in a Brazilian stuttering case.

Speech and language disorders are some of the most common referral reasons to child development centers accounting for approximately 40% of cases. Stuttering is a disorder in which involuntary repetition, prolongation, or cessation of the sound precludes the flow of speech. About 5% of individuals in the general population have a stuttering problem, and about 80% of the affected children recover naturally. The causal factors of stuttering remain uncertain in most cases; studies suggest that genetic factors are responsible for 70% of the variance in liability for stuttering, whereas the remaining 30% is due to environmental effects supporting a complex cause of the disorder. The use of high-resolution genome wide array comparative genomic hybridization has proven to be a powerful strategy to narrow down candidate regions for complex disorders. We report on a case with a complex set of speech and language difficulties including stuttering who presented with a 10 Mb deletion of chromosome region 7q33-35 causing the deletion of several genes and the disruption of CNTNAP2 by deleting the first three exons of the gene. CNTNAP2 is known to be involved in the cause of language and speech disorders and autism spectrum disorder and is in the same pathway as FOXP2, another important language gene, which makes it a candidate gene for causal studies speech and language disorders such as stuttering.

Further characterization of microdeletion syndrome involving 2p15-p16.1.

We report on a patient presenting with cognitive delay, prenatal and postnatal growth deficiency, microcephaly, ptosis of eyelids, high and broad nasal root, and camptodactyly. Analysis of a dense whole genome single-nucleotide polymorphism (SNP) array showed a de novo 3.35 Mb deletion on 2p15-p16.1. In order to study the parental origin of the deletion we analyzed selected SNPs in the deleted area in the proband and her parents showing Mendelian incompatibilities suggesting a de novo deletion on the chromosome of paternal origin. Based on the five cases described previously in the literature, we have narrowed the critical region responsible for the 2p15-p16.1 microdeletion syndrome phenotype. The critical region does not include the VRK2 gene that had been speculated to have a role in cortical dysplasia. However, the association of the VRK2 gene with cortical dysplasia remains to be determined, as MRI imaging of the brain and gene content of the 2p15-16 deletion becomes established in more patients.

A genome-wide association study of cleft lip with and without cleft palate identifies risk variants near MAFB and ABCA4.

Case-parent trios were used in a genome-wide association study of cleft lip with and without cleft palate. SNPs near two genes not previously associated with cleft lip with and without cleft palate (MAFB, most significant SNP rs13041247, with odds ratio (OR) per minor allele = 0.704, 95% CI 0.635-0.778, P = 1.44 x 10(-11); and ABCA4, most significant SNP rs560426, with OR = 1.432, 95% CI 1.292-1.587, P = 5.01 x 10(-12)) and two previously identified regions (at chromosome 8q24 and IRF6) attained genome-wide significance. Stratifying trios into European and Asian ancestry groups revealed differences in statistical significance, although estimated effect sizes remained similar. Replication studies from several populations showed confirming evidence, with families of European ancestry giving stronger evidence for markers in 8q24, whereas Asian families showed stronger evidence for association with MAFB and ABCA4. Expression studies support a role for MAFB in palatal development.