Metabolic characterization of all-trans-retinoic acid (ATRA)-induced craniofacial development of murine embryos using in vivo proton magnetic resonance spectroscopy

LncRNA-NONMMUT100923.1 regulates mouse embryonic palatal shelf adhesion by sponging miR-200a-3p to modulate medial epithelial cell desmosome junction during palatogenesis

Cleft palate (CP) is a common neonatal craniofacial defect caused by the adhesion and fusion dysfunction of bilateral embryonic palatal shelf structures. Long non-coding RNA (lncRNA) is involved in CP formation with regulatory mechanism unknown. In this study, all-trans retinoic acid (ATRA) was used to induced cleft palate in embryonic mice as model group. The RNA-sequencing was performed to screen differentially expressed genes between the normal and model group on embryonic day 16.5, and the expression of LncRNA-NONMMUT100923.1 and miR-200a-3p, Cdsn was confirmed by RT-PCR and western blotting. Colony formation, CCK-8 and EDU assays were performed to measure cell proliferation and apoptosis on mouse embryonic palatal shelf (MEPS) epithelial cells in vitro. Fluorescence in situ hybridization (FISH) and dual luciferase activity assays was used to investigate the regulatory effect of LncRNA-NONMMUT100923.1 on miRNA and its target genes. Up-regulation of LncRNA-NONMMUT100923.1 and Cdsn while downregulation of miR-200a-3p was found in the model group. The sponging effects of LncRNA-NONMMUT100923 on miR-200a-3p and the target gene relations between Cdsn and miR-200a-3p was confirmed. Low expression of miR-200a-3p was related to the increased expressed levels of Cdsn and the proliferation of MEPS epithelial cells. Thus, a potential ceRNA regulatory network in which LncRNA-NONMMUT100923.1 regulates Cdsn expression by competitively binding to endogenous miR-200a-3p during palatogenesis, which may inhibit MEPS adhesion by preventing the disintegration of the desmosome junction in medial edge epithelium cells. These findings indicate the regulatory role of lncRNA and provides a potential direction for target gene therapy of CP.

Computational model for fetal skeletal defects potentially linked to disruption of retinoic acid signaling

All-trans retinoic acid (ATRA) gradients determine skeletal patterning morphogenesis and can be disrupted by diverse genetic or environmental factors during pregnancy, leading to fetal skeleton defects. Adverse Outcome Pathway (AOP) frameworks for ATRA metabolism, signaling, and homeostasis allow for the development of new approach methods (NAMs) for predictive toxicology with less reliance on animal testing. Here, a data-driven model was constructed to identify chemicals associated with both ATRA pathway bioactivity and prenatal skeletal defects. The phenotype data was culled from ToxRefDB prenatal developmental toxicity studies and produced a list of 363 ToxRefDB chemicals with altered skeletal observations. Defects were classified regionally as cranial, post-cranial axial, appendicular, and other (unspecified) features based on ToxRefDB descriptors. To build a multivariate statistical model, high-throughput screening bioactivity data from >8,070 chemicals in ToxCast/Tox21 across 10 in vitro assays relevant to the retinoid signaling system were evaluated and compared to literature-based candidate reference chemicals in the dataset. There were 48 chemicals identified for effects on both in vivo skeletal defects and in vitro ATRA pathway targets for computational modeling. The list included 28 chemicals with prior evidence of skeletal defects linked to retinoid toxicity and 20 chemicals without prior evidence. The combination of thoracic cage defects and DR5 (direct repeats of 5 nucleotides for RAR/RXR transactivation) disruption was the most frequently occurring phenotypic and target disturbance, respectively. This data model provides valuable AOP elucidation and validates current mechanistic understanding. These findings also shed light on potential avenues for new mechanistic discoveries related to ATRA pathway disruption and associated skeletal dysmorphogenesis due to environmental exposures.

A LCMS-based untargeted lipidomics analysis of cleft palate in mouse

BACKGROUND

Recent studies have shown that lipid metabolism was abnormal during the formation of cleft palate. However, the composition of these lipid species remains unclear.

OBJECTIVE

Aims of this study were to identify the lipid species components and reveal the key lipid metabolic disorders in cleft palate formation.

METHODS

