Whole exome sequencing of distant relatives in multiplex families implicates rare variants in candidate genes for oral clefts

Rare variants found in multiplex families with orofacial clefts: Does expanding the phenotype make a difference?

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.

Integrative analysis of transcriptome dynamics during human craniofacial development identifies candidate disease genes

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.

A novel de novo TP63 mutation in whole-exome sequencing of a Syrian family with Oral cleft and ectrodactyly

BACKGROUND

Oral clefts and ectrodactyly are common, heterogeneous birth defects. We performed whole-exome sequencing (WES) analysis in a Syrian family. The proband presented with both orofacial clefting and ectrodactyly but not ectodermal dysplasia as typically seen in ectrodactyly, ectodermal dysplasia, and cleft lip/palate syndrome-3. A paternal uncle with only an oral cleft was deceased and unavailable for analysis.

METHODS

Variant annotation, Mendelian inconsistencies, and novel variants in known cleft genes were examined. Candidate variants were validated using Sanger sequencing, and pathogenicity assessed by knocking out the tp63 gene in zebrafish to evaluate its role during zebrafish development.

RESULTS

Twenty-eight candidate de novo events were identified, one of which is in a known oral cleft and ectrodactyly gene, TP63 (c.956G > T, p.Arg319Leu), and confirmed by Sanger sequencing.

CONCLUSION

TP63 mutations are associated with multiple autosomal dominant orofacial clefting and limb malformation disorders. The p.Arg319Leu mutation seen in this patient is de novo but also novel. Two known mutations in the same codon (c.956G > A, p.(Arg319His; rs121908839, c.955C > T), p.Arg319Cys) cause ectrodactyly, providing evidence that mutating this codon is deleterious. While this TP63 mutation is the best candidate for the patient's clinical presentation, whether it is responsible for the entire phenotype is unclear. Generation and characterization of tp63 knockout zebrafish showed necrosis and rupture of the head at 3 days post-fertilization (dpf). The embryonic phenotype could not be rescued by injection of zebrafish or human messenger RNA (mRNA). Further functional analysis is needed to determine what proportion of the phenotype is due to this mutation.

Genome Analysis Using Whole-Exome Sequencing of Non-Syndromic Cleft Lip and/or Palate from Malagasy Trios Identifies Variants Associated with Cilium-Related Pathways and Asian Genetic Ancestry

Background: Orofacial clefts (OFCs) are common congenital disabilities that can occur as isolated non-syndromic events or as part of Mendelian syndromes. OFC risk factors vary due to differences in regional environmental exposures, genetic variants, and ethnicities. In recent years, significant progress has been made in understanding OFCs, due to advances in sequencing and genotyping technologies. Despite these advances, very little is known about the genetic interplay in the Malagasy population. Methods: Here, we performed high-resolution whole-exome sequencing (WES) on non-syndromic cleft lip with or without palate (nCL/P) trios in the Malagasy population (78 individuals from 26 families (trios)). To integrate the impact of genetic ancestry admixture, we computed both global and local ancestries. Results: Participants demonstrated a high percentage of both African and Asian admixture. We identified damaging variants in primary cilium-mediated pathway genes WNT5B (one family), GPC4 (one family), co-occurrence in MSX1 (five families), WDR11 (one family), and tubulin stabilizer SEPTIN9 (one family). Furthermore, we identified an autosomal homozygous damaging variant in PHGDH (one family) gene that may impact metabiotic activity. Lastly, all variants were predicted to reside on local Asian genetic ancestry admixed alleles. Conclusion: Our results from examining the Malagasy genome provide limited support for the hypothesis that germline variants in primary cilia may be risk factors for nCL/P, and outline the importance of integrating local ancestry components better to understand the multi-ethnic impact on nCL/P.

Identification of Novel Risk Variants of Non-Syndromic Cleft Palate by Targeted Gene Panel Sequencing

Non-syndromic cleft palate (ns-CP) has a genetically heterogeneous aetiology. Numerous studies have suggested a crucial role of rare coding variants in characterizing the unrevealed component of genetic variation in ns-CP called the "missing heritability". Therefore, this study aimed to detect low-frequency variants that are implicated in ns-CP aetiology in the Polish population. For this purpose, coding regions of 423 genes associated with orofacial cleft anomalies and/or involved with facial development were screened in 38 ns-CP patients using the next-generation sequencing technology. After multistage selection and prioritisation, eight novel and four known rare variants that may influence an individual's risk of ns-CP were identified. Among detected alternations, seven were located in novel candidate genes for ns-CP, including COL17A1 (c.2435-1G>A), DLG1 (c.1586G>C, p.Glu562Asp), NHS (c.568G>C, p.Val190Leu-de novo variant), NOTCH2 (c.1997A>G, p.Tyr666Cys), TBX18 (c.647A>T, p.His225Leu), VAX1 (c.400G>A, p.Ala134Thr) and WNT5B (c.716G>T, p.Arg239Leu). The remaining risk variants were identified within genes previously linked to ns-CP, confirming their contribution to this anomaly. This list included ARHGAP29 (c.1706G>A, p.Arg569Gln), FLNB (c.3605A>G, Tyr1202Cys), IRF6 (224A>G, p.Asp75Gly-de novo variant), LRP6 (c.481C>A, p.Pro161Thr) and TP63 (c.353A>T, p.Asn118Ile). In summary, this study provides further insights into the genetic components contributing to ns-CP aetiology and identifies novel susceptibility genes for this craniofacial anomaly.

A Comprehensive Genetic Analysis of Slovenian Families with Multiple Cases of Orofacial Clefts Reveals Novel Variants in the Genes IRF6, GRHL3, and TBX22

Although the aetiology of non-syndromic orofacial clefts (nsOFCs) is usually multifactorial, syndromic OFCs (syOFCs) are often caused by single mutations in known genes. Some syndromes, e.g., Van der Woude syndrome (VWS1; VWS2) and X-linked cleft palate with or without ankyloglossia (CPX), show only minor clinical signs in addition to OFC and are sometimes difficult to differentiate from nsOFCs. We recruited 34 Slovenian multi-case families with apparent nsOFCs (isolated OFCs or OFCs with minor additional facial signs). First, we examined IRF6, GRHL3, and TBX22 by Sanger or whole exome sequencing to identify VWS and CPX families. Next, we examined 72 additional nsOFC genes in the remaining families. Variant validation and co-segregation analysis were performed for each identified variant using Sanger sequencing, real-time quantitative PCR and microarray-based comparative genomic hybridization. We identified six disease-causing variants (three novel) in IRF6, GRHL3, and TBX22 in 21% of families with apparent nsOFCs, suggesting that our sequencing approach is useful for distinguishing syOFCs from nsOFCs. The novel variants, a frameshift variant in exon 7 of IRF6, a splice-altering variant in GRHL3, and a deletion of the coding exons of TBX22, indicate VWS1, VWS2, and CPX, respectively. We also identified five rare variants in nsOFC genes in families without VWS or CPX, but they could not be conclusively linked to nsOFC.

Rare variants found in multiplex families with orofacial clefts: Does expanding the phenotype make a difference?

