<figure><img src="image_1.png" /></figure>

Check for updates

OPEN

# Placental endocrine function shapes cerebellar development and social behavior

Claire-Marie Vacher<sup>1✉</sup>, Helene Lacaille<sup>1</sup>, Jiaqi J. O'Reilly<sup>1</sup>, Jacquelyn Salzbank<sup>1</sup>, Dana Bakalar<sup>2</sup>, Sonia Sebaoui<sup>3</sup>, Philippe Liere<sup>4</sup>, Cheryl Clarkson-Paredes<sup>5</sup>, Toru Sasaki<sup>3</sup>, Aaron Sathyanesan<sup>3</sup>, Panagiotis Kratimenos<sup>3,10</sup>, Jacob Ellegood<sup>1</sup>,<sup>7</sup> Jason P. Lerch<sup>7,8</sup>, Yuka Imamura<sup>9</sup>, Anastas Popratiloff<sup>5,6</sup>, Kazue Hashimoto-Torii<sup>1</sup>,<sup>3,10</sup> Vittorio Gallo<sup>1</sup>,<sup>3,10</sup> Michael Schumacher<sup>4</sup> and Anna A. Penn<sup>1</sup>✉

Compromised placental function or premature loss has been linked to diverse neurodevelopmental disorders. Here we show that placenta allopregnanolone (ALLO), a progesterone-derived GABA-A receptor (GABA<sub>A</sub>R) modulator, reduction alters neurodevelopment in a sex-linked manner. A new conditional mouse model, in which the gene encoding ALLO's synthetic enzyme (akr1c14) is specifically deleted in trophoblasts, directly demonstrated that placental ALLO insufficiency led to cerebellar white matter abnormalities that correlated with autistic-like behavior only in male offspring. A single injection of ALLO or muscimol, a GABA<sub>A</sub>R agonist, during late gestation abolished these alterations. Comparison of male and female human preterm infant cerebellum also showed sex-linked myelination marker alteration, suggesting similarities between mouse placental ALLO insufficiency and human preterm brain development. This study reveals a new role for a placental hormone in shaping brain regions and behaviors in a sex-linked manner. Placental hormone replacement might offer novel therapeutic opportunities to prevent later neurobehavioral disorders.

Neuroplacentalology is a term coined to describe an emerging research area that aims to understand the influence of placental function on the developing brain. Placental dysfunction or pathology has been associated with abnormal neurodevelopment<fnref n="1" />, but the causal mechanisms remain largely unknown. Much research connecting placental compromise to fetal brain injury has focused on gas exchange or nutritional programming, neglecting the placenta's essential neuroendocrine role. Many events, including infection, malnutrition and genetic abnormalities, can disrupt placental function or, as in preterm birth, can abruptly change the hormonal environment of the developing brain. Such changes might alter early brain development or increase the susceptibility of the immature brain to damage.

ALLO, a major GABAergic neurosteroid, might be one key placental hormone in shaping fetal brain. ALLO is synthesized from progesterone and is a positive allosteric modulator of GABA<sub>A</sub>Rs that binds to specific receptor sites and prolongs the opening of the GABA-gated Cl- channel². In the adult brain, ALLO enhances GABAergic inhibition, producing sedative, anxiolytic, anesthetic and anticonvulsant effects³. A large body of evidence suggests that ALLO, through GABA<sub>A</sub>R signaling, is also a potent regulator of many neurodevelopmental processes, including neurogenesis, neuritogenesis, cell survival, synapse stabilization and myelination⁴. In fetal sheep and guinea pigs, brain ALLO is increased in response to immune challenge, hypoxia and other acute prenatal stressors. Blocking this increase pharmacologically with finasteride, a 5α-reductase inhibitor that also blocks ALLO production, leads

to increased apoptosis, excitotoxicity and impaired myelination, particularly in males⁵. These observations support an endogenous neuroprotective role for ALLO and a potential therapeutic efficacy of ALLO analogs as central nervous system drugs but do not directly connect specific placental endocrine functions to changes in neurodevelopment.

Our novel mouse model of specific placental ALLO reduction allowed direct investigation of the neurodevelopmental effect of altered placental endocrine function at a molecular, cellular and functional level. A combination of microscopy, gene and protein expression analysis, neurobehavioral assessments and pharmacological rescue in this model shows that placenta ALLO insufficiency leads to cerebellar myelination abnormalities that correlate with autistic-like behaviors in a sex-linked manner. Understanding how specific placental hormones shape normal brain development and how placental loss or dysfunction contributes to the neurological impairments in those born extremely preterm or after compromised pregnancies lays the groundwork for developing hormone replacement strategies to maintain the normal developmental milieu and protect the brain from further injury.

## Results

**Akr1c14** deletion in *Cyp19*<sup>+</sup> cells impairs placental ALLO production. During fetal development, *akr1c14*, the gene encoding the ALLO-synthesizing enzyme 3α-hydroxysteroid dehydrogenase (3α-HSD), is primarily expressed in the placenta, with up to 60 times higher expression than in mouse brain, and peaks in the

<footnote n="1">¹Department of Pediatrics, Columbia University, New York-Presbyterian Morgan Stanley Children's Hospital, New York, NY, USA. ²National Institutes of Health, Bethesda, MD, USA. ³Center for Neuroscience Research, Children's National Health System, Washington, DC, USA. ⁴U1195 INSERM, Paris-Saclay University, Le Kremlin-Bicêtre Cedex, France. ⁵The George Washington University, Nanofabrication and Imaging Center, Washington, DC, USA.</footnote>

<footnote n="6">⁶The George Washington University, SMHS, Anatomy & Cell Biology, Washington, DC, USA. ⁷Mouse Imaging Centre (MICe), Hospital for Sick Children, Toronto, ON, Canada. ⁸Wellcome Centre for Integrative Neuroimaging (WIN), Nuffield Department of Clinical Neurosciences, University of Oxford, John Radcliffe Hospital, Oxford, UK. ⁹Department of Biochemistry and Molecular Biology, Pennsylvania State University College of Medicine, Pittsburgh, PA, USA. ¹⁰The George Washington University School of Medicine and Health Sciences, Pediatrics, Washington, DC, USA.</footnote>

<footnote>✉e-mail: cv2494@cumc.columbia.edu; aap2216@cumc.columbia.edu</footnote>
a

akr1c14 gene expression

<figure><img src="image_2.png" /></figure>

b

<figure><img src="image_3.png" /></figure>

c

Cyp19-Cre:
R26R-EYFP

akr1c14-
mRNA

Overlap
(+DAPI)

Orthogonal
views

<figure><img src="image_4.png" /></figure>

d

Cyp19-Cre:R26R-EYFP + NeuN

<figure><img src="image_5.png" /></figure>

e

<figure><img src="image_6.png" /></figure>

f

akr1c14-mRNA

<figure><img src="image_7.png" /></figure>

g

akr1c14 gene expression

<figure><img src="image_8.png" /></figure>

h

Allopregnanolone

<figure><img src="image_9.png" /></figure>

**Fig. 1 |** Conditional deletion of *akr1c14* in Cyp19a-expressing trophoblasts results in reduced ALLO levels in the fetal brain. **a**, qRT-PCR for *akr1c14* in wild-type (WT) mice. Data presented as mean fold changes ± s.e.m. E12.5: *n* = 12 C and 12 pIKO; E14.5: *n* = 12 C and 12 pIKO; E15.5: *n* = 10 C and 8 pIKO; E16.5: *n* = 7 C and 7 pIKO; E17.5: *n* = 12 C and 12 pIKO; E19.5: *n* = 12 C and 11 pIKO. One-way ANOVA with Dunnett's multiple comparisons (*P* < 0.0001 compared to E12.5 value). There was no significant difference between males and females. *ywhaz*, tyrosine 3-monooxygenase/tryptophan 5-monooxygenase activation protein zeta. **b**, *akr1c14* genetic locus before and after recombination. Exons 7–9 are conditionally targeted. Cyp19a promoter drives expression of Cre in the placenta only. LoxP sites are indicated by red triangles and exons 5–12 by gray boxes. **c**, In situ hybridization for *akr1c14* in placenta from Cyp19a-Cre:R26R-EYFP mice at E17.5. Arrows show co-localizations. Scale bar, 50 µm. **d**, Sagittal section of a Cyp19a-Cre:R26R-EYFP mouse brain immunostained against YFP (green) and NeuN (red). No Cyp19a promotor activity was evidenced in the brain at P30. Scale bar, 1 mm. **e**, A 164-bp-long PCR product ascertains the presence of the recombined LoxP site in the placenta of pIKO mice. **f**, In situ hybridization for *akr1c14* in the placenta at E17.5. *akr1c14*-mRNA levels are drastically decreased in the spongiotrophoblasts (arrows) of pIKO compared to C mice. Scale bar, 30 µm. **g**, qRT-PCR for *akr1c14* normalized to *pgk1* at E17.5. Data are presented as mean ± s.e.m. Two-tailed unpaired t-test with Welch's correction (*P* < 0.05). Placenta: *n* = 7 C and 8 pIKO (*P* = 0.0247); Brain: *n* = 8 C and 9 pIKO (*P* = 0.122). *Pgk1*, phosphoglycerate kinase 1. **h**, Mass spectrometry ALLO assays at E17.5. Data are presented as mean ± s.e.m. Two-tailed unpaired t-test with Welch's correction (*P* < 0.05). Placenta: *n* = 9 C and 10 pIKO (*P* = 0.0315); Brain: *n* = 8 C and 10 pIKO (*P* = 0.0329).

second half of gestation (Fig. 1a). To assess the neurodevelopmen-
tal role of placental ALLO, we generated *akr1c14*-floxed mice and
crossed them with placenta-specific Cyp19a-Cre mice⁶ (Fig. 1b).
Recombinase activity was confirmed specifically in placenta (Fig. 1c)
but not in brain (Fig. 1d), as evidenced by crossing Cyp19-Cre
with R26R-EYFP reporter mice. *Akr1c14*-mRNA co-localized with
Cyp19a-Cre:R26R-EYFP-positive cells (arrows, Fig. 1c). Polymerase
chain reaction (PCR) confirmed recombination of LoxP sites in pla-
centa, but not brain, of *akr1c14*<sup>Cyp19a</sup> knockout (KO) mice (pIKO
for placental conditional KO) and no recombination in placenta or

brain of control (C) animals (Fig. 1e). In the pIKO mice, *akr1c14*
transcript was significantly reduced in the placenta but not brain
(Fig. 1f,g). Direct steroid measurements revealed that ALLO was
significantly reduced in the placenta and fetal brain in late gesta-
tion (Fig. 1h and Supplementary Table 1). The highest gestational
conversion rate of ALLO (based on progesterone (PROG)/ALLO
ratio) in C mice coincided with the *akr1c14* gene expression peak
at embryonic day 14.5 (E14.5) (Supplementary Table 1 and Fig. 1a).
Thus, *akr1c14* placental gene removal results in greatest placental
ALLO insufficiency in the second half of mouse gestation (E14.5).
The positive placenta-brain correlation of ALLO levels is consistent with placental provision of ALLO to the fetal brain (Extended Data Fig. 1a). Placental ALLO levels were not altered by sex (Extended Data Fig. 1a). ALLO precursors and other neuroactive steroids in the synthesis pathway were not altered (Supplementary Table 1 and Extended Data Fig. 1b–f), with the exception of a transient E14.5 reduction of placental allotetrahydrodeoxycorticosterone (ALLOTHDOC), which is also converted by placental 3α-HSD (Supplementary Table 1). Placental ALLOTHDOC concentrations are negligible when compared to ALLO levels, suggesting that major effects in the plKO are primarily due to placental ALLO insufficiency (Supplementary Table 1). Additionally, the equivalent and specific decrease in both brain and placenta of plKO ALLO levels supports the specificity of this model (Fig. 1h and Extended Data Fig. 1b–f).

**Placental ALLO loss alters postnatal cerebellar myelination.**
Unbiased RNA sequencing (RNA-seq) analysis of cerebral cortex, hippocampus and cerebellum of C and plKO mice at postnatal day 30 (P30) in both sexes was performed to define the most prominent postnatal molecular alterations associated with placental ALLO insufficiency. Threshold criteria used to identify significant differentially expressed genes (DEGs) were validated by RT-PCR and western blots (Extended Data Fig. 2). Male cerebellum was the brain structure affected most by placental ALLO deficiency as measured by the number of DEGs, which was 3–5 times higher than in any other brain region (Extended Data Fig. 3a,b and Supplementary Table 2, sheets 1 and 2). Cerebellar DEGs were equally distributed in up-regulated and down-regulated gene categories (Extended Data Fig. 3a–d). Ingenuity Pathway Analysis (IPA, Qiagen) identified white matter (WM)-associated pathways in the top categories of altered genes in both sexes (Extended Data Fig. 3e,f). Comparison of cerebellar DEGs to previous published oligodendrocyte (OL) and myelin transcriptomes⁷,⁸ revealed an overlap of 204 genes and 68 genes in males and females, respectively (Extended Data Fig. 3g). Striking sexual dimorphism was seen: male and female DEGs were qualitatively different, and most OL-/myelin-related DEGs were up-regulated in plKO males but down-regulated in plKO females (Extended Data Fig. 3h,i).

Investigation of the cerebellar myelin proteins at P30 revealed major sex-linked changes: in plKO males, myelin basic protein (MBP), myelin-associated glycoprotein (MAG) and myelin oligodendrocyte glycoprotein (MOG) levels were all significantly increased (Fig. 2a,b), whereas, in plKO females, MBP cerebellar content was reduced (Fig. 2c,d)—changes confirmed by immunohistochemistry (Fig. 2e–h) and RT-PCR (Extended Data Fig. 2h).

Similar sex-dependent MBP changes were previously described in the cerebellum of guinea pigs delivered preterm⁹. In addition, a migration downshift of MAG on western blots was observed in plKO males but not females (Fig. 2a,c), suggesting increased myelin maturation in males¹⁰. Further structural characterization of the cerebellar WM within lobule VI–VII by scanning electron microscopy (Fig. 2i) revealed thicker myelin (g-ratio measurements; Fig. 2j,k,l,o,s) and fewer unmyelinated axons in male plKOs (Fig. 2p). This increased myelin thickness was not associated with significant changes in axonal inner calibers (Fig. 2t). In contrast, cerebellar axons in plKO females had significantly thinner myelin sheath, and more axons were unmyelinated (Fig. 2m,n,q,r,u), although axonal inner diameters were unchanged (Fig. 2v) compared to C littermates.

A global comparison of regional, sex-linked and genotype-specific brain anatomy was performed by magnetic resonance imaging (MRI) at P30. Most notably, the mean fractional anisotropy (FA) of the cerebellar fiber tracts, measured by diffusion tensor imaging (DTI), was significantly increased in the plKO males, but not females, at P30 (Extended Data Fig. 4a; region effect: *F*<sub>25,1014</sub> = 57.36, *P* < 0.0001; genotype effect: *F*<sub>1,1014</sub> = 16.99, *P* < 0.0001; sex × genotype interaction: *F*<sub>1,1014</sub> = 17.92, *P* < 0.0001; three-way analysis of variance (ANOVA)). In contrast, mean FAs in other WM tracts outside of the cerebellum were not significantly different in plKO mice compared to their C littermates (Extended Data Fig. 4b; region effect: *F*<sub>8,351</sub> = 80.28, *P* < 0.0001; genotype effect: *F*<sub>1,351</sub> = 0.1618, *P* = 0.6877; sex × genotype interaction: *F*<sub>1,351</sub> = 2.589, *P* = 0.1021; three-way ANOVA). The male-specific increase in cerebellar FA is consistent with the enhanced degree of myelination¹¹. OL-/myelin-related DEGs, including *mbp*, were mostly down-regulated in both male and female cerebral cortex in plKO mice (Extended Data Fig. 5). Both the DTI results and RNA-seq transcriptome patterns suggest overall sex-linked and regionally specific myelin abnormalities after placental ALLO loss.

