Skip to main content

ORIGINAL RESEARCH article

Front. Cell Dev. Biol., 27 April 2021
Sec. Molecular and Cellular Pathology

Bhlhe40/Sirt1 Axis-Regulated Mitophagy Is Implicated in All-Trans Retinoic Acid-Induced Spina Bifida Aperta

\r\nLu ZhaoLu ZhaoDan LiuDan LiuWei MaWei MaHui GuHui GuXiaowei WeiXiaowei WeiWenting LuoWenting LuoZhengwei Yuan*Zhengwei Yuan*
  • Key Laboratory of Health Ministry for Congenital Malformation, Shengjing Hospital, China Medical University, Shenyang, China

Neural tube defects (NTDs) are the most severe congenital malformations that result from failure of neural tube closure during early embryonic development, and the underlying molecular mechanisms remain elusive. Mitophagy is the best-known way of mitochondrial quality control. However, the role and regulation of mitophagy in NTDs have not yet been elucidated. In this study, we used an all-trans retinoic acid (ATRA)-induced rat model to investigate mitophagy and its underlying mechanism in spina bifida aperta (SBA). The results of western blot, immunofluorescence and RT-qPCR analyses indicated that mitophagy was impaired and Sirt1 was downregulated in SBA. Administration of resveratrol-a strong specific Sirt1 activator-activated Sirt1, thus attenuating autophagy suppression and ameliorating SBA. RNA-sequencing and bioinformatics analysis results indicated that transcriptional regulation played an important role in NTDs. A luciferase reporter assay was performed to demonstrate that the transcription factor Bhlhe40 directly bound to and negatively regulated Sirt1 expression. Further, we discovered that the Bhlhe40/Sirt1 axis regulated mitophagy in neural stem cells. Collectively, our results for the first time demonstrate that Bhlhe40/Sirt1 axis regulated mitophagy is implicated in ATRA-induced SBA. Our findings provide new insights into pathogenesis of NTDs and a basis for potential therapeutic targets for NTDs.

Introduction

Neural tube defects (NTDs), such as anencephaly and spina bifida, represent the most severe and common congenital malformations of the central nervous system resulting in infant deaths and severe disabilities, bring undue burden to the ill infant and their families as well as society (Greene and Copp, 2014). The global incidence of NTDs is approximately 18.6 in 10,000 births globally (Blencowe et al., 2018). NTDs have a complex etiology and mainly result from incomplete closure of the neural tube during early embryonic development. Failure of neural folding or closure of neuropores is caused by an interaction of genetic and environmental factors, lack of nutrients or other factors (Bhandari and Thada, 2020). Current treatment strategies for NTDs are not satisfactory because the pathogenesis is insufficiently understood (Norkett et al., 2016; Radic et al., 2019). In-depth study of the mechanisms underlying NTD pathogenesis and the development of scientific prevention and treatment methods are of great significance for reducing the malformation and disability rates, thus improving patient quality of life and reducing societal burden.

Neural tube defects are considered polygenic genetic diseases. Extensive studies of vertebrate models have partially elucidated the mechanism of neural tube closure and its failure during NTDs. Studies in various animal models have indicated that the pathogenesis of NTDs involves multiple molecular pathways (Juriloff and Harris, 2018). Neural fold bending, elevation, and fusion is regulated via various molecular events, including the non-canonical wnt/planar cell polarity pathway, sonic hedgehog/bone morphogenetic protein signaling, and several transcription factors regulation (Nikolopoulou et al., 2017; Wang et al., 2019). Recently, it has been found that the dysregulation of apoptosis and autophagy in fetal central nervous tissues is associated with NTD development. Mutations in apoptosis-related genes CASP3, CASP9, and APAF1 reportedly are conductive to human NTDs (Zhou et al., 2018). In another recent study, deletion of the forkhead box O3a transcription factor gene eliminated the inhibitory effect of maternal diabetes on the expression of autophagy-related protein and thus alleviated diabetes-induced cell apoptosis and NTDs in a mouse model (Xu et al., 2019). Further research is still warranted to reveal the mechanisms behind NTDs.

Mitochondria are essential organelles that provide cells with energy and buffering cytoplasmic calcium. In addition, they are the sites of harmful reactive oxygen species (ROS) production and secrete pro-apoptotic factors (Giorgi et al., 2012). A growing body of evidences suggested that neurodevelopment requires well-structured mitochondria to maintain nervous system function (Khacho and Slack, 2018). Notably, mitochondrial damage recently has been proven to be associated with NTDs (Chen X. et al., 2017). Maternal diabetes can significantly increase the number of defective mitochondria in the neuroepithelia of mouse embryos, but does not affect the total number of mitochondria (Xu et al., 2013). Our previous study revealed that compared with normal rat embryos, mitochondrial respiratory efficiency and the adenosine triphosphate content were decreased in rat embryos with all-trans retinoic acid (ATRA)-induced spina bifida aperta (SBA), indicating that mitochondrial dysfunction is associated with ATRA -induced SBA in rat (Xue et al., 2018).

Autophagy is a major cellular pathway that facilitates the elimination of the toxic proteins and damaged organelles by lysosomal machinery. Previous studies have identified a linear association between autophagy dysfunction and NTDs (Xu et al., 2019; Ye et al., 2020). Mitophagy is a selective form of autophagy that mediates the elimination of defective mitochondria and is the only known way to selectively remove entire mitochondria (Lou et al., 2020) and control mitochondrial quality (Liu et al., 2014). Promotion of mitophagy could alleviated the development of neurological deficits due to traumatic brain injury in mice (Lin et al., 2020). Autophagy inducer rapamycin promoted mitophagy and exerted neuroprotective effects in spinal cord ischemia-reperfusion injury mouse model (Li et al., 2018). Therefore, we hypothesized that mitophagy may play an important role in the pathogenesis of NTDs and elucidating its role and underlying molecular mechanisms may reveal potential therapeutic targets.

Master transcription factors regulate autophagy and lysosomal gene expression by coordinating the expansion of autophagy-lysosomal fitness for mitophagy (Fang et al., 2016; Tan and Wong, 2017). These transcription factors include transcription factor EB and forkhead box O, which are activated by silent information regulator factor 2-related enzyme 1 (Sirt1) (Brunet et al., 2004; Ivankovic et al., 2016). Sirt1, best-studied member of the sirtuin family, is a critical regulator of autophagy (Tang, 2016). Recently, studies have shown that Sirt1 is vital for the disposal of defective mitochondria by mitophagy (Fang, 2019; Wang et al., 2020). In mice, environmental stresses reportedly can activate Sirt1 and lead to NTD-like phenotypes (Li et al., 2017). However, a microarray study demonstrated that Sirt1 was downregulated in a valproic acid-induced NTD mouse model (Okada et al., 2005). Further, Sirt1 was decreased in embryonic neural stem cells in diabetic pregnancy due to upregulation of its regulator, miRNA-30b (Ramya et al., 2017). Given these inconsistent results, more in-depth research is urgently needed to elucidate the role of Sirt1 in NTDs.

In the present study, we used an ATRA-induced SBA rat model to investigate mitophagy and Sirt1 expression. RNA-sequencing (RNA-seq) and bioinformatics analysis were performed to find the underlying molecular regulatory mechanism in SBA. Neural stem cells were used to verify the mitophagy and its regulatory mechanism. The aim of our study is to offer new insights into the pathogenesis underlying NTDs.

Materials and Methods

Experimental Animals and Treatments

Outbred female Wistar rats (240–300 g) were obtained from the Animal Center of Shengjing Hospital associated with China Medical University (Shenyang, China) at 10–12 weeks of age. On embryonic day (E) 10, pregnant rats were given a single intragastric dose of ATRA (4% [wt/vol] in olive oil; 140 mg/kg body weight) (Sigma-Aldrich, United States) or control treatment (olive oil only). Resveratrol (RSV; trans-3,5,4′-trihydroxystilbene) (MedChemExpress, NJ, United States; 40 mg/kg body weight), suspended in 50% saline, 5% DMSO, 5% Tween80, and 40% PEG400, was injected intraperitoneally once a day from E8 to E17. From E11 to E19, the posterior spinal cords of embryos were isolated, immediately snap-frozen, and stored at −80 °C until use for detecting gene and protein expression detection. Specimens for histological analysis were isolated, then fixed in 4% paraformaldehyde for 24 h. All animal experiments were approved by the Medical Ethics Committee of Shengjing Hospital of China Medical University.

