Antigen stimulation of peripheral blood mononuclear cells from Mycobacterium bovis infected cattle yields evidence for a novel gene expression program
© Meade et al; licensee BioMed Central Ltd. 2008
Received: 18 March 2008
Accepted: 29 September 2008
Published: 29 September 2008
Bovine tuberculosis (BTB) caused by Mycobacterium bovis continues to cause substantial losses to global agriculture and has significant repercussions for human health. The advent of high throughput genomics has facilitated large scale gene expression analyses that present a novel opportunity for revealing the molecular mechanisms underlying mycobacterial infection. Using this approach, we have previously shown that innate immune genes in peripheral blood mononuclear cells (PBMC) from BTB-infected animals are repressed in vivo in the absence of exogenous antigen stimulation. In the present study, we hypothesized that the PBMC from BTB-infected cattle would display a distinct gene expression program resulting from exposure to M. bovis. A functional genomics approach was used to examine the immune response of BTB-infected (n = 6) and healthy control (n = 6) cattle to stimulation with bovine tuberculin (purified protein derivative – PPD-b) in vitro. PBMC were harvested before, and at 3 h and 12 h post in vitro stimulation with bovine tuberculin. Gene expression changes were catalogued within each group using a reference hybridization design and a targeted immunospecific cDNA microarray platform (BOTL-5) with 4,800 spot features representing 1,391 genes.
250 gene spot features were significantly differentially expressed in BTB-infected animals at 3 h post-stimulation contrasting with only 88 gene spot features in the non-infected control animals (P ≤ 0.05). At 12 h post-stimulation, 56 and 80 gene spot features were differentially expressed in both groups respectively. The results provided evidence of a proinflammatory gene expression profile in PBMC from BTB-infected animals in response to antigen stimulation. Furthermore, a common panel of eighteen genes, including transcription factors were significantly expressed in opposite directions in both groups. Real-time quantitative reverse transcription PCR (qRT-PCR) demonstrated that many innate immune genes, including components of the TLR pathway and cytokines were differentially expressed in BTB-infected (n = 8) versus control animals (n = 8) after stimulation with bovine tuberculin.
The PBMC from BTB-infected animals exhibit different transcriptional profiles compared with PBMC from healthy control animals in response to M. bovis antigen stimulation, providing evidence of a novel gene expression program due to M. bovis exposure.
Mycobacterium bovis infection is the cause of bovine tuberculosis (BTB), an important health problem in cattle that also has zoonotic potential for transmission to humans. The eradication of M. bovis infection in cattle is proving difficult in some developed countries, including the UK and Ireland  due to limitations in the sensitivity of current diagnostics, leading to a failure to detect all infected animals [2, 3]. It is also unclear what role exposure to environmental mycobacterial antigens play in the generation of non-specific immune responses, giving rise to difficulties with test interpretation and reliability. Furthermore, some cattle with advanced disease become anergic, with suppressed cellular immune responses in both the peripheral blood and at the site of the infection , and remain undetected as reservoirs of disease.
The immune response to tuberculosis is a complex process and studies in the bovine model have primarily focused on the adaptive response. The outcome of tuberculosis infection depends on T cell interaction with macrophages , and progression of infection with M. bovis is thought to develop after a shift in the immune system from a protective proinflammatory, cytotoxic T cell response towards a non-protective antibody-mediated response [6, 7]. The timing and potency of the cellular and immunological events that occur immediately post-infection are crucial determinants governing infection  and innate immune responses are considered important for the generation of early, and appropriate adaptive responses to resolve infection [8–11]. Therefore, central to the development of improved or novel diagnostics is increased understanding of the early immune response to tuberculosis in cattle.
At the most basic level, the interplay between the host and pathogen involves changes in gene expression . High-throughput genomics technologies, which offer the ability to survey changes in expression for a large number of genes simultaneously, have been widely used to discern patterns of host gene regulation during infection. Microarray technology has emerged as the method of choice for large-scale gene expression studies that have increased our understanding of host-pathogen interactions [13–20, 12]. These functional genomics studies have also identified new avenues of research for potential control strategies against pathogens . Using this approach we have aimed to gain a better understanding of the molecular regulation of the immune response following M. bovis exposure and infection in cattle, with the expectation of significant benefits in development of new practical tools applicable to disease control .
We have previously used a bovine targeted immunospecific cDNA microarray to study gene expression changes in bovine peripheral blood mononuclear cells (PBMC) from BTB-infected cattle cultured in vitro in the presence of bovine and avian tuberculins . PBMC are an accessible tissue for the development of robust novel diagnostics and previous studies have shown that for bovine tuberculosis, immune responses occurring in the peripheral blood reflect those at the site of disease . In our previous study, antigen stimulation induced significant and specific expression changes in many immune genes that revealed different gene expression patterns in stimulated and non-stimulated PBMC. Although IFNG gene expression was increased in response to antigen stimulation, several other genes were highly expressed earlier in the time course , suggesting that these genes may represent valuable targets for the development of novel diagnostics of M. bovis in cattle early post-infection.
A subsequent comparative study between BTB-infected cattle and healthy control animal PBMC showed that the expression of innate immune genes is repressed in heavily infected cattle in vivo; it also demonstrated that the expression changes for many of these genes represent a BTB-signature of infection . Additional expression profiling using real time qRT-PCR verified that a number of innate immune genes including TLR2 and TLR4 had reduced expression in BTB-infected animals . Failure of M. bovis to activate these immune genes may have resulted from an antigen-induced suppression, or alternatively from a general failure of the immune system in the advanced diseased cattle.
For the present study, we compared PBMC gene expression profiles from natural M. bovis-infected animals with those from control non-infected animals after stimulation with bovine tuberculin (purified protein derivative – PPD-b) antigens. We hypothesized that the PBMC from BTB-infected cattle would display a distinct gene expression program resulting, at least in part, from previous exposure to M. bovis. These data provide evidence for a novel gene expression program in PBMC from infected animals and highlight the value of large-scale genomics approaches to understand immune regulatory mechanisms that in future may form the basis for novel diagnostics and therapeutics.
Distinct microarray gene expression profile in PBMC from BTB-infected cattle at 3 h and 12 h post-stimulation with bovine tuberculin
Inspection of these data revealed that the 250 early differentially expressed spot features included 83 unique genes where two replicate gene spot features were found to be significantly differentially expressed (6% of the 1,391 duplicated genes represented on the array). Fifty-three of these genes were identified as BOTL clones, ESTs derived from genes whose function(s) in cattle are unknown, but were inferred from evolutionary sequence homology with gene orthologs in rodents and humans (see Additional file 1). Similarly, between the later two time points, 18 genes were represented by significantly differentially expressed duplicate significant spot features (P ≤ 0.05).
In contrast to the differential gene expression profile observed in BTB-infected animals in response to bovine tuberculin stimulation, initial gene expression changes in the control PBMC samples were moderate with 88 spot features significantly differentially expressed between 0 h and 3 h (P ≤ 0.05) [see Additional file 2]. These were almost evenly divided between increased (41) and decreased (47) expression at the later time point (Fig. 2b). Between 3 h and 12 h of stimulation, 56 spot features were differentially expressed (see Additional file 2), again almost equally balanced between repression and activation (26 spot features increased in expression and 30 decreased in expression) [Fig. 2b]. The initial 88 spot features and the later 56 spot features significantly differentially expressed represented only 14 and five unique genes represented by duplicate spot features (P ≤ 0.05).
Fold changes for significantly differentially expressed genes on the microarrays ranged from a decrease of 4.38 fold (colony stimulating factor 2 receptor, alpha gene [CSF2RA] in the control animal group after 12 h stimulation with bovine tuberculin) to an increase in expression of 5.58 fold (the major histocompatability complex, class I, A gene [BOLA] in the BTB-infected animal group after 3 h stimulation with bovine tuberculin).
Inversely correlated gene expression changes between BTB-infected and control PBMC after 3 h stimulation with bovine tuberculin
Microarray gene expression fold change values for 18 significant genes within the non-infected control animal group (n = 6) and the BTB-infected animal group (n = 6) between T0 and T3.
BOTL clone ID
Gene ontology (GO) function/s
Fold change infected cattle
Fold change control cattle
Amino acid metabolism
+1.72 ± 0.20
-1.85 ± 0.26
Glioma tumor suppressor candidate region gene 2
+1.98 ± 0.24
-1.92 ± 0.34
Transducin-like enhancer of split 3
Regulation of gene transcription
+1.69 ± 0.13
-1.86 ± 0.30
Dual specificity phosphatase 10
Phosphatase activity and response to stress
-1.44 ± 0.18
+1.57 ± 0.25
v-yes-1 Yamaguchi sarcoma viral related oncogene homolog
Amino acid phosphorylation and receptor signalling
+1.98 ± 0.21
-1.70 ± 0.34
+1.80 ± 0.19
-1.73 ± 0.16
Synaptotagmin binding, cytoplasmic RNA interacting protein
-1.36 ± 0.09
+1.58 ± 0.23
Protein binding and alternative splicing
+1.63 ± 0.21
-1.47 ± 0.52
GTPase, IMAP family member 1
Control of cell survival and nucleotide binding
+1.40 ± 0.18
-1.49 ± 0.13
TAF6 RNA polymerase II, TATA box binding protein (TBP)-associated factor
Regulation of transcription
+1.87 ± 0.31
-1.42 ± 0.17
+1.52 ± 0.13
-1.38 ± 0.12
Similar to Host cell factor C1
Regulation of transcription
+1.78 ± 0.15
-1.76 ± 0.21
GATA binding protein 4
GATA binding protein 4 (GATA4)
Regulation of transcription
+2.49 ± 0.29
-3.56 ± 0.20
Colony stimulating factor 2 receptor, alpha
colony stimulating factor 2 receptor, alpha)
Production, differentiation, and function of granulocytes and macrophages
+3.76 ± 0.41
-2.46 ± 0.26
Splicing factor, arginine/serine-rich 2
NBFGC_BE721857 chemokine (C-C motif) receptor 7 (CCR7)
+2.52 ± 0.24
-1.63 ± 0.19
AXL receptor tyrosine kinase
NBFGC_BE722178 AXL receptor tyrosine kinase
-1.36 ± 0.13
+1.87 ± 0.31
Eukaryotic translation elongation factor 1 gamma
NBFGC_BF230159 EEF1G eukaryotic translation elongation
+2.24 ± 0.15
-1.46 ± 0.24
Chemokine (C-C motif) ligand 1
Small inducible cytokine A1
+5.13 ± 0.34
-2.57 ± 0.36
The expression of 15 of these 18 genes was increased in the BTB-infected animals in response to bovine tuberculin stimulation relative to the unstimulated samples. All 18 genes were expressed in the opposite direction in the control animal group in response to bovine tuberculin stimulation, suggesting a common, but reversed mechanism affecting the regulation of expression or transcriptional trajectory of these genes in the two groups (Fig. 3 and Table 1). Among the 15 genes with significantly increased mRNA levels in PBMC from BTB-infected animals, and significantly decreased expression in PBMC from the control animals post-stimulation, were transcription factor genes (TLE3, TAF6, HCFC1 and GATA4) and genes coding for proteins involved in mRNA processing (SFRS2), nucleotide binding (GIMAP1), amino acid metabolism (SDS), amino acid phosphorylation (LYN), protein binding (PKM2), and translational elongation (EEF1G). Other genes with well-characterised roles in the immune response were also differentially expressed. These included the colony stimulating factor 2 receptor, alpha gene (CSF2RA), involved in the production, differentiation and function of myeloid cells, which was increased by 3.76 fold (P = 0.0057). The expression of the chemokine (C-C motif) ligand 1 gene (CCL1) was increased by 5.13 fold (P = 0.0009). The protein encoded by CCL1 is secreted by activated T cells and displays chemotactic activity for monocytes . Three genes, an RNA binding gene (synaptotagmin binding, cytoplasmic RNA interacting protein–SYNCRIP), a phosphatase activity and stress response gene (dual specificity phosphatase 10–DUSP10) and a signal transduction gene (AXL receptor tyrosine kinase–AXL) were all reduced in their expression in PBMC from BTB-infected animals, but significantly increased in PBMC from control animals in response to bovine tuberculin stimulation (Table 1).
Rapid but transient proinflammatory response to bovine tuberculin in PBMC from BTB-infected animals detected at 3 h post stimulation
Of the 1,391 genes on the BOTL-5 microarray, 250 spot features showed significant differential expression between the BTB-infected animals at 0 h and 3 h post stimulation at the P ≤ 0.05 level (see Additional file 1). Conversely, only 80 spot features were differentially expressed between the later two time points (Fig. 2a). It is clear from the number of genes differentially expressed in PBMC from BTB-infected samples that the majority of differential gene expression occurred early, within 3 h of stimulation with tuberculin antigens. This response was both immediate and transient (Fig. 2a). In contrast, the magnitude of gene expression changes in control animals was significantly different than in the PBMC from BTB-infected samples in response to bovine tuberculin stimulation. Furthermore, the profile of expression indicated that the majority of differential expression occurred between 0 h and 12 h in PBMC from the control animals (Fig. 2b).
A number of genes with well characterised roles in infection and immunity were found to be significantly differentially expressed in the BTB-infected animal group after stimulation with bovine tuberculin. Some of those that were significantly increased in expression by 3 h were the fibroblast growth factor receptor genes (FGFR1) and (FGFR2) [2.27 and 1.91 fold respectively], the colony stimulating factor 2 receptor, alpha gene (CSF2RA) [3.76 fold], the lymphotoxin alpha gene (LTA) [1.76 fold], and a number of kinase genes: AKT1 (2.38 fold), MAP2K1 (1.59 fold), MAP2K7 (2.58 fold), MAPKAPK2 (1.71 fold), MAP4K2 (1.62 fold), MLK3 (3.80 fold), PIK3R5 (1.53 fold) and KSR1 (2.44 fold). The results from the microarray data also highlighted the involvement of a number of genes encoding components of the TLR signalling pathway including CD14 (2.48 fold), IKBKB (1.98 fold) and NFKB1 (2.32 fold), all of which were significantly increased in expression (P ≤ 0.05). A significant increase in a gene involved in antigen presentation (BOLA major histocompatibility complex-class I, A–BOLA-A) was also detected (5.43 fold, P ≤ 0.05).
Among those differentially expressed genes that were expressed at lower levels at 3 h post-stimulation with bovine tuberculin were a number of BOTL clones with homology to known genes. These included PPIA (-1.68 fold) and the protein phosphatase genes, PPP2CA (-1.44 fold), PPP6C (-1.41 fold) and PTPRF (-1.26 fold).
In contrast, 88 gene spot features were significantly differentially expressed in response to bovine tuberculin stimulation in the control animals, 41 of these with increased expression levels and 47 expressed at lower levels (P ≤ 0.05, Fig. 2b). Only 14 genes within this group were represented by duplicate spot features. Among the genes that were significantly increased in expression in response to bovine tuberculin stimulation were AXL (1.87 fold) and MAPK13 (1.84 fold). Genes with significantly decreased expression in response to stimulation were GATA4 (-1.83 fold), MIF (-1.40 fold), and CCR7 (-1.63 fold).
Real-time quantitative reverse transcription PCR (qRT-PCR) supports a differential role for TLR-associated molecules in the early response of BTB-infected animals to bovine tuberculin
Gene expression fold change differences for the BTB-infected animal group (n = 8) between T3 and T0 validated using real time qRT-PCR.
Gene ontology (GO) function/s
Fold change difference
+16.62 ± 8.74
V-akt murine thymoma viral oncogene homolog 1
Protein kinase activity
+4.12 ± 0.46
Inhibitor of kappa light polypeptide gene enhancer in B-cells, kinase beta
Transcription activator activity
+3.64 ± 0.61
Nuclear receptor co-repressor 1
+2.80 ± 0.46
Interleukin 1 receptor antagonist
+2.79 ± 0.60
Toll interacting protein
+2.33 ± 0.34
Lymphotoxin alpha (TNF superfamily, member 1)
+2.27 ± 0.75
Toll-like receptor 4
Bacterial binding and signal transduction
+2.24 ± 0.36
Mitogen-activated protein kinase-activated protein kinase 2
+2.12 ± 0.28
+2.09 ± 0.19
+1.73 ± 0.22
Nuclear factor κβ
+1.40 ± 0.42
Conserved helix-loop-helix ubiquitous kinase
+1.38 ± 0.09
Toll-like receptor adaptor molecule 2
+1.28 ± 0.07
+1.24 ± 0.09
Transforming growth factor, beta 1
+1.19 ± 0.34
Tumor necrosis factor (TNF superfamily, member 2)
+0.94 ± 0.48
Interleukin 1 receptor 2
-1.56 ± 0.63
Gene expression fold change differences for the control non-infected animal group (n = 8) between T3 and T0 validated using real time qRT-PCR.
Gene ontology (GO) function/s
Fold change difference
Toll interacting protein
+4.32 ± 1.87
Interleukin-1 receptor-associated kinase 1
Transcription activator activity
+3.89 ± 1.91
+2.00 ± 0.69
Interleukin 1 receptor, type 1
Protein binding and signal transduction
+1.61 ± 0.55
Interleukin 1 receptor antagonist
Protein binding and signal transduction
+1.53 ± 0.62
Toll-like receptor 2
Bacterial binding and signal transduction
+1.37 ± 0.06
Myeloid differentiation primary response gene (88)
+1.23 ± 0.78
+0.83 ± 0.45
Toll interacting protein (encoded by TOLLIP) is a negative regulator of the TLR pathway and has been characterised in humans and mice . When examined by real time qRT-PCR, it was found to be significantly increased in expression in PBMC of both BTB-infected animal (Fig. 4) and the control samples (Fig. 5) after bovine tuberculin stimulation by 4.32 and 2.33 fold respectively (P = 0.0094 and P = 0.0002). Expression levels of the gene encoding a TLR4-specific molecule, Toll-like receptor adaptor molecule 2 (TICAM2) , were estimated in both groups and although not differentially expressed in the control animal group (data not shown), it was significantly increased by 1.28 fold in the BTB-infected animal group (P = 0.0006, Fig. 4).
The gene encoding interleukin-10 (IL10) was significantly differentially expressed in both groups (Figs. 4 and 5). In the BTB-infected animal group, IL10 expression was increased by 1.24 fold (P = 0.0135). However, in the control animal group, IL10 was significantly increased in expression by 2.00 fold (P = 0.0091). In addition, it was noted that neither NFKB nor IFNG were differentially expressed in the control animal group when tested by real time qRT-PCR (data not shown).
Bovine tuberculosis is the fourth most important livestock disease worldwide . The benefits of developing, applying and maintaining improved control and eradication strategies for BTB are manifold, and directly impact on human and animal health . The specific immune cell signalling pathways involved in the immune response to BTB are highly complex and poorly characterised in cattle. This has obvious limitations for the understanding of, and design of improved diagnostics and effective therapeutics. However, studies on tuberculosis in the human and murine models have highlighted the involvement of cell regulatory signalling pathways in the immune response that are also likely to be relevant in BTB. Whereas traditionally, studies of BTB have focused on adaptive immunity, the findings from these studies are pointing toward a critical role for signalling via the innate immune system, including TLRs in initiating and directing the subsequent immune response and determining the outcome of infection [8–10].
The timing and potency of the cellular and immunological events that occur immediately post-infection are crucial determinants governing infection . Pathogen-induced phenotypic changes in host cells are often accompanied by marked changes in gene expression due to host- and/or pathogen mediated reprogramming of the transcriptome during infection . Previous work by our group compared the gene expression differences between bovine tuberculin-stimulated and non-stimulated PBMC from BTB-infected animals  and demonstrated that stimulation with bovine tuberculin induced significant gene expression changes that can be useful for dissection of the immune response to BTB. Subsequently we identified a novel gene expression profile indicative of innate immune gene repression in heavily infected cattle in vivo . Expression clustering of these data yielded a gene infection signature for disease and highlighted genes and regulatory pathways, including the TLR cell signalling pathway .
In the present study we have shown that after overnight recovery in vitro, PBMC from BTB-infected cattle are significantly more responsive to bovine tuberculin stimulation than control animal PBMC. Gene expression levels were estimated in PBMC from BTB-infected and healthy controls either non-stimulated or at 3 h and 12 h after bovine tuberculin stimulation using a common reference microarray approach (Fig. 1). Significant gene expression changes were observed in response to bovine tuberculin stimulation for both animal groups over a 12 h time course. The microarray data indicated that substantially different gene expression profiles were evident in the BTB-infected animals relative to the control animals (Fig. 2a and Fig. 2b). Following a 12 h incubation with bovine tuberculin, these analyses also showed that the immune response in the infected animal group was both rapid and transient (Fig. 2a). Analysis of gene expression differences across the time course showed that differences between groups were most evident in the period between 0 h and 3 h after bovine tuberculin stimulation–indicating that the early response to bovine tuberculin is substantially different between the two animal groups.
Comparative analysis of the PBMC gene expression program in response to bovine tuberculin identified a panel of 18 genes that were significantly differentially expressed in both animal groups. Interestingly, all of these genes were expressed in opposite directions in the two groups. Expression of 15 of the 18 genes was increased in PBMC from BTB-infected animals including genes encoding proteins involved with transcriptional regulation (TLE3, TAF6, HCFC1 and GATA4), a chemokine (CCL1), and a chemokine receptor(CSF2RA). In contrast, only three genes were decreased in expression in the BTB-infected group. Conversely, the opposite pattern was observed in PBMC from the control animals with 15 genes decreased, and three genes increased in expression relative to the BTB-infected group (Fig. 3 and Table 1). These data suggest a number of important gene targets for further study, as well as identifying cell regulatory pathways that may be differentially regulated in BTB-infected animals.
Of the 250 microarray spot features that were differentially expressed in the BTB-infected animals, those displaying increased expression in response to bovine tuberculin outnumbered those with decreased expression by a factor of two. We have previously shown that the expression of a number of the genes represented by these spot features, including TLR2, TLR4 and NFKB, was significantly repressed in heavily infected cattle in vivo . The transformation associated with gene repression in vivo to gene activation detected in the same animals in vitro was evident in the reversal in direction of expression of a number of well characterised genes including those encoding TLRs, MHC molecules, and cytokines. Furthermore, the 41 genes examined by real time qRT-PCR confirmed the BOTL-5 microarray results and supported a trend towards a proinflammatory immune response to bovine tuberculin in PBMC from BTB-infected cattle. Expression of a gene encoding the key transcription factor and mediator of the immune response (NFKB1) indicated a distinctive proinflammatory gene expression program characterised by a significant two-fold increase in IL8 expression and 16-fold increase in IFNG expression (Fig. 4 and Table 2).
While the differentially expressed genes detected in this study provide evidence of a pre-existing gene expression program most likely caused as a result of M. bovis infection in the diseased cattle, other confounding factors could be responsible for some of the changes detected. Changes or fundamental differences in cell subpopulations between animals and the separation of PBMC from an in vivo immunosuppressive environment  into culture media may have affected some gene expression patterns. Although there was no evidence at post-mortem examination of clinical disease in the cattle caused by other infectious agents, the presence of undetected pathogens may have generated host immune responses in either group of cattle that could have accounted for some of the transcriptional changes detected. In addition, it is important to acknowledge that the control and infected cattle were sampled from different herds, and that this may represent a confounding factor in the analysis of gene expression differences between the groups.
Genes encoding adaptor and mediator molecules of the TLR activation pathway were also profiled by real time qRT-PCR to examine cellular pathways contributing to the differential response between BTB-infected and control animal groups. The gene encoding the Toll-interacting protein (TOLLIP), a negative regulator of TLR signalling  was examined and gene expression was significantly increased in both the control and BTB-infected animal groups (4.32 fold and 2.33 fold, respectively) in response to bovine tuberculin stimulation. TOLLIP has previously been found to impair TLR-2 and TLR-4 activation of NF-κB ; therefore, it might be expected to prevent the downstream signalling from TLR-2 and TLR-4 in healthy control and BTB-infected cattle samples respectively. However, significant increased expression of NFKB1 and genes encoding mediators of NF-κB (CHUK and IKBKB, Fig. 4 and Table 2), coupled with differential expression of over 250 spot features on the microarrays suggested that proinflammatory signalling was not inhibited in the BTB-infected animal samples. Interestingly, another TLR-4-specific accessory molecule, Toll-like receptor adaptor molecule 2 (TICAM2), discovered in studies of M. tuberculosis infection in mice can bypass the inhibitory effects of TOLLIP to transmit cell signals  leading to increased expression of the NFKB1 gene . Expression of the TICAM2 gene was examined using real time qRT-PCR and was found to be increased in the BTB-infected animal samples only (1.28 fold, Fig. 4 and Table 2). The increased expression of TICAM2 may possibly provide a mechanism through which proinflammatory gene activation is regulated in BTB-infected animals (Fig. 2a).
Our results are consistent with proposed mechanisms that drive an ineffective TH2 type response and contribute to the outcome of BTB infection in cattle. It has been suggested that the TH1/TH2 bias of the immune response can be determined by specific TLRs . Furthermore, TLR-2 activation is a less efficient method of proinflammatory gene activation and may play a role in TLR-2 mediated immunosuppression .
In the present study, increased expression of the IL10 and TGFB1 genes was detected in the BTB-infected group in response to bovine tuberculin stimulation (Fig. 4) and elsewhere this has been associated with decreased ability of PBMC and macrophages to restrict mycobacterial growth in both humans [38, 39] and mice [40, 41]. Furthermore, IL-10 has been implicated in the suppression of the proinflammatory immune response and subversion of the host bacteriocidal immune response . The results of the present study suggest that the 16-fold increased expression of the IFNG gene (Table 2) may drive the changes in gene expression detected in the BTB-infected animal group. In addition, although the IL10 gene is significantly increased in expression by 1.23 fold (Table 2), this may be insufficient to compete with the proinflammatory effects of IFN-γ. In this regard, it has been previously shown that the ratio of pro- and anti-inflammatory cytokines will determine the overall outcome of the immune response and subsequent correlated gene activation and/or repression [43, 11].
In many countries the presence of M. bovis-infected wildlife can act as reservoirs of BTB infection for livestock and there is increased information available on the host response to infection in a variety of these species [44, 45]. In our previous study , transcriptional profiling of infected and non-infected control animals in the absence of exogenous antigen stimulation demonstrated decreased expression of MHC class II molecules, and similar findings have been reported in deer in response to natural TB infection .Thacker and colleagues characterised the immunological responses of peripheral blood leukocytes (PBL) from M. bovis-infected and non-infected white-tailed deer to infection by monitoring cytokine gene expression after exogenous antigen stimulation . The infected deer displayed a significant 75-fold increase in the expression of the proinflammatory cytokine, IFN-γ, comparable to the response detected in cattle during the present study. One notable difference was that the increase in IL10 and TGFB1 gene expression in response to bovine tuberculin stimulation of PBMC from infected cattle was not detected in infected deer. However, post-stimulation time points differed between the two studies and this could account for the differences observed.
The results of the present study are consistent with work carried out on Johne's disease in cattle caused by M. avium subsp. paratuberculosis (MAP) [47, 48]. Those studies also detected a novel gene expression profile in PBMC from MAP-infected animals that was both rapid and transient across time [47, 48]. In a separate study, stimulation of PBMC with MAP was shown to suppress the proinflammatory immune response , and the authors also found evidence of a dynamic TH1/TH2 type response, which eventually gives way to a predominantly TH2 like response in MAP-infected animals . Such a temporal expression pattern, where peak production of IL-10 lags behind that of IFN-γ, has also been demonstrated in relation to MAP stimulation of bovine PBMC . An important outcome from these studies involving both MAP-infected and BTB-infected cattle, is that the host gene expression profiles resulting from infection by these two related mycobacteria are remarkably similar.
The outcome of infection with M. bovis is determined by a complex and dynamic interaction between the host immune system and the pathogen . The cellular signalling events governing the immune responses that form the basis of cell-mediated immune BTB diagnostic tests remain incompletely understood. We hypothesized that the PBMC from BTB-infected cattle would display a distinct gene expression program resulting, from previous exposure to M. bovis. The gene expression program observed in PBMC from BTB-infected cattle was substantially different than in PBMC from control animals suggesting that BTB infection can modulate the immune response to bovine tuberculin. To the best of our knowledge, this work is also among the first to report the involvement of TLRs and TLR accessory molecules in this immune response.
The results demonstrate that PBMC from BTB-infected animals are highly responsive to bovine tuberculin stimulation in vitro, due possibly to their removal from an immunosuppressive environment in vivo . Specific genes are commonly involved in the differential response of both control and BTB-infected animals indicating coordinated regulation of the innate immune response to antigen stimulation. The differential responses of a number of genes involved in cell signalling pathways regulating both the innate and adaptive immune response is indicative of a pre-existing gene expression program. The results are in agreement with other human and mouse models using M. tuberculosis [52–55] and bovine studies using MAP. The specific genes activated in response to mycobacterial infection have yet to be fully elucidated, but the overall response profiles seem to be similar among different mammalian hosts (Homo sapiens, Mus musculus and Bos taurus) infected with mycobacterial pathogens (M. tuberculosis, M. bovis and M. avium subsp. paratuberculosis). The results also demonstrate the usefulness of employing the natural host for M. bovis infection as a model to investigate the immune response to tuberculosis using functional genomics technologies.
Future work will concentrate on the elucidation of cell signalling pathways detected during the immune response to bovine tuberculin and their potential roles in the immune repression detected in vivo . It is clear that a pre-existing immune program in PBMC from BTB-infected animals influences the response to bovine tuberculin stimulation and may affect the outcome of infection. Understanding the causes and consequences of this novel gene expression program should ultimately lead to a more complete understanding of the immune response to BTB and inform the development of novel or improved future diagnostics and vaccines.
Experimental animals and infection status
Sixteen cattle were used for this study. The eight infected animals were chosen from herds with a history of M. bovis infection. The animals were selected on the basis of the skin-fold thickness response to bovine and avian tuberculin in the single intradermal comparative tuberculin test (SICTT). The SICTT reactor animals were selected where the skin-fold thickness response to PPD-bovine exceeded that of PPD-avian by at least 12 mm. All of these animals were also measured positive in a whole blood IFN-γ assay . The cattle were confirmed positive for tuberculosis following detailed post-mortem pathological examination and/or culture. There was no evidence at post-mortem examination that any of the BTB-infected animals had clinical diseases caused by other infectious agents. Bronchial, mediastinal, submandibular, retropharyngeal, mesenteric and hepatic lymph nodes and lungs were examined macroscopically for tuberculosis lesions. Suspected lesions were cultured on Stonebrinks and Lowenstein-Jensen media at 37°C for eight weeks to detect M. bovis . The eight non-infected control animals were 2–3 year old unrelated females selected from a Holstein Friesian herd without a recent history of tuberculosis and were SICTT and IFN-γ test negative. This control animal group also underwent comprehensive testing for the following infections: Brucellosis (Brucella abortus), Johne's disease (Mycobacterium avium subsp. paratuberculosis), infectious bovine rhinotracheitis (bovine herpesvirus-1), salmonella (Salmonella typhimurium) and bovine viral diarrhoea (Pestivirus).
Blood sampling and analysis
400 ml of blood was collected from each animal in sterile heparinised bottles. Five ml of blood was used for haematological analysis and leukocyte cell population subsets were compared between infected and control groups as previously described .
PBMC separation, culture, RNA extraction and quality control
PBMC were isolated using the Percoll™ gradient method as previously described . PBMC were seeded at 107 per flask and cultured in RPMI 1640 culture medium supplemented with 5% FBS, 0.1% mercaptoethanol and 0.1% gentamicin. A total of 48 tissue culture flasks represented 16 individual PBMC samples (BTB-infected [n = 8] and healthy control [n = 8]) per time point at 0 h (pre-stimulation), 3 h and 12 h post-stimulation with bovine tuberculin. Previous work in our laboratory using PBMC from BTB-infected cattle indicated that these time points would be the most appropriate for further analysis . All PBMC samples were cultured overnight at 37°C in 5% CO2 to minimize noise in gene expression measurements potentially introduced by the mechanical disruption of cells associated with PBMC isolation. After cells were harvested for time point 0, remaining flasks were stimulated using 50 μg/ml bovine tuberculin and incubated for 3 h and 12 h at 37°C in 5% CO2. Residual cells not seeded for culture in either treatment were immediately dissolved in 3 ml TriReagent® (Molecular Research Centre Inc., Cincinnati, OH) and frozen in 1.5 ml cryotubes at -80°C for use as a common reference RNA (CRR) pool. Total RNA was extracted from PBMC harvested after 3 h and 12 h post-stimulation using a combined TriReagent®, DNase treatment and Qiagen RNeasy® method (Qiagen Ltd., Crawley, UK) according to the manufacturers' instructions. The integrity and stability of RNA samples is crucial for gene expression analyses using microarray technology; therefore, RNA yield and quality were assessed using an Agilent 2100 Bioanalyzer (Agilent Technologies Ireland Ltd., Dublin, Ireland). The two-step method for RNA extraction described above was found to produce RNA of high yield and quality (ratios of 18S to 28S ribosomal RNA averaged > 1.6).
Microarray experimental design
The 3,888 feature BOTL-5 immunogenetic microarray system used has been described previously . The NCBI GEO platform accession for the BOTL-5 microarray is GPL5751. The immunobiology-targeted BOTL-5 array contains 1,391 genes or ESTs spotted in duplicate with multiple additional control features (blank spots, negative spots, housekeeper genes) and is an expanded version of the BOTL-4 array described previously by our group [60, 22]. A reference design was used for microarray hybridizations, such that all RNA samples were labelled using Cy3 and co-hybridized with Cy5 labelled common reference RNA (CRR) pool as described previously . Thirty-six arrays were hybridized in total, representing six individual animal PBMC samples from each treatment group pre-stimulation, and at 3 h and 12 h post-stimulation with bovine tuberculin, as shown in Fig. 1. It was hypothesized that the CRR pool would display similar mRNA expression levels and gene coverage as the target samples, therefore allowing flexible, accurate and consistent comparison of gene expression data across a time course without arbitrarily pairing animals from different groups . The CRR pool contained equal amounts of total RNA from the treated and control animal groups.
cDNA labelling, hybridisation and scanning
cDNA synthesis, Cy3 and Cy5 labelling and microarray hybridizations were performed as previously described . Each labelling reaction contained a total of 8 μg total RNA per sample and 10 μg total RNA from the CRR. Purified labelled samples were combined (either an infected or a control sample combined with a CRR sample) and co-hybridized on the BOTL-5 microarrays using SlideHyb Glass Array Hybridization Buffer #3 (Ambion Ltd.). Microarray hybridizations were performed using a Tecan HS400 hybridisation station (Tecan Ltd.) with the protocol as previously described . Microarrays were scanned immediately using a GenePix 4000B microarray scanner (Molecular Devices Ltd.).
Data processing, normalization and analysis
The working signal intensities were generated using the mean foreground intensity values minus the median background intensity values as outputted from the GenePix Pro 5.0 results file. Two methods of data pre-processing were used to flag unreliable data as previously described .
Median-based normalization, which corrects the data such that all arrays have the same median  was used as previously described . Microarray data analysis (including analyses of the microarray-platform specific false discovery rate) was carried out using class comparisons between experimental groups (parametric t-tests) as implemented in BRB ArrayTools version 3.0  as previously described . Each gene on the BOTL array is represented by duplicate gene spot features. Analyses of the false discovery rate (FDR) according to methods described previously  demonstrated that the differentially expressed gene spot feature lists were robust and reliable. Therefore, to maximise information from these experiments we refer to numbers of differentially expressed gene spot features in the Results section.
Real time qRT-PCR validation of within group differential gene expression profiles
Replicate spot features on the BOTL array were used as a check for the quality control of gene expression data. Each spot was analyzed individually thereby allowing the individual genes to be flagged if expression results from two or more replicates were statistically different. This enabled the identification of differentially expressed genes that had a low probability of being false positives and expedited the choice of target genes for real time qRT-PCR validation of the microarray results.
Real time qRT-PCR primer sequences, optimum primer concentrations and amplicon sizes for all validated genes.
Forward primer (5'-3')
Reverse primer (5'-3')
Amplicon size (bp)
Primer conc. (nM)
We would like to thank the staff at the UCD Lyons Research Farm for assistance with cattle blood sampling and Grace O'Gorman, Mairead Doyle, and Tara Fitzsimons for technical assistance and advice. We also thank Prof. Paul Coussens of the MSU Centre for Animal Functional Genomics for provision of BOTL-5 microarrays. This work was supported by an Investigator Grant from Science Foundation Ireland (Grant No: SFI/01/F.1/B028).
- Phillips CJ, Foster CR, Morris PA, Teverson R: The transmission of Mycobacterium bovis infection to cattle. Res Vet Sci. 2003, 74 (1): 1-15. 10.1016/S0034-5288(02)00145-5.PubMedView Article
- Neill SD, Cassidy J, Hanna J, Mackie DP, Pollock JM, Clements A, Walton E, Bryson DG: Detection of Mycobacterium bovis infection in skin test-negative cattle with an assay for bovine interferon-gamma. Vet Rec. 1994, 135 (6): 134-135.PubMedView Article
- Gormley E, Doyle MB, McGill K, Costello E, Good M, Collins JD: The effect of the tuberculin test and the consequences of a delay in blood culture on the sensitivity of a gamma-interferon assay for the detection of Mycobacterium bovis infection in cattle. Vet Immunol Immunopathol. 2004, 102 (4): 413-420. 10.1016/j.vetimm.2004.08.002.PubMedView Article
- Lepper AW, Pearson CW, Corner LA: Anergy to tuberculin in beef cattle. Aust Vet J. 1977, 53 (5): 214-216. 10.1111/j.1751-0813.1977.tb00188.x.PubMedView Article
- Aldwell FE, Wedlock DN, Slobbe LJ, Griffin JF, Buddle BM, Buchan GS: In vitro control of Mycobacterium bovis by macrophages. Tuberculosis (Edinb). 2001, 81 (1–2): 115-123. 10.1054/tube.2000.0280.View Article
- Pollock JM, Neill SD: Mycobacterium bovis infection and tuberculosis in cattle. Vet J. 2002, 163 (2): 115-127. 10.1053/tvjl.2001.0655.PubMedView Article
- Pollock JM, Welsh MD, McNair J: Immune responses in bovine tuberculosis: towards new strategies for the diagnosis and control of disease. Vet Immunol Immunopathol. 2005, 108 (1–2): 37-43. 10.1016/j.vetimm.2005.08.012.PubMedView Article
- Werling D, Jungi TW: TOLL-like receptors linking innate and adaptive immune response. Vet Immunol Immunopathol. 2003, 91 (1): 1-12. 10.1016/S0165-2427(02)00228-3.PubMedView Article
- Werling D, Coffey TJ: Pattern recognition receptors in companion and farm animals – The key to unlocking the door to animal disease?. Vet J. 2006
- Werling D, Piercy J, Coffey TJ: Expression of TOLL-like receptors (TLR) by bovine antigen-presenting cells-potential role in pathogen discrimination?. Vet Immunol Immunopathol. 2006, 112 (1–2): 2-11. 10.1016/j.vetimm.2006.03.007.PubMedView Article
- Thacker TC, Palmer MV, Waters WR: Associations between cytokine gene expression and pathology in Mycobacterium bovis infected cattle. Vet Immunol Immunopathol. 2007, 119 (3–4): 204-213. 10.1016/j.vetimm.2007.05.009.PubMedView Article
- McGuire K, Glass EJ: The expanding role of microarrays in the investigation of macrophage responses to pathogens. Vet Immunol Immunopathol. 2005, 105 (3–4): 259-275. 10.1016/j.vetimm.2005.02.001.PubMedView Article
- Belcher CE, Drenkow J, Kehoe B, Gingeras TR, McNamara N, Lemjabbar H, Basbaum C, Relman DA: The transcriptional responses of respiratory epithelial cells to Bordetella pertussis reveal host defensive and pathogen counter-defensive strategies. Proc Natl Acad Sci USA. 2000, 97 (25): 13847-13852. 10.1073/pnas.230262797.PubMedPubMed CentralView Article
- Cohen P, Bouaboula M, Bellis M, Baron V, Jbilo O, Poinot-Chazel C, Galiegue S, Hadibi EH, Casellas P: Monitoring cellular responses to Listeria monocytogenes with oligonucleotide arrays. J Biol Chem. 2000, 275 (15): 11181-11190. 10.1074/jbc.275.15.11181.PubMedView Article
- Detweiler CS, Cunanan DB, Falkow S: Host microarray analysis reveals a role for the Salmonella response regulator phoP in human macrophage cell death. Proc Natl Acad Sci USA. 2001, 98 (10): 5850-5855. 10.1073/pnas.091110098.PubMedPubMed CentralView Article
- Kagnoff MF, Eckmann L: Analysis of host responses to microbial infection using gene expression profiling. Curr Opin Microbiol. 2001, 4 (3): 246-250. 10.1016/S1369-5274(00)00198-3.PubMedView Article
- Coussens PM, Colvin CJ, Wiersma K, Abouzied A, Sipkovsky S: Gene expression profiling of peripheral blood mononuclear cells from cattle infected with Mycobacterium paratuberculosis. Infect Immun. 2002, 70 (10): 5494-5502. 10.1128/IAI.70.10.5494-5502.2002.PubMedPubMed CentralView Article
- Wang JP, Rought SE, Corbeil J, Guiney DG: Gene expression profiling detects patterns of human macrophage responses following Mycobacterium tuberculosis infection. FEMS Immunol Med Microbiol. 2003, 39 (2): 163-172. 10.1016/S0928-8244(03)00223-2.PubMedView Article
- Xu Y, Xie J, Li Y, Yue J, Chen J, Chunyu L, Wang H: Using a cDNA microarray to study cellular gene expression altered by Mycobacterium tuberculosis. Chin Med J (Engl). 2003, 116 (7): 1070-1073.
- Blumenthal A, Lauber J, Hoffmann R, Ernst M, Keller C, Buer J, Ehlers S, Reiling N: Common and unique gene expression signatures of human macrophages in response to four strains of Mycobacterium avium that differ in their growth and persistence characteristics. Infect Immun. 2005, 73 (6): 3330-3341. 10.1128/IAI.73.6.3330-3341.2005.PubMedPubMed CentralView Article
- Staudt LM, Brown PO: Genomic views of the immune system. Annu Rev Immunol. 2000, 18: 829-859. 10.1146/annurev.immunol.18.1.829.PubMedView Article
- Meade KG, Gormley E, Park SD, Fitzsimons T, Rosa GJ, Costello E, Keane J, Coussens PM, MacHugh DE: Gene expression profiling of peripheral blood mononuclear cells (PBMC) from Mycobacterium bovis infected cattle after in vitro antigenic stimulation with purified protein derivative of tuberculin (PPD). Vet Immunol Immunopathol. 2006, 113 (1–2): 73-89. 10.1016/j.vetimm.2006.04.012.PubMedView Article
- Rhodes SG, Buddle BM, Hewinson RG, Vordermeier HM: Bovine tuberculosis: immune responses in the peripheral blood and at the site of active disease. Immunology. 2000, 99 (2): 195-202. 10.1046/j.1365-2567.2000.00944.x.PubMedPubMed CentralView Article
- Meade KG, Gormley E, Doyle MB, Fitzsimons T, O' Farrelly C, Costello E, Keane J, Zhao Y, MacHugh DE: Innate gene repression associated with Mycobacterium bovis infection in cattle: toward a gene signature of disease. BMC Genomics. 2007, 8 (1): 400-10.1186/1471-2164-8-400.PubMedPubMed CentralView Article
- Dobbin K, Shih JH, Simon R: Questions and answers on design of dual-label microarrays for identifying differentially expressed genes. J Natl Cancer Inst. 2003, 95 (18): 1362-1369.PubMedView Article
- Simon RM: Design and analysis of DNA microarray investigations. 2003, New York; London: Springer
- Simon R: Microarray-based expression profiling and informatics. Curr Opin Biotechnol. 2008, 19 (1): 26-29. 10.1016/j.copbio.2007.10.008.PubMedPubMed CentralView Article
- Barrett T, Troup DB, Wilhite SE, Ledoux P, Rudnev D, Evangelista C, Kim IF, Soboleva A, Tomashevsky M, Edgar R: NCBI GEO: mining tens of millions of expression profiles–database and tools update. Nucleic acids research. 2007, D760-765. 10.1093/nar/gkl887. 35 Database
- Miller MD, Krangel MS: The human cytokine I-309 is a monocyte chemoattractant. Proc Natl Acad Sci USA. 1992, 89 (7): 2950-2954. 10.1073/pnas.89.7.2950.PubMedPubMed CentralView Article
- Zhang G, Ghosh S: Negative regulation of toll-like receptor-mediated signaling by Tollip. J Biol Chem. 2002, 277 (9): 7059-7065. 10.1074/jbc.M109537200.PubMedView Article
- Yamamoto M, Sato S, Hemmi H, Uematsu S, Hoshino K, Kaisho T, Takeuchi O, Takeda K, Akira S: TRAM is specifically involved in the Toll-like receptor 4-mediated MyD88-independent signaling pathway. Nat Immunol. 2003, 4 (11): 1144-1150. 10.1038/ni986.PubMedView Article
- Perry BD, Randolph TF, McDermott JJ, Sones KR, Thornton PK: Investing in animal health research to alleviate poverty. 2002, Kenya: International Livestock Research Institute (ILRI)
- Gormley E: Diagnosis of Mycobacterium bovis infection in cattle. 27th International Dairy Federation, World Dairy Conference: 2007; Shanghai, China. 2007
- Jenner RG, Young RA: Insights into host responses against pathogens from transcriptional profiling. Nat Rev Microbiol. 2005, 3 (4): 281-294. 10.1038/nrmicro1126.PubMedView Article
- Oshiumi H, Sasai M, Shida K, Fujita T, Matsumoto M, Seya T: TIR-containing adapter molecule (TICAM)-2, a bridging adapter recruiting to toll-like receptor 4 TICAM-1 that induces interferon-beta. J Biol Chem. 2003, 278 (50): 49751-49762. 10.1074/jbc.M305820200.PubMedView Article
- Netea MG, Meer Van der JW, Sutmuller RP, Adema GJ, Kullberg BJ: From the Th1/Th2 paradigm towards a Toll-like receptor/T-helper bias. Antimicrob Agents Chemother. 2005, 49 (10): 3991-3996. 10.1128/AAC.49.10.3991-3996.2005.PubMedPubMed CentralView Article
- Netea MG, Meer Van der JW, Kullberg BJ: Toll-like receptors as an escape mechanism from the host defense. Trends Microbiol. 2004, 12 (11): 484-488. 10.1016/j.tim.2004.09.004.PubMedView Article
- Hirsch CS, Yoneda T, Averill L, Ellner JJ, Toossi Z: Enhancement of intracellular growth of Mycobacterium tuberculosis in human monocytes by transforming growth factor-beta 1. J Infect Dis. 1994, 170 (5): 1229-1237.PubMedView Article
- Boussiotis VA, Tsai EY, Yunis EJ, Thim S, Delgado JC, Dascher CC, Berezovskaya A, Rousset D, Reynes JM, Goldfeld AE: IL-10-producing T cells suppress immune responses in anergic tuberculosis patients. J Clin Invest. 2000, 105 (9): 1317-1325. 10.1172/JCI9918.PubMedPubMed CentralView Article
- Murray PJ: Defining the requirements for immunological control of mycobacterial infections. Trends Microbiol. 1999, 7 (9): 366-372. 10.1016/S0966-842X(99)01567-X.PubMedView Article
- Jacobs M, Brown N, Allie N, Gulert R, Ryffel B: Increased resistance to mycobacterial infection in the absence of interleukin-10. Immunology. 2000, 100 (4): 494-501. 10.1046/j.1365-2567.2000.00053.x.PubMedPubMed CentralView Article
- Weiss DJ, Evanson OA, de Souza C, Abrahamsen MS: A critical role of interleukin-10 in the response of bovine macrophages to infection by Mycobacterium avium subsp paratuberculosis. Am J Vet Res. 2005, 66 (4): 721-726. 10.2460/ajvr.2005.66.721.PubMedView Article
- Welsh MD, Cunningham RT, Corbett DM, Girvin RM, McNair J, Skuce RA, Bryson DG, Pollock JM: Influence of pathological progression on the balance between cellular and humoral immune responses in bovine tuberculosis. Immunology. 2005, 114 (1): 101-111. 10.1111/j.1365-2567.2004.02003.x.PubMedPubMed CentralView Article
- Thacker TC, Palmer MV, Waters WR: Correlation of cytokine gene expression with pathology in white-tailed deer (Odocoileus virginianus) infected with Mycobacterium bovis. Clin Vaccine Immunol. 2006, 13 (6): 640-647. 10.1128/CVI.00024-06.PubMedPubMed CentralView Article
- Lesellier S, Corner L, Costello E, Sleeman P, Lyashchenko K, Greenwald R, Esfandiari J, Singh M, Hewinson RG, Chambers M: Antigen specific immunological responses of badgers (Meles meles) experimentally infected with Mycobacterium bovis. Vet Immunol Immunopathol. 2008, 122 (1–2): 35-45. 10.1016/j.vetimm.2007.11.005.PubMedView Article
- Fernandez de Mera IG, Perez de la Lastra JM, Ayoubi P, Naranjo V, Kocan KM, Gortazar C, de la Fuente J: Differential expression of inflammatory and immune response genes in mesenteric lymph nodes of Iberian red deer (Cervus elaphus hispanicus) naturally infected with Mycobacterium bovis. Dev Comp Immunol. 2008, 32 (2): 85-91.PubMedView Article
- Coussens PM, Colvin CJ, Rosa GJ, Perez Laspiur J, Elftman MD: Evidence for a novel gene expression program in peripheral blood mononuclear cells from Mycobacterium avium subsp. paratuberculosis-infected cattle. Infect Immun. 2003, 71 (11): 6487-6498. 10.1128/IAI.71.11.6487-6498.2003.PubMedPubMed CentralView Article
- Coussens PM, Jeffers A, Colvin C: Rapid and transient activation of gene expression in peripheral blood mononuclear cells from Johne's disease positive cows exposed to Mycobacterium paratuberculosis in vitro. Microb Pathog. 2004, 36 (2): 93-108. 10.1016/j.micpath.2003.09.007.PubMedView Article
- Coussens PM, Verman N, Coussens MA, Elftman MD, McNulty AM: Cytokine gene expression in peripheral blood mononuclear cells and tissues of cattle infected with Mycobacterium avium subsp. paratuberculosis: evidence for an inherent proinflammatory gene expression pattern. Infect Immun. 2004, 72 (3): 1409-1422. 10.1128/IAI.72.3.1409-1422.2004.PubMedPubMed CentralView Article
- Coussens PM: Model for immune responses to Mycobacterium avium subspecies paratuberculosis in cattle. Infect Immun. 2004, 72 (6): 3089-3096. 10.1128/IAI.72.6.3089-3096.2004.PubMedPubMed CentralView Article
- Koul A, Herget T, Klebl B, Ullrich A: Interplay between mycobacteria and host signalling pathways. Nat Rev Microbiol. 2004, 2 (3): 189-202. 10.1038/nrmicro840.PubMedView Article
- Reiling N, Holscher C, Fehrenbach A, Kroger S, Kirschning CJ, Goyert S, Ehlers S: Cutting edge: Toll-like receptor (TLR)2- and TLR4-mediated pathogen recognition in resistance to airborne infection with Mycobacterium tuberculosis. J Immunol. 2002, 169 (7): 3480-3484.PubMedView Article
- Doherty TM, Arditi M: TB, or not TB: that is the question – does TLR signaling hold the answer?. J Clin Invest. 2004, 114 (12): 1699-1703.PubMedPubMed CentralView Article
- Quesniaux V, Fremond C, Jacobs M, Parida S, Nicolle D, Yeremeev V, Bihl F, Erard F, Botha T, Drennan M: Toll-like receptor pathways in the immune responses to mycobacteria. Microbes Infect. 2004, 6 (10): 946-959. 10.1016/j.micinf.2004.04.016.PubMedView Article
- Takeda K, Akira S: Toll-like receptors in innate immunity. Int Immunol. 2005, 17 (1): 1-14. 10.1093/intimm/dxh186.PubMedView Article
- Rothel JS, Jones SL, Corner LA, Cox JC, Wood PR: A sandwich enzyme immunoassay for bovine interferon-gamma and its use for the detection of tuberculosis in cattle. Aust Vet J. 1990, 67 (4): 134-137. 10.1111/j.1751-0813.1990.tb07730.x.PubMedView Article
- Costello E, Quigley F, Flynn O, Gogarty A, McGuirk J, Murphy A, Dolan L: Laboratory examination of suspect tuberculous lesions detected on abattoir postmortem examination of cattle from non-reactor herds. Irish Vet J. 1998, 51 (5): 248-250.
- Ulmer AJ, Scholz W, Ernst M, Brandt E, Flad HD: Isolation and subfractionation of human peripheral blood mononuclear cells (PBMC) by density gradient centrifugation on Percoll. Immunobiology. 1984, 166 (3): 238-250.PubMedView Article
- Coussens PM, Nobis W: Bioinformatics and high throughput approach to create genomic resources for the study of bovine immunobiology. Vet Immunol Immunopathol. 2002, 86 (3–4): 229-244. 10.1016/S0165-2427(02)00005-3.PubMedView Article
- Evans AC, Ireland JL, Winn ME, Lonergan P, Smith GW, Coussens PM, Ireland JJ: Identification of genes involved in apoptosis and dominant follicle development during follicular waves in cattle. Biol Reprod. 2004, 70 (5): 1475-1484. 10.1095/biolreprod.103.025114.PubMedView Article
- Novoradovskaya N, Whitfield ML, Basehore LS, Novoradovsky A, Pesich R, Usary J, Karaca M, Wong WK, Aprelikova O, Fero M: Universal Reference RNA as a standard for microarray experiments. BMC Genomics. 2004, 5 (1): 20-10.1186/1471-2164-5-20.PubMedPubMed CentralView Article
- Quackenbush J: Microarray data normalization and transformation. Nat Genet. 2002, 32 (Suppl): 496-501. 10.1038/ng1032.PubMedView Article
- Simon R, Lam A, Li M-C, Ngan M, Menenzes S, Zhao Y: Analysis of gene expression data using BRB-Array Tools. Cancer Inform. 2007, 2: 11-17.
- Livak KJ, Schmittgen TD: Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method. Methods. 2001, 25 (4): 402-408. 10.1006/meth.2001.1262.PubMedView Article
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.