Skip to main content


Transcription analysis on response of porcine alveolar macrophages to Haemophilus parasuis



Haemophilus parasuis (H. parasuis) is the etiological agent of Glässer's disease in pigs. Currently, the molecular basis of this infection is largely unknown. The innate immune response is the first line of defense against the infectious disease. Systematical analysis on host innate immune response to the infection is important for understanding the pathogenesis of the infectious microorganisms.


A total of 428 differentially expressed (DE) genes were identified in the porcine alveolar macrophages (PAMs) 6 days after H. parasuis infection. These genes were principally related to inflammatory response, immune response, microtubule polymerization, regulation of transcript and signal transduction. Through the pathway analysis, the significant pathways mainly concerned with cell adhesion molecules, cytokine-cytokine receptor interaction, complement and coagulation cascades, toll-like receptor signaling pathway, MAPK signaling pathway, suggesting that the host took different strategies to activate immune and inflammatory response upon H. parasuis infection. The global interactions network and two subnetworks of the proteins encoded by DE genes were analyzed by using STRING. Further immunostimulation analysis indicated that mRNA levels of S100 calcium-binding protein A4 (S100A4) and S100 calcium-binding protein A6 (S100A6) in porcine PK-15 cells increased within 48 h and were sustained after administration of lipopolysaccharide (LPS) and Poly (I:C) respectively. The s100a4 and s100a6 genes were found to be up-regulated significantly in lungs, spleen and lymph nodes in H. parasuis infected pigs. We firstly cloned and sequenced the porcine coronin1a gene. Phylogenetic analysis showed that poCORONIN 1A belonged to the group containing the Bos taurus sequence. Structural analysis indicated that the poCORONIN 1A contained putative domains of Trp-Asp (WD) repeats signature, Trp-Asp (WD) repeats profile and Trp-Asp (WD) repeats circular profile at the N-terminus.


Our present study is the first one focusing on the response of porcine alveolar macrophages to H. parasuis. Our data demonstrate a series of genes are activated upon H. parasuis infection. The observed gene expression profile could help screening the potential host agents for reducing the prevalence of H. parasuis and further understanding the molecular pathogenesis associated with H. parasuis infection in pigs.