Cell Culture and Treatments

Mouse C17.2 neural stem cells were maintained in minimum essential medium (Gibco-Invitrogen, Carlsbad, CA, United States) supplemented with 10% fetal bovine serum (Gibco), 1% penicillin/streptomycin (100 U/mL and 100 μg/mL, respectively; Life Technologies), and 1% minimum essential medium non-essential amino acids (Gibco) in a humidified environment at 37°C with 5% CO2.

C17.2 cells were seeded in 6-well plates (1 × 105 cells/well) and cultured before transfection. Plasmids pcDNA3.1(+)-CMV-Bhlhe40 and pcDNA3.1(+)-CMV-MCS were purchased from SyngenTech (Beijing, China). The cells were transfected with 1.5 μg plasmids or 100 pmol small interfering (si)RNAs in 2 mL of culture medium per well, using LipoFiterTM liposomal transfection reagent (Hanbio Technology, China). After transfection for 48 h, the cells were harvested.

For Sirt1 activation or inhibition, RSV (a specific activator of Sirt1; 5 μM in DMSO) or EX527 (a specific inhibitor of Sirt1; 5 μM in DMSO; MedChemExpress) was added to the cell culture for 12 h.

Transmission Electron Microscopy (TEM)

Mitochondrial structures were examined by TEM. Briefly, embryonic lumbo-sacral neural tube on E12 were collected carefully, immediately fixed in 2.5% glutaraldehyde, and post-fixed in 1% osmium tetraoxide. Neuroepithelia were cut and viewed under a H-7650 transmission electron microscope (Hitachi, Japan) operated at 80 KV.

Mitochondrial Isolation

MinuteTM mitochondrial extraction kits (Invent Biotechnologies, Inc., Plymouth, MN, United States) were used for mitochondrial isolation of C17.2 cells according to the manufacturer’s instructions. The mitochondria were isolated through differential centrifugation and kept in a storage solution with MinuteTM denatured protein solution for western blotting analysis.

Western Blotting Analysis

Total protein was extracted from rat embryonic spinal cords (E11 and E12) or C17.2 cells using radioimmunoprecipitation (RIPA) buffer (Solarbio Science & Technology, Beijing, China). Protein concentrations were determined by a bicinchoninic acid assay (Solarbio Science & Technology). Proteins were separated by 8%–12% sodium dodecyl sulfate-polyacrylamide gel electrophoresis and electrophoretically transferred to polyvinylidene difluoride (PVDF) membranes. After blocking with 5% non-fat milk solution, the membranes were incubated with primary antibodies at 4°C overnight. Following primary antibodies were applied in the experiment: anti-LC3 (Abcam, Cambridge, MA, United States; 1:1000), anti-LAMP2 (ProteinTech, Chicago, IL, United States; 1:1000), anti-Sirt1 (Cell Signaling Technology, Boston, MA, United States; 1:1000), anti-Bhlhe40 (ProteinTech, 1:1000), anti-PINK1 (ProteinTech, 1:1000), anti-VDAC1 (Abcam, 1:2000) and anti-GAPDH (ProteinTech, 1:5000). On the following day, the membranes were washed, and incubated with corresponding secondary antibodies (ProteinTech; 1:5000) at room temperature for 1.5 h. Enhanced chemiluminescent (ECL) reagent (Millipore, Billerica, MA, United States) was applied to detect protein-antibody interactions, which were visualized using a chemiluminescence detection system (C300, Azure Biosystems, United States).

Immunofluorescence Staining

Rat embryos on E12 were collected, then fixed in 4% paraformaldehyde and paraffin embedded. Paraffin-embedded embryos were cut into 4-μm-thick sections and mounted on slides. The sections were deparaffinized in xylene, next hydrated in a graded series of ethanol. C17.2 cells were harvested and fixed with 4% paraformaldehyde, and incubated in 0.5% Triton for 30 min. The sections and cells were incubated with 5% bovine serum albumin (BSA) for 45 min in order to block non-specific sites. Then, they were incubated with normal goat serum, followed by the primary antibodies mouse anti-TOMM20 (Abcam, 1:100) and rabbit anti-LC3 (Abcam, 1:100) in a humidified chamber at 4°C overnight. After sections and cells were incubated with secondary antibodies conjugated with Alexa Fluor 488 (Cell Signaling Technology, 1:100) and Alexa Fluor 555 (Cell Signaling Technology, 1:100) were added and incubated at room temperature in the dark for 2 h. Nuclei were counterstained with 4′,6-diamidino-2-phenylindole. The specimens were observed under an ECLIPSE 80i fluorescence microscope (Nikon, Japan) or a laser scanning confocal microscopy (Ziess, Germany).

Total RNA Isolation and Quantitative Reverse Transcription PCR (RT-qPCR)

Total RNA was extracted from rat embryonic spinal cords (from E11 to E19) or C17.2 cells using TRIzol reagent (Life Technologies, Carlsbad, CA, United States). The RNA was reverse-transcribed into cDNA using a PrimeScript RT kit (TaKaRa Biotech. Co., Ltd., China). RT-qPCR was run in an ABI 7500 Real-time PCR system (Applied Biosystems, United States) using a SYBR green master mix kit (TaKaRa). Primers for rat or mouse genes were designed and synthesized by Sangon Biotech (Sangon Biotech, Shanghai, China) (Table 1). Relative target gene expression was quantified by the 2−ΔΔCt method, using Gapdh as an internal control.

TABLE 1
www.frontiersin.org

Table 1. Sequences of the primers used for RT-qPCR analysis.

RNA-Seq

Total RNA was extracted from E11 and E12 rat embryonic spinal cord tissues. RNA concentrations were determined using an Agilent 2100 Bioanalyzer (Agilent RNA 6000 Nano Kit, United States). RNA purity was evaluated using a NanoDrop instrument. The BGISEQ500 platform (BGI-Shenzhen, China) was used for RNA-seq. The following 12 tissue samples were analyzed: control group of E11 (E11CON1, E11CON2, E11CON3), SBA group of E11 (E11SBA1, E11SBA2, E11SBA3), control group of E12 (E12CON1, E12CON2, E12CON3), SBA group of E12 (E12SBA1, E12SBA2, E12SBA3). Perform quality control on the raw reads to determine whether the sequencing data were fit for follow-up analysis. In total 11.30G RNA-seq data was obtained.

Bioinformatic Analysis

In brief, differentially expressed genes (DEGs) were identified using DESeq2 (v1.4.5) with Q-value ≤ 0.05 (Love et al., 2014). To gain insight into DEG functions, Gene Ontology (GO) and transcription factor enrichment analyses were performed using Phyper, based on a hypergeometric test. Significantly enriched terms were Bonferroni-corrected (Q-value ≤ 0.05) by Bonferroni. The top 20 terms are presented.

Luciferase Reporter Assay

Wild-type and mutant sequences of Sirt1 mRNA sequences containing the putative Bhlhe40-binding sites were synthesized, and cloned into the pGL4.10 vector (Beijing SyngenTech). C17.2 cells were co-transfected with pcDNA3.1(+)-CMV-Bhlhe40 or pcDNA3.1(+)-CMV-MCS and the wild-type or mutant constructs. After incubating for 48 h, the C17.2 cells were harvested and subjected to a dual-luciferase reporter assay system (Promega, Madison, WI, United States) to determine binding between Sirt1 and Bhlhe40.

Chromatin Immunoprecipitation (ChIP) Assay

