Molecular Anatomy of Palate Development

Spatiotemporal Evolution of Developing Palate in Mice

The intricate formation of the palate involves a series of complex events, yet its mechanistic basis remains uncertain. To explore major cell populations in the palate and their roles during development, we constructed a spatiotemporal transcription landscape of palatal cells. Palate samples from C57BL/6 J mice at embryonic days 12.5 (E12.5), 14.5 (E14.5), and 16.5 (E16.5) underwent single-cell RNA sequencing (scRNA-seq) to identify distinct cell subsets. In addition, spatial enhanced resolution omics-sequencing (stereo-seq) was used to characterize the spatial distribution of these subsets. Integrating scRNA-seq and stereo-seq with CellTrek annotated mesenchymal and epithelial cellular components of the palate during development. Furthermore, cellular communication networks between these cell subpopulations were analyzed to discover intercellular signaling during palate development. From the analysis of the middle palate, both mesenchymal and epithelial populations were spatially segregated into 3 domains. The middle palate mesenchymal subpopulations were associated with tooth formation, ossification, and tissue remodeling, with initial state cell populations located proximal to the dental lamina. The nasal epithelium of the palatal shelf exhibited richer humoral immune responses than the oral side. Specific enrichment of Tgfβ3 and Pthlh signals in the midline epithelial seam at E14.5 suggested a role in epithelial-mesenchymal transition. In summary, this study provides high-resolution transcriptomic information, contributing to a deeper mechanistic understanding of palate biology and pathophysiology.

Progress of epigenetic modification of SATB2 gene in the pathogenesis of non-syndromic cleft lip and palate

Non-syndromic Cleft Lip and Palate (NSCLP) is one of the most common congenital craniofacial malformations. However, there is no enough knowledge about its mechanism, even through many relevant studies verify that cleft lip and palate is caused by interactions between environmental and genetic factors. SATB2 gene is one of the most common candidate genes of NSCLP, and the development of epigenetics provides a new direction on pathogenesis of cleft lip and palate. This review summarizes SATB2 gene in the pathogenesis of non-syndromic cleft lip and palate, expecting to provide strategies to prevent and treat cleft and palate in the future.

Single-cell Transcriptomics Reveals Activation of Macrophages in All-trans Retinoic Acid (atRA)-induced Cleft Palate

Cleft palate is among the most common birth defects with an impact on swallowing and speaking and is difficult to diagnose with ultrasound during pregnancy. In this study, we systematically capture the cellular composition of all-trans retinoic acid (atRA)-exposed and normal embryonic gestation 16.5 days mouse palate by the single-cell RNA sequencing technique. The authors identified 14 major cell types with the largest proportion of fibroblasts. The proportion of myeloid cells in atRA-exposed palate was markedly higher than those in the normal palate tissue, especially M1-like macrophages and monocytes. The upregulated genes of the different expression genes between atRA-exposed palate and normal palate tissue were linked to the biological processes of leukocyte chemotaxis and migration. Protein TLR2, CXCR4, THBS1, MRC1, transcription factor encoding genes Cebpb, Fos, Jun, Rela, and signaling pathway IL-17 and phagosome were found to be significantly involved in these processes. Subsequently, cellular communication network analysis suggested that myeloid-centered cell interactions SELL, SELPLG, MIF, CXCL, ANNEXIN, THBS, and NECTIN were significantly more activated in atRA-exposed palate. Overall, we delineate the single-cell landscape of atRA-induced cleft palate, revealing the effects of overexposure to atRA during palate tissue development and providing insights for the diagnosis of cleft palate.

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

Craniofacial disorders arise in early pregnancy and are one of the most common congenital defects. To fully understand how craniofacial disorders arise, it is essential to characterize gene expression during the patterning of the craniofacial region. To address this, we performed bulk and single-cell RNA-seq on human craniofacial tissue from 4-8 weeks post conception. Comparisons to dozens of other human tissues revealed 239 genes most strongly expressed during craniofacial development. Craniofacial-biased developmental enhancers were enriched +/- 400 kb surrounding these craniofacial-biased genes. Gene co-expression analysis revealed that regulatory hubs are enriched for known disease causing genes and are resistant to mutation in the normal healthy population. Combining transcriptomic and epigenomic data we identified 539 genes likely to contribute to craniofacial disorders. While most have not been previously implicated in craniofacial disorders, we demonstrate this set of genes has increased levels of de novo mutations in orofacial clefting patients warranting further study.

