Irf6 directly regulates Klf17 in zebrafish periderm and Klf4 in murine oral epithelium, and dominant-negative KLF4 variants are present in patients with cleft lip and palate

The application of zebrafish model in the study of cleft lip and palate development: A systematic review

Objective

Craniofacial growth and development are more than a scientific curiosity; it is of tremendous interest to clinicians. Insights into the genetic etiology of cleft lip and palate development are essential for improving diagnosis and treatment planning. The purpose of this systematic review was to utilize a zebrafish model to highlight the role of the IRF6 gene in cleft lip and palate development in humans.

Data

This review adhered to the guidelines outlined in the PRISMA statement. Nine studies were included in the analysis.

Sources

This study used major scientific databases such as MEDLINE, EMBASE, Web of Science, and the Zebrafish Information Network and yielded 1275 articles. Two reviewers performed the screening using COVIDENCE™ independently, and a third reviewer resolved any conflicts.

Study selection

After applying the inclusion and exclusion criteria and screening, nine studies were included in the analysis. The Systematic Review Center for Laboratory Animal Experimentation's (SYRCLE's) risk-of-bias tool was used to assess the quality of the included studies.

Results

The main outcome supports the role of the IRF6 gene in zebrafish periderm development and embryogenesis, and IRF6 variations result in cleft lip and palate development. The overall SYRCLE risk of bias was low-medium.

Conclusion

In conclusion, this review indicated the critical role of the IRF6 gene and its downstream genes (GRHL3, KLF17, and ESRP1/2) in the development of cleft lip and palate in zebrafish models. Genetic mutation zebrafish models provide a high level of insights into zebrafish craniofacial development.

Clinical relevance

this review provides a productive avenue for understanding the powerful and conserved zebrafish model for investigating the pathogenesis of human cleft lip and palate.

Periderm Fate during Palatogenesis: TGF-β and Periderm Dedifferentiation

Failure of palatogenesis results in cleft palate, one of the most common congenital disabilities in humans. During the final phases of palatogenesis, the protective function of the peridermal cell layer must be eliminated for the medial edge epithelia to adhere properly, which is a prerequisite for the successful fusion of the secondary palate. However, a deeper understanding of the role and fate of the periderm in palatal adherence and fusion has been hampered due to a lack of appropriate periderm-specific genetic tools to examine this cell type in vivo. Here we used the cytokeratin-6A (Krt-6a) locus to develop both constitutive (Krt6ai-Cre) and inducible (Krt6ai-CreERT2) periderm-specific Cre driver mouse lines. These novel lines allowed us to achieve both the spatial and temporal control needed to dissect the periderm fate on a cellular resolution during palatogenesis. Our studies suggest that, already before the opposing palatal shelves contact each other, at least some palatal periderm cells start to gradually lose their squamous periderm-like phenotype and dedifferentiate into cuboidal cells, reminiscent of the basal epithelial cells seen in the palatal midline seam. Moreover, we show that transforming growth factor-β (TGF-β) signaling plays a critical periderm-specific role in palatogenesis. Thirty-three percent of embryos lacking a gene encoding the TGF-β type I receptor (Tgfbr1) in the periderm display a complete cleft of the secondary palate. Our subsequent experiments demonstrated that Tgfbr1-deficient periderm fails to undergo appropriate dedifferentiation. These studies define the periderm cell fate during palatogenesis and reveal a novel, critical role for TGF-β signaling in periderm dedifferentiation, which is a prerequisite for appropriate palatal epithelial adhesion and fusion.

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.

Consistent downregulation of the cleft lip/palate-associated genes IRF6 and GRHL3 in carcinomas

Interferon Regulatory Factor 6 (IRF6) and Grainyhead Like Transcription Factor 3 (GRHL3) are transcription factors that orchestrate gene regulatory networks required for the balance between keratinocyte differentiation and proliferation. Absence of either protein results in the lack of a normal stratified epidermis with keratinocytes failing to stop proliferating and to terminally differentiate. Numerous pathological variants within IRF6 and GRHL3 have been identified in orofacial cleft-affected individuals and expression of the two transcription factors has been found to be often dysregulated in cancers. However, whether orofacial cleft-associated IRF6 and GRHL3 variants in patients might also affect their cancer risk later in life, is not clear yet. The fact that the role of IRF6 and GRHL3 in cancer remains controversial makes this question even more challenging. Some studies identified IRF6 and GRHL3 as oncogenes, while others could attribute tumor suppressive functions to them. Trying to solve this apparent conundrum, we herein aimed to characterize IRF6 and GRHL3 function in various types of carcinomas. We screened multiple cancer and normal cell lines for their expression, and subsequently proceeded with functional assays in cancer cell lines. Our data uncovered consistent downregulation of IRF6 and GRHL3 in all types of carcinomas analyzed. Reduced levels of IRF6 and GRHL3 were found to be associated with several tumorigenic properties, such as enhanced cell proliferation, epithelial mesenchymal transition, migration and reduced differentiation capacity. Based on our findings, IRF6 and GRHL3 can be considered as tumor suppressor genes in various carcinomas, which makes them potential common etiological factors for cancer and CLP in a fraction of CLP-affected patients.

Zebrafish as model for studies in dentistry

INTRODUCTION

Over the last years, zebrafish has gained prominence in the biomedical community. It is currently considered one of the best vertebrate animal models for various types of studies, such as toxicology and developmental biology.

OBJECTIVE

The aim of this study was to conduct a literature review on the use of zebrafish in dentistry and whether this animal model could be a viable alternative for performing different types of studies in this area.

METHODS

A literature search was performed using the PubMed, Lilacs, Embase, and Dentistry and Oral Sciences Source. The keywords used as search terms were zebrafish and dentistry. The selection criteria were articles published in English that used zebrafish as an animal model in dentistry, oral health, and craniofacial growth/development.

RESULTS

The electronic search of literature yielded 421 articles. After the analysis of the abstracts, 29 articles were selected for an in-depth analysis and reading of the full text.

CONCLUSIONS

All studies included in this review confirm zebrafish's excellence as an animal model for various types of dentistry studies, as well as assisting and complementing other studies involving mammals.

Revisiting the embryogenesis of lip and palate development