**Postnatal OL differentiation rate diverges in plKO males and females.**
Cerebellar myelination begins in late mouse gestation and is primarily a postnatal process¹²⁻¹⁴, so placental endocrine disruption leading to myelination anomalies was unanticipated. The adjustment of myelin sheath thickness notably relies on the axon diameter, whose distribution was unchanged in plKO cerebellum (Fig. 2t,v) and the local density of myelin-competent OLs that are primarily produced before birth as OL progenitor cells (OPCs)¹⁴, at a time when ALLO is highest. ALLO is a potent positive allosteric modulator of GABA<sub>A</sub>R². GABA<sub>A</sub>R signaling is associated with

**Fig. 2 | Placental ALLO insufficiency results in sex-linked cerebellar WM abnormalities at P30.** *a,b,* Western blot analysis of myelin-related proteins in the male cerebellum. Data are presented as mean fold changes ± s.e.m. Multiple unpaired t-tests with Holm-Sidak multiple comparison test (\* *P* < 0.05; \*\*\* *P* < 0.005). MBP: *n* = 12 C and 28 plKO (*P* = 0.038); MAG: *n* = 6 C and 6 plKO (*P* = 0.004); MOG: *n* = 6 C and 6 plKO (*P* < 0.0001). *c,d,* Western blot analysis of myelin-related proteins in the female cerebellum. Data are presented as mean fold changes ± s.e.m. Multiple unpaired t-tests with Holm-Sidak multiple comparison test (\* *P* < 0.05). MBP: *n* = 5 C and 4 plKO (*P* = 0.031); MAG: *n* = 5 C and 4 plKO (*P* = 0.18); MOG: *n* = 5 C and 4 plKO (*P* = 0.54). GAPDH, glyceraldehyde 3-phosphate dehydrogenase. **e–h,** Immunofluorescent staining of MBP in cerebellar lobule VII in males (**e** and **f**) and females (**g** and **h**). Scale bar, 150 µm. **i,** Scanning electron microscopy acquisition of a whole cerebellum ultrathin section showing the different cerebellar lobules (I–X) in a control mouse (×650 low magnification). The region of interest (inter-lobule VI-VII) for high-magnification acquisitions and g-ratio quantifications is indicated by the yellow rectangle. Scale bar, 300 µm. **j,** Schematic representation of an axon and its myelin sheath illustrating g-ratio. **r,** axon inner diameter; **R,** axon outer diameter. **k–n,** Representative scanning electron microscopy acquisitions (×20,000) of myelinated axons in inter-lobule VI-VII WM in males (**k** and **l**) and females (**m** and **n**). Scale bar, 300 nm. **o,** Scatter plot of g-ratios of >400 individual axons in males; data are from one representative control and one representative plKO mouse. Fitted lines are linear regressions. **p,** Percent of unmyelinated axons. Data are presented as mean ± s.e.m. Two-tailed unpaired Student's t-test with Welch's correction (\* *P* < 0.05). *n* = 3 C and 3 plKO (*P* = 0.0275). **q,** Scatter plot of g-ratios of >400 individual axons in females; data are from one representative control and one representative plKO mouse. Fitted lines are linear regressions. **r,** Percent of unmyelinated axons. Data are presented as mean ± s.e.m. Two-tailed unpaired Student's t-test with Welch's correction (\* *P* < 0.05). *n* = 3 C and 3 plKO (*P* = 0.0194). **s,** g-ratios by axon caliber category in males. Data are presented as mean ± s.e.m. Two-way ANOVA with Sidak's multiple comparisons test (\* *P* < 0.05). *n* = 3 C and 3 plKO (<1 µm: *P* = 0.0299; >1 µm: *P* = 0.0165). **t,** Axon caliber distribution in males (*n* = 3 C and 3 plKO). **u,** g-ratios by axon caliber category in females. Data are presented as mean ± s.e.m. Two-way ANOVA with Sidak's multiple comparisons test (\* *P* < 0.05). *n* = 3 C and 3 plKO (<1 µm: *P* = 0.0162; >1 µm: *P* = 0.1185). **v,** Axon caliber distribution in females (*n* = 3 C and 3 plKO).
blockade of OPC proliferation¹⁵,¹⁶, so we tested the hypothesis that prenatal ALLO insufficiency enhances OLs production.

We found that the density of OL lineage cells (Olig2+), and particularly mature OLs (CC1+), were significantly increased in the male plKO cerebellum at P30 compared to controls (Extended Data Fig. 6a,b). Increased OPC proliferation during fetal life is suggested by the results of co-labelling with BrdU at E15.5 and Olig2 expression at P30 (Extended Data Fig. 6c). Additionally, within the OL lineage,

an acceleration of the maturation from PDGFRα+ OPCs to mature CC1+ OLs was observed in the plKO mice at P15 (Extended Data Fig. 6d-i). Interestingly, a sex-divergent progression of the OL lineage was observed at P15. At this age, whereas the density of Olig2+ cells was increased in plKO females in the same range as in plKO males (Extended Data Fig. 6d,j), the OL maturation (given by the CC1+/Olig2+ cell number ratio) was reduced (Extended Data Fig. 6h,k). Overall, these results are consistent with the ultrastructural analysis

<figure><img src="image_10.png" /></figure>

<figure><img src="image_11.png" /></figure>

<figure><img src="image_12.png" /></figure>

<figure><img src="image_13.png" /></figure>

<figure><img src="image_14.png" /></figure>

<figure><img src="image_15.png" /></figure>

<figure><img src="image_16.png" /></figure>

<figure><img src="image_17.png" /></figure>

<figure><img src="image_18.png" /></figure>

<figure><img src="image_19.png" /></figure>

<figure><img src="image_20.png" /></figure>

<figure><img src="image_21.png" /></figure>

$$g = \frac{r}{R}$$

<figure><img src="image_22.png" /></figure>

<figure><img src="image_23.png" /></figure>

MBP

SEM

<figure><img src="image_24.png" /></figure>

<figure><img src="image_25.png" /></figure>

<figure><img src="image_26.png" /></figure>

<figure><img src="image_27.png" /></figure>

<figure><img src="image_28.png" /></figure>

<figure><img src="image_29.png" /></figure>

<figure><img src="image_30.png" /></figure>

<figure><img src="image_31.png" /></figure>
<figure><img src="image_32.png" /></figure>

<figure><img src="image_33.png" /></figure>

<figure><img src="image_34.png" /></figure>

<figure><img src="image_35.png" /></figure>

<figure><img src="image_36.png" /></figure>

<figure><img src="image_37.png" /></figure>

<figure><img src="image_38.png" /></figure>

<figure><img src="image_39.png" /></figure>

<figure><img src="image_40.png" /></figure>

<figure><img src="image_41.png" /></figure>

<figure><img src="image_42.png" /></figure>

**Fig. 3** | Male plKO mice exhibit ASD-like behavior. **a**, Three-chamber sociability test at P30. Data are presented as mean ± s.e.m. Two-way ANOVA with Sidak's multiple comparison test (***P < 0.0001). Males: *n* = 20 C and 28 plKO (*P* = 0.0001); Females: *n* = 16 C and 10 plKO (*P* > 0.9999). **b**, Spontaneous gnawing time over 15 min at P30. Data are presented as mean ± s.e.m. Two-way ANOVA with Sidak's multiple comparison test (*P* < 0.01). Males: *n* = 20 C and 28 plKO (*P* = 0.0094); Females: *n* = 15 C and 11 plKO (*P* = 0.3982). **c**, Spontaneous digging time over 15 min at P30. Data are presented as mean ± s.e.m. Two-way ANOVA with Sidak's multiple comparison test (*P* < 0.05). Males: *n* = 20 C and 28 plKO (*P* = 0.0148); Females: *n* = 15 C and 11 plKO (*P* = 0.4781). **d**, z-standardized single significant behavior readouts that are integrated into the autism composite score in males at P30. Data are presented as mean ± s.e.m. Two-way ANOVA with Sidak's multiple comparisons test (*P* < 0.05; ***P* < 0.0005). *n* = 20 C and 28 plKO. Sociability: *P* = 0.0002; Gnawing: *P* = 0.0239; Digging: *P* = 0.0301. **e**, z-standardized single significant behavior readouts that are integrated into the autism composite score in females at P30. Data are presented as mean ± s.e.m. Two-way ANOVA with Sidak's multiple comparisons test. *n* = 15 C and 10 plKO. Sociability: *P* = 0.9802; Gnawing: *P* = 0.1419; Digging: *P* = 0.2979. **f**, ASD composite severity score at P30. Data are presented as mean ± s.e.m. Two-way ANOVA with Sidak's multiple comparison test (***P* < 0.0001). Males: *n* = 20 C and 28 plKO (*P* < 0.0001); Females: *n* = 16 C and 10 plKO (*P* > 0.9997). **g.h**, Relative frequency distribution of autism composite score bins at P30. Males: *n* = 20 C and 28 plKO; Females: *n* = 16 C and 10 plKO. **i**, Positive correlation between cerebellum MBP levels determined by western blot (normalized with GAPDH) and autism severity score in males at P30 (*n* = 12 C and 28 plKO). Linear regression. Deviation from zero: *P* = 0.0009. **j.k**, USVs at P4, P6, P8 and P11. Individual pups, males (**j**) and females (**k**), were separated from the dam and littermates, and their calls were recorded for 3 min. Data are presented as mean ± s.e.m. Two-way repeated measures ANOVA followed by Sidak's multiple comparisons test (*P* < 0.05). Males: *n* = 13 C and 15 plKO (P4: *P* = 0.0378; P6: *P* = 0.5714; P8: *P* = 0.9299; P11: *P* = 0.9609); Females: *n* = 7 C and 7 plKO (P4: *P* = 0.4587; P6: *P* = 0.6051; P8: *P* > 0.9999; P11: *P* = 0.7135). GAPDH, glyceraldehyde 3-phosphate dehydrogenase.

of myelin and indicate short-term, sex-independent (embryonic OPC proliferation) and long-term, sex-dependent (OL maturation) modifications of cerebellar OL lineage progression and OL maturation in

mice lacking placental ALLO. No volumetric change in the different cerebellar layers (Extended Data Fig. 7a-i) nor Purkinje cell linear density (Extended Data Fig. 7j) were observed in the plKO mice.
<figure><img src="image_43.png" /></figure>

<figure><img src="image_44.png" /></figure>

<figure><img src="image_45.png" /><figcaption>**Fig. 4** | ALLO or muscimol administration during late gestation rescues MBP expression and abnormal behaviors in plKO males at P30. a, Western blot determination of MBP contents in the cerebellum at P30. Normalized data to GAPDH contents are presented as mean fold changes ± s.e.m. Dams received an intraperitoneal injection of ALLO (10 mg kg⁻¹) or muscimol (Musci; 1 mg kg⁻¹) at E15.5. One-way ANOVA with Tukey's multiple comparisons test (\*P < 0.05; \*\*P < 0.01; \*\*\*P < 0.005). n = 8 C, 11 plKO+Veh, 6 plKO+ALLO and 8 plKO+muscimol. Comparisons: C versus plKO+Veh: P = 0.009; C versus plKO+ALLO: P = 0.988; C versus plKO+muscimol: P = 0.607; plKO+Veh versus plKO+ALLO: P = 0.0058; plKO+Veh versus plKO+muscimol: P = 0.027; plKO+ALLO versus plKO+muscimol: P = 0.8436. b,c, ASD composite score and frequency distribution of autism severity score bins at P30. plKO and plKO+Veh mice were combined because there was no significant difference between both groups (named plKO(+Veh)). One-way ANOVA with Tukey's multiple comparisons test (\*\*\*P < 0.005; \*\*\*\*P < 0.0005). n = 20 C, 3 6plKO(+Veh), 17 plKO+ALLO and 9 plKO+muscimol. Comparisons: C versus plKO(+Veh): P < 0.0001; C versus plKO+ALLO: P = 0.6226; C versus plKO+muscimol: P = 0.4933; plKO(+Veh) versus plKO+ALLO: P = 0.0002; plKO(+Veh) versus plKO+muscimol: P = 0.0187; plKO+ALLO versus plKO+muscimol: P = 0.9765. GAPDH, glyceraldehyde 3-phosphate dehydrogenase.</figcaption></figure>

**Placental ALLO insufficiency leads to autism spectrum disorder-like behaviors in males.** Cerebellar circuits provide precise spatiotemporal control of sensorimotor behavior, and cerebellar contribution to cognitive functions has more recently been recognized. Thus, developmental disruption of cerebellar circuits, either through neuronal or glial alterations, might modify motor and cognitive performance. Cerebellar-dependent tests were first selected to assess locomotor and motor learning behaviors. Behavior testing on the Erasmus ladder, an advanced, fully automated system for comprehensive real-time monitoring of cerebellar-mediated motor function, was performed at P30, the earliest validated age¹⁷, to align with the anatomical analysis. At baseline, male plKO mice exhibited a mild alteration of their walking pattern, with higher percentage of long steps (Extended Data Fig. 8a,b). The frequency of missteps, often increased in ataxia models, was unchanged (Extended Data Fig. 8g). The locomotion speed was similar in male plKO and C on the Erasmus ladder, consistent with open field testing results (Extended Data Fig. 8c,d). An associative cerebellar learning task pairing a tone with a ladder obstacle revealed no cerebellar motor learning deficit in male plKO mice (Extended Data Fig. 8e,f). Female plKOs exhibited no locomotor alterations nor learning deficits (Extended Data Fig. 8i–n). Overall, our data suggest that plKO mice do not display major locomotor, ataxia-like or associative motor learning alterations. The only abnormal locomotion feature identified in plKO mice on the Erasmus ladder was an increased frequency of long steps by males. More frequent long steps or longer stride was previously observed in genetic mouse models with impeded intrinsic or synaptic Purkinje cell plasticity¹⁸,¹⁹. However, contrasting outcomes on the Erasmus ladder were also described in mice with related physiological impairments²⁰. Future studies focusing on Purkinje cell physiology might elucidate the precise mechanisms underlying the subtle walking pattern abnormalities seen in plKO males. Interestingly, longer stride was also reported as an autism spectrum disorder (ASD)-associated feature in human and mouse¹⁹,²¹.

To further evaluate balance, motor coordination and learning, mice were tested on an accelerating rotarod. Trained male, but not female, plKO mice displayed longer latency to fall and increased terminal speed (Extended Data Fig. 9a,b). The learning rate of individual mice, given by a linear regression analysis, was also enhanced in male plKO mice (Extended Data Fig. 9c,d). Greater learning in plKO males on the accelerating rotarod might be linked to their increased cerebellar FA, as previously described²². This gain of function, although rare, has been described in ASD mouse models, in which it has been attributed to increased stereotyped motor routine²³. Although the underlying molecular mechanisms and circuits might differ between models, similar repetitive/stereotypical behavior could play a role in the enhanced rotarod learning displayed by plKO males.