Single cell sequencing of the mouse anterior palate reveals mesenchymal heterogeneity

BACKGROUND

Cleft palate is one of the most prevalent birth defects. Mice are useful for studying palate development because of their morphological and genetic similarities to humans. In mice, palate development occurs between embryonic days (E)11.5 to 15.5. Single cell transcriptional profiles of palate cell populations have been a valuable resource for the craniofacial research community, but we lack a single cell transcriptional profile for anterior palate at E13.5, at the transition from proliferation to shelf elevation.

RESULTS

A detailed single cell RNA sequencing analysis reveals heterogeneity in expression profiles of the cell populations of the E13.5 anterior palate. Hybridization chain reaction RNA fluorescent in situ hybridization (HCR RNA FISH) reveals epithelial populations segregate into layers. Mesenchymal populations spatially segregate into four domains. One of these mesenchymal populations expresses ligands and receptors distinct from the rest of the mesenchyme, suggesting that these cells have a unique function. RNA velocity analysis shows two terminal cell states that contribute to either the proximal or distal palatal regions emerge from a single progenitor pool.

CONCLUSION

This single cell resolution expression data and detailed analysis from E13.5 anterior palate provides a powerful resource for mechanistic insight into secondary palate morphogenesis for the craniofacial research community.

Spatiotemporal Expression and Functional Analysis of miRNA-22 in the Developing Secondary Palate

OBJECTIVE

Normal development of the embryonic orofacial region requires precise spatiotemporal coordination between numerous genes. MicroRNAs represent small, single-stranded, non-coding molecules that regulate gene expression. This study examines the role of microRNA-22 (miR-22) in murine orofacial ontogeny.

METHODS

Spatiotemporal and differential expression of miR-22 (mmu-miR-22-3p) within the developing secondary palate was determined by in situ hybridization and quantitative real-time PCR, respectively. Bioinformatic approaches were used to predict potential mRNA targets of miR-22 and analyze their association with cellular functions indispensable for normal orofacial ontogeny. An in vitro palate organ culture system was used to assess the role of miR-22 in secondary palate development.

RESULTS

There was a progressive increase in miR-22 expression from GD12.5 to GD14.5 in palatal processes. On GD12.5 and GD13.5, miR-22 was expressed in the future oral, nasal, and medial edge epithelia. On GD14.5, miR-22 expression was observed in the residual midline epithelial seam (MES), the nasal epithelium and the mesenchyme, but not in the oral epithelium. Inhibition of miR-22 activity in palate organ cultures resulted in failure of MES removal. Bioinformatic analyses revealed potential mRNA targets of miR-22 that may play significant roles in regulating apoptosis, migration, and/or convergence/extrusion, developmental processes that modulate MES removal during palatogenesis.

CONCLUSIONS

Results from the current study suggest a key role for miR-22 in the removal of the MES during palatogenesis and that miR-22 may represent a potential contributor to the etiology of cleft palate.

Essential role of Msx1 in regulating anterior-posterior patterning of the secondary palate in mice

Development of the secondary palate displays molecular heterogeneity along the anterior-posterior axis; however, the underlying molecular mechanism remains largely unknown. MSX1 is an anteriorly expressed transcription repressor required for palate development. Here, we investigate the role of Msx1 in regional patterning of the secondary palate. The Wnt1-Cre-mediated expression of Msx1 (RosaMsx1^Wnt1-Cre) throughout the palatal mesenchyme leads to cleft palate in mice, associated with aberrant cell proliferation and cell death. Osteogenic patterning of the hard palate in RosaMsx1^Wnt1-Cre mice is severely impaired, as revealed by a marked reduction in palatine bone formation and decreased expression of the osteogenic regulator Sp7. Overexpression and knockout of Msx1 in mice show that the transcription repressor promotes the expression of the anterior palate-specific Alx1 but represses the expression of the medial-posterior palate genes Barx1, Meox2, and Tbx22. Furthermore, Tbx22 constitutes a direct Msx1 target gene in the secondary palate, suggesting that Msx1 can directly repress the expression of medial-posterior specific genes. Finally, we determine that Sp7 is downstream of Tbx22 in palatal mesenchymal cells, suggesting that a Msx1/Tbx22/Sp7 axis participates in the regulation of palate development. Our findings unveil a novel role for Msx1 in regulating the anterior-posterior growth and patterning of the secondary palate.

