Maternal diabetes-mediated RORA suppression contributes to gastrointestinal symptoms in autism-like mouse offspring

Background Retinoic acid-related orphan receptor alpha (RORA) has been reported to be suppressed in autistic patients and is associated with autism spectrum disorders (ASD), although the potential role and mechanism of RORA on gastrointestinal (GI) symptoms in ASD patients is still not reported. In this study, we aim to investigate the contribution of RORA to GI symptoms through a maternal diabetes-mediated autism-like mouse model. Results Male offspring of diabetic dams were treated with either superoxide dismutase (SOD) mimetic MnTBAP or RORA agonist SR1078, or were crossbred with intestine epithelial cells (IEC)-specific RORA knockout (RORA−/−) mouse. Gene expression, oxidative stress and inflammation were measured in brain tissues, peripheral blood mononuclear cells (PBMC) and IEC, and GI symptoms were evaluated. Our results showed that SOD mimetic MnTBAP completely, while RORA agonist SR1078 partly, reversed maternal diabetes-mediated oxidative stress and inflammation in the brain, PBMC and IEC, as well as GI symptoms, including intestine permeability and altered gut microbiota compositions. IEC-specific RORA deficiency either mimicked or worsened maternal diabetes-mediated GI symptoms as well as oxidative stress and inflammation in IEC, while there was little effect on maternal diabetes-mediated autism-like behaviors. Conclusions We conclude that RORA suppression contributes to maternal diabetes-mediated GI symptoms in autism-like mouse offspring, this study provides a potential therapeutical target for maternal diabetes-mediated GI symptoms in offspring through RORA activation. Supplementary Information The online version contains supplementary material available at 10.1186/s12868-022-00693-0.