The pregnant mice were divided into experimental group exposed to all-trans retinoic acid (RA-treated group) (n = 12) and control group (n = 12) at embryonic gestation day 10.5 (E0.5). The component of the palatal tissue metabolome was analyzed using a LCMS-based nontargeted lipidomics approach. Multivariate statistical analysis was then carried out to assess the differences between the RA-treated group and the control group.

RESULTS

Twenty-nine lipid species were found to discriminate between RA-treated and control embryos. Among them, 28 lipid species increased and 1 lipid species decreased in the RA-treated group. Among these lipids, 13 were triglycerides, 9 were PEs, 3 were PCs, 2 were PSs, 2 were DGs. Further analysis of the number of carbons and unsaturated bond of triglycerides showed that TGs with high unsaturated bonds constituted a higher fraction in the RA-treated group. A higher amount of triglycerides containing 52, 54, 56, 58, 60 carbons, and 1 to 8 unsaturated bonds. Of note, under RA treatment, TG 50:1, 52:2, 56:6and 60:8 became the most prominent.

CONCLUSION

Lipid metabolism is significantly different in the formation of cleft palate induced by RA, and the unsaturated triglycerides increased in the RA-treated group may play an important role in the formation of cleft palate.

MiR-106a-5p modulates apoptosis and metabonomics changes by TGF-β/Smad signaling pathway in cleft palate

BACKGROUND

The molecular mechanisms of abnormal palatogenesis were investigated in this study. A key regulator, miR-106a-5p, and its target pathway were analyzed.

OBJECTIVES

This research is trying to clarify the underlying mechanism of the modulation of miRNA transcription during the formation of cleft palate by 7T and 9.4T NMR metabolomic platforms.

METHOD

Differentially expressed miRNAs and mRNAs were analyzed by microarray analysis and verified by qRT-PCR. The protein expression in TGFβ signaling pathways were analyzed by Western Blotting. The relationship between miR-106a-5p and TGFβ were analyzed by luciferase reporter assay. Cell apoptosis were analyzed by flow cytometer. And finally, the metabonomics were analyzed by NMR and multivariate data analysis models (MVDA).

RESULTS

The expression of miR-106a-5p increased in cleft palatal tissue and negatively correlated with the protein level of Tgfbr2. The luciferase assay further proved that the tgfbr2 was a direct target of miR-106a-5p. In another aspect, miR-106a-5p increased apoptosis level in palatal mesenchymal cells, possibly because its inhibition of TGFβ signaling pathway. Moreover, low cholesterol and choline levels with high citric acid and lipid levels were observed by 7T and 9.4T NMR metabonomic analysis, which inferred the disorder of cell membrane synthesis in cleft palate formation. Furthermore, transformation from choline to phosphatidylcholine regulated by miR-106a-5p was also disrupted, resulting in phosphatidic choline synthesis disorder and reduced cell membrane synthesis.

CONCLUSIONS

The regulatory mechanism of cleft palate was studied at transcriptional and metabolomics levels, which may provide important information in understanding the primary cause of this abnormality.

Multi‑organ assessment via a 9.4‑Tesla MRS evaluation of metabolites during the embryonic development of cleft palate induced by dexamethasone

The aim of the present study was to determine the association between maternal metabolism and development of the fetal palate, and to suggest a potential non‑invasive prenatal diagnostic method for fetal cleft palate (CP). Dexamethasone (DXM) was used to create a CP mouse model. A 9.4‑Tesla (T) magnetic resonance spectroscopy (MRS) imager was used to measure an array of metabolites in the maternal serum, placental tissue, amniotic fluid and fetal palates. Multivariate statistical analysis was performed using SIMCA‑P 14.1 software. Following DXM treatment, variations were detected in multiple metabolites in the female mice and their fetuses based on 9.4T MRS. It was indicated that in the experimental group during CP formation, leucine, valine, creatine, acetate and citrate levels in the palatal tissue were lower, whereas lactate, alanine, proline/inositol and glutamate‑containing metabolite levels were higher, compared with the levels in the control group. In placental tissue and amniotic fluid, succinate and choline levels were lower in the experimental group. The relative concentrations of cholesterol and lipids in palatal tissues from mice treated with DXM were higher compared with the concentrations in tissues from mice in the control group, with the exception of (CH2)n lipids. In the placental tissue, the alteration in cholesterol level exhibited the opposite trend. Lipid levels for the different lipid forms varied and most of them were unsaturated lipids.

