Systematic analysis of copy number variants of a large cohort of orofacial cleft patients identifies candidate genes for orofacial clefts

Rare variants analyses suggest novel cleft genes in the African population

Non-syndromic orofacial clefts (NSOFCs) are common birth defects with a complex etiology. While over 60 common risk loci have been identified, they explain only a small proportion of the heritability for NSOFCs. Rare variants have been implicated in the missing heritability. Thus, our study aimed to identify genes enriched with nonsynonymous rare coding variants associated with NSOFCs. Our sample included 814 non-syndromic cleft lip with or without palate (NSCL/P), 205 non-syndromic cleft palate only (NSCPO), and 2150 unrelated control children from Nigeria, Ghana, and Ethiopia. We conducted a gene-based analysis separately for each phenotype using three rare-variants collapsing models: (1) protein-altering (PA), (2) missense variants only (MO); and (3) loss of function variants only (LOFO). Subsequently, we utilized relevant transcriptomics data to evaluate associated gene expression and examined their mutation constraint using the gnomeAD database. In total, 13 genes showed suggestive associations (p = E-04). Among them, eight genes (ABCB1, ALKBH8, CENPF, CSAD, EXPH5, PDZD8, SLC16A9, and TTC28) were consistently expressed in relevant mouse and human craniofacial tissues during the formation of the face, and three genes (ABCB1, TTC28, and PDZD8) showed statistically significant mutation constraint. These findings underscore the role of rare variants in identifying candidate genes for NSOFCs.

Rare Variants Analyses Suggest Novel Cleft Genes in the African Population

Non-syndromic orofacial clefts (NSOFCs) are common birth defects with a complex etiology. While over 60 common risk loci have been identified, they explain only a small proportion of the heritability for NSOFC. Rare variants have been implicated in the missing heritability. Thus, our study aimed to identify genes enriched with nonsynonymous rare coding variants associated with NSOFCs. Our sample included 814 non-syndromic cleft lip with or without palate (NSCL/P), 205 non-syndromic cleft palate only (NSCPO), and 2150 unrelated control children from Nigeria, Ghana, and Ethiopia. We conducted a gene-based analysis separately for each phenotype using three rare-variants collapsing models: (1) protein-altering (PA), (2) missense variants only (MO); and (3) loss of function variants only (LOFO). Subsequently, we utilized relevant transcriptomics data to evaluate associated gene expression and examined their mutation constraint using the gnomeAD database. In total, 13 genes showed suggestive associations (p = E-04). Among them, eight genes (ABCB1, ALKBH8, CENPF, CSAD, EXPH5, PDZD8, SLC16A9, and TTC28) were consistently expressed in relevant mouse and human craniofacial tissues during the formation of the face, and three genes (ABCB1, TTC28, and PDZD8) showed statistically significant mutation constraint. These findings underscore the role of rare variants in identifying candidate genes for NSOFCs.

Identification of a Novel Variant of PDGFC Associated with Nonsyndromic Cleft Lip and Palate in a Chinese Family

Nonsyndromic cleft lip with or without cleft palate (NSCL/P) accounts for 70% of the total number of patients with cleft lip with or without cleft palate (CL/P) and is the most common type of congenital deformity of the craniomaxillofacial region. In this study, whole exome sequencing (WES) and Sanger sequencing were performed on affected members of a Han Chinese family, and a missense variant in the platelet-derived growth factor C (PDGFC) gene (NM_016205: c.G93T: p.Q31H) was identified to be associated with NSCL/P. Bioinformatic studies demonstrated that the amino acid corresponding to this variation is highly conserved in many mammals and leads to a glutamine-to-histidine substitution in an evolutionarily conserved DNA-binding domain. It was found that the expression of PDGFC was significantly decreased in the dental pulp stem cells (DPSCs) of NSCL/P cases, compared to the controls, and that the variant (NM_016205: c.G93T) reduced the expression of PDGFC. In addition, the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis showed that Pdgfc deficiency disrupted NSCL/P-related signaling pathways such as the MAPK signaling pathway and cell adhesion molecules. In conclusion, our study identified a missense variant (NM_016205: c.G93T) in exon 1 of PDGFC potentially associated with susceptibility to NSCL/P.

Advances in research of biological functions of Isthmin-1