Whole-exome sequencing (WES) is now a relatively straightforward process to identify causal variants in Mendelian disorders. However, the same is not true for WES 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 whole-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, KLF4, SHROOM3, SMC3, TP63 , and TBX3 in seven 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.

A Retrospective Analysis of Clinically Focused Exome Sequencing Results of 372 Infants with Suspected Monogenic Disorders in China

Objective

The context was designed to optimize the diagnostic utility of clinically focused exome sequencing (CFES) and shorten the diagnostic odyssey among pediatric patients suspected of monogenic disorders (MDs).

Methods

Here, we retrospectively analyzed the clinical notes of 372 patients from different areas in the Jiangxi province that were referred for a diagnostic CFES and analysis from June 2018 to March 2022 with symptoms suggestive of MDs. In our study, preliminary tests using the proband-only clinical exome sequencing as a cost-effective first-tier diagnostic test for pediatric patients with unidentified MDs, supplemented by family segregation studies for targeted variants when indicated.

Results

Probands with confirmed diagnostic (CD) or likely diagnostic (LD) genetic influences accounted for 12% of all cases, whereas those with an uncertain diagnosis accounted for 48%. We also found that systemic primary carnitine deficiency (CDSP) (SLC22A5 gene) and phenylketonuria (PAH gene) were relatively more prevalent, and these patients with CDSP had the most frequent c.1400C > G variant (p.S467C) and c.51C > G variant (p. F17L) in this study. In addition, statistical analysis revealed that the estimates of diagnostic yields varied across certain phenotypic features of patients, and patients with specific phenotypic traits tended to benefit more from CFES.

Conclusion

The CFES may be a first-line genetic test for diagnosing young children with suspected genetic conditions, as it validates the identification of molecular genetics alterations and facilitates comprehensive medical management. Moreover, we found that infants exhibiting metabolism/homeostasis abnormalities, craniofacial /otolaryngology/ ophthalmologic abnormalities, and/or the integument were significantly more likely to receive a genetic diagnosis via CFES than infants without such features. However, due to the current study's low diagnostic yield and inherent limitations, high-quality clinical studies with larger sample sizes are still needed to provide more likely results and confirm our findings.

Using whole exome sequencing to identify susceptibility genes associated with nonsyndromic cleft lip with or without cleft palate

Cleft lip and palate is a common congenital birth defect in humans. Its incidence rate in China is as high as 1.82%, and is now a frequent deformity observed among the Chinese population; moreover, it varies across regions. Although the etiology of nonsyndromic cleft lip with or without cleft palate (NSCL/P) has been widely investigated, the results are inconsistent. The specific genes and mechanisms responsible for NSCL/P have not been fully understood. Whole exome sequencing (WES) is a new strategy for studying pathogenic genes. WES studies on NSCL/P have not been conducted in East China. Therefore, the aim of this study was to screen candidate genes of NSCL/P in East China using WES and analyze the temporal and spatial expressions of the candidate genes during embryonic palatal development. WES was performed in 30 children with NSCL/P from East China to screen candidate genes. A bioinformatics analysis was performed using commercially available software. Variants detected by WES were validated by immunohistochemistry and western blotting. After WES, 506,144 single-nucleotide variant sites were found. The results of database comparison, functional analysis, and mass spectrometry revealed that only the laminin alpha 5 (LAMA5) gene (site: rs145192286) was associated with NSCL/P. Immunohistochemistry results showed that LAMA5 expression in the medial edge epithelium changed with formation, lifting, and contact during palatogenesis. Almost no LAMA5 expression was detected in the palatal mesenchyme or after palatal fusion. Western blotting and immunohistochemistry results showed consistent trends. In conclusion, the WES results shows that the mutation at the site (rs145192286) of LAMA5 is associated with NSCL/P. The temporal and spatial expressions of LAMA5 during palatal development further demonstrate the involvement of this gene. Therefore, we speculate that LAMA5 is a new candidate pathogenic gene of NSCL/P. The identification of new pathogenic genes would help elucidate the pathogenesis of NSCL/P and provide a scientific basis for the prenatal diagnosis, prevention, and treatment of NSCL/P.

Identification of novel susceptibility genes for non-syndromic cleft lip with or without cleft palate using NGS-based multigene panel testing

For non-syndromic cleft lip with or without cleft palate (ns-CL/P), the proportion of heritability explained by the known risk loci is estimated to be about 30% and is captured mainly by common variants identified in genome-wide association studies. To contribute to the explanation of the "missing heritability" problem for orofacial clefts, a candidate gene approach was taken to investigate the potential role of rare and private variants in the ns-CL/P risk. Using the next-generation sequencing technology, the coding sequence of a set of 423 candidate genes was analysed in 135 patients from the Polish population. After stringent multistage filtering, 37 rare coding and splicing variants of 28 genes were identified. 35% of these genetic alternations that may play a role of genetic modifiers influencing an individual's risk were detected in genes not previously associated with the ns-CL/P susceptibility, including COL11A1, COL17A1, DLX1, EFTUD2, FGF4, FGF8, FLNB, JAG1, NOTCH2, SHH, WNT5A and WNT9A. Significant enrichment of rare alleles in ns-CL/P patients compared with controls was also demonstrated for ARHGAP29, CHD7, COL17A1, FGF12, GAD1 and SATB2. In addition, analysis of panoramic radiographs of patients with identified predisposing variants may support the hypothesis of a common genetic link between orofacial clefts and dental abnormalities. In conclusion, our study has confirmed that rare coding variants might contribute to the genetic architecture of ns-CL/P. Since only single predisposing variants were identified in novel cleft susceptibility genes, future research will be required to confirm and fully understand their role in the aetiology of ns-CL/P.

Monogenic Versus Multifactorial Inheritance in the Development of Isolated Cleft Palate: A Whole Genome Sequencing Study

Craniofacial morphogenesis is highly complex, as is the anatomical region involved. Errors during this process, resulting in orofacial clefts, occur in more than 400 genetic syndromes. Some cases of cleft lip and/or palate (CLP) are caused by mutations in single genes; however, complex interactions between genetic and environmental factors are considered to be responsible for the majority of non-syndromic CLP development. The aim of the current study was to identify genetic risk factors in patients with isolated cleft palate (CP) by whole genome sequencing. Patients with isolated CP (n = 30) recruited from the Riga Cleft Lip and Palate Centre, Institute of Stomatology, Riga, were analyzed by whole genome sequencing. Pathogenic or likely pathogenic variants were discovered in genes associated with CP (TBX22, COL2A1, FBN1, PCGF2, and KMT2D) in five patients; hence, rare disease variants were identified in 17% of patients with non-syndromic isolated CP. Our results were relevant to routine genetic counselling practice and genetic testing recommendations. Based on our data, we propose that all newborns with orofacial clefts should be offered genetic testing, at least for a panel of known CLP genes. Only if the results are negative and there is no suggestive family history or additional clinical symptoms (which would support additional exome or genome-wide investigation), should multifactorial empiric recurrence risk prediction tools be applied for families.

Genetic and Epigenetic Studies in Nonsyndromic Oral Clefts

