Volume 5 Supplement 4
Gene expression pattern in swine neutrophils after lipopolysaccharide exposure: a time course comparison
© Sanz-Santos et al; licensee BioMed Central Ltd. 2011
Published: 3 June 2011
Experimental exposure of swine neutrophils to bacterial lipopolysaccharide (LPS) represents a model to study the innate immune response during bacterial infection. Neutrophils can effectively limit the infection by secreting lipid mediators, antimicrobial molecules and a combination of reactive oxygen species (ROS) without new synthesis of proteins. However, it is known that neutrophils can modify the gene expression after LPS exposure. We performed microarray gene expression analysis in order to elucidate the less known transcriptional response of neutrophils during infection.
Blood samples were collected from four healthy Iberian pigs and neutrophils were isolated and incubated during 6, 9 and 18 hrs in presence or absence of lipopolysaccharide (LPS) from Salmonella enterica serovar Typhimurium. RNA was isolated and hybridized to Affymetrix Porcine GeneChip®. Microarray data were normalized using Robust Microarray Analysis (RMA) and then, differential expression was obtained by an analysis of variance (ANOVA).
ANOVA data analysis showed that the number of differentially expressed genes (DEG) after LPS treatment vary with time. The highest transcriptional response occurred at 9 hr post LPS stimulation with 1494 DEG whereas at 6 and 18 hr showed 125 and 108 DEG, respectively. Three different gene expression tendencies were observed: genes in cluster 1 showed a tendency toward up-regulation; cluster 2 genes showing a tendency for down-regulation at 9 hr; and cluster 3 genes were up-regulated at 9 hr post LPS stimulation. Ingenuity Pathway Analysis revealed a delay of neutrophil apoptosis at 9 hr. Many genes controlling biological functions were altered with time including those controlling metabolism and cell organization, ubiquitination, adhesion, movement or inflammatory response.
LPS stimulation alters the transcriptional pattern in neutrophils and the present results show that the robust transcriptional potential of neutrophils under infection conditions, indicating that active regulation of gene expression plays a major role in the neutrophil-mediated- innate immune response.
Neutrophils play a key role in innate immune response. They initiate phagocytosis, degranulation and killing without new synthesis of proteins. However, it has been demonstrated that new gene transcription and protein synthesis are required to maintain full capacity for human neutrophil phagocytosis and associated bactericidal activity [1, 2].
LPS treatment enhances neutrophil bactericidal activity, with an alteration in adhesion, respiratory burst, degranulation and motility [3, 4]. Thus kinetic study of swine neutrophil response to LPS represents an in vitro model to investigate the innate immune response during bacterial infection.
To test the neutrophil transcriptional potential, global gene expression analysis was performed and the results indicated that the LPS-treated neutrophils increase their transcriptional activity by altering genes involved in different cellular processes including transcriptional regulation, cell signalling, cytoskeletal reorganization, etc. Furthermore, inhibition of apoptosis-related functions may indicate an increase of the neutrophil life span at 9 hr post LPS stimulation. Taken together, our results show that (1) relevant changes occurred in the pattern of gene expression of porcine neutrophils after LPS exposure, and (2) the biological functions identified to be meaningful using Ingenuity Pathway Analysis tool (IPA) suggests that these transcriptional changes play an important role in inflammation and immune response associated with LPS stimulation.
Isolation and stimulation of neutrophils
Blood samples were collected from 4 healthy Iberian pigs at the local abattoir. Neutrophils were isolated with dextran and ficoll sedimentation as previously described  and adapted for swine neutrophils. Isolated neutrophils were more than 99% pure and < 1% monocytes were present. Freshly purified neutrophils, ~200 x 106 for each time point, were suspended at 5 x 106 ml-1 in 40 ml RPMI medium (Lonza, Basel, Switzerland) supplemented with 10% heat-inactivated bovine fetal serum. Stimulation was performed in 50 ml Falcon tubes and rotated continuously at 37°C for 6, 9 and 18 hr in presence or absence of 100 ng/ml of LPS from S. Typhimurium. Neutrophil activation was checked by measuring the released IL-8 in the culture medium by ELISA. The manufacturer certified all reagents (serum, buffers, media and containers) as nonpyrogenic.
RNA purification and microarray hybridizations
RNA was isolated with RNeasy columns (Qiagen, Valencia, CA) based on the manufacturer´s protocol. Eluted RNA was treated with RNase-Free DNase Set (Qiagen, Valencia, CA) and the integrity, quality and quantity of RNA were checked in the Agilent Bioanalyser 2100 (Agilent Technologies, Palo Alto, CA) and NanoDrop1000 Spectrophotometer (Thermo Scientific, Wilmington, DE).
Biotinylated cRNA was obtained following the Two-Cycle Eukaryotic Target Labeling System from Affymetrix. The labeled cRNA is then cleaned up, fragmented, and hybridized to Affymetrix Porcine GeneChip according to the manufacturer´s procedures (Expression Analysis Technical Manual, Affymetrix Inc., Santa Clara, CA). Chips were scanned with an Affymetrix GeneChip Scanner 3000 (Affymetrix Inc., Santa Clara, CA).
All data were normalized using the affy library of the Bioconductor package (Gentleman et al., 2004) for R (2.10) and adapted in-house scripts. Features appearing in all replicates and having a signal of 50 over background were normalized using robust multi-array analysis (RMA) . Genes were filtered for based on Affymetrix's present/marginal/absent (PMA) calls using mismatch probe intensity, the ratio of PM to MM, only those probes present in at least 80 % of the slide was used in the analysis of microarray. Differential expression and clustering was performed with Spotfire DecisionSite v 9.0 (SP2) as follows: feature values were obtained by a one-way analysis of variance (ANOVA) with a p-value < 0.05 and hierarchically clustered using UPGMA method and Pearson´s correlation as the similarity measure. Since the Affymetrix Porcine GeneChip® is not fully annotated in all the features, it was re-annotated with Blast2GO  with a minimum E-value of 10-10 and a minimum similarity of 50%.
The networks, functional analysis and canonical pathways were generated through the use of Ingenuity Pathway Analysis (Ingenuity Systems®, http://www.ingenuity.com). IPA provides computational algorithms to identify and generate significant biological networks and pathways that are particularly enriched with our genes of interest. Networks are ranked by a score that takes into account the number of focus genes and the size of the networks, indicating the likelihood of the focus genes in a network being found together by chance. The higher the score (score = -log(p-value)), the lower the probability of finding the observed Network Eligible Molecules in a given network by chance. IPA also gives information on biological functions and canonical pathways. Thus, DEG were grouped into known biological functions, networks and canonical pathways based on human and rodent studies.
Results and discussion
Changes in gene expression profile in neutrophils stimulated with LPS
Biological analysis with Ingenuity Pathway Analysis
IPA identified significant networks, top functions and canonical pathways associated with the DGE for the time course analysis. Networks with score higher than 25 (except for cluster 1) are presented in Additional file 1. Focus genes are shown in bold. Top biofunctions and canonical pathways are listed in Additional file 2.
Cluster 1 contains 8 genes with tendency for up-regulation through the time course, including genes coding for products involved in metabolism (GRHPR, SLC38A1) and immune-related functions, such as adhesion, cellular movement and phagocytosis (SIRPA, SMARCA4) or linked to glucocorticoid signaling (HSPA14), which is known to be implicated in inflammatory response and apoptosis .
Cluster 2 includes 747 down-regulated genes at 9 hr. These molecules are involved in cellular organization, molecular transport and protein trafficking and are related to canonical pathways such as JAK/STAT signaling, Cdc42 signaling and protein ubiquitination. Network 3 (Additional file 3) is related to molecular transport, protein trafficking and cellular development. This network is focused on TGFB, a cytokine implicated in process such as neutrophil locomotion, communication among immune cells or glucocorticoid receptor signaling .
Finally, cluster 3 consists of 335 up-regulated genes. Functions associated with these molecules are related to cellular assembly and reorganization, cellular maintenance and gene expression. Canonical pathways are related to protein ubiquitination signaling, PDGF signaling and IL-3 signaling which is involved in cell survival by activation of JAK/STAT signaling and BCL2 . Network 2 (Additional file 4) highlights NF-κB interactions and covers several canonical pathways such as acute phase response signaling and interferon signaling.
Inhibition of spontaneous apoptosis at 9 hrs
MCL1 encodes myeloid leukemia cell differentiation protein, Mcl-1, that also belongs to the Bcl-2 family. Alternative splicing occurs at this locus and two transcript variants have been identified. The longer gene product (isoform 1) enhances cell survival whereas the shorter gene product (isoform 2) promotes apoptosis . MCL1 is down-regulated at 9 hr suggesting that it could be the shorter isoform of the Mcl-1 protein.
Also, there is a down-regulation of transcripts coding for some anti-bacterial products like NCF1, NCF2 and NCF4. The reactive oxygen system is autocytotoxic for neutrophils and its down-regulation may contribute to the prolongation of life span during inflammation. Therefore, many components of cell survival pathways are regulated by LPS, and these gene products may be involved in maintaining neutrophil viability. Inhibition of neutrophil apoptosis may improve the host defense against infection by increasing functional longevity of the cells at the inflammatory sites. The results of this study indicate that activated neutrophils not only generate cytokines or phagocytic-related response but also changes in the expression of genes coding for receptors, transcription factors, chromatin remodeling, cytoskeleton reorganization and apoptosis-related genes.
Authors acknowledge the computer resources and the technical support provided by the Plataforma Andaluza de Bioinformatica of the University of Malaga. This work was supported by EU funds provided by EADGENE project, and by the project of the Spanish Ministry of Science and Innovation AGL2008-0400.
This article has been published as part of BMC Proceedings Volume 5 Supplement 4, 2011: Proceedings of the International Symposium on Animal Genomics for Animal Health (AGAH 2010). The full contents of the supplement are available online at http://www.biomedcentral.com/1753-6561/5?issue=S4.
- Homburg CH, Ross D: Apoptosis of neutrophils. Curr Opin Hematol. 1996, 3: 94-99. 10.1097/00062752-199603010-00014.View ArticlePubMedGoogle Scholar
- Kasprisin DO, Harris MB: The role of RNA metabolism in polymorphonuclear leukocyte phagocytosis. J Lab Clin Med. 1977, 90: 118-124.PubMedGoogle Scholar
- Cohn ZA, Morse SI: Functional and metabolic properties of polymorphonuclear leucocytes II. The influence of a lipopolysaccharide endotoxin. J Exp Med. 1960, 111: 689-704. 10.1084/jem.111.5.689.PubMed CentralView ArticlePubMedGoogle Scholar
- Dahinden CA, Galanos C, Fehr J: Granulocyte activation by endotoxin I. Correlation between adherence and other granulocyte functions, and role of endotoxin structure on biologic activity. J immunol. 1983, 130: 857-862.PubMedGoogle Scholar
- Heit B, Liu L, Colarusso P, Puri KD, Kubes P: PI3K accelerates, but is not required for, neutrophil chemotaxis to fMLP. J Cell Sci. 2008, 121 (Pt 2): 205-14. 10.1242/jcs.020412.View ArticlePubMedGoogle Scholar
- Irizarry RA, Hobbs B, Collin F, Beazer-Barclay YD, Antonellis KJ, Scherf U, Speed TP: Exploration, normalization, and summaries of high density oligonucleotide array probe level data. Biostatistics. 2003, 4 (2): 249-64. 10.1093/biostatistics/4.2.249.View ArticlePubMedGoogle Scholar
- Conesa A, Götz S, García-Gómez JM, Terol J, Talón M, Robles M: Blast2GO: a universal tool for annotation, visualization and analysis in functional genomics research. Bioinformatics. 2005, 21 (18): 3674-3676. 10.1093/bioinformatics/bti610.View ArticlePubMedGoogle Scholar
- Franchimont D: Overview of the actions of glucocorticoids on the immune response: a good model to characterize new pathways of immunosuppression for new treatment strategies. Ann N Y Acad Sci. 2004, 1024: 124-37. 10.1196/annals.1321.009.View ArticlePubMedGoogle Scholar
- Taylor AW: Review of the activation of TGF-β in immunity. J Leukoc Biol. 2009, 85: 29-33. 10.1189/jlb.0708415.PubMed CentralView ArticlePubMedGoogle Scholar
- Reddy EP, Korapati A, Chaturvedi P, Rane S: IL-3 signalling and the role of Src kinases, JAKs and STATs: a covert liaison unveiled. Oncogene. 2000, 19 (21): 2532-47. 10.1038/sj.onc.1203594.View ArticlePubMedGoogle Scholar
- Savill J: Apoptosis in resolution of inflammation. J Leukoc Biol. 1997, 61: 375-380.PubMedGoogle Scholar
- DeLeo FR: Modulation of phagocyte apoptosis by bacterial pathogens. Apoptosis. 2004, 9: 399-413. 10.1023/B:APPT.0000031448.64969.fa.View ArticlePubMedGoogle Scholar
- Lee A, Whyte MK, Haslett C: Inhibition of apoptosis and prolongation of neutrophil functional longevity by inflammatory mediators. J Leukoc Biol. 1993, 54: 283-288.PubMedGoogle Scholar
- Hockenbery D, Nunez G, Milliman C, Schreiber RD, Korsmeyer SJ: Bcl-2 is an inner mitochondrial membrane protein that blocks programmed cell death. Nature. 1990, 348: 334-336. 10.1038/348334a0.View ArticlePubMedGoogle Scholar
- Bae J, Leo CP, Hsu SY, Hsueh AJW: MCL-1S, a splicing variant of the antiapoptotic BCL-2 family member MCL-1, encodes a proapoptotic protein possessing only the BH3 domain. J. Biol. Chem. 2000, 275: 25255-25261. 10.1074/jbc.M909826199.View ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.