- Research article
- Open Access
Transcriptomics of the late gestation ovine fetal brain: modeling the co-expression of immune marker genes
BMC Genomics volume 15, Article number: 1001 (2014)
Major changes in gene expression occur in the fetal brain to modulate the function of this organ postnatally. Thus, factors can alter the genomics of the fetal brain, predisposing to neurological disorders later in life. We hypothesized that the physiological dynamics of the immune system transcriptome of the fetal brain during the last stage of gestation will reveal patterns of immune function and development in the developing brain. In this study we applied weighted gene co-expression analysis of microarrays performed on ovine fetal brain samples, to model the changes in gene expression throughout the second half of gestation.
Clusters of co-expressed genes that strongly increase in expression toward the first day of extra-uterine life are related to the hematopoietic lineage, while activation of immune pathways is induced after birth. Moreover, the pattern of gene expression suggests induction of tolerance mechanisms, probably necessary to protect highly produced proteins –such as myelin basic protein- from an autoimmune attack.
This study provides insight into the dramatic changes in gene expression that take place in the brain during the fetal life, especially during the last stage of gestation, and suggests that the immune system may have an important role in maturation of the fetal brain, which if disrupted or altered, could have negative consequences in postnatal life.
Normal brain development requires temporal changes in gene expression. In the human brain, the major changes in spatio-temporal gene expression occur during the late gestation fetal life . Little is known about development of immune function in the fetal brain, although there is a growing literature that some chronic diseases of the adult, such as hypertension , schizophrenia , and autoimmune disorders  are at least in part the result of immune dysfunction within the brain. Several studies have analyzed the effect of diverse stimuli (such as nutrition, stress, toxics and so on) on the activity of the immune system in the fetal brain and the repercussions in adulthood (reviewed in [5, 6]). Despite its importance, the ontogeny of immune cells within the fetal brain is not fully understood. The focus of studies on immune development in late gestation has been on development of these cells in somatic sites (i.e., liver, thymus, spleen, etc.). However, a better understanding of the molecular development of the immune system within the brain holds promise for design of novel therapeutics that can redirect the “programmed” brain back to a normal trajectory of development. The purpose of the present study is to model changes and coherence of gene expression in cerebral cortex, brainstem, hippocampus, and hypothalamus in order to obtain a broad view of brain ontogeny throughout the second half of gestation. We found, in addition to expected changes in pathways related to brain structure and metabolism, that there was an over-representation of genes in the hematopoietic pathway. We therefore modeled the changes in genes in the hematopoietic pathway over this period of ontogeny. We expect that the study of genomics related to immune system development in late gestation fetal brain might reveal molecular mechanisms of neuronal protection. Moreover, investigation of the genomics of fetal brain development could provide a better understanding of mechanisms underpinning postnatal neurological disorders.
We used the sheep as an animal model to identify co-expressed genes in the above regions of the ovine fetal brain, from mid-gestation (80 days, n = 4; 100 days, n = 4; 120 days, n = 4; 130 days, n = 4; 145 days, n = 4) to one day of postnatal life (n = 4). In the ewe, gestation length averages 147 days. The ovine fetus is an excellent model to study brain development since the entire gestational equivalent of human brain development occurs in utero. Altricial species (such as rats and mice) are not suitable models for human brain development because of the marked immaturity of the brain at birth in these species. Using a newly-available ovine array and weighted gene co-expression network analysis (WGCNA) on differentially expressed genes (Additional file 1: Figure S1), we identified co-expressed genes in different regions of the ovine fetal brain, from mid-gestation to one day of postnatal life. WGCNA determines pair-wise correlations between gene expression profiles to create modules (clusters) of co-expressed genes (Additional file 2: Figure S2). The resulting modules for each network were related with gestational age to determine gene significance (GS, correlation of i-th gene with the temporal pattern) and module membership (MM; correlation of the i-th gene with respect to its corresponding module) for each gene (Additional file 3: Table S1 and Additional file 4: Figure S3). Top modules were defined as those having the highest positive or negative correlation between GS and MM; i.e., modules with highest increasing or decreasing expression pattern from mid-gestation to 1 day of extra-uterine life (Additional file 5: Figure S4). These modules (or gene clusters) were selected for further analyses.
Results and discussion
Enrichment and network analysis
Significantly enriched Biological Processes (BP; p < 0.01; WebGestalt software ) of the positively correlated top modules from each brain region are shown in Additional file 6: Table S2. Negatively correlated top modules were mainly enriched with cell cycle. In agreement with other studies [1, 9], groups of co-expressed genes that increase in expression toward the end of gestation are related with response to stress and stimulus, myelination, apoptosis, generation of energy, and angiogenesis, while genes that decrease in expression are mainly related to cell cycle.
These events occur because during the last steps of fetal brain development there is a deceleration in cell division (quiescence) simultaneous with an acceleration in neuronal developmental and differentiation [9–11].
Significantly enriched KEGG pathways (p < 0.01; WebGestalt software ) were compared to identify common enriched pathways between all brain regions (Table 1). Some of the common enriched pathways in the positively correlated top modules were expected in the developing fetal brain, such as those related with metabolism. Genes related with endocytosis and lysosomal activity were also expected to be expressed at increasing levels by the end of gestation, since they are particularly relevant in polarized cells like the neurons . Notably, a common enriched KEGG pathway related with the immune function was the hematopoietic cell lineage. We focused further analysis on genes related with the hematopoietic lineage. These genes were visualized as separate networks, according to the brain region, using the Cytoscape software (3.0)  and physical and genetic interactions between genes were inferred by use of the GeneMania plugin  (Figure 1). Then, these networks were merged together to obtain a single network connecting the nodes corresponding to the genes of interest from the four brain regions. The nodes (genes) that comprise this new network represent markers for different immune cells and at different stages of maturation. Therefore, the whole network can be correlated with the hematopoietic lineage (Figure 2), as was indicated by the “hematopoietic cell lineage” as a common enriched KEGG pathway.
These results do not necessarily mean that the brain is a hematopoietic organ similar to liver or spleen during fetal life. Nor can it be concluded from this study that hematopoietic stem cells are present in the ovine fetal brain. Rather, these results strongly suggest the presence of developing myeloid and lymphoid derived cells in different regions of the fetal brain (vide infra).
Detailed analysis of immune marker genes
Marker genes for progenitor immune cells were found in the top positive correlated module of the ovine fetal brainstem. These markers are CD34, CD109 CD44, CD5 and CD9 (Additional file 7: Figure S5). CD109 antigen is expressed by a subpopulation of CD34+ Hematopoietic stem cell (HSC) and progenitor cells in the fetal and adult bone marrow . An isoform of CD44 is expressed on CD34+ cells: HSC and progenitor cells . Progenitors cells of the myeloid lineage, such as colony-forming unit macrophage/dendritic cells, express the colony stimulating factor 1 receptor (CSF1R) for differentiation. The resident myeloid inflammatory cells of the central nervous system (CNS) parenchyma are microglial cells. Thus, expression of CSF1R seems critical for microglia development. In mice, the absence of CSF1R (Csf1r -/-) results in essentially no macrophages/microglia (<99% depletion) in the embryonic and early post-natal brain . In mice, interleukin 34 (IL34) binds to the CSF1R receptor with high affinity to regulate myeloid development and it can substitute for colony stimulating factor 1 (CSF-1) in vivo. IL34 is more highly expressed than CSF-1 mRNA in most of the regions of the developing brain . These genes were part of the top positive correlated module in the cortex network, as well as other genes markers of antigen presenting cells, including CD1D, integrin, alpha M (CD11B) and CD81, (Additional file 8: Figure S6) indicating the development of microglia in the fetal CNS parenchyma.Interestingly, genes related to hematopoiesis increased in expression during fetal life, but genes related to immune system activation and antigen presentation are strongly induced after birth (Figure 3A). These genes are those encoding for the mayor histocompatibility complex class I and II; such as HLA-B and B2M (class I) or CD74, HLA-DMA and HLA-DRB1 (class II), or genes encoding for components of the immunoproteosome necessary to process class I MHC peptides, such as PSMB8 and PSMB10.
Therefore, these results suggest a regulated development of the immune system in the fetal brain. In agreement with this concept, mRNA expression of genes markers for mature T-cells including the three chains of the CD3 marker (CD3G-CD3E-CD3D) were strongly increased after birth in the ovine fetal hippocampus, when measured by RT-PCR (Additional file 9: Figure S7), rising the possibility that T cells are present in the fetal brain.
Circulating lymphoid cells could trespass the relatively permeable blood brain barrier (BBB) during fetal life and reach the brain parenchyma . If any of these cells escape central tolerance, it is possible for them to react against highly antigenic self-proteins which are produced at increasing rate in late gestation, such as the highly antigenic myelin basic protein (MBP) causing the autoimmune demyelinating disease of the CNS known as multiple sclerosis or in the animal model, experimental autoimmune encephalitis (EAE) .
Expression of MBP gene was strongly positively correlated with fetal age for all brain regions. Correlations ranged from 0.87 in hippocampus to 0.92 in cortex, with an average of 0.893. The gene encoding for MBP was one of the top positive correlated modules for all brain regions but also showed a strong increasing mRNA expression measured by RT-PCR in the ovine fetal cortex (Additional file 10: Figure S8A). Consequently, resident cells of the CNS should be able to counteract the action of potential auto-reactive lymphocytes.
Remarkably, an immune marker gene that was part of the top negatively correlated modules for all brain regions was CD24. Correlations of this gene with fetal age ranged from -0.83 in hippocampus to -0.98 in cortex (average -0.895). Thus, the expression trajectory for this gene was the opposite of MBP (Figure 3B). mRNA expression of this gene, as measured by RT-PCR in ovine fetal cortex was dramatically decreased with fetal age (Additional file 10: Figure S8B). CD24 is a surface glycoprotein that is strongly expressed by immature B- and T- lymphocytes, but which disappears from T cell surface as T cell mature. However, activation of T cells resulted in a rapid induction of CD24 expression . CD24 is also expressed on resident cells of the CNS  and it is able to inhibit the proliferation of immature neurons .
Expression of CD24 in both T cells and non-T host cells is required for the development of EAE . Expression of CD24 in T cells actively inhibits clonal deletion (necessary for central tolerance) allowing the expansion of autoreactive T cells . Expression of CD24 in the CNS resident cells, including microglial cells and astrocytes, costimulates myelin antigen specific T cells and promotes their effector function, promoting the autoimmune inflammation in the CNS . Taken together, these data imply that CD24 is involved in the autoimmune reaction to myelin antigens.
Microglia cells are able to acquire a dendritic cell phenotype and mediate peripheral tolerance to avoid progression of EAE . Moreover, protection from EAE development, mediated by dendritic cells, was exclusively associated with low expression of CD24, when compared to EAE controls . Thus, we speculated that the strong downregulation of CD24 with advancing fetal age results in a tolerance mechanism necessary to avoid an autoimmune reaction to myelin antigens, such as MBP, which could be related to the susceptibility to this type of disease later in life.
Accordingly, genes related with immune tolerance followed a pattern of increasing expression in late gestation. The Fc fragment of IgG, low affinity IIb, receptor (FCGRIIB) for example, is an inhibitory receptor on myeloid cells that can regulate T cell tolerance and promotes T-regulatory cell induction by increasing interleukin 10 (IL10) production [29, 30]. Both IL10 and transforming growth factor, beta 1 (TGFb), another cytokine involved in T-regulatory cell induction, showed an increased expression in the cortex with fetal age as measured by microarray and qRT-PCR (Additional file 10: Figure S8C).
This study provides insight into the dramatic changes in gene expression that take place in the brain during the fetal life, especially during the last stage of gestation.
Results from the present study propose a regulated development of the immune system in the ovine fetal brain. This concept implies induction of tolerance mechanisms, to protect fetal myelin from a possible autoimmune attack, and activation of the immune response after birth, which could be part of the physiological adaptation of the newborn for the extrauterine life. These findings could provide the basis for new investigation in this area, in order to understand key regulators of fetal brain development.
Tissues were collected from fetuses at 80 (80 d, n = 4), 96–100 (100 d, n = 4), 120 (120 d, n = 4), 130 (130 d, n = 4), and 142–144 (145 d, n = 4) days of gestation and on the first (1 d, n = 4) day after delivery. Each group included one set of twin fetuses. None of the ewes showed any signs of impending labor. For collection of fetal tissues, ewes were killed with 20 ml of Euthasol solution (7.8 g pentobarbital and 1 g phenytoin sodium; Virbac AH, Fort Worth, TX) administered intravenously, the fetus was quickly removed, and the fetal brain was removed. Brain dissection was performed immediately, using sterile instruments and gloves. First, the head was separated from the backbone cutting cranial to the axis of the vertebrae. Second, the posterior cranial bones were separated using forceps. With a scalpel we cut the pituitary stalk and optic nerves to release the brain, which was then deposited onto a sterile container or surface. Holding the brain upside down, the hypothalamus was dissected, followed by both hippocampi and the brainstem. Finally, a section of the frontal cortex was dissected. Fetal tissues were rapidly frozen in liquid nitrogen and stored at -80°C.
The use of animals in this project was approved by the University of Florida Institutional Animal Care and Use Committee.
RNA extraction and preparation
RNA was extracted from the cortex, hippocampus, hypothalamus and brainstem using Trizol (Invitrogen, Carlsbad, CA) and following the manufacturer’s directions. The RNA was resuspended in RNAsecure, and stored at -80°C in aliquots until use. For microarray analysis, 20 ug of RNA was DNase treated using the Turbo RNase-free DNase kit (Ambion, Foster City, CA), the concentration determined with a Nanodrop spectrophotometer (ND-1000, ThermoFisher, Wilmington DE) and the integrity of the RNA was measured using an Agilent Bioanalyzer, 2100 model. One μg of the DNase-treated RNA was labeled with Cyanine 3 (Cy3) CTP with the Agilent Quick Amp kit (5190–0442, New Castle, DE) according to their methodology, purified with the Qiagen RNeasy kit (Valencia, CA) as according to Agilent’s revision of the Qiagen protocol as shown in the Quick Amp kit protocol except that the microcentrifugation spins were performed at room temperature instead of 4°C. The resulting labeled cRNA was analyzed with the Nanodrop spectrophotometer, and the specific activities and the yields of the cRNAs were calculated. The labeled cRNA was stored at -80°C until use.
This was performed following protocols from Agilent. Briefly, 600 ng of each labeled cRNA was fragmented and then mixed with hybridization buffer using the Agilent gene expression hybridization kit. These were applied to sheep 8 × 15 K array slides (Agilent 019921), containing 8 arrays with 15,208 oligomers with a length of 60 bases and hybridized at 65°C for 17 h at 10 rpm. A total of 3 slides (24 arrays) were employed per each brain region (six gestational ages times four replicates per gestational age). The arrays were washed, dried, stabilized, and scanned with an Agilent G2505B 2 dye scanner at the Interdisciplinary Center for Biotechnology Research at the University of Florida. Features were extracted with Agilent Feature extraction 9.1 software.
The limma package was employed to import the raw data into R (http://www.r-project.org), perform background correction and normalize the data using the quantile normalization method . Control probes and low expressed probes were filtered out, retaining for further analysis the probes that were at least 10% brighter than the negative controls on at least four arrays.
The Bayesian Estimation of Temporal Regulation (BETR) algorithm was used to identify the differentially expressed genes (DEG) for each brain region at a False Discovery Rate (FDR) < 0.05. The first gestational age (80 days) was considered as baseline measurement and compared to the subsequent time points, to correlate the differential expression between time points. The algorithm was applied using the BETR package for R software. A detailed explanation of the mathematical model can be found in Aryee et al. . This method returns the probabilities of differential expression for each gene in the data set. Genes with a probability higher than 99.99% were considered as DEG.
Supervised weighted gene co-expression network analysis (WGCNA)
The DEG for each brain region was subjected to signed-WGCNA. Rows of each data set were collapsed forming an average expression of the genes with the same official symbol, in order to obtain unique identifiers for each gene in the working data set. The automatic method was employed for block-wise network construction and module detection. The co-expression similarity was raised to a soft thresholding power (β) of 12 to calculate adjacency. The adjacency for the signed network is defined as a ij = |(1 + cor(x i ,x j ))/2|β. The resulting modules for each network were related with the gestational age to identify modules, or clusters, of co-expressed genes with increasing or decreasing expression pattern. Gene significance (GS) was defined as the correlation of i-th gene with a temporal pattern. Module membership (MM) was defined as the correlation of the i-th gene respects its corresponding module (the higher is the MM the more connected is the i-th gene with the other genes of the corresponding modules). The correlation coefficient of MM and GS was measure for each module, plotting MM versus GS. Higher correlation between MM and GS means that genes that are highly associated with the temporal pattern are also the central elements of the given module  (Additional file 3: Table S1). Modules that showed the highest correlation between MM and GS, either with increasing pattern (positive GS) or decreasing pattern (negative GS), were selected for Enrichment analysis. The top positively correlated module for each brain region was: the turquoise module for cortex and brainstem, the brown module for hippocampus and the blue module for hypothalamus. The top negatively correlated module for each brain region was: the blue module for cortex and brainstem, the yellow module for hippocampus and the turquoise module for hypothalamus (Additional file 4: Figure S3).
All the analyses were performed with the WGCNA package for R software. More details about the methodology for WGNA for network construction can be found at: http://labs.genetics.ucla.edu/horvath/CoexpressionNetwork/Rpackages/WGCNA/.
The bioinformatics tool applied for this purpose was WebGestalt (WEB-based GEne SeT AnaLysis Toolkit). This program is designed for functional genomic studies that generate large number of gene list. Then, it organizes the large set of genes based on common functional features, like GO categories or biochemical pathways . The list of human official symbols for the genes composing the top modules for each network was submitted for enrichment biological processes and KEGG pathway analysis, selecting H. Sapiens as the organism of interest. The statistical method employed was the hypergeometric test, adjusting the p-values for Benjamini & Hochberg method, selecting the minimun of 5 genes per category and using the human genome as reference set. Significantly enriched Kegg pathways (p < 0.01) were compared to identify common enriched pathways between all brain regions.
A gene network was constructed for each brain region using CytoScape version 3.0  through the GeneMania plugin , which was used to infer network data. The set of functional association data between genes was downloaded from the Homo sapiens database. The list of human official symbols for the genes of interest was input into the GeneMania plugin to retrieve the corresponding association network. The association data employed was protein-protein and protein-DNA interaction. Then, the resulting networks were merged together using the Merge Networks tool of Cytoscape 3.0.
Quantitative Real-Time (qRT)-PCR Validation
The mRNA samples extracted from the four brain regions at the six different gestational ages (4 fetuses/gestational age; 96 samples in total) were converted to cDNA with a High Capacity cDNA Archive kit using the methodology recommended by the kit manufacturer (Applied Biosystems, Foster City, Calif., USA). The newly synthesized cDNA was stored at -20°C until qRT-PCR was performed.
A total of 18 genes were selected for further validation by qRT-PCR. The cDNA from corresponding brain region was used to measure the mRNA of the following genes: Brainstem: CD34, CD109, CD44, CD5 and CD9; Hippocampus: CD3 gamma (CD3G), CD3 delta (C3D) and CD3 epsilon (CD3E); Cortex: CSF-1 CSF1R, IL34, CD11B, CD81, FCGRIIB, IL10, TGFB, CD24 and MBP.
Relative expression of selected genes were determined using primers (Sigma-Aldrich, St Louis, MO) and Sybr Green PCR Master Mix (Applied Biosystems, Foster City, CA). Primers were designed with Primer Express software (Applied Biosystems. Primers for CD109, CD44, CSF1, CSF1R, IL34, CD24 and MBP were designed from the corresponding bovine mRNA. Primers for CD34, CD5, CD9, CD3G, CD3D, CD3E, CD81, CD11B, FCGR2B, IL10 and TGFB were designed from the ovine mRNA. Primers sequences and accession numbers are reported in Additional file 11: Table S3. All primer pairs had efficiencies greater than 95%. The abundance of B-actin mRNA was determined in each brain region sample, using primers and VIC Taqman probes designed from the ovine B-actin sequence and Taqman qRT-PCR master mix (Applied Biosystems, Foster City, CA). All samples were run in triplicate for each gene and for B-actin. Relative mRNA expression of each gene was calculated by determining change in threshold cycle (ΔCt) between the mean Ct for each gene and the mean Ct for B-actin mRNA from the same sample. The effect of gestational age on each gene was analyzed by ANOVA using the ΔCt values. Pairwise comparisons among gestational age means were done by Duncan’s procedure. Data were graphed as the mean fold change in mRNA relative to the 80d group; fold change in each sample was calculated as 2-ΔΔCt, where ΔΔCt is the difference between ΔCt in each sample and the mean ΔCt in the control group. For all statistical analyses, the criterion for achieving statistical significance was P < 0.05.
Kang HJ, Kawasawa YI, Cheng F, Zhu Y, Xu X, Li M, Sousa AM, Pletikos M, Meyer KA, Sedmak G, Guennel T, Shin Y, Johnson MB, Krsnik Z, Mayer S, Fertuzinhos S, Umlauf S, Lisgo SN, Vortmeyer A, Weinberger DR, Mane S, Hyde TM, Huttner A, Reimers M, Kleinman JE, Sestan N: Spatio-temporal transcriptome of the human brain. Nature. 2011, 478: 483-489. 10.1038/nature10523.
Schiffrin EL: T lymphocytes: a role in hypertension?. Curr Opin Nephrol Hypertens. 2010, 19: 181-186. 10.1097/MNH.0b013e3283360a2e.
Brown AS: Prenatal infection as a risk factor for schizophrenia. Schizophr Bull. 2006, 32: 200-202. 10.1093/schbul/sbj052.
Dammann O, Hagberg H, Leviton A: Is periventricular leukomalacia an axonopathy as well as an oligopathy?. Pediatr Res. 2001, 49: 453-457. 10.1203/00006450-200104000-00003.
Bellinger DL, Lubahn C, Lorton D: Maternal and early life stress effects on immune function: relevance to immunotoxicology. J Immunotoxicol. 2008, 5: 419-444. 10.1080/15476910802483415.
Marques AH, O’Connor TG, Roth C, Susser E, Bjorke-Monsen AL: The influence of maternal prenatal and early childhood nutrition and maternal prenatal stress on offspring immune system development and neurodevelopmental disorders. Front Neurosci. 2013, 7: 120-
Ramadoss J, Lunde ER, Chen WJ, West JR, Cudd TA: Temporal vulnerability of fetal cerebellar Purkinje cells to chronic binge alcohol exposure: ovine model. Alcohol Clin Exp Res. 2007, 31: 1738-1745. 10.1111/j.1530-0277.2007.00477.x.
Zhang B, Kirov S, Snoddy J: WebGestalt: an integrated system for exploring gene sets in various biological contexts. Nucleic Acids Res. 2005, 33: W741-W748. 10.1093/nar/gki475.
Colantuoni C, Lipska BK, Ye T, Hyde TM, Tao R, Leek JT, Colantuoni EA, Elkahloun AG, Herman MM, Weinberger DR, Kleinman JE: Temporal dynamics and genetic control of transcription in the human prefrontal cortex. Nature. 2011, 478: 519-523. 10.1038/nature10524.
Johnson MB, Kawasawa YI, Mason CE, Krsnik Z, Coppola G, Bogdanovic D, Geschwind DH, Mane SM, State MW, Sestan N: Functional and evolutionary insights into human brain development through global transcriptome analysis. Neuron. 2009, 62: 494-509. 10.1016/j.neuron.2009.03.027.
Oldham MC, Konopka G, Iwamoto K, Langfelder P, Kato T, Horvath S, Geschwind DH: Functional organization of the transcriptome in human brain. Nat Neurosci. 2008, 11: 1271-1282. 10.1038/nn.2207.
Marin MP, Esteban-Pretel G, Ponsoda X, Romero AM, Ballestin R, Lopez C, Megias L, Timoneda J, Molowny A, Canales JJ, Renau-Piqueras J: Endocytosis in cultured neurons is altered by chronic alcohol exposure. Toxicol Sci. 2010, 115: 202-213. 10.1093/toxsci/kfq040.
Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, Amin N, Schwikowski B, Ideker T: Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 2003, 13: 2498-2504. 10.1101/gr.1239303.
Montojo J, Zuberi K, Rodriguez H, Kazi F, Wright G, Donaldson SL, Morris Q, Bader GD: GeneMANIA cytoscape plugin: fast gene function predictions on the desktop. Bioinformatics. 2010, 26: 2927-2928. 10.1093/bioinformatics/btq562.
Murray LJ, Bruno E, Uchida N, Hoffman R, Nayar R, Yeo EL, Schuh AC, Sutherland DR: CD109 is expressed on a subpopulation of CD34+ cells enriched in hematopoietic stem and progenitor cells. Exp Hematol. 1999, 27: 1282-1294. 10.1016/S0301-472X(99)00071-5.
Neu S, Geiselhart A, Sproll M, Hahn D, Kuci S, Niethammer D, Handgretinger R: Expression of CD44 isoforms by highly enriched CD34-positive cells in cord blood, bone marrow and leukaphereses. Bone Marrow Transplant. 1997, 20: 593-598. 10.1038/sj.bmt.1700940.
Erblich B, Zhu L, Etgen AM, Dobrenis K, Pollard JW: Absence of colony stimulation factor-1 receptor results in loss of microglia, disrupted brain development and olfactory deficits. PLoS One. 2011, 6: e26317-10.1371/journal.pone.0026317.
Wei S, Nandi S, Chitu V, Yeung YG, Yu W, Huang M, Williams LT, Lin H, Stanley ER: Functional overlap but differential expression of CSF-1 and IL-34 in their CSF-1 receptor-mediated regulation of myeloid cells. J Leukoc Biol. 2010, 88: 495-505. 10.1189/jlb.1209822.
Stonestreet BS, Patlak CS, Pettigrew KD, Reilly CB, Cserr HF: Ontogeny of blood–brain barrier function in ovine fetuses, lambs, and adults. Am J Physiol. 1996, 271: R1594-R1601.
Matyszak MK, Perry VH: The potential role of dendritic cells in immune-mediated inflammatory diseases in the central nervous system. Neuroscience. 1996, 74: 599-608. 10.1016/0306-4522(96)00160-1.
Hubbe M, Altevogt P: Heat-stable antigen/CD24 on mouse T lymphocytes: evidence for a costimulatory function. Eur J Immunol. 1994, 24: 731-737. 10.1002/eji.1830240336.
Bai X-F, Li O, Zhou Q, Zhang H, Joshi PS, Zheng X, Liu Y, Wang Y, Zheng P, Liu Y: CD24 controls expansion and persistence of autoreactive T cells in the central nervous system during experimental autoimmune encephalomyelitis. J Exp Med. 2004, 200: 447-458. 10.1084/jem.20040131.
Belvindrah R, Rougon G, Chazal G: Increased neurogenesis in adult mCD24-deficient mice. J Neurosci. 2002, 22: 3594-3607.
Kuchroo VK, Anderson AC, Waldner H, Munder M, Bettelli E, Nicholson LB: T cell response in experimental autoimmune encephalomyelitis (EAE): role of self and cross-reactive antigens in shaping, tuning, and regulating the autopathogenic T cell repertoire. Annu Rev Immunol. 2002, 20: 101-123. 10.1146/annurev.immunol.20.081701.141316.
Carl JW, Liu J-Q, Joshi PS, El-Omrani HY, Yin L, Zheng X, Whitacre CC, Liu Y, Bai X-F: Autoreactive T cells escape clonal deletion in the thymus by a CD24-dependent pathway. J Immunol. 2008, 181: 320-328. 10.4049/jimmunol.181.1.320.
Liu J-Q, Carl JW, Joshi PS, RayChaudhury A, Pu X-A, Shi F-D, Bai X-F: CD24 on the resident cells of the central nervous system enhances experimental autoimmune encephalomyelitis. J Immunol. 2007, 178: 6227-6235. 10.4049/jimmunol.178.10.6227.
Almolda B, Gonzalez B, Castellano B: Activated microglial cells acquire an immature dendritic cell phenotype and may terminate the immune response in an acute model of EAE. J Neuroimmunol. 2010, 223: 39-54. 10.1016/j.jneuroim.2010.03.021.
Liu X, Lindberg R, Xiao BG, Steffensen KR, Leppert D, Link H, Huang YM: CD24 and myosin light polypeptide 2 are involved in prevention of experimental autoimmune encephalomyelitis by myelin basic protein-pulsed dendritic cells. J Neuroimmunol. 2006, 172: 137-144. 10.1016/j.jneuroim.2005.11.013.
Desai DD, Harbers SO, Flores M, Colonna L, Downie MP, Bergtold A, Jung S, Clynes R: Fc gamma receptor IIB on dendritic cells enforces peripheral tolerance by inhibiting effector T cell responses. J Immunol. 2007, 178: 6217-6226. 10.4049/jimmunol.178.10.6217.
Samsom JN, van Berkel LA, van Helvoort JM, Unger WW, Jansen W, Thepen T, Mebius RE, Verbeek SS, Kraal G: Fc gamma RIIB regulates nasal and oral tolerance: a role for dendritic cells. J Immunol. 2005, 174: 5279-5287. 10.4049/jimmunol.174.9.5279.
Smyth GK: Limma: Linear Models for Microarray Data. Bioinformatics and Computational Biology Solutions Using R and Bioconductor, R. Edited by: Gentleman R, Carey V, Dudoit S, Irizarry R, Huber W. 2005, New York: Springer, 397-420.
Aryee MJ, Gutierrez-Pabello JA, Kramnik I, Maiti T, Quackenbush J: An improved empirical bayes approach to estimating differential gene expression in microarray time-course data: BETR (bayesian estimation of temporal regulation). BMC Bioinformatics. 2009, 10: 409-10.1186/1471-2105-10-409.
Zhang B, Horvath S: A general framework for weighted gene co-expression network analysis. Stat Appl Genet Mol Biol. 2005, 4: Article17-
We thank Drs. Nancy Denslow and Elaine Richards for help with the microarray chemistry and for their outstanding advice. We thank Dr. Henry Baker for his advice and consultation. We thank the University of Florida Interdisciplinary Center for Biotechnology Research for technical support. This work was supported by NIH grants HD057561 and HD033053 (to CEW), HD057581 (to MKW) and ARRA supplement to HD057561.
Microarray data have been deposited in the NCBI Gene Expression Omnibus under accession number GSE50460 (http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE50460).
The authors declare that they have no competing interests.
MBR conceived of, designed, and applied the method of gene co-expression analysis to gene array data from fetal brain regions. MBR also provided the basis for this approach by annotating the Agilent 15.5 k ovine array. CEW and MKW collected the tissue and hypothesized that patterns of gene expression will reveal processes critical to the development of the late-gestation fetus and of the biological development that governs readiness for birth. All authors read and approved the final manuscript.
Electronic supplementary material
Additional file 1: Figure S1: Venn Diagram of the number of differentially expressed genes (DEG) in the ovine fetal brain. DEG were selected as those following a temporal profile during the last stage of gestation in different regions of the ovine fetal brain. (PDF 51 KB)
Additional file 2: Figure S2: Weighted gene co-expression network analysis (WGCNA) of the ovine fetal brain. WGCNA identifies modules of co-expressed genes following an increasing or decreasing expression temporal pattern during last stage of gestation in ovine fetal cortex (A), brainstem (B), hippocampus (C) and hypothalamus (D). Modules are colored according the number of genes contained in the module. For example, the module’s color is turquoise if contains the largest number of genes. (PDF 205 KB)
Additional file 3: Table S1: Modules attributes for each brain region. Modules were determined by weighted gene-coexpression analysis on each brain region of the ovine fetal brain. Gene significance (GS) with the temporal pattern (Gestational Age) and module membership (MM) on each module (and respective p-value) was determined for each differentially expressed gene. (XLSX 4 MB)
Additional file 4: Figure S3: Relationship between modules identified in weighted gene co-expression networks and fetal age. The modules were determined in fetal cortex (A), brainstem (B), hippocampus (C) and hypothalamus (D). Positive correlations are showed in red and negative correlations in blue. The stronger the color the higher the correlation index. (PDF 146 KB)
Additional file 5: Figure S4: Module membership versus gene significance plots. Modules composed by highly connected genes with the highest positive or negative correlation with gestational age and the respective dendrogram, identified en network analysis from the ovine fetal cortex (A), brainstem (B), hippocampus (C) and hypothalamus (D). Each plot is colored according the corresponding module. (PDF 956 KB)
Additional file 6: Table S2: Enriched biological processes (BP) in the fetal brain during the last stage of gestation. These BP were significantly enriched with the co-expressed genes forming part of the top modules with highest positive correlation with fetal age in each of the ovine fetal brain regions (i.e., strongly increasing in expression toward the first day of extra-uterine life). (XLSX 12 KB)
Additional file 7: Figure S5: qRT-PCR validation for CD34, CD109, CD44, CD5 and CD9. Gene expression of CD34, CD109, CD44, CD5 and CD9 measured by microarray at 80, 100, 120, 130, 145 days of gestation and 1 day of extra-uterine life and corresponding fold changes in mRNA concentration relative to 80 days, measured by qRT- PCR in samples from ovine fetal brainsteam. Data are fold differences relative to mean expression at 80d. a - different from 80d values; b - different from 100d values; c - different from 120d values; d - different from 130d values; e - different from 145d values; f – different from 1d of extra-uterine life values. For all statistical comparisons, P < 0.05 was used as the criterion for significance. (PDF 402 KB)
Additional file 8: Figure S6: qRT-PCR validation for CSF1R, CSF1, IL34, CD11b and CD81. Gene expression of CSF1R, CSF1, IL34, CD11b and CD81 measured by microarray at 80, 100, 120, 130, 145 days of gestation and 1 day of extra-uterine life and corresponding fold changes in mRNA concentration relative to 80 days, measured by qRT- PCR in samples from ovine fetal cortex. Data are fold differences relative to mean expression at 80d. a - different from 80d values; b - different from 100d values; c - different from 120d values; d - different from 130d values; e - different from 145d values. For all statistical comparisons, P < 0.05 was used as the criterion for significance. (PDF 248 KB)
Additional file 9: Figure S7: qRT-PCR validation for CD3G, CD3D and CD3E. Gene expression of CD3G, CD3D and CD3E measured by microarray at 80, 100, 120, 130, 145 days of gestation and 1 day of extra-uterine life and corresponding fold changes in mRNA concentration relative to 80 days, measured by qRT- PCR in samples from ovine fetal hippocampus. Data are fold differences relative to mean expression at 80d. a - different from 80d values; b - different from 100d values; c - different from 120d values; d - different from 130d values; e - different from 145d values. For all statistical comparisons, P < 0.05 was used as the criterion for significance. (PDF 214 KB)
Additional file 10: Figure S8: qRT-PCR validation for MBP, CD24, FGRIIB, IL10, and TGFB. Gene expression of MBP (A), CD24 (B) and FGRIIB, IL10, TGFB (C) measured by microarray at 80, 100, 120, 130, 145 days of gestation and 1 day of extra-uterine life and corresponding fold changes in mRNA concentration relative to 80 days, measured by qRT- PCR in samples from ovine fetal cortex. Data are fold differences relative to mean expression at 80d. a - different from 80d values; b - different from 100d values; c - different from 120d values; d - different from 130d values; e - different from 145d values. For all statistical comparisons, P < 0.05 was used as the criterion for significance. (PDF 271 KB)
About this article
Cite this article
Rabaglino, M.B., Keller-Wood, M. & Wood, C.E. Transcriptomics of the late gestation ovine fetal brain: modeling the co-expression of immune marker genes. BMC Genomics 15, 1001 (2014). https://doi.org/10.1186/1471-2164-15-1001
- Myelin Basic Protein
- Fetal Brain
- Differentially Express Gene
- Experimental Autoimmune Encephalitis
- Module Membership