The etiology of non-syndromic oral clefts (NSOFC) is complex with genetics, genomics, epigenetics and stochastics factors playing a role. Several approaches have been applied to understand the etiology of non-syndromic oral clefts. These include linkage, candidate gene association studies, genome-wide association studies, whole genome sequencing, copy number variations and epigenetics. In this review we shared these approaches, genes and loci reported in some studies.

To Stick or Not to Stick: Adhesions in Orofacial Clefts

Morphogenesis requires a tight coordination between mechanical forces and biochemical signals to inform individual cellular behavior. For these developmental processes to happen correctly the organism requires precise spatial and temporal coordination of the adhesion, migration, growth, differentiation, and apoptosis of cells originating from the three key embryonic layers, namely the ectoderm, mesoderm, and endoderm. The cytoskeleton and its remodeling are essential to organize and amplify many of the signaling pathways required for proper morphogenesis. In particular, the interaction of the cell junctions with the cytoskeleton functions to amplify the behavior of individual cells into collective events that are critical for development. In this review we summarize the key morphogenic events that occur during the formation of the face and the palate, as well as the protein complexes required for cell-to-cell adhesions. We then integrate the current knowledge into a comprehensive review of how mutations in cell-to-cell adhesion genes lead to abnormal craniofacial development, with a particular focus on cleft lip with or without cleft palate.

Mouse models in palate development and orofacial cleft research: Understanding the crucial role and regulation of epithelial integrity in facial and palate morphogenesis

Cleft lip and cleft palate are common birth defects resulting from genetic and/or environmental perturbations of facial development in utero. Facial morphogenesis commences during early embryogenesis, with cranial neural crest cells interacting with the surface ectoderm to form initially partly separate facial primordia consisting of the medial and lateral nasal prominences, and paired maxillary and mandibular processes. As these facial primordia grow around the primitive oral cavity and merge toward the ventral midline, the surface ectoderm undergoes a critical differentiation step to form an outer layer of flattened and tightly connected periderm cells with a non-stick apical surface that prevents epithelial adhesion. Formation of the upper lip and palate requires spatiotemporally regulated inter-epithelial adhesions and subsequent dissolution of the intervening epithelial seam between the maxillary and medial/lateral nasal processes and between the palatal shelves. Proper regulation of epithelial integrity plays a paramount role during human facial development, as mutations in genes encoding epithelial adhesion molecules and their regulators have been associated with syndromic and non-syndromic orofacial clefts. In this chapter, we summarize mouse genetic studies that have been instrumental in unraveling the mechanisms regulating epithelial integrity and periderm differentiation during facial and palate development. Since proper epithelial integrity also plays crucial roles in wound healing and cancer, understanding the mechanisms regulating epithelial integrity during facial development have direct implications for improvement in clinical care of craniofacial patients.

Genetic models and approaches to study orofacial clefts

INTRODUCTION

Orofacial clefts (OFCs) are common craniofacial birth defects with heterogeneous phenotype and etiology. Geneticists have applied nearly every available method and technology to further our 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 towards 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 distinctions traditional categories used in genetic studies and clinical medicine.

Identification of Novel Genomic Variations in Susceptibility to Nonsyndromic Cleft Lip and Palate Patients

BACKGROUND

Nonsyndromic cleft lip with or without palate (NSCL/P) is a multifactorial and common birth malformation caused by genetic and environmental factors, as well as by teratogens. Genome-wide association studies found genetic variations with modulatory effects of NSCL/P formation in Chinese and Iranian populations. We aimed to identify the susceptibility of single-nucleotide polymorphisms (SNPs) to nonsyndromic cleft lip with or without palate in the Indian population.

MATERIAL AND METHODS

The present study was conducted on NSCL/P cases and controls. Genomic DNA was extracted from peripheral blood and Axiom- Precision Medicine Research Array (PMRA) was performed. The Axiom-PMRA covers 902,527 markers and several thousand novel risk variants. Quality control-passed samples were included for candidate genetic variation identification, gene functional enrichment, and pathway and network analysis.

RESULTS

The genome-wide association study identified fourteen novel candidate gene SNPs that showed the most significant association with the risk of NSCL/P, and eight were predicted to have regulatory sequences.

CONCLUSION

The GWAS study showed novel candidate genetic variations in NSCL/P formations. These findings contribute to the understanding of genetic predisposition to nonsyndromic cleft lip with or without palate.

Identification of Novel Variants in Cleft Palate-Associated Genes in Brazilian Patients With Non-syndromic Cleft Palate Only

The identification of genetic risk factors for non-syndromic oral clefts is of great importance for better understanding the biological processes related to this heterogeneous and complex group of diseases. Herein we applied whole-exome sequencing to identify potential variants related to non-syndromic cleft palate only (NSCPO) in the multiethnic Brazilian population. Thirty NSCPO samples and 30 sex- and genetic ancestry-matched healthy controls were pooled (3 pools with 10 samples for each group) and subjected to whole-exome sequencing. After filtering, the functional affects, individually and through interactions, of the selected variants and genes were assessed by bioinformatic analyses. As a group, 399 variants in 216 genes related to palatogenesis/cleft palate, corresponding to 6.43%, were exclusively identified in the NSCPO pools. Among those genes are 99 associated with syndromes displaying cleft palate in their clinical spectrum and 92 previously related to cleft lip palate. The most significantly biological processes and pathways overrepresented in the NSCPO-identified genes were associated with the folic acid metabolism, highlighting the interaction between LDL receptor-related protein 6 (LRP6) and 5-methyltetrahydrofolate-homocysteine methyltransferase (MTR) that interconnect two large networks. This study yields novel data on characterization of specific variants and complex processes and pathways related to NSCPO, including many variants in genes of the folate/homocysteine pathway, and confirms that variants in genes related to syndromic cleft palate and cleft lip-palate may cause NSCPO.

Cleft lip/palate and educational attainment: cause, consequence or correlation? A Mendelian randomization study

BACKGROUND

Previous studies have found that children born with a non-syndromic orofacial cleft have lower-than-average educational attainment. Differences could be due to a genetic predisposition to low intelligence and academic performance, factors arising due to the cleft phenotype (such as social stigmatization, impaired speech/language development) or confounding by the prenatal environment. A clearer understanding of this mechanism will inform interventions to improve educational attainment in individuals born with a cleft, which could substantially improve their quality of life. We assessed evidence for the hypothesis that common variant genetic liability to non-syndromic cleft lip with or without cleft palate (nsCL/P) influences educational attainment.

METHODS

We performed a genome-wide association study (GWAS) meta-analysis of nsCL/P with 1692 nsCL/P cases and 4259 parental and unrelated controls. Using GWAS summary statistics, we performed Linkage Disequilibrium (LD)-score regression to estimate the genetic correlation between nsCL/P, educational attainment (GWAS n = 766 345) and intelligence (GWAS n = 257 828). We used two-sample Mendelian randomization to evaluate the causal effects of genetic liability to nsCL/P on educational attainment and intelligence.

RESULTS

