Genetic linkage analysis of head and neck cancer in a Spanish family

OBJECTIVES

To describe the genetic variants that may be associated with the development of head and neck cancer (HNC) and functionally validating the molecular implications.

MATERIALS AND METHODS

A prospective observational study was carried out on a family of 3 generations in which 3 members had developed HNC. Peripheral blood sample was taken in a routine procedure for exome sequencing in one relative and genotyping in the remaining twelve relatives. For the functional analysis all-trans retinoic acid (atRA) was extracted from saliva and serum and measured using ultra-performance liquid chromatography-tandem mass spectrometry (UPLC-MS/MS). The presence of HPV-DNA.

RESULTS

None of the patients smoked or consumed alcohol. The presence of HPV DNA was not detected in any of the biopsied samples. A total amount of 6 members out of 13 (46.15%) carried out the same mutation of CYP26B1 (2p13.2; G>T). The mean plasma concentration of atRA was 3.3109 ± 1.4791 pg/mL for the study family and 4.7370 ± 1.5992 pg/mL for the controls (p = 0.042).

CONCLUSION

Lower levels of atRA were confirmed in the study family, which may open the way to the possible relationship between the polymorphism CYP26B1 (2p13.2; G>T) and HNC.

Elongation of the developing spinal cord is driven by Oct4-type transcription factor-mediated regulation of retinoic acid signaling in zebrafish embryos

BACKGROUND

Elongation of the spinal cord is dependent on neural development from neuromesodermal progenitors in the tail bud. We previously showed the involvement of the Oct4-type gene, pou5f3, in this process in zebrafish mainly by dominant-interference gene induction, but, to compensate for the limitation of this transgene approach, mutant analysis was indispensable. pou5f3 involvement in the signaling pathways was another unsolved question.

RESULTS

We examined the phenotypes of pou5f3 mutants and the effects of Pou5f3 activation by the tamoxifen-ERT2 system in the posterior neural tube, together confirming the involvement of pou5f3. The reporter assays using P19 cells implicated tail bud-related transcription factors in pou5f3 expression. Regulation of tail bud development by retinoic acid (RA) signaling was confirmed by treatment of embryos with RA and the synthesis inhibitor, and in vitro reporter assays further showed that RA signaling regulated pou5f3 expression. Importantly, the expression of the RA degradation enzyme gene, cyp26a1, was down-regulated in embryos with disrupted pou5f3 activity.

CONCLUSIONS

The involvement of pou5f3 in spinal cord extension was supported by using mutants and the gain-of-function approach. Our findings further suggest that pou5f3 regulates the RA level, contributing to neurogenesis in the posterior neural tube.

Expression patterns of CYP26A1, FGF8, CDKN1A, and NPVF in the developing rhesus monkey retina

The fovea centralis (fovea) is a specialized region of the primate retina that plays crucial roles in high-resolution visual acuity and color perception. The fovea is characterized by a high density of cone photoreceptors and no rods, and unique anatomical properties that contribute to its remarkable visual capabilities. Early histological analyses identified some of the key events that contribute to foveal development, but the mechanisms that direct the specification of this area are not understood. Recently, the expression of the retinoic acid-metabolizing enzyme CYP26A1 has become a hallmark of some of the retinal specializations found in vertebrates, including the primate fovea and the high-acuity area in avian species. In chickens, the retinoic acid pathway regulates the expression of FGF8 to then direct the development of a rod-free area. Similarly, high levels of CYP26A1, CDKN1A, and NPVF expression have been observed in the primate macula using transcriptomic approaches. However, which retinal cells express these genes and their expression dynamics in the developing primate eye remain unknown. Here, we systematically characterize the expression patterns of CYP26A1, FGF8, CDKN1A, and NPVF during the development of the rhesus monkey retina, from early stages of development in the first trimester until the third trimester (near term). Our data suggest that some of the markers previously proposed to be fovea-specific are not enriched in the progenitors of the rhesus monkey fovea. In contrast, CYP26A1 is expressed at high levels in the progenitors of the fovea, while it localizes in a subpopulation of macular Müller glia cells later in development. Together these data provide invaluable insights into the expression dynamics of several molecules in the nonhuman primate retina and highlight the developmental advancement of the foveal region.

CLMP is a tumor suppressor that determines all-trans retinoic acid response in colorectal cancer

CAR-like membrane protein (CLMP) is a tight junction-associated protein whose mutation is associated with congenital short bowel syndrome (CSBS), but its functions in colorectal cancer (CRC) remain unknown. Here, we demonstrate that CLMP is rarely mutated but significantly decreased in CRC patients, and its deficiency accelerates CRC tumorigenesis, growth, and resistance to all-trans retinoic acid (ATRA). Mechanistically, CLMP recruits β-catenin to cell membrane, independent of cadherin proteins. CLMP-mediated β-catenin translocation inactivates Wnt(Wingless and INT-1)/β-catenin signaling, thereby suppressing CRC tumorigenesis and growth in ApcMin/+, azoxymethane/dextran sodium sulfate (AOM/DSS), and orthotopic CRC mouse models. As a direct target of Wnt/β-catenin, cytochrome P450 hydroxylase A1 (CYP26A1)-an enzyme that degrades ATRA to a less bioactive retinoid-is upregulated by CLMP deficiency, resulting in ATRA-resistant CRC that can be reversed by administering CYP26A1 inhibitor. Collectively, our data identify the anti-CRC role of CLMP and suggest that CYP26A1 inhibitor enable to boost ATRA's therapeutic efficiency.

CYP26B1-related disorder: expanding the ends of the spectrum through clinical and molecular evidence