Clefts of the lip and palate (CLP), the major causes of congenital facial malformation globally, result from failure of fusion of the facial processes during embryogenesis. With a prevalence of 1 in 500-2500 live births, CLP causes major morbidity throughout life as a result of problems with facial appearance, feeding, speaking, obstructive apnoea, hearing and social adjustment and requires complex, multi-disciplinary care at considerable cost to healthcare systems worldwide. Long-term outcomes for affected individuals include increased mortality compared with their unaffected siblings. The frequent occurrence and major healthcare burden imposed by CLP highlight the importance of dissecting the molecular mechanisms driving facial development. Identification of the genetic mutations underlying syndromic forms of CLP, where CLP occurs in association with non-cleft clinical features, allied to developmental studies using appropriate animal models is central to our understanding of the molecular events underlying development of the lip and palate and, ultimately, how these are disturbed in CLP.

Application of zebrafish in the study of craniomaxillofacial developmental anomalies

Craniomaxillofacial developmental anomalies are one of the most prevalent congenital defects worldwide and could result from any disruption of normal development processes, which is generally influenced by interactions between genes and the environment. Currently, with the advances in genetic screening strategies, an increasing number of novel variants and their roles in orofacial diseases have been explored. Zebrafish is recognized as a powerful animal model, and its homologous genes and similar oral structure and development process provide an ideal platform for studying the contributions of genetic and environmental factors to human craniofacial malformations. Here, we reviewed zebrafish models for the study of craniomaxillofacial developmental anomalies, such as human nonsyndromic cleft lip with or without an affected palate and jaw and tooth developmental anomalies. Due to its potential for gene expression and regulation research, zebrafish may provide new perspectives for understanding craniomaxillofacial diseaseand its treatment.

Krüppel-like factors in bone biology

The krüppel-like factor (KLF) family is a group of zinc finger transcription factors and contributes to different cellular processes such as differentiation, proliferation, migration, and apoptosis. While different studies show the roles of this family in skeletal development-specifically in chondrocyte and osteocyte development and bone homeostasis-there are few reviews summarizing their importance. To fill this gap, this review discusses current knowledge on different functions of the KLF family during skeletal development, including their roles in stem cell maintenance and differentiation, cell apoptosis, and cell cycle. To understand the importance of the KLF family, we also review genotype-phenotype correlations in different animal models. We also discuss how KLF proteins function through different signaling pathways and display their paramount importance in skeletal development. To highlight their roles in cartilage- or bone-related cells, we also use single-cell RNA sequencing publicly available data on mouse hindlimb. We also challenge our knowledge of how the KLF family is epigenetically regulated-e.g., using DNA methylation, histone modifications, and noncoding RNAs-during chondrocyte and osteocyte development.

NELL-1 in Genome-Wide Association Studies across Human Diseases

Neural epidermal growth factor-like (EGFL)-like protein (NELL)-1 is a potent and key osteogenic factor in the development and regeneration of skeletal tissues. Intriguingly, accumulative data from genome-wide association studies (GWASs) have started unveiling potential broader roles of NELL-1 beyond its functions in bone and cartilage. With exploration of the genetic variants of the entire genome in large-scale disease cohorts, GWASs have been used for establishing the connection between specific single-nucleotide polymorphisms of NELL1, in addition to osteoporosis, metabolic diseases, inflammatory conditions, neuropsychiatric diseases, neurodegenerative disorders, and malignant tumors. This review summarizes the findings from GWASs on the manifestation, significance level, implications on function, and correlation of specific NELL1 single-nucleotide polymorphisms in various disorders in humans. By offering a unique and comprehensive correlation between genetic variants and plausible functions of NELL1 in GWASs, this review illustrates the wide range of potential effects of a single gene on the pathogenesis of multiple disorders in humans.

Two rare variants reveal the significance of Grainyhead‐like 3 Arginine 391 underlying non‐syndromic cleft palate only

OBJECTIVES

Non-syndromic cleft palate only (NSCPO) is one of the most common craniofacial birth defects with largely undetermined genetic etiology. It has been established that Grainyhead-like 3 (GRHL3) plays an essential role in the pathogenesis of NSCPO. This study aimed to identify and verify the first-reported GRHL3 variant underlying NSCPO among the Chinese cohort.

METHODS

We performed whole-exome sequencing (WES) on a Chinese NSCPO patient and identified a rare variant of GRHL3 (p.Arg391His). A validated deleterious variant p.Arg391Cys was introduced as a positive control. Zebrafish embryos injection, reporter assays, live-cell imaging, and RNA sequencing were conducted to test the pathogenicity of the variants.

RESULTS

Zebrafish embryos microinjection demonstrated that overexpression of the variants could disrupt the normal development of zebrafish embryos. Reporter assays showed that Arg391His disturbed transcriptional activity of GRHL3 and exerted a dominant-negative effect. Interestingly, Arg391His and Arg391Cys displayed distinct nuclear localization patterns from that of wild-type GRHL3 in live-cell imaging. Bulk RNA sequencing suggested that the two variants changed the pattern of gene expression.

CONCLUSIONS

In aggregate, this study identified and characterized a rare GRHL3 variant in NSCPO, revealing the critical role of Arginine 391 in GRHL3. Our findings will help facilitate understanding and genetic counseling of NSCPO.

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.

Aging disrupts circadian gene regulation and function in macrophages

Aging is characterized by an increased vulnerability to infection and the development of inflammatory diseases, such as atherosclerosis, frailty, cancer and neurodegeneration. Here, we find that aging is associated with the loss of diurnally rhythmic innate immune responses, including monocyte trafficking from bone marrow to blood, response to lipopolysaccharide and phagocytosis. This decline in homeostatic immune responses was associated with a striking disappearance of circadian gene transcription in aged compared to young tissue macrophages. Chromatin accessibility was significantly greater in young macrophages than in aged macrophages; however, this difference did not explain the loss of rhythmic gene transcription in aged macrophages. Rather, diurnal expression of Kruppel-like factor 4 (Klf4), a transcription factor (TF) well established in regulating cell differentiation and reprogramming, was selectively diminished in aged macrophages. Ablation of Klf4 expression abolished diurnal rhythms in phagocytic activity, recapitulating the effect of aging on macrophage phagocytosis. Examination of individuals harboring genetic variants of KLF4 revealed an association with age-dependent susceptibility to death caused by bacterial infection. Our results indicate that loss of rhythmic Klf4 expression in aged macrophages is associated with disruption of circadian innate immune homeostasis, a mechanism that may underlie age-associated loss of protective immune responses.

A Novel PAX3 Variant in a Chinese Pedigree with Nonsyndromic Cleft Lip With or Without Palate