The role of lncRNA Meg3 in the proliferation of all-trans retinoic acid-treated mouse embryonic palate mesenchymal cells involves the Smad pathway

Mesenchymal cell proliferation is critical for the growth of the palate shelf. All-trans retinoic acid (atRA), as well as pathways associated with TGF-β/Smad signaling, play crucial roles in the proliferation of mouse embryonic palate mesenchymal (MEPM) cells. We have found that MEPM-cell proliferation was regulated by atRA and exogenous TGF-β3 could significantly antagonize the atRA-mediated suppression of MEPM cell proliferation, which is closely associated with the regulation of TGF-β/Smad signaling pathway. The long non-coding RNA (lncRNA) MEG3 has been reported to activate TGF-β/Smad signaling, thereby regulating cellular proliferation, differentiation, and related processes. Here, we found that Meg3 expression increased significantly in atRA-treated MEPM cells while TGF-β3 treatment markedly inhibited Meg3 expression and antagonized the effect of atRA on Meg3. Moreover, Smad2 was found to interact directly with Meg3, and atRA treatment significantly enriched Meg3 in Smad2-immunoprecipitated samples. After Meg3 deletion, the effects of atRA on the proliferation of MEPM cells and TGF-β3-dependent protein expression were lost. Hence, we speculate that Meg3 has a role in the RA-induced suppression of MEPM cell proliferation by targeting Smad2 and thereby mediating TGF-β/Smad signaling inhibition.

Runx2-Twist1 interaction coordinates cranial neural crest guidance of soft palate myogenesis

Cranial neural crest (CNC) cells give rise to bone, cartilage, tendons, and ligaments of the vertebrate craniofacial musculoskeletal complex, as well as regulate mesoderm-derived craniofacial muscle development through cell-cell interactions. Using the mouse soft palate as a model, we performed an unbiased single-cell RNA-seq analysis to investigate the heterogeneity and lineage commitment of CNC derivatives during craniofacial muscle development. We show that Runx2, a known osteogenic regulator, is expressed in the CNC-derived perimysial and progenitor populations. Loss of Runx2 in CNC-derivatives results in reduced expression of perimysial markers (Aldh1a2 and Hic1) as well as soft palate muscle defects in Osr2-Cre;Runx2fl/fl mice. We further reveal that Runx2 maintains perimysial marker expression through suppressing Twist1, and that myogenesis is restored in Osr2-Cre;Runx2fl/fl;Twist1fl/+ mice. Collectively, our findings highlight the roles of Runx2, Twist1, and their interaction in regulating the fate of CNC-derived cells as they guide craniofacial muscle development through cell-cell interactions.

LncRNA Meg3-mediated regulation of the Smad pathway in atRA-induced cleft palate

Palatal mesenchymal cell proliferation is essential to the process of palatogenesis, and the proliferation of mouse embryonic palate mesenchymal (MEPM) cells is impacted by both all-trans retinoic acid (atRA) and the TGF-β/Smad signaling pathway. The long non-coding RNA (lncRNA) MEG3 has been shown to activate TGF-β/Smad signaling and to thereby regulate cell proliferation, differentiation, and related processes. Herein, we found that atRA treatment (100 mg/kg) promoted Meg3 upregulation in MEPM cells, and that such upregulation was linked to the suppression of MEPM cell proliferation in the context of secondary palate fusion on gestational day (GD) 13 and 14. Moreover, the demethylation of specific CpG sites within the lncRNA Meg3 promoter was detected in atRA-treated MEPM cells, likely explaining the observed upregulation of this lncRNA. Smad signaling was also suppressed by atRA treatment in these cells, and RNA immunoprecipitation analyses revealed that Smad2 can directly interact with Meg3 in MEPM cells following atRA treatment. Therefore, we propose a model wherein Meg3 is involved in the suppression of MEPM cell proliferation, functioning at least in part via interacting with the Smad2 protein and thereby suppressing Smad signaling in the context of atRA-induced cleft palate.

