MSX1 and TGF-beta3 are novel target genes functionally regulated by FOXE1

Clinical and Genetic Correlation in Neurocristopathies: Bridging a Precision Medicine Gap

Neurocristopathies (NCPs) encompass a spectrum of disorders arising from issues during the formation and migration of neural crest cells (NCCs). NCCs undergo epithelial-mesenchymal transition (EMT) and upon key developmental gene deregulation, fetuses and neonates are prone to exhibit diverse manifestations depending on the affected area. These conditions are generally rare and often have a genetic basis, with many following Mendelian inheritance patterns, thus making them perfect candidates for precision medicine. Examples include cranial NCPs, like Goldenhar syndrome and Axenfeld-Rieger syndrome; cardiac-vagal NCPs, such as DiGeorge syndrome; truncal NCPs, like congenital central hypoventilation syndrome and Waardenburg syndrome; and enteric NCPs, such as Hirschsprung disease. Additionally, NCCs' migratory and differentiating nature makes their derivatives prone to tumors, with various cancer types categorized based on their NCC origin. Representative examples include schwannomas and pheochromocytomas. This review summarizes current knowledge of diseases arising from defects in NCCs' specification and highlights the potential of precision medicine to remedy a clinical phenotype by targeting the genotype, particularly important given that those affected are primarily infants and young children.

Canonical Wnt signaling is not required for Tgfb3 expression in the basal medial edge epithelium during palatogenesis

The secondary palate forms from two lateral primordia called the palatal shelves which form a contact in the midline, become adherent at the fusing interface (medial edge epithelia, MEE) and subsequently fuse. The gene encoding transforming growth factor-ß3 (Tgfb3) is strongly and specifically expressed in MEE cells. Our previous study suggested that Tgfb3 expression is controlled via upstream cis-regulatory elements in and around the neighboring Ift43 gene. Another study suggested that the canonical Wnt signaling via ß-Catenin is responsible for the MEE-specific Tgfb3 gene expression, since deletion of the Ctnnb1 gene by a commonly used Keratin 14-Cre (K14Cre) mouse line almost completely abolished Tgfb3 expression in the MEE resulting in cleft palate. Here, we wanted to analyze whether Tcf/Lef consensus binding sites located in the previously identified regions of the Ift43 gene are responsible for the spatiotemporal control of Tgfb3 expression during palatogenesis. We show that contrary to the previous report, deletion of the Ctnnb1 gene in basal MEE cells by the K14Cre driver (the same K14Cre mouse line was used as in the previous study referenced above) does not affect the MEE-specific Tgfb3 expression or TGFß3-dependent palatal epithelial fusion. All mutant embryos showed a lack of palatal rugae accompanied by other craniofacial defects, e.g., a narrow snout and a small upper lip, while only a small subset (<5%) of Ctnnb1 mutants displayed a cleft palate. Moreover, the K14Cre:Ctnnb1 embryos showed reduced levels and altered patterns of Shh expression. Our present data imply that epithelial ß-catenin may not be required for MEE-specific Tgfb3 expression or palatal epithelial fusion.

Disruption of the foxe1 gene in zebrafish reveals conserved functions in development of the craniofacial skeleton and the thyroid

Introduction: Mutations in the FOXE1 gene are implicated in cleft palate and thyroid dysgenesis in humans.

Methods: To investigate whether zebrafish could provide meaningful insights into the etiology of developmental defects in humans related to FOXE1, we generated a zebrafish mutant that has a disruption in the nuclear localization signal in the foxe1 gene, thereby restraining nuclear access of the transcription factor. We characterized skeletal development and thyroidogenesis in these mutants, focusing on embryonic and larval stages.

Results: Mutant larvae showed aberrant skeletal phenotypes in the ceratohyal cartilage and had reduced whole body levels of Ca, Mg and P, indicating a critical role for foxe1 in early skeletal development. Markers of bone and cartilage (precursor) cells were differentially expressed in mutants in post-migratory cranial neural crest cells in the pharyngeal arch at 1 dpf, at induction of chondrogenesis at 3 dpf and at the start of endochondral bone formation at 6 dpf. Foxe1 protein was detected in differentiated thyroid follicles, suggesting a role for the transcription factor in thyroidogenesis, but thyroid follicle morphology or differentiation were unaffected in mutants.

Discussion: Taken together, our findings highlight the conserved role of Foxe1 in skeletal development and thyroidogenesis, and show differential signaling of osteogenic and chondrogenic genes related to foxe1 mutation.

Data on cardiac lncRNA STX18-AS1 expression in developing human hearts and function during in vitro hESC-cardiomyocyte differentiation

This article presents data concerning STX18-AS1, a long noncoding RNA gene identified from a Genome-wide association study of Atrial Septal Defect (ASD). The data describes its expression patterns in human tissues and functions in regulating cardiomyocyte differentiation in vitro. STX18-AS1 is a lncRNA with a higher abundance in developing tissues, including hearts. Its transcription distribution within the embryonic hearts during key heart septation stages supports STX18-AS1's association with risk SNPs for ASD. The CRISPR stem cell pool in which STX18-AS1 was knocked down, showed reduced CM differentiation efficiency and lower expression of key cardiac transcriptional factors. This indicated its regulative role in supporting the lineage specification from cardiac mesoderm into cardiac progenitors and cardiomyocytes. These data can benefit the understanding of human embryonic heart developmental biology, and the time-course changes of cardiac transcriptional factors during in vitro cardiomyocyte differentiation from human embryonic stem cells.

