- Open Access
Gene-expression profiling of calves 6 and 9 months after inoculation with Mycobacterium avium subspecies paratuberculosis
Veterinary Research volume 45, Article number: 96 (2014)
Early detection of Johne's disease (JD) caused by Mycobacterium avium subspecies paratuberculosis (MAP) is essential to reduce transmission; consequently, new diagnostic techniques and approaches to detect MAP or markers of early MAP infection are being explored. The objective was to identify biomarkers associated with MAP infection at 6 and 9 months after oral inoculation. Therefore, gene expression analysis was done using whole blood cells obtained from MAP-infected calves. All MAP-inoculated calves had a cell-mediated immune response (IFN-γ) to Johnin PPD specific antigens, and 60% had an antibody response to MAP antigens. Gene expression analysis at 6 months after inoculation revealed downregulation of chemoattractants, namely neutrophil beta-defensin-9 like peptide (BNBD9-Like), S100 calcium binding protein A9 (s100A9) and G protein coupled receptor 77 (GPR77) or C5a anaphylatoxin chemotactic receptor (C5a2). Furthermore, BOLA/MHC-1 intracellular antigen presentation gene was downregulated 9 months after inoculation. In parallel, qPCR experiments to evaluate the robustness of some differentially expressed genes revealed consistent downregulation of BOLA/MHC-I, BNBD9-Like and upregulation of CD46 at 3, 6, 9, 12, and 15 months after inoculation. In conclusion, measuring the expression of these genes has potential for implementation in a diagnostic tool for the early detection of MAP infection.
Mycobacterium avium subsp. paratuberculosis (MAP) causes Johne's disease (JD) in ruminants, characterized by a chronic granulomatous intestinal infection . Infections with MAP are prevalent worldwide in ruminants, but are also suggested to play a role in Crohn's disease in humans . Both are driving the increasing interest in MAP and JD. The primary route of transmission is fecal-oral, with infected cattle typically remaining subclinical for an extended interval . In the end stages of JD, MAP-infected cows have chronic diarrhoea and emaciation, leading to early culling. Infection with MAP is estimated to cost the Canadian dairy industry $15 to 20'M annually .
Although MAP was discovered as the causative organism of JD early in the 19th century, there are still numerous gaps in our knowledge with regard to JD pathogenesis and diagnostic methods that allow early detection of MAP are lacking. Existing diagnostic tools and strategies are inadequate, due to lack of sensitivity to detect MAP infection during early subclinical stages - when diagnostic test like ELISA and fecal culture fail to detect infection or have poor sensitivity. Hence, novel biomarkers are considered as alternatives for early identification of the infection. Translation genomics and protein arrays have been applied to identify biomarkers for early detection of JD -. Potential biomarkers for diagnosis include genes involved in host stress and immune response to the disease, suggesting that understanding the pathophysiology of MAP infection is crucial in identifying specific biomarkers. Moreover, a practical biomarker needs to be specific, robust and easily measurable without invasive procedures . In a previous study, several putative biomarkers of early (3 months after inoculation) MAP infection were found with particular roles in the immune response . Furthermore, the importance of dose of infection on the discovery of biomarkers was also apparent in that study. The objectives of the current study were to identify potential transcriptional biomarkers for MAP infection at 6 and 9 months stage of inoculation and to confirm the influence of infective dose on these transcripts. Furthermore, we aimed to analyse expression levels of select gene targets over a 15 months period after inoculation to validate their potential use as early diagnostic biomarkers.
Materials and methods
Selection, nutrition, health and husbandry of the animals and the design of the MAP infection experiment have been reported ,. In short, Holstein-Friesian bull calves were procured from Alberta (Canada) dairy farms with zero MAP fecal culture prevalence and a low (< 5%) MAP sero-prevalence. Five calves were orally inoculated with a high dose (HD)  of MAP, 5 with a low dose (LD), and another 5 were kept as non-inoculated controls; the bacteria were a virulent cattle type MAP strain isolated from a clinical Alberta JD case (Cow 69). This isolate has an identical BamHI, PvuII and PstI IS900 - RFLP profile as the reference strain K10 recommended for use in infection trials . Calves were inoculated at 2 weeks of age on 2 consecutive days, with either 5 × 109 CFU (HD) or 5 × 107 CFU (LD). The animal care protocols were approved by the Animal Care Committee of the University of Calgary (M09083).
MAP exposure assays
Both control and MAP-inoculated calves were tested every month for exposure to MAP. Humoral response was determined with a commercial antibody ELISA (Pourquier ELISA™; Institut Pourquier, Montpellier, France) as previously described in  whereas the cell-mediated immune response used an IFN-γ release assay. For the latter, whole blood was incubated in a 24-well plate overnight with Johnin PPD and post-incubation supernatants were used for BOVIGAM® IFN-γ ELISA (Prionics, Schlieren-Zurich, Switzerland). Liquid MAP culture (TREK para-JEM®; TREK Diagnostic Systems, Cleveland, OH, USA) was performed as described previously in  on fecal samples collected weekly, starting in the first month after inoculation and thereafter, monthly until necropsy. All calves were euthanized at 17 months of age. Macroscopic and histological lesions were assessed and bacterial culture was done on numerous tissues .
Sample collection and preparation
Every 3 months, whole blood was collected from HD, LD and control calves in PAXgene® Blood RNA system tubes (PreAnalytix GmbH, Hombrechtikon, Switzerland). Total RNA was extracted from whole blood using the PAXgene® blood miRNA kit (PreAnalytix GmbH) as per kit protocol. The quality of total RNA was determined using an RNA integrity number (Agilent RNA 6000 NanoChip on 2100 Bioanalyzer, Agilent Technologies, Santa Clara, CA, USA). Moreover, 5-10 μg of extracted total RNA was processed using RNeasy Plus Micro kit (Qiagen, Mississauga, ON, Canada) to remove any genomic DNA carryover. The RNA was quantified using NanoDrop ND-1000 (NanoDrop Technologies, Wilmington, DE, USA). Although total RNA was extracted from whole blood collected from trial calves at 3, 6, 9, 12 and 15 months, only total RNA from 6 and 9 months after MAP exposure was used for gene expression analysis by microarray.
Biotin-labelled antisense RNA (aRNA) was prepared from 100 ng of genomic DNA purified from total RNA using a GeneChip® 3'-IVT express kit (Affymetrix, Santa Clara, CA, USA). Furthermore, 12 μg of biotin-labelled aRNA was fragmented (35-200 nt) at 95 °C for 35 min for hybridization. Fragmented biotin-labelled aRNA were hybridized to an Affymetrix® GeneChip® Bovine genome Array at 45 °C for 16-18 h. Arrays were stained with streptavidin-phycoerythrin and washed using Affymetrix® GeneChip Fluidics 450 following manufacturer's protocol and scanned with an Affymetrix® GeneChip Scanner 3000 7G System at the Southern Alberta Microarray Facility (Calgary, AB, Canada).
Gene expression statistical analysis
Gene expression profiles were generated as CEL files using GeneChip® Command Console® Software (AGCC), with GeneSpring™ Gx (Agilent Technologies, Santa Clara, CA, USA) used for statistical analyses. Raw data were normalized and summarized using a robust multichip averaging (RMA) algorithm and Probe Logarithmic Intensity ERror (PLIER) algorithm, as reported . Both RMA and PLIER algorithms employ quantile normalization; the former uses perfect match (PM) to correct the background fluorescence intensity, whereas the latter uses both PM and mismatch probes from the error correction. Gene transcripts were filtered based on probe expression levels (default settings for low and high were 20 and 80%, respectively), with coefficient variation ≤ 50% (part of quality control). In addition, quality control passed probe sets were analysed using ANOVA to check for differentially expressed probe sets among HD, LD and control calves. The Benjamini-Hochberg FDR correction method was used to obtain reliable gene expression data. Finally, 1.5-fold change analysis was done to determine fold change of differentially expressed probes. Hierarchical clustering and Principal Component Analysis (PCA) were performed on the entities designated differentially expressed (by ANOVA) to determine the significant differential gene expression profile between analysed groups.
An additional analysis was performed with ANOVA (p ≤ 0.10) on both RMA and PLIER to obtain entities for systems biology analysis. Ingenuity® Systems Pathway Analysis (IPA; Ingenuity Systems, Redwood City, CA, USA) was used to identify biological relevance, annotate and predict molecular and cellular functions of each gene and identify their involvement in biological processes. Differentially expressed probesets (based on ANOVA) were submitted as Affymetrix® GeneChip® Bovine Genome Array gene set and annotated using bovine genome array databases available on IPA.
Significance of relationship to a functional category or pathway was assigned by calculating the ratio between total molecules in that function by the number of molecules from the dataset. Significance was also estimated by a right-tailed Fisher's exact test (relationships with p ≤ 0.05 were considered significant). In addition, Z-scores were calculated for upstream and downstream analyses and tested for significance to place differentially expressed transcripts in the context of functional categories and pathways.
Reverse transcription for qPCR
Total RNA, extracted from the blood of the calves at 3, 6, 9, 12 and 15 months after infection, was used to prepare cDNA for qPCR analysis. Reverse Transcription (RT) reaction was carried out using The Quantitect® Reverse transcription kit (Qiagen, Mississauga, ON, Canada). For this, 1 μg of genomic DNA purified total RNA was used to prepare cDNA, which involved a genomic DNA elimination reaction and an RT reaction at 42 °C for 10 and 30 min, respectively. The RT mix included oligo-dT and random primers and omniscript and sensiscript reverse transcriptase enzymes to synthesize cDNA. For the post RT-reaction, all cDNA samples was diluted to 100 ng/μL and stored at -20 °C to be used for qPCR validation.
Based on fold change and relevance to the study, 6 genes from 6 months after infection gene expression data and 5 from 9 month gene expression data were chosen for qPCR validation. Intron-spanning primers for these genes were designed using the Primer3 Version 0.4.0 online tool. Primers were synthesized at University of Calgary DNA synthesis lab (Calgary, AB, Canada); primer sequences and information are shown (Tables 1 and 2). Based on previous publications , and geNorm analysis done using qBasePLUS (Biogazalle, Zwijnaarde, Belgium), GAPDH was used as a reference gene for normalization. For these studies, 100 ng of cDNA was used for all samples and experiments were done using Qiagen quantitect SYBR green reagents (Qiagen, Mississauga, ON, Canada) on CFX96™ Real-Time PCR detection system. The MIQE guidelines were adopted for qPCR confirmation assays, as recommended . Prior to gene expression analysis, repeatability and efficiency for each primer was analysed using Bio-Rad CFX Manager™ 2.0 (Bio-Rad, Mississauga, ON, Canada) as part of the primer standardization procedure. Pooled RNA samples were used as non-reverse transcription (NRT) control to measure genomic DNA carry over. Real-time PCR amplification in this study involved initial enzyme activation at 95 °C for 15 min and 45 cycles of 95 °C for 10 s, 15's at 60 °C, and 15's at 72 °C extension step. Melt curve analysis was done at the end of amplification to ensure product specificity. Cycles of threshold values obtained by gene expression qPCR assays were exported in Excel format using Bio-Rad CFX Manager™ 2 for statistical analysis. Calculations for gene expression analysis were done in Excel, using the 2-δδCT method; log2 fold changes obtained were then analysed by ANOVA (p ≤ 0.05) with Tukey post hoc tests done using GraphPad Prism version 5.0c (La Jolla, CA, USA) for Mac OS X.
Additional qPCR reactions were done on selected genes (BOLA, CD46 and BNBD9-like) on samples obtained at 3, 6, 9, 12 and 15 months to characterize expression profiles of these genes over time. Two-way ANOVA (p ≤ 0.05) was performed on the longitudinal gene expression analysis to determine statistical significance for the expression profiles of these genes.
MAP exposure test results
All MAP-inoculated calves (HD and LD) became and remained positive for Johnin PPD-specific IFN-γ response within 3 months after inoculation, whereas all control calves remained negative for the duration of the study. All MAP-inoculated calves, except one LD calf, shed MAP at least once during the 17-month follow-up and had MAP-positive intestinal tissues at necropsy. Furthermore, MAP was detected in fecal samples from one control calf by PCR, but no viable MAP was not cultured.
Control calves remained ELISA-negative, whereas 4 out of 5 (80%) HD MAP-inoculated cattle and 2 out of 5 (40%) LD MAP-inoculated cattle became ELISA-positive.
Gene expression 6 months after inoculation
Gene expression analysis, using RMA and PLIER algorithm summarized entities, identified 101 transcripts that were differentially expressed, based on ANOVA using RMA and 186 transcripts using PLIER algorithm. Likewise, 1.5 fold-change analyses revealed 37 transcripts differentially expressed in RMA and 21 transcripts using the PLIER algorithm. In addition, PCA done on transcripts obtained by ANOVA resulted in clear separation of the samples into HD, LD and control groups (Figure 1A and B). In the overall gene expression profile, both HD and LD calves had 60% of their differentially expressed genes upregulated compared to the control. Furthermore, HD and LD calves had opposite expression profiles for 50% of their differentially expressed genes. Using hierarchical clustering on MAP exposure (HD, LD and control), control calves clustered in a separate branch from HD and LD with a centroid Euclidean distance of 0.32. HD and LD calves clustered as sub clusters under a same branch with a centroid Euclidean distance of 0.17 from each other. Additional ANOVA analysis with p ≤ 0.10 on both RMA and PLIER summarized entities resulted in 1702 entities for systems biology analysis. A complete list of the differentially expressed genes is included (Additional file 1).
Arachidonate 15-lipoxygenase (ALOX15), Arachidonate 5-lipoxygenase activating protein (ALOX5AP), neutrophil beta-defensin-9 like peptide (BNBD9-Like), S100 calcium binding protein A9 (s100A9) and G protein coupled receptor 77 (GPR77) or C5a anaphylatoxin chemotactic receptor C5a2 were downregulated in both HD and LD calves (Table 3).
Systems biology analyses for genes differentially expressed at 6 months after inoculation
Of the 1702 entities submitted to IPA, 1066 were used for systems biology analyses. Downstream function analysis revealed activation of lymphocyte movement and its migration, migration of mononuclear leukocytes and intracellular infection of the cells. In addition, IPA predicted inhibition of phagocytosis by antigen presenting cells, phagocytosis by macrophages, migration of granulocytes, immune response by macrophages, necrosis, apoptosis and downregulation of genes that would inhibit growth. All these predictions had a significant Z-score (either ≥ 2 signifying activation or ≥ 2 indicating downregulation). Pathway analysis indicated that the differentially expressed genes had roles in several biological pathways; the top 5 canonical pathways included the superpathway of Inositol phosphate compounds, D-myo-inositol (1,4,5,6)-Tetrakisphosphate biosynthesis, D-myo-inositol (3,4,5,6)-Tetrakisphosphate biosynthesis, D-myo-inositol-5-phosphate metabolism and 3-phosphoinositide degradation. JAK-STAT signalling pathway, a key regulator of immune response, was among the top canonical pathways (14th rank). The complete list of the canonical pathways influenced by our differentially expressed genes is shown (Additional file 2).
Upstream analysis predicted activation of CD40 ligand (CD40L), CD24, matrix metallopeptidase 14 (membrane-inserted) (MMP14) and TNF family proteins. Transforming growth factor beta 1 (TGF-beta1), several micro RNA's (miRNA) were predicted to be inhibited based on expression levels of molecules differentially expressed among the 3 groups. A list of all the upstream regulators and their direction of activation along with Z-score is shown (Additional file 3).
Gene expression 9 months after inoculation
Analysis of gene expression using RMA identified 31 and 12 transcripts (based on ANOVA and 1.5 fold-change analyses, respectively). Furthermore, ANOVA and 1.5 fold-change analyses using PLIER algorithm detected 49 and 3 transcripts, respectively, that were differentially expressed. Furthermore, PCA analysis on entities differentially expressed by ANOVA (p ≤ 0.05) resulted in clear separation of HD, LD and control groups (Figure 1C and D). In the overall gene expression profile, HD calves had 55% of their transcripts upregulated compared to control calves, whereas the LD group had 86% of transcripts upregulated compared to controls. Hierarchical clustering on conditions (HD, LD and control groups) of differentially expressed entities had HD calves clustered in a separate branch from control and LD groups (centroid Euclidean distance of 0.25). Control and LD calves were subclusters under the same branch (centroid Euclidean distance of 0.20). Additional ANOVA analysis with p ≤ 0.10 performed on both RMA and PLIER algorithm summarized entities resulted in 408 differentially expressed transcripts which were used for systems biology analyses. A complete list of differentially expressed transcripts with their fold-changes is shown (Additional file 4).
Nine months after inoculation, bovine leukocyte antigen (BOLA), interleukin 4 receptor (IL4R), chemokine (C-C motif) receptor 7 and testis expressed 261 (TEX261) genes were downregulated in HD and LD calves. Immunoglobulin super family member 6 (IGSF6), ribosomal protein L15 (RPL15) were upregulated in both LD and HD calves.
Systems biology analyses for genes differentially expressed at 9 months after inoculation
Of the 408 entities submitted to IPA, 275 entities were used for systems biology analyses. Downstream function analysis predicted differentially expressed transcripts to activate autophagy, activation of antigen presenting cells, cytotoxicity of T-lymphocytes, and organisation of the cytoplasm for cell replication. Downstream analysis also predicted inhibition of apoptosis, generation of T-lymphocytes and cell spreading. However, these functions lacked significant Z-scores.
The top 5 canonical pathways included acute myeloid leukaemia signalling, methionine degradation I (to homocysteine), cysteine biosynthesis III (Mammalia), ERK5 signalling and Neurotrophin/TRK signalling. The IL-4 signalling pathway, involved in Th2 response, was the 9th canonical pathway. A complete list of all canonical pathways in which differentially expressed genes were predicted to be involved is shown (Additional file 5). Upstream analysis predicted inhibition of several miRNAs, hypoxia inducible factor 1, alpha subunit (HIF1A), chemokine (C-X-C motif) ligand 12, forkhead box O3 (FOXO3) and myelocytomatosis viral oncogene homolog (Avian) (MYC) and predicted few miRNAs to be activated. A list including complete upstream analysis results is shown (Additional file 6).
Overlap of affected pathways at 3, 6 and 9 months after inoculation
Differentially expressed genes were found in a total of 476 canonical pathways, at either 3 , 6 and 9 months after inoculation, with 55%, 73% and 56% of these pathways affected at the respective time points. A total of 89 (19%) of these pathways were affected at every of the three time points, and 36%, 43% and 24% overlapped at respectively 3 and 6, 6 and 9 and 3 and 9 months after inoculation. At those respective times, only 24%, 17% and 12% of the pathways were unique to these time points.
Real-time qPCR confirmation microarray data
To validate microarray findings, 6 gene targets from the 6 months after inoculation gene expression data and 5 gene targets from 9 months after inoculation gene expression data were chosen to be verified using real-time qPCR. Genes for qPCR validation were chosen based on the relevance of the transcripts to JD and their fold change. Genes selected to confirm gene expression data at 6 months after inoculation were BOLA, BNBD9-Like, ALOX-15, ALOX5AP, GPR77 and s100A9. Similarly, genes selected for qPCR confirmation of gene expression data at 9 months after inoculation were BOLA, TEX261, CCR7, IL4R and IGSF6. Quantitative gene expression data obtained for these genes using qPCR was in agreement with microarray data. Real-time qPCR confirmation of the differentially expressed genes at 6 and 9 months after MAP inoculation are shown (Figures 2 and 3, respectively). A comparison of microarray expression data and qPCR expression data for the confirmed genes obtained 6 and 9 months after MAP inoculation are presented (Tables 3 and 4).
Longitudinal qPCR expression analysis of CD46, BOLA and BNBD9-Like genes
Based on the gene expression studies by microarray analysis at 3, 6 and 9 months after MAP inoculation, CD46, BOLA and BNBD9-Like were differentially expressed at all 3 measured time points. Longitudinal expression profiling of these genes at 3, 6, 9, 12 and 15 months after MAP inoculation by qPCR, demonstrated that CD46 was continuously upregulated in both HD and LD calves from 3 to 15 months after MAP inoculation onwards, except for 3 months after MAP inoculation for LD calves. Both BOLA and BNBD9-Like genes were downregulated in MAP-inoculated calves (HD and LD) at time points 3 to 15 months after inoculation. Two-way ANOVA done on the longitudinal expression values for all 3 genes revealed significant differences (p ≤ 0.01). The relative expression of CD46, BOLA and BNBD9-Like genes across time are shown in Figure 4.
Inoculation with MAP (HD or LD) resulted in a chronic infection in all calves, confirmed by fecal shedding, antibody ELISA  and positive culture of intestinal tissues . However, one LD calf merely had a repeated positive interferon-gamma test as an indication of exposure to MAP.
Early diagnosis of MAP, similar to diagnosis of Mycobacterium tuberculosis, is challenging and largely dependent on detection of cellular immune responses, due to low sensitivity for other tests , at this stage of infection. Regardless, detection of cellular immune responses is undermined by test complexity and specificity, especially during the early stages of the disease ,,. Hence, new biomarkers for early diagnosis of Mycobacterial infections are needed. Unlike previous studies, the current study focused on identifying biomarkers for JD from differentially expressed transcripts in whole blood of HD and LD MAP-inoculated calves at 6 and 9 months after inoculation with MAP.
In this study, key inflammation and immune-related genes (ALOX15, ALOX5AP, BNBD9-Like, GPR77, and S100A9) were downregulated 6 months after inoculation. ALOX15 belongs to lipoxygenases (LOXs) involved in metabolism and production of fatty acid hydroperoxidases . ALOX5AP along with ALOX5 has a role in leukotriene biosynthesis, which is implicated in various inflammatory responses . Enzyme 12/15-LOX or 15-LOX1 is activated by IL-4 in Th2 type response in monocytes and macrophages ; it has a role in resolution of inflammation , and synthesis of lipoxins, which also reduces inflammation by inhibiting chemotaxis, adhesion and superoxide generation . Down-regulation of ALOX15 and ALOX5AP in the early stages of MAP infection correspond to the expected Th1 response-driven IFN-γ production that inhibits expression of LOXs genes. Proteins BNBD9-Like and GPR77 act as chemoattractants to immature dendritic cells ,, which function as professional antigen-presenting cells. S100A9 is a pro-inflammatory damage-associated molecular pattern (DAMPs) molecule secreted by phagocytes, granulocytes, monocytes and early differentiation stage macrophages. DAMPs are recognized by Toll-like receptors (TLRs), leading to their activation, which results in inflammation . Down-regulation of these inflammatory and immune-related genes was consistent with a reduced capacity for an appropriate inflammatory response in MAP-infected animals.
Canonical pathway analysis revealed alterations in phosphoinositol biosynthesis and metabolism, and the JAK-STAT signalling pathway. Phosphatidylinositol signalling pathways are of interest in MAP infection, since these lipids act as chemoattractants and mediate immune cell migration . Furthermore, Mycobacterium spp. is known to control host lipid metabolism to establish an infection . Many genes involved in phosphatidylinositol signalling pathways were downregulated in this study. Conversely, the JAK-STAT pathway is the main signalling cascade involved in activation and regulation of host immune response following activation by cytokines and growth factors . In this study, there were activation of STAT proteins and downregulation of other JAK-STAT association pathways such as Ras and MEK 1/2. Moreover, the negative regulator of the JAK-STAT pathway suppressors of cytokine signalling (SOCS) was upregulated and protein inhibitor of activated STATs (PIAS) was downregulated. The phosphoinositide 3-kinase (PI3K) pathway was also activated, leading to m-TOR mediated activation of STAT3. Therefore, we concluded that the overall JAK-STAT pathway was activated in MAP-inoculated calves 6 months after inoculation.
Lymphocyte and leukocyte migration were activated in the downstream pathway analysis. There was upregulation of chemokine (C-X-C motif) ligand 10 (CXCL10), SH2 domain containing 1A (SH2D1A), TXK tyrosine kinase (TXK); they enable lymphocyte activation and migration -. The latter was consistent with the expected movement of leukocyte and MAP antigen-activated lymphocytes to the site of inflammation or infection.
At 6 months after inoculation, key immune functions such as macrophage response and phagocytosis migration of granulocytes and phagocytosis by other antigen presenting cells were downregulated. This was in agreement with MAP-induced inhibition of macrophage maturation described in other studies ,. Necrosis and apoptosis of cells were concurrently inhibited, as reported at 3 months after inoculation  and also in other studies - as a virulence mechanism of MAP to subvert immune detection and immune escape.
Interestingly, genes that would facilitate failure of growth were downregulated, as reported 3 month after MAP inoculation . This suggests a response counteracting malabsorption, due to early host immune response and inflammation which manifested as histological and gross pathological lesions in MAP-infected tissues , as a homeostatic mechanism to support growth of the animal.
CD40L was activated based on upstream analysis done using differentially expressed transcripts (6 months after inoculation) between HD, LD and control. CD40L is expressed on CD4+ T cells; it interacts with CD40 on the B cells, and activates them to produce a B cell response, immunoglobulin class switch and activation of macrophages to produce IFN-γ . It is noteworthy that 3 of the 5 HD animals in this study had seroconverted 6 months after inoculation. Based on upstream analysis, TGF-beta1 and miRNAs like let-7a-5p and miR-16-5p were suggested to be inhibited at 6 months after inoculation. TGF-beta1 is an anti-inflammatory, pro-fibrotic and macrophage deactivating cytokine, implicated in the pathogenesis of many intracellular pathogens ,. Mycobacterium avium complex (MAC) pathogens induce many cytokines, especially TGF-beta1 ,. TGF-beta1 has an important role in granulomatous infection; its levels are related to intracellular replication of MAC . Future studies analysing TGF-beta1 in subclinical and clinical animals might help to validate it as biomarker for MAP diagnosis and to understand its role in immunopathogenesis.
Micro RNAs regulate expression of many genes by post-translation modification of mRNAs and are associated with cancer and infectious disease -. In Crohn's disease (CD) patients, several miRNAs (e.g.miR-16 and miR-23b) were associated with active and chronically active CD patients ,. In the current study, miR-2277-3p was activated, whereas several other miRNAs were inhibited. The exact role of these miRNAs is not known, but they have potential as markers for JD.
At 9 months after inoculation, gene expression of genes was generally repressed; based on ANOVA, there were only 80 transcripts differentially expressed among HD, LD and control calves (p < 0.05). Similarly, in previous translational studies , on gene expression in macrophages and PBMC's infected with MAP, the number of differentially expressed genes decreased with time after inoculation.
BOLA/MHC-1 was downregulated 9 months after inoculation, in agreement with a previous study . BOLA is a key antigen-presenting protein that carries intracellular proteins to the cell membrane and presents it to cytotoxic T cells for detection and eventual control or killing of intracellular pathogens. In addition, IL4R and CCR7 were also downregulated. IL4R and CCR7 belong to G-protein coupled receptors present on B and T lymphocytes and leukocytes, enabling these cells to move to sites of inflammation and secondary lymphoid organs ,.
Autophagy, antigen presenting cells, cytotoxicity of T-lymphocytes and organisation of the cytoplasm for cell replication were all activated 9 months after inoculation. The role of autophagy in Mycobacterial infections has been reported ,; furthermore, we reported that autophagy was already activated 3 months after inoculation . Perhaps autophagy acts as a compensatory mechanism to present intracellular MAP antigens in the context of reduced antigen presentation by MHC-1/BOLA.
Canonical pathway analysis revealed differentially expressed genes in MAP-infected calves to be involved in IL-4 signalling pathway. IL-4 is a multifunctional cytokine produced by CD4+ Th2 cells, basophils and mast cells. In this study, there was upregulation of son of sevenless (SOS), 1-phosphotidylionositol 3-phosphate, P70 S6 kinase (P70S6K) and nuclear factor of activated T cells (NFAT) belonging to alternate IL-4 signalling pathway. However, IL-4R in the classical IL-4 signalling pathway was downregulated. IL-4 is important for granuloma formation in Mycobacterial infection and defense against Mycobacterial infection . At 9 months after inoculation, not many genes were differentially expressed; consequently, upstream analysis did not lead to many predictions. Activation of miR-22-3p, miR-4283, miR-320b, miR-3175 and few other miRNAs are listed in Additional file 6, according to upstream analysis at 9 months after inoculation.
BOLA and BNBD9-Like genes were downregulated and CD46 gene was upregulated in both LD and HD calves for the duration of the trial (3 to 15 months after inoculation). Identification of genes that are consistently differentially expressed genes warrants a larger cohort study to validate these genes as potential biomarkers for MAP infection. In this study, BNBD9-LIKE and BOLA genes were differentially expressed 6 and 9 months after inoculation. The ubiquitously expressed transmembrane glycoprotein CD46 was upregulated 3 months after inoculation . It prevents unwanted complement-mediated killing of cells by acting as a cofactor in factor-I mediated degradation of C3b and C4b opsonins . Although CD46 has a regulatory role in Th1 responses ,, it also interacts with STE20/SPS1-related proline/alanine-rich kinase (SPAK) kinase and E-cadherin, which have important role sin maintaining epithelial barrier function in intestinal epithelial cells .
A large proportion of pathways, including immune and inflammation pathways, were consistently affected at 3, 6 and 9 months after inoculation, with relatively few pathways being affected at only a single time point. However, no obvious evolution was noticed in the affected pathways over the course of this follow-up.
In conclusion, many genes were differentially expressed between MAP-infected animals and control animals at 6 and 9 months after inoculation. BOLA, BNBD9-Like and CD46 were differentially expressed in both LD and HD calves as early as 3 months after inoculation until at least 15 months after inoculation. Dose of infection had an influence on the fold change of differentially expressed genes. Migration and trafficking of leukocytes and lymphocytes were activated, but there was more downregulation of immune response by inhibition of phagocytosis, downregulation of antigen presentation, macrophage phagocytosis and response. Furthermore, MHC-1/BOLA was downregulated and autophagy was activated, potentially compensating for any loss in antigen presentation by PRRs on the account of downregulation of MHC-1/BOLA. Finally, we propose that CD46, BNBD9-Like and BOLA are further considered as potential biomarkers for diagnosis of JD.
Manning EJ, Collins MT:Mycobacterium avium subsp. paratuberculosis: pathogen, pathogenesis and diagnosis. Rev Sci Tech. 2001, 20: 133-150.
Barkema HW, De Buck JM, Ghosh S, Kaplan GG, Rioux KP: Crohn's disease in humans and Johne's disease in cattle - Linked diseases?. Zoonatic Pathogens In The Food Chain. Edited by: Krause D, Hendrick S. 2010, CABI, Wallingford, Oxfordshire; Cambridge, MA, xii, 242 p
Whittington RJ, Sergeant ES: Progress towards understanding the spread, detection and control of Mycobacterium avium subsp paratuberculosis in animal populations. Aust Vet J. 2001, 79: 267-278. 10.1111/j.1751-0813.2001.tb11980.x.
Tiwari A, VanLeeuwen JA, McKenna SL, Keefe GP, Barkema HW: Johne's disease in Canada Part I: clinical symptoms, pathophysiology, diagnosis, and prevalence in dairy herds. Can Vet J. 2006, 47: 874-882.
Dargatz DA, Byrum BA, Barber LK, Sweeney RW, Whitlock RH, Shulaw WP, Jacobson RH, Stabel JR: Evaluation of a commercial ELISA for diagnosis of paratuberculosis in cattle. J Am Vet Med Assoc. 2001, 218: 1163-1166. 10.2460/javma.2001.218.1163.
Alinovi CA, Ward MP, Lin TL, Moore GE, Wu CC: Real-time PCR, compared to liquid and solid culture media and ELISA, for the detection of Mycobacterium avium ssp. paratuberculosis. Vet Microbiol. 2009, 136: 177-179. 10.1016/j.vetmic.2008.10.012.
Wells SJ, Collins MT, Faaberg KS, Wees C, Tavornpanich S, Petrini KR, Collins JE, Cernicchiaro N, Whitlock RH: Evaluation of a rapid fecal PCR test for detection of Mycobacterium avium subsp. paratuberculosis in dairy cattle. Clin Vaccine Immunol. 2006, 13: 1125-1130. 10.1128/CVI.00236-06.
Stabel JR: Production of gamma-interferon by peripheral blood mononuclear cells: an important diagnostic tool for detection of subclinical paratuberculosis. J Vet Diagn Invest. 1996, 8: 345-350. 10.1177/104063879600800311.
Mattos AM, Almeida CS, Franken KL, Alves CC, Abramo C, de Souza MA, L’Hotellier M, Alves MJ, Ferreira AP, Oliveira SC, Ottenhoff TH, Teixeira HC: Increased IgG1, IFN-gamma, TNF-alpha and IL-6 responses to Mycobacterium tuberculosis antigens in patients with tuberculosis are lower after chemotherapy. Int Immunol. 2010, 22: 775-782. 10.1093/intimm/dxq429.
Seth M, Lamont EA, Janagama HK, Widdel A, Vulchanova L, Stabel JR, Waters WR, Palmer MV, Sreevatsan S: Biomarker discovery in subclinical mycobacterial infections of cattle.PLoS One 2009, 4:e5478.,
Ruhwald M, Aabye MG, Ravn P: IP-10 release assays in the diagnosis of tuberculosis infection: current status and future directions. Expert Rev Mol Diagn. 2012, 12: 175-187. 10.1586/erm.11.97.
David J, Barkema HW, Mortier RAR, Ghosh S, Le Luo G, De Buck J: Gene expression profiling and putative biomarkers of calves 3 months after infection with Mycobacterium avium subspecies paratuberculosis. Vet Immunol Immunopathol. 2014, 160: 107-117. 10.1016/j.vetimm.2014.04.006.
Mortier RA, Barkema HW, Bystrom JM, Illanes O, Orsel K, Wolf R, Atkins G, De Buck J: Evaluation of age-dependent susceptibility in calves infected with two doses ofMycobacterium aviumsubspeciesparatuberculosisusing pathology and tissue culture.Vet Res 2013, 44:94.,
Kuehnel MP, Goethe R, Habermann A, Mueller E, Rohde M, Griffiths G, Valentin-Weigand P: Characterization of the intracellular survival of Mycobacterium avium ssp. paratuberculosis: phagosomal pH and fusogenicity in J774 macrophages compared with other mycobacteria. Cell Microbiol. 2001, 3: 551-566. 10.1046/j.1462-5822.2001.00139.x.
Hines ME, Stabel JR, Sweeney RW, Griffin F, Talaat AM, Bakker D, Benedictus G, Davis WC, de Lisle GW, Gardner IA, Juste RA, Kapur V, Koets A, McNair J, Pruitt G, Whitlock RH: Experimental challenge models for Johne's disease: a review and proposed international guidelines. Vet Microbiol. 2007, 122: 197-222. 10.1016/j.vetmic.2007.03.009.
Mortier RAR, Barkema HW, Negron ME, Orsel K, Wolf R, De Buck J: Antibody responses early after experimental infection with Mycobacterium avium subspecies paratuberculosis in dairy calves. Diary Sci. 2014, 97: 5558-5565. 10.3168/jds.2014-8139.
Mortier RAR, Barkema HW, Orsel K, Wolf R, De Buck J: Shedding patterns of dairy calves experimentally infected with Mycobacterium avium subspecies paratuberculosis.Vet Res 2014, 45:71.,
Robinson TL, Sutherland IA, Sutherland J: Validation of candidate bovine reference genes for use with real-time PCR. Vet Immunol Immunopathol. 2007, 115: 160-165. 10.1016/j.vetimm.2006.09.012.
Jin W, Olson EN, Moore SS, Basarab JA, Basu U, Guan LL: Transcriptome analysis of subcutaneous adipose tissues in beef cattle using 3' digital gene expression-tag profiling. J Anim Sci. 2012, 90: 171-183. 10.2527/jas.2011-4229.
Bustin SA, Benes V, Garson JA, Hellemans J, Huggett J, Kubista M, Mueller R, Nolan T, Pfaffl MW, Shipley GL, Vandesompele J, Wittwer CT: The MIQE guidelines: minimum information for publication of quantitative real-time PCR experiments. Clin Chem. 2009, 55: 611-622. 10.1373/clinchem.2008.112797.
Paolicchii FA, Zumarraga MJ, Gioffre A, Zamorano P, Morsella C, Verna A, Cataldi A, Alito A, Romano M: Application of different methods for the diagnosis of paratuberculosis in a dairy cattle herd in Argentina. J Vet Med B Infect Dis Vet Public Health. 2003, 50: 20-26. 10.1046/j.1439-0450.2003.00606.x.
Detjen AK, Keil T, Roll S, Hauer B, Mauch H, Wahn U, Magdorf K: Interferon-gamma release assays improve the diagnosis of tuberculosis and nontuberculous mycobacterial disease in children in a country with a low incidence of tuberculosis. Clin Infect Dis. 2007, 45: 322-328. 10.1086/519266.
Gwozdz JM, Thompson KG, Murray A, Reichel MP, Manktelow BW, West DM: Comparison of three serological tests and an interferon-gamma assay for the diagnosis of paratuberculosis in experimentally infected sheep. Aust Vet J. 2000, 78: 779-783. 10.1111/j.1751-0813.2000.tb10452.x.
Wallis RS, Pai M, Menzies D, Doherty TM, Walzl G, Perkins MD, Zumla A: Biomarkers and diagnostics for tuberculosis: progress, needs, and translation into practice. Lancet. 2010, 375: 1920-1937. 10.1016/S0140-6736(10)60359-5.
Conrad DJ: The arachidonate 12/15 lipoxygenases. A review of tissue expression and biologic function. Clin Rev Allergy Immunol. 1999, 17: 71-89. 10.1007/BF02737598.
Strom JO, Strid T, Hammarstrom S: Disruption of the alox5ap gene ameliorates focal ischemic stroke: possible consequence of impaired leukotriene biosynthesis.BMC Neurosci 2012, 13:146.,
Conrad DJ, Kuhn H, Mulkins M, Highland E, Sigal E: Specific inflammatory cytokines regulate the expression of human monocyte 15-lipoxygenase. Proc Natl Acad Sci U S A. 1992, 89: 217-221. 10.1073/pnas.89.1.217.
Serhan CN, Jain A, Marleau S, Clish C, Kantarci A, Behbehani B, Colgan SP, Stahl GL, Merched A, Petasis NA, Chan L, Van Dyke TE: Reduced inflammation and tissue damage in transgenic rabbits overexpressing 15-lipoxygenase and endogenous anti-inflammatory lipid mediators. J Immunol. 2003, 171: 6856-6865. 10.4049/jimmunol.171.12.6856.
Gronert K, Maheshwari N, Khan N, Hassan IR, Dunn M, Laniado Schwartzman M: A role for the mouse 12/15-lipoxygenase pathway in promoting epithelial wound healing and host defense. J Biol Chem. 2005, 280: 15267-15278. 10.1074/jbc.M410638200.
Serhan CN: Lipoxins and aspirin-triggered 15-epi-lipoxins are the first lipid mediators of endogenous anti-inflammation and resolution. Prostaglandins Leukot Essent Fatty Acids. 2005, 73: 141-162. 10.1016/j.plefa.2005.05.002.
Mackenzie-Dyck S, Attah-Poku S, Juillard V, Babiuk LA, van Drunen Littel-van den Hurk S: The synthetic peptides bovine enteric beta-defensin (EBD), bovine neutrophil beta-defensin (BNBD) 9 and BNBD 3 are chemotactic for immature bovine dendritic cells. Vet Immunol Immunopathol. 2011, 143: 87-107. 10.1016/j.vetimm.2011.06.028.
Ohno M, Hirata T, Enomoto M, Araki T, Ishimaru H, Takahashi TA: A putative chemoattractant receptor, C5L2, is expressed in granulocyte and immature dendritic cells, but not in mature dendritic cells. Mol Immunol. 2000, 37: 407-412. 10.1016/S0161-5890(00)00067-5.
Foell D, Wittkowski H, Vogl T, Roth J: S100 proteins expressed in phagocytes: a novel group of damage-associated molecular pattern molecules. J Leukoc Biol. 2007, 81: 28-37. 10.1189/jlb.0306170.
Pike MC, Bruck ME, Arndt C, Lee CS: Chemoattractants stimulate phosphatidylinositol-4-phosphate kinase in human polymorphonuclear leukocytes. J Biol Chem. 1990, 265: 1866-1873.
Kusner DJ, Hall CF, Schlesinger LS: Activation of phospholipase D is tightly coupled to the phagocytosis of Mycobacterium tuberculosis or opsonized zymosan by human macrophages. J Exp Med. 1996, 184: 585-595. 10.1084/jem.184.2.585.
Aaronson DS, Horvath CM: A road map for those who don't know JAK-STAT. Science. 2002, 296: 1653-1655. 10.1126/science.1071545.
Stanford MM, Issekutz TB: The relative activity of CXCR3 and CCR5 ligands in T lymphocyte migration: concordant and disparate activities in vitro and in vivo. J Leukoc Biol. 2003, 74: 791-799. 10.1189/jlb.1102547.
Koyasu S: The role of PI3K in immune cells. Nat Immunol. 2003, 4: 313-319. 10.1038/ni0403-313.
Suzuki N, Nara K, Suzuki T: Skewed Th1 responses caused by excessive expression of Txk, a member of the Tec family of tyrosine kinases, in patients with Behcet's disease. Clin Med Res. 2006, 4: 147-151. 10.3121/cmr.4.2.147.
Smythies LE, Sellers M, Clements RH, Mosteller-Barnum M, Meng G, Benjamin WH, Orenstein JM, Smith PD: Human intestinal macrophages display profound inflammatory anergy despite avid phagocytic and bacteriocidal activity. J Clin Invest. 2005, 115: 66-75. 10.1172/JCI200519229.
Kabara E, Coussens PM: Infection of primary bovine macrophages withMycobacterium aviumsubspeciesparatuberculosissuppresses host cell apoptosis.Front Microbiol 2012, 3:215.,
Behar SM, Divangahi M, Remold HG: Evasion of innate immunity by Mycobacterium tuberculosis: is death an exit strategy?. Nat Rev Microbiol. 2010, 8: 668-674.
Schaible UE, Winau F, Sieling PA, Fischer K, Collins HL, Hagens K, Modlin RL, Brinkmann V, Kaufmann SH: Apoptosis facilitates antigen presentation to T lymphocytes through MHC-I and CD1 in tuberculosis. Nat Med. 2003, 9: 1039-1046. 10.1038/nm906.
van Kooten C, Banchereau J: CD40-CD40 ligand. J Leukoc Biol. 2000, 67: 2-17.
Border WA, Ruoslahti E: Transforming growth factor-beta in disease: the dark side of tissue repair. J Clin Invest. 1992, 90: 1-7. 10.1172/JCI115821.
Barral A, Barral-Netto M, Yong EC, Brownell CE, Twardzik DR, Reed SG: Transforming growth factor beta as a virulence mechanism for Leishmania braziliensis. Proc Natl Acad Sci U S A. 1993, 90: 3442-3446. 10.1073/pnas.90.8.3442.
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: 1229-1237. 10.1093/infdis/170.5.1229.
Horsburgh CR, Mason UG, Farhi DC, Iseman MD: Disseminated infection with Mycobacterium avium-intracellulare. A report of 13 cases and a review of the literature. Medicine. 1985, 64: 36-48. 10.1097/00005792-198501000-00003.
Bermudez LE, Young LS, Martinelli J, Petrofsky M: Exposure to ethanol up-regulates the expression of Mycobacterium avium complex proteins associated with bacterial virulence. J Infect Dis. 1993, 168: 961-968. 10.1093/infdis/168.4.961.
Schultz NA, Dehlendorff C, Jensen BV, Bjerregaard JK, Nielsen KR, Bojesen SE, Calatayud D, Nielsen SE, Yilmaz M, Hollander NH, Andersen KK, Johansen JS: MicroRNA biomarkers in whole blood for detection of pancreatic cancer. JAMA. 2014, 311: 392-404. 10.1001/jama.2013.284664.
Hogfeldt T, Johnsson P, Grander D, Bahnassy AA, Porwit A, Eid S, Osterborg A, Zekri AR, Lundahl J, Khaled MH, Mellstedt H, Moshfegh A: Expression of microRNA-1234 related signal transducer and activator of transcription 3 in patients with diffuse large B-cell lymphoma of activated B-cell like type from high and low infectious disease areas. Leuk Lymphoma. 2013, 55: 1158-1165. 10.3109/10428194.2013.824077.
Chen X, Zhou JY, Zhou JY: MicroRNA-34a: role in cancer and cardiovascular disease. Curr Drug Targets. 2014, 15: 361-373. 10.2174/1389450115666140120102935.
Wu F, Zhang S, Dassopoulos T, Harris ML, Bayless TM, Meltzer SJ, Brant SR, Kwon JH: Identification of microRNAs associated with ileal and colonic Crohn's disease. Inflamm Bowel Dis. 2010, 16: 1729-1738. 10.1002/ibd.21267.
Wu F, Guo NJ, Tian H, Marohn M, Gearhart S, Bayless TM, Brant SR, Kwon JH: Peripheral blood microRNAs distinguish active ulcerative colitis and Crohn's disease. Inflamm Bowel Dis. 2011, 17: 241-250. 10.1002/ibd.21450.
Machugh DE, Taraktsoglou M, Killick KE, Nalpas NC, Browne JA, DE Park S, Hokamp K, Gormley E, Magee DA: Pan-genomic analysis of bovine monocyte-derived macrophage gene expression in response to in vitro infection with Mycobacterium avium subspecies paratuberculosis.Vet Res 2012, 43:25.,
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 Immunol. 2003, 71: 6487-6498. 10.1128/IAI.71.11.6487-6498.2003.
Purdie AC, Plain KM, Begg DJ, de Silva K, Whittington RJ: Expression of genes associated with the antigen presentation and processing pathway are consistently regulated in early Mycobacterium avium subsp. paratuberculosis infection. Comp Immunol Microbiol Infect Dis. 2012, 35: 151-162. 10.1016/j.cimid.2011.12.007.
Schaeuble K, Hauser MA, Singer E, Groettrup M, Legler DF: Cross-talk between TCR and CCR7 signaling sets a temporal threshold for enhanced T lymphocyte migration. J Immunol. 2011, 187: 5645-5652. 10.4049/jimmunol.1101850.
Johnson LA, Jackson DG: Control of dendritic cell trafficking in lymphatics by chemokines. Angiogenesis. 2014, 17: 335-345. 10.1007/s10456-013-9407-0.
Delgado MA, Deretic V: Toll-like receptors in control of immunological autophagy. Cell Death Differ. 2009, 16: 976-983. 10.1038/cdd.2009.40.
Sugawara I, Yamada H, Mizuno S, Iwakura Y: IL-4 is required for defense against mycobacterial infection. Microbiol Immunol. 2000, 44: 971-979. 10.1111/j.1348-0421.2000.tb02592.x.
Liszewski MK, Post TW, Atkinson JP: Membrane cofactor protein (MCP or CD46): newest member of the regulators of complement activation gene cluster. Ann Rev Immunol. 1991, 9: 431-455. 10.1146/annurev.iy.09.040191.002243.
Astier A, Trescol-Biemont MC, Azocar O, Lamouille B, Rabourdin-Combe C: Cutting edge: CD46, a new costimulatory molecule for T cells, that induces p120CBL and LAT phosphorylation. J Immunol. 2000, 164: 6091-6095. 10.4049/jimmunol.164.12.6091.
Kemper C, Chan AC, Green JM, Brett KA, Murphy KM, Atkinson JP: Activation of human CD4+ cells with CD3 and CD46 induces a T-regulatory cell 1 phenotype. Nature. 2003, 421: 388-392. 10.1038/nature01315.
Cardone J, Al-Shouli S, Kemper C: A novel role for CD46 in wound repair.Front Immunol 2011, 2:28.,
Alberta Innovates - Technology Futures, Alberta Innovates - Bio Solutions, the Alberta Livestock and Meat Agency, Alberta Milk, Dairy Farmers of Canada and the Natural Sciences and Engineering Research Council of Canada supported this work. We thank Dr Robert Wolf for sample collection, Xiuling Wang from the Health Science Centre Microarray facility for assistance with the microarray analysis, and Dr John Kastelic for editing the manuscript.
The authors declare that they have no competing interests.
JD carried out the studies and prepared the manuscript. LG contributed to development of methods, and acquisition and analysis of the data. HB contributed to study design and helped write the manuscript. JDB conceived of the study and participated in its design and coordination and helped to draft the manuscript. All authors read and approved the final manuscript.
Electronic supplementary material
About this article
Cite this article
David, J., Barkema, H.W., Guan, L.L. et al. Gene-expression profiling of calves 6 and 9 months after inoculation with Mycobacterium avium subspecies paratuberculosis. Vet Res 45, 96 (2014) doi:10.1186/s13567-014-0096-5
- Canonical Pathway
- Mycobacterium Avium Subspecies Paratuberculosis
- Control Calf
- Mycobacterium Avium Subspecies Paratuberculosis Infection
- System Biology Analysis