ChIP assays were carried out essentially following the protocol from the SimpleChIP Enzymatic Chromatin IP Kit (Cell Signaling Technology, 9002). Briefly, C17.2 cells were treated with 1% formaldehyde for 10 min at room temperature to crosslink proteins to DNA. The chromatin was harvested and fragmented using enzymatic digestion. An aliquot of each sample was set aside as input control, while the remaining portion was subjected to immunoprecipitation with anti-Bhlhe40 antibodies (Novus Biologicals Co., Ltd., Shanghai, China) overnight at 4°C. IgG was used as a control for non-specific immunoprecipitation of DNA. Immunoprecipitated complex was treated with protease and the fragment of the Sirt1 promoter containing the E-box was amplified by RT-qPCR. The sequences of the primers were mouse Sirt1-promoter-F: 5′-GCGCCATCGCAAACTTGA-3′ and mouse Sirt1-promoter-R: 5′-CACAACACGCCGGGTCA-3′.

Statistical Analysis

All results were analyzed with SPSS 17.0 software (SPSS Inc., Chicago, IL, United States). Data are expressed as the mean ± standard deviations (SD). Means of two groups were compared using two-tailed Student’s t-tests. The significant differences in Sirt1 expression among the different treatment groups were performed by one-way analysis of variance (ANOVA) test. p < 0.05 was considered statistically significant.

Results

Mitophagy in the Spinal Cord Is Inhibited in Embryonic Rats With SBA

Mitochondrial morphology in the neuroepithelium of spinal cords of rat embryos at E12 was detected by TEM. As shown in Figure 1A, mitochondrial outer membrane wrapped the mitochondria, which were horizontally arranged in ridges. In the SBA group, the mitochondrial ridges were destroyed or even disappeared, and the matrix electron density was reduced. We found that compared to the control group, the number of defective mitochondria was significantly increased in the SBA group (Figure 1B).

FIGURE 1
www.frontiersin.org

Figure 1. Autophagy and mitophagy in SBA rat embryos. (A) TEM detection of mitochondria in the neuroepithelium of spinal cords of rat embryos of the CON and SBA groups. Arrows indicate mitochondria. Scale bar = 500 nm. (B) Defective mitochondria rates in neuroepithelial cells from E12 embryos. Defective mitochondria rate = number of defective mitochondria divided by total number of mitochondria. (C) Western blot showing LC3 expression in rat embryonic spinal cord tissues in the CON and SBA groups. (D) Quantitation of LC3 II expression (n = 6). (E) Western blot showing LAMP2 expression in rat embryonic spinal cord tissues in the CON and SBA groups. (F) Quantitation of LAMP2 expression (n = 6). (G) Double immunofluorescence staining with LC3-labeled autophagosomes (green) and TOMM20-labeled mitochondria (red) in rat embryonic spinal cords. Scale bar = 10 μm. (H) Relative fluorescence intensities of LC3 and TOMM20 in the SBA group versus the CON group (n = 4). Results are expressed as the mean ± SD. *p < 0.05 vs. control group.

Next, the expression of the autophagy marker microtubule associated protein 1 light chain 3 (LC3) and the lysosomal-specific marker lysosome-associated membrane protein 2 (LAMP2) was investigated in ATRA-induced SBA and control rats by western blotting. LC3 II expression was significantly decreased in the SBA group as compared to the CON group on E12 (Figures 1C,D). LAMP2 protein expression was notably decreased by 44.85% in the SBA group as compared to the CON group (Figures 1E,F).

To determine the occurrence of mitophagy in the ATRA-induced SBA rat model, colocalization of LC3-stained autophagosomes and translocase of outer mitochondrial membrane 20 (TOMM20)-stained mitochondria in the neuroepithelium of spinal cords was analyzed. As shown in Figures 1G,H, consistent with the western blotting results, LC3 staining was markedly weaker in the SBA group than in the CON group. TOMM20 antibody staining were relatively uniform, and no significant difference was noted between the CON and SBA group. Further, there was a lower level of LC3 puncta with TOMM20-stained mitochondria were observed in ATRA-induced SBA tissues, indicating mitophagy inhibition.

Sirt1 Expression in the Spinal Cord Is Downregulated in Embryonic Rats With SBA

Next, we investigated Sirt1 expression in the spinal cord in SBA model and control rats. RT-qPCR results indicated that Sirt1 mRNA was decreased in the SBA group as compared to the CON group at each time point between E11 and E19 (Figure 2A). To explore whether Sirt1 protein levels were consistent with the mRNA levels, samples collected during the critical period of neural tube closure (E11 and E12) were selected for western blot analysis. Sirt1 protein was significantly downregulated in the SBA group as compared to the CON group on E11 (Figures 2B,C) and E12 (Figures 2D,E). There results suggested that Sirt1 may be involved in NTDs formation.

FIGURE 2
www.frontiersin.org

Figure 2. Sirt1 expression in ATRA-induced SBA. (A) Relative Sirt1 mRNA expression in rat embryonic spinal cord tissues of the CON and SBA groups from E11 to E19 as assessed by RT-qPCR (n = 5). Results are expressed as the mean ± SD. *p < 0.05 vs. control group. (B) Western blot of Sirt1 in rat embryonic spinal cord tissues on E11 in the CON and SBA groups. (C) Quantitation of Sirt1 protein on E11 (n = 6). (D) Western blot of Sirt1 in rat embryonic spinal cord tissues on E12 in the CON and SBA groups. (E) Quantitation of Sirt1 protein on E12 (n = 6). Results are expressed as the mean ± SD. *p < 0.05 vs. control group.

RSV Increases Sirt1 and Ameliorates Autophagy Suppression in vivo

To confirm the function of Sirt1 in ATRA-induced SBA, pregnant dams were treated with RSV (a Sirt1 agonist) from E8 to E17. As shown in Figures 3A,B, Sirt1 expression was decreased in ATRA-treated rats as compared to CON rats on E18. Sirt1 levels in embryonic spinal cord tissues were significantly elevated after treatment with RSV. The protein level of the mitophagic marker PINK1 was also decreased in ATRA-treated rats as compared to CON rats on E18, and RSV significantly enhanced PINK1 expression as compared to ATRA treatment alone (Figures 3C,D). Further, LC3 II expression in the embryonic spinal cord tissues was significantly reduced in ATRA-treated group as compared to the CON group on E18, whereas RSV significantly increased the levels of LC3 II as compared to the ATRA treatment alone on E18 (Figures 3E,F). LAMP2 protein levels in the embryos were significantly decreased by approximately 42.98% in the ATRA-treated group as compared to the CON group on E18, whereas RSV significantly increased these levels as compared to ATRA treatment alone on E18 (Figures 3G,H).

FIGURE 3
www.frontiersin.org

Figure 3. RSV activates Sirt1 and ameliorates autophagy suppression in rat embryos with ATRA-induced SBA. (A) Western blot of Sirt1 protein expression in rat embryonic spinal cord tissues in the CON, ATRA, and ATRA + RSV groups. (B) Quantitation of Sirt1 protein (n = 6). (C) Western blot of PINK1 protein expression in rat embryonic spinal cord tissues in the CON, ATRA, and ATRA + RSV groups. (D) Quantitation of PINK1 protein (n = 6). (E) Western blot of LC3 protein expression in rat embryonic spinal cord tissues in the CON, ATRA, and ATRA + RSV groups. (F) Quantitation of LC3 II protein (n = 6). (G) Western blot of LAMP2 protein in rat embryonic spinal cord tissues in the CON, ATRA, and ATRA + RSV groups. (H) Quantitation of LAMP2 protein (n = 6). The CON group tissues were taken from normal embryos, ATRA group and ATRA + SRV group tissues were taken from embryos with SBA phenotypes. Results are expressed as the mean ± SD. p < 0.05 vs. CON group; #p < 0.05 vs. ATRA group.

Importantly, RSV treatment reduced the ATRA-induced SBA incidence rate in rat embryos (Table 2). These results implied that RSV treatment enhanced autophagy/mitophagy and reduced the SBA rate in ATRA-induced SBA rats.

TABLE 2
www.frontiersin.org

Table 2. RSV treatment reduces the SBA rate.

Transcriptome Alterations in ATRA-Induced SBA Rats

To explore the mechanism of NTDs and the genes regulating Sirt1 in NTDs, RNA-seq was used to identify DEGs between CON and SBA rat embryos on E11 and E12. As shown in Figure 4A, we identified DEGs, including 638 upregulated and 377 downregulated genes, in the E11SBA group compared to the E11CON group, and 872 DEGs, including 507 upregulated and 365 downregulated genes, in the E12SBA group compared to the E12CON group. The Venn diagram in Figure 4B shows that 324 DEGs were in common on E11 and E12 (Figure 4B).

