Interferon regulatory factor 6 promotes differentiation of the periderm by activating expression of Grainyhead-like 3

Neural crest and periderm-specific requirements of Irf6 during neural tube and craniofacial development

IRF6 is a key genetic determinant of cleft lip and palate. The ability to interrogate post-embryonic requirements of Irf6 has been hindered, as global Irf6 ablation in the mouse causes neonatal lethality. Prior work analyzing Irf6 in mice defined its role in the embryonic surface epithelium and periderm, where it regulates cell proliferation and differentiation. Several reports have also described Irf6 expression in other cell types, such as muscle, and neuroectoderm. However, analysis of a functional role in non-epithelial cells has been incomplete due to the severity and lethality of the Irf6 knockout model and the paucity of work with a conditional Irf6 allele. Here we describe the generation and characterization of a new Irf6 floxed mouse model and analysis of Irf6 ablation in periderm and neural crest lineages. This work found that loss of Irf6 in periderm recapitulates a mild Irf6 null phenotype, suggesting that Irf6-mediated signaling in periderm plays a crucial role in regulating embryonic development. Further, conditional ablation of Irf6 in neural crest cells resulted in an anterior neural tube defect of variable penetrance. The generation of this conditional Irf6 allele allows for new insights into craniofacial development and new exploration into the post-natal role of Irf6.

Glucose modulates IRF6 transcription factor dimerization to enable epidermal differentiation

Non-energetic roles for glucose are largely unclear, as is the interplay between transcription factors (TFs) and ubiquitous biomolecules. Metabolomic analyses uncovered elevation of intracellular glucose during differentiation of diverse cell types. Human and mouse tissue engineered with glucose sensors detected a glucose gradient that peaked in the outermost differentiated layers of the epidermis. Free glucose accumulation was essential for epidermal differentiation and required the SGLT1 glucose transporter. Glucose affinity chromatography uncovered glucose binding to diverse regulatory proteins, including the IRF6 TF. Direct glucose binding enabled IRF6 dimerization, DNA binding, genomic localization, and induction of IRF6 target genes, including essential pro-differentiation TFs GRHL1, GRHL3, HOPX, and PRDM1. These data identify a role for glucose as a gradient morphogen that modulates protein multimerization in cellular differentiation.

ARHGAP29 promotes keratinocyte proliferation and migration in vitro and is dispensable for in vivo wound healing

BACKGROUND

RhoA GTPases play critical roles in actin cytoskeletal remodeling required for controlling a diverse range of cellular functions including cell proliferation, adhesion, migration and changes in cell shape, all required for cutaneous wound healing. RhoA cycles between an active GTP-bound and an inactive GDP-bound form, a process regulated by guanine nucleotide exchange factors (GEFs) and GTPase-activating proteins (GAPs). ARHGAP29 is a GAP expressed in skin keratinocytes and is decreased in the absence of interferon regulator factor 6, a critical regulator of cell proliferation, migration, and wound healing. However, the role for ARHGAP29 in keratinocyte biology is unknown.

RESULTS

We generated ARHGAP29 knockdown keratinocyte cell lines and show they displayed increased filamentous actin, phospho-myosin regulatory light chain, cell area and population doubling time. Furthermore, we found that ARHGAP29 knockdown keratinocytes displayed significant delays in scratch wound closure in both single and collective cell migration conditions; these delays were rescued by both adding back ARHGAP29 or adding a ROCK inhibitor to ARHGAP29 knockdown cells. In vivo, however, Arhgap29 heterozygotes or keratinocyte-specific knockouts showed on-time wound healing.

CONCLUSIONS

These data demonstrate that ARHGAP29 is required for keratinocyte morphology, proliferation and migration in vitro but is dispensable during wound healing in vivo.

Rare variants in PRKCI cause Van der Woude syndrome and other features of peridermopathy

Van der Woude syndrome (VWS) is an autosomal dominant disorder characterized by lower lip pits and orofacial clefts (OFCs). With a prevalence of approximately 1 in 35,000 live births, it is the most common form of syndromic clefting and may account for ~2% of all OFCs. The majority of VWS is attributed to genetic variants in IRF6 (~70%) or GRHL3 (~5%), leaving up to 25% of individuals with VWS without a molecular diagnosis. Both IRF6 and GRHL3 function in a transcriptional regulatory network governing differentiation of periderm, a single layer of epithelial cells that prevents pathological adhesions during palatogenesis. Disruption of this layer results in a spectrum of phenotypes ranging from lip pits and OFCs to severe pterygia and other congenital anomalies that can be incompatible with life. Understanding the mechanisms of peridermopathies is vital in improving health outcomes for affected individuals. We reasoned that genes encoding additional members of the periderm gene regulatory network, including kinases acting upstream of IRF6 (i.e., atypical protein kinase C family members, RIPK4, and CHUK), are candidates to harbor variants resulting in VWS. Consistent with this prediction, we identified 6 de novo variants (DNs) and 11 rare variants in PRKCI, an atypical protein kinase C, in 17 individuals with clinical features consistent with syndromic OFCs and peridermopathies. Of the identified DNs, 4 were identical p.(Asn383Ser) variants in unrelated individuals with syndromic OFCs, indicating a likely hotspot mutation. We also performed functional validation of 12 variants using the enveloping layer in zebrafish embryos, a structure analogous to the periderm. Three patient-specific alleles (p.Arg130His, p.(Asn383Ser), and p.Leu385Phe) were confirmed to be loss-of-function variants. In summary, we identified PRKCI as a novel causal gene for VWS and syndromic OFC with other features of peridermopathies.

A scATAC-seq atlas of stasis zone in rat skin burn injury wound process

Burn injuries often leave behind a "stasis zone", a region of tissue critically important for determining both the severity of the injury and the potential for recovery. To understand the intricate cellular and epigenetic changes occurring within this critical zone, we utilized single-cell assay for transposase-accessible chromatin sequencing (scATAC-seq) to profile over 31,500 cells from both healthy rat skin and the stasis zone at nine different time points after a burn injury. This comprehensive approach revealed 26 distinct cell types and the dynamic shifts in the proportions of these cell types over time. We observed distinct gene activation patterns in different cell types at various stages post-burn, highlighting key players in immune activation, tissue regeneration, and blood vessel repair. Importantly, our analysis uncovered the regulatory networks governing these genes, offering valuable insights into the intricate mechanisms orchestrating burn wound healing. This comprehensive cellular and molecular atlas of the stasis zone provides a powerful resource for developing targeted therapies aimed at improving burn injury recovery and minimizing long-term consequences.

