Skip to main content
Download PDF

Fetal hypoxia results in sex- and cell type-specific alterations in neonatal transcription in rat oligodendrocyte precursor cells, microglia, neurons, and oligodendrocytes

Cell & Bioscience volume 13, Article number: 58 (2023) Cite this article

Abstract

Background

Fetal hypoxia causes vital, systemic, developmental malformations in the fetus, particularly in the brain, and increases the risk of diseases in later life. We previously demonstrated that fetal hypoxia exposure increases the susceptibility of the neonatal brain to hypoxic-ischemic insult. Herein, we investigate the effect of fetal hypoxia on programming of cell-specific transcriptomes in the brain of neonatal rats.

Results

We obtained RNA sequencing (RNA-seq) data from neurons, microglia, oligodendrocytes, A2B5+ oligodendrocyte precursor cells, and astrocytes from male and female neonatal rats subjected either to fetal hypoxia or control conditions. Substantial transcriptomic responses to fetal hypoxia occurred in neurons, microglia, oligodendrocytes, and A2B5+ cells. Not only were the transcriptomic responses unique to each cell type, but they also occurred with a great deal of sexual dimorphism. We validated differential expression of several genes related to inflammation and cell death by Real-time Quantitative Polymerase Chain Reaction (qRT-PCR). Pathway and transcription factor motif analyses suggested that the NF-kappa B (NFκB) signaling pathway was enriched in the neonatal male brain due to fetal hypoxia, and we verified this result by transcription factor assay of NFκB-p65 in whole brain.

Conclusions

Our study reveals a significant impact of fetal hypoxia on the transcriptomes of neonatal brains in a cell-specific and sex-dependent manner, and provides mechanistic insights that may help explain the development of hypoxic-ischemic sensitive phenotypes in the neonatal brain.

Background

Large clinical and preclinical studies have revealed that an adverse fetal environment increases the risk of cardiovascular, metabolic, and neurobehavioral diseases; these observations have led to what is known as the “fetal origins of adult disease” hypothesis [1,2,3,4]. Fetal hypoxia is a common form of prenatal stress that can occur under a variety of conditions, including pregnancy at high altitude, maternal smoking, maternal anemia, placental insufficiency, and cord compression. Fetal hypoxia can cause developmental malformations, particularly in the brain [5,6,7,8]. This can lead to multiple neurodevelopmental disorders in offspring, such as cerebral palsy, autism, and cognitive dysfunction [8, 9]. Further, our previous studies in rodents demonstrate that fetal hypoxia promotes the development of neurological disease-sensitive phenotypes in offspring, including neonatal hypoxic-ischemic encephalopathy, ischemic stroke, and age-related cognitive decline [10,11,12,13].

Due to the protracted nature of brain development, the time window of fetal hypoxia exposure is a critical factor that affects the ultimate consequences of fetal hypoxia. The processes at each developmental stage affect a specific set of developmental processes in the brain, and interfering with each of these can have a specific set of consequences on specific cell types. Neurons (NRs) are the dominant cell type and basic functional unit in the brain; a series of events related to NR growth proceed during the perinatal period, such as neurogenesis and synaptogenesis, to form the mature nervous system [14, 15]. Glial cells, including microglia (MGs), astrocytes (ACs), and oligodendrocytes (ODCs), comprise at least 50% of all cells in the brain [16, 17]. The prenatal period is a critical stage for glial proliferation, maturation and plasticity. For example, the second wave of ODC progenitor cell specification and migration to the central nervous system (CNS) occurs at this time [18,19,20]. The migration of MGs into the brain happens around the mid-gestation stage (E9) in mouse, and MGs start to achieve ramified morphology and distribute in the cerebral cortex prenatally [21, 22]. Astrogenesis initiates at the end of gestation (E18) and lasts until neonatal day 7 in mice [23]. The interaction of glial cells with NRs is also indispensable for neurodevelopment and later physiological maintenance. Mounting evidence suggests that prenatal hypoxia affects the number of NRs, synaptic density, MG morphology and function, AC differentiation and migration, and ODC maturation [24,25,26].

Although the deleterious effect of fetal hypoxia on neurodevelopment and progression of neurological disorders in later life have been phenotypically characterized, the underlying cellular and molecular mechanisms remain unclear. Several mechanisms have been proposed, however, including increased apoptotic cell death and loss of stem or progenitor cells due to energy deprivation [27], increased reactive oxygen species (ROS) and reactive nitrogen species (RNS), and maturation stage-specific vulnerability of ODC progenitors in the periventricular leukomalacia region [28, 29]. Hypoxia-induced fetal brain injury is also associated with an inflammatory response via overexpression of pro-inflammatory cytokines and pro-apoptotic proteins [9, 30].

Male infants are known to have a larger inflammatory response to ischemic insults than females [31], and fetal hypoxia has been observed to have a sexually dimorphic effect on the mitochondria of fetal guinea pigs [32]. The precise mechanisms that give rise to sexually dimorphic outcomes, and their precise relationship to more general, sex-independent outcomes, are unknown. Understanding the nature of sexually dimorphic responses to fetal hypoxia therefore has the potential to help understand the precise mechanisms of action of fetal hypoxia.

In order to gain insights into the specific brain cell types affected by fetal hypoxia, as well as the degree of sexual dimorphism, we performed next-generation RNA sequencing (RNA-seq) of multiple brain cell types in neonatal rats subjected to either fetal hypoxia or normoxia (controls): NRs, MGs, ACs, ODCs, and ODC precursor cells (A2B5+). We observed that fetal hypoxia induced substantial transcriptomic changes in the neonatal brain in a cell-type specific and sex-dependent manner in NRs, MGs, ODCs, and A2B5+ cells. With pathway and transcription factor motif analyses, we found that immune- and inflammation-related genes were highly enriched in each cell type. A functional assay revealed that the NF-kappa B (NFκB) signaling pathway is specifically activated in the male brain. Thus, this study demonstrates that fetal hypoxia has unique effects on neonatal transcription across multiple brain cell types, and that many of these alterations are sexually dimorphic. The diversity of cell- and sex-specific responses to fetal hypoxia may underly the diversity of disease outcomes that can occur after fetal hypoxia exposure.

Results

Sex- and cell type-dependent transcriptional changes in the neonatal rat brain due to fetal hypoxia

First, we measured the brain and body masses of neonatal rats subjected to fetal hypoxia versus controls. Both brain and body masses were decreased by fetal hypoxia, but the brain to body mass ratio was increased (Additional file 1: Fig. S1A), confirming a typical asymmetric growth restriction resulting from fetal hypoxia [12, 13].

We therefore performed RNA-seq on A2B5+ cells, ODCs, ACs, NRs, and MGs in both male and female rat pups subject to fetal hypoxia or control conditions. Four biological replicates from each treatment group were sequenced, for a total of 80 RNA-seq samples. Prior to sequencing, all RNA samples were subject to and passed quality control measures (Additional file 1: Fig. S1B). We sequenced to an average depth of 24.3 million reads per sample, across all samples. All samples showed uniformly high sequencing quality across both read pairs and passed the FastQC [33] quality control measures (Additional file 1: Fig. S1C). Fragments per kilobase per million (FPKM) normalized values of gene expression were highly correlated across replicates within the same treatment group (Additional file 1: Fig. S1D).

We then determined the set of differentially expressed genes (DEGs) of male vs. female, as well as fetal hypoxia treatment versus controls. 2,195 fetal hypoxia-induced DEGs (hDEGs) were identified, with the largest number of hDEGs occurring in ODCs, NRs, and MGs (Fig. 1A). In terms of the number of hDEGs, the cell type that was most affected was male ODC, which had 1,114 hDEGs (Additional file 2: Table S1). Interestingly, fetal hypoxia affected the expression of more genes than were altered by sex (Additional file 2: Table S2). A small number of genes showed interaction effects, whereby changes in both sex and hypoxia treatment together affected gene expression more than the sum of the individual effects (Additional file 2: Table S3). Strikingly, the number of hDEGs was highly sex-dependent in all brain cell types except ACs (Figs. 1B–G). While some of these sex-specific hDEGs were sexually dimorphic in normoxic control animals, many were not (Fig. 1A). All cell types contained sex-dependant hDEGs that weren’t sexually dimorphic in normoxic controls, with the largest number of these occurring in ODCs. As an example, Nfe2l2 expression was affected by hypoxia to a much greater extent in ODC males than females, and had similar expression in male and female controls (Additional file 1: Fig. S1E).

Fig. 1
figure 1