FIGURE 4
www.frontiersin.org

Figure 4. Transcriptome alterations in the ATRA-induced SBA rat model (E11 and E12). (A) Number of DEGs between SBA and CON rat embryos on E11 and E12. DESeq2 | log2FC| ≥ 0.3, Q ≤ 0.05. (B) Venn diagram of DEGs in E11CON vs. E11SBA and E12CON vs. E12SBA. (C) Top 20 enriched BP terms in GO analysis. (D) Top 20 enriched CC terms in GO analysis. (E) Top 20 enriched MF terms in GO analysis. (F) Transcription factor enrichment analysis of DEGs. Terms were considered significantly enriched when Q ≤ 0.05.

Gene Ontology analysis was conducted to speculate the potential functions of those 324 DEGs. GO analysis describes genes in terms of biological process (BP) (Supplementary Table 1), cellular component (CC) (Supplementary Table 2) and molecular function (MF) (Supplementary Table 3). The top 20 GO processes ranked by enrichment score [−log10 (Q-value)] are presented. The most enriched GO terms included positive regulation of transcription by RNA polymerase II, regulation of transcription, DNA-templated, multicellular organism development and somitogenesis in BP (Figure 4C); extracellular space, extracellular matrix, receptor complex and transcription factor complex in CC (Figure 4D); sequence-specific DNA binding, DNA-binding transcription activator activity, RNA polymerase II-specific, RNA polymerase II proximal promoter sequence-specific DNA binding and DNA-binding transcription factor activity, RNA polymerase II-specific in MF (Figure 4E). Then, transcription factor enrichment analysis of the 324 DEGs further revealed that differential genes are mainly enriched in the basic helix-loop-helix (bHLH) transcription factor family was strongly enriched (Figure 4F) (Supplementary Tables 4, 5). Together, these results suggested that transcription factors are strongly involved in the regulation of ATRA-induced SBA.

Bhlhe40 Expression Is Increased in the Spinal Cord in Embryonic Rats With SBA

We selected bHLH transcription factor family genes that regulate neurodevelopment (Bhlhe40, Atoh1, Neurod1, Nhlh1, Olig3, Mesp2 and Msgn1) for RT-qPCR verification. Promoter region prediction using the JASPAR database1 indicated that these genes harbor potential Sirt1-binding sites. The average fragments per kilobase of transcript per million fragments mapped (FPKM) values of these bHLH transcription factor genes in SBA rat embryos are presented in Figure 5A.

FIGURE 5
www.frontiersin.org

Figure 5. Expression of bHLH transcription factors in spinal cords of rat embryos with SBA-induced ATRA. (A) Average fragments per kilobase of transcript per million fragments mapped values of transcription factors in SBA rat embryos (n = 3). (B–H) Relative mRNA expression of Bhlhe40, Atoh1, Neurod1, Nhlh1, Olig3, Mesp2, and Msgn1 in rat embryonic spinal cord tissues in the CON and SBA groups on E11 and E12 as assessed by RT-qPCR (n = 6). Results are expressed as the mean ± SD. p < 0.05 vs. control group. (I) Western blot of Bhlhe40 expression. (J) Quantitation of Bhlhe40 expression (n = 6). Results are expressed as the mean ± SD. p < 0.05 vs. control group. (K) Relative Bhlhe40 mRNA expression in rat embryonic spinal cord tissues of the CON and SBA groups from E11 to E19 as assessed by RT-qPCR (n = 5). Results are expressed as the mean ± SD. p < 0.05 vs. control group.

RT-qPCR results showed that among the upregulated genes, the mRNA level of Bhlhe40 was increased 1.97-folds and 1.95-folds in the SBA group as compared to the CON group on E11 and E12, respectively (Figure 5B). The mRNA expression of Mesp2 and Msgn1 was downregulated in SBA, which was consistent with the RNA-seq results (Figures 5G,H). However, there were no significant changes in Nhlh1, Olig3, Atoh1, and Neurod1 expression in the SBA group as compared to the CON group on both E11 and E12 (Figures 5C–F). As Bhlhe40 (encoding basic helix-loop-helix family member e40) was significantly upregulated in the spinal cord in SBA embryos at an early stage and it was predicted to binding to Sirt1 with the highest, we also detected the protein expression of Bhlhe40. Compared to the CON group, protein expression of Bhlhe40 was significantly increased in the SBA group as compared to CON group (Figures 5I,J). RT-qPCR results indicated that Bhlhe40 mRNA was significantly increased in the SBA group as compared to the CON group at each time point between E11 and E17. The mRNA level of Bhlhe40 was increased, but has no statistically significant in the SBA group as compared to the CON group at E18 and E19 (Figure 5K).

Bhlhe40 Binds to the Sirt1 Promoter and Suppresses Its Transcription

Promoter region prediction showed that Bhlhe40 bind on the conserved E-box motifs of the Sirt1 promoter (Figure 6A). Luciferase reporter assays confirmed binding between Sirt1 and Bhlhe40 in C17.2 cells. C17.2 cells transfected with Sirt1-promoter-WT-luc and Bhlhe40 showed significant lower luciferase activity than Sirt1-promoter-WT-luc and pcDNA transfected C17.2 cells (Figure 6B). We next performed a ChIP assay to confirm the binding of Bhlhe40 to the Sirt1 promoter in C17.2 cells. The ChIP assay results revealed that Bhlhe40 interacted with the binding sites of the Sirt1 promoter in C17.2 cells (Figure 6C).

FIGURE 6
www.frontiersin.org

Figure 6. Bhlhe40 negatively regulates Sirt1. (A) Schematic representation of Sirt1 mRNA depicting Bhlhe40-binding sites in its promoter. (B) Luciferase reporter assays were conducted using a Sirt1 promoter with deleted E-box motifs after C17.2 cells were transfection with pcDNA or Bhlhe40 (n = 3). Results are expressed as the mean ± SD. *p < 0.05 vs. Sirt1-promoter-WT-luc + pcDNA group. (C) Binding of Bhlhe40 to Sirt1 promoters was examined in C17.2 by ChIP assay. Results are expressed as the mean ± SD. *p < 0.05 vs. IgG group. (D) RT-qPCR analysis of the Bhlhe40 mRNA level in C17.2 cells transfected with pcDNA or Bhlhe40 (n = 3). (E) RT-qPCR analysis of the Sirt1 mRNA level in C17.2 cells transfected with pcDNA or Bhlhe40 (n = 3). (F) Western blot analysis of Bhlhe40 and Sirt1 protein levels in C17.2 cells transfected with pcDNA or Bhlhe40. (G) Quantitation of Bhlhe40 (n = 3). (H) Quantitation of Sirt1 (n = 3). Results are expressed as the mean ± SD. *p < 0.05 vs. pcDNA group. (I) RT-qPCR analysis of the Bhlhe40 mRNA level in C17.2 cells transfected with NC or si-Bhlhe40 (n = 3). (J) RT-qPCR analysis of the Sirt1 mRNA level in C17.2 cells transfected with NC or si-Bhlhe40 (n = 3). (K) Western blot analysis of Bhlhe40 and Sirt1 protein levels in C17.2 cells transfected with NC or si-Bhlhe40. (L) Quantitation of Bhlhe40 (n = 3). (M) Quantitation of Sirt1 (n = 3). Results are expressed as the mean ± SD. p < 0.05 vs. NC group.

We next transfected C17.2 cells with Bhlhe40 overexpression plasmid or Bhlhe40 siRNA to study the effect of Bhlhe40 on Sirt1. Bhlhe40 overexpression remarkably increased Bhlhe40 mRNA (Figure 6D) and Bhlhe40 protein (Figures 6H,J) expression. Overexpressed Bhlhe40 notably decreased Sirt1 mRNA (Figure 6E) and Sirt1 protein (Figures 6H,K) expression levels. Moreover, Bhlhe40 siRNA significantly decreased the mRNA (Figure 6F) and protein (Figures 6I,L) levels of Bhlhe40, and significantly increased the mRNA (Figure 6G) and protein (Figures 6I,M) levels of Sirt1. These findings indicated that Bhlhe40 regulates Sirt1 transcription by directly binding to its promoter and that Bhlhe40 and Sirt1 expression is negatively correlated.