High incidence and geographic distribution of cleft palate in Finland are associated with the IRF6 gene

In Finland, the frequency of isolated cleft palate (CP) is higher than that of isolated cleft lip with or without cleft palate (CL/P). This trend contrasts to that in other European countries but its genetic underpinnings are unknown. We conducted a genome-wide association study in the Finnish population and identified rs570516915, a single nucleotide polymorphism highly enriched in Finns, as strongly associated with CP (P = 5.25 × 10-34, OR = 8.65, 95% CI 6.11-12.25), but not with CL/P (P = 7.2 × 10-5), with genome-wide significance. The risk allele frequency of rs570516915 parallels the regional variation of CP prevalence in Finland, and the association was replicated in independent cohorts of CP cases from Finland (P = 8.82 × 10-28) and Estonia (P = 1.25 × 10-5). The risk allele of rs570516915 alters a conserved binding site for the transcription factor IRF6 within an enhancer (MCS-9.7) upstream of the IRF6 gene and diminishes the enhancer activity. Oral epithelial cells derived from CRISPR-Cas9 edited induced pluripotent stem cells demonstrate that the CP-associated allele of rs570516915 concomitantly decreases the binding of IRF6 and the expression level of IRF6, suggesting impaired IRF6 autoregulation as a molecular mechanism underlying the risk for CP.

Combining genetic and single-cell expression data reveals cell types and novel candidate genes for orofacial clefting

Non-syndromic cleft lip with/without cleft palate (nsCL/P) is one of the most common birth defects and has a multifactorial etiology. To date, over 45 loci harboring common risk variants have been identified. However, the effector genes at these loci, and the cell types that are affected by risk alleles, remain largely unknown. To address this, we combined genetic data from an nsCL/P genome-wide association study (GWAS) with single-cell RNA sequencing data obtained from the heads of unaffected human embryos. Using the recently developed single-cell disease relevance score (scDRS) approach, we identified two major cell types involved in nsCL/P development, namely the epithelium and the HAND2+ pharyngeal arches (PA). Combining scDRS with co-expression networks and differential gene expression analysis, we prioritized nsCL/P candidate genes, some of which were additionally supported by GWAS data (e.g., CTNND1, PRTG, RPL35A, RAB11FIP1, KRT19). Our results suggest that specific epithelial and PA sub-cell types are involved in nsCL/P development, and harbor a substantial fraction of the genetic risk for nsCL/P.

Single-cell transcriptome and chromatin accessibility mapping of upper lip and primary palate fusion

Cleft lip and/or primary palate (CL/P) represent a prevalent congenital malformation, the aetiology of which is highly intricate. Although it is generally accepted that the condition arises from failed fusion between the upper lip and primary palate, the precise mechanism underlying this fusion process remains enigmatic. In this study, we utilized transposase-accessible chromatin sequencing (scATAC-seq) and single-cell RNA sequencing (scRNA-seq) to interrogate lambdoidal junction tissue derived from C57BL/6J mouse embryos at critical stages of embryogenesis (10.5, 11.5 and 12.5 embryonic days). We successfully identified distinct subgroups of mesenchymal and ectodermal cells involved in the fusion process and characterized their unique transcriptional profiles. Furthermore, we conducted cell differentiation trajectory analysis, revealing a dynamic repertoire of genes that are sequentially activated or repressed during pseudotime, facilitating the transition of relevant cell types. Additionally, we employed scATAC data to identify key genes associated with the fusion process and demonstrated differential chromatin accessibility across major cell types. Finally, we constructed a dynamic intercellular communication network and predicted upstream transcriptional regulators of critical genes involved in important signalling pathways. Our findings provide a valuable resource for future studies on upper lip and primary palate development, as well as congenital defects.

Dissecting the regulatory logic of specification and differentiation during vertebrate embryogenesis

The interplay between transcription factors and chromatin accessibility regulates cell type diversification during vertebrate embryogenesis. To systematically decipher the gene regulatory logic guiding this process, we generated a single-cell multi-omics atlas of RNA expression and chromatin accessibility during early zebrafish embryogenesis. We developed a deep learning model to predict chromatin accessibility based on DNA sequence and found that a small number of transcription factors underlie cell-type-specific chromatin landscapes. While Nanog is well-established in promoting pluripotency, we discovered a new function in priming the enhancer accessibility of mesendodermal genes. In addition to the classical stepwise mode of differentiation, we describe instant differentiation, where pluripotent cells skip intermediate fate transitions and terminally differentiate. Reconstruction of gene regulatory interactions reveals that this process is driven by a shallow network in which maternally deposited regulators activate a small set of transcription factors that co-regulate hundreds of differentiation genes. Notably, misexpression of these transcription factors in pluripotent cells is sufficient to ectopically activate their targets. This study provides a rich resource for analyzing embryonic gene regulation and reveals the regulatory logic of instant differentiation.

Functional analysis of ESRP1/2 gene variants and CTNND1 isoforms in orofacial cleft pathogenesis

Orofacial cleft (OFC) is a common human congenital anomaly. Epithelial-specific RNA splicing regulators ESRP1 and ESRP2 regulate craniofacial morphogenesis and their disruption result in OFC in zebrafish, mouse and humans. Using esrp1/2 mutant zebrafish and murine Py2T cell line models, we functionally tested the pathogenicity of human ESRP1/2 gene variants. We found that many variants predicted by in silico methods to be pathogenic were functionally benign. Esrp1 also regulates the alternative splicing of Ctnnd1 and these genes are co-expressed in the embryonic and oral epithelium. In fact, over-expression of ctnnd1 is sufficient to rescue morphogenesis of epithelial-derived structures in esrp1/2 zebrafish mutants. Additionally, we identified 13 CTNND1 variants from genome sequencing of OFC cohorts, confirming CTNND1 as a key gene in human OFC. This work highlights the importance of functional assessment of human gene variants and demonstrates the critical requirement of Esrp-Ctnnd1 acting in the embryonic epithelium to regulate palatogenesis.

High incidence and geographic distribution of cleft palate cases in Finland are associated with a regulatory variant in IRF6