Objectives: Nonsyndromic cleft lip with or without palate (NSCL/P) is a common congenital orofacial defect, which is associated with severe disruption of orofacial development. The present study was designed to identify potential underlying gene variants in a Chinese pedigree with NSCL/P, in which the proband and the proband's father were affected. Methods: DNA was extracted from the participants' peripheral venous blood, and whole-exome sequencing was performed on the proband and the proband's parents. Results: After filtering, a paired box gene 3 (PAX3) missense variant (c.92C>G_p.Thr31Ser) was identified, which was verified by Sanger sequencing. This variant, which was not present in 113 unrelated healthy individuals or in a Chinese public database, may affect the transcription inhibition domain of the PAX3 protein. Conservation analysis and in silico predictions suggested that this variant may be evolutionarily conserved and potentially deleterious. In addition, it was reported that mice with PAX3 variants show cleft palates. Thus, the PAX3 missense variant (c.92C>G_p.Thr31Ser) is a candidate causative variant in this family. Conclusions: To the best of our knowledge, the present study is the first to report on a PAX3 variant in a pedigree with NSCL/P. The present study further suggests that PAX3 may be associated with CL/P etiology.

Lack of IRF6 Disrupts Human Epithelial Homeostasis by Altering Colony Morphology, Migration Pattern, and Differentiation Potential of Keratinocytes

Variants within the gene encoding for the transcription factor Interferon Regulatory Factor 6 (IRF6) are associated with syndromic and non-syndromic Cleft Lip/Palate (CLP) cases. IRF6 plays a vital role in the regulation of the proliferation/differentiation balance in keratinocytes and is involved in wound healing and migration. Since a fraction of CLP patients undergoing corrective cleft surgery experience wound healing complications, IRF6 represents an interesting candidate gene linking the two processes. However, Irf6 function has been mainly studied in mice and knowledge on IRF6 in human cells remains sparse. Here, we aimed to elucidate the role of IRF6 in human postnatal skin- and oral mucosa-derived keratinocytes. To do so, we applied CRISPR/Cas9 to ablate IRF6 in two TERT-immortalized keratinocyte cultures, which we used as model cell lines. We show that IRF6 controls the appearance of single cells and colonies, with the latter being less cohesive in its absence. Consequently, IRF6 knockout keratinocytes often moved as single cells instead of a collective epithelial sheet migration but maintained their epithelial character. Lack of IRF6 triggered severe keratinocyte differentiation defects, which were already apparent in the stratum spinosum and extended to the stratum corneum in 3D organotypic skin cultures, while it did not alter their growth rate. Finally, proteomics revealed that most of the differentially expressed proteins in the absence of IRF6 could be associated with differentiation, cell-cell adhesion as well as immune response. Our data expand the knowledge on IRF6 in human postnatal keratinocytes, which will help to better understand IRF6-related pathologies.

An Irf6-Esrp1/2 regulatory axis controls midface morphogenesis in vertebrates

Irf6 and Esrp1 are important for palate development across vertebrates. In zebrafish, we found that irf6 regulates the expression of esrp1 We detailed overlapping Irf6 and Esrp1/2 expression in mouse orofacial epithelium. In zebrafish, irf6 and esrp1/2 share expression in periderm, frontonasal ectoderm and oral epithelium. Genetic disruption of irf6 and esrp1/2 in zebrafish resulted in cleft of the anterior neurocranium. The esrp1/2 mutant also developed cleft of the mouth opening. Lineage tracing of cranial neural crest cells revealed that the cleft resulted not from migration defect, but from impaired chondrogenesis. Analysis of aberrant cells within the cleft revealed expression of sox10, col1a1 and irf6, and these cells were adjacent to krt4+ and krt5+ cells. Breeding of mouse Irf6; Esrp1; Esrp2 compound mutants suggested genetic interaction, as the triple homozygote and the Irf6; Esrp1 double homozygote were not observed. Further, Irf6 heterozygosity reduced Esrp1/2 cleft severity. These studies highlight the complementary analysis of Irf6 and Esrp1/2 in mouse and zebrafish, and identify a unique aberrant cell population in zebrafish expressing sox10, col1a1 and irf6 Future work characterizing this cell population will yield additional insight into cleft pathogenesis.

Cellular and molecular mechanisms in the development of a cleft lip and/or cleft palate; insights from zebrafish (Danio rerio)

Human cleft lip and/or palate (CLP) are immediately recognizable congenital abnormalities of the face. Lip and palate develop from facial primordia through the coordinated activities of ectodermal epithelium and neural crest cells (NCCs) derived from ectomesenchyme tissue. Subtle changes in the regulatory mechanisms of NCC or ectodermal epithelial cells can result in CLP. Genetic and environmental contributions or a combination of both play a significant role in the progression of CLP. Model organisms provide us with a wealth of information in understanding the pathophysiology and genetic etiology of this complex disease. Small teleost, zebrafish (Danio rerio) is one of the popular model in craniofacial developmental biology. The short generation time and large number of optically transparent, easily manipulated embryos increase the value of zebrafish to identify novel candidate genes and gene regulatory networks underlying craniofacial development. In addition, it is widely used to identify the mechanisms of environmental teratogens and in therapeutic drug screening. Here, we discuss the value of zebrafish as a model to understand epithelial and NCC induced ectomesenchymal cell activities during early palate morphogenesis and robustness of the zebrafish in modern research on identifying the genetic and environmental etiological factors of CLP.

A Novel Van der Woude Syndrome-Causing IRF6 Variant Is Subject to Incomplete Non-sense-Mediated mRNA Decay Affecting the Phenotype of Keratinocytes

Van der Woude syndrome (VWS) is a genetic syndrome that leads to typical phenotypic traits, including lower lip pits and cleft lip/palate (CLP). The majority of VWS-affected patients harbor a pathogenic variant in the gene encoding for the transcription factor interferon regulatory factor 6 (IRF6), a crucial regulator of orofacial development, epidermal differentiation and tissue repair. However, most of the underlying mechanisms leading from pathogenic IRF6 gene variants to phenotypes observed in VWS remain poorly understood and elusive. The availability of one VWS individual within our cohort of CLP patients allowed us to identify a novel VWS-causing IRF6 variant and to functionally characterize it. Using VWS patient-derived keratinocytes, we reveal that most of the mutated IRF6_VWS transcripts are subject to a non-sense-mediated mRNA decay mechanism, resulting in IRF6 haploinsufficiency. While moderate levels of IRF6_VWS remain detectable in the VWS keratinocytes, our data illustrate that the IRF6_VWS protein, which lacks part of its protein-binding domain and its whole C-terminus, is noticeably less stable than its wild-type counterpart. Still, it maintains transcription factor function. As we report and characterize a so far undescribed VWS-causing IRF6 variant, our results shed light on the physiological as well as pathological role of IRF6 in keratinocytes. This acquired knowledge is essential for a better understanding of the molecular mechanisms leading to VWS and CLP.