The Bhlhe40/Sirt1 Axis Regulates Mitophagy in vitro

To investigate the regulatory role of the Bhlhe40/Sirt1 axis in mitophagy function, we used western blotting to examine autophagy-related proteins and double immunofluorescence staining to detect mitophagy in C17.2 cells. Bhlhe40 overexpression significantly decreased the protein levels of Sirt1 (Figures 7A–C), the autophagy marker LC3 II (Figures 7A,D), and the lysosomal-specific marker LAMP2 (Figures 7A,E), which were rescued by RSV treatment. Inversely, Bhlhe40 siRNA significantly increased the protein levels of Sirt1 (Figures 7F–H), LC3 II (Figures 7F,I), and LAMP2 (Figures 7F,J), which were abolished by the Sirt1 inhibitor EX527.

FIGURE 7
www.frontiersin.org

Figure 7. The Bhlhe40/Sirt1 axis regulates autophagy/mitophagy-related proteins in C17.2 cells. (A) Western blot analysis of Bhlhe40, Sirt1, LC3, and LAMP2 protein levels in C17.2 cells treated with different combinations of pcDNA, Bhlhe40, DMSO, and RSV. Quantitation of Bhlhe40 (B), Sirt1 (C), LC3 II (D), and LAMP2 (E) protein (n = 3). Results are expressed as the mean ± SD. p < 0.05 vs. pcDNA + DMSO group; #p < 0.05 vs. Bhlhe40 + DMSO group. (F) Western blot analysis of Bhlhe40, Sirt1, LC3 II, and LAMP2 protein in C17.2 cells treated with different combinations of NC, si-Bhlhe40, DMSO, and EX527. Quantitation of Bhlhe40 (G), Sirt1 (H), LC3 (I), and LAMP2 (J) protein (n = 3). Results are expressed as the mean ± SD. p < 0.05 vs. NC + DMSO group; #p < 0.05 vs. si-Bhlhe40 + DMSO group.

Consistent with the western blotting results, LC3 staining was weakened in the Bhlhe40 overexpression group compared with the pcDNA group, whereas it was increased by RSV treatment (Figures 8A,B). Further, we observed a lower level of colocalization of stained LC3 puncta with TOMM20-stained mitochondria, indicating mitophagy inhibition. The Sirt1 activator RSV attenuated mitophagy inhibition caused by Bhlhe40 in C17.2 cells. PINK1 is characterized by expressing on damaged mitochondria, which regulates Parkin trafficking to phagosomes and drives dysfunctional mitochondria removal via mitophagy (Gelmetti et al., 2017). Western blotting analysis also showed PINK1 was decreased by Bhlhe40 overexpression, whereas increased after RSV treatment in mitochondria of C17.2 (Figures 8C,D). Together, these data implied that Bhlhe40 regulates mitophagy in C17.2 cells by targeting Sirt1.

FIGURE 8
www.frontiersin.org

Figure 8. The Sirt1 activator RSV attenuated mitophagy inhibition caused by Bhlhe40 in C17.2 cells. (A) C17.2 cells were double immunofluorescence-stained, with LC3-labeled autophagosomes (green) and TOMM20-labeled mitochondria (red). Scale bar = 25 μm. (B) Relative fluorescence intensities of LC3 and TOMM20 in each group (n = 3). Results are expressed as the mean ± SD. *p < 0.05 vs. pcDNA + DMSO group; # p < 0.05 vs. Bhlhe40 + DMSO group. (C) Western blot analysis of PINK1 in mitochondria isolated from C17.2 cells treated with different combinations of pcDNA, Bhlhe40, DMSO, and RSV. (D) Quantitation of PINK1 protein (n = 3). Results are expressed as the mean ± SD. p < 0.05 vs. pcDNA + DMSO group; #p < 0.05 vs. Bhlhe40 + DMSO group.

Discussion

Despite advances in the elucidation of the pathogenesis and in the treatment of NTDs, current treatments are still limited and unsatisfactory. In this study, we made three major findings. First, we observed mitophagy impairment in ATRA-induced rat embryonic SBA. Second, administration of RSV activated Sirt1 to attenuate autophagy suppression and ameliorate SBA in vivo. Third, we demonstrated for the first time that the Bhlhe40/Sirt1 axis regulates mitophagy function in vitro.

The maintenance of a healthy and functioning mitochondrial network is critical during development (Pickles et al., 2018). Mitochondria are complex organelles that serve as the cellular powerhouse and are involved in a wide range of key cellular functions, including redox signaling, and the regulation of metabolic homeostasis, cells survival, and cell death (Lou et al., 2020). Emerging evidence suggests that mitochondria are critical regulatory organelles in determining neural stem cell fate, both in initial neural development and in adult neurogenesis (Khacho et al., 2019). Our previous study have showed that mitochondrial respiratory efficiency, adenosine triphosphate content, as well as antioxidant enzyme activity were decreased in spinal cords of early embryos with SBA, which may have caused high ROS levels, leading to mitochondrial dysfunction (Xue et al., 2018). In this study, we observed mitochondrial crest blurring and mitochondrial swelling in early embryos with SBA by TEM, which suggested that damaged mitochondria may not be effectively eliminated.

Autophagy, a lysosome dependent pathway that leads to the degradation of dysfunctional organelles and proteins, plays a crucial role in early embryonic development (Tsukamoto et al., 2008). LC3 is a soluble protein widely distributed in mammalian tissues as well as in cultured cells. LC3 II is formed by conjugation of LC3 I to phosphatidylethanolamine and is defined as a marker of autophagosomes (Tanida et al., 2008). In the present study, LC3 II was observed to be downregulated in the spinal cord in ATRA-induced SBA rats compared with control rats. We also detected the lysosomal-specific marker LAMP2, as lysosomes play a key degrading role during autophagy (Zhou et al., 2013). LAMP2 expression was significantly decreased in the spinal cords of ATRA-induced SBA rats compared with control rats, which suggested that autophagy was inhibited. In addition to maintaining cellular homeostasis, autophagy eliminates damaged organelles for cellular quality control purposes. Notably, damaged mitochondria can be eliminated through the mitophagy pathway to maintain mitochondrial homeostasis, energy production, and neuronal survival (Wang and Wang, 2019; Wang et al., 2020). In other words, mitophagy dysfunction can result in the accumulation of damaged mitochondria, which in turn is associated with a various neurological pathological conditions (Yang et al., 2018). To investigate whether mitophagy is involved in the pathogenesis of NTDs, we examined the expression of double-labeling for LC3 and TOMM20 by double-labeling immunofluorescence analysis. TOMM20 plays an essential role as a receptor for proteins targeted to mitochondria (Park et al., 2019). In our study, results showed the integration of autophagosomes and damaged mitochondria was inhibited in the SBA group, which suggested that mitophagy was blocked. These results implied that mitophagy was inhibited, which may be the reason of mitochondrial dysfunction in the spinal cords of embryonic rats with SBA.

Sirt1 is important for the turnover of defective mitochondria by mitophagy, and it upregulates autophagy/mitophagy through multiple known pathways (Tang, 2016; Fang, 2019). Mice lacking the gene Sirt1 were small, exhibited notable developmental defects of the retina and heart, and died a few hours after birth (Cheng et al., 2003; Wang et al., 2008; Zhang et al., 2020), suggesting that Sirt1 plays a vital role in mouse embryonic stem cell maintenance and embryonic development (Tang et al., 2017). Indeed, Sirt1 plays protective roles in some neurodegenerative diseases, including motor neuron diseases, and Parkinson’s disease (Herskovits and Guarente, 2014). In our study, Sirt1 expression was significantly decreased in ATRA-induced SBA compared to normal rat embryonic spinal cord tissues. RSV, a natural plant polyphenol, has also been identified as a Sirt1 activator, and a number of studies have linked the beneficial effects of RSV to modulation of Sirt1’s activity (Galiniak et al., 2019). Further, RSV has shown neural protective effects of RSV has been shown in the central nervous system in experimental models (Jardim et al., 2018). In addition, RSV protected against oxidative stress and cell apoptosis in embryos in diabetic embryopathy model rats (Singh et al., 2011). In our rat model, we found RSV attenuated autophagy/mitophagy-associated proteins (Sirt1, PINK1, LC3 II and LAMP2) suppression by ATRA-induced SBA. Moreover, RSV administration reduced the incidence rate of SBA. These data indicated that RSV centrally activates Sirt1, enhancing autophagy/mitophagy and improving the embryonic outcome.