There was limited evidence for shared genetic aetiology or causal relationships between nsCL/P and educational attainment [genetic correlation (rg) -0.05, 95% confidence interval (CI) -0.12 to 0.01, P 0.13; MR estimate (βMR) -0.002, 95% CI -0.009 to 0.006, P 0.679) or intelligence (rg -0.04, 95% CI -0.13 to 0.04, P 0.34; βMR -0.009, 95% CI -0.02 to 0.002, P 0.11).

CONCLUSIONS

Common variants are unlikely to predispose individuals born with nsCL/P to low educational attainment or intelligence. This is an important first step towards understanding the aetiology of low educational attainment in this group.

The PAX1 locus at 20p11 is a potential genetic modifier for bilateral cleft lip

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.

Co-occurrence of orofacial clefts and clubfoot phenotypes in a sub-Saharan African cohort: Whole-exome sequencing implicates multiple syndromes and genes

Background

Orofacial clefts (OFCs) are congenital malformations of the face and palate, with an incidence of 1 per 700 live births. Clubfoot or congenital talipes equinovarus (CTEV) is a three‐dimensional abnormality of the leg, ankle, and feet that leads to the anomalous positioning of foot and ankle joints and has an incidence of 1 per 1000 live births. OFCs and CTEV may occur together or separately in certain genetic syndromes in addition to other congenital abnormalities. Here, we sought to decipher the genetic etiology of OFC and CTEV that occurred together in six probands.

Methods

At the time of recruitment, the most clinically obvious congenital anomalies in these individuals were the OFC and CTEV. We carried out whole‐exome sequencing (WES) on DNA samples from probands and available parents employing the Agilent SureSelect XT kit and Illumina HiSeq2500 platform, followed by bioinformatics analyses. WES variants were validated by clinical Sanger Sequencing.

Results

Of the six probands, we observed probable pathogenic genetic variants in four. In three probands with probable pathogenic genetic variants, each individual had variants in three different genes, whereas one proband had probable pathogenic variant in just one gene. In one proband, we observed variants in DIS3L2, a gene associated with Perlman syndrome. A second proband had variants in EPG5 (associated with Vici Syndrome), BARX1 and MKI67, while another proband had potentially etiologic variants in FRAS1 (associated with Fraser Syndrome 1), TCOF1 (associated with Treacher Collins Syndrome 1) and MKI67. The last proband had variants in FRAS1, PRDM16 (associated with Cardiomyopathy, dilated, 1LL/Left ventricular noncompaction 8) and CHD7 (associated with CHARGE syndrome/Hypogonadotropic hypogonadism 5 with or without anosmia).

Conclusion

Our results suggest that clubfoot and OFCs are two congenital abnormalities that can co‐occur in certain individuals with varying genetic causes and expressivity, warranting the need for deep phenotyping.

A Novel IRF6 Variant Detected in a Family With Nonsyndromic Cleft Lip and Palate by Whole Exome Sequencing

ABSTRACT

Nonsyndromic cleft lip with or without cleft palate (NSCL/P) is the most common congenital craniofacial malformation, and its harmful effects on affected individuals and families are apparent. The causative genes and their mechanisms are not completely clear, although several studies have been conducted. Accordingly, in the present study, we recruited a Han Chinese family with hereditary NSCL/P to explore the possible causative variants of this disease using whole exome sequencing. Bioinformatics screening and analysis, mutation function prediction, species conservation analysis, and homology protein modeling were used to identify the variants and evaluate their influence. A mutation in the interferon regulatory factor 6 (IRF6) gene (c.961C>T; p.Val321Met) was detected as a candidate causative variant and predicted to be deleterious. The codon was found to be conserved in many species, and the residue change caused by this mutation changed the structure of IRF6 to a certain degree. The findings suggest that this IRF6 variant is probably the pathogenic cause of NSCL/P in this family. Our results further provide evidence that IRF6 variants play a role in the etiology of NSCL/P.

A novel splicing mutation of ARHGAP29 is associated with nonsyndromic cleft lip with or without cleft palate

BACKGROUND

Nonsyndromic cleft lip with or without cleft palate (NSCL/P) is one of the most common birth defects, and occurs in approximately 1/700 live births worldwide. The correlation between the ABCA4-ARHGAP29 region and NSCL/P was first identified by genome-wide association studies (GWAS), but few reports have examined NSCL/P caused by ARHGAP29 mutations in the Chinese population.

METHODS

We performed chromosome microarray analysis (CMA) for two consecutive abnormal fetuses and whole exome sequencing (WES) for the family, including 3 patients and 2 normal family members, Sanger sequencing and RT-PCR were used to confirm the mutation.

RESULTS

We identified a novel splice donor mutation (ARHGAP29 c.1920 + 1G > A) in two consecutive NSCL/P fetuses, and the variant was inherited from the mother and grandfather. The mutation caused abnormal skipping of exon 17, and the mRNA level of ARHGAP29 was significantly decreased compared to the wild type.

CONCLUSIONS

In this study, we successfully diagnosed the genetic cause of NSCL/P in a family and first report that the c.1920 + 1G > A mutation in ARHGAP29 is associated with NSCL/P. Our study enriches the genetic landscape of NSCL/P, extends the mutation spectrum of ARHGAP29, and provides a new direction for the diagnosis of NSCL/P in patients and its prenatal diagnosis in fetuses.

Genome-wide Enrichment of De Novo Coding Mutations in Orofacial Cleft Trios

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.

Detection of rare disease variants in extended pedigrees using RVS

Summary

Family-based sequencing studies enable researchers to identify highly penetrant genetic variants too rare to be tested in conventional case-control studies, by studying co-segregation of variant and disease phenotypes. When multiple affected subjects in a family are sequenced, the probability that a variant or a set of variants is shared identical-by-descent by some or all affected relatives provides evidence against the null hypothesis of complete absence of linkage and association. The Rare Variant Sharing software package RVS implements a suite of tools to assess association and linkage between rare genetic variants and a dichotomous disease indicator in family pedigrees.

Availability and Implementation

RVS is available as open source software from the Bioconductor webpage at https://bioconductor.org/packages/release/bioc/html/RVS.html.

Supplementary information

Supplementary data are available at Bioinformatics online.

CDH1 Mutation Distribution and Type Suggests Genetic Differences between the Etiology of Orofacial Clefting and Gastric Cancer

Pathogenic variants in CDH1, encoding epithelial cadherin (E-cadherin), have been implicated in hereditary diffuse gastric cancer (HDGC), lobular breast cancer, and both syndromic and non-syndromic cleft lip/palate (CL/P). Despite the large number of CDH1 mutations described, the nature of the phenotypic consequence of such mutations is currently not able to be predicted, creating significant challenges for genetic counselling. This study collates the phenotype and molecular data for available CDH1 variants that have been classified, using the American College of Medical Genetics and Genomics criteria, as at least ‘likely pathogenic’, and correlates their molecular and structural characteristics to phenotype. We demonstrate that CDH1 variant type and location differ between HDGC and CL/P, and that there is clustering of CL/P variants within linker regions between the extracellular domains of the cadherin protein. While these differences do not provide for exact prediction of the phenotype for a given mutation, they may contribute to more accurate assessments of risk for HDGC or CL/P for individuals with specific CDH1 variants.

