Skip to main page content
U.S. flag

An official website of the United States government

Dot gov

The .gov means it’s official.
Federal government websites often end in .gov or .mil. Before sharing sensitive information, make sure you’re on a federal government site.

Https

The site is secure.
The https:// ensures that you are connecting to the official website and that any information you provide is encrypted and transmitted securely.

Access keys NCBI Homepage MyNCBI Homepage Main Content Main Navigation
. 2020 Oct;375(1):210-222.
doi: 10.1124/jpet.120.266171. Epub 2020 Jul 13.

Sex- and Region-Specific Differences in the Transcriptomes of Rat Microglia from the Brainstem and Cervical Spinal Cord

Affiliations

Sex- and Region-Specific Differences in the Transcriptomes of Rat Microglia from the Brainstem and Cervical Spinal Cord

Andrea C Ewald et al. J Pharmacol Exp Ther. 2020 Oct.

Abstract

The neural control system underlying breathing is sexually dimorphic with males being more vulnerable to dysfunction. Microglia also display sex differences, and their role in the architecture of brainstem respiratory rhythm circuitry and modulation of cervical spinal cord respiratory plasticity is becoming better appreciated. To further understand the molecular underpinnings of these sex differences, we performed RNA sequencing of immunomagnetically isolated microglia from brainstem and cervical spinal cord of adult male and female rats. We used various bioinformatics tools (Gene Ontology, Kyoto Encyclopedia of Genes and Genomes, Reactome, STRING, MAGICTRICKS) to functionally categorize identified gene sets, as well as to pinpoint common transcriptional gene drivers that may be responsible for the observed transcriptomic differences. We found few sex differences in the microglial transcriptomes derived from the brainstem, but several hundred genes were differentially expressed by sex in cervical spinal microglia. Comparing brainstem and spinal microglia within and between sexes, we found that the major factor guiding transcriptomic differences was central nervous system (CNS) location rather than sex. We further identified key transcriptional drivers that may be responsible for the transcriptomic differences observed between sexes and CNS regions; enhancer of zeste homolog 2 emerged as the predominant driver of the differentially downregulated genes. We suggest that functional gene alterations identified in metabolism, transcription, and intercellular communication underlie critical microglial heterogeneity and sex differences in CNS regions that contribute to respiratory disorders categorized by dysfunction in neural control. These data will also serve as an important resource data base to advance our understanding of innate immune cell contributions to sex differences and the field of respiratory neural control. SIGNIFICANCE STATEMENT: The contributions of central nervous system (CNS) innate immune cells to sexually dimorphic differences in the neural circuitry controlling breathing are poorly understood. We identify key transcriptomic differences, and their transcriptional drivers, in microglia derived from the brainstem and the C3-C6 cervical spinal cord of healthy adult male and female rats. Gene alterations identified in metabolism, gene transcription, and intercellular communication likely underlie critical microglial heterogeneity and sex differences in these key CNS regions that contribute to the neural control of breathing.

PubMed Disclaimer

Figures