In Finland the frequency of isolated cleft palate (CP) is higher than that of isolated cleft lip with or without cleft palate (CL/P). This trend contrasts to that in other European countries but its genetic underpinnings are unknown. We performed a genome-wide association study for orofacial clefts, which include CL/P and CP, in the Finnish population. We identified rs570516915, a single nucleotide polymorphism that is highly enriched in Finns and Estonians, as being strongly associated with CP ( P = 5.25 × 10 -34 , OR = 8.65, 95% CI 6.11-12.25), but not with CL/P ( P = 7.2 × 10 -5 ), with genome-wide significance. The risk allele frequency of rs570516915 parallels the regional variation of CP prevalence in Finland, and the association was replicated in independent cohorts of CP cases from Finland ( P = 8.82 × 10 -28 ) and Estonia ( P = 1.25 × 10 -5 ). The risk allele of rs570516915 disrupts a conserved binding site for the transcription factor IRF6 within a previously characterized enhancer upstream of the IRF6 gene. Through reporter assay experiments we found that the risk allele of rs570516915 diminishes the enhancer activity. Oral epithelial cells derived from CRISPR-Cas9 edited induced pluripotent stem cells demonstrate that the CP-associated allele of rs570516915 concomitantly decreases the binding of IRF6 and the expression level of IRF6 , suggesting impaired IRF6 autoregulation as a molecular mechanism underlying the risk for CP.

Functional analysis of ESRP1/2 gene variants and CTNND1 isoforms in orofacial cleft pathogenesis

Orofacial cleft (OFC) is a common human congenital anomaly. Epithelial-specific RNA splicing regulators ESRP1 and ESRP2 regulate craniofacial morphogenesis and their disruption result in OFC in zebrafish, mouse and humans. Using esrp1/2 mutant zebrafish and murine Py2T cell line models, we functionally tested the pathogenicity of human ESRP1/2 gene variants. We found that many variants predicted by in silico methods to be pathogenic were functionally benign. Esrp1 also regulates the alternative splicing of Ctnnd1 and these genes are co-expressed in the embryonic and oral epithelium. In fact, over-expression of ctnnd1 is sufficient to rescue morphogenesis of epithelial-derived structures in esrp1/2 zebrafish mutants. Additionally, we identified 13 CTNND1 variants from genome sequencing of OFC cohorts, confirming CTNND1 as a key gene in human OFC. This work highlights the importance of functional assessment of human gene variants and demonstrates the critical requirement of Esrp-Ctnnd1 acting in the embryonic epithelium to regulate palatogenesis.

Neural crest and periderm-specific requirements of Irf6 during neural tube and craniofacial development

IRF6 is a key genetic determinant of syndromic and non-syndromic cleft lip and palate. The ability to interrogate post-embryonic requirements of Irf6 has been hindered, as global Irf6 ablation in the mouse causes neonatal lethality. Prior work analyzing Irf6 in mouse models defined its role in the embryonic surface epithelium and periderm where it is required to regulate cell proliferation and differentiation. Several reports have also described Irf6 gene expression in other cell types, such as muscle, and neuroectoderm. However, analysis of a functional role in non-epithelial cell lineages has been incomplete due to the severity and lethality of the Irf6 knockout model and the paucity of work with a conditional Irf6 allele. Here we describe the generation and characterization of a new Irf6 floxed mouse model and analysis of Irf6 ablation in periderm and neural crest lineages. This work found that loss of Irf6 in periderm recapitulates a mild Irf6 null phenotype, suggesting that Irf6-mediated signaling in periderm plays a crucial role in regulating embryonic development. Further, conditional ablation of Irf6 in neural crest cells resulted in an anterior neural tube defect of variable penetrance. The generation of this conditional Irf6 allele allows for new insights into craniofacial development and new exploration into the post-natal role of Irf6.

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.

Molecular and spatial landmarks of early mouse skin development

A wealth of specialized cell populations within the skin facilitates its hair-producing, protective, sensory, and thermoregulatory functions. How the vast cell-type diversity and tissue architecture develops is largely unexplored. Here, with single-cell transcriptomics, spatial cell-type assignment, and cell-lineage tracing, we deconstruct early embryonic mouse skin during the key transitions from seemingly uniform developmental precursor states to a multilayered, multilineage epithelium, and complex dermal identity. We identify the spatiotemporal emergence of hair-follicle-inducing, muscle-supportive, and fascia-forming fibroblasts. We also demonstrate the formation of the panniculus carnosus muscle (PCM), sprouting blood vessels without pericyte coverage, and the earliest residence of mast and dendritic immune cells in skin. Finally, we identify an unexpected epithelial heterogeneity within the early single-layered epidermis and a signaling-rich periderm layer. Overall, this cellular and molecular blueprint of early skin development-which can be explored at https://kasperlab.org/tools-establishes histological landmarks and highlights unprecedented dynamic interactions among skin cells.

Functions of the RIP kinase family members in the skin

The receptor interacting protein kinases (RIPK) are a family of serine/threonine kinases that are involved in the integration of various stress signals. In response to several extracellular and/or intracellular stimuli, RIP kinases engage signaling cascades leading to the activation of NF-κB and mitogen-activated protein kinases, cell death, inflammation, differentiation and Wnt signaling and can have kinase-dependent and kinase-independent functions. Although it was previously suggested that seven RIPKs are part of the RIPK family, phylogenetic analysis indicates that there are only five genuine RIPKs. RIPK1 and RIPK3 are mainly involved in controlling and executing necroptosis in keratinocytes, while RIPK4 controls proliferation and differentiation of keratinocytes and thereby can act as a tumor suppressor in skin. Therefore, in this review we summarize and discuss the functions of RIPKs in skin homeostasis as well as the signaling pathways involved.

Simultaneous Occurrence of Hypospadias and Bilateral Cleft Lip and Jaw in a Crossbred Calf: Clinical, Computer Tomographic, and Genomic Characterization

Congenital abnormalities in animals, including abnormalities of the cleft lip and jaw and hypospadias have been reported in all domesticated species. They are a major concern for breeders due to the increased economic loss they entail. In this article, we described a congenital bilateral cheilognathoschisis (cleft lip and jaw) with campylognathia in association with penile hypospadias and preputial hypoplasia with failure of preputial fusion in a Bos taurus crossbred Piedmontese × Wagyu calf. Clinical examination, computed tomography, and whole genome sequencing were performed to describe and identify a possible cause of the abnormalities. Clinical examination revealed a bilateral cheilognathoschisis of approximately 4 cm in length and 3 cm in width in the widest part, with computer tomography analyses confirming the bilateral absence of the processus nasalis of the incisive bone and the lateral deviation of the processus palatinus towards the left side. Genomic data analyses identified 13 mutations with a high impact on the products of the following overlapped genes: ACVR1, ADGRA2, BHMT2, BMPR1B, CCDC8, CDH1, EGF, F13A1, GSTP1, IRF6, MMP14, MYBPHL, and PHC2 with ADGRA2, EGF, F13A1, GSTP1, and IRF6 having mutations in a homozygous state. The whole genome investigation indicates the involvement of multiple genes in the birth defects observed in this case.