The pig is an important agricultural animal and is an excellent mammalian model for biomedical research [1, 2]. H. parasuis is the etiological agent of porcine polyserositis and arthritis (Glässer's disease) characterized by fibrinous polyserositis, meningitis and polyarthritis, causing severe economic losses to the swine industry [3]. To date, 15 serovars of H. parasuis have been identified [4]. H. parasuis infection can be acute or chronic, depending on the immunological status of the herd [3]. The infection by H. parasuis has become an increasing threat to early-weaned pigs and in pig herds of high health status [5, 6].

The innate immune response in vertebrates is the first defense line against invading microorganisms. The main players in innate immunity are phagocytes such as neutrophils, dendritic cells and macrophages [7]. As a major component of the host innate immunity, macrophages have essential roles in host defense to infection, because they often mediate the killing of microbes as well as initiate, maintain and resolve host inflammatory responses by releasing cytokines and chemokines [811]. Bacterial pathogens that overcome host defenses ensure their ability to survive and propagate [10, 12]. The diversity of bacteria and the differences in their pathogenesis may lead to pathogen-specific responses of macrophages [10]. A greater understanding of the complex interactions, which occur between the macrophages and pathogen, could lead to the identification of the host defense strategies and the complementary pathogen evasion strategies [8, 10]. The interactions between H. parasuis with porcine alveolar macrophages have been studied [13], but the detailed mechanisms of how porcine alveolar macrophages response to H. parasuis infection are not well elucidated. The high throughput cDNA microarray represents a powerful tool for analyzing the molecular events in bacteria-host cell interactions [14]. This technology has been useful in identifying changes in gene expression both in cultured cells and in whole organisms infected with pathogens [12, 15, 16]. In this study, we applied this high throughput cDNA microarray assay to improve our understanding of the innate immune response of macrophages to H. parasuis infection.


Clinical evaluation of infected pigs

In the challenge group, at 144 h post-infection, all three pigs had a rectal temperature of over 40.6°C and displayed lateral recumbency or labored breathing. At 144 h post-infection, all pigs including three control pigs were euthanized by intravenous administration of an overdose of sodium pentobarbital. In the challenge group, severe fibrinous polyserositis, arthritis and meningitis were observed at necropsy. On the other hand, in the control group, three pigs remained clinically normal throughout the experiment and did not have lesions at necropsy. The detection of H. parasuis by bacterial isolation, nested PCR and LAMP in different samples are shown in Table 1. The results indicated that the H. parasuis could be detected in the lymph nodes, lungs and spleen in all of the three pigs that challenged with H. parasuis serovar 5 SH0165 strain. In contrast, H. parasuis could not be detected by the three methods in the control group.

Table 1 Results of culture, nested PCR and LAMP analysis for three pigs challenged with H. parasuis serovar 5, presented as the number of pigs positive/pigs samples

Overview of differential expressed genes in PAMs

To investigate the dynamic gene transcriptional profiles of PAMs in response to H. parasuis infection, six microarrays were used in this experiment, corresponding to the RNAs from PAMs of three H. parasuis infected piglets and three controls. The total RNA samples were hybridized with Affymetrix GeneChip Porcine Genome Array, and the microarray data were analyzed using Significance Analysis of Microarrays (SAM) [17]. Hybridization results indicated that 14,228 and 13,813 probes sets, corresponding to 58.9% and 57.3% of all probe sets, were detected in H. parasuis serovar 5 and mock-infected PAMs (Additional file 1). After quantile normalization and statistical analysis, 623 transcripts were identified at SAM |Score(d)| ≥ 2. Furthermore, Genes whose relative transcription levels showed a fold change FC ≥ 1.33 and SAM |Score(d)| ≥ 2 were considered to be up-regulated, and those with an FC ≤ 0.75 and SAM |Score(d)| ≥ 2 were considered to be down-regulated. In this study, 575 transcripts showed a level of expression that differed significantly from that of the control group with H. parasuis serovar 5 SH0165 strain infected group. These included 428 genes annotated with DAVID or by searching against the GenBank database (Additional file 2). Among these, 338 genes were up-regulated and 90 genes were down-regulated (Figure 1).

Figure 1

A hierarchical cluster of 428 transcripts following H. parasuis infection in PAMs. Each row represents a separate transcript and each column represents a separate piglet. Color legend is on the left, red indicates increased transcript expression levels, whereas green indicates decreased levels compared with normal samples. Up-/down-regulated response transcripts are highlighted on the right with the No. in each cluster in parentheses. Fold change is calculated based on the mean intensity value from 3 donors within each group.

The functions of the DE genes were analyzed by Molecular Annotation System 3.0 software [18]. In the MAS 3.0 tool, the GO terms and KEGG pathways are ranked with statistical significance by calculating their p-values based on hypergeometric distribution [19]. Go terms and KEGG pathways with p- values less than 0.05 are considered statistically significant [20]. Among the 428 annotated genes, a total of 229 genes were grouped into 156 categories based on biological process Gene Ontology (GO) terms with the p-values less than 0.05 (Additional file 3). Several GO terms were associated with the immune system. These were "inflammatory response" (GO:0006954, p = 2.06E-04), "immune response" (GO:0006955, p = 1.75E-03), "complement activation, classical pathway" (GO:0006958, p = 9.48E-04) and "leukocyte migration" (GO:0050900, p = 9.55E-03). Particularly, the DE genes associated with immune and inflammatory response suggested that they play roles in host defense response to H. parasuis infection (Table 2). To gain insight into the different biological processes of H. parasuis infection, a pathway analysis by KEGG database was performed on DE genes (Additional file 4) [18]. The significant pathways mainly contained: cell adhesion molecules (p = 3.74E-12), cytokine-cytokine receptor interaction (p = 2.18E-10), complement and coagulation cascades (p = 1.16E-07), toll-like receptor signaling pathway (p = 1.93E-05), MAPK signaling pathway (p = 7.63E-04), which suggested that the host took different strategies to activate immune and inflammatory response upon H. parasuis infection.

Table 2 The DE genes associated with immune and inflammatory response in PAM after H. parasuis serovar 5 infection 6 days

Validation of microarray data by quantitative real-time PCR (qPCR)

In order to confirm the statistical significance of our findings, we performed quantitative real-time PCR (qPCR) analysis of the relevant genes in our original samples used in microarray study. Eleven genes were selected for qPCR analysis. Ten selected genes that were up-regulated in microarray also showed significantly higher expression in H. parasuis serovar 5 infected samples than in the control samples determined by qPCR analysis. The ppp1r13l gene that was down-regulated in microarray data also showed significantly lower expression in H. parasuis serovar 5 infected samples than in the control samples by qPCR (Table 3).

Table 3 Validation of microarray results by qPCR

STRING analysis of the relationships between DE genes

STRING is a web-based interface that could predict protein associations which can mean direct physical binding and can also mean indirect interaction such as participation in the same metabolic pathway or cellular process on the basis of genomic context, high-throughput experiments, co-expression and data from the literature [21, 22]. DE genes were analyzed using STRING for predicting network of proteins encoded by DE genes. Among the 428 annotated DE genes, 236 genes containing 181 up-regulated genes and 55 down-regulated genes were eligible to STRING analyses when the Sus Scrofa database was chosen. In order to seek the possibility of the associations between DE genes, the combined score of 0.15 was chosen. The network of predicted associations for all of the DE genes encoded proteins are shown in Additional file 5. Some molecules are the key molecules that link to other proteins according to the STRING analysis. However, many proteins do not link to others, indicating that their functions are unrelated or unknown. As shown in Figure 2A, a total of 12 DE genes encoded proteins are associated with IL-1β according to the textmining evidence, and they form the IL-1β network. Furthermore, the CD14 and SOD2 are associated with IL-1β according to the co-expression evidence. Two of phagocytosis-related genes (cd14, fcgr2β) are associated with a total of 11 DE genes encoded proteins (Figure 2B) according to the textmining evidence, and they form the phagocytosis network. Furthermore, the IL-1β, FGL-2, CCL2 and FCGR2β are associated with CD14 according to the co-expression evidence.

Figure 2

STRING analysis of the relationship between DE genes. The DE genes were analyzed using the Sus Scrofa STRING database. The network nodes represent the proteins encoded by the DE genes. Seven different colored link a number of nodes and represent seven types of evidence used in predicting associations. A red line indicates the presence of fusion evidence; a green line represents neighborhood evidence, a blue line represents coocurrence evidence; a purple line represents experimental evidence; a yellow line represents textmining evidence; a light blue line represents database evidence and a black line represents coexpression evidence. A: The network of DE genes related to IL-1β. B: The network of DE genes related to phagocytosis.

Identification of novel infection-related DE genes

In order to identify novel candidates for disease-related DE genes, we evaluated the DE genes that were not highlighted in the KEGG or STRING analysis. Although the s100a4, s100a6 and coronin 1a were not highlighted in the KEGG or STRING analysis, they were found to play roles in the immune response [2334]. These observations suggest that the three genes may be novel candidates for disease-related DE genes.

Among the three genes, coronin 1a has not been identified in pigs before. So we cloned and sequenced the porcine coronin 1a gene according to description of Liu et al [35] and the sequence was submitted to the GenBank [GenBank: JN092377]. The full-length cDNA of porcine coronin 1a contains 1386 bp and 461 amino acid residues. Multiple sequence alignment with the identified coronin 1a of cattle, human, mouse, rat and the predicted coronin 1a of other species showed that the nucleotide sequence of the poCORONIN 1A ORF is 93.58%, 93.58%, 92.42%, 92.35%, 92.28%, 91.70%, 87.81% and 86.87% identical to that of panda, cattle, human, chimpanzee, northern white-cheeked gibbon, common marmoset, mouse and rat coronin 1a, respectively. At the amino acid level, the corresponding identities were 96.96%, 97.61%, 95.66%, 95.66%, 95.88%, 96.10%, 93.71% and 92.84%, respectively. To define the molecular evolutionary history of poCORONIN 1A, protein sequences from 9 vertebrates were obtained to construct a phylogenetic tree. Phylogenetic analysis showed that poCORONIN 1A belongs to the group containing the Bos taurus sequence (Additional file 6). Structural analysis with the ExPASy server indicated that the poCORONIN 1A contains putative domains of Trp-Asp (WD) repeats signature, Trp-Asp (WD) repeats profile and Trp-Asp (WD) repeats circular profile at the N-terminus (Additional file 7).

Expression analyses of S100A4, S100A6 in PK15 cells stimulated with LPS and Poly (I:C)

In order to investigate the expression patterns of s100a4 and s100a6 under general conditions that mimic bacterial and viral infection, the immunostimulation assay was carried out in PK-15 cells by using the LPS and Poly (I:C) as the stimulators.

Overnight cultures of PK-15 cells were treated with 1 μg/ml LPS or 10 μg/ml Poly (I:C) for 0, 2, 6, 12, 24 and 48 h. LPS and Poly (I:C) stimulation did not induce expression of porcine s100a4 until 48 h (Figure 3A, C). LPS stimulation induced expression of s100a6 at 2 h and 12 h, after which s100a6 expression dropped and plateaued for 24-48 h (Figure 3B). After Poly (I:C) stimulation, the expression of s100a6 reached the peak at 12 h, after which s100a6 expression dropped at 24 h, and the up-regulation of s100a6 was again observed at 48 h (Figure 3D). These observations indicate that both LPS and Poly (I:C) can induce the expression of porcine s100a4 and s100a6 in vitro.

Figure 3

Kinetic immune stimuli analyses challenged by LPS and Poly (I:C) in PK-15 cells. A, B: LPS-induced expression of porcine s100a4 and s100a6 in PK-15 cells respectively. PK-15 cells were cultured with 1 μg/mL LPS for 48 h. C, D: Poly (I:C) induced expression of porcine s100a4 and s100a6 in PK-15 cells respectively. PK-15 cells were cultured with 10 μg/mL Poly (I:C) for 48 h. Relative expression of s100a4 and s100a6 were detected by qPCR and normalized to the expression of GAPDH. The fold increase is expressed as the mean of three replicates with SEM by comparison with the control (0 h). QPCR was performed using primers described in Table 3. The significance of difference for the expression compared to the untreated control (0 h) was calculated using two-directional paired Student's T-test. ** p ≤ 0.01; * p ≤ 0.05.

In vivo expression of s100a4 and s100a6 in pigs with systemic infection of H. parasuis

In order to understand the expression of the s100a4 and s100a6 in pigs with systemic infection of H. parasuis, the different tissues obtained from the H. parasuis infected pigs and the controls were selected for the qPCR analysis. Our qPCR examination demonstrated that the increasing expression of s100a4 was observed in the lungs, spleen and lymph nodes of pigs infected with H. parasuis for 6 days (Figure 4A). The expression of s100a6 in the lungs, spleen and lymph nodes had the same expression tendencies (Figure 4B). However, in brain and heart of H. parasuis infected pigs, the expression of s100a4 and s100a6 did not show significant changes compared to the controls.

Figure 4

Quantitative expression of s100a4 and s100a6 in five tissues from pigs with Glässer's disease. A: Increased in vivo gene expression of s100a4 in lungs, spleen, lymph nodes of pigs with Glässer's disease. B: Increased in vivo gene expression of s100a6 in lungs, spleen, lymph nodes of pigs with Glässer's disease. Relative expression of s100a4 and s100a6 were detected by qPCR and normalized to the expression of GAPDH. The fold increase is expressed as the mean of three replicates with SEM by comparison with the control. The significance of difference for the expression compared to the control was calculated using two-directional paired Student's T-test. ** p ≤ 0.01; * p ≤ 0.05.


During infection, H. parasuis has to reach the lung and survive the host pulmonary defenses before invading the blood stream [13]. In the lung, bacteria have to confront alveolar macrophages, whose main roles include: ingestion of bacteria by phagocytosis, destruction of bacteria within phagolysosomes and recruitment of inflammatory cells to the site of infection via chemokines and acute-phase proteins [36].

Phagocytosis is a cytoskeleton-dependent process of engulfment of large particles, and macrophages could present a restricted number of phagocytic receptors that induce rearrangements in the actin cytoskeleton that lead to the internalization of the particle [36]. Phagocytosis is a key mechanism used by macrophages to control virulent Pasteurellaceae, such as Pasteurella multocida, Haemophilus parasuis, Haemophilus influenzae, Actinobacillus pleuropneumoniae [13, 3740]. In this study, the cd14 [4144], hmox1 [45], fcgr2β [46] and abca1 [47] genes, which were identified as DE genes, were also found to be involved in the phagocytosis. Meanwhile, the STRING analysis indicated that many DE gene encoded proteins could interact with CD14 and FCGR2β molecules, suggesting that PAMs may upregulate these genes to facilitate the phagocytosis of H. parasuis or other cells to play their immunological roles.

One of the important mechanisms used by macrophage to play its immunological functions is to kill bacteria by the activation and recruitment of antibacterial effectors to the phagolysosome [36]. The fusion of phagosomes with lysosomes results in the formation of phagolysosomes [11, 36]. In our study, we found two DE genes that were related to the formation of phagolysosome, that is, smpd1 [48] and coronin 1a [2931]. Interestingly, many groups have reported that the CORONIN 1A could prevent lysosomal delivery and allow the bacteria to survive intracellularly [30, 3234]. In the course of infection, H. parasuis has to survive from the host pulmonary defense, such as alveolar macrophages, to produce disease. In this way, the up-regulation of coronin 1a gene may facilitate the H. parasuis in producing the disease.

Interleukin-1 beta (IL-1β) is an important inflammation-associated gene that is up-regulated in many microarray experiments [11, 14, 16, 49]. Interestingly, Wilkinson et al reported that an increase in IL-1β gene expression is observed in H. parasuis- infected lungs [50]. In our study, IL-1β was also up-regulated in H. parasuis- infected PAMs. Unsurprisingly, STRING analysis also revealed that many molecules encoded by up-regulated genes interact with IL-1β and form the IL-1β network. Meanwhile, the pathway analysis indicated that IL-1β is in some pathways, such as cytokine-cytokine receptor interaction (p = 2.18E-10), MAPK signaling pathway (p = 7.63E-04), and toll-like receptor signaling pathway (p = 1.93E-05). CCL5/RANTES plays an important role in regulating the movements of inflammatory cells to the infection sites [51, 52]. Many viruses, such as Japanese encephalitis virus (JEV) [53], respiratory syncytial virus (RSV) [54], influenza virus A [55] and porcine reproductive and respiratory syndrome virus (PRRSV) [51] have been shown to induce CCL5. In addition, some papers have reported that the CCL5 could be induced in macrophages by bacterial infections, such as Salmonella typhimurium infection, Streptococcus pyogenes infection and Lactobacillus rhamnosus infection [56, 57]. Interestingly, in our study, the up-regulation of CCL5 was observed in H. parasuis-infected PAMs, which suggested that CCL5 plays a role in the host response against H. parasuis infection. Thus, during the H. parasuis infection, the PAMs mount a powerful inflammatory response in an effort to clear this pathogen. Alternatively, the influx of inflammatory cells to the site of infection may provide additional host cells for H. parasuis infection. However, sustained or excessive production of inflammatory cytokines can have damaging consequences. To counterbalance inflammatory cytokines, anti-inflammatory cytokines are produced. Anti-inflammatory cytokines include interleukin 10 (IL-10), transforming growth factor β (TGF-β), and IL-1 receptor antagonist (IL-1RA) [11, 58, 59]. Wilkinson et al reported that the IL-1β and its antagonist, IL-1RA are both more highly expressed in "susceptible" animals challenged with H. parasuis [50]. In our study, TGF-β, an anti-inflammatory cytokine, was increased in H. parasuis infection group. During H. parasuis infection, anti-inflammatory signals may decrease the potentially damaging effects of proinflammatory cytokines on host tissue.

Macrophage also effectively controls bacterial infection by producing of reactive species such as oxygen species and nitric oxide (NO). Sustained production of NO endows macrophages with cytostatic or cytotoxic activity against viruses, bacteria, fungi, protozoa, helminths and tumor cells. Unsurprisingly, H. parasuis infection could cause up-regulated expression of a large set of genes involved in the nitric oxide production. These genes were: spr, rora, klrk1, sod2 and il-1β [6066]. The up-regulated genes related to the nitric oxide production may contribute to the PAM for confronting H. parasuis infection.

The DE genes that are related to phagocytosis, formation of phagolysosome, chemokines production, and nitric oxide production may help us to better understand the complicated mechanisms by which PAMs play their functions. Another highlight of our study is the new identified candidate genes that may be implicated in the pathogenesis of Glässer's disease. These genes could help to screen the potential host agents for reducing the prevalence of H. parasuis and further understand the molecular pathogenesis associated with H. parasuis infection in pigs. These genes are: s100a4, s100a6, caveolin 2 and ppp1r13l.

S100 A4 and S100 A6 belong to the S100 family that contained 2 EF-hand calcium-binding motifs [23, 27]. Two of S100 family genes (S100 calcium-binding protein A8 and A9) are dramatically up-regulated in spleen and lung following H. parasuis infection [67]. Meanwhile, many other S100 family genes are up-regulated following different bacterial and viral infection [16, 18, 67, 68], suggesting that the S100 family genes play roles in the immune response to infections. In our study, the S100 calcium-binding protein A4 and A6 were up-regulated after H. parasuis infection when determined by microarray and qPCR. Further immunostimulation analysis indicated that the mRNA levels of S100 calcium-binding protein A4 (S100A4) and S100 calcium-binding protein A6 (S100A6) in porcine PK-15 cells increased within 48 h and were sustained after administration of LPS and Poly (I:C) respectively. We also found that the s100a4 and s100a6 genes were up-regulated in lungs, spleen and lymph nodes in H. parasuis infected pigs. Interestingly, the kidney fibrosis [2426, 69], liver fibrosis [70, 71], lung fibrosis [7274], cardiac fibrosis [23, 75, 76] and peritoneal fibrosis [77] are found to be related to the expression of s100a4. Glässer's disease is characterized mainly by fibrinous polyserositis, meningitis, and arthritis. In this way, we hypothesized that the increase expression of s100a4 may underlie fibrosis during H. parasuis infection in pigs. Meanwhile, some reports indicated that s100a6 plays roles in cell proliferation and signaling transduction [27, 28]. Therefore, the s100a4 and s100a6 genes could be two novel genes related to H. parasuis infection.

Caveolins are the major components and protein markers of caveolae that are 50-100 nm invaginations of membrane. The caveolin gene family includes three members in vertebrates, caveolin-1, caveolin-2 and caveolin-3, of which caveolin 1 and caveolin 2 have been detected in mouse macrophages [35, 78]. Caveolin 1 molecule is related to H. parasuis infection [35]. Caveolin 2, which localizes to the Golgi complex but redistributes to plasma membrane, caveolae and rafts when co-expressed with caveolin 1, is a potential key molecule related to the Pseudomonas infection causing pneumonia in patients with cystic fibrosis and other immunocompromising conditions [79, 80]. In our study, the caveolin 2 gene was highly expressed in PAM isolated from the H. parasuis serovar 5 challenged group. Therefore, in addition to the caveolin 1 gene, the caveolin 2 gene may be a novel candidate gene related to H. parasuis infection.

The NF-kappa B (NF-κB) signaling pathway is important in signal transduction during the innate immune response [36]. NF-κB signaling relies on the targeting of IκB (inhibitor of NF-κB) subunit to the proteasome to allow NF-κB to translocate from the cytosol to the nucleus where it activates gene transcription [81]. The PPP1R13L is mentioned as a novel inhibitor of NF-κB [82]. In our study, microarray and qPCR analysis indicated that the mRNA of PPP1R13L was down-regulated significantly compared to control. The IPA network indicated that the PPP1R13L could directly or indirectly interacts with many molecules, such as micro RNAs, transcriptions, enzymes, and cytokines (Additional file 8), suggesting that ppp1r13l gene is an innate immune related gene that plays a role in PAM during H. parasuis infection. The detailed mechanism of ppp1r13l gene in NF-κB signaling pathway in H. parasuis infected PAM needs further studies.


This is the first study focusing on response of porcine alveolar macrophages to Haemophilus parasuis by using the Affymetrix GeneChip Porcine Genome Array. Although great efforts have been made to understand the molecular basis of H. parasuis infection, the cellular response to H. parasuis infection is still largely unknown. The high-density cDNA array technology to analysis of H. parasuis-infected PAM could improve our understanding of the H. parasuis infection. Our data show that a series of genes are activated upon H. parasuis infection. These genes are involved in inflammatory response, immune response, microtubule polymerization, regulation of transcript and signal transduction. Particularly, some genes related to phagocytosis, formation of phagolysosome, chemokines production and nitric oxide production could contribute to explain the complicated mechanisms by which PAM played its functions. Some new identified genes may also provide implication on the pathogenesis of Glässer's disease caused by H. parasuis.


Animals for Microarray experiment and porcine alveolar macrophages isolation

All animals' tissue collection procedures were performed according to protocols approved by the Hubei Province PR China for Biological Studies Animal Care and Use Committee. Six piglets which were obtained from a commercial herd free of Glässer's disease were weaned at 27 days, shipped to the Animal Disease Center of Huazhong Agricultural University, and raised with isolation facilities. Three piglets were randomly allocated to the non infected group and three to the infected group. The three piglets were intratracheally challenged with H. parasuis strain 0165 (serovar 5) at a dose of 6 × 109 colony-forming units (CFU). The noninfected group piglets were treated similarly with identical volume of PBS served as control. All piglets were determined to the HPS-free by serum indirect haemagglutination (IHA) test before artificial bacterial challenges. Clinical signs and lesions of Glässer's disease were apparent in the challenged group at 6 days post-infection (dpi). All piglets were slaughtered at 6 dpi. Bacterial isolation, nested PCR and LAMP were performed after the piglets were killed at 6 dpi. PAMs were isolated according to Olvera's description [13]. Briefly, Bronchoalveolar lavage of the lungs was performed with 100 mL aliquots of sterile PBS containing gentamicin at 70 μg/mL (Sigma-Aldrich). To collect the porcine alveolar macrophages (PAM), lavage fluids were centrifuged at 230 g for 15 min, and then cells were washed twice with Dulbecco's Modified Eagle's Medium (DMEM) with gentamicin (50 μg/mL). PAM isolation was confirmed by detection of macrophage markers (SWC3, CD169 and SLAII) in the cells by flow cytometry.

RNA preparation for Microarray experiment

Total RNA were extracted from PAM of each group with Trizol (Invitrogen) then quantified using the NanoDrop 1000 Spectrophotometer (Thermo Fisher Scientific Inc., USA). The quality of the RNA was checked by formaldehyde denaturing gel electrophoresis in 1.2% agarose gels, which showed dispersed bands (28S and 18S) without any obvious smearing patterns that would indicate degradation.

Microarray hybridization and data analyses

Affymetrix GeneChip Porcine Genome Array, which contains 24,123 probe sets to interrogate 23,256 transcripts in pig, represents 20,201 genes, was used in microarray analysis. Hybridization, data capture and analysis were performed by CapitalBio Corporation (Beijing, China), a service provider authorized by Affymetrix Inc. (Santa Clara, CA). Briefly, a total of 1 μg RNA was used for cDNA synthesis and to produce biotin-tagged cRNA with GeneChip IVT Labeling kit (Affymetrix). A total of 15 μg fragmented cRNA, with contol oligo B2 and eukaryotic hybridization controls (bioB, bioC, bioD, cre) was hybridized to each GeneChip array at 45°C for 16 hours (Affymetrix Gene Chip Hybridization Oven 640) according to manufacturer's instructions. After hybridization, the GeneChip arrays were washed and stained with streptavidin phycoerythrin onan (SAPE) with Affymetrix Fluidics Station 450 followed by scanning with the Affymetrix GeneChip Scanner 3000. Six microarrays were used in the experiment, corresponding to the RNAs from PAMs of three H. parasuis infected piglets and three controls.

The hybridization data were analyzed using GeneChip Operating Software (GCOS, version 1.4), which uses statistical criteria to generate a 'present' or 'absent' call for genes represented by each probe set on the array. The scanned images were first assessed by visual inspection and then analyzed to generate raw data files saved as CEL files using the default setting of GCOS 1.4. Microarray data were normalized using the robust multi-array average (RMA) method [83], which consists of three steps: background correction, quantile normalization (each performed at the individual probe level), and robust linear model fit using log-transformed intensities (at the probe set level). Significance Analysis of Microarrays (SAM) add-in to Microsoft Excel was used for comparisons of replicate array experiments. SAM identifies genes with statistically significant changes in expression by assimilating a set of gene-specific t-tests, and provides an estimate of the false discovery rate (FDR) from randomly generated data. Genes with scores higher than a threshold value or genes with FDR value lower than the threshold value were deemed potentially significant. Furthermore, fold-change analysis which calculates the ratios of geometric means of expression intensities of H. parasuis-infected PAMs relative to controls was performed. These ratios were reported as the up-or down-fold change. In this study, genes were considered statistically significant if they had SAM |Score(d)| ≥ 2 [84, 85] and exhibited a fold change ≥ 1.33 and ≤ 0.75. DE genes performed for hierarchical cluster (Ver.3.0) and TreeView (Ver.1.1.1) analyses. Genes with significant similarities to the transcripts in nr database based on BLASTX searches were selected for GO analysis, performed by MAS 3.0 software which was based on DAVID database (CapitalBio, Beijing, China) [16]. Annotation results were obtained by inputting the list of gene symbol as identifier [18]. The Pathway analysis was done using the MAS 3.0 software which was based on the Kyoto Encyclopedia of Genes and Genomes database (KEGG) (CapitalBio, Beijing, China). All microarray results from this study were deposited in NCBI'S Gene Expression Omnibus (GEO) database, accession numbers are: Platform, GPL 3533, Samples, GSM 747145, GSM 747146, GSM 747147, GSM, 747148, GSM 747149, GSM 747150 with the series accession number GSE 30172.

QPCR analysis

Total RNA were extracted from the PAMs of each group with Trizol (Invitrogen) and 5 μg of total RNA were used for first strand cDNA synthesis by using Superscript II cDNA amplification System (Invitrogen) following manufacturer's instructions. Real-time PCR was performed using LightCycler 480 (Roche Applied Science) and Quantitect SYBR Green PCR kit (Roche) following the companies' instructions. Briefly, PCR assay was performed under the following conditions: 95°C for 15 sec, 55°C for 15 sec and 72°C for 15 sec. Real-time PCR primers for each gene were indicated in Table 3. All the primers were originally designed using Primer 3 software (Rozen & Skaletsky, 2000) or according to the published papers. Results were calculated by minus delta delta threshold cycle (-ddCt) method. Briefly, the threshold cycle Ct1 of each sample reaction were deducted with the threshold cycle Ct2 of GAPDH reaction for normalization, then deducted from the threshold cycle Ct3 of calibration control (40 Cycles in this experiment); thus, the final result was represented by the formula: Ct3-(Ct1-Ct2).

Expression of S100A4, S100A6 in PK-15 cells stimulated with LPS and Poly (I:C)

PK-15 cells have been shown especially useful for the study of infectious disease processes in swine [86, 87]. In this study, 12 groups (with three repeats in each group, ~1 × 105 cells/samples) of PK-15 cells were grown in culture medium (DMEM) supplemented with 10% heat-inactivated fetal bovine serum at 37°C with 5% CO2. Adherent PK-15 cells were obtained by washing off nonadherent cells with warm culture medium and PBS twice, respectively. Adherent cells were further cultured in DMEM (control samples) or treated with 1 μg/mL LPS (Sigma-Aldrich, E.coli 0127:B8) or 10 μg/mL Poly (I:C) (Sigma-Aldrich) respectively (stimulated samples) for 0 h, 2 h, 6 h, 12 h, 24 h and 48 h. Cells were harvested and total RNA were extracted as described above.

DNA preparation from bacterial isolates and clinical samples

Bacterial cultures were harvested from trypticase soy agar (TSA) using an inoculation loop and were placed into a 1.5 mL tube to which was added with 500 μL of phosphate buffered saline (PBS). One milliliter of the fluid and 0.5 g of the tissue samples were respectively placed in sterile tubes containing 5 mL of trypticase soy broth (TSB), 5 μL nicotinamide adenine dinucleotide (NAD) and 500 μL sterilized fetal bovine serum and then incubated for 8 h at 37°C with agitation. Five hundred microliters of the suspension was removed to a new 1.5 mL tube. Tubes containing bacteria, tissue and fluid suspensions were centrifuged at 13,400 g for 5 min. After centrifugation, the supernatant was discarded and the remaining pellet was suspended in 200 μL PBS, boiled for 10 min. After boiling, tubes were centrifuged at 13,400 g for 5 min. Fifty microliters of supernatant from each sample containing extracted DNA were mixed with 50 μL of Tris-EDTA buffer and stored at 4°C. This final solution was used as DNA template in nested PCR and LAMP reaction. The primers for nested PCR and LAMP were listed in Additional file 9. The procedure of bacterial isolation, nested PCR and LAMP were carried out according to description of Wang et al [6].

Detection of s100a4 and s100a6 expression in different tissues

Three pigs in H. parasuis infection group and control group were selected for the analysis of s100a4 and s100a6 expression in different tissues. Total RNA from 5 porcine organs (inguinal lymph node, heart, spleen, lung, brain) was isolated with RNAprep pure Tissue Kit (TianGen Biotech (Beijing) Co., Ltd). Total RNA was then quantified by NanoDrop 1000 Spectrophotometer (Thermo Fisher Scientific Inc., USA). The quality of the RNA was checked by formaldehyde denaturing gel electrophoresis in 1.2% agarose gels, which showed dispersed bands (28S and 18S) without any obvious smearing patterns that would indicate degradation. Two microgram of total RNA was used for reverse transcription polymerase chain reaction, using the TransSript First Strand cDNA Synthesis SuperMix according to the manufacturer's instructions (TianGen Biotech (Beijing) Co., Ltd). The qPCR assays were performed and analyzed as described above, with primers listed in Table 3.

Data for STRING and IPA analysis

Differentially expressed (DE) genes were analyzed using STRING, a database of known and predicted protein interaction for DE gene encoded proteins. The results were obtained by inputting the list of gene symbol as identifier (organism = sus scrofa, combined score = 0.15). Ppp1r13l gene was selected for network exploration using Ingenuity Pathway Analysis (Ingenuity® Systems, The data set containing gene identifier and corresponding expression value was uploaded into in the application. The identifier was mapped to its corresponding object in Ingenuity's Knowledge Base; Network Eligible molecules were then overlaid onto a global molecular network so that network of Network Eligible Molecules could be algorithmically generated based on their connectivity.


H. parasuis :

Haemophilus parasuis


differentially expressed


fold change


Gene Ontology


quantitative real-time PCR


porcine alveolar macrophages



Poly (I:C):

polyinosinic acid-polycytidylic acid


pig kidney-15


loop-mediated isothermal amplification. TSA: trypticase soy agar


trypticase soy broth


nicotinamide adenine dinucleotide. MAS: Molecular Annotation System. SAM: Significance Analysis of Microarrays. FDR: false discovery rate. IPA: Ingenuity Pathway Analysis. DAVID: Database for Annotation, Visualization and Integrated Discovery. ORF: open reading frame.


  1. 1.

    Lunney JK: Advances in swine biomedical model genomics. Int J Biol Sci. 2007, 3: 179-184.

  2. 2.

    Yang W, Pin C, Haibing G, Hui L, Qigai H: Loop mediated isothermal amplification targeting the apxIVA gene for detection of Actinobacillus pleuropneumoniae. FEMS microbiol. lett. 2009, 300: 83-89. 10.1111/j.1574-6968.2009.01779.x.

  3. 3.

    Oliveira S, Pijoan C: Haemophilus parasuis: new trends on diagnosis, epidemiology and control. Vet Microbiol. 2004, 99: 1-12. 10.1016/j.vetmic.2003.12.001.

  4. 4.

    Angen O, Oliver S, Ahrens P, Svensmark B, Leser TD: Development of an improved species specific PCR test for detection of Haemophilus parasuis. Vet Microbiol. 2007, 119: 266-276. 10.1016/j.vetmic.2006.10.008.

  5. 5.

    Amano H, Shibata M, Kajio N, Morozumi T: Pathologic observations of pigs intranasally inoculated with serovar 1, 4 and 5 of Haemophilus parasuis using immunoperoxidase method. J Vet Med Sci. 1994, 56: 639-644. 10.1292/jvms.56.639.

  6. 6.

    Yang W, Ying F, Yingyu L, Pin C, Wentao L, Shuqing L, Haoyong Z, Qigai H: Development and evaluation of loop-mediated isothermal amplification for rapid detection of Haemophilus parasuis. FEMS Microbiol. Lett. 2010, 313: 54-60. 10.1111/j.1574-6968.2010.02126.x.

  7. 7.

    Aderem A, Ulevitch RJ: Toll-like receptors in the induction of the innate immune response. Nature. 2000, 406: 782-787. 10.1038/35021228.

  8. 8.

    McGuire K, Glass EJ: The expanding role of microarrays in the investigation of macrophage responses to pathogens. Vet Immunol Immunop. 2005, 105: 259-275. 10.1016/j.vetimm.2005.02.001.

  9. 9.

    Andersson H, Hartmanova B, Ryden P, Noppa L, Naslund L, Sjostedt A: A microarray analysis of the murine macrophage response to infection with Francisella tularensis LVS. J Med Microbiol. 2006, 55: 1023-1033. 10.1099/jmm.0.46553-0.

  10. 10.

    Nau GJ, Richomond JFL, Schlesinger A, Jennings EG, Lander ES, Young RA: Human macrophage activation programs induced by bacterial pathogens. P Natl Acad Sci USA. 2002, 99: 1503-1508. 10.1073/pnas.022649799.

  11. 11.

    Eskra L, Mathison A, Splitter G: Microarray analysis of mRNA levels from RAW2647 macrophages infected with Brucella abortus. Infect Immun. 2003, 71: 1125-1133. 10.1128/IAI.71.3.1125-1133.2003.

  12. 12.

    Detweiler CS, Cunanan DB, Falkow S: Host microarray analysis reveals a role for the Salmonella response regulator phoP in human macrophage cell death. P Natl Acad Sci USA. 2001, 98: 5850-5855. 10.1073/pnas.091110098.

  13. 13.

    Olvera A, Ballester M, Norfrarias M, Sibila M, Aragon V: Differences in phagocytosis susceptibility in Haemophilus parasuis strains. Vet Res. 2009, 40: 24-10.1051/vetres/2009007.

  14. 14.

    Auger E, Deslandes V, Ramjeet M, Contreras I, Nash JHE, Harel J, Gottschalk M, Olivier M, Jacques M: Host pathogen interactions of Actinobacillus pleuropneumoniae with porcine lung and tracheal epithelial cells. Infect Immun. 2009, 77: 1426-1441. 10.1128/IAI.00297-08.

  15. 15.

    Ichikawa JK, Bangera ANMG, Geiss GK, Wout AB, Bumgarner RE, Lory S: Interaction of Pseudomonas aeruginosa with epithelial cells: identification of differentially regulated genes by expression microarray analysis of human cDNAs. P Natl Acad Sci USA. 2000, 97: 9659-9664.

  16. 16.

    Ran L, Zhang A, Chen B, Teng L, Wang Y, Chen H, Jin M: Response of swine spleen to Streptococcus suis infection revealed by transcription analysis. BMC genomics. 2010, 11: 556-10.1186/1471-2164-11-556.

  17. 17.

    Tusher VG, Tibshirani R, Chu G: Significance analysis of microarrays applied to the ionizing radiation response. P Natl Acad Sci USA. 2001, 98: 5116-5121. 10.1073/pnas.091062498.

  18. 18.

    Li Y, Zhou H, Wen Z, Wu S, Huang C, Jia G, Chen H, Jin Me: Transcription analysis on response of swine lung to H1N1 swine influenza virus. BMC genomics. 2011, 12: 398-10.1186/1471-2164-12-398.

  19. 19.

    Wang D, Pan Y, Zhao X, Zhu L, Fu B, Li Z: Genome-wide temporal-spatial gene expression profiling of drought responsiveness in rice. BMC genomics. 2011, 12: 149-10.1186/1471-2164-12-149.

  20. 20.

    Zou X, Zou L, Foldager C, Bendtsen M, Feng W, Bunger C: Different mechanisms of spinal fusion using equine bone protein extract, rhBMP-2 and autograft during the process of anterior lumbar interbody fusion. Biomaterials. 2009, 30: 991-1004. 10.1016/j.biomaterials.2008.10.061.

  21. 21.

    Bonetta L: Protein-protein interactions: Interactome under construction. Nature. 2010, 468: 851-854. 10.1038/468851a.

  22. 22.

    Von Mering C, Lars J, Jensen , Berend , Snel , Sean D, Hooper , Markus , Krupp , Mathilide , Foglierini , Nelly , Jouffre , Martijin A, Huynen , Peer , Bork : STRING: known and predicted protein Cprotein associations, integrated and transferred across organisms. Nucleic Acids Res. 2005, 33: 433-437.

  23. 23.

    Schneider M, Kostin S, Strom C, Aplin M, Lyngbak S, Theilade J, Grigorian M, Andersen CB, Lukanidin E, Hansen JL, Sheikh SP: S100A4 is upregulated in injured myocardium and promotes growth and survival of cardiac myocytes. Cardiovasc Res. 2007, 75: 40-50. 10.1016/j.cardiores.2007.03.027.

  24. 24.

    Iwano M, Plieth D, Danoff TM, Xue C, Okada H, Neilson EG: Evidence that fibroblasts derive from epithelium during tissue fibrosis. J Clin Invest. 2002, 110: 341-350.

  25. 25.

    Nishitani Y, Iwano M, Yanaguchi Y, Nakatani K, Akai Y, Nishino T, Shiiki H, Kanauchi M, Saito Y, Neilson EG: Fibroblast-specific protein 1 is a specific prognostic marker for renal survival in patients with IgAN1. Kidney Int. 2005, 68: 1078-1085. 10.1111/j.1523-1755.2005.00500.x.

  26. 26.

    Rivard CJ, Brown LM, Almeida NE, Maunsbach AB, Pihakashi-Maunsach K, Andres-Hernando , Capasso JM, Berl T: Expression of the calcium-binding protein S100A4 is markedly up-regulated by osmotic stress and is involved in the renal osmoadaptive response. J Biol Chem. 2007, 282: 6644-6652.

  27. 27.

    Breen EC, Tang K: Calcyclin (S100A6) regulates pulmonary fibroblast proliferation, morphology, and cytoskeletal organization in vitro. J Cell Biochem. 2003, 88: 848-854. 10.1002/jcb.10398.

  28. 28.

    Slomnicki LP, Nawrot B, Lesniak W: S100A6 binds p53 and affects its activity. Int J Biochem Cell B. 2009, 41: 784-790. 10.1016/j.biocel.2008.08.007.

  29. 29.

    Suzuki K, Nishihata J, Arai Y, Honma N, Yamamoto K, Irimura T, Toyoshima S: Molecular cloning of a novel actin-binding protein, p57, with a WD repeat and a leucine zipper motif. FEBS lett. 1995, 364: 283-288. 10.1016/0014-5793(95)00393-N.

  30. 30.

    Ferrari G, Langen H, Naito M, Pieters J: A coat protein on phagosomes involved in the intracellular survival of mycobacteria. Cell. 1999, 97: 435-447. 10.1016/S0092-8674(00)80754-0.

  31. 31.

    De Hostos EL: The coronin family of actin-associated proteins. Trends Cell Biol. 1999, 9: 345-350. 10.1016/S0962-8924(99)01620-7.

  32. 32.

    Gatfield J, Albrecht I, Zanolari B, Steinmetz MO, Pieters J: Association of the leukocyte plasma membrane with the actin cytoskeleton through coiled coil-mediated trimeric coronin 1 molecules. Mol Biol Cell. 2005, 16: 2786-2798. 10.1091/mbc.E05-01-0042.

  33. 33.

    Gatfield J, Pieters J: Essential role for cholesterol in entry of mycobacteria into macrophages. Science. 2000, 288: 1647-1651. 10.1126/science.288.5471.1647.

  34. 34.

    Tanigawa K, Suzuki K, Kimura H, Takeshita F, Wu H, Akama T, Kawashima A, Ishii N: Tryptophan aspartate-containing coat protein (CORO1A) suppresses Toll-like receptor signalling in Mycobacterium leprae infection. Clin Exp Immunol. 2009, 156: 495-501. 10.1111/j.1365-2249.2009.03930.x.

  35. 35.

    Liu XD, Chen HB, Tong Q, Li XY, Zhu MJ, Wu ZF, Zhou R, Zhao SH: Molecular Characterization of Caveolin-1 in Pigs Infected with Haemophilus parasuis. J Immunol. 2011, 186: 3031-3046. 10.4049/jimmunol.0902687.

  36. 36.

    Rosenberger CM, Finlay BB: Phagocyte sabotage: disruption of macrophage signalling by bacterial pathogens. Nat Rev Mol Cell Bio. 2003, 4: 385-396. 10.1038/nrm1104.

  37. 37.

    Harper M, Boyce JD, Adler B: Pasteurella multocida pathogenesis: 125 years after Pasteur. FEMS microbiol. lett. 2006, 265: 1-10. 10.1111/j.1574-6968.2006.00442.x.

  38. 38.

    Inzana TJ: Capsules and virulence in the HAP group of bacteria. Can J Vet Res. 1990, 54 (Suppl): S22-27.

  39. 39.

    Inzana T, Ma J, Workman T, Gogolewski RP, Anderson P: Virulence properties and protective efficacy of the capsular polymer of Haemophilus (Actinobacillus) pleuropneumoniae serotype 5. Infect Immun. 1988, 56: 1880-1889.

  40. 40.

    Noel GJ, Hoiseth SK, Edelson PJ: Type b capsule inhibits ingestion of Haemophilus influenzae by murine macrophages: studies with isogenic encapsulated and unencapsulated strains. J Infect Dis. 1992, 166: 178-182. 10.1093/infdis/166.1.178.

  41. 41.

    Jiang Z, Georgel P, Du X, Shamel L, Sovath S, Mudd S, Huber M, Kalis C, Keck S, Galanos C: CD14 is required for MyD88-independent LPS signaling. Nature immunology. 2005, 6: 565-570. 10.1038/ni1207.

  42. 42.

    Schiff DE, Kline L, Soldau K, Lee JD, Pugin J, Tobias PS, Ulevitch RJ: Phagocytosis of gram-negative bacteria by a unique CD14-dependent mechanism. J Leukoc Biol. 1997, 62: 786-794.

  43. 43.

    Schlegel RA, Krahling S, Callahan MK, Williamson P: CD14 is a component of multiple recognition systems used by macrophages to phagocytose apoptotic lymphocytes. Cell death and differ. 1999, 6: 583-592. 10.1038/sj.cdd.4400529.

  44. 44.

    Wright SD, Ramos RA, Tobias PS, Ulevitch RJ, Mathison JC: CD14, a receptor for complexes of lipopolysaccharide (LPS) and LPS binding protein. Science. 1990, 249: 1431-1433. 10.1126/science.1698311.

  45. 45.

    Zhao C, Zhang H, Wong WC, Sem X, Han H, Ong SM, Tan YC, Yeap WH, Gan CS, Ng KQ: Identification of novel functional differences in monocyte subsets using proteomic and transcriptomic methods. J Proteome Res. 2009, 8: 4028-4038. 10.1021/pr900364p.

  46. 46.

    Bliska JB, Black DS: Inhibition of the Fc receptor-mediated oxidative burst in macrophages by the Yersinia pseudotuberculosis tyrosine phosphatase. Infect Immun. 1995, 63: 681-685.

  47. 47.

    Hamon Y, Chambenoit O, Luciani M, Toti F, Chaslin S, Freyssinet J, Devaux PF, McNeish J, Marguet D: ABC1 promotes engulfment of apoptotic cells and transbilayer redistribution of phosphatidylserine. Nat Cell Biol. 2000, 2: 399-406. 10.1038/35017029.

  48. 48.

    Uterm hlen O, Herz J, Schramm M, Kr nke M, et al: Fusogenicity of membranes: the impact of acid sphingomyelinase on innate immune responses. Immunobiology. 2008, 213: 307-314. 10.1016/j.imbio.2007.10.016.

  49. 49.

    Gladue DP, Zhu J, Holinka LG, Sainz IF, Carrillo C, Prarat MV, O'Donnell V, Borca MV: Patterns of gene expression in swine macrophages infected with classical swine fever virus detected by microarray. Virus Res. 2010, 151: 10-18. 10.1016/j.virusres.2010.03.007.

  50. 50.

    Wilkinson JM, Sargent CA, Pantoja LG, Tucher AW: Gene expression profiling in the lungs of pigs with different susceptibilities to Glässer's disease. BMC genomics. 2010, 11: 455-10.1186/1471-2164-11-455.

  51. 51.

    Wang YW, Luo R, Fang L, Wang D, Bi J, Chen H, Xiao S: Porcine reproductive and respiratory syndrome virus (PRRSV) infection activates chemokine RANTES in MARC-145 cells. Mol Immunol. 2010, 48: 586-591.

  52. 52.

    Zlotnik A, Yoshie O: Chemokines: A New Classification Review System and Their Role in Immunity. Immunity. 2000, 12: 121-127. 10.1016/S1074-7613(00)80165-X.

  53. 53.

    Chen CJ, Chen JH, Chen SY, Liao SL, Raung SL: Upregulation of RANTES gene expression in neuroglia by Japanese encephalitis virus infection. J Virol. 2004, 78: 12107-12119. 10.1128/JVI.78.22.12107-12119.2004.

  54. 54.

    Casola A, Garofalo RP, Haeberle H, Elliott TF, Lin R, Jamaluddin M, Brasier AR: Multiple cis regulatory elements control RANTES promoter activity in alveolar epithelial cells infected with respiratory syncytial virus. J Virol. 2001, 75: 6428-6439. 10.1128/JVI.75.14.6428-6439.2001.

  55. 55.

    Kujime K, Hashimoto S, Gon Y, Shimizu K, Horie T: p38 mitogen-activated protein kinase and c-jun-NH2-terminal kinase regulate RANTES production by influenza virus-infected human bronchial epithelial cells. J Immunol. 2000, 164: 3222-3228.

  56. 56.

    Veckman V, Miettinen M, Matikainen S, Lande R, Giacomini E, Coccia M, Julkunen I: Lactobacilli and streptococci induce inflammatory chemokine production in human macrophages that stimulates Th1 cell chemotaxis. J Leukocyte Biol. 2003, 74: 395-402. 10.1189/jlb.0402212.

  57. 57.

    Pietilä TE, Veckman V, Kyllonen P, Lahteenmaki K, Korhonen TK, Julknen I: Activation, cytokine production, and intracellular survival of bacteria in Salmonella-infected human monocyte-derived macrophages and dendritic cells. J Leukocyte Biol. 2005, 78: 909-920. 10.1189/jlb.1204721.

  58. 58.

    Hamilton TA, Ohmori Y, Tebo JM, Kishore R: Regulation of macrophage gene expression by pro-and anti-inflammatory cytokines. Pathobiology. 2000, 67: 241-244.

  59. 59.

    Opal SM, Depalo VA: Anti-Inflammatory Cytokines. Chest. 2000, 117: 1162-1172. 10.1378/chest.117.4.1162.

  60. 60.

    Franscini N, Bachli EB, Blau N, Fischler M, Walter R, Schlaffner A, Schoedon G: Functional tetrahydrobiopterin synthesis in human platelets. Circulation. 2004, 110: 186-192. 10.1161/01.CIR.0000134281.82972.57.

  61. 61.

    Ikeda M, Matsui K, Ishihara Y, Morita I, Murota S, Yuasa T, Miyatake T: Cerebellar nitric oxide synthase, cGMP and motor function in two lines of cerebellar mutant mice, Staggerer and Wriggle Mouse Sagami. Neurosci lett. 1994, 168: 65-68. 10.1016/0304-3940(94)90417-0.

  62. 62.

    Diefenbach A, Hsia JK, Hsiung MYB, Raulet DH: A novel ligand for the NKG2D receptor activates NK cells and macrophages and induces tumor immunity. Eur J Immunol. 2003, 33: 381-391. 10.1002/immu.200310012.

  63. 63.

    Wu Z, Qiu L: Effect of nitric oxide on iron-mediated cytotoxicity in primary cultured renal proximal tubules. Cell Biochem Funct. 2001, 19: 237-247. 10.1002/cbf.921.

  64. 64.

    Alderton WK, Cooper CE, Knowles RG: Nitric oxide synthases: structure, function and inhibition. Biochem J. 2001, 357: 593-615. 10.1042/0264-6021:3570593.

  65. 65.

    Bakalian A, Kopmels B, Messer A, Fradelizi D, Delhaye-Bouchaud N, Wollman E, Mariani J: Peripheral macrophage abnormalities in mutant mice with spinocerebellar degeneration. Res Immunol. 1992, 143: 129-139. 10.1016/0923-2494(92)80090-8.

  66. 66.

    Jamieson AM, Diefenbach A, McMahon CW, Xiong N, Carlyle JR, Raulet DH: The role of the NKG2D immunoreceptor in immune cell activation and natural killing. Immunity. 2002, 17: 19-29. 10.1016/S1074-7613(02)00333-3.

  67. 67.

    Chen H, Lunney JK, Cheng L, Li X, Cao J, Zhu M, Zhao : S. Porcine S100A8 and S100A9: molecular characterizations and crucial functions in response to Haemophilus parasuis infection. Dev Comp Immunol. 35: 490-500.

  68. 68.

    Foell D, Frosch M, Sorg C, Roth J: Phagocyte-specific calcium-binding S100 proteins as clinical laboratory markers of inflammation. Clinica Chim Acta. 2004, 344: 37-51. 10.1016/j.cccn.2004.02.023.

  69. 69.

    Xu Y, Wan J, Jiang D, Wu X: BMP-7 counteracts TGF-beta1-induced epithelial-to-mesenchymal transition in human renal proximal tubular epithelial cells. J Nephrol. 2009, 22: 403-410.

  70. 70.

    Zeisberg M, Yang C, Martino M, Duncan M, Rieder F, Tanjore H, Kalluri R: Fibroblasts derive from hepatocytes in liver fibrosis via epithelial to mesenchymal transition. J Biol Chem. 2007, 282: 23337-23347. 10.1074/jbc.M700194200.

  71. 71.

    Robertson H, Kirby JA, Yip WW, Jones DE, Burt AD: Biliary epithelial-mesenchymal transition in posttransplantation recurrence of primary biliary cirrhosis. Hepatology. 2007, 45: 977-981. 10.1002/hep.21624.

  72. 72.

    Greenway S, van Suylen RJ, Du Marchie Sarvaas G, Kwan E: S100A4/Mts1 produces murine pulmonary artery changes resembling plexogenic arteriopathy and is increased in human plexogenic arteriopathy. Am J Pathol. 2004, 164: 253-262. 10.1016/S0002-9440(10)63115-X.

  73. 73.

    Kwapiszewska G, Wilhelm J, Wolff S, Laumanns I, Koenig IR, Ziegler A, Seeger W, Bohle RM, Weissmann N, Fink L: Expression profiling of laser-microdissected intrapulmonary arteries in hypoxia-induced pulmonary hypertension. Respir Res. 2005, 6: 109-10.1186/1465-9921-6-109.

  74. 74.

    Lawson WE, Polosukhin VV, Zoia O, Stathopoulos GT, Han W, Plieth D, Loyd JE, Neilson EG, Blackwell TS: Characterization of fibroblast-specific protein 1 in pulmonary fibrosis. Am J Respir Crit Care Med. 2005, 171: 899-907. 10.1164/rccm.200311-1535OC.

  75. 75.

    Inamoto S, Murao S, Yokoyama M, Kitazawa S, Maeda S: Isoproterenol-induced myocardial injury resulting in altered S100A4 and S100A11 protein expression in the rat. Pathol Int. 2000, 50: 480-485. 10.1046/j.1440-1827.2000.01069.x.

  76. 76.

    Zeisberg EM, Tarnavski O, Zeisberg M, Dorfman AL, McMullen JR, Gustafsson E, Chandraker A, Yuan X, Pu WT, Roberts AB, Neilson EG, Sayegh MH, Izumo S, Kalluri R: Endothelial-to-mesenchymal transition contributes to cardiac fibrosis. Nat Med. 2007, 13: 952-961. 10.1038/nm1613.

  77. 77.

    Okada H, Inoue T, Kanno Y, Kobayashi T, Watanabe Y, Ban S, Neilson EG, Suzuki H: Selective depletion of fibroblasts preserves morphology and the functional integrity of peritoneum in transgenic mice with peritoneal fibrosing syndrome1. Kidney international. 2003, 64: 1722-1732. 10.1046/j.1523-1755.2003.00290.x.

  78. 78.

    Gargalovic P, Dory L: Caveolin-1 and Caveolin-2 expression in mouse macrophages. J Biol Chem. 2001, 276: 26164-26170. 10.1074/jbc.M011291200.

  79. 79.

    Zaas DW, Duncan MJ, Li G, Wright JR, Abraham SN: Pseudomonas invasion of type I pneumocytes is dependent on the expression and phosphorylation of caveolin-2. J Biol Chem. 2005, 280: 4864-4872.

  80. 80.

    Mora R, Bonilha VL, Marmorstein A, Scherer PE, Brown D, Lisanti MP, Rodriguez-Boulan E: Caveolin-2 localizes to the golgi complex but redistributes to plasma membrane, caveolae, and rafts when co-expressed with caveolin-1. J Biol Chem. 1999, 274: 25708-25717. 10.1074/jbc.274.36.25708.

  81. 81.

    Yang JP, Hori M, Sanda T, Okamoto T: Identification of a novel inhibitor of nuclear factor-κB, RelA-associated inhibitor. J Biol Chem. 1999, 274: 15662-15670. 10.1074/jbc.274.22.15662.

  82. 82.

    Takada N, Sanda T, Okamoto H, Yang JP, Asamitsu K, Sarol L, Kimura G, Uranishi H, Tetsuka T, Okamoto T: RelA-associated inhibitor blocks transcription of human immunodeficiency virus type 1 by inhibiting NF-kappa B and Sp1 actions. J Virol. 2002, 76: 8019-8030. 10.1128/JVI.76.16.8019-8030.2002.

  83. 83.

    Irizarry RA, Hobbs B, Collin F, Beazer-Barclay YD, Antonellis KJ, Schert U, Speed TP: Exploration, normalization, and summaries probe level data. Biostatistics. 2003, 4: 249-264. 10.1093/biostatistics/4.2.249.

  84. 84.

    Kasler HG, Young BD, Mottet D, Lim HW, Collins AM, Olson EN, Verdin E: Histone deacetylase 7 regulates cell survival and TCR signaling in CD4/CD8 double-positive thymocytes. J Immunol. 2011, 186: 4782-4793. 10.4049/jimmunol.1001179.

  85. 85.

    Wu Z, Zhang W, Lu Y, Lu C: Transcriptome profiling of zebrafish infected with Streptococcus suis. Microbial pathogenesis. 2010, 48: 178-187. 10.1016/j.micpath.2010.02.007.

  86. 86.

    Boury BN, Stabel TJ, Kehrli JM, Taylor M: Comparison of the PK (15)-and WEHI 164 (clone 13)-based bioassays for detection of porcine tumor necrosis factor. Am J Vet Res. 1997, 58: 1115-1119.

  87. 87.

    Chen H, Li C, Fang M, Zhu M, Li X, Zhou R, Li K, Zhao S: Understanding Haemophilus parasuis infection in porcine spleen through a transcriptomics approach. BMC genomics. 2009, 10: 64-10.1186/1471-2164-10-64.

Download references


We thank Dr. Liangjun Zhao (Wayne state university, USA), Xiangdong Liu, Hongbo Chen (Huazhong agricultural university, Wuhan, China) for the technical assistance.

Author information

Correspondence to Qigai He.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

YW and CL carried out all works in the lab and drafted the manuscript. YF, XL, WL, SL, YL, YL, and CC participated in the animal challenge experiment and immunoassays. CC and JCA participated in the experiment design and coordination. QH conceived the study, and participated in its coordination and helped to draft the manuscript. All authors read and approved the final manuscript.

Electronic supplementary material

Additional file 1: PAMs transcriptome analysis following H. parasuis infection using the Affymetrix Porcine Genechip. (XLS 3 MB)

Additional file 2: 575 transcripts that are differentially expressed in PAM following H. parasuis infection. (XLS 124 KB)

Additional file 3: Categories of annotated DE genes based on biological process GO term. Many categories shared the same transcripts. (XLS 38 KB)

Additional file 4: KEGG Pathway analysis of annotated DE genes. (XLS 26 KB)

Additional file 5: STRING analysis of all annotated DE genes. (DOC 348 KB)

Additional file 6: Phylogenetic tree of coronin 1a. (DOC 120 KB)

Additional file 7: Multiple alignment of the porcine coronin 1a protein with other 8 known coronin 1a proteins. (DOC 742 KB)

Additional file 8: Ppp1r13l related genes by IPA analysis. (DOC 156 KB)

Additional file 9: Primers for nested PCR and LAMP. (DOC 33 KB)

Authors’ original submitted files for images

Rights and permissions

Reprints and Permissions

About this article

Cite this article

Wang, Y., Liu, C., Fang, Y. et al. Transcription analysis on response of porcine alveolar macrophages to Haemophilus parasuis. BMC Genomics 13, 68 (2012).

Download citation


  • Differentially Express Gene
  • Japanese Encephalitis Virus
  • S100a6 Expression
  • Porcine Alveolar Macrophage
  • S100 Family Gene