Unique and sexually dimorphic transcriptomic changes in the neonatal brain due to fetal hypoxia. A. Heatmap of log2 Fold Change (FC) values of hDEGs from the indicated samples, along with normalized gene expression (log2 (FPKM + 1)) of male and female normoxic controls. Genes are ordered by hierarchical clustering of the log2 FC values. In the log2 FC columns (red-blue), M indicates hypoxia vs. control in males, F indicates hypoxia vs. control in females, and N indicates male vs. female in normoxic controls. In the gene expression columns (white-black), M indicates the pooled gene expression of normoxic males, and F indicates the pooled gene expression of normoxic females. B. Bar plots of the number of hDEGs in samples of the indicated cell type and sex. Bars on the positive y-axis indicate the number of DEGs upregulated due to fetal hypoxia, while bars on the negative y-axis indicate the number of DEGs downregulated due to fetal hypoxia. C–G. Venn diagrams showing the number of hDEGs distinct and overlapping between male (M) and female (F). P-values displayed are for the null hypothesis that the number of male and female hDEGs are the same, by Fisher’s Exact Test. P-value cutoffs are as follows: * p < .01; ** p < .001; *** p < 10–5; **** p < 10–10

Full size image

Considering all hDEGs equally as above does not provide any information about which genes have the largest changes in absolute expression and therefore could have a larger influence on certain phenotypes. We therefore performed Principal Component Analysis (PCA) on the RNA-seq data in a way that weights genes with high expression variability more heavily than those with low expression variability (see Methods). When including all samples in a single PCA plot, glial cells (MG, AC, and ODC) clustered closely together and were far from NRs (Additional file 1: Fig. S1F), consistent with the fact that glial cells are more similar to each other than to NRs in terms of neonatal phenotypes and function. The effect of sex and hypoxia on the PCA were both small relative to differences between cell types, again consistent with expectation, since large transcriptome variation across cell types is well known, whereas the number of hDEGs is small relative to the transcriptome sizes of the cells. Our PCA is therefore capable of recapitulating large-scale phenotypic differences between samples.

Given that the results of our global PCA analysis were consistent with expected overall phenotypic variation, which is primarily driven by differences in cell type transcriptomes, we performed a PCA on the samples of each cell type separately in order to more clearly assess the effects of sex and fetal hypoxia on highly changed transcripts. Interestingly, and consistent with the relative number of hDEGs in each sample (Fig. 1), the relative strength and degree of sexual dimorphism of the hypoxia effect was unique to each cell type. Based on PCA, hypoxia was the major effector in A2B5+ cells (Fig. 2A and Additional file 1: Fig. S2A). Although there was a significant difference in the number of hDEGs between A2B5+ males and females (Fig. 1C), the PCA of A2B5+ cells suggests a much smaller degree of sexual dimorphism in genes with the highest variability in A2B5+ cells. On the other hand, hypoxia only had a weak effect in ACs when compared to sex differences in the PCA (Fig. 2B and Additional file 1: Fig. S2B), consistent with the relatively low number of AC hDEGs (Fig. 1B). MGs were substantially affected by hypoxia, but the effect on the PCA was stronger on average in males than in females (Fig. 2C and Additional file 1: Fig. S2C), consistent with a greater number of MG hDEGs in males than females (Fig. 1E). A weak effect of hypoxia was evident in the PCA for NRs (Fig. 2D and Additional file 1: Fig. S2D) despite a large number of hDEGs in NRs (Fig. 1F), suggesting that in NRs, hypoxia induces relatively small changes in absolute gene expression as compared to sex. In the PCA for ODCs, hypoxia had a larger effect on females on average than males (Fig. 2E and Additional file 1: Fig. S2E) despite the fact that there were more hDEGs in males than females (Fig. 1G). Thus, highly changed transcripts are affected more in female ODCs than males, despite there being many more male than female hDEGs. Interestingly, the sex difference in normoxic MGs was largely ablated by fetal hypoxia, suggesting that, for transcripts with the largest changes in MG, the sexual dimorphism is driven by the convergence of sexually distinct normal phenotypes to a common hypoxia phenotype. On the other hand, the PCA of ODCs suggests that under normal conditions, the most highly variable ODC genes are expressed similarly in males and females, but these genes respond to fetal hypoxia uniquely in each sex.

Fig. 2
figure 2

Principal component analysis of RNA-seq data. The first two principal components in cell type specific PCAs for A. A2B5 + cells; B. astrocytes; C. microglia; D. neurons; E. oligodendrocytes

Full size image

Taken together, these results demonstrate that fetal hypoxia induces unique and sex-dependent transcriptomic alterations in neonatal A2B5+ cells, MGs, NRs, and ODCs, and suggest that these transcriptomic alterations may result in significant and sexually dimorphic alterations in the phenotypes of each of these cell types.

Pathway analysis of the hDEGs

To get insights into the biological consequences of fetal hypoxia on the neonatal brain, we performed pathway analyses on the hDEGs independently for each cell type and sex. Each cell type and sex had a unique set of enriched pathways (Fig. 3). Phagocytosis and cell death were primarily enriched in NRs and ODCs, while all cell types had enrichment of at least one pathway related to immunity and inflammation. For example, the NF-kappa B signaling pathway was significantly enriched in male A2B5+ cells and MGs, and both male and female NRs and ODCs. Transcriptional changes in both the canonical and non-canonical NF-kappa B pathways occurred in male MGs (Additional file 1: Fig. S3A). Diagrams of additional pathways are included for reference (Additional file 1: Fig. S3B–K).

Fig. 3
figure 3

Enriched pathways from fetal hypoxia induced DEGs. Plot indicating the significantly enriched pathways (P < .05) in each of the cell types and sexes indicated

Full size image

Transcription factor binding motifs associated with hDEGs in brain cells

An important driver of environmentally-induced changes in gene expression is changes in transcription factor (TF) binding at enhancers. We therefore performed a motif enrichment analysis at enhancers near hDEGs for all brain cell types that had enhancer annotations available in the rat (Fig. 4A–D). Only TFs that are expressed according to the RNA-seq data for the given cell type and sex were shown. Male MGs manifested the highest number of hypoxia-enriched, putative TF binding sites (TFBSs). Examples illustrating the expression of some of these TFs, and changes in gene expression of a few example genes with nearby TF motifs are shown in Additional file 1: Fig. S4A–L. Interestingly, the NF Kappa B1 subunit (NFκB1) motif was significantly enriched in male MGs and ODCs, consistent with the enrichment of the NF Kappa B pathway amongst hDEGs in both cell types (Fig. 3).

Fig. 4
figure 4

Enriched TF binding motifs near fetal hypoxia-induced, neonatal DEGs. Bar plots indicating the number of putative TF binding sites (TFBSs) at enhancers, based on sequence motif hits near fetal hypoxia DEGs in A. A2B5 + cells; B. microglia; C. neurons; D. oligodendrocytes. Indicated P-values are by Fisher’s Exact Test

Full size image

Validation of fetal hypoxia-induced DE of genes related to inflammation and cell death by qRT-PCR

The results above based on RNA-seq demonstrate putative alterations in a number of genes involved in both inflammation and cell death. To validate this, we performed qRT-PCR to confirm select DEGs in male and female MGs, NRs, ODCs, and A2B5+ cells. The following hDEGs were validated in the indicated cell types: Il10, Tnfsf13b, Nfkbib, and Ikbke in male MGs; Abcg2 and Hspa1a in female MGs; Alox5, Steap3, Tnf, and Lcn2 in male NRs; Alox5 and Tnf in female NRs; Tnf, Lox, and Lcn2 in male ODCs; Il1b and Alox15 in female ODCs; Il1a, Il1b, and Lox in male A2B5+ cells; Abcg2 and Hspa1a in female A2B5+ cells (Fig. 5A). The qRT-PCR results were strongly correlated to the corresponding results in RNA-seq (Fig. 5B), indicating a strong concordance between the DEG results of RNA-seq and of qRT-PCR analysis.

Fig. 5
figure 5

Validation of inflammation and cell death related fetal hypoxia-induced, neonatal DEGs. A. Quantification of mRNA abundance of genes by qRT-PCR in the indicated cell types sorted from neonatal brains with fetal hypoxia or normoxia exposure. NR, neurons; MG, microglia; ODC, oligodendrocytes; A2B5, A2B5 + cells; M, male; F, female. Error bars indicate mean ± SEM. P-values are by t-test, with n = 5 pups per group. *P < .05. B. Comparison of hypoxia-induced log2 Fold Changes according to RNA-seq vs. qRT-PCR of hDEGs validated by qRT-PCR

Full size image

Fetal hypoxia enhanced the activity of NFκB-p65 in the brain