CYP26B1 metabolizes retinoic acid in the developing embryo to regulate its levels. A limited number of individuals with pathogenic variants in CYP26B1 have been documented with a varied phenotypic spectrum, spanning from a severe manifestation involving skull anomalies, craniosynostosis, encephalocele, radio-humeral fusion, oligodactyly, and a narrow thorax, to a milder presentation characterized by craniosynostosis, restricted radio-humeral joint mobility, hearing loss, and intellectual disability. Here, we report two families with CYP26B1-related phenotypes and describe the data obtained from functional studies of the variants. Exome and Sanger sequencing were used for variant identification in family 1 and family 2, respectively. Family 1 reflects a mild phenotype, which includes craniofacial dysmorphism with brachycephaly (without craniosynostosis), arachnodactyly, reduced radioulnar joint movement, conductive hearing loss, learning disability-and compound heterozygous CYP26B1 variants: (p.[(Pro118Leu)];[(Arg234Gln)]) were found. In family 2, a stillborn fetus presented a lethal phenotype with spina bifida occulta, hydrocephalus, poor skeletal mineralization, synostosis, limb defects, and a synonymous homozygous variant in CYP26B1: c.1083C > A. A minigene assay revealed that the synonymous variant created a new splice site, removing part of exon 5 (p.Val361_Asp382del). Enzymatic activity was assessed using a luciferase assay, demonstrating a notable reduction in exogenous retinoic acid metabolism for the variant p.Val361_Asp382del. (~ 3.5 × decrease compared to wild-type); comparatively, the variants p.(Pro118Leu) and p.(Arg234Gln) demonstrated a partial loss of metabolism (1.7× and 2.3× reduction, respectively). A proximity-dependent biotin identification assay reaffirmed previously reported ER-resident protein interactions. Additional work into these interactions is critical to determine if CYP26B1 is involved with other biological events on the ER. Immunofluorescence assay suggests that mutant CYP26B1 is still localized in the endoplasmic reticulum. These results indicate that novel pathogenic variants in CYP26B1 result in varying levels of enzymatic activity that impact retinoic acid metabolism and relate to the distinct phenotypes observed.

Hyperglycemia alters retinoic acid catabolism in embryos exposed to a maternal diabetic milieu

Pregestational diabetes is highly associated with increased risk of birth defects. We previously reported that the expression of Cyp26a1, the major catabolizing enzyme for controlling retinoic acid (RA) homeostasis, is significantly down-regulated in embryos of diabetic mice, thereby increasing the embryo's susceptibility to malformations caused by RA dysregulation. However, the underlying mechanism for the down-regulation of Cyp26a1 remains unclear. This study aimed to investigate whether elevated maternal blood glucose in the diabetic milieu is a critical factor for the altered Cyp26a1 expression. Streptozotozin-induced diabetic pregnant mice were treated with phlorizin (PHZ) to reduce blood glucose concentrations via induction of renal glucosuria. Embryonic Cyp26a1 expression level, RA catabolic activity and susceptibility to various RA-induced abnormalities were examined. To test the dose-dependent effect of glucose on Cyp26a1 level, early head-fold stage rat embryos of normal pregnancy were cultured in vitro with varying concentrations of D-glucose, followed by quantification of Cyp26a1 transcripts. We found that Cyp26a1 expression, which was down-regulated in diabetic pregnancy, could be normalized under reduced maternal blood glucose level, concomitant with an increase in RA catabolic activity in embryonic tissues. Such normalization could successfully reduce the susceptibility to different RA-induced malformations including caudal regression, cleft palate and renal malformations. The expression level of Cyp26a1 in the embryo was inversely correlated with D-glucose concentrations. Diabetic patients suffer from retinopathy, dermopathy, male infertility and increased cancer risk. Coincidentally, RA dysregulation is also associated with these health problems. Our results provided evidence that elevated glucose can down-regulate Cyp26a1 expression level and disturb RA homeostasis, shedding light on the possibility of affecting the health of diabetic patients via a similar mechanism.

Cyp26a1 supports postnatal retinoic acid homeostasis and glucoregulatory control

Considerable evidence confirms the importance of Cyp26a1 to all-trans-retinoic acid (RA) homeostasis during embryogenesis. In contrast, despite its presence in postnatal liver as a potential major RA catabolizing enzyme and its acute sensitivity to induction by RA, some data suggested that Cyp26a1 contributes only marginally to endogenous RA homeostasis postnatally. We report reevaluation of a conditional Cyp26a1 knockdown in the postnatal mouse. The current results show that Cyp26a1 mRNA in WT mouse liver increases 16-fold upon refeeding after a fast, accompanied by an increased rate of RA elimination and a 41% decrease in the RA concentration. In contrast, Cyp26a1 mRNA in the refed homozygotic knockdown reached only 2% of its extent in WT during refeeding, accompanied by a slower rate of RA catabolism and no decrease in liver RA, relative to fasting. Refed homozygous knockdown mice also had decreased Akt1 and 2 phosphorylation and pyruvate dehydrogenase kinase 4 (Pdk4) mRNA and increased glucokinase (Gck) mRNA, glycogen phosphorylase (Pygl) phosphorylation, and serum glucose, relative to WT. Fasted homozygous knockdown mice had increased glucagon/insulin relative to WT. These data indicate that Cyp26a1 participates prominently in moderating the postnatal liver concentration of endogenous RA and contributes essentially to glucoregulatory control.

A Functional Polymorphism Downstream of Vitamin A Regulator Gene CYP26B1 Is Associated with Hand Osteoarthritis

While genetic analyses have revealed ~100 risk loci associated with osteoarthritis (OA), only eight have been linked to hand OA. Besides, these studies were performed in predominantly European and Caucasian ancestries. Here, we conducted a genome-wide association study in the Han Chinese population to identify genetic variations associated with the disease. We recruited a total of 1136 individuals (n = 420 hand OA-affected; n = 716 unaffected control subjects) of Han Chinese ancestry. We carried out genotyping using Axiom Asia Precisi on Medicine Research Array, and we employed the RegulomeDB database and RoadMap DNase I Hypersensitivity Sites annotations to further narrow down our potential candidate variants. Genetic variants identified were tested in the Geisinger's hand OA cohort selected from the Geisinger MyCode community health initiative (MyCode^®). We also performed a luciferase reporter assay to confirm the potential impact of top candidate single-nucleotide polymorphisms (SNPs) on hand OA. We identified six associated SNPs (p-value = 6.76 × 10^-7-7.31 × 10^-6) clustered at 2p13.2 downstream of the CYP26B1 gene. The strongest association signal identified was rs883313 (p-value = 6.76 × 10^-7, odds ratio (OR) = 1.76), followed by rs12713768 (p-value = 1.36 × 10^-6, OR = 1.74), near or within the enhancer region closest to the CYP26B1 gene. Our findings showed that the major risk-conferring CC haplotype of SNPs rs12713768 and rs10208040 [strong linkage disequilibrium (LD); D' = 1, r² = 0.651] drives 18.9% of enhancer expression activity. Our findings highlight that the SNP rs12713768 is associated with susceptibility to and severity of hand OA in the Han Chinese population and that the suggested retinoic acid signaling pathway may play an important role in its pathogenesis.