Inactivation of Zeb1 in GRHL2-deficient mouse embryos rescues mid-gestation viability and secondary palate closure

Cleft lip and palate are common birth defects resulting from failure of the facial processes to fuse during development. The mammalian grainyhead-like (Grhl1-3) genes play key roles in a number of tissue fusion processes including neurulation, epidermal wound healing and eyelid fusion. One family member, Grhl2, is expressed in the epithelial lining of the first pharyngeal arch in mice at embryonic day (E)10.5, prompting analysis of the role of this factor in palatogenesis. Grhl2-null mice die at E11.5 with neural tube defects and a cleft face phenotype, precluding analysis of palatal fusion at a later stage of development. However, in the first pharyngeal arch of Grhl2-null embryos, dysregulation of transcription factors that drive epithelial-mesenchymal transition (EMT) occurs. The aberrant expression of these genes is associated with a shift in RNA-splicing patterns that favours the generation of mesenchymal isoforms of numerous regulators. Driving the EMT perturbation is loss of expression of the EMT-suppressing transcription factors Ovol1 and Ovol2, which are direct GRHL2 targets. The expression of the miR-200 family of microRNAs, also GRHL2 targets, is similarly reduced, resulting in a 56-fold upregulation of Zeb1 expression, a major driver of mesenchymal cellular identity. The critical role of GRHL2 in mediating cleft palate in Zeb1^−/− mice is evident, with rescue of both palatal and facial fusion seen in Grhl2^−/−;Zeb1^−/− embryos. These findings highlight the delicate balance between GRHL2/ZEB1 and epithelial/mesenchymal cellular identity that is essential for normal closure of the palate and face. Perturbation of this pathway may underlie cleft palate in some patients.

Summary: Epithelial transcription factor GRHL2 is required for face closure while mesenchymal transcription factor ZEB1 is required for palate closure. Surprisingly, animals lacking both factors close their face and secondary palate.

Analysis of zebrafish periderm enhancers facilitates identification of a regulatory variant near human KRT8/18

Genome-wide association studies for non-syndromic orofacial clefting (OFC) have identified single nucleotide polymorphisms (SNPs) at loci where the presumed risk-relevant gene is expressed in oral periderm. The functional subsets of such SNPs are difficult to predict because the sequence underpinnings of periderm enhancers are unknown. We applied ATAC-seq to models of human palate periderm, including zebrafish periderm, mouse embryonic palate epithelia, and a human oral epithelium cell line, and to complementary mesenchymal cell types. We identified sets of enhancers specific to the epithelial cells and trained gapped-kmer support-vector-machine classifiers on these sets. We used the classifiers to predict the effects of 14 OFC-associated SNPs at 12q13 near KRT18. All the classifiers picked the same SNP as having the strongest effect, but the significance was highest with the classifier trained on zebrafish periderm. Reporter and deletion analyses support this SNP as lying within a periderm enhancer regulating KRT18/KRT8 expression.

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.

Exploring GRHL3 polymorphisms and SNP-SNP interactions in the risk of non-syndromic oral clefts in the Brazilian population

OBJECTIVE

To investigate the association of single-nucleotide polymorphisms (SNP) in grainyhead-like 3 (GRHL3) and to verify its possible interactions with others genes responsible for craniofacial development in the risk of non-syndromic oral cleft (NSOC).

METHODS

Applying TaqMan allelic discrimination assays, we evaluated GRHL3 SNPs (rs10903078, rs41268753, and rs4648975) in an ancestry-structured case-control sample composed of 1,127 Brazilian participants [272 non-syndromic cleft palate only (NSCPO), 242 non-syndromic cleft lip only (NSCLO), 319 non-syndromic cleft lip and palate (NSCLP), and 294 healthy controls]. Additionally, SNP-SNP interactions of GRHL3 and previously reported variants in FAM49A, FOXE1, NTN1, and VAX1 were verified in non-syndromic cleft lip with or without cleft palate (NSCL ± P). To eliminate false-positive associations, Bonferroni correction or 1,000 permutation method was applied.

RESULTS

The multiple logistic regression analysis showed that the CC genotype of rs10903078 (p = .03) and the haplotype C-C formed by the SNPs rs10903078 and rs41268753 (p = .04) were associated with NSCLO, but the p-values did not withstand Bonferroni correction. However, SNP-SNP test revealed significant interactions between GRHL3 SNPs and FAM49A (rs7552), FOXE1 (rs3758249), VAX1 (rs7078160 and rs751231), and NTN1 (rs9891446).

CONCLUSIONS

Our results confirm the importance of GRHL3 and its interactions with previously NSOC-associated genes, including FAM49A, FOXE1, NTN1, and VAX1, in the pathogenesis of NSOC in the Brazilian population.

Tissue-Specific Transcriptomes Reveal Gene Expression Trajectories in Two Maturing Skin Epithelial Layers in Zebrafish Embryos

Epithelial cells are the building blocks of many organs, including skin. The vertebrate skin initially consists of two epithelial layers, the outer periderm and inner basal cell layers, which have distinct properties, functions, and fates. The embryonic periderm ultimately disappears during development, whereas basal cells proliferate to form the mature, stratified epidermis. Although much is known about mechanisms of homeostasis in mature skin, relatively little is known about the two cell types in pre-stratification skin. To define the similarities and distinctions between periderm and basal skin epithelial cells, we purified them from zebrafish at early development stages and deeply profiled their gene expression. These analyses identified groups of genes whose tissue enrichment changed at each stage, defining gene flow dynamics of maturing vertebrate epithelia. At each of 52 and 72 hr post-fertilization (hpf), more than 60% of genes enriched in skin cells were similarly expressed in both layers, indicating that they were common epithelial genes, but many others were enriched in one layer or the other. Both expected and novel genes were enriched in periderm and basal cell layers. Genes encoding extracellular matrix, junctional, cytoskeletal, and signaling proteins were prominent among those distinguishing the two epithelial cell types. In situ hybridization and BAC transgenes confirmed our expression data and provided new tools to study zebrafish skin. Collectively, these data provide a resource for studying common and distinguishing features of maturing epithelia.

The RIPK4-IRF6 signalling axis safeguards epidermal differentiation and barrier function