ARHGAP29 is required for keratinocyte proliferation and migration

BACKGROUND

RhoA GTPase plays critical roles in actin cytoskeletal remodeling required for controlling a diverse range of cellular functions including cell proliferation, cell adhesions, migration and changes in cell shape. RhoA cycles between an active GTP-bound and an inactive GDP-bound form, a process that is regulated by guanine nucleotide exchange factors (GEFs), and GTPase-activating proteins (GAPs). ARHGAP29 is a GAP expressed in keratinocytes of the skin and is decreased in the absence of Interferon Regulator Factor 6, a critical regulator of cell proliferation and migration. However, the role for ARHGAP29 in keratinocyte biology is unknown.

RESULTS

Novel ARHGAP29 knockdown keratinocyte cell lines were generated using both CRISPR/Cas9 and shRNA technologies. Knockdown cells exhibited significant reduction of ARHGAP29 protein (50-80%) and displayed increased filamentous actin (stress fibers), phospho-myosin light chain (contractility), cell area and population doubling time. Furthermore, we found that ARHGAP29 knockdown keratinocytes displayed significant delays in scratch wound closure in both single cell and collective cell migration conditions. Particularly, our results show a reduction in path lengths, speed, directionality and persistence in keratinocytes with reduced ARHGAP29. The delay in scratch closure was rescued by both adding back ARHGAP29 or adding a ROCK inhibitor to ARHGAP29 knockdown cells.

CONCLUSIONS

These data demonstrate that ARHGAP29 is required for keratinocyte morphology, proliferation and migration mediated through the RhoA pathway.

The IRF6 AP-2α binding site polymorphism relate to the severity of non-syndromic orofacial cleft of Indonesian patients

BACKGROUND

IRF6 AP-2α binding site polymorphism is known as IRF6 rs642961. It has been associated with a nonsyndromic orofacial cleft (NS OFC). This study aimed to determine the IRF6 rs642961 as a risk factor associated with NS OFC and its phenotypes.

METHODS

The case-control design used for 264 subjects consists of 158 NS CLP subjects (42 CU CLP, 34 CB CLP, 33 CLO, 49 CPOs) and 106 healthy controls. The DNA is extracted from venous blood. The segment of IRF6 rs642961 amplified by polymerase chain reaction (PCR) followed by restriction fragment length of polymorphisms (RFLPs) used the MspI digestion enzyme. The qPCR method to identify the mRNA expression levels of the IRF6 gene rs642961 was analyzed by the Livak method.

RESULTS

The study results show that in NS CB CLP phenotype as the most severe phenotype of NS OFC, the Odds Ratio (OR) of A mutant allele was 5.094 (CI=1.456-17.820; P=0.011) and the OR of AA homozygous mutant genotype was 13.481 (CI=2.648-68.635; P=0.001). There are different levels of mRNA expression changes from NS OFC and its phenotypes. It is substantial among the 2^-ΔΔCt and the group of AA, GA, and GG genotypes (P<0.05); in the NS CPO phenotype, it shows IRF6 mRNA under-expression in GA, AA genotypes while in other phenotypes it shows IRF6 mRNA overexpression.

CONCLUSIONS

The IRF6 AP-2α binding site polymorphism is strongly associated with the severity of NS OFC, and this polymorphism has a functional role in affecting IRF6 mRNA expression that is variable in each phenotype.

Negative regulation of interferon regulatory factor 6 (IRF6) in interferon and NF-κB signalling pathways of common carp (Cyprinus carpio L.)

BACKGROUND

Interferon (IFN) regulatory factors (IRFs) is a kind of transcription factors, which play an important role in regulating the expression of type I IFN and related genes. In mammals, IRF6 is not relevant with IFN expression, while zebrafish IRF6 was reported to be a positive regulator of IFN expression and could be phosphorylated by both MyD88 and TBK1. However, the role of IRF6 in the immune response and IFN transcription of common carp is unknown.

RESULTS

In the present study, the cDNA of IRF6 gene (CcIRF6) was cloned from common carp using RACE technique, with a total length of 1905 bp, encoding 471 amino acid residues, which possesses two functional domains of DBD and IAD. Similarity analysis showed that CcIRF6 had more than 50% similarity with IRFs of other vertebrates, and had the highest similarity with grass carp and zebrafish, among which the DBD domain was much more conserved. The phylogenetic analysis showed that CcIRF6 is in the branch of Osteichthyes and has the closest relationship with grass carp. In healthy common carp, the CcIRF6 was expressed in all the examined tissues, with the highest level in the oral epithelium, and the lowest level in the head kidney. After intraperitoneal injection of poly(I:C) or Aeromonas hydrophila, the expression of CcIRF6 increased in spleen, head kidney, foregut and hindgut of common carp. Moreover, poly(I:C), LPS, PGN and flagellin induced the expression of CcIRF6 in peripheral leukocytes and head kidney leukocytes of common carp in vitro. In EPC cells, CcIRF6 inhibited the expression of some IFN-related genes and pro-inflammatory cytokines, and dual luciferase reporter assay showed that CcIRF6 reduced the activity of IFN and NF-κB reporter genes.

CONCLUSIONS

The present study suggests that CcIRF6 is involved in the antiviral and antibacterial immune response of common carp, and negatively regulate the expression of IFN and NF-κB signalling pathways, which provides a theoretical basis for the study and prevention of fish disease pathogenesis.

Origin, form and function of extraembryonic structures in teleost fishes

Teleost eggs have evolved a highly derived early developmental pattern within vertebrates as a result of the meroblastic cleavage pattern, giving rise to a polar stratified architecture containing a large acellular yolk and a small cellular blastoderm on top. Besides the acellular yolk, the teleost-specific yolk syncytial layer (YSL) and the superficial epithelial enveloping layer are recognized as extraembryonic structures that play critical roles throughout embryonic development. They provide enriched microenvironments in which molecular feedback loops, cellular interactions and mechanical signals emerge to sculpt, among other things, embryonic patterning along the dorsoventral and left-right axes, mesendodermal specification and the execution of morphogenetic movements in the early embryo and during organogenesis. An emerging concept points to a critical role of extraembryonic structures in reinforcing early genetic and morphogenetic programmes in reciprocal coordination with the embryonic blastoderm, providing the necessary boundary conditions for development to proceed. In addition, the role of the enveloping cell layer in providing mechanical, osmotic and immunological protection during early stages of development, and the autonomous nutritional support provided by the yolk and YSL, have probably been key aspects that have enabled the massive radiation of teleosts to colonize every ecological niche on the Earth. This article is part of the theme issue 'Extraembryonic tissues: exploring concepts, definitions and functions across the animal kingdom'.

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.

