- Research article
- Open Access
Characterization of paucibacillary ileal lesions in sheep with subclinical active infection by Mycobacterium avium subsp. paratuberculosis
- Salvatore Pisanu†1,
- Tiziana Cubeddu†2,
- Carla Cacciotto1,
- Ylenia Pilicchi2,
- Daniela Pagnozzi1,
- Sergio Uzzau1, 3,
- Stefano Rocca2 and
- Maria Filippa Addis1, 4Email authorView ORCID ID profile
© The Author(s) 2018
- Received: 7 September 2018
- Accepted: 13 November 2018
- Published: 4 December 2018
Paratuberculosis (PTB) or Johne’s disease is a contagious enteritis of ruminants caused by Mycobacterium avium subsp. paratuberculosis (MAP). Ovine PTB is less understood than bovine PTB, especially concerning paucibacillary infection and its evolution into clinical disease. We combined shotgun proteomics, histopathology and immunohistochemistry for the characterization of ileal tissues collected from seven asymptomatic sheep negative to serum ELISA, positive to feces and tissue MAP IS900 and F57 PCR, histologically classified as paucibacillary, actively infected, together with 3 MAP-free controls (K). Following shotgun proteomics with label-free quantitation and differential analysis, 96 proteins were significantly changed in PTB vs K, and were mostly involved in immune defense processes and in the macrophage-MAP interaction. Principal component analysis (PCA) of protein abundances highlighted two PTB sample clusters, PTB1 and PTB2, indicating a dichotomy in their proteomic profiles. This was in line with the PCA of histopathology data and was related to features of type 2 (PTB1) and type 3a (PTB2) lesions, respectively. PTB2 proteomes differed more than PTB1 proteomes from K: 43 proteins changed significantly only in PTB2 and 11 only in PTB1. The differential proteins cathelicidin, haptoglobin, S100A8 and S100A9 were evaluated by immunohistochemistry. K tissues were negative to cathelicidin and haptoglobin and sparsely positive to S100A8 and S100A9. PTB tissues were positive to all four proteins, with significantly more cells in PTB2 than in PTB1. In conclusion, we described several pathways altered in paucibacillary PTB, highlighted some proteomic differences among paucibacillary PTB cases, and identified potential markers for disease understanding, staging, and detection.
Paratuberculosis (PTB) or Johne’s disease is a chronic, contagious enteritis of ruminants caused by Mycobacterium avium subspecies paratuberculosis (MAP). PTB is especially relevant in farmed ruminants [1–4] for the economic consequences caused by increase in mortality, decrease in milk production, and weight loss [5, 6]. In addition, viable MAP can be found in pasteurized milk and milk products with a potential risk of zoonotic transmission [7–9].
MAP and its role in PTB have been the subject of numerous studies on disease progression and evolution in cattle. Our knowledge in this respect has increased considerably, but some aspects remain unclear also in this species, and more effective tools for disease diagnosis and control are still needed . MAP transmission can occur by the fecal–oral route, in utero, and by ingestion of contaminated colostrum or milk [10–15]. Once the bacterium reaches the intestine, it is taken up by M cells and translocated across the intestinal mucosa, where it is internalized by naive macrophages and can lead to persistent infection . In cattle, the disease typically goes from an early subclinical, paucibacillary phase to a later multibacillary phase with severe clinical manifestations of the disease that include wasting and profuse, watery diarrhea . In sheep, PTB is more insidious and classification is less clearly defined [2, 17], leading to a great underestimation of its worldwide diffusion. According to several authors, the largest number of sheep in an affected flock are infected but asymptomatic [18, 19]. These can either have a latent infection, encompassing the presence of MAP in tissues without signs of disease (subclinical infection), or an active infection, in which histological signs of the disease can be observed (subclinical disease) [20, 21]. A histological scoring system  classifies these signs into mild lesions represented by small focal granulomata of epithelioid cells limited to the Peyer patches (type 1) or extending to the adjacent mucosa (type 2), and more severe lesions with a multifocal cellular infiltration in mucosal areas not associated with lymphoid tissues and extending into the submucosa, with thickening of the mucosa and atrophy of villi (type 3a). According to Pérez and coworkers , symptomatic animals with the end-point disease (clinical disease) present two different types of lesion, lepromatous (type 3b) or tuberculoid (type 3c). Based on the degree of colonization, lesions can also be classified into paucibacillary, with few or no acid-fast bacilli (AFB), and multibacillary, with abundant AFB. Type 1, 2 and 3c lesions are paucibacillary, while type 3b are multibacillary. Type 3a lesions are mainly paucibacillary, but multibacillary patterns can also be present, indicating that a crucial “transition stage” may occur in animals with this type of lesions . Actively infected sheep with the subclinical disease can either remain such for their whole life, acting as MAP reservoirs and shedders, or develop clinical disease by showing either paucibacillary (3c) or multibacillary (3b) lesions. At present, factors, pathways, stages and dynamics of disease progression are not completely understood , and the best diagnostic approach remains post-mortem evaluation with histopathology image analysis, which represents the best indicator to confirm PTB and define its stage [1, 2, 23, 24].
The availability of detailed information on the proteomic alterations introduced in ovine ileal tissues by the subclinical, paucibacillary, active MAP infection would provide novel information for understanding disease mechanisms, discover protein markers with potential for disease staging and, possibly, indicate host proteins that may be shed and detected in feces in vivo. This would also help experimental infection studies aimed at understanding PTB progression to clinical disease. Currently, however, proteomic information on MAP-infected tissues is limited to a preliminary study recently published by our research group and describing the changes occurring in the ileum of sheep with multibacillary, clinical PTB .
In this study, we combined histopathology, proteomics, and immunohistochemistry for the characterization of paucibacillary MAP lesions found in actively infected, subclinical sheep. Ileal tissues of asymptomatic, serum ELISA-negative, feces and tissue PCR-positive sheep were subjected to histopathological analysis and those classified as paucibacillary PTB were analyzed by shotgun proteomics in comparison to MAP-free controls. Differential proteins were investigated for host pathways activated during infection. Finally, proteins of interest were evaluated by immunohistochemistry in MAP-infected tissues and in MAP-free controls.
Animals and tissues
Sheep ileal tissues used for the study belonged to a MAP-positive flock of Sarda sheep (N = 174). The flock was monitored by the farm veterinarian under an owner voluntary basis by clinical examination, evaluation of PTB symptoms, serological screening by ELISA for presence of anti-MAP antibodies in serum (IDEXX Laboratories, Inc., Westbrook, MA, USA), and IS900 and F57 PCR for presence of MAP in feces, as described previously . Out of 174 sheep, 21 were both ELISA and PCR-positive, 20 were ELISA-positive and PCR-negative, and 15 were ELISA-negative and PCR-positive. The latter 15 were identified as those most probably affected by paucibacillary paraTBC. At routine slaughtering, the intestinal packages belonging to these 15 sheep (all females between 3 and 4 years of age) were retrieved at the slaughterhouse and brought to the Department of Veterinary Medicine at the University of Sassari for gross pathological anatomy examination. Separate aliquots of tissue from the proximal, intermediate, and distal ileum (N = 45 tissue samples) were collected from each sheep. Matched tissue aliquots were frozen at −80 °C or formalin-fixed and paraffin embedded for downstream molecular and histopathological characterization.
Histopathological analysis and molecular characterization of ileal lesions
Three micrometre sections from paraffin-embedded tissue blocks from each ileal tract were subjected to haematoxylin–eosin and Ziehl–Neelsen (ZN) staining and examined to confirm presence of MAP and associated lesions. PCR was carried out on all 45 frozen tissues samples for confirmation of MAP infection, as described previously . Briefly, DNA was isolated with the DNeasy® Blood and Tissue kit (Qiagen, Hilden, Germany), and tested for presence of the MAP specific sequences IS900 and F57 by qualitative PCR. The MAP type was also characterized by IS1311, RFLP, and sequencing, as described previously [25, 26]. All 15 tested animals carried the sheep MAP strain (S). FFPE tissue sections obtained from the proximal, intermediate, and distal ileum of the 15 ELISA-negative, PCR-positive sheep (N = 45 tissue samples) were evaluated by histopathological analysis. Out of these 15 sheep, 7 showed histopathological lesions compatible with paucibacillary paraTBC and were selected for proteomic analysis (Additional file 1). Negative control tissues (K) were obtained from three sheep of a certified MAP-free flock. K sheep were negative to serum ELISA and feces PCR. Their intestinal tissues were negative to PCR, microscopical evaluation and Ziehl–Neelsen staining, and showed normal anatomy and histology.
Tissue processing for protein extraction and shotgun proteomic analysis
Tissue sample preparation was carried out as previously described . Briefly, 100 mg of tissue from the distal ileum of the 7 PTB and 3 K sheep (N = 10 tissue samples) was minced with a sterile knife, placed in Eppendorf safe-lock tubes (Eppendorf, Hamburg, Germany), immersed in extraction buffer and subjected to homogenization in a TissueLyser mechanical homogenizer (Qiagen, Hilden, Germany) followed by three cycles of sonication and freeze-thawing. Then, the extract was clarified by centrifugation, subjected to nucleic acid digestion, and protein concentration was determined. Extracts and residual pellets (N = 20) were processed by filter-aided sample preparation (FASP). Briefly, protein samples were subjected to reduction, alkylation, and trypsin digestion on Amicon Ultra-0.5 centrifugal filter units with Ultracel-10 membrane (Millipore, Billerica, MA, USA) to obtain sample digests for shotgun proteomic analysis. Peptide concentration of digests was determined by measuring absorbance at 280 nm with a NanoDrop 2000 spectrophotometer (Thermo Scientific, San Jose, CA, USA) using MassPREP E. coli Digest Standard (Waters, Milford, MA, USA) to create a calibration curve.
Shotgun proteomic analysis of peptides
All peptide mixtures were analyzed by liquid chromatography-tandem mass spectrometry (LC–MS/MS) in duplicate runs (two technical replicates for each peptide mixture) on a Q-Exactive interfaced with an UltiMate 3000 RSLCnanoLC system (Thermo Scientific, San Jose, CA, USA), as described previously . Four microgram of each peptide mixture were concentrated and washed onto a trapping precolumn (Acclaim PepMap C18, 75 µm × 2 cm nanoViper, 3 µm, 100 Å, Thermo Scientific) and fractionated on a C18 RP column (Acclaim PepMap RSLC C18, 75 µm × 50 cm nanoViper, 2 µm, 100 Å, Thermo Scientific) at a flow rate of 250 nL/min using a linear gradient of 245 min from 5 to 37.5% eluent B (0.1% formic acid in 80% acetonitrile) in eluent A (0.1% formic acid). Fragmentation occurred by higher energy collisional dissociation (HCD) and nitrogen as the collision gas. Proteome Discoverer (version 1.4; Thermo Scientific, Bremen, Germany) was used for protein identification using Sequest-HT as search engine. Each MS/MS spectrum was analyzed as follows. Database: Bos taurus, Ovis aries, and Mycobacterium avium downloaded from UniProtKB/Swiss-Prot (release 2017 06); enzyme: trypsin, with two missed cleavages allowed; precursor mass tolerance: 10 ppm; MS/MS tolerance: 0.02 Da; charge states: + 2, + 3, and + 4; cysteine carbamidomethylation as static modification and methionine oxidation as dynamic modifications. The percolator algorithm [28, 29] was used for assessing the significance of protein identification (P < 0.01) and for peptide validation (false discovery rate, FDR, < 0.01%) [28, 30]. Only rank 1 peptides and only proteins identified with at least two peptides and two spectral counts were considered. Peptide and protein grouping according to the Proteome Discoverer’s algorithm were allowed, applying the strict maximum parsimony principle. The protein list for each intestinal sample was built up by merging the data from the LC–MS/MS runs of the protein extract and of the residual pellet, generating 10 protein identification lists (1 for each sheep sample). For proteins having more than one entry, only those with the highest number of unique peptides and spectral counts (SpC) were selected for downstream analyses.
Principal component analysis and hierarchical clustering of samples
Principal component analysis (PCA) and hierarchical clustering of samples according to proteomic results were carried out based on the normalized spectral abundance factor (NSAF) values of all identified proteins, obtained according to Old et al. , as an indicator of protein abundance . PCA and hierarchical clustering of samples according to histopathological results were carried out based on the scores generated upon detailed histological assessment by two observers grading from 0 (absence of feature) to 3 (maximum degree of feature) and reported in Additional file 1. Data analysis was carried out with Perseus (v.22.214.171.124).
Differential proteomic analysis
To estimate protein abundance and to compare the levels of the same protein among sample groups, a spectral counting approach was applied as described by Old et al.  and Zybailov et al. . The RSC, that is the log2 ratio of protein abundances between two experimental groups, was calculated as described previously  and expressed as fold change (FC) . Statistical significance of differences in protein abundance was assessed by the beta-binomial test with FDR correction by Benjamini-Hochberg [35, 36]. Only proteins with FC ≥ 2 or ≤ −2, with a p-value ≤ 0.05 and present in at least two samples of each group being compared were considered significant.
To investigate the biological role of differential proteins, their functional association at a system level  was assessed by performing a knowledge-based protein–protein interaction network analysis by means of the biological interface and web resource STRING (version 10.5) after replacing all Ovis aries UniProt IDs with the corresponding Bos taurus UniProt IDs. Specifically, KEGG enrichment analysis was aimed at mapping the differential proteins in biological pathways, while gene ontology (GO) analysis was aimed at mapping them into biological processes (BP), molecular functions (MF), and cellular components (CC).
Immunohistochemistry (IHC) was carried out as described previously . Antibodies and dilutions were as follows: rabbit anti-CAMP (Sigma-Aldrich), 1:1500; rabbit anti-S100A8 (Sigma-Aldrich), 1:1000; rabbit anti-haptoglobin (Thermo Fisher Scientific), 1:100; rabbit anti-S100A9 (Sigma-Aldrich), 1:1000. Signal evaluation was carried out by counting positive cells in 10 random fields at 200× magnification (PC). Statistical analysis of IHC results was carried out with GraphPad Prism version 5.03 for Windows (GraphPad Software, La Jolla, CA). According to the Shapiro–Wilk normality test the data followed a non-normal distribution, and a nonparametric Mann–Whitney U-test was therefore applied.
Shotgun proteomic results, principal component analysis, and hierarchical clustering
Differential proteomic analysis
Differential proteins in all PTB, PTB1 and PTB2 vs K samples, respectively
All PTB vs K
PTB1 vs K
PTB2 vs K
Histone H2B type 1
Ubiquitin-60S ribosomal protein L40
Myosin light chain 1/3, skeletal muscle isoform
Small nuclear ribonucleoprotein-associated protein N
Keratin, type I cytoskeletal 17
Chaperone protein HtpG
Primary amine oxidase, liver isozyme
Fatty acid-binding protein, intestinal
Protein transport protein Sec61 subunit gamma
ORM1-like protein 2
Leukocyte surface antigen CD47
Ras-related protein Rab-5A
Leukocyte elastase inhibitor
Scavenger receptor cysteine-rich type 1 protein M130
Ras-related protein Rab-13
Allograft inflammatory factor 1
Platelet-activating factor acetylhydrolase 2
Serine beta-lactamase-like protein LACTB
Neutral amino acid transporter B (0)
Microsomal triglyceride transfer protein large subunit
Ribosome production factor 2 homolog
Glutathione S-transferase A1
Myeloid differentiation primary response protein MyD88
Small ubiquitin-related modifier 2
Signal transducer and activator of transcription 5A
Ras-related protein Rab-28
Histone H2B type 1-N
Keratin, type II cytoskeletal 73
Sodium/potassium-transporting ATPase subunit alpha-2
Proteasome subunit beta type-7
Chitinase-3-like protein 1
Bifunctional methylenetetrahydrofolate dehydrogenase
Neutrophil cytosol factor 2
Cation-independent mannose-6-phosphate receptor
Procollagen-lysine,2-oxoglutarate 5-dioxygenase 1
Vesicle-associated membrane protein-associated protein B
V-type proton ATPase subunit B
Biogenesis of lysosome-related organelles complex 1 sub 3
ERO1-like protein alpha
Placenta-specific protein 9
V-type proton ATPase catalytic subunit A
Four and a half LIM domains protein 3
Ribosome biogenesis protein BRX1 homolog
V-type proton ATPase subunit H
Serum amyloid A protein
Cell cycle control protein 50A
Notchless protein homolog 1
Mast cell protease 1A
Carbonic anhydrase 3
Intestinal-type alkaline phosphatase
Trefoil factor 3
Dynein light chain 2
Proteasome subunit beta type-6
Neural cell adhesion molecule 1
Histone H2B type 1-K
Inter-alpha-trypsin inhibitor heavy chain H4
ORM1-like protein 1
Translation machinery-associated protein 7
Elongation factor 1-alpha 2
Prostaglandin G/H synthase 1
Small nuclear ribonucleoprotein-associated protein B’
Tubulin alpha-3 chain
Tubulin alpha-1D chain
A total of 96 proteins showed statistically significant differences (Table 1). Of these, 69 were higher and 27 were lower in PTB vs K, respectively. Of the 69 increased proteins, 33 were statistically significant in all PTB samples, 6 only in PTB1, and 30 only in PTB2. Of the 27 decreased proteins, 9 were statistically significant in all PTB samples, 5 only in PTB1, and 13 only in PTB2 (Table 1). The higher number of differential proteins detected in PTB2 samples indicates that these differ more than PTB1 samples from K, and suggests the presence of two different levels of tissue involvement in the PTB sample set. Eleven proteins were identified in the MAP database and included in the differential analysis (Additional file 3). Among them, only chaperone protein HtpG passed the thresholds for FC value and statistical significance (Table 1).
KEGG enrichment and gene ontology analysis of differential proteins
Differential proteins were evaluated for protein networks by STRING. KEGG enrichment and GO results are summarized in Additional file 4.
IHC evaluation of selected differential proteins
Differential shotgun proteomics enables to investigate the changes occurring in the protein makeup of a tissue in a condition of interest, such as in a pathological state against the physiological state. The proteomic approach presents advantages and disadvantages in comparison to gene expression strategies. One advantage is the ability to truly evaluate the extent of protein abundance, going beyond the estimates based on gene expression. In fact, several proteins do not follow a strict relationship with gene transcription but are regulated at translational or post-translational level, as in the case of cathelicidins [38, 43]. Moreover, adding to changes in protein abundance, changes in protein localization can often be revealed including protein capture or release phenomena, provided that histological evaluations are carried out for further investigation or validation. On the other hand, the number of samples in a high-performance shotgun proteomics study is typically limited by analytical throughput and cost issues. In addition, the sensitivity of proteomics is not as high as transcriptomics, and some relevant low-abundance, transient mediators may not be detected, including cytokines.
In this study, the application of shotgun proteomics to paucibacillary MAP-infected sheep ileal tissues (PTB) and matched MAP-free controls (K) highlighted several host proteins and pathways that are altered by MAP upon infection. In addition, proteomic results and histopathological grading revealed the presence of two PTB sample clusters, one with a higher (PTB2) and one with a lower proteomic divergence (PTB1) from K. That is, several proteins associated with MAP infection and with its pathological process changed more intensely in the PTB2 than in the PTB1 sample cluster. This was also reflected by the different abundance of the MAP chaperone protein HtpG in the two clusters, with FC values of 5.32 in PTB1 and 11.19 in PTB2, respectively. Based on histological observations, the PTB1 and PTB2 proteomic profiles might relate to paucibacillary type 2 and type 3a lesions , indicating that the different levels of tissue alteration and lesion severity found in sheep with subclinical, active, paucibacillary MAP infection might be associated to specific changes in their proteomic profiles.
According to KEGG enrichment analysis, the differential proteins seen in intestinal tissues of PTB sheep were mainly involved in phagosome formation, lysosome function and tuberculosis. This is due to the MAP strategy of survival within macrophages: impairment of phagosome and lysosome fusion. The mechanism is shared with other mycobacteria including the more intensively investigated M. tuberculosis . Interestingly, most of these proteins changed significantly only in PTB2. Three out of eight proteins of the phagosome pathway were subunits of the V-type proton ATPase complex (V-ATPase). V-ATPase is responsible for phagolysosome acidification. By specifically excluding V-ATPase from this compartment, mycobacteria can survive and multiply within macrophages [10, 44]. Different authors have reported a higher expression of V-ATPase in MAP-infected macrophages as compared to macrophages infected with non-pathogenic mycobacteria [45–47]. This finding is in line with our previous work on multibacillary PTB, where we observed an increase in most V-ATPase subunits. In that work, consistently with a higher MAP load and lesion severity, fold changes were very high for all subunits and peaked up to a FC of 18.09 for subunit H, the specific target for mycobacterium-mediated exclusion [25, 44].
Changes in Rab GTPases were also observed. It is known that Rab5 stimulates fusion of early endosomes while Rab7 promotes fusion of mature phagosomes with endosomes and lysosomes . By retaining Rab5, MAP impairs maturation of endosomes into functional mycobactericidal compartments . In line with this, we observed that Rab5 was significantly increased in the PTB2 sample cluster (but not in PTB1). On the other hand, Rab7 was unchanged in both clusters, highlighting the role of Rab5 in the mycobacterial pathogenicity mechanism. The increase seen in Rab13 does also relate to the relevance of the phagocytotic pathway in the host defense from MAP. Tubulins are involved in phagosome maturation, and their massive decrease seen in infected tissues might be related to the impairment of this process in MAP-infected macrophages.
According to GO analysis, most proteins undergoing significant changes in MAP-infected tissues were involved in defense response, inflammatory response and acute phase response. Two isoforms of cathelicidin were increased in all PTB samples, with higher FC values in PTB2 vs PTB1. This result was also in line with our previous study on multibacillary PTB, (FCs of 24.53 for cathelicidin 1 and of 15.79 for cathelicidin 2) . Cathelicidins are well-known antimicrobial peptides with a crucial role in the intracellular killing of mycobacteria in macrophages [39, 50]. Besides their direct antimicrobial functions, cathelicidins play multiple roles as inflammation mediators , and their importance in the innate immune defense of ruminant is highlighted by the unusually large number of genes compared to a single copy in most other mammalians, such as humans and mice . Consistently with the proteomic data, IHC confirmed the absence of cathelicidin in K tissues and its presence in PTB tissues, with a significantly higher abundance in PTB2 vs PTB1. However, it was not possible to specifically identify the cathelicidin isoform(s) by IHC due to the polyclonal anti-cathelicidin antibodies used. Apparently, DCs were responsible for most of the cathelicidin signal; however, such DC localization might be due to protein capture rather than to protein expression by DCs themselves. Based on studies in humans, cathelicidin (LL-37) is internalized by DCs with subsequent localization primarily in the cytoplasmic compartment and then in the nucleus. This eventually leads to a suppression of their response to toll-like receptor ligands and renders DCs less capable to activate T lymphocytes. In other words, cathelicidins inhibit DC function [51, 53, 54]. Nevertheless, in human tuberculosis this is inverted by presence of vitamin D, in a finely regulated equilibrium between inhibition and stimulation of cathelicidin production, DC regulation and Th1 differentiation . This is only part of the puzzle governing the combined action of genetic background, farming variables and environmental conditions (such as nutritional balance and exposure to light) in resistance of ruminants to mycobacterial diseases, in which the role of cathelicidin regulation of DCs might deserve further investigation.
Several proteins involved in the acute phase response were significantly changed in PTB vs K. Interestingly, haptoglobin, S100A8 and serum amyloid protein A were seen only in PTB2 samples and were not detected in PTB1 or K. This might underscore a potential role for these proteins as markers of disease progression, possibly indicating the activation of an acute inflammatory process when MAP infection is less controlled by the host. Of note, haptoglobin showed the highest increase among all the proteins significant only in PTB2 vs K, with a FC of 13.84. In our previous study on multibacillary MAP infection, haptoglobin showed a FC of 22.93 . Here, upon IHC validation, haptoglobin was abundantly present in tissues of the PTB2 cluster while it was not detected in tissues of the PTB1 cluster. Notably, haptoglobin-positive cells appeared to be almost exclusively macrophages. This might be related to its scavenging functions; in fact, haptoglobin captures hemoglobin by forming a complex that is promptly internalized by macrophages. Its endocytosis and targeting to the lysosome are mediated by the macrophage receptor CD163 (scavenger receptor cysteine-rich type 1 protein M130) . Consistently with this observation, CD163 was significantly increased only in the PTB2 Cluster. Once again, the observed increase in protein abundance within the affected tissue is probably due to selective capture, rather than increased expression, by the cell. An increase in the CD163 receptor has been observed also in ileal tissues of cattle infected by MAP . It would be interesting to investigate how an increased availability of iron within the macrophage phagosomes is reflected on MAP pathogenicity.
S100A8 plays a prominent role in the regulation of inflammation and immune response. It can induce neutrophil chemotaxis and adhesion, and it is predominantly found as the calprotectin heterodimer in association with S100A9, with a wide plethora of intra- and extracellular functions . Fecal calprotectin is currently used in human medicine to differentiate organic from functional bowel disorders . Upon IHC validation, S100A8 was detected in all PTB samples, but few, scattered positive cells were also present in K tissues. The number of S100A8-positive cells was about three times higher in PTB2 than PTB1. In consideration of its known association with S100A9, it was surprising that the latter protein was not detected by proteomics in any of the experimental samples. However, S100A9 was readily detected by IHC in all PTB samples. The number of S100A9-positive cells was higher than that of S100A8-positive cells, and the cell type was also different. A possible explanation is that S100 proteins can also exist as homodimers , and a differential expression or a selective capture of S100A9 by a specific cell type might be occurring.
Another protein significantly increased in PTB was protein transport protein Sec61 (SEC61). As recently discovered, SEC61 is present in antigen-containing endosomes of DCs after stimulation with microbial substances and is essential for endosome-to-cytosol translocation and for antigen cross-presentation to T cells . This pathway enables DCs to present extracellular antigens onto MHC I molecules, thereby stimulating naive cytotoxic CD8+ T cells into activated cytotoxic CD8+ T cells. This is only part of the complex mechanisms responsible for the predominance of the Th1 or Th2 response, a crucial step in the control of mycobacterial infections.
H2B was the most increased protein in all PTB samples. This might relate to the well-known antimicrobial role of histones. On the other hand, ubiquitin was the second highest protein, and this might indicate the occurrence of an extensive histone ubiquitination . Nevertheless, it should be reminded that shotgun proteomics results as analyzed and presented here can only provide information on relative protein abundance. These data can give useful hints on protein networks and pathways activated upon infection, but indications on protein structure, protein interactions or post-translational modifications need dedicated experimental investigations or data analysis algorithms.
Adding to proteins increased in PTB2, proteins significantly increased only in samples of the PTB1 cluster (such as fatty acid-binding protein, glutathione S-transferase, Na/K-transporting ATPase, argininosuccinate synthase) or detected only in the PTB1 cluster (such as Rab-28), could also be worth investigating further; these proteins might represent potential indicators of a better ability to contrast disease and, should this be demonstrated, have potential as targets in genetic selection for resistance to mycobacterial diseases.
Finally, as anticipated, it should be reminded that this study was carried out on a limited number of animals due to methodological constraints. Stringent statistic procedures were applied in differential proteomic analysis and IHC validation confirmed the main study findings, but observations such as the relationship of PTB1 and PTB2 clusters with type 2 and type 3a lesions or any hypothesis on their connection with disease progression or increased pathogenicity and associated pathways will need to be investigated further with studies on a larger, well-characterized sample cohort by means of higher throughput techniques, such as IHC or other molecular tests, for a thorough validation on the entire range of lesions including also multibacillary cases. It should also be highlighted that low-abundance mediators such as cytokines were not identified in this study, as in most shotgun proteomics studies. However, these have been investigated in detail by other authors in paucibacillary PTB [17, 20, 61, 62].
In conclusion, this work provided the first detailed proteomic characterization of paucibacillary MAP-infected sheep ileal tissues in comparison to MAP-free tissues. The differential proteomic analysis combined with immunohistochemical validation highlighted several changes occurring in PTB tissues and provided novel molecular information on the distinct levels of tissue involvement that can be found within the asymptomatic, paucibacillary condition. Adding to useful insights into disease processes and pathways, we identified several proteins that are absent in healthy tissues and are present in MAP-infected tissues, with differences in abundance levels or localization possibly reflecting different degrees of lesion severity. These host protein markers can assist future investigations aimed at understanding PTB evolution and progression, either with classical post-mortem analysis of tissues or, hopefully, by enabling the development of immunoassay tools for their in vivo detection in feces.
The authors declare that they have no competing interests.
Design of study and experiments: SP, TC, SU, SR, MFA. Carrying out of laboratory and field activities: SP, TC, CC, YP, SR. Analysis of results and interpretation of data: SP, TC, SR, MFA. Manuscript drafting: SP, TC, MFA. Manuscript revision: CC, DP, SR, SU, MFA. All authors read and approved the final manuscript.
The authors wish to thank the farm veterinarian Francesco Salis for his valuable contribution to this work.
Ethics approval and consent to participate
The farmer was informed of the study and intestinal samples were collected at the slaughterhouse during routine procedures. Animals were not sacrificed for the purposes of the study but as a result of the food production procedures. The study did not influence in any way the killing of animals.
The study was financially supported by Sardegna Ricerche with the Grant Art. 26 2013, and with funds granted by the University of Sassari. The funding bodies had no role in the design of the study and collection, analysis, and interpretation of data and in writing the manuscript.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Hines SA, Buergelt CD, Wilson JH, Bliss EL (1987) Disseminated Mycobacterium paratuberculosis infection in a cow. J Am Vet Med Assoc 190:681–683PubMedGoogle Scholar
- Windsor PA (2015) Paratuberculosis in sheep and goats. Vet Microbiol 181:161–169View ArticlePubMedGoogle Scholar
- Carrigan MJ, Seaman JT (1990) The pathology of Johne’s disease in sheep. Aust Vet J 67:47–50View ArticlePubMedGoogle Scholar
- Corpa JM, Garrido J, García Marín JF, Pérez V (2000) Classification of lesions observed in natural cases of paratuberculosis in goats. J Comp Pathol 122:255–265View ArticlePubMedGoogle Scholar
- McGregor H, Abbott KA, Whittington RJ (2015) Effects of Mycobacterium avium subsp. paratuberculosis infection on serum biochemistry, body weight and wool growth in Merino sheep: a longitudinal study. Small Rumin Res 125:146–153View ArticleGoogle Scholar
- Bush RD, Windsor PA, Toribio JA (2006) Losses of adult sheep due to ovine Johne’s disease in 12 infected flocks over a 3-year period. Aust Vet J 84:246–253View ArticlePubMedGoogle Scholar
- Millar D, Ford J, Sanderson J, Withey S, Tizard M, Doran T, Hermon-Taylor J (1996) IS900 PCR to detect Mycobacterium paratuberculosis in retail supplies of whole pasteurized cows’ milk in England and Wales. Appl Environ Microbiol 62:3446–3452PubMedPubMed CentralGoogle Scholar
- Grant IR, Ball HJ, Neill SD, Rowe MT (1996) Inactivation of Mycobacterium paratuberculosis in cows’ milk at pasteurization temperatures. Appl Environ Microbiol 62:631–636PubMedPubMed CentralGoogle Scholar
- Monif GRG (2016) The Mycobacterium avium subspecies paratuberculosis dilemma. Biol Med (Aligarh) 8:270Google Scholar
- Arsenault RJ, Maattanen P, Daigle J, Potter A, Griebel P, Napper S (2014) From mouth to macrophage: mechanisms of innate immune subversion by Mycobacterium avium subsp. paratuberculosis. Vet Res 45:54View ArticlePubMedPubMed CentralGoogle Scholar
- Chiodini RJ, Van Kruiningen HJ, Merkal RS (1984) Ruminant paratuberculosis (Johne’s disease): the current status and future prospects. Cornell Vet 74:218–262PubMedGoogle Scholar
- Stabel JR (2000) Transitions in immune responses to Mycobacterium paratuberculosis. Vet Microbiol 77:465–473View ArticlePubMedGoogle Scholar
- Sweeney RW (1996) Transmission of paratuberculosis. Vet Clin North Am Food Anim Pract 12:305–312View ArticlePubMedGoogle Scholar
- Lambeth C, Reddacliff LA, Windsor P, Abbott KA, McGregor H, Whittington RJ (2004) Intrauterine and transmammary transmission of Mycobacterium avium subsp. paratuberculosis in sheep. Aust Vet J 82:504–508View ArticlePubMedGoogle Scholar
- Nebbia P, Robino P, Zoppi S, De Meneghi D (2006) Detection and excretion pattern of Mycobacterium avium subspecies paratuberculosis in milk of asymptomatic sheep and goats by nested-PCR. Small Rumin Res 66:116–120View ArticleGoogle Scholar
- Ponnusamy D, Periasamy S, Tripathi BN, Pal A (2013) Mycobacterium avium subsp. paratuberculosis invades through M cells and enterocytes across ileal and jejunal mucosa of lambs. Res Vet Sci 94:306–312View ArticlePubMedGoogle Scholar
- Gillan S, O’Brien R, Hughes AD, Griffin JFT (2010) Identification of immune parameters to differentiate disease states among sheep infected with Mycobacterium avium subsp. paratuberculosis. Clin Vaccine Immunol 17:108–117View ArticlePubMedGoogle Scholar
- Begara-Mcgorum I, Wildblood LA, Clarke CJ, Connor KM, Stevenson K, McInnes CJ, Sharp JM, Jones DG (1998) Early immunopathological events in experimental ovine paratuberculosis. Vet Immunol Immunopathol 63:265–287View ArticlePubMedGoogle Scholar
- Sergeant ESG, Marshall DJ, Eamens GJ, Kearns C, Whittington RJ (2003) Evaluation of an absorbed ELISA and an agar-gel immuno-diffusion test for ovine paratuberculosis in sheep in Australia. Prev Vet Med 61:235–248View ArticlePubMedGoogle Scholar
- Marquetoux N, Mitchell R, Ridler A, Heuer C, Wilson P (2018) A synthesis of the patho-physiology of Mycobacterium avium subspecies paratuberculosis infection in sheep to inform mathematical modelling of ovine paratuberculosis. Vet Res 49:27View ArticlePubMedPubMed CentralGoogle Scholar
- Whittington RJ, Begg DJ, de Silva K, Purdie AC, Dhand NK, Plain KM (2017) Case definition terminology for paratuberculosis (Johne’s disease). BMC Vet Res 13:328View ArticlePubMedPubMed CentralGoogle Scholar
- Pérez V, García Marín JF, Badiola JJ (1996) Description and classification of different types of lesion associated with natural paratuberculosis infection in sheep. J Comp Pathol 114:107–122View ArticlePubMedGoogle Scholar
- Kurade NP, Tripathi BN, Rajukumar K, Parihar NS (2004) Sequential development of histologic lesions and their relationship with bacterial isolation, fecal shedding, and immune responses during progressive stages of experimental infection of lambs with Mycobacterium avium subsp. paratuberculosis. Vet Pathol 41:378–387View ArticlePubMedGoogle Scholar
- Hailat NQ, Hananeh W, Metekia AS, Stabel JR, Al-Majali A, Lafi S (2010) Pathology of subclinical paratuberculosis (Johne’s Disease) in Awassi sheep with reference to its occurrence in Jordan. Vet Med (Praha) 55:590–602View ArticleGoogle Scholar
- Pisanu S, Cubeddu T, Uzzau S, Rocca S, Addis MF (2017) Proteomic changes in the ileum of sheep infected with Mycobacterium avium subspecies paratuberculosis. Vet J 219:1–3View ArticlePubMedGoogle Scholar
- Collins DM, De Zoete M, Cavaignac SM (2002) Mycobacterium avium subsp. paratuberculosis strains from cattle and sheep can be distinguished by a PCR test based on a novel DNA sequence difference. J Clin Microbiol 40:4760–4762View ArticlePubMedPubMed CentralGoogle Scholar
- Pisanu S, Biosa G, Carcangiu L, Uzzau S, Pagnozzi D (2018) Comparative evaluation of seven commercial products for human serum enrichment/depletion by shotgun proteomics. Talanta 185:213–220View ArticlePubMedGoogle Scholar
- The M, MacCoss MJ, Noble WS, Käll L (2016) Fast and accurate protein false discovery rates on large-scale proteomics data sets with Percolator 3.0. J Am Soc Mass Spectrom 27:1719–1727View ArticlePubMedPubMed CentralGoogle Scholar
- Käll L, Canterbury JD, Weston J, Noble WS, MacCoss MJ (2007) Semi-supervised learning for peptide identification from shotgun proteomics datasets. Nat Methods 4:923–925View ArticlePubMedGoogle Scholar
- The M, Tasnim A, Käll L (2016) How to talk about protein-level false discovery rates in shotgun proteomics. Proteomics 16:2461–2469View ArticlePubMedPubMed CentralGoogle Scholar
- Old WM, Meyer-Arendt K, Aveline-Wolf L, Pierce KG, Mendoza A, Sevinski JR, Resing KS, Ahn NG (2005) Comparison of label-free methods for quantifying human proteins by shotgun proteomics. Mol Cell Proteomics 4:1487–1502View ArticlePubMedGoogle Scholar
- Pagnozzi D, Biosa G, Addis MF, Mastrandrea S, Masala G, Uzzau S (2014) An easy and efficient method for native and immunoreactive Echinococcus granulosus antigen 5 enrichment from hydatid cyst fluid. PLoS One 9:e104962View ArticlePubMedPubMed CentralGoogle Scholar
- Zybailov B, Mosley AL, Sardiu ME, Coleman MK, Florens L, Washburn MP (2006) Statistical analysis of membrane proteome expression changes in Saccharomyces cerevisiae. J Proteome Res 5:2339–2347View ArticlePubMedGoogle Scholar
- Tanca A, Pagnozzi D, Burrai G, Polinas M, Uzzau S, Antuofermo E, Addis MF (2012) Comparability of differential proteomics data generated from paired archival fresh-frozen and formalin-fixed samples by GeLC-MS/MS and spectral counting. J Proteomics 77:561–576View ArticlePubMedGoogle Scholar
- Pham TV, Piersma SR, Warmoes M, Jimenez CR (2010) On the beta-binomial model for analysis of spectral count data in label-free tandem mass spectrometry-based proteomics. Bioinformatics 26:363–369View ArticlePubMedGoogle Scholar
- Kammers K, Cole RN, Tiengwe C, Ruczinski I (2015) Detecting significant changes in protein abundance. EuPA Open Proteomics 7:11–19View ArticlePubMedPubMed CentralGoogle Scholar
- Szklarczyk D, Franceschini A, Wyder S, Forslund K, Heller D, Huerta-Cepas J, Simonovic M, Roth A, Santos A, Tsafou KP, Kuhn M, Bork P, Jensen LJ, von Mering C (2015) STRING v10: protein-protein interaction networks, integrated over the tree of life. Nucleic Acids Res 43:D447–D452View ArticleGoogle Scholar
- Cubeddu T, Cacciotto C, Pisanu S, Tedde V, Alberti A, Pittau M, Dore S, Cannas A, Uzzau S, Rocca S, Addis MF (2017) Cathelicidin production and release by mammary epithelial cells during infectious mastitis. Vet Immunol Immunopathol 189:66–70View ArticleGoogle Scholar
- Sonawane A, Santos JC, Mishra BB, Jena P, Progida C, Sorensen OE, Gallo R, Appelberg R, Griffiths G (2011) Cathelicidin is involved in the intracellular killing of mycobacteria in macrophages. Cell Microbiol 13:1601–1617View ArticlePubMedPubMed CentralGoogle Scholar
- Zanetti M (2005) The role of cathelicidins in the innate host defenses of mammals. Curr Issues Mol Biol 7:179–196PubMedGoogle Scholar
- Gill N, Wlodarska M, Finlay BB (2011) Roadblocks in the gut: barriers to enteric infection. Cell Microbiol 13:660–669View ArticlePubMedGoogle Scholar
- Burri E, Beglinger C (2014) The use of fecal calprotectin as a biomarker in gastrointestinal disease. Expert Rev Gastroenterol Hepatol 8:197–210View ArticlePubMedGoogle Scholar
- Chromek M, Slamová Z, Bergman P, Kovács L, Podracká L, Ehrén I, Hökfelt T, Gudmundsson GH, Gallo RL, Agerberth B, Brauner A (2006) The antimicrobial peptide cathelicidin protects the urinary tract against invasive bacterial infection. Nat Med 12:636–641View ArticlePubMedGoogle Scholar
- Wong D, Bach H, Sun J, Hmama Z, Av-Gay Y (2011) Mycobacterium tuberculosis protein tyrosine phosphatase (PtpA) excludes host vacuolar-H+-ATPase to inhibit phagosome acidification. Proc Natl Acad Sci U S A 108:19371–19376View ArticlePubMedPubMed CentralGoogle Scholar
- Weiss DJ, Evanson OA, Moritz A, Deng MQ, Abrahamsen MS (2002) Differential responses of bovine macrophages to Mycobacterium avium subsp. paratuberculosis and Mycobacterium avium subsp. avium. Infect Immun 70:5556–5561View ArticlePubMedPubMed CentralGoogle Scholar
- Murphy JT, Sommer S, Kabara EA, Verman N, Kuelbs MA, Saama P, Halgren R, Coussens PM (2006) Gene expression profiling of monocyte-derived macrophages following infection with Mycobacterium avium subspecies avium and Mycobacterium avium subspecies paratuberculosis. Physiol Genomics 28:67–75View ArticlePubMedGoogle Scholar
- Weiss DJ, Evanson OA, de Souza C, Abrahamsen MS (2005) A critical role of interleukin-10 in the response of bovine macrophages to infection by Mycobacterium avium subsp. paratuberculosis. Am J Vet Res 66:721–726View ArticlePubMedGoogle Scholar
- Chua CEL, Gan BQ, Tang BL (2011) Involvement of members of the Rab family and related small GTPases in autophagosome formation and maturation. Cell Mol Life Sci 68:3349–3358View ArticlePubMedGoogle Scholar
- Fratti RA, Backer JM, Gruenberg J, Corvera S, Deretic V (2001) Role of phosphatidylinositol 3-kinase and Rab5 effectors in phagosomal biogenesis and mycobacterial phagosome maturation arrest. J Cell Biol 154:631–644View ArticlePubMedPubMed CentralGoogle Scholar
- Liu PT, Stenger S, Tang DH, Modlin RL (2007) Cutting edge: vitamin D-mediated human antimicrobial activity against Mycobacterium tuberculosis is dependent on the induction of cathelicidin. J Immunol 179:2060–2063View ArticlePubMedGoogle Scholar
- Agier J, Efenberger M, Brzezińska-Blaszczyk E (2015) Cathelicidin impact on inflammatory cells. Cent Eur J Immunol 40:225–235View ArticlePubMedPubMed CentralGoogle Scholar
- Bovine Genome Sequencing and Analysis Consortium, Elsik CG, Tellam RL, Worley KC, Gibbs RA, Muzny DM, Weinstock GM, Adelson DL, Eichler EE, Elnitski L, Guigó R, Hamernik DL, Kappes SM, Lewin HA, Lynn DJ, Nicholas FW, Reymond A, Rijnkels M, Skow LC, Zdobnov EM, Schook L, Womack J, Alioto T, Antonarakis SE, Astashyn A, Chapple CE, Chen HC, Chrast J, Câmara F, Ermolaeva O, Henrichsen CN et al (2009) The genome sequence of taurine cattle: a window to ruminant biology and evolution. Science 324:522–528View ArticleGoogle Scholar
- Bandholtz L, Ekman GJ, Vilhelmsson M, Buentke E, Agerberth B, Scheynius A, Gudmundsson GH (2006) Antimicrobial peptide LL-37 internalized by immature human dendritic cells alters their phenotype. Scand J Immunol 63:410–419View ArticlePubMedGoogle Scholar
- Kandler K, Shaykhiev R, Kleemann P, Klescz F, Lohoff M, Vogelmeier C, Bals R (2006) The anti-microbial peptide LL-37 inhibits the activation of dendritic cells by TLR ligands. Int Immunol 18:1729–1736View ArticlePubMedGoogle Scholar
- Rode AKO, Kongsbak M, Hansen MM, Lopez DV, Levring TB, Woetmann A, Ødum N, Bonefeld CM, Geisler C (2017) Vitamin D counteracts Mycobacterium tuberculosis-induced cathelicidin downregulation in dendritic cells and allows Th1 differentiation and IFNγ secretion. Front Immunol 8:656View ArticlePubMedPubMed CentralGoogle Scholar
- Graversen JH, Madsen M, Moestrup SK (2002) CD163: a signal receptor scavenging haptoglobin-hemoglobin complexes from plasma. Int J Biochem Cell Biol 34:309–314View ArticlePubMedGoogle Scholar
- Fernández M, Benavides J, Castaño P, Elguezabal N, Fuertes M, Muñoz M, Royo M, Ferreras MC, Pérez V (2017) Macrophage subsets within granulomatous intestinal lesions in bovine paratuberculosis. Vet Pathol 54:82–93View ArticlePubMedGoogle Scholar
- Ryckman C, Vandal K, Rouleau P, Talbot M, Tessier PA (2003) Proinflammatory activities of S100: proteins S100A8, S100A9, and S100A8/A9 induce neutrophil chemotaxis and adhesion. J Immunol 170:3233–3242View ArticlePubMedGoogle Scholar
- Zehner M, Marschall AL, Bos E, Schloetel JG, Kreer C, Fehrenschild D, Limmer A, Ossendorp F, Lang T, Koster AJ, Dübel S, Burgdorf S (2015) The translocon protein Sec61 mediates antigen transport from endosomes in the cytosol for cross-presentation to CD8+ T cells. Immunity 42:850–863View ArticlePubMedGoogle Scholar
- Chandrasekharan MB, Huang F, Sun ZW (2010) Histone H2B ubiquitination and beyond: regulation of nucleosome stability, chromatin dynamics and the trans-histone H3 methylation. Epigenetics 5:460–468View ArticlePubMedPubMed CentralGoogle Scholar
- Smeed JA, Watkins CA, Rhind SM, Hopkins J (2007) Differential cytokine gene expression profiles in the three pathological forms of sheep paratuberculosis. BMC Vet Res 3:18View ArticlePubMedPubMed CentralGoogle Scholar
- de Silva K, Plain K, Purdie A, Begg D, Whittington R (2018) Defining resilience to mycobacterial disease: characteristics of survivors of ovine paratuberculosis. Vet Immunol Immunopathol 195:56–64View ArticlePubMedGoogle Scholar