To validate the activation of the NFκB pathway due to fetal hypoxia in the neonatal rat brain, we conducted a protein assay of nuclear translocation of a key NFκB subunit, p65, in the neonatal brain. Importantly, we found that NFκB-p65 activity was significantly increased by fetal hypoxia in male, but not female, neonatal brains (Fig. 6).

Fig. 6
figure 6

Effect of fetal hypoxia on NFκB-p65 activity in the brains of neonatal rats. Plots showing the average (tops of bars), standard error (error bars), and distribution (individual points) of the fetal hypoxia-induced increase of NFκB-p65 nuclear translocation in whole brains of the indicated samples. P-values are determined by two-way ANOVA. *P < 0.05, hypoxia versus control

Full size image

Discussion

In this study, we used RNA-seq of sorted brain cells of rat pups to analyze neonatal transcriptomic changes associated with fetal hypoxia. To the best of our knowledge, this is the first study to measure the transcriptional profile of cell populations in the developing brain of rats using RNA-seq. The results indicate that fetal hypoxia induces substantial transcriptomic alterations in neonatal NRs and glial cells, in a cell type- and sex-specific manner.

Our previous studies demonstrate that fetal hypoxia exposure increases the susceptibility of the neonatal brain to hypoxic-ischemic insult [12, 34], which mainly occurs during the perinatal period in humans. Brain development of humans at birth is roughly equivalent to rodents at postnatal days 7–10 [35,36,37,38,39]. The main goal of this study was to explore the transcriptomic response of brain cell populations with hypoxic-ischemic sensitive phenotypes, so we sampled rat brains at postnatal day 7. Due to the heterogeneity of brain development and the ability of self-repair or -renewal, the transcriptomic response of each brain cell type may be different at different life stages following hypoxia exposure in the fetal stage. For example, our previous study determined that gestational hypoxia significantly increases the expression of genes related to astrogliosis and microgliosis, but doesn`t affect the expression of synaptic proteins, in the brain of 2-month-old mouse offspring [13].

Our pathway analyses suggest that immune and inflammatory pathways may be highly affected by fetal hypoxia in neonatal MGs, NRs, and ODCs. This could contribute to the activation of cell death pathways in NRs and ODCs, thus leading to increased susceptibility of the brain to later insults. In addition, we demonstrated an overall reduction in neonatal brain mass due to fetal hypoxia. An increase in inflammation and cell death could be the cause of this reduction in brain mass, although further studies are needed to directly link the two. We performed qRT-PCR in MGs, NRs, ODCs, and A2B5+ cells, and directly demonstrated that fetal hypoxia induces neonatal transcriptomic changes in genes with known involvement in inflammatory and cell death pathways.

Among the pathways and TF motifs we analyzed bioinformatically, NFκB signaling was commonly enriched in multiple brain cell types, and this culminated in verification of an increased nuclear translocation and activity of NFκB-p65 in hypoxia-exposed male brains. Sexually dimorphic responses to fetal hypoxia have been observed in the past, but the precise mechanisms are unknown. The finding that NFκB activity is increased in the neonatal brain in a sex-dependent manner may provide insight into the mechanisms involved in sexually dimorphic responses to fetal hypoxia. The NFκB signaling pathway plays a key role in inflammation by regulating immune response-related gene expression.

Our recent study showed that NFκB-p65 activation increases the production of an important inflammatory response gene, IL-1β, after neonatal hypoxic-ischemic brain injury in mice [40]. Further, it has been reported that fetal stress causes priming of MG activation [41,42,43], which is regulated by the NFκB signaling pathway [44]. We validated immune response related hDEGs in male MGs using qRT-PCR and found that Il10, Ikbke, and Nfkbib were downregulated and Tnfsf13b was upregulated. Il10 exerts anti-inflammatory actions in glial cells and can inhibit activation of the NFκB pathway [45, 46]. Nfkbib and Ikbke encode proteins of NFκB endogenous inhibitors [47]. Tnfsf13b is a member of the TNF superfamily, which is a sensome gene of MGs [48], and is active, along with its receptor, in inflammation [49]. Thus the results of our study suggest that fetal hypoxia may prime male MGs through the activation of the NFκB signaling pathway.

Interestingly, immune and inflammation pathways were also activated by fetal hypoxia in ODCs and A2B5+ cells. For example, the key cytokines Il1b and Il1a were upregulated in the ODC lineage, which can result in excessive inflammation and ODC and A2B5+ cell death [50]. IL-1β is consistently upregulated by perinatal encephalopathy in both humans and experimental models [51,52,53], and inhibition of IL-1β reduces ODC loss following prenatal inflammation in near-term fetal sheep [54]. These findings suggest that increased inflammation after hypoxia exposure is associated with a greater risk of white matter injury after birth. As a pleiotropic cytokine, IL-1β can induce multiple effects in the injured brain. In additional to its detrimental effect, IL-1β also participates in CNS repair in the adult brain after brain injury [55]. Moreover, IL-1β regulates ODC development by blocking proliferation at the late progenitor/pro-ODC (O4+) stage, but does not affect proliferation of early progenitors (A2B5+) [56], suggesting that the upregulation of IL-1β and other immune-related pathways in A2B5+ cells may be a tissue repair mechanism to rescue fetal stress-induced ODC death. Our pathway analysis also showed that cell death pathways were enriched in NRs and ODCs, but not A2B5+ cells. Therefore, the potential double-faced actions of pro-inflammatory cytokines in ODC and A2B5+ cells are worth further investigation.

Phagocytosis and cell death pathways were enriched in NRs and ODCs after fetal hypoxia, which may also be ascribed to immune and inflammation responses. Phagocytic NRs have been observed, with proposed involvements in debris clearance, neuronal death, and cell-to-cell spreading of disease [57]. It is still unclear which types of hypoxia-induced cell death occur in early neonatal brains. For the first time, we identified necroptosis and ferroptosis as cell death pathways that may be upregulated in neonatal NRs and ODCs exposed to fetal hypoxia. Necroptosis is unregulated necrotic death, stimulated by secretion of cytokines/chemokines, resulting in inflammation, which is initiated by TNF superfamily receptors, toll-like receptors, and interferon receptors [58, 59]. All these pathways were enriched in neonatal NRs and ODCs by RNA-seq pathway enrichment analysis.

Ferroptosis is another cell death pathway enriched in NRs and ODCs, and is characterized by a large amount of iron accumulation and lipid ROS during the cell death process [60]. In the present study, necroptosis and ferroptosis pathway alterations were supported by qRT-PCR of Tnf, Steap3, Alox5, and Lcn2 in NRs, and Il1β, Alox15, and Pygl in ODCs. Steap3 is a ferroptosis-related gene encoding an iron transporter, and dysregulation of which leads to impaired iron homeostasis and neurodegeneration [61]. Neuronal derived lipocalin-2 (Lcn2), another hDEG in NR, is a signature of several CNS anomalies and participates in inflammation, NR death and neurodegeneration [62,63,64,65]. Alox5 and Alox15 were also hDEGs in NR; these genes regulate lipid peroxidation and cause ferroptosis and inflammation [66, 67]. Moreover, Alox15 expression is increased in the brains of periventricular leukomalacia (PVL) patients [68]. The upregulation of cell death-related genes in ODCs agrees with previous findings that fetal hypoxia induces ODC loss and white matter injury. On the other hand, cell death-related genes in NRs, such as Tnf, Lcn2 and Alox5, were downregulated, which has been reported to be associated with anti-inflammation and neuroprotection [65, 69, 70]. This result may suggest a recovery process in neuronal cells in the early neonatal stage from prenatal hypoxia exposure.

The proper development and function of the CNS is particularly sensitive to oxygen concentration, and hypoxia is associated with developmental malfunctions in offspring. The impact of fetal hypoxia on the neonatal brain at the cellular level is tightly regulated by transcriptomic changes at the molecular level. Our data suggest that ODCs have the greatest transcriptomic response to fetal hypoxia, followed by NRs, MGs, and A2B5+ cells. These data indicate that the oligodendroglial lineage is highly vulnerable to fetal hypoxia in the developing brain. Our findings are in line with the notion that PVL is mainly caused by the dysregulation of the oligodendroglial lineage in perinatal white matter [71,72,73]. Along with white matter injury, cortical and subcortical grey matter is also significantly affected by fetal hypoxia, resulting in neuronal damage in the premature brain [27, 74]. Indeed, fetal hypoxia can result in both subcortical white matter and gray matter injury in preterm sheep [75,76,77]. PVL white matter lesions are usually surrounded by astrogliosis and microgliosis [78]. The perinatal period is a critical window for microglia distribution, ramification, and maturation in the brain. Thus, MGs during this period are particularly sensitive to environmental influences, which may result in long-term changes in MG cell number and/or phenotypes [13, 79]. We have previously demonstrated that fetal hypoxia increases microgliosis and MG activation in adult mouse offspring [13].

AC was the cell type with the smallest number of hDEGs. ACs originate from neural cells at the end of the neurogenic process (E18) in rodents. Thus, a possible explanation for this is that fetal hypoxia may affect neural differentiation into ACs, leading to a reduced ratio of NRs to ACs in the immature brain [25, 80], instead of a direct effect on the AC transcriptome. Indeed, our recent study showed that prenatal hypoxia exposure significantly increases astrogliosis in the brain of adult offspring [13]. More research is needed to understand the mechanistic alterations of neural differentiation in response to fetal hypoxia.

Conclusions

In the present study, we performed RNA-seq on sorted brain cells to analyze transcriptomic changes in rat pups after fetal hypoxia exposure. We generated global gene expression profiles for five major brain cell populations: NRs, MGs, ODCs, ACs and A2B5+ cells. Our data indicate that immune- and inflammation-related pathways are activated in neonatal NRs, MGs, ODCs, and A2B5+ cells. We validated a sexually dimorphic, male-specific activation of the NFκB signaling pathway, and RNA-seq signal of numerous genes is affected in a sexually dimorphic manner within NRs, MG, ODCs, and A2B5+ cells. NFκB and other inflammatory and immune related pathways identified in our study may be major mediators of fetal hypoxia-induced immature brain malformation and damage, which is associated with increased risk for neurological disorders in later life. The genes and pathways identified in this study are therefore potential therapeutic targets for attenuating fetal hypoxia-induced immature brain malformation and damage.

Methods

Experimental animals

Pregnant Sprague–Dawley (SD) rats (Charles River Laboratories, Portage, MI, USA) were randomly divided into two groups: (i) the normoxic control and (ii) hypoxic exposure (10.5% O2) from Day 15 to Day 21 of gestation, as we previously described [11, 34, 81,82,83]. The pups were sacrificed at neonatal day 7, and the brains were collected for cell sorting. All experimental procedures and protocols were approved by the Institutional Animal Care and Use Committee of Loma Linda University and followed the guidelines in the National Institutes of Health Guide for the Care and Use of Laboratory Animals.

Cell sorting

Five types of brain cells (microglia, neurons, astrocytes, oligodendrocytes, A2B5+ glia precursor cells) were sorted from P7 pups using magnetic-activated cell sorting (MACS) technology. Briefly, brains were harvested from neonatal rats and then dissociated to single-cell suspensions using a Neural Tissue Dissociation Kit-T (Miltenyi Biotec) according to the manufacturer’s instructions. Myelin debris in single-cell suspensions was depleted by using Myelin Removal Beads II (Miltenyi Biotec) according to the manufacturer’s instructions. The unlabeled cell fraction after myelin removal was subjected to subsequent isolation of different neural cell populations, with the following MicroBead kits (Miltenyi Biotec): for microglia, CD11b/c MicroBeads; for astrocytes (precursor/immature cells), Anti-GLAST (ACSA-1) MicroBeads; and for oligodendrocytes (pre/immature/mature cells), Anti-O4 MicroBeads. For neurons (progenitor/immature cell), Anti-PSA-NCAM MicroBeads were used for sorting after depletion of A2B5+ glial progenitor cells using Anti-A2B5 MicroBeads.

RNA isolation

Total RNA was isolated from each cell type using the Qiagen miRNeasy kit (Cat # 217004) following manufacturer’s protocol. In brief, QIAzol Lysis Reagent was used to disrupt and homogenize the cells and the lysate was incubated at room temp for 5 min. Then chloroform was added to the lysate, and shaken vigorously for 15 s. The mixture was re-incubated for 5 min and centrifuged for 15 min at 14000 RPM at 4 °C. The upper aqueous phase was collected to a fresh tube and 1.5 volumes of 100% ethanol was added. Then the sample was processed through a RNeasy Mini spin column. On-column DNase digestion was performed to effectively remove DNA during RNA purification by using the RNase-Free DNase set (Qiagen, cat# 79254) according to the manufacturers’ instructions. After final column wash, total RNA was eluted out of mini column membranes using 30 µl nuclease free water (Qiagen). Integrity of RNA was confirmed by determining RIN number. Total RNA was delivered to the Center for Genomics at Loma Linda University for library construction and RNA sequencing.

RNA-Seq library construction

The Ovation® Rat RNA-Seq System (NuGEN Technologies, San Carlos, CA) was used per manufacturer’s instructions to construct RNA-seq libraries. 100 ng of total RNA was used as input. First and second strands of cDNA were synthesized from total RNA (100 ng) spiked with 1 µl of 1:500 diluted ERCC ExFold RNA Spike-In Mix 2 (Life Technologies, Carlsbad, CA) at a final concentration of 1%. Following primer annealing and cDNA synthesis, end-repair, adaptor index ligation, and strand selection was conducted. Barcodes with unique indices was used per sample for multiplexing. Ribosomal RNA depletion was performed by using custom InDA-C primer mixture SS5 V8 for rat. Finally, libraries were amplified for 13 cycles (Mastercycler® pro, Eppendorf, Hamburg, Germany), and purified with Agencourt XP beads (Beckman Coulter, Indianapolis, IN).

Library quantification, quality control (QC) and sequencing

The final amplified libraries were purified using Agencourt XP beads and quantified using the Qubit dsDNA HS Kit on Qubit 3.0 Fluorometer (Life Technologies, Carlsbad, CA). Quality and peak size were determined with the D1000 ScreenTape on Agilent 2200 TapeStation (Agilent Technologies, Santa Clara, CA). Final libraries were pooled and sequnced on an Illumina NextSeq 550 or HiSeq 4000 (Illumina, Inc., San Diego, CA). Paired-end reads of 75 bp were generated.

RNA-seq data analysis

Raw fastq files were trimmed to remove adaptors and low quality bases using trim_galore [84] version 0.6.6., a wrapper for cutadapt [85] (version 2.8). Trimmed fastq files were then mapped to the rn7 reference genome using tophat (version 2.1.1) [86] using the refFlat gene annotation for rn7, obtained from the UCSC table browser [87]. A custom script was used to select only reads that mapped uniquely to one location on the genome for downstream analysis by taking only reads marked with “NH:i:1” in the tophat2 output.bam files. Normalized gene expression (FPKM) was quantified for each sample using cufflinks version 2.2.1 [88]. Differential gene expression analyses were performed using cuffdiff version 2.2.1 [89] using q < 0.05 and |log2FC|> 1 as the cutoffs for differential expression. To determine genes with interaction effects between sex and hypoxia, a two-way analysis of variance (ANOVA) was performed on FPKM values of genes output by cufflinks, using sex and hypoxia as the factors. The analysis was done on each cell type independently in R version 4.0.3, using the function “aov()”, with the input argument “scale.unit = F”. Interaction P-values from the ANOVA were input into the R function “qvalue()”, and genes with q < 0.05 for the interaction term were used in the final list of interaction genes.

Statistical comparison of the number of male vs. female hDEGs

For each cell type, a 2X2 contingency matrix,

$$\begin{array}{ll} {\text{X}}_{11} & {\text{X}}_{12} \\ {\text{X}}_{21} & {\text{X}}_{22} ,\end{array}$$

was constructed, where X11 = # of female hDEGs, X12 = # of non-hDEGs in female, X21 = # of male hDEGs, X22 = # of non-hDEGs in male. Non-hDEGs are the total number of genes analylzed minus the number of hDEGs in the given sex and cell type. Fisher’s Exact Test was then calculated on each 2X2 contingency matrix in order to obtain a P-value for the null hypothesis that the proportion of genes that were hDEGs in male is the same as the proportion in female.

Principal component analyses

FPKM values output by cufflinks, as described above, were used as inputs for the PCA. The analysis was done in R version 4.0.3. The R function “aov()” was used to perform a linear analysis of variance (ANOVA) at each gene, using each sample as a factor. The R function “qvalue()” was then used to calculate ANOVA q-values for all genes, based on the P-values output by “aov()”. For genes with q < 0.05 in the ANOVA, the FPKM of all replicates were then input into the R function “PCA()”, with the input argument “scale.unit = F”, in order to generate PCA plots.

Pathway analyses

DEG symbols for each cell type were input separately into the Enrichr pathway analysis online tool [90]. Enriched pathways in this paper are from the KEGG 2019 Mouse pathway analysis in the “legacy” tab of the Enrichr website. Pathway shown in this page are rat-specific diagrams (species code rno) from the Kyoto Encyclopedia of Genes and Genomes website [91].

Transcription factor motif enrichment analyses

Annotated enhancers for rat ODCs, ODC precursors, and NRs were obtained from EnhancerAtlas2.0 [92]. EnhancerAtlas2.0 does not contain rat MG enhancers, so instead we used published DNAse-I Hypersensitive Sites (DHSSs) taken from rat MG under a variety of conditions associated with inflammation [93]. Both the DHSSs and enhancers from EnhancerAtlas2.0 were in rn5 coordinates, and so all were converted to rn7 coordinates using UCSC liftOver [94].

For each brain cell type, we compiled a list of expressed genes by taking the set of genes with FPKM > 1 in at least one treatment for the cell type (hypoxia or control), with replicates of each sample pooled before calculating FPKM values in cufflinks as described above. All regions within 500 base pairs from any transcription start site (TSS) of an expressed gene were removed from subsequent analysis, as it would be ambiguous in those cases whether the sites are acting as enhancers or promoters. We also removed enhancers/DHSSs that were > 10 kilobases from the nearest TSS of an expressed gene.

We extracted the DNA sequences at the remaining enhancers/DHSSs, using bedtools getfasta (version 2.27.1. TF motif scanning was then done at these sequences using fimo [95] version 5.3.0, using a combination of the jolma2013.meme, JASPAR_CORE_2014_vertebrates.meme, and uniprobe_mouse.meme TF motif data files provided by fimo for motif scanning. Then, in each cell type, for each TF that was expressed in the cell type, we used Fisher’s exact test to calculate the P-value for the enrichment of the number of enhancers containing a motif for the TF at enhancers near DEGs relative to enhancers/DHSSs near all expressed genes (with proximity defined as within 10 kilobases of the TSS).

qRT-PCR assay

Total RNA was converted to cDNA using Maxima First Strand cDNA Synthesis kit with dsDNase (Thermo Fisher) following manufacture’s protocol. Briefly, 1 µg total RNA in 10 µl containing 1X DNAse buffer and 1 µl DNase was incubated for 2 min at 37 °C, followed by incubation on ice for 5 min. 10 µl of 1X Maxima Enzyme Mix was next added to it and incubated for 10 min at 25 °C followed by 30 min at 50 °C, 85 °C for 5 min and then hold at 4 °C until ready for QPCR. The specific qPCR primers for target genes were designed using Primer Blast tool. https://www.ncbi.nlm.nih.gov/tools/primer-blast/index.cgi?LINK_LOC=BlastHome. All primers were synthesized at the Integrated DNA Technologies. The cDNA template (10 ng), forward and reverse primers each at 0.5 µM, IQ SYBR Green Supermix (BIO-RAD) were mixed in a total of 20 µl volume. β-Actin primers were used as an internal control. Each qPCR was done for 30 cycles. The relative expression of gene of interest was calculated by the 2−ΔΔCT method [96] and is presented as the fold modulation in hypoxia relative to the control. PCR primer details are in the supplemental materials (Additional file 2: Table S4).

Nuclear extraction

The Nuclear Extract Kit (Active Motif, Cat # 40010) was used to isolated nuclear from the brain of rat pups according to manufacturer’s instruction. The brain hemisphere was transferred to a pre-chilled glass Dounce homogenizer with 1 ml of Hypotonic Buffer supplemented with protease inhibitors, phosphatase inhibitors, DTT, and detergent. The tissue was homogenized 30 times on ice bath until dispersed to a homogeneous cell suspension. Then the cell suspension was incubated 15 min on ice, vortexed and re-incubated on ice for 15 min. The cell suspension was first centrifuged at 2500 RPM for 10 min to get rid of tissue debris. The collected supernatant was centrifuged at 14000 RPM for 30 min again. The pellet was resuspended in 100 µl Complete Lysis buffer containing DTT, AM1, protease inhibitor, detergent and vortexed at maximum strength for 20 s followed by incubation on ice. This procedure was repeated 10 times. Finally, the content was centrifuged for 20 min at 14000RPM at 4 °C. The supernatant was the Nuclear Extract (NE), which was stored at −80 °C. The protein concentration of each NE was determined by Bradford Method before TransAM NFκBp65 assay.

TransAM NFκBp65 assay

The TransAM NFκBp65 kit (Active Motif, Cat # 40096) was used to measure the activation of NFκB pathway according to the manufacture’s instruction. Complete Lysis Buffer (CLB) and Complete Binding Buffer (CBB) were prepared using the Kit reagents. CBB (30 µl) was first added to each well followed by NE (20 µg protein in 20 µl CLB) prepared from each tissue. The assay was done in triplicate. For negative controls, only 20 µl CLB was added to 30 µl of CBB. For positive control wells, 2.5 µg of stimulated Jurkat cell nuclear extract (kit component) was used. When all samples were added to the wells, the adhesive cover was added to seal the plate. The plate was incubated at room temperature for 1 h with gentle rocking at around 100 RPM. Then, each well was washed 3 times with 200 µl wash buffer followed by 100 µl of NFκB antibody (1:1000 dilution in antibody binding buffer). The plate was sealed and incubated at room temperature for 1 h with no agitation. After wash, 100 µl of HRP conjugated secondary antibody (1:1000 dilution in antibody binding buffer) was added and incubated at room temperature for another 1 h with no agitation. After wash, 100 µl of Developing solution was added to each well, incubated for 5 min at room temp protected from light, until the blue color develops, and the colorimetric reaction was stopped by adding stop solution. The absorbance was read in a 96 well plate reader at 450 nm with a reference wave-length of 655 nm.

Statistical analysis

For non-bioinformatic data analysis, data were presented as mean ± SEM, and were analyzed using Student’s t-test or two-way ANOVA followed by Tukey’s post hoc test. P-value < 0.05 was considered statistically significant. Experimental number (n) represents brains of animals from different dams.

Availability of data and materials

RNA-seq data that support the findings of this study have been deposited in NCBI GEO with the accession code GSE216946. All custom scripts used to process and analyze the data and to obtain the results described in this paper have been deposited in the GitHub repository (https://github.com/ikremsky/Scripts-for-Fetal-hypoxia-results-in-sex--and-cell-type-specific-alterations-in-neonatal-transcript).

References

  1. de Boo HA, Harding JE. The developmental origins of adult disease (Barker) hypothesis. Aust N Z J Obstet Gynaecol. 2006;46:4–14.

    Article  PubMed  Google Scholar 

  2. Warner MJ, Ozanne SE. Mechanisms involved in the developmental programming of adulthood disease. Biochem J. 2010;427:333–47.

    Article  CAS  PubMed  Google Scholar 

  3. Barker DJ. Developmental origins of adult health and disease. J Epidemiol Community Health. 2004;58:114–5.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  4. Langley-Evans SC, McMullen S. Developmental origins of adult disease. Med Princ Pract. 2010;19:87–98.

    Article  PubMed  Google Scholar 

  5. Katz JM, Pacia SV, Devinsky O. Current management of epilepsy and pregnancy: fetal outcome, congenital malformations, and developmental delay. Epilepsy Behav. 2001;2:119–23.

    Article  PubMed  Google Scholar 

  6. Miller SL, Huppi PS, Mallard C. The consequences of fetal growth restriction on brain structure and neurodevelopmental outcome. J Physiol. 2016;594:807–23.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  7. Phillips TJ, Scott H, Menassa DA, Bignell AL, Sood A, Morton JS, Akagi T, Azuma K, Rogers MF, Gilmore CE, et al. Treating the placenta to prevent adverse effects of gestational hypoxia on fetal brain development. Sci Rep. 2017;7:9079.

    Article  PubMed  PubMed Central  Google Scholar 

  8. Li Y, Gonzalez P, Zhang L. Fetal stress and programming of hypoxic/ischemic-sensitive phenotype in the neonatal brain: mechanisms and possible interventions. Prog Neurobiol. 2012;98:145–65.

    Article  PubMed  PubMed Central  Google Scholar 

  9. Ma Q, Xiong F, Zhang L. Gestational hypoxia and epigenetic programming of brain development disorders. Drug Discov Today. 2014;19:1883–96.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  10. Li Y, Ma Q, Dasgupta C, Halavi S, Hartman RE, Xiao D, Zhang L. Inhibition of DNA methylation in the developing rat brain disrupts sexually dimorphic neurobehavioral phenotypes in adulthood. Mol Neurobiol. 2017;54:3988–99.

    Article  CAS  PubMed  Google Scholar 

  11. Li Y, Xiao D, Yang S, Zhang L. Promoter methylation represses AT2R gene and increases brain hypoxic-ischemic injury in neonatal rats. Neurobiol Dis. 2013;60:32–8.

    Article  CAS  PubMed  Google Scholar 

  12. Gonzalez-Rodriguez PJ, Xiong F, Li Y, Zhou J, Zhang L. Fetal hypoxia increases vulnerability of hypoxic-ischemic brain injury in neonatal rats: role of glucocorticoid receptors. Neurobiol Dis. 2014;65:172–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  13. Shen G, Hu S, Zhao Z, Zhang L, Ma Q. Antenatal hypoxia accelerates the onset of Alzheimer’s disease pathology in 5xFAD mouse model. Front Aging Neurosci. 2020;12:251.

    Article  PubMed  PubMed Central  Google Scholar 

  14. Angevine JB Jr, Sidman RL. Autoradiographic study of cell migration during histogenesis of cerebral cortex in the mouse. Nature. 1961;192:766–8.

    Article  PubMed  Google Scholar 

  15. Rodier PM. Chronology of neuron development: animal studies and their clinical implications. Dev Med Child Neurol. 1980;22:525–45.

    Article  CAS  PubMed  Google Scholar 

  16. Azevedo FA, Carvalho LR, Grinberg LT, Farfel JM, Ferretti RE, Leite RE, Jacob Filho W, Lent R, Herculano-Houzel S. Equal numbers of neuronal and nonneuronal cells make the human brain an isometrically scaled-up primate brain. J Comp Neurol. 2009;513:532–41.

    Article  PubMed  Google Scholar 

  17. Pelvig DP, Pakkenberg H, Stark AK, Pakkenberg B. Neocortical glial cell numbers in human brains. Neurobiol Aging. 2008;29:1754–62.

    Article  CAS  PubMed  Google Scholar 

  18. Chen VS, Morrison JP, Southwell MF, Foley JF, Bolon B, Elmore SA. Histology atlas of the developing prenatal and postnatal mouse central nervous system, with emphasis on prenatal days E7.5 to E18.5. Toxicol Pathol. 2017;45:705–44.

    Article  PubMed  PubMed Central  Google Scholar 

  19. Rowitch DH, Kriegstein AR. Developmental genetics of vertebrate glial-cell specification. Nature. 2010;468:214–22.

    Article  CAS  PubMed  Google Scholar 

  20. Kessaris N, Pringle N, Richardson WD. Specification of CNS glia from neural stem cells in the embryonic neuroepithelium. Philos Trans R Soc Lond B Biol Sci. 2008;363:71–85.

    Article  CAS  PubMed  Google Scholar 

  21. Menassa DA, Gomez-Nicola D. Microglial dynamics during human brain development. Front Immunol. 2018;9:1014.

    Article  PubMed  PubMed Central  Google Scholar 

  22. Schulz C, Gomez Perdiguero E, Chorro L, Szabo-Rogers H, Cagnard N, Kierdorf K, Prinz M, Wu B, Jacobsen SE, Pollard JW, et al. A lineage of myeloid cells independent of Myb and hematopoietic stem cells. Science. 2012;336:86–90.

    Article  CAS  PubMed  Google Scholar 

  23. Reemst K, Noctor SC, Lucassen PJ, Hol EM. The indispensable roles of microglia and astrocytes during brain development. Front Hum Neurosci. 2016;10:566.

    Article  PubMed  PubMed Central  Google Scholar 

  24. Lawrence KM, McGovern PE, Mejaddam A, Rossidis AC, Baumgarten H, Kim AG, Grinspan JB, Licht DJ, Radaelli E, Rychik J, et al. Prenatal hypoxemia alters microglial morphology in fetal sheep. J Thorac Cardiovasc Surg. 2019;1:112.

    Google Scholar 

  25. Nalivaeva NN, Turner AJ, Zhuravin IA. Role of Prenatal Hypoxia in Brain Development, Cognitive Functions, and Neurodegeneration. Front Neurosci. 2018;12:825.

    Article  PubMed  PubMed Central  Google Scholar 

  26. Wang B, Zeng H, Liu J, Sun M. Effects of prenatal hypoxia on nervous system development and related diseases. Front Neurosci. 2021;15: 755554.

    Article  PubMed  PubMed Central  Google Scholar 

  27. Rees S, Harding R, Walker D. The biological basis of injury and neuroprotection in the fetal and neonatal brain. Int J Dev Neurosci. 2011;29:551–63.

    Article  PubMed  PubMed Central  Google Scholar 

  28. Riddle A, Luo NL, Manese M, Beardsley DJ, Green L, Rorvik DA, Kelly KA, Barlow CH, Kelly JJ, Hohimer AR, Back SA. Spatial heterogeneity in oligodendrocyte lineage maturation and not cerebral blood flow predicts fetal ovine periventricular white matter injury. J Neurosci. 2006;26:3045–55.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  29. Back SA, Han BH, Luo NL, Chricton CA, Xanthoudakis S, Tam J, Arvin KL, Holtzman DM. Selective vulnerability of late oligodendrocyte progenitors to hypoxia-ischemia. J Neurosci. 2002;22:455–63.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  30. Guo R, Hou W, Dong Y, Yu Z, Stites J, Weiner CP. Brain injury caused by chronic fetal hypoxemia is mediated by inflammatory cascade activation. Reprod Sci. 2010;17:540–8.

    Article  CAS  PubMed  Google Scholar 

  31. Mirza MA, Ritzel R, Xu Y, McCullough LD, Liu F. Sexually dimorphic outcomes and inflammatory responses in hypoxic-ischemic encephalopathy. J Neuroinflammation. 2015;12:32.

    Article  PubMed  PubMed Central  Google Scholar 

  32. Song H, Telugu BP, Thompson LP. Sexual dimorphism of mitochondrial function in the hypoxic guinea pig placenta. Biol Reprod. 2019;100:208–16.

    Article  PubMed  Google Scholar 

  33. Andrews S. FastQC. https://www.bioinformatics.babraham.ac.uk/projects/fastqc/. Accessed 23 Sep 2021.

  34. Li Y, Ma Q, Halavi S, Concepcion K, Hartman RE, Obenaus A, Xiao D, Zhang L. Fetal stress-mediated hypomethylation increases the brain susceptibility to hypoxic-ischemic injury in neonatal rats. Exp Neurol. 2016;275(Pt 1):1–10.

    Article  CAS  PubMed  Google Scholar 

  35. Semple BD, Blomgren K, Gimlin K, Ferriero DM, Noble-Haeusslein LJ. Brain development in rodents and humans: Identifying benchmarks of maturation and vulnerability to injury across species. Prog Neurobiol. 2013;106–107:1–16.

    Article  PubMed  Google Scholar 

  36. Dobbing J, Sands J. Comparative aspects of the brain growth spurt. Early Hum Dev. 1979;3:79–83.

    Article  CAS  PubMed  Google Scholar 

  37. Bockhorst KH, Narayana PA, Liu R, Ahobila-Vijjula P, Ramu J, Kamel M, Wosik J, Bockhorst T, Hahn K, Hasan KM, Perez-Polo JR. Early postnatal development of rat brain: in vivo diffusion tensor imaging. J Neurosci Res. 2008;86:1520–8.

    Article  CAS  PubMed  Google Scholar 

  38. Catalani A, Sabbatini M, Consoli C, Cinque C, Tomassoni D, Azmitia E, Angelucci L, Amenta F. Glial fibrillary acidic protein immunoreactive astrocytes in developing rat hippocampus. Mech Ageing Dev. 2002;123:481–90.

    Article  CAS  PubMed  Google Scholar 

  39. Craig A, Ling Luo N, Beardsley DJ, Wingate-Pearse N, Walker DW, Hohimer AR, Back SA. Quantitative analysis of perinatal rodent oligodendrocyte lineage progression and its correlation with human. Exp Neurol. 2003;181:231–40.

    Article  PubMed  Google Scholar 

  40. Ma Q, Dasgupta C, Shen G, Li Y, Zhang L. MicroRNA-210 downregulates TET2 and contributes to inflammatory response in neonatal hypoxic-ischemic brain injury. J Neuroinflammation. 2021;18:6.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  41. Diz-Chaves Y, Astiz M, Bellini MJ, Garcia-Segura LM. Prenatal stress increases the expression of proinflammatory cytokines and exacerbates the inflammatory response to LPS in the hippocampal formation of adult male mice. Brain Behav Immun. 2013;28:196–206.

    Article  CAS  PubMed  Google Scholar 

  42. Slusarczyk J, Trojan E, Glombik K, Budziszewska B, Kubera M, Lason W, Popiolek-Barczyk K, Mika J, Wedzony K, Basta-Kaim A. Prenatal stress is a vulnerability factor for altered morphology and biological activity of microglia cells. Front Cell Neurosci. 2015;9:82.

    Article  PubMed  PubMed Central  Google Scholar 

  43. Desplats P, Gutierrez AM, Antonelli MC, Frasch MG. Microglial memory of early life stress and inflammation: susceptibility to neurodegeneration in adulthood. Neurosci Biobehav Rev. 2020;117:232–42.

    Article  PubMed  Google Scholar 

  44. Taetzsch T, Benusa S, Levesque S, Mumaw CL, Block ML. Loss of NF-kappaB p50 function synergistically augments microglial priming in the middle-aged brain. J Neuroinflammation. 2019;16:60.

    Article  PubMed  PubMed Central  Google Scholar 

  45. Landskron G, De la Fuente M, Thuwajit P, Thuwajit C, Hermoso MA. Chronic inflammation and cytokines in the tumor microenvironment. J Immunol Res. 2014;2014: 149185.

    Article  PubMed  PubMed Central  Google Scholar 

  46. Wang P, Wu P, Siegel MI, Egan RW, Billah MM. Interleukin (IL)-10 inhibits nuclear factor kappa B (NF kappa B) activation in human monocytes. IL-10 and IL-4 suppress cytokine synthesis by different mechanisms. J Biol Chem. 1995;1:9558–63.

    Article  Google Scholar 

  47. Hsueh YS, Chang HH, Shan YS, Sun HS, Fletcher JA, Li CF, Chen LT. Nuclear KIT induces a NFKBIB-RELA-KIT autoregulatory loop in imatinib-resistant gastrointestinal stromal tumors. Oncogene. 2019;38:6550–65.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  48. Crotti A, Ransohoff RM. Microglial physiology and pathophysiology: insights from genome-wide transcriptional profiling. Immunity. 2016;44:505–15.

    Article  CAS  PubMed  Google Scholar 

  49. Fan H, Lu B, Cao C, Li H, Yang D, Huang L, Ding T, Wu M, Lu G. Plasma TNFSF13B and TNFSF14 function as inflammatory indicators of severe adenovirus pneumonia in pediatric patients. Front Immunol. 2020;11: 614781.

    Article  CAS  PubMed  Google Scholar 

  50. Patel JR, Klein RS. Mediators of oligodendrocyte differentiation during remyelination. FEBS Lett. 2011;585:3730–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  51. Girard S, Sebire G, Kadhim H. Proinflammatory orientation of the interleukin 1 system and downstream induction of matrix metalloproteinase 9 in the pathophysiology of human perinatal white matter damage. J Neuropathol Exp Neurol. 2010;69:1116–29.

    Article  CAS  PubMed  Google Scholar 

  52. Girard S, Kadhim H, Larouche A, Roy M, Gobeil F, Sebire G. Pro-inflammatory disequilibrium of the IL-1 beta/IL-1ra ratio in an experimental model of perinatal brain damages induced by lipopolysaccharide and hypoxia-ischemia. Cytokine. 2008;43:54–62.

    Article  CAS  PubMed  Google Scholar 

  53. Favrais G, van de Looij Y, Fleiss B, Ramanantsoa N, Bonnin P, Stoltenburg-Didinger G, Lacaud A, Saliba E, Dammann O, Gallego J, et al. Systemic inflammation disrupts the developmental program of white matter. Ann Neurol. 2011;70:550–65.

    Article  CAS  PubMed  Google Scholar 

  54. Kelly SB, Stojanovska V, Zahra VA, Moxham A, Miller SL, Moss TJM, Hooper SB, Nold MF, Nold-Petry CA, Dean JM, et al. Interleukin-1 blockade attenuates white matter inflammation and oligodendrocyte loss after progressive systemic lipopolysaccharide exposure in near-term fetal sheep. J Neuroinflammation. 2021;18:189.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  55. Mason JL, Suzuki K, Chaplin DD, Matsushima GK. Interleukin-1beta promotes repair of the CNS. J Neurosci. 2001;21:7046–52.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  56. Vela JM, Molina-Holgado E, Arevalo-Martin A, Almazan G, Guaza C. Interleukin-1 regulates proliferation and differentiation of oligodendrocyte progenitor cells. Mol Cell Neurosci. 2002;20:489–502.

    Article  CAS  PubMed  Google Scholar 

  57. Bowen S, Ateh DD, Deinhardt K, Bird MM, Price KM, Baker CS, Robson JC, Swash M, Shamsuddin W, Kawar S, et al. The phagocytic capacity of neurones. Eur J Neurosci. 2007;25:2947–55.

    Article  PubMed  Google Scholar 

  58. Dhuriya YK, Sharma D. Necroptosis: a regulated inflammatory mode of cell death. J Neuroinflammation. 2018;15:199.

    Article  PubMed  PubMed Central  Google Scholar 

  59. Degterev A, Hitomi J, Germscheid M, Ch’en IL, Korkina O, Teng X, Abbott D, Cuny GD, Yuan C, Wagner G, et al. Identification of RIP1 kinase as a specific cellular target of necrostatins. Nat Chem Biol. 2008;4:313–21.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  60. Li J, Cao F, Yin HL, Huang ZJ, Lin ZT, Mao N, Sun B, Wang G. Ferroptosis: past, present and future. Cell Death Dis. 2020;11:88.

    Article  PubMed  PubMed Central  Google Scholar 

  61. Lumsden AL, Rogers JT, Majd S, Newman M, Sutherland GT, Verdile G, Lardelli M. Dysregulation of neuronal iron homeostasis as an alternative unifying effect of mutations causing familial Alzheimer’s Disease. Front Neurosci. 2018;12:533.

    Article  PubMed  PubMed Central  Google Scholar 

  62. Xing C, Wang X, Cheng C, Montaner J, Mandeville E, Leung W, van Leyen K, Lok J, Wang X, Lo EH. Neuronal production of lipocalin-2 as a help-me signal for glial activation. Stroke. 2014;45:2085–92.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  63. Jha MK, Lee S, Park DH, Kook H, Park KG, Lee IK, Suk K. Diverse functional roles of lipocalin-2 in the central nervous system. Neurosci Biobehav Rev. 2015;49:135–56.

    Article  CAS  PubMed  Google Scholar 

  64. Zhao N, Xu X, Jiang Y, Gao J, Wang F, Xu X, Wen Z, Xie Y, Li J, Li R, et al. Lipocalin-2 may produce damaging effect after cerebral ischemia by inducing astrocytes classical activation. J Neuroinflammation. 2019;16:168.

    Article  PubMed  PubMed Central  Google Scholar 

  65. Lee S, Lee WH, Lee MS, Mori K, Suk K. Regulation by lipocalin-2 of neuronal cell death, migration, and morphology. J Neurosci Res. 2012;90:540–50.

    Article  CAS  PubMed  Google Scholar 

  66. Tang D, Chen X, Kang R, Kroemer G. Ferroptosis: molecular mechanisms and health implications. Cell Res. 2021;31:107–25.

    Article  CAS  PubMed  Google Scholar 

  67. Ivanov I, Kuhn H, Heydeck D. Structural and functional biology of arachidonic acid 15-lipoxygenase-1 (ALOX15). Gene. 2015;573:1–32.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  68. Haynes RL, van Leyen K. 12/15-lipoxygenase expression is increased in oligodendrocytes and microglia of periventricular leukomalacia. Dev Neurosci. 2013;35:140–54.

    Article  CAS  PubMed  Google Scholar 

  69. Joshi YB, Giannopoulos PF, Chu J, Sperow M, Kirby LG, Abood ME, Pratico D. Absence of ALOX5 gene prevents stress-induced memory deficits, synaptic dysfunction and tauopathy in a mouse model of Alzheimer’s disease. Hum Mol Genet. 2014;23:6894–902.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  70. Rudhard Y, Sengupta Ghosh A, Lippert B, Bocker A, Pedaran M, Kramer J, Ngu H, Foreman O, Liu Y, Lewcock JW. Identification of 12/15-lipoxygenase as a regulator of axon degeneration through high-content screening. J Neurosci. 2015;35:2927–41.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  71. Blumenthal I. Periventricular leucomalacia: a review. Eur J Pediatr. 2004;163:435–42.

    Article  PubMed  Google Scholar 

  72. Leviton A, Gilles F. Ventriculomegaly, delayed myelination, white matter hypoplasia, and “periventricular” leukomalacia: how are they related? Pediatr Neurol. 1996;15:127–36.

    Article  CAS  PubMed  Google Scholar 

  73. Volpe JJ. Neurobiology of periventricular leukomalacia in the premature infant. Pediatr Res. 2001;50:553–62.

    Article  CAS  PubMed  Google Scholar 

  74. Ma Q, Zhang L. Epigenetic programming of hypoxic-ischemic encephalopathy in response to fetal hypoxia. Prog Neurobiol. 2015;124:28–48.

    Article  CAS  PubMed  Google Scholar 

  75. George S, Scotter J, Dean JM, Bennet L, Waldvogel HJ, Guan J, Faull RL, Gunn AJ. Induced cerebral hypothermia reduces post-hypoxic loss of phenotypic striatal neurons in preterm fetal sheep. Exp Neurol. 2007;203:137–47.

    Article  CAS  PubMed  Google Scholar 

  76. Bennet L, Roelfsema V, George S, Dean JM, Emerald BS, Gunn AJ. The effect of cerebral hypothermia on white and grey matter injury induced by severe hypoxia in preterm fetal sheep. J Physiol. 2007;578:491–506.

    Article  CAS  PubMed  Google Scholar 

  77. Dean JM, Gunn AJ, Wassink G, George S, Bennet L. Endogenous alpha2-adrenergic receptor-mediated neuroprotection after severe hypoxia in preterm fetal sheep. Neuroscience. 2006;142:615–28.

    Article  CAS  PubMed  Google Scholar 

  78. Zaghloul N, Ahmed M. Pathophysiology of periventricular leukomalacia: what we learned from animal models. Neural Regen Res. 2017;12:1795–6.

    Article  PubMed  PubMed Central  Google Scholar 

  79. Hanamsagar R, Alter MD, Block CS, Sullivan H, Bolton JL, Bilbo SD. Generation of a microglial developmental index in mice and in humans reveals a sex difference in maturation and immune reactivity. Glia. 2017;65:1504–20.

    Article  PubMed  PubMed Central  Google Scholar 

  80. Blutstein T, Castello MA, Viechweg SS, Hadjimarkou MM, McQuail JA, Holder M, Thompson LP, Mong JA. Differential responses of hippocampal neurons and astrocytes to nicotine and hypoxia in the fetal guinea pig. Neurotox Res. 2013;24:80–93.

    Article  CAS  PubMed  Google Scholar 

  81. Patterson AJ, Chen M, Xue Q, Xiao D, Zhang L. Chronic prenatal hypoxia induces epigenetic programming of PKC{epsilon} gene repression in rat hearts. Circ Res. 2010;107:365–73.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  82. Lv J, Ma Q, Dasgupta C, Xu Z, Zhang L. Antenatal hypoxia and programming of glucocorticoid receptor expression in the adult rat heart. Front Physiol. 2019;10:323.

    Article  PubMed  PubMed Central  Google Scholar 

  83. Xiong F, Lin T, Song M, Ma Q, Martinez SR, Lv J, MataGreenwood E, Xiao D, Xu Z, Zhang L. Antenatal hypoxia induces epigenetic repression of glucocorticoid receptor and promotes ischemic-sensitive phenotype in the developing heart. J Mol Cell Cardiol. 2016;91:160–71.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  84. Krueger F. Trim Galore. https://www.bioinformatics.babraham.ac.uk/projects/trim_galore/. Accessed 1 Apr 2021.

  85. Martin M. Cutadapt removes adapter sequences from high-throughput sequencing reads. EMBnet journal. 2011;17:10–2.

    Article  Google Scholar 

  86. Kim D, Pertea G, Trapnell C, Pimentel H, Kelley R, Salzberg SL. TopHat2: accurate alignment of transcriptomes in the presence of insertions, deletions and gene fusions. Genome Biol. 2013;14:1123.

    Article  Google Scholar 

  87. Karolchik D, Hinrichs AS, Furey TS, Roskin KM, Sugnet CW, Haussler D, Kent WJ. The UCSC table browser data retrieval tool. Nucleic Acids Res. 2004;32:D493–6.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  88. Trapnell C, Williams BA, Pertea G, Mortazavi A, Kwan G, van Baren MJ, Salzberg SL, Wold BJ, Pachter L. Transcript assembly and quantification by RNA-Seq reveals unannotated transcripts and isoform switching during cell differentiation. Nat Biotechnol. 2010;28:511–5.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  89. Trapnell C, Roberts A, Goff L, Pertea G, Kim D, Kelley DR, Pimentel H, Salzberg SL, Rinn JL, Pachter L. Differential gene and transcript expression analysis of RNA-seq experiments with tophat and cufflinks. Nat Protoc. 2012;7:562–78.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  90. Chen EY, Tan CM, Kou Y, Duan Q, Wang Z, Meirelles GV, Clark NR, Ma’ayan A. Enrichr: interactive and collaborative HTML5 gene list enrichment analysis tool. BMC Bioinformatics. 2013;14:128.

    Article  PubMed  PubMed Central  Google Scholar 

  91. Kanehisa M, Goto S. KEGG: kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 2000;28:27–30.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  92. Gao T, Qian J. EnhancerAtlas 2.0: an updated resource with enhancer annotation in 586 tissue/cell types across nine species. Nucleic Acids Res. 2020;48:58–64.

    CAS  Google Scholar 

  93. Przanowski P, Mondal SS, Cabaj A, Debski KJ, Wojtas B, Gielniewski B, Kaza B, Kaminska B, Dabrowski M. Open chromatin landscape of rat microglia upon proinvasive or inflammatory polarization. Glia. 2019;67:2312–28.

    Article  PubMed  Google Scholar 

  94. Kent WJ, Sugnet CW, Furey TS, Roskin KM, Pringle TH, Zahler AM, Haussler D. The human genome browser at UCSC. Genome Res. 2002;12:996–1006.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  95. Grant CE, Bailey TL, Noble WS. FIMO: scanning for occurrences of a given motif. Bioinformatics. 2011;27:1017–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  96. Huang L, Chen X, Dasgupta C, Chen W, Song R, Wang C, Zhang L. Foetal hypoxia impacts methylome and transcriptome in developmental programming of heart disease. Cardiovasc Res. 2019;115:1306–19.

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

None.

Funding

This work was supported in part by the National Institutes of Health grants NS103017 to LZ and NS112404 to QM.

Author information

Author notes

  1. Isaac Kremsky and Qingyi Ma contributed equally

Authors and Affiliations

  1. Department of Basic Sciences, Loma Linda University School of Medicine, Loma Linda, CA, USA

    Isaac Kremsky, Qingyi Ma, Xin Chen, Samir Ali, Shawnee Angeloni, Charles Wang & Lubo Zhang

  2. Center for Genomics, Loma Linda University School of Medicine, Loma Linda, CA, USA

    Isaac Kremsky, Xin Chen & Charles Wang

  3. Lawrence D. Longo MD Center for Perinatal Biology, Department of Basic Sciences, Loma Linda University School of Medicine, Loma Linda, CA, 92350, USA

    Qingyi Ma, Bo Li, Chiranjib Dasgupta & Lubo Zhang

Authors
  1. Isaac Kremsky
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Qingyi Ma
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Bo Li
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Chiranjib Dasgupta
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. Xin Chen
    View author publications

    You can also search for this author in PubMed Google Scholar

  6. Samir Ali
    View author publications

    You can also search for this author in PubMed Google Scholar

  7. Shawnee Angeloni
    View author publications

    You can also search for this author in PubMed Google Scholar

  8. Charles Wang
    View author publications

    You can also search for this author in PubMed Google Scholar

  9. Lubo Zhang
    View author publications

    You can also search for this author in PubMed Google Scholar

Contributions

LZ conceived and designed the study. BL and CD performed experiments. CW oversaw transcriptomic sequencing. IK managed Bioinformatics analysis. S Ali, S Angeloni, XC and IK performed Bioinformatics analysis. IK, QM and LZ wrote the manuscript. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Lubo Zhang.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1: Figure S1.

Quality control of RNA-seq data. Figure S2. Related to Fig. 2. Figure S3. Related to Fig. 3. Figure S4. Related to Fig. 4.

Additional file 2: Table S1.

Hypoxia vs. Normoxia DEG (hDEG) counts. Table S2. Male vs. Female DEG counts. Table S3. Interaction DEG counts. Table S4. PCR primer details.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Kremsky, I., Ma, Q., Li, B. et al. Fetal hypoxia results in sex- and cell type-specific alterations in neonatal transcription in rat oligodendrocyte precursor cells, microglia, neurons, and oligodendrocytes. Cell Biosci 13, 58 (2023). https://doi.org/10.1186/s13578-023-01012-8

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13578-023-01012-8

Keyword

Cell & Bioscience

ISSN: 2045-3701

Contact us