In Vitro Transfection of Up-Regulated Genes Identified in Favorable-Outcome Neuroblastoma into Cell Lines

We previously used microarrays to show that high expression of DHRS3, NROB1, and CYP26A1 predicts favorable NB outcomes. Here, we investigated whether expression of these genes was associated with suppression of NB cell (SK-N-SH, NB12, and TGW) growth. We assessed morphology and performed growth, colony-formation, and migration assays, as well as RNA sequencing. The effects of the transient expression of these genes were also assessed with a tetracycline-controlled expression (Tet-On) system. Gene overexpression reduced cell growth and induced morphological senescence. Gene-expression analysis identified pathways involving cellular senescence and cell adhesion. In these cells, transduced gene dropout occurred during passage, making long-term stable gene transfer difficult. Tet-On-induced gene expression caused more pronounced cell-morphology changes. Specifically, DHRS3 and NROB1 led to rapid inhibition and arrest of cell growth, though CYP26A1 did not affect cell-growth rate or cell cycle. DHRS3 arrested the cell cycle by interacting with the all-trans-retinol pathway and drove differentiation and senescence in tumors. Overexpression of these genes reduced the malignant grade of these cells. A new therapeutic strategy might be the induction of these genes, as they suppress the growth of high-risk neuroblastoma and lead to differentiation and senescence.

Distinctive Metabolism-Associated Gene Clusters That Are Also Prognostic in Intrahepatic Cholangiocarcinoma and Hepatocellular Carcinoma

Objective

To offer new prognostic evaluations by exploring potentially distinctive genetic features of hepatocellular carcinoma (HCC) and intrahepatic cholangiocarcinoma (ICC).

Methods

There were 12 samples for gene expression profiling processes in this study. These included three HCC lesion samples and their matched adjacent nontumor liver tissues obtained from patients with HCC, as well as three ICC samples and their controls collected similarly. In addition to the expression matrix generated on our own, profiles of other cohorts from The Cancer Genome Atlas (TCGA) program and the Gene Expression Omnibus (GEO) were also employed in later bioinformatical analyses. Differential analyses, functional analyses, protein interaction network analyses, and gene set variation analyses were used to identify key genes. To establish the prognostic models, univariate/multivariate Cox analyses and subsequent stepwise regression were applied, with the Akaike information criterion evaluating the goodness of fitness.

Results

The top three pathways enriched in HCC were all metabolism-related; they were fatty acid degradation, retinol metabolism, and arachidonic acid metabolism. In ICC, on the other hand, additional pathways related to fat digestion and absorption and cholesterol metabolism were identified. Consistent characteristics of such a metabolic landscape were observed across different cohorts. A prognostic risk score model for calculating HCC risk was constructed, consisting of ADH4, ADH6, CYP2C9, CYP4F2, and RDH16. This signature predicts the 3-year survival with an AUC area of 0.708 (95%CI = 0.644 to 0.772). For calculating the risk of ICC, a prognostic risk score model was built upon the expression levels of CYP26A1, NAT2, and UGT2B10. This signature predicts the 3-year survival with an AUC area of 0.806 (95% CI = 0.664 to 0.947).

Conclusion

HCC and ICC share commonly abrupted pathways associated with the metabolism of fatty acids, retinol, arachidonic acids, and drugs, indicating similarities in their pathogenesis as primary liver cancers. On the flip side, these two types of cancer possess distinctive promising biomarkers for predicting overall survival or potential targeted therapies.

CYP26A1 Is a Novel Cancer Biomarker of Pancreatic Carcinoma: Evidence from Integration Analysis and In Vitro Experiments

Background

CYP26A1 has been reported in multiple cancers. However, the role of CYP26A1 in pancreatic cancer (PC) has not been explored.

Method

The public data used for this study was obtained from The Cancer Genome Atlas (TCGA), Gene Expression Omnibus (GEO), and Cancer Cell Line Encyclopedia (CCLE) cell lines. CCK8, colony formation, and EdU assay were used to assess the proliferation ability of cancer cells. Transwell and wound healing assays were used to evaluate the invasion and migration ability of cancer cells. qRT-PCR and western blot assays were used to analyze the RNA and protein level of genes. Survival package was used for prognosis analysis. Gene Set Enrichment Analysis (GSEA) was used to identify biological pathway differences between two groups. ssGSEA analysis was used to quantify the immune microenvironment in PC tissue. GDSC and TIDE analyses were used for sensitivity analysis of chemotherapy and immunotherapy.

Results

Our results showed that CYP26A1 was overexpressed in PC tissue and cell lines. Meanwhile, metastatic PC cell lines tend to have a higher CYP26A1 level compared with the primary PC cell lines based on CCLE data. Moreover, CYP26A1 was associated with worse clinical features. Also, we found that CYP26A1 had a satisfactory efficiency in predicting overall survival, disease-specific survival, and progression-free interval of PC patients, independent of other clinical features. In vitro experiments indicated that CYP26A1 could significantly facilitate the proliferation, invasion, and migration ability of PC cells. GSEA showed that the pathways of angiogenesis, E2F target, MYC target, mTORC signaling, G2M checkpoint, and epithelial-mesenchymal transition were activated in high CYP26A1 patients. Immune infiltration analysis showed that CYP26A1 was positively correlated with macrophages, Th1 cells, and Treg cells, but negatively correlated with Th17 cells. TIDE analysis showed that non_responder patients had a higher CYP26A1 level compared with predicted responder patients of immunotherapy. Drug sensitivity analysis and assay showed that CYP26A1 could increase the chemotherapy sensitivity of gemcitabine.

Conclusions

In summary, CYP26A1 promotes PC progression and is a novel biomarker of PC, with potential for clinical application.

Multilayer omics analysis reveals a non-classical retinoic acid signaling axis that regulates hematopoietic stem cell identity

Hematopoietic stem cells (HSCs) rely on complex regulatory networks to preserve stemness. Due to the scarcity of HSCs, technical challenges have limited our insights into the interplay between metabolites, transcription, and the epigenome. In this study, we generated low-input metabolomics, transcriptomics, chromatin accessibility, and chromatin immunoprecipitation data, revealing distinct metabolic hubs that are enriched in HSCs and their downstream multipotent progenitors. Mechanistically, we uncover a non-classical retinoic acid (RA) signaling axis that regulates HSC function. We show that HSCs rely on Cyp26b1, an enzyme conventionally considered to limit RA effects in the cell. In contrast to the traditional view, we demonstrate that Cyp26b1 is indispensable for production of the active metabolite 4-oxo-RA. Further, RA receptor beta (Rarb) is required for complete transmission of 4-oxo-RA-mediated signaling to maintain stem cells. Our findings emphasize that a single metabolite controls stem cell fate by instructing epigenetic and transcriptional attributes.