Understanding the development of oral epithelial organs through single cell transcriptomic analysis

During craniofacial development, the oral epithelium begins as a morphologically homogeneous tissue that gives rise to locally complex structures, including the teeth, salivary glands and taste buds. How the epithelium is initially patterned and specified to generate diverse cell types remains largely unknown. To elucidate the genetic programs that direct the formation of distinct oral epithelial populations, we mapped the transcriptional landscape of embryonic day 12 mouse mandibular epithelia at single cell resolution. Our analysis identified key transcription factors and gene regulatory networks that define different epithelial cell types. By examining the spatiotemporal patterning process along the oral-aboral axis, our results propose a model in which the dental field is progressively confined to its position by the formation of the aboral epithelium anteriorly and the non-dental oral epithelium posteriorly. Using our data, we also identified Ntrk2 as a proliferation driver in the forming incisor, contributing to its invagination. Together, our results provide a detailed transcriptional atlas of the embryonic mandibular epithelium, and unveil new genetic markers and regulators that are present during the specification of various oral epithelial structures.

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.

Grainyhead-like (Grhl) Target Genes in Development and Cancer

Grainyhead-like (GRHL) factors are essential, highly conserved transcription factors (TFs) that regulate processes common to both natural cellular behaviours during embryogenesis, and de-regulation of growth and survival pathways in cancer. Serving to drive the transcription, and therefore activation of multiple co-ordinating pathways, the three GRHL family members (GRHL1-3) are a critical conduit for modulating the molecular landscape that guides cellular decision-making processes during proliferation, epithelial-mesenchymal transition (EMT) and migration. Animal models and in vitro approaches harbouring GRHL loss or gain-of-function are key research tools to understanding gene function, which gives confidence that resultant phenotypes and cellular behaviours may be translatable to humans. Critically, identifying and characterising the target genes to which these factors bind is also essential, as they allow us to discover and understand novel genetic pathways that could ultimately be used as targets for disease diagnosis, drug discovery and therapeutic strategies. GRHL1-3 and their transcriptional targets have been shown to drive comparable cellular processes in Drosophila, C. elegans, zebrafish and mice, and have recently also been implicated in the aetiology and/or progression of a number of human congenital disorders and cancers of epithelial origin. In this review, we will summarise the state of knowledge pertaining to the role of the GRHL family target genes in both development and cancer, primarily through understanding the genetic pathways transcriptionally regulated by these factors across disparate disease contexts.

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

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

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.

Uncovering the Pathogenesis of Orofacial Clefts Using Bioinformatics Analysis

OBJECTIVE

Many genes have been found to be associated with the occurrence of the orofacial clefts (OFC). The links between these pathogenic genes are rarely studied. In this study, bioinformatics analysis were performed in order to find associations between OFC-related genes and provide new ideas for etiology study of OFCs.

METHODS

Orofacial clefts-related genes were searched and identified from the Online Mendelian Inheritance of Man (OMIM.org). These genes were then analyzed by bioinformatics methods, including protein-protein interaction network, functional enrichment analysis, module analysis, and hub genes analysis.

RESULTS

After searching the database of OMIM.org and removing duplicate results, 279 genes were finally obtained. These genes were involved to 369 pathways in biological process, 56 in cell component, 64 in molecular function, and 45 in the Kyoto Encyclopedia of Genes and Genomes. Most identified genes were significantly enriched in embryonic appendage morphogenesis (29.17%), embryonic limb morphogenesis (6.06%), and limb development (4.33%) for biological process (Fig. 5A); ciliary tip (42.86%), MKS complex (28.57%), ciliary basal body (14.29%), and ciliary membrane (14.29%) for cell component. The top 10 hub genes were identified, including SHH, GLI2, PTCH1, SMAD4, FGFR1, BMP4, SOX9, SOX2, RUNX2, and CDH1.

CONCLUSIONS

Bioinformatics methods were used to analyze OFC-related genes in this study, including hub gene identifying and analysis, protein-protein interaction network construction, and functional enrichment analysis. Several potential mechanisms related to occurrence of OFCs were also discussed. These results may be helpful for further studies of the etiology of OFC.

Single-cell atlas of early chick development reveals gradual segregation of neural crest lineage from the neural plate border during neurulation

The epiblast of vertebrate embryos is comprised of neural and non-neural ectoderm, with the border territory at their intersection harboring neural crest and cranial placode progenitors. Here, we a generate single-cell atlas of the developing chick epiblast from late gastrulation through early neurulation stages to define transcriptional changes in the emerging ‘neural plate border’ as well as other regions of the epiblast. Focusing on the border territory, the results reveal gradual establishment of heterogeneous neural plate border signatures, including novel genes that we validate by fluorescent in situ hybridization. Developmental trajectory analysis infers that segregation of neural plate border lineages only commences at early neurulation, rather than at gastrulation as previously predicted. We find that cells expressing the prospective neural crest marker Pax7 contribute to multiple lineages, and a subset of premigratory neural crest cells shares a transcriptional signature with their border precursors. Together, our results suggest that cells at the neural plate border remain heterogeneous until early neurulation, at which time progenitors become progressively allocated toward defined neural crest and placode lineages. The data also can be mined to reveal changes throughout the developing epiblast.

Wholemount In-Situ Hybridization (WISH) in Zebrafish Embryos to Analyze Craniofacial Development

Wholemount in-situ hybridization in zebrafish is a powerful technique for visualizing spatiotemporal gene expression during development. Here we describe a technique to detect endogenous mRNA expression in zebrafish that can be adapted to use on embryos from the single-cell stage until 5 days postfertilization.

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.

Grhl3 promotes retention of epidermal cells under endocytic stress to maintain epidermal architecture in zebrafish

Epithelia such as epidermis cover large surfaces and are crucial for survival. Maintenance of tissue homeostasis by balancing cell proliferation, cell size, and cell extrusion ensures epidermal integrity. Although the mechanisms of cell extrusion are better understood, how epithelial cells that round up under developmental or perturbed genetic conditions are reintegrated in the epithelium to maintain homeostasis remains unclear. Here, we performed live imaging in zebrafish embryos to show that epidermal cells that round up due to membrane homeostasis defects in the absence of goosepimples/myosinVb (myoVb) function, are reintegrated into the epithelium. Transcriptome analysis and genetic interaction studies suggest that the transcription factor Grainyhead-like 3 (Grhl3) induces the retention of rounded cells by regulating E-cadherin levels. Moreover, Grhl3 facilitates the survival of MyoVb deficient embryos by regulating cell adhesion, cell retention, and epidermal architecture. Our analyses have unraveled a mechanism of retention of rounded cells and its importance in epithelial homeostasis.