A new FOXE1 homozygous frameshift variant expands the genotypic and phenotypic spectrum of Bamforth-Lazarus syndrome

Bamforth-Lazarus syndrome is a rare autosomal recessive disease caused by biallelic loss-of-function variants in the FOXE1 gene. The condition is characterized by congenital hypothyroidism due to thyroid agenesis or thyroid hypoplasia, cleft palate, spiky hair, with or without choanal atresia, and bifid epiglottis. To date, seven pathogenic variants have been reported in the FOXE1 gene causing Bamforth-Lazarus syndrome. Here we report a novel homozygous loss-of-function variant in the FOXE1 gene NM_004473.4:c.141dupC:p.(Leu49Profs*75) leading to congenital hypothyroidism due to thyroid agenesis, scalp hair abnormalities, cleft palate, small areola, cafe-au-lait spots, mild bilateral hearing loss, skin abnormalities, and facial dysmorphism. We describe the evolving phenotype in the patient with age and review previous variants reported in FOXE1. This report further expands the clinical and molecular spectrum of Bamforth-Lazarus syndrome.

Novel insights from a multiomics dissection of the hayflick limit

The process wherein dividing cells exhaust proliferative capacity and enter into replicative senescence has become a prominent model for cellular aging in vitro. Despite decades of study, this cellular state is not fully understood in culture and even much less so during aging. Here, we revisit Leonard Hayflick’s original observation of replicative senescence in WI-38 human lung fibroblasts equipped with a battery of modern techniques including RNA-seq, single-cell RNA-seq, proteomics, metabolomics, and ATAC-seq. We find evidence that the transition to a senescent state manifests early, increases gradually, and corresponds to a concomitant global increase in DNA accessibility in nucleolar and lamin associated domains. Furthermore, we demonstrate that senescent WI-38 cells acquire a striking resemblance to myofibroblasts in a process similar to the epithelial to mesenchymal transition (EMT) that is regulated by t YAP1/TEAD1 and TGF-β2. Lastly, we show that verteporfin inhibition of YAP1/TEAD1 activity in aged WI-38 cells robustly attenuates this gene expression program.

Expression Analysis of FGF/FGFR and FOX Family Proteins in Mucosal Tissue Obtained from Orofacial Cleft-Affected Children

Orofacial clefts affect hundreds of thousands of children worldwide annually and are usually corrected by a series of surgeries extending to childhood. The underlying mechanisms that lead to clefts are still unknown, mainly because of the multifactorial etiology and the myriad of interactions between genes and environmental factors. In the present study, we investigated the role and expression of candidate genes belonging to the FGF/FGFR signaling pathway and FOX family in tissue material obtained from 12 pediatric patients undergoing cleft correction surgery. The expression was investigated using immunohistochemistry (IHC) and chromogenic in-situ hybridization (CISH) in three cell/tissue types-epithelial cells, connective tissue, and endothelial cells. We found elevated expression of FGFR1 in epithelial cells while no expression was observed in endothelial cells. Further, our results elucidate the potential pathogenetic role of FGFR1 in cellular proliferation, local site inflammation, and fibrosis in cleft patients. Along with bFGF (also called FGF2), FGFR1 could play a pro-inflammatory role in clefts. Over-amplification of FGFR2 in some patients, along with bFGF, could potentially suggest roles for these genes in angiogenesis. Additionally, increased expression of FOXE1 (also called TTF2) contributes to local site inflammation. Finally, zero to low amplification of FOXO1 could suggest its potential role in inducing oxidative stress in the endothelium along with reduced epithelial apoptosis.

Cleft Candidate Genes and Their Products in Human Unilateral Cleft Lip Tissue

Cleft lip and palate are common congenital pathologies that affect the human population worldwide. The formation of cleft lip is associated with multiple genes and their coded proteins, which regulate the development of craniofacial region, but the exact role of these factors is not always clear. The use of morphological studies for evaluation of human cleft-affected tissue has been limited because of insufficiency of available pathological material. The aim of this study was to detect and compare the immunohistochemical expression of cleft candidate gene coded proteins (DLX4, MSX2, HOXB3, SHH, PAX7, SOX3, WNT3A, and FOXE1) in the non-syndromic unilateral cleft lip patient tissue and control group tissue. A semiquantitative counting method was used to evaluate the tissue in biotin-streptavidin-stained slides. Statistically significant differences between the patient and control groups were found for the number of immunoreactive structures for SHH (p = 0.019) and FOXE1 (p = 0.011) in the connective tissue and SOX3 (p = 0.012) in the epithelium. Multiple statistically significant very strong and strong correlations were found between the immunoreactives in cleft-affected tissue. These significant differences and various correlations indicate that multiple morphopathogenetic pathways are possibly involved in unilateral cleft lip pathogenesis. Therefore, we further discuss these possible interactions.

Orofacial clefts embryology, classification, epidemiology, and genetics

