Altered brain gene expression but not steroid biochemistry in a genetic mouse model of neurodevelopmental disorder
Molecular Autism volume 5, Article number: 21 (2014)
The 39,XY*O mouse, which lacks the orthologues of the ADHD and autism candidate genes STS (steroid sulphatase) and ASMT (acetylserotonin O-methyltransferase), exhibits behavioural phenotypes relevant to developmental disorders. The neurobiology underlying these phenotypes is unclear, although there is evidence for serotonergic abnormalities in the striatum and hippocampus.
Using microarray and quantitative gene expression analyses, and gas chromatography–mass spectrometry, we compared brain gene expression and steroid biochemistry in wildtype (40,XY) and 39,XY*O adult mice to identify non-obvious genetic and endocrine candidates for between-group differences in behaviour and neurochemistry. We also tested whether acute STS inhibition by COUMATE in wildtype (40,XY) adult male mice recapitulated any significant gene expression or biochemical findings from the genetic comparison. Data were analysed by unpaired t-test or Mann Whitney U-test depending on normality, with a single factor of KARYOTYPE.
Microarray analysis indicated seven robust gene expression differences between the two groups (Vmn2r86, Sfi1, Pisd-ps1, Tagap1, C1qc, Metap1d, Erdr1); Erdr1 and C1qc expression was significantly reduced in the 39,XY*O striatum and hippocampus, whilst the expression of Dhcr7 (encoding 7-dehydrocholesterol reductase, a modulator of serotonin system development), was only reduced in the 39,XY*O hippocampus. None of the confirmed gene expression changes could be recapitulated by COUMATE administration. We detected ten free, and two sulphated steroids in 40,XY and 39,XY*O brain; surprisingly, the concentrations of all of these were equivalent between groups.
Our data demonstrate that the mutation in 39,XY*O mice: i) directly disrupts expression of the adjacent Erdr1 gene, ii) induces a remarkably limited suite of downstream gene expression changes developmentally, with several of relevance to associated neurobehavioural phenotypes and iii) does not elicit large changes in brain steroid biochemistry. It is possible that individuals with STS/ASMT deficiency exhibit a similarly specific pattern of gene expression changes to the 39,XY*O mouse, and that these contribute towards their abnormal neurobiology. Future work may focus on whether complement pathway function, mitochondrial metabolism and cholesterol biosynthesis pathways are perturbed in such subjects.
Autism spectrum disorders (ASDs) and Attention Deficit Hyperactivity Disorder (ADHD) exhibit partially overlapping neurobehavioural symptoms, frequent comorbidity, altered monoaminergic function, and shared genetic aetiology . Both types of disorder are significantly more frequently diagnosed in males than in females, suggesting a potential role for sex-linked genetic risk variants , and both can have long-term adverse consequences for example, increased risk for alcohol dependence in ADHD  or lack of independence, close social ties and employment in ASDs .
Both cytogenetic deletions at Xp22.32 encompassing the X-linked STS gene (encoding the enzyme steroid sulphatase) and its immediate neighbours, and inactivating point mutations within STS, appear to predispose to ADHD (particularly the inattentive subtype) ; larger cytogenetic deletions encompassing STS and more distant contiguous genes (notably NLGN4X) seem to predispose to autism and related disorders . Polymorphisms within STS are associated with ADHD risk [6, 7] and cognitive function in individuals with ADHD , whilst the gene is expressed in regions of the developing brain whose structure is known to be altered in ADHD cases . Steroid sulphatase cleaves sulphate groups from a variety of steroid hormones (for example, dehydroepiandrosterone sulphate, DHEAS) thereby altering their activity and/or specificity, and subsequent developmental and physiological effects . As sulphated and non-sulphated steroid hormones can act as modulators at key neurotransmitter receptors, including N-methyl-D-aspartic acid (NMDA) and γ-aminobutyric acid type A (GABAA) receptors , lack of STS developmentally could potentially elicit important effects on neuronal organisation processes mediated by these neurotransmitters .
Inactivating mutations within the ASMT gene, located within the pseudoautosomal region of the human X chromosome and encoding the enzyme acetylserotonin O-methyltransferase that catalyses the final step in melatonin biosynthesis, have been suggested as being potentially pathogenic in a variety of psychiatric and developmental conditions, including ASDs [11–17]. Such mutations may act to reduce systemic melatonin levels, a reported feature of individuals with ASDs . Alternatively, or additionally, they could affect upstream substrate levels in the brain or blood platelets, for example, of the growth factor serotonin (5-hydroxytryptamine, 5-HT) or blood cell function ; elevated platelet serotonin levels are a consistent finding in ASD cases .
The 39,XY*O mouse lacks both the pseudoautosomal Sts and Asmt genes (and hence their expression in all tissues) as a consequence of an end-to-end fusion of the X and Y chromosomes ; as such, it has some degree of construct validity as a genetic mouse model for neurodevelopmental disorders. On an MF1 outbred albino strain background, this mouse also exhibits considerable face validity for such disorders: it is inattentive , hyperactive, emotionally hyper-reactive (showing increased indices of stress in novel or arousing environments), occasionally aggressive , and perseverative (showing persistent responding in the absence of reinforcement) [19, 22] and exhibits reduced systemic DHEA levels . Whilst melatonin levels in wildtype and 39,XY*O MF1 male mice remain to be determined, other outbred albino strains are known to produce significant quantities of the hormone .
Currently, the neurobiology of the 39,XY*O mouse is poorly defined, although we have previously shown that it exhibits altered monoaminergic chemistry, notably elevated hippocampal and striatal serotonin levels and reduced 5-HT turnover in these regions [19, 22]. Interestingly, however, the 39,XY*O mouse, in contrast to individuals with ADHD, exhibits enhanced behavioural inhibition relative to 40,XY male controls as indexed by performance on murine versions of the 5-choice Serial Reaction Time Task and the Stop Signal Reaction Time Task [20, S.T., O.A.O. and W.D., unpublished observations]. We have previously shown that acute administration of one dose of the specific steroid sulphatase inhibitor COUMATE to wildtype male mice also results in inattention  and enhanced behavioural inhibition [S.T., O.A.O. and W.D., unpublished observations], suggesting that these phenotypes in the 39,XY*O mouse are due to the ongoing activity of the enzyme. Other phenotypes in the 39,XY*O mouse (for example, hyperactivity and anxiety) cannot be recapitulated by acute inhibition of steroid sulphatase , suggesting that they may arise from the developmental effects of deficiency for the enzyme, or alternatively from neuroendocrinological abnormalities as a consequence of ASMT deficiency.
Here, we further investigated the neurobiology of the 39,XY*O mouse using two methods, with a view to identifying biological correlates of the behaviours mentioned above. First, we compared gene expression in adult 40,XY and 39,XY*O whole brain tissue by microarray to identify non-obvious genetic changes between the two groups, that is changes that could not be predicted a priori on the basis of known biology; we assayed whole brain tissue given that Sts is widely expressed throughout the mouse brain, and because disruptions to multiple brain regions were likely to underpin the 39,XY*O behavioural phenotypes. We also tested to see whether any of the significant changes were seen in the 39,XY*O hippocampus or striatum (given the known changes in 5-HT levels), or in COUMATE-treated male mouse whole brain (and hence whether they could explain the effects of acute STS deficiency on attention and behavioural inhibition). Second, we characterised brain steroids in the adult 40,XY mouse for the first time, and compared this profile to that of the 39,XY*O mouse. We anticipated that these combined approaches might reveal new pathophysiological mechanisms underlying phenotypes associated with STS/ASMT deficiency specifically and neurodevelopmental disorders more generally, and might highlight novel factors protecting against behavioural disinhibition. We report that 39,XY*O mice show a highly specific pattern of gene expression changes, with several being of relevance to developmental disorders, but surprisingly, no large changes in brain steroid levels.
All experimental research involving animals was performed in accordance with the Animal (Scientific Procedures) Act 1986 (United Kingdom) and was approved by the UK Home Office. 39,XY*O and 40,XY mice, matched for genetic background (predominantly MF1 strain), were bred and genotyped for the presence or absence of Sts at the Medical Research Council National Institute for Medical Research (MRC NIMR) as described previously ; at six weeks of age, behaviourally naïve, group-housed 40,XY (n = 10) and 39,XY*O (n = 10) mice were culled by cervical dislocation between 1200 and 1500 hours, and their brains frozen immediately at −80°C for subsequent RNA extraction. Some 40,XY and 39,XY*O mice were subsequently transferred to Cardiff University, where they were treated with antibiotics to treat a Pasteurella pneumotropica infection prior to release onto the open racks and behavioural testing; these mice were then culled by cervical dislocation between 1100 and 1300 hours (aged 10 to 12 months), and their hippocampi and striata dissected from the left hemisphere and frozen immediately at −80°C (n = 6 per group). For the pharmacological study, behaviourally naïve, six-week old, male MF1 mice (Harlan, UK) were administered 7-O-sulphamoyl-4-methylcoumarin (COUMATE) at a dose previously shown to influence behaviour  (10 mg/kg in 0.5% methylcellulose 0.9% saline vehicle, per os (po), n = 12) or vehicle alone (n = 12); 24 hours later they were culled by cervical dislocation between 1600 and 1830 hours, their brains dissected and immediately frozen at −80°C. For the brain steroid measurements, behaviourally-naïve, group-housed mice (n = 8 per group) from MRC NIMR were culled by cervical dislocation between 1100 and 1500 hours (aged five to seven months), their brains dissected and frozen immediately at −80°C.
RNA extraction and cDNA biosynthesis
High quality total RNA (absorbance ratios of 1.7 to 2.1) for the microarray study was extracted from 40,XY and 39,XY*O hemibrain tissue using RNeasy Lipid Tissue Kit (Qiagen, Manchester, UK), and DNase-treated using RNase-free DNase kit (Qiagen, Manchester, UK). For other tissue samples, high quality RNA was extracted by homogenising the tissue in TRIReagent (Sigma-Aldrich, Gillingham, UK), precipitating with isopropanol, washing in 75% ethanol, and resuspending in RNase-free distilled water; DNase treatment was then carried out using the TURBO DNA-free kit (Life Technologies, Paisley, UK). cDNA was synthesised from 1 to 2 μg RNA per sample using EcoDry Premix with random hexamers (Clontech, Mountain View, CA, USA), and diluted 25-fold.
Total RNA was supplied to Source Bioscience Ltd. (Nottingham, UK) where it underwent quality control checks prior to conversion to cRNA and hybridisation to the Mouse Gene 1.0 ST array (Affymetrix, Santa Clara, USA) according to their standard protocol. To reduce the possibility of spurious expression differences, we carried out four separate hybridisations: 40,XY1 (five pooled hemibrains), 40,XY2 (five pooled hemibrains), 39,XY*O1 (five pooled hemibrains), and 39,XY*O2 (five pooled hemibrains). Array data were subject to normalisation using the Robust Multipoint Average (RMA) procedure (background adjustment, quantile normalisation and summarisation) within the Affymetrix Expression Console prior to being loaded in the CLC Bio Genomics Workbench to examine differential expression across samples. Microarray data are available via ArrayExpress  (accession number E-MEXP-3943).
Quantitative polymerase chain (qPCR) analysis
Quantitative PCR (qPCR) analysis of gene expression was performed using a Rotorgene 6000 coupled with a CAS1200 automated set-up, and utilizing standard consumables (Qiagen, Manchester, UK). PCR reactions were performed using 5 μl cDNA mix and 200 nM custom-designed primers (Additional file 1) and Quantace SensiMix (Bioline, London, UK). qPCR data were analysed using ΔCt methods as described previously . Briefly, individual PCR reaction data were normalised to the mean of at least three ‘housekeeping gene standards’ (Hprt, Gapdh, Actb and Rn18s, whose expression was significantly correlated within the samples), giving a value known as ΔCt. The data were processed further to show how the 39,XY*O values vary with respect to those of the 40,XY (ΔΔCt) and allow simpler graphical representation after transformation (2−ΔΔC t).
Brain steroid analysis
Pooled samples (40,XY: n = 4, 39,XY*O: n = 4) were generated for each group (two whole brains per pool). Free steroids and their sulphate esters were extracted from brain tissue and fractionated prior to derivatisation with methoxyamine (MO) and trimethylsilylimidazole (TMSI) for gas chromatography-mass spectrometry (GC-MS) as previously described  except that steroid sulphates were deconjugated by incubation with Helix sulphatase rather than by solvolysis . Briefly, after drying down and redissolving in sodium acetate (0.5 M, pH 5.0), the sulphates were deconjugated by incubation with the arylsulphatase enzyme from Helix pomatia Type H1 (Sigma-Aldrich, Gillingham, UK) at 1 mg/ml overnight at 40°C, then at 55°C for two hours. The deconjugated steroids could then be recovered by loading the incubate onto a 60 mg hydrophobic-lipophilic balance (HLB) cartridge (Waters, Milford, MA, USA). After a wash with 5 ml water, the steroids were eluted from this cartridge with 4 ml ethyl acetate. These steroids liberated by deconjugation of the sulphate fraction were then dried down under nitrogen and derivatised with MO and TMSI for GC-MS alongside the steroids from the free steroid fraction of mouse brain. Increasing amounts (1.4 to 23.1 ng) of standard reference steroids were derivatised at the same time for calibration and each sample and calibration reference included a fixed amount (50 ng) of the internal standards 16-dehydropregnenolone and 6α-methyl-17-hydroxyprogesterone. Column conditions for the GC were as previously described, with elution of the MO-TMS-steroids detected by selective ion monitoring as follows: androsterone, 270 and 360; dehydroepiandrosterone, 268 and 358; epiandrosterone, 270 and 360; androstenediol, 129 and 239; 17β-oestradiol, 129, 285 and 416; testosterone, 129 and 389; allopregnanolone, 298 and 388; 16-dehydropregnenolone, 400 and 415; allopregnanediol, 269 and 346; pregnenolone, 288 and 386; epiallopregnanolone, 243 and 388; 5α-dihydroprogesterone, 288 and 343; 20β-dihydropregnenolone, 372 and 462; 20α-dihydropregnenolone, 372 and 462; progesterone, 341 and 372; 5α,20α-tetrahydroprogesterone, 289 and 303; 6α-methyl-17-hydroxyprogesterone, 443 and 474.
For the microarray data, statistical significance of the changes in mean expression levels per group was examined by the unpaired two-tailed t-test with a single overall factor of KARYOTYPE (that is 40,XY or 39,XY*O) within CLC Bio Genomics Workbench; P-values taking into account correction for False Discovery Rate were also calculated. For the qPCR analyses, samples > two standard deviations from the group mean were excluded as outliers; statistical significance of the changes in expression levels was examined using 2-∆Ct values by the unpaired t-test (normally distributed data according to Shapiro-Wilks test) or Mann-Whitney U-test (non-normally distributed data according to Shapiro-Wilks test) with a single factor of KARYOTYPE. Where variance within the two groups was unequal, Greenhouse-Geisser corrected degrees of freedom are indicated. Where the direction of effects could be predicted, that is in confirming microarray findings, a one-tailed test was used; in other cases, a two-tailed test was used. For the brain steroid hormone analyses, data were tested for normality as above, and analysed by the unpaired t-test with a single factor of KARYOTYPE. P-values ≤ 0.05 were regarded as significant, and all data are presented as mean values ± standard error of the mean.
Filtering for nominal P-values < 0.05, and >2-fold change identified just 13 genes of interest in addition to Sts, which was, as expected, found to be significantly downregulated; of these genes, eight appeared to be upregulated, and five downregulated (Table 1). Asmt was not represented on the array. Just one gene, Erdr1 (NM_133362.1) survived correction for False Discovery Rate (P = 0.044). Of particular note was the finding that two genes that lie adjacent to one another on chromosome 11 (Sfi1 and Pisd-ps1) were both upregulated, a result highly unlikely to occur by chance. A further 1,536 transcripts (out of 26,166 RefSeq transcripts represented on the array) were potentially more subtly differentially expressed (that is P < 0.05, but <2-fold expression change); of these, 47 had previously been implicated in autism and/or ADHD according to relevant online databases [28–30]: Abat, Arhgap15, Cdh8, Chrna7, Cldn5, Cntn3, Crygc, Cttnbp2, Ctnna2, Dctn5, Dhcr7, Dhrs3, Ehmt1, Elovl6, Eif4e1b, Eln, Fbxo40, Gclc, Gdi1, Gfod1, Gpd2, Grid2, Hepacam, Hes1, Hes6, Irak1, Kcnj10, Lphn3, Lrp2, Lrrn3, Mcph1, Mctp1, Mid1, Nlgn1, Npy5r, Nucb1, Park2, Pdzd4, Plcb1, Prodh, Prrt2, Prune2, Scn2a1, Slc25a12, Srsf3, Syngap1 and Ythdc2.
Confirmation of microarray results by quantitative PCR (qPCR)
Initially, we sought to confirm microarray changes > 2-fold and with a P-value < 0.05 by qPCR using individual samples (n = 10 per group). The direction of effect was consistent between the microarray and the qPCR analysis for 12 of the 13 genes in Table 1. However, the expression of only seven of the 13 genes was confirmed as being significantly different between the two experimental groups by qPCR (Figure 1A): six upregulated: Vmn2r86: t8.00 = −6.11, P < 0.001 (expression higher than 40,XY mean in 10/10 39,XY*O samples), Sfi1: t16 = −8.62, P < 0.001 (expression higher than 40,XY mean in 10/10 39,XY*O samples), Pisd-ps1: t11.84 = −1.84, P < 0.05 (expression higher than 40,XY mean in 7/10 39,XY*O samples), Tagap1: t17 = −3.32, P < 0.005 (expression higher than 40,XY mean in 8/10 39,XY*O samples), C1qc: t9.07 = −11.47, P < 0.001 (expression higher than 40,XY mean in 10/10 39,XY*O samples), Metap1d: t9.29 = −12.02, P < 0.001 (expression higher than 40,XY mean in 10/10 39,XY*O samples), and one downregulated: Erdr1: t9.00 = 5.96, P < 0.001 (expression lower than 40,XY mean in 10/10 39,XY*O samples). Prima facie, these genes do not appear to function within common known pathways, being involved in signal transduction (Vmn2r86), phosphatase activity regulation (Sfi1), complement activation in the immune response (C1qc), mitochondrial proteolysis (Metap1d) and regulation of cell proliferation and migration (Erdr1) . The expression of Fam177a showed a trend towards being more highly expressed in 39,XY*O brain than 40,XY brain (t17 = −1.51, P = 0.08). The expression of the Gm16432, Trdn, Gm12696, Fut11 and AA388235 genes did not differ significantly as a function of KARYOTYPE (t18 = −0.67, P = 0.26, t18 = 0.45, P = 0.33, t18 = −0.18, P = 0.43, t17 = 1.08, P = 0.15 and P = 0.94 respectively).
The Erdr1 gene is located adjacent to the fusion point of the X and Y chromosomes in the 39,XY*O mouse; therefore we tested whether two additional genes not represented on the microarray, which lie directly adjacent to Erdr1 (G530011O06Rik and LOC100861696)  were also differentially expressed in 39,XY*O and 40,XY brain; this did not appear to be the case: G530011O06Rik: t13.4 = −0.037, P = 0.97, LOC100861696: P = 0.33 (see Additional file 2).
Subsequently, we examined the expression of a selection of genes with known, or suspected, roles in neurodevelopmental disorders with < 2-fold change (1.04 to 1.25-fold change) but P-values < 0.05 using qPCR. Of the 11 genes examined, none were confirmed as being significantly differentially expressed between 40,XY and 39,XY*O groups (Figure 1B, Chrna7: t17 = −0.15, P = 0.44, Abat: P = 0.33, Grm1: t17 = 1.26, P = 0.12, Lrrn3: t18 = 0.77, P = 0.23, Ppp2r5c: t17 = −0.62, P = 0.28, Ercc4: t18 = 0.32, P = 0.38, Plcb1: t16 = 0.41, P = 0.34, Ctnna2: t17 = 1.71, P = 0.06, Dhcr7: t18 = 1.04, P = 0.16, Nrg3: t18 = −0.15, P = 0.45, Plekhb1 t18 = −0.27, P = 0.40).
Gene expression changes in the hippocampus and striatum of 39,XY*O mice
Of the seven genes whose expression levels differed in 40,XY and 39,XY*O whole brain tissue, six were expressed in the adult hippocampus, and all showed reduced mean levels in 39,XY*O tissue (Figure 2A). However, for only two of these (Erdr1 and C1qc), was this reduction statistically significant: Sfi1: t10 = 1.81, P = 0.10, Pisd-ps1: t10 = 1.20, P = 0.26, Tagap1: t10 = 0.78, P = 0.45, Metap1d: t10 = 1.04, P = 0.33, Erdr1: t5.00 = 3.17, P = 0.025 (expression lower than 40,XY mean in 6/6 39,XY*O samples), C1qc: P = 0.004 (expression lower than 40,XY mean in 6/6 39,XY*O samples)). We also examined hippocampal expression of Dhcr7 (encoding 7-dehydrocholesterol reductase, P < 0.05 in the microarray analysis) given that this gene had previously been implicated in neurodevelopmental disorders, and that its deficiency had been implicated in altered serotonergic function . We found that in 39,XY*O hippocampus, Dhcr7 expression was significantly downregulated relative to 40,XY expression: t5.97 = 3.47, P = 0.013 (expression lower than 40,XY mean in 6/6 39,XY*O samples) (Figure 2A).
Of the seven genes whose expression levels differed in 40,XY and 39,XY*O whole brain tissue, six were expressed in the adult striatum; we also examined Dhcr7 expression in this tissue. Both Erdr1 and C1qc genes were found to be significantly downregulated in 39,XY*O striatum: t5.19 = 3.12, P = 0.025 and t10 = 8.69, P < 0.001 respectively (expression lower than 40,XY mean in 6/6 39,XY*O samples for both genes). The expression of the remaining five genes was equivalent in 40,XY and 39,XY*O striatum (Figure 2B, Sfi1: t10 = 1.35, P = 0.21, Pisd-ps1: t10 = −0.58, P = 0.58, Tagap1: t6.59 = 1.29, P = 0.24, Metap1d: t10 = 1.04, P = 0.32, Dhcr7: t10 = 0.33, P = 0.75).
Brain gene expression changes in wildtype COUMATE-treated male mice
To determine whether the behavioural effects arising from acute steroid sulphatase inhibition could be due to altered expression of the Vmn2r86, Sfi1, Pisd-ps1, Tagap1, C1qc, or Map1d genes, we compared their expression in COUMATE and vehicle-treated MF1 male hemibrain. Erdr1 expression was not examined in COUMATE-treated brain given that its reduced expression in the 39,XY*O brain was due to a direct effect of the genetic mutation and not a downstream effect of Sts or Asmt deficiency. None of the six genes assayed were significantly differentially expressed in these tissues: Vmn2r86 (t22 = −0.99, P = 0.33), Sfi1 (t21 = 0.01, P = 0.99), Pisd-ps1 (t20 = 0.84, P = 0.41), Tagap1 (t22 = −0.82, P = 0.42), C1qc (t22 = −0.42, P = 0.68), Metap1d (t21 = 0.64, P = 0.53) (Additional file 3).
Steroid hormone analysis
Of the 15 steroids whose identity could reliably be determined from their elution and ion profile in our GC-MS analyses, ten had concentrations above the limit of detection in the free steroid fraction: (dehydroepiandrosterone (DHEA), epiandrosterone, testosterone, allopregnanolone, pregnenolone, epiallopregnanolone, 5α-dihydroprogesterone, 20α-dihydropregnenolone, progesterone and 5α,20α-tetrahydroprogesterone). The concentrations of these ten steroids did not differ significantly between 40,XY and 39,XY*O brains, although there was considerable intra-group variability, particularly in the mutant group (Table 2).
In the steroid sulphate fraction, two steroids had concentrations above the limit of detection (DHEAS and epiandrosterone sulphate); the concentrations of both of these steroids did not differ significantly in 40,XY and 39,XY*O brain (Table 2).
The 39,XY*O mouse model exhibits face and construct validity for neurodevelopmental disorders, showing behavioural endophenotypes associated with such conditions, and lacking two genes (Sts and Asmt) whose human orthologues have been implicated in ADHD and autism pathogenesis respectively. In this study, we investigated the neurobiology of this model by two methods with a view to identifying mechanisms by which loss of function of these genes might contribute towards behavioural pathology; this work is important given the current lack of availability of single gene knockout models for Sts and Asmt.
Microarray and quantitative PCR analyses comparing 40,XY and 39,XY*O brain tissue identified a surprisingly small number of robust gene expression differences between the groups. Our inability to verify differential microarray expression calls at P < 0.05 for multiple transcripts of relevance to autism and ADHD suggests that we are likely to have successfully identified all the genes that truly differ in their expression across the whole brain in 40,XY and 39,XY*O mice. However, it is possible that there are further group differences in gene expression within specific brain regions; indeed, our microarray analysis did not indicate altered expression of Htr2c (encoding the serotonin 2c receptor) which we have previously shown to be upregulated in the 39,XY*O hippocampus .
Erdr1 gene expression was significantly reduced in 39,XY*O whole brain, and in hippocampal and striatal dissections; this gene is retained in the 39,XY*O mouse and is located adjacent to the fusion point of the X and Y chromosomes . Our current results suggest that in addition to deleting of the Sts and Asmt genes, the lesion in 39,XY*O mice disrupts a genetic element that enhances Erdr1 expression. Theoretically, in concert with loss of Sts and Asmt, reduced expression of the widely-expressed, but poorly-characterised, erythroid differentiation regulator 1 protein encoded by Erdr1 could contribute towards downstream gene expression changes, abnormal striatal and hippocampal monoamine neurochemistry, and behavioural phenotypes in the 39,XY*O model; brain and behavioural investigations in mice in which the function of this gene alone is disrupted will help clarify the extent and specificity of this contribution.
Our microarray and qPCR analyses also showed that the expression of C1qc (encoding the protein complement component 1, q subcomponent, c chain) was upregulated in the brains of young, behaviourally-naïve 39,XY*O mice, but downregulated in the striata and hippocampi of older, behaviourally-trained mutant animals. These data suggest a potential basis for the neurochemical and behavioural abnormalities seen in the 39,XY*O mouse that may be sensitive to spatiotemporal or environmental regulation, and further suggest the possibility that the neurobehavioural pathology in individuals lacking functional STS and/or ASMT proteins may be due, in part, to altered C1QC levels. Previous animal and clinical studies have implicated aberrant expression of C1q family members in developmental and behavioural phenotypes. In rodents, C1q deletion results in altered synaptic elimination [33, 34], C1qc expression levels are altered in a model of developmental hippocampal pathology , and C1qc expression is associated with behavioural phenotypes (notably the consumption of ethanol relative to water) . In man, individuals with autism can exhibit elevated C1QC serum levels  and altered gastrointestinal C1q deposition [38, 39]. Whether the 39,XY*O mouse exhibits alterations in synaptic structure/function or hippocampal structure, or heightened alcohol preference remains to be investigated. It also remains to be seen whether individuals with elevated C1QC levels and autism possess genetic mutations in either STS or ASMT, and whether individuals lacking functional STS and/or ASMT genes are at increased risk of alcohol dependence.
By the same logic, the present findings further indicate that disrupted expression of Metap1d and/or Sfi1 could play a role in 39,XY*O phenotypes, and in developmental phenotypes associated with Xp22.3 mutations, but probably not through influencing striatal or hippocampal physiology. There is some evidence for an association between a linkage block at 2q31.1 containing METAP1D and autism , whilst copy number variants encompassing SFI1 have previously been identified in autism and related developmental disorders [41–45].
We also showed that the genetic mutation in 39,XY*O mice resulted in reduced hippocampal expression of the Dhcr7 gene. DHCR7 is a known modulator of serotonergic system development in mammals ; therefore, its reduced expression represents a strong candidate mechanism for abnormal serotonin levels in the 39,XY*O hippocampus and associated behavioural phenotypes [19, 22]. In man, defects in the DHCR7 enzyme underlie Smith-Lemli-Opitz syndrome (SLOS). Individuals with SLOS exhibit a range of behavioural symptoms with some overlap with autism, including: hyperactivity, aggression, insomnia, self-injurious behaviour, sensory hypersensitivity and repetitive behaviours ; interestingly, several of these behavioural abnormalities are also observed in the 39,XY*O mouse [19, 21, 22], indicating that hippocampal loss of DHCR7 function may underlie key SLOS phenotypes, and suggesting the 39,XY*O mouse as a potential novel model for aspects of the syndrome.
Genetic and functional work in mice has indicated a link between steroid sulphatase and aggressive behaviour [46, 47]. Consistent with this, 39,XY*O mice , and mice co-administered COUMATE and DHEAS , exhibit elevated levels of aggression. Previous rodent studies have demonstrated that major urinary proteins may elicit aggressive behaviour through their actions at sensory neurons expressing Vmn2r putative pheromone receptors . Of the 124 genes within the Vmn2r family, the expression of just one, Vmn2r86, was significantly altered in 39,XY*O brain; this increased expression thus represents an excellent candidate mechanism underlying aggression in these mutant mice. The fact that there is no human orthologue of Vmn2r86 may explain why individuals with Xp22.3 deletions encompassing STS and/or ASMT do not consistently show obvious aggressive tendencies.
Acute administration of COUMATE, a specific steroid sulphatase inhibitor, given at a dose known to induce behavioural changes, did not recapitulate any of the whole brain gene expression changes seen in the 39,XY*O mouse. There are two obvious possibilities why this might be the case: i) the expression changes in the 39,XY*O mouse result from abnormal developmental expression of STS and/or ii) the gene expression changes in the 39,XY*O mouse are the result of loss of function of ASMT, or reduced expression of Erdr1. A previous study found no effect of acute administration of COUMATE on the concentrations of endogenous DHEAS or DHEA in whole mouse brain, although the drug did reduce entry of systemic DHEAS into the brain . Thus, the molecular basis of COUMATE-induced behavioural changes remains obscure; it is plausible that the drug induces brain region-specific gene expression changes that we were unable to detect, and this possibility remains to be investigated.
Our current analyses provide, for the first time, a systematic profile of the steroid milieu in the mouse brain. There was substantial overlap between the free and sulphated steroids that were detectable in the adult male mouse brain (present data), and those that were most readily detectable in the adult male rat brain  consistent with a degree of cross-species homology, although, interestingly, concentrations of most compounds tended to be higher in mouse brain. We found no significant differences in the concentrations of the detectable compounds between 40,XY and 39,XY*O brains, consistent with an absence of large between-group differences in steroid brain biochemistry. This finding, taken together with our previous data showing reduced serum DHEA levels in 39,XY*O mice  suggests the possibility that where STS is absent developmentally, as in 39,XY*O mice, a compensatory mechanism is recruited to cleave sulphated steroid esters in brain, but not in peripheral tissues. Due to the difficulty of generating 39,XY*O mice and precisely genetically-matched controls, and the apparent variability in brain steroid levels in the mutant group, our study had limited power, with several steroids below the limit of detection. As such, we cannot completely exclude the possibility that there are subtle differences in levels of one or more steroids within 40,XY and 39,XY*O brain tissue.
Here, we have shown that the genetic mutation in 39,XY*O mice, in addition to deleting the Sts and Asmt genes, also results in significant downregulation of the adjacent Erdr1 gene. In the absence of Sts or Asmt single gene knockout mice, the present study gives the first clues as to possible downstream gene expression changes that might result from the loss of one (or both) of these genes; of the limited number of robust resultant gene expression changes, several may be pertinent to 39,XY*O neurochemical and behavioural phenotypes, and hence, to similar phenotypes in individuals with loss of function mutations within STS or ASMT. Importantly, it should be noted that just because gene expression differences across genotypes are large, they might not necessarily be biologically significant.
Future functional validation and pathway analysis studies in the 39,XY*O mouse throughout development, incorporating an examination of the spatiotemporal dynamics of protein changes indicated by our present study, should further elucidate the neurobiological pathways by which the 39,XY*O mutation gives rise to behavioural phenotypes analogous to those seen in neurodevelopmental disorders. Future work might examine whether the gene expression changes seen here are recapitulated in accessible tissues from individuals lacking STS and/or ASMT, and could test for abnormalities in complement pathway function, mitochondrial metabolism and cholesterol biosynthesis. Should this be the case, these physiological abnormalities could modulate either the risk of developing ADHD or autism, and/or to the clinical course of these disorders.
Our analyses indicated no large differences in brain steroid concentration between 40,XY and 39,XY*O adult mice, and hence, suggest that altered steroid biochemistry may not be a significant contributor to abnormal brain and behavioural phenotypes in this mouse model, nor to similar phenotypes in individuals with mutations in STS and/or ASMT. However, a developmental difference in brain steroid levels between Sts/Asmt-deficient and wildtype subjects may plausibly exist and contribute towards between-group adult behavioural phenotypes.
Attention Deficit Hyperactivity Disorder
autism spectrum disorders
gas chromatography-mass spectrometry
- MRC NIMR:
Medical Research Council National Institute for Medical Research
polymerase chain reaction
quantitative polymerase chain reaction
Robust Multipoint Average
Rommelse NN, Geurts HM, Franke B, Buitelaar JK, Hartman CA: A review on cognitive and brain endophenotypes that may be common in autism spectrum disorder and attention-deficit/hyperactivity disorder and facilitate the search for pleiotropic genes. Neurosci Biobehav Rev. 2011, 35: 1363-1396. 10.1016/j.neubiorev.2011.02.015.
Baron-Cohen S, Lombardo MV, Auyeung B, Ashwin E, Chakrabarti B, Knickmeyer R: Why are autism spectrum conditions more prevalent in males?. PLoS Biol. 2011, 9: e1001081-10.1371/journal.pbio.1001081.
Groenman AP, Oosterlaan J, Rommelse N, Franke B, Roeyers H, Oades RD, Sergeant JA, Buitelaar JK, Faraone SV: Substance use disorders in adolescents with attention deficit hyperactivity disorder: a 4-year follow-up study. Addiction. 2013, 108: 1503-1511. 10.1111/add.12188.
Howlin P, Goode S, Hutton J, Rutter M: Adult outcome for children with autism. J Child Psychol Psychiatry. 2004, 45: 212-229. 10.1111/j.1469-7610.2004.00215.x.
Kent L, Emerton J, Bhadravathi V, Weisblatt E, Pasco G, Willatt LR, McMahon R, Yates JR: X-linked ichthyosis (steroid sulphatase deficiency) is associated with increased risk of attention deficit hyperactivity disorder, autism and social communication deficits. J Med Genet. 2008, 45: 519-524. 10.1136/jmg.2008.057729.
Brookes KJ, Hawi Z, Kirley A, Barry E, Gill M, Kent L: Association of the steroid sulphatase (STS) gene with attention deficit hyperactivity disorder. Am J Med Genet B Neuropsychiatr Genet. 2008, 147B: 1531-1535. 10.1002/ajmg.b.30873.
Brookes KJ, Hawi Z, Park J, Scott S, Gill M, Kent L: Polymorphisms of the steroid sulphatase (STS) gene are associated with attention deficit hyperactivity disorder and influence brain tissue mRNA expression. Am J Med Genet B Neuropsychiatr Genet. 2010, 153B: 1417-1424. 10.1002/ajmg.b.31120.
Stergiakouli E, Langley K, Williams H, Walters J, Williams NM, Suren S, Giegling I, Wilkinson LS, Owen MJ, O’Donovan MC, Thapar A, Davies W: Steroid sulfatase is a potential modifier of cognition in attention deficit hyperactivity disorder. Genes Brain Behav. 2011, 10: 334-344. 10.1111/j.1601-183X.2010.00672.x.
Reed MJ, Purohit A, Woo LW, Newman SP, Potter BV: Steroid sulphatase: molecular biology, regulation, and inhibition. Endocr Rev. 2005, 26: 171-202. 10.1210/er.2004-0003.
Penn AA: Early brain wiring: activity-dependent processes. Schizophr Bull. 2001, 27: 337-347. 10.1093/oxfordjournals.schbul.a006880.
Chaste P, Clement N, Botros HG, Guillaume JL, Konyukh M, Pagan C, Scheid I, Nygren G, Anckarsater H, Rastam M, Ståhlberg O, Gillberg IC, Melke J, Delorme R, Leblond C, Toro R, Huguet G, Fauchereau F, Durand C, Boudarene L, Serrano E, Lemière N, Launay JM, Leboyer M, Jockers R, Gillberg C, Bourgeron T: Genetic variations of the melatonin pathway in patients with attention-deficit and hyperactivity disorders. J Pineal Res. 2011, 51: 394-399. 10.1111/j.1600-079X.2011.00902.x.
Etain B, Dumaine A, Bellivier F, Pagan C, Francelle L, Goubran-Botros H, Moreno S, Deshommes J, Moustafa K, Le Dudal K, Mathieu F, Henry C, Kahn JP, Launay JM, Mühleisen TW, Cichon S, Bourgeron T, Leboyer M, Jamain S: Genetic and functional abnormalities of the melatonin biosynthesis pathway in patients with bipolar disorder. Hum Mol Genet. 2012, 21: 4030-4037. 10.1093/hmg/dds227.
Galecki P, Szemraj J, Bartosz G, Bienkiewicz M, Galecka E, Florkowski A, Lewinski A, Karbownik-Lewinska M: Single-nucleotide polymorphisms and mRNA expression for melatonin synthesis rate-limiting enzyme in recurrent depressive disorder. J Pineal Res. 2010, 48: 311-317. 10.1111/j.1600-079X.2010.00754.x.
Jonsson L, Ljunggren E, Bremer A, Pedersen C, Landen M, Thuresson K, Giacobini M, Melke J: Mutation screening of melatonin-related genes in patients with autism spectrum disorders. BMC Med Genomics. 2010, 3: 10-10.1186/1755-8794-3-10.
Melke J, Goubran Botros H, Chaste P, Betancur C, Nygren G, Anckarsater H, Rastam M, Stahlberg O, Gillberg IC, Delorme R, Chabane N, Mouren-Simeoni MC, Fauchereau F, Durand CM, Chevalier F, Drouot X, Collet C, Launay JM, Leboyer M, Gillberg C, Bourgeron T: Abnormal melatonin synthesis in autism spectrum disorders. Mol Psychiatry. 2008, 13: 90-98. 10.1038/sj.mp.4002016.
Toma C, Rossi M, Sousa I, Blasi F, Bacchelli E, Alen R, Vanhala R, Monaco AP, Jarvela I, Maestrini E: Is ASMT a susceptibility gene for autism spectrum disorders? A replication study in European populations. Mol Psychiatry. 2007, 12: 977-979. 10.1038/sj.mp.4002069.
Wang L, Li J, Ruan Y, Lu T, Liu C, Jia M, Yue W, Liu J, Bourgeron T, Zhang D: Sequencing ASMT identifies rare mutations in Chinese Han patients with autism. PLoS One. 2013, 8: e53727-10.1371/journal.pone.0053727.
Veenstra-VanderWeele J, Blakely RD: Networking in autism: leveraging genetic, biomarker and model system findings in the search for new treatments. Neuropsychopharmacology. 2012, 37: 196-212. 10.1038/npp.2011.185.
Trent S, Dean R, Veit B, Cassano T, Bedse G, Ojarikre OA, Humby T, Davies W: Biological mechanisms associated with increased perseveration and hyperactivity in a genetic mouse model of neurodevelopmental disorder. Psychoneuroendocrinology. 2013, 38: 1370-1380. 10.1016/j.psyneuen.2012.12.002.
Davies W, Humby T, Kong W, Otter T, Burgoyne PS, Wilkinson LS: Converging pharmacological and genetic evidence indicates a role for steroid sulphatase in attention. Biol Psychiatry. 2009, 66: 360-367. 10.1016/j.biopsych.2009.01.001.
Trent S, Dennehy A, Richardson H, Ojarikre OA, Burgoyne PS, Humby T, Davies W: Steroid sulfatase-deficient mice exhibit endophenotypes relevant to attention deficit hyperactivity disorder. Psychoneuroendocrinology. 2012, 37: 221-229. 10.1016/j.psyneuen.2011.06.006.
Trent S, Cassano T, Bedse G, Ojarikre OA, Humby T, Davies W: Altered serotonergic function may partially account for behavioral endophenotypes in steroid sulfatase-deficient mice. Neuropsychopharmacology. 2012, 37: 1267-1274. 10.1038/npp.2011.314.
Vivien-Roels B, Malan A, Rettori MC, Delagrange P, Jeanniot JP, Pévet P: Daily variations in pineal melatonin concentrations in inbred and outbred mice. J Biol Rhythms. 1998, 13: 403-409. 10.1177/074873098129000228.
Isles AR, Davies W, Burrmann D, Burgoyne PS, Wilkinson LS: Effects on fear reactivity in XO mice are due to haploinsufficiency of a non-PAR X gene: implications for emotional function in Turner’s syndrome. Hum Mol Genet. 2004, 13: 1849-1855. 10.1093/hmg/ddh203.
Ebner MJ, Corol DI, Havlikova H, Honour JW, Fry JP: Identification of neuroactive steroids and their precursors and metabolites in adult male rat brain. Endocrinology. 2006, 147: 179-190. 10.1210/en.2005-1065.
Sujkovic E, Mileusnic R, Fry JP: Metabolism of neuroactive steroids in day-old chick brain. J Neurochem. 2009, 109: 348-359. 10.1111/j.1471-4159.2009.05965.x.
SFARI Gene: a modular database for autism research. [https://gene.sfari.org/autdb/Welcome.do]
AutDB Autism Database. [http://autism.mindspec.org/autdb/Welcome.do]
ADHDgene: a genetic database for Attention Deficit Hyperactivity Disorder. [http://adhd.psych.ac.cn]
Mouse Genome Informatics. [http://www.informatics.jax.org]
Waage-Baudet H, Lauder JM, Dehart DB, Kluckman K, Hiller S, Tint GS, Sulik KK: Abnormal serotonergic development in a mouse model for the Smith-Lemli-Opitz syndrome: implications for autism. Int J Dev Neurosci. 2003, 21: 451-459. 10.1016/j.ijdevneu.2003.09.002.
Boulanger LM: Immune proteins in brain development and synaptic plasticity. Neuron. 2009, 64: 93-109. 10.1016/j.neuron.2009.09.001.
Stevens B, Allen NJ, Vazquez LE, Howell GR, Christopherson KS, Nouri N, Micheva KD, Mehalow AK, Huberman AD, Stafford B, Sher A, Litke AM, Lambris JD, Smith SJ, John SW, Barres BA: The classical complement cascade mediates CNS synapse elimination. Cell. 2007, 131: 1164-1178. 10.1016/j.cell.2007.10.036.
Hagihara H, Takao K, Walton NM, Matsumoto M, Miyakawa T: Immature dentate gyrus: an endophenotype of neuropsychiatric disorders. Neural Plast. 2013, 2013: 318596-
Tabakoff B, Saba L, Kechris K, Hu W, Bhave SV, Finn DA, Grahame NJ, Hoffman PL: The genomic determinants of alcohol preference in mice. Mamm Genome. 2008, 19: 352-365. 10.1007/s00335-008-9115-z.
Corbett BA, Kantor AB, Schulman H, Walker WL, Lit L, Ashwood P, Rocke DM, Sharp FR: A proteomic study of serum from children with autism showing differential expression of apolipoproteins and complement proteins. Mol Psychiatry. 2007, 12: 292-306.
Torrente F, Anthony A, Heuschkel RB, Thomson MA, Ashwood P, Murch SH: Focal-enhanced gastritis in regressive autism with features distinct from Crohn’s and Helicobacter pylori gastritis. Am J Gastroenterol. 2004, 99: 598-605. 10.1111/j.1572-0241.2004.04142.x.
Torrente F, Ashwood P, Day R, Machado N, Furlano RI, Anthony A, Davies SE, Wakefield AJ, Thomson MA, Walker-Smith JA, Murch SH: Small intestinal enteropathy with epithelial IgG and complement deposition in children with regressive autism. Mol Psychiatry. 2002, 7: 375-382. 10.1038/sj.mp.4001077.
Segurado R, Conroy J, Meally E, Fitzgerald M, Gill M, Gallagher L: Confirmation of association between autism and the mitochondrial aspartate/glutamate carrier SLC25A12 gene on chromosome 2q31. Am J Psychiatry. 2005, 162: 2182-2184. 10.1176/appi.ajp.162.11.2182.
Rosenfeld JA, Ballif BC, Torchia BS, Sahoo T, Ravnan JB, Schultz R, Lamb A, Bejjani BA, Shaffer LG: Copy number variations associated with autism spectrum disorders contribute to a spectrum of neurodevelopmental disorders. Genet Med. 2010, 12: 694-702. 10.1097/GIM.0b013e3181f0c5f3.
van der Zwaag B, Franke L, Poot M, Hochstenbach R, Spierenburg HA, Vorstman JA, van Daalen E, de Jonge MV, Verbeek NE, Brilstra EH, van’ t Slot R, Ophoff RA, van Es MA, Blauw HM, Veldink JH, Buizer-Voskamp JE, Beemer FA, van den Berg LH, Wijmenga C, van Amstel HK, van Engeland H, Burbach JP, Staal WG: Gene-network analysis identifies susceptibility genes related to glycobiology in autism. PLoS One. 2009, 4: e5324-10.1371/journal.pone.0005324.
Griswold AJ, Ma D, Cukier HN, Nations LD, Schmidt MA, Chung RH, Jaworski JM, Salyakina D, Konidari I, Whitehead PL, Wright HH, Abramson RK, Williams SM, Menon R, Martin ER, Haines JL, Gilbert JR, Cuccaro ML, Pericak-Vance MA: Evaluation of copy number variations reveals novel candidate genes in autism spectrum disorder-associated pathways. Hum Mol Genet. 2012, 21: 3513-3523. 10.1093/hmg/dds164.
Kaminsky EB, Kaul V, Paschall J, Church DM, Bunke B, Kunig D, Moreno-De-Luca D, Moreno-De-Luca A, Mulle JG, Warren ST, Richard G, Compton JG, Fuller AE, Gliem TJ, Huang S, Collinson MN, Beal SJ, Ackley T, Pickering DL, Golden DM, Aston E, Whitby H, Shetty S, Rossi MR, Rudd MK, South ST, Brothman AR, Sanger WG, Iyer RK, Crolla JA, Thorland EC, Aradhya S, Ledbetter DH, Martin CL: An evidence-based approach to establish the functional and clinical significance of copy number variants in intellectual and developmental disabilities. Genet Med. 2011, 13: 777-784. 10.1097/GIM.0b013e31822c79f9.
Prasad A, Merico D, Thiruvahindrapuram B, Wei J, Lionel AC, Sato D, Rickaby J, Lu C, Szatmari P, Roberts W, Fernandez BA, Marshall CR, Hatchwell E, Eis PS, Scherer SW: A discovery resource of rare copy number variations in individuals with autism spectrum disorder. G3 (Bethesda). 2012, 2: 1665-1685. 2012.
Roubertoux PL, Carlier M, Degrelle H, Haas-Dupertuis MC, Phillips J, Moutier R: Co-segregation of intermale aggression with the pseudoautosomal region of the Y chromosome in mice. Genetics. 1994, 136: 225-230.
Le Roy I, Mortaud S, Tordjman S, Donsez-Darcel E, Carlier M, Degrelle H, Roubertoux PL: Genetic correlation between steroid sulfatase concentration and initiation of attack behavior in mice. Behav Genet. 1999, 29: 131-136. 10.1023/A:1021664607131.
Nicolas LB, Pinoteau W, Papot S, Routier S, Guillaumet G, Mortaud S: Aggressive behavior induced by the steroid sulphatase inhibitor COUMATE and by DHEAS in CBA/H mice. Brain Res. 2001, 922: 216-222. 10.1016/S0006-8993(01)03171-7.
Chamero P, Marton TF, Logan DW, Flanagan K, Cruz JR, Saghatelian A, Cravatt BF, Stowers L: Identification of protein pheromones that promote aggressive behaviour. Nature. 2007, 450: 899-902. 10.1038/nature05997.
Nicolas LB, Fry JP: The steroid sulphatase inhibitor COUMATE attenuates rather than enhances access of dehydroepiandrosterone sulphate to the brain in the mouse. Brain Res. 2007, 1174: 92-96.
The work was supported by a Medical Research Council United Kingdom (MRC UK) New Investigator Research Grant to WD (G0900636), by a Research Councils UK Fellowship to WD, and by MRC UK funding to OAO (U117532009). The funding bodies had no role in the design, collection, analysis or interpretation of data, nor in the writing of the manuscript, nor in the decision to submit the manuscript for publication. COUMATE was a kind gift to JPF from Dr Laurent Nicolas.
The authors declare that they have no competing interests.
ST, JPF and WD conceived and designed the experiments. ST, JPF and OAO performed the experiments. ST, JPF and WD analysed the data. ST, JFP, OAO and WD wrote the paper. All authors read and approved the final manuscript.
Electronic supplementary material
About this article
Cite this article
Trent, S., Fry, J.P., Ojarikre, O.A. et al. Altered brain gene expression but not steroid biochemistry in a genetic mouse model of neurodevelopmental disorder. Molecular Autism 5, 21 (2014). https://doi.org/10.1186/2040-2392-5-21
- Acetylserotonin O-methyltransferase
- Steroid sulphatase