Skip to main content

Proteomic analysis of protein purified derivative of Mycobacterium bovis



Tuberculin skin test based on in vivo intradermal inoculation of purified protein derivative from Mycobacterium bovis (bPPD) is the diagnostic test for the control and surveillance of bovine tuberculosis (bTB).


Proteomic analysis was performed on different bPPD preparations from M. bovis, strain AN5. Proteins were precipitated from bPPD solutions by TCA precipitation. The proteome of bPPD preparations was investigated by bottom-up proteomics, which consisted in protein digestion and nano-LC–MS/MS analysis. Mass spectrometry analysis was performed on a Q-exactive hybrid quadrupole-Orbitrap mass spectrometer coupled online to an Easy nano-LC1000 system.


Three hundred and fifty-six proteins were identified and quantified by at least 2 peptides (99% confidence per peptide). One hundred and ninety-eight proteins, which had not been previously described, were detected; furthermore, the proteomic profile shared 80 proteins with previous proteomes from bPPDs from the United Kingdom and Brazil and 139 protein components from bPPD from Korea. Locus name of M. bovis (Mb) with orthologs from M. tuberculosis H37Rv, comparative gene and protein length, molecular mass, functional categories, gene name and function of each protein were reported. Ninety-two T cell mycobacterial antigens responsible for delayed-type hypersensitivity were detected, fifty-two of which were not previously reported in any bPPD proteome. Data are available via ProteomeXchange with identifier PXD005920.


This study represents the highest proteome coverage of bPPD preparations to date. Since proteins perform cellular functions essential to health and/or disease, obtaining knowledge of their presence and variance is of great importance in understanding disease states and for advancing translational studies. Therefore, to better understand Mycobacterium tuberculosis complex biology during infection, survival, and persistence, the reproducible evaluation of the proteins that catalyze and control these processes is critically important. More active and more specific tuberculins would be desirable. Indeed, many antigens contained within bPPD are currently responsible for the cross-reactivity resulting in false-positive results as they are shared between non-tuberculous and tuberculous mycobacteria.


Tuberculosis (TB), a zoonotic disease, is a major global human health problem, with 10.4 million new cases of active disease and nearly 1.8 million deaths estimated for 2015 [1]. The disease has similarly heavy consequences for a broad range of animal species thus being a recognized public veterinary health problem in many countries [2, 3].

Tuberculosis in bovines (bTB), caused predominantly by Mycobacterium bovis (M. bovis) a member of the Mycobacterium tuberculosis complex, is a disease still endemic in many countries [4]. bTB is the cause of significant economic hardship to the livestock industry with estimates of >50 million cattle infected worldwide [5] and is of zoonotic importance [6]. Indeed, although M. bovis is known to be the most important infectious agent responsible for bTB, however, it may cause tuberculosis in humans (hTB) both in developing and developed countries [79].

Moreover, bTB is subject to comprehensive control measures in order to limit both zoonotic transmission and economic losses. Such control is typically based on test-and-slaughter schemes, which require the accurate diagnosis of infected animals [10].

The diagnostic test for the control and surveillance of bTB used worldwide is the Tuberculin Skin Test (TST), which is based on in vivo intradermal inoculation of purified protein derivative from M. bovis (bPPD) alone or in combination with M. avium (aPPD). Those animals that react to PPD are isolated and slaughtered [11, 12]. Despite intensive eradication efforts over decades, bTB persists as a costly disease with adverse impacts on animal health and welfare, trade of animals and animal products, and livelihoods of producers, and continues to be a problem with global perspectives [4, 13]. It has been suggested that TST is a good herd test but a poor test for identifying individual infected animals [4]. Furthermore limitations in specificity and sensitivity of bPPD are additional factors contributing to the persistence of bTB [14]. However, TST is the gold standard for determining whether an individual animal is infected with bTB.

bPPD is a poorly characterized and ill-defined mix of proteins, lipids, and carbohydrates [11, 13] and little is known regarding what compounds are responsible for the delayed-type hypersensitivity (DTH) response [15, 16]. More defined knowledge on PPD composition and contribution of individual antigens in TST would give a better insight into the molecular mechanism behind the complex would therefore allow a better selection of proteins specific to M. tuberculosis [17]. Therefore, the identification of the molecular composition would facilitate the development of a more refined reagent [15].

A few proteomic studies have been performed on bPPD composition. Borsuk et al. [18] reported the first proteomic study from bPPD from the United Kingdom (UK) and from Brazil (BR). Cho et al. [19] described proteome profiles of bPPD from Korea (KR). More recently, Gcebe et al. [20] carried out a proteomic analysis of bPPD obtained from Prionics at The Netherlands.

The aim of the present paper is to report proteomic profiles detected on four bPPD preparations used for control and surveillance of bTB in most European countries.


bPPD samples