Isthmin-1 (ISM1) was initially thought to be a brain secretory factor, but with the development of technical means of research and the refinement of animal models, numerous studies have shown that this molecule is expressed in multiple tissues, suggesting that it may have multiple biological functions. As a factor that regulates growth and development, ISM1 is expressed in different animals with spatial and temporal variability and can coordinate the normal development of multiple organs. Recent studies have found that under the dependence of a non-insulin pathway, ISM1 can lower blood glucose, inhibit insulin-regulated lipid synthesis, promote protein synthesis, and affect the body's glucolipid and protein metabolism. In addition, ISM1 plays an important role in cancer development by promoting apoptosis and anti-angiogenesis, and by regulating multiple inflammatory pathways to influence the body's immune response. The purpose of this paper is to summarize relevant research results from recent years and to describe the key features of the biological functions of ISM1. We aimed to provide a theoretical basis for the study of ISM1 related diseases, and potential therapeutic strategies. The main biological functions of ISM1. Current studies on the biological functions of ISM1 focus on growth and development, metabolism, and anticancer treatment. During embryonic development, ISM1 is dynamically expressed in the zebrafish, African clawed frog, chick, mouse, and human, is associated with craniofacial malformations, abnormal heart localization, and hematopoietic dysfunction. ISM1 plays an important role in regulating glucose metabolism, lipid metabolism, and protein metabolism in the body. ISM1 affects cancer development by regulating cellular autophagy, angiogenesis, and the immune microenvironment.

Café-au-lait Spots and Cleft Palate: Not a Chance Association

The recognition of syndromic forms of cleft palate is important for condition-specific management. Here, we report a patient with cleft palate, congenital heart disease, intellectual disability, and café-au-lait spots who had a deletion of chromosome 15q14. The identification of the precise breakpoints using a Nanopore-based long-read sequencer showed that the deletion spanned MEIS2 and SPRED1 loci. Cleft palate and café-au-lait spots can be ascribed to MEIS2 and SPRED1, respectively. Patients with cleft palate and café-au-lait spots should be encouraged to undergo a detailed genomic evaluation, including screening for a 15q14 deletion, to enable appropriate anticipatory medico-surgical management and genetic counseling.

Genetics and orofacial clefts: a clinical perspective

Orofacial clefts (OFCs) are the most common congenital craniofacial anomaly seen in humans. Most OFCs are sporadic and isolated - these are thought to be multifactorial in origin. Chromosomal and monogenic variants account for the syndromic forms and for some of the non-syndromic inherited forms. This review discusses the importance of genetic testing and the current clinical strategy to deliver a genomics service that is of direct benefit to patients and their families.

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.

Alternative polyadenylation transcriptome-wide association study identifies APA-linked susceptibility genes in brain disorders

