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

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

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

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

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

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

Genome-wide association studies in dogs and humans identify ADAMTS20 as a risk variant for cleft lip and palate.

Cleft lip with or without cleft palate (CL/P) is the most commonly occurring craniofacial birth defect. We provide insight into the genetic etiology of this birth defect by performing genome-wide association studies in two species: dogs and humans. In the dog, a genome-wide association study of 7 CL/P cases and 112 controls from the Nova Scotia Duck Tolling Retriever (NSDTR) breed identified a significantly associated region on canine chromosome 27 (unadjusted p=1.1 x 10(-13); adjusted p= 2.2 x 10(-3)). Further analysis in NSDTR families and additional full sibling cases identified a 1.44 Mb homozygous haplotype (chromosome 27: 9.29 - 10.73 Mb) segregating with a more complex phenotype of cleft lip, cleft palate, and syndactyly (CLPS) in 13 cases. Whole-genome sequencing of 3 CLPS cases and 4 controls at 15X coverage led to the discovery of a frameshift mutation within ADAMTS20 (c.1360_1361delAA (p.Lys453Ilefs*3)), which segregated concordant with the phenotype. In a parallel study in humans, a family-based association analysis (DFAM) of 125 CL/P cases, 420 unaffected relatives, and 392 controls from a Guatemalan cohort, identified a suggestive association (rs10785430; p =2.67 x 10-6) with the same gene, ADAMTS20. Sequencing of cases from the Guatemalan cohort was unable to identify a causative mutation within the coding region of ADAMTS20, but four coding variants were found in additional cases of CL/P. In summary, this study provides genetic evidence for a role of ADAMTS20 in CL/P development in dogs and as a candidate gene for CL/P development in humans.

Interaction between IRF6 and TGFA genes contribute to the risk of nonsyndromic cleft lip/palate.

Previous evidence from tooth agenesis studies suggested IRF6 and TGFA interact. Since tooth agenesis is commonly found in individuals with cleft lip/palate (CL/P), we used four large cohorts to evaluate if IRF6 and TGFA interaction contributes to CL/P. Markers within and flanking IRF6 and TGFA genes were tested using Taqman or SYBR green chemistries for case-control analyses in 1,000 Brazilian individuals. We looked for evidence of gene-gene interaction between IRF6 and TGFA by testing if markers associated with CL/P were overtransmitted together in the case-control Brazilian dataset and in the additional family datasets. Genotypes for an additional 142 case-parent trios from South America drawn from the Latin American Collaborative Study of Congenital Malformations (ECLAMC), 154 cases from Latvia, and 8,717 individuals from several cohorts were available for replication of tests for interaction. Tgfa and Irf6 expression at critical stages during palatogenesis was analyzed in wild type and Irf6 knockout mice. Markers in and near IRF6 and TGFA were associated with CL/P in the Brazilian cohort (p<10(-6)). IRF6 was also associated with cleft palate (CP) with impaction of permanent teeth (p<10(-6)). Statistical evidence of interaction between IRF6 and TGFA was found in all data sets (p = 0.013 for Brazilians; p = 0.046 for ECLAMC; p = 10(-6) for Latvians, and p = 0.003 for the 8,717 individuals). Tgfa was not expressed in the palatal tissues of Irf6 knockout mice. IRF6 and TGFA contribute to subsets of CL/P with specific dental anomalies. Moreover, this potential IRF6-TGFA interaction may account for as much as 1% to 10% of CL/P cases. The Irf6-knockout model further supports the evidence of IRF6-TGFA interaction found in humans.

Follow-up association studies of chromosome region 9q and nonsyndromic cleft lip/palate.

Cleft lip/palate comprises a large fraction of all human birth defects, and is notable for its significant lifelong morbidity and complex etiology. Several studies have shown that genetic factors appear to play a significant role in the etiology of cleft lip/palate. Human chromosomal region 9q21 has been suggested in previous reports to contain putative cleft loci. Moreover, a specific region (9q22.3-34.1) was suggested to present a approximately 45% probability of harboring a cleft susceptibility gene. Fine mapping of 50 SNPs across the 9q22.3-34.11 region was performed to test for association with cleft lip/palate in families from United States, Spain, Turkey, Guatemala, and China. We performed family-based analyses and found evidence of association of cleft lip/palate with STOM (rs306796) in Guatemalan families (P = 0.004) and in all multiplex families pooled together (P = 0.002). This same SNP also showed borderline association in the US families (P = 0.04). Under a nominal value of 0.05, other SNPs also showed association with cleft lip/palate and cleft subgroups. SNPs in STOM and PTCH genes and nearby FOXE1 were further associated with cleft phenotypes in Guatemalan and Chinese families. Gene prioritization analysis revealed PTCH and STOM ranking among the top fourteen candidates for cleft lip/palate among 339 genes present in the region. Our results support the hypothesis that the 9q22.32-34.1 region harbors cleft susceptibility genes. Additional studies with other populations should focus on these loci to further investigate the participation of these genes in human clefting.