To explore the molecular mechanisms underlying mitophagy in NTDs, further, we compared fetal spinal cord RNA profiles of SBA embryos and normal embryos in the early embryonic stage using RNA-seq. We identified 324 DEGs between the CON and ATRA-induced SBA groups on both E11 and E12. Importantly, GO BP and MF term analysis showed that transcription factors were strongly enriched among the DEGs. Moreover, transcription factor enrichment analysis of the DEGs indicated that the bHLH transcription factor family was significantly enriched. Transcription factors containing the bHLH domain participate in several stages of neurodevelopment, including neurite growth and neuron differentiation (Kabos et al., 2002). Among the candidate DEGs, Bhlhe40 was selected for further study because it was significantly upregulated in the spinal cords of SBA embryos at an early stage and it was predicted to binding to Sirt1 with the highest.

Bhlhe40 is extensively and strongly expressed in nervous tissues (Shi et al., 2010). As a transcriptional repressor, Bhlhe40 binds to target gene promoters to suppress gene expression. Bhlhe40 regulates different targets in different organs, and may have differential roles in various diseases. For example, Bhlhe40 facilitates PI3K/Akt/mTOR activation and triggers tumor progression in breast cancer, whereas it represses STAT1 expression and activates tumor suppression in clear cell carcinoma (Kiss et al., 2020). Downregulation of Bhlhe40 in the prefrontal cortex of spontaneously hypertensive rats improved their hyperactive behavior. Overexpressed Bhlhe40 inhibited the differentiation of PC12 cells in vitro, whereas Bhlhe40 silencing promoted neurite axon and dendrite outgrowth (Wu et al., 2017). Remarkably, we found that Bhlhe40 was significantly upregulated in ATRA-induced SBA rat embryonic spinal cords. This suggested that Bhlhe40 may play a regulatory role in the onset of NTDs. We identified two possible Bhlhe40-response elements (consensus E-box: CACGTG) in the promoter region of Sirt1. As expected, forced expression of Bhlhe40 significantly repressed Sirt1 promoter activity. Importantly, Bhlhe40 overexpression repressed Sirt1 mRNA and protein expression, and Bhlhe40 siRNA significantly enhanced Sirt1 mRNA and protein expression. Taken together, these results suggest that Bhlhe40 negatively regulates Sirt1 transcription by directly binding to its promoter.

As a transcriptional inhibitor, Bhlhe40 is involved various cellular events, including cell proliferation and differentiation, apoptosis, neurogenesis and tumorigenesis (Kabos et al., 2002; Zheng et al., 2009; Kiss et al., 2020). However, the potential role of Bhlhe40 in autophagy/mitophagy had not been elucidated. In our study, overexpressed Bhlhe40 significantly decreased the expression of PINK1, LC3 II and LAMP2, and LC3 puncta colocalized less with TOMM20-stained mitochondria in Bhlhe40-overexpressing cells, indicating mitophagy inhibition. RSV is capable of stimulates autophagy/mitophagy by activating Sirt1, and therefore has shown great potential in the treatment of several diseases. Pre-treatment with RSV reduced the benzo (a) pyrene-induced mitochondrial dysfunction and apoptosis via mitophagy regulation in neurons (Kang et al., 2020). In support of this, our study showed that RSV treatment strongly ameliorated Bhlhe40-induced LC3 II and LAMP2 inhibition in C17.2 cells. Consistent herewith, the Sirt1 inhibitor EX527 (Chen G. et al., 2017), reversed the effects of Bhlhe40 siRNA on LC3 II and LAMP2 expression. PINK1-mediated mitophagy is currently the most well-established way mediate mitophagy in mammalian cells (Youle and Narendra, 2011). Depletion of PINK1, a key mitophagic regulator, impaired mitophagy and blocked the protective role of Sirt1 against compression induced senescence in nucleus pulposus cells (Wang et al., 2020). In our study, we showed that increased level of PINK1 protein after Sirt1 activator RSV treatment in isolated mitochondrial protein. Moreover, suppression of mitophagosomes formation by Bhlhe40 was rescued by RSV in C17.2 cells. Together, these data suggested that the Sirt1/Bhlhe40 axis regulates mitophagy. Further studies on the function of the Bhlhe40/Sirt1 axis in vivo are needed.

In summary, the present study provided evidence of abnormal mitophagy and differential Sirt1 expression in a rat model of ATRA-induced SBA. Moreover, we demonstrated that maternal RSV treatment activated Sirt1 to attenuate autophagy suppression and ameliorate SBA in vivo. Transcriptome analysis showed that transcriptional regulation plays an important role in NTDs. Furthermore, we found that the Bhlhe40/Sirt1 axis regulates mitophagy in neural stem cells. Our findings provide new insights in the pathogenesis of NTDs and may pave the way for potential therapeutic targets for NTDs.

Data Availability Statement

The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author.

Ethics Statement

The animal study was reviewed and approved by the Animal Ethics Committee of Shengjing Hospital of China Medical University.

Author Contributions

ZY and LZ contributed to the study concepts. LZ contributed to the experiment performance. LZ, DL, WM, HG, XW, and WL contributed to the data analysis. LZ wrote the manuscript. ZY made manuscript revisions. All authors contributed to the article and approved the submitted version.

Funding

This work was supported by the National Key R&D Program of China (2016YFC1000505); the National Natural Foundation of China (81871219); the China Postdoctoral Science Foundation (2019M661165); and the LiaoNing Revitalization Talents Program (No. XLYC1902099).

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

We thank all the participants for helpful discussions and suggestions of this study.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fcell.2021.644346/full#supplementary-material

Footnotes

  1. ^ http://jaspar.genereg.net/

References

Bhandari, J., and Thada, P. K. (2020). Neural Tube Disorders. StatPearls [Internet]. Treasure Island, FL: StatPearls Publishing.

Google Scholar

Blencowe, H., Kancherla, V., Moorthie, S., Darlison, M. W., and Modell, B. (2018). Estimates of global and regional prevalence of neural tube defects for 2015: a systematic analysis. Ann. N. Y. Acad. Sci. 1414, 31–46. doi: 10.1111/nyas.13548

PubMed Abstract | CrossRef Full Text | Google Scholar

Brunet, A., Sweeney, L. B., Sturgill, J. F., Chua, K. F., Greer, P. L., Lin, Y., et al. (2004). Stress-dependent regulation of FOXO transcription factors by the SIRT1 deacetylase. Science 303, 2011–2015. doi: 10.1126/science.1094637

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, G., Zhang, B., Xu, H., Sun, Y., Shi, Y., Luo, Y., et al. (2017). Suppression of Sirt1 sensitizes lung cancer cells to WEE1 inhibitor MK-1775-induced DNA damage and apoptosis. Oncogene 36, 6863–6872. doi: 10.1038/onc.2017.297

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, X., Zhong, J., Dong, D., Liu, G., and Yang, P. (2017). Endoplasmic reticulum stress-induced CHOP inhibits PGC-1α and causes mitochondrial dysfunction in diabetic embryopathy. Toxicol. Sci. 158, 275–285. doi: 10.1093/toxsci/kfx096

PubMed Abstract | CrossRef Full Text | Google Scholar

Cheng, H. L., Mostoslavsky, R., Saito, S., Manis, J. P., Gu, Y., Patel, P., et al. (2003). Developmental defects and p53 hyperacetylation in Sir2 homolog (SIRT1)-deficient mice. Proc. Natl. Acad. Sci. U S A. 100, 10794–10799. doi: 10.1073/pnas.1934713100

PubMed Abstract | CrossRef Full Text | Google Scholar