The integrity of the mammalian epidermis depends on a balance of proliferation and differentiation in the resident population of stem cells¹. The kinase RIPK4 and the transcription factor IRF6 are mutated in severe developmental syndromes in humans, and mice lacking these genes display epidermal hyperproliferation and soft-tissue fusions that result in neonatal lethality^2-5. Our understanding of how these genes control epidermal differentiation is incomplete. Here we show that the role of RIPK4 in mouse development requires its kinase activity; that RIPK4 and IRF6 expressed in the epidermis regulate the same biological processes; and that the phosphorylation of IRF6 at Ser413 and Ser424 primes IRF6 for activation. Using RNA sequencing (RNA-seq), histone chromatin immunoprecipitation followed by sequencing (ChIP-seq) and assay for transposase-accessible chromatin using sequencing (ATAC-seq) of skin in wild-type and IRF6-deficient mouse embryos, we define the transcriptional programs that are regulated by IRF6 during epidermal differentiation. IRF6 was enriched at bivalent promoters, and IRF6 deficiency caused defective expression of genes that are involved in the metabolism of lipids and the formation of tight junctions. Accordingly, the lipid composition of the stratum corneum of Irf6^-/- skin was abnormal, culminating in a severe defect in the function of the epidermal barrier. Collectively, our results explain how RIPK4 and IRF6 function to ensure the integrity of the epidermis and provide mechanistic insights into why developmental syndromes that are characterized by orofacial, skin and genital abnormalities result when this axis goes awry.

Targeted sequencing reveals KLF4 gene associated with NSCL/P in Western Han Chinese

OBJECTIVE

Non-syndromic cleft lip with or without palate (NSCL/P) is one of the most common congenital birth defects affected by both gene and environmental modifiers. Recent studies revealed that KLF4 played an important role during craniofacial development. This study aimed to evaluate the association between KLF4 gene variations and NSCL/P in a Western Han Chinese population.

METHODS

In this study, 159 unrelated cases with NSCL/P were recruited. We used targeted region sequencing around KLF4 gene, and single variation association analysis and gene-based burden analysis were then performed.

RESULTS

116 SNPs and 30 indels were found. Burden analysis showed no statistical significance. Association analysis results showed rs6477574 (P = 0.014 and OR = 1.58), rs112488950 (P = 0.0024 and OR = 1.98) and rs111357138 (P = 0.0025 and OR = 1.96) were associated with NSCL/P.

CONCLUSIONS

The results revealed that KLF4 variations were associated with NSCL/P in a Western Han Chinese population, providing a new reference for future genetic study of NSCL/P.

Zebrafish and medaka as models for biomedical research of bone diseases

The identification of disease-causing mutations has in recent years progressed immensely due to whole genome sequencing approaches using patient material. The task accordingly is shifting from gene identification to functional analysis of putative disease-causing genes, preferably in an in vivo setting which also allows testing of drug candidates or biotherapeutics in whole animal disease models. In this review, we highlight the advances made in the field of bone diseases using small laboratory fish, focusing on zebrafish and medaka. We particularly highlight those human conditions where teleost models are available.

Identification of regulatory elements recapitulating early expression of L-plastin in the zebrafish enveloping layer and embryonic periderm

We have cloned and characterized an intronic fragment of zebrafish lymphocyte cytosolic protein 1 (lcp1, also called L-plastin) that drives expression to the zebrafish enveloping layer (EVL). L-plastin is a calcium-dependent actin-bundling protein belonging to the plastin/fimbrin family of proteins, and is necessary for the proper migration and attachment of several adult cell types, including leukocytes and osteoclasts. However, in zebrafish lcp1 is abundantly expressed much earlier, during differentiation of the EVL. The cells of this epithelial layer migrate collectively, spreading vegetally over the yolk. L-plastin expression persists into the larval periderm, a transient epithelial tissue that forms the first larval skin. This finding establishes that L-plastin is activated in two different embryonic waves, with a distinct regulatory switch between the early EVL and the later leukocyte. To better study L-plastin expressing cells we attempted CRISPR/Cas9 homology-driven recombination (HDR) to insert a self-cleaving peptide (Cre-P2A-EGFP-CAAX) downstream of the native lcp1 promoter. This produced a stable zebrafish line expressing Cre recombinase in EVL nuclei and green fluorescence in EVL cell membranes. In vivo tracking of these labeled cells provided enhanced views of EVL migration behavior, membrane extensions, and mitotic events. Finally, we experimentally dissected key elements of the targeted lcp1 locus, discovering a ∼300 bp intronic sequence sufficient to drive EVL expression. The lcp1: Cre-P2A-EGFP-CAAX zebrafish should be useful for studying enveloping layer specification, gastrulation movements and periderm development in this widely used vertebrate model. In addition, the conserved regulatory sequences we have isolated predict that L-plastin orthologs may have a similar early expression pattern in other vertebrate embryos.

Investigation of candidate genes of non-syndromic cleft lip with or without cleft palate, using both case-control and family-based association studies

OBJECTIVE

Non-syndromic cleft of the lip and/or palate (NSCL/P) is one of the most common polygenic diseases. In this study, both case-control and family-based association study were used to confirm whether the Single Nucleotide Polymorphisms (SNPs) were associated with NSCL/P.

METHODS

A total of 37 nuclear families and 189 controls were recruited, whose blood DNA was extracted and subjected to genotyping of SNPs of 27 candidate genes by polymerase chain reaction-improved multiple ligase detection reaction technology (PCR-iMLDR). Case-control statistical analysis was performed using the SPSS 19.0. Haplotype Relative Risk (HRR), transmission disequilibrium test (TDT), and Family-Based Association Test (FBAT) were used to test for over-transmission of the target alleles in case-parent trios. The gene-gene interactions on NSCL/P were analyzed by Unphased-3.1.4.

RESULTS

In case-control statistical analysis, only C14orf49 chr14_95932477 had statistically significant on genotype model (P = .03) and allele model (P = .03). Seven SNPs had statistically significant on TDT. None of 26 alleles has association with NSCL/P on FBAT. Some SNPs had haplotype-haplotype interactions and genotype-genotype interactions.

CONCLUSION

C14orf49 chr14_95932477 was significantly different between cases and controls on genotype model and allele model by case-control design. Seven SNPs were significantly different on HRR. Four SNPs were significantly different on TDT.

2,3,7,8-Tetrachlorodibenzo-p-dioxin and TGFβ3-Mediated Mouse Embryonic Palatal Mesenchymal Cells