Identification of circular RNA-associated competing endogenous RNA network in the development of cleft palate

Circular RNAs (circRNAs) serve as competing endogenous RNAs (ceRNAs) and indirectly regulate gene expression through shared microRNAs (miRNAs). However, the regulatory mechanisms of circRNA as ceRNA associated with the fusion of palatal shelves in palatogenesis are yet unclear. This study aimed to explore the potential mechanism underlying the role of circRNA as ceRNA in cleft palate (CP). First, we systematically analyzed RNA-seq and miRNA-seq data after high-throughput sequencing for embryonic palatal shelf tissues from a mouse CP model induced by maternal exposure to all-trans retinoic acid on embryonic gestation day 14.5 (E14.5). Thirty-nine circRNAs, 18 miRNAs, and 936 messenger RNAs (mRNAs) were significantly dysregulated (log₂ [fold change {FC}] > 1; P < 0.05). Thereafter, we constructed a circRNA-associated ceRNA network. Finally, we determined the circRNA_0954-miRNA-881-3p-PRKAR1α ceRNA network as a hub involved in palatogenesis. Gene Ontology analysis revealed that ceRNA-related genes were associated with facial morphogenesis and developmental gene silencing. Kyoto Encyclopedia of Genes and Genomes pathway analysis indicated that ceRNA-related genes are involved in apoptosis (P < 0.05, fold enrichment >1). Quantitative reverse transcription polymerase chain reaction was performed to verify the results of ceRNA analysis. We found that the circRNA-miRNA-mRNA ceRNA network is involved in palatogenesis. The present results imply that circRNA_0954-miRNA-881-3p-PRKAR1α ceRNA network may cause dysfunctional palatal fusion and might facilitate the development of novel epigenetic biomarkers to treat CP in the future.

AMBRA1-mediated autophagy and apoptosis associated with an epithelial-mesenchymal transition in the development of cleft palate induced by all-trans retinoic acid

BACKGROUND

Autophagy and apoptosis are involved in embryogenesis. However, little is known about the regulatory mechanism of AMBRA1-mediated autophagy and apoptosis associated with epithelial-mesenchymal transition (EMT) in the development of cleft palate (CP). This study is aimed to elucidate a novel regulatory mechanism by which AMBRA1 regulates autophagy and apoptosis associated with EMT during palatal fusion.

METHODS

We performed lncRNA and mRNA co-expression profile analysis on embryonic gestation day 14.5 (E14.5) mouse embryos from control (n=3) and all-trans retinoic acid-treated (to induce cleft palate, n=3) C57BL/6J mice. Functional prediction for transcription factor (TF)-target gene relationship, which was obtained using Gene Ontology/Kyoto Encyclopedia of Genes and Genomes analyses (GO/KEGG) pathway analysis, identified the regulatory "lncRNA-TF-target gene" using the trans model.

RESULTS

The trans analysis revealed that some TFs (e.g., LEF1, SMAD4, and FOXD3) regulate lncRNA and gene expression. Finally, we identified a NONMMUT034790.2-LEF1-AMBRA1 trans-regulatory network associated with CP. Our results indicate that AMBRA1 might be a novel epigenetic biomarker in palatogenesis.

CONCLUSIONS

AMBRA1-mediated autophagy and apoptosis associated with EMT by a NONMMUT034790.2-LEF1-AMBRA1 trans-regulatory network might be an important mechanism underlying dysfunctional palatal fusion.

Genome-Wide mRNA-Seq Profiling Reveals that LEF1 and SMAD3 Regulate Epithelial-Mesenchymal Transition Through the Hippo Signaling Pathway During Palatal Fusion

BACKGROUND

Epithelial-mesenchymal transition (EMT) of the medial edge epithelium (MEE) occurs through fusion of the palatal shelves and is a crucial step in palatogenesis. The key genes, however, and the related signaling pathway of EMT are not yet fully understood. Therefore, the aim of this study was to reveal the key genes and the related signaling pathway of EMT during palatal fusion.

MATERIALS AND METHODS

C57BL/6J mice at embryonic gestation day 14.5 (E14.5; n = 6) were used to establish the cleft palate model for mRNA-Seq (HiSeq X Ten). The Gene Ontology (GO) and the Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses were performed for functional annotations of the differentially expressed genes. Quantitative polymerase chain reaction (qPCR) assays were used to validate the RNAseq data.