Fang, E. F. (2019). Mitophagy and NAD+ inhibit alzheimer disease. Autophagy 15, 1112–1114. doi: 10.1080/15548627.2019.1596497

PubMed Abstract | CrossRef Full Text | Google Scholar

Fang, E. F., Kassahun, H., Croteau, D. L., Scheibye-Knudsen, M., Marosi, K., Lu, H., et al. (2016). NAD+ replenishment improves lifespan and healthspan in ataxia telangiectasia models via mitophagy and DNA repair. Cell Metab. 24, 566–581. doi: 10.1016/j.cmet.2016.09.004

PubMed Abstract | CrossRef Full Text | Google Scholar

Galiniak, S., Aebisher, D., and Bartusik-Aebisher, D. (2019). Health benefits of resveratrol administration. Acta. Biochim. Pol. 66, 13–21. doi: 10.18388/abp.2018_2749

CrossRef Full Text | Google Scholar

Gelmetti, V., De Rosa, P., Torosantucci, L., Marini, E. S., Romagnoli, A., Di Rienzo, M., et al. (2017). PINK1 and BECN1 relocalize at mitochondria-associated membranes during mitophagy and promote ER-mitochondria tethering and autophagosome formation. Autophagy 13, 654–669. doi: 10.1080/15548627.2016.1277309

PubMed Abstract | CrossRef Full Text | Google Scholar

Giorgi, C., Baldassari, F., Bononi, A., Bonora, M., De Marchi, E., Marchi, S., et al. (2012). Mitochondrial Ca(2+) and apoptosis. Cell Calcium 52, 36–43. doi: 10.1016/j.ceca.2012.02.008

PubMed Abstract | CrossRef Full Text | Google Scholar

Greene, N. D., and Copp, A. J. (2014). Neural tube defects. Annu. Rev. Neurosci. 37, 221–242. doi: 10.1146/annurev-neuro-062012-170354

PubMed Abstract | CrossRef Full Text | Google Scholar

Herskovits, A. Z., and Guarente, L. (2014). SIRT1 in neurodevelopment and brain senescence. Neuron 81, 471–483. doi: 10.1016/j.neuron.2014.01.028

PubMed Abstract | CrossRef Full Text | Google Scholar

Ivankovic, D., Chau, K. Y., Schapira, A. H., and Gegg, M. E. (2016). Mitochondrial and lysosomal biogenesis are activated following PINK1/parkin-mediated mitophagy. J. Neurochem. 136, 388–402. doi: 10.1111/jnc.13412

PubMed Abstract | CrossRef Full Text | Google Scholar

Jardim, F. R., de Rossi, F. T., Nascimento, M. X., da Silva Barros, R. G., Borges, P. A., Prescilio, I. C., et al. (2018). Resveratrol and brain mitochondria: a review. Mol. Neurobiol. 55, 2085–2101. doi: 10.1007/s12035-017-0448-z

PubMed Abstract | CrossRef Full Text | Google Scholar

Juriloff, D. M., and Harris, M. J. (2018). Insights into the etiology of mammalian neural tube closure defects from developmental, genetic and evolutionary studies. J. Dev. Biol. 6:22. doi: 10.3390/jdb6030022

PubMed Abstract | CrossRef Full Text | Google Scholar

Kabos, P., Kabosova, A., and Neuman, T. (2002). Neuronal injury affects expression of helix-loop-helix transcription factors. Neuroreport 13, 2385–2388. doi: 10.1097/00001756-200212200-00002

PubMed Abstract | CrossRef Full Text | Google Scholar

Kang, R. R., Sun, Q., Chen, K. G., Cao, Q. T., Liu, C., Liu, K., et al. (2020). Resveratrol prevents benzo(a)pyrene-induced disruption of mitochondrial homeostasis via the AMPK signaling pathway in primary cultured neurons. Environ. Pollut. 261:114207. doi: 10.1016/j.envpol.2020.114207

PubMed Abstract | CrossRef Full Text | Google Scholar

Khacho, M., Harris, R., and Slack, R. S. (2019). Mitochondria as central regulators of neural stem cell fate and cognitive function. Nat. Rev. Neurosci. 20, 34–48. doi: 10.1038/s41583-018-0091-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Khacho, M., and Slack, R. S. (2018). Mitochondrial dynamics in the regulation of neurogenesis: from development to the adult brain. Dev. Dyn. 247, 47–53. doi: 10.1002/dvdy.24538

PubMed Abstract | CrossRef Full Text | Google Scholar

Kiss, Z., Mudryj, M., and Ghosh, P. M. (2020). Non-circadian aspects of BHLHE40 cellular function in cancer. Genes Cancer 11, 1–19. doi: 10.18632/genesandcancer.201

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, G., Jiapaer, Z., Weng, R., Hui, Y., Jia, W., Xi, J., et al. (2017). Dysregulation of the SIRT1/OCT6 axis contributes to environmental stress-induced neural induction defects. Stem Cell Rep. 8, 1270–1286. doi: 10.1016/j.stemcr.2017.03.017

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, Q., Gao, S., Kang, Z., Zhang, M., Zhao, X., Zhai, Y., et al. (2018). Rapamycin enhances mitophagy and attenuates apoptosis after spinal ischemia-reperfusion injury. Front. Neurosci. 12:865. doi: 10.3389/fnins.2018.00865

PubMed Abstract | CrossRef Full Text | Google Scholar

Lin, C., Li, N., Chang, H., Shen, Y., Li, Z., Wei, W., et al. (2020). Dual effects of thyroid hormone on neurons and neurogenesis in traumatic brain injury. Cell Death Dis. 11:671. doi: 10.1038/s41419-020-02836-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, L., Sakakibara, K., Chen, Q., and Okamoto, K. (2014). Receptor mediated mitophagy in yeast and mammalian systems. Cell Res. 24, 787–795. doi: 10.1038/cr.2014.75

PubMed Abstract | CrossRef Full Text | Google Scholar

Lou, G., Palikaras, K., Lautrup, S., Scheibye-Knudsen, M., Tavernarakis, N., and Fang, E. F. (2020). Mitophagy and neuroprotection. Trends. Mol. Med. 26, 8–20. doi: 10.1016/j.molmed.2019.07.002

PubMed Abstract | CrossRef Full Text | Google Scholar

Love, M. I., Huber, W., and Anders, S. (2014). Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 15:550. doi: 10.1186/s13059-014-0550-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Nikolopoulou, E., Galea, G. L., Rolo, A., Greene, N. D., and Copp, A. J. (2017). Neural tube closure: cellular, molecular and biomechanical mechanisms. Development 144, 552–566. doi: 10.1242/dev.145904

PubMed Abstract | CrossRef Full Text | Google Scholar

Norkett, W., McLone, D. G., and Bowman, R. (2016). Current management strategies of hydrocephalus in the child with open spina bifida. Top Spinal Cord Inj. Rehabil. 22, 241–246. doi: 10.1310/sci2204-241

PubMed Abstract | CrossRef Full Text | Google Scholar

Okada, A., Kushima, K., Aoki, Y., Bialer, M., and Fujiwara, M. (2005). Identification of early-responsive genes correlated to valproic acid-induced neural tube defects in mice. Birth Defects Res. A Clin. Mol. Teratol. 73, 229–238. doi: 10.1002/bdra.20131

PubMed Abstract | CrossRef Full Text | Google Scholar

Park, S. H., Lee, A. R., Choi, K., Joung, S., Yoon, J. B., and Kim, S. (2019). TOMM20 as a potential therapeutic target of colorectal cancer. BMB Rep. 52, 712–717. doi: 10.5483/BMBRep.2019.52.12.249

PubMed Abstract | CrossRef Full Text | Google Scholar

Pickles, S., Vigié, P., and Youle, R. J. (2018). Mitophagy and quality control mechanisms in mitochondrial maintenance. Curr. Biol. 28, R170–R185. doi: 10.1016/j.cub.2018.01.004

PubMed Abstract | CrossRef Full Text | Google Scholar

Radic, J. A. E., Illes, J., and McDonald, P. J. (2019). Fetal repair of open neural tube defects: ethical, legal, and social issues. Camb. Q. Healthc. Ethics 28, 476–487. doi: 10.1017/S0963180119000409