2,3,7,8-Tetrachlorodibenzo-p-dioxin (TCDD) is a well-known environmental teratogenic effector for cleft palate. Transforming growth factor 3 (TGF-β3) is an essential growth factor for palatogenesis. The objective of this study is to clarify the effects of TCDD and TGF-β3 in mouse embryonic palatal mesenchymal (MEPM) cells. The effects of 10 nM TCDD, 10 ng/mL TGF-β3, or a combination of 10 nM TCDD and 10 ng/mL TGF-β3 on MEPM cells were revealed by cell and biological methods. With the increase in TCDD (0.5-10 nM), the expression of TGF-β3 increased, but at TCDD concentrations greater than 10 nM, the expression of TGF-β3 reduced. The viabilities of MEPM cells decreased in the 10 nM TCDD-treated group. But the viabilities increased in the 10 ng/mL TGF-β3-treated group, and the viabilities were intermediate in the group treated with a combination of 10 nM TCDD and 10 ng/mL TGF-β3. This phenomenon was the same as that of the motilities. In addition, we found that the expression of p-Smad2, p-Smad3,and Smad7 were increased by TCDD, TGF-β3, combination of TCDD and TGF-β3, but the expression of Smad4 were decreased by TCDD, TGF-β3, combination of TCDD and TGF-β3. These data revealed that TCDD and TGF-β3 interacted and affected MEPM cells.

Keratinocytes Isolated From Individual Cleft Lip/Palate Patients Display Variations in Their Differentiation Potential in vitro

To gain more understanding of the complex molecular processes underlying cleft lip/palate (CLP), we established a unique human cell bank, consisting of keratinocytes and corresponding fibroblasts from individual CLP patients as a new study tool. After their careful characterization, we used such patient-derived cell cultures as well as control keratinocytes for in vitro differentiation and proliferation assays. Foreskin-derived control cells as a group showed significant higher induction of the late differentiation markers Loricrin and Filaggrin than the group of CLP patients-derived keratinocytes. Additionally, we detected great variations between individual CLP keratinocyte cell cultures in regard to their potential to terminally differentiate as assessed by the induction of Loricrin and Filaggrin. Primary patient cell cultures that did not properly differentiate, exhibited high proliferation rates. Moreover, we could correlate the expression levels of transcription factor IRF6 to the ability of individual cell cultures to terminally differentiate. Using clinically relevant, patient-derived cells, our results suggest that some of the genetic predispositions causing CLP might also lead to deficiencies in keratinocyte differentiation manifested in in vitro assays.

2,3,7,8-Tetrachlorodibenzo-p-Dioxin and TGF-β3 Mediated-Mouse Embryonic Palatal Mesenchymal Cells

2,3,7,8-Tetrachlorodibenzo-p-dioxin (TCDD) is a well-known environmental teratogenic agent for cleft palate. But transforming growth factor β3 (TGF-β3) is an essential growth factor for palatogenesis. This study is to clarify effects of TCDD and TGF-β3 in mouse embryonic palatal mesenchymal (MEPM) cells. The result showed that with increase of TCDD (0.5 nM-10 nM), the expression of TGF-β3 increased, but after 10 nM TCDD, the expression of TGF-β3 reduced. The viabilities of MEPM cells decreased in 10 nM TCDD-treated group. But the viabilities increased in 10 ng/mL TGF-β3-treated group, or the viabilities were between that of them in combination of 10 nM TCDD and 10 ng/mL TGF-β3-treated group. This phenomenon was the same as the motilities. In addition, we found that the expression of phosphorylated Smad2/3 and Smad7 was increased by 10 nM TCDD, 10 ng/mL TGF-β3, or combination of 10 nM TCDD and 10 ng/mL TGF-β3 induced, but the expression of Smad4 was decreased. These data revealed that the TGF-β/Smad signaling pathway affected TCDD and TGF-β3 in MEPM cells.

Knockdown of Crispld2 in zebrafish identifies a novel network for nonsyndromic cleft lip with or without cleft palate candidate genes

Orofacial development is a multifaceted process involving tightly regulated genetic signaling networks, that when perturbed, lead to orofacial abnormalities including cleft lip and/or cleft palate. We and others have shown an association between the cysteine-rich secretory protein LCCL domain containing 2 (CRISPLD2) gene and nonsyndromic cleft lip with or without cleft palate (NSCLP). Further, we demonstrated that knockdown of Crispld2 in zebrafish alters neural crest cell migration patterns resulting in abnormal jaw and palate development. In this study, we performed RNA profiling in zebrafish embryos and identified 249 differentially expressed genes following knockdown of Crispld2. In silico pathway analysis identified a network of seven genes previously implicated in orofacial development for which differential expression was validated in three of the seven genes (CASP8, FOS, and MMP2). Single nucleotide variant (SNV) genotyping of these three genes revealed significant associations between NSCLP and FOS/rs1046117 (GRCh38 chr14:g.75746690 T > C, p = 0.0005) in our nonHispanic white (NHW) families and MMP2/rs243836 (GRCh38 chr16:g.55534236 G > A; p = 0.002) in our Hispanic families. Nominal association was found between NSCLP and CASP8/rs3769825 (GRCh38 chr2:g.202111380 C > A; p < 0.007). Overtransmission of MMP2 haplotypes were identified in the Hispanic families (p < 0.002). Significant gene-gene interactions were identified for FOS-MMP2 in the NHW families and for CASP8-FOS in the NHW simplex family subgroup (p < 0.004). Additional in silico analysis revealed a novel gene regulatory network including five of these newly identified and 23 previously reported NSCLP genes. Our results demonstrate that animal models of orofacial clefting can be powerful tools to identify novel candidate genes and gene regulatory networks underlying NSCLP.

Zebrafish models of orofacial clefts

Zebrafish is a model organism that affords experimental advantages toward investigating the normal function of genes associated with congenital birth defects. Here we summarize zebrafish studies of genes implicated in orofacial cleft (OFC). The most common use of zebrafish in this context has been to explore the normal function an OFC-associated gene product in craniofacial morphogenesis by inhibiting expression of its zebrafish ortholog. The most frequently deployed method has been to inject embryos with antisense morpholino oligonucleotides targeting the desired transcript. However, improvements in targeted mutagenesis strategies have led to widespread adoption of CRISPR/Cas9 technology. A second application of zebrafish has been for functional assays of gene variants found in OFC patients; such in vivo assays are valuable because the success of in silico methods for testing allele severity has been mixed. Finally, zebrafish have been used to test the tissue specificity of enhancers that harbor single nucleotide polymorphisms associated with risk for OFC. We review examples of each of these approaches in the context of genes that are implicated in syndromic and non-syndromic OFC. Developmental Dynamics 246:897-914, 2017. © 2017 Wiley Periodicals, Inc.

Rapid functional analysis of computationally complex rare human IRF6 gene variants using a novel zebrafish model