Four bPPD preparations from M. bovis, strain AN5, were examined. They were from Spain (S), manufactured by CZV company [CZ Veterinaria S.A., Porriño, (Pontevedra) Spain] and commercialized in Spain, Portugal, France, Germany, Ireland, Greece, United Kingdom, Belgium, Hungary, Bulgaria, Italy, Romania; two from Italy (one prepared by Istituto Zooprofilattico dell’Umbria e delle Marche, Perugia (Ip) and one by Istituto Zooprofilattico Sperimentale dell’Abruzzo e del Molise “G. Caporale”, Teramo (It), respectively) and commercialized in Italy and one from Netherlands (NL), manufactured by Prionics Lelystad BV, Lelystad, Holland and commercialized in Germany, Denmark, Sweden, Norway and Netherlands.


All chemicals were from Sigma (St. Louis, MO), unless otherwise specified.

Protein digestion and peptide dimethylation labelling

Proteins were precipitated from bPPD solutions by TCA precipitation. Briefly, proteins were precipitated in 10% TCA overnight, then pelleted at 12,000 g for 30 min at 4 °C. Pellets were washed with (1) ethyl ether (2) acetone. Proteins were resuspended in 100 mM triethyl ammonium bicarbonate buffer (TEAB) containing 0.2% SDS. Protein content was determined by the BCA protein assay using BSA as standard for the calibration curve (BCA protein assay, Thermo Scientific, Rockford, USA). For each sample, a 100 μg aliquot was subjected to protein reduction (10 mM DTT, 1 h at 37 °C) and alkylation (24 mM iodoacetamide, 1 h at 37 °C). The excess of iodoacetamide was quenched by additional 2 mM DTT. SDS concentration was brought to 0.05% by addition of HPLC water; then, 1 μg of proteomics grade trypsin was added, and digestion was allowed to proceed overnight at 37 °C.

Samples were labelled by either “light” (L), “medium” (M) or “heavy” (H) dimethyl labelling [21]. An aliquot of each sample, containing 25 μg of protein, was transferred to a separate Eppendorf vial. Then, 4 μL of either 4% (v/v) regular formaldehyde (L) or CD2O (M), or 13CD2O (H) plus 4 μL of either 0.6 M NaBH3CN (L, M) or NaBD3CN (H) were added. Reductive amination was allowed to proceed at room temperature for 1 h with shaking. To quench the reaction, 16 μL of 1% (v/v) ammonia solution and 8 μL of 5% formic acid were added to the samples. Samples were labelled in duplicates with label swapping, as follows: sample S (L, M, H), sample It (M, L), sample NL (H, L), sample Ip (H, M). Triplets were obtained by combining 5 μg of each labelled sample as follows: (1) S(L):It(M):NL(H); (2) It(L):S(M):Ip(H); (3) NL(L):Ip(M):S(H). Sample S was used as reference sample: thus, it was present in all triplets. Each mix was fractionated by SCX StageTips [22] using Empore SCX Extraction Disks. Briefly, samples were diluted 20-fold in solution A (0.5% formic acid, FA in 80% ACN), and then loaded onto a column prepared from a 10 μL micropipette tip stacked with two layers of the SCX resin, previously conditioned with 20 µl of solution A. The plugs were washed twice with 20 μL of solution A. Then, peptides were stepwise eluted by delivering six 14 μL aliquots of eluent of increasing ionic strength. The first five eluent solutions contained 20% acetonitrile and 0.5% formic acid (v/v) plus the following amount of ammonium acetate: (1) 50 mM, (2) 75 mM, (3) 100 mM, (4) 150 mM, (5) 250 mM. The sixth eluent solution was composed of 20% acetonitrile and 500 mM ammonium acetate. Eluates were evaporated to dryness and resuspended in 12 μL of mobile phase A (see below).

Nano-LC–MS/MS and data analysis

Mass spectrometry analysis was performed on a Q-Exactive Hybrid Quadrupole-Orbitrap Mass Spectrometer coupled online to an Easy nano-LC1000 system (Thermo Fisher Scientific, Germany). Peptide mixtures (3 µL) were loaded at a flow rate of 500 nL/min onto a silica capillary tip (75 µm i.d., length 10 cm) packed in house with 3 µm C18 silica particles (Dr. Maisch, Germany). Gradient elution was from 8% B (0,1% FA in 80% ACN) to 35% in 55 min, then from 35% B to 100% B in 5 min. Column equilibration (20 min) was at 2% B. MS acquisition was performed in positive ion mode with a nanoelectrospray voltage of 1800 V. Mass spectra were obtained in data-dependent acquisition (DDA) mode, using a top 12 method consisting in a survey full scan across the m/z range 350–1800, followed by MS/MS scans on the twelve most intense precursor ions accumulated for a maximum of 60 ms. Full scan acquisition parameters were: 70,000 FWHM, AGC target 1e6, maximum IT 50 ms, scan range 350 to 1800 m/z. Instead, dd-MS2 acquisition parameters were set as following: 17 500 FWHM, AGC target 1e5, maximum IT 60 ms, isolation window 1.6 m/z, scan range 200 to 2000 m/z. Collision energy was set at 25%.