Background
Autism spectrum disorders (ASD) are characterized as a group of neurodevelopment disorders with impaired social interaction and communication as well as restricted behaviors that may be triggered by multistage prenatal risk factor exposure [1][2][3]. Over the past few decades, the prevalence of ASD has increased to around 1:59, with a male/female ratio of 4/1 [2,3]. A variety of risk factors, including epigenetics and environmental exposures, have been reported to be associated with ASD [4,5]. We have previously found that prenatal exposure to either maternal diabetes [6][7][8] or hormones (e.g. progestin and androgen) [9][10][11][12] contributes to ASD development through oxidative stress, epigenetic changes and gene suppression [13,14], while detailed mechanism remains largely unknown.
Accumulating evidence shows that gastrointestinal (GI) symptoms are concurrent and potentially associated with ASD [15]. Gut microbiota may participate in ASD pathogenesis [16], and clinical study shows that microbiota transfer therapy may significantly improve GI and ASD symptoms in autistic patients [17]. Gut microbiota may contribute to brain dysfunction and play a role in ASD development through the gut-brain axis, while the related mechanisms and cause-effect relationship between ASD and GI problems are still unclear [15,[17][18][19][20].
Retinoic acid-related orphan receptor alpha (RORA) is involved in many pathophysiological processes [21,22] and has been shown to be associated with ASD development, with suppressed expression of RORA being identified in autistic patients [23]. Aromatase, encoded by gene cytochrome P450, family 19 (CYP19A1), is the rate-limiting enzyme that is responsible for the conversion of androgen into estradiol (E2) and has been reported to be suppressed in ASD patients and is regulated by RORA either directly or indirectly [24][25][26][27][28]. It has been also reported that RORA regulates the expression of superoxide dismutase 2 (SOD2) through retinoic acidrelated orphan receptor elements (RORE) located on the SOD2 promoter [29]. Furthermore, we have early found that maternal diabetes-mediated autism-like mouse offspring has reduced RORA expression in some brain regions as well as in intestinal epithelial cells (IEC). We then hypothesize that RORA may play a role in GI symptoms in maternal diabetes-mediated autism-like mouse offspring through pathways involving the suppression of RORA and its target genes CYP19A1 and SOD2 in the gut system [30,31].
We have previously found that maternal diabetes-mediated mouse offspring show autism-like behaviors with several core features of ASD, which including decreased ultrasonic vocalizations as well as suppressed social recognition and interaction [8,32]. In this study, we aim to investigate the potential role of RORA on GI symptoms in maternal diabetes-mediated autism-like mouse models. Our results showed that maternal diabetes-mediated offspring showed significantly suppressed expression of RORA and its target genes CYP19A1 and SOD2 in the brain, peripheral blood mononuclear cells (PBMC) and IEC. The offspring also showed GI symptoms, including oxidative stress, pro-inflammatory cytokine release, increased intestinal permeability and altered gut microbiota composition. Treatment of SOD mimetic MnTBAP [8] completely, while RORA agonist SR1078 [33] partly, reversed this effect. Further experiments using the RORA knockout (RORA −/− ) mouse model showed that intestine-specific RORA deficiency mimicked or worsened maternal diabetes-mediated oxidative stress and inflammation in IEC as well as GI symptoms. We conclude that RORA suppression contributes to maternal diabetesmediated GI symptoms in mouse offspring.

Results
RORA agonist and SOD mimetic reverse maternal diabetes-mediated RORA suppression and oxidative stress in the brain Dams were treated by control (CTL/VEH), diabetes (STZ/VEH), diabetes plus SOD mimetic MnT-BAP (STZ/MnTBAP), or diabetes plus RORA agonist SR1078 (STZ/SR1078) during pregnancy. The amygdala neurons and brain tissues, including the amygdala, hypothalamus and hippocampus, were isolated from subsequent male offspring for biomedical analysis. We first evaluated the gene expression of RORA and its target genes from amygdala tissues. The results showed that maternal diabetes (STZ/VEH) treatment significantly deceased the mRNA levels of RORA, CYP19A1 and SOD2 (P < 0.0001), respectively, compared to the CTL/ VEH group; SOD mimetic MnTBAP treatment (STZ/ MnTBAP) completely reversed this effect; while RORA agonist (STZ/SR1078) completely reversed the expression of CYP19A1, partly reversed on SOD2 (P < 0.01), but showed no effect on RORA (see Fig. 1a). We also measured the protein levels by western blotting, and an expression pattern similar with that of mRNA levels was observed (see Fig. 1b, c, Additional file 1: Fig. S1a). We then evaluated oxidative stress, and the results showed that maternal diabetes (STZ/VEH) significantly increased ROS formation (see Fig. 1d) and 3-nitrotyrosin formation (see Fig. 1e) (P < 0.0001), respectively, compared to CTL/VEH group. STZ/MnTBAP treatment completely, while RORA agonist (STZ/SR1078) partly, reversed this effect. We also evaluated SOD2 activity, and the results showed that maternal diabetes (STZ/VEH) significantly decreased SOD2 activity (P < 0.0001) compared to the CTL/VEH group; treatments of either STZ/MnTBAP or STZ/SR1078 completely reversed this effect (see Fig. 1f ). Finally, we evaluated 8-oxo-dG formation from amygdala neurons. The results showed that maternal diabetes (STZ/VEH) significantly increased 8-oxo-dG formation (P < 0.0001) compared to the CTL/VEH group; treatments of either STZ/MnTBAP or STZ/SR1078 had no significant effect (see Fig. 1g, h). In addition, we evaluated gene expression in the hypothalamus and hippocampus, and the results showed that maternal diabetes (STZ/ VEH) significantly decreased expression of RORA and CYP19A1 in both the hypothalamus and hippocampus compared to the CTL/VEH group. STZ/MnTBAP treatment completely reversed this effect on both genes, while STZ/MnTBAP treatment completely reversed this effect on CYP19A1 expression but had no effect on RORA. Additionally, none of the treatments had any effect on SOD2 expression in tissues of the hypothalamus and hippocampus (see Additional file 1: Fig. S2). We conclude that RORA agonist and SOD mimetic reverse maternal diabetes-mediated RORA suppression and oxidative stress in the brain.

RORA agonist and SOD mimetic ameliorate maternal diabetes-mediated autism-like behaviors in male offspring
We evaluated the effects of RORA agonist and SOD mimetic MnTBAP on maternal diabetes-mediated autism-like behaviors. Our result showed that offspring from the maternal diabetes (STZ/VEH) treatment group displayed significantly decreased ultrasonic vocalizations (P < 0.0001) compared to CTL/VEH treatment; both SOD mimetic (STZ/MnTBAP) or RORA agonist (STZ/ SR1078) partly reversed this effect (P < 0.01, see Fig. 2a). Additionally, mice in the maternal diabetes (STZ/VEH) group spent significantly less time sniffing (P < 0.0001), mounting (P < 0.01), and interacting in total (P < 0.0001) during the social interaction tests compared to the CTL/ VEH group; both STZ/MnTBAP and STZ/SR1078 treatment partly reversed the effect on sniffing (P < 0.01) and total interaction time (P < 0.01) while completely reversing the effect on mounting interaction time. None of the treatments appeared to have a significant effect on the amount of time spent grooming partners (see Fig. 2b). Furthermore, mice from the maternal diabetic (STZ/ VEH) treatment group spent significantly less time in The amygdala neurons were isolated on embryonic day (E18) from the above treatment for immunostaining. g Quantitation of 8-oxox-dG staining, n = 5. h Representative pictures for 8-oxo-dG staining (green) and DAPI staining for nuclei (blue). The one-way ANOVA analysis was performed to determine statistical significance of different groups. *P < 0.05, vs. CTL/VEH group; ¶ P < 0.05, vs. STZ/VEH group. Data were expressed as mean ± SD the stranger 1 side (P < 0.01) and more time in the empty side (P < 0.01) during the sociability test (see Fig. 2c); while spent significant more time in the stranger 1 side (P < 0.01) and less time in the stranger 2 side (P < 0.01) in the social novelty test compared to the CTL/VEH treatment group (see Fig. 2d); both STZ/MnTBAP and STZ/ SR1078 treatment completely reversed the effect. We conclude that RORA agonist and SOD mimetic ameliorate maternal diabetes-mediated autism-like behaviors in male offspring.

RORA agonist and SOD mimetic reverse maternal diabetes-mediated pre-inflammatory cytokine release in PBMC
We evaluated the potential effect of RORA agonist on pre-inflammatory cytokine release in PBMC cells. Our results showed that maternal diabetes (STZ/VEH) treatment significantly increased mRNA levels of IL1β, IL6, MCP1 and IL17A (P < 0.0001), respectively, compared to the CTL/VEH treatment; SOD mimetic (STZ/MnTBAP) treatment completely reversed the effects on the expression of IL6, MCP1 and IL17A, while partly reversed the effects on IL1β expression (P < 0.001). RORA agonist (STZ/ SR1078) completely reversed the effects on the expression of IL1β and IL6 and partly reversed the effects on MCP1 (P < 0.001) and IL17A (P < 0.01) expression (see Fig. 3a). We also measured the cytokine secretion in cultured PBMC, which included IL1β (see Fig. 3b), IL6 (see Fig. 3c), MCP1 (see Fig. 3d) and IL17A (see Fig. 3e). The results showed that maternal diabetes (STZ/VEH) treatment significantly increased pro-inflammatory cytokine secretion compared to the CTL/VEH treatment (P < 0.0001), and treatments of STZ/MnTBAP and STZ/SR1078 either completely or partly reversed the effects of maternal diabetes. We further evaluated the gene expression of RORA and its target genes in PBMC, and the results showed that maternal diabetes (STZ/VEH) treatment significantly decreased mRNA levels of RORA and its target genes CYP19A1 and SOD2 Dams were divided into the following treatment groups: control (CTL/VEH), diabetes (STZ/VEH), diabetes plus SOD mimetic MnTBAP (STZ/MnTBAP), or diabetes plus RORA agonist SR1078 (STZ/SR1078) during pregnancy, and the subsequent male offspring were used for evaluation of autism-like behaviors. a Ultrasonic vocalization, n = 9. b Social interaction (SI) test, with total interaction time and amount of time spent following, mounting, grooming, and sniffing any body parts of the other mouse were calculated, n = 9. c, d Three-chambered social tests, n = 9. c Time spent in chamber for sociability. d Time spent in chamber for social novelty. The one-way ANOVA analysis was performed to determine statistical significance of different groups. *P < 0.05, vs. CTL/VEH group; ¶ P < 0.05, vs. STZ/VEH group. Data were expressed as mean ± SD (P < 0.0001) compared to CTL/VEH treatment, and STZ/ MnTBAP treatment completely, while STZ/SR1078 treatment partly, reversed this effect (see Fig. 3f). In addition, we evaluated the potential effect of maternal diabetes on epigenetic changes on the RORA promoter in PBMC cells, and the results showed that the different treatments had no effect on DNA methylation (see Additional file 1: Fig.  S3), Histone 4 (H4) methylation (see Additional file 1: Fig.  S4a), or histone acetylation (see Additional file 1: Fig. S4b). Finally, we measured histone 3 methylation on the RORA promoter, and the results showed that the maternal diabetes (STZ/VEH) treatment group displayed significantly increased H3K9me3 modification (P < 0.0001) compared to the CTL/VEH treatment. STZ/MnTBAP treatment completely, while STZ/SR1078 treatment partly reversed this effect (P < 0.01). The treatments had no effect on H3K9me2, H3K27me2 or H3K27me3 modifications (see Fig. 3g). We conclude that RORA agonist and SOD mimetic reverse maternal diabetes-mediated pre-inflammatory cytokine release in PBMC cells.

RORA agonist and SOD mimetic reverse maternal diabetes-mediated RORA suppression and oxidative stress in IEC
We evaluated the potential effect of maternal diabetes and RORA agonist on gene expression and oxidative stress in IEC cells. The results showed that maternal diabetes (STZ/VEH) significantly decreased the mRNA levels of RORA, CYP19A1 and SOD2 (P < 0.0001) compared to the CTL/VEH group; MnTBAP treatment (STZ/MnTBAP) completely, while STZ/SR1078 treatment partly (P < 0.01), reversed this effect (see Fig. 4a). We also measured the protein levels by western blotting, and an expression pattern similar with that of mRNA levels was observed (see Fig. 4b, c, Additional file 1: Fig. S1b). We then evaluated oxidative stress, and the results showed that maternal diabetes (STZ/VEH) significantly increased ROS formation (see Fig. 4d) and 3-nitrotyrosin formation (see Fig. 4e), respectively, compared to the CTL/VEH group (P < 0.0001); STZ/ MnTBAP treatment completely, while RORA agonist (STZ/SR1078) partly, reversed this effect. We finally evaluated the SOD2 activity, and the results showed that offspring in the maternal diabetes (STZ/VEH) group showed decreased SOD2 activity (P < 0.0001) compared to the CTL/VEH group; again, STZ/MnT-BAP treatment completely, while RORA agonist (STZ/ SR1078) partly (P < 0.01), reversed this effect (see Fig. 4f ). We conclude that RORA agonist and SOD mimetic reverse maternal diabetes-mediated RORA suppression and oxidative stress in IEC. RORA agonist and SOD mimetic reverse maternal diabetes-mediated pre-inflammatory cytokine release in PBMC. Dams were divided into the following treatment groups: control (CTL/VEH), diabetes (STZ/VEH), diabetes plus SOD mimetic MnTBAP (STZ/MnTBAP), or diabetes plus RORA agonist SR1078 (STZ/SR1078) during pregnancy, and the PBMC were isolated from male offspring for biomedical analysis. a mRNA levels of pre-inflammatory cytokines, n = 4. b IL-1β secretion, n = 9. c IL-6 secretion, n = 9. d MCP1 secretion, n = 9. e IL17A secretion, n = 9. f mRNA levels of RORA and target genes, n = 4. g ChIP analysis on the RORA promoter, n = 4. The one-way ANOVA analysis was performed to determine statistical significance of different groups. *P < 0.05, vs. CTL/VEH group; ¶ P < 0.05, vs. STZ/VEH group. Data were expressed as mean ± SD

RORA agonist and SOD mimetic reverse maternal diabetes-mediated gastrointestinal dysfunction
We evaluated the potential effect of maternal diabetes and RORA agonist on gastrointestinal function. First, we measured cytokine secretion from isolated IEC cells, and the results showed that maternal diabetes (STZ/VEH) significantly increased cytokine secretion of IL1β (see Fig. 5a), IL6 (see Fig. 5b), MCP1 (see Fig. 5c) and IL17A (see Fig. 5d) compared to the CTL/ VEH group (P < 0.0001); treatments of STZ/MnTBAP and STZ/SR1078 partly reversed this effect. We then evaluated GI function by intestinal permeability assay using FITC-dextran, and the results showed that maternal diabetes (STZ/VEH) significantly increased intestinal permeability (P < 0.0001) compared to the CTL/VEH group; STZ/MnTBAP treatment completely, while STZ/ SR1078 treatment partly (P < 0.01), reversed this effect (see Fig. 5e). We then evaluated the microbial populations of the treated mice using 16S rRNA gene sequencing, and the results showed that maternal diabetes (STZ/ VEH) resulted in significant differences in gut microbial composition compared to CTL/VEH group. In general, the phyla Firmicutes and Proteobacteria dominate the microbiome of CTL/VEH mice, while a shift towards Firmicutes and Actinobacteria occurred in STZ/VEH mice; and treatments of both STZ/MnTBAP and STZ/ SR1078 partly reversed this effect (see Fig. 5f ). Additionally, there was no difference in microbial species richness (see Fig. 5g) and diversity (see Fig. 5h) across the treatment groups. Furthermore, we evaluated the relative abundance of different gut organisms. The genus level of Mucispirillum (g_Mucispirillum) was significantly decreased in the STZ/VEH treatment group (P < 0.0001) compared to the CTL/VEH group; STZ/MnTBAP treatment completely, while STZ/SR1078 treatment partly (P < 0.01), reversed this effect (see Fig. 5i). In addition, a genus belonging to the phylum Deferribacteres (p_Deferribacteres) was significantly decreased in STZ/VEH treatment compared to the CTL/VEH group, and treatments of STZ/MnTBAP and STZ/SR1078 completely reversed this effect. On the other hand, genus belonging to the phylum Proteobacteria (p_Proteobacteria) (P < 0.001) and Tenericutes (p_Tenericutes) (P < 0.001) were significantly increased in the STZ/VEH treatment group compared to the CTL/VEH group; treatments of STZ/MnTBAP and STZ/SR1078 partly reversed this effect (see Fig. 5j). In general, our results indicate that RORA agonist and SOD mimetic reverse maternal diabetes-mediated

Intestine-specific RORA deficiency mimics or worsens maternal diabetes-mediated RORA suppression and oxidative stress in IEC
Dams in either the control (CTL) or diabetes (STZ) groups were crossbred with either wild type (WT) or intestinal epithelial-specific RORA knockout (RORA −/− ) mice, and the IEC cells were isolated from subsequent male offspring for biomedical analysis. We first evaluated gene expression from IEC, and the results showed that maternal diabetes (STZ/WT) significantly decreased mRNA levels of RORA, CYP19A1 and SOD2, respectively, compared to the CTL/WT group. RORA deficiency (RORA −/− ), including treatments of CTL/RORA −/− and STZ/RORA −/− , either mimicked or worsened this effect (see Fig. 6a). We also measured the protein levels by western blotting, and an expression pattern similar to that of mRNA levels was observed (see Fig. 6b, c, Additional file 1: Fig.  S1c). We then evaluated oxidative stress, and the results showed that maternal diabetes (STZ/WT) significantly increased ROS formation (see Fig. 6d) and 3-nitrotyrosin formation (see Fig. 6e), respectively, compared to the CTL/WT group; RORA deficiency (RORA −/− ) either mimicked or worsened this effect. We then evaluated SOD2 activity, and the results showed that maternal diabetes (STZ/WT) significantly decreased SOD2 activity compared to the CTL/WT group (see Fig. 6f ); we then evaluated 8-oxo-dG formation by immunostaining in IEC, and the results showed that maternal diabetes (STZ/WT) significantly increased 8-oxo-dG formation compared to the CTL/WT group (see Fig. 6g, h); again, RORA deficiency (RORA −/− ) either mimicked or worsened these effects. We conclude that intestine-specific RORA deficiency mimics secretion; c MCP1 secretion; d IL17A secretion. e Intestinal permeability assay by FITC-dextran, n = 5. f-j Gut microbiota analysis, n = 9; f overview of the identified relative frequencies of different phyla found in treated mice; g species richness; h species diversity; i relative abundance of Mucispirillumn at genus level; j relative abundance of different bacteria in phylum level. The one-way ANOVA analysis was performed to determine statistical significance of different groups. *P < 0.05, vs. CTL/VEH group; ¶ P < 0.05, vs. STZ/VEH group. Data were expressed as mean ± SD or worsens maternal diabetes-mediated RORA suppression and oxidative stress in IEC.

Intestine-specific RORA deficiency does not affect maternal diabetes-mediated autism-like behaviors
We evaluated the potential effect of intestine-specific RORA deficiency on maternal diabetes-mediated gene expression and autism-like behaviors. Our results showed that maternal diabetes (STZ/WT) treatment significantly decreased the mRNA levels of RORA, CYP19A1 and SOD2 in the amygdala compared to the CTL/WT group, and intestine-specific RORA deficiency (RORA −/− ) showed no effect on gene expression (see Additional file 1: Fig. S5a). We then evaluated autism-like behaviors in those animals, and the results showed that maternal diabetes (STZ/WT) treatment significantly decreased ultrasonic vocalizations (see Additional file 1: Fig. S5b).
Moreover, offspring from the STZ/WT group spent significantly less time interacting with stranger mice during the social interaction tests (see Additional file 1: Fig. S5c) and showed significant less interest in sociability (see Additional file 1: Fig. S5d) and social novelty (see Additional file 1: Fig. S5e) compared to the CTL/WT group. RORA deficiency (RORA −/− ), including treatments of CTL/RORA −/− and STZ/RORA −/− , showed no effect on these behaviors. We conclude that intestine-specific RORA deficiency does not affect maternal diabetesmediated autism-like behaviors.

Intestine-specific RORA deficiency mimics or worsens maternal diabetes-mediated gastrointestinal dysfunction
We evaluated the potential effect of RORA deficiency on maternal diabetes-mediated gastrointestinal dysfunction. We first measured cytokine secretion from isolated Intestine-specific RORA deficiency mimics or worsens maternal diabetes-mediated RORA suppression and oxidative stress in IEC. Dams in either the control (CTL) or diabetes (STZ) groups were crossbred with either wild type (WT) or intestinal epithelial-specific RORA knockout (RORA −/− ) mice, and the IEC cells were isolated from subsequent male offspring for biomedical analysis. a mRNA levels by qPCR, n = 4. b Quantitation of protein levels, n = 5. c Representative pictures of western blotting for b, and full-length blots are presented in Additional file 1: Fig. S1c. d ROS formation, n = 5. e 3-Nitrotyrosine formation, n = 5. f SOD2 activity, n = 5. g, h The IEC were isolated for 8-oxo-dG immunostaining; g representative pictures for 8-oxo-dG staining (green) and DAPI staining for nuclei (blue); h quantitation of 8-oxox-dG staining, n = 5; The two-way ANOVA analysis was performed to determine statistical significance of different groups. ***P < 0.0001; **P < 0.001; *P < 0.01; N/S, no significance. Data were expressed as mean ± SD IEC cells, and the results showed that maternal diabetes (STZ/WT) significantly increased cytokine secretion of IL1β (see Fig. 7a), IL6 (see Fig. 7b), MCP1 (see Fig. 7c) and IL17A (see Fig. 7d) compared to CTL/WT group; RORA deficiency (RORA −/− ), including treatments of CTL/RORA −/− and STZ/RORA −/− , either mimicked or worsened this effect. We then evaluated GI function by intestinal permeability assay, and the results showed that maternal diabetes (STZ/WT) significantly increased intestinal permeability compared to CTL/WT group; RORA deficiency (RORA −/− ) either mimicked or worsened this effect (see Fig. 7e). We then evaluated the microbial populations of the treated mice using 16S rRNA gene sequencing, and the results showed that maternal diabetes (STZ/WT) treatment achieved significant differences in gut microbial composition compared to the CTL/WT group; in general, the phyla Firmicutes and Proteobacteria dominate the microbiome of CTL/WT mice, while a shift towards Firmicutes and Actinobacteria occurred in STZ/WT mice; RORA deficiency (RORA −/− ) mimicked the effect of STZ/WT treatment (see Fig. 7f ). Again, there was no difference in the microbial species richness (see Fig. 7g) and diversity (see Fig. 7h) among those treatments. Furthermore, we evaluated the relative abundance of different gut organisms. The genus level of Mucispirillum (g_Mucispirillum) was significantly decreased in the STZ/WT treatment group compared to the CTL/WT group; RORA deficiency (RORA −/− ) either mimicked or worsened this effect (see Fig. 7i). In addition, there were significantly fewer organisms of a genus belonging to the phylum Deferribacteres (p_Deferribacteres) in STZ/WT treatment compared to the CTL/WT group; RORA deficiency (RORA −/− ) either mimicked or worsened this effect. On the other hand, the number organisms of a genus belonging to the phylum Proteobacteria (p_Proteobacteria) and Tenericutes (p_Tenericutes) was significantly increased in the STZ/ WT treatment group compared to CTL/WT group; Fig. 7 Intestine-specific RORA deficiency mimics or worsens maternal diabetes-mediated gastrointestinal dysfunction. Dams from either the control (CTL) or diabetes (STZ) groups were crossbred with either wild type (WT) or intestinal epithelial-specific RORA knockout (RORA −/− ), and the gastrointestinal function of male offspring was evaluated. a-d The IEC cells were isolated to evaluate cytokine secretion, n = 9; a IL-1β secretion; b IL-6 secretion; c MCP1 secretion; d IL17A secretion. e Intestinal permeability assay by FITC-dextran, n = 5. f-j Gut microbiota analysis, n = 9; f overview of the identified relative frequencies of different phyla found in treated mice; g species richness; h species diversity; i relative abundance of Mucispirillumn at genus level; j relative abundance of different bacteria in phylum level. The two-way ANOVA analysis was performed to determine statistical significance of different groups. ***P < 0.0001; **P < 0.001; *P < 0.01; N/S, no significance. Data were expressed as mean ± SD and RORA deficiency (RORA −/− ) either mimicked or worsened this effect (see Fig. 7j). In general, our results indicate that intestine-specific RORA deficiency mimics or worsens maternal diabetes-mediated gastrointestinal dysfunction, including cytokine release, increased intestinal permeability and changes of the gut microbial composition.

Schematic model for maternal diabetes-mediated GI symptoms through RORA suppression in autism-like offspring
Maternal diabetes-mediated oxidative stress [34] triggers suppression of RORA and its target genes (CYP19A1 and SOD2) in the brain, leading to ASD development. Simultaneously, maternal diabetes-induces RORA suppression by epigenetic changes on the RORA promoter in PBMC, resulting in pro-inflammatory cytokine release, which including IL1β, IL6, MCP1 and IL17A. Furthermore, maternal diabetes-mediated RORA suppression in IEC leads to inflammation, increased intestinal permeability and altered gut microbial composition, subsequently resulting in GI symptoms. On the other hand, the cause-effect relationship between ASD and GI symptoms in maternal diabetes-mediated mouse model remains unclear (see Fig. 8).

Discussion
In this study, we found that maternal diabetes-mediated autism-like mouse offspring have suppressed expression of RORA and its target gene in addition to oxidative stress and inflammation in brain tissues, PBMC and IEC, as well GI symptoms. RORA agonist SR1078 partly reversed this effect, and intestine-specific RORA deficiency mimics or worsens maternal diabetes-mediated GI symptoms.

RORA suppression in autism-like behaviors
It has previously been reported that RORA expression is reduced in the brains of ASD patients [23], while its potential role in ASD development is still largely unknown. Our results showed that RORA expression was significantly suppressed in brain tissues of maternal diabetes-mediated autism-like mouse models. SOD mimetic MnTBAP completely, while RORA agonist SR1078 partly, reversed maternal diabetes-mediated autism-like behaviors, indicating that maternal diabetes-mediated oxidative stress plays a major role [8]. This effect is at least partly mediated through RORA and its target genes, including CYP19A1 and SOD2, since CYP19A1 converts androgen into estradiol, which directly or indirectly regulates many antioxidant genes [35,36]. RORA may play an important role in sex bias in ASD development [37] by modulating the expression of CYP19A1 and SOD2 as well as oxidative stress in neurons [8,28].

RORA suppression in PBMC and its contribution to GI symptoms
Our results showed that maternal diabetes-mediated autism-like offspring had significant GI symptoms, together with suppressed RORA expression in PBMC [38] and increased secretion of proinflammatory cytokines, including IL1β, IL6 [39], MCP1 and IL17A [40]. It has been reported that altered gut microbiota composition in ASD is associated with systemic immune dysfunction, such as chronic inflammation [38], and RORA regulates the inflammatory signaling pathway by modulation of macrophage development and activation [41,42]. We suppose that maternal diabetes-mediated RORA suppression in PBMC may contribute to GI symptoms in autism-like offspring by modulation of systemic inflammation in PBMC cells.

RORA suppression in IEC and its contribution to GI symptoms
We found that maternal diabetes-mediated autism-like offspring had RORA suppression in IEC with significant GI symptoms [43][44][45], including oxidative stress and inflammation in IEC [46], increased intestine permeability and altered microbiota compositions. SOD mimetic MnTBAP completely, while RORA agonist SR1078 partly, reversed this effect. This indicates that RORA suppression may play an important role in maternal diabetes-mediated GI symptoms, although this does not exclude other potential factors. Furthermore, IEC-specific RORA deficiency mimicked or worsened maternal diabetes-mediated GI symptoms, while surprisingly, it showed little effect on maternal diabetes-mediated autism-like behaviors. It has been widely reported that gut microbiota may modulate animal behaviors through the gut-brain axis and fecal microbiota transplantation (FMT) [47] rescues social behaviors in ASD model [48,49]. Our results indicate that IEC may not play a major role in gut microbiotamediated animal behaviors, while it has been reported that gut microbiota may modulate autism-like behaviors through the vagus nerve or enteric nervous system [50,51], and this may be the potential explanation for our findings.

Limitations
We have previously showed that the amygdala plays a major role in maternal diabetes-mediated autism-like behavior [6][7][8]; therefore, gene expression in the amygdala was a key focus in this study, while this does not exclude the potential contribution and effect from other brain area on animal behaviors. Furthermore, only male offspring were studied here in order to avoid potential interference from estrogen/estrogen receptor-mediated effects in female offspring [6][7][8], and this may potentially affect the conclusions about the role of the estrous cycle in the assessment of animal behaviors and GI symptoms. Additionally, the lack of IEC-specific RORA expression mouse is a potential limitation to investigate the mechanism of RORA and its contribution to GI symptoms.

Conclusions
Maternal diabetes-mediated autism-like mouse offspring display RORA suppression in the brain, PBMC and IEC with oxidative stress and inflammation. They also show a variety of GI symptoms, including increased intestine permeability and altered gut microbiota composition. Treatment with either SOD mimetic or RORA agonist partly reverses this effect. Furthermore, IEC-specific RORA deficiency either mimics or worsens maternal diabetes-mediated GI symptoms. We conclude that RORA suppression contributes to maternal diabetes-mediated GI symptoms in autism-like mouse offspring, this study provides a potential therapeutical target for maternal diabetes-mediated GI symptoms in offspring through RORA activation.

Methods
A detailed description can be found in Additional file 1: Data S1, and the related primers used in this study are shown in Additional file 1: Table S1.

In vivo mouse experiments
The animal protocol conformed to US NIH guidelines (Guide for the Care and Use of Laboratory Animals, No. 85-23, revised 1996), and was reviewed and approved by the Institutional Animal Care and Use Committee from Foshan Maternity and Child Health Care Hospital and Hainan Women and Children's Medical Center, and this study is in accordance with the ARRIVE guidelines. In this study, 10 dams were assigned for each exposure, and one representative offspring was picked up randomly from each dam for experiments and analysis, and 9 representative offspring were usually picked up from 10 in total in case animal may die by accident during experimental process. Each animal was firstly used for behavior analysis, then was sacrificed for tissue collection and subsequent biological analysis. Each experimental offspring was in separate cage since they came from different dams.

Generation of intestine-specific RORA knockout mice
The RORA fl/fl mouse, which has loxP flanking sites targeting exon 3 of the RORA gene, was generated by in vitro fertilization and was obtained for the study as a generous gift from Dr. Haimou Zhang from Hubei University. The Villin-cre (Vil1-cre) transgenic mice (#021504) have the mouse villin 1 promoter directing expression of Cre recombinase to villus and crypt epithelial cells of the small and large intestines, was obtained from Jackson Laboratories. To generate intestine-specific RORA −/− null mouse (Vil1-cre-RORA fl/fl ), the RORA fl/fl mice were cross-bred with Vil1-cre mice for over 4 generations on the C57BL/6J background. Positive offspring were confirmed by genotyping through PCR using specific primers (see Additional file 1: Table S1) for the presence of both loxP sites within RORA alleles and Cre recombinase [22].

Generation of diabetic mice
All the experimental mice were either RORA wild type (WT) or RORA null (RORA −/− ) mice with a C57BL/6J mixed genetic background. Chronic diabetic female mice were induced by injection of 35 mg/kg streptozotocin (STZ), and the diabetic mice were confirmed with blood glucose > 250 mg/dl.

Mouse protocol 1: prenatal treatment of diabetes or chemicals
Verified pregnant dams (n = 9 for each group) were randomly assigned to the following 4 groups: Group 1: CTL group mice receiving only vehicle treatment (CTL/VEH); Group 2: STZ mice receiving only vehicle treatment (STZ/VEH); Group 3: CTL group mice receiving 10 mg/ kg/day of MnTBAP (dissolved in DMSO) injection (CTL/ MnTBAP); Group 4. STZ mice receiving 10 mg/kg/day of RORA agonist SR1078 (dissolved in DMSO) injection (STZ/SR1078). The injection was conducted on days 1, 4, 7, 10, 13, 16 and 19 of pregnancy, respectively. Neurons from the amygdala were isolated on embryonic day 18 (E18) as described below. The male offspring (1 male offspring from each dam was randomly selected for experiments) were separated from the dams on day 21, fed with normal chow, and underwent behavior tests until 7-8 weeks of age. GI symptoms were evaluated, the whole blood was collected by heart puncture and the peripheral blood mononuclear cells (PBMC) were isolated, and the intestine epithelial cells (IEC) were isolated as described below for further biomedical analysis,

Animal behavior test
The animal behavior test of offspring was usually carried out at 7-8 weeks of age unless otherwise indicated. Autism-like behavior was evaluated using ultrasonic vocalization (USV), social interaction (SI) tests and a three-chambered social test [52][53][54], and detailed methods are described in Additional file 1.

Isolation of mouse intestine epithelial cells (IEC)
The protocol for isolation of IEC cells was based on the previously described method with minor modifications. In brief, the small and large intestines were harvested individually from treated mice and rinsed extensively with RPMI-1640 media (from Lonza) after Peyer's patches were removed (for small intestine). The rinsed intestines were opened longitudinally and macerated; the tissue pieces were shaken gently in RPMI-1640 containing 2 mM EDTA and 10% fetal calf serum. The tissue preparations were passed through 70-µm mesh filters, and the resulting single-cell suspensions were applied to Percoll (from Sigma) density gradients of 25%, 40%, and 75%. After centrifugation at 2000×g for 20 min, the interface between the 25% and 40% layers was collected to obtain IECs. The cells were stained using antibodies for either epithelial cell adhesion molecule (EpCAM, from Biolegend) or CD45 (from Biolegend) and nucleic acid dye (Via-Probe, from BD Biosiences). The Via-Probe − /CD45 − /EpCAM + IEC were sorted using BD FACSMelody TM Cell Sorter (BD Biosciences) for further biomedical analysis [55,56].

Intestinal permeability assay
The protocol was followed based on the previously described method with minor modifications. In brief, treated mice were fasted for 4 h before the experiment and then the FITC-dextran (50 mg/mL, Cat# 46944 from Sigma) was gavaged to mice (600 mg/kg). After 4 h, the whole blood was collected by cardiac puncture and placed at room temperature for 1 h before being centrifuged at a speed of 3000 rpm for 10 min. The supernatant was then transferred to a new tube for further centrifugation at a speed of 12,000 rpm for 10 min at 4 °C. The subsequent supernatant (serum) was diluted with equal volume of PBS and 100 µL diluted serum was added to a 96-cell microplate. The concentration of FITC in serum was determined at excitation/emission wavelengths of 485/530nm using a FLx800 microplate fluorescence reader (Bio-Tek). The serial diluted FITC-dextran (0, 0.5, 1, 2, 4, 6, 8, 10 µg/µL) was used as standards. Serum of mice administered with PBS was used as negative controls [17,57].

Fecal microbiome analysis
Fecal samples of the experimental mice were collected and stored at − 80 °C before being processed. The microbial DNA was extracted using a QIAamp Fast DNA Stool Mini Kit (from Qiagen) according to the manufacturer's protocol [58]. The purity and concentration of the extracted DNA were detected using agarose gel electrophoresis. Fecal microbiota was studied by performing V3-V4 16S rDNA amplicon sequencing in order to obtain the operational taxonomic units (OTU) defining the bacterial communities [59]. Sequencing samples from frozen fecal pellets were prepared, sequenced and subsequently processed using the MiSeq Pe300 Sequencing