A developmental stage-specific network approach for studying dynamic co-regulation of transcription factors and microRNAs during craniofacial development

Craniofacial development is regulated through dynamic and complex mechanisms that involve various signaling cascades and gene regulations. Disruption of such regulations can result in craniofacial birth defects. Here, we propose the first developmental stage-specific network approach by integrating two crucial regulators, transcription factors (TFs) and microRNAs (miRNAs), to study their co-regulation during craniofacial development. Specifically, we used TFs, miRNAs and non-TF genes to form feed-forward loops (FFLs) using genomic data covering mouse embryonic days E10.5 to E14.5. We identified key novel regulators (TFs Foxm1, Hif1a, Zbtb16, Myog, Myod1 and Tcf7, and miRNAs miR-340-5p and miR-129-5p) and target genes (Col1a1, Sgms2 and Slc8a3) expression of which changed in a developmental stage-dependent manner. We found that the Wnt-FoxO-Hippo pathway (from E10.5 to E11.5), tissue remodeling (from E12.5 to E13.5) and miR-129-5p-mediated Col1a1 regulation (from E10.5 to E14.5) might play crucial roles in craniofacial development. Enrichment analyses further suggested their functions. Our experiments validated the regulatory roles of miR-340-5p and Foxm1 in the Wnt-FoxO-Hippo subnetwork, as well as the role of miR-129-5p in the miR-129-5p-Col1a1 subnetwork. Thus, our study helps understand the comprehensive regulatory mechanisms for craniofacial development.

Transcriptional analysis of cleft palate in TGFβ3 mutant mice

Cleft palate (CP) is one of the most common craniofacial birth defects, impacting about 1 in 800 births in the USA. Tgf-β3 plays a critical role in regulating murine palate development, and Tgf-β3 null mutants develop cleft palate with 100% penetrance. In this study, we compared global palatal transcriptomes of wild type (WT) and Tgf-β3 −/− homozygous (HM) mouse embryos at the crucial palatogenesis stages of E14.5, and E16.5, using RNA-seq data. We found 1,809 and 2,127 differentially expressed genes at E16.5 vs. E14.5 in the WT and HM groups, respectively (adjusted p < 0.05; |fold change|> 2.0). We focused on the genes that were uniquely up/downregulated in WT or HM at E16.5 vs. E14.5 to identify genes associated with CP. Systems biology analysis relating to cell behaviors and function of WT and HM specific genes identified functional non-Smad pathways and preference of apoptosis to epithelial-mesenchymal transition. We identified 24 HM specific and 11 WT specific genes that are CP-related and/or involved in Tgf-β3 signaling. We validated the expression of 29 of the 35 genes using qRT-PCR and the trend of mRNA expression is similar to that of RNA-seq data . Our results enrich our understanding of genes associated with CP that are directly or indirectly regulated via TGF-β.

CXCL12-CXCR4 Interplay Facilitates Palatal Osteogenesis in Mice

