<!DOCTYPE article PUBLIC "-//NLM//DTD JATS (Z39.96) Journal Archiving and Interchange DTD v1.1 20151215//EN" "JATS-archivearticle1.dtd"> <article xmlns:ali="http://www.niso.org/schemas/ali/1.0/" xmlns:xlink="http://www.w3.org/1999/xlink" article-type="research-article" dtd-version="1.1"><front><journal-meta><journal-id journal-id-type="nlm-ta">elife</journal-id><journal-id journal-id-type="publisher-id">eLife</journal-id><journal-title-group><journal-title>eLife</journal-title></journal-title-group><issn pub-type="epub" publication-format="electronic">2050-084X</issn><publisher><publisher-name>eLife Sciences Publications, Ltd</publisher-name></publisher></journal-meta><article-meta><article-id pub-id-type="publisher-id">60100</article-id><article-id pub-id-type="doi">10.7554/eLife.60100</article-id><article-categories><subj-group subj-group-type="display-channel"><subject>Research Article</subject></subj-group><subj-group subj-group-type="heading"><subject>Genetics and Genomics</subject></subj-group><subj-group subj-group-type="heading"><subject>Neuroscience</subject></subj-group></article-categories><title-group><article-title>Sequential perturbations to mouse corticogenesis following in utero maternal immune activation</article-title></title-group><contrib-group><contrib contrib-type="author" equal-contrib="yes" id="author-193179"><name><surname>Canales</surname><given-names>Cesar P</given-names></name><contrib-id authenticated="true" contrib-id-type="orcid">https://orcid.org/0000-0003-2505-8367</contrib-id><xref ref-type="aff" rid="aff1">1</xref><xref ref-type="fn" rid="equal-contrib1">†</xref><xref ref-type="other" rid="fund1"/><xref ref-type="fn" rid="con1"/><xref ref-type="fn" rid="conf1"/></contrib><contrib contrib-type="author" equal-contrib="yes" id="author-194423"><name><surname>Estes</surname><given-names>Myka L</given-names></name><xref ref-type="aff" rid="aff1">1</xref><xref ref-type="fn" rid="equal-contrib1">†</xref><xref ref-type="fn" rid="con2"/><xref ref-type="fn" rid="conf1"/></contrib><contrib contrib-type="author" equal-contrib="yes" id="author-194424"><name><surname>Cichewicz</surname><given-names>Karol</given-names></name><contrib-id authenticated="true" contrib-id-type="orcid">https://orcid.org/0000-0001-5926-3663</contrib-id><xref ref-type="aff" rid="aff1">1</xref><xref ref-type="fn" rid="equal-contrib1">†</xref><xref ref-type="fn" rid="con3"/><xref ref-type="fn" rid="conf1"/></contrib><contrib contrib-type="author" id="author-194425"><name><surname>Angara</surname><given-names>Kartik</given-names></name><xref ref-type="aff" rid="aff2">2</xref><xref ref-type="fn" rid="con4"/><xref ref-type="fn" rid="conf1"/></contrib><contrib contrib-type="author" id="author-194426"><name><surname>Aboubechara</surname><given-names>John Paul</given-names></name><xref ref-type="aff" rid="aff1">1</xref><xref ref-type="fn" rid="con5"/><xref ref-type="fn" rid="conf1"/></contrib><contrib contrib-type="author" id="author-194427"><name><surname>Cameron</surname><given-names>Scott</given-names></name><xref ref-type="aff" rid="aff1">1</xref><xref ref-type="fn" rid="con6"/><xref ref-type="fn" rid="conf1"/></contrib><contrib contrib-type="author" id="author-194428"><name><surname>Prendergast</surname><given-names>Kathryn</given-names></name><xref ref-type="aff" rid="aff1">1</xref><xref ref-type="other" rid="fund8"/><xref ref-type="fn" rid="con7"/><xref ref-type="fn" rid="conf1"/></contrib><contrib contrib-type="author" id="author-194429"><name><surname>Su-Feher</surname><given-names>Linda</given-names></name><xref ref-type="aff" rid="aff1">1</xref><xref ref-type="fn" rid="con8"/><xref ref-type="fn" rid="conf1"/></contrib><contrib contrib-type="author" id="author-194430"><name><surname>Zdilar</surname><given-names>Iva</given-names></name><xref ref-type="aff" rid="aff1">1</xref><xref ref-type="fn" rid="con9"/><xref ref-type="fn" rid="conf1"/></contrib><contrib contrib-type="author" id="author-194431"><name><surname>Kreun</surname><given-names>Ellie J</given-names></name><xref ref-type="aff" rid="aff1">1</xref><xref ref-type="fn" rid="con10"/><xref ref-type="fn" rid="conf1"/></contrib><contrib contrib-type="author" id="author-194432"><name><surname>Connolly</surname><given-names>Emma C</given-names></name><xref ref-type="aff" rid="aff1">1</xref><xref ref-type="fn" rid="con11"/><xref ref-type="fn" rid="conf1"/></contrib><contrib contrib-type="author" id="author-194433"><name><surname>Seo</surname><given-names>Jin Myeong</given-names></name><xref ref-type="aff" rid="aff1">1</xref><xref ref-type="fn" rid="con12"/><xref ref-type="fn" rid="conf1"/></contrib><contrib contrib-type="author" id="author-194434"><name><surname>Goon</surname><given-names>Jack B</given-names></name><xref ref-type="aff" rid="aff1">1</xref><xref ref-type="fn" rid="con13"/><xref ref-type="fn" rid="conf1"/></contrib><contrib contrib-type="author" id="author-194435"><name><surname>Farrelly</surname><given-names>Kathleen</given-names></name><xref ref-type="aff" rid="aff1">1</xref><xref ref-type="fn" rid="con14"/><xref ref-type="fn" rid="conf1"/></contrib><contrib contrib-type="author" id="author-194436"><name><surname>Stradleigh</surname><given-names>Tyler W</given-names></name><xref ref-type="aff" rid="aff1">1</xref><xref ref-type="fn" rid="con15"/><xref ref-type="fn" rid="conf1"/></contrib><contrib contrib-type="author" id="author-126525"><name><surname>van der List</surname><given-names>Deborah</given-names></name><xref ref-type="aff" rid="aff1">1</xref><xref ref-type="fn" rid="con16"/><xref ref-type="fn" rid="conf1"/></contrib><contrib contrib-type="author" id="author-194437"><name><surname>Haapanen</surname><given-names>Lori</given-names></name><xref ref-type="aff" rid="aff3">3</xref><xref ref-type="fn" rid="con17"/><xref ref-type="fn" rid="conf1"/></contrib><contrib contrib-type="author" id="author-194438"><name><surname>Van de Water</surname><given-names>Judy</given-names></name><xref ref-type="aff" rid="aff3">3</xref><xref ref-type="fn" rid="con18"/><xref ref-type="fn" rid="conf1"/></contrib><contrib contrib-type="author" id="author-194439"><name><surname>Vogt</surname><given-names>Daniel</given-names></name><xref ref-type="aff" rid="aff2">2</xref><xref ref-type="fn" rid="con19"/><xref ref-type="fn" rid="conf1"/></contrib><contrib contrib-type="author" corresp="yes" equal-contrib="yes" id="author-1986"><name><surname>McAllister</surname><given-names>A Kimberley</given-names></name><contrib-id authenticated="true" contrib-id-type="orcid">https://orcid.org/0000-0001-9177-9889</contrib-id><email>kmcallister@ucdavis.edu</email><xref ref-type="aff" rid="aff1">1</xref><xref ref-type="fn" rid="equal-contrib2">‡</xref><xref ref-type="other" rid="fund2"/><xref ref-type="other" rid="fund6"/><xref ref-type="other" rid="fund4"/><xref ref-type="fn" rid="con20"/><xref ref-type="fn" rid="conf1"/></contrib><contrib contrib-type="author" corresp="yes" equal-contrib="yes" id="author-112175"><name><surname>Nord</surname><given-names>Alex S</given-names></name><contrib-id authenticated="true" contrib-id-type="orcid">https://orcid.org/0000-0003-4259-7514</contrib-id><email>asnord@ucdavis.edu</email><xref ref-type="aff" rid="aff1">1</xref><xref ref-type="fn" rid="equal-contrib2">‡</xref><xref ref-type="other" rid="fund2"/><xref ref-type="other" rid="fund3"/><xref ref-type="other" rid="fund4"/><xref ref-type="other" rid="fund5"/><xref ref-type="fn" rid="con21"/><xref ref-type="fn" rid="conf1"/></contrib><aff id="aff1"><label>1</label><institution>Center for Neuroscience, UC Davis</institution><addr-line><named-content content-type="city">Davis</named-content></addr-line><country>United States</country></aff><aff id="aff2"><label>2</label><institution>Department of Pediatrics & Human Development, Michigan State University</institution><addr-line><named-content content-type="city">East Lansing</named-content></addr-line><country>United States</country></aff><aff id="aff3"><label>3</label><institution>Division of Rheumatology, Allergy and Clinical Immunology, UC Davis</institution><addr-line><named-content content-type="city">Davis</named-content></addr-line><country>United States</country></aff></contrib-group><contrib-group content-type="section"><contrib contrib-type="senior_editor"><name><surname>Cheah</surname><given-names>Kathryn Song Eng</given-names></name><role>Senior Editor</role><aff><institution>The University of Hong Kong</institution><country>Hong Kong</country></aff></contrib><contrib contrib-type="editor"><name><surname>Bhattacharyya</surname><given-names>Anita</given-names></name><role>Reviewing Editor</role><aff><institution>University of Wisconsin, Madison</institution><country>United States</country></aff></contrib></contrib-group><author-notes><fn fn-type="con" id="equal-contrib1"><label>†</label><p>These authors contributed equally to this work</p></fn><fn fn-type="con" id="equal-contrib2"><label>‡</label><p>These authors also contributed equally to this work</p></fn></author-notes><pub-date date-type="publication" publication-format="electronic"><day>05</day><month>03</month><year>2021</year></pub-date><pub-date pub-type="collection"><year>2021</year></pub-date><volume>10</volume><elocation-id>e60100</elocation-id><history><date date-type="received" iso-8601-date="2020-06-16"><day>16</day><month>06</month><year>2020</year></date><date date-type="accepted" iso-8601-date="2021-03-02"><day>02</day><month>03</month><year>2021</year></date></history><permissions><copyright-statement>© 2021, Canales et al</copyright-statement><copyright-year>2021</copyright-year><copyright-holder>Canales et al</copyright-holder><ali:free_to_read/><license xlink:href="http://creativecommons.org/licenses/by/4.0/"><ali:license_ref>http://creativecommons.org/licenses/by/4.0/</ali:license_ref><license-p>This article is distributed under the terms of the <ext-link ext-link-type="uri" xlink:href="http://creativecommons.org/licenses/by/4.0/">Creative Commons Attribution License</ext-link>, which permits unrestricted use and redistribution provided that the original author and source are credited.</license-p></license></permissions><self-uri content-type="pdf" xlink:href="elife-60100-v2.pdf"/><abstract><p>In utero exposure to maternal immune activation (MIA) is an environmental risk factor for neurodevelopmental and neuropsychiatric disorders. Animal models provide an opportunity to identify mechanisms driving neuropathology associated with MIA. We performed time-course transcriptional profiling of mouse cortical development following induced MIA via poly(I:C) injection at E12.5. MIA-driven transcriptional changes were validated via protein analysis, and parallel perturbations to cortical neuroanatomy were identified via imaging. MIA-induced acute upregulation of genes associated with hypoxia, immune signaling, and angiogenesis, by 6 hr following exposure. This acute response was followed by changes in proliferation, neuronal and glial specification, and cortical lamination that emerged at E14.5 and peaked at E17.5. Decreased numbers of proliferative cells in germinal zones and alterations in neuronal and glial populations were identified in the MIA-exposed cortex. Overall, paired transcriptomic and neuroanatomical characterization revealed a sequence of perturbations to corticogenesis driven by mid-gestational MIA.</p></abstract><kwd-group kwd-group-type="author-keywords"><kwd>MIA</kwd><kwd>RNA-Seq</kwd><kwd>transcriptomic</kwd><kwd>brain development</kwd><kwd>WGCNA</kwd></kwd-group><kwd-group kwd-group-type="research-organism"><title>Research organism</title><kwd>Mouse</kwd></kwd-group><funding-group><award-group id="fund1"><funding-source><institution-wrap><institution-id institution-id-type="FundRef">http://dx.doi.org/10.13039/100000025</institution-id><institution>National Institute of Mental Health</institution></institution-wrap></funding-source><award-id>2T32MH073124-16</award-id><principal-award-recipient><name><surname>Canales</surname><given-names>Cesar P</given-names></name></principal-award-recipient></award-group><award-group id="fund2"><funding-source><institution-wrap><institution-id institution-id-type="FundRef">http://dx.doi.org/10.13039/100000025</institution-id><institution>National Institute of Mental Health</institution></institution-wrap></funding-source><award-id>5R21MH116681-02</award-id><principal-award-recipient><name><surname>McAllister</surname><given-names>Kim</given-names></name><name><surname>Nord</surname><given-names>Alex S</given-names></name></principal-award-recipient></award-group><award-group id="fund3"><funding-source><institution-wrap><institution-id institution-id-type="FundRef">http://dx.doi.org/10.13039/100000882</institution-id><institution>Brain Research Foundation</institution></institution-wrap></funding-source><principal-award-recipient><name><surname>Nord</surname><given-names>Alex S</given-names></name></principal-award-recipient></award-group><award-group id="fund4"><funding-source><institution-wrap><institution-id institution-id-type="FundRef">http://dx.doi.org/10.13039/501100001631</institution-id><institution>UCD</institution></institution-wrap></funding-source><principal-award-recipient><name><surname>McAllister</surname><given-names>Kim</given-names></name><name><surname>Nord</surname><given-names>Alex S</given-names></name></principal-award-recipient></award-group><award-group id="fund5"><funding-source><institution-wrap><institution-id institution-id-type="FundRef">http://dx.doi.org/10.13039/100000057</institution-id><institution>National Institute of General Medical Sciences</institution></institution-wrap></funding-source><award-id>GM119831</award-id><principal-award-recipient><name><surname>Nord</surname><given-names>Alex S</given-names></name></principal-award-recipient></award-group><award-group id="fund6"><funding-source><institution-wrap><institution-id institution-id-type="FundRef">http://dx.doi.org/10.13039/100000025</institution-id><institution>National Institute of Mental Health</institution></institution-wrap></funding-source><award-id>P50 MH106438</award-id><principal-award-recipient><name><surname>McAllister</surname><given-names>A Kimberley</given-names></name></principal-award-recipient></award-group><award-group id="fund8"><funding-source><institution-wrap><institution-id institution-id-type="FundRef">http://dx.doi.org/10.13039/100000025</institution-id><institution>National Institute of Mental Health</institution></institution-wrap></funding-source><award-id>T32MH112507</award-id><principal-award-recipient><name><surname>Prendergast</surname><given-names>Kathryn</given-names></name></principal-award-recipient></award-group><funding-statement>The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.</funding-statement></funding-group><custom-meta-group><custom-meta specific-use="meta-only"><meta-name>Author impact statement</meta-name><meta-value>Mid-gestational exposure to maternal immune activation drives a sequence of transcriptional signatures and developmental pathology in embryonic mouse cortex, culminating in altered lamination and cellular lineage specification.</meta-value></custom-meta></custom-meta-group></article-meta></front><body><sec id="s1" sec-type="intro"><title>Introduction</title><p>Epidemiological association between maternal infection and neurodevelopmental disorders (NDDs) has been found for autism spectrum disorder (ASD), schizophrenia (SZ), bipolar disorder (BPD), anxiety, and major depressive disorder (MDD) (<xref ref-type="bibr" rid="bib10">Brown, 2011</xref>; <xref ref-type="bibr" rid="bib26">Estes and McAllister, 2015</xref>; <xref ref-type="bibr" rid="bib61">Parboosing et al., 2013</xref>; <xref ref-type="bibr" rid="bib62">Patterson, 2009</xref>). Indeed, maternal immune activation (MIA) itself is sufficient to produce NDD-relevant outcomes in offspring in animal models (<xref ref-type="bibr" rid="bib12">Canetta et al., 2014</xref>; <xref ref-type="bibr" rid="bib47">Machado et al., 2015</xref>). Offspring from female mice exposed to the toll-like receptor 3 (TLR3) agonist and viral mimic polyinosinic:polycytidylic acid [poly(I:C)], and offspring in maternal immune activation models using the bacterial mimic lipopolysaccharide (LPS) (<xref ref-type="bibr" rid="bib17">Depino, 2015</xref>; <xref ref-type="bibr" rid="bib59">O'Loughlin et al., 2017</xref>; <xref ref-type="bibr" rid="bib71">Simões et al., 2018</xref>), recapitulate neuropathologies and aberrant behaviors seen in offspring born to flu-infected dams (<xref ref-type="bibr" rid="bib69">Shi et al., 2003</xref>; <xref ref-type="bibr" rid="bib72">Smith et al., 2007</xref>). MIA in mice produces behavioral and cognitive abnormalities with relevance across multiple NDD diagnostic boundaries (<xref ref-type="bibr" rid="bib27">Estes and McAllister, 2016</xref>; <xref ref-type="bibr" rid="bib49">Malkova et al., 2012</xref>; <xref ref-type="bibr" rid="bib64">Richetto et al., 2014</xref>).</p><p>MIA-induced outcomes in adult offspring have been well characterized; however, mechanisms of initiation and progression of NDD-related brain pathology remain poorly defined. Poly(I:C) increases proinflammatory cytokines such as IL-1β, IL-6, CXCL1, and TNF-α in the maternal circulation (<xref ref-type="bibr" rid="bib6">Ballendine et al., 2015</xref>; <xref ref-type="bibr" rid="bib16">Dahlgren et al., 2006</xref>; <xref ref-type="bibr" rid="bib32">Hsiao and Patterson, 2011</xref>; <xref ref-type="bibr" rid="bib53">Meyer et al., 2006b</xref>; <xref ref-type="bibr" rid="bib52">Meyer et al., 2006a</xref>; <xref ref-type="bibr" rid="bib72">Smith et al., 2007</xref>). Downstream pathological changes in the brains of offspring have been reported across cell types, from neurons to microglia, and these changes have been linked to perturbations in proliferation, migration, and maturation in the developing brain (<xref ref-type="bibr" rid="bib29">Haida et al., 2019</xref>; <xref ref-type="bibr" rid="bib70">Shin Yim et al., 2017</xref>; <xref ref-type="bibr" rid="bib73">Soumiya et al., 2011</xref>; <xref ref-type="bibr" rid="bib89">Zhang et al., 2019</xref>; <xref ref-type="bibr" rid="bib91">Zhao et al., 2019</xref>). Recent studies have identified acute transcriptomic changes in fetal brain in mice and rats at single time-points (3–4 hr) following MIA, with significant overlap between transcriptome changes in the cortices of MIA-exposed offspring and those altered in the brains of children with ASD (<xref ref-type="bibr" rid="bib46">Lombardo et al., 2018</xref>). Although MIA could cause a dynamic sequence of transcriptional changes in the fetal brain following acute exposure, no integrated picture of MIA-associated transcriptomic pathology exists that maps changes spanning from the initial acute response to MIA through birth. Here, we combined time-course transcriptomics with neuroanatomy to map changes in prenatal mouse cortical development following MIA induced by poly(I:C) injection at E12.5.</p></sec><sec id="s2" sec-type="results"><title>Results</title><p>We used paired transcriptomic and neuroanatomical analyses to map changes in the fetal brain following MIA in a poly(I:C) mouse model (<xref ref-type="fig" rid="fig1">Figure 1</xref>, <xref ref-type="fig" rid="fig1s1">Figure 1—figure supplement 1</xref>, <xref ref-type="fig" rid="fig1s2">Figure 1—figure supplement 2</xref> and <xref ref-type="fig" rid="fig1s3">Figure 1—figure supplement 3</xref>). Induction of MIA via poly(I:C) injection at E12.5 produces relevant phenotypes in mice (<xref ref-type="bibr" rid="bib14">Choi et al., 2016</xref>; <xref ref-type="bibr" rid="bib70">Shin Yim et al., 2017</xref>; <xref ref-type="bibr" rid="bib72">Smith et al., 2007</xref>) and the specific implementation of this model was recently validated in our hands to produce sufficient levels of MIA to mimic viral infection and cause aberrant behavioral outcomes in offspring (<xref ref-type="bibr" rid="bib25">Estes et al., 2020</xref>). The dosage of 30 mg/ml poly(I:C) elicited substantial elevations in maternal serum IL-6 (on average 2200 pg/ml) 4 hr following injection (<xref ref-type="fig" rid="fig1s1">Figure 1—figure supplement 1</xref>). Pregnant mice were injected with saline or poly(I:C) at E12.5, and dorsal telencephalon (E12.5 + 6 hr, E14.5, E17.5, and birth (P0)) was microdissected (<xref ref-type="fig" rid="fig1">Figure 1a</xref>). RNA-seq datasets were generated from male and female embryos from 28 independent litters across control and MIA groups, typically with one to three embryos represented per independent litter (<xref ref-type="supplementary-material" rid="supp20">Supplementary file 20</xref>, <xref ref-type="supplementary-material" rid="supp5">Supplementary file 5</xref>, File source data 1). For control and MIA, 7 <italic>vs.</italic> 7 samples at E12.5, 13 <italic>vs.</italic> 13 samples at E14.5, 6 <italic>vs.</italic> 6 samples at E17.5, and 10 <italic>vs.</italic> 12 samples at P0, were compared in a differential expression (DE) analysis. For transcriptomic analysis, we used a strategy of first defining DE genes at each time-point using the general linear model (GLM) approach implemented in edgeR, followed by mapping DE signatures to systems-level expression patterns via module assignment using weighted gene co-expression network analysis (WGCNA).</p><fig-group><fig id="fig1" position="float"><label>Figure 1.</label><caption><title>Differential expression across embryonic brain development following maternal immune activation via poly(I:C) injection at E12.5.</title><p>(<bold>a</bold>) Schematic representation of MIA model and study design. E: embryonic day, P: postnatal day, IHC: immunohistochemistry, WB: western blot. Poly(I:C) was injected at E12.5. Samples for RNA-seq were collected at E12.5 + 6 hr, E14.5, E17.5, and at birth (P0). IHC and WB analysis were conducted on a separate animal cohort at E17.5 and P0. (<bold>b, c</bold>) Principal component analysis (PCA) of RNA-seq data indicates that developmental age accounts for the majority of variance across samples. Age (DPC) or sex are represented as colored symbols in (<bold>b</bold>) and (<bold>c</bold>) respectively, and poly(I:C) or saline treatment indicated by circles or triangles, respectively (<bold>d</bold>) Heatmap representing relative gene expression changes between control and MIA samples across time-points. Hierarchical clustering by relative fold changes shows stage-specific differential expression signatures (DE genes shown have FDR < 0.05 and log<sub>2</sub> fold change (log<sub>2</sub>FC)>1 or < −1). (<bold>e</bold>) Volcano plots of DE effect size and significance show stage-specific differences in number of DE genes and DE effect size, with the strongest dysregulation at E17.5. Colors represent directionality and statistical significance. (<bold>f</bold>) Numbers of upregulated and downregulated DE genes at FDR < 0.05 and p<0.05 thresholds again show varying DE genes numbers with a peak of dysregulated genes at E17.5. (<bold>g</bold>) Intersection of stage-specific DE genes with the 82 SFARI autism-associated mouse gene orthologs that were expressed at measurable levels in the RNA-seq data. E17.5 DE genes (FDR < 0.05) were enriched for SFARI genes (p=3.9e-04, hypergeometric test). Concentric circles represent developmental time-points. Light red, upregulated DE genes (p<0.05); dark red, upregulated DE genes (FDR < 0.05); light blue, downregulated DE genes (p<0.05); dark blue, downregulated DE genes (FDR < 0.05).</p></caption><graphic mime-subtype="jpeg" mimetype="image" xlink:href="elife-60100.xml.media/fig1.jpg"/></fig><fig id="fig1s1" position="float" specific-use="child-fig"><label>Figure 1—figure supplement 1.</label><caption><title>Validation of responsiveness to poly(I:C) based on dam serum IL-6 levels.</title><p>Injection of poly(I:C) at E12.5 consistently elevates serum IL-6 in dams 4 hr following injection, compared with little response in dams injected with saline.</p></caption><graphic mime-subtype="jpeg" mimetype="image" xlink:href="elife-60100.xml.media/fig1-figsupp1.jpg"/></fig><fig id="fig1s2" position="float" specific-use="child-fig"><label>Figure 1—figure supplement 2.</label><caption><title>Males and females show largely concordant DE signatures following MIA.</title><p>(<bold>a</bold>) Numbers of upregulated and downregulated DE genes with p < 0.05 and FDR < 0.05, in the DE analysis lacking the sex covariate, are similar to the DE model including the sex covariate (<xref ref-type="fig" rid="fig2">Figure 2c</xref>). (<bold>b</bold>) Numbers of upregulated and downregulated DE genes with p < 0.05 and FDR < 0.05, in a DE analysis testing for differential expression between sexes identifies few DE genes, suggesting limited sex dimorphism between samples. MIA conditions and sequencing lanes were set as covariates for the DE model. (<bold>c, d</bold>) Numbers of upregulated and downregulated DE genes with p < 0.05 and FDR < 0.05, for the DE analysis performed on (<bold>c</bold>) males and (<bold>b</bold>) females demonstrate robust DE at E17.5 in both sexes.</p></caption><graphic mime-subtype="jpeg" mimetype="image" xlink:href="elife-60100.xml.media/fig1-figsupp2.jpg"/></fig><fig id="fig1s3" position="float" specific-use="child-fig"><label>Figure 1—figure supplement 3.</label><caption><title>Sex-stratified comparison of DE genes and effect sizes suggest increased DE effect sizes in females.</title><p>(<bold>a</bold>) Correlation of the relative effect sizes (log<sub>2</sub>FC) for genes passing FDR < 0.05 identifies high degree of concordance between MIA responses between sexes. Genes passing FDR threshold in either or both sexes are color coded. (<bold>b</bold>) Venn diagrams illustrating overlap between DE gene sets passing FDR < 0.05 at E14.5 and E17.5, in DE models containing samples from both sexes, males, or females. The DE model comprising individuals of both sexes includes most of the DE signature of each sex. Female DE analysis identifies limited and unique DE signature. (<bold>c</bold>) Examples of E14.5 and E17.5 genes significant in females but not males. At E14.5 genes were selected to pass FDR < 0.05 and logFC > 0.5 or logFC < −0.5 in females, but not logFC > 0.5 or logFC < −0.5 in males. The gene-set selection process was analogous for E17.5 with the addition of FDR < 0.05 condition in males. After gene set filtering, twelve genes with the lowest FDR in females were plotted.</p></caption><graphic mime-subtype="jpeg" mimetype="image" xlink:href="elife-60100.xml.media/fig1-figsupp3.jpg"/></fig></fig-group><p>Principal component analysis (PCA) of MIA and control RNA-seq samples showed developmental age accounts for the majority of variance across samples, with some additional separation between MIA and control groups (<xref ref-type="fig" rid="fig1">Figure 1b–c</xref>). We performed stage-specific differential expression analysis for male and female samples and for both sexes with sex as a covariate. For each time-point, DE genes were defined using a stringent false discovery rate (FDR) < 0.05 threshold and a more inclusive p < 0.05 threshold (<xref ref-type="fig" rid="fig1">Figure 1d,e</xref>, <xref ref-type="supplementary-material" rid="supp1">Supplementary file 1</xref>–<xref ref-type="supplementary-material" rid="supp5">5</xref>). Sex-stratified analysis indicated generally concordant DE differences between male and female samples following MIA, although DE effect sizes appeared stronger in females at E14.5 and E17.5, leading to a larger number of DE genes identified independently in females compared to males (<xref ref-type="fig" rid="fig1s2">Figure 1—figure supplement 2</xref>, <xref ref-type="supplementary-material" rid="supp6">Supplementary files 6</xref>–<xref ref-type="supplementary-material" rid="supp15">15</xref>). Overall, sex-stratified differences in DE genes were generally subtle, as demonstrated by high DE correlation and shared DE gene sets between sexes, with consistent findings across sexes among key DE genes (<xref ref-type="fig" rid="fig1s3">Figure 1—figure supplement 3</xref>). Samples from both male and female offspring were used for overall DE analysis models, with sex included as a covariate.</p><p>For the DE gene set passing the stringent FDR < 0.05 threshold, there were varying numbers of DE genes across time-points (<xref ref-type="fig" rid="fig1">Figure 1d</xref>). The strongest transcriptional signature was observed at E17.5 (2621 up and 3058 down), suggesting dramatic impact of MIA on cortical development at this time-point. P0 represented the subtlest transcriptomic signature, with no genes passing the stringent FDR < 0.05 threshold, and with the p < 0.05 DE genes showing a strong upregulation bias. We tested for overlap between DE RNA-seq genes and high confidence autism-associated genes in the Simons Foundation Autism Research Initiative (SFARI) gene database (<xref ref-type="bibr" rid="bib8">Basu et al., 2009</xref>). The 82 mouse orthologs of SFARI ASD genes that were expressed at measurable levels were significantly enriched among E17.5 DE genes, with 25 upregulated and 19 downregulated DE genes at FDR < 0.05 (p=3.9e-04, hypergeometric test), and 28 upregulated and 28 downregulated DE genes at p < 0.05 (p = 2.8e-06, hypergeometric test) (<xref ref-type="fig" rid="fig1">Figure 1e</xref>, <xref ref-type="supplementary-material" rid="supp16">Supplementary files 16</xref>–<xref ref-type="supplementary-material" rid="supp17">17</xref>). This enrichment during peak transcriptomic dysregulation suggests involvement of ASD-relevant pathways in MIA etiology here and demonstrates general NDD relevance of our model. Detailed results of DE analysis are reported in <xref ref-type="supplementary-material" rid="supp1">Supplementary files 1</xref>–<xref ref-type="supplementary-material" rid="supp15">15</xref> and can be visualized using our interactive online browser.</p><p>We next sought to capture MIA-induced gene regulatory changes at the systems and network level across embryonic cortex development using WGCNA (<xref ref-type="bibr" rid="bib39">Langfelder and Horvath, 2008</xref>; <xref ref-type="bibr" rid="bib90">Zhang and Horvath, 2005</xref>; <xref ref-type="fig" rid="fig2">Figure 2</xref>, <xref ref-type="fig" rid="fig2s1">Figure 2—figure supplement 1</xref>, <xref ref-type="fig" rid="fig2s2">Figure 2—figure supplement 2</xref>). In this approach, genes with correlated expression patterns are assigned into modules, enabling identification of gene sets with shared expression and function in neurodevelopment. WGCNA co-expression modules were arbitrarily named after colors, with Grey reserved for genes with no strongly correlated expression patterns (<xref ref-type="fig" rid="fig2">Figure 2a</xref>). Our initial analysis identified 10 modules (<xref ref-type="fig" rid="fig2s1">Figure 2—figure supplement 1a</xref>), of which correlated modules were combined to produce the final set of five modules and the unassigned Grey set. Modules were tested for association with sample age, MIA versus saline condition, and sex (<xref ref-type="fig" rid="fig2">Figure 2b</xref>). The modules captured dynamic trajectories of neurodevelopmental gene expression (<xref ref-type="fig" rid="fig2">Figure 2c</xref>). The Blue and Turquoise modules represented genes that increased or decreased in expression respectively, during neurodevelopment (Blue: p=9e-34; Turquoise: p=1e-29, <xref ref-type="fig" rid="fig2s1">Figure 2—figure supplement 1b</xref>). The other modules showed more complex patterns of expression. Each WGCNA module included DE genes, with differences in what time-point had the most extensive DE, indicating module-specific timing of perturbation associated with MIA (<xref ref-type="fig" rid="fig2s1">Figure 2—figure supplement 1c</xref>). DE genes were enriched for module-specific GO terms, with GO enrichment findings similar using all module genes and using an alternative rank-based gene set enrichment that is not dependent on FDR or p-value (<xref ref-type="bibr" rid="bib56">Mootha et al., 2003</xref>; <xref ref-type="bibr" rid="bib75">Subramanian et al., 2005</xref>; <xref ref-type="fig" rid="fig2">Figure 2d</xref>, <xref ref-type="supplementary-material" rid="supp18">Supplementary file 18</xref>; <xref ref-type="fig" rid="fig2s2">Figure 2—figure supplement 2</xref>, <xref ref-type="supplementary-material" rid="supp19">Supplementary file 19</xref>).</p><fig-group><fig id="fig2" position="float"><label>Figure 2.</label><caption><title>Weighted gene co-expression network analysis (WGCNA) reveals discrete gene networks perturbed by MIA across embryonic cortical neurodevelopment.</title><p>(<bold>a</bold>) WGCNA cluster dendrogram of time series samples identifies genes that were assigned into co-expression modules. Based on similarities in module gene expression, six of the original 10 modules were grouped into two larger modules, BrRePi: Brown, Red and Pink; YeMaBl: Yellow, Magenta and Black. (<bold>b</bold>) Heatmap of correlation between gene expression modules and experimental traits; age, condition (saline vs poly(I:C)), and sex. Blue and Turquoise modules are strongly associated with age; Green and Grey modules are significantly associated with MIA condition. Numerical values represent signed Pearson’s correlation coefficients, with Student asymptotic p values in brackets. Green represents negative and red represents positive correlation. Color intensity signifies the strength of the correlation. (<bold>c</bold>) Module eigengene expression for MIA and control groups plotted by time-point illustrates expression trajectories across developmental stages, capturing module- and stage-specific differences between MIA and control groups. (<bold>d</bold>) Heatmap of enrichment of DE genes for representative gene ontology biological processes (GO BP) by developmental time-point shows stage- and module-specific transcriptional pathology. Y-axis rows show enrichment for GO BP terms among module-specific DE genes (FDR < 0.05). Heatmap color scale represents relative fold enrichment for the GO terms among DE genes. P0 not shown due to insufficient DE gene numbers and absence of GO BP terms passing enrichment criteria.</p></caption><graphic mime-subtype="jpeg" mimetype="image" xlink:href="elife-60100.xml.media/fig2.jpg"/></fig><fig id="fig2s1" position="float" specific-use="child-fig"><label>Figure 2—figure supplement 1.</label><caption><title>Summary of WGCNA module eigengene expression, association with experimental variables, and overlap with DE genes.</title><p>(<bold>a</bold>) Eigengene trajectories across neurodevelopment of the ten original co-expression modules. (<bold>b</bold>) DE (p < 0.05) gene set enrichment analysis in WGCNA modules across the developmental time-points of the study. Values representing the number of upregulated and downregulated genes are plotted above and below 0, respectively. Asterisks indicate statistically significant enrichment of DE genes (p < 0.05) in a module (p < 0.05; hypergeometric test).</p></caption><graphic mime-subtype="jpeg" mimetype="image" xlink:href="elife-60100.xml.media/fig2-figsupp1.jpg"/></fig><fig id="fig2s2" position="float" specific-use="child-fig"><label>Figure 2—figure supplement 2.</label><caption><title>Stage-specific heatmaps of GO BP enrichment analysis of upregulated and downregulated differentially expressed genes in WGCNA modules.</title><p>(<bold>a</bold>) E12.5 upregulated DE genes (FDR < 0.05). (<bold>b</bold>) E12.5 downregulated DE genes (FDR < 0.05). (<bold>c</bold>) E14.5 upregulated DE genes (FDR < 0.05). (<bold>d</bold>) E14.5 downregulated DE genes (FDR < 0.05). (<bold>e</bold>) E17.5 upregulated DE genes (FDR < 0.05). (<bold>f</bold>) E17.5 downregulated DE genes (FDR < 0.05). (<bold>g</bold>) P0 upregulated DE genes (FDR < 0.05). (<bold>h</bold>) P0 downregulated DE genes (FDR < 0.05). Only modules with more than five DE genes (FDR < 0.05) are shown in the heatmaps.</p></caption><graphic mime-subtype="jpeg" mimetype="image" xlink:href="elife-60100.xml.media/fig2-figsupp2.jpg"/></fig></fig-group><p>The WGCNA-resolved modules enabled mapping of DE signatures to neurodevelopmental processes that were acutely induced by MIA (<xref ref-type="fig" rid="fig3">Figure 3</xref>, <xref ref-type="fig" rid="fig3s1">Figure 3—figure supplement 1</xref>). There was significant association between genes in Grey (p = 0.007) and Green (p = 0.03) modules and MIA treatment, and marginal significance with MIA for the YeMaBl (p = 0.05) module (<xref ref-type="fig" rid="fig2s1">Figure 2—figure supplement 1b</xref>). The acute, initiating signaling pathways, found 6 hr after poly(I:C) injection captured by the Green module and by a subset of genes within the Grey module, show a strong signature of upregulated DE genes from E12.5 to E14.5 (<xref ref-type="fig" rid="fig3">Figure 3a</xref>). Analysis of the GO BP enrichment of upregulated DE genes in these modules identified activation of immune-related pathways, including defense response to virus, as well as angiogenesis and Vascular Endothelial Growth Factor A (VEGFA) signaling (<xref ref-type="fig" rid="fig2">Figures 2c</xref> and <xref ref-type="fig" rid="fig3">3b</xref>, <xref ref-type="fig" rid="fig2s2">Figure 2—figure supplement 2a,c</xref>). We tested for protein-protein association networks among the acute E12.5 DE (FDR < 0.05) genes using STRING (<xref ref-type="bibr" rid="bib77">Szklarczyk et al., 2019</xref>), and found significantly more interactions than expected by chance (observed edges = 161, expected edges = 46, enrichment = 3.5, STRING p < 1.0e-16) (<xref ref-type="fig" rid="fig3">Figure 3c</xref>). This intersecting protein network was associated with metabolism, hypoxia, and stress (<xref ref-type="fig" rid="fig3">Figure 3b</xref>). The network includes a highly interacting core gene set with hub nodes, such as <italic>Vegfa</italic>, that have many spokes and may be signaling factors that direct the MIA response (<xref ref-type="fig" rid="fig3">Figure 3b,c</xref>). Developmental expression profiles for DE genes that were associated with this interaction network, are shown in <xref ref-type="fig" rid="fig3">Figure 3d</xref>, and stratified by sex in <xref ref-type="fig" rid="fig3s1">Figure 3—figure supplement 1a</xref>. These signatures capture a complex but coherent transcriptional response involving hypoxia, immune, metabolic, and angiogenesis pathways that are strongly and transiently induced in fetal brain within 6 hr following mid-gestational MIA.</p><fig-group><fig id="fig3" position="float"><label>Figure 3.</label><caption><title>Green and Grey co-expression modules capture acute transcriptional response at 6 hr and 2 days following MIA exposure.</title><p>(<bold>a</bold>) Violin plots visualize distribution of log<sub>2</sub> fold changes of Green and Grey module gene expression between control and MIA animals. At E12.5, an initial set of genes in the Green and Grey modules are induced and DE in MIA samples. By E14.5, generalized module expression exhibits induction in the MIA samples, with particularly strong change for the Green module, where nearly all genes are upregulated. Genes with expression trajectories shown in (<bold>d</bold>) are labeled. (<bold>b</bold>) Gene set enrichment analysis of GO BP terms significantly enriched among DE genes (FDR < 0.05) at E12.5 and E14.5 in Green and Grey modules showing upregulation of angiogenesis and immune pathways. Representative enriched GO BP terms with p < 0.05 colored by enrichment; gray represents enrichment p > 0.05 (Fisher’s exact test). (<bold>c</bold>) STRING protein–protein interaction network of E12.5 DE genes (FDR < 0.05) colored by annotation to GO BP terms. There were significantly more interactions than expected by chance among these genes. (<bold>d</bold>) RPKM expression plots of genes that are associated with this network show acute DE at E12.5. Stars represent FDR < 0.05.</p></caption><graphic mime-subtype="jpeg" mimetype="image" xlink:href="elife-60100.xml.media/fig3.jpg"/></fig><fig id="fig3s1" position="float" specific-use="child-fig"><label>Figure 3—figure supplement 1.</label><caption><title>Sex-stratified RPKM gene trajectories of genes presented in the manuscript showing concordant effects between males and females.</title><p>(<bold>a</bold>) Vegfa, Flt1, Pdk1, Ldha, Bnip3, Il20rb. (<bold>b</bold>) Mki67, Eomes, Pax6, Sox9, Tbr1, Bcl11b, Satb2, Cux1. (<bold>c</bold>) Gfap, Dlx2, Olig2, Gad1.</p></caption><graphic mime-subtype="jpeg" mimetype="image" xlink:href="elife-60100.xml.media/fig3-figsupp1.jpg"/></fig></fig-group><p>The most pervasive WGCNA-resolved DE signatures following MIA impacted a large proportion of genes in the Blue and Turquoise modules by E14.5 and peaked at E17.5 (<xref ref-type="fig" rid="fig1">Figure 1c</xref>). GO terms related to proliferation and cell cycle were enriched in the Blue module among downregulated DE genes, and terms related to neuronal differentiation and synapses were enriched in the Turquoise module among upregulated DE genes (<xref ref-type="fig" rid="fig2">Figure 2c</xref>). Considering the changes in proliferative and lamination markers associated with Blue and Turquoise genes, we sought to validate such changes at the protein level in independent samples (<xref ref-type="fig" rid="fig4">Figure 4</xref>, <xref ref-type="fig" rid="fig4s1">Figure 4—figure supplement 1</xref>). Four out of five DE genes tested validated with concordant changes at the protein level as determined by western blot (WB) analysis at E17.5, showing either statistical significance or trends in the same direction as observed in our transcriptomic analysis (<xref ref-type="fig" rid="fig4">Figure 4a–b</xref>, <xref ref-type="fig" rid="fig4s1">Figure 4—figure supplement 1</xref>). The DE signatures associated with Blue and Turquoise modules were evidenced by separation between MIA and saline samples in PCA plots at E14.5 and E17.5, with E17.5 MIA samples clustered closer to the P0 samples than saline age-matched counterparts (<xref ref-type="fig" rid="fig1">Figure 1b</xref>). To test for MIA-associated perturbation to temporal full transcriptome expression patterns, we generated a linear regression model using RNA-seq data principal components 1–5 and RNA-seq covariates to model age in saline samples. We then used this model to predict age in MIA samples (<xref ref-type="fig" rid="fig4">Figure 4c</xref>). E17.5 MIA samples exhibited older predicted age values relative to E17.5 saline controls (p = 4.5e-06, T-test), and E14.5 MIA samples showed a similar trend (p = 0.07, T-test). These findings suggest that mid-gestational MIA at E12.5 perturbs the timing of major downstream neurodevelopmental processes at E14.5 and E17.5 in embryonic cerebral cortex.</p><fig-group><fig id="fig4" position="float"><label>Figure 4.</label><caption><title>Altered cortical proliferation and lamination dynamics in fetal offspring five days after induction of MIA.</title><p>(<bold>a</bold>) RPKM trajectories of proliferative and cortical lamination markers that were DE and validated at protein level. (<bold>b</bold>) Quantified protein levels from WB analysis validated DE changes at E17.5. Protein expression data are relative to GAPDH expression (n = 4 control, n = 5 MIA; PAX6 p=0.039, TBR1 p=0.615, CTIP2 p=0.0090, SATB2 p=0.1160, CUX1 p=0.2112; two tailed Student’s t-test). Individual blots are in supplementary data (<xref ref-type="fig" rid="fig4s1">Figure 4—figure supplement 1a–c,h–i</xref>). (<bold>c</bold>) Temporal modeling of the RNA-seq data suggests acceleration of the neurodevelopmental program in MIA animals at E17.5 (p =4.5e-06, Student’s t-test) and similar trend at E14.5 (p = 0.07, Student’s t-test). Actual age (X-axis) vs predicted age (Y-axis) calculated by the linear model. Control samples were used for training the model. Lines connect average values, points depict individual samples and are jittered along the X-axis. (<bold>d</bold>) GO BP categories highlighting module-specific enrichment of E17.5 DE genes in the proliferation (Blue) and neuronal (Turquoise) WGCNA modules. (<bold>e</bold>) Representative images of coronal brain sections from E17.5 saline and MIA animals show reduced progenitor populations (neural stem cells: PAX6, SOX9, SOX2, intermediate progenitors: TBR2, and general proliferative markers: Ki67, PH3) and reduced active proliferation (Ki67:PH3 ratio). Staining for all cell populations is qualitatively decreased in the cortices from MIA offspring. Yellow box indicates the region of the neocortex utilized for Ki67 and Ph3 cell counts, and white boxes indicate regions of the VZ and upper layers of the nCtx utilized for SOX9 and SOX2 cell counts reported in (f-h). Scale bars = 100 µm, (labels: nCtx: neocortex, VZ: ventricular zone, SVZ: sub-ventricular zone, dDG: developing dentate gyrus, GE: ganglionic eminence). (<bold>f</bold>) SOX9 positive cell density in the VZ as well as in the upper layers of the nCtx shows reduced and ectopic SOX9+ cells in the MIA group (n = 3 per condition, SOX9 VZ *p=0.0198, SOX9 nCtx **p=0.0058, two tailed Student’s t-test) (<bold>g</bold>) SOX2 positive cell density counts show a specific reduction in the VZ zone in the MIA group (n = 5 control group, n = 7 MIA group, SOX2 VZ **p=0.0006). (<bold>h</bold>) Ki67 and PH3 cell density quantification along the entire length of the VZ confirms this reduction in proliferative cells in MIA samples (Ki67 n = 5 control group, n = 4 MIA group, Ki67 *p=0.0198; PH3 n = 8 control group, n = 6 MIA group, PH3 **p=0.0003, two tailed Student’s t-test) as well as a reduction in the Ki67/PH3 ratio (n = 5 control group, n = 3 MIA group, *p=0.0005 Chi-Squared).</p></caption><graphic mime-subtype="jpeg" mimetype="image" xlink:href="elife-60100.xml.media/fig4.jpg"/></fig><fig id="fig4s1" position="float" specific-use="child-fig"><label>Figure 4—figure supplement 1.</label><caption><title>Raw western blot and RNA-seq supporting data to proliferation, lamination, and cell-specificity analyses.</title><p>(<bold>a</bold>) Western blot demonstrates reduced PAX6 protein level in tissue extracts from poly(I:C) treated animals at E17.5, (<bold>b</bold>) but not at P0. (<bold>c</bold>) Quantified PAX6 protein levels from western blot analysis at P0 (n = 6 for each control and poly(I:C) groups; p=0.766, two tailed Student’s t-test). (<bold>d–g</bold>) RNA-seq RPKM expression trajectories throughout development of (<bold>d</bold>) Ki67, (<bold>e</bold>) Tbr2, (<bold>f</bold>) Pax6, and (<bold>g</bold>) Sox9. (<bold>h</bold>) Individual western blots detecting TBR1, CTIP2, SATB2 and CUX1 from whole forebrain extracts at E17.5 and (<bold>i</bold>) CTIP2 at P0. (<bold>j</bold>) E17.5 RNA-seq RPKM expression levels of Pax6, Neurog2, Notch1, and (<bold>k</bold>) Gfap and Olig2.</p></caption><graphic mime-subtype="jpeg" mimetype="image" xlink:href="elife-60100.xml.media/fig4-figsupp1.jpg"/></fig></fig-group><p>To characterize the neuroanatomical specificity of Blue and Turquoise dysregulation transcriptional signatures (<xref ref-type="fig" rid="fig4">Figure 4d</xref>), we performed a static immunofluorescent labeling comparison of late neurogenesis, cortical lamination and cell specification at E17.5 in an independent MIA cohort. E17.5 represents the end of the peak of neurogenesis during cortical development, before the neuronal-glial transition that occurs around E18.5. At this stage, it is expected that progenitor markers, such as PAX6, SOX9 and more specifically SOX2, will identify a small population of progenitors that are maintained and retain the ability to produce neurons, as well as subpopulations of proliferative cells commencing gliogenesis (<xref ref-type="bibr" rid="bib50">Martynoga et al., 2012</xref>). PAX6+, SOX9+ and SOX2+ cell distribution patterns showed that progenitors localized correctly around the cortical ventricular zone in coronal brain sections; all analyzed cell populations were significantly reduced in MIA brains compared to controls (<xref ref-type="fig" rid="fig4">Figure 4e–f</xref>). The PAX6 finding in particular is consistent with our WB analysis at E17.5 (PAX6, p=0.039, two tailed Student’s t-test) (<xref ref-type="fig" rid="fig4">Figure 4b</xref>, <xref ref-type="fig" rid="fig4s1">Figure 4—figure supplement 1a</xref>), with protein levels that either resolved or were too subtle to be distinguishable by birth (p=0.766, two tailed Student’s t-test) (<xref ref-type="fig" rid="fig4">Figure 4a</xref>, <xref ref-type="fig" rid="fig4s1">Figure 4—figure supplement 1b–c,i</xref>). MIA brains showed a decrease in SOX9+ cells in the ventricular zone (VZ), which is expressed in progenitors (<xref ref-type="bibr" rid="bib68">Scott et al., 2010</xref>). We also observed increased SOX9+ cells above the VZ in the neocortex, which can remain expressed in a subset of cells, including oligodendrocyte progenitors and maturing astrocytes (<xref ref-type="bibr" rid="bib37">Klum et al., 2018</xref>; <xref ref-type="bibr" rid="bib76">Sun et al., 2017</xref>) (SOX9 VZ, p=0.0198; nCtx, p=0.0058, two tailed Student’s t-test). Reduced expression of SOX2, a progenitor marker that is restricted to the VZ, was also observed in MIA brains (SOX2 VZ, p=0.0006 two tailed Student’s t-test) (<xref ref-type="fig" rid="fig4">Figure 4f</xref>). Overall, these data are consistent with the transcriptomic downregulation signature of proliferating cells in MIA brains.</p><p>To determine if the decrease in progenitor markers might coincide with decreased proliferation, we next labeled cells undergoing proliferation, using the marker KI67 (any stage of cell proliferation) and the Phospho Histone 3 (PH3) marker, which labels a subset of cycling cells in mitosis (M) phase (<xref ref-type="bibr" rid="bib30">Hendzel et al., 1997</xref>; <xref ref-type="bibr" rid="bib84">Walton et al., 2019</xref>). We also examined the post-mitotic intermediate progenitor marker TBR2. Analysis of KI67+ and PH3+ cells showed reduced numbers of labeled cells along with a qualitative reduction in numbers and ectopic positioning of postmitotic intermediate progenitor cells in the MIA group (TBR2+). Reduction in Ki67+ cells was observed in several regions of the developing cortex, including zones defined as S1 and S1DZ (<xref ref-type="bibr" rid="bib14">Choi et al., 2016</xref>; <xref ref-type="bibr" rid="bib70">Shin Yim et al., 2017</xref>), as well as the ganglionic eminences (<xref ref-type="fig" rid="fig4">Figure 4e</xref>). Quantification of Ki67+ cells and PH3+ cells in the neocortex confirmed the MIA-induced reduction in progenitor and mitotic populations at E17.5 (Ki67, p=0.0198; PH3, p=0.0003, two tailed Student’s t-test). The ratio of PH3+ (mitosis) cells to Ki67+ (all stages of cell division) cells was then used to determine whether cell cycle kinetics or the number of actively dividing cells in M phase was different in the MIA-exposed group. Ki67/PH3 ratios indicated a reduced proportion of progenitor cells in M phase in MIA brains at this specific time-point (p=0.0005 Chi-squared) (<xref ref-type="fig" rid="fig4">Figure 4f</xref>). These results provide evidence that the reduced proliferation signatures following MIA that were observed in the transcriptomic profiles are associated with changes at the protein level. Moreover, within the reduced progenitor cell populations at E17.5, the number of cells in the mitosis stage of the cell cycle are decreased to a greater extent.</p><p>MIA has been associated with altered cortical thickness and reported to cause cortical dysplasia in young adult MIA offspring (<xref ref-type="bibr" rid="bib72">Smith et al., 2007</xref>; <xref ref-type="bibr" rid="bib73">Soumiya et al., 2011</xref>; <xref ref-type="bibr" rid="bib82">Tsukada et al., 2015</xref>), with stereotypic dysplasia impacting lamination localized to the somatosensory cortex reported in one MIA model (<xref ref-type="bibr" rid="bib14">Choi et al., 2016</xref>; <xref ref-type="bibr" rid="bib70">Shin Yim et al., 2017</xref>). A number of cortical layer markers were DE at E17.5 in RNA, with validated changes at the protein level (<xref ref-type="fig" rid="fig4">Figure 4a–b</xref>, <xref ref-type="fig" rid="fig3s1">Figure 3—figure supplement 1b</xref>). To examine cortical lamination, distribution of cortical layer-specific markers was analyzed via IHC in the same samples as above (<xref ref-type="fig" rid="fig5">Figure 5</xref>, <xref ref-type="fig" rid="fig4s1">Figure 4—figure supplement figure 1 h-i</xref>).</p><fig id="fig5" position="float"><label>Figure 5.</label><caption><title>MIA impacts cortical lamination during development.</title><p>(<bold>a</bold>) Coronal fetal brain sections from E17.5 saline (control) and MIA offspring show altered lamination pattern. Immunostaining using antibodies to identify deep cortical lamination layers using CTIP2 (green), and TBR1 (red) markers as well as markers of postmitotic neurons, SATB2 and CUX1, in more superficial cortical layers in brains derived from saline (control) or MIA injected dams. Representative images are shown. CTIP2/TBR1 overlap, represented by yellow color, was generally present in control brains, compared to fully laminated deep layers in MIA brains (no overlap illustrated by adjacent green and red rectangles). Boxes indicate magnified areas used for quantification and CTIP2/TBR1 overlap analyses. Scale bars = 100 µm. (<bold>b</bold>) TBR1/CTIP2 lamination overlap occurrence is present in all control brains but absent in the majority of MIA samples at E17.5. Brains from three MIA litters and two control litters were co-stained with TBR1 and CTIP2 cortical markers (as shown in a-b) and analyzed for lamination overlap occurrence (n = 4 control, n = 6 MIA). While all control brains showed lamination overlap between TBR1/CTIP2 layers, absence of TBR1/CTIP2 lamination overlap was observed in four out of six analyzed MIA brains. (<bold>c</bold>) Representative images of SATB2 and TBR1 staining at higher resolution show a change in cell morphology consistent with a difference in migration between MIA and control groups at E17.5. (<bold>d</bold>) Schematic representation illustrates the oval-radially migrating cell body shapes in control brains and the more rounded cell body shape in the MIA brains typical of static cells after migration. (<bold>e</bold>) MIA-induced reduction in TBR1+, CTIP2+, but not in SATB2+ cell counts (n = 4–5 brains per condition, TBR1 p=0.0138; CTIP2 p=0.0119; SATB2 p=0.306). (<bold>f</bold>) Significant decrease in CTIP2 and increase in SATB2 thickness indicates relative expansion of more superficial layers in MIA brains at E17.5 (Layer-specific thickness analysis: n = 4–8 per condition, TBR1 p=0.081; CTIP2 p=0.024; SATB2 p=0.039; CUX1 p=0.430, two tailed Student’s t-test). (<bold>g</bold>) Representative caudal and rostral coronal sections of E17.5 control and poly(I:C) brains stained with DAPI show no obvious cortical dysplasia or gross anatomy alterations. Yellow shading represents neocortical areas considered for thickness measurements. nCtx: neocortex, Th: thalamus, St: striatum, HC: hippocampus, Cp: caudate-putamen. (<bold>h</bold>) Quantification of cortical thickness measured at 30% and 70% distance from the dorsal midline and cortical hemispheric circumference indicates no differences in overall thickness of these neocortical regions (n = 5 control, n = 8 poly(I:C); caudal 30% p=0.5252; caudal 70% p=0.2481; rostral 30% p=0.4459; rostral 70% p=0.4583, two tailed Student’s t-test).</p></caption><graphic mime-subtype="jpeg" mimetype="image" xlink:href="elife-60100.xml.media/fig5.jpg"/></fig><p>Examination of the distribution of the lamination markers TBR1 (layer 6), CTIP2 (layer five and a subset of cells in layer 6), SATB2 (layers 2–4, and a subset of cells in layer 5), and CUX1 (layers 2 and 3), revealed MIA-induced cortical layer-specific alterations in lamination, thickness, and cell density at E17.5 (<xref ref-type="fig" rid="fig5">Figure 5a–b,e–f</xref>). Overall, MIA brains appeared to be at a more mature stage of lamination at E17.5, showing a full separation of cortical deep layers more typical of later ages and in contrast to the partial layer overlap found in control offspring at E17.5 (<xref ref-type="fig" rid="fig5">Figure 5a–b</xref>). Similarly, SATB2 and TBR1 co-staining suggested differences in lamination and migration state between the two groups. Oval cell body shapes were observed, which resembled radially migrating cells (<xref ref-type="bibr" rid="bib15">Cooper, 2013</xref>) were observed in both SATB2 and TBR1 cells in control tissue, whereas these cells showed a more rounded/post-migratory cell body shape in the MIA brains (<xref ref-type="fig" rid="fig5">Figure 5c–d</xref>). The number of TBR1+ and CTIP2+ cells, but not SATB2+ cells, were also significantly reduced in the neocortex from MIA offspring (TBR1+ cells, p=0.0138; CTIP2+ cells, p=0.0119; SATB2+ cells, p=0.306, two tailed Student t-test) (<xref ref-type="fig" rid="fig5">Figure 5e</xref>). The thickness of these layers and also a more external layer represented by CUX1, were measured by comparing the radial spread of stained cells within the neocortex. We found a significant reduction in CTIP2 layer thickness and a similar trend in TBR1 layer thickness in the MIA brains. In contrast, SATB2 showed an MIA-induced significant increase in layer thickness with no change in the more external layer represented by CUX1 (TBR1 thickness, p=0.081; CTIP2 thickness, 0.024; SATB2 thickness, p=0.039; CUX1, p=0.430, two tailed Student’s t-test). (<xref ref-type="fig" rid="fig5">Figure 5f</xref>). These findings are consistent with individual trajectories of DE genes and WB analyses (<xref ref-type="fig" rid="fig4">Figure 4a–b</xref>, <xref ref-type="fig" rid="fig4s1">Figure 4—figure supplement 1h</xref>). Finally, to assess cortical gross anatomy, we looked for MIA-associated dysplasic alterations across the brain (<xref ref-type="bibr" rid="bib14">Choi et al., 2016</xref>) and measured the thickness of the developing somatosensory cortex at two anteroposterior regions at E17.5. In contrast to previous reports (<xref ref-type="bibr" rid="bib14">Choi et al., 2016</xref>), we did not find obvious cortical dysplasia or gross alterations in brain anatomy in the analyzed rostral and caudal regions, including somatosensory cortex and S1DZ (<xref ref-type="fig" rid="fig5">Figure 5e–f</xref>).</p><p>To investigate other cell populations that may be impacted by MIA during development, we analyzed astrocytes, oligodendrocytes and GABAergic interneuron cell populations (<xref ref-type="fig" rid="fig6">Figure 6</xref>, <xref ref-type="fig" rid="fig4s1">Figure 4—figure supplement 1k</xref>). Concordant with RNA-seq results, the number of cells expressing GFAP, a protein found in radial glia and a subset of astrocytes (<xref ref-type="bibr" rid="bib23">Eng et al., 2000</xref>), was increased in MIA offspring, mapping specifically to cells that seemed to emerge ventral from the neocortical VZ (<xref ref-type="fig" rid="fig6">Figure 6a–b,q</xref>, <xref ref-type="fig" rid="fig4s1">Figure 4—figure supplement 1k</xref>). Since astrocytes also express SOX9 (<xref ref-type="bibr" rid="bib37">Klum et al., 2018</xref>; <xref ref-type="bibr" rid="bib76">Sun et al., 2017</xref>) and radial glia also express PAX6, which is decreased by MIA, this finding is consistent with the possibility that the increased SOX9+ cells in the dorsal regions of the neocortex, described earlier, may be indicative of early differentiated astrocytes (<xref ref-type="fig" rid="fig4">Figure 4f</xref>). OLIG2, which is expressed in oligodendrocyte progenitors and mature cells as well as transiently in astrocyte progenitors (<xref ref-type="bibr" rid="bib18">Dimou et al., 2008</xref>), was also upregulated at E17.5 in the MIA group. OLIG2+ cells were increased in MIA brains compared to controls (n = 3 per group, p=0.0134, Student’s t-test) (<xref ref-type="fig" rid="fig6">Figure 6c–d,k</xref>, <xref ref-type="fig" rid="fig4s1">Figure 4—figure supplement 1k</xref>). PDGFRα, a marker of early developing oligodendrocytes, was also increased in MIA brains compared to controls (n = 7 per group, zone 1 p<0.0001, zone 2 p=0.0063, zone 3 p=0.0125, Student’s t-test) (<xref ref-type="fig" rid="fig6">Figure 6e–f,l</xref>). Thus, these experiments provide evidence that MIA leads to increased numbers of additional cell types (putative astrocytes and oligodendrocytes) that are typically born at later ages during corticogenesis, reinforcing findings of accelerated cell development that is supported by an increased representation of later-born cell types at E17.5 in brains exposed to MIA.</p><fig id="fig6" position="float"><label>Figure 6.</label><caption><title>MIA alters development of neuronal and glial cell types in fetal offspring at E17.5.</title><p>(<bold>a–j</bold>) Coronal fetal brain sections from three independent E17.5 saline litters (n = 4 control) and four independent MIA litters (n = 6 MIA) were immunostained for cell-type markers. (<bold>a,b</bold>) Altered distribution and numbers of GFAP+ cells (radial glia and astrocytes) in MIA samples. (<bold>c-j</bold>) Reduced sub-population of interneurons and oligodendrocytes using OLIG2, PDGFRα, DLX2, and GAD67. (<bold>k-p</bold>) Quantification of the cell density in (<bold>c-j</bold>). Overall MIA-associated increase in OLIG2+ cells (<bold>k</bold>), with further specific increased PDGFRα+ cells representing oligodendrocytes (<bold>l</bold>), and no change in overall DLX2+ interneurons (<bold>m,o</bold>), but increased overall number of DLX2+/GAD67+ interneurons (<bold>n,p</bold>) are observed in MIA neocortex. Student’s t-test: OLIG2 p=0.0134; DLX2+/GAD67+ cells p=0.0004. Bin quantification analyses of PDGFRα+, DLX2+, and GAD67+ immunoreactive cells in zones 1–3 (Z1 to Z3) indicate that changes are not restricted to a particular zones of the neocortex. Arrows denote co-labeled cells. Scale bars in (a-h) = 100 µm and in (i-j) = 30 µm. Abbreviations: neocortex (nCtx); marginal zone (mz); ventricular zone (VZ). (<bold>q</bold>) RPKM trajectories of the cell identity markers tested in (<bold>a-k</bold>) across the developmental time-points of the study show parallel differences in mRNA expression of cell-type markers following MIA induction. <italic>Gfap</italic> and <italic>Olig2</italic> mRNA was significantly upregulated in MIA samples at E17.5. <italic>Dlx2</italic> was significantly downregulated at E17.5, with no significant changes observed for <italic>Gad1</italic> (encodes GAD67) or for <italic>Pdgfra</italic> (not shown).</p></caption><graphic mime-subtype="jpeg" mimetype="image" xlink:href="elife-60100.xml.media/fig6.jpg"/></fig><p>Finally, we examined DLX2 and GAD67, markers of committed and maturing cortical interneurons. Although <italic>Dlx2</italic> expression was downregulated in the MIA group from bulk tissue dissections at E17.5, the number of cells positive for DLX2 protein in the neocortex were comparable when assessed via IHC in the MIA and control groups (n = 3 per group, total DLX2 p=0.5494; by zones DLX2, zone 1 P=0.9848, zone 2 P=0.6903, zone 3 P=0.2714, Student’s t-test) (<xref ref-type="fig" rid="fig6">Figure 6g–h and m,o</xref>). In contrast, GAD protein in this DLX2+ subpopulation, assessed by the pan-GABAergic marker GAD67+, was increased across the neocortex of MIA brains, and was not restricted to any specific zones of the neocortex (n = 3 per group, total DLX2+/GAD67 p=0.0004; by zones DLX2+/GAD67, zone 1 P=0.2277, zone 2 P=0.1352, zone 3 P=0.2862, Student’s t-test) (<xref ref-type="fig" rid="fig6">Figure 6i–j and n,p</xref>). Since GAD67 expression increases in developing cortical interneurons, the decrease in <italic>Dlx2</italic> transcript expression but no difference in cells expressing DLX2 protein could indicate increased interneuron maturation but no change in interneuron numbers at this age. While neuroanatomical characterization here was limited to E17.5 and interrogated a subset of genes and processes of interest, our IHC findings correlate with transcriptional profiles, providing validation of MIA-associated DE in developing cortex.</p></sec><sec id="s3" sec-type="discussion"><title>Discussion</title><p>Studies leveraging MIA models in mice have demonstrated the role of maternal immune signaling as a risk factor for cellular, anatomical, and behavioral outcomes in offspring. Nonetheless, a systems-level perspective of the early neurodevelopmental trajectory of pathology associated with MIA has been lacking. Our results offer a time course resolved map of molecular targets aberrantly transcribed in developing cerebral cortex following MIA. MIA-associated transcriptional signatures were robust, changed greatly across the surveyed time- points, and revealed novel pathways and processes involved in pathology. Among MIA-associated DE up-regulated and down-regulated genes, there was strong enrichment for ASD-relevant loci, indicating the general translational relevance of our MIA model to ASD. Both male and female samples were examined, and we did not find evidence for strong sex-specific DE signatures or impacts on neuroanatomy. Elucidating the trajectory of transcriptional changes across embryonic cortical development following mid-gestational MIA represents a valuable step toward understanding the progression of MIA-induced pathology in the developing brain of offspring.</p><p>Our results suggest two waves of transcriptional pathology following MIA. First, we found initial induction of immune, metabolic, and angiogenesis gene expression modules specific to MIA exposure and largely restricted to the days following MIA. There is overlap in the acute induction signatures identified here and acute changes reported in other mid-gestational MIA models, including those using a range of injection gestational ages and rat models of LPS (<xref ref-type="bibr" rid="bib13">Carpentier et al., 2013</xref>; <xref ref-type="bibr" rid="bib54">Meyer, 2014</xref>; <xref ref-type="bibr" rid="bib52">Meyer et al., 2006a</xref>; <xref ref-type="bibr" rid="bib59">O'Loughlin et al., 2017</xref>; <xref ref-type="bibr" rid="bib64">Richetto et al., 2014</xref>; <xref ref-type="bibr" rid="bib72">Smith et al., 2007</xref>; <xref ref-type="bibr" rid="bib87">Wu et al., 2017</xref>). Previous studies linked maternal immune challenge, prenatal brain hypoxia, and impacts on proliferation (<xref ref-type="bibr" rid="bib74">Stolp et al., 2011</xref>), neurogenesis, and cortical lamination (<xref ref-type="bibr" rid="bib13">Carpentier et al., 2013</xref>). Thus, it is possible that the acute signatures identified here and mapped to the Green and YeMaBl modules as well as a subset of genes in the unassigned Grey module, capture a generalized response in fetal brain across mid-gestational MIA models. Notably, induction of hypoxia, cellular metabolism and neuroimmune transcriptional signatures, particularly captured here by the Green module, persist at reduced levels at least through P0, and may indicate lasting perturbation to these pathways. Significant association of the Grey unassigned module with the MIA phenotype was driven by a subset of genes that showed strong acute induction at E12.5 and E14.5. By altering WGCNA parameters, we could assign the Grey module genes to co-expression modules but those alterations increased splitting, resulting in a high number of modules and reduced module coherence without clarifying additional MIA-specific modules made up by the relevant Grey module MIA-associated genes. Thus, we opted for reduced splitting to capture major co-expression trends.</p><p>Identifying the early biological processes triggered by MIA are of critical interest for understanding how MIA initiates pathology. Our analysis identified potential drivers of these induced signaling pathways, including targets of substantial interest such as <italic>Vegfa,</italic> that link hypoxia, neurogenesis and cortical angiogenesis (<xref ref-type="bibr" rid="bib4">Argaw et al., 2012</xref>; <xref ref-type="bibr" rid="bib78">Takahashi and Shibuya, 2005</xref>). <italic>Vegfa</italic> was upregulated in MIA-exposed animal at all time-points and was among the strongest DE genes at E12.5 and P0. <italic>Vegfa</italic> is critical for cortical development promoting neurogenesis and neurite outgrowth, as well as angiogenesis (<xref ref-type="bibr" rid="bib9">Bayless et al., 2016</xref>; <xref ref-type="bibr" rid="bib48">Mackenzie and Ruhrberg, 2012</xref>; <xref ref-type="bibr" rid="bib66">Rosenstein et al., 2010</xref>). Upregulation of <italic>Vegfa</italic> is found in several NDDs and psychiatric disorders, including SZ and depression, and was also reported in a different MIA mouse model (<xref ref-type="bibr" rid="bib5">Arrode-Brusés and Brusés, 2012</xref>; <xref ref-type="bibr" rid="bib31">Hohman et al., 2015</xref>; <xref ref-type="bibr" rid="bib33">Iga et al., 2007</xref>; <xref ref-type="bibr" rid="bib34">Isung et al., 2012</xref>; <xref ref-type="bibr" rid="bib36">Kahl et al., 2009</xref>; <xref ref-type="bibr" rid="bib41">Lee and Kim, 2012</xref>; <xref ref-type="bibr" rid="bib45">Lizano et al., 2018</xref>; <xref ref-type="bibr" rid="bib88">Xie et al., 2017</xref>). Induction of angiogenesis in the developing cortex also independently causes a transition from proliferative state to neuronal maturation (<xref ref-type="bibr" rid="bib9">Bayless et al., 2016</xref>; <xref ref-type="bibr" rid="bib35">Javaherian and Kriegstein, 2009</xref>; <xref ref-type="bibr" rid="bib79">Tan et al., 2016</xref>). Other genes of interest include <italic>Pdk1</italic> (pyruvate dehydrogenase kinase 1), which regulates glucose metabolism (<xref ref-type="bibr" rid="bib58">Newington et al., 2012</xref>), <italic>Ldha</italic> (lactate dehydrogenase A), which facilitates glycolysis and participates in brain angiogenesis (<xref ref-type="bibr" rid="bib43">Lin et al., 2018</xref>), <italic>Bnip3</italic> (BCL2 interacting protein 3), which is involved in cell death and mitochondrial autophagy (<xref ref-type="bibr" rid="bib44">Liu and Frazier, 2015</xref>), and <italic>Il20rb</italic> (interleukin 20 receptor subunit beta), a receptor subunit for pro-inflammatory interleukins IL-19, IL-20 and IL-24 (<xref ref-type="bibr" rid="bib1">Akdis et al., 2016</xref>). The early perturbed pathways observable within hours of MIA induction may be critical initiating steps leading to the altered general timing of fetal brain development described here. Our findings expand understanding of the molecular basis of these MIA-induced initial perturbations and reveal potentially key drivers of these changes.</p><p>The second component of the neurodevelopmental MIA signature defined here was perturbation to normal patterns of both proliferative and maturational neurodevelopmental modules that manifest several days after MIA induction. Initially, at six hours after MIA, there is a shift in expression of differentiation and neuronal maturation processes (e.g. synaptogenesis) to increased proliferative processes (e.g. positive regulation of cell cycle). This pattern parallels results described for the acute MIA response from other mouse and rat MIA models, where similar pathways were down- or up-regulated in fetal brain in the hours following MIA (<xref ref-type="bibr" rid="bib55">Money et al., 2018</xref>; <xref ref-type="bibr" rid="bib60">Oskvig et al., 2012</xref>). At this developmental stage, our data suggest that MIA induces an acute expenditure of the neocortical, and potentially other, progenitor pools, which in turn leads to an expedited maturation of cells that would normally be generated later during development. Specifically, by E17.5, multiple progenitor and cell cycle markers were significantly decreased in the MIA brains while concordantly, we observed an increase in maturational markers and laminar position that resembled a more developed cortex. Finally, astrocytes, oligodendrocytes, GABAergic interneurons and radially migrating excitatory neuron precursors were effected, suggesting an unbiased impact on progenitors of multiple brain cell types. Although these findings are indicative of overall decreased neurogenesis by E17.5, specific birth-dating and time-course experiments are needed to test whether progenitor cells exit the cell cycle differently and if overall neurogenesis is perturbed beyond the E17.5 time-point reported here. Nonetheless, our time course RNA-seq design revealed that the transcriptional changes following MIA can be separated into discrete biological and temporal expression modules that support altered proliferation and maturation of cell types following MIA.</p><p>Our findings of disrupted cortical lamination following the acute MIA signature are consistent with, but not identical to previous studies showing cortical dysplasia associated with MIA-induced aberrant behavior in offspring (<xref ref-type="bibr" rid="bib14">Choi et al., 2016</xref>; <xref ref-type="bibr" rid="bib70">Shin Yim et al., 2017</xref>). In those studies, dysplasia resolved at later postnatal stages but defects in lamination persisted into adulthood. Although we did not find evidence of aberrant cortical patches in the somatosensory cortex, secondary motor cortex or S1 and S1DZ, as previously reported, our results do show altered cortical lamination at E17.5. Our findings suggest that perturbation to neuronal migration, potentially due to differences in an early exit of pyramidal neuron precursors from the cell cycle that led to cells further along in their radial migration/lamination, developmental progression, at least in part contribute to MIA-associated lamination differences identified at E17.5 here. Thus, results from multiple studies using MIA models to date are consistent in showing a detrimental impact of MIA on cortical patterning; however, we also acknowledge that our and other studies may have differential variables, including onset of MIA, dose, mouse background, and endpoint measurements, which will be important factors to assess in future studies. Finally, bulk tissue RNA-seq is limited in cell type-specific resolution and signatures could additionally be in part driven by cortex-adjacent structures, such as primordial hippocampus, that may have also been captured at the youngest age. Future work using single-cell approaches is needed to resolve the impacted cell types and compare transcriptional pathology following MIA.</p><p>While significantly subtler in overall effects compared to earlier time-points, we also see evidence for altered neuronal transcription and lasting immune activation at P0. By P0, these signatures were reduced, which suggests a transient or attenuating effect on cortical patterning. It is also possible that biological variability in neonatal mice or technical variability impacted detection of gene expression changes in these P0 samples. Nonetheless, the subtle changes in gene expression at P0 are consistent with deficits in synapse formation that require changes in expression of immune molecules in neurons from newborn MIA offspring (<xref ref-type="bibr" rid="bib22">Elmer et al., 2013</xref>).</p><p>Consistent with links between NDDs and developmental timing, atypical neurodevelopmental gene expression networks have been reported in a cortical neurodevelopmental model using autism-patient-derived induced pluripotent stem cells (<xref ref-type="bibr" rid="bib9">Bayless et al., 2016</xref>; <xref ref-type="bibr" rid="bib67">Schafer et al., 2019</xref>). The majority of the ASD-relevant genes that were DE following MIA (<xref ref-type="fig" rid="fig1">Figure 1e</xref>) were associated with perturbation to the modules associated with neurodevelopmental processes. These ASD-relevant genes map to functional and biological pathways generally associated with ASD and NDDs, including gene regulation, neurogenesis, and synapse development, structure, and function. The broad pathway overlap and presence of up- and down-regulation among DE ASD-relevant genes supports their generalized relevance. These overlapping genes represent targets of interest for future studies investigating parallel mechanisms between genetic and MIA-associated NDD risk. Our findings suggest some aspects of the transition from neurogenesis to astrogenesis appear to occur at an earlier stage following MIA based on similar signatures of SOX9+ and GFAP+ cells reported in literature (<xref ref-type="bibr" rid="bib7">Bansod et al., 2017</xref>), and increased OLIG2+ cells and more mature GAD67+ DLX2+ interneurons in the cortex following MIA. These findings are additionally consistent with recent reports using poly(I:C) MIA models that identified perturbed GABAergic interneuron development (<xref ref-type="bibr" rid="bib81">Thion et al., 2019</xref>), and maturation of GAD+ neuroblast after birth (<xref ref-type="bibr" rid="bib83">Vasistha et al., 2020</xref>). Future work is needed to determine if and how MIA-associated dynamic transcriptional changes manifest postnatally and whether perturbation to embryonic development impacts cortical cytoarchitecture and function at later postnatal stages.</p><p>Altogether, our analyses indicate that MIA drives acute and lasting transcriptional changes that impact specific cell populations during forebrain development, disrupting normal developmental timing and causing increased numbers of astrocytes, oligodendrocytes, and altered molecular properties of GABAergic interneurons. The timing of MIA during gestation is known to cause a range of phenotypes in offspring, and thus the findings here regarding mid-gestational MIA at E12.5 may not generalize to other models with different timing of MIA induction, especially for later gestational MIA after neurogenesis is largely complete. Further, the strong embryonic signatures we observe were generally consistent across sexes, raising the question of whether secondary insults or sex-specific protective effects contribute to some reports of sex differences in pathology following MIA (<xref ref-type="bibr" rid="bib29">Haida et al., 2019</xref>; <xref ref-type="bibr" rid="bib57">Naviaux et al., 2013</xref>; <xref ref-type="bibr" rid="bib72">Smith et al., 2007</xref>). Nonetheless, the overlapping results between our study and published single time-point studies of the acute MIA fetal brain transcriptional response suggest some consistency across bacterial and viral mimics, and rat and mouse models. Our study substantially expands previous work, revealing new molecular pathways and providing a temporal map of MIA-induced changes across fetal neurodevelopment that will be relevant to understanding the pathology of, and informing future treatments for, NDDs.</p></sec><sec id="s4" sec-type="materials|methods"><title>Materials and methods</title><table-wrap id="keyresource" position="anchor"><label>Key resources table</label><table frame="hsides" rules="groups"><thead><tr><th valign="top">Reagent type (species) or resource</th><th valign="top">Designation</th><th valign="top">Source or reference</th><th valign="top">Identifiers</th><th valign="top">Additional information</th></tr></thead><tbody><tr><td valign="top">Chemical compound, drug</td><td valign="top">poly(I:C) dsRNA</td><td valign="top">Sigma Aldrich</td><td valign="top">P0913</td><td valign="top">Lot # 016M1451V</td></tr><tr><td valign="top">Antibody</td><td valign="top">Rabbit polyclonal anti PAX-6</td><td valign="top">Covance</td><td valign="top">PRB-278P-100</td><td valign="top">Dilution: (1:250) for IHC; <break/>(1:3000) for WB</td></tr><tr><td valign="top">Antibody</td><td valign="top">Rabbit polyclonal anti-TBR1</td><td valign="top">Abcam</td><td valign="top">ab31940</td><td valign="top">Dilution: (1:500) for IHC; <break/>(1:2000) WB</td></tr><tr><td valign="top">Antibody</td><td valign="top">Rat monoclonal anti-CTIP2</td><td valign="top">Abcam</td><td valign="top">ab18465</td><td valign="top">Dilution (1:250) for IHC; <break/>(1:1000) for WB</td></tr><tr><td valign="top">Antibody</td><td valign="top">Rabbit polyclonal anti-CUX1</td><td valign="top">Abclonal</td><td valign="top">A2213</td><td valign="top">Dilution (1:200) for IHC; <break/>(1:1000) for WB</td></tr><tr><td valign="top">Antibody</td><td valign="top">Mouse monoclonal anti-SATB2</td><td valign="top">Abcam</td><td valign="top">Ab51502</td><td valign="top">Dilution (1:500) for IHC; <break/>(1:2000) for WB</td></tr><tr><td valign="top">Antibody</td><td valign="top">Rabbit monoclonal anti-KI67</td><td valign="top">Cell Signaling</td><td valign="top">12202</td><td valign="top">Dilution (1:500) for IHC</td></tr><tr><td valign="top">Antibody</td><td valign="top">Goat polyclonal anti-SOX9</td><td valign="top">R and D Systems</td><td valign="top">AF3075</td><td valign="top">Dilution (1:500) for IHC</td></tr><tr><td valign="top">Antibody</td><td valign="top">Rabbit polyclonal anti-PH3</td><td valign="top">Cell Signaling</td><td valign="top">9701</td><td valign="top">Dilution (1:500) for IHC</td></tr><tr><td valign="top">Antibody</td><td valign="top">Rat monoclonal anti-TBR2</td><td valign="top">Thermo Fisher Scientific</td><td valign="top">14-4875-82</td><td valign="top">Dilution (1:500) for IHC</td></tr><tr><td valign="top">Antibody</td><td valign="top">anti-DLX2</td><td valign="top">John Rubenstein Lab</td><td valign="top">N/A</td><td valign="top">Dilution (1:200) for IHC</td></tr><tr><td valign="top">Antibody</td><td valign="top">Rabbit polyclonal anti-GFAP</td><td valign="top">Agilent Dako</td><td valign="top">Z0334</td><td valign="top">Dilution (1:250) for IHC</td></tr><tr><td valign="top">Antibody</td><td valign="top">Rabbit polyclonal anti-OLIG2</td><td valign="top">Millipore Sigma</td><td valign="top">AB9610</td><td valign="top">Dilution (1:500) for IHC</td></tr></tbody></table></table-wrap><sec id="s4-1"><title>Animal care and use</title><p>All studies were conducted in compliance with NIH guidelines and approved protocols from the University of California Davis Animal Care and Use Committee. C57BL/6N females (Charles River, Kingston, NY), were bred in house and maintained on a 12:12 hr light dark cycle at 20 ± 1°C, food and water available ad libitum. Mice utilized in this study do harbor Segmented Filamentous Bacteria (SFB) (<xref ref-type="bibr" rid="bib24">Estes et al., 2019</xref>). Males and females embryos were analyzed across experiments, sex was determined as previously described (<xref ref-type="bibr" rid="bib51">McFarlane et al., 2013</xref>).</p></sec><sec id="s4-2"><title>Maternal immune activation</title><p>Assessment of females baseline immunoreactivity before pregnancy and maternal immune activation were performed as previously described (<xref ref-type="bibr" rid="bib25">Estes et al., 2020</xref>). Since female mice display a wide range of baseline immunoreactivity to poly(I:C) that dictates susceptibility and resilience of offspring from later pregnancies to MIA, virgin females in the lowest 25th percentile of immune responsiveness to poly(I:C) (IL-6 serum levels) were excluded from the study to reduce variability and ensure sufficient immune responsiveness of all included dams. E12.5 pregnancies were determined by visualizing vaginal seminal plug (noted as E0.5) and by body weight increase. IL-6 serum levels were measured at 2.5 hr post-injection. Temperature, weight, and sickness behavior were tracked as previously described (<xref ref-type="bibr" rid="bib25">Estes et al., 2020</xref>).</p></sec><sec id="s4-3"><title>Tissue dissections for RNA-seq and immunoblotting</title><p>Dissections of dorsal telencephalon/pallium or developing cortex were performed at E12.5 + 6 hr, E14.5, E17.5, and P0 from control (saline) and poly(I:C)-exposed MIA litters. At E12.5, due to the smaller embryo size, dissections were made via a diagonal cut that maximized collection of dorsal telencephalon, but could have also included some subpallium/ventral telencephalon and superficial tissues. For later time-points, dissection was of pallium/developing cortex. Potential variance in the individual dissections should be mitigated by the inclusion of multiple replicates across both groups. Dissections included both hemispheres from males and females, with exact numbers and sample details reported in <xref ref-type="supplementary-material" rid="supp20">Supplementary file 20</xref>. Independent litters were used for RNA-seq and immunoblotting.</p></sec><sec id="s4-4"><title>RNA-seq</title><p>Following tissue collection as described above, total RNA was isolated using Ambion RNAqueous total RNA Isolation Kit and assayed using Agilent RNA 6000 Nano Bioanalyzer kit/instrument. Stranded mRNA libraries were prepared using TruSeq Stranded mRNA kits. Eight to 12 samples per lane were pooled and sequenced on an Illumina HiSeq 4000 instrument using a single-end 50 bp protocol. Reads were aligned to mouse genome (mm9) using STAR (version 2.4.2a) (<xref ref-type="bibr" rid="bib19">Dobin et al., 2013</xref>) and gene counts produced using featureCounts (<xref ref-type="bibr" rid="bib42">Liao et al., 2014</xref>) and mm9 knownGenes. Quality assayed using FastQC (<xref ref-type="bibr" rid="bib3">Andrews, 2010</xref>) and RSeQC (<xref ref-type="bibr" rid="bib85">Wang et al., 2012</xref>), and samples that exhibited 3′ bias or poor exon distribution were discarded. Raw RNA-seq fastq files and a gene count matrix is available on GEO <ext-link ext-link-type="uri" xlink:href="https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE166376">GSE166376</ext-link>.</p></sec><sec id="s4-5"><title>Immunoblotting</title><p>Tissue dissections from male and female littermates from at least three litters per condition were performed as described above and collected in HBSS, immediately frozen on dry ice, and stored at −80°C until processed. Samples were lysed in 50 mM Tris HCl, pH 8, 140 nM NaCl, 1 mM EDTA, 10% glycerol, 0.5% NP40 and 0.25% Triton with protease inhibitor cocktail (Roche). After sonication, samples were spun down and the supernatant was used for a BCA Bradford assay using the Spectramax 190 plate reader to assess protein concentration using a standard curve. Twelve or 18 μg of protein were run on a 10 or 12% Tris acetate gel using the Mini-PROTEAN system (BioRad). Membranes were blocked in Odyssey blocking buffer (TBS; LiCor) and probed with the indicated primary antibodies overnight at 4°C. Resolved proteins were visualized in two channels using fluorescent secondary antibodies at 680 and 800 nm on an Odyssey Clx infrared imaging system (LiCor). Specific band intensities for all detected antibodies were quantified using the manufacturer’s software ImageSoft (LiCor) and normalized to GAPDH loading control. Antibodies used were anti-PAX6 (1:250; cat #PRB-278P-100; Covance, Princeton NJ.), anti-TBR1 (1:500; cat#ab31940; Abcam, Cambridge, MA), anti-CTIP2 (1:250; cat# ab18465; Abcam, Cambridge, MA), anti-CUX1 (1:200; cat#A2213; Abclonal, Woburn, MA) and anti-SATB2 (1:500; cat#ab51502; Abcam, Cambridge, MA).</p></sec><sec id="s4-6"><title>Gross anatomy, immunohistochemistry, and image quantification analyses</title><p>Histology was performed in triplicate sections in brains using male and female pups from at least two litters. Sections were selected to focus on the developing somatosensory cortex. Details on sample individual sample usage for each IHC experiment is detailed in <xref ref-type="supplementary-material" rid="supp21">Supplementary files 21</xref> and <xref ref-type="supplementary-material" rid="supp22">22</xref>. Embryos used in IHC experiments came from independent litters versus RNA-seq and immunoblotting. Morphological parameters were measured using FIJI ImageJ (NIH), by an experimenter blinded to group and sex. Brain regions were identified based on anatomical landmarks as previously described (<xref ref-type="bibr" rid="bib28">Gompers et al., 2017</xref>). No data points were excluded. For IHC staining, fetal brains were dissected and fixed with 4% paraformaldehyde/PBS solution overnight at 4°C. Tissues were then equilibrated in approximately 15 mL of 30% sucrose/PBS solution until they sank to the bottom of a conical tube. Equilibrated brains were embedded in Optimum Cutting Temperature (OCT) compound (Tissue-Tek, Torrance, CA) and frozen on dry-ice. OCT-embedded brain blocks were cryo-sectioned on a coronal plane (30 µm). Immunostaining was performed in free-floating sections with agitation. Sections were washed five times in PBST (PBS with 0.05% Triton X-100, 5 min each) and antigen retrieval was performed using 1x Citrate buffer pH6.0 antigen retriever solution (cat# C9999; Millipore-Sigma, Burlington, MA), at 60°C for 1 hr. Sections were washed five times in PBST (5 min each), permeabilized in PBS containing 0.5% Triton for 20 min and blocked for 1 hr at room temperature in 5% milk/PBST. Primary antibodies were incubated overnight at 4°C with orbital agitation (40–50 rpm). All antibodies used for this study were validated and their specificity validated. The following primary antibodies were used: anti-PAX6 (1:250; cat #PRB-278P-100; Covance, Princeton NJ.), anti-KI67 (1:500; cat#12202; Cell Signaling, Danvers, MA), anti-SOX9 (1:500 cat# AF3075, R and D systems, Minneapolis, MN), anti-PH3 (1:500 cat# 9701, Cell Signaling, Danvers, MA), anti-TBR2 (1:500; cat#14-4875-82; Thermo Fisher Scientific, Waltham, MA.), anti-TBR1 (1:500; cat#ab31940; Abcam, Cambridge, MA), anti-CTIP2 (1:250; cat# ab18465; Abcam, Cambridge, MA), anti-CUX1 (1:200; cat#A2213; Abclonal, Woburn, MA.), anti-SATB2 (1:500; cat#ab51502; Abcam, Cambridge, MA.), anti-DLX2 (1:200; generous gift from John Rubenstein, UCSF), anti-GFAP (1:250, cat#Z0334; Agilent Dako, Santa Clara, CA), anti-OLIG2 (1:500, cat#AB9610, MilliporeSigma, Burlington, MA). After primary antibody incubation, free-floating sections were washed five times in PBST (5 min each). Species-specific fluorophores-conjugated IgG (1:1000; Invitrogen-Thermo Fisher Scientific, Waltham, MA) were used as secondary antibodies (45 min, RT). 40,6-Diamidino-2-phenylindole (DAPI) (1:10000; Millipore-Sigma, Burlington, MA) was used for nuclear staining (20 min, RT). For cell counting, boxes were drawn on specific areas of the neocortex for progenitor and mitosis markers; SOX2, SOX9, and PH3. This box was 300 × 300 pixels on an image of the neocortex obtained using a ×20 objective. To be consistent, the ventral aspect of the box was drawn along the VZ and one corner touched the medial aspect of the VZ. The same strategy was used to draw boxes from Ki67 images except the box was increased in size to 500 × 500 pixels to account for the additional Ki67+ cells seen in neocortical areas above the SVZ. SOX9 boxes were drawn along the dorsal border edge (600 pixels in width) of the neocortex and extended ventrally two thirds the length of the neocortex in order to avoid the VZ/SVZ region. For laminar zone counts, a box spanning 300 pixels was drawn along the VZ and the dorsal/ventral proportions were equally divided into three zones. These box sizes allowed for a sufficient number of cells to be counted from each tissue. Laminar thickness analyses were performed by measuring the spread of all positive cells within the cortical plate per individual marker. For all quantification analyses, one to three sections per brain utilizing the developing hippocampus as a landmark for positioning were used.</p></sec><sec id="s4-7"><title>Bioinformatics analysis</title><p>Bioinformatic analysis was performed using R programming language version 3.5.1 (<xref ref-type="bibr" rid="bib63">R Development Core Team, 2015</xref>) run in RStudio integrated development environment version 1.2.1269 (<xref ref-type="bibr" rid="bib80">Team R, 2018</xref>). Plots were generated using ggplot2 R package version 3.1.0 (<xref ref-type="bibr" rid="bib86">Wickham, 2009</xref>). Heatmaps were generated using pheatmap R package 1.0.10 (<xref ref-type="bibr" rid="bib38">Kolde, 2018</xref>). The analysis scripts are available at: <ext-link ext-link-type="uri" xlink:href="https://github.com/NordNeurogenomicsLab/Publications/tree/master/Canales_eLife_2021">https://github.com/NordNeurogenomicsLab/Publications/tree/master/Canales_eLife_2021</ext-link>; <xref ref-type="bibr" rid="bib11">Canales, 2021</xref>; copy archived at <ext-link ext-link-type="uri" xlink:href="https://archive.softwareheritage.org/swh:1:dir:5f2ff57077946f69fd4fb99a2df240c5b5411a1e;origin=https://github.com/NordNeurogenomicsLab/Publications/;visit=swh:1:snp:afd3044f9dbd07b9cb7620b8ea08b8928c8fde6f;anchor=swh:1:rev:28836c8758908e130f1f6d8bfb3b6a112c6cbd1b/">swh:1:rev:28836c8758908e130f1f6d8bfb3b6a112c6cbd1b</ext-link>.</p></sec><sec id="s4-8"><title>Differential expression analysis</title><p>Raw count data for all samples were used as input along with sample information for differential expression analysis using edgeR (<xref ref-type="bibr" rid="bib65">Robinson et al., 2010</xref>). Genes with minimum log<sub>2</sub><underline>r</underline>eads <underline>p</underline>er <underline>k</underline>ilobase per <underline>m</underline>illion (RPKM) expression of −2 in at least two samples were included for analysis, resulting in a final set of 17,195 genes for differential testing. For time-point differential expression analysis genes with count per million (CPM) >0.1 in at least two samples were considered, resulting in E12.5, 16,396; E14.5, 16,658; E17.5, 16,229; P0, 16,613 genes in differential expression analyses. Principal component analysis indicated that the strongest driver of variance across samples was developmental age. Tagwise dispersion estimates were generated, and differential expression analysis was performed with edgeR using a generalized linear model, separately for each developmental time-point, including sex, sequencing run factor, and treatment as the variable for testing. Normalized expression levels were generated using the edgeR rpkm function. Normalized log<sub>2</sub> RPKM values were used for plotting of summary heatmaps and of expression data for individual genes.</p></sec><sec id="s4-9"><title>SFARI gene set enrichment analysis</title><p>The autism risk gene-set was downloaded from <ext-link ext-link-type="uri" xlink:href="https://gene.sfari.org/">https://gene.sfari.org/</ext-link> on 02-26-2019 and is included in this manuscript as <xref ref-type="supplementary-material" rid="supp17">Supplementary file 17</xref>. High confidence risk genes, annotated as ‘gene-score’ 1 and 2 were selected, and their orthologs were found using getLDS function from the biomaRt R package (<xref ref-type="bibr" rid="bib21">Durinck et al., 2009</xref>; <xref ref-type="bibr" rid="bib20">Durinck et al., 2005</xref>). Overlap of up- and downregulated DE genes with SFARI ortholog genes was calculated and plotted using a custom R script. Statistical significance of overlap was tested with hypergeometric test using the following R script: sum(dhyper(t:b, a, n - a, b)):</p><list list-type="simple"><list-item><p><italic>t</italic> = number of overlapping DE genes with SFARI gene orthologs with gene_score 1 or 2</p></list-item><list-item><p><italic>b</italic> = number of SFARI gene orthologs with gene_score 1 or 2</p></list-item><list-item><p><italic>a</italic> = number of DE genes at a developmental stage</p></list-item><list-item><p><italic>n</italic> = number of genes at a developmental stage <italic>sum</italic> returns the sum of all the values present in its arguments. <italic>dhyper</italic> calculates density for the hypergeometric distribution.</p></list-item></list></sec><sec id="s4-10"><title>WGCNA</title><p>We used the WGCNA R package, version 1.66 (<xref ref-type="bibr" rid="bib40">Langfelder and Horvath, 2012</xref>; <xref ref-type="bibr" rid="bib39">Langfelder and Horvath, 2008</xref>) to construct signed co-expression networks using the entire dataset containing 24,015 genes. After the network construction, the gene set was filtered for minimal gene expression at an RPKM value of 0.25 or higher in at least two sample, resulting in a dataset consisting 17,195 genes. A correlation matrix using the biweight midcorrelation between all genes was computed for all relevant samples. The soft thresholding power was estimated and used to derive an adjacency matrix exhibiting approximate scale-free topology (<italic>R</italic><sup>2</sup> >0.85). The adjacency matrix was transformed to a topological overlap matrix (TOM). The matrix 1-TOM was used as the input to calculate co-expression modules using hierarchical clustering. Modules were branches of the hierarchical cluster tree base, with minimum module size set to 10 genes. Pearson’s correlation coefficients were used to calculate correlation between sample traits (e.g., sex, treatment) and modules. The expression profile of a given module was summarized by the module eigengene (ME). Modules with highly correlated MEs (correlation >0.80) were merged together. The module connectivity (kME) of each gene was calculated by correlating the gene expression profile with module eigengenes. The module connectivity (kME) of each gene was calculated by correlating the gene expression profile with module eigengenes. Genes with no network correlation were placed into the module Grey. Following manual data inspection further highly correlated modules were merged: BrRePi: Brown, Red, and Pink, YeMaBl: Yellow, Magenta, and Black.</p></sec><sec id="s4-11"><title>Gene Ontology enrichment analysis</title><p>Mouse Gene Ontology (GO) data was downloaded from Bioconductor (org.Mm.eg.db). We used the TopGO R package version 2.34.0 (<xref ref-type="bibr" rid="bib2">Alexa and Rahnenfuhrer, 2019</xref>) to test for enrichment of GO terms. For the analysis presented here, we restricted our testing to GO Biological Process annotations and required a minimal node size (number of genes annotated to GO terms) of 20. We used the internal ‘weight01’ testing framework and the Fisher test, a strategy recommended for gene set analysis that generally accounts for multiple testing comparisons. For GO BP analysis, we reported terms with p-value<0.05. For all enrichment analysis, the test set of DE genes was compared against the background set of genes expressed in our study based on minimum read-count cutoffs described above. Heatmaps showing positive log<sub>2</sub> (expected/observed) values were plotted for GO terms of interest.</p></sec><sec id="s4-12"><title>Protein–protein interaction</title><p>Protein–protein interaction enrichment and network generation for E12.5 DE (FDR < 0.05) gene sets was performed using STRING (<xref ref-type="bibr" rid="bib77">Szklarczyk et al., 2019</xref>), version 11, considering only experimentally and text mining interactions, with at least medium interaction confidence score. Disconnected nodes were removed from the network.</p></sec><sec id="s4-13"><title>Data analysis</title><p>For all experiments, data was collected from at least two (usually three, and sometimes more, as indicated in each analysis) experiments and is presented as the mean ± SEM. Protein validation data were analyzed by Student’s <italic>t</italic>-test or one-way ANOVA, followed where appropriate by Tukey’s honestly significant difference <italic>post hoc</italic> test (Graphpad Prism v.7). Significance was defined as: *p<0.05, **p<0.01, ***p<0.001, ****p<0.0001. Samples for RNA-seq were randomly collected across litters and processed blind to experimental condition. For differential gene expression analysis using edgeR, differences were considered statistically significant at high stringency for FDR < 0.05, or with reduced stringency for p values < 0.05.</p></sec></sec></body><back><ack id="ack"><title>Acknowledgements</title><p>This work was supported by NIH NIMH (project # 5R21MH116681-02), NIH NIGMS (R35 GM119831), NARSAD Young Investigator Grant from Brain and Behavior Research Foundation and The UCD Clinical Translational Science Center. CPC is recipient of NIH NIMH Institutional National Research Service Award Autism Research Training Program fellowship (project # 2T32MH073124-16).</p></ack><sec id="s5" sec-type="additional-information"><title>Additional information</title><fn-group content-type="competing-interest"><title>Competing interests</title><fn fn-type="COI-statement" id="conf1"><p>No competing interests declared</p></fn></fn-group><fn-group content-type="author-contribution"><title>Author contributions</title><fn fn-type="con" id="con1"><p>Conceptualization, Formal analysis, Validation, Visualization, Methodology, Writing - original draft, Project administration, Writing - review and editing</p></fn><fn fn-type="con" id="con2"><p>Conceptualization, Visualization, Methodology, Writing - original draft, Writing - review and editing</p></fn><fn fn-type="con" id="con3"><p>Conceptualization, Data curation, Software, Formal analysis, Visualization, Writing - original draft, Writing - review and editing</p></fn><fn fn-type="con" id="con4"><p>Formal analysis, Methodology</p></fn><fn fn-type="con" id="con5"><p>Methodology</p></fn><fn fn-type="con" id="con6"><p>Methodology</p></fn><fn fn-type="con" id="con7"><p>Methodology</p></fn><fn fn-type="con" id="con8"><p>Formal analysis, Methodology</p></fn><fn fn-type="con" id="con9"><p>Methodology</p></fn><fn fn-type="con" id="con10"><p>Methodology</p></fn><fn fn-type="con" id="con11"><p>Methodology</p></fn><fn fn-type="con" id="con12"><p>Methodology</p></fn><fn fn-type="con" id="con13"><p>Data curation, Methodology</p></fn><fn fn-type="con" id="con14"><p>Methodology</p></fn><fn fn-type="con" id="con15"><p>Methodology</p></fn><fn fn-type="con" id="con16"><p>Methodology</p></fn><fn fn-type="con" id="con17"><p>Methodology</p></fn><fn fn-type="con" id="con18"><p>Conceptualization, Resources, Supervision, Investigation, Methodology, Writing - review and editing</p></fn><fn fn-type="con" id="con19"><p>Formal analysis, Supervision, Investigation, Writing - review and editing</p></fn><fn fn-type="con" id="con20"><p>Conceptualization, Resources, Supervision, Funding acquisition, Investigation, Visualization, Project administration, Writing - review and editing</p></fn><fn fn-type="con" id="con21"><p>Conceptualization, Resources, Data curation, Software, Formal analysis, Supervision, Funding acquisition, Investigation, Visualization, Methodology, Writing - original draft, Project administration, Writing - review and editing</p></fn></fn-group><fn-group content-type="ethics-information"><title>Ethics</title><fn fn-type="other" id="fn1"><p>Animal experimentation: This study was conducted in compliance with NIH guidelines and approved protocols from the University of California, Davis Animal Care and Use Committee (IACUC) protocol #20229.</p></fn></fn-group></sec><sec id="s6" sec-type="supplementary-material"><title>Additional files</title><supplementary-material id="supp1"><label>Supplementary file 1.</label><caption><title>Differential gene expression data of control vs MIA samples, and WGCNA gene module assignment at E12.5.</title><p>Columns indicate gene symbols (gene_name), log<sub>2</sub> fold change (LogFC), log<sub>2</sub> counts per million (LogCPM), likelihood ratio statistics (LR), p value (PValue), false discovery rate adjusted p value (FDR), and WGCNA module assignment (moduleColors and moduleColors_merged). A total of 16,396 genes were tested.</p></caption><media mime-subtype="zip" mimetype="application" xlink:href="elife-60100-supp1-v2.zip"/></supplementary-material><supplementary-material id="supp2"><label>Supplementary file 2.</label><caption><title>Differential gene expression data of control vs MIA samples, and WGCNA gene module assignment at E14.5.</title><p>Columns indicate gene symbols (gene_name), log<sub>2</sub> fold change (LogFC), log<sub>2</sub> counts per million (LogCPM), likelihood ratio statistics (LR), p value (PValue), false discovery rate adjusted p value (FDR), and WGCNA module assignment (moduleColors and moduleColors_merged). A total of 16,658 genes were tested.</p></caption><media mime-subtype="zip" mimetype="application" xlink:href="elife-60100-supp2-v2.zip"/></supplementary-material><supplementary-material id="supp3"><label>Supplementary file 3.</label><caption><title>Differential gene expression data of control vs MIA samples, and WGCNA gene module assignment at E17.5.</title><p>Columns indicate gene symbols (gene_name), log<sub>2</sub> fold change (LogFC), log<sub>2</sub> counts per million (LogCPM), likelihood ratio statistics (LR), p value (PValue), false discovery rate adjusted p value (FDR), and WGCNA module assignment (moduleColors and moduleColors_merged). A total of 16,229 genes were tested.</p></caption><media mime-subtype="zip" mimetype="application" xlink:href="elife-60100-supp3-v2.zip"/></supplementary-material><supplementary-material id="supp4"><label>Supplementary file 4.</label><caption><title>Differential gene expression data of control vs MIA samples, and WGCNA gene module assignment at E12.5.</title><p>Columns indicate gene symbols (gene_name), log<sub>2</sub> fold change (LogFC), log<sub>2</sub> counts per million (LogCPM), likelihood ratio statistics (LR), p value (PValue), false discovery rate adjusted p value (FDR), and WGCNA module assignment (moduleColors and moduleColors_merged). A total of 16,613 genes were tested.</p></caption><media mime-subtype="zip" mimetype="application" xlink:href="elife-60100-supp4-v2.zip"/></supplementary-material><supplementary-material id="supp5"><label>Supplementary file 5.</label><caption><title>RNA-seq sample metadata including Sample ID, age, treatment (Condition), classification as per their IL-6 prenatal levels (Response), sequencer lane batch parameters (HiSeq Lane, Library barcodes), animal and litter parameters (Animal ID, Litter ID, Liter Size), animal sex inferred from <italic>Xist</italic> expression (Sex by <italic>Xist</italic> counts), fastq file name, and fastq file md5sum hash value (md5sum).</title><p>Out of 74 samples, 73 were tracked for the Litter ID (explained in the Comments).</p></caption><media mime-subtype="zip" mimetype="application" xlink:href="elife-60100-supp5-v2.zip"/></supplementary-material><supplementary-material id="supp6"><label>Supplementary file 6.</label><caption><title>Differential gene expression data of control vs MIA male samples at E12.5.</title><p>Columns indicate gene symbols (gene_name), log<sub>2</sub> fold change (LogFC), log<sub>2</sub> counts per million (LogCPM), likelihood ratio statistics (LR), p value (PValue), and false discovery rate adjusted p value (FDR). A total of 16,171 genes were tested.</p></caption><media mime-subtype="zip" mimetype="application" xlink:href="elife-60100-supp6-v2.zip"/></supplementary-material><supplementary-material id="supp7"><label>Supplementary file 7.</label><caption><title>Differential gene expression data of control vs MIA male samples at E14.5.</title><p>Columns indicate gene symbols (gene_name), log<sub>2</sub> fold change (LogFC), log<sub>2</sub> counts per million (LogCPM), likelihood ratio statistics (LR), p value (PValue), and false discovery rate adjusted p value (FDR). A total of 16,568 genes were tested.</p></caption><media mime-subtype="zip" mimetype="application" xlink:href="elife-60100-supp7-v2.zip"/></supplementary-material><supplementary-material id="supp8"><label>Supplementary file 8.</label><caption><title>Differential gene expression data of control vs MIA male samples at E17.5.</title><p>Columns indicate gene symbols (gene_name), log<sub>2</sub> fold change (LogFC), log<sub>2</sub> counts per million (LogCPM), likelihood ratio statistics (LR), p value (PValue), and false discovery rate adjusted p value (FDR). A total of 15,883 genes were tested.</p></caption><media mime-subtype="zip" mimetype="application" xlink:href="elife-60100-supp8-v2.zip"/></supplementary-material><supplementary-material id="supp9"><label>Supplementary file 9.</label><caption><title>Differential gene expression data of control vs MIA male samples at P0.</title><p>Columns indicate gene symbols (gene_name), log<sub>2</sub> fold change (LogFC), log<sub>2</sub> counts per million (LogCPM), likelihood ratio statistics (LR), p value (PValue), and false discovery rate adjusted p value (FDR). A total of 16,394 genes were tested.</p></caption><media mime-subtype="zip" mimetype="application" xlink:href="elife-60100-supp9-v2.zip"/></supplementary-material><supplementary-material id="supp10"><label>Supplementary file 10.</label><caption><title>Differential gene expression data of control vs MIA female samples at E12.5.</title><p>Columns indicate gene symbols (gene_name), log<sub>2</sub> fold change (LogFC), log<sub>2</sub> counts per million (LogCPM), likelihood ratio statistics (LR), p value (PValue), and false discovery rate adjusted p value (FDR). A total of 16,161 genes were tested.</p></caption><media mime-subtype="zip" mimetype="application" xlink:href="elife-60100-supp10-v2.zip"/></supplementary-material><supplementary-material id="supp11"><label>Supplementary file 11.</label><caption><title>Differential gene expression data of control vs MIA female samples at E14.5.</title><p>Columns indicate gene symbols (gene_name), log<sub>2</sub> fold change (LogFC), log<sub>2</sub> counts per million (LogCPM), likelihood ratio statistics (LR), p value (PValue), and false discovery rate adjusted p value (FDR). A total of 16,350 genes were tested.</p></caption><media mime-subtype="zip" mimetype="application" xlink:href="elife-60100-supp11-v2.zip"/></supplementary-material><supplementary-material id="supp12"><label>Supplementary file 12.</label><caption><title>Differential gene expression data of control vs MIA female samples at E17.5.</title><p>Columns indicate gene symbols (gene_name), log<sub>2</sub> fold change (LogFC), log<sub>2</sub> counts per million (LogCPM), likelihood ratio statistics (LR), p value (PValue), and false discovery rate adjusted p value (FDR). A total of 16,093 genes were tested.</p></caption><media mime-subtype="zip" mimetype="application" xlink:href="elife-60100-supp12-v2.zip"/></supplementary-material><supplementary-material id="supp13"><label>Supplementary file 13.</label><caption><title>Differential gene expression data of control vs MIA female samples at P0.</title><p>Columns indicate gene symbols (gene_name), log<sub>2</sub> fold change (LogFC), log<sub>2</sub> counts per million (LogCPM), likelihood ratio statistics (LR), p value (PValue), and false discovery rate adjusted p value (FDR). A total of 16,487 genes were tested.</p></caption><media mime-subtype="zip" mimetype="application" xlink:href="elife-60100-supp13-v2.zip"/></supplementary-material><supplementary-material id="supp14"><label>Supplementary file 14.</label><caption><title>Differential gene expression data for all four developmental ages (Age) analyzed without sex covariate.</title><p>Columns indicate gene symbols (gene_name), log<sub>2</sub> fold change (LogFC), log<sub>2</sub> counts per million (LogCPM), likelihood ratio statistics (LR), p value (PValue), and false discovery rate adjusted p value (FDR).</p></caption><media mime-subtype="zip" mimetype="application" xlink:href="elife-60100-supp14-v2.zip"/></supplementary-material><supplementary-material id="supp15"><label>Supplementary file 15.</label><caption><title>Differential gene expression data contrasting male vs female samples at each developmental timepoint (Age).</title><p>Columns indicate gene symbols (gene_name), log<sub>2</sub> fold change (LogFC), log<sub>2</sub> counts per million (LogCPM), likelihood ratio statistics (LR), p value (PValue), and false discovery rate adjusted p value (FDR).</p></caption><media mime-subtype="zip" mimetype="application" xlink:href="elife-60100-supp15-v2.zip"/></supplementary-material><supplementary-material id="supp16"><label>Supplementary file 16.</label><caption><title>High confidence SFARI mouse orthologs genes used in the analysis.</title><p>Columns indicate gene symbols and names, chromosomal location, genetic category, gene and syndromic scores, number of reports and Mouse Genome Informatics (MGI) symbols for mouse orthologs. 89 mouse orthologs with high confidence SFARI gene scores of one or two were included in the analysis.</p></caption><media mime-subtype="zip" mimetype="application" xlink:href="elife-60100-supp16-v2.zip"/></supplementary-material><supplementary-material id="supp17"><label>Supplementary file 17.</label><caption><title>Complete SFARI gene database available at the time of conducting the analysis.</title></caption><media mime-subtype="zip" mimetype="application" xlink:href="elife-60100-supp17-v2.zip"/></supplementary-material><supplementary-material id="supp18"><label>Supplementary file 18.</label><caption><title>Gene ontology biological process enrichment analysis for genes passing FDR < 0.05, for each developmental timepoint and direction of differential expression, summarized in <xref ref-type="fig" rid="fig2">Figure 2d</xref> heatmap.</title><p>Columns signify gene ontology ID number (GO.ID), gene ontology term (Term), the number of annotated genes in the term (Annotated), the number of genes passing FDR < 0.05 for up- or downregulated genes (Significant), the expected number of genes under the null hypothesis, Fisher’s p value (classicFisher), developmental age (Age), direction of differential expression relative to saline (DE_direction), enrichment.</p></caption><media mime-subtype="zip" mimetype="application" xlink:href="elife-60100-supp18-v2.zip"/></supplementary-material><supplementary-material id="supp19"><label>Supplementary file 19.</label><caption><title>Gene ontology biological process enrichment analysis for genes passing FDR < 0.05, for each developmental timepoint, direction of differential expression, and WGCNA module membership.</title><p>Columns signify gene ontology ID number (GO.ID), gene ontology term (Term), the number of annotated genes in the term (Annotated), the number of genes passing FDR < 0.05 for up- or downregulated genes and WGCNA module (Significant), the expected number of genes under the null hypothesis, Fisher’s p value (classicFisher), WGCNA module membership (module_color), enrichment, direction of differential expression (Direction), and developmental age (Age).</p></caption><media mime-subtype="zip" mimetype="application" xlink:href="elife-60100-supp19-v2.zip"/></supplementary-material><supplementary-material id="supp20"><label>Supplementary file 20.</label><caption><title>Numbers of RNA-seq samples divided by sex and treatment condition.</title></caption><media mime-subtype="zip" mimetype="application" xlink:href="elife-60100-supp20-v2.zip"/></supplementary-material><supplementary-material id="supp21"><label>Supplementary file 21.</label><caption><title>Summary of sample numbers (N) used for neuroanatomical IHC validation at E17.5.</title><p>Numbers represent samples analyzed for at least one of the individual listed markers or co-stains. Details of specific use of individual samples can be found in <xref ref-type="supplementary-material" rid="supp22">supplementary file 22</xref>.</p></caption><media mime-subtype="xlsx" mimetype="application" xlink:href="elife-60100-supp21-v2.xlsx"/></supplementary-material><supplementary-material id="supp22"><label>Supplementary file 22.</label><caption><title>List of all MIA and control samples analyzed through neuroanatomical IHC validation at E17.5, and its specific use for individual experiments.</title></caption><media mime-subtype="xlsx" mimetype="application" xlink:href="elife-60100-supp22-v2.xlsx"/></supplementary-material><supplementary-material id="transrepform"><label>Transparent reporting form</label><media mime-subtype="docx" mimetype="application" xlink:href="elife-60100-transrepform-v2.docx"/></supplementary-material></sec><sec id="s7" sec-type="data-availability"><title>Data availability</title><p>Raw and gene-count RNA-seq sequencing data are deposited in GEO under accession number GSE166376. Data analysis code is available at: <ext-link ext-link-type="uri" xlink:href="https://github.com/NordNeurogenomicsLab/Publications/tree/master/Canales_eLife_2021">https://github.com/NordNeurogenomicsLab/Publications/tree/master/Canales_eLife_2021</ext-link>, copy archived at <ext-link ext-link-type="uri" xlink:href="https://archive.softwareheritage.org/swh:1:rev:28836c8758908e130f1f6d8bfb3b6a112c6cbd1b/">https://archive.softwareheritage.org/swh:1:rev:28836c8758908e130f1f6d8bfb3b6a112c6cbd1b/</ext-link> Detailed results of differential gene expression analysis are included in the manuscript and supporting files (Supplementary files 1-15), and can be visualized using our interactive online browser at: <ext-link ext-link-type="uri" xlink:href="https://nordlab.shinyapps.io/MIA_RPKM_plots/">https://nordlab.shinyapps.io/MIA_RPKM_plots/</ext-link>.</p><p>The following datasets were generated:</p><p><element-citation id="dataset1" publication-type="data" specific-use="isSupplementedBy"><person-group person-group-type="author"><name><surname>Cichewicz</surname><given-names>K</given-names></name><name><surname>Nord</surname><given-names>AS</given-names></name></person-group><year iso-8601-date="2021">2021</year><data-title>A temporal map of maternal immune activation-induced changes reveals a shift in neurodevelopmental timing and perturbed cortical development in mice</data-title><source>NCBI Gene Expression Omnibus</source><pub-id assigning-authority="NCBI" pub-id-type="accession" xlink:href="https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE166376">GSE166376</pub-id></element-citation></p><p><element-citation id="dataset2" publication-type="data" specific-use="isSupplementedBy"><person-group person-group-type="author"><name><surname>Cichewicz</surname><given-names>K</given-names></name><name><surname>Nord</surname><given-names>AS</given-names></name></person-group><year iso-8601-date="2021">2021</year><data-title>A temporal map of maternal immune activation-induced changes reveals a shift in neurodevelopmental timing and perturbed cortical development in mice</data-title><source>NordNeurogenomicsLab</source><pub-id assigning-authority="other" pub-id-type="accession" xlink:href="https://github.com/NordNeurogenomicsLab/Publications/tree/master/Canales_eLife_2021">eLife_2021</pub-id></element-citation></p></sec><ref-list><title>References</title><ref id="bib1"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Akdis</surname> <given-names>M</given-names></name><name><surname>Aab</surname> <given-names>A</given-names></name><name><surname>Altunbulakli</surname> <given-names>C</given-names></name><name><surname>Azkur</surname> <given-names>K</given-names></name><name><surname>Costa</surname> <given-names>RA</given-names></name><name><surname>Crameri</surname> <given-names>R</given-names></name><name><surname>Duan</surname> <given-names>S</given-names></name><name><surname>Eiwegger</surname> <given-names>T</given-names></name><name><surname>Eljaszewicz</surname> <given-names>A</given-names></name><name><surname>Ferstl</surname> <given-names>R</given-names></name><name><surname>Frei</surname> <given-names>R</given-names></name><name><surname>Garbani</surname> <given-names>M</given-names></name><name><surname>Globinska</surname> <given-names>A</given-names></name><name><surname>Hess</surname> <given-names>L</given-names></name><name><surname>Huitema</surname> <given-names>C</given-names></name><name><surname>Kubo</surname> <given-names>T</given-names></name><name><surname>Komlosi</surname> <given-names>Z</given-names></name><name><surname>Konieczna</surname> <given-names>P</given-names></name><name><surname>Kovacs</surname> <given-names>N</given-names></name><name><surname>Kucuksezer</surname> <given-names>UC</given-names></name><name><surname>Meyer</surname> <given-names>N</given-names></name><name><surname>Morita</surname> <given-names>H</given-names></name><name><surname>Olzhausen</surname> <given-names>J</given-names></name><name><surname>O'Mahony</surname> <given-names>L</given-names></name><name><surname>Pezer</surname> <given-names>M</given-names></name><name><surname>Prati</surname> <given-names>M</given-names></name><name><surname>Rebane</surname> <given-names>A</given-names></name><name><surname>Rhyner</surname> <given-names>C</given-names></name><name><surname>Rinaldi</surname> <given-names>A</given-names></name><name><surname>Sokolowska</surname> <given-names>M</given-names></name><name><surname>Stanic</surname> <given-names>B</given-names></name><name><surname>Sugita</surname> <given-names>K</given-names></name><name><surname>Treis</surname> <given-names>A</given-names></name><name><surname>van de Veen</surname> <given-names>W</given-names></name><name><surname>Wanke</surname> <given-names>K</given-names></name><name><surname>Wawrzyniak</surname> <given-names>M</given-names></name><name><surname>Wawrzyniak</surname> <given-names>P</given-names></name><name><surname>Wirz</surname> <given-names>OF</given-names></name><name><surname>Zakzuk</surname> <given-names>JS</given-names></name><name><surname>Akdis</surname> <given-names>CA</given-names></name></person-group><year iso-8601-date="2016">2016</year><article-title>Interleukins (from IL-1 to IL-38), interferons, transforming growth factor β, and TNF-α: receptors, functions, and roles in diseases</article-title><source>Journal of Allergy and Clinical Immunology</source><volume>138</volume><fpage>984</fpage><lpage>1010</lpage><pub-id pub-id-type="doi">10.1016/j.jaci.2016.06.033</pub-id><pub-id pub-id-type="pmid">27577879</pub-id></element-citation></ref><ref id="bib2"><element-citation publication-type="software"><person-group person-group-type="author"><name><surname>Alexa</surname> <given-names>A</given-names></name><name><surname>Rahnenfuhrer</surname> <given-names>J</given-names></name></person-group><year iso-8601-date="2019">2019</year><data-title>topGO: Enrichment Analysis for Gene Ontology</data-title><version designator="2.38.1">2.38.1</version><publisher-name>R package</publisher-name><ext-link ext-link-type="uri" xlink:href="https://bioconductor.org/packages/release/bioc/html/topGO.html">https://bioconductor.org/packages/release/bioc/html/topGO.html</ext-link></element-citation></ref><ref id="bib3"><element-citation publication-type="software"><person-group person-group-type="author"><name><surname>Andrews</surname> <given-names>S</given-names></name></person-group><year iso-8601-date="2010">2010</year><data-title>FastQC A Quality Control tool for High Throughput Sequence Data</data-title><source>Babraham Bioinformatics</source><version designator="3"> GPL v3</version><ext-link ext-link-type="uri" xlink:href="https://www.bioinformatics.babraham.ac.uk/projects/fastqc/">https://www.bioinformatics.babraham.ac.uk/projects/fastqc/</ext-link></element-citation></ref><ref id="bib4"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Argaw</surname> <given-names>AT</given-names></name><name><surname>Asp</surname> <given-names>L</given-names></name><name><surname>Zhang</surname> <given-names>J</given-names></name><name><surname>Navrazhina</surname> <given-names>K</given-names></name><name><surname>Pham</surname> <given-names>T</given-names></name><name><surname>Mariani</surname> <given-names>JN</given-names></name><name><surname>Mahase</surname> <given-names>S</given-names></name><name><surname>Dutta</surname> <given-names>DJ</given-names></name><name><surname>Seto</surname> <given-names>J</given-names></name><name><surname>Kramer</surname> <given-names>EG</given-names></name><name><surname>Ferrara</surname> <given-names>N</given-names></name><name><surname>Sofroniew</surname> <given-names>MV</given-names></name><name><surname>John</surname> <given-names>GR</given-names></name></person-group><year iso-8601-date="2012">2012</year><article-title>Astrocyte-derived VEGF-A drives blood-brain barrier disruption in CNS inflammatory disease</article-title><source>Journal of Clinical Investigation</source><volume>122</volume><fpage>2454</fpage><lpage>2468</lpage><pub-id pub-id-type="doi">10.1172/JCI60842</pub-id><pub-id pub-id-type="pmid">22653056</pub-id></element-citation></ref><ref id="bib5"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Arrode-Brusés</surname> <given-names>G</given-names></name><name><surname>Brusés</surname> <given-names>JL</given-names></name></person-group><year iso-8601-date="2012">2012</year><article-title>Maternal immune activation by poly I:c induces expression of cytokines IL-1β and IL-13, chemokine MCP-1 and colony stimulating factor VEGF in fetal mouse brain</article-title><source>Journal of Neuroinflammation</source><volume>9</volume><elocation-id>83</elocation-id><pub-id pub-id-type="doi">10.1186/1742-2094-9-83</pub-id><pub-id pub-id-type="pmid">22546005</pub-id></element-citation></ref><ref id="bib6"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Ballendine</surname> <given-names>SA</given-names></name><name><surname>Greba</surname> <given-names>Q</given-names></name><name><surname>Dawicki</surname> <given-names>W</given-names></name><name><surname>Zhang</surname> <given-names>X</given-names></name><name><surname>Gordon</surname> <given-names>JR</given-names></name><name><surname>Howland</surname> <given-names>JG</given-names></name></person-group><year iso-8601-date="2015">2015</year><article-title>Behavioral alterations in rat offspring following maternal immune activation and ELR-CXC chemokine receptor antagonism during pregnancy: implications for neurodevelopmental psychiatric disorders</article-title><source>Progress in Neuro-Psychopharmacology and Biological Psychiatry</source><volume>57</volume><fpage>155</fpage><lpage>165</lpage><pub-id pub-id-type="doi">10.1016/j.pnpbp.2014.11.002</pub-id></element-citation></ref><ref id="bib7"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Bansod</surname> <given-names>S</given-names></name><name><surname>Kageyama</surname> <given-names>R</given-names></name><name><surname>Ohtsuka</surname> <given-names>T</given-names></name></person-group><year iso-8601-date="2017">2017</year><article-title>Hes5 regulates the transition timing of neurogenesis and gliogenesis in mammalian neocortical development</article-title><source>Development</source><volume>144</volume><fpage>3156</fpage><lpage>3167</lpage><pub-id pub-id-type="doi">10.1242/dev.147256</pub-id></element-citation></ref><ref id="bib8"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Basu</surname> <given-names>SN</given-names></name><name><surname>Kollu</surname> <given-names>R</given-names></name><name><surname>Banerjee-Basu</surname> <given-names>S</given-names></name></person-group><year iso-8601-date="2009">2009</year><article-title>AutDB: a gene reference resource for autism research</article-title><source>Nucleic Acids Research</source><volume>37</volume><fpage>D832</fpage><lpage>D836</lpage><pub-id pub-id-type="doi">10.1093/nar/gkn835</pub-id></element-citation></ref><ref id="bib9"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Bayless</surname> <given-names>NL</given-names></name><name><surname>Greenberg</surname> <given-names>RS</given-names></name><name><surname>Swigut</surname> <given-names>T</given-names></name><name><surname>Wysocka</surname> <given-names>J</given-names></name><name><surname>Blish</surname> <given-names>CA</given-names></name></person-group><year iso-8601-date="2016">2016</year><article-title>Zika virus infection induces cranial neural crest cells to produce cytokines at levels detrimental for neurogenesis</article-title><source>Cell Host & Microbe</source><volume>20</volume><fpage>423</fpage><lpage>428</lpage><pub-id pub-id-type="doi">10.1016/j.chom.2016.09.006</pub-id><pub-id pub-id-type="pmid">27693308</pub-id></element-citation></ref><ref id="bib10"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Brown</surname> <given-names>AS</given-names></name></person-group><year iso-8601-date="2011">2011</year><article-title>The environment and susceptibility to schizophrenia</article-title><source>Progress in Neurobiology</source><volume>93</volume><fpage>23</fpage><lpage>58</lpage><pub-id pub-id-type="doi">10.1016/j.pneurobio.2010.09.003</pub-id><pub-id pub-id-type="pmid">20955757</pub-id></element-citation></ref><ref id="bib11"><element-citation publication-type="software"><person-group person-group-type="author"><name><surname>Canales</surname> <given-names>C</given-names></name></person-group><year iso-8601-date="2021">2021</year><data-title>Publications</data-title><source>Software Heritage</source><version designator="swh:1:rev:28836c8758908e130f1f6d8bfb3b6a112c6cbd1b">swh:1:rev:28836c8758908e130f1f6d8bfb3b6a112c6cbd1b</version><ext-link ext-link-type="uri" xlink:href="https://archive.softwareheritage.org/swh:1:dir:5f2ff57077946f69fd4fb99a2df240c5b5411a1e;origin=https://github.com/NordNeurogenomicsLab/Publications/;visit=swh:1:snp:afd3044f9dbd07b9cb7620b8ea08b8928c8fde6f;anchor=swh:1:rev:28836c8758908e130f1f6d8bfb3b6a112c6cbd1b/">https://archive.softwareheritage.org/swh:1:dir:5f2ff57077946f69fd4fb99a2df240c5b5411a1e;origin=https://github.com/NordNeurogenomicsLab/Publications/;visit=swh:1:snp:afd3044f9dbd07b9cb7620b8ea08b8928c8fde6f;anchor=swh:1:rev:28836c8758908e130f1f6d8bfb3b6a112c6cbd1b/</ext-link></element-citation></ref><ref id="bib12"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Canetta</surname> <given-names>SE</given-names></name><name><surname>Bao</surname> <given-names>Y</given-names></name><name><surname>Co</surname> <given-names>MD</given-names></name><name><surname>Ennis</surname> <given-names>FA</given-names></name><name><surname>Cruz</surname> <given-names>J</given-names></name><name><surname>Terajima</surname> <given-names>M</given-names></name><name><surname>Shen</surname> <given-names>L</given-names></name><name><surname>Kellendonk</surname> <given-names>C</given-names></name><name><surname>Schaefer</surname> <given-names>CA</given-names></name><name><surname>Brown</surname> <given-names>AS</given-names></name></person-group><year iso-8601-date="2014">2014</year><article-title>Serological documentation of maternal influenza exposure and bipolar disorder in adult offspring</article-title><source>American Journal of Psychiatry</source><volume>171</volume><fpage>557</fpage><lpage>563</lpage><pub-id pub-id-type="doi">10.1176/appi.ajp.2013.13070943</pub-id><pub-id pub-id-type="pmid">24480930</pub-id></element-citation></ref><ref id="bib13"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Carpentier</surname> <given-names>PA</given-names></name><name><surname>Haditsch</surname> <given-names>U</given-names></name><name><surname>Braun</surname> <given-names>AE</given-names></name><name><surname>Cantu</surname> <given-names>AV</given-names></name><name><surname>Moon</surname> <given-names>HM</given-names></name><name><surname>Price</surname> <given-names>RO</given-names></name><name><surname>Anderson</surname> <given-names>MP</given-names></name><name><surname>Saravanapandian</surname> <given-names>V</given-names></name><name><surname>Ismail</surname> <given-names>K</given-names></name><name><surname>Rivera</surname> <given-names>M</given-names></name><name><surname>Weimann</surname> <given-names>JM</given-names></name><name><surname>Palmer</surname> <given-names>TD</given-names></name></person-group><year iso-8601-date="2013">2013</year><article-title>Stereotypical alterations in cortical patterning are associated with maternal illness-induced placental dysfunction</article-title><source>The Journal of Neuroscience</source><volume>33</volume><fpage>16874</fpage><lpage>16888</lpage><pub-id pub-id-type="doi">10.1523/JNEUROSCI.4654-12.2013</pub-id><pub-id pub-id-type="pmid">24155294</pub-id></element-citation></ref><ref id="bib14"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Choi</surname> <given-names>GB</given-names></name><name><surname>Yim</surname> <given-names>YS</given-names></name><name><surname>Wong</surname> <given-names>H</given-names></name><name><surname>Kim</surname> <given-names>S</given-names></name><name><surname>Kim</surname> <given-names>H</given-names></name><name><surname>Kim</surname> <given-names>SV</given-names></name><name><surname>Hoeffer</surname> <given-names>CA</given-names></name><name><surname>Littman</surname> <given-names>DR</given-names></name><name><surname>Huh</surname> <given-names>JR</given-names></name></person-group><year iso-8601-date="2016">2016</year><article-title>The maternal interleukin-17a pathway in mice promotes autism-like phenotypes in offspring</article-title><source>Science</source><volume>351</volume><fpage>933</fpage><lpage>939</lpage><pub-id pub-id-type="doi">10.1126/science.aad0314</pub-id><pub-id pub-id-type="pmid">26822608</pub-id></element-citation></ref><ref id="bib15"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Cooper</surname> <given-names>JA</given-names></name></person-group><year iso-8601-date="2013">2013</year><article-title>Mechanisms of cell migration in the nervous system</article-title><source>Journal of Cell Biology</source><volume>202</volume><fpage>725</fpage><lpage>734</lpage><pub-id pub-id-type="doi">10.1083/jcb.201305021</pub-id></element-citation></ref><ref id="bib16"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Dahlgren</surname> <given-names>J</given-names></name><name><surname>Samuelsson</surname> <given-names>AM</given-names></name><name><surname>Jansson</surname> <given-names>T</given-names></name><name><surname>Holmäng</surname> <given-names>A</given-names></name></person-group><year iso-8601-date="2006">2006</year><article-title>Interleukin-6 in the maternal circulation reaches the rat fetus in mid-gestation</article-title><source>Pediatric Research</source><volume>60</volume><fpage>147</fpage><lpage>151</lpage><pub-id pub-id-type="doi">10.1203/01.pdr.0000230026.74139.18</pub-id><pub-id pub-id-type="pmid">16864694</pub-id></element-citation></ref><ref id="bib17"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Depino</surname> <given-names>AM</given-names></name></person-group><year iso-8601-date="2015">2015</year><article-title>Early prenatal exposure to LPS results in anxiety- and depression-related behaviors in adulthood</article-title><source>Neuroscience</source><volume>299</volume><fpage>56</fpage><lpage>65</lpage><pub-id pub-id-type="doi">10.1016/j.neuroscience.2015.04.065</pub-id><pub-id pub-id-type="pmid">25943476</pub-id></element-citation></ref><ref id="bib18"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Dimou</surname> <given-names>L</given-names></name><name><surname>Simon</surname> <given-names>C</given-names></name><name><surname>Kirchhoff</surname> <given-names>F</given-names></name><name><surname>Takebayashi</surname> <given-names>H</given-names></name><name><surname>Götz</surname> <given-names>M</given-names></name></person-group><year iso-8601-date="2008">2008</year><article-title>Progeny of Olig2-expressing progenitors in the gray and white matter of the adult mouse cerebral cortex</article-title><source>Journal of Neuroscience</source><volume>28</volume><fpage>10434</fpage><lpage>10442</lpage><pub-id pub-id-type="doi">10.1523/JNEUROSCI.2831-08.2008</pub-id><pub-id pub-id-type="pmid">18842903</pub-id></element-citation></ref><ref id="bib19"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Dobin</surname> <given-names>A</given-names></name><name><surname>Davis</surname> <given-names>CA</given-names></name><name><surname>Schlesinger</surname> <given-names>F</given-names></name><name><surname>Drenkow</surname> <given-names>J</given-names></name><name><surname>Zaleski</surname> <given-names>C</given-names></name><name><surname>Jha</surname> <given-names>S</given-names></name><name><surname>Batut</surname> <given-names>P</given-names></name><name><surname>Chaisson</surname> <given-names>M</given-names></name><name><surname>Gingeras</surname> <given-names>TR</given-names></name></person-group><year iso-8601-date="2013">2013</year><article-title>STAR: ultrafast universal RNA-seq aligner</article-title><source>Bioinformatics</source><volume>29</volume><fpage>15</fpage><lpage>21</lpage><pub-id pub-id-type="doi">10.1093/bioinformatics/bts635</pub-id><pub-id pub-id-type="pmid">23104886</pub-id></element-citation></ref><ref id="bib20"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Durinck</surname> <given-names>S</given-names></name><name><surname>Moreau</surname> <given-names>Y</given-names></name><name><surname>Kasprzyk</surname> <given-names>A</given-names></name><name><surname>Davis</surname> <given-names>S</given-names></name><name><surname>De Moor</surname> <given-names>B</given-names></name><name><surname>Brazma</surname> <given-names>A</given-names></name><name><surname>Huber</surname> <given-names>W</given-names></name></person-group><year iso-8601-date="2005">2005</year><article-title>BioMart and bioconductor: a powerful link between biological databases and microarray data analysis</article-title><source>Bioinformatics</source><volume>21</volume><fpage>3439</fpage><lpage>3440</lpage><pub-id pub-id-type="doi">10.1093/bioinformatics/bti525</pub-id></element-citation></ref><ref id="bib21"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Durinck</surname> <given-names>S</given-names></name><name><surname>Spellman</surname> <given-names>PT</given-names></name><name><surname>Birney</surname> <given-names>E</given-names></name><name><surname>Huber</surname> <given-names>W</given-names></name></person-group><year iso-8601-date="2009">2009</year><article-title>Mapping identifiers for the integration of genomic datasets with the R/Bioconductor package biomaRt</article-title><source>Nature Protocols</source><volume>4</volume><fpage>1184</fpage><lpage>1191</lpage><pub-id pub-id-type="doi">10.1038/nprot.2009.97</pub-id><pub-id pub-id-type="pmid">19617889</pub-id></element-citation></ref><ref id="bib22"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Elmer</surname> <given-names>BM</given-names></name><name><surname>Estes</surname> <given-names>ML</given-names></name><name><surname>Barrow</surname> <given-names>SL</given-names></name><name><surname>McAllister</surname> <given-names>AK</given-names></name></person-group><year iso-8601-date="2013">2013</year><article-title>MHCI requires MEF2 transcription factors to negatively regulate synapse density during development and in disease</article-title><source>Journal of Neuroscience</source><volume>33</volume><fpage>13791</fpage><lpage>13804</lpage><pub-id pub-id-type="doi">10.1523/JNEUROSCI.2366-13.2013</pub-id><pub-id pub-id-type="pmid">23966700</pub-id></element-citation></ref><ref id="bib23"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Eng</surname> <given-names>LF</given-names></name><name><surname>Ghirnikar</surname> <given-names>RS</given-names></name><name><surname>Lee</surname> <given-names>YL</given-names></name></person-group><year iso-8601-date="2000">2000</year><article-title>Glial fibrillary acidic protein: gfap-thirty-one years (1969-2000)</article-title><source>Neurochemical Research</source><volume>25</volume><fpage>1439</fpage><lpage>1451</lpage><pub-id pub-id-type="doi">10.1023/a:1007677003387</pub-id><pub-id pub-id-type="pmid">11059815</pub-id></element-citation></ref><ref id="bib24"><element-citation publication-type="preprint"><person-group person-group-type="author"><name><surname>Estes</surname> <given-names>ML</given-names></name><name><surname>Farrelly</surname> <given-names>K</given-names></name><name><surname>Cameron</surname> <given-names>S</given-names></name><name><surname>Aboubechara</surname> <given-names>JP</given-names></name><name><surname>Haapanen</surname> <given-names>L</given-names></name><name><surname>Schauer</surname> <given-names>JD</given-names></name><name><surname>Horta</surname> <given-names>A</given-names></name><name><surname>Prendergast</surname> <given-names>K</given-names></name><name><surname>MacMahon</surname> <given-names>JA</given-names></name><name><surname>Shaffer</surname> <given-names>CI</given-names></name><name><surname>Ct</surname> <given-names>L</given-names></name><name><surname>Kincheloe</surname> <given-names>GN</given-names></name><name><surname>Tan</surname> <given-names>DJ</given-names></name><name><surname>Van der</surname> <given-names>LD</given-names></name><name><surname>Bauman</surname> <given-names>MD</given-names></name><name><surname>Carter</surname> <given-names>CS</given-names></name><name><surname>de</surname> <given-names>WJV</given-names></name><name><surname>McAllister</surname> <given-names>AK</given-names></name></person-group><year iso-8601-date="2019">2019</year><article-title>Enhancing rigor and reproducibility in maternal immune activation models: practical considerations and predicting resilience and susceptibility using baseline immune responsiveness before pregnancy</article-title><source>bioRxiv</source><pub-id pub-id-type="doi">10.1101/699983</pub-id></element-citation></ref><ref id="bib25"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Estes</surname> <given-names>ML</given-names></name><name><surname>Prendergast</surname> <given-names>K</given-names></name><name><surname>MacMahon</surname> <given-names>JA</given-names></name><name><surname>Cameron</surname> <given-names>S</given-names></name><name><surname>Aboubechara</surname> <given-names>JP</given-names></name><name><surname>Farrelly</surname> <given-names>K</given-names></name><name><surname>Sell</surname> <given-names>GL</given-names></name><name><surname>Haapanen</surname> <given-names>L</given-names></name><name><surname>Schauer</surname> <given-names>JD</given-names></name><name><surname>Horta</surname> <given-names>A</given-names></name><name><surname>Shaffer</surname> <given-names>IC</given-names></name><name><surname>Le</surname> <given-names>CT</given-names></name><name><surname>Kincheloe</surname> <given-names>GN</given-names></name><name><surname>Tan</surname> <given-names>DJ</given-names></name><name><surname>van der List</surname> <given-names>D</given-names></name><name><surname>Bauman</surname> <given-names>MD</given-names></name><name><surname>Carter</surname> <given-names>CS</given-names></name><name><surname>Van de Water</surname> <given-names>J</given-names></name><name><surname>McAllister</surname> <given-names>AK</given-names></name></person-group><year iso-8601-date="2020">2020</year><article-title>Baseline immunoreactivity before pregnancy and poly(I:c) dose combine to dictate susceptibility and resilience of offspring to maternal immune activation</article-title><source>Brain, Behavior, and Immunity</source><volume>88</volume><fpage>619</fpage><lpage>630</lpage><pub-id pub-id-type="doi">10.1016/j.bbi.2020.04.061</pub-id><pub-id pub-id-type="pmid">32335198</pub-id></element-citation></ref><ref id="bib26"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Estes</surname> <given-names>ML</given-names></name><name><surname>McAllister</surname> <given-names>AK</given-names></name></person-group><year iso-8601-date="2015">2015</year><article-title>Immune mediators in the brain and peripheral tissues in autism spectrum disorder</article-title><source>Nature Reviews Neuroscience</source><volume>16</volume><fpage>469</fpage><lpage>486</lpage><pub-id pub-id-type="doi">10.1038/nrn3978</pub-id><pub-id pub-id-type="pmid">26189694</pub-id></element-citation></ref><ref id="bib27"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Estes</surname> <given-names>ML</given-names></name><name><surname>McAllister</surname> <given-names>AK</given-names></name></person-group><year iso-8601-date="2016">2016</year><article-title>Maternal immune activation: implications for neuropsychiatric disorders</article-title><source>Science</source><volume>353</volume><fpage>772</fpage><lpage>777</lpage><pub-id pub-id-type="doi">10.1126/science.aag3194</pub-id><pub-id pub-id-type="pmid">27540164</pub-id></element-citation></ref><ref id="bib28"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Gompers</surname> <given-names>AL</given-names></name><name><surname>Su-Feher</surname> <given-names>L</given-names></name><name><surname>Ellegood</surname> <given-names>J</given-names></name><name><surname>Copping</surname> <given-names>NA</given-names></name><name><surname>Riyadh</surname> <given-names>MA</given-names></name><name><surname>Stradleigh</surname> <given-names>TW</given-names></name><name><surname>Pride</surname> <given-names>MC</given-names></name><name><surname>Schaffler</surname> <given-names>MD</given-names></name><name><surname>Wade</surname> <given-names>AA</given-names></name><name><surname>Catta-Preta</surname> <given-names>R</given-names></name><name><surname>Zdilar</surname> <given-names>I</given-names></name><name><surname>Louis</surname> <given-names>S</given-names></name><name><surname>Kaushik</surname> <given-names>G</given-names></name><name><surname>Mannion</surname> <given-names>BJ</given-names></name><name><surname>Plajzer-Frick</surname> <given-names>I</given-names></name><name><surname>Afzal</surname> <given-names>V</given-names></name><name><surname>Visel</surname> <given-names>A</given-names></name><name><surname>Pennacchio</surname> <given-names>LA</given-names></name><name><surname>Dickel</surname> <given-names>DE</given-names></name><name><surname>Lerch</surname> <given-names>JP</given-names></name><name><surname>Crawley</surname> <given-names>JN</given-names></name><name><surname>Zarbalis</surname> <given-names>KS</given-names></name><name><surname>Silverman</surname> <given-names>JL</given-names></name><name><surname>Nord</surname> <given-names>AS</given-names></name></person-group><year iso-8601-date="2017">2017</year><article-title>Germline Chd8 haploinsufficiency alters brain development in mouse</article-title><source>Nature Neuroscience</source><volume>20</volume><fpage>1062</fpage><lpage>1073</lpage><pub-id pub-id-type="doi">10.1038/nn.4592</pub-id><pub-id pub-id-type="pmid">28671691</pub-id></element-citation></ref><ref id="bib29"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Haida</surname> <given-names>O</given-names></name><name><surname>Al Sagheer</surname> <given-names>T</given-names></name><name><surname>Balbous</surname> <given-names>A</given-names></name><name><surname>Francheteau</surname> <given-names>M</given-names></name><name><surname>Matas</surname> <given-names>E</given-names></name><name><surname>Soria</surname> <given-names>F</given-names></name><name><surname>Fernagut</surname> <given-names>PO</given-names></name><name><surname>Jaber</surname> <given-names>M</given-names></name></person-group><year iso-8601-date="2019">2019</year><article-title>Sex-dependent behavioral deficits and neuropathology in a maternal immune activation model of autism</article-title><source>Translational Psychiatry</source><volume>9</volume><elocation-id>e0457</elocation-id><pub-id pub-id-type="doi">10.1038/s41398-019-0457-y</pub-id></element-citation></ref><ref id="bib30"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Hendzel</surname> <given-names>MJ</given-names></name><name><surname>Wei</surname> <given-names>Y</given-names></name><name><surname>Mancini</surname> <given-names>MA</given-names></name><name><surname>Van Hooser</surname> <given-names>A</given-names></name><name><surname>Ranalli</surname> <given-names>T</given-names></name><name><surname>Brinkley</surname> <given-names>BR</given-names></name><name><surname>Bazett-Jones</surname> <given-names>DP</given-names></name><name><surname>Allis</surname> <given-names>CD</given-names></name></person-group><year iso-8601-date="1997">1997</year><article-title>Mitosis-specific phosphorylation of histone H3 initiates primarily within pericentromeric heterochromatin during G2 and spreads in an ordered fashion coincident with mitotic chromosome condensation</article-title><source>Chromosoma</source><volume>106</volume><fpage>348</fpage><lpage>360</lpage><pub-id pub-id-type="doi">10.1007/s004120050256</pub-id><pub-id pub-id-type="pmid">9362543</pub-id></element-citation></ref><ref id="bib31"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Hohman</surname> <given-names>TJ</given-names></name><name><surname>Bell</surname> <given-names>SP</given-names></name><name><surname>Jefferson</surname> <given-names>AL</given-names></name><collab>Alzheimer’s Disease Neuroimaging Initiative</collab></person-group><year iso-8601-date="2015">2015</year><article-title>The role of vascular endothelial growth factor in neurodegeneration and cognitive decline: exploring interactions with biomarkers of alzheimer disease</article-title><source>JAMA Neurology</source><volume>72</volume><fpage>520</fpage><lpage>529</lpage><pub-id pub-id-type="doi">10.1001/jamaneurol.2014.4761</pub-id><pub-id pub-id-type="pmid">25751166</pub-id></element-citation></ref><ref id="bib32"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Hsiao</surname> <given-names>EY</given-names></name><name><surname>Patterson</surname> <given-names>PH</given-names></name></person-group><year iso-8601-date="2011">2011</year><article-title>Activation of the maternal immune system induces endocrine changes in the placenta via IL-6</article-title><source>Brain, Behavior, and Immunity</source><volume>25</volume><fpage>604</fpage><lpage>615</lpage><pub-id pub-id-type="doi">10.1016/j.bbi.2010.12.017</pub-id><pub-id pub-id-type="pmid">21195166</pub-id></element-citation></ref><ref id="bib33"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Iga</surname> <given-names>J</given-names></name><name><surname>Ueno</surname> <given-names>S</given-names></name><name><surname>Yamauchi</surname> <given-names>K</given-names></name><name><surname>Numata</surname> <given-names>S</given-names></name><name><surname>Tayoshi-Shibuya</surname> <given-names>S</given-names></name><name><surname>Kinouchi</surname> <given-names>S</given-names></name><name><surname>Nakataki</surname> <given-names>M</given-names></name><name><surname>Song</surname> <given-names>H</given-names></name><name><surname>Hokoishi</surname> <given-names>K</given-names></name><name><surname>Tanabe</surname> <given-names>H</given-names></name><name><surname>Sano</surname> <given-names>A</given-names></name><name><surname>Ohmori</surname> <given-names>T</given-names></name></person-group><year iso-8601-date="2007">2007</year><article-title>Gene expression and association analysis of vascular endothelial growth factor in major depressive disorder</article-title><source>Progress in Neuro-Psychopharmacology & Biological Psychiatry</source><volume>31</volume><fpage>658</fpage><lpage>663</lpage><pub-id pub-id-type="doi">10.1016/j.pnpbp.2006.12.011</pub-id><pub-id pub-id-type="pmid">17257729</pub-id></element-citation></ref><ref id="bib34"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Isung</surname> <given-names>J</given-names></name><name><surname>Aeinehband</surname> <given-names>S</given-names></name><name><surname>Mobarrez</surname> <given-names>F</given-names></name><name><surname>Mårtensson</surname> <given-names>B</given-names></name><name><surname>Nordström</surname> <given-names>P</given-names></name><name><surname>Asberg</surname> <given-names>M</given-names></name><name><surname>Piehl</surname> <given-names>F</given-names></name><name><surname>Jokinen</surname> <given-names>J</given-names></name></person-group><year iso-8601-date="2012">2012</year><article-title>Low vascular endothelial growth factor and interleukin-8 in cerebrospinal fluid of suicide attempters</article-title><source>Translational Psychiatry</source><volume>2</volume><elocation-id>e196</elocation-id><pub-id pub-id-type="doi">10.1038/tp.2012.123</pub-id><pub-id pub-id-type="pmid">23168996</pub-id></element-citation></ref><ref id="bib35"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Javaherian</surname> <given-names>A</given-names></name><name><surname>Kriegstein</surname> <given-names>A</given-names></name></person-group><year iso-8601-date="2009">2009</year><article-title>A stem cell niche for intermediate progenitor cells of the embryonic cortex</article-title><source>Cerebral Cortex</source><volume>19</volume><fpage>i70</fpage><lpage>i77</lpage><pub-id pub-id-type="doi">10.1093/cercor/bhp029</pub-id></element-citation></ref><ref id="bib36"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Kahl</surname> <given-names>KG</given-names></name><name><surname>Bens</surname> <given-names>S</given-names></name><name><surname>Ziegler</surname> <given-names>K</given-names></name><name><surname>Rudolf</surname> <given-names>S</given-names></name><name><surname>Kordon</surname> <given-names>A</given-names></name><name><surname>Dibbelt</surname> <given-names>L</given-names></name><name><surname>Schweiger</surname> <given-names>U</given-names></name></person-group><year iso-8601-date="2009">2009</year><article-title>Angiogenic factors in patients with current major depressive disorder comorbid with borderline personality disorder</article-title><source>Psychoneuroendocrinology</source><volume>34</volume><fpage>353</fpage><lpage>357</lpage><pub-id pub-id-type="doi">10.1016/j.psyneuen.2008.09.016</pub-id><pub-id pub-id-type="pmid">19062198</pub-id></element-citation></ref><ref id="bib37"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Klum</surname> <given-names>S</given-names></name><name><surname>Zaouter</surname> <given-names>C</given-names></name><name><surname>Alekseenko</surname> <given-names>Z</given-names></name><name><surname>Björklund</surname> <given-names>AK</given-names></name><name><surname>Hagey</surname> <given-names>DW</given-names></name><name><surname>Ericson</surname> <given-names>J</given-names></name><name><surname>Muhr</surname> <given-names>J</given-names></name><name><surname>Bergsland</surname> <given-names>M</given-names></name></person-group><year iso-8601-date="2018">2018</year><article-title>Sequentially acting SOX proteins orchestrate astrocyte‐ and oligodendrocyte‐specific gene expression</article-title><source>EMBO Reports</source><volume>19</volume><elocation-id>e46635</elocation-id><pub-id pub-id-type="doi">10.15252/embr.201846635</pub-id></element-citation></ref><ref id="bib38"><element-citation publication-type="software"><person-group person-group-type="author"><name><surname>Kolde</surname> <given-names>R</given-names></name></person-group><year iso-8601-date="2018">2018</year><source>pheatmap</source><version designator="1.0.12">Pretty Heatmaps version 1.0.12</version></element-citation></ref><ref id="bib39"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Langfelder</surname> <given-names>P</given-names></name><name><surname>Horvath</surname> <given-names>S</given-names></name></person-group><year iso-8601-date="2008">2008</year><article-title>WGCNA: an R package for weighted correlation network analysis</article-title><source>BMC Bioinformatics</source><volume>9</volume><elocation-id>559</elocation-id><pub-id pub-id-type="doi">10.1186/1471-2105-9-559</pub-id><pub-id pub-id-type="pmid">19114008</pub-id></element-citation></ref><ref id="bib40"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Langfelder</surname> <given-names>P</given-names></name><name><surname>Horvath</surname> <given-names>S</given-names></name></person-group><year iso-8601-date="2012">2012</year><article-title>Fast R functions for robust correlations and hierarchical clustering</article-title><source>Journal of Statistical Software</source><volume>46</volume><fpage>1</fpage><lpage>17</lpage><pub-id pub-id-type="doi">10.18637/jss.v046.i11</pub-id><pub-id pub-id-type="pmid">23050260</pub-id></element-citation></ref><ref id="bib41"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Lee</surname> <given-names>BH</given-names></name><name><surname>Kim</surname> <given-names>YK</given-names></name></person-group><year iso-8601-date="2012">2012</year><article-title>Increased plasma VEGF levels in major depressive or manic episodes in patients with mood disorders</article-title><source>Journal of Affective Disorders</source><volume>136</volume><fpage>181</fpage><lpage>184</lpage><pub-id pub-id-type="doi">10.1016/j.jad.2011.07.021</pub-id><pub-id pub-id-type="pmid">21862441</pub-id></element-citation></ref><ref id="bib42"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Liao</surname> <given-names>Y</given-names></name><name><surname>Smyth</surname> <given-names>GK</given-names></name><name><surname>Shi</surname> <given-names>W</given-names></name></person-group><year iso-8601-date="2014">2014</year><article-title>featureCounts: an efficient general purpose program for assigning sequence reads to genomic features</article-title><source>Bioinformatics</source><volume>30</volume><fpage>923</fpage><lpage>930</lpage><pub-id pub-id-type="doi">10.1093/bioinformatics/btt656</pub-id><pub-id pub-id-type="pmid">24227677</pub-id></element-citation></ref><ref id="bib43"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Lin</surname> <given-names>H</given-names></name><name><surname>Muramatsu</surname> <given-names>R</given-names></name><name><surname>Maedera</surname> <given-names>N</given-names></name><name><surname>Tsunematsu</surname> <given-names>H</given-names></name><name><surname>Hamaguchi</surname> <given-names>M</given-names></name><name><surname>Koyama</surname> <given-names>Y</given-names></name><name><surname>Kuroda</surname> <given-names>M</given-names></name><name><surname>Ono</surname> <given-names>K</given-names></name><name><surname>Sawada</surname> <given-names>M</given-names></name><name><surname>Yamashita</surname> <given-names>T</given-names></name></person-group><year iso-8601-date="2018">2018</year><article-title>Extracellular lactate dehydrogenase A release from damaged neurons drives central nervous system angiogenesis</article-title><source>EBioMedicine</source><volume>27</volume><fpage>71</fpage><lpage>85</lpage><pub-id pub-id-type="doi">10.1016/j.ebiom.2017.10.033</pub-id><pub-id pub-id-type="pmid">29248508</pub-id></element-citation></ref><ref id="bib44"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Liu</surname> <given-names>KE</given-names></name><name><surname>Frazier</surname> <given-names>WA</given-names></name></person-group><year iso-8601-date="2015">2015</year><article-title>Phosphorylation of the BNIP3 C-Terminus inhibits mitochondrial damage and cell death without blocking autophagy</article-title><source>PLOS ONE</source><volume>10</volume><elocation-id>e0129667</elocation-id><pub-id pub-id-type="doi">10.1371/journal.pone.0129667</pub-id><pub-id pub-id-type="pmid">26102349</pub-id></element-citation></ref><ref id="bib45"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Lizano</surname> <given-names>P</given-names></name><name><surname>Lutz</surname> <given-names>O</given-names></name><name><surname>Ling</surname> <given-names>G</given-names></name><name><surname>Padmanabhan</surname> <given-names>J</given-names></name><name><surname>Tandon</surname> <given-names>N</given-names></name><name><surname>Sweeney</surname> <given-names>J</given-names></name><name><surname>Tamminga</surname> <given-names>C</given-names></name><name><surname>Pearlson</surname> <given-names>G</given-names></name><name><surname>Ruaño</surname> <given-names>G</given-names></name><name><surname>Kocherla</surname> <given-names>M</given-names></name><name><surname>Windemuth</surname> <given-names>A</given-names></name><name><surname>Clementz</surname> <given-names>B</given-names></name><name><surname>Gershon</surname> <given-names>E</given-names></name><name><surname>Keshavan</surname> <given-names>M</given-names></name></person-group><year iso-8601-date="2018">2018</year><article-title>VEGFA GENE variation influences hallucinations and frontotemporal morphology in psychotic disorders: a B-SNIP study</article-title><source>Translational Psychiatry</source><volume>8</volume><elocation-id>215</elocation-id><pub-id pub-id-type="doi">10.1038/s41398-018-0271-y</pub-id><pub-id pub-id-type="pmid">30310054</pub-id></element-citation></ref><ref id="bib46"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Lombardo</surname> <given-names>MV</given-names></name><name><surname>Moon</surname> <given-names>HM</given-names></name><name><surname>Su</surname> <given-names>J</given-names></name><name><surname>Palmer</surname> <given-names>TD</given-names></name><name><surname>Courchesne</surname> <given-names>E</given-names></name><name><surname>Pramparo</surname> <given-names>T</given-names></name></person-group><year iso-8601-date="2018">2018</year><article-title>Maternal immune activation dysregulation of the fetal brain transcriptome and relevance to the pathophysiology of autism spectrum disorder</article-title><source>Molecular Psychiatry</source><volume>23</volume><fpage>1001</fpage><lpage>1013</lpage><pub-id pub-id-type="doi">10.1038/mp.2017.15</pub-id><pub-id pub-id-type="pmid">28322282</pub-id></element-citation></ref><ref id="bib47"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Machado</surname> <given-names>CJ</given-names></name><name><surname>Whitaker</surname> <given-names>AM</given-names></name><name><surname>Smith</surname> <given-names>SE</given-names></name><name><surname>Patterson</surname> <given-names>PH</given-names></name><name><surname>Bauman</surname> <given-names>MD</given-names></name></person-group><year iso-8601-date="2015">2015</year><article-title>Maternal immune activation in nonhuman primates alters social attention in juvenile offspring</article-title><source>Biological Psychiatry</source><volume>77</volume><fpage>823</fpage><lpage>832</lpage><pub-id pub-id-type="doi">10.1016/j.biopsych.2014.07.035</pub-id><pub-id pub-id-type="pmid">25442006</pub-id></element-citation></ref><ref id="bib48"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Mackenzie</surname> <given-names>F</given-names></name><name><surname>Ruhrberg</surname> <given-names>C</given-names></name></person-group><year iso-8601-date="2012">2012</year><article-title>Diverse roles for VEGF-A in the nervous system</article-title><source>Development</source><volume>139</volume><fpage>1371</fpage><lpage>1380</lpage><pub-id pub-id-type="doi">10.1242/dev.072348</pub-id><pub-id pub-id-type="pmid">22434866</pub-id></element-citation></ref><ref id="bib49"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Malkova</surname> <given-names>NV</given-names></name><name><surname>Yu</surname> <given-names>CZ</given-names></name><name><surname>Hsiao</surname> <given-names>EY</given-names></name><name><surname>Moore</surname> <given-names>MJ</given-names></name><name><surname>Patterson</surname> <given-names>PH</given-names></name></person-group><year iso-8601-date="2012">2012</year><article-title>Maternal immune activation yields offspring displaying mouse versions of the three core symptoms of autism</article-title><source>Brain, Behavior, and Immunity</source><volume>26</volume><fpage>607</fpage><lpage>616</lpage><pub-id pub-id-type="doi">10.1016/j.bbi.2012.01.011</pub-id><pub-id pub-id-type="pmid">22310922</pub-id></element-citation></ref><ref id="bib50"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Martynoga</surname> <given-names>B</given-names></name><name><surname>Drechsel</surname> <given-names>D</given-names></name><name><surname>Guillemot</surname> <given-names>F</given-names></name></person-group><year iso-8601-date="2012">2012</year><article-title>Molecular control of neurogenesis: a view from the mammalian cerebral cortex</article-title><source>Cold Spring Harbor Perspectives in Biology</source><volume>4</volume><elocation-id>a008359</elocation-id><pub-id pub-id-type="doi">10.1101/cshperspect.a008359</pub-id><pub-id pub-id-type="pmid">23028117</pub-id></element-citation></ref><ref id="bib51"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>McFarlane</surname> <given-names>L</given-names></name><name><surname>Truong</surname> <given-names>V</given-names></name><name><surname>Palmer</surname> <given-names>JS</given-names></name><name><surname>Wilhelm</surname> <given-names>D</given-names></name></person-group><year iso-8601-date="2013">2013</year><article-title>Novel PCR assay for determining the genetic sex of mice</article-title><source>Sexual Development</source><volume>7</volume><fpage>207</fpage><lpage>211</lpage><pub-id pub-id-type="doi">10.1159/000348677</pub-id><pub-id pub-id-type="pmid">23571295</pub-id></element-citation></ref><ref id="bib52"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Meyer</surname> <given-names>U</given-names></name><name><surname>Feldon</surname> <given-names>J</given-names></name><name><surname>Schedlowski</surname> <given-names>M</given-names></name><name><surname>Yee</surname> <given-names>BK</given-names></name></person-group><year iso-8601-date="2006">2006a</year><article-title>Immunological stress at the maternal-foetal interface: a link between neurodevelopment and adult psychopathology</article-title><source>Brain, Behavior, and Immunity</source><volume>20</volume><fpage>378</fpage><lpage>388</lpage><pub-id pub-id-type="doi">10.1016/j.bbi.2005.11.003</pub-id><pub-id pub-id-type="pmid">16378711</pub-id></element-citation></ref><ref id="bib53"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Meyer</surname> <given-names>U</given-names></name><name><surname>Nyffeler</surname> <given-names>M</given-names></name><name><surname>Engler</surname> <given-names>A</given-names></name><name><surname>Urwyler</surname> <given-names>A</given-names></name><name><surname>Schedlowski</surname> <given-names>M</given-names></name><name><surname>Knuesel</surname> <given-names>I</given-names></name><name><surname>Yee</surname> <given-names>BK</given-names></name><name><surname>Feldon</surname> <given-names>J</given-names></name></person-group><year iso-8601-date="2006">2006b</year><article-title>The time of prenatal immune challenge determines the specificity of inflammation-mediated brain and behavioral pathology</article-title><source>Journal of Neuroscience</source><volume>26</volume><fpage>4752</fpage><lpage>4762</lpage><pub-id pub-id-type="doi">10.1523/JNEUROSCI.0099-06.2006</pub-id><pub-id pub-id-type="pmid">16672647</pub-id></element-citation></ref><ref id="bib54"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Meyer</surname> <given-names>U</given-names></name></person-group><year iso-8601-date="2014">2014</year><article-title>Prenatal poly(i:c) exposure and other developmental immune activation models in rodent systems</article-title><source>Biological Psychiatry</source><volume>75</volume><fpage>307</fpage><lpage>315</lpage><pub-id pub-id-type="doi">10.1016/j.biopsych.2013.07.011</pub-id><pub-id pub-id-type="pmid">23938317</pub-id></element-citation></ref><ref id="bib55"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Money</surname> <given-names>KM</given-names></name><name><surname>Barke</surname> <given-names>TL</given-names></name><name><surname>Serezani</surname> <given-names>A</given-names></name><name><surname>Gannon</surname> <given-names>M</given-names></name><name><surname>Garbett</surname> <given-names>KA</given-names></name><name><surname>Aronoff</surname> <given-names>DM</given-names></name><name><surname>Mirnics</surname> <given-names>K</given-names></name></person-group><year iso-8601-date="2018">2018</year><article-title>Gestational diabetes exacerbates maternal immune activation effects in the developing brain</article-title><source>Molecular Psychiatry</source><volume>23</volume><fpage>1920</fpage><lpage>1928</lpage><pub-id pub-id-type="doi">10.1038/mp.2017.191</pub-id><pub-id pub-id-type="pmid">28948973</pub-id></element-citation></ref><ref id="bib56"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Mootha</surname> <given-names>VK</given-names></name><name><surname>Lindgren</surname> <given-names>CM</given-names></name><name><surname>Eriksson</surname> <given-names>KF</given-names></name><name><surname>Subramanian</surname> <given-names>A</given-names></name><name><surname>Sihag</surname> <given-names>S</given-names></name><name><surname>Lehar</surname> <given-names>J</given-names></name><name><surname>Puigserver</surname> <given-names>P</given-names></name><name><surname>Carlsson</surname> <given-names>E</given-names></name><name><surname>Ridderstråle</surname> <given-names>M</given-names></name><name><surname>Laurila</surname> <given-names>E</given-names></name><name><surname>Houstis</surname> <given-names>N</given-names></name><name><surname>Daly</surname> <given-names>MJ</given-names></name><name><surname>Patterson</surname> <given-names>N</given-names></name><name><surname>Mesirov</surname> <given-names>JP</given-names></name><name><surname>Golub</surname> <given-names>TR</given-names></name><name><surname>Tamayo</surname> <given-names>P</given-names></name><name><surname>Spiegelman</surname> <given-names>B</given-names></name><name><surname>Lander</surname> <given-names>ES</given-names></name><name><surname>Hirschhorn</surname> <given-names>JN</given-names></name><name><surname>Altshuler</surname> <given-names>D</given-names></name><name><surname>Groop</surname> <given-names>LC</given-names></name></person-group><year iso-8601-date="2003">2003</year><article-title>PGC-1alpha-responsive genes involved in oxidative phosphorylation are coordinately downregulated in human diabetes</article-title><source>Nature Genetics</source><volume>34</volume><fpage>267</fpage><lpage>273</lpage><pub-id pub-id-type="doi">10.1038/ng1180</pub-id><pub-id pub-id-type="pmid">12808457</pub-id></element-citation></ref><ref id="bib57"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Naviaux</surname> <given-names>RK</given-names></name><name><surname>Zolkipli</surname> <given-names>Z</given-names></name><name><surname>Wang</surname> <given-names>L</given-names></name><name><surname>Nakayama</surname> <given-names>T</given-names></name><name><surname>Naviaux</surname> <given-names>JC</given-names></name><name><surname>Le</surname> <given-names>TP</given-names></name><name><surname>Schuchbauer</surname> <given-names>MA</given-names></name><name><surname>Rogac</surname> <given-names>M</given-names></name><name><surname>Tang</surname> <given-names>Q</given-names></name><name><surname>Dugan</surname> <given-names>LL</given-names></name><name><surname>Powell</surname> <given-names>SB</given-names></name></person-group><year iso-8601-date="2013">2013</year><article-title>Antipurinergic therapy corrects the autism-like features in the poly(IC) mouse model</article-title><source>PLOS ONE</source><volume>8</volume><elocation-id>e57380</elocation-id><pub-id pub-id-type="doi">10.1371/journal.pone.0057380</pub-id><pub-id pub-id-type="pmid">23516405</pub-id></element-citation></ref><ref id="bib58"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Newington</surname> <given-names>JT</given-names></name><name><surname>Rappon</surname> <given-names>T</given-names></name><name><surname>Albers</surname> <given-names>S</given-names></name><name><surname>Wong</surname> <given-names>DY</given-names></name><name><surname>Rylett</surname> <given-names>RJ</given-names></name><name><surname>Cumming</surname> <given-names>RC</given-names></name></person-group><year iso-8601-date="2012">2012</year><article-title>Overexpression of pyruvate dehydrogenase kinase 1 and lactate dehydrogenase A in nerve cells confers resistance to amyloid β and other toxins by decreasing mitochondrial respiration and reactive oxygen species production</article-title><source>Journal of Biological Chemistry</source><volume>287</volume><fpage>37245</fpage><lpage>37258</lpage><pub-id pub-id-type="doi">10.1074/jbc.M112.366195</pub-id><pub-id pub-id-type="pmid">22948140</pub-id></element-citation></ref><ref id="bib59"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>O'Loughlin</surname> <given-names>E</given-names></name><name><surname>Pakan</surname> <given-names>JMP</given-names></name><name><surname>Yilmazer-Hanke</surname> <given-names>D</given-names></name><name><surname>McDermott</surname> <given-names>KW</given-names></name></person-group><year iso-8601-date="2017">2017</year><article-title>Acute in utero exposure to lipopolysaccharide induces inflammation in the pre- and postnatal brain and alters the glial cytoarchitecture in the developing amygdala</article-title><source>Journal of Neuroinflammation</source><volume>14</volume><elocation-id>212</elocation-id><pub-id pub-id-type="doi">10.1186/s12974-017-0981-8</pub-id><pub-id pub-id-type="pmid">29096641</pub-id></element-citation></ref><ref id="bib60"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Oskvig</surname> <given-names>DB</given-names></name><name><surname>Elkahloun</surname> <given-names>AG</given-names></name><name><surname>Johnson</surname> <given-names>KR</given-names></name><name><surname>Phillips</surname> <given-names>TM</given-names></name><name><surname>Herkenham</surname> <given-names>M</given-names></name></person-group><year iso-8601-date="2012">2012</year><article-title>Maternal immune activation by LPS selectively alters specific gene expression profiles of interneuron migration and oxidative stress in the fetus without triggering a fetal immune response</article-title><source>Brain, Behavior, and Immunity</source><volume>26</volume><fpage>623</fpage><lpage>634</lpage><pub-id pub-id-type="doi">10.1016/j.bbi.2012.01.015</pub-id><pub-id pub-id-type="pmid">22310921</pub-id></element-citation></ref><ref id="bib61"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Parboosing</surname> <given-names>R</given-names></name><name><surname>Bao</surname> <given-names>Y</given-names></name><name><surname>Shen</surname> <given-names>L</given-names></name><name><surname>Schaefer</surname> <given-names>CA</given-names></name><name><surname>Brown</surname> <given-names>AS</given-names></name></person-group><year iso-8601-date="2013">2013</year><article-title>Gestational influenza and bipolar disorder in adult offspring</article-title><source>JAMA Psychiatry</source><volume>70</volume><fpage>677</fpage><lpage>685</lpage><pub-id pub-id-type="doi">10.1001/jamapsychiatry.2013.896</pub-id><pub-id pub-id-type="pmid">23699867</pub-id></element-citation></ref><ref id="bib62"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Patterson</surname> <given-names>PH</given-names></name></person-group><year iso-8601-date="2009">2009</year><article-title>Immune involvement in schizophrenia and autism: etiology, pathology and animal models</article-title><source>Behavioural Brain Research</source><volume>204</volume><fpage>313</fpage><lpage>321</lpage><pub-id pub-id-type="doi">10.1016/j.bbr.2008.12.016</pub-id><pub-id pub-id-type="pmid">19136031</pub-id></element-citation></ref><ref id="bib63"><element-citation publication-type="software"><person-group person-group-type="author"><collab>R Development Core Team</collab></person-group><year iso-8601-date="2015">2015</year><data-title>R: a language and environment for statistical computing</data-title><publisher-loc>Vienna, Austria</publisher-loc><publisher-name>R Foundation for Statistical Computing</publisher-name><ext-link ext-link-type="uri" xlink:href="http://www:r-project.org">http://www:r-project.org</ext-link></element-citation></ref><ref id="bib64"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Richetto</surname> <given-names>J</given-names></name><name><surname>Calabrese</surname> <given-names>F</given-names></name><name><surname>Riva</surname> <given-names>MA</given-names></name><name><surname>Meyer</surname> <given-names>U</given-names></name></person-group><year iso-8601-date="2014">2014</year><article-title>Prenatal immune activation induces maturation-dependent alterations in the prefrontal GABAergic transcriptome</article-title><source>Schizophrenia Bulletin</source><volume>40</volume><fpage>351</fpage><lpage>361</lpage><pub-id pub-id-type="doi">10.1093/schbul/sbs195</pub-id><pub-id pub-id-type="pmid">23328159</pub-id></element-citation></ref><ref id="bib65"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Robinson</surname> <given-names>MD</given-names></name><name><surname>McCarthy</surname> <given-names>DJ</given-names></name><name><surname>Smyth</surname> <given-names>GK</given-names></name></person-group><year iso-8601-date="2010">2010</year><article-title>edgeR: a bioconductor package for differential expression analysis of digital gene expression data</article-title><source>Bioinformatics</source><volume>26</volume><fpage>139</fpage><lpage>140</lpage><pub-id pub-id-type="doi">10.1093/bioinformatics/btp616</pub-id><pub-id pub-id-type="pmid">19910308</pub-id></element-citation></ref><ref id="bib66"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Rosenstein</surname> <given-names>JM</given-names></name><name><surname>Krum</surname> <given-names>JM</given-names></name><name><surname>Ruhrberg</surname> <given-names>C</given-names></name></person-group><year iso-8601-date="2010">2010</year><article-title>VEGF in the nervous system</article-title><source>Organogenesis</source><volume>6</volume><fpage>107</fpage><lpage>114</lpage><pub-id pub-id-type="doi">10.4161/org.6.2.11687</pub-id><pub-id pub-id-type="pmid">20885857</pub-id></element-citation></ref><ref id="bib67"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Schafer</surname> <given-names>ST</given-names></name><name><surname>Paquola</surname> <given-names>ACM</given-names></name><name><surname>Stern</surname> <given-names>S</given-names></name><name><surname>Gosselin</surname> <given-names>D</given-names></name><name><surname>Ku</surname> <given-names>M</given-names></name><name><surname>Pena</surname> <given-names>M</given-names></name><name><surname>Kuret</surname> <given-names>TJM</given-names></name><name><surname>Liyanage</surname> <given-names>M</given-names></name><name><surname>Mansour</surname> <given-names>AA</given-names></name><name><surname>Jaeger</surname> <given-names>BN</given-names></name><name><surname>Marchetto</surname> <given-names>MC</given-names></name><name><surname>Glass</surname> <given-names>CK</given-names></name><name><surname>Mertens</surname> <given-names>J</given-names></name><name><surname>Gage</surname> <given-names>FH</given-names></name></person-group><year iso-8601-date="2019">2019</year><article-title>Pathological priming causes developmental gene network heterochronicity in autistic subject-derived neurons</article-title><source>Nature Neuroscience</source><volume>22</volume><fpage>243</fpage><lpage>255</lpage><pub-id pub-id-type="doi">10.1038/s41593-018-0295-x</pub-id><pub-id pub-id-type="pmid">30617258</pub-id></element-citation></ref><ref id="bib68"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Scott</surname> <given-names>CE</given-names></name><name><surname>Wynn</surname> <given-names>SL</given-names></name><name><surname>Sesay</surname> <given-names>A</given-names></name><name><surname>Cruz</surname> <given-names>C</given-names></name><name><surname>Cheung</surname> <given-names>M</given-names></name><name><surname>Gomez Gaviro</surname> <given-names>MV</given-names></name><name><surname>Booth</surname> <given-names>S</given-names></name><name><surname>Gao</surname> <given-names>B</given-names></name><name><surname>Cheah</surname> <given-names>KS</given-names></name><name><surname>Lovell-Badge</surname> <given-names>R</given-names></name><name><surname>Briscoe</surname> <given-names>J</given-names></name></person-group><year iso-8601-date="2010">2010</year><article-title>SOX9 induces and maintains neural stem cells</article-title><source>Nature Neuroscience</source><volume>13</volume><fpage>1181</fpage><lpage>1189</lpage><pub-id pub-id-type="doi">10.1038/nn.2646</pub-id><pub-id pub-id-type="pmid">20871603</pub-id></element-citation></ref><ref id="bib69"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Shi</surname> <given-names>L</given-names></name><name><surname>Fatemi</surname> <given-names>SH</given-names></name><name><surname>Sidwell</surname> <given-names>RW</given-names></name><name><surname>Patterson</surname> <given-names>PH</given-names></name></person-group><year iso-8601-date="2003">2003</year><article-title>Maternal influenza infection causes marked behavioral and pharmacological changes in the offspring</article-title><source>The Journal of Neuroscience</source><volume>23</volume><fpage>297</fpage><lpage>302</lpage><pub-id pub-id-type="doi">10.1523/JNEUROSCI.23-01-00297.2003</pub-id><pub-id pub-id-type="pmid">12514227</pub-id></element-citation></ref><ref id="bib70"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Shin Yim</surname> <given-names>Y</given-names></name><name><surname>Park</surname> <given-names>A</given-names></name><name><surname>Berrios</surname> <given-names>J</given-names></name><name><surname>Lafourcade</surname> <given-names>M</given-names></name><name><surname>Pascual</surname> <given-names>LM</given-names></name><name><surname>Soares</surname> <given-names>N</given-names></name><name><surname>Yeon Kim</surname> <given-names>J</given-names></name><name><surname>Kim</surname> <given-names>S</given-names></name><name><surname>Kim</surname> <given-names>H</given-names></name><name><surname>Waisman</surname> <given-names>A</given-names></name><name><surname>Littman</surname> <given-names>DR</given-names></name><name><surname>Wickersham</surname> <given-names>IR</given-names></name><name><surname>Harnett</surname> <given-names>MT</given-names></name><name><surname>Huh</surname> <given-names>JR</given-names></name><name><surname>Choi</surname> <given-names>GB</given-names></name></person-group><year iso-8601-date="2017">2017</year><article-title>Reversing behavioural abnormalities in mice exposed to maternal inflammation</article-title><source>Nature</source><volume>549</volume><fpage>482</fpage><lpage>487</lpage><pub-id pub-id-type="doi">10.1038/nature23909</pub-id><pub-id pub-id-type="pmid">28902835</pub-id></element-citation></ref><ref id="bib71"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Simões</surname> <given-names>LR</given-names></name><name><surname>Sangiogo</surname> <given-names>G</given-names></name><name><surname>Tashiro</surname> <given-names>MH</given-names></name><name><surname>Generoso</surname> <given-names>JS</given-names></name><name><surname>Faller</surname> <given-names>CJ</given-names></name><name><surname>Dominguini</surname> <given-names>D</given-names></name><name><surname>Mastella</surname> <given-names>GA</given-names></name><name><surname>Scaini</surname> <given-names>G</given-names></name><name><surname>Giridharan</surname> <given-names>VV</given-names></name><name><surname>Michels</surname> <given-names>M</given-names></name><name><surname>Florentino</surname> <given-names>D</given-names></name><name><surname>Petronilho</surname> <given-names>F</given-names></name><name><surname>Réus</surname> <given-names>GZ</given-names></name><name><surname>Dal-Pizzol</surname> <given-names>F</given-names></name><name><surname>Zugno</surname> <given-names>AI</given-names></name><name><surname>Barichello</surname> <given-names>T</given-names></name></person-group><year iso-8601-date="2018">2018</year><article-title>Maternal immune activation induced by lipopolysaccharide triggers immune response in pregnant mother and fetus, and induces behavioral impairment in adult rats</article-title><source>Journal of Psychiatric Research</source><volume>100</volume><fpage>71</fpage><lpage>83</lpage><pub-id pub-id-type="doi">10.1016/j.jpsychires.2018.02.007</pub-id><pub-id pub-id-type="pmid">29494891</pub-id></element-citation></ref><ref id="bib72"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Smith</surname> <given-names>SE</given-names></name><name><surname>Li</surname> <given-names>J</given-names></name><name><surname>Garbett</surname> <given-names>K</given-names></name><name><surname>Mirnics</surname> <given-names>K</given-names></name><name><surname>Patterson</surname> <given-names>PH</given-names></name></person-group><year iso-8601-date="2007">2007</year><article-title>Maternal immune activation alters fetal brain development through interleukin-6</article-title><source>Journal of Neuroscience</source><volume>27</volume><fpage>10695</fpage><lpage>10702</lpage><pub-id pub-id-type="doi">10.1523/JNEUROSCI.2178-07.2007</pub-id><pub-id pub-id-type="pmid">17913903</pub-id></element-citation></ref><ref id="bib73"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Soumiya</surname> <given-names>H</given-names></name><name><surname>Fukumitsu</surname> <given-names>H</given-names></name><name><surname>Furukawa</surname> <given-names>S</given-names></name></person-group><year iso-8601-date="2011">2011</year><article-title>Prenatal immune challenge compromises the normal course of neurogenesis during development of the mouse cerebral cortex</article-title><source>Journal of Neuroscience Research</source><volume>89</volume><fpage>1575</fpage><lpage>1585</lpage><pub-id pub-id-type="doi">10.1002/jnr.22704</pub-id><pub-id pub-id-type="pmid">21732402</pub-id></element-citation></ref><ref id="bib74"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Stolp</surname> <given-names>HB</given-names></name><name><surname>Turnquist</surname> <given-names>C</given-names></name><name><surname>Dziegielewska</surname> <given-names>KM</given-names></name><name><surname>Saunders</surname> <given-names>NR</given-names></name><name><surname>Anthony</surname> <given-names>DC</given-names></name><name><surname>Molnár</surname> <given-names>Z</given-names></name></person-group><year iso-8601-date="2011">2011</year><article-title>Reduced ventricular proliferation in the foetal cortex following maternal inflammation in the mouse</article-title><source>Brain</source><volume>134</volume><fpage>3236</fpage><lpage>3248</lpage><pub-id pub-id-type="doi">10.1093/brain/awr237</pub-id><pub-id pub-id-type="pmid">21964917</pub-id></element-citation></ref><ref id="bib75"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Subramanian</surname> <given-names>A</given-names></name><name><surname>Tamayo</surname> <given-names>P</given-names></name><name><surname>Mootha</surname> <given-names>VK</given-names></name><name><surname>Mukherjee</surname> <given-names>S</given-names></name><name><surname>Ebert</surname> <given-names>BL</given-names></name><name><surname>Gillette</surname> <given-names>MA</given-names></name><name><surname>Paulovich</surname> <given-names>A</given-names></name><name><surname>Pomeroy</surname> <given-names>SL</given-names></name><name><surname>Golub</surname> <given-names>TR</given-names></name><name><surname>Lander</surname> <given-names>ES</given-names></name><name><surname>Mesirov</surname> <given-names>JP</given-names></name></person-group><year iso-8601-date="2005">2005</year><article-title>Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles</article-title><source>PNAS</source><volume>102</volume><fpage>15545</fpage><lpage>15550</lpage><pub-id pub-id-type="doi">10.1073/pnas.0506580102</pub-id><pub-id pub-id-type="pmid">16199517</pub-id></element-citation></ref><ref id="bib76"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Sun</surname> <given-names>W</given-names></name><name><surname>Cornwell</surname> <given-names>A</given-names></name><name><surname>Li</surname> <given-names>J</given-names></name><name><surname>Peng</surname> <given-names>S</given-names></name><name><surname>Osorio</surname> <given-names>MJ</given-names></name><name><surname>Aalling</surname> <given-names>N</given-names></name><name><surname>Wang</surname> <given-names>S</given-names></name><name><surname>Benraiss</surname> <given-names>A</given-names></name><name><surname>Lou</surname> <given-names>N</given-names></name><name><surname>Goldman</surname> <given-names>SA</given-names></name><name><surname>Nedergaard</surname> <given-names>M</given-names></name></person-group><year iso-8601-date="2017">2017</year><article-title>SOX9 is an Astrocyte-Specific nuclear marker in the adult brain outside the neurogenic regions</article-title><source>The Journal of Neuroscience</source><volume>37</volume><fpage>4493</fpage><lpage>4507</lpage><pub-id pub-id-type="doi">10.1523/JNEUROSCI.3199-16.2017</pub-id><pub-id pub-id-type="pmid">28336567</pub-id></element-citation></ref><ref id="bib77"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Szklarczyk</surname> <given-names>D</given-names></name><name><surname>Gable</surname> <given-names>AL</given-names></name><name><surname>Lyon</surname> <given-names>D</given-names></name><name><surname>Junge</surname> <given-names>A</given-names></name><name><surname>Wyder</surname> <given-names>S</given-names></name><name><surname>Huerta-Cepas</surname> <given-names>J</given-names></name><name><surname>Simonovic</surname> <given-names>M</given-names></name><name><surname>Doncheva</surname> <given-names>NT</given-names></name><name><surname>Morris</surname> <given-names>JH</given-names></name><name><surname>Bork</surname> <given-names>P</given-names></name><name><surname>Jensen</surname> <given-names>LJ</given-names></name><name><surname>Mering</surname> <given-names>CV</given-names></name></person-group><year iso-8601-date="2019">2019</year><article-title>STRING v11: protein-protein association networks with increased coverage, supporting functional discovery in genome-wide experimental datasets</article-title><source>Nucleic Acids Research</source><volume>47</volume><fpage>D607</fpage><lpage>D613</lpage><pub-id pub-id-type="doi">10.1093/nar/gky1131</pub-id><pub-id pub-id-type="pmid">30476243</pub-id></element-citation></ref><ref id="bib78"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Takahashi</surname> <given-names>H</given-names></name><name><surname>Shibuya</surname> <given-names>M</given-names></name></person-group><year iso-8601-date="2005">2005</year><article-title>The vascular endothelial growth factor (VEGF)/VEGF receptor system and its role under physiological and pathological conditions</article-title><source>Clinical Science</source><volume>109</volume><fpage>227</fpage><lpage>241</lpage><pub-id pub-id-type="doi">10.1042/CS20040370</pub-id><pub-id pub-id-type="pmid">16104843</pub-id></element-citation></ref><ref id="bib79"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Tan</surname> <given-names>X</given-names></name><name><surname>Liu</surname> <given-names>WA</given-names></name><name><surname>Zhang</surname> <given-names>XJ</given-names></name><name><surname>Shi</surname> <given-names>W</given-names></name><name><surname>Ren</surname> <given-names>SQ</given-names></name><name><surname>Li</surname> <given-names>Z</given-names></name><name><surname>Brown</surname> <given-names>KN</given-names></name><name><surname>Shi</surname> <given-names>SH</given-names></name></person-group><year iso-8601-date="2016">2016</year><article-title>Vascular influence on ventral telencephalic progenitors and neocortical interneuron production</article-title><source>Developmental Cell</source><volume>36</volume><fpage>624</fpage><lpage>638</lpage><pub-id pub-id-type="doi">10.1016/j.devcel.2016.02.023</pub-id><pub-id pub-id-type="pmid">27003936</pub-id></element-citation></ref><ref id="bib80"><element-citation publication-type="software"><person-group person-group-type="author"><collab>Team R</collab></person-group><year iso-8601-date="2018">2018</year><source>RStudio: Integrated Development Environment for R</source><version designator="1.4">1.4</version><ext-link ext-link-type="uri" xlink:href="https://rstudio.com/products/rstudio/">https://rstudio.com/products/rstudio/</ext-link></element-citation></ref><ref id="bib81"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Thion</surname> <given-names>MS</given-names></name><name><surname>Mosser</surname> <given-names>CA</given-names></name><name><surname>Férézou</surname> <given-names>I</given-names></name><name><surname>Grisel</surname> <given-names>P</given-names></name><name><surname>Baptista</surname> <given-names>S</given-names></name><name><surname>Low</surname> <given-names>D</given-names></name><name><surname>Ginhoux</surname> <given-names>F</given-names></name><name><surname>Garel</surname> <given-names>S</given-names></name><name><surname>Audinat</surname> <given-names>E</given-names></name></person-group><year iso-8601-date="2019">2019</year><article-title>Biphasic impact of prenatal inflammation and macrophage depletion on the wiring of neocortical inhibitory circuits</article-title><source>Cell Reports</source><volume>28</volume><fpage>1119</fpage><lpage>1126</lpage><pub-id pub-id-type="doi">10.1016/j.celrep.2019.06.086</pub-id><pub-id pub-id-type="pmid">31365857</pub-id></element-citation></ref><ref id="bib82"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Tsukada</surname> <given-names>T</given-names></name><name><surname>Simamura</surname> <given-names>E</given-names></name><name><surname>Shimada</surname> <given-names>H</given-names></name><name><surname>Arai</surname> <given-names>T</given-names></name><name><surname>Higashi</surname> <given-names>N</given-names></name><name><surname>Akai</surname> <given-names>T</given-names></name><name><surname>Iizuka</surname> <given-names>H</given-names></name><name><surname>Hatta</surname> <given-names>T</given-names></name></person-group><year iso-8601-date="2015">2015</year><article-title>The suppression of maternal-fetal leukemia inhibitory factor signal relay pathway by maternal immune activation impairs brain development in mice</article-title><source>PLOS ONE</source><volume>10</volume><elocation-id>e0129011</elocation-id><pub-id pub-id-type="doi">10.1371/journal.pone.0129011</pub-id><pub-id pub-id-type="pmid">26043040</pub-id></element-citation></ref><ref id="bib83"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Vasistha</surname> <given-names>NA</given-names></name><name><surname>Pardo-Navarro</surname> <given-names>M</given-names></name><name><surname>Gasthaus</surname> <given-names>J</given-names></name><name><surname>Weijers</surname> <given-names>D</given-names></name><name><surname>Müller</surname> <given-names>MK</given-names></name><name><surname>García-González</surname> <given-names>D</given-names></name><name><surname>Malwade</surname> <given-names>S</given-names></name><name><surname>Korshunova</surname> <given-names>I</given-names></name><name><surname>Pfisterer</surname> <given-names>U</given-names></name><name><surname>von Engelhardt</surname> <given-names>J</given-names></name><name><surname>Hougaard</surname> <given-names>KS</given-names></name><name><surname>Khodosevich</surname> <given-names>K</given-names></name></person-group><year iso-8601-date="2020">2020</year><article-title>Maternal inflammation has a profound effect on cortical interneuron development in a stage and subtype-specific manner</article-title><source>Molecular Psychiatry</source><volume>25</volume><fpage>2313</fpage><lpage>2329</lpage><pub-id pub-id-type="doi">10.1038/s41380-019-0539-5</pub-id><pub-id pub-id-type="pmid">31595033</pub-id></element-citation></ref><ref id="bib84"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Walton</surname> <given-names>CC</given-names></name><name><surname>Zhang</surname> <given-names>W</given-names></name><name><surname>Patiño-Parrado</surname> <given-names>I</given-names></name><name><surname>Barrio-Alonso</surname> <given-names>E</given-names></name><name><surname>Garrido</surname> <given-names>JJ</given-names></name><name><surname>Frade</surname> <given-names>JM</given-names></name></person-group><year iso-8601-date="2019">2019</year><article-title>Primary neurons can enter M-phase</article-title><source>Scientific Reports</source><volume>9</volume><elocation-id>4594</elocation-id><pub-id pub-id-type="doi">10.1038/s41598-019-40462-4</pub-id><pub-id pub-id-type="pmid">30872646</pub-id></element-citation></ref><ref id="bib85"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Wang</surname> <given-names>L</given-names></name><name><surname>Wang</surname> <given-names>S</given-names></name><name><surname>Li</surname> <given-names>W</given-names></name></person-group><year iso-8601-date="2012">2012</year><article-title>RSeQC: quality control of RNA-seq experiments</article-title><source>Bioinformatics</source><volume>28</volume><fpage>2184</fpage><lpage>2185</lpage><pub-id pub-id-type="doi">10.1093/bioinformatics/bts356</pub-id><pub-id pub-id-type="pmid">22743226</pub-id></element-citation></ref><ref id="bib86"><element-citation publication-type="book"><person-group person-group-type="author"><name><surname>Wickham</surname> <given-names>H</given-names></name></person-group><year iso-8601-date="2009">2009</year><source>Ggplot2, Elegant Graphics for Data Analysis, Ggplot2</source><publisher-name>Springer</publisher-name><pub-id pub-id-type="doi">10.1007/978-0-387-98141-3</pub-id></element-citation></ref><ref id="bib87"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Wu</surname> <given-names>WL</given-names></name><name><surname>Hsiao</surname> <given-names>EY</given-names></name><name><surname>Yan</surname> <given-names>Z</given-names></name><name><surname>Mazmanian</surname> <given-names>SK</given-names></name><name><surname>Patterson</surname> <given-names>PH</given-names></name></person-group><year iso-8601-date="2017">2017</year><article-title>The placental interleukin-6 signaling controls fetal brain development and behavior</article-title><source>Brain, Behavior, and Immunity</source><volume>62</volume><fpage>11</fpage><lpage>23</lpage><pub-id pub-id-type="doi">10.1016/j.bbi.2016.11.007</pub-id><pub-id pub-id-type="pmid">27838335</pub-id></element-citation></ref><ref id="bib88"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Xie</surname> <given-names>T</given-names></name><name><surname>Stathopoulou</surname> <given-names>MG</given-names></name><name><surname>de Andrés</surname> <given-names>F</given-names></name><name><surname>Siest</surname> <given-names>G</given-names></name><name><surname>Murray</surname> <given-names>H</given-names></name><name><surname>Martin</surname> <given-names>M</given-names></name><name><surname>Cobaleda</surname> <given-names>J</given-names></name><name><surname>Delgado</surname> <given-names>A</given-names></name><name><surname>Lamont</surname> <given-names>J</given-names></name><name><surname>Peñas-LIedó</surname> <given-names>E</given-names></name><name><surname>LLerena</surname> <given-names>A</given-names></name><name><surname>Visvikis-Siest</surname> <given-names>S</given-names></name></person-group><year iso-8601-date="2017">2017</year><article-title>VEGF-related polymorphisms identified by GWAS and risk for major depression</article-title><source>Translational Psychiatry</source><volume>7</volume><elocation-id>e1055</elocation-id><pub-id pub-id-type="doi">10.1038/tp.2017.36</pub-id><pub-id pub-id-type="pmid">28267147</pub-id></element-citation></ref><ref id="bib89"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Zhang</surname> <given-names>J</given-names></name><name><surname>Jing</surname> <given-names>Y</given-names></name><name><surname>Zhang</surname> <given-names>H</given-names></name><name><surname>Bilkey</surname> <given-names>DK</given-names></name><name><surname>Liu</surname> <given-names>P</given-names></name></person-group><year iso-8601-date="2019">2019</year><article-title>Maternal immune activation altered microglial immunoreactivity in the brain of postnatal day 2 rat offspring</article-title><source>Synapse</source><volume>73</volume><elocation-id>e22072</elocation-id><pub-id pub-id-type="doi">10.1002/syn.22072</pub-id><pub-id pub-id-type="pmid">30256454</pub-id></element-citation></ref><ref id="bib90"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Zhang</surname> <given-names>B</given-names></name><name><surname>Horvath</surname> <given-names>S</given-names></name></person-group><year iso-8601-date="2005">2005</year><article-title>A general framework for weighted gene co-expression network analysis</article-title><source>Statistical Applications in Genetics and Molecular Biology</source><volume>4</volume><elocation-id>Article17</elocation-id><pub-id pub-id-type="doi">10.2202/1544-6115.1128</pub-id><pub-id pub-id-type="pmid">16646834</pub-id></element-citation></ref><ref id="bib91"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Zhao</surname> <given-names>Q</given-names></name><name><surname>Wang</surname> <given-names>Q</given-names></name><name><surname>Wang</surname> <given-names>J</given-names></name><name><surname>Tang</surname> <given-names>M</given-names></name><name><surname>Huang</surname> <given-names>S</given-names></name><name><surname>Peng</surname> <given-names>K</given-names></name><name><surname>Han</surname> <given-names>Y</given-names></name><name><surname>Zhang</surname> <given-names>J</given-names></name><name><surname>Liu</surname> <given-names>G</given-names></name><name><surname>Fang</surname> <given-names>Q</given-names></name><name><surname>You</surname> <given-names>Z</given-names></name></person-group><year iso-8601-date="2019">2019</year><article-title>Maternal immune activation-induced PPARγ-dependent dysfunction of microglia associated with neurogenic impairment and aberrant postnatal behaviors in offspring</article-title><source>Neurobiology of Disease</source><volume>125</volume><fpage>1</fpage><lpage>13</lpage><pub-id pub-id-type="doi">10.1016/j.nbd.2019.01.005</pub-id><pub-id pub-id-type="pmid">30659984</pub-id></element-citation></ref></ref-list></back><sub-article article-type="decision-letter" id="sa1"><front-stub><article-id pub-id-type="doi">10.7554/eLife.60100.sa1</article-id><title-group><article-title>Decision letter</article-title></title-group><contrib-group><contrib contrib-type="editor"><name><surname>Bhattacharyya</surname><given-names>Anita</given-names></name><role>Reviewing Editor</role><aff><institution>University of Wisconsin, Madison</institution><country>United States</country></aff></contrib></contrib-group><contrib-group><contrib contrib-type="reviewer"><name><surname>Bhattacharyya</surname><given-names>Anita</given-names> </name><role>Reviewer</role><aff><institution>University of Wisconsin, Madison</institution><country>United States</country></aff></contrib><contrib contrib-type="reviewer"><name><surname>Stevens</surname><given-names>Hanna</given-names> </name><role>Reviewer</role><aff><country>United States</country></aff></contrib></contrib-group></front-stub><body><boxed-text><p>In the interests of transparency, eLife publishes the most substantive revision requests and the accompanying author responses.</p></boxed-text><p><bold>Acceptance summary:</bold></p><p>The manuscript by Canales provides a comprehensive study of the effects of maternal immune activation (MIA) using a rigorous and carefully validated MIA mouse model. A major strength of the manuscript is the investigation of changes across time. The results are of broad interest to the fields of developmental neuroscience and neurodevelopmental disorders, where MIA is a risk factor for schizophrenia and autism.</p><p><bold>Decision letter after peer review:</bold></p><p>Thank you for submitting your article "Sequential perturbations to mouse corticogenesis following in utero maternal immune activation" for consideration by <italic>eLife</italic>. Your article has been reviewed by three peer reviewers, including Anita Bhattacharyya as the Reviewing Editor and Reviewer #1, and the evaluation has been overseen by Kathryn Cheah as the Senior Editor. The following individual involved in review of your submission has agreed to reveal their identity: Hanna Stevens (Reviewer #2).</p><p>The reviewers have discussed the reviews with one another and the Reviewing Editor has drafted this decision to help you prepare a revised submission.</p><p>As the editors have judged that your manuscript is of interest, but as described below that additional experiments are required before it is published, we would like to draw your attention to changes in our revision policy that we have made in response to COVID-19 (https://elifesciences.org/articles/57162). First, because many researchers have temporarily lost access to the labs, we will give authors as much time as they need to submit revised manuscripts. We are also offering, if you choose, to post the manuscript to bioRxiv (if it is not already there) along with this decision letter and a formal designation that the manuscript is "in revision at <italic>eLife</italic>". Please let us know if you would like to pursue this option. (If your work is more suitable for medRxiv, you will need to post the preprint yourself, as the mechanisms for us to do so are still in development.)</p><p>Your manuscript provides a comprehensive study of the effects of maternal immune activation (MIA) using a rigorous and carefully validated MIA model. A major strength of the manuscript is the investigation of changes across time. The results are of broad interest to the fields of developmental neuroscience and neurodevelopmental disorders, where MIA is a risk factor for schizophrenia and autism. Despite these strengths, there are significant problems that limit the potential impact of the results. There are areas in which additional support is needed for the conclusions and the authors should be careful to state conclusions that are supported by specific experimental results.</p><p>Essential revisions:</p><p>(1a) The data in Figure 4 are compelling but do not yet definitively support a conclusion of reduced proliferation at E17.5 following MIA.</p><p>(1b) In addition, the authors should avoid hyperbole when describing these results ("MIA profoundly alters cortical proliferation.")</p><p>(1c) Additional quantification of cell types would bolster the results: The authors should quantify Ki67+ cells in the same manner that they quantified PH3+ cells and analyze expression of the progenitor marker SOX2, which tends to be more restricted to the VZ, in addition to PAX6 and SOX9.</p><p>(1d) In addition to measuring cortical thickness and counting SOX9+ and PH3+ cells, birth-dating can help support the under-proliferation/precocious neurogenesis hypothesis to test whether neural progenitor cells are exiting the cell cycle and differentiating early, migrating out of the VZ early, or both. Alternatively, an analysis of cell cycle kinetics would be useful; delayed mitosis and extended cell cycle time could contribute to the reduction in early-born cell types.</p><p>(1e) In the Discussion, the authors describe transcriptional signatures as "robust". What is the basis of this claim?</p><p>(2a) Additional analyses are needed to support the results of layer specific differences induced by MIA in in Figure 5.</p><p>(2b) Tbr1 thickness is not significantly different, as suggested in the text.</p><p>(2c) The authors should clarify whether Satb2 thickness measured using all Satb2+ layers or just upper layers.</p><p>(2d) A subset of layer 6 neurons expresses Ctip2, so it is hard to determine whether the phenotype represents a reduction in layer 6 or layer 5 Ctip2+ neurons, or a reduction in Ctip2+ neurons in general. The authors should co-label with Tbr1 and additional deep layer markers, such as Satb2 (which is expressed in a subset of layer 5 neurons), Fezf2, Sox5, and/or Rbp4 to bet a better idea of layer 6-layer 5 separation. The results do not rule out possibility that migration may be affected by MIA and the authors should address this point.</p><p>(3a) The Olig2 result suggests that later born cells are over-represented in the MIA cortex at the expense of younger-born cells, which could have a secondary effect on interneuron development. The authors should be careful to interpret this as a difference in progenitor maturation, not neuronal maturation.</p><p>(3b) In addition, other later-born cell type markers, such as Olig1 or Pdgfra, should be analyzed.</p><p>(4) Transcriptional analysis</p><p>(4a) Since the majority of variation is coming from developmental age, more subtle expression changes due to MIA might be missed by PCA and WGCNA.</p><p>(4b) The authors should discuss the definition of the "grey module" more deeply and discuss the fact that it is the only one that correlates with MIA.</p><p>(4c) The results that no DE genes with an FDR < 0.05 at P0 is surprising. The authors should address explanations for why expression appears to normalize suddenly, both biological explanations and technical (e.g., result of noisy or stochastic genes being falsely identified as DEGs).</p><p>(5a) A more nuanced acknowledgment of the different sources of tissue, in particular when it may have a bearing on the interpretations, is important.</p><p>(5b) The focus on cortical development is warranted and borne out in analysis at later timepoints but early timepoint tissue includes more than cortex.</p><p>(5c) RNA sequencing analysis could include tissue type as a covariate (total forebrain for E12.5 and dorsal forebrain for E14.5, E17.5, and P0).</p><p>(5d) Given that protein was taken from total forebrain at E17.5 for western analysis, the legends of the figures that show these results should not call this "cortical extracts" and the high expression of ctip2 in ventral forebrain at E17.5 should inform the reporting/interpretation of the results.</p><p>(5e) The differential development of the hippocampus should be considered. The authors should address whether the primordial hippocampus may be the prominent source of RNA difference and whether pax6 or Sox9 is elevated in that region, determining whether radial glia or newly differentiated astrocytes may be a source in this particular domain. This may reveal an interesting component of the developmental advancement, given the slightly different timing of hippocampal development relative to neocortex.</p><p>(6) The discussion of the enrichment of DE ASD related genes and the relationship to MIA is underdeveloped. What specific conclusions/parallels can be made?</p><p>(7) Additional details are needed to address rigor of analyses:</p><p>(7a) The authors should state explicitly how many biological replicates were used per condition per stage.</p><p>(7b) The authors should add a clustering analysis in Figure 1 (or move the PCA to Figure 1) to help make it more explicit.</p><p>(7c) Perhaps a uMAP or tSNE rather than simple PCA?</p><p>(7d) For immunostaining the authors should state whether the samples taken from the same litters as those samples used for RNA sequencing.</p><p>(7e) Parameters for cellular data and image analysis are not defined. No details are provided about number of sections and cells that were counted with ImageJ and how chosen fields were randomized.</p></body></sub-article><sub-article article-type="reply" id="sa2"><front-stub><article-id pub-id-type="doi">10.7554/eLife.60100.sa2</article-id><title-group><article-title>Author response</article-title></title-group></front-stub><body><disp-quote content-type="editor-comment"><p>Essential revisions:</p><p>(1a) The data in Figure 4 are compelling but do not yet definitively support a conclusion of reduced proliferation at E17.5 following MIA.</p></disp-quote><p>We thank the reviewers for suggesting experiments to further elucidate proliferative phenotypes. As described below, we performed new IHC experiments for SOX2 and Ki67/PH3 labeling from the same animals and re-analyzed our initial results as recommended. The revised version of the manuscript better characterizes progenitor cells and provides stronger evidence of reduced mitosis at E17.5 following MIA (See response 1c below).</p><disp-quote content-type="editor-comment"><p>(1b) In addition, the authors should avoid hyperbole when describing these results ("MIA profoundly alters cortical proliferation.")</p></disp-quote><p>We removed the word “profoundly”. The title of the figure now reads:</p><p>“Figure 4. Altered cortical proliferation and lamination dynamics in fetal offspring five days after induction of MIA.”</p><disp-quote content-type="editor-comment"><p>(1c) Additional quantification of cell types would bolster the results: The authors should quantify Ki67+ cells in the same manner that they quantified PH3+ cells and analyze expression of the progenitor marker SOX2, which tends to be more restricted to the VZ, in addition to PAX6 and SOX9.</p></disp-quote><p>Following the review request, we quantified Ki67+ cells and performed new SOX2 IHC. The revised results are consistent with reduced progenitor cell findings we reported in the initial submission, with the new data improving the strength and rigor of this analysis. Specifically, we found significantly reduced Ki67+ cell counts and reduced VZ-specific Sox2<sup>+</sup> cell counts in in MIA-exposed group. The new data has been incorporated into Figure 4 and is described in the revised text below:</p><p>“To characterize the neuroanatomical specificity of Blue and Turquoise dysregulation transcriptional signatures (Figure 4D), we performed a static immuno-fluorescent labeling comparison of late neurogenesis, cortical lamination and cell specification at E17.5 in an independent MIA cohort. E17.5 represents the end of the peak of neurogenesis during cortical development, before the neuronal-glial transition that occurs around E18.5. […] Reduced expression of SOX2, a progenitor marker that is restricted to the VZ, was also observed in MIA brains (SOX2 VZ, P=0.0006 two tailed Student’s t-test) (Figure 4F)”.</p><p>“To determine if the decrease in progenitor markers might coincide with decreased proliferation, we next labeled cells undergoing proliferation, using the marker KI67 (any stage of cell proliferation) and the Phospho Histone 3 (PH3) marker, which labels a subset of cycling cells in mitosis (M) phase (Hendzel et al., 1997; Walton et al., 2019). […] Moreover, within the reduced progenitor cell populations at E17.5, the number of cells in the mitosis stage of the cell cycle are decreased to a greater extent.”</p><disp-quote content-type="editor-comment"><p>1d) In addition to measuring cortical thickness and counting SOX9+ and PH3+ cells, birth-dating can help support the under-proliferation/precocious neurogenesis hypothesis to test whether neural progenitor cells are exiting the cell cycle and differentiating early, migrating out of the VZ early, or both. Alternatively, an analysis of cell cycle kinetics would be useful; delayed mitosis and extended cell cycle time could contribute to the reduction in early-born cell types.</p></disp-quote><p>We agree that experiments such as birth dating, analysis of cell cycle timing and exit, and analysis of other time-points across the window of neurogenesis will be needed to ultimately resolve the impacts of MIA across neurogenesis and gliogenesis. Regrettably, such in depth analysis is beyond the scope of this manuscript. Nonetheless, we were able to follow the alternative strategy proposed in the review summary, providing an initial perspective on cell cycle kinetics via comparing the ratio of overall cycling cell numbers (Ki67+) with cells undergoing mitosis and labeled by Phospho Histone 3 (PH3), which is present in M phase. In addition to an overall decrease in VZ progenitor cell populations, the reduced Ki67/PH3 ratio revealed a significantly decreased proportion of progenitor cells undergoing mitosis in MIA brains at E17.5. This demonstrates that cell cycle kinetics and mitosis are altered beyond a simple reduction in progenitor cells at E17.5. The altered ratio is consistent with reduced time in M phase and/or a decreased number of cycling progenitors. This result is exciting and extends our findings. In addition to describing the new finding, we additionally acknowledge in the Discussion section that birth-dating, cell cycle kinetics, and time course experiments are needed to comprehensively determine altered proliferative and neurogenesis dynamics following MIA. The new experiments and analysis are included in Figure 4. The following text has been modified in the revised manuscript. To avoid reproduction, see text in [1c] as well as below.</p><p>“Although these findings are indicative of overall decreased neurogenesis by E17.5, specific birth-dating and time course experiments are needed to test whether progenitor cells exit the cell cycle differently and if overall neurogenesis is perturbed beyond the E17.5 time-point reported here.”</p><disp-quote content-type="editor-comment"><p>(1e) In the Discussion, the authors describe transcriptional signatures as "robust". What is the basis of this claim?</p></disp-quote><p>We describe the transcriptional signatures as “robust” because at E14.5 and E17.5 a large number of genes in our differential expression analysis surpass stringent statistical significance of FDR < 0.05, and substantial fold change (FC) difference, with many genes exceeding 2 FC (log2 = 1) or even 4 FC (log2 = 2) (Figure 1B-D). We acknowledge that the term robust is non-specific, but thought it captured the overall results. We are happy to substitute for another word or drop if the reviewers prefer, but we left this phrasing for now.</p><disp-quote content-type="editor-comment"><p>(2a) Additional analyses are needed to support the results of layer specific differences induced by MIA in in Figure 5.</p></disp-quote><p>We thank the reviewers for the suggested revisions to better define lamination phenotypes. In addition to clarifying the methods and more accurately describing the results, we performed new IHC experiments with TBR1 and SATB2 co-staining, as requested, and re-analyzed cell count and layer thickness data for consistency. The new IHC results further define differences in in MIA-exposed fetal cortex and provide evidence that the altered lamination is associated with differences in neuronal migration.</p><disp-quote content-type="editor-comment"><p>(2b) Tbr1 thickness is not significantly different, as suggested in the text.</p></disp-quote><p>The thickness of the Tbr1+ layer is not significantly decreased but the number of Tbr1+ cells is different. We apologize for this error. We made the following changes to the text:</p><p>“The number of TBR1+ and CTIP2+ cells, but not SATB2+ cells, were also significantly reduced in the neocortex from MIA offspring (TBR1+ cells, P = 0.0138; CTIP2+ cells, P = 0.0119; SATB2+ cells, P = 0.306, two tailed Student t-test) (Figure 5E)”</p><p>And:</p><p>“We found a significant reduction in CTIP2 layer thickness and a similar trend in TBR1 layer thickness in the MIA brains.”</p><disp-quote content-type="editor-comment"><p>(2c) The authors should clarify whether Satb2 thickness measured using all Satb2+ layers or just upper layers.</p></disp-quote><p>For consistency, all thickness analyses (including SATB2) were performed by measuring all positive cells for each marker. We clarify how these layer specific thickness analyses were performed when introducing these assays in the Results and added details to the Materials and methods section under “Gross anatomy, immunohistochemistry and image quantification analysis”:</p><p>“Examination of the distribution of the lamination markers TBR1 (layer 6), CTIP2 (layer 5 and a subset of cells in layer 6), SATB2 (layers 2-4, and a subset of cells in layer 5), and CUX1 (layers 2 and 3), revealed MIA-induced cortical layer-specific alterations in lamination, thickness, and cell density at E17.5 (Figure 5A-B, E-F).”</p><p>“Laminar thickness analyses were performed by measuring the spread of all positive cells within the cortical plate per individual marker.”</p><p>In Figure 5 legend: “(F) laminar thickness analysis, based on the spread of all labeled cells, are shown.”</p><disp-quote content-type="editor-comment"><p>(2d) A subset of layer 6 neurons expresses Ctip2, so it is hard to determine whether the phenotype represents a reduction in layer 6 or layer 5 Ctip2+ neurons, or a reduction in Ctip2+ neurons in general. The authors should co-label with Tbr1 and additional deep layer markers, such as Satb2 (which is expressed in a subset of layer 5 neurons), Fezf2, Sox5, and/or Rbp4 to bet a better idea of layer 6-layer 5 separation. The results do not rule out possibility that migration may be affected by MIA and the authors should address this point.</p></disp-quote><p>We are grateful for the reviewer’s point that Ctip2 is not a definitive marker distinguishing layer 5 or 6, and we have clarified this in the Figure 5 legend. As requested, we performed new experiments co-staining for TBR1 and SATB2 to further investigate the layer 5/6 lamination phenotype via labeling SATB2+ putative upper layer cells that form layers 2-4, as well as some Layer 5 cells. These experiments show that SATB2+ cells exhibit similar differences in positioning and differential overlap with Tbr1+ cells as Ctip2+ cells. These findings indicate generalized perturbation to lamination, rather than a specific impact on Ctip2+ layer 5/6 cells. Our new experiments show that SATB2+ cells in the MIA group are already in their more superficial position and show a more rounded morphology, indicating a more complete stage of lamination and further progression of migration of SATB2+ cells. These findings provide further evidence of lamination defects, and are consistent with our general interpretation that MIA alters developing cortical anatomy.</p><p>Of interest and relevant to the reviewer comment on migration, the new co-staining experiments also provide evidence of alterations in migration. By analyzing higher magnification images, we identified differences in cell body shape of both TBR1 and SATB2 positive cells between the MIA and control groups, indicative of an active radial migration state observable in control animals versus a more static rounded cell body shape in the MIA brains (shown in Figure 5D). This finding is consistent with our hypothesis that timing of lamination is altered following MIA exposure, with greater number of cells that were still migrating evident in the control versus MIA group.</p><p>These new TBR1+/SATB2+ co-labeling experiments are presented in Figure 5C-D and described in the Results and Discussion sections as follows:</p><p>“Examination of the distribution of the lamination markers TBR1 (layer 6), CTIP2 (layer 5 and a subset of cells in layer 6), SATB2 (layers 2-4, and a subset of cells in layer 5), and CUX1 (layers 2 and 3), revealed MIA-induced cortical layer-specific alterations in lamination, thickness, and cell density at E17.5 (Figure 5A-B, E-F). […] In contrast to previous reports (Choi et al., 2016), we did not find obvious cortical dysplasia or gross alterations in brain anatomy in the analyzed rostral and caudal regions, including somatosensory cortex and S1DZ (Figure 5E-F).”</p><disp-quote content-type="editor-comment"><p>(3a) The Olig2 result suggests that later born cells are over-represented in the MIA cortex at the expense of younger-born cells, which could have a secondary effect on interneuron development. The authors should be careful to interpret this as a difference in progenitor maturation, not neuronal maturation.</p></disp-quote><p>We agree that the cell-type specific impacts could be either primary or secondary, and are likely to be driven by complex interactions between cell-autonomous and non-cell autonomous signaling. Examining the impact of MIA on specific neuronal and glial populations (including interneurons) has been a focus of other published studies that we reference. We hope our findings spur further future work. We also agree that this should be referred to as progenitor maturation, and updated the text accordingly. The text below from the Discussion now addresses this point:</p><p>“At this developmental stage, our data suggest that MIA induces an acute expenditure of the neocortical, and potentially other, progenitor pools, which in turn leads to an expedited maturation of cells that would normally be generated later during development. […] Although these findings are indicative of overall decreased neurogenesis by E17.5, specific birth-dating and time course experiments are needed to test whether progenitor cells exit the cell cycle differently and if overall neurogenesis is perturbed beyond the E17.5 time-point reported here.”</p><disp-quote content-type="editor-comment"><p>(3b) In addition, other later-born cell type markers, such as Olig1 or Pdgfra, should be analyzed.</p></disp-quote><p>Following the reviewer’s recommendation, we added new IHC experiments staining for PDGRFa. We include quantification analysis by zones for this new marker, following the same strategy we used for DLX2 in the initial submission. PDGFRa is specific to immature oligodendrocytes and showed similar over-representation as seen in our original OLIG2 IHC analysis. These new experiments are presented as Figure 6E, F, Ll. The paragraph in the text that describes these new findings is below:</p><p>“To investigate other cell populations that may be impacted by MIA during development, we analyzed astrocytes, oligodendrocytes and GABAergic interneuron cell populations. […] Thus, these experiments provide evidence that MIA leads to increased numbers of additional cell types (putative astrocytes and oligodendrocytes) that are typically born at later ages during corticogenesis, reinforcing findings of accelerated cell development that is supported by an increased representation of later-born cell types at E17.5 in brains exposed to MIA.”</p><disp-quote content-type="editor-comment"><p>(4) Transcriptional analysis</p><p>(4a) Since the majority of variation is coming from developmental age, more subtle expression changes due to MIA might be missed by PCA and WGCNA.</p></disp-quote><p>This concern seems to be based on some misunderstanding of our approach to identify and characterize MIA transcriptional signatures. We fully agree that MIA-associated differential expression signatures will not be adequately captured by PCA and WGCNA. While PCA and WGCNA are powerful approaches, they are not used as primary methods for differential expression analysis. Indeed, we use these tools as secondary in the context of contextualizing MIA-associated differential expression (DE) signatures at a systems level. For primary DE analysis and identification of DE genes, we used a well-established method that maximizes power via controlling for co-variates, the generalized linear model framework implemented in edgeR. We have added text to clarify this strategy:</p><p>“For transcriptomic analysis, we used a strategy of first defining DE genes at each time-point using the general linear model (GLM) approach implemented in edgeR, followed by mapping DE signatures to systems-level expression patterns via module assignment using weighted gene co-expression network analysis (WGCNA)”</p><p>As described in the manuscript, we applied a pairwise DE testing strategy, examining RNA-seq each time-point independently and including technical batch and sex as covariates in the GLM. We also tested a full model, including all samples with time-point as a covariate, but feel this model inflated significance values for some genes that were DE across the time course and missed stage-specific changes. We also performed DE analysis on sex-stratified samples, i.e. independently on males and females, as reported in the manuscript. This DE analysis, gene ontology enrichment analysis, and modeling of the developmental timing are reported in the Results and figures (Figure 1D-F, Figure 4C).</p><p>Our use of PCA was for initial insights into overall variance across all samples, as is common in RNA-seq studies. This dimensionality reduction approach is useful for detecting technical irregularities (quality control) and for visualizing major components of variation across samples. In this analysis, PC1 and PC2 captured variation associated with developmental age, with some further stratification between MIA and saline samples. We now include these PCA plots as Figure 1B, as requested. We later took advantage of the utility of PC1 and PC2 for capturing developmental age to build an age-predictive model to statistically test the observation that MIA samples at E14.5 and E17.5 exhibit differences in predicted age, as reported in the Results and presented in Figure 4. PCA itself was not used to determine or test differential expression.</p><p>Finally, we used weighted gene co-expression network analysis (WGCNA) as a method used for identifying module patterns and contextualizing and grouping DE genes at the systems/module level. We describe how DE genes map to WGCNA co-expression modules in Figure 2 and focus on DE genes and specific modules in Figures 3-5. This approach enabled association of DE genes to major acute and secondary transcriptional changes. Overall, our paired use of DE (here via GLM) and WGCNA is a standard approach in the field to link gene-level DE changes with systems-level signatures in transcriptomics data.</p><disp-quote content-type="editor-comment"><p>(4b) The authors should discuss the definition of the "grey module" more deeply and discuss the fact that it is the only one that correlates with MIA.</p></disp-quote><p>We first want to clarify that Grey and Green modules are both significantly associated with the MIA phenotype, and YeMaBl as well demonstrates marginal significance. Thus, the review comment that the Grey module was the only MIA-associated module is incorrect. This was described in the following section of the Results in the initial manuscript: “The WGCNA-resolved modules enabled mapping of DE signatures to neurodevelopmental processes that were acutely induced by MIA. There was significant association between genes in Grey (P = 0.007) and Green (P = 0.03) modules and MIA treatment, and marginal significance with MIA for the YeMaBl (P = 0.05) module (Figure 2—figure supplement 1B).” We go into depth interrogating these three MIA associated modules in the Results text, in Figure 3, and in the Discussion.</p><p>Nonetheless, we acknowledge the general point that it is interesting that the Grey module was associated with MIA at the trait level considering this module is made up of the unassigned genes that did not exhibit strong co-expression with the other modules. We explored this finding during our initial analysis to ensure that this was not due to WGCNA parameters. In other words, we tested if the MIA-associated genes in Grey would be assigned to a specific new module if we changed the module assignment parameters. When we increased the number of modules, the composition of Grey shifted, but the association with MIA remained. Further, increasing the number of modules resulted in over-splitting and some reduced biological coherence of individual modules. We thus elected to use WGCNA parameters that generated modules that effectively captured general trends of developmental expression as well as MIA-associated signatures, as reported in the manuscript. In response to the reviewer’s point, we added text to the Discussion regarding the Grey-MIA association:</p><p>“Significant association of the Grey unassigned module with the MIA phenotype was driven by a subset of genes that showed strong acute induction at E12.5 and E14.5. […] Thus, we opted for reduced splitting to capture major co-expression trends.”</p><disp-quote content-type="editor-comment"><p>(4c) The results that no DE genes with an FDR < 0.05 at P0 is surprising. The authors should address explanations for why expression appears to normalize suddenly, both biological explanations and technical (e.g., result of noisy or stochastic genes being falsely identified as DEGs).</p></disp-quote><p>We agree that it was somewhat surprising that the P0 DE signatures were weaker at the gene level compared to the other time-points. This result does not appear to be due to our experimental approach, since the P0 time-point was based on similar sample numbers, sex distribution, batch complexity similar, sequencing depth, and general quality metrics as for the other time-points in our study. In our original manuscript, we speculated that the reduced signature at P0 supports that the significant DE signatures at the earlier time-points reflect transient perturbation at the transcriptional level. As we found strong reproducible relationships between DE signatures at E12.5 and E14.5 and between E14.5 and E17.5, we believe that the MIA signatures of these time-points are unlikely to be false positives. As noted in the review, we agree with the possibility that some unknown technical issue or that general increased biological variability could both reduce the power to detect P0 DE in our experiments. Pursuing the time course approach beyond P0 in future work will be useful to understand the continuing effects, but such further investigation is beyond the scope of this manuscript. We updated the relevant sentence in the Discussion to incorporate suggestions made in the reviewer critique:</p><p>“While significantly subtler in overall effects compared to earlier time-points, we also see evidence for altered neuronal transcription and lasting immune activation at P0. […] Nonetheless, the subtle changes in gene expression at P0 are consistent with deficits in synapse formation that require changes in expression of immune molecules in neurons from newborn MIA offspring (Elmer et al., 2013).”</p><disp-quote content-type="editor-comment"><p>(5a) A more nuanced acknowledgment of the different sources of tissue, in particular when it may have a bearing on the interpretations, is important.</p></disp-quote><p>We acknowledge the importance of communicating our dissection strategy and source of tissue in order to clearly interpret RNA-seq results. We attempted to clarify the methods and add discussion to address the relevant points below.</p><disp-quote content-type="editor-comment"><p>(5b) The focus on cortical development is warranted and borne out in analysis at later timepoints but early timepoint tissue includes more than cortex.</p></disp-quote><p>First, we apologize that the description of our dissection methods was not specific enough and have clarified the Materials and methods text in the revised manuscript (see below). We have also made clear in the Discussion that some subpallial tissues may have been collected in addition to cortex at the youngest age examined, and that this difference in approach may drive some aspects of the differential expression signature (see text below). Finally, we clarify that our strategy at E12.5 was to collect mainly dorsal telencephalon/pallium, but that our approach likely also captured some ventral telencephalon/sub-pallial structures as well as superficial tissue.</p><p>“Tissue dissections for RNA-seq and Immunoblotting. Dissections of dorsal telencephalon/pallium or developing cortex were performed at E12.5 + 6 hrs, E14.5, E17.5, and P0 from control (saline) and poly(I:C)-exposed MIA litters. At E12.5, due to the smaller embryo size, dissections were made via a diagonal cut that maximized collection of dorsal telencephalon, but could have also included some subpallium/ventral telencephalon and superficial tissues. […] Independent litters were used for RNA-seq and immunoblotting.”</p><p>“Finally, bulk tissue RNA-seq is limited in cell type-specific resolution and signatures could additionally be in part driven by cortex-adjacent structures, such as primordial hippocampus, that may have also been captured at the youngest age. Future work using single cell approaches is needed to resolve the impacted cell types and compare transcriptional pathology following MIA.”</p><disp-quote content-type="editor-comment"><p>(5c) RNA sequencing analysis could include tissue type as a covariate (total forebrain for E12.5 and dorsal forebrain for E14.5, E17.5, and P0).</p></disp-quote><p>All differential expression analysis was performed independently for each developmental stage via comparison of matched MIA and control samples. Thus, all samples considered in each DE GLM analysis and the same tissue make up and same developmental stage, precluding the inclusion of any covariate. Within stage, dissections were performed blinded to condition and variation across samples should be random and thus averaged out across replicates. We hope this clarifies the method and addresses this comment. To clarify this in the manuscript, we describe these methods as follows:</p><p>“Tagwise dispersion estimates were generated, and differential expression analysis was performed with edgeR using a generalized linear model, separately for each developmental time-point, including sex, sequencing run factor, and treatment as the variable for testing.”</p><disp-quote content-type="editor-comment"><p>(5d) Given that protein was taken from total forebrain at E17.5 for western analysis, the legends of the figures that show these results should not call this "cortical extracts" and the high expression of ctip2 in ventral forebrain at E17.5 should inform the reporting/interpretation of the results.</p></disp-quote><p>We apologize for this error in the methods. The dissection strategy for RNA-seq and immunoblotting sample collection at E17.5 was identical and only includes developing cortex. Thus RNA and protein quantitative analyses were performed on tissues collected using the same methods. We reorganized the Materials and methods section to highlight this by changing the subsection to the subtitle “Dissection methods for RNA-seq and immunoblotting.” We also removed the term “cortical extracts” from the legend for Figure 4—figure supplement 1A. See response text in [5b] for specific wording.</p><disp-quote content-type="editor-comment"><p>(5e) The differential development of the hippocampus should be considered. The authors should address whether the primordial hippocampus may be the prominent source of RNA difference and whether pax6 or Sox9 is elevated in that region, determining whether radial glia or newly differentiated astrocytes may be a source in this particular domain. This may reveal an interesting component of the developmental advancement, given the slightly different timing of hippocampal development relative to neocortex.</p></disp-quote><p>These questions regarding hippocampal contributions are interesting. Considering the dissection strategy of targeting pallium (due to challenge of isolating cortex), the primordial hippocampus was included and differences in this structure could contribute to our differential expression findings. While there may be some impact from structures such as early hippocampus (or alternatively ventral telencephalon as raised earlier), we believe the general findings do strongly implicate the pallium/cortex as the major source of DE signatures and importantly, our IHC results verified phenotypes in the developing cortex. Future studies examining region-specific differences across non-cortical structures are certainly of interest, though we believe they are outside the scope of our study. Nonetheless, we recognize the limits of manual dissection and bulk tissue strategies, and have added a clarifying sentence into the Discussion – see [5b] response text.</p><disp-quote content-type="editor-comment"><p>(6) The discussion of the enrichment of DE ASD related genes and the relationship to MIA is underdeveloped. What specific conclusions/parallels can be made?</p></disp-quote><p>We added text to the Discussion to further discuss this relationship. We note that the overlap was general, included both up- and down-regulated DE genes, and that overlapping ASD risk genes were annotated across biological pathways rather than specific to any particular function or process. Thus, while we think this is indicative of general translational relevance, we don’t see a clear conclusion beyond the general overlap implicating shared processes across corticogenesis. Thus, we tried to offer some further thoughts, but want to avoid speculation too far beyond the evidence. If the reviewers had further suggestions for how to summarize, we are open to amending.</p><p>“Among MIA-associated DE up-regulated and down-regulated genes, there was strong enrichment for ASD-relevant loci, indicating the general translational relevance of our MIA model to ASD.”</p><p>“The majority of the ASD-relevant genes that were DE following MIA (Figure 1E) were associated with perturbation to the modules associated with neurodevelopmental processes. […] These overlapping genes represent targets of interest for future studies investigating parallel mechanisms between genetic and MIA-associated NDD risk.”</p><disp-quote content-type="editor-comment"><p>(7) Additional details are needed to address rigor of analyses:</p><p>(7a) The authors should state explicitly how many biological replicates were used per condition per stage.</p></disp-quote><p>In the Results section we added the following references to the supplementary files, which provide details of the sample size split by sex (Supplementary file 5), as well as a complete RNA-seq metadata spreadsheet (Supplementary file 20). The following underlined sentence detailing our sample size was also added to the Results section:</p><p>“RNA-seq datasets were generated from male and female embryos from 28 independent litters across control and MIA groups, typically with 1-3 embryos represented per independent litter (Supplementary file 20, Supplementary file 5). […] For transcriptomic analysis, we used a strategy of first defining DE genes at each time-point using the general linear model (GLM) approach implemented in edgeR, followed by mapping DE signatures to systems-level expression patterns via module assignment using weighted gene co-expression network analysis (WGCNA).”</p><disp-quote content-type="editor-comment"><p>(7b) The authors should add a clustering analysis in Figure 1 (or move the PCA to Figure 1) to help make it more explicit.</p></disp-quote><p>PCA plots are now included in Figure 1.</p><disp-quote content-type="editor-comment"><p>(7c) Perhaps a uMAP or tSNE rather than simple PCA?</p></disp-quote><p>We respectfully note that we report results of a bulk RNA-seq experiment and use the PCA solely to capture major components of variation. uMAP and tSNE are more sophisticated dimensionality reduction techniques used primarily to deal with visualization of complex, and most commonly single cell, relationships. With our simpler data, we are unclear what the benefit would be of using these methods for further exploration of sample relationships. Perhaps incorporating PCA plots into Figure 1, as requested, should clarify how the PCA results were used and interpreted and address this point.</p><disp-quote content-type="editor-comment"><p>(7d) For immunostaining the authors should state whether the samples taken from the same litters as those samples used for RNA sequencing.</p><p>(7e) Parameters for cellular data and image analysis are not defined. No details are provided about number of sections and cells that were counted with ImageJ and how chosen fields were randomized.</p></disp-quote><p>We thank the reviewers for bringing this to our attention, we now realize that although all numbers were included for all experiments in our original submission, the presentation of this information lacked consistency and could be confusing. The following modifications were made to address this concern in our revised manuscript:</p><p>First, we further clarified the samples used for IHC in two new supplementary files. We note that since we attempted to use similar tissue sections for consistency, there was a limit on the number of analyses that could be done on individual samples. Further, we used some new litters/embryos for the experiments in the revision. Thus, we provide a simple summary in Supplementary file 21 including sample numbers (N) used for neuroanatomical IHC validation at E17.5, where numbers represent samples analyzed for at least one of the individual listed markers or co-stains. Details of specific use of individual samples is reported in Supplementary file 22.</p><p>Second, we added missing information regarding the imaging quantification methods. We apologize for this omission. We have now added details of parameters used for cellular counting data to the Materials and methods section and clarified that the analysis was performed in triplicate sections. Please see below our modified text, which is now part of the revised Materials and methods subsection “Gross anatomy, immunohistochemistry and image quantification analyses”:</p><p>“Histology was performed in triplicate sections in brains using male and female pups from at least two litters. Sections were selected to focus on the developing somatosensory cortex. Details on sample individual sample usage for each IHC experiment is detailed in Supplementary files 21 and 22.”</p><p>“For cell counting, boxes were drawn on specific areas of the neocortex for progenitor and mitosis markers; SOX2, SOX9 and PH3. This box was 300 x 300 pixels on an image of the neocortex obtained using a 20x objective. To be consistent, the ventral aspect of the box was drawn along the VZ and one corner touched the medial aspect of the VZ. The same strategy was used to draw boxes from Ki67 images except the box was increased in size to 500 x 500 pixels to account for the additional Ki67+ cells seen in neocortical areas above the SVZ. SOX9 boxes were drawn along the dorsal border edge (600 pixels in width) of the neocortex and extended ventrally two thirds the length of the neocortex in order to avoid the VZ/SVZ region. For laminar zone counts, a box spanning 300 pixels was drawn along the VZ and the dorsal/ventral proportions were equally divided into 3 zones. These box sizes allowed for a sufficient number of cells to be counted from each tissue. Laminar thickness analyses were performed by measuring the spread of all positive cells within the cortical plate per individual marker. For all quantification analyses, 1-3 sections per brain utilizing the developing hippocampus as a landmark for positioning were used.”</p></body></sub-article></article>