RESULTS

A total of 936 differentially expressed genes, including 558 upregulated and 378 downregulated genes were identified in cases versus controls, respectively. Among these genes, the GO analysis showed that Lymphoid Enhancer-Binding Factor 1 (LEF1) and SMAD Family Member 3 (SMAD3) significantly enriched biological processes, which were EMT related. The KEGG analysis showed that these genes regulated EMT through the Hippo signaling pathway. LEF1 and SMAD3 were downregulated, and the qPCR results corroborated the RNA-seq data.

CONCLUSIONS

These results demonstrate that LEF1 and SMAD3 inhibits EMT at the MEE through the Hippo signaling pathway; and that this could contribute to cleft palate formation in embryonic palatal fusion at E 14.5.

A novel lncRNA-mediated trans-regulatory mechanism in the development of cleft palate in mouse

BACKGROUND

Increasing evidence indicates that long non-coding RNAs (lncRNAs) play crucial regulatory roles in epithelial-mesenchymal transition (EMT). However, the regulatory mechanisms during EMT of the medial edge epithelium (MEE) remain elusive. The aim of this work is to reveal a novel lncRNA-regulated dysfunction of EMT involved in the development of cleft palate (CP).

METHODS

C57BL/6 J mice at embryonic gestation day 14.5 (n = 6, 3 case samples vs. 3 control samples) were used to establish the CP model for lncRNA-mRNA co-expression profile analysis after high-throughput sequencing. Functional predictions for the differentially expressed lncRNA-mRNA co-expression with transcription factor (TF)-target gene relationship Gene Ontology/Kyoto Encyclopedia of Genes and Genomes pathway (GO/KEGG) analyses identified the regulatory "lncRNA-TF-target gene" trans model.

RESULTS

A total of 583 differentially expressed lncRNAs and 703 differentially expressed mRNAs were identified. The results of trans analysis revealed that some TFs (LEF1, SMAD4, and FOXD3) regulate lncRNAs and gene expression. Finally, we identified the NONMMUT034790.2-LEF1-SMAD7 co-expression trans-regulatory network that might be associated with CP.

CONCLUSIONS

Our results revealed that NONMMUT034790.2 might be a novel epigenetic biomarker in CP. The integration of lncRNA modulators into trans-regulatory networks will further enhance our understanding of lncRNA functions and regulatory mechanisms during palatal fusion in ATRA-induced mouse CP.

DNA hypermethylation of Fgf16 and Tbx22 associated with cleft palate during palatal fusion

OBJECTIVE

Cleft palate (CP) is a congenital birth defect caused by the failure of palatal fusion. Little is known about the potential role of DNA methylation in the pathogenesis of CP. This study aimed to explore the potential role of DNA methylation in the mechanism of CP.

METHODOLOGY

We established an all-trans retinoic acid (ATRA)-induced CP model in C57BL/6J mice and used methylation-dependent restriction enzymes (MethylRAD, FspEI) combined with high-throughput sequencing (HiSeq X Ten) to compare genome-wide DNA methylation profiles of embryonic mouse palatal tissues, between embryos from ATRA-treated vs. untreated mice, at embryonic gestation day 14.5 (E14.5) (n=3 per group). To confirm differentially methylated levels of susceptible genes, real-time quantitative PCR (qPCR) was used to correlate expression of differentially methylated genes related to CP.

RESULTS

We identified 196 differentially methylated genes, including 17,298 differentially methylated CCGG sites between ATRA-treated vs. untreated embryonic mouse palatal tissues (P<0.05, log2FC>1). The CP-related genes Fgf16 (P=0.008, log2FC=1.13) and Tbx22 (P=0.011, log2FC=1.64,) were hypermethylated. Analysis of Fgf16 and Tbx22, using Gene Ontology (GO) and the Kyoto Encyclopedia of Genes and Genomes (KEGG), identified 3 GO terms and 1 KEGG pathway functionally related to palatal fusion. The qPCR showed that changes in expression level negatively correlated with methylation levels.

CONCLUSIONS

Taken together, these results suggest that hypermethylation of Fgf16 and Tbx22 is associated with decreased gene expression, which might be responsible for developmental failure of palatal fusion, eventually resulting in the formation of CP.

Temporal transcriptomic analysis of metabolic genes in maternal organs and placenta during murine pregnancy