Orofacial clefts (OFCs) rank as the second most common congenital birth defect in the United States after Down syndrome and are the most common head and neck congenital malformations. They are classified as cleft lip with or without cleft palate (CL/P) and cleft palate only (CPO). OFCs have significant psychological and socio-economic impact on patients and their families and require a multidisciplinary approach for management and counseling. A complex interaction between genetic and environmental factors contributes to the incidence and clinical presentation of OFCs. In this comprehensive review, the embryology, classification, epidemiology and etiology of clefts are thoroughly discussed and a "state-of-the-art" snapshot of the recent advances in the genetics of OFCs is presented.

Unraveling the Complex Interplay Between Transcription Factors and Signaling Molecules in Thyroid Differentiation and Function, From Embryos to Adults

Thyroid differentiation of progenitor cells occurs during embryonic development and in the adult thyroid gland, and the molecular bases of these complex and finely regulated processes are becoming ever more clear. In this Review, we describe the most recent advances in the study of transcription factors, signaling molecules and regulatory pathways controlling thyroid differentiation and development in the mammalian embryo. We also discuss the maintenance of the adult differentiated phenotype to ensure the biosynthesis of thyroid hormones. We will focus on endoderm-derived thyroid epithelial cells, which are responsible for the formation of the thyroid follicle, the functional unit of the thyroid gland. The use of animal models and pluripotent stem cells has greatly aided in providing clues to the complicated puzzle of thyroid development and function in adults. The so-called thyroid transcription factors - Nkx2-1, Foxe1, Pax8 and Hhex - were the first pieces of the puzzle identified in mice. Other transcription factors, either acting upstream of or directly with the thyroid transcription factors, were subsequently identified to, almost, complete the puzzle. Among them, the transcription factors Glis3, Sox9 and the cofactor of the Hippo pathway Taz, have emerged as important players in thyroid differentiation and development. The involvement of signaling molecules increases the complexity of the puzzle. In this context, the importance of Bmps, Fgfs and Shh signaling at the onset of development, and of TSH, IGF1 and TGFβ both at the end of terminal differentiation in embryos and in the adult thyroid, are well recognized. All of these aspects are covered herein. Thus, readers will be able to visualize the puzzle of thyroid differentiation with most - if not all - of the pieces in place.

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.

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.

Transcriptional landscape of the embryonic chicken Müllerian duct

Background

Müllerian ducts are paired embryonic tubes that give rise to the female reproductive tract in vertebrates. Many disorders of female reproduction can be attributed to anomalies of Müllerian duct development. However, the molecular genetics of Müllerian duct formation is poorly understood and most disorders of duct development have unknown etiology. In this study, we describe for the first time the transcriptional landscape of the embryonic Müllerian duct, using the chicken embryo as a model system. RNA sequencing was conducted at 1 day intervals during duct formation to identify developmentally-regulated genes, validated by in situ hybridization.

Results

This analysis detected hundreds of genes specifically up-regulated during duct morphogenesis. Gene ontology and pathway analysis revealed enrichment for developmental pathways associated with cell adhesion, cell migration and proliferation, ERK and WNT signaling, and, interestingly, axonal guidance. The latter included factors linked to neuronal cell migration or axonal outgrowth, such as Ephrin B2, netrin receptor, SLIT1 and class A semaphorins. A number of transcriptional modules were identified that centred around key hub genes specifying matrix-associated signaling factors; SPOCK1, HTRA3 and ADGRD1. Several novel regulators of the WNT and TFG-β signaling pathway were identified in Müllerian ducts, including APCDD1 and DKK1, BMP3 and TGFBI. A number of novel transcription factors were also identified, including OSR1, FOXE1, PRICKLE1, TSHZ3 and SMARCA2. In addition, over 100 long non-coding RNAs (lncRNAs) were expressed during duct formation.

Conclusions

This study provides a rich resource of new candidate genes for Müllerian duct development and its disorders. It also sheds light on the molecular pathways engaged during tubulogenesis, a fundamental process in embryonic development.

Association between forkhead box E1 polymorphisms and risk of non-syndromic cleft lip with or without cleft palate: A meta-analysis

OBJECTIVE

The purpose of the present work was to investigate the association between forkhead box E1 (FOXE1) and the risk of non-syndromic cleft lip with or without cleft palate (NSCL/P).

MATERIALS AND METHODS

Relevant studies were searched in several professional databases up to 31 July 2019. The pooled odds ratios (ORs) and 95% confidence intervals (95% CIs) were calculated using a fixed-effect model or a random-effect model to analyse the relationship between FOXE1 polymorphisms and NSCL/P.

RESULTS

A total of four single nucleotide polymorphisms (SNPs), including rs3758249, rs4460498, rs1443434 and rs10217225, were analysed. The overall findings showed that FOXE1 rs4460498 was statistically associated with NSCL/P (including cleft lip with or without cleft palate (CL/P) and cleft palate only (CPO)). Genotypes CC and CT of rs4460498 were significantly more closely correlated with NSCL/P (including CL/P and CPO) than genotype TT (NSCL/P: TT vs CC, OR = 0.630, P = .000; TT vs TC + CC, OR = 0.775, P = .020; CL/P: TT vs CC, OR = 0.664, P = .000; TT vs TC + CC, OR = 0.738, P = .006. CPO: TT vs CC, OR = 0.761, P = .027; TT vs TC + CC, OR = 0.792, P = .045). For rs10217225, only the TT genotype might have contributed to the elevated risk of CL/P (TT vs CC OR = 2.236, P = .000). The other FOXE1 polymorphisms were not associated with NSCLP, CL/P or CPO.