The functional GRHL3-filaggrin axis maintains a tumor differentiation potential and influences drug sensitivity

Current therapies for treating heterogeneous cancers such as head and neck squamous cell carcinoma (HNSCC) are non-selective and are administered independent of response biomarkers. Therapy resistance subsequently emerges, resulting in increased cellular proliferation that is associated with loss of differentiation. Whether a cancer cell differentiation potential can dictate therapy responsiveness is still currently unknown. A multi-omic approach integrating whole-genome and whole-transcriptome sequencing with drug sensitivity was employed in a HNSCC mouse model, primary patients’ data, and human cell lines to assess the potential of functional differentiation in predicting therapy response. Interestingly, a subset of HNSCC with effective GRHL3-dependent differentiation was the most sensitive to inhibitors of PI3K/mTOR, c-Myc, and STAT3 signaling. Furthermore, we identified the GRHL3-differentiation target gene Filaggrin (FLG) as a response biomarker and more importantly, stratified HNSCC subsets as treatment resistant based on their FLG mutational profile. The loss of FLG in sensitive HNSCC resulted in a dramatic resistance to targeted therapies while the GRHL3-FLG signature predicted a favorable patient prognosis. This study provides evidence for a functional GRHL3-FLG tumor-specific differentiation axis that regulates targeted therapy response in HNSCC and establishes a rationale for clinical investigation of differentiation-paired targeted therapy in heterogeneous cancers.

Association analysis of SNPs in GRHL3, FAF1, and KCNJ2 with NSCPO sub-phenotypes in Han Chinese

OBJECTIVES

Non-syndromic cleft palate only (NSCPO) is a common congenital deformity with complex etiologies. GRHL3, FAF1, and KCNJ2 have been reported to be involved in the pathogenesis of NSCPO. Up till now, there have been no replication studies based on large Han Chinese. Therefore, this study aimed to investigate associations between GRHL3, FAF1, KCNJ2, and NSCPO sub-phenotypes patients in Han Chinese.

MATERIALS AND METHODS

Firstly, we selected 2 SNPs based on previous literatures: FAF1 (rs3827730) and GRHL3 (rs41268753). Also, we selected 8 tagSNPs in GRHL3 (rs557811, rs609352, rs10903078, rs6659209, rs12401714, rs12568599, rs3887581, rs12024148) and 2 tagSNPs in KCNJ2 (rs75855040 and rs236514). Afterward, we evaluated these SNPs among 1668 NSCPO patients and 1811 normal controls from Han Chinese. Following data were analyzed by PLINK and Haploview program.

RESULTS

Association analysis under additive model showed that allele A at rs12568599 in GRHL3 gene is significantly associated with NSCPO (p = 0.0034, OR = 1.38 and 95%CI: 1.11-1.72) and its sub-phenotype incomplete cleft palate (ICP) (p = 0.0039, OR = 1.4 and 95%CI: 1.11-1.75), and it could increase the risk of both NSCPO and ICP.

CONCLUSIONS

This study firstly found that rs12568599 in GRHL3 is associated with NSCPO and ICP in Han Chinese, indicating that sub-phenotypes of NSCPO have different etiologies.

Grainyhead-like transcription factors: guardians of the skin barrier

There has been selective pressure to maintain a skin barrier since terrestrial animals evolved 360 million years ago. These animals acquired an unique integumentary system with a keratinized, stratified, squamous epithelium surface barrier. The barrier protects against dehydration and entry of microbes and toxins. The skin barrier centres on the stratum corneum layer of the epidermis and consists of cornified envelopes cemented by the intercorneocyte lipid matrix. Multiple components of the barrier undergo cross-linking by transglutaminase (TGM) enzymes, while keratins provide additional mechanical strength. Cellular tight junctions also are crucial for barrier integrity. The grainyhead-like (GRHL) transcription factors regulate the formation and maintenance of the integument in diverse species. GRHL3 is essential for formation of the skin barrier during embryonic development, whereas GRHL1 maintains the skin barrier postnatally. This is achieved by transactivation of Tgm1 and Tgm5, respectively. In addition to its barrier function, GRHL3 plays key roles in wound repair and as an epidermal tumour suppressor. In its former role, GRHL3 activates the planar cell polarity signalling pathway to mediate wound healing by providing directional migration cues. In squamous epithelium, GRHL3 regulates the balance between proliferation and differentiation, and its loss induces squamous cell carcinoma (SCC). In the skin, this is mediated through increased expression of MIR21, which reduces the expression levels of GRHL3 and its direct target, PTEN, leading to activation of the PI3K-AKT signalling pathway. These data position the GRHL family as master regulators of epidermal homeostasis across a vast gulf of evolutionary history.

The power of zebrafish models for understanding the co-occurrence of craniofacial and limb disorders

Craniofacial and limb defects are two of the most common congenital anomalies in the general population. Interestingly, these defects are not mutually exclusive. Many patients with craniofacial phenotypes, such as orofacial clefting and craniosynostosis, also present with limb defects, including polydactyly, syndactyly, brachydactyly, or ectrodactyly. The gene regulatory networks governing craniofacial and limb development initially seem distinct from one another, and yet these birth defects frequently occur together. Both developmental processes are highly conserved among vertebrates, and zebrafish have emerged as an advantageous model due to their high fecundity, relative ease of genetic manipulation, and transparency during development. Here we summarize studies that have used zebrafish models to study human syndromes that present with both craniofacial and limb phenotypes. We discuss the highly conserved processes of craniofacial and limb/fin development and describe recent zebrafish studies that have explored the function of genes associated with human syndromes with phenotypes in both structures. We attempt to identify commonalities between the two to help explain why craniofacial and limb anomalies often occur together.

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.

Regulatory mechanisms governing epidermal stem cell function during development and homeostasis

Cell divisions and cell-fate decisions require stringent regulation for proper tissue development and homeostasis. The mammalian epidermis is a highly organized tissue structure that is sustained by epidermal stem cells (ESCs) that balance self-renewal and cell-fate decisions to establish a protective barrier, while replacing dying cells during homeostasis and in response to injury. Extensive work over past decades has provided insights into the regulatory mechanisms that control ESC specification, self-renewal and maintenance during different stages of the lifetime of an organism. In this Review, we discuss recent findings that have furthered our understanding of key regulatory features that allow ESCs to establish a functional barrier during development and to maintain tissue homeostasis in adults.