Polymorphisms in Vitamin A-Related Genes and Their Functions in Autoimmune Thyroid Disease

Background: Vitamin A is a factor that suppresses immune responses, including T helper (Th)1 and Th17 responses. However, there has been no report showing the association between vitamin A-related genes (CYP26B1, RARB, and RARG) and the prognosis of autoimmune thyroid disease (AITD). The objective of this study was to clarify the association between vitamin A-related genes and the susceptibility and prognosis of AITD. Methods: We genotyped polymorphisms in genes encoding vitamin A-related molecules using the polymerase chain reaction-restriction fragment length polymorphism method. The proportion of T helper cells was analyzed by flow cytometry. Serum interleukin (IL)-17 and interferon (IFN)-γ were examined by enzyme-linked immunosorbent assay. Results:CYP26B1 rs3768641 GG genotype and G allele were significantly more frequent in patients with mild Hashimoto's thyroiditis (HT) than in those with severe HT (p = 0.0013 and 0.0024, respectively). The RARB rs1997352 CC genotype was significantly more frequent in HT patients than in controls (p = 0.0207). The proportion of Th17 cells was significantly higher in CYP26B1 rs2241057 TT genotype than C carrier (CC+CT genotypes) (p = 0.0385), in RARB rs1997352 A carrier (AA+AC genotypes) than those with CC genotype (p = 0.0246), and in RARG rs7398676 G carrier (GG+GT genotypes) than in TT genotype (p = 0.0249). In the RARB rs1997352 polymorphism, HT patients with a high concentration of IFN-γ (≥150 ng/mL) were more frequent in the CC genotype than in A carriers (AA+AC genotypes) (p = 0.0226). Serum levels of IL-17 were significantly elevated in subjects with the TT genotype of the CYP26B1 rs2241057 single nucleotide polymorphism (SNP) (p = 0.0026) and in subjects with the GG genotype of the CYP26B1 rs3798641 SNP (p = 0.030). Subjects with a high concentration of IL-17 (≥0.71 pg/mL) were more frequent in RARG 7398676 G carriers (GG+GT genotypes) than in TT genotype (p = 0.0218). Conclusions: Polymorphisms in the CYP26B1 gene were related to the proportion of Th17 cells, the level of IL-17 and the severity of HT. Polymorphisms in RAR were related to the proportion of Th17 cells, concentrations of IFN-γ and IL-17, and susceptibility to HT.

A novel Cytochrome P450 26A1 expressing NK cell subset at the mouse maternal-foetal interface

Cyp26a1 had important roles in mouse embryo implantation and was highly expressed in some of NK cells at the human maternal-foetal interface in early pregnancy. However, the regulatory effect of Cyp26a1 on NK cells remains poorly understood. Through qPCR and flow cytometric assays, we found that Cyp26a1 was expressed by mouse uterine NK cells but not spleen NK cells during the peri-implantation period and there was a group of NK cells that highly expressed Cyp26a1, that is Cyp26a1+ NK cell subset. single cell-population transcriptome sequencing on Cyp26a1+ NK and Cyp26a1- NK cell subsets was performed. We found that there were 3957 differentially expressed genes in the Cyp26a1+ NK cell subset with a cut-off of fold change ≥2 and FDR < 0.01, 2509 genes were up-regulated and 1448 genes were down-regulated in Cyp26a1+ NK cell subset. Moreover, cytokine-cytokine receptor interaction signalling pathway and natural killer cell-mediated cytotoxicity signalling pathway were enriched according to KEGG pathway enrichment analysis. We further found that the expression of Gzma and Klrg1 was significantly increased and Fcgr4 was significantly decreased when inhibiting Cyp26a1. Our experimental results show that there is a novel NK cell subset of Cyp26a1+ NK cells in mouse uterus and Cyp26a1 can regulate the gene expression of Gzma, Klrg1 and Fcgr4 in the Cyp26a1+ NK cells.

Liposomal delivery of hydrophobic RAMBAs provides good bioavailability and significant enhancement of retinoic acid signalling in neuroblastoma tumour cells

Retinoid treatment is employed during residual disease treatment in neuroblastoma, where the aim is to induce neural differentiation or death in tumour cells. However, although therapeutically effective, retinoids have only modest benefits and suffer from poor pharmacokinetic properties. In vivo, retinoids induce CYP26 enzyme production in the liver, enhancing their own rapid metabolic clearance, while retinoid resistance in tumour cells themselves is considered to be due in part to increased CYP26 production. Retinoic acid metabolism blocking agents (RAMBAs), which inhibit CYP26 enzymes, can improve retinoic acid (RA) pharmacokinetics in pre-clinical neuroblastoma models. Here, we demonstrate that in cultured neuroblastoma tumour cells, RAMBAs enhance RA action as seen by morphological differentiation, AKT signalling and suppression of MYCN protein. Although active as retinoid enhancers, these RAMBAs are highly hydrophobic and their effective delivery in humans will be very challenging. Here, we demonstrate that such RAMBAs can be loaded efficiently into cationic liposomal particles, where the RAMBAs achieve good bioavailability and activity in cultured tumour cells. This demonstrates the efficacy of RAMBAs in enhancing retinoid signalling in neuroblastoma cells and shows for the first time that liposomal delivery of hydrophobic RAMBAs is a viable approach, providing novel opportunities for their delivery and application in humans.

Cyp26b1 is an essential regulator of distal airway epithelial differentiation during lung development

Proper organ development depends on coordinated communication between multiple cell types. Retinoic acid (RA) is an autocrine and paracrine signaling molecule essential for the development of most organs, including the lung. Despite extensive work detailing effects of RA deficiency in early lung morphogenesis, little is known about how RA regulates late gestational lung maturation. Here, we investigate the role of the RA catabolizing protein Cyp26b1 in the lung. Cyp26b1 is highly enriched in lung endothelial cells (ECs) throughout development. We find that loss of Cyp26b1 leads to reduction of alveolar type 1 cells, failure of alveolar inflation and early postnatal lethality in mouse. Furthermore, we observe expansion of distal epithelial progenitors, but no appreciable changes in proximal airways, ECs or stromal populations. Exogenous administration of RA during late gestation partially mimics these defects; however, transcriptional analyses comparing Cyp26b1-/- with RA-treated lungs reveal overlapping, but distinct, responses. These data suggest that defects observed in Cyp26b1-/- lungs are caused by both RA-dependent and RA-independent mechanisms. This work reports crucial cellular crosstalk during lung development involving Cyp26b1-expressing endothelium and identifies a novel RA modulator in lung development.