Cranial neural crest cells (CNCCs), identified by expression of transcription factor Sox9, migrate to the first branchial arch and undergo proliferation and differentiation to form the cartilage and bone structures of the orofacial region, including the palatal bone. Sox9 promotes osteogenic differentiation and stimulates CXCL12-CXCR4 chemokine-receptor signaling, which elevates alkaline phosphatase (ALP)-activity in osteoblasts to initiate bone mineralization. Disintegration of the midline epithelial seam (MES) is crucial for palatal fusion. Since we earlier demonstrated chemokine-receptor mediated signaling by the MES, we hypothesized that chemokine CXCL12 is expressed by the disintegrating MES to promote the formation of an osteogenic center by CXCR4-positive osteoblasts. Disturbed migration of CNCCs by excess oxidative and inflammatory stress is associated with increased risk of cleft lip and palate (CLP). The cytoprotective heme oxygenase (HO) enzymes are powerful guardians harnessing injurious oxidative and inflammatory stressors and enhances osteogenic ALP-activity. By contrast, abrogation of HO-1 or HO-2 expression promotes pregnancy pathologies. We postulate that Sox9, CXCR4, and HO-1 are expressed in the ALP-activity positive osteogenic regions within the CNCCs-derived palatal mesenchyme. To investigate these hypotheses, we studied expression of Sox9, CXCL12, CXCR4, and HO-1 in relation to palatal osteogenesis between E15 and E16 using (immuno)histochemical staining of coronal palatal sections in wild-type (wt) mice. In addition, the effects of abrogated HO-2 expression in HO-2 KO mice and inhibited HO-1 and HO-2 activity by administrating HO-enzyme activity inhibitor SnMP at E11 in wt mice were investigated at E15 or E16 following palatal fusion. Overexpression of Sox9, CXCL12, CXCR4, and HO-1 was detected in the ALP-activity positive osteogenic regions within the palatal mesenchyme. Overexpression of Sox9 and CXCL12 by the disintegrating MES was detected. Neither palatal fusion nor MES disintegration seemed affected by either HO-2 abrogation or inhibition of HO-activity. Sox9 progenitors seem important to maintain the CXCR4-positive osteoblast pool to drive osteogenesis. Sox9 expression may facilitate MES disintegration and palatal fusion by promoting epithelial-to-mesenchymal transformation (EMT). CXCL12 expression by the MES and the palatal mesenchyme may promote osteogenic differentiation to create osteogenic centers. This study provides novel evidence that CXCL12-CXCR4 interplay facilitates palatal osteogenesis and palatal fusion in mice.

Transforming Growth Factor-Beta and Sonic Hedgehog Signaling in Palatal Epithelium Regulate Tenascin-C Expression in Palatal Mesenchyme During Soft Palate Development

During palatogenesis, the palatal shelves first grow vertically on either side of the tongue before changing their direction of growth to horizontal. The extracellular matrix (ECM) plays an important role in these dynamic changes in palatal shelf morphology. Tenascin-C (TNC) is an ECM glycoprotein that shows unique expression in the posterior part of the palatal shelf, but little is known about the regulation of TNC expression. Since transforming growth factor-beta-3 (TGF-β3) and sonic hedgehog (SHH) signaling are known to play important roles in palatogenesis, we investigated whether TGF-β3 and SHH are involved in the regulation of TNC expression in the developing palate. TGF-β3 increased the expression of TNC mRNA and protein in primary mouse embryonic palatal mesenchymal cells (MEPM) obtained from palatal mesenchyme dissected at embryonic day 13.5-14.0. Interestingly, immunohistochemistry experiments revealed that TNC expression was diminished in K14-cre;Tgfbr2 ^fl/fl mice that lack the TGF-β type II receptor in palatal epithelial cells and exhibit cleft soft palate, whereas TNC expression was maintained in Wnt1-cre;Tgfbr2 ^fl/fl mice that lack the TGF-β type II receptor in palatal mesenchymal cells and exhibit a complete cleft palate. SHH also increased the expression of TNC mRNA and protein in MEPM cells. However, although TGF-β3 up-regulated TNC mRNA and protein expression in O9-1 cells (a cranial neural crest cell line), SHH did not. Furthermore, TGF-β inhibited the expression of osteoblastic differentiation markers (osterix and alkaline phosphatase) and induced the expression of fibroblastic markers (fibronectin and periostin) in O9-1 cells, whereas SHH did not affect the expression of osteoblastic and fibroblastic markers in O9-1 cells. However, immunohistochemistry experiments showed that TNC expression was diminished in the posterior palatal shelves of Shh^-/+ ;MFCS4 ^+/- mice, which have deficient SHH signaling in the posterior palatal epithelium. Taken together, our findings support the proposal that TGF-β and SHH signaling in palatal epithelium co-ordinate the expression of TNC in the posterior palatal mesenchyme through a paracrine mechanism. This signal cascade may work in the later stage of palatogenesis when cranial neural crest cells have differentiated into fibroblast-like cells. The spatiotemporal regulation of ECM-related proteins by TGF-β and SHH signaling may contribute not only to tissue construction but also to cell differentiation or determination along the anterior-posterior axis of the palatal shelves.

