Bo Zhou1,2,3†
Xiaoli Zheng1,2,3†Yunhua Chen4Xuehui Yan1,2,3Jinggang Peng1,2,3Yibu Liu1,2,3Yi Zhang1,2,3Lei Tang1,2,3*
Min Wen1,2,3*
- 1State Key Laboratory of Functions and Applications of Medicinal Plants, Guizhou Medical University, Guiyang, China
- 2Guizhou Provincial Engineering Technology Research Center for Chemical Drug R&D, Guizhou Medical University, Guiyang, China
- 3College of Pharmacy, Guizhou Medical University, Guiyang, China
- 4College of Basic Medical, Guizhou Medical University, Guiyang, China
Background: Some studies have shown that arginine vasopressin (AVP) can significantly improve the social interaction disorder of autism, but the mechanism remains unclear.
Methods: Female Wistar rats were intraperitoneally injected with VPA or normal saline at embryonic day 12.5 to establish an autism model or normal control in their offspring. Male offspring prenatally exposed to VPA were randomly assigned to two groups: the VPA-induced autism model group and the AVP group. The rats in the AVP group were treated with intranasal AVP at postnatal day (PND) 21 and for 3 weeks. The VPA-induced autism model group was given the same dose of normal saline in the same way. Behavioral responses were evaluated in the open field and three-chambered social test apparatus; the expression levels of AVP in serum were detected by enzyme-linked immunosorbent assay kit, and the gene expression levels on the amygdala were measured by RNA-seq at PND42.
Results: Intranasal administration of AVP can significantly improve the social interaction disorder and elevate the levels of AVP in serum. Transcriptome sequencing results showed that 518 differently expressed genes (DEGs) were identified in the VPA-induced autism model group compared with the control in this study. Gene Ontology biological process enrichment analysis of DEGs showed that the VPA-induced autism model group had significant nervous system developmental impairments compared with the normal group, particularly in gliogenesis, glial cell differentiation, and oligodendrocyte differentiation. Gene Set Enrichment Analysis (GSEA) enrichment analysis also showed that biological process of oligodendrocyte differentiation, axoneme assembly, and axon ensheathment were inhibited in the VPA-induced autism model group. Pathway enrichment analysis of DEGs between the control and VPA-induced autism model group showed that the PI3K/AKT and Wnt pathways were significantly dysregulated in the VPA-induced autism model group. Few DEGs were found when compared with the transcriptome between the VPA-induced autism model group and the AVP treatment group. GSEA enrichment analysis showed deficits in oligodendrocyte development and function were significantly improved after AVP treatment; the pathways were mainly enriched in the NOTCH, mitogen-activated protein kinase, and focal adhesion signaling pathways, but not in the PI3K/AKT and Wnt pathways. The expression patterns analysis also showed the same results.
Conclusion: AVP can significantly improve the social interaction disorder of VPA-induced autism model, and AVP may target behavioral symptoms in autism by modulating the vasopressin pathways, rather than primary disease mechanisms.
Introduction
Autism spectrum disorder (ASD) is a complex neurodevelopmental disorder primarily characterized by deficits in social interaction and communication as well as repetitive and stereotypic behavior. The current prevalence rate of ASD is 14.7 per 1,000 children in the United States (Christensen et al., 2016) and 7 per 1,000 in China (Zhou et al., 2020). It is more prevalent in males than in females, with reported ratios ranging from 2:1 to 5:1 (Lord et al., 2020). Although many researches have been carried out, the precise molecular mechanisms of ASD are still not fully elucidated, and no medications are currently approved for ameliorating ASD’s core social behavior deficits.
Arginine vasopressin (AVP) is a neuropeptide hormone synthesized and secreted by separate neuronal populations in the paraventricular nucleus (PVN) of the hypothalamic (Benarroch, 2013; Sparapani et al., 2021). It has been known for several decades that AVP plays a critical role in promoting the formation of complex mammalian social behavior such as social cognition, social recognition, social exploration, and preference by acting on AVP receptors which mainly located in hippocampus, amygdala, striatum, and other “social brains.” In recent years, a growing number of scientific reports have shown that AVP is a crucial factor in neurodevelopmental disorders, including the ASD. Oztan et al. (2018, 2020) (Parker et al., 2018) find that children with ASD have lower cerebrospinal fluid (CSF) AVP concentrations compared with control children and that patients with ASD with the lowest CSF AVP concentrations have the most severe symptoms. Zhang H. F. et al. (2016) find that the plasma AVP levels were associated with repetitive behavior, and ASD children with higher plasma AVP levels tended to have lower levels of repetitive. In addition, lots of researches have shown that correction of vasopressin deficit may ameliorate social deficits of autism (Borie et al., 2021). Parker et al. (2019) presented the safety and efficacy of 4 weeks’ intranasal AVP administration to improve social abilities in children with ASD using a double-blind, randomized, placebo-controlled trial design. Wu et al. (2021) studies have shown that infant rats exposed to VPA showed obviously impaired communication and repetitive behaviors with reduced number of AVP-immunoreactive cells in PVN and content of AVP in CSF. The postnatal subcutaneous injection of AVP can alleviate social preference deficits and stereotyped behaviors, accompanied with the increase in AVP concentrations in the CSF. Clearly, AVP signaling pathway may be a promising therapeutic target for autism. Today, the AVP has already entered phase II clinical trials (NCT01962870, NCT03204786) for the treatment of autism. Although intranasally administered AVP will observably achieve behavioral effects, the precise central mechanisms remain elusive.
To explore the potential for AVP as a treatment for autism and to further explore their mechanisms of action, we compared the amygdala transcriptome changes before and after AVP treatment in VPA-induced autism rat model.
Materials and Methods
Materials
Chemicals
Valproic acid sodium salt (P4543-10G) was purchased from Sigma-Aldrich Co. (St. Louis, MO, United States). Argipressin (AVP, HY-P0049) was purchased from MedChemExpress. RAT AVP enzyme-linked immunosorbent assay (ELISA) KIT (55R-1951) was purchased from Fitzgerald (United States).
Animals
Male and female Wistar rats weighing 270–290 g were obtained from the Department of Experimental Animal Center of Guizhou Medical University. Animals were housed individually with water and chow freely available under a regulated environment (23°C ± 2°C; 50% ± 10% humidity) with a 12/12-h light–dark cycle. All experiments were approved by the Guizhou Medical University Animal Care and Use Committee.
Methods
Animal Model
As previously described (Schneider and Przewlocki, 2005; Dai et al., 2018), female and male rats were allowed to mate overnight, and the morning when spermatozoa were found was designated as embryonic day 0.5 (E0.5). The pregnant rats were randomly distributed into two groups: VPA group (n = 10) and saline group (n = 5). On E12.5, rats in the VPA group were intraperitoneally injected with sodium valproate (dose of 600 mg/kg, 250 mg/mL dissolved in physiological saline); the saline groups received the same volume of normal saline at the same time. The day of birth of the offspring was marked as postnatal day 1 (PND1). After weaning at PND21, offspring of the same sex were housed separately, with four to five per cage. To assess the treatment effects of AVP, the offspring of VPA group were randomly divided into two groups: VPA-induced autism model groups (n = 10) and AVP treatment group (n = 10). AVP treatment group received a daily intranasal of AVP (dose of 400 μg/kg, 2.5 mg/mL dissolved in NS) from PND21 to PND42. The offspring of saline groups were marked as control group. The control groups and the VPA-induced autism model group were given the same amount of saline. All experiments were carried out on male offspring. The experimental procedure is shown in Figure 1.
Figure 1. Schematic representation of the experimental procedure.
Development and Behavioral Tests
Open Eyes Test
To examine the maturation process in the pups, we monitored eye opening status from PND12–16 once daily and scored as follows: 0 = both eyes closed, 1 = one eye open, and 2 = both eyes open.
Swimming Performance
The swimming test measures motor development and integration of coordinated series of reflex responses (Schneider and Przewlocki, 2005). It was conducted on PND8, PND10, PND12, PND14, and PND16. Each animal was put at the center of an container filled with water (28–29°C) and was observed for 5–10 s. Swimming performance was scored as follows: 0 = head and nose below the surface; 1 = nose below the surface; 2 = nose and top of head at or above the surface, but ears still below the surface; 3 = the same as in 2 except that water line was at mid-ear level; and 4 = the same as in 3 except that water line was at the bottom of ears.
Open Field Test
To test for generalized anxiety, the open field test was used at PND42. The tested rat was placed in the center of a squared open field box (80 cm L × 80 cm W × 40 cm H) and allowed to explore freely for 5 min. Locomotion trajectory of the rat was automatically detected and was recorded, and the percentage of the time spent in the central zone in 5 min was calculated.
Three-Chamber Test
The three-chamber test was performed at PND42 as previously described (Rein et al., 2020). A polyvinyl chloride box was divided in three compartments (each chamber is 40 cm L × 40 cm W × 40 cm H), and both side compartments contained an empty perforated cup. First, the tested mouse was allowed to explore freely the whole setting, with all doors open for 10 min (phase 1). After this habituation period, the mouse was restricted in the central compartment, whereas an unfamiliar rat of the same sex (stranger 1) was placed under one of the cups (sides alternated between each rat). The tested mouse was then allowed to explore the whole apparatus for 10 min (phase 2). After that, it was restricted to the central compartment, whereas another unfamiliar rat of the same sex (stranger 2) was placed under the other cup. The tested mouse could then again freely explore the whole apparatus for 10 min (phase 3). In all three phases, time spent in each compartment was manually recorded, and the social preference indexes (Rein et al., 2020) were calculated as follows: social preference indexes = (TS − TNS)/(TS + TNS); TS, social stimulus, the interaction time with strange 1 in phase 2 and the interaction time with strange 2 in phase 3; TNS, non-social stimulus, the interaction time in empty in phase 2 and the interaction time with strange 1 in phase 3. To minimize the impact from residual rat odors, the entire apparatus was thoroughly cleaned with 70% ethanol at the beginning of each trial.
Self-Grooming Test
The procedure of self-grooming was acquired as described previously, to assess the repetitive behavior at PND42 (Mahmood et al., 2020). In brief, a standard cage of the rat was used to measure repetitive behavior in the self-grooming test. The dimensions of cages were 25 cm wide, 45 cm long, and 20 cm high. After a 5-min habituation in the cage, the cumulative grooming time spent for all body parts of each mouse was calculated by using a stopwatch for 5 min.
The Expression Levels of Arginine Vasopressin in Serum
Arginine vasopressin content in the serum was assayed using a commercially available ELISA kit (Fitzgerald, 55R-1951).
RNA-seq
Sampling and RNA Isolation
Rats were decapitated, and the amygdala was quickly dissected and frozen in liquid nitrogen for 2 h. The frozen amygdalae were stored at −80°C until used. Total RNA was extracted and purified using TRIzol reagent (Invitrogen, Carlsbad, CA, United States) following the manufacturer’s procedure. The RNA concentration and purity of each sample were quantified using NanoDrop ND-1000 (NanoDrop, Wilmington, DE, United States). RNA quality was measured by Bioanalyzer 2100 (Agilent, Santa Clara, CA, United States). All RNA samples included in the expression analysis had an A260/A280 absorbance ratio greater than 1.8 and RNA integrity number > 7.0.
RNA-seq
An RNA-seq library of each sample was prepared, and 2 × 150-bp paired-end sequencing (PE150) was performed on an Illumina Novaseq™ 6000 (LC-Bio Technology Co., Ltd., Hangzhou, China) following the vendor’s recommended protocol.
Sequence Analysis
Raw data files in FASTQ format were generated from the Illumina Novaseq™ 6000. The reads after quality control and preprocessing were aligned to the reference genome (rattus_norvegicus6.0, v101) using the HISAT2 (Kim et al., 2019) and gene expression quantified using HTSEQ (Anders et al., 2015). The differentially expressed genes (DEGs) were selected with fold changes (FCs) > 1.3 and with statistical significance [false discovery rate (FDR) < 0.05] by DESeq2 (Love et al., 2014); the Gene Ontology (GO) enrichment and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment of the DEGs were analyzed using clusterProfiler (Yu et al., 2012). Gene Set Enrichment Analysis (GSEA) was performed to identify the sets of related genes that might be systematically altered in each group by GSEA V4.1.0 (Mootha et al., 2003; Subramanian et al., 2005). The biological process (c5.go.bp.v7.4) and KEGG term (c2.cp.kegg.v7.4) were annotated by Molecular Signatures Database v7.4 (Liberzon et al., 2011, 2015). The significantly enriched gene sets were selected with |NES|>1 (Normalized enrichment score) and with NOM p < 0.05 (Nominal p-value). Soft clustering was performed using Mfuzz (Kumar and Futschik, 2007) to mine the expression patterns of genes in three groups. In the clustering analysis, the optimal cluster number was calculated using default parameters. Meanwhile, the min score (membership) threshold was set to 0.6. Each cluster was subsequent functional annotation analysis using the online tool g:Profiler1 (Raudvere et al., 2019).
Quantitative Reverse Transcription–Polymerase Chain Reaction Confirmation
Total RNA from amygdala was extracted utilizing TRIzol reagent (Invitrogen) and reverse transcribed to cDNA using the Thermo Scientific Revert Aid First Strand cDNA Synthesis Kit (Thermo Fisher Scientific Inc., Waltham, MA, United States). Relative expression level of target mRNA was normalized to GAPDH, and relative expression ratio of a target gene was calculated using the mRNA by the 2–ΔΔCT method (Livak and Schmittgen, 2001). The primer sequences used are listed in Table 1.
Table 1. Primer sequences.
Statistical Analysis
All data are represented as the mean ± SD. Statistical analysis of multiple-group comparisons was performed by one-way analysis of variance, and two-group comparisons were analyzed by the two-tailed Student t test using STATA 14.0 (StataCorp, College Station, TX, United States). Statistical significance was set at p < 0.05.
Results
Development and Behavioral Tests
Eye Opening
The eye opening scores between the two groups were significantly different on PND13 (p < 0.05) and PND14 (p < 0.05) with significantly lower eye opening score in the VPA group. On PND12, PND15, and PND16, the eye opening score was similar, and no statistically significant difference was seen between the two groups (Figure 2). The results indicate that VPA may significantly delay the maturation of rats.
Figure 2. Eye opening test in each group (n = 6 in each group). The eye opening score was significantly delayed in the VPA group on PNDs 13 and 14. Compared with the control group: *p < 0.05.
Swimming Performance
Compared with the control group, the ontogeny of swimming behavior was significantly delayed in VPA rats on PND9 (p < 0.01) and PND11 (p < 0.01), and no significant differences were found on PND13 and PND15 (Figure 3). The results indicate that VPA may significantly delay motor development and attenuate integration of coordinated series of reflexes.
Figure 3. Swimming test in each group (n = 6 in each group). The swimming score was significantly delayed in the VPA group on PND9 and PND11. Compared with the control group: **p < 0.01.
Open Field Test
When compared with the control group, the anxiety-like behavior was significantly enhanced (it spent more time staying in the edge and corner and less time in the central zone than did the control group) in the VPA-induced autism model group (p < 0.01), and the anxiety-like behavior was significantly decreased after the treatment with AVP (p < 0.01; Figure 4).
Figure 4. Open field test (n = 6 in each group). (A) locomotion trajectory of control group, (B) locomotion trajectory of VPA-induced autism model group, (C) locomotion trajectory of AVP group, (D) compare the time spend in center in each group. Compared with the control group: **p < 0.01; compared with the VPA-induced autism model group: ##p < 0.01.
Three-Chamber Test
In social preference test (phase2), each group rats spent more time in the side with stranger 1 than in the empty cage (p < 0.01). The time in the side with stranger 1 were significantly lower in the VPA-induced autism model group than control (p < 0.01) and AVP group (p < 0.01; Figure 5A). The social preference index in social preference test was significantly decreased in the VPA-induced autism model group compared with the control group (p < 0.01) and AVP group (p < 0.01; Figure 5B). In the social novelty test (phase 3), the control and AVP groups spent more time in the side with stranger 2 than in the side with stranger 1 (p < 0.01). The VPA-induced autism rats seemed to lose interest in social novelty with conspecifics, spending a comparable time exploring the cage containing stranger 1 rat and stranger 2 (p > 0.01; Figure 5C). Furthermore, the social preference index in social novelty test was significantly decreased in the VPA-induced autism model group compared with the control group and AVP group (p < 0.01; Figure 5D). The results suggest that the social interaction was impaired in the VPA-induced autism model group, and the impairment was significantly improved following AVP treatment.
Figure 5. The results of three-chamber test in each group (n = 6 in each group). (A) The results of social preference test in each group. Comparison of the time spent in the side with stranger 1 vs. the side with empty: *p < 0.05, **p < 0.01. (B) Social preference index in social preference test. Compared with the control group: **p < 0.01; compared with the autism group: ##p < 0.01. (C) The results of social novelty test in each group. Comparison of the time spent in the side with stranger 1 vs. the side with stranger 2: *p < 0.05, **p < 0.01. (D) Social preference index in social novelty test. Compared with the control group: **p < 0.01; compared with the autism group: ##p < 0.01.
Self-Grooming Test
To test for repetitive/stereotypic behavior, the self-grooming test was used. When compared with the control group, the cumulative self-grooming time was significantly prolonged in the VPA-induced autism model group (p < 0.01), and the cumulative time was significantly shorter in the AVP treatment group (p < 0.01) when compared with the VPA-induced autism rats (Figure 6).
Figure 6. Repetitive/stereotypic behavior test in rats (n = 6 in each group). Compared with the control group: **p < 0.01; compared with the VPA-induced autism model group: ##p < 0.01.
The Expression Levels of Arginine Vasopressin in Serum
Compared with the control group, the AVP level in the serum significantly decreased in the VPA-induced autism model group (p < 0.01). After AVP treatment, the AVP levels in the serum were significantly increased (p < 0.01; Figure 7).
Figure 7. The expression levels of serum AVP in each group (n = 5 in each group). Compared with the control group: **p < 0.01; compared with the VPA-induced autism model group: ##p < 0.01.
RNA-seq
Quality Control and Read Mapping
The number of total reads, clean reads, the sequencing error rate, and the percentage of Q20, Q30, and the GC content of all samples met the quality control requirements of sequencing. More than 95% of the clean reads were mapped to the rat genome, and more than 75% of them were mapped to a single genome location.
Differentially Expressed Analysis
There are 518 genes (216 up and 302 down) whose expression was significantly different with an adjusted FDR p < 0.05 and FC > 1.3 (only 20 upregulated and 7 downregulated genes were found in FC > 2) in the VPA-induced autism model group compared with the control (the heatmap of the top 50 DEGs are shown in Figure 8A). We performed GO biology process (BP) and KEGG pathway functional enrichment analyses to explore the potential biological functions of these differential expression genes. The results of GO BP were significantly enriched in nervous system development, tissue development and remodeling, and extracellular structure organization (Figures 8B,C). Clearly, the VPA-induced autism model group had significant developmental impairments compared with the control group, particularly in the nervous system [gliogenesis, glial cell differentiation, and oligodendrocyte (OL) differentiation]. Further analysis showed that neurodevelopmental disorders and tissue developmental disorders were associated with down-regulated genes, whereas extracellular matrix remodeling was associated with up-regulated genes (Figures 8D,E). The KEGG pathway enrichment analysis showed that these differential expression genes were significantly enriched in six pathways, such as PI3K-Akt signaling pathway, Wnt signaling pathway, protein digestion and absorption, and so on (Figure 8F). We also analyzed the intersection between our differential expression gene list and the list of 1,010 genes that have evidence of genetic association with ASD from the SFARI human gene list2. Of 518 differential expression genes, 27 were found in the SFARI database, which are mainly involved in nervous system development, synapse maturation, synapse organization, protein localization to synapse, glial cell migration, and so on.
Figure 8. Differentially expressed analysis in autism group compared with the control group (n = 5 in each group). (A) Hierarchical clustering heatmap of the top 50 DEGs in the VPA-exposed group compared with the control group. (B) GO biology process enrichment analysis of the DEGs. The x axis indicates the enrichment ratio, and the y axis indicates the GO terms. The size of the circle represents the gene number, and the color of the circle indicates the value (adjusted p value). (C) Directed acyclic graph of GO biology process in DEGs, the depth of the color represents the degree of enrichment. (D,E) GO biology process enrichment analysis of down-regulated and up-regulated genes. The size of the circle represents the gene number, and the color of the line indicates clustering of biological processes (F) KEGG enrichment analysis of DEGs. The x axis indicates the enrichment ratio, and the y axis indicates the KEGG terms. The size of the circle represents the gene number, and the color of the circle indicates the value (adjusted p value).
There are nine genes (seven up and two down) whose expression was significantly different with an adjusted FDR p < 0.05 and FC > 1.3 (no significantly different genes were found in FC > 2) in the VPA-induced autism model group compared with the AVP group. The number of DEGs between the VPA-induced autism model and the AVP group was limited and too few for confident GO and KEGG annotation.
Gene Set Enrichment Analysis
Traditional KEGG analysis strategies usually focus on a handful of gene that exhibit differences between two states of interest. Although useful, they are easily affected by the filter threshold (logFC and FDR) and lead to miss some genes with moderate differential expression but important biological significance. To overcome many genes with moderate but meaningful expression changes are discarded by the strict cutoff value, which leads to a reduction in statistical power, we conducted a GSEA (Subramanian et al., 2005). From a biological perspective, GSEA methods are promising because functionally related genes often display a coordinated expression to accomplish their roles in the cell. In the VPA-induced autism model group compared with the control, the results of GSEA GO BP enrichment analysis showed that 239 gene sets are significantly upregulated in the VPA-induced autism model group (NOM p < 0.05, |NES|>1). The top 5 terms were significantly enriched in mitochondrial electron transport NADH to ubiquinone, ATP synthesis coupled electron transport, positive regulation of telomerase RNA localization to Cajal body, regulation of cellular amino acid metabolic process, and oxidative phosphorylation. Two hundred eighty-seven gene sets are significantly downregulated in the VPA-induced autism model group (NOM p < 0.05, |NES|>1). The top 5 GO terms were significantly enriched in OL differentiation, axoneme assembly, axon ensheathment in central nervous system (CNS), negative regulation of glial cell differentiation, and microtubule bundle formation. The results of GSEA KEGG enrichment analysis showed that 18 gene sets are significantly upregulated in the VPA-induced autism model group (NOM p < 0.05, |NES|>1). The top 5 terms were significantly enriched in proteasome, Parkinson’s disease, oxidative phosphorylation, Alzheimer’s disease, and Huntington disease. Four gene sets (base excision repair, ether lipid metabolism, notch signaling pathway, and peroxisome proliferator–activated receptor [PPAR] signaling pathway) were significantly downregulated in the VPA-induced autism model group (NOM p < 0.05, |NES|>1).
In the VPA-induced autism model group compared with the AVP group, the results of GSEA GO BP enrichment analysis showed that 96 gene sets are significantly upregulated in the VPA-induced autism model group (NOM p < 0.05, |NES|>1). The top 5 GO terms were significantly enriched in neuropeptide signaling pathway, γ-aminobutyric acid (GABA) signaling pathway, synaptic transmission GABAergic, chromatin remodeling at centromere, and establishment of mitochondrion localization. Six hundred sixty-six gene sets are significantly downregulated in the VPA-induced autism model group (NOM p < 0.05, |NES| > 1). The top 5 GO terms were significantly enriched in positive regulation of response to cytokine stimulus, positive regulation of interleukin 4 production, cellular extravasation, cell differentiation involved in metanephros development, and regulation of lymphocyte homeostasis. The results of GSEA KEGG enrichment analysis showed that nine gene sets are significantly upregulated in the VPA-induced autism model group (NOM p < 0.05, |NES|>1). The top 5 terms were significantly enriched in neuroactive ligand–receptor interaction, O glycan biosynthesis, aminoacyl tRNA biosynthesis, proteasome, and citrate cycle (tricarboxylic acid cycle). Twenty-seven gene sets are significantly downregulated in the VPA-induced autism model group (NOM p < 0.05, |NES|>1). The top 5 pathways were significantly enriched in complement and coagulation cascades, Leishmania infection, cytokine–cytokine receptor interaction, prostate cancer, and basal cell carcinoma.
The intersection of the gene expression variation between down in the VPA-induced autism model group (con vs. mod) and up in AVP (mod vs. AVP) was analyzed by Evenn3. The results showed that 102 BPs were remarkably downregulated in the VPA-induced autism model group and upregulated in AVP treatment, such as OL differentiation, glial cell differentiation, gliogenesis, and so on (Figures 9A–F). Furthermore, some neurodevelopment-related BPs such as glial cell proliferation, glial cell migration, neural tube development, and so on were not changed in the VPA-induced autism model group but remarkably upregulated after AVP treatment. Only notch signaling pathway was found remarkably downregulated in the VPA-induced autism model group and upregulated in AVP treatment. Some neurodevelopment-related pathways, such as mitogen-activated protein kinase (MAPK) signaling pathway, focal adhesion, and so on, were not changed in the VPA-induced autism model group but remarkably upregulated after AVP treatment (Figures 10A–F).
Figure 9. Enrichment of related biological process by GSEA. (A) Gliogenesis in control (con) vs. autism (mod), (B) glial cell differentiation in control (con) vs. autism (mod), (C) oligodendrocyte differentiation in control (con) vs. autism (mod), (D) GLIOGENESIS in AVP vs. autism (mod), (E) glial cell differentiation in AVP vs. autism (mod), and (F) oligodendrocyte differentiation in AVP vs. autism (mod). The significantly enriched gene set was selected with |NES|>1 and with NOM p < 0.05.
Figure 10. Enrichment of related signaling pathway by GSEA. (A) Notch signaling pathway in control (con) vs. autism (mod). (B) MAPK signaling pathway in control (con) vs. autism (mod). (C) Focal adhesion signaling pathway in control (con) vs. autism (mod). (D) Notch signaling pathway in AVP vs. autism (mod). (E) MAPK signaling pathway in AVP vs. autism (mod). (F) Focal adhesion signaling pathway in AVP vs. autism (mod). The significantly enriched gene set was selected with |NES|>1 and with NOM p < 0.05.
Dynamic Expression Patterns of Genes
The expression trends of genes at each group were analyzed using Mfuzz clustering analysis, and eight clusters with different change trends were screened out (Figure 11). Genes in clusters 1 and 8 presented a trend of increasing in the VPA-induced autism model group and decreasing in the AVP group. Genes in clusters 5 and 7 presented a trend of decreasing in the VPA-induced autism model group and increasing in the AVP group. Genes in clusters 2 and 6 presented a trend of gradually declining in the VPA-induced autism model group and AVP group. Genes in clusters 3 and 4 presented a trend of continuous increase in the VPA-induced autism model group and AVP group.
Figure 11. The series of diagrams illustrates the patterns of dynamic changes of gene in each group. Clusters 1 and 8 presented a trend of increasing in the VPA-induced autism model group and decreasing in the AVP group. Clusters 5 and 7 presented a trend of decreasing in the VPA-induced autism model group and increasing in the AVP group. Clusters 2 and 6 presented a trend of gradually declining in the VPA-induced autism model group and AVP group, and clusters 3 and 4 presented a trend of continuous increase in the VPA-induced autism model group and AVP group.
We further analyzed the gene biology function in clusters 5 and 7. A large number of development-related biological processes are enriched in cluster 5, such as neurodevelopment, neurogenesis, gliogenesis, glial cell differentiation, and migration, and these genes all point together to notch signal pathway. In cluster 7, lots of genes were enriched in the cellular developmental process, cell differentiation, cell adhesion, cell migration, and so on, and these genes were mainly involved in the focal adhesion pathway and MAPK pathway. The results of expression trend analysis were basically consistent with those of analysis of GSEA.
The mRNA Levels of Oligodendrocyte and Myelin Development-Related Genes
Compared with the control group, the mRNA levels of OL and myelin development-related genes, such as Sox10 (a transcription factor that directs neural stem cells toward the glial lineage), Olig2 (a transcription factors necessary for OL development), and MBP (a structural component of myelin, expressed exclusively by myelinating glia), were significantly decreased in the VPA-induced autism model group (p < 0.05), and the mRNA levels were improved after AVP treatment (p < 0.05) (Figure 12).
Figure 12. The mRNA levels of oligodendrocyte and myelin development-related genes (n = 5 in each group). Compared with the control group: **p < 0.01; compared with the VPA-induced autism model group: ##p < 0.01.
Discussion
Our study showed that rats exposed to VPA at E12.5 had significantly delayed growth development and impaired social behaviors, and the social behaviors were improved significantly by AVP treatment. The results of this study are consistent with previous reports (Wu et al., 2021).
As a hub in the “social brain,” the amygdala, which was located within the anterior medial portion of the brain’s temporal lobe, maintains diverse connections with the neocortex, basal ganglia, hippocampus, thalamus, and other brain cortical areas via a network of white matter tracts, responsible for information processing that subserve emotional learning, social cognition, and social interaction (Baron-Cohen et al., 2000; Zalla and Sperduti, 2013; Ferrara et al., 2021; Raam and Hong, 2021; Xu et al., 2021). A large number of studies have demonstrated that amygdala structural and/or functional abnormalities are associated with social malfunctioning in ASD (Shen et al., 2016; Avino et al., 2018; Herrero et al., 2020; Xu et al., 2020; Seguin et al., 2021; Vitor-Vieira et al., 2021). So, we speculated that the prosocial behavior of AVP may be related to the amygdala. In order to further explore its mechanism of action, we focused on transcriptome analysis in the amygdala after AVP treatment.
Deficits in Oligodendrocyte Development and Function Potentially Represent a Common Molecular Pathway Disrupted in Autism Spectrum Disorder
Five hundred eighteen DEGs were identified in the VPA-induced autism model group compared with the control in this study. GO biological process enrichment analysis of DEGs showed that the VPA-induced autism model group had significant nervous system developmental impairments compared with the normal group, particularly in gliogenesis, glial cell differentiation, and OL differentiation. GSEA enrichment analysis also showed that the biological process of OL differentiation, axoneme assembly, and axon ensheathment were inhibited in the VPA-induced autism model group. Further analysis showed that neurodevelopmental disorders were associated with the transcription inhibition of OL development and myelination-related genes. The results were consistent with the research of Zhang X. F. et al. (2020) who found the transcript levels of myelin-related genes were reduced in poly(I:C)-exposed mouse offspring in the MIA model. Obviously, deficits in OL development and function may be one of the core mechanisms of autism. This may be related to VPA directly inhibits histone deacetylase (Kataoka et al., 2013), causing transient hyperacetylation in the brain, inhibiting development-related genes transcription (Kataoka et al., 2013; Taleb et al., 2021), causes a delay in OL differentiation and hypomyelination. In addition, VPA exposure may impairs repair of DNA damage (Servadio et al., 2018), modifies cholesterol/isoprenoid metabolism, and reduces the number of OLs leading to lower myelin and cholesterol levels (Cartocci et al., 2018).
OLs are glial cells in the CNS, responsible for myelin sheath formation, which allows fast signal transmission, provides metabolic support to axons (Cherchi et al., 2021), and contributes to the optimal information processing in complex neural networks to maintaining the functional connectome of the brain. Deficits in OL development and function will lead to demyelination and axonal dysregulation, and disruption in neuron–glia interactions promotes autistic-like features (Bronzuoli et al., 2018; Duncan et al., 2021). Ample evidence has revealed that OL development disorders are closely related to autism. Graciarena et al. (2018) found that reduction in the proliferation of oligodendroglial cells and low levels of myelin basic protein has also been implicated in ASD pathogenesis. The analyses of DEGs highlighted OL dysregulation, which we confirmed in two additional mouse models of syndromic ASD (Phan et al., 2020). Zhang X. F. et al. (2020) found oligodendroglia-related gene mRNA, which are critically involved in myelination and OL progenitor differentiation, is significantly decreased in amygdala in maternal immune activation ASD model. Astrocytes and OLs may contribute to neurochemical imbalances described in the autistic brain by disrupting neurotransmission or modifying axonal conduction (Voineagu et al., 2011; Gupta et al., 2014).
PI3K/AKT and Wnt Pathway May Be the Core Mechanism of Autism
Pathway enrichment analysis of DEGs between the control group and VPA-induced autism model group showed that the PI3K/AKT and Wnt pathways were significantly dysregulated in the VPA-induced autism model group.
PI3K/AKT Signaling Pathway Represents an Essential Signaling Mechanism for Mammalian Enzyme-Related Receptors
PI3K/AKT signaling pathway represents an essential signaling mechanism for mammalian enzyme-related receptors in transducing signals or biological processes such as cell development, differentiation, cell survival, protein synthesis, and metabolism (Sharma and Mehan, 2021). Up to now, there is still not a unified conclusion on the expression changes of PI3K/Akt pathway in autism. Nicolini et al. (2015) find that PI3K–AKT–mTOR signaling pathways were downregulated in idiopathic autism fusiform gyrus and in neocortex of valproic acid–induced rat model. Luo et al. (2020) find that the PI3K/Akt pathway was inactivated in BTBR mice; overexpressed moesin significantly restored the activity of PI3K/Akt pathway and improved the autistic-like behaviors. Tylee et al. (2017) found that PI3K–AKT–mTOR and RAS–MAPK signaling cascades were diminished, and ribosomal translation and natural killer cell–related activity increased in ASD blood transcriptome. Conversely, Zhang J. et al. (2016) find that p-PI3K and p-AKT significantly increased in the hippocampus in autism rat model, and inhibition of mTOR by NVP-BEZ235 significantly reduced the activity of PI3K/Akt and improved social interaction in VPA-induced ASD. Zhang W. et al. (2020) find that small GTPase gene RAB39b mutation promotes PI3K–AKT–mTOR activity and alters cortical neurogenesis, leading to macrocephaly and autistic-like behaviors. The inconsistent results may be due to the sample size, different brain regions, different developmental time, and different research methods.
Wnt/β-Catenin Signaling Pathway Plays a Central Role in Neurodevelopment
Wnt/β-catenin signaling pathway plays a central role in neurodevelopment, and perturbation Wnt signaling may trigger the advent of disorders related to the structures and functions of the CNS (Inestrosa and Varela-Nallar, 2015; Kumar et al., 2019). Zhang et al. (2012) (Qin et al., 2016; Kuo and Liu, 2018) find that prenatal administration of VPA causes the upregulation of the Wnt/β-catenin signaling pathway, which facilitates susceptibility to autism. Sulindac treatment ameliorates autism-like behavioral phenotypes probably at least in part via the downregulation of the canonical Wnt/b-catenin (Zhang et al., 2015). Vallee et al. (2019) found Wnt/β-catenin pathway was upregulated in autism, and PPARγ agonists may play potential treatment for ASD by inhibiting the canonical Wnt/β-catenin pathway. Emery (2010) found the Wnt/β-catenin pathway plays an important role in OL development; mutant mice with elevated Wnt/β-catenin signaling in the OL lineage display blocked differentiation and hypomyelination. It is obvious that the PI3K/AKT and Wnt pathways may be the core mechanism of autism; targeting these pathways can be beneficial for the improvement of autism (Sharma and Mehan, 2021).
Arginine Vasopressin Improves the Autism-Like Behavior Through Non-core Pathways
Although AVP treatment significantly improved social interaction disorders in VPA-induced autism model rats, the mechanism remains unclear. Few DEGs were found when compared with the transcriptome between the VPA-induced autism model group and AVP treatment group. GSEA enrichment analysis showed that deficits in OL development and function were significantly improved after AVP treatment, the pathways were mainly enriched in the NOTCH, MAPK, focal adhesion signaling pathways, but not in the PI3K/AKT and Wnt pathway. The expression patterns analysis also showed the same results. Thus, we speculate that AVP can compensate for the dysregulation of PI3K/AKT and Wnt pathways by regulating the NOTCH, MAPK, and focal adhesion signaling pathways, thus promoting OL development and myelin formation and improving autism-like behavior. A large number of studies have confirmed that activation of NOTCH (Kim et al., 2020; Li et al., 2021), MAPK (Ishii et al., 2016; Lorenzati et al., 2021), and Focal adhesion (Forrest et al., 2009; Lafrenaye and Fuss, 2010) signaling pathway can regulate the proliferation, differentiation, and migration of OL and myelin formation. In addition, crosstalk between these pathways and PI3K/AKT and Wnt pathway has been reported in many neurological or non-neurological disorders (Bai et al., 2019; Ishii et al., 2019, 2021; Worthmuller and Ruegg, 2020; Acar et al., 2021).
Conclusion
Collectively, these results demonstrate that AVP can significantly improve the social interaction disorder of VPA-induced autism model, and AVP may target behavioral symptoms in autism by modulating the vasopressin pathways, rather than primary disease mechanisms.
Data Availability Statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/supplementary material.
Ethics Statement
The animal study was reviewed and approved by Guizhou Medical University.
Author Contributions
BZ, LT, and MW: conceptualization, writing, supervision, project administration, and funding acquisition. XZ, JP, and YZ: animal model. XZ, XY, and YL: behavioral tests. BZ and YC: RNA-seq data analysis. All authors read and approved the final manuscript.
Funding
This study was supported by Ph.D. research startup foundation of Guizhou Medical University [(2020)008 and (2020)010].
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.
Publisher’s Note
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.
Acknowledgments
We are grateful to Guizhou Provincial Engineering Technology Research Center for Chemical Drug R&D.
Footnotes
- ^ https://biit.cs.ut.ee/gprofiler/gost
- ^ https://gene.sfari.org/database/human-gene/
- ^ http://www.ehbio.com/test/venn/#/
References
Acar, A., Hidalgo-Sastre, A., Leverentz, M. K., Mills, C. G., Woodcock, S., Baron, M., et al. (2021). Inhibition of Wnt signalling by Notch via two distinct mechanisms. Sci. Rep. 11:9096. doi: 10.1038/s41598-021-88618-5
Anders, S., Pyl, P. T., and Huber, W. (2015). HTSeq–a Python framework to work with high-throughput sequencing data. Bioinformatics 31, 166–169. doi: 10.1093/bioinformatics/btu638
Avino, T. A., Barger, N., Vargas, M. V., Carlson, E. L., Amaral, D. G., Bauman, M. D., et al. (2018). Neuron numbers increase in the human amygdala from birth to adulthood, but not in autism. Proc. Natl. Acad. Sci. U.S.A. 115, 3710–3715. doi: 10.1073/pnas.1801912115
Bai, C., Zhang, H., Zhang, X., Yang, W., Li, X., and Gao, Y. (2019). MiR-15/16 mediate crosstalk between the MAPK and Wnt/beta-catenin pathways during hepatocyte differentiation from amniotic epithelial cells. Biochim. Biophys. Acta Gene Regul. Mech. 1862, 567–581. doi: 10.1016/j.bbagrm.2019.02.003
Baron-Cohen, S., Ring, H. A., Bullmore, E. T., Wheelwright, S., Ashwin, C., and Williams, S. C. (2000). The amygdala theory of autism. Neurosci. Biobehav. Rev. 24, 355–364. doi: 10.1016/s0149-7634(00)00011-7
Benarroch, E. E. (2013). Oxytocin and vasopressin: social neuropeptides with complex neuromodulatory functions. Neurology 80, 1521–1528. doi: 10.1212/WNL.0b013e31828cfb15
Borie, A. M., Dromard, Y., Guillon, G., Olma, A., Manning, M., Muscatelli, F., et al. (2021). Correction of vasopressin deficit in the lateral septum ameliorates social deficits of mouse autism model. J. Clin. Invest. 131:e144450. doi: 10.1172/JCI144450
Bronzuoli, M. R., Facchinetti, R., Ingrassia, D., Sarvadio, M., Schiavi, S., Steardo, L., et al. (2018). Neuroglia in the autistic brain: evidence from a preclinical model. Mol. Autism 9:66. doi: 10.1186/s13229-018-0254-0
Cartocci, V., Catallo, M., Tempestilli, M., Segatto, M., Pfrieger, F. W., Bronzuoli, M. R., et al. (2018). Altered brain cholesterol/isoprenoid metabolism in a rat model of autism spectrum disorders. Neuroscience 372, 27–37. doi: 10.1016/j.neuroscience.2017.12.053
Cherchi, F., Pugliese, A. M., and Coppi, E. (2021). Oligodendrocyte precursor cell maturation: role of adenosine receptors. Neural Regen. Res. 16, 1686–1692. doi: 10.4103/1673-5374.306058
Christensen, D. L., Baio, J., Van Naarden Braun, K., Bilder, D., Charles, J., Constantino, J. N., et al. (2016). Prevalence and characteristics of autism spectrum disorder among children aged 8 years–autism and developmental disabilities monitoring network, 11 Sites, United States, 2012. MMWR Surveill. Summ. 65, 1–23. doi: 10.15585/mmwr.ss6802a1
Dai, Y. C., Zhang, H. F., Schon, M., Bockers, T. M., Han, S. P., Han, J. S., et al. (2018). Neonatal oxytocin treatment ameliorates autistic-like behaviors and oxytocin deficiency in valproic acid-induced rat model of autism. Front. Cell Neurosci. 12:355. doi: 10.3389/fncel.2018.00355
Duncan, G. J., Simkins, T. J., and Emery, B. (2021). Neuron-oligodendrocyte interactions in the structure and integrity of axons. Front. Cell Dev. Biol. 9:653101. doi: 10.3389/fcell.2021.653101
Emery, B. (2010). Regulation of oligodendrocyte differentiation and myelination. Science 330, 779–782. doi: 10.1126/science.1190927
Ferrara, N. C., Trask, S., Avonts, B., Loh, M. K., Padival, M., and Rosenkranz, J. A. (2021). Developmental shifts in amygdala activity during a high social drive state. J. Neurosci. 41, 9308–9325. doi: 10.1523/JNEUROSCI.1414-21.2021
Forrest, A. D., Beggs, H. E., Reichardt, L. F., Dupree, J. L., Colello, R. J., and Fuss, B. (2009). Focal adhesion kinase (FAK): a regulator of CNS myelination. J. Neurosci. Res. 87, 3456–3464. doi: 10.1002/jnr.22022
Graciarena, M., Seiffe, A., Nait-Oumesmar, B., and Depino, A. M. (2018). Hypomyelination and oligodendroglial alterations in a mouse model of autism spectrum disorder. Front. Cell Neurosci. 12:517. doi: 10.3389/fncel.2018.00517
Gupta, S., Ellis, S. E., Ashar, F. N., Moes, A., Bader, J. S., Zhan, J., et al. (2014). Transcriptome analysis reveals dysregulation of innate immune response genes and neuronal activity-dependent genes in autism. Nat. Commun. 5:5748. doi: 10.1038/ncomms6748
Herrero, M. J., Velmeshev, D., Hernandez-Pineda, D., Sethi, S., Sorrells, S., Banerjee, P., et al. (2020). Identification of amygdala-expressed genes associated with autism spectrum disorder. Mol. Autism 11:39. doi: 10.1186/s13229-020-00346-1
Inestrosa, N. C., and Varela-Nallar, L. (2015). Wnt signalling in neuronal differentiation and development. Cell Tissue Res. 359, 215–223. doi: 10.1007/s00441-014-1996-4
Ishii, A., Furusho, M., and Bansal, R. (2021). Mek/ERK1/2-MAPK and PI3K/Akt/mTOR signaling plays both independent and cooperative roles in Schwann cell differentiation, myelination and dysmyelination. Glia 69, 2429–2446. doi: 10.1002/glia.24049
Ishii, A., Furusho, M., Dupree, J. L., and Bansal, R. (2016). Strength of ERK1/2 MAPK activation determines its effect on myelin and axonal integrity in the adult CNS. J. Neurosci. 36, 6471–6487. doi: 10.1523/JNEUROSCI.0299-16.2016
Ishii, A., Furusho, M., Macklin, W., and Bansal, R. (2019). Independent and cooperative roles of the Mek/ERK1/2-MAPK and PI3K/Akt/mTOR pathways during developmental myelination and in adulthood. Glia 67, 1277–1295. doi: 10.1002/glia.23602
Kataoka, S., Takuma, K., Hara, Y., Maeda, Y., Ago, Y., and Matsuda, T. (2013). Autism-like behaviours with transient histone hyperacetylation in mice treated prenatally with valproic acid. Int. J. Neuropsychopharmacol. 16, 91–103. doi: 10.1017/S1461145711001714
Kim, D., Paggi, J. M., Park, C., Bennett, C., and Salzberg, S. L. (2019). Graph-based genome alignment and genotyping with HISAT2 and HISAT-genotype. Nat. Biotechnol. 37, 907–915. doi: 10.1038/s41587-019-0201-4
Kim, H. K., Lee, D. W., Kim, E., Jeong, I., Kim, S., Kim, B. J., et al. (2020). Notch signaling controls oligodendrocyte regeneration in the injured telencephalon of adult zebrafish. Exp. Neurobiol. 29, 417–424. doi: 10.5607/en20050
Kumar, L., and Futschik, M. E. (2007). Mfuzz: a software package for soft clustering of microarray data. Bioinformation 2, 5–7. doi: 10.6026/97320630002005
Kumar, S., Reynolds, K., Ji, Y., Gu, R., Rai, S., and Zhou, C. J. (2019). Impaired neurodevelopmental pathways in autism spectrum disorder: a review of signaling mechanisms and crosstalk. J. Neurodev. Disord. 11:10. doi: 10.1186/s11689-019-9268-y
Kuo, H. Y., and Liu, F. C. (2018). Molecular pathology and pharmacological treatment of autism spectrum disorder-like phenotypes using rodent models. Front. Cell Neurosci. 12:422. doi: 10.3389/fncel.2018.00422
Lafrenaye, A. D., and Fuss, B. (2010). Focal adhesion kinase can play unique and opposing roles in regulating the morphology of differentiating oligodendrocytes. J. Neurochem. 115, 269–282. doi: 10.1111/j.1471-4159.2010.06926.x
Li, C., Xie, Z., Xing, Z., Zhu, H., Zhou, W., Xie, S., et al. (2021). The notch signaling pathway regulates differentiation of ng2 cells into oligodendrocytes in demyelinating diseases. Cell Mol. Neurobiol. doi: 10.1007/s10571-021-01089-0 [Epub ahead of print].
Liberzon, A., Birger, C., Thorvaldsdottir, H., Ghandi, M., Mesirov, J. P., and Tamayo, P. (2015). The molecular signatures database (MSigDB) hallmark gene set collection. Cell Syst. 1, 417–425. doi: 10.1016/j.cels.2015.12.004
Liberzon, A., Subramanian, A., Pinchback, R., Thorvaldsdottir, H., Tamayo, P., and Mesirov, J. P. (2011). Molecular signatures database (MSigDB) 3.0. Bioinformatics 27, 1739–1740. doi: 10.1093/bioinformatics/btr260
Livak, K. J., and Schmittgen, T. D. (2001). Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method. Methods 25, 402–408. doi: 10.1006/meth.2001.1262
Lord, C., Brugha, T. S., Charman, T., Cusack, J., Dumas, G., Frazier, T., et al. (2020). Autism spectrum disorder. Nat. Rev. Dis. Primers 6:5.
Lorenzati, M., Boda, E., Parolisi, R., Bonato, M., Borsello, T., Herdegen, T., et al. (2021). c-Jun N-terminal kinase 1 (JNK1) modulates oligodendrocyte progenitor cell architecture, proliferation and myelination. Sci. Rep. 11:7264. doi: 10.1038/s41598-021-86673-6
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
Luo, T., Ou, J. N., Cao, L. F., Peng, X. Q., Li, Y. M., and Tian, Y. Q. (2020). The autism-related lncRNA MSNP1AS regulates moesin protein to influence the RhoA, Rac1, and PI3K/Akt pathways and regulate the structure and survival of neurons. Autism Res. 13, 2073–2082. doi: 10.1002/aur.2413
Mahmood, H. M., Aldhalaan, H. M., Alshammari, T. K., Alqasem, M. A., Alshammari, M. A., Albekairi, N. A., et al. (2020). The role of nicotinic receptors in the attenuation of autism-related behaviors in a murine BTBR T + tf/J Autistic model. Autism Res. 13, 1311–1334. doi: 10.1002/aur.2342
Mootha, V. K., Lindgren, C. M., Eriksson, K. F., Subramanian, A., Sihag, S., Lehar, J., et al. (2003). PGC-1alpha-responsive genes involved in oxidative phosphorylation are coordinately downregulated in human diabetes. Nat. Genet. 34, 267–273. doi: 10.1038/ng1180
Nicolini, C., Ahn, Y., Michalski, B., Rho, J. M., and Fahnestock, M. (2015). Decreased mTOR signaling pathway in human idiopathic autism and in rats exposed to valproic acid. Acta Neuropathol. Commun. 3:3. doi: 10.1186/s40478-015-0184-4
Oztan, O., Garner, J. P., Constantino, J. N., and Parker, K. J. (2020). Neonatal CSF vasopressin concentration predicts later medical record diagnoses of autism spectrum disorder. Proc. Natl. Acad. Sci. U.S.A. 117, 10609–10613. doi: 10.1073/pnas.1919050117
Oztan, O., Garner, J. P., Partap, S., Sherr, E. H., Hardan, A. Y., Farmer, C., et al. (2018). Cerebrospinal fluid vasopressin and symptom severity in children with autism. Ann. Neurol. 84, 611–615. doi: 10.1002/ana.25314
Parker, K. J., Garner, J. P., Oztan, O., Tarara, E. R., Li, J., Sclafani, V., et al. (2018). Arginine vasopressin in cerebrospinal fluid is a marker of sociality in nonhuman primates. Sci. Transl. Med. 10:eaam9100. doi: 10.1126/scitranslmed.aam9100
Parker, K. J., Oztan, O., Libove, R. A., Mohsin, N., Karhson, D. S., Sumiyoshi, R. D., et al. (2019). A randomized placebo-controlled pilot trial shows that intranasal vasopressin improves social deficits in children with autism. Sci. Transl. Med. 11:eaau7356. doi: 10.1126/scitranslmed.aau7356
Phan, B. N., Bohlen, J. F., Davis, B. A., Ye, Z., Chen, H. Y., Mayfield, B., et al. (2020). A myelin-related transcriptomic profile is shared by Pitt-Hopkins syndrome models and human autism spectrum disorder. Nat. Neurosci. 23, 375–385. doi: 10.1038/s41593-019-0578-x
Qin, L., Dai, X., and Yin, Y. (2016). Valproic acid exposure sequentially activates Wnt and mTOR pathways in rats. Mol. Cell Neurosci. 75, 27–35. doi: 10.1016/j.mcn.2016.06.004
Raam, T., and Hong, W. (2021). Organization of neural circuits underlying social behavior: a consideration of the medial amygdala. Curr. Opin. Neurobiol. 68, 124–136. doi: 10.1016/j.conb.2021.02.008
Raudvere, U., Kolberg, L., Kuzmin, I., Arak, T., Adler, P., Peterson, H., et al. (2019). G:Profiler: a web server for functional enrichment analysis and conversions of gene lists (2019 update). Nucleic Acids Res. 47, W191–W198. doi: 10.1093/nar/gkz369
Rein, B., Ma, K., and Yan, Z. (2020). A standardized social preference protocol for measuring social deficits in mouse models of autism. Nat. Protoc. 15, 3464–3477. doi: 10.1038/s41596-020-0382-9
Schneider, T., and Przewlocki, R. (2005). Behavioral alterations in rats prenatally exposed to valproic acid: animal model of autism. Neuropsychopharmacology 30, 80–89. doi: 10.1038/sj.npp.1300518
Seguin, D., Pac, S., Wang, J., Nicolson, R., Martinez-Trujillo, J., and Duerden, E. G. (2021). Amygdala subnuclei development in adolescents with autism spectrum disorder: association with social communication and repetitive behaviors. Brain Behav. 11:e2299. doi: 10.1002/brb3.2299
Servadio, M., Manduca, A., Melancia, F., Leboffe, L., Schiavi, S., Campolongo, P., et al. (2018). Impaired repair of DNA damage is associated with autistic-like traits in rats prenatally exposed to valproic acid. Eur. Neuropsychopharmacol. 28, 85–96. doi: 10.1016/j.euroneuro.2017.11.014
Sharma, A., and Mehan, S. (2021). Targeting PI3K-AKT/mTOR signaling in the prevention of autism. Neurochem. Int. 147:105067. doi: 10.1016/j.neuint.2021.105067
Shen, M. D., Li, D. D., Keown, C. L., Lee, A., Johnson, R. T., Angkustsiri, K., et al. (2016). Functional connectivity of the amygdala is disrupted in preschool-aged children with autism spectrum disorder. J. Am. Acad. Child Adolesc. Psychiatry 55, 817–824. doi: 10.1016/j.jaac.2016.05.020
Sparapani, S., Millet-Boureima, C., Oliver, J., Mu, K., Hadavi, P., Kalostian, T., et al. (2021). The biology of vasopressin. Biomedicines 9:89.
Subramanian, A., Tamayo, P., Mootha, V. K., Mukherjee, S., Ebert, B. L., Gillette, M. A., et al. (2005). Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc. Natl. Acad. Sci. U.S.A. 102, 15545–15550. doi: 10.1073/pnas.0506580102
Taleb, A., Lin, W., Xu, X., Zhang, G., Zhou, Q. G., Naveed, M., et al. (2021). Emerging mechanisms of valproic acid-induced neurotoxic events in autism and its implications for pharmacological treatment. Biomed. Pharmacother. 137:111322. doi: 10.1016/j.biopha.2021.111322
Tylee, D. S., Hess, J. L., Quinn, T. P., Barve, R., Huang, H., Zhang-James, Y., et al. (2017). Blood transcriptomic comparison of individuals with and without autism spectrum disorder: a combined-samples mega-analysis. Am. J. Med. Genet. B Neuropsychiatr. Genet. 174, 181–201. doi: 10.1002/ajmg.b.32511
Vallee, A., Vallee, J. N., and Lecarpentier, Y. (2019). PPARgamma agonists: potential treatment for autism spectrum disorder by inhibiting the canonical WNT/beta-catenin pathway. Mol. Psychiatry 24, 643–652. doi: 10.1038/s41380-018-0131-4
Vitor-Vieira, F., Vilela, F. C., and Giusti-Paiva, A. (2021). Hyperactivation of the amygdala correlates with impaired social play behavior of prepubertal male rats in a maternal immune activation model. Behav. Brain Res. 414:113503. doi: 10.1016/j.bbr.2021.113503
Voineagu, I., Wang, X., Johnston, P., Lowe, J. K., Tian, Y., Horvath, S., et al. (2011). Transcriptomic analysis of autistic brain reveals convergent molecular pathology. Nature 474, 380–384. doi: 10.1038/nature10110
Worthmuller, J., and Ruegg, C. (2020). The crosstalk between FAK and Wnt signaling pathways in cancer and its therapeutic implication. Int. J. Mol. Sci. 21:9107. doi: 10.3390/ijms21239107
Wu, J., Dai, Y. C., Lan, X. Y., Zhang, H. F., Bai, S. Z., Hu, Y., et al. (2021). Postnatal AVP treatments prevent social deficit in adolescence of valproic acid-induced rat autism model. Peptides 137:170493. doi: 10.1016/j.peptides.2021.170493
Xu, Q., Zuo, C., Liao, S., Long, Y., and Wang, Y. (2020). Abnormal development pattern of the amygdala and hippocampus from childhood to adulthood with autism. J. Clin. Neurosci. 78, 327–332. doi: 10.1016/j.jocn.2020.03.049
Xu, S., Jiang, M., Liu, X., Sun, Y., Yang, L., Yang, Q., et al. (2021). Neural circuits for social interactions: from microcircuits to input-output circuits. Front. Neural Circuits 15:768294. doi: 10.3389/fncir.2021.768294
Yu, G., Wang, L. G., Han, Y., and He, Q. Y. (2012). clusterProfiler: an R package for comparing biological themes among gene clusters. OMICS 16, 284–287. doi: 10.1089/omi.2011.0118
Zalla, T., and Sperduti, M. (2013). The amygdala and the relevance detection theory of autism: an evolutionary perspective. Front. Hum. Neurosci. 7:894. doi: 10.3389/fnhum.2013.00894
Zhang, H. F., Dai, Y. C., Wu, J., Jia, M. X., Zhang, J. S., Shou, X. J., et al. (2016). Plasma oxytocin and arginine-vasopressin levels in children with autism spectrum disorder in china: associations with symptoms. Neurosci. Bull. 32, 423–432. doi: 10.1007/s12264-016-0046-5
Zhang, J., Zhang, J. X., and Zhang, Q. L. (2016). PI3K/AKT/mTOR-mediated autophagy in the development of autism spectrum disorder. Brain Res. Bull. 125, 152–158. doi: 10.1016/j.brainresbull.2016.06.007
Zhang, W., Ma, L., Yang, M., Shao, Q., Xu, J., Lu, Z., et al. (2020). Cerebral organoid and mouse models reveal a RAB39b-PI3K-mTOR pathway-dependent dysregulation of cortical development leading to macrocephaly/autism phenotypes. Genes Dev. 34, 580–597. doi: 10.1101/gad.332494.119
Zhang, X. F., Chen, T., Yan, A., Xiao, J., Xie, Y. L., Yuan, J., et al. (2020). Poly(I:C) challenge alters brain expression of oligodendroglia-related genes of adult progeny in a mouse model of maternal immune activation. Front. Mol. Neurosci. 13:115. doi: 10.3389/fnmol.2020.00115
Zhang, Y., Sun, Y., Wang, F., Wang, Z., Peng, Y., and Li, R. (2012). Downregulating the canonical Wnt/beta-catenin signaling pathway attenuates the susceptibility to autism-like phenotypes by decreasing oxidative stress. Neurochem. Res. 37, 1409–1419. doi: 10.1007/s11064-012-0724-2
Zhang, Y., Yang, C., Yuan, G., Wang, Z., Cui, W., and Li, R. (2015). Sulindac attenuates valproic acid-induced oxidative stress levels in primary cultured cortical neurons and ameliorates repetitive/stereotypic-like movement disorders in Wistar rats prenatally exposed to valproic acid. Int. J. Mol. Med. 35, 263–270. doi: 10.3892/ijmm.2014.1996
Keywords: autism spectrum disorder, neurodevelopmental, arginine vasopressin, amygdala, oligodendrocyte
Citation: Zhou B, Zheng X, Chen Y, Yan X, Peng J, Liu Y, Zhang Y, Tang L and Wen M (2022) The Changes of Amygdala Transcriptome in Autism Rat Model After Arginine Vasopressin Treatment. Front. Neurosci. 16:838942. doi: 10.3389/fnins.2022.838942
Received: 18 December 2021; Accepted: 31 January 2022;
Published: 25 March 2022.
Edited by:
Atsushi Yoshiki, RIKEN BioResource Research Center (BRC), JapanReviewed by:
Hidenori Tabata, Aichi Human Service Center, JapanMladen-Roko Rasin, Rutgers, The State University of New Jersey, United States
Copyright © 2022 Zhou, Zheng, Chen, Yan, Peng, Liu, Zhang, Tang and Wen. 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: Lei Tang, dGxlaTE5NzRAaG90bWFpbC5jb20=; Min Wen, NjcwNDE4MDc3QHFxLmNvbQ==
†These authors have contributed equally to this work