Craniofacial Analysis May Indicate Co-Occurrence of Skeletal Malocclusions and Associated Risks in Development of Cleft Lip and Palate

Non-syndromic orofacial clefts encompass a range of morphological changes affecting the oral cavity and the craniofacial skeleton, of which the genetic and epigenetic etiologic factors remain largely unknown. The objective of this study is to explore the contribution of underlying dentofacial deformities (also known as skeletal malocclusions) in the craniofacial morphology of non-syndromic cleft lip and palate patients (nsCLP). For that purpose, geometric morphometric analysis was performed using full skull cone beam computed tomography (CBCT) images of patients with nsCLP (n = 30), normocephalic controls (n = 60), as well as to sex- and ethnicity- matched patients with an equivalent dentofacial deformity (n = 30). Our outcome measures were shape differences among the groups quantified via principal component analysis and associated principal component loadings, as well as mean shape differences quantified via a Procrustes distance among groups. According to our results, despite the shape differences among all three groups, the nsCLP group shares many morphological similarities in the maxilla and mandible with the dentofacial deformity group. Therefore, the dentoskeletal phenotype in nsCLP could be the result of the cleft and the coexisting dentofacial deformity and not simply the impact of the cleft.

PBX-WNT-P63-IRF6 pathway in nonsyndromic cleft lip and palate

Nonsyndromic cleft lip and palate (NSCLP) is one of the most common craniofacial anomalies in humans, affecting more than 135,000 newborns worldwide. NSCLP has a multifactorial etiology with more than 50 genes postulated to play an etiologic role. The genetic pathway comprised of Pbx-Wnt-p63-Irf6 genes was shown to control facial morphogenesis in mice and proposed as a regulatory pathway for NSCLP. Based on these findings, we investigated whether variation in PBX1, PBX2, and TP63, and their proposed interactions were associated with NSCLP. Fourteen single nucleotide variants (SNVs) in/nearby PBX1, PBX2, and TP63 were genotyped in 780 NSCLP families of nonHispanic white (NHW) and Hispanic ethnicities. Family-based association tests were performed for individual SNVs stratified by ethnicity and family history of NSCLP. Gene-gene interactions were also tested. A significant association was found for PBX2 rs3131300 and NSCLP in combined Hispanic families (p = .003) while nominal association was found for TP63 rs9332461 in multiplex Hispanic families (p = .005). Significant haplotype associations were observed for PBX2 in NHW (p = .0002) and Hispanic families (p = .003), and for TP63 in multiplex Hispanic families (.003). An independent case-control group was used to validate findings, and significant associations were found with PBX1 rs6426870 (p = .007) and TP63 rs9332461 (p = .03). Gene-gene interactions were detected between PBX1/PBX2/TP63 with IRF6 in NHW families, and between PBX1 with WNT9B in both NHW and Hispanic families (p < .0018). This study provides the first evidence for a role of PBX1 and PBX2, additional evidence for the role of TP63, and support for the proposed PBX-WNT-TP63-IRF6 regulatory pathway in the etiology of NSCLP.

Non-Syndromic Cleft Lip with or without Cleft Palate: Genome-Wide Association Study in Europeans Identifies a Suggestive Risk Locus at 16p12.1 and Supports SH3PXD2A as a Clefting Susceptibility Gene

Non-syndromic cleft lip with or without cleft palate (nsCL/P) ranks among the most common human congenital malformations, and has a multifactorial background in which both exogenous and genetic risk factors act in concert. The present report describes a genome-wide association study (GWAS) involving a total of 285 nsCL/P patients and 1212 controls from the Netherlands and Belgium. Twenty of the 40 previously reported nsC/LP susceptibility loci were replicated, which underlined the validity of this sample. SNV-based analysis of the data identified an as yet unreported suggestive locus at chromosome 16p12.1 (p-value of the lead SNV: 4.17 × 10^-7). This association was replicated in two of three patient/control replication series (Central European and Yemeni). Gene analysis of the GWAS data prioritized SH3PXD2A at chromosome 10q24.33 as a candidate gene for nsCL/P. To date, support for this gene as a cleft gene has been restricted to data from zebrafish and a knockout mouse model. The present GWAS was the first to implicate SH3PXD2A in non-syndromic cleft formation in humans. In summary, although performed in a relatively small sample, the present GWAS generated novel insights into nsCL/P etiology.

Likely Pathogenic Variants in One Third of Non-Syndromic Discontinuous Cleft Lip and Palate Patients

Oral clefts are composed of cleft of the lip, cleft of the lip and palate, or cleft of the palate, and they are associated with a wide range of expression and severity. When cleft of the palate is associated with cleft of the lip with preservation of the primary palate, it defines an atypical phenotype called discontinuous cleft. Although this phenotype may represent 5% of clefts of the lip and/or palate (CLP), it is rarely specifically referred to and its pathophysiology is unknown. We conducted whole exome sequencing (WES) and apply a candidate gene approach to non-syndromic discontinuous CLP individuals in order to identify genes and deleterious variants that could underlie this phenotype. We discovered loss-of-function variants in two out of the seven individuals, implicating FGFR1 and DLG1 genes, which represents almost one third of this cohort. Whole exome sequencing of clinically well-defined subgroups of CLP, such as discontinuous cleft, is a relevant approach to study CLP etiopathogenesis. It could facilitate more accurate clinical, epidemiological and fundamental research, ultimately resulting in better diagnosis and care of CLP patients. Non-syndromic discontinuous cleft lip and palate seems to have a strong genetic basis.

A Novel CDH1 Mutation Causing Reduced E-Cadherin Dimerization Is Associated with Nonsyndromic Cleft Lip With or Without Cleft Palate

Aims: Cleft lip with or without cleft palate (CL/P) is a common birth defect with an average prevalence of 1/700 to 1/1000. Almost 70% of CL/P cases are nonsyndromic CL/P (NSCL/P). The aim of this study was to identify the underlying cause of a four-generation Chinese family with autosomal dominant NSCL/P. Methods: Genomic DNA was extracted from peripheral blood leukocytes, and whole-exome sequencing was carried out to identify the underlying genetic cause of the disorder. The mutation was confirmed by Sanger sequencing and polymerase chain reaction-restriction fragment length polymorphism methods. Western blotting and coimmunoprecipitation were used to analyze the protein expression level and adhesive dimerization of the CDH1 mutants. Slow aggregation assays were conducted to investigate the cell-cell adhesion ability. Results: A novel missense mutation (c.468G>C/p.Trp156Cys) of CDH1 was identified in the proband and the mutation was shown to cosegregate with the phenotype in the family. Furthermore, we found that the p.Trp156Cys mutation led to decreased E-cadherin dimerization and cell-cell adhesion ability. Conclusions: Our findings identified a novel CDH1 variant (c.468G>C/p.Trp156Cys) responsible for NSCL/P in a Chinese family, which expanded the mutational spectrum of the CDH1 gene and may contribute to understanding the molecular basis of NSCL/P.