Genetics and signaling mechanisms of orofacial clefts

Craniofacial development involves several complex tissue movements including several fusion processes to form the frontonasal and maxillary structures, including the upper lip and palate. Each of these movements are controlled by many different factors that are tightly regulated by several integral morphogenetic signaling pathways. Subject to both genetic and environmental influences, interruption at nearly any stage can disrupt lip, nasal, or palate fusion and result in a cleft. Here, we discuss many of the genetic risk factors that may contribute to the presentation of orofacial clefts in patients, and several of the key signaling pathways and underlying cellular mechanisms that control lip and palate formation, as identified primarily through investigating equivalent processes in animal models, are examined.

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.

Non-syndromic Cleft Palate: An Overview on Human Genetic and Environmental Risk Factors

The epithelial and mesenchymal cells involved in early embryonic facial development are guided by complex regulatory mechanisms. Any factor perturbing the growth, approach and fusion of the frontonasal and maxillary processes could result in orofacial clefts that represent the most common craniofacial malformations in humans. The rarest and, probably for this reason, the least studied form of cleft involves only the secondary palate, which is posterior to the incisive foramen. The etiology of cleft palate only is multifactorial and involves both genetic and environmental risk factors. The intention of this review is to give the reader an overview of the efforts made by researchers to shed light on the underlying causes of this birth defect. Most of the scientific papers suggesting potential environmental and genetic causes of non-syndromic cleft palate are summarized in this review, including genome-wide association and gene–environment interaction studies.

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.

Insight Into the Function of RIPK4 in Keratinocyte Differentiation and Carcinogenesis

The receptor-interacting protein kinase 4 (RIPK4), a member of the RIPK family, was originally described as an interaction partner of protein kinase C (PKC) β and PKCδ. RIPK4 is identified as a key regulator of keratinocyte differentiation, cutaneous inflammation, and cutaneous wound repair. The mechanism by which RIPK4 integrates upstream signals to initiate specific responses remains elusive. Previous studies have indicated that RIPK4 can regulate several signaling pathways, including the NF-κB, Wnt/β-catenin, and RAF/MEK/ERK pathways. Furthermore, RIPK4-related biological signaling pathways interact with each other to form a complex network. Mounting evidence suggests that RIPK4 is aberrantly expressed in various kinds of cancers. In several types of squamous cell carcinoma (SCC), the mutations that drive aggressive SCC have been found in RIPK4. In addition, the function of RIPK4 in carcinogenesis is probably tissue-specific, since RIPK4 can play a dual role as both a tumor promoter and a tumor suppressor in different tumor types. Therefore, RIPK4 may represent as an independent prognostic factor and a promising novel therapeutic target, which can be used to identify the risks of patients and guide personalized treatments. In future, RIPK4-interacting pathways and precise molecular targets need to be investigated in order to further elucidate the mechanisms underlying epidermal differentiation and carcinogenesis.

Machine learning in prediction of genetic risk of nonsyndromic oral clefts in the Brazilian population

OBJECTIVES

Genetic variants in multiple genes and loci have been associated with the risk of nonsyndromic cleft lip with or without cleft palate (NSCL ± P). However, the estimation of risk remains challenge, because most of these variants are population-specific rendering the identification of the underlying genetic risk difficult. Herein we examined the use of machine learning network in previously reported single nucleotide polymorphisms (SNPs) to predict risk of NSCL ± P in the Brazilian population.

MATERIALS AND METHODS

Random forest and neural network methods were applied in 72 SNPs in a case-control sample composed by 722 NSCL ± P and 866 controls for discrimination of NSCL ± P risk. SNP-SNP interactions and functional annotation biological processes associated with the identified NSCL ± P risk genes were verified.

RESULTS

Supervised random forest decision trees revealed high scores of importance for the SNPs rs11717284 and rs1875735 in FGF12, rs41268753 in GRHL3, rs2236225 in MTHFD1, rs2274976 in MTHFR, rs2235371 and rs642961 in IRF6, rs17085106 in RHPN2, rs28372960 in TCOF1, rs7078160 in VAX1, rs10762573 and rs2131960 in VCL, and rs227731 in 17q22, with an accuracy of 99% and an error rate of approximately 3% to predict the risk of NSCL ± P. Those same 13 SNPs were considered the most important for the neural network to effectively predict NSCL ± P risk, with an overall accuracy of 94%. Multivariate regression model revealed significant interactions among all SNPs, with an exception of those in FGF12 and MTHFD1. The most significantly biological processes for selected genes were those involved in tissue and epithelium development; neural tube closure; and metabolism of methionine, folate, and homocysteine.

CONCLUSIONS

Our results provide novel clues for genetic mechanism studies of NSCL ± P and point out for a machine learning model composed by 13 SNPs that is capable of predicting NSCL ± P risk.

CLINICAL RELEVANCE

Although validation is necessary, this genetic panel can be useful in the near future to assist in NSCL ± P genetic counseling.

Functional Characterization of a Novel IRF6 Frameshift Mutation From a Van Der Woude Syndrome Family

Background

Loss-of-function mutations in interferon regulatory factor-6 (IRF6) are responsible for about 70% of cases of Van Der Woude Syndrome (VWS), an autosomal dominant developmental disorder characterized by pits and/or sinuses of the lower lip and cleft lip, cleft palate, or both.

Methods

We collected a Chinese Han VWS pedigree, performed sequencing and screening for the causal gene mutant. Initially, species conservation analysis and homology protein modeling were used to predict the potential pathogenicity of mutations. To test whether a VWS family-derived mutant variant of IRF6 retained function, we carried out rescue assays in irf6 maternal-null mutant zebrafish embryos. To assess protein stability, we overexpressed reference and family-variants of IRF6 in vitro.

Results

We focused on a VWS family that includes a son with bilateral lip pits, uvula fissa and his father with bilateral cleft lip and palate. After sequencing and screening, a frameshift mutation of IRF6 was identified as the potential causal variant (NM.006147.3, c.1088-1091delTCTA; p.Ile363ArgfsTer33). The residues in this position are strongly conserved among species and homology modeling suggests the variant alters the protein structure. In irf6 maternal-null mutant zebrafish embryos the periderm differentiates abnormally and the embryos rupture and die during gastrulation. Injection of mRNA encoding the reference variant of human IRF6, but not of the frame-shift variant, rescued such embryos through gastrulation. Upon overexpression in HEK293FT cells, the IRF6 frame-shift mutant was relatively unstable and was preferentially targeted to the proteasome in comparison to the reference variant.