CONCLUSION

The meta-analysis provided confirmation that the polymorphism of FOXE1 rs10217225 was correlated with an increased risk of CL/P, and the polymorphism of FOXE1 rs4460498 was a protective factor for NSCL/P, including CLP and CPO.

A systematic review and meta-analysis on protective role of forkhead box E1 (FOXE1) polymorphisms in susceptibility to non-syndromic cleft lip/palate

INTRODUCTION

Several environmental and genetic factors have a role in the aetiology of non-syndromic cleft lip/palate (NSCL/P). This meta-analysis evaluated the association of rs3758249 and rs4460498 forkhead box E1 (FOXE1) polymorphisms with the NSCL/P risk.

MATERIALS AND METHODS

The Scopus, Cochrane Library, Web of Science, and PubMed databases were searched for articles published until March 2019. The analyses were performed by Review Manager 5.3 using the crude odds ratio (OR) and 95% confidence interval (CI) for a strong association between FOXE1 polymorphisms and the risk of NSCL/P.

RESULTS

Out of 161 articles retrieved from the databases, four case-control articles were involved in the meta-analysis. The pooled ORs of rs4460498 polymorphism based on allelic, homozygous, heterozygous, dominant, and recessive models were 0.74 (95% CI: 0.69, 0.80; P<0.00001), 0.43 (95% CI: 0.30, 0.61; P<0.00001), 0.66 (95% CI: 0.55, 0.80; P<0.0001), 0.66 (95% CI: 0.59, 0.73; P<0.00001), and 0.70 (95% CI: 0.60, 0.82; P<0.0001), respectively; whereas, the pooled OR of rs3758249 polymorphism were 0.86 (95% CI: 0.71, 1.04; P=0.12), 0.68 (95% CI: 0.57, 0.82; P<0.0001), 0.79 (95% CI: 0.57, 1.09; P=0.15), 0.79 (95% CI: 0.58, 1.08; P=0.14), and 0.80 (95% CI: 0.68, 0.95; P=0.010) for the afore-mentioned models, respectively.

CONCLUSIONS

The results showed that the T allele, TT, and CT genotypes of rs4460498 polymorphism were significantly associated with a decreased risk of NSCL/P; whereas, for rs3758249 polymorphism, only the AA genotype had a significant protective role in NSCL/P. Thus, FOXE1 is strongly associated with NSCL/P in the populations.

Development of the thyroid gland

Thyroid hormones are crucial for organismal development and homeostasis. In humans, untreated congenital hypothyroidism due to thyroid agenesis inevitably leads to cretinism, which comprises irreversible brain dysfunction and dwarfism. Elucidating how the thyroid gland - the only source of thyroid hormones in the body - develops is thus key for understanding and treating thyroid dysgenesis, and for generating thyroid cells in vitro that might be used for cell-based therapies. Here, we review the principal mechanisms involved in thyroid organogenesis and functional differentiation, highlighting how the thyroid forerunner evolved from the endostyle in protochordates to the endocrine gland found in vertebrates. New findings on the specification and fate decisions of thyroid progenitors, and the morphogenesis of precursor cells into hormone-producing follicular units, are also discussed.

MYH9 binds to lncRNA gene PTCSC2 and regulates FOXE1 in the 9q22 thyroid cancer risk locus

A locus on chromosome 9q22 harbors a SNP (rs965513) firmly associated with risk of papillary thyroid carcinoma (PTC). The locus also comprises the forkhead box E1 (FOXE1) gene, which is implicated in thyroid development, and a long noncoding RNA (lncRNA) gene, papillary thyroid cancer susceptibility candidate 2 (PTCSC2). How these might interact is not known. Here we report that PTCSC2 binds myosin-9 (MYH9). In a bidirectional promoter shared by FOXE1 and PTCSC2, MYH9 inhibits the promoter activity in both directions. This inhibition can be reversed by PTCSC2, which acts as a suppressor. RNA knockdown of FOXE1 in primary thyroid cells profoundly interferes with the p53 pathway. We propose that the interaction between the lncRNA, its binding protein MYH9, and the coding gene FOXE1 underlies the predisposition to PTC triggered by rs965513.

Determinants of orofacial clefting I: Effects of 5-Aza-2'-deoxycytidine on cellular processes and gene expression during development of the first branchial arch

In this study, we identify gene targets and cellular events mediating the teratogenic action(s) of 5-Aza-2'-deoxycytidine (AzaD), an inhibitor of DNA methylation, on secondary palate development. Exposure of pregnant mice (on gestation day (GD) 9.5) to AzaD for 12h resulted in the complete penetrance of cleft palate (CP) in fetuses. Analysis of cells of the embryonic first branchial arch (1-BA), in fetuses exposed to AzaD, revealed: 1) significant alteration in expression of genes encoding several morphogenetic factors, cell cycle inhibitors and regulators of apoptosis; 2) a decrease in cell proliferation; and, 3) an increase in apoptosis. Pyrosequencing of selected genes, displaying pronounced differential expression in AzaD-exposed 1-BAs, failed to reveal significant alterations in CpG methylation levels in their putative promoters or gene bodies. CpG methylation analysis suggested that the effects of AzaD on gene expression were likely indirect.

Loss of polarity protein AF6 promotes pancreatic cancer metastasis by inducing Snail expression

