- 1Department of Obstetrics and Gynaecology, Faculty of Medicine, CHU Ste-Justine Research Centre, Université de Montréal, Montréal, QC, Canada
- 2Department of Neurosciences, Faculty of Medicine, CHU Ste-Justine Research Centre, Université de Montréal, Montréal, QC, Canada
- 3Faculty of Veterinary Medicine, Animal Reproduction Research Centre, Université de Montréal, Montréal, QC, Canada
- 4Neuroimmunology Unit, Montréal Neurological Institute, McGill University, Montréal, QC, Canada
- 5Department of Clinical Sciences, Faculty of Veterinary Medicine, Université de Montréal, QC, Canada
- 6Département de Pathologie, University Hospital Ste-Justine, Université de Montréal, QC, Canada
- 7Department of Pediatrics, Faculty of Medicine, CHU Ste-Justine Research Centre, Université de Montréal, Montréal, QC, Canada
Objective: Neuroinflammation in utero may result in life-long neurological disabilities. The molecular mechanisms whereby microglia contribute to this response remain incompletely understood.
Methods: Lipopolysaccharide (LPS) or saline were administered intravenously to non-anesthetized chronically instrumented near-term fetal sheep to model fetal inflammation in vivo. Microglia were then isolated from in vivo LPS and saline (naïve) exposed animals. To mimic the second hit of neuroinflammation, these microglia were then re-exposed to LPS in vitro. Cytokine responses were measured in vivo and subsequently in vitro in the primary microglia cultures derived from these animals. We sequenced the whole transcriptome of naïve and second hit microglia and profiled their genetic expression to define molecular pathways disrupted during neuroinflammation.
Results: In vivo LPS exposure resulted in IL-6 increase in fetal plasma 3 h post LPS exposure. Even though not histologically apparent, microglia acquired a pro-inflammatory phenotype in vivo that was sustained and amplified in vitro upon second hit LPS exposure as measured by IL-1β response in vitro and RNAseq analyses. While NFKB and Jak-Stat inflammatory pathways were up regulated in naïve microglia, heme oxygenase 1 (HMOX1) and Fructose-1,6-bisphosphatase (FBP) genes were uniquely differentially expressed in the second hit microglia. Compared to the microglia exposed to LPS in vitro only, the transcriptome of the in vivo LPS pre-exposed microglia showed a diminished differential gene expression in inflammatory and metabolic pathways prior and upon re-exposure to LPS in vitro. Notably, this desensitization response was also observed in histone deacetylases (HDAC) 1, 2, 4, and 6. Microglial calreticulin/LRP genes implicated in microglia-neuronal communication relevant for the neuronal development were up regulated in second hit microglia.
Discussion: We identified a unique HMOX1down and FBPup phenotype of microglia exposed to the double-hit suggesting interplay of inflammatory and metabolic pathways. Our findings suggest that epigenetic mechanisms mediate this immunological and metabolic memory of the prior inflammatory insult relevant to neuronal development and provide new therapeutic targets for early postnatal intervention to prevent brain injury.
Introduction
Brain injury acquired antenatally remains a major cause of long-term neurodevelopmental sequelae (Saigal and Doyle, 2008). There is growing clinical and experimental evidence for maternal and fetal infection acting via systemic and neuroinflammation to cause fetal brain injury or contributing to in utero asphyxial brain injury with consequences for postnatal health (Hagberg et al., 2002; Rees and Inder, 2005; Wang et al., 2006; Gotsch et al., 2007; Murthy and Kennea, 2007; Fahey, 2008).
In humans, the main cause of fetal inflammation is chorioamnionitis, a frequent condition affecting 10% of all pregnancies and up to 40% of preterm births. Chorioamnionitis is associated with ~nine-fold increased risk for cerebral palsy spectrum disorders with life lasting neurological deficits and an increased risk for acute or life-long morbidity and mortality (Fahey, 2008; Agrawal and Hirsch, 2012; Fishman and Gelber, 2012).
In addition to short-term brain damage, neuroimmune responses to in utero infection may also have long-term health consequences, the “second hit” hypothesis: In adults, exposure to inflammatory stimuli can activate microglia (glial priming, reviewed in Billiards et al., 2006; Karrow, 2006; Bilbo and Schwarz, 2009; Bilbo and Tsang, 2010; Ajmone-Cat et al., 2013; Bolton et al., 2014).
We hypothesized that an inflammatory response induced by lipopolysaccharide (LPS) will result in microglial activation reflecting neuroinflammation. To test the “second hit” hypothesis, we developed a protocol to culture fetal sheep microglia and re-expose them to LPS under in vitro conditions allowing a more mechanistic study of their phenotype.
Materials and Methods
Ethics Statement
This study was carried out in strict accordance with the recommendations in the Guide for the Care and Use of Laboratory Animals of the National Institutes of Health. The respective in vivo and in vitro protocols were approved by the Committee on the Ethics of Animal Experiments of the Université de Montréal (Permit Number: 10-Rech-1560).
Anesthesia and Surgical Procedure
We instrumented pregnant time-dated ewes at 126 days of gestation (dGA, ~0.86 gestation) with arterial, venous and amniotic catheters and ECG electrodes (Frasch et al., 2007). Ovine singleton fetuses of mixed breed were surgically instrumented with sterile technique under general anesthesia (both ewe and fetus). In case of twin pregnancy the larger fetus was chosen based on palpating and estimating the intertemporal diameter. The total duration of the procedure was approximately 2 h. Antibiotics were administered to the mother intravenously (trimethoprim sulfadoxine 5 mg/kg body weight) as well as to the fetus intravenously and into the amniotic cavity (ampicillin 250 mg). Amniotic fluid lost during surgery was replaced with warm saline. The catheters exteriorized through the maternal flank were secured to the back of the ewe in a plastic pouch. For the duration of the experiment the ewe was returned to a metabolic cage, where she could stand, lie and eat ad libitum while we monitored the non-anesthetized fetus without sedating the mother. During postoperative recovery antibiotic administration was continued for 3 days. Arterial blood was sampled for evaluation of maternal and fetal condition and catheters were flushed with heparinized saline to maintain patency.
In vivo Experimental Protocol
Postoperatively, all animals were allowed 3 days to recover before starting the experiments. On these 3 days, at 9.00 am 3 mL arterial plasma sample were taken for blood gasses and cytokine analysis. Each experiment commenced at 9.00 am with a 1 h baseline measurement followed by the respective intervention as outlined below. FHR and arterial blood pressure was monitored continuously (CED, Cambridge, UK, and NeuroLog, Digitimer, Hertfordshire, UK). Blood samples (3 mL) were taken for arterial blood gasess, lactate, glucose, and base excess (ABL800Flex, Radiometer) and cytokines at the time points 0 (baseline), +1 (i.e., after LPS administration), +3, +6, +24, +48, and +54 h (i.e., before sacrifice at day 3). For the cytokine analysis, plasma was spun at 4°C (4 min, 4000 g, Eppendorf 5804R, Mississauga, ON), decanted and stored at −80°C for subsequent ELISAs. After the +54 h (Day 3) sampling, the animals were sacrificed with an overdose of barbiturate (30 mg pentobarbital sodium, Fatal-Plus; Vortech Pharmaceuticals, Dearborn, MI) and a post mortem was carried out during which fetal gender and weight were determined. The fetal brain was then perfusion-fixed with 250 mL of cold saline followed by 250 mL of 4% paraformaldehyde and processed for histochemical analysis or dissected for cell culture (details see in vitro microglia culture paragraph). Fetal growth was assessed by body, brain, liver, and maternal weights.
Nine fetuses were used as controls receiving NaCl 0.9%. Twelve fetuses received LPS (400 ng/fetus/day) derived from E. coli (Sigma L5293, from E. coli O111:B4, ready-made solution containing 1 mg/ml of LPS) were administered intravenously to fetuses on days 1 and 2 at 10.00 am to mimic high levels of endotoxin in fetal circulation over several days as it may occur in chorioamnionitis.
In vitro Microglia Culture Protocol
Fetal sheep brain tissues were obtained during sheep autopsy after completion of the experiment for in vitro study (Figure 1A). The non-instrumented, untreated twins were designated “naïve” (NC, no LPS exposure in vivo) and NL when exposed to LPS in vitro for the first time. Instrumented animals that received LPS in vivo (SHC) were used for 2nd hit LPS exposure in vitro (SHL). Fetal sheep microglia culture protocol was adapted from an established human adult and fetal microglia culture protocol that was modified to include a myelin removal step following high-speed centrifugation (Durafourt et al., 2013). Briefly, fetal sheep cells were plated on poly-L-lysine (PLL)-coated tissue culture flasks at a concentration of 2 × 106 cells /ml in DMEM with 5% heat-inactivated fetal bovine serum (Gibco, Canada Origin), 1% penicillin/ streptomycin, and 1% glutamine (5% DMEM), in which microglia are preferable to grow (Durafourt et al., 2013). Cells were allowed to incubate for seven days at 37°C, 5% CO2, followed by media change by centrifugation and addition of re-suspended cells back to the culture flask. Cells were continued to incubate for seven more days with 5% DMEM at 37°C, 5% CO2, before floating cells were collected. Carefully collecting the floating microglia to avoid contamination with astrocytes and oligodendrocytes, the cells were incubated in 24-well plates at 1 × 105 cells/1.82 cm2 surface area with 1 mL of 5% DMEM for another 4–5 days, and then treated with or w/o LPS (100 ng/ml, Sigma L5024, from E. coli O127, B8) for 6 h. Cell conditioned media were collected for cytokine analysis, 0.5 ml TriZol were added per well for RNA extraction.
Figure 1. Fetal sheep in vivo LPS exposure causes a systemic inflammation. (A) Experimental design. In vivo, in vitro and RNAseq experiments are illustrated. In vivo study, Control (saline) and LPS; in vitro study, cultured cells derived from in vivo Control animal, named as Naïve, whereas cells derived from LPS-exposed animal named as 2nd hit (second hit, SH), there are four experimental groups: naïve Control (NC), naïve LPS (NL), 2nd hit control (SHC), and 2nd hit LPS (SHL), respectively. For RNAseq data comparisons, we first compared pair 1 (n = 6) Control (NC) vs. LPS-exposed naïve microglia (NL); then pair 2 (n = 4) naïve control (NC) vs. 2nd hit Control (SHC); and finally pair 3 (n = 4) naïve LPS-exposed microglia (NL) vs. 2nd hit LPS-exposed microglia (SHL). (B) IL-6 levels peaked at 3 h in fetal sheep plasma following LPS administration in vivo (***P < 0.001). In vivo Control group, n = 9, in vivo LPS group, n = 12. Blood samples were collected in heparinized syringe from fetal arterial catheter, plasma was obtained by centrifugation. A sheep specific IL-6 ELISA was performed to measure the cytokine levels. (C) No evidence of Iba1+ fetal brain microglia inflammatory response to LPS exposure in vivo. Normalized Iba1 + signal (microglia) fluorescence in six randomly chosen high power fields per brain region is shown in hippocampus [CA1, CA234, and dentate gyrus (DG) subregions], cortical gray matter (GM) layers GM13, GM46, white matter (WM), and thalamus. GEE model for prediction of Iba1+ normalized signal intensity: group main effect p = 0.62; brain region main effect p < 0.001; brain region*group interaction p = 0.13.
To verify microglia purity, a portion of floating cells were cultured in 24-well plates at above conditions for flow cytometry analysis, cell morphology was documented with light microscopy (see Supplementary Material). Another portion of floating cells were plated into Lab-Tek eight well chamber glass slide (Thermo Scientific) and treated with or w/o LPS for immunocytochemistry analysis, in this experiment, some wells of astrocytes cultured at DMEM with 10% FCS were included for comparison.
Measurements of Inflammatory Responses
Measurement of Cytokines in Plasma and Cell Culture Media
Cytokine concentrations in plasma (IL-6) and cell culture media (IL-1β) were determined by using an ovine-specific sandwich ELISA. Briefly, 96-well plates (Nunc Maxisorp, high capacity microtitre wells) were pre-coated with the capture antibody, the mouse anti sheep monoclonal antibodies (IL-6, MCA1659; IL-1β, MCA1658, Bio Rad AbD Serotec) at a concentration 4 μg/ml on ELISA plates at 4°C for overnight, after 3 times wash with washing buffer (0.05% Tween 20 in PBS, PBST), plates were then blocked for 1 h with 1% BSA in PBST for plasma samples or 10% FBS for cell culture media. Recombinant sheep proteins (IL-6, Protein Express Cat. no 968-305; IL-1 β, Cat. no 968-405) were used as ELISA standard. All standards and samples (50 μl per well) were run in duplicate. Rabbit anti sheep polyclonal antibodies (detection antibody IL-6, AHP424; IL-1β, AHP423, Bio Rad AbD Serotec) at a concentration of 4 μg/ml were applied in wells and incubated for 30 min at room temperature. Plates were washed with washing buffer for 5–7 times between each step. Detection was accomplished by assessing the conjugated enzyme activity (goat anti-rabbit IgG-HRP, dilution 1:5000, Jackson ImmunoResearch, Cat. No 111-035-144) via incubation with TMB substrate solution (BD OptEIA TMB substrate Reagent Set, BD Biosciences Cat. No 555214), color development reaction was stopped with 2 N sulphuric acid. Plates were read on ELISA plate reader at 450 nm, with 570 nm wavelength correction (EnVision 2104 Multilabel Reader, Perkin Elmer). The sensitivity of IL-6 ELISA for plasma was 16 pg/ml, the sensitivity of IL-1b ELISA for media was 41.3 pg/ml, respectively. For all assays, the intra-assay and inter-assay coefficients of variance was <5 and <10%, respectively.
Immunofluorescence Imaging Analysis
Complete brain was taken from the fetus during necropsy after perfusion and immediately immersed in 4% PFA for 48–72 h. The tissue sample was then washed and stored in 1 × PBS buffer changed daily for 3 days. Finally, the brain was stored in 70% ethanol until further processing. All the brain tissue samples were kept at 4°C when they were in liquid. The fetal brains were cut into two equal halves of left and right hemispheres, and then sliced coronally and placed into cassettes to be processed with Leica TP 1020 Automatic Tissue Processor (Leica Instruments, Mussloch, Germany). The tissues were embedded in paraffin with Leica EG 1160 Paraffin Embedding Center (Leica Instruments, Mussloch, Germany). Five-micrometer slices were obtained from slicing the embedded tissue samples with the Leica RM2145 Rotary Microtome (Leica Instruments, Mussloch, Germany), and mounted on the Fisherbrand Colorfrost Plus microscope slides (Fischer Scientific). The sectioned brain tissue samples went through de-paraffinization with CitroSolv (Fischer Scientific), 100, 95, 70, and 50% ethanol at room temperature, and antigen retrieval with 10 mM citrate buffer at pH 6 before being washed with water and 1 × PBS, and blocked by Background Sniper Blocking Reagent (Biocare Medical, Cat. No BS966JJ). Then the sections were incubated with the primary antibody (Iba1, rabbit polyclonal antibody 4, 1:250 dilution, Wako, Cat No. 019-19741) for 1 h, followed by washing with 1× PBS and incubation with secondary antibody (Alexa Fluor 568 goat anti-rabbit IgG, 1:400 dilution, Life Technologies, Cat no A-11011) for 30 min in the dark. After that, the sections were washed again with 1× PBS, and the nuclei were counterstained with DAPI (1:4000 dilution, Sigma D-9564). Finally, the sections were cover-slipped with Fisherfinest Premium Cover glass (22 × 50-1, Fisher Scientific) and Fluoromount-G (SouthernBiotech, Cat no 0100-01) mounting medium, and viewed after 24 h of drying. Widefield fluorescence microscopy was performed on the stained brain tissue samples with a Zeiss Axiovert 200 M inverted microscope (Jena, Germany), at the magnification of 40× using a HBO100 mercury-arc lamp as a light source. The images were captured using a Zeiss Axiocam HRm (high-resolution monochrome) CCD (charged-coupled device) camera. Six high power field (HPF) images at 40× magnification were obtained for each animal. Multichannel imaging was used with the Iba1 channel and the DAPI channel for obtaining the pictures used for macrophage quantification. Appropriate ranges of color were selected showing positive contiguous cytoplasmic staining as a criterion for microglia cell count scoring which were then applied uniformly to calibrated images for all brain regions (Figure 2). Scoring was performed in a blinded fashion to experimental groups. To normalize for cell density Iba1+ signal over the whole area measured (100 sq micron) was divided by the respective optical intensity values for each HPF according to Lin et al. (2000).
Figure 2. Purity validation of fetal sheep brain primary microglia in cultures and LPS second hit. (A) Photomicrographs (ICC) confirming cell purity. Iba1+ staining in microglia vs. undetectable GFAP signal in microglia indicating no contamination with astrocytes in the culture. Microglia were cultured in eight-well chamber slide with DMEM +5% FCS for 4–7 days. Microglia were collected from the floating fraction and stained for Iba1 and GFAP. Scale bar = 50 μm. Magnification 40× for both images. (B) Purity of fetal sheep brain primary microglia cultures was verified by flowcytometry. After several days in culture, fetal sheep microglia were scraped from the wells using a cell scraper and blocked for 30 min using normal mouse IgG and 10% human serum. Cells were then stained using a FITC-conjugated monoclonal bovine anti-CD11b (1:40, Bio-Rad) on ice for 20 min. Cells were washed in FACS buffer and analyzed using a BD FACSCalibur and FlowJo software. (C) Microglia from in vivo LPS exposed brain appear more aggregated in vitro than microglia derived from in vivo controls. Microglia were cultured in 24-well plates with DMEM +5% FCS for 4–7 days, when images were taken. Cells were extracted from a twin control fetal brain and an in vivo LPS exposed fetal brain (Magnification 20× for both images). (D) Effect of “second-hit” in vitro LPS treatment on microglial phenotype. BOTTOM: IL-1β concentration in conditioned media of microglia derived from fetal sheep brain that were exposed to LPS vs. saline in vivo (***P < 0.0001). Cultures from in vivo LPS-exposed (SH), n = 4, cultures from in vivo Control (Naïve), n = 10. Cell culture media supernatant was obtained by centrifugation upon cell culture termination. A sheep specific IL-1β ELISA was performed to measure cytokine levels in cell culture media. In vitro, at baseline, microglia secreted more IL-1β in the in vivo LPS group (SHC) than in naïve Control (NC). LPS re-exposure (SHL) further increased IL-1β vs. naïve LPS (NL) by ~4.6-fold.
RNAseq Approach
To extract and quantify RNA, total RNA was extracted from cultured microglia using TRIzol Reagent (Life Technologies). To obtain enough RNA, same treatment cells were pooled in one replicate. RNA quantity and quality (RNA integrity number, RIN) was determined by using a RNA Nano Chips (Agilent RNA 6000 Nano Chips) with Agilent 2100 BioAnalyzer. All samples had a RIN-value ranging from 6 to 8.5, except for one sample having RIN = 5.5 but an acceptable 84% of transcripts mapped, which did not affect the read count for this sample.
A total of eight samples from four set of replicates were selected for RNA sequencing at high throughput, of which three replicates were derived from in vivo control fetal sheep and one replicate was from in vivo LPS-exposed (second hit) fetal sheep. RNAseq libraries were prepared by using Illumina TruSeq RNA Sample Preparation v2 kit (Illumina) and quality control was performed on the BioAnalyzer. Single-end 50-bp sequencing was performed at high throughput on an Illumina HiSeq2500 at the CHU Ste-Justine Core Facility Sequencing Platform. Raw data and RNAseq data discussed in this publication were deposited on NCBI and are accessible online with the GEO accession number GSE71037.
RNAseq Data Analysis
Reads alignment to the reference genome
To maximize the amount of genes covered, raw data were mapped to the reference genome of the sheep Ovis aris v3.1 from NCBI and Ensembl (GCA_000298735.1) as transcriptome reference. Index of the reference fasta file were built with Bowtie2 (Langmead and Salzberg, 2012), we then trimmed the adaptor of the fastQ files with TrimGalore, and mapped reads to the reference with Tophat2 (Kim et al., 2013). From the aligned reads from Tophat2, the number of reads per gene were counted with HTseq and assembled into a matrix containing the read count of each gene per sample (Anders et al., 2015).
Normalization and transcriptome analysis
Among packages available to test for differential expression, DESeq2 provides methods suited for the use of replicates; it uses negative binomial generalized linear models and estimates dispersion and logarithmic fold changes. We used DESeq2 to normalize the dataset, generated log2-fold changes and adjacent P-values (padj) and therefore, to find differentially expressed (DE) genes in microglia (Love et al., 2014). We first compared NC to NL to understand gene expression in naïve microglial cells after in vitro exposure to LPS. Then, due to the lack of replicates, we were not able to compare second hit microglial cells to their respective second hit control. Instead, we compared the genetic expression difference in response to a pre-exposure in vitro to LPS in NC vs. SHC. Finally, we assessed genetic expression in SHL compared to NL. A gene was considered differentially expressed (DE) if its adjacent p-value was strictly lower than 0.1. Pools of DE up and down regulated genes were clustered and visualized in heat maps, generated in R using the log2 normalized counts and the heatmap 0.2 method of the gplots library (Warnes et al., 2009).
Gene selection and Gene Ontology (GO)
The sheep genome is not yet supported by most gene ontology platforms, therefore, downstream analyses were performed with orthologs in the human genome Homo sapiens. ToppGenes and ToppCluster (Chen et al., 2009; Kaimal et al., 2010) were used to test for functional annotation enrichment analyses of biologic process and pathway with a false discovery rate correction of 0.05 (Franceschini et al., 2013). Gene Ontology was then performed with Gorilla and significant networks (P < 0.03) were selected for further discussion (Bauer et al., 2008; Eden et al., 2009).
Quantitative Real-Time PCR Analysis
The expression profiles of candidate genes were validated by real-time qRT-PCR. Total RNAs were subjected to cDNA synthesis using a QuantiTech Rev. Transcription Kit (Qiagen). HMOX1 and FBP mRNA were quantified by qRT-PCR using a QuantiFast SYBR Green PCR Kit (Qiagen) with STRATAGENE 3000 P, mRNA relative expression was calculated by the 2−ΔΔCt method over housekeeping gene GAPDH compared to baseline (Livak and Schmittgen, 2001). Sheep specific HMOX1 primers were designed with primer3 (Untergasser et al., 2012) and FBP primers were designed using Integrated DNA Technologies online tool and listed in Table 1.
Statistical Analyses
Generalized estimating equations (GEE) modeling was used to assess the effects of LPS while accounting for repeated measurements on fetal blood gasses and acid-base status, plasma and in vitro cytokines, cardiovascular responses [AR(1) correlation matrix to account for temporal structure] and in vivo Iba1+ fluorescence (independent correlation matrix to deal with the spatial distribution of Iba1+ fluorescence across the brain regions). We used a linear scale response model with LPS and time or brain regions as predicting factors to assess their interactions using maximum likelihood estimate and Type III analysis with Wald Chi-square statistics. Correlation analysis was performed using Spearman correlation coefficient. SPSS Version 21 was used for these analyses (IBM SPSS Statistics, IBM Corporation, Armonk, NY). Significance was assumed for p < 0.05. Results are provided as means ± SEM. Not all measurements were obtained for each animal studied.
Results
In vivo Studies
Cohorts' Characteristics
Maternal venous blood gasses, pH, and lactate did not significantly change during the experiments and were within physiological range throughout the experiment for both groups. Maternal and fetal cohort's characteristics are summarized in Table 2. Gestational age at time of the experimental day 1 averaged 130 days ± 1.3 dGA (term 145 dGA). Overall, mother and fetus were considered healthy based upon a physical examination and laboratory data collected.
Clinical-chemical Data
Clinical-chemical data, reported elsewhere, (Durosier et al., 2015) are summarized in Table 2 and were within physiological range for both groups. We found significant time*LPS interactions for pH (P = 0.03), pO2, pCO2, lactate, and BE (all P < 0.001).
Cardiovascular Analysis
As reported (Durosier et al., 2015), we found time-LPS interactions for mean arterial blood pressure and fetal heart rate responses (P = 0.015 and P < 0.001, respectively).
Plasma Cytokines Response to LPS
In vivo LPS exposure resulted in a peak of IL-6 at 3 h compared to baseline. We detected time-LPS interaction for fetal plasma IL-6 (P < 0.001, Figure 1B).
In vivo LPS Effect on Neuroinflammation
To assess the effect of in vivo LPS exposure on neuroinflammation in situ we measured microglial activation as Iba1+ immunofluorescence signal. We found a significant brain region main effect (p < 0.001), but no group main effect (p = 0.62) and no significant brain region*group interaction (p = 0.13) (Figure 1C). Thus, in vivo LPS exposure did not cause any measurable neuroinflammation as can be seen with higher doses of LPS using the same microglia marker (Keogh et al., 2012; Kuypers et al., 2013).
Overall, fetuses responded to the in vivo LPS exposure with signs of moderate sepsis as evident by the changes observed with arterial blood gas, pH, lactate, but with no signs of cardiovascular decompensation. Despite the systemic response to the LPS challenge, we observed no signs of neuroinflammation in vivo.
In vitro Studies
Having established a moderate LPS-induced in vivo fetal systemic inflammation paradigm without overt neuroinflammation in situ, we next aimed to test the functional properties of microglia exposed to LPS in vivo in an in vitro setting allowing characterization of microglial cytokine secretion and transcriptome profiles in response to LPS.
Primary Fetal Sheep Microglia Culture
In vitro studies were conducted in primary cultures derived from six controls (naïve) and from two in vivo LPS-exposed animals (SH). We were able to perform 1–2 in vitro replicates per each animal depending on cell numbers obtained.
We identified oligodendrocytes and neurons in the initial cell isolation in addition to microglia and astrocytes (data not shown). To enrich for microglia we subjected the cells to a second step as detailed in Methods. To verify cell culture purity, we performed immunofluorescence staining with a microglia marker confirming that the isolated primary microglia was very high (Figure 2A), whereas an astrocyte marker, GFAP, was absent from the cell population. To further verify cell purity, flow cytometry CD11b-FITC antibody was performed resulting in 96% of cultured cells are CD11b+ (Figure 2B), further indicating that a highly pure microglia population was obtained.
We used the purified highly enriched microglia cultures to pursue the second hit paradigm, i.e., how these cells behave in vitro in dependence on previous in vivo LPS exposure (Figure 1A). We found that microglia from in vivo LPS exposed fetal brain differ in morphology, showing more aggregation or clumping compared to naïve microglia (Figure 2C). This finding indicates that these microglia might have already been activated by LPS exposure in vivo.
Next, we investigated cytokine secretion properties of these cells in the absence or presence of LPS. For IL-1β, we found that in vitro LPS administration resulted in increased IL-1β in microglia compared to control cell cultures; this IL-1β response was potentiated by 4.6-fold in cells derived from animals with in vivo LPS exposure: 1884 ± 481 pg/ml vs. 406.14 ± 193 pg/ml (all p < 0.001, Figure 2D). Moreover, even in the absence of LPS in vitro, at baseline, microglia from the in vivo LPS group secreted more IL-1β (208.1 ± 16.63 pg/ml vs. 44.97 ± 59.21 pg/ml) with the fold increase being concordant with the level of gene expression increase (all P < 0.001, Table 3). Other pro-inflammatory cytokines of interest such as IL-6 and TNF-α were undetectable in cell-conditioned media (ELISA data not shown). Our findings suggest that a pro-inflammatory microglial phenotype acquired during in vivo exposure to LPS is sustained in vitro (second hit paradigm).
Table 3. Gene expression summary in naïve (one time exposure to LPS in vitro) and second hit (exposed once in vivo and second time in vitro) microglia.
RNAseq Approach
General overview of the whole transcriptome sequencing
To explore the genomic landscape of fetal sheep microglia, we sequenced the transcriptome of naïve and “second hit” microglia. As a quality control, we tested the expression levels of GFAP and TNFα across our three comparisons, and confirmed that all cells in our platform shared the same gene expression characteristics of microglia (Table 3). As a control measure for cell purity, our data confirmed the presence of TGF-β1 in each sample, as previously reported (Butovsky et al., 2014). To further confirm cell purity and the findings on protein level (ELISA), our transcriptome analysis showed a 1.654-fold increase of IL-1β (log2 = 0.726) between naïve and second hit LPS-exposed microglia (Table 3, respectively NL and SHL).
Firstly, we compared gene expression between the naïve controls and naïve LPS-exposed microglial cells. We found 258 differentially expressed genes (padj < 0.1), among which, 205 genes were up regulated and 53 were down regulated. We selected relevant differentially expressed genes with ToppCluster (logP>4.00) based on their role in the immune response (Figures 3A,B).
Figure 3. Heat maps of the gene expression in microglia cells exposed to LPS. Selected up regulated (red) and down regulated (blue) genes are listed; genes were selected with ToppCluster based on their significance in the immune response (logP>4.00) (A) Heat map of 205 differentially expressed (padj < 0.1) up regulated genes (red) in NL microglia. (B) Heat map of 53 differentially expressed down regulated genes in NL microglia, among selected genes indicated in blue, HMOX1 was strongly differentially expressed (log2 = −2.686 and ). In both up and down regulated genes, we observed a different behavior for NL3 that did not affect our differential analysis. (C) Heat map of the 4112 most differentially expressed and up regulated genes in SHC compared to NC cells. Selected genes with ToppCluster (logP>4.00) include GABA receptor genes and genes related to the transport of ion, Calcium, Sodium, and Potassium. IL1B was up regulated with log2 = 1.766 () corresponding to a 3.40-fold increase (D) Among the 2530 most differentially expressed and down regulated genes identified, HMOX1 was significantly down regulated (log2 = −4.462 and ). The reported genes were selected with ToppCluster (logP > 4.00) for their implication in the inflammatory response. (E) Differentially expressed up and down regulated genes in SHL compared to NL microglia. HMOX1 was the only down regulated gene (log2 = −4.303, ). When comparing common genes in SHC and SHL cells, the gene FBP was unique to SHL (red rectangle). NC, Naïve control microglia; NL, Naïve LPS-exposed microglia; SHC, Second hit control microglia; SHL, Second hit LPS-exposed microglia.
Then, to better understand the effect of an in vivo pre-exposure to LPS on biological processes, we compared gene expression between the naïve and second hit controls, i.e., NC and SHC, respectively. We found 6642 differentially expressed genes, among which, we identified 4112 up regulated and 2530 down regulated genes. Selection of relevant genes with ToppCluster (logP>4.00) showed that up regulated genes in SHC are mainly composed of GABA genes and genes responsible for calcium, potassium, and second messengers transport. Differentially expressed down regulated genes comprised the genes of the NF-κB signaling pathway and the HMOX1 gene, responsible for iron metabolism (log2 = −4.462 and ).
Finally, in an effort to discover the differences in response between NL and SHL, we compared gene expression of the NL set of three replicates and SHL (n = 1). We identified a total of six differentially expressed genes: five were up regulated and one gene, HMOX1, was strongly down regulated (HMOX1down, log2 = −4.303 and ). Among the five differentially up regulated genes identified, Fructose-1,6-bisphosphatase (FBP) was uniquely differentially expressed in second hit LPS-exposed microglia (FBPup, log2 = 4.057 and ). The expression profile of HMOX1 and FBP were assessed by quantitative real-time PCR (qRT-PCR). The results showed that the expressions of HMOX1 and FBP were consistent with the expressions from the transcriptome analyses (Figures 4A,B). The roles of these genes are discussed below.
Figure 4. Validation of RNAseq findings by qRT-PCR analysis of HMOX1 and FBP. (A) Quantification of HMOX1 in each NL microglia compared to averaged-NC. (B) Quantification of HMOX1 and FBP in SHL compared to averaged-NL microglia. Avg, Averaged; NC, Naïve control microglia; NL, Naïve LPS-exposed microglia; SHL, Second hit LPS-exposed microglia.
Discussion
We established for the first time an in vivo–in vitro endotoxin double-hit mammalian microglia experimental model to mimic multiple perinatal neuroinflammation episodes. The isolation of viable and highly purified microglia populations from in vivo LPS-exposed brain allowed an in vitro characterization of this cell type. Our most striking discovery was that the fetal inflammatory microglial phenotype acquired during in vivo exposure to LPS, even if not histologically apparent, is sustained and potentiated in vitro upon re-exposure to LPS. The subsequent RNA sequencing of the microglial genome revealed a unique HMOX1down and FBPup phenotype of microglia exposed to the double-hit, suggesting interplay of inflammatory and metabolic pathways.
In vivo–In vitro Model of Perinatal Inflammation Double-hit
Intrauterine exposure to inflammatory stimuli may switch innate immunity cells such as macrophages and microglia to a reactive phenotype (“priming”). Confronted with renewed inflammatory stimuli during labor or postnatally (especially in preterm neonates in the intensive care unit), such sensitized cells can sustain a chronic or exaggerated production of proinflammatory cytokines associated with neurodevelopmental deficits persisting into adulthood (double-hit hypothesis) (Larouche et al., 2005; Spencer et al., 2006; Wang et al., 2006).
Experimentally induced inflammation in chronically instrumented non-anesthetized fetal sheep is a well-established in vivo model of fetal physiology (Prout et al., 2010, 2012). Primary microglia cultures in different species have been reported for decades (Stansley et al., 2012). We integrated both in vivo and in vitro models into a new, hybrid system adding the layer of the whole transcriptome analysis using RNAseq analyses. The chief advantage of the new in vivo–in vitro model presented here is that it allows us to examine microglia responses to LPS-induced double-hit inflammation in situ and in vitro on integrative physiological, protein and genomic levels, and in a physiologically and clinically meaningful context. This approach has the potential to uncover hitherto unseen relationships between brain and immune system on different scales of organization in the perinatal stage of development, which might accelerate discovery of new treatment strategies.
In vivo, our experimental cohort's morphometric, arterial blood gasses, acid-base status and cardiovascular characteristics were within physiological range and representative for late-gestation fetal sheep as a model of human fetal development near term (Frasch et al., 2007; Rurak and Bessette, 2013). As reported elsewhere (Durosier et al., 2015), the effect of the low LPS dose we administered on the arterial blood gasses, acid-base status and cardiovascular responses is compatible with a mild septicemia (mild compensated metabolic acidemia and hypoxia) evidenced by a transient rise of IL-6 at 3 h without overt shock, i.e., without cardiovascular decompensation. Similar levels of systemic IL-6 have been reported (Prout et al., 2010).
In vitro, we developed a new microglia isolation protocol that combines the human adult and fetal brain microglia isolation protocols (Durafourt et al., 2013) and successfully collected a highly enriched microglia population. The use of a modified cell isolation approach from fetal brain tissue is mainly due to the higher degree of myelination in the adult human brain compared to that of a near-term fetus (Durafourt et al., 2013). We were able to attain high purity of microglia, which we validated by flow cytometry and immunocytochemistry. Moreover, RNAseq showed a consistent and constant low level of expression of the astrocyte marker GFAP further confirming cell purity.
In this study, the morphology of microglia from in vivo LPS exposed fetal brain was distinguished by more aggregation or clumping and less ramification compared to naïve microglia (Suzuki et al., 2006; Henkel et al., 2009). This suggests that microglia exposed to LPS in vivo may have been already activated before plating in cultures.
Single-hit LPS Exposure In Vitro Results in up Regulation of Inflammatory Pathways JAK-STAT and NFKB and Down Regulation of Metabolic Pathways
Gene ontology analyses of DE genes in NL microglia revealed an up regulation of inflammatory pathways NFκB, PIK3-Akt, and Jak-STAT. Interestingly, this was accompanied by a down regulation of metabolic pathways in LPS-induced inflammatory response (Figures 3A,B and Table 3). These findings may be explained, at least in part, by the emerging role of energy-sensing AMPK signaling in microglia, which links the inflammatory and metabolic regulatory networks (Frasch, 2014), We will return to this observation in Section Double-hit LPS Exposure of Microglia In vivo and In vitro Is Uniquely Characterized By a HMOX1down/FBPup Phenotype.
Among differentially expressed genes selected, NFKB (log2 = 2.676 and ) and JAK3 (log2 = 2.965 ) were up regulated in NL microglia. We then investigated the expression of genes involved in the NFκB and JAK-Stat pathways; our data showed that IL1B (log2 = 7.578), TNF (log2 = 4.990), NFKBIA (log2 = 2.578 and ), and RELB (log2 = 1.503) were up regulated in NL microglia. Gene ontology analysis revealed down regulation of the energy consuming processes and up regulation of energy conserving processes, as evidenced by the down regulation of genes related to glycolysis (GPI) and up regulation of gluconeogenesis (FBP) and the insulin signaling pathway. Furthermore, Gene Ontology of up regulated genes revealed that the GO term “immune system process” clustered key genes of inflammatory pathways, such as, JAK3, NFKBIA, and NFKBIB (GO:0002376 and P = 9.56 × 10−8). Differentially expressed down regulated genes also affected “the immune system process” (GO:0002376 and P = 3.24 × 10−4) and cellular response to metal ion (GO: 0071248 and P = 6.05 × 10−6). HMOX1 and FOS clustered in both GO terms underlying the potential role of HMOX1 in the immune system in relation with FOS. Analysis of all down regulated genes showed that the “metabolic process” (GO:0008152 and P = 2.38 × 10−8) was also globally affected (data not shown).
Double-hit LPS Exposure of Microglia In vivo and In vitro is Uniquely Characterized By a HMOX1down/FBPup Phenotype
Interestingly, HMOX1 gene expression showed a strong down regulation in SHL and SHC by four-fold (Table 3, Figures 3D,E). The level of expression of HMOX1 was higher in SHC than in SHL (log2 = −4.462 and ; log2 = −4.303 and , respectively). Such differences in response patterns were observed in other genes as well suggesting a memory of inflammation induced by pre-exposure to LPS in vivo.
HMOX1 role in microglia is yet to be fully understood. Across the three group comparisons (Figure 1A), HMOX1 was significantly down regulated and in vivo pre-exposure to LPS seemed to further enhance the down regulation of HMOX1 in response to second-hit in vitro LPS stimulation. We confirmed by RT-PCR that transcript amounts of HMOX1 are low in NL and SHL microglia (Figures 4A,B). HMOX1 was suggested to play an anti-inflammatory role in LPS-induced murine adult cell line macrophages via the activation of the Nrf2/ARE pathway (Ye et al., 2014). Pre-treatment with Oroxylin A, an inhibitor of LPS-induced mRNA, substantially increased the levels of NRF-2 and heme oxygenase 1. The response of SHC and SHL compared to their single-hit microglia counterpart NL showed that HMOX1 and NRF-2 had a greater down regulation after pre-exposure to LPS in vivo (Table 3), supporting the potential role of HMOX1 in the inflammatory response and as a determinant of microglial phenotype.
While the role of FBP in inflammation is unclear, its neuroprotective effect in brain injury models was suggested through various mechanisms. During hypoxia, FBP supports ATP production via stimulation of glycolysis which results in maintenance of normal intracellular calcium levels via the phospholipase-C dependent MAP kinase signaling pathway (Bickler and Kelleher, 1992; Fahlman et al., 2002). When comparing gene expression in SHC (Figures 3C,D), we observed that genes responsible for the transport of calcium, potassium and second messengers were also differentially expressed and up regulated.
It was previously reported that FBP dose-dependently suppressed LPS-induced nitric oxide (NO) production, and higher FBP doses were also associated with decreased levels of the transcription factor activator protein AP-1 in primary neonatal murine microglia cultures (Kim et al., 2012). We confirmed this observation in our comparison of SHL to SHC, wherein we observed that up regulation of FBP was concordant with lower expression of NOS1AP (log2 = −1.393) suggesting lower production of NO. We also observed that lower FBP transcripts amount in NL was accompanied by higher expression level of NOS production related genes (Table 3). The authors observed that FBP had an effect on the binding of transcription factors to DNA: FBP diminished the binding of AP-1 to DNA, but NFKB and CREB did not seem affected. We found down regulation of AP-1 (log2 = −2.385) and a slight down regulation of NFKB and CEBP, though CREB1 remained unaffected. In the SHL, our results confirmed that DNA binding nuclear factors were not strongly down regulated upon higher transcript level of FBP. The authors suggested that FBP inhibits iNOS expression by blocking the JNK/p38 MAPK pathway. We confirmed that JNK related genes may have lower expression level, however we did not observe any marked difference for P38.
A common theme within the newly found HMOX1down/FBPup phenotype appears to be its memory of the “energy restoring direction” following in vivo exposure to LPS. This metabolic effect is evidenced for example by up regulation of AMPK, insulin, growth arrest processes, mitochondrial biogenesis signaling pathways and down regulation of mTOR signaling pathway and such energy consuming processes as cell growth and protein synthesis (cf. Table 3).
Does In vivo Endotoxin Exposure Induce Transcriptome Memory of Inflammation in Fetal Microglia Mediated By Epigenetic Mechanisms?
Pre-exposure to LPS in vivo affected globally the transcriptome of microglia (Table 3). We observed that SHC microglia had a diminished response in gene expression of inflammatory pathways NF-κB, JAK-Stat, and PIK3-Akt compared to the behavior of the NL microglia; this phenomenon was sustained in SHL microglia. As mentioned in Section Double-hit LPS Exposure of Microglia in vivo and In vitro Is Uniquely Characterized By a HMOX1down/FBPup Phenotype, this was also true for the metabolic pathways. This desensitization was also observed in histone deacetylase 1 (HDAC1 and 6), which was DE up regulated by two-fold in NL (log2 = 2.271 and ), and up regulated by less than one-fold in SHC and SHL (log2 = 0.676 and log2 = 0.145, respectively, Table 3). This HDAC1 profile was accompanied by a less than one-fold down regulation of HDAC6 (log2 = −0.688 and ) followed again by desensitization in microglia exposed to LPS in vivo. HDAC4 was 2.5-fold up regulated (log2 = 1.292 and padj = 0.133) in NL vs. NC microglia followed also by desensitization in the comparison to the in vivo pre-exposed microglia; meanwhile, HDAC2 showed a less than one-fold down regulation (log2 = −2.746 and ) in NC vs. SHC microglia, with no detectable change in microglia exposed to LPS in vitro only or upon double-hit exposure. In parallel, HMOX1 was down regulated by four-fold in SHL and SHC, and by two-fold in NL microglia, and did not seem to have a diminished response in SHL. These findings underscore the potential role of HDAC1, 2, 4, and 6 in the memory of the in vivo exposure to inflammation in line with the histone code hypothesis (Jenuwein and Allis, 2001).
In light of the putative epigenetic mechanisms underlying our findings of single- and double-hit LPS signatures in microglial transcriptomes, it remains to be tested whether these signatures are indeed unique to LPS or apply more widely for perinatal exposures to other stressors, such as the psychosocial stress, e.g., caused by fear (Shapiro et al., 2013; Monteleone et al., 2014; Metz et al., 2015). Forced-swim stress applied over 4 days in adult male mice induced changes lasting at least two following weeks in neuronal acetylcholine esterase (AChE) expression via an epigenetic mechanism of hypoacetylation, with near-exclusive enrichment of HDAC4, and hypermethylation of histone H3K9 at a specific promoter of AChE (mP1c) with resulting suppression of the mE1c exon expression levels (Meshorer et al., 2002; Sailaja et al., 2012). Interestingly, a non-exclusive increase of HDAC-1, 2, and 7 was also detected. Animals showed anxiety-like behavior and this behavior as well as the AChE chromatine structure and the entire HDAC enrichment profile were reversed by NaBu, an HDAC inhibitor; the restoration of mE1c expression level was however due to HDAC4 inhibition entirely. AChE-R is the alternative splicing soluble variant of AChE-S in neurons; AChE-R production increases under various stress influences (Soreq and Seidman, 2001). This splicing switch can be induced by short-lasting (minutes) stress exposures, but can then last for weeks as shown in adult neuronal and hippocampal slice cultures (Meshorer et al., 2002; Sailaja et al., 2012). NaBu restored this splicing switch with regard to reduction in AChE-R, although the renewed increase of AChE-S variant was incomplete compared to the non-stressed animals. This finding is particularly interesting, as it sheds a new light on how stress may modulate inflammation via epigenetic mechanisms impacting the pro-inflammatory AChE. AChE inhibition restricts inflammation not only in the peripheral organs, but also in the brain (Pollak et al., 2005). The incomplete restoration of AChE-S suggests a complex regulatory network controlling AChE-S/AChE-R ratio in response to stress. Ultimately, such shifts in AChE presence in intercellular space may have long-lasting effects on cholinergic transmission with regard to cognition (cf. Section Microglial LRP-mediated Neuronal Phagocytosis May Be Enhanced By In utero Exposure to Inflammation) and neuroinflammation. Adding to the complexity of epigenetic regulation of cholinergic signaling and neuroinflammation, microRNA (miRNA)-132 has been shown in adult murine model and cell lines to potentiate cholinergic anti-inflammatory signaling in the periphery, myeloid cells in particular, and in the brain by inhibiting AChE expression (Shaked et al., 2009). The role of miRNA-132 in microglia is not yet known, but the evidence is growing for the overall importance of miRNA signaling in determining the polarization and phenotype of microglia and myeloid cells in general (Ponomarev et al., 2013).
Fear represents a model system to study chronic impact of stress on epigenome and cardiovascular system (Shenhar-Tsarfaty et al., 2015). As noted above, this approach may relate conceptually to our current findings bringing together the effects of in vivo endotoxin exposure as a stressor on the brain's microglial transcriptome and the cardiovascular system. Interestingly, changes in miRNA-608 activity on AChE binding sites in the brain (e.g., due to single nucleotide polymorphisms, SNPs) concomittantly raise levels of anxiety and blood pressure in adult mice and humans by decreasing the inhibition of AChE expression, while reducing CDC42 and IL-6 levels, important pro-inflammatory mediators (Hanin et al., 2014). This link between epigenetic signaling mechanisms, stress, and cardiovascular system is further strengthened by the recent study, in adult humans showing synergistic effects of fear as a stressor on heart rate and inflammation with cholinergic signaling playing a central role in modulating both systems (Shenhar-Tsarfaty et al., 2014, 2015). We found an increase of heart rate and a slight drop of blood pressure within the time frame of the IL-6 peak following LPS injection to the ovine fetus (Durosier et al., 2015). However, this effect appeared to dissipate at 54 h following the initial LPS exposure. Still, our experimental design does not allow drawing conclusions whether such intrauterine exposure to low-dose endotoxin concentrations may induce lasting cardiovascular changes along with alterations in innate immune responses upon repeated exposure to inflammatory stimuli. This remains subject of future studies. Interestingly, BCHE, but not ACHE, showed DE and less than one-fold down regulation (log2 = −3.197 and ) in NC vs. SHC microglia, with no detectable change in microglia exposed to LPS in vitro only or upon double-hit exposure; both BCHE and ACHE were also less than one-fold down regulated, but not differentially expressed in all other comparisons. In this regard, the potential role of serum cholinesterases as easily accessible biomarkers of neuroimmune function, along with heart rate variability monitoring, present an attractive opportunity to translate these insights into bedside applications to improve perinatal health outcomes (Durosier et al., 2015; Lake et al., 2014; Shenhar-Tsarfaty et al., 2014).
In summary, microglia pre-exposed to inflammation in vivo seem to acquire a memory of inflammation that reflects on the transcriptome by an overall decreased response in inflammatory pathways while the production of the pro-inflammatory cytokine IL-1β is up regulated. In light of the above discussion, our findings lend support to the notion of an inflammation memory in SHC sustained in SHL microglia that may be mediated by epigenetic regulatory processes involving histone acetylation and miRNA signaling. The intriguing link to the metabolic processes and cardiovascular system also deserves attention in future studies. Additional mechanistic studies (knockout, knockdown, or overexpression) are needed to validate these observations.
Microglial LRP-Mediated Neuronal Phagocytosis May be Enhanced By in Utero Exposure to Inflammation
Calreticulin (CRT) exposure on the surface of viable or apoptotic neurons is required for their phagocytosis via low-density lipoprotein receptor-related protein (LRP) receptors on LPS-stimulated primary culture rat microglia (Fricker et al., 2012). We found that the gene LRP6 is significantly up regulated after LPS exposure in vitro in NL microglia (log2 = 1.052 and ) and the activation of LRP6 is sustained in vitro in SHC microglia (log2 = 1.850 and ), i.e., after the LPS exposure in vivo. LRP1B (log2 = 6.380 and ) and LRP2 (log2 = 6.571 and ) were also strongly up regulated in SHC microglia. LRP1B and LRP2 showed a four-fold up regulation in NL microglia, however, adjacent p-values were not consistent to support this observation.
We show that a single LPS exposure in vivo or in vitro suffices to up regulate LRP genes suggesting that in utero exposure to inflammation may alter microglial—neuronal communication making CRT expressing neurons vulnerable to LRP-mediated phagocytosis. Our data does not allow validating the idea that double hit exposure to an inflammatory stimulus enhances up regulation of microglial LRP, because we could not test directly SHL vs. SHC (cf. Section Methodological Considerations, Discussion on limitations of RNAseq approach). Future studies, should estimate genetic expression profile of SHL compared to SHC.
Methodological Considerations
We could not detect any in situ neuroinflammation using Iba1, a well-established myeloid cell marker in sheep and other species. Despite the lack of overt neuroinflammation seen in situ we demonstrated a pattern of LPS-induced systemic IL-6 cytokine production in vivo and microglial IL-1β cytokine secretion in vitro. This further supports the notion that even subtle LPS exposures in utero in vivo may polarize microglia toward a neuroinflammatory phenotype without or with secondary re-exposure to an inflammatory stimulus. The LPS-triggered rise of IL-6 in plasma is in line with animal and human studies at this developmental stage (Duncombe et al., 2010; Chan et al., 2013). Microglia in vitro have been shown to secrete IL-1β preferentially when challenged with LPS, while IL-6 secretion is a hallmark of cultured astrocytes in rat (Gottschall et al., 1994). Our findings are consistent with literature and further support the cell culture purity.
In parallel to our team, the feasibility of creating a mixed primary fetal ovine brain culture has been recently, demonstrated (Weaver-Mikaere et al., 2012). We have advanced this work by focusing on late rather than mid-gestation fetuses and creating primary pure microglial culture rather than mixed culture. This allowed us to then study the microglia-specific effects of the double hit in vivo/in vitro LPS exposure on the secretion profile of the inflammatory cytokine IL-1β and the high-throughput transcriptome.
In this study, we did not discriminate between the various phenotypes of the endogenous microglia as well as the microglia recruited to the brain during the inflammatory process via the blood brain barrier, whose permeability increases under conditions of hypoxia/ischemia and fetal inflammatory response (Hutton et al., 2007, 2008; Butovsky et al., 2014; Yamasaki et al., 2014; Greter et al., 2015; Sadowska et al., 2015). Considering the mild, low-dose LPS exposure, we speculate that no recruitment of peripheral monocytes was triggered. However, we cannot state with certainty whether the microglial memory of inflammation was entirely newly established upon in vivo LPS exposure, or certain pre-existing sub-populations of microglia responded differentially to the endotoxin; another possibility needing validation remains that progenitor cells from the periphery differentiated accordingly. Hence, future studies, perhaps using single cell RNAseq, will elucidate whether the “memory” is entirely newly established, carried by a subpopulation of endogenous or periphery-recruited microglia. Isolating single cells and expanding them in culture may be another approach to test these hypotheses.
In our approach, we used DESeq2 to normalize read counts and identify differentially expressed genes. DESeq2 was specifically designed to estimate differential expression in a dataset containing replicates for both control and treatment samples. Our method used a large number of animals allowing us to have replicates for naïve control and LPS-exposed microglia. However, a limitation of our RNAseq analysis is the lack of replicates for SHC and SHL preventing us from comparing SHC to SHL directly. Other platforms meant to analyze samples without replicate could have been used here. However, we chose not to disrupt the analytical pipeline and keep the statistical analysis consistent throughout the analysis. Despite quality control measures prior to sequencing, the sample NL3 had a different expression pattern than the two other NL samples. We believe it is not related to RNA quality, and may have been due to environmental or other physiological conditions of the animal that we were not aware of at the time of the experiment. In interrogating the differential gene expression, we have ensured that the partially deviating pattern observed in sample NL3 did not confound our findings (Figures 3A,B).
Conclusions
Inflammatory microglial phenotype acquired during in vivo exposure to LPS is sustained and potentiated in vitro upon re-exposure to LPS. We identified a unique HMOX1down and FBPup phenotype of microglia exposed to the double-hit. Our results also suggest that microglia may have acquired in vivo a memory of inflammation regulated by an epigenetic process that should be confirmed by further epigenomic studies. This model allows studying mechanisms of fetal neuroinflammation in utero in vivo and in vitro to identify potential therapeutic targets for early postnatal intervention to prevent brain injury.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Acknowledgments
The authors thank Dora Siontas, Manon Blain for cell culture and ICC, Lamia Naouel Hachehouche for help with setting up the sheep specific cytokine ELISA assay, Vania Yotova for RNAseq library preparation, Jean-Christopher Grenier for alignment to the reference genome and read count, St-Hyacinthe CHUV team and M. Michel-Robinson for technical assistance. Supported by grants from the Canadian Institute of Health Research (CIHR) (MF); Fonds de la recherche en santé du Québec (FRSQ) (MF) and Molly Towell Perinatal Research Foundation (MF); QTNPR (by CIHR) (LD).
Supplementary Material
The Supplementary Material for this article can be found online at: http://journal.frontiersin.org/article/10.3389/fncel.2015.00294
References
Agrawal, V., and Hirsch, E. (2012). Intrauterine infection and preterm labor. Semin. Fetal Neonatal Med. 17, 12–19. doi: 10.1016/j.siny.2011.09.001
Ajmone-Cat, M. A., Mancini, M., De Simone, R., Cilli, P., and Minghetti, L. (2013). Microglial polarization and plasticity: evidence from organotypic hippocampal slice cultures. Glia 61, 1698–1711. doi: 10.1002/glia.22550
Anders, S., Pyl, T. P., and Huber, W. (2015). HTSeq—A Python framework to work with high-throughput sequencing data. Bioinformatics 31, 166–169. doi: 10.1093/bioinformatics/btu638
Bauer, S., Grossmann, S., Vingron, M., and Robinson, P. N. (2008). Ontologizer 2.0–a multifunctional tool for GO term enrichment analysis and data exploration. Bioinformatics 24, 1650–1651. doi: 10.1093/bioinformatics/btn250
Bickler, P. E., and Kelleher, J. A. (1992). Fructose-1,6-bisphosphate stabilizes brain intracellular calcium during hypoxia in rats. Stroke 23, 1617–1622. doi: 10.1161/01.STR.23.11.1617
Bilbo, S. D., and Schwarz, J. M. (2009). Early-life programming of later-life brain and behavior: a critical role for the immune system. Front. Behav. Neurosci. 3:14. doi: 10.3389/neuro.08.014.2009
Bilbo, S. D., and Tsang, V. (2010). Enduring consequences of maternal obesity for brain inflammation and behavior of offspring. FASEB J. 24, 2104–2115. doi: 10.1096/fj.09-144014
Billiards, S. S., Haynes, R. L., Folkerth, R. D., Trachtenberg, F. L., Liu, L. G., Volpe, J. J., et al. (2006). Development of microglia in the cerebral white matter of the human fetus and infant. J. Comp. Neurol. 497, 199–208. doi: 10.1002/cne.20991
Bolton, J. L., Auten, R. L., and Bilbo, S. D. (2014). Prenatal air pollution exposure induces sexually dimorphic fetal programming of metabolic and neuroinflammatory outcomes in adult offspring. Brain Behav. Immun. 37, 30–44. doi: 10.1016/j.bbi.2013.10.029
Butovsky, O., Jedrychowski, M. P., Moore, C. S., Cialic, R., Lanser, A. J., Gabriely, G., et al. (2014). Identification of a unique TGF-beta-dependent molecular and functional signature in microglia. Nat. Neurosci. 17, 131–143. doi: 10.1038/nn.3599
Chan, C. J., Summers, K. L., Chan, N. G., Hardy, D. B., and Richardson, B. S. (2013). Cytokines in umbilical cord blood and the impact of labor events in low-risk term pregnancies. Early Hum. Dev. 89, 1005–1010. doi: 10.1016/j.earlhumdev.2013.08.017
Chen, J., Bardes, E. E., Aronow, B. J., and Jegga, A. G. (2009). ToppGene Suite for gene list enrichment analysis and candidate gene prioritization. Nucleic Acids Res. 37, W305–W311. doi: 10.1093/nar/gkp427
Duncombe, G., Veldhuizen, R. A., Gratton, R. J., Han, V. K., and Richardson, B. S. (2010). IL-6 and TNFalpha across the umbilical circulation in term pregnancies: relationship with labour events. Early Hum. Dev. 86, 113–117. doi: 10.1016/j.earlhumdev.2010.01.027
Durafourt, B. A., Moore, C. S., Blain, M., and Antel, J. P. (2013). Isolating, culturing, and polarizing primary human adult and fetal microglia. Methods Mol. Biol. 1041, 199–211. doi: 10.1007/978-1-62703-520-0_19
Durosier, L. D., Herry, C., Cortes, M., Cao, M., Burns, P., Desrochers, A., et al. (2015). Does heart rate variability provide a signature of fetal systemic inflammatory response in a fetal sheep model of lipopolysaccharide-induced sepsis? Physiol. Meas. (in press).
Eden, E., Navon, R., Steinfeld, I., Lipson, D., and Yakhini, Z. (2009). GOrilla: a tool for discovery and visualization of enriched GO terms in ranked gene lists. BMC Bioinformatics 10:48. doi: 10.1186/1471-2105-10-48
Fahey, J. O. (2008). Clinical management of intra-amniotic infection and chorioamnionitis: a review of the literature. J. Midwifery Womens Health 53, 227–235. doi: 10.1016/j.jmwh.2008.01.001
Fahlman, C. S., Bickler, P. E., Sullivan, B., and Gregory, G. A. (2002). Activation of the neuroprotective ERK signaling pathway by fructose-1,6-bisphosphate during hypoxia involves intracellular Ca2+ and phospholipase C. Brain Res. 958, 43–51. doi: 10.1016/S0006-8993(02)03433-9
Fishman, S. G., and Gelber, S. E. (2012). Evidence for the clinical management of chorioamnionitis. Semin. Fetal Neonatal Med. 17, 46–50. doi: 10.1016/j.siny.2011.09.002
Franceschini, A., Szklarczyk, D., Frankild, S., Kuhn, M., Simonovic, M., Roth, A., et al. (2013). STRING v9.1: protein-protein interaction networks, with increased coverage and integration. Nucleic Acids Res. 41, D808–D815. doi: 10.1093/nar/gks1094
Frasch, M. G. (2014). Putative role of AMPK in fetal adaptive brain shut-down: linking metabolism and inflammation in the brain. Front. Neurol. 5:150. doi: 10.3389/fneur.2014.00150
Frasch, M. G., Müller, T., Wicher, C., Weiss, C., Löhle, M., Schwab, K., et al. (2007). Fetal body weight and the development of the control of the cardiovascular system in fetal sheep. J. Physiol. 579, 893–907. doi: 10.1113/jphysiol.2006.124800
Fricker, M., Olíva-Martin, M. J., and Brown, G. C. (2012). Primary phagocytosis of viable neurons by microglia activated with LPS or Abeta is dependent on calreticulin/LRP phagocytic signalling. J. Neuroinflammation 9, 196. doi: 10.1186/1742-2094-9-196
Gotsch, F., Romero, R., Kusanovic, J. P., Mazaki-Tovi, S., Pineles, B. L., Erez, O., et al. (2007). The fetal inflammatory response syndrome. Clin. Obstet. Gynecol. 50, 652–683. doi: 10.1097/GRF.0b013e31811ebef6
Gottschall, P. E., Tatsuno, I., and Arimura, A. (1994). Regulation of interleukin-6 (IL-6) secretion in primary cultured rat astrocytes: synergism of interleukin-1 (IL-1) and pituitary adenylate cyclase activating polypeptide (PACAP). Brain Res. 637, 197–203. doi: 10.1016/0006-8993(94)91233-5
Greter, M., Lelios, I., and Croxford, A. L. (2015). Microglia versus myeloid cell nomenclature during brain inflammation. Front. Immunol. 6:249. doi: 10.3389/fimmu.2015.00249
Hagberg, H., Peebles, D., and Mallard, C. (2002). Models of white matter injury: comparison of infectious, hypoxic-ischemic, and excitotoxic insults. Ment. Retard. Dev. Disabil. Res. Rev. 8, 30–38. doi: 10.1002/mrdd.10007
Hanin, G., Shenhar-Tsarfaty, S., Yayon, N., Yau, Y. H., Bennett, E. R., Sklan, E. H., et al. (2014). Competing targets of microRNA-608 affect anxiety and hypertension. Hum. Mol. Genet. 23, 4569–4580. doi: 10.1093/hmg/ddu170
Henkel, J. S., Beers, D. R., Zhao, W., and Appel, S. H. (2009). Microglia in ALS: the good, the bad, and the resting. J. Neuroimmune Pharmacol. 4, 389–398. doi: 10.1007/s11481-009-9171-5
Hutton, L. C., Castillo-Melendez, M., Smythe, G. A., and Walker, D. W. (2008). Microglial activation, macrophage infiltration, and evidence of cell death in the fetal brain after uteroplacental administration of lipopolysaccharide in sheep in late gestation. Am. J. Obstet. Gynecol. 198, 117.e1–117.e11. doi: 10.1016/j.ajog.2007.06.035
Hutton, L. C., Castillo-Melendez, M., and Walker, D. W. (2007). Uteroplacental inflammation results in blood brain barrier breakdown, increased activated caspase 3 and lipid peroxidation in the late gestation ovine fetal cerebellum. Dev. Neurosci. 29, 341–354. doi: 10.1159/000105475
Jenuwein, T., and Allis, C. D. (2001). Translating the histone code. Science 293, 1074–1080. doi: 10.1126/science.1063127
Kaimal, V., Bardes, E. E., Tabar, S. C., Jegga, A. G., and Aronow, B. J. (2010). ToppCluster: a multiple gene list feature analyzer for comparative enrichment clustering and network-based dissection of biological systems. Nucleic Acids Res. 38, W96–W102. doi: 10.1093/nar/gkq418
Karrow, N. A. (2006). Activation of the hypothalamic-pituitary-adrenal axis and autonomic nervous system during inflammation and altered programming of the neuroendocrine-immune axis during fetal and neonatal development: lessons learned from the model inflammagen, lipopolysaccharide. Brain Behav. Immun. 20, 144–158. doi: 10.1016/j.bbi.2005.05.003
Keogh, M. J., Bennet, L., Drury, P. P., Booth, L. C., Mathai, S., Naylor, A. S., et al. (2012). Subclinical exposure to low-dose endotoxin impairs EEG maturation in preterm fetal sheep. Am. J. Physiol. Regul. Integr. Comp. Physiol. 303, R270–R278. doi: 10.1152/ajpregu.00216.2012
Kim, D., Pertea, G., Trapnell, C., Pimentel, H., Kelley, R., and Salzberg, S. L. (2013). TopHat2: accurate alignment of transcriptomes in the presence of insertions, deletions and gene fusions. Genome Biol. 14, R36. doi: 10.1186/gb-2013-14-4-r36
Kim, Y. C., Park, T. Y., Baik, E., and Lee, S. H. (2012). Fructose-1,6-bisphosphate attenuates induction of nitric oxide synthase in microglia stimulated with lipopolysaccharide. Life Sci. 90, 365–372. doi: 10.1016/j.lfs.2011.12.011
Kuypers, E., Jellema, R. K., Ophelders, D. R., Dudink, J., Nikiforou, M., Wolfs, T. G., et al. (2013). Effects of intra-amniotic lipopolysaccharide and maternal betamethasone on brain inflammation in fetal sheep. PLoS ONE 8:e81644. doi: 10.1371/journal.pone.0081644
Lake, D. E., Fairchild, K. D., and Moorman, J. R. (2014). Complex signals bioinformatics: evaluation of heart rate characteristics monitoring as a novel risk marker for neonatal sepsis. J. Clin. Monit. Comput. 28, 329–339. doi: 10.1007/s10877-013-9530-x
Langmead, B., and Salzberg, S. L. (2012). Fast gapped-read alignment with Bowtie 2. Nat. Methods 9, 357–359. doi: 10.1038/nmeth.1923
Larouche, A., Roy, M., Kadhim, H., Tsanaclis, A. M., Fortin, D., and Sébire, G. (2005). Neuronal injuries induced by perinatal hypoxic-ischemic insults are potentiated by prenatal exposure to lipopolysaccharide: animal model for perinatally acquired encephalopathy. Dev. Neurosci. 27, 134–142. doi: 10.1159/000085985
Lin, J. W., Ju, W., Foster, K., Lee, S. H., Ahmadian, G., Wyszynski, M., et al. (2000). Distinct molecular mechanisms and divergent endocytotic pathways of AMPA receptor internalization. Nat. Neurosci. 3, 1282–1290. doi: 10.1038/81814
Livak, K. J., and Schmittgen, T. D. (2001). Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method. Methods 25, 402–408. doi: 10.1006/meth.2001.1262
Love, M. I., Huber, W., and Anders, S. (2014). Moderated estimation of fold change and dispersion for RNA-Seq data with DESeq2. Genome Biol. 15, 550. doi: 10.1186/s13059-014-0550-8
Meshorer, E., Erb, C., Gazit, R., Pavlovsky, L., Kaufer, D., Friedman, A., et al. (2002). Alternative splicing and neuritic mRNA translocation under long-term neuronal hypersensitivity. Science 295, 508–512. doi: 10.1126/science.1066752
Metz, G. A., Ng, J. W., Kovalchuk, I., and Olson, D. M. (2015). Ancestral experience as a game changer in stress vulnerability and disease outcomes. Bioessays 37, 602–611. doi: 10.1002/bies.201400217
Monteleone, M. C., Adrover, E., Pallarés, M. E., Antonelli, M. C., Frasch, A. C., and Brocco, M. A. (2014). Prenatal stress changes the glycoprotein GPM6A gene expression and induces epigenetic changes in rat offspring brain. Epigenetics 9, 152–160. doi: 10.4161/epi.25925
Murthy, V., and Kennea, N. L. (2007). Antenatal infection/inflammation and fetal tissue injury. Best Pract. Res. Clin. Obstet. Gynaecol. 21, 479–489. doi: 10.1016/j.bpobgyn.2007.01.010
Pollak, Y., Gilboa, A., Ben-Menachem, O., Ben-Hur, T., Soreq, H., and Yirmiya, R. (2005). Acetylcholinesterase inhibitors reduce brain and blood interleukin-1beta production. Ann. Neurol. 57, 741–745. doi: 10.1002/ana.20454
Ponomarev, E. D., Veremeyko, T., and Weiner, H. L. (2013). MicroRNAs are universal regulators of differentiation, activation, and polarization of microglia and macrophages in normal and diseased CNS. Glia 61, 91–103. doi: 10.1002/glia.22363
Prout, A. P., Frasch, M. G., Veldhuizen, R. A., Hammond, R., Ross, M. G., and Richardson, B. S. (2010). Systemic and cerebral inflammatory response to umbilical cord occlusions with worsening acidosis in the ovine fetus. Am. J. Obstet. Gynecol. 202, 82.e81–82.e89. doi: 10.1016/j.ajog.2009.08.020
Prout, A. P., Frasch, M. G., Veldhuizen, R., Hammond, R., Matushewski, B., and Richardson, B. S. (2012). The impact of intermittent umbilical cord occlusions on the inflammatory response in pre-term fetal sheep. PLoS ONE 7:e39043. doi: 10.1371/journal.pone.0039043
Rees, S., and Inder, T. (2005). Fetal and neonatal origins of altered brain development. Early Hum. Dev. 81, 753–761. doi: 10.1016/j.earlhumdev.2005.07.004
Rurak, D., and Bessette, N. W. (2013). Changes in fetal lamb arterial blood gas and acid-base status with advancing gestation. Am. J. Physiol. Regul. Integr. Comp. Physiol. 304, R908–R916. doi: 10.1152/ajpregu.00430.2012
Sadowska, G. B., Chen, X., Zhang, J., Lim, Y. P., Cummings, E. E., Makeyev, O., et al. (2015). Interleukin-1β transfer across the blood-brain barrier in the ovine fetus. J. Cereb. Blood Flow Metab. doi: 10.1038/jcbfm.2015.134. [Epub ahead of print].
Saigal, S., and Doyle, L. W. (2008). An overview of mortality and sequelae of preterm birth from infancy to adulthood. Lancet 371, 261–269. doi: 10.1016/S0140-6736(08)60136-1
Sailaja, B. S., Cohen-Carmon, D., Zimmerman, G., Soreq, H., and Meshorer, E. (2012). Stress-induced epigenetic transcriptional memory of acetylcholinesterase by HDAC4. Proc. Natl. Acad. Sci. U.S.A. 109, E3687–E3695. doi: 10.1073/pnas.1209990110
Shaked, I., Meerson, A., Wolf, Y., Avni, R., Greenberg, D., Gilboa-Geffen, A., et al. (2009). MicroRNA-132 potentiates cholinergic anti-inflammatory signaling by targeting acetylcholinesterase. Immunity 31, 965–973. doi: 10.1016/j.immuni.2009.09.019
Shapiro, G. D., Fraser, W. D., Frasch, M. G., and Séguin, J. R. (2013). Psychosocial stress in pregnancy and preterm birth: associations and mechanisms. J. Perinat. Med. 41, 631–645. doi: 10.1515/jpm-2012-0295
Shenhar-Tsarfaty, S., Berliner, S., Bornstein, N. M., and Soreq, H. (2014). Cholinesterases as biomarkers for parasympathetic dysfunction and inflammation-related disease. J. Mol. Neurosci. 53, 298–305. doi: 10.1007/s12031-013-0176-4
Shenhar-Tsarfaty, S., Yayon, N., Waiskopf, N., Shapira, I., Toker, S., Zaltser, D., et al. (2015). Fear and C-reactive protein cosynergize annual pulse increases in healthy adults. Proc. Natl. Acad. Sci. U.S.A. 112, E467–E471. doi: 10.1073/pnas.1418264112
Soreq, H., and Seidman, S. (2001). Acetylcholinesterase–new roles for an old actor. Nat. Rev. Neurosci. 2, 294–302. doi: 10.1038/35067589
Spencer, S. J., Auer, R. N., and Pittman, Q. J. (2006). Rat neonatal immune challenge alters adult responses to cerebral ischaemia. J. Cereb. Blood Flow Metab. 26, 456–467. doi: 10.1038/sj.jcbfm.9600206
Stansley, B., Post, J., and Hensley, K. (2012). A comparative review of cell culture systems for the study of microglial biology in Alzheimer's disease. J. Neuroinflammation 9, 115. doi: 10.1186/1742-2094-9-115
Suzuki, T., Hide, I., Matsubara, A., Hama, C., Harada, K., Miyano, K., et al. (2006). Microglial alpha7 nicotinic acetylcholine receptors drive a phospholipase C/IP3 pathway and modulate the cell activation toward a neuroprotective role. J. Neurosci. Res. 83, 1461–1470. doi: 10.1002/jnr.20850
Untergasser, A., Cutcutache, I., Koressaar, T., Ye, J., Faircloth, B. C., Remm, M., et al. (2012). Primer3–new capabilities and interfaces. Nucleic Acids Res. 40, e115. doi: 10.1093/nar/gks596
Wang, X., Rousset, C. I., Hagberg, H., and Mallard, C. (2006). Lipopolysaccharide-induced inflammation and perinatal brain injury. Semin. Fetal Neonatal Med. 11, 343–353. doi: 10.1016/j.siny.2006.04.002
Warnes, G. R., Bolker, B., Bonebakker, L., Gentleman, R., Huber, W., Liaw, A., et al. (2009). Gplots: Various R Programming Tools for Plotting Data. R package version 2.
Weaver-Mikaere, L., Gibbons, H. M., De Silva, D., and Fraser, M. (2012). Primary mixed glial cultures from fetal ovine forebrain are a valid model of inflammation-mediated white matter injury. Dev. Neurosci. 34, 30–42. doi: 10.1159/000338039
Yamasaki, R., Lu, H., Butovsky, O., Ohno, N., Rietsch, A. M., Cialic, R., et al. (2014). Differential roles of microglia and monocytes in the inflamed central nervous system. J. Exp. Med. 211, 1533–1549. doi: 10.1084/jem.20132477
Keywords: brain, neuroinflammation, bioinformatics, RNAseq, sheep, metabolism, cytokines, epigenetics
Citation: Cao M, Cortes M, Moore CS, Leong SY, Durosier LD, Burns P, Fecteau G, Desrochers A, Auer RN, Barreiro LB, Antel JP and Frasch MG (2015) Fetal microglial phenotype in vitro carries memory of prior in vivo exposure to inflammation. Front. Cell. Neurosci. 9:294. doi: 10.3389/fncel.2015.00294
Received: 03 May 2015; Accepted: 16 July 2015;
Published: 04 August 2015.
Edited by:
Carlos Barcia, Universitat Autònoma de Barcelona, SpainReviewed by:
Hermona Soreq, The Hebrew University of Jerusalem, IsraelGunnar P. H. Dietz, Schwabe Pharma Deutschland, Germany
Copyright © 2015 Cao, Cortes, Moore, Leong, Durosier, Burns, Fecteau, Desrochers, Auer, Barreiro, Antel and Frasch. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Martin G. Frasch, Départment Obstétrique-Gynécologie, CHU Ste-Justine Research Centre, Université de Montréal, 3175, Côte-Sainte-Catherine, Montréal, QC H3T 1C5, Canada, martin.frasch@recherche-ste-justine.qc.ca
†These authors have contributed equally to this work.