Retinoic acid degradation shapes zonal development of vestibular organs and sensitivity to transient linear accelerations

Each vestibular sensory epithelium in the inner ear is divided morphologically and physiologically into two zones, called the striola and extrastriola in otolith organ maculae, and the central and peripheral zones in semicircular canal cristae. We found that formation of striolar/central zones during embryogenesis requires Cytochrome P450 26b1 (Cyp26b1)-mediated degradation of retinoic acid (RA). In Cyp26b1 conditional knockout mice, formation of striolar/central zones is compromised, such that they resemble extrastriolar/peripheral zones in multiple features. Mutants have deficient vestibular evoked potential (VsEP) responses to jerk stimuli, head tremor and deficits in balance beam tests that are consistent with abnormal vestibular input, but normal vestibulo-ocular reflexes and apparently normal motor performance during swimming. Thus, degradation of RA during embryogenesis is required for formation of highly specialized regions of the vestibular sensory epithelia with specific functions in detecting head motions.

CYP26B1 and its implications in lymphangiogenesis: Literature review and study of rare variants in two families

CYP26B1 is a member of the cytochrome P450 family and is responsible for the break-down of retinoic acid for which appropriate levels are important for normal development of the cardiovascular and lymphatic systems. In a cohort of 235 patients with lymphatic malformations, we performed genetic testing for the CYP26B1 gene. These probands had previously tested negative for known lymphedema genes. We identified two heterozygous missense CY-P26B1 variants in two patients. Our bioinformatic study suggested that alterations caused by these variants have no major effect on the overall stability of CYP26B1 protein structure. Balanced levels of retinoic acid maintained by CYP26B1 are crucial for the lymphatic system. We identified that CYP26B1 could be involved in predisposition for lymphedema. We propose that CYP26B1 be further explored as a new candidate gene for genetic testing of lymphedema patients.

Manipulation of retinoic acid signaling in the nucleus accumbens shell alters rat emotional behavior

Novel targets for depression and anxiety disorders are necessary for the development of more effective pharmacotherapeutics. Our previous study found that the retinoic acid (RA) signaling pathway is the signaling pathway most enhanced in the nucleus accumbens (NAc) shell, a region important for depression, anxiety, and addiction. Genetic manipulations of RA signaling in the NAc affecting addiction-related behavior prompted our study of the role of retinoic acid signaling in depression-related and anxiety-related behavior using in vivo RNA interference. Knockdown of the retinoic acid degradation enzyme cytochrome p450 family 26 subfamily b member 1 (Cyp26b1) in the nucleus accumbens shell increased depression-related behavior while decreasing anxiety-like behavior. Knockdown of the retinoic acid binding protein, cellular RA binding protein 2 (Crabp2), also increased depression-related behavior. Knockdown of another RA binding partner fatty acid binding protein 5 (Fabp5), did not alter these behaviors. These results further support the contention that RA signaling in the NAc shell can affect emotional behavior and that targeting some components of this pathway could be a promising avenue for developing novel treatments for depression and anxiety.

Gene Expression in the Supporting Cells at the Onset of Meiosis in Rabbit Gonads

Subsequent to somatic gonadal sexual differentiation, germ cells enter meiosis or mitotic arrest in the ovary or testis, respectively. Among mice, these processes occur almost synchronically in fetal gonads and depend, among other factors, on the levels of retinoic acid (RA). In contrast to those in mice, rabbit germ cells enter meiosis or mitotic arrest after birth and coexist with proliferating germ cells. Here, we studied the somatic cell context in which germ cells enter meiosis or mitotic arrest in the rabbit. Using confocal immunofluorescence and real-time PCR, we studied the expression profiles of ALDH1A1 and ALDH1A2 and, comprising 2 genes required for RA synthesis, 2 meiosis markers STRA8 and SYCP3 as well as 2 genes involved in meiosis inhibition, CYP26B1 and NANOS2. We found that although both meiosis and mitotic arrest initiate after birth, these 2 processes are regulated in a way similar to the human fetal gonad. Current results reinforce the value of the neonatal rabbit gonad as an alternative experimental model for analyzing the direct effect of environmental factors during critical stages of germ cell establishment.

CYP26A1 gene promoter is a useful tool for reporting RAR-mediated retinoid activity

Of numerous genes regulated by retinoic acid (RA), CYP26A1 is the most inducible gene by RA. In this study, we have used a shortened construct form, E4, of the CYP26A1 gene promoter, in a promoter-less vector with either luciferase or red fluorescent protein (RFP) as the reporter gene and have tested its responses to retinoids in transfected HepG2 and HEK293T cells. The promoter responded linearly to a wide concentration range of RA in cells cotransfected with retinoic acid receptors. It also responded quantitatively to retinol and other retinoids. An isolated clonal line of HEK293T cells permanently transfected with the promoter driving the expression of RFP responded to both RA and retinol, and the responses could be measured by fluorescence microscopy and flow cytometry. The promoter was used to assess the retinoid activity of 3 novel synthetic retinoid analogues, as well as of the intact serum samples of rats. Among the synthetic retinoid analogues tested, EC23 is more potent than RA at lower concentrations and was more stable than RA. The retinoid activities could be measured in control rat serum samples and were increased in the serum of RA-treated rats. This system offers a biologically-based alternative to mass-based retinoid analysis.

Ethanol promotes differentiation of embryonic stem cells through retinoic acid receptor-γ