Pancreatic cancer (PC) is a particularly lethal form of cancer with high potential for metastasis to distant organs. Disruption of cell polarity is a hallmark of advanced epithelial tumours. Here we show that the polarity protein AF6 (afadin and MLLT4) is expressed at low levels in PC. We demonstrate that depletion of AF6 markedly promotes proliferation and metastasis of PC cells through upregulation of the expression of Snail protein, and this requires the nuclear localization of AF6. Furthermore, AF6 deficiency in PC cells leads to increased formation of a Dishevelled 2 (Dvl2)-FOXE1 complex on the promoter region of Snail gene, and activation of Snail expression. Altogether, our data established AF6 as a potential inhibitor of metastasis in PC cells. Targeting the Dvl2-FOXE1-Snail signalling axis may thus represent a promising therapeutic strategy.

Thyroid transcription factors in development, differentiation and disease

Identification of the thyroid transcription factors (TTFs), NKX2-1, FOXE1, PAX8 and HHEX, has considerably advanced our understanding of thyroid development, congenital thyroid disorders and thyroid cancer. The TTFs are fundamental to proper formation of the thyroid gland and for maintaining the functional differentiated state of the adult thyroid; however, they are not individually required for precursor cell commitment to a thyroid fate. Although knowledge of the mechanisms involved in thyroid development has increased, the full complement of genes involved in thyroid gland specification and the signals that trigger expression of the genes that encode the TTFs remain unknown. The mechanisms involved in thyroid organogenesis and differentiation have provided clues to identifying the genes that are involved in human congenital thyroid disorders and thyroid cancer. Mutations in the genes that encode the TTFs, as well as polymorphisms and epigenetic modifications, have been associated with thyroid pathologies. Here, we summarize the roles of the TTFs in thyroid development and the mechanisms by which they regulate expression of the genes involved in thyroid differentiation. We also address the implications of mutations in TTFs in thyroid diseases and in diseases not related to the thyroid gland.

Palatogenesis and cutaneous repair: A two-headed coin

BACKGROUND

The reparative mechanism that operates following post-natal cutaneous injury is a fundamental survival function that requires a well-orchestrated series of molecular and cellular events. At the end, the body will have closed the hole using processes like cellular proliferation, migration, differentiation and fusion.

RESULTS

These processes are similar to those occurring during embryogenesis and tissue morphogenesis. Palatogenesis, the formation of the palate from two independent palatal shelves growing towards each other and fusing, intuitively, shares many similarities with the closure of a cutaneous wound from the two migrating epithelial fronts.

CONCLUSIONS

In this review, we summarize the current information on cutaneous development, wound healing, palatogenesis and orofacial clefting and propose that orofacial clefting and wound healing are conserved processes that share common pathways and gene regulatory networks.

Role of angiogenesis-related genes in cleft lip/palate: review of the literature

OBJECTIVES

Cleft lip and cleft palate (CLP) are the most common congenital craniofacial anomalies. They have a multifactorial etiology and result from an incomplete fusion of the facial buds. Two main mechanisms, acting alone or interacting with each other, were evidenced in this fusion defect responsible for CLP: defective tissue development and/or defective apoptosis in normal or defective tissues. The objective of this work was to study the implication and role of angiogenesis-related genes in the etiology of CL/P.

METHODS

Our methodological approach included a systematic and thorough analysis of the genes involved in CL/P (syndromic and non-syndromic forms) including previously identified genes but also genes that could potentially be angiogenesis-related (OMIM, Pub Med).We studied the interactions of these different genes and their relationships with potential environmental factors.

RESULTS

TGFβ, FGA, PDGFc, PDGFRa, FGF, FGFR1, FGFR2 growth factors as well as MMP and TIMP2 proteolytic enzymes are involved in the genesis of CLP (P>L). Furthermore, 18 genes involved in CLP also interact with angiogenesis-related genes.

DISCUSSION

Even if the main angiogenesis-related genes involved in CLP formation are genes participating in several biological activities and their implication might not be always related to angiogenesis defects, they nevertheless remain an undeniably important research pathway. Furthermore, their interactions with environmental factors make them good candidates in the field of CLP prevention.

Control elements targeting Tgfb3 expression to the palatal epithelium are located intergenically and in introns of the upstream Ift43 gene

Tgfb3 is strongly and specifically expressed in the epithelial tips of pre-fusion palatal shelves where it plays a critical non-redundant role in palatal fusion in both medial edge epithelial (MEE) cells and in a thin layer of flattened peridermal cells that covers the MEE. It is not known how Tgfb3 expression is regulated in these specific cell types. Using comparative genomics and transgenic reporter assays, we have identified cis-regulatory elements that could control Tgfb3 expression during palatogenesis. Our results show that a 61-kb genomic fragment encompassing the Tgfb3 gene drives remarkably specific reporter expression in the MEE and adjacent periderm. Within this fragment, we identified two small, non-coding, evolutionarily conserved regions in intron 2 of the neighboring Ift43 gene, and a larger region in the intervening sequence between the Ift43 and Tgfb3 genes, each of which could target reporter activity to the tips of pre-fusion/fusing palatal shelves. Identification of the cis-regulatory sequences controlling spatio-temporal Tgfb3 expression in palatal shelves is a key step toward understanding upstream regulation of Tgfb3 expression during palatogenesis and should enable the development of improved tools to investigate palatal epithelial fusion.