Alternative polyadenylation (APA) plays an essential role in brain development; however, current transcriptome-wide association studies (TWAS) largely overlook APA in nominating susceptibility genes. Here, we performed a 3' untranslated region (3'UTR) APA TWAS (3'aTWAS) for 11 brain disorders by combining their genome-wide association studies data with 17,300 RNA-seq samples across 2,937 individuals. We identified 354 3'aTWAS-significant genes, including known APA-linked risk genes, such as SNCA in Parkinson's disease. Among these 354 genes, ~57% are not significant in traditional expression- and splicing-TWAS studies, since APA may regulate the translation, localization and protein-protein interaction of the target genes independent of mRNA level expression or splicing. Furthermore, we discovered ATXN3 as a 3'aTWAS-significant gene for amyotrophic lateral sclerosis, and its modulation substantially impacted pathological hallmarks of amyotrophic lateral sclerosis in vitro. Together, 3'aTWAS is a powerful strategy to nominate important APA-linked brain disorder susceptibility genes, most of which are largely overlooked by conventional expression and splicing analyses.

Genome-wide analysis of copy-number variation in humans with cleft lip and/or cleft palate identifies COBLL1, RIC1, and ARHGEF38 as clefting genes

Cleft lip with or without cleft palate (CL/P) is a common birth defect with a complex, heterogeneous etiology. It is well established that common and rare sequence variants contribute to the formation of CL/P, but the contribution of copy-number variants (CNVs) to cleft formation remains relatively understudied. To fill this knowledge gap, we conducted a large-scale comparative analysis of genome-wide CNV profiles of 869 individuals from the Philippines and 233 individuals of European ancestry with CL/P with three primary goals: first, to evaluate whether differences in CNV number, amount of genomic content, or amount of coding genomic content existed within clefting subtypes; second, to assess whether CNVs in our cohort overlapped with known Mendelian clefting loci; and third, to identify unestablished Mendelian clefting genes. Significant differences in CNVs across cleft types or in individuals with non-syndromic versus syndromic clefts were not observed; however, several CNVs in our cohort overlapped with known syndromic and non-syndromic Mendelian clefting loci. Moreover, employing a filtering strategy relying on population genetics data that rare variants are on the whole more deleterious than common variants, we identify several CNV-associated gene losses likely driving non-syndromic clefting phenotypes. By prioritizing genes deleted at a rare frequency across multiple individuals with clefts yet enriched in our cohort of individuals with clefts compared to control subjects, we identify COBLL1, RIC1, and ARHGEF38 as clefting genes. CRISPR-Cas9 mutagenesis of these genes in Xenopus laevis and Danio rerio yielded craniofacial dysmorphologies, including clefts analogous to those seen in human clefting disorders.

Low-level complex mosaic with multiple cell lines affecting the 18q21.31q21.32 region in a patient with de novo 18q terminal deletion

We describe a 5-year-old girl who was diagnosed at birth with 18q de novo homogeneous deletion at G-banding karyotype. Her clinical condition, characterized by hypotonia, psychomotor retardation, short stature, deafness secondary to bilateral atresia of the external auditory canals, was in agreement with the 18q deletion syndrome though presence of coloboma of a single eye only suggested a mosaic condition as an unusual sign. By combining multiple technologies including array-CGH, FISH, and WGS, we found that the terminal deletion 18q21.32q23 (21 Mb) was in segmental mosaicism of the proximal region 18q21.31q21.32 (2.7 Mb), which showed a variable number of copies: one, two, or three, in 7, 41 and 55% of the cells respectively. Breakpoint junction analysis demonstrated the presence of an inv-dup del (18q) with a disomic segment of 4.7 kb between the inverted and non-inverted copies of the duplicated region 18q21.31q21.32. From these results, we propose that all three types of abnormal chr18 (the inv-dup del and the two 18q terminal deletions of different sizes) arisen from breaks in a dicentric mirror chromosome 18q, either in more than one embryo cell or from subsequent breaking-fusion-bridge cycles. The duplication region was with identical polymorphisms as in all non-recurrent inv-dup del rearrangements though, in contrast with most of them, the 18q abnormality was of maternal origin. Taking into account that distal 18q deletions are not rarely associated with inv-dup del(18q) cell lines, and that the non-disjunction of chromosome 18 takes place especially at maternal meiosis II rather than meiosis I, multiple rescue events starting from trisomic zygotes could be considered alternative to the postmitotic ones. From the clinical point of view, our case, as well as those of del(18q) in mosaic with the dic(18q), shows that the final phenotype is the sum of the different cell lines that acted on embryonic development with signs typical of both the 18q deletion syndrome and trisomy 18. Asymmetrical malformations, such as coloboma of the iris only in the right eye, confirm the underlying mosaicism regardless of whether it is still detectable in the blood.

The Evaluation of FGFR1, FGFR2 and FOXO1 in Orofacial Cleft Tissue

Although cleft lip with or without cleft palate (CL/P) is one of the most common congenital anomalies worldwide, the morphopathogenesis of non-syndromic orofacial clefts is still unclear. Many candidate genes have been proposed to play a causal role; however, only a few have been confirmed, leaving many still to be assessed. Taking into account the significance of FGFR1, FGFR2 and FOXO1 in embryogenesis, the aim of this work was to detect and compare the three candidate genes in cleft-affected lip and palatine tissue. Ten soft tissue samples were taken during cheiloplasty and veloplasty. The signals of the candidate genes were visualized using chromogenic in situ hybridization and analyzed using a semi-quantitative method. No statistically important difference in the distribution of FGFR1, FGFR2 and FOXO1 between neither the patients’ lip and vomer mucosa nor the control group was observed. Statistically significant very strong and strong correlations were found between genes in the lip and palatine tissue. The expression of FGFR1, FGFR2 and FOXO1 in cleft-affected lip and palatine tissue seems to be highly individual. Numerous intercorrelations between the genes do not exclude their role in the possible complex morphopathogenesis of orofacial clefts.

MEIS2 (15q14) gene deletions in siblings with mild developmental phenotypes and bifid uvula: documentation of mosaicism in an unaffected parent

MEIS2 (Meis homeobox 2) encodes a homeobox protein in the three amino acid loop extension (TALE) family of highly conserved homeodomain-containing transcription regulators important for development. MEIS2 deletions/mutations have been associated with cleft lip/palate, dysmorphic facial features, cardiac defects, as well as intellectual disability at a variable severity. Here we report on one familial case that two affected siblings carry the same non-mosaic ~ 423 kb genomic deletion at 15q14 encompassing the entirety of CDIN1 and the last three exons (ex. 10, 11, 12) of the MEIS2 gene, while their unaffected father is mosaic for the same deletion in about 10% lymphocytes. Both siblings presented with mild developmental delay and bifid uvula, while no congenital cardiac abnormalities were identified. The elder sister also showed syncopal episodes and mild speech delay and the father had atrial septal defects. This is the first report showing multiple family members inherit a genomic deletion resulting in a MEIS2 partial truncation from a mosaic parent. Taken all together, this study has important implications for genetic counseling regarding recurrence risk and also points to the importance of offering MEIS2 gene tests covering both point mutations and microdeletions to individuals with milder bifid uvula and developmental delay.

Clinical Practice Guidelines on the Treatment of Patients with Cleft Lip, Alveolus, and Palate: An Executive Summary

Significant treatment variation exists in the Netherlands between teams treating patients with cleft lip, alveolus, and/or palate, resulting in a confusing and undesirable situation for patients, parents, and practitioners. Therefore, to optimize cleft care, clinical practice guidelines (CPGs) were developed. The aim of this report is to describe CPG development, share the main recommendations, and indicate knowledge gaps regarding cleft care. Together with patients and parents, a multidisciplinary working group of representatives from all relevant disciplines assisted by two experienced epidemiologists identified the topics to be addressed in the CPGs. Searching the Medline, Embase, and Cochrane Library databases identified 5157 articles, 60 of which remained after applying inclusion and exclusion criteria. We rated the quality of the evidence from moderate to very low. The working group formulated 71 recommendations regarding genetic testing, feeding, lip and palate closure, hearing, hypernasality, bone grafting, orthodontics, psychosocial guidance, dentistry, osteotomy versus distraction, and rhinoplasty. The final CPGs were obtained after review by all stakeholders and allow cleft teams to base their treatment on current knowledge. With high-quality evidence lacking, the need for additional high-quality studies has become apparent.

Identification of Copy Number Variation Among Nonsyndromic Cleft Lip and or Without Cleft Palate With Hypodontia: A Genome-Wide Association Study

Nonsyndromic cleft lip and or without cleft palate (NSCL/P) with the hypodontia is a common developmental abnormality in humans and animals. This study identified the genetic aberration involved in both NSCL/P and hypodontia pathogenesis. A cross-sectional study using genome-wide study copy number variation-targeted CytoScan 750K array carried out on salivary samples from 61 NSCL/P and 20 noncleft with and without hypodontia Malay subjects aged 7-13 years old. Copy number variations (CNVs) of SKI and fragile histidine triad (FHIT) were identified in NSCL/P and noncleft children using quantitative polymerase chain reaction (qPCR) as a validation analysis. Copy number calculated (CNC) for each gene determined with Applied Biosystems CopyCaller Software v2.0. The six significant CNVs included gains (12q14.3, 15q26.3, 1p36.32, and 1p36.33) and losses (3p14.2 and 4q13.2) in NSCL/P with hypodontia patients compared with the NSCL/P only. The genes located in these regions encoded LEMD3, IGF1R, TP73, SKI, FHIT, and UGT2β15. There were a significant gain and loss of both SKI and FHIT copy number in NSCL/P with hypodontia compared with the noncleft group (p < 0.05). The results supported that CNVs significantly furnish to the development of NSCL/P with hypodontia.

A novel MEIS2 mutation explains the complex phenotype in a boy with a typical NF1 microdeletion syndrome

Concurrence of distinct genetic conditions in the same patient is not rare. Several cases involving neurofibromatosis type 1 (NF1) have recently been reported, indicating the need for more extensive molecular analysis when phenotypic features cannot be explained by a single gene mutation. Here, we describe the clinical presentation of a boy with a typical NF1 microdeletion syndrome complicated by cleft palate and other dysmorphic features, hypoplasia of corpus callosum, and partial bicoronal craniosynostosis caused by a novel 2bp deletion in exon 2 of Meis homeobox 2 gene (MEIS2) inherited from the mildly affected father. This is only the second case of an inherited MEIS2 intragenic mutation reported to date. MEIS2 is known to be associated with cleft palate, intellectual disability, heart defects, and dysmorphic features. Our clinical report suggests that this gene may also have a role in cranial morphogenesis in humans, as previously observed in animal models.

Clinical application of chromosomal microarray analysis for fetuses with craniofacial malformations

Background

The potential correlations between chromosomal abnormalities and craniofacial malformations (CFMs) remain a challenge in prenatal diagnosis. This study aimed to evaluate 118 fetuses with CFMs by applying chromosomal microarray analysis (CMA) and G-banded chromosome analysis.

Results

Of the 118 cases in this study, 39.8% were isolated CFMs (47/118) whereas 60.2% were non-isolated CFMs (71/118). The detection rate of chromosomal abnormalities in non-isolated CFM fetuses was significantly higher than that in isolated CFM fetuses (26/71 vs. 7/47, p = 0.01). Compared to the 16 fetuses (16/104; 15.4%) with pathogenic chromosomal abnormalities detected by karyotype analysis, CMA identified a total of 33 fetuses (33/118; 28.0%) with clinically significant findings. These 33 fetuses included cases with aneuploidy abnormalities (14/118; 11.9%), microdeletion/microduplication syndromes (9/118; 7.6%), and other pathogenic copy number variations (CNVs) only (10/118; 8.5%).We further explored the CNV/phenotype correlation and found a series of clear or suspected dosage-sensitive CFM genes including TBX1, MAPK1, PCYT1A, DLG1, LHX1, SHH, SF3B4, FOXC1, ZIC2, CREBBP, SNRPB, and CSNK2A1.

Conclusion

These findings enrich our understanding of the potential causative CNVs and genes in CFMs. Identification of the genetic basis of CFMs contributes to our understanding of their pathogenesis and allows detailed genetic counselling.

Expanded Prader-Willi Syndrome due to an Unbalanced de novo Translocation t(14;15): Report and Review of the Literature

In the present study, we report a case of a female infant with a de novo unbalanced t(14;15) translocation resulting in a 14-Mb deletion of the 15q11.1q14 region. The deletion includes the 15q11.2q13 Prader-Willi syndrome (PWS) critical region, while no known deleted genes are found in the 14qter region. According to literature review, patients with similar or larger deletions in the 15q region exhibit an expanded phenotype of PWS with case-specific atypical features such as severe retardation, absence of speech, microcephaly, retrognathia, bifid uvula, ear malformations, and heart defects in addition to typical features of PWS. Our proband exhibited increased deep tendon reflexes, an atypical feature which is not reported in the reviewed literature. The severity of the phenotype is not directly associated with the size of the deletion; however, using a combination of methods, the identification of breakpoints and the deleted genes can be helpful for the prognostication in patients with atypical PWS deletions.

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.

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.

Craniofacial, orofacial and dental disorders: the role of the RAS/ERK pathway

Deviations from the precisely coordinated programme of human head development can lead to craniofacial and orofacial malformations often including a variety of dental abnormalities too. Although the aetiology is still unknown in many cases, during the last decades different intracellular signalling pathways have been genetically linked to specific disorders. Among these pathways, the RAS/extracellular signal-regulated kinase (ERK) signalling cascade is the focus of this review since it encompasses a large group of genes that when mutated cause some of the most common and severe developmental anomalies in humans. We present the components of the RAS/ERK pathway implicated in craniofacial and orodental disorders through a series of human and animal studies. We attempt to unravel the specific molecular targets downstream of ERK that act on particular cell types and regulate key steps in the associated developmental processes. Finally we point to ambiguities in our current knowledge that need to be clarified before RAS/ERK-targeting therapeutic approaches can be implemented.

Deletions and loss-of-function variants in TP63 associated with orofacial clefting

We aimed to identify novel deletions and variants of TP63 associated with orofacial clefting (OFC). Copy number variants were assessed in three OFC families using microarray analysis. Subsequently, we analyzed TP63 in a cohort of 1072 individuals affected with OFC and 706 population-based controls using molecular inversion probes (MIPs). We identified partial deletions of TP63 in individuals from three families affected with OFC. In the OFC cohort, we identified several TP63 variants predicting to cause loss-of-function alleles, including a frameshift variant c.569_576del (p.(Ala190Aspfs*5)) and a nonsense variant c.997C>T (p.(Gln333*)) that introduces a premature stop codon in the DNA-binding domain. In addition, we identified the first missense variants in the oligomerization domain c.1213G>A (p.(Val405Met)), which occurred in individuals with OFC. This variant was shown to abrogate oligomerization of mutant p63 protein into oligomeric complexes, and therefore likely represents a loss-of-function allele rather than a dominant-negative. All of these variants were inherited from an unaffected parent, suggesting reduced penetrance of such loss-of-function alleles. Our data indicate that loss-of-function alleles in TP63 can also give rise to OFC as the main phenotype. We have uncovered the dosage-dependent functions of p63, which were previously rejected.

Unmasking familial CPX by WES and identification of novel clinical signs

Mutations in the T-Box transcription factor gene TBX22 are found in X-linked Cleft Palate with or without Ankyloglossia syndrome (CPX syndrome). In addition to X-linked inheritance, ankyloglossia, present in the majority of CPX patients, is an important diagnostic marker, but it is frequently missed or unreported, as it is a "minor" feature. Other described anomalies include cleft lip, micro and/or hypodontia, and features of CHARGE syndrome. We conducted whole exome sequencing (WES) on 22 individuals from 17 "a priori" non-syndromic cleft lip and/or cleft palate (CL/P) families. We filtered the data for heterozygous pathogenic variants within a set of predefined candidate genes. Two canonical splice-site mutations were found in TBX22. Detailed re-phenotyping of the two probands and their families unravelled orofacial features previously not associated with the CPX phenotypic spectrum: choanal atresia, Pierre-Robin sequence, and overgrowths on the posterior edge of the hard palate, on each side of the palatal midline. This study emphasizes the importance of WES analysis in familial CLP cases, combined with deep (reverse) phenotyping in "a priori" non-syndromic clefts.

GWAS reveals loci associated with velopharyngeal dysfunction

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

Application of high-resolution array platform for genome-wide copy number variation analysis in patients with nonsyndromic cleft lip and palate

BACKGROUND

Although more than 14 loci may be involved in the development of nonsyndromic cleft lip and palate (NSCLP), the etiology has not been fully elucidated due to genetic and environmental risk factor interactions. Despite advances in identifying genes associated with the NSCLP development using traditional genetic mapping strategies of candidate genes, genome-wide studies, and epidemiologic and linkage analysis, microarray techniques have become important complementary tools in the search for potential causative oral clefts genes in genetic studies. Microarray hybridization enables scanning of the whole genome and detecting copy number variants (CNVs). Although common benign CNVs are often smaller, with sizes smaller than 20 kb, here we reveal small exonic CNVs based on the importance of the encompassed genes in cleft lip and palate phenotype.

METHODS

Microarray hybridization analysis was performed in 15 individuals with NSCLP.

RESULTS

We identified 11 exonic CNVs affecting at least one exon of the candidate genes. Thirteen candidate genes (COL11A1-1p21; IRF6-1q32.3; MSX1-4p16.2; TERT-5p15.33; MIR4457-5p15.33; CLPTM1L-5p15.33; ESR1-6q25.1; GLI3-7p13; FGFR-8p11.23; TBX1-22q11.21; OFD-Xp22; PHF8-Xp11.22; and FLNA-Xq28) overlapped with the CNVs identified.

CONCLUSIONS

Considering the importance to NSCLP, the microdeletions that encompass MSX1, microduplications over TERT, MIR4457, CLPTM1L, and microduplication of PHF8 have been identified as small CNVs related to sequence variants associated with oral clefts susceptibility. Our findings represent a preliminary study on the clinical significance of small CNVs and their relationship with genes implicated in NSCLP.

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

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

Identification of novel candidate disease genes from de novo exonic copy number variants

BACKGROUND

Exon-targeted microarrays can detect small (<1000 bp) intragenic copy number variants (CNVs), including those that affect only a single exon. This genome-wide high-sensitivity approach increases the molecular diagnosis for conditions with known disease-associated genes, enables better genotype-phenotype correlations, and facilitates variant allele detection allowing novel disease gene discovery.

METHODS

We retrospectively analyzed data from 63,127 patients referred for clinical chromosomal microarray analysis (CMA) at Baylor Genetics laboratories, including 46,755 individuals tested using exon-targeted arrays, from 2007 to 2017. Small CNVs harboring a single gene or two to five non-disease-associated genes were identified; the genes involved were evaluated for a potential disease association.

RESULTS

In this clinical population, among rare CNVs involving any single gene reported in 7200 patients (11%), we identified 145 de novo autosomal CNVs (117 losses and 28 intragenic gains), 257 X-linked deletion CNVs in males, and 1049 inherited autosomal CNVs (878 losses and 171 intragenic gains); 111 known disease genes were potentially disrupted by de novo autosomal or X-linked (in males) single-gene CNVs. Ninety-one genes, either recently proposed as candidate disease genes or not yet associated with diseases, were disrupted by 147 single-gene CNVs, including 37 de novo deletions and ten de novo intragenic duplications on autosomes and 100 X-linked CNVs in males. Clinical features in individuals with de novo or X-linked CNVs encompassing at most five genes (224 bp to 1.6 Mb in size) were compared to those in individuals with larger-sized deletions (up to 5 Mb in size) in the internal CMA database or loss-of-function single nucleotide variants (SNVs) detected by clinical or research whole-exome sequencing (WES). This enabled the identification of recently published genes (BPTF, NONO, PSMD12, TANGO2, and TRIP12), novel candidate disease genes (ARGLU1 and STK3), and further confirmation of disease association for two recently proposed disease genes (MEIS2 and PTCHD1). Notably, exon-targeted CMA detected several pathogenic single-exon CNVs missed by clinical WES analyses.

CONCLUSIONS

Together, these data document the efficacy of exon-targeted CMA for detection of genic and exonic CNVs, complementing and extending WES in clinical diagnostics, and the potential for discovery of novel disease genes by genome-wide assay.

Tubby domain superfamily protein is required for the formation of the 7S SNARE complex in Drosophila

Tubby domain superfamily protein (TUSP) is a distant member of the Tubby-like protein (TULP) family. Although other TULPs play important roles in sensation, metabolism, and development, the molecular functions of TUSP are completely unknown. Here, we explore the function of TUSP in the Drosophila nervous system where it is expressed in all neurons. Tusp mutant flies exhibit a temperature-sensitive paralysis. This paralysis can be rescued by tissue-specific expression of Tusp in the giant fibers and peripherally synapsing interneurons of the giant fiber system, a well-characterized neuronal circuit that mediates rapid escape behavior in flies. Consistent with this paralytic phenotype, we observed a profound reduction in the assembly of the ternary 7S SNARE complex that is required for neurotransmitter release despite seeing no changes in the expression of each individual SNARE complex component. Together, these data suggest TUSP is a novel regulator of SNARE assembly and, therefore, of neurotransmitter release.

Array comparative genomic hybridization and genomic sequencing in the diagnostics of the causes of congenital anomalies

The aim of this review is to provide the current state of knowledge about the usefulness of modern genetic technologies in uncovering the causality of isolated and multiple congenital anomalies. Array comparative genomic hybridization and next-generation sequencing have revolutionized the clinical approach to patients with these phenotypes. Both technologies enable early diagnosis, especially in clinically challenging newborn populations, and help to uncover genetic defects associated with various phenotypes. The application of both complementary methods could assist in identifying many variants that may simultaneously be involved in the development of a number of isolated or multiple congenital anomalies. Both technologies carry serious variant misinterpretation risks as well. Therefore, the methods of variant classification and accessible variant databases are mentioned. A useful strategy of clinical genetic testing with the application of both methodologies is presented. Finally, future directions and challenges are briefly commented on in this review.

SATB2-associated syndrome: Mechanisms, phenotype, and practical recommendations

The SATB2‐associated syndrome is a recently described syndrome characterized by developmental delay/intellectual disability with absent or limited speech development, craniofacial abnormalities, behavioral problems, dysmorphic features, and palatal and dental abnormalities. Alterations of the SATB2 gene can result from a variety of different mechanisms that include contiguous deletions, intragenic deletions and duplications, translocations with secondary gene disruption, and point mutations. The multisystemic nature of this syndrome demands a multisystemic approach and we propose evaluation and management guidelines. The SATB2‐associated syndrome registry has now been started and that will allow gathering further clinical information and refining the provided surveillance recommendations. © 2016 The Authors. American Journal of Medical Genetics Part A Published by Wiley Periodicals, Inc.