Functional annotation of RNAseq results
We performed a functional group enrichment of differentially expressed genes (DEGs) in mPFC using the DAVID Annotation Tool (http://david.abcc.ncifcrf.gov)(Jiao et al., 2012) considering GO (Gene Ontology) biological processes and KEGG pathways (Kyoto Encyclopedia of Genes and Genomes). Then, we searched for over‐represented transcription factor‐binding sites (TFBS) using the information of MsigDB (https://www.gsea-msigdb.org/gsea/msigdb) integrated on WebGestalt2019 (http://www.webgestalt.org/)(Liao et al., 2019), and the default parameters, applying the weighted set cover method to reduce redundancy. In both analyses, the Benjamin–Hochberg procedure was performed for multiple testing. Finally, we investigated the existence of gene networks with Ingenuity Pathway Analysis 8.8 software (IPA, http://www.ingenuity.com/products/ipa; Ingenuity Systems, Redwood City, CA, USA) (Krämer et al., 2014) after selecting genes with fold-change >|1.2|.