Exploratory genotype-phenotype correlations of facial form and asymmetry in unaffected relatives of children with non-syndromic cleft lip and/or palate

Family relatives of children with nonsyndromic cleft lip with or without cleft palate (NSCL/P) who presumably carry a genetic risk yet do not manifest overt oral clefts, often present with distinct facial morphology of unknown genetic etiology. This study investigates distinct facial morphology among unaffected relatives and examines whether candidate genes previously associated with overt NSCL/P and left-right body patterning are correlated with such facial morphology. Cases were unaffected relatives of individuals with NSCL/P (n = 188) and controls (n = 194) were individuals without family history of NSCL/P. Cases and controls were genotyped for 20 SNPs across 13 candidate genes for NSCL/P (PAX7, ABCA4-ARHGAP29, IRF6, MSX1, PITX2, 8q24, FOXE1, TGFB3 and MAFB) and left-right body patterning (LEFTY1, LEFTY2, ISL1 and SNAI1). Facial shape and asymmetry phenotypes were obtained via principal component analyses and Procrustes analysis of variance from 32 coordinate landmarks, digitized on 3D facial images. Case-control comparisons of phenotypes obtained were performed via multivariate regression adjusting for age and gender. Phenotypes that differed significantly (P < 0.05) between cases and controls were regressed on the SNPs one at a time. Cases had significantly (P < 0.05) more profile concavity with upper face retrusion, upturned noses with obtuse nasolabial angles, more protrusive chins, increased lower facial heights, thinner and more retrusive lips and more protrusive foreheads. Furthermore, cases showed significantly more directional asymmetry compared to controls. Several of these phenotypes were significantly associated with genetic variants (P < 0.05). Facial height and width were associated with SNAI1. Midface antero-posterior (AP) projection was associated with LEFTY1. The AP position of the chin was related to SNAI1, IRF6, MSX1 and MAFB. The AP position of the forehead and the width of the mouth were associated with ABCA4-ARHGAP29 and MAFB. Lastly, facial asymmetry was related to LEFTY1, LEFTY2 and SNAI1. This study demonstrates that, genes underlying lip and palate formation and left-right patterning also contribute to facial features characteristic of the NSCL/P spectrum.

Signaling networks in palate development

Palatogenesis, the formation of the palate, is a dynamic process regulated by a complex series of context-dependent morphogenetic signaling events. Many genes involved in palatogenesis have been discovered through the use of genetically manipulated mouse models as well as from human genetic studies, but the roles of these genes and their products in signaling networks regulating palatogenesis are still poorly known. In this review, we give a brief overview on palatogenesis and introduce key signaling cascades leading to formation of the intact palate. Moreover, we review conceptual differences between pathway biology and network biology and discuss how some of the recent technological advances in conjunction with mouse genetic models have contributed to our understanding of signaling networks regulating palate growth and fusion. For further resources related to this article, please visit the WIREs website.

CONFLICT OF INTEREST

The authors have declared no conflicts of interest for this article.

Strong association of variants around FOXE1 and orofacial clefting

Nonsyndromic orofacial clefting (nsOFC) is a common, complex congenital disorder. The most frequent forms are nonsyndromic cleft lip with or without cleft palate (nsCL/P) and nonsyndromic cleft palate only (nsCPO). Although they are generally considered distinct entities, a recent study has implicated a region around the FOXE1 gene in both nsCL/P and nsCPO. To investigate this hypothesis, we analyzed the 2 most strongly associated markers (rs3758249 and rs4460498) in 2 independent samples of differing ethnicities: Central European (949 nsCL/P cases, 155 nsCPO cases, 1163 controls) and Mayan Mesoamerican (156 nsCL/P cases, 10 nsCPO cases, 338 controls). While highly significant associations for both single-nucleotide polymorphisms were obtained in nsCL/P (rs4460498: p Europe = 6.50 × 10(-06), p Mayan = .0151; rs3758249: p Europe = 2.41 × 10(-05), p Mayan = .0299), no association was found in nsCPO (p > .05). Genotyping of rs4460498 in 472 independent European trios revealed significant associations for nsCL/P (p = .016) and nsCPO (p = .043). A meta-analysis of all data revealed a genomewide significant result for nsCL/P (p = 1.31 × 10(-08)), which became more significant when nsCPO cases were added (p nsOFC = 1.56 × 10(-09)). These results strongly support the FOXE1 locus as a risk factor for nsOFC. With the data of the initial study, there is now considerable evidence that this locus is the first conclusive risk factor shared between nsCL/P and nsCPO.

Search for genetic modifiers of IRF6 and genotype-phenotype correlations in Van der Woude and popliteal pterygium syndromes

Van der Woude syndrome is the most common form of syndromic orofacial clefting, accounting for 1-2% of all patients with cleft lip and/or cleft palate. Van der Woude and popliteal pterygium syndromes are caused by mutations in IRF6, but phenotypic variability within and among families with either syndrome suggests that other genetic factors contribute to the phenotypes. The aim of this study was to identify common variants acting as genetic modifiers of IRF6 as well as genotype-phenotype correlations based on mutation type and location. We identified an association between mutations in the DNA-binding domain of IRF6 and limb defects (including pterygia). Although we did not detect formally significant associations with the genes tested, borderline associations suggest several genes that could modify the VWS phenotype, including FOXE1, TGFB3, and TFAP2A. Some of these genes are hypothesized to be part of the IRF6 gene regulatory network and may suggest additional genes for future study when larger sample sizes are also available. We also show that families with the Van de Woude phenotype but in whom no mutations have been identified have a lower frequency of cleft lip, suggesting there may be locus and/or mutation class differences in Van de Woude syndrome.