Conclusion

In this VWS pedigree, a novel frameshift of IRF6 was identified as the likely causative gene variant. It is a lost function mutation which could not rescue abnormal periderm phenotype in irf6 maternal-null zebrafish and which causes the protein be unstable through proteasome-dependent degradation.

SPECC1L regulates palate development downstream of IRF6

SPECC1L mutations have been identified in patients with rare atypical orofacial clefts and with syndromic cleft lip and/or palate (CL/P). These mutations cluster in the second coiled-coil and calponin homology domains of SPECC1L and severely affect the ability of SPECC1L to associate with microtubules. We previously showed that gene-trap knockout of Specc1l in mouse results in early embryonic lethality. We now present a truncation mutant mouse allele, Specc1lΔC510, that results in perinatal lethality. Specc1lΔC510/ΔC510 homozygotes showed abnormal palate rugae but did not show cleft palate. However, when crossed with a gene-trap allele, Specc1lcGT/ΔC510 compound heterozygotes showed a palate elevation delay with incompletely penetrant cleft palate. Specc1lcGT/ΔC510 embryos exhibit transient oral epithelial adhesions at E13.5, which may delay shelf elevation. Consistent with oral adhesions, we show periderm layer abnormalities, including ectopic apical expression of adherens junction markers, similar to Irf6 hypomorphic mutants and Arhgap29 heterozygotes. Indeed, SPECC1L expression is drastically reduced in Irf6 mutant palatal shelves. Finally, we wanted to determine if SPECC1L deficiency also contributed to non-syndromic (ns) CL/P. We sequenced 62 Caucasian, 89 Filipino, 90 Ethiopian, 90 Nigerian and 95 Japanese patients with nsCL/P and identified three rare coding variants (p.Ala86Thr, p.Met91Iso and p.Arg546Gln) in six individuals. These variants reside outside of SPECC1L coiled-coil domains and result in milder functional defects than variants associated with syndromic clefting. Together, our data indicate that palate elevation is sensitive to deficiency of SPECC1L dosage and function and that SPECC1L cytoskeletal protein functions downstream of IRF6 in palatogenesis.

Interrogating the Grainyhead-like 2 (Grhl2) genomic locus identifies an enhancer element that regulates palatogenesis in mouse

The highly-conserved Grainyhead-like (Grhl) transcription factors are critical regulators of embryogenesis that regulate cellular survival, proliferation, migration and epithelial integrity, especially during the formation of the craniofacial skeleton. Family member Grhl2 is expressed throughout epithelial tissues during development, and loss of Grhl2 function leads to significant defects in neurulation, abdominal wall closure, formation of the face and fusion of the maxilla/palate. Whereas numerous downstream target genes of Grhl2 have been identified, very little is known about how this crucial developmental transcription factor itself is regulated. Here, using in silico and in utero expression analyses and functional deletion in mice, we have identified a novel 2.4 ​kb enhancer element (mm1286) that drives reporter gene expression in a pattern that strongly recapitulates endogenous Grhl2 in the craniofacial primordia, modulates Grhl2 expression in these tissues, and augments Grhl2-mediated closure of the secondary palate. Deletion of this genomic element, in the context of inactivation of one allele of Grhl2 (through generation of double heterozygous Grhl2^+/-;mm1286^+/- mice), results in a significant predisposition to palatal clefting at birth. Moreover, we found that a highly conserved 325 bp region of mm1286 is both necessary and sufficient for mediating the craniofacial-specific enhancer activity of this region, and that an extremely well-conserved 12-bp sequence within this element (CTGTCAAACAGGT) substantially determines full enhancer function. Together, these data provide valuable new insights into the upstream genomic regulatory landscape responsible for transcriptional control of Grhl2 during palatal closure.

Epithelial Migration and Non-adhesive Periderm Are Required for Digit Separation during Mammalian Development

The fusion of digits or toes, syndactyly, can be part of complex syndromes, including van der Woude syndrome. A subset of van der Woude cases is caused by dominant-negative mutations in the epithelial transcription factor Grainyhead like-3 (GRHL3), and Grhl3^-/-mice have soft-tissue syndactyly. Although impaired interdigital cell death of mesenchymal cells causes syndactyly in multiple genetic mutants, Grhl3^-/- embryos had normal interdigital cell death, suggesting alternative mechanisms for syndactyly. We found that in digit separation, the overlying epidermis forms a migrating interdigital epithelial tongue (IET) when the epithelium invaginates to separate the digits. Normally, the non-adhesive surface periderm allows the IET to bifurcate as the digits separate. In contrast, in Grhl3^-/- embryos, the IET moves normally between the digits but fails to bifurcate because of abnormal adhesion of the periderm. Our study identifies epidermal developmental processes required for digit separation.

A Novel Combined Scientific and Artistic Approach for the Advanced Characterization of Interactomes: The Akirin/Subolesin Model

The main objective of this study was to propose a novel methodology to approach challenges in molecular biology. Akirin/Subolesin (AKR/SUB) are vaccine protective antigens and are a model for the study of the interactome due to its conserved function in the regulation of different biological processes such as immunity and development throughout the metazoan. Herein, three visual artists and a music professor collaborated with scientists for the functional characterization of the AKR2 interactome in the regulation of the NF-κB pathway in human placenta cells. The results served as a methodological proof-of-concept to advance this research area. The results showed new perspectives on unexplored characteristics of AKR2 with functional implications. These results included protein dimerization, the physical interactions with different proteins simultaneously to regulate various biological processes defined by cell type-specific AKR–protein interactions, and how these interactions positively or negatively regulate the nuclear factor kappa-light-chain-enhancer of activated B cells (NF-κB) signaling pathway in a biological context-dependent manner. These results suggested that AKR2-interacting proteins might constitute suitable secondary transcription factors for cell- and stimulus-specific regulation of NF-κB. Musical perspective supported AKR/SUB evolutionary conservation in different species and provided new mechanistic insights into the AKR2 interactome. The combined scientific and artistic perspectives resulted in a multidisciplinary approach, advancing our knowledge on AKR/SUB interactome, and provided new insights into the function of AKR2–protein interactions in the regulation of the NF-κB pathway. Additionally, herein we proposed an algorithm for quantum vaccinomics by focusing on the model proteins AKR/SUB.