Large-scale sequencing efforts have captured a rapidly growing catalogue of genetic variations. However, the accurate establishment of gene variant pathogenicity remains a central challenge in translating personal genomics information to clinical decisions. Interferon Regulatory Factor 6 (IRF6) gene variants are significant genetic contributors to orofacial clefts. Although approximately three hundred IRF6 gene variants have been documented, their effects on protein functions remain difficult to interpret. Here, we demonstrate the protein functions of human IRF6 missense gene variants could be rapidly assessed in detail by their abilities to rescue the irf6^-/- phenotype in zebrafish through variant mRNA microinjections at the one-cell stage. The results revealed many missense variants previously predicted by traditional statistical and computational tools to be loss-of-function and pathogenic retained partial or full protein function and rescued the zebrafish irf6^-/- periderm rupture phenotype. Through mRNA dosage titration and analysis of the Exome Aggregation Consortium (ExAC) database, IRF6 missense variants were grouped by their abilities to rescue at various dosages into three functional categories: wild type function, reduced function, and complete loss-of-function. This sensitive and specific biological assay was able to address the nuanced functional significances of IRF6 missense gene variants and overcome many limitations faced by current statistical and computational tools in assigning variant protein function and pathogenicity. Furthermore, it unlocked the possibility for characterizing yet undiscovered human IRF6 missense gene variants from orofacial cleft patients, and illustrated a generalizable functional genomics paradigm in personalized medicine.

Advances in sequencing technologies have led to rapid increases in personalized genetics information. Millions of differences exist when comparing the genomes of two individuals, accounting for the diversity of humans and occasionally disease pathogenicity. Accurate determination of the functional consequences of these differences is critical for translating personal genetic information into clinical decision. Various methods have been devised to meet this challenge. However, they are imperfect, especially in their abilities to interpret rare variants. Rare variants must be evaluated by multi-pronged approaches to establish disease pathogenicity, one crucial approach being functional testing through experimental models. Here, we utilized a zebrafish model to rapidly evaluate the protein functions of rare gene variants of IRF6, a gene of established importance in orofacial cleft pathogenesis. IRF6 functions are well-conserved across species, and allowed us to test the functions of human IRF6 variants by their abilities to rescue zebrafish embryos depleted of irf6. Many variants previously labeled pathogenic and loss-of-function retained full wild type-level protein activity, suggesting they could potentially represent benign, rather than disease-causing genetic variations. This paradigm is applicable to other genes and will allow researchers to rapidly triage gene variants for further study and clinicians to provide more accurate genetic diagnoses.

ARHGAP29 Mutation Is Associated with Abnormal Oral Epithelial Adhesions

Nonsyndromic cleft lip and/or palate (NSCL/P) is a prevalent birth defect of complex etiology. Previous studies identified mutations in ARHGAP29 associated with an increased risk for NSCL/P. To investigate the effects of ARHGAP29 in vivo, we generated a novel murine allele by inserting a point mutation identified in a patient with NSCL/P. This single-nucleotide variation of ARHGAP29 translates to an early nonsense mutation (K326X), presumably resulting in loss-of-function (LoF). Embryos from Arhgap29^K326X/+ intercrosses were harvested at various time points. No homozygous Arhgap29^K326X animals were found in the 45 analyzed litters, assessed as early as embryonic day 8.5 (e8.5). Coronal sectioning of e13.5 and e14.5 heads revealed that 59% of Arhgap29^K326X/+ mice ( n = 37) exhibited improper epithelial contact between developing oral structures, while none were observed in wild types ( n = 10). In addition, Arhgap29^K326X/+ embryos exhibited a significantly higher percentage of maxillary epithelium in contact with mandibular epithelium. Immunofluorescent analyses of the periderm and oral adhesions revealed the presence of Arhgap29 in periderm cells. These cells were p63 negative, keratin 17 positive, and keratin 6 positive and present at sites of adhesion, although occasionally disorganized. Oral adhesions did not appear to impair palatogenesis, as all analyzed Arhgap29^K326X/+ embryos showed confluent palatal mesenchyme and epithelium at e18.5 ( n = 16), and no mice were found with a cleft at birth. Collectively, our data demonstrate that ARHGAP29 is required for embryonic survival and that heterozygosity for LoF variants of Arhgap29 increases the incidence and length of oral adhesions at a critical time point during orofacial development. In conclusion, we validate the LoF nature of the human K326X mutation in vivo and reveal a previously unknown effect of Arhgap29 in murine craniofacial development.

Periderm: Life-cycle and function during orofacial and epidermal development

Development of the secondary palate involves a complex series of embryonic events which, if disrupted, result in the common congenital anomaly cleft palate. The secondary palate forms from paired palatal shelves which grow initially vertically before elevating to a horizontal position above the tongue and fusing together in the midline via the medial edge epithelia. As the epithelia of the vertical palatal shelves are in contact with the mandibular and lingual epithelia, pathological fusions between the palate and the mandible and/or the tongue must be prevented. This function is mediated by the single cell layered periderm which forms in a distinct and reproducible pattern early in embryogenesis, exhibits highly polarised expression of adhesion complexes, and is shed from the outer surface as the epidermis acquires its barrier function. Disruption of periderm formation and/or function underlies a series of birth defects that exhibit multiple inter-epithelial adhesions including the autosomal dominant popliteal pterygium syndrome and the autosomal recessive cocoon syndrome and Bartsocas Papas syndrome. Genetic analyses of these conditions have shown that IRF6, IKKA, SFN, RIPK4 and GRHL3, all of which are under the transcriptional control of p63, play a key role in periderm formation. Despite these observations, the medial edge epithelia must rapidly acquire the capability to fuse if the palatal shelves are not to remain cleft. This process is driven by TGFβ3-mediated, down-regulation of p63 in the medial edge epithelia which allows periderm migration out of the midline epithelial seam and reduces the proliferative potential of the midline epithelial seam thereby preventing cleft palate. Together, these findings indicate that periderm plays a transient but fundamental role during embryogenesis in preventing pathological adhesion between intimately apposed, adhesion-competent epithelia.

Alternative splicing and start sites: Lessons from the Grainyhead-like family

The two main mechanisms that expand the proteomic output of eukaryotic genes are alternative splicing and alternative translation initiation signals. Despite being essential to generate isoforms of gene products that create functional diversity during development, the impact of these mechanisms on fine-tuning regulatory gene networks is still underappreciated. In this review, we use the Grainyhead-like (Grhl) family as a case study to illustrate the importance of isoforms when investigating transcription factor family function during development and disease, and highlight the potential for differential modulation of downstream target genes. We provide insights into the importance of considering alternative gene products when designing, undertaking, and analysing primary research, and the effect that isoforms may have on development. This review also covers known mutations in Grhl family members, and postulates how genetic changes may dictate transcriptional specificity between the Grhl family members. It also contrasts and compares the available literature on the function and importance of the Grhl isoforms, and highlights current gaps in our understanding of their regulatory gene networks in development and disease.