High-Resolution Epigenomic Atlas of Human Embryonic Craniofacial Development

Defects in patterning during human embryonic development frequently result in craniofacial abnormalities. The gene regulatory programs that build the craniofacial complex are likely controlled by information located between genes and within intronic sequences. However, systematic identification of regulatory sequences important for forming the human face has not been performed. Here, we describe comprehensive epigenomic annotations from human embryonic craniofacial tissues and systematic comparisons with multiple tissues and cell types. We identified thousands of tissue-specific craniofacial regulatory sequences and likely causal regions for rare craniofacial abnormalities. We demonstrate significant enrichment of common variants associated with orofacial clefting in enhancers active early in embryonic development, while those associated with normal facial variation are enriched near the end of the embryonic period. These data are provided in easily accessible formats for both craniofacial researchers and clinicians to aid future experimental design and interpretation of noncoding variation in those affected by craniofacial abnormalities.

Gene-gene interaction among cell adhesion genes and risk of nonsyndromic cleft lip with or without cleft palate in Chinese case-parent trios

Background

Nonsyndromic cleft lip with or without cleft palate (NSCL/P) is a common birth defect with complex etiology. One strategy for studying the genetic risk factors of NSCL/P is to consider gene–gene interaction (G × G) among gene pathways having a role in craniofacial development. The present study aimed to investigate the G × G among cell adhesion gene pathway.

Methods

We carried out an interaction analysis of eight genes involved in cell adherens junctions among 806 NSCL/P Chinese case‐parent trios originally recruited for a genome‐wide association study (GWAS). Regression‐based approach was used to test for two‐way G × G interaction, while machine learning algorithm was run for exploring both two‐way and multi‐way interaction that may affect the risk of NSCL/P.

Results

A two‐way ACTN1 × CTNNB1 interaction reached the adjusted significance level. The single nucleotide polymorphisms pair composed of rs17252114 (CTNNB1) and rs1274944 (ACTN1) yielded a p value of .0002, and this interaction was also supported by the logic regression algorithm. Higher order interactions involving ACTN1, CTNNB1, and CDH1 were picked out by logic regression, suggesting a potential role in NSCL/P risk.

Conclusion

This study suggests for the first time evidence of both two‐way and multi‐way G × G interactions among cell adhesion genes contributing to the NSCL/P risk.

Rediscovering the value of families for psychiatric genetics research

As it is likely that both common and rare genetic variation are important for complex disease risk, studies that examine the full range of the allelic frequency distribution should be utilized to dissect the genetic influences on mental illness. The rate limiting factor for inferring an association between a variant and a phenotype is inevitably the total number of copies of the minor allele captured in the studied sample. For rare variation, with minor allele frequencies of 0.5% or less, very large samples of unrelated individuals are necessary to unambiguously associate a locus with an illness. Unfortunately, such large samples are often cost prohibitive. However, by using alternative analytic strategies and studying related individuals, particularly those from large multiplex families, it is possible to reduce the required sample size while maintaining statistical power. We contend that using whole genome sequence (WGS) in extended pedigrees provides a cost-effective strategy for psychiatric gene mapping that complements common variant approaches and WGS in unrelated individuals. This was our impetus for forming the "Pedigree-Based Whole Genome Sequencing of Affective and Psychotic Disorders" consortium. In this review, we provide a rationale for the use of WGS with pedigrees in modern psychiatric genetics research. We begin with a focused review of the current literature, followed by a short history of family-based research in psychiatry. Next, we describe several advantages of pedigrees for WGS research, including power estimates, methods for studying the environment, and endophenotypes. We conclude with a brief description of our consortium and its goals.

Genomic analyses in African populations identify novel risk loci for cleft palate

Orofacial clefts are common developmental disorders that pose significant clinical, economical and psychological problems. We conducted genome-wide association analyses for cleft palate only (CPO) and cleft lip with or without palate (CL/P) with ~17 million markers in sub-Saharan Africans. After replication and combined analyses, we identified novel loci for CPO at or near genome-wide significance on chromosomes 2 (near CTNNA2) and 19 (near SULT2A1). In situ hybridization of Sult2a1 in mice showed expression of SULT2A1 in mesenchymal cells in palate, palatal rugae and palatal epithelium in the fused palate. The previously reported 8q24 was the most significant locus for CL/P in our study, and we replicated several previously reported loci including PAX7 and VAX1.

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

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

Wnt signaling in orofacial clefts: crosstalk, pathogenesis and models

Diverse signaling cues and attendant proteins work together during organogenesis, including craniofacial development. Lip and palate formation starts as early as the fourth week of gestation in humans or embryonic day 9.5 in mice. Disruptions in these early events may cause serious consequences, such as orofacial clefts, mainly cleft lip and/or cleft palate. Morphogenetic Wnt signaling, along with other signaling pathways and transcription regulation mechanisms, plays crucial roles during embryonic development, yet the signaling mechanisms and interactions in lip and palate formation and fusion remain poorly understood. Various Wnt signaling and related genes have been associated with orofacial clefts. This Review discusses the role of Wnt signaling and its crosstalk with cell adhesion molecules, transcription factors, epigenetic regulators and other morphogenetic signaling pathways, including the Bmp, Fgf, Tgfβ, Shh and retinoic acid pathways, in orofacial clefts in humans and animal models, which may provide a better understanding of these disorders and could be applied towards prevention and treatments.

Inferring disease risk genes from sequencing data in multiplex pedigrees through sharing of rare variants

We previously demonstrated how sharing of rare variants (RVs) in distant affected relatives can be used to identify variants causing a complex and heterogeneous disease. This approach tested whether single RVs were shared by all sequenced affected family members. However, as with other study designs, joint analysis of several RVs (e.g., within genes) is sometimes required to obtain sufficient statistical power. Further, phenocopies can lead to false negatives for some causal RVs if complete sharing among affected is required. Here, we extend our methodology (Rare Variant Sharing, RVS) to address these issues. Specifically, we introduce gene-based analyses, a partial sharing test based on RV sharing probabilities for subsets of affected relatives and a haplotype-based RV definition. RVS also has the desirable feature of not requiring external estimates of variant frequency or control samples, provides functionality to assess and address violations of key assumptions, and is available as open source software for genome-wide analysis. Simulations including phenocopies, based on the families of an oral cleft study, revealed the partial and complete sharing versions of RVS achieved similar statistical power compared with alternative methods (RareIBD and the Gene-Based Segregation Test), and had superior power compared with the pedigree Variant Annotation, Analysis, and Search Tool (pVAAST) linkage statistic. In studies of multiplex cleft families, analysis of rare single nucleotide variants in the exome of 151 affected relatives from 54 families revealed no significant excess sharing in any one gene, but highlighted different patterns of sharing revealed by the complete and partial sharing tests.

Clinical spectrum and pleiotropic nature of CDH1 germline mutations