Fig. 1.
Fig. 1.
Female and male brainstem microglia share similar transcriptomes. (A) Principle component analysis of male and female brainstem (BS) microglia (n = 5/treatment). (B) Venn diagram showing the total numbers of shared and unique genes between males and females. There were 13,517 shared genes that did not significantly differ in either sex. There were 28 differentially expressed genes that significantly differed (FDR < 0.05) between females and males. Upregulated genes are shown in red, and downregulated genes are shown in blue. (C) Significantly upregulated genes (12) in female brainstem microglia vs. males. Genes highlighted in pink are X chromosome–linked. (D) Significantly downregulated genes (16) in female brainstem microglia vs. males. Genes highlighted in pink are X chromosome–linked, and those in blue are Y chromosome–linked.
Fig. 2.
Fig. 2.
Cervical spinal cord microglial transcriptomes vary in females compared with males. (A) Principle component analysis of male and female cervical spinal cord microglia (SC; n = 5/treatment). (B) Volcano plot of differentially expressed genes in female vs. male cervical spinal cord microglia. Red dots represent significantly (FDR < 0.05) upregulated genes, whereas blue dots represent significantly downregulated genes. FC, fold change. (C) Functional enrichment analyses of differentially expressed genes from Biologic Process (GO; left), KEGG pathway (middle), and Reactome pathway (right). Shown are the top five most significant functional enrichments in upregulated genes (red) and in downregulated genes (blue). Enrichment categories are displayed on the y-axis and graphed by FDR, with lower FDR values indicating greater significance. Parentheses show the number of genes differentially expressed compared with the number of total background genes expressed in that category. Biologic Process (GO) categories: translation, peptide metabolic process (Pep Met Proc), cellular amide metabolic process (Cell Am Met Proc), organonitrogen compound biosynthetic process (Org Comp Bio Proc), organonitrogen compound metabolic process (Org Comp Met Proc). KEGG pathway categories: mammalian target of rapamycin signaling pathway (mTOR), breast cancer, ribosome, proteasome, oxidative phosphorylation (Ox Phos), Huntington’s disease (Huntington’s), Alzheimer’s disease (Alzheimer’s). Reactome pathway categories: nonsense mediated decay enhanced by the exon junction complex (NMD En EJC); antigen processing–cross presentation (Antigen Proc – Cross Pres); metabolism of RNA (RNA Met); formation of the ternary complex and, subsequently, the 43S complex (Ternary Complex). (D) STRING network analyses of the 395 upregulated genes. (E) STRING network analyses of the 511 downregulated genes. (F) MAGICTRICKS summary output for upregulated genes in female cervical spinal microglia (top 30 factors). (G) MAGICTRICKS summary output for downregulated genes in female cervical spinal microglia (top 30 factors). TF, transcription factor.
Fig. 3.
Fig. 3.
Transcriptomic divergence between male cervical spinal and brainstem microglia. (A) Principle component analysis of male cervical spinal cord (SC) and brainstem (BS) microglia (n = 5/treatment). (B) Volcano plot of differentially expressed genes in male cervical spinal cord microglia. Red dots represent significantly (FDR < 0.05) upregulated genes, blue dots represent significantly downregulated genes. FC, fold change. (C) Functional enrichment analyses of differentially expressed genes from Biologic Process (GO; left), KEGG pathway (middle), and Reactome pathway (right). Shown are the top five most significant functional enrichments in upregulated genes (red) and in downregulated genes (blue). Enrichment categories are displayed on the y-axis and graphed by FDR, with lower FDR values indicating greater significance. Parentheses show the number of genes differentially expressed in that category compared with the number of total background genes expressed in the category. Biologic Process (GO) categories: positive regulation of nucleic acid–templated transcription (Pos Reg NA Trans); positive regulation of RNA metabolic process (Pos Reg RNA Met Proc); positive regulation of transcription by RNA polymerase II (Pos Reg RNA Poly II Trans); positive regulation of transcription, DNA templated (Pos Reg DNA Trans); positive regulation of nucleobase-containing metabolic compound (Pos Reg Nuc Met Proc); cell communication (Cell Comm); ion transmembrane transport (Ion Transm Transport); signaling; response to stimulus (Resp to Stim); signal transduction (Sig Transd). KEGG pathway categories: tumor necrosis factor signaling pathway (TNF), interleukin-17 signaling pathway (IL-17), osteoclast differentiation (Osteo Diff), NF-κB signaling pathway, hepatitis B, focal adhesion, extracellular matrix–receptor interaction (ECM Int), regulation of actin cytoskeleton (Reg Act Cyto), arrhythmogenic right ventricular cardiomyopathy (ARVC), hypertrophic cardiomyopathy (HCM). Reactome pathway categories: extracellular matrix organization (Extra Mat Org), signaling by receptor tyrosine kinases (Sig Rec Tyr Kin), transport of small molecules (Trans Sm Molec). (D) STRING network analyses of the 826 upregulated genes. (E) STRING network analyses of the 636 downregulated genes. (F) MAGICTRICKS summary output for upregulated genes in male cervical spinal microglia (top 30 factors). (G) MAGICTRICKS summary output for downregulated genes in male cervical spinal microglia (all 22 factors). TF, transcription factor.
Fig. 4.
Fig. 4.
Transcriptomic divergence between female cervical spinal and brainstem microglia. (A) Principle component analysis of female cervical spinal cord (CSC) and brainstem (BS) microglia (n = 5/treatment). (B) Volcano plot of differential gene expression changes in female C3–C6 cervical spinal cord vs. brainstem microglia. Red dots represent significantly (FDR < 0.05) upregulated genes, and blue dots represent significantly downregulated genes. FC, fold change. (C) Functional enrichment analyses of differentially expressed genes from Biologic Process (GO; left), KEGG pathway (middle), and Reactome pathway (right). Shown are the top five most significant functional enrichments in upregulated genes (red) and in downregulated genes (blue). Enrichment categories are displayed on the y-axis and graphed by FDR, with lower FDR values indicating greater significance. Parentheses show the number of genes differentially expressed in that category compared with the number of total background genes expressed in the category. Biologic Process (GO) categories: positive regulation of transcription by RNA polymerase II (Pos Reg RNA Poly II Trans), regulation of RNA metabolic process (Reg RNA Met Proc), regulation of gene expression (Reg Gene Exp), regulation of macromolecule metabolic process (Reg Macro Met Proc), positive regulation of macromolecule metabolic process (Pos Reg Macro Met Proc), organonitrogen compound biosynthetic process (Org Comp Bio Proc), translation, cellular process, cation transport, peptide metabolic process (Pep Met Proc). KEGG pathway categories: tumor necrosis factor signaling pathway (TNF), Epstein-Barr virus signaling pathway (EBV), interleukin-17 signaling pathway (IL-17), osteoclast differentiation (Osteo Diff), spliceosome, oxidative phosphorylation (Ox Phos), Parkinson’s disease (Parkinson’s), ribosome, Alzheimer’s disease (Alzheimer’s), Huntington’s disease (Huntington’s). Reactome pathway categories: mRNA 3′-end processing (mRNA 3′ Proc); chromatin modifying enzymes (Chrom Mod Enz); cleavage of the growing transcript in the termination region (Transcript cleavage); respiratory electron transport, ATP synthesis by chemiosmotic coupling, and heat production by uncoupling proteins (Elect. Transp. Chain); the citric acid cycle and respiratory electron transport (TCA and Resp ET); respiratory electron transport (Resp ET); complex I biogenesis (Comp Bio); nonsense mediated decay enhanced by the exon junction complex (NMD Enh EJC). (D) STRING network analyses of the 1167 upregulated genes. (E) STRING network analyses of the 786 downregulated genes. (F) MAGICTRICKS summary output for cervical spinal cord upregulated genes (top 30 factors). (G) EZH2 is the sole factor identified from MAGICTRICKS summary output for downregulated genes in cervical spinal microglia. TF, transcription factor.
Fig. 5.
Fig. 5.
Functional gene categories have little overlap between males and females. (A) Principle component analyses of microglia from male and female brainstem (BS) and cervical spinal cord (SC) (n = 5/treatment) show that CNS region, not sex, is the stronger driver of differential gene expression. (B) Venn diagram showing the total number of differentially upregulated genes in male and female spinal microglia compared with brainstem microglia. Of the differentially upregulated genes in spinal compared with brainstem microglia, 304 were shared by both sexes, 522 genes were unique to males, and 863 were unique to females. Genes expressed in males are shown in green, and genes expressed in females are shown in purple. (C) Venn diagram showing the total number of differentially upregulated genes in male and female spinal microglia compared with brainstem microglia. Of the differentially upregulated genes in spinal compared with brainstem microglia, 514 were shared by both sexes, 122 genes were unique to males, and 272 were unique to females. Genes expressed in males are shown in green, and genes expressed in females are shown in purple. Tables Biologic Process (GO class) (D), KEGG pathway (E), and Reactome pathway (F) summarize the top five functional gene categories into which the differentially expressed genes in spinal microglia fall, among both both sexes. Categories that are comparatively upregulated (↑ red shading) and downregulated (↓ blue shading) in male and female spinal microglia (relative to the expression of those genes in male and female brainstem microglia) are shown separately and demonstrate that very few functional gene categories overlap between the sexes. ECM, extracellular matrix; IL-17, interleukin-17; NMD, nonsense mediated decay; TCA, tricarboxylic acid; TNF, tumor necrosis factor.

Similar articles

Cited by

References

    1. Agosto-Marlin IM, Nichols NL, Mitchell GS. (2018) Systemic inflammation inhibits serotonin receptor 2-induced phrenic motor facilitation upstream from BDNF/TrkB signaling. J Neurophysiol 119:2176–2185. - PMC - PubMed
    1. Ahuja D, Mateika JH, Diamond MP, Badr MS. (2007) Ventilatory sensitivity to carbon dioxide before and after episodic hypoxia in women treated with testosterone. J Appl Physiol (1985) 102:1832–1838. - PubMed
    1. Armoskus C, Moreira D, Bollinger K, Jimenez O, Taniguchi S, Tsai H-W. (2014) Identification of sexually dimorphic genes in the neonatal mouse cortex and hippocampus. Brain Res 1562:23–38. - PMC - PubMed
    1. Baldy C, Fournier S, Boisjoly-Villeneuve S, Tremblay ME, Kinkead R. (2018) The influence of sex and neonatal stress on medullary microglia in rat pups. Exp Physiol 103:1192–1199. - PubMed
    1. Becklake MR, Kauffmann F. (1999) Gender differences in airway behaviour over the human life span. Thorax 54:1119–1138. - PMC - PubMed

Publication types