Patterns of FOXE1 expression in papillary thyroid carcinoma by immunohistochemistry

BACKGROUND

FOXE1, a thyroid-specific transcription factor also known as TTF-2, was recently identified as a major genetic risk factor for papillary thyroid carcinoma (PTC). Its role in thyroid carcinogenesis, however, remains unknown. The purpose of the present study was to assess the relationship between the FOXE1 immunohistochemical features and the clinical and genetic characteristics of PTC.

METHODS

Immunohistochemical staining of FOXE1 was performed in 48 PTC cases. Two single nucleotide polymorphisms immediately inside (rs1867277) or in the vicinity (rs965513) of the FOXE1 gene were genotyped by direct sequencing. Histopathological, clinical, and genetic data were included in statistical analyses.

RESULTS

FOXE1 exhibited cytoplasmic overexpression in tumor tissue compared to the normal counterpart (p<0.001). Both cancer and normal thyroid cells demonstrated the highest FOXE1 scores in the areas closest to the tumor border (<300 μm) compared with more distant areas (p<0.001). No differences in FOXE1 staining distributions were found between microcarcinomas and PTC of larger size, between different histopathological variants of PTC, and encapsulated and nonencapsulated tumors. Multivariate regression analysis revealed that nuclear FOXE1 expression in neoplastic cells in the vicinity of the tumor border independently associated with the genotype at rs1867277 (the dominant model of inheritance, p=0.037) and tumor multifocality (p=0.032), and with marginal significance with capsular invasion (p=0.051).

CONCLUSIONS

FOXE1 overexpression and translocation to the cytoplasm are phenotypic hallmarks of tumor cells suggesting that FOXE1 is involved in the pathogenesis of PTC. Nuclear FOXE1 expression in tumor cells in the vicinity of the PTC border is associated with the presence of a risk allele of rs1867277 (c.-238G>A) in the 5' untranslated region of the FOXE1 gene, as well as with pathological characteristics of PTC, suggesting possible FOXE1 involvement in the facilitation of tumor development beginning at an early stage.

New insights into FoxE1 functions: identification of direct FoxE1 targets in thyroid cells

BACKGROUND

FoxE1 is a thyroid-specific forkhead transcription factor essential for thyroid gland development, as well as for the maintenance of the thyroid differentiated state in adults. FoxE1 recognizes and binds to a short DNA sequence present in thyroglobulin (Tg) and thyroperoxidase (Tpo) promoters, but FoxE1 binding to regulatory regions other than Tg and Tpo promoters remains almost unexplored. Improving knowledge of the regulatory functions of FoxE1 is necessary to clarify its role in endocrine syndromes and cancer susceptibility.

METHODOLOGY/PRINCIPAL FINDING

In order to further investigate downstream FoxE1 targets, we performed a genome-wide expression screening after knocking-down FoxE1 and obtained new insights into FoxE1 transcriptional networks in thyroid follicular cells. After validation, we confirmed Adamts9, Cdh1, Duox2 and S100a4 as upregulated genes and Casp4, Creld2, Dusp5, Etv5, Hsp5a, Nr4a2 and Tm4sf1 as downregulated genes when FoxE1 was silenced. In promoter regions of putative FoxE1-regulated genes and also in the promoters of the classical thyroid genes Nis, Pax8 and Titf1, we performed an in silico search of the FoxE1 binding motif that was in close proximity to the NF1/CTF binding sequence, as previously described for other forkhead factors. Using chromatin immunoprecipitation we detected specific in vivo FoxE1 binding to novel regulatory regions in two relevant thyroid genes, Nis and Duox2. Moreover, we demonstrated simultaneous binding of FoxE1 and NF1/CTF to the Nis upstream enhancer region, as well as a clear functional activation of the Nis promoter by both transcription factors.

CONCLUSIONS/SIGNIFICANCE

In search for potential downstream mediators of FoxE1 function in thyroid cells, we identified two novel direct FoxE1 target genes. To our knowledge, this is the first evidence regarding the implication of Nis and Duox2 in executing the transcriptional program triggered by FoxE1. Furthermore, this study points out the important role of FoxE1 in the regulation of a large number of genes in thyroid cells.

Molecular signaling along the anterior-posterior axis of early palate development

Cleft palate is a common congenital birth defect in humans. In mammals, the palatal tissue can be distinguished into anterior bony hard palate and posterior muscular soft palate that have specialized functions in occlusion, speech or swallowing. Regulation of palate development appears to be the result of distinct signaling and genetic networks in the anterior and posterior regions of the palate. Development and maintenance of expression of these region-specific genes is crucial for normal palate development. Numerous transcription factors and signaling pathways are now recognized as either anterior- (e.g., Msx1, Bmp4, Bmp2, Shh, Spry2, Fgf10, Fgf7, and Shox2) or posterior-specific (e.g., Meox2, Tbx22, and Barx1). Localized expression and function clearly highlight the importance of regional patterning and differentiation within the palate at the molecular level. Here, we review how these molecular pathways and networks regulate the anterior-posterior patterning and development of secondary palate. We hypothesize that the anterior palate acts as a signaling center in setting up development of the secondary palate.