Integrating gene regulatory pathways into differential network analysis of gene expression data

The advent of next-generation sequencing has introduced new opportunities in analyzing gene expression data. Research in systems biology has taken advantage of these opportunities by gleaning insights into gene regulatory networks through the analysis of gene association networks. Contrasting networks from different populations can reveal the many different roles genes fill, which can lead to new discoveries in gene function. Pathologies can also arise from aberrations in these gene-gene interactions. Exposing these network irregularities provides a new avenue for understanding and treating diseases. A general framework for integrating known gene regulatory pathways into a differential network analysis between two populations is proposed. The framework importantly allows for any gene-gene association measure to be used, and inference is carried out through permutation testing. A simulation study investigates the performance in identifying differentially connected genes when incorporating known pathways, even if the pathway knowledge is partially inaccurate. Another simulation study compares the general framework with four state-of-the-art methods. Two RNA-seq datasets are analyzed to illustrate the use of this framework in practice. In both examples, the analysis reveals genes and pathways that are known to be biologically significant along with potentially novel findings that may be used to motivate future research.

Spatio-Temporal Expression and Functional Analysis of miR-206 in Developing Orofacial Tissue

BACKGROUND

Development of the mammalian palate is dependent on precise, spatiotemporal expression of a panoply of genes. MicroRNAs (miRNAs), the largest family of noncoding RNAs, function as crucial modulators of cell and tissue differentiation, regulating expression of key downstream genes.

OBSERVATIONS

Our laboratory has previously identified several developmentally regulated miRNAs, including miR-206, during critical stages of palatal morphogenesis. The current study reports spatiotemporal distribution of miR-206 during development of the murine secondary palate (gestational days 12.5-14.5).

RESULT AND CONCLUSION

Potential cellular functions and downstream gene targets of miR-206 were investigated using functional assays and expression profiling, respectively. Functional analyses highlighted potential roles of miR-206 in governing TGFß- and Wnt signaling in mesenchymal cells of the developing secondary palate. In addition, altered expression of miR-206 within developing palatal tissue of TGFß3-/- fetuses reinforced the premise that crosstalk between this miRNA and TGFß3 is crucial for secondary palate development.

Laser Capture Microdissection of Murine Embryonic Neural Crest Cells

The purpose of this chapter is to provide a step-by-step protocol to enable performance of laser capture microdissection (LCM) on tissue sections from mammalian embryos or postnatal organism stages in order to collect pure populations of neural crest cells from which sufficient amounts of nucleic acids and/or protein can be obtained for quantitative analysis. The methods (1) define a strategy to genetically and indelibly label mammalian neural crest-derived cells with a fluorescent marker, thus enabling their isolation throughout the pre- and postnatal life span of the organism, and (2) describe subsequent isolation by LCM of the labeled neural crest cells, or their derivatives, from embryonic/postnatal tissue cryosections. Details are provided for using the Arcturus PixCell^®IIe Laser Capture Microdissection System (Arcturus) and CapSure LCM Caps (Thermo Fisher Scientific), to which the selected cells adhere upon laser-mediated capture. The protocol outlined herein can be applied in any situation wherein limited cellular samples are available for isolation by LCM. Nucleic acids or proteins can be extracted from LCM-isolated cells and processed for high-density gene expression profiling analyses (microarrays or RNA sequencing), Real-Time PCR (q-PCR) for specific candidate gene expression, investigation of DNA methylation, as well as for varied protein analyses.

Associations Between Genetic Data and Quantitative Assessment of Normal Facial Asymmetry