Beside its role in motor learning and coordination, the cerebellum contributes to cognitive processing and emotional control. Cerebellar impairments are often associated with ASD symptoms in humans and genetic mouse models²⁴,²⁵. Specifically, cerebellar WM abnormalities, including transient hyperplasia or dysregulation of myelin-related genes, were reported in patients with ASD²⁶,²⁷. We next asked whether mice lacking placental ALLO show autistic-like behavioral features by evaluating adult sociability and motor stereotypies, as well as isolation calling response at young ages. The three-chamber social behavior test revealed significant social interaction deficits in plKO males but not females (Fig. 3a). Stereotyped motor behavior was then investigated in spontaneously behaving mice. Male, but not female, plKO mice exhibited specific motor stereotypies, particularly increased gnawing (Fig. 3b). In addition, male plKO mice spent less time digging spontaneously or during a marble burying test (Fig. 3c and Extended Data Fig. 9e). Difficulties in social interactions and motor stereotypies are two ASD hallmarks,
<figure><img src="image_46.png" /><figcaption>**Fig. 5** | Myelin proteins are dysregulated in a sex-linked manner in the cerebellar vermis of preterm infants. **a**, Schematic of the age terminology during the perinatal period used for the human study. **b**, Distribution of donors' gestational age at birth. **c**, Distribution of donors' corrected age at death. Graphs in **b** and **c** show box and whisker plots (including minima, maxima and median values) with single values. **d,e**, Western blot analysis of myelin-related proteins in the cerebellar vermis of term and preterm male infants at an average corrected age of 6 weeks. Data are presented as mean fold changes ± s.e.m. (n = 6 T and 4 PT). Multiple unpaired t-tests with Holm-Sidak multiple comparison test (*P* < 0.05). MBP: P = 0.014; MAG: P = 0.044; MOG: P = 0.024. **f,g**, Western blot analysis of myelin-related proteins in the cerebellar vermis of term and preterm female infants at an average corrected age of 6 weeks. Data are presented as mean fold changes ± s.e.m. Multiple unpaired t-tests with Holm-Sidak multiple comparison test (*P* < 0.05). MBP: n = 3 T and 6 PT (P = 0.015); MAG: n = 3 T and 6 PT (P = 0.169); MOG: n = 3 T and 5 PT (P = 0.179). **h-j**, qRT-PCR for myelin- and OL-related genes normalized to *cyc1* in male infants. Data are presented as mean fold changes ± s.e.m. Two-tailed unpaired Student's t-test with Welch's correction (*P* < 0.05, ****P < 0.0001). **mbp**: n = 7 T and 5 PT (P = 0.046); **pdgfra**: n = 9 T and 5 PT (P = 0.0482); **cspg4**: n = 9 T and 5 PT (P < 0.0001). **k,l**, qRT-PCR for myelin- and OL-related genes normalized to *cyc1* in female infants. Data are presented as mean fold changes ± s.e.m. Two-tailed unpaired Student's t-test with Welch's correction. **mbp**: n = 4 T and 6 PT (P = 0.9923); **pdgfra**: n = 4 T and 6 PT (P = 0.9122). **Cspg4**, chondroitin sulfate proteoglycan 4 gene; **Cyc1**, cytochrome C1 gene; **GAPDH**, glyceraldehyde 3-phosphate dehydrogenase; **Pdgfra**, platelet-derived growth factor receptor alpha gene. **PT**, preterm; **T**, term.</figcaption></figure>

especially in males²⁸. In addition, both increased and decreased digging and marble burying were noted as secondary ASD-like features related to emotional state and anxiety in several genetic ASD mouse models¹⁹,²¹,²⁹⁻³¹.

Behavioral sex-linked differences between C and pIKO mice can be detected even earlier in postnatal development. The ultrasonic vocalizations (USVs) emitted by the pups when separated from the dam (referred to as isolation calling response) were evaluated as a measure of the degree of aversive affective state and early communicative behavior. We found that pIKO males produced fewer USVs than their littermate controls when separated from the dam (Fig. 3j). By contrast, no significant change was evidenced in the

pIKO females (Fig. 3k). Deficient isolation calling response has been detected in different mouse models of ASD³², indicating that insufficient provision of placental ALLO to the fetal brain can lead to a wide array of autistic-like behaviors in the male progeny.

**Cerebellar MBP correlates with ASD-like symptom severity.** We calculated an autism severity composite index based on previously established scoring systems³³. Behaviors were z-standardized with higher values representing higher symptom severity. Individual behavioral scores were significantly higher in pIKO males compared to C but not in females (Fig. 3d,e). Thus, the combined z-scores yielded an elevated autism severity score in males (Fig. 3f). The
relative frequency distribution of individual autism severity scores followed Gaussian distribution curves that differed in plKO males, which shifted to the right along the *x* axis (Fig. 3g,h). Interestingly, some of the plKO mice exhibited scores similar to control mice, suggesting a spectrum of behavioral changes in male plKOs. This variability might arise from natural variations in placental progesterone and ALLO levels, variable deletion of *akr1c14*, mixed genetic background of the mice or subtle changes in hormone exposure due to in utero position of the fetus.

A significant positive correlation between autism symptom severity scores and cerebellar MBP levels was seen in plKO males (Fig. 3i). This finding is of particular interest knowing that transient WM hyperplasia, myelin thickening and/or over-expression of OL lineage and myelin markers such as *olig2*, *mbp* and *mag* have been documented in the cerebellum of patients with ASD<sup>26,27</sup> and several genetic ASD-like mouse models<sup>34,35</sup>. Together with the enhanced FA in the cerebellum-related tracts described above (Extended Data Fig. 4a), this correlation could reflect a short-range hyperconnectivity and/or a transient WM hyperplasia as seen in patients with autism<sup>26,36</sup>. The abnormal maturation speed of the OL lineage in plKO mice suggests dynamic, rather than uniform, changes of WM development, similar to the pattern seen across the lifespan of patients with ASD<sup>26</sup>. Overall, our placental mouse model recapitulates several characteristics of human studies, including the male greater vulnerability to perinatal brain injuries and ASD<sup>19,23,27,37</sup>.

**ALLO or GABA<sub>A</sub> agonist rescues cerebellar and behavioral impairments.** The therapeutic potential of ALLO to rescue cerebellar WM and behavioral abnormalities was then tested by administering ALLO (10 mg kg<sup>-1</sup>) or vehicle (sesame oil) to dams carrying mixed litters at E15.5, the period of maximal difference of ALLO levels in the plKO mice compared to C. A single injection of ALLO during gestation prevented the cerebellar MBP up-regulation (Fig. 4a), normalized OL density in the cerebellar WM (Extended Data Fig. 10a) and reduced the autism severity score (Fig. 4b,c). To assess the specific contribution of GABA<sub>A</sub>R signaling to this rescue, muscimol (1 mg kg<sup>-1</sup>), a selective GABA<sub>A</sub> agonist, was administered instead of ALLO. Muscimol injection resulted in similar molecular and behavioral rescue in plKO males (Fig. 4a-c), although the contribution of other minor, GABA<sub>A</sub>-independent ALLO actions cannot be completely excluded<sup>38-40</sup>. Interestingly, a fraction of C mice exposed to ALLO treatment exhibited decreased social preference and subsequent higher autism score (Extended Data Fig. 10b,c), suggesting that excess ALLO exposure during fetal life might be detrimental. Interestingly, increased PROG exposure, along with other steroids, was previously linked to ASD risk in boys<sup>41</sup>. Our results support the potential therapeutic utility of ALLO administration during gestation if ALLO or its precursors are determined to be low (as might occur with chronic placental insufficiency) but, additionally, suggest the need to maintain fetal ALLO exposure within an appropriate physiological window.

**Myelin genes are similarly altered in preterm infants and plKO mice.** To further evaluate the translational potential of our findings, we asked whether preterm birth, a condition characterized by premature loss of placenta and a risk factor for ASD<sup>42</sup>, was associated with cerebellar myelin alterations similar to those seen in plKO mice. We examined 24 postmortem cerebellar vermises from term (&gt;37 weeks of gestation) and preterm (&lt;34 weeks of gestation) infants who died at a mean corrected age of 6 weeks (Fig. 5a-c and Supplementary Table 3). Compelling similarities between our mouse model and human pathology were found in cerebellar MBP expression in preterm males and females (Fig. 5d-h,k) compared to corrected age-matched specimens. Furthermore, the migration downshift of MAG on western blot gels was observed in human preterm male cerebellum, similar to our findings in plKO males (Figs. 2a and 5d), again suggesting accelerated maturation of myelin<sup>10</sup>. The higher levels of *pdgfra*- and *cspg4*-mRNA in the cerebellum of preterm males also suggest increased OPC proliferation (Fig. 5i,j,l), as temporarily seen in the plKO males (Extended Data Fig. 6e). This striking correlation between mice exposed to ALLO insufficiency and infants born preterm suggests that altered ALLO exposure might contribute to abnormal preterm brain development.

Discussion

Our study empirically defines a critical role for a specific placental hormone that can alter fetal brain development in late gestation with postnatal developmental consequences. Many neuroactive hormones are produced in the placenta in late gestation and might affect specific stages of neurogenesis and gliogenesis, neuronal migration and circuit formation. We focused on ALLO because this placental steroid hormone normally peaks in the second half of gestation, when placental insufficiency becomes evident and when preterm birth occurs. Interestingly, in severe pre-eclampsia, a pregnancy condition associated with both preterm birth and ASD, plasma ALLO levels in pregnant women are reduced, as compared to control group<sup>43</sup>, strengthening the potential link among ALLO, prematurity and ASD risk.

Many brain regions have been linked to autism, but structural and functional cerebellar abnormalities are among the most consistent findings<sup>44-46</sup>. Early cerebellar injury or dysfunction carries a high relative risk (up to 36×) for ASD, exceeded only by genetic factors<sup>45</sup>. Both neuronal dysfunction, especially of cerebellar Purkinje cells, and WM changes have been linked to ASD phenotypes. Cerebellar transcriptome analysis highlighted WM-associated pathways as a top category of altered genes. Anatomical myelination changes, consistent with these transcriptomic findings, suggest that oligodendrogenesis and myelination significantly contribute to the observed outcomes. In our mouse model, Purkinje cell number appears normal (Extended Data Fig. 7j), but the contribution of Purkinje cell physiology to the myelin and behavioral phenotypes will be assessed in future studies. Additionally, the local specificity of changes will need to be further assessed. Motor and non-motor functions of the cerebellum are controlled by specific subdivisions—that is, lobules and modules—that are characterized by specific inputs and outputs, Purkinje cell firing frequency, degree of myelination and neurochemical signatures (for example, zebrin expression)<sup>47</sup>. Future studies addressing whether the increased myelin thickness in the cerebellum of plKO mice is global or restricted to functional zones, and whether it correlates with Purkinje cell firing frequency, might elucidate how placental ALLO insufficiency alters specific cerebellar functions that are linked to ASD.

Our experiments demonstrate a treatment-reversible link between specific placental ALLO insufficiency and sex-specific cerebellar myelination associated with an ASD-like phenotype. It is well documented that female and male brains might take different developmental trajectories in physiological and pathological conditions<sup>48</sup>, but there is still limited mechanistic information on the underlying causes of sex differences<sup>49</sup>. The results we describe here suggest that, in the case of ALLO insufficiency and cerebellar myelination, this sex-linked divergence occurs primarily in the postnatal period. Furthermore, the sex-linked phenotypes of our mouse model parallel human pathological features, strengthening the translational value of our findings. These experiments suggest that the cerebellum is the primary brain region affected by the lack of placental ALLO in males. Other brain regions might be altered specifically in females, as our transcriptomic assessment suggests (Extended Data Fig. 5). Undoubtedly, such complex sexual dimorphism builds upon the interplay of a myriad of factors throughout prenatal and postnatal development. Further investigation of brain region-specific, sex-linked trajectories under normal and pathological conditions can pave the way for personalized risk
assessment, diagnosis and treatment informed by patient sex and hormone levels.

Increased risk of ASD diagnosis after extremely preterm birth⁵⁰, particularly in males, suggests that loss of placental hormones might contribute to specific human neurobehavioral outcomes in previously unanticipated ways. Our findings lay the ground-work for developing hormone replacement strategies to normalize the developmental milieu and prevent long-term impairments in neurobehavior.

## Online content

Any methods, additional references, Nature Research report-ing summaries, source data, extended data, supplementary infor-mation, acknowledgements, peer review information; details of author contributions and competing interests; and statements of data and code availability are available at https://doi.org/10.1038/s41593-021-00896-4.

Received: 3 March 2021; Accepted: 23 June 2021;

Published online: 16 August 2021

## References

1. Kratimenos, P. & Penn, A. A. Placental programming of neuropsychiatric disease. *Pediatr. Res.* **86**, 157–164 (2019).

2. Hosie, A. M., Wilkins, M. E., da Silva, H. M. & Smart, T. G. Endogenous steroids regulate GABA<sub>A</sub> receptors through two discrete transmembrane sites. *Nature* **444**, 486–489 (2006).

3. Beledi, D., Hogenkamp, D., Gee, K. W. & Lambert, J. J. Realising the therapeutic potential of neuroactive steroid modulators of the GABA<sub>A</sub> receptor. *Neurobiol. Stress* **12**, 100207 (2020).

4. Schumacher, M. et al. Revisiting the roles of progesterone and allopregnanolone in the nervous system: resurgence of the progesterone receptors. *Prog. Neurobiol.* **113**, 6–39 (2014).

5. Brunton, P. J., Russell, J. A. & Hirst, J. J. Allopregnanolone in the brain: protecting pregnancy and birth outcomes. *Prog. Neurobiol.* **113**, 106–136 (2014).

6. Wenzel, P. L. & Leone, G. Expression of Cre recombinase in early diploid trophoblast cells of the mouse placenta. *Genesis* **45**, 129–134 (2007).

7. Marques, S. et al. Oligodendrocyte heterogeneity in the mouse juvenile and adult central nervous system. *Science* **352**, 1326–1329 (2016).

8. Thakurela, S. et al. The transcriptome of mouse central nervous system myelin. *Sci. Rep.* **6**, 25828 (2016).

9. Shaw, J. C., Palliser, H. K., Palazzi, K. & Hirst, J. J. Administration of progesterone throughout pregnancy increases maternal steroids without adverse effect on mature oligodendrocyte immunostaining in the guinea pig. *Reprod. Sci.* **25**, 395–405 (2018).

10. Marton, L. S. & Stefansson, K. Developmental alterations in molecular weights of proteins in the human central nervous system that react with antibodies against myelin-associated glycoprotein. *J. Cell Biol.* **99**, 1642–1646 (1984).

11. Chang, E. H. et al. The role of myelination in measures of white matter integrity: combination of diffusion tensor imaging and two-photon microscopy of CLARITY intact brains. *Neuroimage* **147**, 253–261 (2017).

12. Baumann, N. & Pham-Dinh, D. Biology of oligodendrocyte and myelin in the mammalian central nervous system. *Physiol. Rev.* **81**, 871–927 (2001).

13. Chen, V. S. et al. Histology atlas of the developing prenatal and postnatal mouse central nervous system, with emphasis on prenatal days E7.5 to E18.5. *Toxicol. Pathol.* **45**, 705–744 (2017).

14. van Tilborg, E. et al. Origin and dynamics of oligodendrocytes in the developing brain: Implications for perinatal white matter injury. *Glia* **66**, 221–238 (2018).

15. Hamilton, N. B. et al. Endogenous GABA controls oligodendrocyte lineage cell number, myelination, and CNS internode length. *Glia* **65**, 309–321 (2017).

16. Zonouzi, M. et al. GABAergic regulation of cerebellar NG2 cell development is altered in perinatal white matter injury. *Nat. Neurosci.* **18**, 674–682 (2015).

17. Sathyanesan, A., Kundu, S., Abbah, J. & Gallo, V. Neonatal brain injury causes cerebellar learning deficits and Purkinje cell dysfunction. *Nat. Commun.* **9**, 3235 (2018).

18. Grasselli, G. et al. SK2 channels in cerebellar Purkinje cells contribute to excitability modulation in motor-learning-specific memory traces. *PLoS Biol.* **18**, e3000596 (2020).

19. Piochon, C. et al. Cerebellar plasticity and motor learning deficits in a copy-number variation mouse model of autism. *Nat. Commun.* **5**, 5586 (2014).