Mechanisms of thyroid development and dysgenesis: an analysis based on developmental stages and concurrent embryonic anatomy

Thyroid dysgenesis is the most common cause of congenital hypothyroidism that affects 1 in 3000 newborns. Although a number of pathogenetic mutations in thyroid developmental genes have been identified, the molecular mechanism of disease is unknown in most cases. This chapter summarizes the current knowledge of normal thyroid development and puts the different developmental stages in perspective, from the time of foregut endoderm patterning to the final shaping of pharyngeal anatomy, for understanding how specific malformations may arise. At the cellular level, we will also discuss fate determination of follicular and C-cell progenitors and their subsequent embryonic growth, migration, and differentiation as the different thyroid primordia evolve and merge to establish the final size and shape of the gland.

A Bmp reporter with ultrasensitive characteristics reveals that high Bmp signaling is not required for cortical hem fate

Insights into Bone morphogenetic protein (Bmp) functions during forebrain development have been limited by a lack of Bmp signaling readouts. Here we used a novel Bmp signaling reporter ("BRE-gal" mice) to study Bmp signaling in the dorsal telencephalon. At early stages, BRE-gal expression was restricted to the dorsal telencephalic midline. At later stages, strong BRE-gal expression occurred in neurons of the marginal zone and dentate gyrus. Comparisons to nuclear phospho-Smad1/5/8 (pSmad) and Msx1 indicated that BRE-gal expression occurred exclusively in neural cells with high-level Bmp signaling. BRE-gal responsiveness to Bmps was confirmed in reporter-negative cortical cells cultured with Bmp4, and both in vivo and in vitro, BRE-gal expression was switch-like, or ultrasensitive. In the early dorsal telencephalon, BRE-gal expression negatively correlated with the cortical selector gene Lhx2, indicating a BRE-gal expression border that coincides with the cortex-hem boundary. However, in Lhx2 null chimeras, neither BRE-gal nor nuclear pSmad increases were observed in ectopic hem cells. These findings establish BRE-gal as an ultrasensitive reporter of Bmp signaling in the dorsal telencephalon, imply that hem fate can be specified at different Bmp signaling intensities, and suggest that Lhx2 primarily regulates the responses to--rather than the intensity of--Bmp signaling in dorsal telencephalic cells.

Downregulation of Foxe1 by HR suppresses Msx1 expression in the hair follicles of Hr(Hp) mice

Hairless (HR), a transcriptional cofactor, is highly expressed in the skin and brain. To characterize the effects of HR expression in the skin, we examined its capacity for transcriptional regulation of its target genes in mouse skin and keratinocytes. We found that Foxe1 mRNA expression was suppressed in HR-overexpressing skin, as well as in HR-expressing keratinocytes. In turn, Msx1 expression was downregulated contingent on Foxe1 downregulation in skin and keratinocytes. We also found that expression of Sfrp1 was also correlated with that of Foxe1. Further investigation of the mechanisms involved in the transcriptional regulation of these genes will facilitate our understanding of the relationship among genes involved in hair follicle morphogenesis and cycling.

Accelerated evolution of 3'avian FOXE1 genes, and thyroid and feather specific expression of chicken FoxE1

Background

The forkhead transcription factor gene E1 (FOXE1) plays an important role in regulation of thyroid development, palate formation and hair morphogenesis in mammals. However, avian FOXE1 genes have not been characterized and as such, codon evolution of FOXE1 orthologs in a broader evolutionary context of mammals and birds is not known.

Results

In this study we identified the avian FOXE1 gene in chicken, turkey and zebra finch, all of which consist of a single exon. Chicken and zebra finch FOXE1 are uniquely located on the sex-determining Z chromosome. In situ hybridization shows that chicken FOXE1 is specifically expressed in the developing thyroid. Its expression is initiated at the placode stage and is maintained during the stages of vesicle formation and follicle primordia. Based on this expression pattern, we propose that avian FOXE1 may be involved in regulating the evagination and morphogenesis of thyroid. Chicken FOXE1 is also expressed in growing feathers. Sequence analysis identified two microdeletions in the avian FOXE1 genes, corresponding to the loss of a transferable repression domain and an engrailed homology motif 1 (Eh1) C-terminal to the forkhead domain. The avian FOXE1 proteins exhibit a significant sequence divergence of the C-terminus compared to those of amphibian and mammalian FOXE1. The codon evolution analysis (dN/dS) of FOXE1 shows a significantly increased dN/dS ratio in the avian lineages, consistent with either a relaxed purifying selection or positive selection on a few residues in avian FOXE1 evolution. Further site specific analysis indicates that while relaxed purifying selection is likely to be a predominant cause of accelerated evolution at the 3'-region of avian FOXE1, a few residues might have evolved under positive selection.

Conclusions

We have identified three avian FOXE1 genes based on synteny and sequence similarity as well as characterized the expression pattern of the chicken FOXE1 gene during development. Our evolutionary analyses suggest that while a relaxed purifying selection is likely to be the dominant force driving accelerated evolution of avian FOXE1 genes, a few residues may have evolved adaptively. This study provides a basis for future genetic and comparative biochemical studies of FOXE1.