Ethanol (EtOH) is a teratogen, but its teratogenic mechanisms are not fully understood. The alcohol form of vitamin A (retinol/ROL) can be oxidized to all-trans-retinoic acid (RA), which plays a critical role in stem cell differentiation and development. Using an embryonic stem cell (ESC) model to analyze EtOH's effects on differentiation, we show here that EtOH and acetaldehyde, but not acetate, increase differentiation-associated mRNA levels, and that EtOH decreases pluripotency-related mRNAs. Using reporter assays, ChIP assays, and retinoic acid receptor-γ (RARγ) knockout ESC lines generated by CRISPR/Cas9 and homologous recombination, we demonstrate that EtOH signals via RARγ binding to RA response elements (RAREs) in differentiation-associated gene promoters or enhancers. We also report that EtOH-mediated increases in homeobox A1 (Hoxa1) and cytochrome P450 family 26 subfamily A member 1 (Cyp26a1) transcripts, direct RA target genes, require the expression of the RA-synthesizing enzyme, aldehyde dehydrogenase 1 family member A2 (Aldh1a2), suggesting that EtOH-mediated induction of Hoxa1 and Cyp26a1 requires ROL from the serum. As shown with CRISPR/Cas9 knockout lines, the retinol dehydrogenase gene Rdh10 and a functional RARE in the ROL transporter stimulated by retinoic acid 6 (Stra6) gene are required for EtOH induction of Hoxa1 and Cyp26a1 We conclude that EtOH stimulates stem cell differentiation by increasing the influx and metabolism of ROL for downstream RARγ-dependent transcription. In stem cells, EtOH may shift cell fate decisions to alter developmental outcomes by increasing endogenous ROL/RA signaling via increased Stra6 expression and ROL oxidation.

Altered Expression of Retinol Metabolism-Related Genes in an ANIT-Induced Cholestasis Rat Model

Cholestasis is defined as a reduction of bile secretion caused by a dysfunction of bile formation. Insufficient bile secretion into the intestine undermines the formation of micelles, which may result in the reduced absorption of lipids and fat-soluble vitamins. Here, we investigated the retinol homeostasis and the alterations of retinol metabolism-related genes, including β-carotene 15,15′ monooxygenase (BCMO), lecithin:retinol acyltransferase (LRAT), aldehyde dehydrogenase (ALDH), cytochrome P450 26A1 (CYP26A1), and retinoic acid receptors (RAR) β, in a α-naphthyl isothiocyanate (ANIT)-induced cholestasis rat model. Moreover, we examined the expression of the farnesoid X receptor (FXR) target genes. Our results showed that plasma retinol levels were decreased in ANIT rats compared to control rats. On the contrary, hepatic retinol levels were not different between the two groups. The expression of FXR target genes in the liver and intestine of cholestasis model rats was repressed. The BCMO expression was decreased in the liver and increased in the intestine of ANIT rats compared to control rats. Finally, the hepatic expression of LRAT, RARβ, and ALDH1A1 in cholestatic rats was decreased compared to the control rats, while the CYP26A1 expression of the liver was not altered. The increased expression of intestinal BCMO in cholestasis model rats might compensate for decreased circulatory retinol levels. The BCMO expression might be regulated in a tissue-specific manner to maintain the homeostasis of retinol.

RARγ is required for mesodermal gene expression prior to gastrulation in Xenopus

The developing vertebrate embryo is exquisitely sensitive to retinoic acid (RA) concentration, particularly during anteroposterior patterning. In contrast to Nodal and Wnt signaling, RA was not previously considered to be an instructive signal in mesoderm formation during gastrulation. Here, we show in Xenopus that RARγ is indispensable for the expression of early mesoderm markers and is, therefore, an obligatory factor in mesodermal competence and/or maintenance. We identified several novel targets upregulated by RA receptor signaling in the early gastrula that are expressed in the circumblastoporal ring and linked to mesodermal development. Despite overlapping expression patterns of the genes encoding the RA-synthesizing enzyme Aldh1a2 and the RA-degrading enzyme Cyp26a1, RARγ1 functions as a transcriptional activator in early mesoderm development, suggesting that RA ligand is available to the embryo earlier than previously appreciated. RARγ1 is required for cellular adhesion, as revealed by spontaneous dissociation and depletion of ncam1 mRNA in animal caps harvested from RARγ1 knockdown embryos. RARγ1 knockdown obliterates somite boundaries, and causes loss of Myod protein in the presomitic mesoderm, but ectopic, persistent expression of Myod protein in the trunk. Thus, RARγ1 is required for stabilizing the mesodermal fate, myogenic commitment, somite boundary formation, and terminal skeletal muscle differentiation.

RARγ is required for mesodermal gene expression prior to gastrulation in Xenopus

The developing vertebrate embryo is exquisitely sensitive to retinoic acid (RA) concentration, particularly during anteroposterior patterning. In contrast to Nodal and Wnt signaling, RA was not previously considered to be an instructive signal in mesoderm formation during gastrulation. Here, we show in Xenopus that RARγ is indispensable for the expression of early mesoderm markers and is, therefore, an obligatory factor in mesodermal competence and/or maintenance. We identified several novel targets upregulated by RA receptor signaling in the early gastrula that are expressed in the circumblastoporal ring and linked to mesodermal development. Despite overlapping expression patterns of the genes encoding the RA-synthesizing enzyme Aldh1a2 and the RA-degrading enzyme Cyp26a1, RARγ1 functions as a transcriptional activator in early mesoderm development, suggesting that RA ligand is available to the embryo earlier than previously appreciated. RARγ1 is required for cellular adhesion, as revealed by spontaneous dissociation and depletion of ncam1 mRNA in animal caps harvested from RARγ1 knockdown embryos. RARγ1 knockdown obliterates somite boundaries, and causes loss of Myod protein in the presomitic mesoderm, but ectopic, persistent expression of Myod protein in the trunk. Thus, RARγ1 is required for stabilizing the mesodermal fate, myogenic commitment, somite boundary formation, and terminal skeletal muscle differentiation.

Uncovering the genetic lesions underlying the most severe form of Hirschsprung disease by whole-genome sequencing

Hirschsprung disease (HSCR) is a complex birth defect characterized by the lack of ganglion cells along a variable length of the distal intestine. A large proportion of HSCR patients remain genetically unexplained. We applied whole-genome sequencing (WGS) on 9 trios where the probands are sporadically affected with the most severe form of the disorder and harbor no coding sequence variants affecting the function of known HSCR genes. We found de novo protein-altering variants in three intolerant to change genes-CCT2, VASH1, and CYP26A1-for which a plausible link with the enteric nervous system (ENS) exists. De novo single-nucleotide and indel variants were present in introns and non-coding neighboring regions of ENS-related genes, including NRG1 and ERBB4. Joint analysis with those inherited rare variants found under recessive and/or digenic models revealed both patient-unique and shared genetic features where rare variants were found to be enriched in the extracellular matrix-receptor (ECM-receptor) pathway (p = 3.4 × 10-11). Delineation of the genetic profile of each patient might help finding common grounds that could lead to the discovery of shared molecules that could be used as drug targets for the currently ongoing cell therapy effort which aims at providing an alternative to the surgical treatment.