Genome scan, fine-mapping, and candidate gene analysis of non-syndromic cleft lip with or without cleft palate reveals phenotype-specific differences in linkage and association results.

OBJECTIVES: Non-syndromic orofacial clefts, i.e. cleft lip (CL) and cleft palate (CP), are among the most common birth defects. The goal of this study was to identify genomic regions and genes for CL with or without CP (CL/P). METHODS: We performed linkage analyses of a 10 cM genome scan in 820 multiplex CL/P families (6,565 individuals). Significant linkage results were followed by association analyses of 1,476 SNPs in candidate genes and regions, utilizing a weighted false discovery rate (wFDR) approach to control for multiple testing and incorporate the genome scan results. RESULTS: Significant (multipoint HLOD >or=3.2) or genome-wide-significant (HLOD >or=4.02) linkage results were found for regions 1q32, 2p13, 3q27-28, 9q21, 12p11, 14q21-24 and 16q24. SNPs in IRF6 (1q32) and in or near FOXE1 (9q21) reached formal genome-wide wFDR-adjusted significance. Further, results were phenotype dependent in that the IRF6 region results were most significant for families in which affected individuals have CL alone, and the FOXE1 region results were most significant in families in which some or all of the affected individuals have CL with CP. CONCLUSIONS: These results highlight the importance of careful phenotypic delineation in large samples of families for genetic analyses of complex, heterogeneous traits such as CL/P.

AXIS inhibition protein 2, orofacial clefts and a family history of cancer.

BACKGROUND: Cancer and congenital malformations occasionally may have a common etiology. The authors investigated whether families with one or more members affected by orofacial clefts (that is, families segregating orofacial clefts) had an increased cancer incidence when compared with control families. METHODS: The authors assessed 75 white families with nonsyndromic cleft lip with or without cleft palate (CL/P) and 93 white control families regarding a history of cancer. They used chi(2) and Fisher exact tests to determine significant differences. They then performed molecular studies with genes in which mutations have been independently associated with both cancer and craniofacial anomalies in a total of 111 families with CL/P. RESULTS: The families with CL/P reported a family history of cancer more often than did control families (P <.001), and they had higher rates of specific cancer types: colon (P <.001), brain (P = .003), leukemia (P = .005), breast (P = .009), prostate (P = .01), skin (P = .01), lung (P = .02) and liver (P = .02). The authors detected overtransmission of AXIS inhibition protein 2 (AXIN2) in CL/P probands (P = .003). CONCLUSION: Families segregating CL/P may have an increased susceptibility to cancer, notably colon cancer. Furthermore, AXIN2, a gene that when mutated increases susceptibility to colon cancer, also is associated with CL/P. CLINICAL IMPLICATIONS: People who are at a higher risk of developing disease need to adopt a healthier lifestyle, including avoiding exposure to risk factors that may interact with their genotypes.

The PDGF-C regulatory region SNP rs28999109 decreases promoter transcriptional activity and is associated with CL/P.

Human linkage and association studies suggest a gene(s) for nonsyndromic cleft lip with or without cleft palate (CL/P) on chromosome 4q31-q32 at or near the platelet-derived growth factor-C (PDGF-C) locus. The mouse Pdgfc(-/-) knockout shows that PDGF-C is essential for palatogenesis. To evaluate the role of PDGF-C in human clefting, we performed sequence analysis and SNP genotyping using 1048 multiplex CL/P families and 1000 case-control samples from multiple geographic origins. No coding region mutations were identified, but a novel -986 C>T SNP (rs28999109) was significantly associated with CL/P (P=0.01) in cases from Chinese families yielding evidence of linkage to 4q31-q32. Significant or near-significant association was also seen for this and several other PDGF-C SNPs in families from the United States, Spain, India, Turkey, China, and Colombia, whereas no association was seen in families from the Philippines, and Guatemala, and case-controls from Brazil. The -986T allele abolished six overlapping potential transcription regulatory motifs. Transfection assays of PDGF-C promoter reporter constructs show that the -986T allele is associated with a significant decrease (up to 80%) of PDGF-C gene promoter activity. This functional polymorphism acting on a susceptible genetic background may represent a component of human CL/P etiology.

Candidate gene/loci studies in cleft lip/palate and dental anomalies finds novel susceptibility genes for clefts.