Maternal pregnancy adaptation is crucial for fetal development and long-term health. Complex interactions occur between maternal digestive and excretory systems as they interface with the developing fetus through the placenta, and transcriptomic regulation in these organs throughout pregnancy is poorly understood. Our objective is to characterize transcriptomic changes across gestation in maternal organs and placenta. Gene expression was quantified in the kidney, liver, and small intestine harvested from nonpregnant and pregnant FVB mice at four time points and placenta at three time points (N = 5/time point) using Affymetrix Mouse Gene 1.0 ST arrays. In maternal organs, we identified 476 genes in the liver, 207 genes in the kidney, and 27 genes in the small intestine that were differentially expressed across gestation (False Discovery Rate [FDR] adjusted q < 0.05). The placenta had a total of 1576 differentially expressed genes between the placenta at either/gd15 or gd19 compared to gd10. We identified a number of pathways enriched for genes differentially expressed across gestation, including 5 pathways in the placenta, 9 pathways in the kidney, and 28 pathways in the liver, including the citrate cycle, retinol metabolism, bile acid synthesis, and steroid bile synthesis, which play functional roles in fetal development and pregnancy maintenance. Characterization of normal longitudinal changes that occur in pregnancy provides context to understand how perturbations in these biochemical pathways and perturbations in nutrient signaling may impact pregnancy.

Correlation between HDAC4 enhancer DNA methylation and mRNA expression during palatal fusion induced by all-trans retinoic acid

Epithelial-mesenchymal transformation of the medial edge epithelium is the most crucial process in embryonic palatal fusion. This study aimed to explore the relationship and potential mechanism between enhancer DNA methylation and mRNA expression of histone deacetylase 4 (HDAC4) during palatal fusion induced by maternal exposure to all-trans retinoic acid (ATRA). Pregnant mice were administered ATRA (70 mg/kg) by gavage at embryonic gestation day 10.5 (E10.5) to establish a cleft palate (CP) model in C57BL/6J mice. Control groups were given an equivalent volume of corn oil. Pregnant mice were dissected at E14.5 (n = 6) to obtain embryonic palates. HDAC4 enhancer DNA methylation data were obtained from a previous MethylRAD-seq. Methylation-specific polymerase chain reaction (MSP) and real-time quantitative PCR were used to quantify enhancer methylation and the mRNA expression level of HDAC4. Enhancer DNA methylation at a non-CpG site within the HDAC4 gene was hyper-methylated at E14.5 (P: 0.011, log₂ FC:1.67). The MSP results indicated a similar trend, in agreement with the MethylRAD-seq results. The change in the HDAC4 expression level was negatively correlated with its enhancer DNA methylation level, at the non-CpG site, during palatal fusion induced by ATRA. Enhancer DNA methylation of HDAC4 might play an important regulatory role during palatogenesis, especially in embryonic palatal fusion at E 14.5, and may facilitate the development of novel epigenetic biomarkers in the treatment of CP.

Prenatal retinoic acid exposure reveals candidate genes for craniofacial disorders

Syndromes that display craniofacial anomalies comprise a major class of birth defects. Both genetic and environmental factors, including prenatal retinoic acid (RA) exposure, have been associated with these syndromes. While next generation sequencing has allowed the discovery of new genes implicated in these syndromes, some are still poorly characterized such as Oculo-Auriculo-Vertebral Spectrum (OAVS). Due to the lack of clear diagnosis for patients, developing new strategies to identify novel genes involved in these syndromes is warranted. Thus, our study aimed to explore the link between genetic and environmental factors. Owing to a similar phenotype of OAVS reported after gestational RA exposures in humans and animals, we explored RA targets in a craniofacial developmental context to reveal new candidate genes for these related disorders. Using a proteomics approach, we detected 553 dysregulated proteins in the head region of mouse embryos following their exposure to prenatal RA treatment. This novel proteomic approach implicates changes in proteins that are critical for cell survival/apoptosis and cellular metabolism which could ultimately lead to the observed phenotype. We also identified potential molecular links between three major environmental factors known to contribute to craniofacial defects including maternal diabetes, prenatal hypoxia and RA exposure. Understanding these links could help reveal common key pathogenic mechanisms leading to craniofacial disorders. Using both in vitro and in vivo approaches, this work identified two new RA targets, Gnai3 and Eftud2, proteins known to be involved in craniofacial disorders, highlighting the power of this proteomic approach to uncover new genes whose dysregulation leads to craniofacial defects.