p63 exerts spatio-temporal control of palatal epithelial cell fate to prevent cleft palate

Cleft palate is a common congenital disorder that affects up to 1 in 2500 live births and results in considerable morbidity to affected individuals and their families. The aetiology of cleft palate is complex with both genetic and environmental factors implicated. Mutations in the transcription factor p63 are one of the major individual causes of cleft palate; however, the gene regulatory networks in which p63 functions remain only partially characterized. Our findings demonstrate that p63 functions as an essential regulatory molecule in the spatio-temporal control of palatal epithelial cell fate to ensure appropriate fusion of the palatal shelves. Initially, p63 induces periderm formation and controls its subsequent maintenance to prevent premature adhesion between adhesion-competent, intra-oral epithelia. Subsequently, TGFβ3-induced down-regulation of p63 in the medial edge epithelia of the palatal shelves is a pre-requisite for palatal fusion by facilitating periderm migration from, and reducing the proliferative potential of, the midline epithelial seam thereby preventing cleft palate.

Cleft palate is a serious congenital condition which affects approximately 1 in every 2500 births. Cleft palate occurs when the palatal shelves fail to grow, adhere or fuse during development. Mutations in the gene encoding the transcription factor p63 result in cleft palate in humans and mice. However, the role of p63 and how it controls the network of genes to regulate palate development is not well understood.In this study, we demonstrate that p63 controls the spatio-temporal regulation of palatal epithelial cell fate to ensure appropriate palatal adhesion: p63 induces the formation of a flattened layer of epithelial (periderm) cells and controls its subsequent maintenance. We also demonstrate that TGFβ3-induced, down-regulation of p63 in the medial edge epithelial cells, through which the palatal shelves adhere and fuse, controls Jag2-induced periderm migration to the oral and nasal epithelial triangles. In addition, p63 plays a central role in maintaining the proliferative potential of the basal layer of the medial edge epithelia. Our study provides significant new insights into the mechanisms that regulate development of the palate by establishing the role of p63 in governing the fate of the midline epithelial cells.

Association studies of low-frequency coding variants in nonsyndromic cleft lip with or without cleft palate

Nonsyndromic cleft lip with or without cleft palate (NSCL/P) is a group of common human birth defects with complex etiology. Although genome-wide association studies have successfully identified a number of risk loci, these loci only account for about 20% of the heritability of orofacial clefts. The "missing" heritability may be found in rare variants, copy number variants, or interactions. In this study, we investigated the role of low-frequency variants genotyped in 1995 cases and 1626 controls on the Illumina HumanCore + Exome chip. We performed two statistical tests, Sequence Kernel Association Test (SKAT) and Combined Multivariate and Collapsing (CMC) method using two minor allele frequency cutoffs (1% and 5%). We found that a burden of low-frequency coding variants in N4BP2, CDSN, PRTG, and AHRR were associated with increased risk of NSCL/P. Low-frequency variants in other genes were associated with decreased risk of NSCL/P. These results demonstrate that low-frequency variants contribute to the genetic etiology of NSCL/P.

Identification of common non-coding variants at 1p22 that are functional for non-syndromic orofacial clefting

Genome-wide association studies (GWAS) do not distinguish between single nucleotide polymorphisms (SNPs) that are causal and those that are merely in linkage-disequilibrium with causal mutations. Here we describe a versatile, functional pipeline and apply it to SNPs at 1p22, a locus identified in several GWAS for non-syndromic cleft lip with or without cleft palate (NS CL/P). First we amplified DNA elements containing the ten most-highly risk-associated SNPs and tested their enhancer activity in vitro, identifying three SNPs with allele-dependent effects on such activity. We then used in vivo reporter assays to test the tissue-specificity of these enhancers, chromatin configuration capture to test enhancer-promoter interactions, and genome editing in vitro to show allele-specific effects on ARHGAP29 expression and cell migration. Our results further indicate that two SNPs affect binding of CL/P-associated transcription factors, and one affects chromatin configuration. These results translate risk into potential mechanisms of pathogenesis.

Nonsyndromic cleft lip with or without cleft palate and cancer: Evaluation of a possible common genetic background through the analysis of GWAS data

Previous research suggests a genetic overlap between nonsyndromic cleft lip with or without cleft palate (NSCL/P) and cancer. The aim of the present study was to identify common genetic risk loci for NSCL/P and cancer entities that have been reported to co-occur with orofacial clefting. This was achieved through the investigation of large genome-wide association study datasets. Investigations of 12 NSCL/P single nucleotide polymorphisms (SNPs) in 32 cancer datasets, and 204 cancer SNPs in two NSCL/P datasets, were performed. The SNPs rs13041247 (20q12) and rs6457327 (6p21.33) showed suggestive evidence for an association with both NSCL/P and a specific cancer entity. These loci harbor genes of biological relevance to oncogenesis (MAFB and OCT4, respectively). This study is the first to characterize possible pleiotropic risk loci for NSCL/P and cancer in a systematic manner. The data represent a starting point for future research by identifying a genetic link between NSCL/P and cancer.

Small teleost fish provide new insights into human skeletal diseases

Small teleost fish such as zebrafish and medaka are increasingly studied as models for human skeletal diseases. Efficient new genome editing tools combined with advances in the analysis of skeletal phenotypes provide new insights into fundamental processes of skeletal development. The skeleton among vertebrates is a highly conserved organ system, but teleost fish and mammals have evolved unique traits or have lost particular skeletal elements in each lineage. Several unique features of the skeleton relate to the extremely small size of early fish embryos and the small size of adult fish used as models. A detailed analysis of the plethora of interesting skeletal phenotypes in zebrafish and medaka pushes available skeletal imaging techniques to their respective limits and promotes the development of new imaging techniques. Impressive numbers of zebrafish and medaka mutants with interesting skeletal phenotypes have been characterized, complemented by transgenic zebrafish and medaka lines. The advent of efficient genome editing tools, such as TALEN and CRISPR/Cas9, allows to introduce targeted deficiencies in genes of model teleosts to generate skeletal phenotypes that resemble human skeletal diseases. This review will also discuss other attractive aspects of the teleost skeleton. This includes the capacity for lifelong tooth replacement and for the regeneration of dermal skeletal elements, such as scales and fin rays, which further increases the value of zebrafish and medaka models for skeletal research.