CDH1 encodes E-cadherin, a key protein in adherens junctions. Given that E-cadherin is involved in major cellular processes such as embryogenesis and maintenance of tissue architecture, it is no surprise that deleterious effects arise from its loss of function. E-cadherin is recognised as a tumour suppressor gene, and it is well established that CDH1 genetic alterations cause diffuse gastric cancer and lobular breast cancer—the foremost manifestations of the hereditary diffuse gastric cancer syndrome. However, in the last decade, evidence has emerged demonstrating that CDH1 mutations can be associated with lobular breast cancer and/or several congenital abnormalities, without any personal or family history of diffuse gastric cancer. To date, no genotype–phenotype correlations have been observed. Remarkably, there are reports of mutations affecting the same nucleotide but inducing distinct clinical outcomes. In this review, we bring together a comprehensive analysis of CDH1-associated disorders and germline alterations found in each trait, providing important insights into the biological mechanisms underlying E-cadherin’s pleiotropic effects. Ultimately, this knowledge will impact genetic counselling and will be relevant to the assessment of risk of cancer development or congenital malformations in CDH1 mutation carriers.

Genetic investigation into an increased susceptibility to biliary atresia in an extended New Zealand Māori family

Background

Biliary atresia (BA), a fibrosing disorder of the developing biliary tract leading to liver failure in infancy, has an elevated incidence in indigenous New Zealand (NZ) Māori. We investigated a high rate of BA in a group of children (n = 12) belonging to a single Māori iwi (or ‘tribe’, related through a remote ancestor).

Methods

Population and geographical data was used to estimate the rate of BA in Māori sub-groups, and a pedigree linking most of the affected children was constructed from oral and documented history. Array genotyping was used to examine hypotheses about the inheritance of a possible genetic risk factor, and the history of the affected population, and Exome Sequencing to search for candidate genes.

Results

Most of these affected children (n = 7) link to a self-reported pedigree and carry a 50-fold increase in BA risk over unrelated Māori (χ² = 296P < 0.001, 95% CI 23–111). Genetic analysis using FEstim and SNP array genotypes revealed no evidence for elevated consanguinity between parents of affected children (FEstim: F (2,21) = 0.469, P > 0.63). Genome-wide quantitation of intervals of contiguous, homozygous-by-state markers reached a similar conclusion (F (2,399) = 1.99, P = 0.138). Principal component analysis and investigation with STRUCTURE found no evidence of increased allele frequency of either a recessive variant, or additive, low-risk variants due to reproductive isolation. To identify candidate causal factors, Exome Sequencing datasets were scrutinised for shared rare coding variants across 8 affected individuals. No rare, non-synonymous, phylogenetically conserved variants were common to 6 or more affected children.

Conclusion

The substantially elevated risk for development of BA in this subgroup could be mediated by genetic factors, but the iwi exhibits no properties indicative of recent or remote reproductive isolation. Resolution of any risk loci may rely on extensive genomic sequencing studies in this iwi or investigation of other mechnaisms such as copy number variation.

Electronic supplementary material

The online version of this article (10.1186/s12920-018-0440-0) contains supplementary material, which is available to authorized users.

Simulating pedigrees ascertained for multiple disease-affected relatives

Background

Studies that ascertain families containing multiple relatives affected by disease can be useful for identification of causal, rare variants from next-generation sequencing data.

Results

We present the R package SimRVPedigree, which allows researchers to simulate pedigrees ascertained on the basis of multiple, affected relatives. By incorporating the ascertainment process in the simulation, SimRVPedigree allows researchers to better understand the within-family patterns of relationship amongst affected individuals and ages of disease onset.

Conclusions

Through simulation, we show that affected members of a family segregating a rare disease variant tend to be more numerous and cluster in relationships more closely than those for sporadic disease. We also show that the family ascertainment process can lead to apparent anticipation in the age of onset. Finally, we use simulation to gain insight into the limit on the proportion of ascertained families segregating a causal variant. SimRVPedigree should be useful to investigators seeking insight into the family-based study design through simulation.

Electronic supplementary material

The online version of this article (10.1186/s13029-018-0069-6) contains supplementary material, which is available to authorized users.

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.

A novel non-coding RNA within an intron of CDH2 and association of its SNP with non-syndromic cleft lip and palate

BACKGROUND

Genome-wide linkage analysis and whole genome sequencing in a Van der Woude syndrome (VWS) family revealed that the SNP, rs539075, within intron 2 of the cadherin 2 gene (CDH2) co-segregated with the disease phenotype.

RESULTS

A study with nonsyndromic cleft lip with or without cleft palate (NSCL ± P) cases (N = 292) and controls (N = 287) established association of this SNP with NSCL ± P as a risk factor. RT-PCR based expression analysis of the SNP-harbouring region of intron 2 of CDH2 in the clefted lip and/or palate tissues of 16 patients revealed that the mutant allele expressed in all those individuals having it (hetero-/homozygous), whereas the wild type allele expressed in <50% of the samples in which it was present. The intronic transcript was also present in the prospective lip and palate region of 13.5 dpc mouse embryo, detected by RNA in situ hybridization and RT-PCR.

CONCLUSIONS

These results including the in silico, characterization of the ~200 nt-intronic transcript showed that conformationally it fits best with noncoding small RNA, possibly a precursor of miRNA. Its function in the orofacial organogenesis remains to be elucidated which will enable us to define the role of this mutant ncRNA in the clefting of lip and palate.

Variants in members of the cadherin-catenin complex, CDH1 and CTNND1, cause blepharocheilodontic syndrome

Blepharocheilodontic syndrome (BCDS) consists of lagophthalmia, ectropion of the lower eyelids, distichiasis, euryblepharon, cleft lip/palate and dental anomalies and has autosomal dominant inheritance with variable expression. We identified heterozygous variants in two genes of the cadherin-catenin complex, CDH1, encoding E-cadherin, and CTNND1, encoding p120 catenin delta1 in 15 of 17 BCDS index patients, as was recently described in a different publication. CDH1 plays an essential role in epithelial cell adherence; CTNND1 binds to CDH1 and controls the stability of the complex. Functional experiments in zebrafish and human cells showed that the CDH1 variants impair the cell adhesion function of the cadherin-catenin complex in a dominant-negative manner. Variants in CDH1 have been linked to familial hereditary diffuse gastric cancer and invasive lobular breast cancer; however, no cases of gastric or breast cancer have been reported in our BCDS cases. Functional experiments reported here indicated the BCDS variants comprise a distinct class of CDH1 variants. Altogether, we identified the genetic cause of BCDS enabling DNA diagnostics and counseling, in addition we describe a novel class of dominant negative CDH1 variants.

Analysis of sequence data to identify potential risk variants for oral clefts in multiplex families

BACKGROUND

Nonsyndromic oral clefts are craniofacial malformations, which include cleft lip with or without cleft palate. The etiology for oral clefts is complex with both genetic and environmental factors contributing to risk. Previous genome-wide association (GWAS) studies have identified multiple loci with small effects; however, many causal variants remain elusive.

METHODS

In this study, we address this by specifically looking for rare, potentially damaging variants in family-based data. We analyzed both whole exome sequence (WES) data and whole genome sequence (WGS) data in multiplex cleft families to identify variants shared by affected individuals.

RESULTS