PubMed Abstract | CrossRef Full Text | Google Scholar

Ramya, S., Shyamasundar, S., Bay, B. H., and Dheen, S. T. (2017). Maternal diabetes alters expression of microRNAs that regulate genes critical for neural tube development. Front. Mol. Neurosci. 10:237. doi: 10.3389/fnmol.2017.00237

PubMed Abstract | CrossRef Full Text | Google Scholar

Shi, X., Zheng, Y., Ma, W., and Wang, Y. (2010). Possible involvement of DEC1 on the adverse effects of quinolone antibiotics. Toxicology 271, 1–4. doi: 10.1016/j.tox.2010.03.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Singh, C. K., Kumar, A., Hitchcock, D. B., Fan, D., Goodwin, R., LaVoie, H. A., et al. (2011). Resveratrol prevents embryonic oxidative stress and apoptosis associated with diabetic embryopathy, and improves glucose and lipid profile of diabetic dam. Mol. Nutr. Food Res. 55, 1186–1196. doi: 10.1002/mnfr.201000457

PubMed Abstract | CrossRef Full Text | Google Scholar

Tan, S., and Wong, E. (2017). Mitophagy transcriptome: mechanistic insights into polyphenol-mediated mitophagy. Oxid. Med. Cell. Longev. 2017:9028435. doi: 10.1155/2017/9028435

PubMed Abstract | CrossRef Full Text | Google Scholar

Tang, B. L. (2016). Sirt1 and the mitochondria. Mol. Cells 39, 87–95. doi: 10.14348/molcells.2016.2318

PubMed Abstract | CrossRef Full Text | Google Scholar

Tang, S., Fang, Y., Huang, G., Xu, X., Padilla-Banks, E., Fan, W., et al. (2017). Methionine metabolism is essential for SIRT1-regulated mouse embryonic stem cell maintenance and embryonic development. EMBO J. 36, 3175–3193. doi: 10.15252/embj.201796708

PubMed Abstract | CrossRef Full Text | Google Scholar

Tanida, I., Ueno, T., and Kominami, E. (2008). LC3 and autophagy. Methods Mol. Biol. 445, 77–88. doi: 10.1007/978-1-59745-157-4_4

CrossRef Full Text | Google Scholar

Tsukamoto, S., Kuma, A., Murakami, M., Kishi, C., Yamamoto, A., and Mizushima, N. (2008). Autophagy is essential for preimplantation development of mouse embryos. Science 321, 117–120. doi: 10.1126/science.1154822

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, M., Marco, P., Capra, V., and Kibar, Z. (2019). Update on the role of the non-canonical Wnt/Planar cell polarity pathway in neural tube defects. Cells 8:1198. doi: 10.3390/cells8101198

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, R., and Wang, G. (2019). Autophagy in mitochondrial quality control. Adv. Exp. Med. Biol. 1206, 421–434. doi: 10.1007/978-981-15-0602-4_19

CrossRef Full Text | Google Scholar

Wang, R. H., Sengupta, K., Li, C., Kim, H. S., Cao, L., Xiao, C., et al. (2008). Impaired DNA damage response, genome instability, and tumorigenesis in SIRT1 mutant mice. Cancer Cell 14, 312–323. doi: 10.1016/j.ccr.2008.09.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, Y., Wang, H., Zhuo, Y., Hu, Y., Zhang, Z., Ye, J., et al. (2020). SIRT1 alleviates high-magnitude compression-induced senescence in nucleus pulposus cells via PINK1-dependent mitophagy. Aging (Albany NY) 12, 16126–16141. doi: 10.18632/aging.103587

PubMed Abstract | CrossRef Full Text | Google Scholar

Wu, L. H., Cheng, W., Yu, M., He, B. M., Sun, H., Chen, Q., et al. (2017). Nr3C1-Bhlhb2 axis dysregulation is involved in the development of attention deficit hyperactivity. Mol. Neurobiol. 54, 1196–1212. doi: 10.1007/s12035-015-9679-z

PubMed Abstract | CrossRef Full Text | Google Scholar

Xu, C., Chen, X., Reece, E. A., Lu, W., and Yang, P. (2019). The increased activity of a transcription factor inhibits autophagy in diabetic embryopathy. Am. J. Obstet. Gynecol. 220, 108.e1–108.e12. doi: 10.1016/j.ajog.2018.10.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Xu, C., Li, X., Wang, F., Weng, H., and Yang, P. (2013). Trehalose prevents neural tube defects by correcting maternal diabetes-suppressed autophagy and neurogenesis. Am. J. Physiol. Endocrinol. Metab. 305, E667–E678. doi: 10.1152/ajpendo.00185.2013

PubMed Abstract | CrossRef Full Text | Google Scholar

Xue, J., Gu, H., Liu, D., Ma, W., Wei, X., Zhao, L., et al. (2018). Mitochondrial dysfunction is implicated in retinoic acid-induced spina bifida aperta in rat fetuses. Int. J. Dev. Neurosci. 68, 39–44. doi: 10.1016/j.ijdevneu.2018.04.003

PubMed Abstract | CrossRef Full Text | Google Scholar

Yang, X., Yan, X., Yang, D., Zhou, J., Song, J., and Yang, D. (2018). Rapamycin attenuates mitochondrial injury and renal tubular cell apoptosis in experimental contrast-induced acute kidney injury in rats. Biosci. Rep. 38:BSR20180876. doi: 10.1042/BSR20180876

PubMed Abstract | CrossRef Full Text | Google Scholar

Ye, J., Tong, Y., Lv, J., Peng, R., Chen, S., Kuang, L., et al. (2020). Rare mutations in the autophagy-regulating gene AMBRA1 contribute to human neural tube defects. Hum. Mutat. 41, 1383–1393. doi: 10.1002/humu.24028

PubMed Abstract | CrossRef Full Text | Google Scholar

Youle, R. J., and Narendra, D. P. (2011). Mechanisms of mitophagy. Nat. Rev. Mol. Cell Biol. 12, 9–14. doi: 10.1038/nrm3028

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, W., Feng, Y., Guo, Q., Guo, W., Xu, H., Li, X., et al. (2020). SIRT1 modulates cell cycle progression by regulating CHK2 acetylation-phosphorylation. Cell Death Differ. 27, 482–496. doi: 10.1038/s41418-019-0369-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Zheng, Y., Jia, Y., Wang, Y., Wang, M., Li, B., Shi, X., et al. (2009). The hypoxia-regulated transcription factor DEC1 (Stra13, SHARP-2) and its expression in human tissues and tumours. OMICS 13, 301–306. doi: 10.1089/omi.2009.0014

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhou, J., Tan, S. H., Nicolas, V., Bauvy, C., Yang, N. D., Zhang, J., et al. (2013). Activation of lysosomal function in the course of autophagy via mTORC1 suppression and autophagosome-lysosome fusion. Cell Res. 23, 508–523. doi: 10.1038/cr.2013.11

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhou, X., Zeng, W., Li, H., Chen, H., Wei, G., Yang, X., et al. (2018). Rare mutations in apoptosis related genes APAF1, CASP9, and CASP3 contribute to human neural tube defects. Cell Death Dis. 9:43. doi: 10.1038/s41419-017-0096-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: spina bifida aperta, mitophagy, Bhlhe40, Sirt1, autophagy, resveratrol

Citation: Zhao L, Liu D, Ma W, Gu H, Wei X, Luo W and Yuan Z (2021) Bhlhe40/Sirt1 Axis-Regulated Mitophagy Is Implicated in All-Trans Retinoic Acid-Induced Spina Bifida Aperta. Front. Cell Dev. Biol. 9:644346. doi: 10.3389/fcell.2021.644346

Received: 21 December 2020; Accepted: 06 April 2021;
Published: 27 April 2021.

Edited by:

Roland Wohlgemuth, Lodz University of Technology, Poland

Reviewed by:

Richard H. Finnell, Baylor College of Medicine, United States
Silvia Pomella, Bambino Gesù Children Hospital (IRCCS), Italy

Copyright © 2021 Zhao, Liu, Ma, Gu, Wei, Luo and Yuan. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Zhengwei Yuan, eXVhbnp3QGhvdG1haWwuY29t

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.