20. Vinieza Veloz, M. F. et al. Cerebellar control of gait and interlimb coordination. *Brain Struct. Funct.* **220**, 3513–3536 (2015).

21. Kindregan, D., Gallagher, L. & Gormley, J. Gait deviations in children with autism spectrum disorders: a review. *Autism Res Treat.* **2015**, 741480 (2015).

22. Scholz, J., Niibori, Y., Frankland, P. W. & Lerch, J. P. Rotarod training in mice is associated with changes in brain structure observable with multimodal MRI. *Neuroimage* **107**, 182–189 (2015).

23. Rothwell, P. E. et al. Autism-associated neuroligin-3 mutations commonly impair striatal circuits to boost repetitive behaviors. *Cell* **158**, 198–212 (2014).

24. Stoodley, C. J. et al. Altered cerebellar connectivity in autism and cerebellar-mediated rescue of autism-related behaviors in mice. *Nat. Neurosci.* **20**, 1744–1751 (2017).

25. Wang, S. S., Kloth, A. D. & Badura, A. The cerebellum, sensitive periods, and autism. *Neuron* **83**, 518–532 (2014).

26. Courchesne, E. et al. Unusual brain growth patterns in early life in patients with autistic disorder: an MRI study. *Neurology* **57**, 245–254 (2001).

27. Zeidan-Chulia, F. et al. Up-regulation of oligodendrocyte lineage markers in the cerebellum of autistic patients: evidence from network analysis of gene expression. *Mol. Neurobiol.* **53**, 4019–4025 (2016).

28. Silverman, J. L., Yang, M., Lord, C. & Crawley, J. N. Behavioural phenotyping assays for mouse models of autism. *Nat. Rev. Neurosci.* **11**, 490–502 (2010).

29. Gallo, I. et al. Formyl peptide receptor as a novel therapeutic target for anxiety-related disorders. *PLoS ONE* **9**, e114626 (2014).

30. Kim, H., Lim, C. S. & Kaang, B. K. Neuronal mechanisms and circuits underlying repetitive behaviors in mouse models of autism spectrum disorder. *Behav. Brain Funct.* **12**, 3 (2016).

31. Lugo, J. N. et al. Deletion of PTEN produces autism-like behavioral deficits and alterations in synaptic proteins. *Front. Mol. Neurosci.* **7**, 27 (2014).

32. Premoli, M. et al. Specific profile of ultrasonic communication in a mouse model of neurodevelopmental disorders. *Sci. Rep.* **9**, 15912 (2019).

33. El-Kordi, A. et al. Development of an autism severity score for mice using *Nlgn4* null mutants as a construct-valid model of heritable monogenic autism. *Behav. Brain Res.* **251**, 41–49 (2013).

34. Fraser, M. M., Bayazitov, I. T., Zakharenko, S. S. & Baker, S. J. Phosphatase and tensin homolog, deleted on chromosome 10 deficiency in brain causes defects in synaptic structure, transmission and plasticity, and myelination abnormalities. *Neuroscience* **151**, 476–488 (2008).

35. Pacey, L. K. et al. Delayed myelination in a mouse model of fragile X syndrome. *Hum. Mol. Genet* **22**, 3920–3930 (2013).

36. Uddin, L. Q., Supekar, K. & Menon, V. Reconceptualizing functional brain connectivity in autism from a developmental perspective. *Front Hum. Neurosci.* **7**, 458 (2013).

37. Dachew, B. A., Mamun, A., Maravilla, J. C. & Alati, R. Pre-eclampsia and the risk of autism-spectrum disorder in offspring: meta-analysis. *Br. J. Psychiatry* **212**, 142–147 (2018).

38. Guennoun, R. Progesterone in the brain: hormone, neurosteroid and neuroprotectant. *Int. J. Mol. Sci.* **21**, 5271 (2020).

39. Pang, Y., Dong, J. & Thomas, P. Characterization, neurosteroid binding and brain distribution of human membrane progesterone receptors δ and ε (mPRδ and mPRε) and mPRδ involvement in neurosteroid inhibition of apoptosis. *Endocrinology* **154**, 283–295 (2013).

40. Thomas, P. & Pang, Y. Anti-apoptotic actions of allopregnanolone and ganaxolone mediated through membrane progesterone receptors (PAQRs) in neuronal cells. *Front. Endocrinol. (Lausanne)* **11**, 417 (2020).

41. Baron-Cohen, S. et al. Foetal oestrogens and autism. *Mol. Psychiatry* **25**, 2970–2978 (2020).

42. Agrawal, S., Rao, S. C., Bulsara, M. K. & Patole, S. K. Prevalence of autism spectrum disorder in preterm infants: a meta-analysis. *Pediatrics* **142**, e20180134 (2018).

43. Hertig, A. et al. Steroid profiling in preeclamptic women: evidence for aromatase deficiency. *Am. J. Obstet. Gynecol.* **203**, 477 e471–477 (2010).

44. Becker, E. B. & Stoodley, C. J. Autism spectrum disorder and the cerebellum. *Int. Rev. Neurobiol.* **113**, 1–34 (2013).

45. Fatemi, S. H. et al. Consensus paper: pathological role of the cerebellum in autism. *Cerebellum* **11**, 777–807 (2012).

46. Rogers, T. D. et al. Is autism a disease of the cerebellum? An integration of clinical and pre-clinical research. *Front. Syst. Neurosci.* **7**, 15 (2013).

47. De Zeeuw, C. I. & Ten Brinke, M. M. Motor learning and the cerebellum. *Cold Spring Harb. Perspect. Biol.* **7**, a021683 (2015).

48. McCarthy, M. M. Sex differences in the developing brain as a source of inherent risk. *Dialogues Clin. Neurosci.* **18**, 361–372 (2016).

49. McCarthy, M. M., Arnold, A. P., Ball, G. F., Blaustein, J. D. & De Vries, G. J. Sex differences in the brain: the not so inconvenient truth. *J. Neurosci.* **32**, 2241–2247 (2012).

50. Limperopoulos, C. Autism spectrum disorders in survivors of extreme prematurity. *Clin. Perinatol.* **36**, 791–805 (2009).

**Publisher's note** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
<figure><img src="image_47.png" /></figure>

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons license,

unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2021
## Methods

Mice. All the procedures on experimental animals were performed in accordance with the protocol approved by the Institutional Animal Care and Use Committee at Children's National Medical Center (protocol no. 30534 (PI: Penn)) and at Columbia University Medical Center (protocol nos. AC-AABE 6553 (PI: Penn) and AC-AABF 5550 (PI: Yang)). Mice were housed on a 12-h light/dark cycle and had ad libitum access to food and water.

The *akr1c14*-floxed mouse was constructed with inGenious Targeting Laboratories. The construct targeted a 3.15-kb region that includes exons 7–9 of *akr1c14*, with a 5' homology arm extending to the position of a LoxP-FRT flanked Neo cassette and a 3' homology arm extending to a single LoxP site just 5' of exon 10 (Fig. 1b). The targeting vector was electroporated into embryonic stem cells, which were cloned and screened. Confirmed clones were micro-injected into C57Bl/6 blastocysts. The Neo cassette was later removed using flipase. The resulting mice are designated as the *akr1c14*-floxed line. *akr1c14*<sup>fl/fl</sup> and *akr1c14*<sup>fl/fl</sup> mice appear phenotypically normal. The final placental 3α-HSD knockout was generated by successive crosses of progenitor *akr1c14*-floxed mice and Cyp19a-Cre mice (generous gift of G. Leone), which express Cre exclusively in placental trophoblast cells<sup>5</sup>. Mouse, unlike human, does not express endogenous placental Cyp19a, which codes for aromatase, but the human Cyp19a promotor effectively drives placenta-specific Cre expression<sup>6</sup>. Cyp19a-Cre females were crossed to homozygous *akr1c14*<sup>fl/fl</sup> males. Female offspring of these crosses that were positive for Cre with two floxed alleles were used as the dams of experimental mice. These *akr1c14*<sup>fl/fl</sup> females carrying Cyp19a-Cre were crossed to *akr1c14*<sup>fl/fl</sup> males, so that half of resulting offspring are both homozygous *akr1c14*<sup>fl/fl</sup> and positive for placental Cre. In tests of recombination specificity (Fig. 1e), the presence of both LoxP sites (indicating un-recombined DNA) was ascertained using a primer set, including a forward primer homologous to a region 5' of exon 7, and a reverse primer incorporating the 3' LoxP site between exons 9 and 10. This primer pair was specific to the floxed allele due to the inclusion of the LoxP site in the reverse primer sequence. It did not amplify the wild-type allele. The presence of the recombined LoxP site was ascertained using a primer set, which flanks exon 10 such that the recombined allele (164 bp in length) amplifies, whereas the intact allele is too long to amplify (3,396 bp). These mice in which *akr1c14* is deleted specifically in the placenta are *akr1c14*<sup>Cyp19a-/-</sup>, designated as *akr1c14*<sup>Cyp19aKO</sup> (pIKO) for simplicity. An additional half of offspring are negative for placental Cre and designated as controls (C).

Cyp19a-Cre transgenic mice were crossed with ROSA26-YFP mice (The Jackson Laboratory, no. 006148) to obtain Cyp19a-Cre:R26R-EYFP littermates.

**Human samples.** Human cerebellar vermes were obtained from the National Institutes of Health (NIH) NeuroBioBank (request ID no. 709). Donors consisted of 6-week-old term infants (average age) and corrected age-matched preterm infants, excluding those with major congenital anomalies or known genetic diagnoses and those with meningitis or stroke as cause of death. Sex, gestational weeks, absolute age (postnatal age) and corrected age are shown in Supplementary Table 3. Frozen tissues were preserved at –80 °C.

**Drug injections.** Pregnant dams received one intraperitoneal injection of allopregnanolone (Tocris Bioscience, no. 3653) diluted in sesame oil (vehicle; Sigma-Aldrich, no. S3547) at E15.5 at a dose of 10 mg kg⁻¹ of body weight. Muscimol (Tocris, no. 0289) was diluted in saline solution and injected intraperitoneally in pregnant mice at E15.5 (1 mg kg⁻¹). Allopregnanolone and muscimol injections were done in the morning, during sleep time, to limit the effect of the drug-induced sedation on dams' activity. Doses were chosen based on prior pharmacological studies in mice<sup>31,52</sup>, and the timing was an empirical choice based on *akr1c14* gene expression peak at E14.5 (Fig. 1a). BrdU (50 mg kg⁻¹) dissolved in saline solution was injected intraperitoneally in dams at E15.5.

**Mass spectrometry.** *Steroid extraction.* Pregnenolone (PREG), PROG, ALLO, epiallopregnanolone (EPIALLO), ALLOTHDOC and 3α5α-tetrahydrotestosterone (3α5α-THT) were identified and quantified simultaneously in individual tissues by gas chromatography–tandem mass spectrometry (GC–MS/MS) as previously described<sup>53</sup>. Placental (65–185 mg) and fetal brain (61–130 mg) tissues of male and female wild-type and pIKO mice at E17.5 were weighed and stored at –20 °C until GC–MS/MS analysis. Briefly, steroids were first extracted from placentas, and brains with 10 volumes of methanol (MeOH) and the following internal standards were added to the extracts for steroid quantification: 2 ng of epietocholanolone for PREG, ALLO, EPIALLO, ALLOTHDOC and 3α5α-THT and 2 ng of ¹³C₂-PROG for PROG. Samples were purified and fractionated by solid-phase extraction with the recycling procedure<sup>54</sup>. The extracts were dissolved in 1 ml of MeOH and applied to the C18 cartridge (500 mg, 6 ml, International Sorbent Technology), followed by 5 ml of MeOH/H₂O (85/15, vol/vol). The flow-through, containing the free steroids, was collected and dried. After a previous re-conditioning of the same cartridge with 5 ml of H₂O, the dried samples were dissolved in MeOH/H₂O (2/8, vol/vol) and re-applied. The cartridge was then washed with 5 ml of H₂O and 5 ml of MeOH/H₂O (1/1, vol/vol), and unconjugated steroids were eluted with 5 ml of MeOH/H₂O (9/1, vol/vol). The fraction containing the unconjugated steroids was then filtered and further purified by high-performance liquid chromatography

(HPLC). The HPLC system is composed of a WPS-3000SL analytical autosampler and an LPG-3400SD quaternary pump gradient coupled with a SR-3000 fraction collector (Thermo Fisher Scientific). The HPLC separation was achieved with a LiChrosorb Diol column (25 cm, 4.6 mm, 5 µm) in a thermostated block at 30 °C. The column was equilibrated in a solvent system of 90% heptane and 10% of a mixture composed of heptane/isopropanol (85/15, vol/vol). Elution was performed at a flow rate of 1 ml min⁻¹, first 90% heptane and 10% of heptane/isopropanol (85/15, vol/vol) for 15 min, and then with a linear gradient to 100% acetone in 2 min. The column was washed with acetone for 15 min. Steroids were collected in the time range of 15–29 min and were derivatized with 25 µl of HFBA and 25 µl of anhydrous acetone for 1 h at 20 °C. Samples were dried under a stream of nitrogen and resuspended in heptane.

**GC–MS/MS analysis.** GC–MS/MS analysis of the biological extracts was performed using an A1310 autosampler, a Trace 1310 gas chromatograph and a TSQ 8000 mass spectrometer (Thermo Fisher Scientific). Injection was performed in the splitless mode at 250 °C (1 min of splitless time), and the temperature of the gas chromatograph oven was initially maintained at 80 °C for 1 min and ramped between 80 °C and 200 °C at 20 °C min⁻¹ and then ramped to 300 °C at 5 °C min⁻¹ and finally ramped to 350 °C at 30 °C min⁻¹. The helium carrier gas flow was maintained constant at 1 ml min⁻¹ during the analysis. The transfer line and ionization chamber temperatures were 300 °C and 200 °C, respectively. Electron impact ionization was used for mass spectrometry with ionization energy of 70 eV, and GC–MS/MS analysis was performed in multiple reaction monitoring mode with argon as the collision gas. GC–MS/MS signals were evaluated using a computer workstation by means of the software Excalibur, release 3.0 (Thermo Fisher Scientific). Identification of steroids was supported by their retention time and two or three transitions. Quantification was performed according to the more abundant transition with a previously established calibration curve. The range of the limit of detection was roughly 0.5–20 pg according to the steroid structure. The GC–MS/MS analytical procedure was fully validated in terms of accuracy, reproducibility and linearity in mouse brain<sup>53</sup>.

**3' mRNA sequencing.** Tissue collection and RNA extraction: Cerebellum, bilateral cerebral cortices and bilateral hippocampi from C and pIKO mice (males and females) at P30 (three animals per group) were dissected and flash-frozen. Tissue homogenization and total RNA extraction were performed using the *mirVana Isolation Kit* (Thermo Fisher Scientific, no. AM1560) according to the manufacturer's instructions.

**Library preparation and sequencing for mRNA.** The cDNA libraries were prepared using the QuantSeq 3' mRNA-Seq Library Prep Kit FWD for Illumina (Lexogen) as per the manufacturer's instructions. Briefly, total RNA was reverse transcribed using oligo (dT) primers. The second cDNA strand was synthesized by random priming, in a manner that DNA polymerase is efficiently stopped when reaching the next hybridized random primer; therefore, only the fragment closed to the 3' end was captured for indexed adapter ligation and PCR amplification. The processed libraries were assessed for its size distribution and concentration using BioAnalyzer High Sensitivity DNA Kit (Agilent Technologies, no. 5067-4626). Pooled libraries were diluted to 2 nM in EB buffer (Qiagen, no. 19086) and then denatured using the Illumina protocol. The libraries were pooled and diluted to 3 nM using 10 mM Tris-HCl, pH 8.5, and then denatured using the Illumina protocol. The denatured libraries were diluted to 10 pM by pre-chilled hybridization buffer and loaded onto a TruSeq v2 Rapid flow cell on an Illumina HiSeq 2500 and run for 50 cycles using a single-read recipe according to the manufacturer's instructions. De-multiplexed sequencing reads were generated using Illumina bcl2fastq (release version 2.18.0.12) allowing no mismatches in the index read.