Genetic contribution of retinoid-related genes to neural tube defects

Rare variants are considered underlying causes of complex diseases. The complex and severe group of disorders called neural tube defects (NTDs) results from failure of the neural tube to close during early embryogenesis. Neural tube closure requires the coordination of numerous signaling pathways, including the precise regulation of retinoic acid (RA) concentration, which is controlled by enzymes involved in RA synthesis and degradation. Here, we used a case-control mutation screen study to reveal rare variants in retinoid-related genes in a Han Chinese NTD population by sequencing six genes in 355 NTD cases and 225 controls. More specific rare variants were found in exonic and upstream regions in NTD cases. The RA-responsive genes CYP26A1, CRABP1, and ALDH1A2 harbored NTD-specific rare variants in their upstream regions. Unexpectedly, the majority of missense variants in NTD cases were found in CYP26B1, which encodes a RA degradation enzyme, whereas no missense variants in this gene were found in controls. Functional analysis indicated that the CYP26B1 NTD variants were inefficient in the degradation of RA using assays of RA-induced transcription and RA-initiated neuronal differentiation. Our study supports the contribution of rare variants in RA-related genes to the etiology of human NTDs.

The cycling gonad: retinoic-acid synthesis and degradation patterns during adult zebrafish Danio rerio oogenesis

The expression pattern of genes coding for enzymes of the retinoic acid (RA) synthetic and degradation pathways was characterized in adult female zebrafish Danio rerio. Females were conditioned until maturation and post-spawn expression dynamics were determined. A striking upregulation of cyp26b1, but not cyp26a1, was observed following egg deposition, decreasing to initial levels during recovery. A similar, yet lower, fluctuation was observed for aldh1a2 and rdh10a, the enzymes participating in the two-step RA biosynthesis cascade. The present work highlights the dynamics of the adult D. rerio oogenesis and uncovers novel, yet elusive, metabolic contributors. Possible compartmentalized roles for the different gene paralogue isoforms are discussed.

In vivo zebrafish morphogenesis shows Cyp26b1 promotes tendon condensation and musculoskeletal patterning in the embryonic jaw

Integrated development of diverse tissues gives rise to a functional, mobile vertebrate musculoskeletal system. However, the genetics and cellular interactions that drive the integration of muscle, tendon, and skeleton are poorly understood. In the vertebrate head, neural crest cells, from which cranial tendons derive, pattern developing muscles just as tendons have been shown to in limb and trunk tissue, yet the mechanisms of this patterning are unknown. From a forward genetic screen, we determined that cyp26b1 is critical for musculoskeletal integration in the ventral pharyngeal arches, particularly in the mandibulohyoid junction where first and second arch muscles interconnect. Using time-lapse confocal analyses, we detail musculoskeletal integration in wild-type and cyp26b1 mutant zebrafish. In wild-type fish, tenoblasts are present in apposition to elongating muscles and condense in discrete muscle attachment sites. In the absence of cyp26b1, tenoblasts are generated in normal numbers but fail to condense into nascent tendons within the ventral arches and, subsequently, muscles project into ectopic locales. These ectopic muscle fibers eventually associate with ectopic tendon marker expression. Genetic mosaic analysis demonstrates that neural crest cells require Cyp26b1 function for proper musculoskeletal development. Using an inhibitor, we find that Cyp26 function is required in a short time window that overlaps the dynamic window of tenoblast condensation. However, cyp26b1 expression is largely restricted to regions between tenoblast condensations during this time. Our results suggest that degradation of RA by this previously undescribed population of neural crest cells is critical to promote condensation of adjacent scxa-expressing tenoblasts and that these condensations are subsequently required for proper musculoskeletal integration.

A novel de novo mutation in MYT1, the unique OAVS gene identified so far

Oculo-auriculo-vertebral spectrum (OAVS) is a developmental disorder characterized by hemifacial microsomia associated with ear, eyes and vertebrae malformations showing highly variable expressivity. Recently, MYT1, encoding the myelin transcription factor 1, was reported as the first gene involved in OAVS, within the retinoic acid (RA) pathway. Fifty-seven OAVS patients originating from Brazil were screened for MYT1 variants. A novel de novo missense variant affecting function, c.323C>T (p.(Ser108Leu)), was identified in MYT1, in a patient presenting with a severe form of OAVS. Functional studies showed that MYT1 overexpression downregulated all RA receptors genes (RARA, RARB, RARG), involved in RA-mediated transcription, whereas no effect was observed on CYP26A1 expression, the major enzyme involved in RA degradation, Moreover, MYT1 variants impacted significantly the expression of these genes, further supporting their pathogenicity. In conclusion, a third variant affecting function in MYT1 was identified as a cause of OAVS. Furthermore, we confirmed MYT1 connection to RA signaling pathway.

Toxicity of triphenyltin on the development of retinal axons in zebrafish at low dose

The impacts of triphenyltin (TPT) on ecological health are of particular concern due to the unexpectedly high levels found in wild fish around the world. Here, zebrafish embryos were exposed to TPT via in ovo nano-injection to study its toxicity on the development of retinal axons in fish. Lipophilic dye labeling revealed obvious defects in retinal axon development in larvae with normally shaped eyes, with incidences of 0, 1.08%, 2.66%, 4.26%, and 6.85% observed in the control, 0.8, 4.0, 20.0, and 100ng TPT-Cl/g wet weight (ww) exposure groups, respectively, showing a dose-dependent increase. Since the lowest observable effective concentration of TPT to induce retinal axon development defects was 0.8ng TPT-Cl/g ww, which is lower than the concentrations in wild fish eggs, this defect would occur in wild fish larvae. Alterations in the expressions of pax6 and ephrinBs, which regulate the establishment of retinal polarity, were correlated with defect incidence. Expression levels of the CYP26A1 gene and protein were significantly up-regulated in all exposure groups compared with the control, which may lead to significant decreases in concentrations of all-trans retinoic acid (atRA). Such a disruption of RA metabolism would, at least partly, contribute to the incidence of developmental defects in retinal axons. This study is the first to report that TPT can interfere with development of retinal axons in fish at low dose.

Diverse Regulation of Vitamin D Receptor Gene Expression by 1,25-Dihydroxyvitamin D and ATRA in Murine and Human Blood Cells at Early Stages of Their Differentiation

