Genome-wide immunity studies in the rabbit: transcriptome variations in peripheral blood mononuclear cells after in vitro stimulation by LPS or PMA-Ionomycin
© Jacquier et al.; licensee Biomed Central. 2015
Received: 25 July 2014
Accepted: 2 January 2015
Published: 23 January 2015
Our purpose was to obtain genome-wide expression data for the rabbit species on the responses of peripheral blood mononuclear cells (PBMCs) after in vitro stimulation by lipopolysaccharide (LPS) or phorbol myristate acetate (PMA) and ionomycin. This transcriptome profiling was carried out using microarrays enriched with immunity-related genes, and annotated with the most recent data available for the rabbit genome.
The LPS affected 15 to 20 times fewer genes than PMA-Ionomycin after both 4 hours (T4) and 24 hours (T24), of in vitro stimulation, in comparison with mock-stimulated PBMCs. LPS induced an inflammatory response as shown by a significant up-regulation of IL12A and CXCL11 at T4, followed by an increased transcription of IL6, IL1B, IL1A, IL36, IL37, TNF, and CCL4 at T24. Surprisingly, we could not find an up-regulation of IL8 either at T4 or at T24, and detected a down-regulation of DEFB1 and BPI at T24. A concerted up-regulation of SAA1, S100A12 and F3 was found upon stimulation by LPS.
PMA-Ionomycin induced a very early expression of Th1, Th2, Treg, and Th17 responses by PBMCs at T4. The Th1 response increased at T24 as shown by the increase of the transcription of IFNG and by contrast to other cytokines which significantly decreased from T4 to T24 (IL2, IL4, IL10, IL13, IL17A, CD69) by comparison to mock-stimulation. The granulocyte-macrophage colony-stimulating factor (CSF2) was by far the most over-expressed gene at both T4 and T24 by comparison to mock-stimulated cells, confirming a major impact of PMA-Ionomycin on cell growth and proliferation. A significant down-regulation of IL16 was observed at T4 and T24, in agreement with a role of IL16 in PBMC apoptosis.
We report new data on the responses of PBMCs to LPS and PMA-Ionomycin in the rabbit species, thus enlarging the set of mammalian species for which such reports exist. The availability of the rabbit genome assembly together with high throughput genomic tools should pave the way for more intense genomic studies for this species, which is known to be a very relevant biomedical model in immunology and physiology.
KeywordsRabbit Transcriptome Immune response Microarray Peripheral blood mononuclear cells In vitro stimulation LPS PMA-Ionomycin
The rabbit (Oryctolagus cuniculus) is a wild and domesticated species that is used for many purposes including production for meat and fur, biotechnological applications (e.g. antibody production, cloning), and as a tool for the study of nutrition, reproduction, toxicology, pharmaceutical research, and pathologies [1-4]. In fact, the rabbit has been extensively used as a model to study infectious diseases such as tuberculosis, syphilis, anthrax, tularemia, poxvirus diseases, and hepatitis E. Rabbits are also used for research on autoimmune diseases such as systemic lupus erythematosus . The rabbit immune system generates antibody diversity and optimizes affinity through mechanisms that are more efficient than those of mice and other rodents. Due to their high specificity and affinity, rabbit antibodies are commercially available for a large number of target antigens.
The rabbit has been included in the Mammalian Genome Project  and a 2X followed by a 7X coverage version of the rabbit genome sequence has been released in 2005 and 2009, respectively. The reference genome assembly OryCun2.0 (AAGW00000000.2, http://www.ncbi.nlm.nih.gov/genome/?term=oryctolagus%20cuniculus) has provided insights to decipher the impact of the domestication process on shaping genomes . Genomic tools are growing for the rabbit species, complementary to the onset of more and more cost-effective high throughput sequencing methodologies, paving the way to large scale genome-wide expression studies . However, very few genome-wide expression studies have been reported to date in the rabbit and most of which have focused on early embryogenesis [8-10] pluripotent stem cells , implantation during gestation , ocular research [13,14] and response of different organs to Eimeria infections [15-17]. Although rabbits have proven very useful in immunological and infectious disease research, genome-wide expression studies targeting the immune response to various challenges are still limited. Our aim was to create genome-wide gene expression datasets representative of various in vitro stimulations of immune cells in rabbits in order to increase knowledge on the rabbit immune responses, to provide new rabbit-specific data for comparative immunology, and to promote that transcriptomics-based approaches are efficient to provide molecular phenotype information to analyze immunity in the rabbit species. We chose to study variations of the transcriptome in rabbit peripheral blood mononuclear cells (PBMCs) using two different stimulations expected to induce either an innate immune response (lipopolysaccharide (LPS stimulation)) or a lymphocyte proliferation (mixture of phorbol myristate acetate (PMA) and ionomycin (PMA-Ionomycin stimulation)). LPS is part of the outermost layer of gram-negative bacteria and is a pathogen-associated molecular pattern (PAMP) used for in vitro studies of the innate immune response after bacterial infection. PMA, in conjunction with ionomycin, act preferentially on growth and proliferation of immune cells, and stimulate the intracellular production of the cytokines IL2 and IL4. Both stimulations were chosen as they are classically used to induce distinct types of immune responses for further characterization. This study stands as the rabbit counterpart of a similar study that was conducted in pig , and that has provided much new data on the main gene pathways that were modified following each stimulation. For this purpose, we have designed a rabbit long oligonucleotide-based DNA microarray, by enriching a previously customized gene expression design (GPL16482) with a set of well-annotated genes known to be involved in immune and inflammatory responses. We report the striking features of gene expression changes in response to both stimulations over time, relying on gene expression microarrays that are informative, cost-effective and easy to use for genome-wide expression studies in the rabbit. To our knowledge, our study reports for the first time the modifications of transcriptome profiling of PBMCs upon LPS of PMA-Ionomycin stimulation in the rabbit species.
A well-annotated gene expression platform enriched in immunity-related genes for the rabbit species
The rabbit gene expression microarray used in this study combines the customized microarray design (GPL16482) with a set of 453 specific rabbit genes involved in the immune system process (GO: 0002376), and identified to be missing on the previous array. Gene set information is available in the Additional file 1. This custom microarray has been updated with the latest annotation generated by the rabbit sequencing consortium .
Number of genes found up- or down-regulated (FDR < 0.05) after stimulation, and percentage of probe annotation
Total DE genes
Percentage of annotation
LPS at T4
LPS at T24
PMA-Ionomycin at T4
PMA-Ionomycin at T24
Clustering of samples and comparative effect of LPS and PMA-Ionomycin stimulations
Top-ten DE genes after LPS or PMA-Ionomycin stimulation, T4 and T24
Up - regulation
Down - regulation
PMA - Ionomycin
Up - regulation
Down - regulation
Genes and biological functions affected by LPS stimulation
At T24, the LPS activation revealed a switch towards a majority of genes being up-regulated (212 up vs. 165 down-regulated genes; see Table 1 and Figure 2), with higher FC values than at T4, ranging from −9.64 to 50.10 (Table 2). The most up-regulated gene was the coagulation factor III (thromboplastin, tissue factor) (F3). Among the top ten up-regulated genes (Table 2), an important subset of genes was connected to the inflammation process, and included genes that encode the proinflammatory cytokines IL6, IL1B, and IL36G, the Serum Amyloid A1 protein (SAA1), the chemokine (C-C motif) ligand 4 (CCL4), and the interleukin 10 (IL10) known to have pleiotropic effects on immunoregulation and inflammation. TNF was also found to be up-regulated (Additional file 4) together with the gene encoding the TNFAIP3 Interacting Protein 3 (TNIP3). In addition, the heparin-binding EGF-like growth factor (HBEGF) and the S100 calcium binding protein A12 (S100A12) were in the top ten up-regulated genes. The HBEGF gene has been reported to encode a cell surface binding protein potentially involved in macrophage-mediated cellular proliferation, and the S100A12 gene has been shown to be involved in cell cycle progression and differentiation. The most down-regulated gene was the Palladin, Cytoskeletal Associated Protein gene (PALLD). The top-ten down-regulated genes (Table 2) also included the bactericidal/permeability-increasing protein gene (BPI) and the defensin beta 1 gene (DEFB1). BPI encodes a LPS binding protein and has bactericidal activity on gram-negative organisms. DEFB1 belongs to the defensins reported to form a family of microbicidal and cytotoxic peptides made by neutrophils, and encodes an antimicrobial peptide involved in the resistance of epithelial surfaces to microbial colonization.
List of the 10 genes found DE at T4 and T24 upon LPS stimulation
For the global analysis of biological functions affected by LPS stimulation, the DE gene lists were submitted to the Ingenuity Pathway Analysis (IPA) system. From this analysis 262 (82.9% of DE genes) and 322 (85.4% of DE genes) genes were identified by IPA at T4 and T24, respectively. The global picture provided by Figure 2 highlights that most genes were found down-regulated at T4 conversely to T24. The foremost over-represented down-regulated biological functions (P < 0.05) at T4 were related to Molecular and Cellular Functions (Cellular Function and Maintenance, Cell Death and Survival, Cellular Growth and Proliferation). However, a number of genes involved in biological functions related to inflammatory disease or response were repressed at T4. By contrast, at T24, there was a significant enrichment of up-regulated genes associated with biological functions related to Diseases and Disorders such as inflammatory response, Inflammatory Disease and Infectious Disease. For the Molecular and Cellular Functions, the most represented functions were Cell Death and Survival, Cellular Growth and Proliferation, Cellular Movement, Cellular Development and Cell-to-Cell Signaling and Interactions. For the Physiological System Development and Functions, the most represented function was the Hematological System Development and Function.
Transcription profile of cytokines in PMBCs stimulated by LPS at T4 and T24
Genes and biological functions affected by PMA-Ionomycin stimulation
Two large sets of genes were found DE at T4 and T24 (Table 1), with a slight tendency toward more up-regulated genes than down-regulated genes at both times (2273 down- versus 2477 up-regulated probes at T4; 3414 down- versus 4039 up-regulated probes at T24). As shown in Table 2 (and Additional files 5 (T4) and Additional file 6 (T24)), there was a decrease over time in the FC values for up-regulated genes (maximum FC was 227.98 and 57.58 at T4 and T24, respectively) but an increase over time in the FC of down-regulated genes (minimum FC was −48.36 and −146.61 at T4 and T24, respectively). A set of 3120 DE probes was shared at both times, and could be divided into a subset of 2378 probes that were continually DE in the same directions at T4 and T24 (1091 probes down-regulated and 1287 probes up-regulated), a subset of 515 probes that were down-regulated at T4 and up-regulated at T24, and a third subset of 227 probes that were up-regulated at T4 but down-regulated at T24.
The genes CSF2, IL2, and RGCC (gene regulator of cell cycle) were among the top ten genes found over-expressed at T4 and T24. The granulocyte-macrophage colony-stimulating factor (CSF2) was by far the most over-expressed gene with a positive FC of 227.98 at T4 and 57.58 at T24 (Table 2). CSF2 encodes a cytokine that controls the production, differentiation, and function of granulocytes and macrophages, and the RGCC gene encodes a protein thought to regulate cell cycle progression. All top ten over-expressed genes at T4 were still up-regulated at T24. Conversely, the genes BCO2, CD70, CCNB1, and NCAPG found in the top ten over-expressed genes at T24 were found to be down-regulated at T4, showing a switch in gene regulation between T4 and T24. The beta-carotene oxygenase 2 gene (BCO2) encodes an enzyme, which oxidizes beta-carotene during the biosynthesis of vitamin A. CD70 molecule gene (CD70) encodes a cytokine that belongs to the TNF ligand family. CD70 is a ligand for TNFRSF27/CD27 and a surface antigen on activated B and T lymphocytes, and contributes to T cell activation. The cyclin B1 gene (CCNB1) encodes a regulatory protein involved in mitosis. The non-SMC condensin I complex, subunit G gene (NCAPG) encodes a subunit of the condensin complex, which is responsible for the condensation and stabilization of chromosomes during mitosis and meiosis.
The genes LTB and FPR1 were both found in the top ten most repressed genes at T4 and T24 post stimulation. The lymphotoxin beta (TNF superfamily, member 3) (LTB) encodes a protein which is an inducer of the inflammatory response system and is involved in normal development of lymphoid tissue. In addition, there was an intensification of the repression of the two genes S100A8 (S100 calcium binding protein A8) and S100A12 (S100 calcium binding protein A12), which encode proteins of the S100 family containing 2 EF-hand calcium-binding motifs, involved in cell cycle progression and differentiation. The defensin beta 1 (DEFB1) gene was repressed only at T24.
Transcription profiles of cytokines and transcription factors in PMBCs stimulated by PMA-Ionomycin at T4 and T24
The transcription factors TBX21, STAT5A, and STAT5B were up-regulated at T4 and T24, whereas GATA1 was slightly repressed at both times (Table 5). TBX21 is the T-box 21 gene that stands as the human ortholog of the mouse Tbx21/Tbet gene. This gene was reported to act as a Th1 cell-specific transcription factor and controls the expression of IFNG . STAT5A and STAT5B are both members of the signal transducer and activator of transcription (STAT) family of transcription factors that respond to cytokines and growth factors. STAT5A encodes a protein that is activated by IL2 and CSF2, and STAT5B encodes a protein reported to mediate the signal transduction triggered by various cell ligands, such as IL2 and IL4. At T24, two additional members of the STAT family of transcription factors, STAT4 and STAT1, were found up-regulated and two other members of the GATA families were slightly repressed (GATA2 and GATAD2B). FOXP3 was not found affected by the PMA-Ionomycin stimulation. We checked the differential expression of the thrombospondin gene (THBS1) that was found to be the most repressed gene in porcine PBMCs stimulated for 24 hours by PMA-Ionomycin . In rabbit PBMCs, the THBS1 gene was repressed only at T24 and with a limited FC of −4.75.
The genes encoding class II molecules of the major histocompatibility complex (MHC) were found to be repressed. The most repressed genes belong to the DR series followed by DM, DQ and DP. This repression was initiated early after stimulation since DMA, DMB and DRB1 were already found down-regulated at T4 (Additional file 7). In addition, few probes corresponding to MHC class I transcripts were found down-regulated. Interestingly, and in agreement with data obtained in pig , we observed an up-regulation of genes encoding heat shock proteins (HSPs) and many genes involved in the cascade of peptide processing before loading to the MHC molecule binding groove (components of the proteasome (PSMB gene series), calnexin (CANX), tapasin (TAP1)), despite an overall down-regulation of genes encoding class I and II molecules (Additional file 7).
At T4 and T24, subsets of 3575 (75.3%) and 5153 (69.1%) DE genes were mapped into the IPA system. The relative representation of biological functions between T4 and T24 was quite similar (Figure 2). By excluding all genes associated to Cancer, IPA revealed an enrichment of genes principally related to Molecular and Cellular Functions (cell death and survival, cellular growth and proliferation, cellular development, cellular function and maintenance), followed by an enrichment of genes associated with Diseases and Disorders (gastrointestinal disease, infectious disease, inflammatory response), and Physiological System Development and Function (organismal survival, organismal development, hematological system development and function). These observations were reinforced by the analysis of the set of 3123 genes found to be DE at both T4 and T24, and showed a predominance of genes associated with Molecular and Cellular Functions. The 515 genes found down-regulated at T4 and up-regulated at T24 were predominantly associated to the biological functions Cell Cycle, DNA Replication Recombination and Repair, and Cellular Assembly and Organization. For the 227 genes found up-regulated at T4 and down-regulated at T24, there was an enrichment in biological functions relating to Cell Death and Survival, Cell Growth and Proliferation, Cell Development, Tissue Morphology, and Hematological System Development and Function. For these two gene sets, the enriched biological functions aligned well with those found for the genes that revealed no inverted regulation compared to mock-stimulation at T4 and T24.
Validation of differentially expressed genes by qRT-PCR
Cytokine profile induced by LPS stimulation in rabbit PBMCs
LPS is an outer membrane component of Gram-negative bacteria and the crystal structure of the trans-envelope protein complex inserting LPS in the outer leaflet of the outer envelope has been very recently described [20,21]. Lipid A is the main PAMP of LPS. Interactions between LPS and TLR4, a sensor for LPS  at the surface of macrophages and dendritic cells, trigger the expression of proinflammatory cytokines and activate the innate immune system by MyD88 (Myeloid differentiation primary response gene) dependent or MyD88 independent pathways . In our study, an early up-regulation of IL12A and CXCL11 together with a slight repression of CD14 was found at T4. This LPS-induced CD14 repression observed at T4 was comparable to in vivo human studies , highlighting the immunodepressive effect of LPS on rabbit PBMCs. Up-regulation of a wider range of inflammatory cytokines was observed at T24 with a differential expression between mock-stimulation and LPS-stimulation for IL6, IL1A, IL1B, and TNF. In LPS-stimulated human whole blood , it has been reported that gene activations of both IL6 and TNF were induced early, after only 2–4 hours and then rapidly declined. Several studies have shown the synergy between IL1B and TNF [26,27]. When both cytokines were administered to healthy humans, they induced fever, inflammation, tissue destruction, or resulted in the development of sickness behavior, and depressed mood when combined with IL6 [28,29]. Moreover, IL1B and TNF have been shown to stimulate the expression of IL36G in keratinocytes [30,31], confirming results at T24. Our results suggest a significant triggering of the MyD88 dependent pathway in dendritic cells that is characterized by IL12 production and followed by a T cell differentiation toward Th1 phenotype [32,33], as represented in Figure 5A. The onset of the inflammation process was detected at T4 but clearly amplified at T24 as represented in Figure 5B, showing that the data obtained at T4 did not provide an accurate overall picture of gene expression picture upon LPS stimulation. In addition, gene transcription levels did not return to basal levels at T24, as already shown in pigs . This late increase in inflammatory cytokine transcription contrasts with results on bovine PBMCs for which the inflammation process was triggered much earlier with a maximal expression of inflammation gene expression at T4 . For both T4 and T24 after LPS activation, the interferon regulatory factor 5 (IRF5) was found to be under-expressed. In 2011, Krausgruber et al.  highlighted a critical role for the transcription factor IRF5 in human M1 macrophage polarization and defined it as a transcriptional repressor. Moreover, an activation of the Th1 response by IRF5 has been reported.
Putative specific features of rabbit PBMCs upon LPS stimulation
IL8 (CXCL8 in human) was not DE in our experimental conditions at the transcript level at both T4 and T24, whereas nucleotide probes specific for this gene are available on the microarrays used for the transcriptome analyses. Indeed IL8 is a chemokine known to act as one of the major mediators of the inflammatory response and has been reported to be up-regulated after LPS stimulation in many mammals including cattle [34,36] and pig [37,18]. As opposed to other species, the possibility that T4 might be too late for detecting an up-regulation of IL8 in rabbits cannot be excluded. Recent reports have shown that bacterial-host interactions affect histone acetylation, phosphorylation and methylation state at the TLR4 and IL8 promoter in host cells [38,39]. In cattle, it has been hypothesized that host individual variability in response to LPS stimulation could be due to epigenetic variations modulating individual IL8 expression in PBMCs  and in dermal fibroblasts .
In addition, DEFB1 and BPI were down-regulated at T24 compared to mock-stimulation, whereas they were expected to be up-stimulated due to their bactericidal activity [40-42]. Defensins form a family of microbicidal and cytotoxic peptides, and DEFB1 encodes an antimicrobial peptide implicated in the resistance of epithelial surfaces to microbial colonization . Beta defensins have thus an antimicrobial activity that defend epithelial surfaces including the skin, gastrointestinal, and respiratory tracts. BPI encodes a LPS binding protein associated with neutrophil granules in humans, and has bactericidal activity on gram-negative organisms . Upon LPS activation, BPI may have several functions, as inhibition of endothelial cell growth and inhibition of dendritic cell maturation, or as an anti-angiogenic, chemoattractant or opsonization agent . To our knowledge, no study has shown down-regulation of DEFB1 or BPI in rabbit PBMCs upon LPS activation. However, a down-regulated expression of DEFB1 has been shown following infections by several pathogens (Shigella and Cryptosporidium parvum) in epithelial cells of the digestive and respiratory tracts in humans or mice [46-48]. In our experimental model, we do not have a clear explanation for the down-regulation of DEFB1 and BPI, but our results suggest the onset of an immune escape mechanism upon LPS stimulation that needs further investigation. We cannot rule out that both genes are induced early after stimulation, before the time window chosen for this study. Additional data are required to better assess why no increase of IL8 transcription was detected during the global inflammatory process induced by LPS stimulation in rabbit PBMCs between 4 and 24 hours post stimulation, and why a down-regulation of bactericidal activity-related genes was found at T24. It could be interesting to carry out new kinetics studies that include earlier time points (e.g. one hour post-stimulation), as well as a mid-time point (e.g. 12 hours post stimulation) and a late time point at 48 hours post-stimulation, thus providing new data to validate or not the hypothesis that rabbits differ from mice and humans for the in vitro production of IL8, and up-regulation of DEFB1 and BPI by PBMCs upon LPS stimulation.
Sets of genes commonly affected by LPS stimulation in various mammalian species
Our results have confirmed a wide range of genes, other than cytokines and chemokines, which were affected by LPS stimulation in rabbit PBMCs. This has already been reported in other species, thus illustrating shared innate immunity responses among mammals.
The SAA1 gene belongs to the family encoding the serum amyloid A (SAA) proteins. As reported in pigs , the SAA1 gene was found up-regulated upon LPS stimulation in PBMCs at T24. The association of increased levels of SAA proteins with inflammatory states has been largely reported , but the precise functions of SAA proteins still need clarification. SAA has been reported to be involved in the establishment and maintenance of inflammation, and as an opsonin it facilitates phagocytosis of Gram-negative bacteria . SAA acts as an antiapoptotic agent for neutrophils , and directly drives a regulatory program in neutrophils by stimulating them to produce IL10 that has anti-inflammatory effects . The strong up-regulation of IL10 detected in rabbit PBMCs at T24 might be related to this role of SAA. In addition, SAA was recently shown to induce mitogenic signals in regulatory T cells, by enhancing the secretion of IL1B and IL6 by monocytes, and driving them toward a phenotype resembling immature dendritic cells . Indeed, IL6 and IL1B genes were both strongly up-regulated together with SAA1. Our results suggest that the rabbit is an interesting species for further study of SAA-associated regulatory mechanisms, which could be shared across mammals.
The S100A12 gene was strongly up-regulated in rabbit LPS-stimulated PBMCs at T24, as reported in pig PBMCs under similar LPS stimulation conditions . S100A12 gene belongs to the family that encodes the S100 proteins. Approximately 25 different S100 proteins have been characterized so far . A common feature shared by all the S100 proteins is the presence of a pair of calcium-binding helix-loop-helix domains referred to as EF-hand calcium-binding regions toward either end of the protein and separated by a hinge region . The broad diversity of S100 proteins is related to a multitude of biological functions such as cell growth, cell differentiation, and cell survival in numerous physiological and pathological conditions in all cells of the body. Furthermore, all S100 proteins are likely to induce some changes in the actin cytoskeleton organization . S100A12 is expressed in the myeloid cell lineage, and is found in abundance in granulocytes, monocytes, and lymphocytes in human. Extracellular S100A12 can induce directional migration and chemotactic responsiveness of monocytes and neutrophils in vitro .
At T24 post LPS stimulation, we detected a very strong transcription increase of the gene F3 that encodes the coagulation factor III or thromboplastin, a cell surface glycoprotein that enables cells to initiate the blood coagulation cascades. It has been long reported that blood cells generate thromboplastin activity upon incubation with LPS or Gram-negative bacteria . It has been further established that monocytes are the blood cells that possess this property , and that collaboration with lymphocytes is required for a rapid induction of procoagulant activity . Induction of pro-coagulant activity of hepatic macrophages by LPS has already been demonstrated in rabbits . In this paper, we demonstrate by a genome-wide and a non-focused approach that a pro-coagulant activity is also strongly induced by LPS in PBMCs in the rabbit.
A concerted action of S100A12, SAA, and thrombotic functions during inflammatory diseases has been reported in humans . Here, we clearly show that in rabbit PBMCs stimulated in vitro by LPS, a strong combined increase in the expression of the corresponding rabbit genes (S100A12, SAA1, F3) is detected.
PMA-Ionomycin stimulation affects many more genes than LPS stimulation
PMA is a diester of phorbol and a potent tumor promoter often used in biomedical research to activate the signal transduction enzyme protein kinase C (PKC). For in vitro cell stimulation, PMA is used in conjunction with ionomycin, a Ca2+ ionophore produced by the bacterium Streptomyces conglobatus. The number of DE genes after PMA-Ionomycin stimulation was found to be 15 and 20 times higher at T4 and T24, respectively, by comparing to the number of DE genes after LPS stimulation. These findings are in agreement with a similar study in swine in which we reported that ten times more genes were DE after PMA-Ionomycin stimulation than after LPS stimulation . This difference is consistent with the dynamics of the responses to LPS and PMA-Ionomycin, since it was observed at T4 and maintained at T24. As previously discussed , this observation might reflect that LPS targets only a limited set of monocytes and macrophages expressing CD14 , while PMA-Ionomycin has a much wider spectrum of target cells. We found 515 probes that were down- and up-regulated at T4 and T24, respectively, and vice versa for another 227 probe set. All these genes relate to the same biological functions than those associated to genes found only up- or down-regulated at both T4 and T24, thus likely reflecting a dynamic process that includes combined actions of genes toward a concerted immune response.
PMA-Ionomycin stimulation induced early Th2, iTreg (induced Tregs), and Th17 responses at T4 followed by an increase in the Th1 response at T24
At T4, the cytokine profile revealed a very strong up-regulation of IL13 and IL4, which are produced by Th2 population subsets . Both cytokines were still up-regulated at T24 in comparison to mock-stimulated cells but with a much lower fold-change than at T4 (Table 5). We could not detect up-regulation of transcription factors that induce cells to commit to a Th2 phenotype such as STAT6 and GATA3. Similarly, IL10 could be detected at T4 but the up-regulation was decreased at T24, as previously described in pigs . This cytokine is identified as a hallmark of iTreg lymphocytes . The FOXP3 gene known to be involved in iTreg lymphocyte differentiation was not found to be affected by PMA-Ionomycin stimulation. In addition, IL17F was moderately up-regulated at T4 (Table 5) in comparison to mock-stimulated PBMCs, and IL17A was found to be DE at both T4 and T24 with a decrease in the FC between stimulated and non-stimulated PBMCs at T24 in comparison to T4. All these results suggest a very early expression of Th2, iTreg, and Th17-associated cytokines with a progressive slow-down after T4. Conversely, we detected an increase in Th1 phenotypes between T4 and T24. Indeed, the two lineage-defining transcription factors STAT4 and TBX21 (or T-bet) were significantly up-regulated at T4, and at T4 and T24, respectively. STAT4 and TBX21 genes drive the Th1 cell commitment and IFNG is the main Th1 cytokine. These data show that PMA-Ionomycin stimulation induces a Th1 response that is not only maintained, but was found to increase at T24. This coincides with results we have previously reported, where a predominance of a Th1 response after PMA-Ionomycin treatment of PBMCs was seen in swine for T24. Our results in rabbits suggest that the Th1 response induced by PMA-Ionomycin might be maintained on a longer term than other responses, thus confirming a predominance of a Th1 response after 24 h stimulation as previously shown in pigs . Recent reports have demonstrated a plasticity of T-helper cell differentiation , thus suggesting that T cell differentiation is not always the result of a direct differentiation from naïve T cells but that new cell commitments may also occur after a reprogramming of differentiated cells . From our data, we cannot determine whether the observed Th1 phenotype increase is due to the reprogramming of differentiated T cells or to the direct differentiation of naïve T cells.
A strong up-regulation of CSF2 and a down-regulation of IL16 induced by PMA-Ionomycin stimulation in favor of cell proliferation
The most up-regulated gene at T4 and T24 after PMA-Ionomycin stimulation was CSF2, a cytokine that functions as a white blood cell growth factor. CSF2 stimulates stem cells to produce granulocytes (neutrophils, eosinophils, and basophils) and monocytes . By contrast, the cytokine IL16 was down-regulated at both T4 and T24, with a high negative fold-change.
Several studies have shown that IL16 is constitutively expressed in peripheral blood monocytes and that there is a direct link between IL16 expression and apoptosis [64,65], and with autoimmune and allergic diseases . Whereas some genes (including IL10, IL16, and TLR6) have been shown to be down-regulated during Fas-mediated or camptothecin-induced apoptosis in human polymorphonuclear leukocytes , our results confirmed previous studies where down-regulation of IL16 is associated with T cell activation and proliferation [68,69]. IL16 has also been described as an important growth-promoting factor in multiple myeloma .
In this report, we have characterized the gene profiling modification of rabbit PBMCs upon in vitro stimulation by either LPS or a mixture of PMA and ionomycin, thus providing genome-wide expression data on immune responses for the rabbit species. Our results were consistent with data previously cited for other mammalian species, thus confirming that the rabbit is a relevant model for immunity studies.
Presently, the rabbit is no longer a species lacking genomic tools. The availability of a 7X coverage of the genome sequence , together with the progress of related information on annotation and SNP variations, constitute highly valuable tools for launching large-scale genome-wide studies in this species . For transcriptomic studies, it is now widely acknowledged that RNAseq-approaches are well adapted to exhaustively detect all transcription variations between different conditions together with SNP variations among individuals. The rabbit is a very useful model in biomedical research [4,71]. As a livestock species, there is a need for large scale studies on resistance to disease and reduction of the use of antibiotics in production systems. The rabbit species does not have the same financial means as other species such as cattle. This study shows that a commercial well-annotated microarray is now available for rabbits, and can be used in combination with RNAseq approaches. Due to their efficiency in the study of immune responses, as well as their cost-effectiveness, the use of microarrays can be easily expanded to large cohorts of animals either for biomedical research or for improving knowledge on immunity and resistance to diseases to promote safe and efficient farm systems.
Enrichment and annotation steps
The first step was to compare the existing Agilent customized microarray design 8x60K (GPL16482) with the set of genes involved in the pig immune response  and then using the Gene Ontology term GO:0002376 (immune system process) on rabbit, pig, mouse and human species. Immunity-related rabbit sequences were retrieved by GeneID and RefSeq search or by analysis for sequence similarity by BLAST. All 60-mer length oligonucleotide probes were designed using eArray (Agilent technologies). All new probes were spotted in triplicates. The platform was registered under the GEO accession number GPL16709.
The second step was the annotation of the probes present on the microarray. Preliminary work based on chromosomal coordinates was made. Then we performed a BLAST search of the rabbit probe sequences versus a new RNAseq annotated database (Broad Institute, Cambridge, MA, USA).
PBMCs: isolation and stimulation
All experiments involving animals complied with the institutional guidelines of the French Ethical Committee. All animals were bred and slaughtered in an approved animal facility (UCEA, approval number A-78-322-4). Experimental protocols were conducted under the supervision of the animal welfare committee in charge of this animal facility (Protocol registered under number 272). V. Duranthon holds an authorization certificate issued by the French governmental administration to experiment on animals (Authorization number B78-101). PBMCs were obtained from 13 mL of heparinized whole blood of four adult female New Zealand White rabbits (numbered 1, 2, 4, and 5) by density gradient separation using ficoll-histopaque. A total of 5 × 105 cells were suspended into RPMI 1640 medium (BioWhittaker, Belgium) supplemented with 200 mM/L of L-glutamine (Sigma-Aldrich, Saint Louis, MO, USA), 10% heat inactivated fetal bovine serum, 1% penicillin and streptomycin (Life Technologies, Carlsbad, CA, USA). PBMCs were stimulated with 1 μg/mL of LPS (Sigma, France) or a mixture of PMA (Sigma, France) at 10 ng/mL and ionomycin (Sigma, France) at 1 μg/mL. For mock-stimulation, a buffer without stimulation agents was added to the cells similarly to stimulation conditions. Cells were incubated at 37°C with 5% CO2 and harvested 4 and 24 hours post-stimulation for further RNA extraction.
Total RNA was extracted from PBMCs using the RNeasy Mini Kit (Qiagen, Valencia, CA, USA) and residual contaminating genomic DNA was cleaned by using RNase-free DNase I (Qiagen, Valencia, CA, USA). RNA samples were quantified using a NanoDrop ND-1000 spectrophotometer (Thermo Fisher Scientific, Waltham, MA, USA). The RNA integrity was then assessed on a Bioanalyzer 2100 using RNA 6000 Nano chips (Agilent Technologies, Santa Clara, CA, USA).
RNA labeling and microarray processing
Transcriptional profiling was performed using the previously described custom array. A total of 24 samples was processed. All steps were performed by the CRB GADIE facility (INRA Jouy-en-Josas, France, http://crb-gadie.inra.fr/). Cyanine-3 (Cy3) labeled cRNA were prepared using 200 ng of total RNA using the One-Color Low Input Quick Amp Labeling kit (Agilent Technologies, Santa Clara, CA, USA) following the recommended protocol. Specific activities and cRNA yields were determined using the NanoDrop ND-1000 (Thermo Fisher Scientific, Waltham, MA, USA).
For each sample, 600 ng of Cy3-labeled cRNA (specific activity > 6.0 pmol Cy3/μg of cRNA) were fragmented at 60°C for 30 minutes in a reaction volume of 25 μl containing 25× Agilent Fragmentation Buffer and 10× Agilent Blocking Agent, following the manufacturer’s instructions. Subsequently, 25 μl of 2× Agilent Hybridization Buffer were added to the fragmentation mixture and hybridized to the Rabbit Custom Gene Expression Microarrays (Agilent Technologies, design ID: 042421) for 17 hours at 65°C in a rotating Agilent hybridization oven (Agilent Technologies). After hybridization, microarrays were washed 1 minute at room temperature with the GE Wash Buffer 1 (Agilent Technologies) and 1 minute at 37°C using the GE Wash Buffer 2 (Agilent Technologies), then dried instantly.
Immediately after washing, the slides were scanned using a G2565CA Scanner System (Agilent Technologies), which used a scan protocol with a resolution of 3 μm and a dynamic range of 20 bits. The resulting .tiff images were analyzed with the Feature Extraction Software v10.7.3.1 (Agilent Technologies), using the GE2_107_Sep09 protocol. The microarray data were submitted to the GEO and received the accession number GSE59263.
Microarray data were analyzed using the R/Bioconductor software package Limma (Linear Models for Microarray Data) . The median signals of all probes were first Log 2-transformed, and then these data were normalized by a median normalization. After averaging the normalized data of the probes targeting the same genes, a linear model was fitted for all probes using lmFit function, with type of stimulation and time as factors. As we had a low number of biological replicates, we also used the empirical Bayes approach to compute moderated t-statistics and log-odds of differential expression. The differential analysis was done with averaged and single probes, and DE genes were identified based on contrasts for the (stimulation × time levels) interactions. The P-values were corrected for multiple testing using a false discovery rate method (q-value < 0.05), which provides an estimate of the fraction of false discoveries among the significant terms. A PCA and a HCA were performed. In the results and discussion, the FC values were transformed into linear values for a better biological understanding. Venn diagrams were produced using the Limma package . Finally, the significant DE genes were analyzed with IPA (Ingenuity Systems, Mountain View, CA; http://www.ingenuity.com) to obtain the top biological functions and relevant networks.
Real time quantitative reverse transcription PCR
For real-time PCR analyses, 100 ng of DNase I treated total RNA were reverse-transcribed using Superscript III with random hexamers (Life Technologies, Carlsbad, CA, USA). Triplicate reactions were performed in a final volume of 20 μL, mixing 20 μL cDNA, 120 nM primers and 10 μL SYBR Green PCR Master Mix (Life Technologies, Carlsbad, CA, USA), using a QuantStudio 12 K Flex system (Life Technologies, Carlsbad, CA, USA). Primers were defined either using the Primer3 Software (Life Technologies, Carlsbad, CA, USA) or manually designed (Additional file 9). The genes B2M and GAPDH were chosen as housekeeping genes and the 2-∆∆Ct method was used to calculate the fold change in gene expression. For comparison with microarray data, Spearman correlation was calculated between the Log2 (FC) provided by the microarray study and the Log2 (Ct) provided by the RT-qPCR study .
Peripheral blood mononuclear cell
Phorbol myristate acetate
False discovery rate
Quantitative Real-Time PCR
The authors are grateful to Federica Di Palma, Jessica Alföldi, and Jeremy Johnson (Broad Institute, Cambridge, MA, https://www.broadinstitute.org) for their scientific support on the annotation of the rabbit genome. We also thank the RGB-Net COST action TD1101 for efficient networking. We would like to acknowledge the SIGENAE platform (http://www.sigenae.org/) for their support on the microarray annotation and the CRB GADIE for their help in submitting the data. The project was funded by the Carnot Institute ICSA (AGIS project). VJ’s PhD fellowship is jointly funded by the three INRA divisions: Animal Genetics, Animal Health, and Animal Physiology. VJ was also supported by Agreenium (EIR-A program).
- Fortun-Lamothe L., Gidenne T. Recent advances in the digestive physiology of the growing rabbit. In Recent Advances in Rabbit Sciences. Edited by Maertens L., Coudert P. ILVO: Melle, Belgium; 2006:201–210.
- Rogel-Gaillard C, Ferrand N, Hayes H. Rabbit. In Genome Mapping and Genomics in Domestic Animals. Edited by C. Kole and N. Cockett. Springer Berlin Heidelberg. 2009:165–230.
- Fischer B, Chavatte-Palmer P, Viebahn C, Navarrete Santos A, Duranthon V. Rabbit as a reproductive model for human health. Reproduction. 2012;144:1–10.View ArticlePubMedGoogle Scholar
- Miller I, Rogel-Gaillard C, Spina D, Fontanesi L, De Almeida AM. The rabbit as an experimental and production animal: from genomics to proteomics. Curr Protein Pept Sci. 2014;15:134–45.View ArticlePubMedGoogle Scholar
- Rai G, Ray S, Milton J, Yang J, Ren P, Lempicki R, et al. Gene expression profiles in a rabbit model of systemic lupus erythematosus autoantibody production. J Immunol. 2010;185:4446–56.View ArticlePubMed CentralPubMedGoogle Scholar
- Lindblad-Toh K, Garber M, Zuk O, Lin MF, Parker BJ, Washietl S, et al. A high-resolution map of human evolutionary constraint using 29 mammals. Nature. 2011;478:476–82.View ArticlePubMed CentralPubMedGoogle Scholar
- Carneiro M, Rubin CJ, Di Palma F, Albert FW, Alföldi J, Martinez-Barrio A, et al. The rabbit genome provides new insight into the genetics of animal domestication. Science, in press.
- Léandri RD, Archilla C, Bui LC, Peynot N, Liu Z, Cabau C, et al. Revealing the dynamics of gene expression during embryonic genome activation and first differentiation in the rabbit embryo with a dedicated array screening. Physiol Genomics. 2009;36:98–113.View ArticlePubMedGoogle Scholar
- Naturil-Alfonso C, Vicente JS, Peñaranda DS, Marco-Jiménez F. Up-regulation of insulin-like growth factor I and uteroglobin in in vivo-developed parthenogenetic embryos. Reprod Domest Anim. 2013;48:126–30.View ArticlePubMedGoogle Scholar
- Saenz-de-Juano MD, Marco-Jiménez F, Peñaranda DS, Joly T, Vicente JS. Effects of slow freezing procedure on late blastocyst gene expression and survival rate in rabbit. Biol Reprod. 2012;87:91.View ArticlePubMedGoogle Scholar
- Osteil P, Tapponnier Y, Markossian S, Godet M, Schmaltz-Panneau B, Jouneau L, et al. Induced pluripotent stem cells derived from rabbits exhibit some characteristics of naïve pluripotency. Biol Open. 2013;2:613–28.View ArticlePubMed CentralPubMedGoogle Scholar
- Saenz-de-Juano MD, Marco-Jimenez F, Schmaltz-Panneau B, Jimenez-Trigos E, Viudes-de-Castro MP, Peñaranda DS, et al. Vitrification alters rabbit foetal placenta at transcriptomic and proteomic level. Reproduction. 2014;147:789–801.View ArticlePubMedGoogle Scholar
- Popp MP, Liu L, Timmers A, Esson DW, Shiroma L, Meyers C, et al. Development of a microarray chip for gene expression in rabbit ocular research. Mol Vis. 2007;13:164–73.PubMed CentralPubMedGoogle Scholar
- Ker-Woon C, Abd Ghafar N, Hui CK, Mohd Yusof YA. Effect of acacia honey on cultured rabbit corneal keratocytes. BMC Cell Biol. 2014;15:19.View ArticlePubMed CentralPubMedGoogle Scholar
- Dkhil MA, Abdel-Maksoud MA, Al-Quraishy S, Abdel-Baki AA, Wunderlich F. Gene expression in rabbit appendices infected with Eimeria coecicola. Vet Parasitol. 2012;186:222–8.View ArticlePubMedGoogle Scholar
- Al-Quraishy S. Exogenous and endogenous stages of Eimeria perforans naturally infected domestic rabbit (Oryctolagus cuniculus) in Saudi Arabia: Light microscopic study. Saudi J Biol Sci. 2012;19:31–4.View ArticlePubMed CentralPubMedGoogle Scholar
- Dkhil MA, Al-Quraishy S, Abdel-Baki AA, Delic D, Wunderlich F. Eimeria coecicola: spleen response of Oryctolagus cuniculus. Exp Parasitol. 2013;133:137–43.View ArticlePubMedGoogle Scholar
- Gao Y, Flori L, Lecardonnel J, Esquerré D, Hu ZL, Teillaud A, et al. Transcriptome analysis of porcine PBMCs after in vitro stimulation by LPS or PMA-ionomycin using an expression array targeting the pig immune response. BMC Genomics. 2010;11:292.View ArticlePubMed CentralPubMedGoogle Scholar
- O'Shea JJ, Paul WE. Mechanisms underlying lineage commitment and plasticity of helper CD4+ T cells. Science. 2010;327:1098–102.View ArticlePubMed CentralPubMedGoogle Scholar
- Qiao S, Luo Q, Zhao Y, Zhang XC, Huang Y. Structural basis for lipopolysaccharide insertion in the bacterial outer membrane. Nature. 2014;511:108–11.View ArticlePubMedGoogle Scholar
- Dong H, Xiang Q, Gu Y, Wang Z, Paterson NG, Stansfeld PJ, et al. Structural basis for outer membrane lipopolysaccharide insertion. Nature. 2014;511:52–6.View ArticlePubMedGoogle Scholar
- Poltorak A, He X, Smirnova I, Liu MY, Van Huffel C, Du X, et al. Defective LPS signaling in C3H/HeJ and C57BL/10ScCr mice: mutations in Tlr4 gene. Science. 1998;282:2085–8.View ArticlePubMedGoogle Scholar
- Lu J, Bai L, Sun H, Nikolovska-Coleska Z, McEachern D, Qiu S, et al. SM-164: a novel, bivalent Smac mimetic that induces apoptosis and tumor regression by concurrent removal of the blockade of cIAP-1/2 and XIAP. Cancer Res. 2008;68:9384–93.View ArticlePubMed CentralPubMedGoogle Scholar
- Hojman H, Lounsbury D, Harris H, Horn JK. Immunodepressive effects of LPS on monocyte CD14 in vivo. J Surg Res. 1997;69:7–10.View ArticlePubMedGoogle Scholar
- DeForge LE, Remick DG. Kinetics of TNF, IL-6, and IL-8 gene expression in LPS-stimulated human whole blood. Biochem Biophys Res Commun. 1991;174:18–24.View ArticlePubMedGoogle Scholar
- Dinarello CA. Proinflammatory cytokines. Chest. 2000;118:503–8.View ArticlePubMedGoogle Scholar
- Raices RM, Kannan Y, Sarkar A, Bellamkonda-Athmaram V, Wewers MD. A synergistic role for IL-1beta and TNFalpha in monocyte-derived IFNgamma inducing activity. Cytokine. 2008;44:234–41.View ArticlePubMed CentralPubMedGoogle Scholar
- Dantzer R, O'Connor JC, Freund GG, Johnson RW, Kelley KW. From inflammation to sickness and depression: when the immune system subjugates the brain. Nat Rev Neurosci. 2008;9:46–56.View ArticlePubMed CentralPubMedGoogle Scholar
- Hannestad J, Subramanyam K, Dellagioia N, Planeta-Wilson B, Weinzimmer D, Pittman B, et al. Glucose metabolism in the insula and cingulate is affected by systemic inflammation in humans. J Nucl Med. 2012;53:601–7.View ArticlePubMed CentralPubMedGoogle Scholar
- Carrier Y, Ma HL, Ramon HE, Napierata L, Small C, O'Toole M, et al. Inter-regulation of Th17 cytokines and the IL-36 cytokines in vitro and in vivo: implications in psoriasis pathogenesis. J Invest Dermatol. 2011;131:2428–37.View ArticlePubMedGoogle Scholar
- Lian LH, Milora KA, Manupipatpong KK, Jensen LE. The double-stranded RNA analogue polyinosinic-polycytidylic acid induces keratinocyte pyroptosis and release of IL-36γ. J Invest Dermatol. 2012;132:1346–53.View ArticlePubMed CentralPubMedGoogle Scholar
- Kaisho T, Hoshino K, Iwabe T, Takeuchi O, Yasui T, Akira S. Endotoxin can induce MyD88-deficient dendritic cells to support T(h)2 cell differentiation. Int Immunol. 2002;14:695–700.View ArticlePubMedGoogle Scholar
- Cekic C, Casella CR, Sag D, Antignano F, Kolb J, Suttles J, et al. MyD88-dependent SHIP1 regulates proinflammatory signaling pathways in dendritic cells after monophosphoryl lipid A stimulation of TLR4. J Immunol. 2011;186:3858–65.View ArticlePubMed CentralPubMedGoogle Scholar
- Doherty R, O'Farrelly C, Meade KG. Epigenetic regulation of the innate immune response to LPS in bovine peripheral blood mononuclear cells (PBMC). Vet Immunol Immunopathol. 2013;154:102–10.View ArticlePubMedGoogle Scholar
- Krausgruber T, Blazek K, Smallie T, Alzabin S, Lockstone H, Sahgal N, et al. IRF5 promotes inflammatory macrophage polarization and TH1-TH17 responses. Nat Immunol. 2011;12:231–8.View ArticlePubMedGoogle Scholar
- Green BB, Kerr DE. Epigenetic contribution to individual variation in response to lipopolysaccharide in bovine dermal fibroblasts. Vet Immunol Immunopathol. 2014;157:49–58.View ArticlePubMed CentralPubMedGoogle Scholar
- Lin G, Pearson AE, Scamurra RW, Zhou Y, Baarsch MJ, Weiss DJ, et al. Regulation of interleukin-8 expression in porcine alveolar macrophages by bacterial lipopolysaccharide. J Biol Chem. 1994;269:77–85.PubMedGoogle Scholar
- Schmeck B, Beermann W, van Laak V, Zahlten J, Opitz B, Witzenrath M, et al. Intracellular bacteria differentially regulated endothelial cytokine release by MAPK-dependent histone modification. J Immunol. 2005;175:2843–50.View ArticlePubMedGoogle Scholar
- Raymond B, Batsche E, Boutillon F, Wu YZ, Leduc D, Balloy V, et al. Anthrax lethal toxin impairs IL-8 expression in epithelial cells through inhibition of histone H3 modification. PLoS Pathog. 2009;5:e1000359.View ArticlePubMed CentralPubMedGoogle Scholar
- Chadebech P, Goidin D, Jacquet C, Viac J, Schmitt D, Staquet MJ. Use of human reconstructed epidermis to analyze the regulation of beta-defensin hBD-1, hBD-2, and hBD-3 expression in response to LPS. Cell BiolToxicol. 2003;19:313–24.Google Scholar
- Chen QX, Lv C, Huang LX, Cheng BL, Xie GH, Wu SJ, et al. Genomic variations within DEFB1 are associated with the susceptibility to and the fatal outcome ofsevere sepsis in Chinese Han population. Genes Immun. 2007;8:439–43.View ArticlePubMedGoogle Scholar
- Shuyi Y, Feng W, Jing T, Hongzhang H, Haiyan W, Pingping M, et al. Human beta-defensin-3 (hBD-3) upregulated by LPS via epidermal growth factor receptor (EGFR) signaling pathways to enhance lymphatic invasion of oral squamous cell carcinoma. Oral Surg Oral Med Oral Pathol Oral Radiol Endod. 2011;112:616–25.View ArticlePubMedGoogle Scholar
- Diamond DL, Kimball JR, Krisanaprakornkit S, Ganz T, Dale BA. Detection of beta-defensins secreted by human oral epithelial cells. J Immunol Methods. 2001;256:65–76.View ArticlePubMedGoogle Scholar
- Alexander S, Bramson J, Foley R, Xing Z. Protection from endotoxemia by adenoviral-mediated gene transfer of human bactericidal/permeability-increasing protein. Blood. 2004;103:93–9.View ArticlePubMedGoogle Scholar
- Balakrishnan A, Marathe SA, Joglekar M, Chakravortty D. Bactericidal/permeability increasing protein: a multifaceted protein with functions beyond LPS neutralization. Innate Immun. 2013;19:339–47.View ArticlePubMedGoogle Scholar
- Islam D, Bandholtz L, Nilsson J, Wigzell H, Christensson B, Agerberth B, et al. Downregulation of bactericidal peptides in enteric infections: a novel immune escape mechanism with bacterial DNA as a potential regulator. Nat Med. 2001;7:180–5.View ArticlePubMedGoogle Scholar
- Zaalouk TK, Bajaj-Elliott M, George JT, McDonald V. Differential regulation of beta-defensin gene expression during Cryptosporidium parvum infection. Infect Immun. 2004;72:2772–9.View ArticlePubMed CentralPubMedGoogle Scholar
- Sperandio B, Regnault B, Guo J, Zhang Z, Stanley Jr SL, Sansonetti PJ, et al. Virulent Shigella flexneri subverts the host innate immune response through manipulation of antimicrobial peptide gene expression. J Exp Med. 2008;205:1121–32.View ArticlePubMed CentralPubMedGoogle Scholar
- Obici L, Merlini G. Amyloidosis in autoinflammatory syndromes. Autoimmun Rev. 2012;12:14–7.View ArticlePubMedGoogle Scholar
- Shah C, Hari-Dass R, Raynes JG. Serum amyloid A is an innate immune opsonin for Gram-negative bacteria. Blood. 2006;108:1751–7.View ArticlePubMedGoogle Scholar
- Christenson K, Björkman L, Tängemo C, Bylund J. Serum amyloid A inhibits apoptosis of human neutrophils via a P2X7-sensitive pathway independent of formyl peptide receptor-like 1. J Leukoc Biol. 2008;83:139–48.View ArticlePubMedGoogle Scholar
- De Santo C, Arscott R, Booth S, Karydis I, Jones M, Asher R, et al. Invariant NKT cells modulate the suppressive activity of IL-10-secreting neutrophils differentiated with serum amyloid A. Nat Immunol. 2010;11:1039–46.View ArticlePubMed CentralPubMedGoogle Scholar
- Nguyen KD, Macaubas C, Truong P, Wang N, Hou T, Yoon T, et al. Serum amyloid A induces mitogenic signals in regulatory T cells via monocyte activation. Mol Immunol. 2014;59:172–9.View ArticlePubMedGoogle Scholar
- Gross SR, Sin CG, Barraclough R, Rudland PS. Joining S100 proteins and migration: for better or for worse, in sickness and in health. Cell Mol Life Sci. 2014;71:1551–79.View ArticlePubMedGoogle Scholar
- Barraclough R, Savin J, Dube SK, Rudland PS. Molecular cloning and sequence of the gene for p9Ka. A cultured myoepithelial cell protein with strong homology to S-100, a calcium-binding protein. J Mol Biol. 1987;198:13–20.View ArticlePubMedGoogle Scholar
- Lerner KG, Glick B, McDuffie FC. Role of the bursa of Fabricius in IgG and IgM production in the chicken: evidence for the role of a non-bursal site in the development of humoral immunity. J Immunol. 1971;107:493–503.PubMedGoogle Scholar
- Rivers RP, Hathaway WE, Weston WL. The endotoxin-induced coagulant activity of human monocytes. Br J Haematol. 1975;30:311–6.View ArticlePubMedGoogle Scholar
- Levy GA, Leibowitz JL, Edgington TS. Induction of monocyte procoagulant activity by murine hepatitis virus type 3 parallels disease susceptibility in mice. J Exp Med. 1981;154:1150–63.View ArticlePubMedGoogle Scholar
- Maier RV, Ulevitch RJ. The induction of a unique procoagulant activity in rabbit hepatic macrophages by bacterial lipopolysaccharides. J Immunol. 1981;127:1596–600.PubMedGoogle Scholar
- Chung YM, Goyette J, Tedla N, Hsu K, Geczy CL. S100A12 suppresses pro-inflammatory, but not pro-thrombotic functions of serum amyloid A. PLoS One. 2013;8:e62372.View ArticlePubMed CentralPubMedGoogle Scholar
- Piriou-Guzylack L, Salmon H. Membrane markers of the immune cells in swine: an update. Vet Res. 2008;39:54.View ArticlePubMedGoogle Scholar
- Ledger TN, Pinton P, Bourges D, Roumi P, Salmon H, Oswald IP. Development of a macroarray to specifically analyze immunological gene expression in swine. Clin Diagn Lab Immunol. 2004;11:691–8.PubMed CentralPubMedGoogle Scholar
- Lehtonen A, Matikainen S, Miettinen M, Julkunen I. Granulocyte-macrophage colony-stimulating factor (GM-CSF)-induced STAT5 activation and target-gene expression during human monocyte/macrophage differentiation. J Leukoc Biol. 2002;71:511–9.PubMedGoogle Scholar
- Elssner A, Doseff AI, Duncan M, Kotur M, Wewers MD. IL-16 is constitutively present in peripheral blood monocytes and spontaneously released during apoptosis. J Immunol. 2004;172:7721–5.View ArticlePubMedGoogle Scholar
- Benoit M, Fenollar F, Raoult D, Mege JL. Increased levels of circulating IL-16 and apoptosis markers are related to the activity of Whipple's disease. PLoS One. 2007;2:e494.View ArticlePubMed CentralPubMedGoogle Scholar
- Bowler RP, Bahr TM, Hughes G, Lutz S, Kim YI, Coldren CD, et al. Integrative omics approach identifies interleukin-16 as a biomarker of emphysema. OMICS. 2013;17:619–26.View ArticlePubMed CentralPubMedGoogle Scholar
- Kobayashi SD, Voyich JM, Braughton KR, DeLeo FR. Down-regulation of proinflammatory capacity during apoptosis in human polymorphonuclear leukocytes. J Immunol. 2003;170:3357–68.View ArticlePubMedGoogle Scholar
- Feske S, Giltnane J, Dolmetsch R, Staudt LM, Rao A. Gene regulation mediated by calcium signals in T lymphocytes. Nat Immunol. 2001;2:316–24.View ArticlePubMedGoogle Scholar
- Wilson KC, Center DM, Cruikshank W. The effect of interleukin-16 and its precursor on T lymphocyte activation and growth. Growth Factors. 2004;22:97–104.View ArticlePubMedGoogle Scholar
- Atanackovic D, Hildebrandt Y, Templin J, Cao Y, Keller C, Panse J, et al. Role of interleukin 16 in multiple myeloma. J Natl Cancer Inst. 2012;104:1005–20.View ArticlePubMedGoogle Scholar
- Fontanesi L, Martelli PL, Scotti E, Russo V, Rogel-Gaillard C, Casadio R, et al. Exploring copy number variation in the rabbit (Oryctolagus cuniculus) genome by array comparative genome hybridization. Genomics. 2012;100:245–51.View ArticlePubMedGoogle Scholar
- Smyth GK, Michaud J, Scott HS. Use of within-array replicate spots for assessing differential expression in microarray experiments. Bioinformatics. 2005;21:2067–75.View ArticlePubMedGoogle Scholar
- Morey JS, Ryan JC, Van Dolah FM. Microarray validation: factors influencing correlation between oligonucleotide microarrays and real-time PCR. Biol Proced Online. 2006;8:175–93.View ArticlePubMed CentralPubMedGoogle Scholar
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.