**Data analysis.** After the quality and polyA trimming by BBDuk<sup>55</sup> and alignment by HISAT2 (version 2.1.0)<sup>56</sup>, read counts were calculated using HTSeq<sup>57</sup> by supplementing Ensembl gene annotation (GRCm38.78). DE analysis and MA plots were done using TCC R package (version 1.12.1, https://bioconductor.rikен.jp/packages/3.3/bioc/manuals/TCC/man/TCC.pdf)<sup>58</sup>. After getting DEGs, hierarchical clustering of DEGs as a heat map was performed using Partek Genomics Suite (version 6.6, http://www.partek.com/partek-genomics-suite/) (Partek). No outlier was identified in our biological replicates, as indicated by principal component analysis and scatter matrices. DEGs met the following guidelines: *P* < 0.05; fragments per kilobase of transcript per million mapped reads values >1; and fold change > 1.5. The limits of these cutoffs were validated by RT-PCR on samples from other mouse cohorts. IPA was used to identify the top biological functions and disease processes that were differentially regulated. This software is based on the Ingenuity Pathway Knowledge Base for genetic interaction, which derives from the scientific literature, each network connection being supported by previous publications<sup>59</sup>. To determine and visualize the degree of gene overlaps in datasets, Venn analysis was performed using Venny 2.1 (http://bioinfogp.cnb.csic.es/tools/venny/).

**RT-PCR.** Tissues were homogenized in TRIzol Reagent (Thermo Fisher Scientific); total RNA was extracted with the RNeasy Mini Kit (Qiagen, no. 74104).
Next, 1 µg of RNA was used to make cDNA with the iScript cDNA Synthesis Kit (Bio-Rad, no. 1708891). All primer pairs were designed and validated in-house for efficiency and specificity. RT-PCR experiments were performed on cDNA samples in the presence of PowerUp SYBR Green Master Mix (Thermo Fisher Scientific, no. 1725271) with specific primers at 100 nM using the CFX96 Touch Real-Time PCR Detection System (Bio-Rad). The cDNA-generated signals for target genes were internally corrected with phosphoglycerate kinase 1 (pgk1) in mouse postnatal brains, tyrosine 3-monooxygenase/tryptophan 5-monooxygenase activation protein zeta (ywhaz) in mouse placentas and embryonic brains or cytochrome C1 (cyc1) in human tissues. The regulation fold changes were determined with the 2⁻⁵ΔCq method⁶⁰.

**MRI. Brain preparation.** Brains left within the skull were prepared for MRI using the following method. Mice were anesthetized and intracardially perfused at a rate of 1 ml min⁻¹, with 30 ml of 0.1 M PBS containing 10 U/100 ml heparin (Sigma-Aldrich) and 2 mM ProHance (gadolinium contrast agent, Bracco Diagnostics) followed by 30 ml of 4% paraformaldehyde (PFA) and 2 mM ProHance. After perfusion, the mice were decapitated, and the skin and the lower jaw were removed. Brains in skulls were post-fixed in 4% PFA and 2 mM ProHance at 4 °C overnight and then kept in 0.1 M PBS containing 2 mM ProHance and 0.02% sodium azide for at least 1 month before MRI scanning⁶¹,⁶².

**DTI acquisition.** A multi-channel 7-Tesla MRI scanner (Agilent) was used to image the brains. For DTI Imaging, 12 brains were scanned in parallel using a 3D diffusion-weighted fast spin echo scan⁶³,⁶⁴. Parameters for the DTI sequence were as follows: TR of 270 ms, echo train length of 6, first TE of 32 ms and a TE of 10 ms for the remaining five echoes, 1 average, field of view of 14 mm × 14 mm × 25 mm and a matrix size of 180 × 180 × 324, resulting in an image with 78-µm isotropic voxels. Five *b*=0 s mm⁻² images and 30 high *b*-value (*b* = 2,147 s mm⁻²) images in 30 different directions were acquired, using the Jones30 scheme⁶⁵. Total imaging time was approximately 12 h. The FSL software package (FMRIB) was used to generate FA and mean diffusivity (MD) maps for each brain sample.

**MRI registration and analysis.** To assess any changes to the mouse brains due to genotype and sex, the *b*=0 s mm⁻² images were registered linearly (6 followed by 12 parameter) and non-linearly together. A combination of mni_autoreg tools⁶⁶ and ANTs (advanced normalization tools)⁶⁷ were used to perform the registrations. A population atlas, representing the average anatomy of the study samples, was created when all scans were resampled with an appropriate transform. The final registration results are the individual images deformed into alignment with one another in an unbiased manner. The analysis of the individual DTI parameters compare intensity differences of FA and MD between genotype and sex. Warping an available classified MRI mouse brain atlas onto the population average allows for significant differences in the diffusion measures of segmented structures to be calculated⁶⁸⁻⁷¹.

**Histology. Section preparation and staining.** Animals were anesthetized using isoflurane (Isothix, Henry Schein Animal Health) and transcardially perfused with 20 ml of 1× PBS followed by 30 ml of 4% PFA, as previously described. Brains were post-fixed for 24 h in 4% PFA and cryoprotected in 20% sucrose in PBS. Serial 40-µm-thick sagittal brain/cerebellar sections were then collected using a sliding microtome. Placentas were collected at E17.5 after deep anaesthesia of the dams, drop-fixed in 4% PFA and cryoprotected in 30% sucrose. Frozen 20-µm cross-sections were obtained with a cryostat. Immunohistochemistry on brain or placenta sections was performed using the following antibodies: rabbit anti-MBP (Abcam, no. ab40390, 1:500), mouse anti-APC clone CC1 (EMD Millipore, no. MABC20, 1:500), rabbit anti-Olig2 (Abcam, no. ab9610, 1:500), mouse anti-neuronal nuclei (NeuN) (Millipore, no. MAB377, 1:500), goat anti-NeuroD1 (RD Systems, no. AF2746, 1:500), rabbit anti-calbindin (Swant, no. CB38, 1:1,000), rat anti-mouse PDGFRα (CD140a; BD Biosciences, 17-1401-81, 1:500) and chicken anti-GFP (Abcam, no. ab13970, 1:500). Sections were incubated with primary antibodies overnight at 4 °C containing 0.3% Triton and 10% normal donkey serum. Sections were then incubated with secondary antibodies (1:500) together with DAPI (Invitrogen, no. D1306, 1:1,000) for 2 h at room temperature. FluoroMyelin Green (Invitrogen, F34651, 1:500) was also added at this time for those sections. Secondary antibodies were used as follows: donkey anti-mouse Alexa-488 (Invitrogen, A-21202), donkey anti-rabbit Alexa-488 (Invitrogen, A-21206), donkey anti-mouse Alexa-555 (Invitrogen, A-31570), donkey anti-rabbit Alexa-555 (Invitrogen, A-31572), donkey anti-mouse Alexa-647 (Invitrogen, A-31571) and donkey anti-rabbit Alexa-647 (Invitrogen, A-31573). Finally, floating sections were mounted and cover-slipped using ProLong Gold Antifade Mountant (Thermo Fisher Scientific, P36930) before imaging.

**In situ hybridization.** For analysis of *akr1c14*-mRNA, E17.5 placental tissue was drop-fixed overnight in 4% PFA and then transferred through a sucrose gradient of 10%, 20% and 30% over 3 d. Tissue was cross-sectioned at 12 µm and assayed for *akr1c14*-mRNA using the Advanced Cell Diagnostics RNAScope system and an ACD probe designed to target the 621–1,989 region of *akr1c14* (NM_134072.1).

**Fluorescence imaging and cell counting.** Placenta cross-sections and whole brain or cerebellar sagittal sections were imaged using a virtual slide microscope using cellSens version 2.3 (VS 120, Olympus Life Science) under a ×20 objective. High-magnification images were taken under a confocal microscope with LAS X Life Science Microscope software version 2.7 (TCS SP8, Leica Microsystems). z-stack images were acquired with a step size of 1.55 µm and viewed using NIH ImageJ 1.53c (http://imagej.nih.gov/ij). The different cerebellar layers were delineated manually using the freehand selection tool in ImageJ. The density of Olig2-, CC1-, PDGFRα- and BrdU-positive cells was estimated using Stereo Investigator (MBF Bioscience). After delineating the cerebellar WM using the freehand selection tool on corresponding DAPI images, immunopositive cells were manually counted under the ×40 objective (Zeiss Axio Imager M2) using the optical fractional as a probe in grids randomly distributed and over covering 30% of the cerebellar WM. Cell densities were then calculated by using the individual WM area values and the mean measured section thickness. The linear density of Purkinje cells (Purkinje cells per millimeter) was determined in the lobule VI–VII by drawing a freehand line through the center of the cell bodies.

**Electron microscopy. Sample preparation.** P30 mice were anesthetized with a 10:1 mixture of ketamine/xylazine in saline vehicle such that each mouse receives 100 mg kg⁻¹ of ketamine and 10 mg kg⁻¹ of xylazine, respectively. Under deep anaesthesia, mice were perfused with 20 ml of 0.12 M cacodylate buffer followed by 30 ml of fixative made of 2.5% glutaraldehyde and 1% PFA (electron microscopy grade, EMS nos. 16221 and 15713) in cacodylate buffer (pH 7.4). Brains were removed and post-fixed for 1 h in the same fixative. Next, 300-µm-thick sagittal slices were made using a vibratome (VT 1000S, Leica Biosystems). Slices were post-fixed with 1% osmium tetroxide in 0.12 M cacodylate buffer for 2 h at room temperature and ‘en block’ stained with 1% uranyl acetate in 0.1 M acetate buffer overnight at room temperature. After several washes in acetate buffer, slices were dehydrated by passing them through increasing concentrations of ethyl alcohol (up to 100%). Slices were then progressively infiltrated in epoxy resin and placed for 48 h at 60 °C for resin polymerization. Then, 100-nm-thick sagittal ultrasects were performed using an ultramicrotome (Ultraprot UC7, Leica Biosystems) through the whole cerebellum (+ brainstem) or the whole anterior brain (for corpus callosum analyses).

**Image acquisition.** Large ultrathin sections absorbed on silicon were observed with a FEI Helios NanoLab 660 FIBSEM field emission scanning electron microscope (FEI, Thermo Fisher Scientific) using high-resolution immersion mode and equipped with a solid-state concentric ring (insertable) back-scattering electron detector. Image acquisition was done using MAPS Software version 3.7 (Thermo Fisher Scientific). Low-magnification images (×600) were first taken to delineate our regions of interest (vermal inter-lobule 6–7 of the cerebellum). Then, high magnifications of two-dimensional image registration (×20,000–50,000) were taken using 4 kV/0.2 nAmp landing electron beam.

**Myelin thickness measurement.** The outer diameter (including myelin sheath) and inner diameter of at least 400 randomly selected myelinated axons were measured on one ultrathin section per animal (three animals per group). The g-ratio (equal to the ratio of the inner-to-outer diameter of a myelinated axon) of each axon was calculated. Scatter plots of the axon caliber and g-ratios were then analyzed. The slopes and elevations of the linear regressions from C and pIKO mice were compared.

**Western blots.** Whole mouse cerebellums and human vermal samples were homogenized in 250 µl of radioimmunoprecipitation assay lysis buffer consisting of (in mM) 50 Tris-HCl, pH 7.4, 150 NaCl, 2 EDTA, 50 NaF, 1 Na₂VO₃, 1% Triton X-100, 0.1% SDS, 0.5% Na-deoxycholate and a Protease/Phosphatase Inhibitor Cocktail (Santa Cruz Biotechnology). After centrifugation at 14,000g for 10 min, protein concentration was determined using Bradford protein assay kit (Bio-Rad). Sample total proteins were resolved by sodium dodecyl sulfate-polyacrylamide gel electrophoresis using 10% Bis-Tris precast gel (Thermo Fisher Scientific) and transferred to polyvinylidene fluoride membranes. Membranes were incubated with blocking buffer consisting of 4% non-fat milk in 1% Tween-20 in Tris-buffered saline (TBS-T) for 1 h, followed by overnight incubation at 4 °C with one of the following primary antibodies diluted in 3% BSA TBST-T: rabbit anti-MBP (Abcam, no. ab40390, 1:500), rabbit anti-MAG (Thermo Fisher Scientific, no. PA5-79620, 1:500), rabbit anti-MOG (Abcam, 1:500), mouse anti-cofilin-1 (Cfl1) (Santa Cruz Biotechnology, no. sc-53934, 1:500), mouse anti-prosaposin (Santa Cruz Biotechnology, no. sc-390184, 1:500), rabbit anti-hnRNPK (R332) (Cell Signaling Technology, no. 4675, 1:500) and rabbit anti-GAPDH (Cell Signaling Technology, no. 5174, 1:2,000). After three washes with TBS-T, membranes were incubated with horseradish peroxide-conjugated secondary antibodies (Jackson ImmunoResearch, 1:2,000), and protein bands were visualized using the chemiluminescent ECL detection system (Bio-Rad) according to the manufacturer’s instruction. Signal intensities of protein bands were quantified using Image J software (http://rsb.info.nih.gov/ij/) and normalized with GAPDH as an internal control.
**Behavioral testing.** Tests were performed on 30- and 60-d-old mice. No more than three behavioral tests were done on each mouse, accordingly to our Institutional Animal Care and Use Committee recommendations. A minimal resting period of 3 d was allowed between two consecutive tests.

**Pup ultrasonic vocalizations.** Mouse pups emit USVs when separated from the dam. This innate behavior is essential for eliciting maternal care behaviors that the pups need to survive. Separation-induced vocalizations were tested on postnatal days 4, 6, 8 and 11, with vocalizations declining drastically after P12. USVs were recorded for 3 min, and the number of calls was tracked using Avisoft-SASLab Pro software (Avisoft Bioacoustics).

**Erasmus ladder.** The Erasmus ladder (Noldus) is a fully automated system allowing for the assessment of baseline gait and locomotor coordination and associative cerebellar learning. The instrument consists of a horizontal ladder with 2 × 37 pressure-sensitive rungs for both the left and right sides, stationed between two shelter boxes. Each shelter is equipped with a white LED spotlight and pressurized air outlet serving as cues to signal time to departure from shelters. As the mouse is prompted to move, pressure sensors continuously monitor and analyze the walking pattern of the mouse in real time, sending data to the computer system, which, in turn, adjusts air pressure, calculates interventions and predicts future steps. Through this prediction, the software can also be programmed to move rungs mid-trial by high-speed pneumatic slides to create an obstacle. Mice were trained with 42 runs per day for four successive days without obstacles. Missteps, trial times and jumps were measured for each run. After training, cerebellar associative motor learning was assessed on days 4–8 with the introduction of a tone as a conditioning stimulus (CS) and a rising rung (obstacle) as the unconditioned stimulus (US). These perturbations were used to evaluate mouse ability to adapt motor behavior in real time. Repeated perturbations were used to assess associative motor learning throughout successive trials. With the inter-stimulus interval fixed at 285 ms, mice that associated the tone with the obstacle learned to increase walking speed to avoid being hit. In this sense, they decreased their 'pre-perturbation step time' (step before US) and their 'post-perturbation step time' (step just after the US). The Erasmus ladder is particularly valuable for detailed phenotyping of deficits associated with cerebellar pathology¹⁷.

**Rotarod.** The accelerating rotarod performance test is a standard rodent assay for motor-associated functions such as coordination and balance. A 4-d testing paradigm (three sessions per day) was performed as previously described³⁻⁵ with a continuously accelerating rod that forces animals to modulate their balance. For the first 2 d, animals were placed on the rod, which accelerated from 4 to 40 r.p.m. in 180 s. The test was stopped when the animal fell off the rod, was unable to maintain coordinated movements or 180 s passed. Additionally, if animals displayed 'cartwheeling', or clutching the rod and circling without walking, the test was stopped. To increase the challenge, for days 3–4 the acceleration of the rod increased from 8 to 80 r.p.m. in 180 s. The latency to fall from the rod and terminal velocity were then recorded. Linear regression analysis on data from each individual mouse was used to estimate initial motor coordination (given by the intercept) and learning rate (given by the slope).

**Spontaneous behavior/stereotypies.** Assessment of repetitive behavior was done to test autism-like behavior. Mice were allowed to habituate to the testing room before placement in a new, empty cage with bedding materials. A video camera was set with a lateral view to record spontaneous behavior for 15 min. Videos were manually scored for time spent engaging in behaviors, including digging, gnawing, grooming, rearing and jumping.

**Marble burying.** The anxiety and stereotypic and/or obsessive-compulsive-like behavior components of ASD-like behavior can be assessed in mice by their behavioral responses to new 'diggable' media⁷,²³ in the marble burying test. This test was performed in a clean, large (26 × 16 cm) cage filled with 5 cm of bedding and 12 glass marbles evenly spaced on the bedding surface. The animal was left to explore the cage for 30 min undisturbed. Marbles covered 2/3 or more with bedding were counted as buried.

**Socialization test.** The three-chamber test (Crawley's paradigm⁴) was used to assess social behavior as previously described³⁰. A three-chamber Plexiglas box was used with small openings in each dividing wall to allow free access to each chamber. The center chamber was kept empty, whereas the flanking chambers were each equipped with an identical wire cup. After a 10-min habituation period to allow for exploration of the three-chamber equipment (session 1), an unfamiliar adult mouse of similar weight and coloration was placed within the wire cup in one of the side chambers and an 'object' within the wire cup in the opposite chamber. The dividers were raised to allow the test subject to move freely throughout all three chambers of the apparatus over a 10-min test session (session 2). The second session was analyzed for social preference given by the social preference index (SPI). SPI was calculated as follows: if time spent with the mouse (social) is S and the object (non-social) is NS, then $SPI = S/(S + NS)$.

**Determination of the autism composite severity score.** We used the same procedure as previously described in ref. ³³. The score calculation allows for the combination of individual discrete symptoms of the autistic syndrome, in a continuous manner, so that higher values reflect higher severity of autistic-like behavior. For each mouse, we recorded their spontaneous behaviors (grooming, digging, rearing, jumping and gnawing) and determined their SPI in the three-chamber test. The significant readouts (digging time, gnawing time and SPI) were then z-standardized and averaged to establish the autism composite severity score.

**Statistics and reproducibility.** Mouse assignment to groups was not randomized, because it relies on their known genotype (C versus pKO), but mice from at least three different litters were used in each group for all experiments. No litter-associated effect was seen. All data analyses were conducted blinded to group allocation. Statistical analysis was carried out with GraphPad Prism 7 software. Detection of EYFP in Cyp19-cre:R26R-EYFP placenta and brain sections, and in situ hybridizations in placenta sections, were performed in at least three animals. Cell and area quantification on cerebellar immunofluorescent-labelled sections was done on six sections per animal, homogeneously distributed throughout the cerebellar vermis. Myelinated axon g-ratios were measured from more than 400 axons on one electron microscopy scan per animal. For western blot analyses, at least two technical replicates were performed. When more than ten samples, and, thus, several membranes, needed to be compared at once, two common samples were resolved in the different electrophoresis gels and used as reference for quantification. The sample sizes are shown in the legends and chosen to meet or exceed sample sizes typically used in the field. Robust regression and outlier removal methodology with Q value equal to 1 was used to determine outliers. Before conducting all analyses, variable distributions were evaluated for normality. If distribution was not normal, non-parametric methods were used. Differential gene expression across development in C and pKO mice was analyzed using one-way ANOVA with Dunnett's multiple comparisons test. Comparisons between two groups were analyzed using parametric test (unpaired t-test with Welch's correction) or non-parametric test (Mann-Whitney test). The Holm-Sidak method was used for multiple unpaired t-test comparison. Analyses involving data from three or more groups while considering only one independent variable were performed using one-way ANOVA with Tukey's multiple comparisons test. Two-way ANOVA was used to compare means across two or more dependent variables. Two-way repeated-measures ANOVA with Sidak's multiple comparisons tests were used for the Erasmus ladder and the rotarod. Summary data are presented in the text as mean ± s.e.m. from *n* animals. Differences were considered significant at *P* < 0.05. A supplementary statistics methods checklist is available.

**Reporting Summary.** Further information on research design is available in the Nature Research Reporting Summary linked to this article.

## Data availability

Raw and processed data are available from the corresponding authors upon reasonable request. The RNA-seq data discussed in this publication have been deposited in NCBI's Gene Expression Omnibus (Vacher et al., 2021) and are accessible through GEO Series accession number GSE173440. Source data are provided with this paper.

## References

51. Irwin, R. W. & Brinton, R. D. Allopregnanolone as regenerative therapeutic for Alzheimer's disease: translational development and clinical promise. *Prog. Neurobiol.* **113**, 40–55 (2014).

52. Ngo, D. H. & Vo, T. S. An updated review on pharmaceutical properties of gamma-aminobutyric acid. *Molecules* **24**, 2678 (2019).

53. Zhu, X. et al. A role of endogenous progesterone in stroke cerebroprotection revealed by the neural-specific deletion of its intracellular receptors. *J. Neurosci.* **37**, 10998–11020 (2017).

54. Liere, P. et al. Novel lipoidal derivatives of pregnenolone and dehydroepiandrosterone and absence of their sulfated counterparts in rodent brain. *J. Lipid Res.* **45**, 2287–2302 (2004).

55. Kechin, A., Boyarskikh, U., Kel, A. & Filipenko, M. cutPrimers: a new tool for accurate cutting of primers from reads of targeted next generation sequencing. *J. Comput Biol.* **24**, 1138–1143 (2017).

56. Goldstein, L. D. et al. Prediction and quantification of splice events from RNA-seq data. *PLoS ONE* **11**, e0156132 (2016).

57. Anders, S., Pyl, P. T. & Huber, W. HTSeq—a Python framework to work with high-throughput sequencing data. *Bioinformatics* **31**, 166–169 (2015).

58. Sun, J., Nishiyama, T., Shimizu, K. & Kadota, K. TCC: an R package for comparing tag count data with robust normalization strategies. *BMC Bioinformatics* **14**, 219 (2013).

59. Kramer, A., Green, J., Pollard, J. Jr. & Tugendreich, S. Causal analysis approaches in ingenuity pathway analysis. *Bioinformatics* **30**, 523–530 (2014).

60. Livak, K. J. & Schmittgen, T. D. Analysis of relative gene expression data using real-time quantitative PCR and the $2^{-\Delta\Delta C_T}$ method. *Methods* **25**, 402–408 (2001).
61. Cahill, L. S. et al. Preparation of fixed mouse brains for MRI. *Neuroimage* **60**, 933–939 (2012).

62. de Guzman, A. E., Wong, M. D., Gleave, J. A. & Nieman, B. J. Variations in post-perfusion immersion fixation and storage alter MRI measurements of mouse brain morphometry. *Neuroimage* **142**, 687–695 (2016).

63. Lerch, J. P., Sled, J. G. & Henkelman, R. M. MRI phenotyping of genetically altered mice. *Methods Mol. Biol.* **711**, 349–361 (2011).

64. Nieman, B. J., Flenniken, A. M., Adamson, S. L., Henkelman, R. M. & Sled, J. G. Anatomical phenotyping in the brain and skull of a mutant mouse by magnetic resonance imaging and computed tomography. *Physiol. Genomics* **24**, 154–162 (2006).

65. Jones, D. K., Horsfield, M. A. & Simmons, A. Optimal strategies for measuring diffusion in anisotropic systems by magnetic resonance imaging. *Magn. Reson. Med.* **42**, 515–525 (1999).

66. Collins, D. L., Neelin, P., Peters, T. M. & Evans, A. C. Automatic 3D intersubject registration of MR volumetric data in standardized Talairach space. *J. Comput. Assist. Tomogr.* **18**, 192–205 (1994).

67. Avants, B. B. et al. A reproducible evaluation of ANTs similarity metric performance in brain image registration. *Neuroimage* **54**, 2033–2044 (2011).

68. Dorr, A. E., Lerch, J. P., Spring, S., Kabani, N. & Henkelman, R. M. High resolution three-dimensional brain atlas using an average magnetic resonance image of 40 adult C57BL/6J mice. *Neuroimage* **42**, 60–69 (2008).

69. Richards, K. et al. Segmentation of the mouse hippocampal formation in magnetic resonance images. *Neuroimage* **58**, 732–740 (2011).

70. Steadman, P. E. et al. Genetic effects on cerebellar structure across mouse models of autism using a magnetic resonance imaging atlas. *Autism Res.* **7**, 124–137 (2014).

71. Ullmann, J. F., Watson, C., Janke, A. L., Kurniawan, N. D. & Reutens, D. C. A segmentation protocol and MRI atlas of the C57BL/6J mouse neocortex. *Neuroimage* **78**, 196–203 (2013).

72. Egashira, N. et al. Impaired social interaction and reduced anxiety-related behavior in vasopressin V1a receptor knockout mice. *Behav. Brain Res.* **178**, 123–127 (2007).

73. Gyertyan, I. Analysis of the marble burying response: marbles serve to measure digging rather than evoke burying. *Behav. Pharmacol.* **6**, 24–31 (1995).

74. Nadler, J. J. et al. Automated apparatus for quantitation of social approach behaviors in mice. *Genes Brain Behav.* **3**, 303–314 (2004).

## Acknowledgements

We thank G. Leone for providing us with Cyp19a-Cre mice. We also would like to express our gratitude to J. Scafidi, L.-J. Chew and I. Papazoglou for their comments on

the manuscript. We thank M. Yang, the head of the Neurobehavioral Core at Columbia University Medical Center, for giving us access to the equipment. We thank E. Rafikian for scoring the ultrasonic vocalizations. We also thank the contribution of families to the NIH NeuroBioBank, as well as J. Cottrell (NIH NeuroBioBank), who assisted us throughout the process of human tissue selection. We gratefully acknowledge the Children's Research Animal Facility, directed by J. Bradford, for taking care of our mouse colonies and the DC-IDDRC Animal Neurobehavioral Core, headed by J. Corbin, for providing advice, space and equipment to conduct the behavioral phenotyping of our mice. Light and confocal microscopic analyses were carried out at the DC-IDDRC's Cell and Tissue Microscopy Core. The DC-IDDRC cores are supported by National Institutes of Health NICHD U54HD090257 (V.G.). This work was funded by National Institutes of Health grants R01HD092593 and 3R01HD092593-S1 (A.A.P.), R37NS109478 (V.G.) and F31HD098886 (J.J.O.), the Simons Foundation (SFARI Pilot Award no. 572832; A.A.P.), the Children's National Board of Visitors (A.A.P. and V.G.) and the Research Foundation of Cerebral Palsy Alliance (no. 3720; A.A.P.).

## Author contributions

C.-M.V. and A.A.P. conceived the project and designed the experiments. C.-M.V. and J.S. performed the pilot experiments. J.J.O., H.L. and D.B. performed the genetic validation of the model. Y.I., K.H.-T., T.S. and C.-M.V. performed the RNA sequencing analysis. H.L. and S.S. performed the qRT-PCRs. C.-M.V. and J.S. performed the histological experiments. P.K. counted Purkinje cells. C.-M.V., J.S. and S.S. performed the behavioral experiments. A.S. and V.G. contributed to the Erasmus ladder data acquisition and analysis. C.C.-P., A.P. and C.-M.V. performed the electron microscopy experiments. S.S. and H.L. performed the molecular and biochemical experiments on human tissues. P.L. and M.S. performed the mass spectrometry experiments. J.E. and J.P.L. performed the brain MRI. S.S. and J.S. performed the genotyping. C.-M.V. and A.A.P. wrote the manuscript. All authors reviewed and revised the manuscript.

## Competing interests

The authors declare no competing financial interests.

## Additional information

**Extended data** is available for this paper at https://doi.org/10.1038/s41593-021-00896-4.

**Supplementary information** The online version contains supplementary material available at https://doi.org/10.1038/s41593-021-00896-4.

**Correspondence and requests for materials** should be addressed to C.-M.V. or A.A.P.

**Peer review information** *Nature Neuroscience* thanks Margaret McCarthy and the other, anonymous, reviewer(s) for their contribution to the peer review of this work.

**Reprints and permissions information** is available at www.nature.com/reprints.
<figure><img src="image_48.png" /></figure>

<figure><img src="image_49.png" /></figure>

<figure><img src="image_50.png" /></figure>

<figure><img src="image_51.png" /></figure>

<figure><img src="image_52.png" /></figure>

<figure><img src="image_53.png" /></figure>

**Extended Data Fig. 1** | Chemical validation of the plKO model by mass spectrometry. a, Linear regression of ALLO levels in placenta and brain, all sexes and genotypes included. ALLO levels in placenta and brain show a positive correlation (p = 0.0008; n = 18). b–f, Mass spectrometry data in placenta and brain at E17.5. Data is presented as means ± SEM (n = 9 C and 10 plKO). Two-way ANOVA with Sidak's multiple comparison test. Sex mixed as no sex effect was evidenced. ALLO, allopregnanolone; EPIALLO, epiallopregnanolone; PREG, pregnenolone; PROG, progesterone; ALLOTHDOC, allotetrahydrodeoxycorticosterone.
a

efnb3

<figure><img src="image_54.png" /></figure>

b

sbf1

<figure><img src="image_55.png" /></figure>

c

fabp7

<figure><img src="image_56.png" /></figure>

d

tcf7l2

<figure><img src="image_57.png" /></figure>

e

nf1

<figure><img src="image_58.png" /></figure>

f

plp1

<figure><img src="image_59.png" /></figure>

g

<figure><img src="image_60.png" /></figure>

h

<figure><img src="image_61.png" /></figure>

i

<figure><img src="image_62.png" /></figure>

Extended Data Fig. 2 | See next page for caption.
**Extended Data Fig. 2 | Molecular validation of the RNA sequencing cutoffs.** Changes in mRNA (a-f) or protein (g-h) contents of selected cerebellar dysregulated genes after normalization. **a**, *efnb3*, ephrinB3 (RNAseq rank #26/1192 in males; up-regulated; q = 0.018; p = 3.41E-05). **b**, *sbf1*, SET Binding Factor 1 (RNAseq rank #82/1192 in males; down-regulated; q = 0.067; p = 0.00039). **c**, *fabp7*, Fatty Acid Binding Protein 7 (RNAseq rank #84/1192 in males; up-regulated; q = 0.0725; p = 0.00043). **d**, *tcf7l2*, transcription factor 7 like 2 (RNAseq rank #795/1192; up-regulated; q = 0.514; p = 0.029). **e**, *nfl*, Neurofibromin 1 (RNAseq rank #1189/1192 in males; down-regulated; q = 0.59; p = 0.0499). **f**, *plp1*, Proteolipid Protein 1 (RNAseq rank #64/389 in females; down-regulated; q = 1; p = 0.0054). Normalized qRT-PCR data (to *pgk1*) is presented as mean fold changes ± SEM, n = 10 samples/group.). Two-way ANOVA with Sidak's multiple comparison test (a-f). \*p < 0.05, \*\*p < 0.01, \*\*\*p < 0.0001. **g-h**, Cerebellum protein contents determined by Western blot in males, and normalized to GAPDH levels. hnRNPκ/Cfl1:n = 7 C and 7pIKO; PSAP: n = 6 C and 6 pIKO. Multiple unpaired t tests with Holm-Sidak multiple comparison test (\*\*p < 0.01, \*\*\*p < 0.005, \*\*\*\*p < 0.0001). PSAP, prosaposin (RNAseq rank #950/1192 in males; up-regulated; q = 0.5329; p = 0.0359). Cfl1, cofilin-1 (RNAseq rank #27/1192 in males; up-regulated; q = 0.0192; p = 4.03E-05). hnRNPκ, heterogeneous nuclear ribonucleoprotein κ (rank #787; up-regulated; q = 0.509; p = 0.0284). **i**, *mbp*, myelin basic protein. Normalized qRT-PCR data (to *ywhaz*, Tyrosine 3-Monooxygenase/Tryptophan 5-Monooxygenase Activation Protein Zeta). Two-way ANOVA with Sidak's multiple comparison test (n = 11 C and 11 pIKO; \*p < 0.05). See Extended Data Table 1 sheets 1 and 2 for list of ranked DEGs.
<figure><img src="image_63.png" /><figcaption>b</figcaption></figure>

<figure><img src="image_64.png" /><figcaption>c</figcaption></figure>

<figure><img src="image_65.png" /><figcaption>d</figcaption></figure>

<figure><img src="image_66.png" /><figcaption>e</figcaption></figure>

<table><caption>f</caption><thead><tr><th>Diseases and Functions Annotation in Males</th><th>p value</th></tr></thead><tbody><tr><td>atrophy of white matter</td><td>5.34E-03</td></tr><tr><td>abnormal morphology of white matter</td><td>9.47E-03</td></tr><tr><td>long-term potentiation</td><td>1.69E-02</td></tr><tr><td>migration of Purkinje cells</td><td>1.76E-02</td></tr><tr><td>cell-cell contact</td><td>2.62E-02</td></tr><tr><td>apoptosis</td><td>3.15E-02</td></tr></tbody></table>

<table><caption>i</caption><thead><tr><th>Diseases or Functions Annotation in Females</th><th>p value</th></tr></thead><tbody><tr><td>microtubule dynamics</td><td>1.23E-02</td></tr><tr><td>cell death of granule cells</td><td>1.81E-02</td></tr><tr><td>synaptic transmission</td><td>1.94E-02</td></tr><tr><td>neurotransmission</td><td>2.82E-02</td></tr><tr><td>abnormal morphology of white matter</td><td>4.18E-02</td></tr><tr><td>apoptosis</td><td>4.44E-02</td></tr></tbody></table>

<figure><img src="image_67.png" /><figcaption>g</figcaption></figure>

<figure><img src="image_68.png" /></figure>

$$M_{\text{Cereb}} \cap (\text{OL} \cup \text{Myelin}) = 204 \text{ genes}$$

$$F_{\text{Cereb}} \cap (\text{OL} \cup \text{Myelin}) = 68 \text{ genes}$$

<figure><img src="image_69.png" /><figcaption>h</figcaption></figure>
**Extended Data Fig. 3** | Long-term impact of placental ALLO insufficiency on brain transcriptome: highlight on the myelin-associated genes in the cerebellum. **a-b**, Number of differentially expressed genes (DEGs) that are up- or down-regulated in the cerebral cortex, hippocampus and cerebellum of pIKO vs C mice at P30 in males (blue) and females (pink). The majority of DEGs were found in the male cerebellum. **c-d**, MA plots for the cerebellar RNAseq analysis of pIKO vs C mice at P30. M, log ratio; A, mean average. DEGs are represented as colored dots (blue for males and pink for females). **e-f**, Top significant diseases and functions IPA annotations. **g**, Venn diagrams showing the overlaps between cerebellar DEGs in the pIKO males (M_Cereb) and females (F_Cereb), and oligodendrocyte (OL) and myelin transcriptomes³,⁴. **h-i**, Heatmaps of OL-and myelin-related cerebellar DEGs in C and pIKO samples. Genes and samples were hierarchically clustered based on Pearson distance of z-score data and average linkage. n = 3 independent samples/group⁷,⁸.
a Mean FA in cerebellar tracts

<figure><img src="image_70.png" /></figure>

b Mean FA in extra-cerebellar tracts

<figure><img src="image_71.png" /></figure>

<figure><img src="image_72.png" /></figure>
**Extended Data Fig. 4** | Diffusion Tensor Image (DTI) analysis at P30. **a**, Heatmap representing mean fractional anisotropies (FA) in the cerebellum related tracts. Mean FA is increased in the cerebellar WM of plKO males compared to C littermates. Data is presented in a hierarchical tree [region effect: F(25, 1014)=57.36, p < 0.0001; sex effect: F(1,1014)=2.097, p = 0.1479; genotype effect: F(1,1014)=16.99, p < 0.0001; sex × genotype interaction: F(1,1014)=17.92, p < 0.0001; 3-way ANOVA]. **b**, Heatmap representing mean FA in extra-cerebellar tracts. No significant change was evidenced between groups when compared for genotype or sex [region effect: F(8,351)=80.28, p < 0.0001; sex effect: F(1,351)=1.011, p = 0.3152; genotype effect: F(1,351)=0.1618, p = 0.6877; sex × genotype interaction: F(1,351)=2.589, p = 0.1021; 3-way ANOVA]. Males: n = 11 C and 13plKO; Females: n = 8 C and 11 plKO.
<figure><img src="image_73.png" /><figcaption>a</figcaption></figure>

$M\_Cortex \cap (OL \cup Myelin) = 46$ genes

$F\_Cortex \cap (OL \cup Myelin) = 90$ genes

<figure><img src="image_74.png" /><figcaption>b</figcaption></figure>

<figure><img src="image_75.png" /></figure>

OL/Myelin-related DEGs
(M_Cortex)

Extended Data Fig. 5 | See next page for caption.
**Extended Data Fig. 5** | Long-term impact of placental ALLO insufficiency on cerebral cortex transcriptome at P30: highlight on the OL- and myelin-related genes. **a**, Venn diagrams showing the overlaps between cortical DEGs in the plKO males (M_Cortex) and females (F_Cortex), oligodendrocyte (OL) and myelin transcriptomes³,⁴. **b-c**, Heatmaps of OL-and myelin-related DEGs in the cerebral cortex of C and plKO samples. Genes and samples were hierarchically clustered based on Pearson distance of z-score data and average linkage. n = 3 independent samples/group.
<figure><img src="image_76.png" /></figure>

<figure><img src="image_77.png" /></figure>

<figure><img src="image_78.png" /></figure>

<figure><img src="image_79.png" /></figure>

<figure><img src="image_80.png" /></figure>

<figure><img src="image_81.png" /></figure>

<figure><img src="image_82.png" /></figure>

<figure><img src="image_83.png" /></figure>

<figure><img src="image_84.png" /></figure>

<figure><img src="image_85.png" /></figure>

<figure><img src="image_86.png" /></figure>

Extended Data Fig. 6 | See next page for caption.
**Extended Data Fig. 6** | Oligodendrocyte lineage progression during the postnatal period. **a**, Immunofluorescent staining of Olig2 and CC1 in inter-lobule VI-VII WM at P30. Scale bar, 60 µm. **b**, Quantification of Olig2+ and CC1+ cell density in the male cerebellar WM at P30. Data presented as means ± SEM. Two-way ANOVA with Sidak's multiple comparisons test (***p < 0.005, ****p < 0.0001). Olig2: n = 7 C and 8 plKO (p < 0.0001); CC1: n = 6 C and 7plKO (p = 0.0005). **c**, Quantification of BrdU+/Olig2+ cell density in the male cerebellar WM at P30. BrdU (50 mg/kg) was injected to dams at E15.5. Data presented as means ± SEM. Two-tailed unpaired Student's t test with Welch's correction (*p* < 0.05). n = 5 C and 5plKO (p = 0.0107). **d**, Quantification of Olig2+ cell density in the male cerebellar WM at P7 and P15. Data is presented as means ± SEM. Two-way ANOVA with Sidak's multiple comparisons test (***p < 0.005, ****p < 0.00010). P7: n = 3 mice/group (p = 0.0026); P15: n = 6 mice/group (p < 0.0001). The density of Olig2+ cells in the cerebellar WM is significantly increased in plKO males as compared with C littermates. **e**, OPC (PDGFRα + cell) density at P7, P15 and P30. Data is presented as means ± SEM. Two-way ANOVA with Sidak's multiple comparisons test (*p* < 0.05). P7: n = 3 mice/group (p = 0.0435); P15: n = 6 mice/group (p = 0.3644); P30: n = 6 C and 8 plKO (p = 0.8558). The density of OPCs in the cerebellar WM is transiently higher in plKO males as compared with C littermates at P7. **f**, Quantification of mature OL (CC1 + cell) density within the male cerebellar WM at P7 and P15. Data is presented as means ± SEM. Two-way ANOVA with Sidak's multiple comparisons test (***p < 0.005). P7: n = 3 mice/group (p = 0.4951); P15: n = 6 mice/group (p < 0.0004). The density of mature OLs in the cerebellar WM is significantly increased in plKO males as compared with C littermates at P15. **g-i**, Ratios of OPCs (PDGFRα + cells) or mature OLs (CC1 + cells) within the total oligodendrocyte lineage (Olig2 + cells) in the male cerebellar WM across postnatal period. Data is presented as means ± SEM. Two-way ANOVA with Sidak's multiple comparisons test (**p < 0.005). P7: n = 3 mice/group [p = 0.9741 (PDGFRα), p = 0.9393 (CC1)]; P15: n = 6 mice/group [p = 0.002 (PDGFRα), p = 0.001 (CC1)]; P30: n = 7 C and 8 plKO [p = 0.052 (PDGFRα), p = 0.17 (CC1)]. At P15, a transient acceleration of the OL maturation is observed in the male plKO. **j**, Quantification of total OL lineage (Olig2 +) and mature OL (CC1 +) cell densities in the female cerebellar WM at P15. Data is presented as means ± SEM. Two-way ANOVA with Sidak's multiple comparisons test (*p* < 0.01). Olig2 +: n = 7 C and 5 plKO (p = 0.0058); CC1 +: n = 6 C and 5 plKO (p = 0.9984). The density of total OL is increased in the cerebellar WM of plKO females compared to C littermates. **k**, Ratio of mature OLs (CC1 + cells) within the total OL lineage (Olig2 + cells) in the female cerebellar WM at P15. Data is presented as means ± SEM. Two-tailed unpaired Student's t test with Welch's correction (*p* < 0.05). n = 6 C and 5plKO (p = 0.0225). A significantly reduction of OL maturation is observed in the female plKO as compared with littermate C.
<figure><img src="image_87.png" /></figure>

**Extended Data Fig. 7** | Anatomical analysis of cerebellar grey and white matters in C and pIKO at P30. **a–h**, Cerebellar sagittal 40-µm-thick sections through the vermis in males. **a–b**, FluoroMyelin (FM) green staining. **c–d**, NeuN-immunofluorescence (marker for neuronal nuclei). **e–f**, NeuroD1-immunofluorescence (marker for granule cells). **g–h**, Calbindin-immunofluorescence (Purkinje cell marker; in magenta) and DAPI (in blue). **i**, Cerebellar layer volumes were unchanged in pIKO when compared with C mice. Data is presented as means ± SEM from 6 sections/animal. Two-way ANOVA with Sidak's multiple comparisons test. n = 5 C and 7pIKO (WM: p = 0.48; GCL: p = 0.98; ML: p = 0.22). **j**, Purkinje cell linear density in the lobule VI-VII of C and pIKO males at P30. Data is presented as means ± SEM from 6 sections/animal. Two-tailed unpaired Student's t test with Welch's correction. n = 6 C and 7pIKO (p = 0.556). Calb, calbindin; GCL, granule cell layer; ML, molecular layer; WM, white matter. Scale bar, 200 µm.
a

Long steps

<figure><img src="image_88.png" /></figure>

b

Short steps

<figure><img src="image_89.png" /></figure>

c

Speed

<figure><img src="image_90.png" /></figure>

d

Open field

<figure><img src="image_91.png" /></figure>

e

Pre-perturbation time

<figure><img src="image_92.png" /></figure>

f

Post-perturbation time

<figure><img src="image_93.png" /></figure>

g

Missteps

<figure><img src="image_94.png" /></figure>

h

Jumps

<figure><img src="image_95.png" /></figure>

i

Long steps

<figure><img src="image_96.png" /></figure>

j

Short steps

<figure><img src="image_97.png" /></figure>

k

Pre-perturbation time

<figure><img src="image_98.png" /></figure>

l

Post-perturbation time

<figure><img src="image_99.png" /></figure>

m

Missteps

<figure><img src="image_100.png" /></figure>

n

Jumps

<figure><img src="image_101.png" /></figure>

Extended Data Fig. 8 | See next page for caption.
**Extended Data Fig. 8** | Locomotor assessment in the Erasmus ladder at P30. **a-h**, Performance on the Erasmus ladder. **a-b**, The analysis of baseline gait and coordination on days 1-4 showed significantly longer steps in adult male plKO mice than in their control littermates. Data is presented as means ± SEM. Two-way repeated-measures ANOVA with Sidak's multiple comparisons test (\*p < 0.05). n= 8 C and 6 plKO. Long steps=Day 1: p = 0.64; Day 2: p = 0.031; Day 3: p = 0.039; Day 4: p = 0.011. Short steps=Day 1: p = 0.61; Day 2: p = 0.029; Day 3: p = 0.039; Day 4: p = 0.015. **c**, Speed on the ladder was unchanged between C and plKO mice. Data is presented as means ± SEM. Two-way repeated-measures ANOVA with Sidak's multiple comparisons test (n= 8 C and 6 plKO). **d**, Locomotor speed was unchanged in the open-field. Data is presented as means ± SEM. Two-tailed unpaired Student's t test with Welch's correction (n = 6 C and 8 plKO). **e-h**, The associative cerebellar learning was tested with the introduction of a conditioned tone and obstacle on days 5-8. The plKO mice exhibited no motor learning deficit. Data is presented as means ± SEM. Two-way repeated-measures ANOVA with Sidak's multiple comparisons test (n = 8 C and 6 plKO). **i-j**, The analysis of baseline gait and coordination on days 1-4 showed no difference in female plKO mice vs their control littermates. Data is presented as means ± SEM. Two-way repeated-measures ANOVA with Sidak's multiple comparisons test (n = 7 C and 8 plKO). **k-n**, The associative cerebellar learning was tested with the introduction of a conditioned tone and obstacle on days 5-8. The plKO mice exhibited no motor learning deficit. Data is presented as means ± SEM. Two-way repeated-measures ANOVA with Sidak's multiple comparisons test (n = 7 C and 8 plKO).
<figure><img src="image_102.png" /><figcaption>b Rotarod (latency to fall)</figcaption></figure>

<figure><img src="image_103.png" /><figcaption>c Rotarod (learning curve)</figcaption></figure>

<figure><img src="image_104.png" /><figcaption>d Rotarod (learning curve)</figcaption></figure>

<figure><img src="image_105.png" /><figcaption>e Marble burying test</figcaption></figure>

<figure><img src="image_106.png" /></figure>

Extended Data Fig. 9 | See next page for caption.
**Extended Data Fig. 9** | Locomotor and balance behavioral testing at P30. a–d, Accelerating rotarod performance at 4 to 40 rpm (trials 1-6) and 8 to 80 rpm (trials 7-12). a–b, Time to fall off is presented at 4 to 40 rpm and 8 to 80 rpm in males (a) and females (b). Data is presented as means ± SEM. Two-way repeated-measures ANOVA with Sidak's multiple comparisons test (*p* < 0.05, **p** < 0.01, ***p** < 0.005). Males: n = 9 C and 7 pIKO (Trial1: p = 0.4346; Trial2: p = 0.9961; Trial3: p = 0.9918; Trial4: p > 0.9999; Trial5: p = 0.0386; Trial6: p = 0.0016; Trial7: p = 0.9420; Trial8: p = 0.9995; Trial9: p = 0.4957; Trial10: p = 0.9118; Trial11: p = 0.9377; Trial12: p = 0.0495); Females: n = 6 C and 4 pIKO (Any trial: p > 0.98). c–d, Terminal speed of rotation on the accelerating rotarod, presented at 4 to 40 rpm (trials 1-6) and 8 to 80 rpm (trials 7-12), was used to calculate learning rate (linear regression). The slopes from C and pIKO are significantly different in males (c) (F = 11.32, DFd=20, p = 0.0031) but not in females (d). Data is presented as means ± SEM. Two-way repeated-measures ANOVA with Sidak's multiple comparisons test (*p* < 0.05, **p** < 0.01, ***p** < 0.005). Males: n = 9 C and 7 pIKO (Trial1: p > 0.99; Trial2: p = 0.4176; Trial3: p = 0.0328; Trial4: p = 0.6343; Trial5: p = 0.9589; Trial6: p = 0.0328; Trial7: p = 0.6343; Trial8: p = 0.9589; Trial9: p = 0.0372; Trial10: p = 0.0328; Trial11: p = 0.4022; Trial12: p = 0.0045); Females: n = 6 C and 4 pIKO (Any trial: p > 0.9). e, Marble burying test. Graph represents the number of marbles at least 2/3rd buried over 30-min. Data is presented as means ± SEM. Two-tailed unpaired Student's t test with Welch's correction (*p* < 0.05). n = 7 mice/group (p = 0.023).
<figure><img src="image_107.png" /><figcaption>b</figcaption></figure>

<figure><img src="image_108.png" /><figcaption>c</figcaption></figure>

<figure><img src="image_109.png" /></figure>

**Extended Data Fig. 10 | Cerebellar WM and behavioral impairments in male plKO mice can be prevented by pharmacological treatments at E15.5.**

a, Quantification of Olig2+ and CC1+ cell densities in the male cerebellar WM (n=7 C, 8 plKO+Veh and 4 plKO+ALLO). Data presented as means ± SEM. Two-way ANOVA with Sidak's multiple comparisons test (*p<0.05; **p<0.01). Olig2+ = C vs plKO+Veh: p=0.0024; C vs plKO+ALLO: p=0.9982; plKO+Veh vs plKO+ALLO: p=0.041. CC1+ : C vs plKO+Veh: p=0.0014; C vs plKO+ALLO: p=0.9995; plKO+Veh vs plKO+ALLO: p=0.0354.
b, 3-chamber sociability test. Data is presented as means ± SEM. Two-tailed unpaired Student's t test with Welch's correction (***p<0.001).
n=29 C(+Veh) and 14 C + ALLO (p<0.0001).
c, ASD composite severity score. Data presented as means ± SEM. Two-tailed unpaired Student's t test
with Welch's correction (*p<0.05). n=28 C(+Veh) and 13 C + ALLO (p=0.0261).
# Reporting Summary

Nature Research wishes to improve the reproducibility of the work that we publish. This form provides structure for consistency and transparency in reporting. For further information on Nature Research policies, see our Editorial Policies and the Editorial Policy Checklist.

## Statistics

For all statistical analyses, confirm that the following items are present in the figure legend, table legend, main text, or Methods section.

<figure><img src="image_110.png" /></figure>

Our web collection on statistics for biologists contains articles on many of the points above.

## Software and code

Policy information about availability of computer code

**Data collection** Image acquisition on the scanning electron microscope was done using MAPS Software version 3.7 (ThermoFisher). Image acquisition on the confocal microscope was done with LAS X Life Science Microscope software version 2.7 (Leica Microsystems). Image acquisition on the scanning fluorescent microscope was done using cellSens version 2.3 (Olympus). FSL software package 5.0.11 (FMRIB, Oxford, UK) was used to generate fractional anisotropy (FA) and mean diffusivity (MD) maps.

**Data analysis** Data was analyzed using ImageJ 1.53c (NIH), Microsoft Excel 2016 and Graphpad Prism 7. Sequencing data was analyzed using TCC R package (version 1.12.1), Partek Genomics Suite (version 6.6), Ingenuity Pathway Analysis (IPA, Qiagen) and Venny 2.1.

For manuscripts utilizing custom algorithms or software that are central to the research but not yet described in published literature, software must be made available to editors and reviewers. We strongly encourage code deposition in a community repository (e.g. GitHub). See the Nature Research guidelines for submitting code & software for further information.

## Data

Policy information about availability of data

All manuscripts must include a data availability statement. This statement should provide the following information, where applicable:

- Accession codes, unique identifiers, or web links for publicly available datasets

- A list of figures that have associated raw data

- A description of any restrictions on data availability

Raw and processed data are available from the corresponding authors upon request.
# Field-specific reporting

Please select the one below that is the best fit for your research. If you are not sure, read the appropriate sections before making your selection.

<figure><img src="image_111.png" /></figure>

Life sciences

<figure><img src="image_112.png" /></figure>

Behavioural & social sciences

<figure><img src="image_113.png" /></figure>

Ecological, evolutionary & environmental sciences

For a reference copy of the document with all sections, see [nature.com/documents/nr-reporting-summary-flat.pdf](https://nature.com/documents/nr-reporting-summary-flat.pdf)

# Life sciences study design

All studies must disclose on these points even when the disclosure is negative.

**Sample size**

The sample sizes are shown in the legends and chosen to meet or exceed sample sizes typically used in the field (El Khordi et al., 2013; Lerch et al., 2011; Zonouzi et al., 2015). Preliminary experiments were performed when possible to determine requirements for sample size, taking into consideration resources available and ethical, reductionist animal use. Sample size was somewhat limited by the litter sizes and Mendelian ratios of complex genotypes but mice from at least 3 litters were used to overcome this limitation.

**Data exclusions**

- Effect of muscimol treatment on social behavior: one mouse was excluded because it refused to get out of the middle chamber in the 3-chamber equipment.

- Erasmus Ladder: 5 mice were removed from the experiment because they refused to perform the 42 runs in a row per day required for the 2-way Repeated Measures ANOVA analysis.

- Western blot on human tissue: One sample was excluded due to poor tissue quality.

- ROUT methodology was utilized to determine outliers.

**Replication**

Every experiment was performed on animals from at least 3 different litters. Two technical replicates were performed for RT-PCR and Western blot analyses on mouse and human samples. All attempts at replication were successful.

**Randomization**

Sample assignment to groups was not randomized, since it relies on their known genotype (C vs pIKO mice) or gestational age (T and PT infants). Internal controls (housekeeping proteins and genes) were used in the molecular and biochemical experiments.

**Blinding**

All data collection and analyses were conducted blind to group allocation.

# Reporting for specific materials, systems and methods

We require information from authors about some types of materials, experimental systems and methods used in many studies. Here, indicate whether each material, system or method listed is relevant to your study. If you are not sure if a list item applies to your research, read the appropriate section before selecting a response.

## Materials & experimental systems

<table><thead><tr><th>n/a</th><th>Involved in the study</th></tr></thead><tbody><tr><td>☐</td><td>Antibodies</td></tr><tr><td>☒</td><td>Eukaryotic cell lines</td></tr><tr><td>☒</td><td>Palaeontology and archaeology</td></tr><tr><td>☐</td><td>Animals and other organisms</td></tr><tr><td>☐</td><td>Human research participants</td></tr><tr><td>☒</td><td>Clinical data</td></tr><tr><td>☒</td><td>Dual use research of concern</td></tr></tbody></table>

## Methods

<table><thead><tr><th>n/a</th><th>Involved in the study</th></tr></thead><tbody><tr><td>☒</td><td>ChIP-seq</td></tr><tr><td>☒</td><td>Flow cytometry</td></tr><tr><td>☐</td><td>MRI-based neuroimaging</td></tr></tbody></table>

## Antibodies

**Antibodies used**

Antibodies were used as follows for immunohistochemistry: : rabbit anti-MBP (Abcam, Cambridge, MA, USA; #ab40390; 1:500), mouse anti-APC clone CC1 (EMD Millipore, Burlington, NA, USA; #MABC20; 1:500), rabbit anti-Olig2 (Abcam; #ab9610, 1:500), mouse anti-neuronal nuclei (NeuN) (Millipore; #MAB377; 1:500), goat anti-NeuroD1 (RD Systems, Minneapolis, MN, USA; #AF2746; 1:500), rabbit anti-calbindin (Swant, Marly, Switzerland; #CB38, 1:1000), rat anti-mouse PDGFRα (CD140a; BD Biosciences, San Jose, CA, USA; #17-1401-81; 1:500), chicken anti-GFP (Abcam; #ab13970; 1:500), donkey anti-mouse Alexa-488 (Invitrogen; #A-21202), donkey anti-rabbit Alexa-488 (Invitrogen; #A-21206), donkey anti-mouse Alexa-555 (Invitrogen; #A-31570), donkey anti-rabbit Alexa-555 (Invitrogen; #A-31572), donkey anti-mouse Alexa-647 (Invitrogen; #A-31571), and donkey anti-rabbit Alexa-647 (Invitrogen; #A-31573).

Antibodies were used as follows for Western blots: rabbit anti-MBP (Abcam; #ab40390; 1:500), rabbit anti-MAG (Thermo Fisher Scientific; #PA5-79620; 1:500), rabbit anti-MOG (Abcam; #ab32760; 1:500), mouse anti-cofilin-1 (Cfl1) (Santa Cruz Biotechnology, Dallas, TX, USA; #sc-53934; 1:500), mouse anti-prosaposin (PSAP) (Santa Cruz Biotechnology; #sc-390184; 1:500), rabbit anti-hnRNPk (R332) (Cell Signaling Technology; #4675; 1:500) and rabbit anti-GAPDH (Cell Signaling Technology, Danvers, MA, USA; #5174; 1:2000).

**Validation**

These antibodies were commercially available and validated as noted on the manufacturers website. Additional verifications were
Validation

done as follows. The specificity of the secondary fluorescent antibodies was checked by omitting the primary antibodies. The validity of the primary antibodies was tested with the molecular weight of the bands detected by Western blot in mouse and human samples. The antibodies (or clones) are referenced in the Antibodyregistry.org.

## Animals and other organisms

Policy information about studies involving animals; ARRIVE guidelines recommended for reporting animal research

Laboratory animals

The strains used in this study were: Cyp-19a-Cre, Akr1c14-floxed and ROSA26-YFP mice. Males and females were collected at embryonic days: E12.5, E14.5, E15.5, E16.5, E17.5, E19.5 and postnatal days P30.
Mice were housed in a climate controlled facility with a fixed 12/12 hour light cycle with ad libitum access to food and water. All behavioral experiments were conducted between 8 am and 5 pm.

Wild animals

The study did not involve wild animals.

Field-collected samples

The study did not involve samples collected from fields.

Ethics oversight

All the procedures on experimental animals were performed in accordance with the protocols approved by the Institutional Animal Care and Use Committee at Children's National Medical Center (protocol # 30534 (PI: Penn)) and at Columbia University Medical Center (protocols # AC-AABE 6553 (PI: Penn) and AC-AABF 5550 (PI: Yang)).

Note that full information on the approval of the study protocol must also be provided in the manuscript.

## Human research participants

Policy information about studies involving human research participants

Population characteristics

Describe the covariate-relevant population characteristics of the human research participants (e.g. age, gender, genotypic information, past and current diagnosis and treatment categories). If you filled out the behavioural & social sciences study design questions and have nothing to add here, write "See above."

Recruitment

Describe how participants were recruited. Outline any potential self-selection bias or other biases that may be present and how these are likely to impact results.

Ethics oversight

Identify the organization(s) that approved the study protocol.

Note that full information on the approval of the study protocol must also be provided in the manuscript.

## Magnetic resonance imaging

### Experimental design

Design type

MRI was used to examine the structural integrity of white matter tracts in the brains of 4% paraformaldehyde fixed brains.

Design specifications

4% paraformaldehyde perfused brains

Behavioral performance measures

no tasks performed during

### Acquisition

Imaging type(s)

structural/diffusion

Field strength

7T

Sequence & imaging parameters

Parameters for the DTI sequence were as follows: TR of 270 ms, echo train length of 6, first TE of 32 ms, and a TE of 10 ms for the remaining 5 echoes, 1 average, FOV of 14 mm x 14 mm x 25 mm and a matrix size of 180 x 180 x 324 resulting in an image with 78 µm isotropic voxels. Five b=0 s/mm2 images and 30 high b value (b=2147 s/mm2) in 30 different directions were acquired, using the Jones scheme. Total imaging time was approximately 12 hours.

Area of acquisition

whole brain

Diffusion MRI

<figure><img src="image_114.png" /></figure>

Used

<figure><img src="image_115.png" /></figure>

Not used

Parameters

Parameters for the DTI sequence were as follows: TR of 270 ms, echo train length of 6, first TE of 32 ms, and a TE of 10 ms for the remaining 5 echoes, 1 average, FOV of 14 mm x 14 mm x 25 mm and a matrix size of 180 x 180 x 324 resulting in an image with 78 µm isotropic voxels. Five b=0 s/mm2 images and 30 high b value (b=2147 s/mm2) in 30 different directions were acquired, using the Jones scheme. Total imaging time was approximately 12 hours.

### Preprocessing

Preprocessing software

FSL software package (FMRIB, Oxford, UK)
Normalization

To assess any changes to the mouse brains due to genotype and sex, the b=0 s/mm2 images were registered linearly (6 followed by 12 parameter) and non linearly together. A combination of mni_autoreg tools and ANTS (advanced normalization tools) were used to perform the registrations.

Normalization template

A population atlas, representing the average anatomy of the study samples, was created when all scans were resampled with an appropriate transform. The final registration results are the individual images deformed into alignment with one another in an unbiased manner. Warping an available classified MRI mouse brain atlas onto the population average allows for significant differences in the diffusion measures of segmented structures to be calculated.

Noise and artifact removal

N/A

Volume censoring

N/A

Statistical modeling & inference

Model type and settings

This study examined fixed brains, without any behavioral tasks, no statistical modeling was involved.

Effect(s) tested

The analysis of the individual DTI parameters compare intensity differences of FA between genotype and sex, in different ROIs, using 3-way ANOVA.

Specify type of analysis:

<figure><img src="image_116.png" /></figure>

Whole brain

<figure><img src="image_117.png" /></figure>

ROI-based

<figure><img src="image_118.png" /></figure>

Both

Anatomical location(s)
- Cerebellum-related tracts
- Other major cerebral tracts: anterior commissure, posterior commissure, corpus callosum, internal capsule, fornix, mammillothalamic tract, corticospinal tract, cerebral peduncle.

Statistic type for inference
(See Eklund et al. 2016)

Statistical analysis was based on brain region clusters: cerebellar WM tracts and other major (extra-cerebellar) WM tracts

Correction

N/A

Models & analysis

n/a

Involved in the study

<figure><img src="image_119.png" /></figure>

<figure><img src="image_120.png" /></figure>

Functional and/or effective connectivity

<figure><img src="image_121.png" /></figure>

<figure><img src="image_122.png" /></figure>

Graph analysis

<figure><img src="image_123.png" /></figure>

<figure><img src="image_124.png" /></figure>

Multivariate modeling or predictive analysis