PURPOSE: We revisited 42 families with two or more cleft-affected siblings who participated in previous studies. Complete dental information was collected to test the hypothesis that dental anomalies are part of the cleft phenotype spectrum, and can provide new opportunities for identification of cleft susceptibility genes. METHODS: Genotypes from 1489 single nucleotide polymorphism markers located in 150 candidate genes/loci were reanalyzed. Two sets of association analyses were carried out. First, we ran the analysis solely on the cleft status. Second, we assigned affection to any cleft or dental anomaly (tooth agenesis, supernumerary teeth, and microdontia) and repeated the analysis. RESULTS: Significant over-transmission was seen for a single nucleotide polymorphism in ankyrin repeat and sterile alpha motif domain containing 6 (rs4742741, 9q22.33; P = 0.0004) when a dental anomaly phenotype was included in the analysis. Significant over-transmission was also seen for a single nucleotide polymorphism in ERBB2 (rs1810132, 17q21.1; P = 0.0006). In the clefts only data, the most significant result was also for ERBB2 (P = 0.0006). Other markers with suggestive P values included interferon regulatory factor 6 and 6q21-q23 loci. In contrast to the above results, suggestive over-transmission of markers in GART, DPF3, and neurexin 3 were seen only when the dental anomaly phenotype was included in the analysis. CONCLUSIONS: These findings support the hypothesis that some loci may contribute to both clefts and congenital dental anomalies. Thus, including dental anomalies information in the genetics analysis of cleft lip and palate will provide new opportunities to map susceptibility loci for clefts.

Rethinking isolated cleft palate: evidence of occult lip defects in a subset of cases.

Emerging research suggests that subepithelial defects of the upper lip musculature are part of the phenotypic spectrum of cleft lip and/or palate (CL/P) and may represent an occult, subclinical manifestation of the anomaly. The present study investigates whether similar occult lip defects are present in individuals affected with isolated cleft palate (CP). To this end, upper lip ultrasounds of 33 CP cases (12 males, 21 females) were evaluated retrospectively for the presence of discontinuities (i.e., breaks) within the orbicularis oris muscle (OOM). In four CP cases (2 males, 2 females), distinct discontinuities of the OOM were identified. Of the remaining CP individuals, 23 demonstrated normal lip morphology on ultrasound (7 males, 16 females), while, in 6 cases (3 males, 3 females), a definitive evaluation was not possible. As CP and CL/P are traditionally thought to be etiologically distinct, these findings raise the possibility that some CP cases may be misclassified. Such diagnostic errors could have important implications for recurrence risk estimation and studies aimed at discovering etiology. (c) 2008 Wiley-Liss, Inc.

A genome wide linkage scan for cleft lip and palate and dental anomalies.

We revisited 46 families with two or more siblings affected with an orofacial cleft that participated in previous genome wide studies and collected complete dental information. Genotypes from 392 microsatellite markers at 10 cM intervals were reanalyzed. We carried out four sets of genome wide analyses. First, we ran the analysis solely on the cleft status. Second, we assigned to any dental anomaly (tooth agenesis, supernumerary teeth, and microdontia) an affection status, and repeated the analysis. Third, we ran only the 19 families where the proband had a cleft with no dental anomalies. Finally, we ran only the 27 families that had a proband with cleft and additional dental anomalies outside the cleft area. Chromosomes (1, 2, 6, 8, 16, and 19) presented regions with LOD scores >2.0. Chromosome 19 has the most compelling results in our study. The LOD scores increased from 3.11 (in the scan of all 46 families with clefts as the only assigned affection status) to 3.91 when the 19 families whose probands present with no additional dental anomalies were studied, suggesting the interval 19p13.12-19q12 may contain a gene that contributes to clefts but not to dental anomalies. On the other hand, we found a LOD score of 3.00 in the 2q22.3 region when dental anomalies data were added to the analysis to define affection status. Our preliminary results support the hypothesis that some loci may contribute to both clefts and congenital dental anomalies. Also, adding dental anomalies information will provide new opportunities to map susceptibility loci for clefts.

An ordered subset approach to including covariates in the transmission disequilibrium test.

Clinical heterogeneity of a disease may reflect an underlying genetic heterogeneity, which may hinder the detection of trait loci. Consequently, many statistical methods have been developed that allow for the detection of linkage and/or association signals in the presence of heterogeneity.This report describes the work of two parallel investigations into similar approaches to ordered subset analysis, based on an observed covariate, in the framework of family-based association analysis using Genetic Analysis Workshop 15 simulated data.With an appropriate choice of covariate, both approaches allow detection of two loci that are undetectable by the classical transmission-disequilibrium test. For a third locus, detectable by the classical transmission-disequilibrium test, a substantial increase of power of detection is shown.