Genome-Wide DNA Methylation Profile of Gene cis-Acting Element Methylations in All-trans Retinoic Acid-Induced Mouse Cleft Palate

DNA methylation epigenetically regulates gene expression. This study is aimed to investigate genome-wide DNA methylations involved in the regulation of palatal fusion in the all-trans retinoic acid-induced mouse cleft palate model. There were 4,718,556 differentially CCGG methylated sites and 367,504 CCWGG methylated sites for 1497 genes between case and control embryonic mouse palatal tissues. The enhancers (HDAC4 and SMAD3) and promoter (MID1) of these three genes had cis-acting element methylation. HDAC4 is localized within the CCWGG, while MID1 and SMAD3 are localized within the CCGG of the gene intron. The methylation-specific polymerase chain reaction data confirmed the MethylRAD-seq results, while the quantitative reverse transcriptase-polymerase chain reaction result showed that changes in gene expression inversely were associated with the cis-acting element methylation of the gene during retinoic acid-induced palatal fusion. The GO and KEGG data showed that these three genes could regulate cell proliferation, skeletal muscle fiber development, and development-related gene signaling or activity. The cis-acting element methylation of HDAC4, SMAD3, and MID1 may play a regulatory role during palatal fusion. Further research is needed to verify these novel epigenetic biomarkers for cleft palate.

Genome-Wide DNA Methylation Analysis During Palatal Fusion Reveals the Potential Mechanism of Enhancer Methylation Regulating Epithelial Mesenchyme Transformation

Epithelial mesenchyme transformation (EMT) of the medial edge epithelium (MEE) is the crucial process during palatal fusion. This work is aimed to elucidate the enhancer regulatory mechanism by genome-wide DNA methylation analysis of EMT during palatal fusion. Over 800 million clean reads, 325 million enzyme reads, and 234 million mapping reads were generated. The mapping rate was 68.85-74.32%, which included differentially methylated 17299 CCGG sites and 2363 CCWGG sites (p < 0.05, log₂FC >1). Methylated sites in intron and intergenic regions were more compared to other regions of all DNA elements. GO and KEGG analysis indicated that differential methylation sites related to embryonic palatal fusion genes (HDAC4, TCF7L2, and PDGFRB) at the enhancer were located on CCWGG region of non-CpG islands. In addition, the results showed that the enhancer for HDAC4 was hypermethylated, whereas the enhancers for TCF7L2 and PDGFRB were hypomethylated. The methylation status of enhancer regions of HDAC4, PDGFRB, and TCF7L2, involved in the regulation of the EMT during palatal fusion, may enlighten the development of novel epigenetic biomarkers in the treatment of cleft palate.

All-Trans Retinoic Acid-Induced Craniofacial Malformation Model: A Prenatal and Postnatal Morphological Analysis

Objective

To characterize the prenatal and postnatal craniofacial bone development in mouse model of all-trans retinoic acid (ATRA) exposure at different ages by a quantitative and morphological analysis of skull morphology.

Methods

Pregnant mice were exposed to ATRA at embryonic day 10 (E10) and 13 (E13) by oral gavage. Skulls of mice embryos at E19.5 and adult mice at postnatal day 35 (P35) were collected for high-resolution microcomputed tomography (microCT) imaging scanning and section HE staining. Reconstruction and measurement of mouse skulls were performed for prenatal and postnatal analysis of the control and ATRA-exposed mice.

Results

Craniofacial malformations in mouse models caused by ATRA exposure were age dependent. ATRA exposure at E10 induced cleft palate in 81.8% of the fetuses, whereas the palatine bone of E13-exposed mice was intact. Inhibitions of maxilla and mandible development with craniofacial asymmetry induced were observed at E19.5 and P35. Compared with control and E13-exposed mice, the palatine bones of E10-exposed mice were not elevated and were smaller in dimension. Some E10-exposed mice exhibited other craniofacial abnormalities, including premature fusion of mandibular symphysis with a missing mandibular incisor and a smaller mandible. Severe deviated snouts and amorphous craniofacial suture were detected in E13-exposed mice at P35.

Conclusion

These morphological variations in E10- and E13-exposed mice suggested that ATRA was teratogenic in craniofacial bone development in mice and the effect was age dependent.