Vitamin D receptor (VDR) is present in multiple blood cells, and the hormonal form of vitamin D, 1,25-dihydroxyvitamin D (1,25D) is essential for the proper functioning of the immune system. The role of retinoic acid receptor α (RARα) in hematopoiesis is very important, as the fusion of RARα gene with PML gene initiates acute promyelocytic leukemia where differentiation of the myeloid lineage is blocked, followed by an uncontrolled proliferation of leukemic blasts. RARα takes part in regulation of VDR transcription, and unliganded RARα acts as a transcriptional repressor to VDR gene in acute myeloid leukemia (AML) cells. This is why we decided to examine the effects of the combination of 1,25D and all-trans-retinoic acid (ATRA) on VDR gene expression in normal human and murine blood cells at various steps of their development. We tested the expression of VDR and regulation of this gene in response to 1,25D or ATRA, as well as transcriptional activities of nuclear receptors VDR and RARs in human and murine blood cells. We discovered that regulation of VDR expression in humans is different from in mice. In human blood cells at early stages of their differentiation ATRA, but not 1,25D, upregulates the expression of VDR. In contrast, in murine blood cells 1,25D, but not ATRA, upregulates the expression of VDR. VDR and RAR receptors are present and transcriptionally active in blood cells of both species, especially at early steps of blood development.

Perturbation of Retinoid Homeostasis Increases Malformation Risk in Embryos Exposed to Pregestational Diabetes

Pregestational diabetes is highly associated with an increased risk of birth defects. However, factors that can increase or reduce the expressivity and penetrance of malformations in pregnancies in women with diabetes remain poorly identified. All-trans retinoic acid (RA) plays crucial roles in embryogenesis. Here, we find that Cyp26a1, which encodes a key enzyme for catabolic inactivation of RA required for tight control of local RA concentrations, is significantly downregulated in embryos of diabetic mice. Embryonic tissues expressing Cyp26a1 show reduced efficiency of RA clearance. Embryos exposed to diabetes are thus sensitized to RA and more vulnerable to the deleterious effects of increased RA signaling. Susceptibility to RA teratogenesis is further potentiated in embryos with a preexisting genetic defect of RA metabolism. Increasing RA clearance efficiency using a preconditioning approach can counteract the increased susceptibility to RA teratogenesis in embryos of diabetic mice. Our findings provide new insight into gene-environment interactions that influence individual risk in the manifestation of diabetes-related birth defects and shed light on environmental risk factors and genetic variants for a stratified medicine approach to screening women with diabetes who are of childbearing age and assessing the risk of birth defects during pregnancy.

Retinoic acid catabolizing enzyme CYP26C1 is a genetic modifier in SHOX deficiency

Mutations in the homeobox gene SHOX cause SHOX deficiency, a condition with clinical manifestations ranging from short stature without dysmorphic signs to severe mesomelic skeletal dysplasia. In rare cases, individuals with SHOX deficiency are asymptomatic. To elucidate the factors that modify disease severity/penetrance, we studied a three-generation family with SHOX deficiency. The variant p.Phe508Cys of the retinoic acid catabolizing enzyme CYP26C1 co-segregated with the SHOX variant p.Val161Ala in the affected individuals, while the SHOX mutant alone was present in asymptomatic individuals. Two further cases with SHOX deficiency and damaging CYP26C1 variants were identified in a cohort of 68 individuals with LWD The identified CYP26C1 variants affected its catabolic activity, leading to an increased level of retinoic acid. High levels of retinoic acid significantly decrease SHOX expression in human primary chondrocytes and zebrafish embryos. Individual morpholino knockdown of either gene shortens the pectoral fins, whereas depletion of both genes leads to a more severe phenotype. Together, our findings describe CYP26C1 as the first genetic modifier for SHOX deficiency.

Accelerated Skeletal Maturation in Disorders of Retinoic Acid Metabolism: A Case Report and Focused Review of the Literature

Nutritional excess of vitamin A, a precursor for retinoic acid (RA), causes premature epiphyseal fusion, craniosynostosis, and light-dependent retinopathy. Similarly, homozygous loss-of-function mutations in CYP26B1, one of the major RA-metabolizing enzymes, cause advanced bone age, premature epiphyseal fusion, and craniosynostosis. In this paper, a patient with markedly accelerated skeletal and dental development, retinal scarring, and autism-spectrum disease is presented and the role of retinoic acid in longitudinal bone growth and skeletal maturation is reviewed. Genetic studies were carried out using SNP array and exome sequencing. RA isomers were measured in the patient, family members, and in 18 age-matched healthy children using high-performance liquid chromatography coupled to tandem mass spectrometry. A genomic SNP array identified a novel 8.3 megabase microdeletion on chromosome 10q23.2-23.33. The 79 deleted genes included CYP26A1 and C1, both major RA-metabolizing enzymes. Exome sequencing did not detect any variants that were predicted to be deleterious in the remaining alleles of these genes or other known retinoic acid-metabolizing enzymes. The patient exhibited elevated plasma total RA (16.5 vs. 12.6±1.5 nM, mean±SD, subject vs. controls) and 13-cisRA (10.7 nM vs. 6.1±1.1). The findings support the hypothesis that elevated RA concentrations accelerate bone and dental maturation in humans. CYP26A1 and C1 haploinsufficiency may contribute to the elevated retinoic acid concentrations and clinical findings of the patient, although this phenotype has not been reported in other patients with similar deletions, suggesting that other unknown genetic or environmental factors may also contribute.

Gambogic acid causes fin developmental defect in zebrafish embryo partially via retinoic acid signaling

Gambogic acid (GA), the major active ingredient of gamboge, has been approved by the Chinese Food and Drug Administration for clinical trials in cancer patients due to its strong anticancer activity. However, our previous research showed that GA was teratogenic against zebrafish fin development. To explore the teratogenicity and the underlying mechanisms, zebrafish (Danio rerio) embryos were used. The morphological observations revealed that GA caused fin defects in zebrafish embryos in a concentration-dependent manner. The critical exposure time of GA to reveal teratogenicity was before 8 hpf (hours post fertilization). LC/MS/MS analysis revealed that a maximum bioconcentration of GA was occurred at 4 hpf. Q-PCR data showed that GA treatment resulted in significant inactivation of RA signaling which could be partially rescued by the exogenous supply of RA. These results indicate the potential teratogenicity of GA and provide evidence for a caution in its future clinic use.