Human facial asymmetry is due to a complex interaction of genetic and environmental factors. To identify genetic influences on facial asymmetry, we developed a method for automated scoring that summarizes local morphology features and their spatial distribution. A genome-wide association study using asymmetry scores from two local symmetry features was conducted and significant genetic associations were identified for one asymmetry feature, including genes thought to play a role in craniofacial disorders and development: NFATC1, SOX5, NBAS, and TCF7L1. These results provide evidence that normal variation in facial asymmetry may be impacted by common genetic variants and further motivate the development of automated summaries of complex phenotypes.

A Combined PLS and Negative Binomial Regression Model for Inferring Association Networks from Next-Generation Sequencing Count Data

A major challenge of genomics data is to detect interactions displaying functional associations from large-scale observations. In this study, a new cPLS-algorithm combining partial least squares approach with negative binomial regression is suggested to reconstruct a genomic association network for high-dimensional next-generation sequencing count data. The suggested approach is applicable to the raw counts data, without requiring any further pre-processing steps. In the settings investigated, the cPLS-algorithm outperformed the two widely used comparative methods, graphical lasso, and weighted correlation network analysis. In addition, cPLS is able to estimate the full network for thousands of genes without major computational load. Finally, we demonstrate that cPLS is capable of finding biologically meaningful associations by analyzing an example data set from a previously published study to examine the molecular anatomy of the craniofacial development.

Systems biology of facial development: contributions of ectoderm and mesenchyme

The rapid increase in gene-centric biological knowledge coupled with analytic approaches for genomewide data integration provides an opportunity to develop systems-level understanding of facial development. Experimental analyses have demonstrated the importance of signaling between the surface ectoderm and the underlying mesenchyme are coordinating facial patterning. However, current transcriptome data from the developing vertebrate face is dominated by the mesenchymal component, and the contributions of the ectoderm are not easily identified. We have generated transcriptome datasets from critical periods of mouse face formation that enable gene expression to be analyzed with respect to time, prominence, and tissue layer. Notably, by separating the ectoderm and mesenchyme we considerably improved the sensitivity compared to data obtained from whole prominences, with more genes detected over a wider dynamic range. From these data we generated a detailed description of ectoderm-specific developmental programs, including pan-ectodermal programs, prominence- specific programs and their temporal dynamics. The genes and pathways represented in these programs provide mechanistic insights into several aspects of ectodermal development. We also used these data to identify co-expression modules specific to facial development. We then used 14 co-expression modules enriched for genes involved in orofacial clefts to make specific mechanistic predictions about genes involved in tongue specification, in nasal process patterning and in jaw development. Our multidimensional gene expression dataset is a unique resource for systems analysis of the developing face; our co-expression modules are a resource for predicting functions of poorly annotated genes, or for predicting roles for genes that have yet to be studied in the context of facial development; and our analytic approaches provide a paradigm for analysis of other complex developmental programs.

Computational Model of Secondary Palate Fusion and Disruption

Morphogenetic events are driven by cell-generated physical forces and complex cellular dynamics. To improve our capacity to predict developmental effects from chemical-induced cellular alterations, we built a multicellular agent-based model in CompuCell3D that recapitulates the cellular networks and collective cell behavior underlying growth and fusion of the mammalian secondary palate. The model incorporated multiple signaling pathways (TGFβ, BMP, FGF, EGF, and SHH) in a biological framework to recapitulate morphogenetic events from palatal outgrowth through midline fusion. It effectively simulated higher-level phenotypes (e.g., midline contact, medial edge seam (MES) breakdown, mesenchymal confluence, and fusion defects) in response to genetic or environmental perturbations. Perturbation analysis of various control features revealed model functionality with respect to cell signaling systems and feedback loops for growth and fusion, diverse individual cell behaviors and collective cellular behavior leading to physical contact and midline fusion, and quantitative analysis of the TGF/EGF switch that controls MES breakdown-a key event in morphogenetic fusion. The virtual palate model was then executed with theoretical chemical perturbation scenarios to simulate switch behavior leading to a disruption of fusion following chronic (e.g., dioxin) and acute (e.g., retinoic acid) chemical exposures. This computer model adds to similar systems models toward an integrative "virtual embryo" for simulation and quantitative prediction of adverse developmental outcomes following genetic perturbation and/or environmental disruption.