Here we present the results from these analyses. Our most interesting finding was from a single Syrian family, which showed enrichment of nonsynonymous and potentially damaging rare variants in two genes: CASP9 and FAT4.

CONCLUSION

Neither of these candidate genes has previously been associated with oral clefts and, if confirmed as contributing to disease risk, may indicate novel biological pathways in the genetic etiology for oral clefts.

Candidate Genes for Nonsyndromic Cleft Palate Detected by Exome Sequencing

Nonsyndromic cleft palate only (nsCPO) is a facial malformation that has a livebirth prevalence of 1 in 2,500. Research suggests that the etiology of nsCPO is multifactorial, with a clear genetic component. To date, genome-wide association studies have identified only 1 conclusive common variant for nsCPO, that is, a missense variant in the gene grainyhead-like-3 ( GRHL3). Thus, the underlying genetic causes of nsCPO remain largely unknown. The present study aimed at identifying rare variants that might contribute to nsCPO risk, via whole-exome sequencing (WES), in multiply affected Central European nsCPO pedigrees. WES was performed in 2 affected first-degree relatives from each family. Variants shared between both individuals were analyzed for their potential deleterious nature and a low frequency in the general population. Genes carrying promising variants were annotated for 1) reported associations with facial development, 2) multiple occurrence of variants, and 3) expression in mouse embryonic palatal shelves. This strategy resulted in the identification of a set of 26 candidate genes that were resequenced in 132 independent nsCPO cases and 623 independent controls of 2 different ethnicities, using molecular inversion probes. No rare loss-of-function mutation was identified in either WES or resequencing step. However, we identified 2 or more missense variants predicted to be deleterious in each of 3 genes ( ACACB, PTPRS, MIB1) in individuals from independent families. In addition, the analyses identified a novel variant in GRHL3 in 1 patient and a variant in CREBBP in 2 siblings. Both genes underlie different syndromic forms of CPO. A plausible hypothesis is that the apparently nonsyndromic clefts in these 3 patients might represent hypomorphic forms of the respective syndromes. In summary, the present study identified rare variants that might contribute to nsCPO risk and suggests candidate genes for further investigation.

Copy Number Changes Identified Using Whole Exome Sequencing in Nonsyndromic Cleft Lip and Palate in a Honduran Population

BACKGROUND

The majority of cleft lip with or without cleft palate cases appear as an isolated, nonsyndromic entity (NSCLP). With the advent of next generation sequencing, whole exome sequencing (WES) has been used to identify single nucleotide variants and insertion/deletions which cause or increase the risk of NSCLP. However, to our knowledge, there are no published studies using WES in NSCLP to investigate copy number changes (CNCs), which are a major component of human genetic variation. Our study aimed to identify CNCs associated with NSCLP in a Honduran population using WES.

METHODS

WES was performed on two to four members of 27 multiplex Honduran families. CNCs were identified using two algorithms, CoNIFER and XHMM. Priority was given to CNCs that were identified in more than one patient and had variant frequencies of less than 5% in reference data sets.

RESULTS

WES completion was defined as >90% of the WES target at >8 × coverage and >80% of the WES target at >20 × coverage. Twenty-four CNCs that met our inclusion criteria were identified by both CoNIFER and XHMM. These CNCs were confirmed using quantitative PCR. Pedigree analysis produced three CNCs corresponding to ADH7, AHR, and CRYZ segregating with NSCLP. Two of the three CNCs implicate genes, AHR and ADH7, whose known biological functions could plausibly play a role in NSCLP.

CONCLUSION

WES can be used to detect candidate CNCs that may be involved in the pathophysiology of NSCLP. Birth Defects Research 109:1257-1267, 2017. © 2017 Wiley Periodicals, Inc.

ACSS2 gene variant associated with cleft lip and palate in two independent Hispanic populations

OBJECTIVES/HYPOTHESIS

A candidate variant (p.Val496Ala) of the ACSS2 gene (T > C missense, rs59088485 variant at chr20: bp37 33509608) was previously found to consistently segregate with nonsyndromic cleft lip and/or palate (NSCLP) in three Honduran families. Objectives of this study were 1) to investigate the frequency of this ACSS2 variant in Honduran unrelated NSCLP patients and unrelated unaffected controls and 2) to investigate the frequency of this variant in Colombian unrelated affected NSCLP patients and unrelated unaffected controls.

STUDY DESIGN

Case-control studies.

METHODS

Sanger sequencing of 99 unrelated Honduran NSCLP patients and 215 unrelated unaffected controls for the p.Val496Ala ACSS2 variant was used to determine the carrier frequency in NSCLP patients and controls. Sanger sequencing of 230 unrelated Colombian NSCLP patients and 146 unrelated unaffected controls for the p.Val496Ala ACSS2 variant was used to determine the carrier frequency in NSCLP patients and controls.

RESULTS

In the Honduran population, the odds ratio of having NSCLP among carriers of the p.Val496Ala ACSS2 variant was 4.0 (P = .03), with a carrier frequency of seven of 99 (7.1%) in unrelated affected and four of 215 (1.9%) in unrelated unaffected individuals. In the Colombian population, the odds ratio of having NSCLP among carriers of the p.Val496Ala ACSS2 variant was 2.6 (P = .04), with a carrier frequency of 23 of 230 (10.0%) in unrelated affected and six of 146 (4.1%) in unrelated unaffected individuals.

CONCLUSIONS

These findings support the role of ACSS2 in NSCLP in two independent Hispanic populations from Honduras and Colombia.

LEVEL OF EVIDENCE

NA Laryngoscope, 127:E336-E339, 2017.

Intraflagellar transport 88 (IFT88) is crucial for craniofacial development in mice and is a candidate gene for human cleft lip and palate

Ciliopathies are pleiotropic human diseases resulting from defects of the primary cilium, and these patients often have cleft lip and palate. IFT88 is required for the assembly and function of the primary cilia, which mediate the activity of key developmental signaling pathways. Through whole exome sequencing of a family of three affected siblings with isolated cleft lip and palate, we discovered that they share a novel missense mutation in IFT88 (c.915G > C, p.E305D), suggesting this gene should be considered a candidate for isolated orofacial clefting. In order to evaluate the function of IFT88 in regulating craniofacial development, we generated Wnt1-Cre;Ift88fl/fl mice to eliminate Ift88 specifically in cranial neural crest (CNC) cells. Wnt1-Cre;Ift88fl/flpups died at birth due to severe craniofacial defects including bilateral cleft lip and palate and tongue agenesis, following the loss of the primary cilia in the CNC-derived palatal mesenchyme. Loss of Ift88 also resulted in a decrease in neural crest cell proliferation during early stages of palatogenesis as well as a downregulation of the Shh signaling pathway in the palatal mesenchyme. Importantly, Osr2KI-Cre;Ift88fl/flmice, in which Ift88 is lost specifically in the palatal mesenchyme, exhibit isolated cleft palate. Taken together, our results demonstrate that IFT88 has a highly conserved function within the primary cilia of the CNC-derived mesenchyme in the lip and palate region in mice and is a strong candidate as an orofacial clefting gene in humans.