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|.