Raw data were processed by Proteome Discoverer 1.4 using the Sequest algorithm and searched against the protein sequence database of Mycobacterium tuberculosis Complex (June 2015, 6415 proteins). Search criteria were set as follows: enzyme trypsin, maximum two missed cleavages, Dimethyl (Any N-Terminus), Dimethyl (K), Carbamidomethyl (C) as static modifications (dimethyl modifications were set to either “light”, “medium” or “heavy” in three parallel searches), Oxidation (Met) as dynamic modification, MS tolerance 10 ppm, MS/MS tolerance 0.02 Da. Search results were filtered by q values using Percolator [23], integrated in Proteome Discoverer (q value <0.01). As negative examples for the classifier, Percolator peptide hits derived from searching a decoy database composed of reversed protein sequences were used. Protein hits based on two successful peptide identifications in at least 2 out 3 LC-MS/MS data sets were considered valid. Initial quantification was performed in Proteome discoverer using default parameters for triplex dimethyl labelling. Advanced quantification, such as data normalization and permutation statistical analysis of the peptide ratios was performed using quantitative proteomics p value calculator (QPPC) [24]. Number of permutations was set to 10,000. Minimum number of observations was 2; thus, proteins identified by two peptides but quantified by a single unique peptide were not assigned a fold change. For each binary comparison (It versus S, NL versus S, Ip versus S) proteins whose fold change was either >2.0 or <0.5, with an associated p value <0.05 in both replicate analysis were considered differentially abundant. Finally, BoviList and TubercuList databases ( provided information on annotated M. bovis and M. tuberculosis H37Rv genes and proteins, including molecular weights and functional annotation as well as orthologous genes of these two highly related strains. Uniprot ( and KEGG ( databases have also been utilized in order to obtain further detailed information about proteins of M. bovis, strain AN5 and M. tuberculosis, strain H37Rv.

The mass spectrometry proteomics data have been deposited to the ProteomeXchange Consortium via the PRIDE [25] partner repository with the dataset identifier PXD005920.


Identification of bPPD proteins by list-based nano-LC-MS/MS

The proteome of four bPPD preparations was investigated by bottom-up proteomics, which consisted in protein digestion and nano-LC-MS/MS analysis. Altogether, we identified 356 proteins by 2 or more peptides (99% confidence per peptide) in at least 2 out 3 LC-MS/MS data sets (Fig. 1). A vast majority (~75%) of these proteins was found to be shared among four bPPDs. Our proteomic profile was composed of 198 not previously described proteins; furthermore, nineteen proteins were found to be shared with bPPD UK and bPPD BR proteome [18] and seventy-eight with bPPD KR [19]. Sixty one proteins were found in common both with bPPD UK as well as BR and bPPD KR proteomes (Fig. 2). Additional file 1: Table S1 reports locus names of M. bovis (Mb) and their orthologs from M. tuberculosis H37Rv (Rv) including comparative gene and protein lengths, molecular mass, functional categories, gene name and protein function of both mycobacterium strains. Finally, Additional file 1: Table S1 shows a comparative protein expression profile with two previous bPPD proteomic analysis [18, 19]. Furthermore, our proteomic profile shared 158 proteins with the proteome of M. tuberculosis PPD (MtbPPD) investigated by Prasad et al. [17] who were able to identify 265 proteins, the largest number of proteins found in MtbPPD so far.

Fig. 1
figure 1

Total number of proteins identified within four bPPD preparations

Fig. 2
figure 2

Venn diagram showing the unique and shared proteins detected in this study compared with previous bPPD proteomic profiles

BoviList functional classification codes

The proteins of bPPDs were grouped into eight functional categories, with the most representative (33%) being the functional category 7 (intermediary metabolism and respiration). The other proteins were grouped into the functional categories that include 0 (virulence, detoxification, adaptation—3%), 2 (lipid metabolism—5%), 3 (information pathways—20%), 4 (cell wall and cell processes—13%), 9 (regulatory proteins—4%), 10 (conserved hypotheticals—19%) and 16 (conserved hypotheticals with an orthologue in M. tuberculosis—2%). From a comparative point of view, all proteins we have identified in bPPDs have protein equivalents from M. tuberculosis H37Rv (, but some of them have a different structure and appear to belong to different categories (see Additional file 1: Table S1—red and blue colors, respectively). Accordingly, protein equivalents from M. tuberculosis H37Rv were grouped in seven functional categories.

Mycobacterial proteins containing possible T-cell antigens

Ninety-two mycobacterial proteins containing possible T-cell antigens and responsible for delayed-type hypersensitivity (DTH) have been identified. M. tuberculosis, strain H37Rv, and M. bovis, strain AN5, are highly related. Accordingly, some proteins including Mb0448, Mb1918c, Mb2002c, Mb3789, Mb3834c, Mb3904, Mb3905, just like their orthologs from M. tuberculosis H37RV, could have more than 20 different T cell epitopes since they are known to be evolutionarily hyperconserved [26, 27]. Numerous proteins such as Mb0584, Mb1767, Mb2056c, Mb2057c, Mb2058, Mb 2493c, Mb2656, Mb2657, Mb2659c, Mb3155, Mb3157c, referred to as latency antigens involved in the latent infection, appear to be as a part of the so-called dormancy (DosR) regulon, the expression of which is observed as part of adaptive response of M. tuberculosis complex to hypoxia [28, 29].

Fifty-two T cell antigens including Mb1301c (gene name: lprA), Mb1868c (gene name: glcB), Mb1950 (gene name: aceA), Mb3646c (gene name: espA), Mb3911c (gene name: espB) recently identified from M. tuberculosis, strain H37RV, via throughput proteome screening [30], have not been reported previously in bPPDs. Fourteen of them belong to a novel series of in vivo expressed M. tuberculosis (IVE-TB) T cell antigens [31, 32]. Two reactivation associated antigens, namely Mb0391c (gene name: clpB) and Mb2492c (gene name: rplB) and three lipoprotein T cell antigens, namely Mb0959 (gene name: pstS1), Mb0956c (gene name: pstS2) and Mb0951 (gene name: pstS3), known to generate very high levels of cytokine secretion [33], were also detected. Nine antigens were in common with bPPD KR only [19]. Thirty-one antigens were found to be shared both with bPPDs from UK as well as BR [18] and with bPPD from KR [19]. Our results are summarized in Fig. 3.

Fig. 3
figure 3

Venn diagram showing the unique and shared T cell mycobacterial antigens seen in this study compared to T cell antigens found in other bPPD proteomic analysis

Semi-quantitative analysis of overall protein abundance and quantitative comparison of PPDs by dimethyl labelling

Following database searches, qualitative (Additional file 2: Table S2) and relative quantitative information (Additional file 3: Table S3) was obtained. In Additional file 2: Table S2, proteins are sorted by the parameter “% total weighted spectra”. This is a semi-quantitative measure of protein abundance within a mixture (expressed in estimated % weight), which is based on normalized spectral counting [34]. It is worthwhile noting that the top ten proteins in the list comprise two proteins, namely: ESAT-6-like protein EsxN and Immunogenic protein MPB63, not previously reported. These ten most abundant proteins, whose list is reported in Table 1, constituted over 67% of total weighted spectra.

Table 1 The most abundant proteins in PPD preparations

Besides estimating average protein abundance by spectral counting, accurate quantification by mass spectrometry based on isotopic labelling was used to compare protein levels in each mixture. Proteins whose relative abundance changed, in both replicates, significantly (p value <0.05) by twofold or more in samples NL, It, Ip with respect to the reference sample S, are listed in Additional file 4: Table S4. Comparative analysis by dimethyl labelling showed that PPD preparations did not generally differ significantly in terms of protein composition, especially for what concerns the most abundant proteins. Two notable exceptions were: (1) protein 14 kDa antigen, whose levels in samples NL and It were comparable to sample S (fold change close to 1), but was found at 0.4-fold change in sample Ip, and (2) chaperonin GroES, whose levels were comparable in samples Ip and It relative to sample S (fold change close to 1 for Ip, fold change = 0.3 for sample It, but not found statistically significant in both replicates), but was found at relative fold change 0.3 in sample NL (Table 1). In total, with respect to reference sample S, 44, 34 and 30 proteins were found at significantly different levels in samples It, NL and Ip, respectively. These proteins were mostly present at very low levels in all preparations.


Comparative proteomics was performed on four bPPD preparations. Proteomic profile of all bPPDs was composed of 356 proteins thus representing the highest proteome coverage of bPPD preparations to date. We deciphered 198 new, never previously reported proteins in the protein expression profile. Altogether, 512 proteins of M. bovis PPDs, strain AN5, have been identified so far. Recently, 132 protein components were also revealed from a commercial bPPD preparation [20]. Unfortunately, this last study did not report the complete list of proteins identified by Mb number; therefore, it was not possible to carry out any comparative studies.

The worldwide used diagnostic test for the control and surveillance of bTB is the TST, based on the detection of cell mediated immunity under the exposure to bPPD, the composition of which is highly complex and remains ill-defined [4, 12].

Since proteins perform cellular functions essential to health and/or disease, obtaining knowledge of their presence and variance is of great importance in understanding disease states and for advancing translational studies [35]. Therefore, to better understand M. tuberculosis complex biology during infection, survival, and persistence, the reproducible evaluation of the proteins that catalyze and control these processes is critically important [36].

The emerging field of proteomics has contributed greatly to improving our understanding of the M. tuberculosis complex quite recently. Proteomics is currently in transition from pure basic research to medical application [37].

All proteins we detected in bPPDs have human equivalents. All bovine and human proteins were characterized by a very strong similarity and a remarkable identity. Indeed, all of them shared characteristics such as gene and protein length as well as molecular mass and function. However, we showed that some bovine and human proteins encoded by ortholog genes belong to different functional categories. Accordingly, it has been shown that some genes of members of M. tuberculosis complex can have functional polymorphisms and encoded proteins responsible for some phenotypic differences between M. bovis and the other members of M. tuberculosis complex [38].

The proteomic content of the four PPD preparations was characterized by a remarkable presence of chaperone proteins such as HspX, DnaK, GroEs, GroEl. These proteins are known to share a high homology (upwards of 70%) and are conserved amongst most mycobacterial species [18, 19]. They are believed to be the main proteins of the current diagnostic test responsible for high level of false positive responses [15].

A quantitative but not qualitative difference in protein content was seen in our bPPD proteomes. It is worthwhile noting that subtle differences in culture conditions, sterilization methods, protein precipitation methods, peptide fractionation process, trypsin efficiency may result in differences in proteomic profiles of bPPDs [19].

It has been suggested that proteomic analysis of different bPPD preparations could improve current diagnostic tests and gain insights into the immune response seen in TB disease [19].

We detected 92 mycobacterial proteins potentially involved in DTH, thus deciphering 52 new antigens not previously reported in bPPDs. Thirty six proteins were identified as T cell mycobacterial antigens in bPPDs from UK and from BR [18]; bPPD KR proteomes revealed 46 proteins playing a role in cellular immunity [19]. Altogether, ninety-eight mycobacterial antigens which play a central role in DTH have been identified in all examined bPPDs. Furthermore, from a comparative point of view, a very large number of T cell antigens was found to be shared between bPPDs of this study and MtbPPDs [26, 2830, 32, 3941].

The number of novel T cell mycobacterial antigens is increasingly detected in bPPD. There is a need to establish a better and more detailed understanding of T-cell biology through comparative investigations to decipher immune mechanisms that control mycobacterial infections which appear to rely heavily on the cellular immune system [27, 42]. Accordingly, more active and more specific tuberculins would be desirable. To date, limited progress has been achieved in this field, mainly because of the ill-defined nature of the antigens present in tuberculins as well as the complexity of PPD production.

We believe that developing more proper and defined antigens will be crucial to increase specificity and sensitivity of PPDs. It is worthwhile remembering that there is the urgency to improve specificity of PPD since as cross-reactive responses to bPPD (false-positive results) may occur as many antigens contained within bPPD are shared between non-tuberculous and tuberculous mycobacterial [4, 11].


Improvement and/or replacement of classical PPD composition with novel more specific reagents remain demanding [15].

To control TB disease, new strategies are needed to prepare better tools based on identification of novel antigens important for developing diagnostic tests which could be not only more accurate and sensitive but also capable of differentiating infected and uninfected vaccinated animals (“DIVA” tests), thus offering perspectives to introduce potential vaccination within existing eradication programs [43, 44].

Indeed, understanding immunity to M. bovis is a continuing challenge and one that is of interest to the fields of human and animal medicine alike [6]. It is important to note that ruminant also have greater similarity to humans than do the often-used mouse model with regard to several aspects of immune physiology, increasing the utility of cattle as a model [45]. Therefore, comparative immunology studies will continue to provide mutual benefit to TB research in both man and animals [46].

Accordingly, the One Health approach, believed to be a critical necessity to address zoonotic diseases [2], is clearly warranted for tuberculosis. The One Health concept is a worldwide strategy for expanding interdisciplinary collaborations and communications in all aspect of health care for humans, animals and the environment which, as far as tuberculosis disease is concerned, can speed the development of new diagnostic tests for humans and livestock thus improving tuberculosis surveillance, control, and eradication programs.





tuberculosis in humans


bovine TB


bovine purified protein derivative


avium purified protein derivative


bPPD from Brazil


bPPD from United Kingdom


bPPD from Korea


delayed type hypersensitivity


M. tuberculosis PPD


tuberculin Skin Test

Ip :

Istituto Zooprofilattico dell’Umbria e delle Marche, Perugia

It :

Istituto Zooprofilattico Sperimentale dell’Abruzzo e del Molise, Teramo


liquid chromatography tandem mass spectrometry






  1. World Health Organization (WHO). Global tuberculosis report 2016. Geneva: World Health Organization; 2016. p. 5.

    Google Scholar 

  2. Kaneene JB, Miller RA, Kaplan B, Steel JH, Thoen CH. Preventing and controlling zoonotic tuberculosis: a One Health approach. Vet Ital. 2014;50:7–22.

    PubMed  Google Scholar 

  3. Lewerin SS. Tuberculosis and One Health—what is in a Name? Front Vet Sci. 2015;2:54. doi:10.3389/fvets.2015.00054.

    Google Scholar 

  4. Schiller I, Oesch B, Vordermeier HM, Palmer MV, Harris BN, Orloski KA, Buddle BM, Thacker TC, Lyashchenko KP, Waters WR. Bovine tuberculosis: a review of current and emerging diagnostic techniques in view of their relevance for disease control and eradication. Transbound Emerg Dis. 2010;57:205–20.

    CAS  PubMed  Google Scholar 

  5. Tuggle CK, Waters WR. Tuberculosis-resistant transgenic cattle. Proc Natl Acad Sci USA. 2015;112:3854–5.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  6. McGill JL, Sacco RE, Baldwin CL, Telfer JC, Palmer MV, Waters WR. The role of gamma delta T cells in immunity to Mycobacterium bovis infection in cattle. Vet Immunol Immunopathol. 2015;159:133–43.

    Article  Google Scholar 

  7. Müller B, Dürr S, Alonso S, Hattendorf J, Laisse CJM, Parsons SDC, van Helden PD, Zinsstag J. Zoonotic Mycobacterium bovis-induced tuberculosis in humans. Emerg Infect Dis. 2013;19:899–908.

    Article  PubMed  PubMed Central  Google Scholar 

  8. El-Sayed A, El-Shannat S, Kamel M, Castañeda-Vazquez MA, Castañeda-Vazquez H. Molecular epidemiology of Mycobacterium bovis in humans and cattle. Zoonoses Public Health. 2016;63:251–64.

    Article  CAS  PubMed  Google Scholar 

  9. Scott C, Cavanaugh JS, Pratt R, Silk BJ, LoBue P, Moonan PK. Human tuberculosis caused by Mycobacterium bovis in the United States, 2006–2013. Clin Infect Dis. 2016;63:594–601.

    Article  PubMed  Google Scholar 

  10. Parsons SDC, McGill K, Doyle MB, Goosen WJ, van Heldem PD, Gormley E. Antigen-specific IP-10 release is a sensitive biomarker of Mycobacterium bovis infection in cattle. PLoS ONE. 2016;11(5):e0155440. doi:10.1371/journal.pone.0155440.

    Article  PubMed  PubMed Central  Google Scholar 

  11. Chen S, Parlane NA, Lee J, Wedlock N, Buddle BM, Rehm BHA. New skin test for detection of bovine tuberculosis on the basis of antigen-displaying polyester inclusions produced by recombinant Escherichia coli. Appl Environ Microbiol. 2014;80:2526–35.

    Article  PubMed  PubMed Central  Google Scholar 

  12. Eirin ME, Macias A, Magnano G, Morsella C, Mendez L, Blanco FC, Bianco MV, Severina W, Alito A, de los AngelesPando M, Singh M, Spallek R, Paolicchi FA, Bigi F, Cataldi AA. Identification and evaluation of new Mycobacterium bovis antigens in the in vitro interferon gamma release assay for bovine tuberculosis diagnosis. Tuberculosis. 2015;95:795–801.

    Article  CAS  PubMed  Google Scholar 

  13. Waters WR, Pallmer MV, Stafne MR, Bass KE, Maggioli MF, Thacker TC, Linscott R, Lawrence JC, Nelson JT, Esfandiari J, Greenwald R, Lyashchenko KP. Effects of serial skin testing with purified protein derivative on the level and quality of antibodies to complex and defined antigens in Mycobacterium bovis-infected cattle. Clin Vaccine Immunol. 2015;22:641–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  14. Karolemeas K, de la Rua-Domenech R, Cooper R, Goodchild AV, Clifton-Hadley RS, Conlan AJ. Estimation of the relative sensitivity of the comparative tuberculin skin test in tuberculous cattle herds subjected to depopulation. PLoS ONE. 2012;7(8):e43217. doi:10.1371/journal.pone.0043217.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  15. Yang H, Kruh-Garcia NA, Dobos KM. Purified protein derivatives of tuberculin: past, present, and future. FEMS Immunol Med Microbiol. 2012;66:273–80.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  16. Flores-Villalva S, Rogríguez-Hernández E, Rubio-Venegas Y, Cantó-Alarcón JG, Milián-Suazo F. What can proteomics tell us about tuberculosis? J Microbiol Biotechnol. 2015;25:1181–94.

    Article  CAS  PubMed  Google Scholar 

  17. Prasad TSK, Verma R, Kumar S, Nirujogi RS, Sathe GJ, Madugundu AK, Sharma J, Puttamallesh VN, Ganjiwale A, Myneedu VP, Chatterjee A, Pandey A, Harsha HC, Narayana J. Proteomic analysis of purified protein derivative of Mycobacterium tuberculosis. Clin Proteom. 2013;10:8–16.

    Article  Google Scholar 

  18. Borsuk S, Newcombe J, Mendum TA, Dellagostin OA, McFadden J. Identification of proteins from tuberculin purified protein derivative (PPD) by LC-MS/MS. Tuberculosis. 2009;89:423–30.

    Article  CAS  PubMed  Google Scholar 

  19. Cho YS, Jang YB, Lee SE, Cho JY, Ahn JM, Hwang I, Heo E, Nam HM, Cho D, Her M, Jean YH, Jung SC, Kim JM, Lee HS, Lee K, Belisle JT. Short communication: proteomic characterization of tuberculin purified protein derivative from Mycobacterium bovis. Res Vet Sci. 2015;101:117–9.

    Article  CAS  PubMed  Google Scholar 

  20. Gcebe N, Michel A, van Pittius NCG, Rutten V. Comparative genomics and proteomic analysis of four non-tuberculous Mycobacterium species and Mycobacterium tuberculosis complex: occurrence of shared immunogenic proteins. Front Microbiol. 2016;7:795. doi:10.3389/fmicb.2016.00795.

    Article  PubMed  PubMed Central  Google Scholar 

  21. Boersema PJ, Raijmakers R, Lemeer S, Mohammed S, Heck AJR. Multiplex peptide stable isotope dimethyl labeling for quantitative proteomics. Nat Protoc. 2009;4:484–94.

    Article  CAS  PubMed  Google Scholar 

  22. Rappsilber J, Mann M, Ishihama Y. Protocol for micro-purification, enrichment, pre-fractionation and storage of peptides for proteomics using StageTips. Nat Protoc. 2007;2:1896–906.

    Article  CAS  PubMed  Google Scholar 

  23. Spivak M, Weston J, Bottou L, Käll L, Noble WS. Improvements to the percolator algorithm for peptide identification from shotgun proteomics data sets. J Proteome Res. 2009;8:3737–45.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  24. Chen D, Shah A, Nguyen H, Loo D, Inder KL, Hill MM. Online quantitative proteomics p-value calculator for permutation-based statistical testing of peptides ratios. J Proteome Res. 2014;13:4184–91.

    Article  CAS  PubMed  Google Scholar 

  25. Vizcaíno JA, Csordas A, del-Toro N, Dianes JA, Griss J, Lavidas I, Mayer G, Perez-Riverol Y, Reisinger F, Ternent T, Xu QW, Wang R, Hermjakob H. 2016 update of the PRIDE database and related tools. Nucleic Acids Res. 2016;44(D1):D447–56.

    Article  PubMed  Google Scholar 

  26. Comas I, Chakravartti J, Small PM, Galagan J, Niemann S, Kremer K, Ernst JD, Gagneux S. Human T-cell epitopes of Mycobacterium tuberculosis are evolutionarily hyperconserved. Nat Genet. 2010;42:498–503.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  27. Axelsson-Robertson R, Magalhaes I, Parida SK, Zumla A, Maeure M. The immunological footprint of Mycobacterium tuberculosis T-cell epitope recognition. J Infect Dis. 2012;205:S301–15.

    Article  CAS  PubMed  Google Scholar 

  28. Leyten EMS, Lin MY, Franken KLMC, Friggen AH, Prins C, van Meijgaarden KE, Voskuil MI, Welding K, Andersen P, Schoolnik GK, Arend SM, Ottenhoff THM, Klein MR. Human T-cell responses to 25 novel antigens encoded by genes of the dormancy regulon of Mycobacterium tuberculosis. Microbes Infect. 2006;8:2052–60.

    Article  CAS  PubMed  Google Scholar 

  29. Gideon HP, Wilkinson KA, Rustad TR, Oni T, Guio H, Sherman DR, Vordermeier HM, Robertson BD, Young DB, Wilkinson RJ. Bioinformatic and empirical analysis of novel hypoxia-inducible targets of the human antituberculosis T cell response. J Immunol. 2012;189:5867–76.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  30. Nayak K, Jing L, Russell RM, Davies DH, Hermanson G, Molina DM, Liang X, Sherman DR, Kwok WW, Yang J, Kenneth J, Ahamed SF, Chandele A, Murali-Krishna K, Koelle DM. Identification of novel Mycobacterium tuberculosis CD4 T-cell antigens via throughput proteome screening. Tuberculosis. 2015;95:275–87.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  31. Serra-Vidal MM, Latorre I, Franken KLCM, Diaz J, de Souza-Galvão ML, Casas I, Maldonado J, Milà C, Solsona MC, Jimenez-Fuentes MA, Altet N, Lacoma A, Ruiz-Manzano J, Ausina B, Prat C, Ottenhoff TH, Dominguez J. Immunogenicity of 60 novel latency-related antigens of Mycobacterium tuberculosis. Front Microbiol. 2014;5:517. doi:10.3389/fmicb.2014.00517.

    Article  PubMed  PubMed Central  Google Scholar 

  32. Coppola M, van Meijgaarden KE, Franken KLMC, Commandeur S, Dolganov G, Kramnik I, Schoolnik GK, Comas I, Lund O, Prins C, van den Eeden SJF, Korsvold GE, Oftung F, Geluk A, Ottenhoff THM. New genome-wide algorithm identifies novel in-vivo expressed antigens inducing human T-cell responses with classical and unconventional cytokine profiles. Sci Rep. 2016;6:37792. doi:10.1038/srep37793.

    Article  Google Scholar 

  33. Becker K, Sander P. Mycobacterium tuberculosis lipoproteins in virulence and immunity—fighting with a double-edged sword. FEBS Lett. 2016. doi:10.1002/1873-3468.12273.

    Google Scholar 

  34. Paoletti AC, Parmely TJ, Tomomori-Sato C, Sato S, Zhu D, Conaway RC, Conaway JW, Florens L, Washburn MP. Quantitative proteomic analysis of distinct mammalian Mediator complexes using normalized spectral abundance factors. Proc Natl Acad Sci USA. 2006;103:18928–33.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  35. Baker ES, Liu T, Petyuk VA, Burnum-Johnson KE, Ibrahim YM, Anderson GA, Smith RD. Mass spectrometry for translational proteomics: progress and clinical implications. Genome Med. 2012;4:63. doi:10.1186/gm364.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  36. Schubert OT, Mouritsen J, Ludwig C, Röst HL, Rosenberger G, Arthur PK, Claassen M, Campbell DS, Sun Z, Farrah T, Gengenbacher M, Maiolica A, Kaufmann SHE, Moritz RL, Aebersold R. The Mtb proteome library: a resource of assays to quantify the complete proteome of Mycobacterium tuberculosis. Cell Host Microbe. 2013;13:602–12.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  37. Gengenbacher M, Mouritsen J, Schubert OT, Aebersold R, Kaufmann SH. Mycobacterium tuberculosis in the Proteomics Era. Microbiol Spectr. 2014. doi:10.1128/microbiolspec.MGM2-0020-2013.

    PubMed  Google Scholar 

  38. Peirs P, Parmentier B, De Wit L, Content J. The Mycobacterium bovis homologous protein of the Mycobacterium tuberculosis serine/threonine protein kinase Mbk (PknD) is truncated. FEMS Microbiol Lett. 2000;188:135139.

    Article  Google Scholar 

  39. Cho YS, Dobos KM, Prenni J, Yang H, Hess A, Rosenkrands I, Andersen P, Ryoo SW, Bai GH, Brennan MJ, Izzo A, Bielefeldt-Ohmann H, Belisle JT. Deciphering the proteome of the in vivo diagnostic reagent “purified protein derivative” from Mycobacterium tuberculosis. Proteomics. 2012;12:979–91.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  40. Choi HG, Choi S, Back YW, Park HS, Bae HS, Choi CH, Kim HJ. Mycobacterium tuberculosis Rv2882c protein induces activation of macrophages through TLR4 and exhibits vaccine potential. PLoS ONE. 2016;11(10):e0164458. doi:10.1371/journal.pone.0164458.

    Article  PubMed  PubMed Central  Google Scholar 

  41. Vita R, Overton JA, Greenbaum JA, Ponomarenko J, Clark JD, Cantrell R, Wheeler DK, Gabbard JL, Sette A, Peters B. The immune epitope database (IEDB) 3.0. Nucleic Acids Res. 2015;43:D405–12. doi:10.1093/nar/gku938.

    Article  PubMed  Google Scholar 

  42. Lindestam Arlehamn CL, Lewinsohn D, Sette A, Lewinsohn D. Antigens for CD4 and CD8 T cells in tuberculosis. Cold Spring Harb Perspect Med. 2014;4:a018465.

    Article  PubMed  PubMed Central  Google Scholar 

  43. Jones GJ, Whelan A, Clifford D, Coad M, Vordermeier HM. Improved skin test for differential diagnosis of bovine tuberculosis by the addition of Rv3020c-derived peptides. Vaccine Clin Immunol. 2012;19:620–2.

    Article  CAS  Google Scholar 

  44. Conlan AJ, Brooks Pollock E, McKinley TJ, Mitchell AP, Jones GJ, Vordermeier M, Wood JL. Potential benefits of cattle vaccination as a supplementary control for bovine tuberculosis. PLoS Comput Biol. 2015;11:e1004038. doi:10.1371/journal.pcbi.1004038.

    Article  PubMed  PubMed Central  Google Scholar 

  45. Baldwin CL, Telfer JC. The bovine model for elucidating the role of γδ T cells in controlling infectious diseases of importance to cattle and humans. Mol Immunol. 2015;66:35–47.

    Article  CAS  PubMed  Google Scholar 

  46. Waters WR, Maggioli MF, McGill JL, Lyashchenko KP, Palmer MV. Relevance of bovine tuberculosis research to the understanding of human disease: historical perspectives, approaches, and immunologic mechanisms. Vet Immunol Immunopathol. 2014;159:113–32.

    Article  CAS  PubMed  Google Scholar 

Download references

Authors’ contributions

SR coordinated the study and interpreted the results; RL and DMC helped to interpret the results; MC and VR participated in data analyses; MV prepared samples and performed proteomic analysis on LC/MS–MS; MG and GC supervised the laboratory work; SR and FR designed the study and wrote the paper. All authors read and approved the final manuscript.


We wish to thank Dr. Maria Saveria Campo, professor Emeritus of Glasgow University for reading the manuscript and for her useful comments.

Competing interests

The authors declare that they have no competing interests.

Availability of data and materials

All data produced for the purpose of this project are available for consultation and for the public, upon request to the corresponding author (SR) at

Author information

Authors and Affiliations


Corresponding author

Correspondence to Sante Roperto.

Additional files


Additional file 1: Table S1. This table shows: the list of the identified proteins within all four bPPD preparations, ordered locus names of M. bovis (Mb) and M. tuberculosis, strain H37Rv, (Rv), functional categories of identified proteins, their molecular mass, gene and protein lengths and gene name and protein function. A comparative protein expression profile with two previous bPPD proteomic analysis is also shown. Differences in functional category among identified proteins are in blue; differences in molecular mass, gene and protein lengths are in red. These findings were from TubercuList, BoviList, Uniprot and KEGG databases (See text for further details).


Additional file 2: Table S2. The list of identified proteins in bPPDs. Sequence coverage, number of unique peptides, as well as total weighted spectra as a semi-quantitative measure of protein abundance, are reported.


Additional file 3: Table S3. Dimethyl labeling-based quantification of proteins present in samples It, NLand Ip relative to sample S. Result of single replicates as well as average H:L and M:L ratios, their associated p-values and number of observations (unique peptides) are reported.


Additional file 4: Table S4. Proteins found differentially abundant in samples It, NL and Ip compared to sample S. The file contains three data sheets reporting 44, 34 and 30 proteins found at significantly different levels in samples It, NL and Ip, respectively. These proteins were mostly present at very low levels in all preparations. Fold changes >2 and <0.5 are reported in two separate lists.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, 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 ( applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Roperto, S., Varano, M., Russo, V. et al. Proteomic analysis of protein purified derivative of Mycobacterium bovis . J Transl Med 15, 68 (2017).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: