Skip to main content

Second generation multiple reaction monitoring assays for enhanced detection of ultra-low abundance Mycobacterium tuberculosis peptides in human serum

Abstract

Background

Mycobacterium tuberculosis (Mtb) is the causative agent of Tuberculosis (TB), the number one cause of death due to an infectious disease. TB diagnosis is performed by microscopy, culture or PCR amplification of bacterial DNA, all of which require patient sputum or the biopsy of infected tissue. Detection of mycobacterial products in serum, as biomarkers of diagnosis or disease status would provide an improvement over current methods. Due to the low-abundanceĀ of mycobacterial products in serum, we have explored exosome enrichment to improve sensitivity. Mtb resides intracellularly where its secreted proteins have been shown to be packaged into host exosomes and released into the bloodstream. Exosomes can be readily purified assuring an enrichment of mycobacterial analytes from the complex mix of host serum proteins.

Methods

Multiple reaction monitoring assays were optimized for the enhanced detection of 41 Mtb peptides in exosomes purified from the serum of individuals with TB. Exosomes isolated from the serum of healthy individuals was used to create and validate a unique data analysis algorithm and identify filters to reduce the rate of false positives, attributed to host m/z interference. The final optimized method was tested in 40 exosome samples from TB positive patients.

Results

Our enhanced methods provide limit of detection and quantification averaging in the low femtomolar range for detection of mycobacterial products in serum. At least one mycobacterial peptide was identified in 92.5% of the TB positive patients. Four peptides from the Mtb proteins, Cfp2, Mpt32, Mpt64 and BfrB, show normalized total peak areas significantly higher in individuals with active TB as compared to healthy controls; three of the peptides from these proteins have not previously been associated with serum exosomes from individuals with active TB disease. Some of the detected peptides were significantly associated with specific geographical locations, highlighting potential markers that can be linked to the Mtb strains circulating within each given region.

Conclusions

An enhanced MRM method to detect ultra-low abundance Mtb peptides in human serum exosomes is demonstrated, highlighting the potential of this methodology for TB diagnostic biomarker development.

Background

The World Health Organization estimates that 2 billion people globally are infected with Mycobacterium tuberculosis (Mtb) [1]. One of the biggest hurdles for the global control of tuberculosis (TB) is the lack of point-of-care, rapid and accurate diagnostic tools. Current diagnostics rely on sputum availability to confirm the presence of bacteria by microscopy, live bacilli by culture or using molecular tools to detect pathogen DNA. However, microscopy only detects 20ā€“80% of all active TB cases [2] and although diagnosis by sputum culture is highly sensitive, it can take 4ā€“6Ā weeks to yield results; this results in delays in treatment and continued transmission of disease in the community. To reduce rates of TB worldwide, simpler more easily scalable accurate diagnostics are needed, prompting research on novel biomarker discovery and new diagnostic assays, with an emphasis on alternative non-sputum diagnostics targeting blood, urine, and breath.

We previously demonstrated that during cellular infection with Mtb, mycobacterial products can be incorporated into host cell exosomes [3]; the presence of Mtb proteins were confirmed in exosomes isolated from a variety of fluids from infected animal models [4, 5]. Exosomes are 100Ā nm vesicles generated by all nucleated cell types and released into biofluids, including blood, urine and sputum [6,7,8,9,10]. Purification of exosomes from serum is a facile means to reduce the complexity of the fluid and concentrate the encapsulated Mtb proteins. Despite the simplification, the exosomes of interest, that is those derived from the Mtb-infected cells, are still a minor component of the entire exosome population. Therefore coupling this purification with a sensitive downstream detection platform is critical for the detection of mycobacterial proteins as potential biomarker candidates.

Using targeted mass spectrometry we previously conducted a pilot study to determine if 33 mycobacterial proteins that we previously identified in cell culture [3] and animal studies were also present in exosomes purified from the serum of individuals diagnosed with active disease or known to have latent TB infection [11]. We found at least one of the 76 peptides from 27 of the 33 Mtb proteins in at least a single individual of a cohort or 57 subjects; the remaining 6 Mtb proteins were not identified in any of the 57 samples. This discovery experiment provided us with a preliminary list of Mtb protein candidate biomarkers and a rapid method for triaging peptide candidates to proceed to assay refinement and larger verification studies.

The main aim of this study was to enhance an MRM method to detect Mtb peptides in serum exosomes from TB patients. The aforementioned study by our group used 17 unrefined targeted MRM assays to screen for the presence of peptides from 33 Mtb proteins [11]. The goal of this study is to validate these findings using refined MRM methods to reduce the number of channels being monitored and improve sensitivity and specificity by spiking isotopically labeled peptide standards and including control samples from healthy donors. Novel data analysis methodology was applied to confirm the presence of 19 low-abundance mycobacterial proteins in samples dominated by host proteins. In addition, in this study we have analyzed samples from both HIV+ and HIVāˆ’Ā TB patients. Individuals with HIV/TB co-infection are often difficult to diagnose using current techniques due to paucibacillary presentation of the disease which can compromise detection by microscopy, and while detection by culture has higher positivity rates, the collection of several samples is often needed to achieve a conclusive diagnosis in these patients [12, 13]. HIV+ infected patients also have a lower rate of positivity using the TB skin test or Interferon Gamma Release Assays due to their inherent immunocompromised status [13,14,15]. Thus, a diagnostic method that can efficiently detect both HIV+ and HIVāˆ’Ā TB disease and/or TB infection at similar rates is needed. The ultimate goal is to define biomarkers of active disease (regardless of HIV status), discovered and verified by MRM-MS, and to translate the promising candidates to a point-of-care platform for use in the field that is independent of complex sample processing and high-end instrumentation.

Methods

Study design

Serum samples from TB negative (suspect) and TB positive patients were obtained from the Foundation for Innovative Novel Diagnostics (FIND) specimen repository (Geneva, Switzerland). Serum samples from healthy donors, with no history of tuberculosis were purchased from Bioreclamation IVT (Westbury, NY). All samples were stored at āˆ’80Ā Ā°C upon arrival until processed. For initial method development and optimization we used 16Ā TB negative and 20Ā TB positive serum samples from 4 different geographic locations (TableĀ 1); these samples were pooled and used as background matrix. For MRM method validation we used 20 individual healthy controls and 40 individual TB positive serum samples all with culture confirmed pulmonary tuberculosis (TableĀ 2).

TableĀ 1 Patient breakdown of samples included in the pooled matrix sample
TableĀ 2 Patient breakdown of samples included in the assay verification set

Sample processing

Serum samples (250Ā ĀµL) were centrifuged to remove whole cells/cellular debris. Exoquick (System Biosciences, Palo Alto, CA) was added to the cleaned serum at a 4:1 ratio (sample:reagent), incubated at 4Ā Ā°C for 30Ā min, and exosomes were pelleted by centrifugation at 1.5Ā k x g for 30Ā min, as per manufacturer recommendation. The exosomes were suspended in 250 ĀµL of PBS and micro bicinchoninic acid assay was performed to quantify protein content. To construct the pooled matrix stock, 50Ā Āµg (protein) of purified exosomes from the 36 samples listed in TableĀ 1, were mixed. 20Ā Āµg of pooled matrix or 50Ā Āµg (1.18Ā Ā±Ā 0.41Ā ĀµL) of individual exosome sample were run into a NuPAGE Novex 4-12% Bisā€“Tris Gel 1.0Ā mm gel for 5Ā min in NuPAGE MES SDS Running Buffer (Life Technologies, Carlsbad, CA) to trap the Exoquick polymer prior to peptide extraction. In-gel digest with sequencing-grade trypsin (Roche, Switzerland) was performed at a 1:20 (enzyme:substrate) ratio, overnight at 37Ā Ā°C, as previously described [3]. The extracted peptides were dried and suspended at a final concentration of 1Ā Āµg/ĀµL in 3% acetonitrile (ACN), 0.5% formic acid (FA) in water.

Peptide standards

Isotope-labelled standards for each of our peptides of interest were purchased from New England Peptide (Gardner, MA). All peptides were prepared to a minimum of 95% purity and confirmed by HPLC. QC of chromatography was monitored by the addition of indexed retention time (iRT) standards (Biognosys AG, Switzerland) [16]. 10Ā nM mixes of the isotope-labelled peptide standards and iRT mix at a 0.4Ɨ final concentration were spiked into each sample.

Multiple reaction monitoring

Daily Skyline (64-bit) was used to build and optimize the multiple reaction monitoring (MRM) methods for the relative quantification of peptides [17]. For this study, two MRM methods were built in Skyline. The first (MRM-1), included 20 peptides from 9 Mtb proteins (TableĀ 3). The second method (MRM-2), included 10 Mtb proteins and 21 peptides (TableĀ 3). Briefly, FASTA-formatted sequences of all 19 proteins were used for in silico tryptic (KR|P) digestion with peptides being selected based on previous discovery studies [4, 18]. Both double and triple charge precursor ions were empirically tested and selected based on their performance. ā€œyā€ ions for each transition were selected based on a library built from LCā€“MS/MS data acquired in an Orbitrap Velos (Thermo Scientific) (Additional file 1). The resultant methods were exported to Masslynx (Waters Corporation, Milford, MA). All method development was performed using a 10Ā nM mix of all 41 heavy labeled peptides (K^Ā =Ā Lysine, 13C6, 15N2 or R^Ā =Ā Arginine, 13C6, 15N4) (New England Peptide, Gardner, MA) in a 1Ā Āµg/ĀµL matrix background (as described above). One and a half microliters (1.5Ā ĀµL) of the heavy labeled peptide mix in background matrix were then injected into the LCā€“MS/MS system consisting of a Waters nanoACQUITY UPLC coupled to a Waters Xevo TQ-S mass spectrometer fitted with a Trizaic source. The instrument was operated with MassLynx V4.1 SCN810 (Waters Corporation, Milford, MA). Chromatography was performed on a 150Ā ĀµmĀ Ć—Ā 50Ā mm Ion key packed with BEH C18 130Ā Ć…, 1.7Ā um. The chromatography length and gradient were optimized to separate all peptides as much as possible so that at least 12 points per each transition were acquired. Peptides were separated using gradient elution with a stable flow rate of 3.06 ĀµL/min. A linear method consisting of 2Ā min of equilibration in 97% buffer A (99.9% water with 0.1% formic acid) and 3% buffer B (99.9% ACN with 0.1% formic acid), followed by a 45Ā min linear gradient to 22% buffer B. The method finished with 5Ā min wash at 97% buffer B and final equilibration at 3% buffer B for 5Ā min. The column was maintained at 45Ā Ā°C during analysis, and the samples were kept at 4Ā Ā°C at all times. The mass spectrometer was operated in selective reaction monitoring mode using electrospray ionization in nanospray positive ion mode, with a capillary voltage of 3.6Ā kV and a source temperature of 100Ā Ā°C. Cone voltage was static at 35Ā V and the collision energies were in silico predicted by Skyline for each compound individually (Additional file 1). The final methods included the selection of the 5 most abundant transitions per peptide and the precursor ion with the best peak shape and overall signal.

TableĀ 3 Peptides/proteins included in MRM assays 1 and 2

MRM of clinical samples

Trypsin digested clinical samples were resuspended at a final concentration of 1Ā Āµg/ĀµL in 5% Acetonitrile, 0.1% Formic acid containing 10Ā nM internal standard mix. After resuspension each sample was centrifuged 10Ā min to pellet minor impurities. Supernatant was then transferred to a MS vial and placed in the Xevo-TQS autosampler. One and a half microliters (1.5Ā ĀµL) were injected into the instrument and data was acquired as described above monitoring for both heavy and light forms for a total of 10 transitions per peptide. Blank runs were run every 4 samples. After MRM analysis of each clinical sample, raw files were imported into Skyline where prior to data processing in Excel (see section below), each sample was manually validated for quality (i.e. retention time; peak shape; and intensity and peak boundaries of internal standards). The majority of samples were run only once. However in cases were the sample did not pass the manual QC (usually due to retention time drift), that sample was re-injected at least once. In some cases, even after multiple injections, data did not pass manual QC. In those case that particular sample/peptide was removed from the results.

Determination of limit of detection (LOD) and limit of quantification (LOQ)

Serial dilutions (10ā€“0.01Ā nM) of the heavy labeled peptides were prepared in 1Ā Āµg/ĀµL of pooled matrix. 1.5Ā ĀµL of each dilution were injected into the LCā€“MS/MS system described above using a 2 or 4Ā min window for acquisition of precursor and transition of each peptide. All dilution points were run in triplicate. Raw data resulting from each of the dilution points were exported into Skyline. Peaks were manually validated and boundaries adjusted if necessary. Only the final set of transitions after validation with clinical samples was used for LOD/LOQ purposes. Area under the peak for each of the validated transitions was exported into excel and the LOD/LOQ was calculated using the LINEST method. Briefly, a linear regression was performed and the LOD and LOQ were calculated as shown in Eqs.Ā 1 and 2 respectively, where m is the slope of the curve and s(y) is the standard deviation of the y values.

$${\text{LOD}} = 3*\frac{m}{s\left( y \right)}$$
(1)
$${\text{LOQ}} = 10*\frac{m}{s\left( y \right)}$$
(2)

All LOD and LOQ values are expressed in fmol/ĀµL (Additional file 2).

Liquid chromatography-multiple reaction monitoring mass spectrometry analysis

Data processing was performed using Skyline software. Manual inspection and border adjustment of isotopic standard peaks was performed for each peptide. While five transitions were monitored for each peptide, the final analysis used 3ā€“5 transitions selected based on low background/noise level. Any standard peptide which showed either an aberrant transition ratio or retention time (RT) was used to disqualify the inclusion of the native peptide from sample analysis.

Data post-processing

After initial screening in Skyline, peak areas for individual transitions and RTs for both native and standard peaks were exported to Excel. In total 19 methods for processing the raw data were trialed (data not shown). The final method, summarized in Fig.Ā 1, utilized the sum of the transition peak areas or total peak area (TPA) for each native peptide, followed by normalization by the TPA of the isotope-labelled standard (peak area ratio). In order for a peak to be qualified for inclusion, three minimal qualifications were set: (1) the TPA of the native peptide must be composed of values from 3 or more transitions, (2) the RT of native transition peaks must be within 0.1Ā min of the standard transition peak, and (3) the final normalized TPA must exceed the cut-off determined by the healthy donor samples. The post data-processing protocol, including the first two qualifications for inclusion stated above, were applied to the raw values for the healthy data to formulate the normalized TPAs. The healthy cut-off value was established as the mean normalized TPA plus three times the standard deviation of all 20 samples and are summarized in Additional file 3. False discovery rates, calculated as the percent of healthy samples with nTPA values above the cut-off is included in Additional file 3; the maximal number of controls with a positive signal in any sample is 1 (5% FDR). The final Skyline data files have been deposited to the Panorama Repository (https://panoramaweb.org/) [19]. Unpaired, two-tailed t-tests were performed using GraphPad Prism version 6.00 for Windows (La Jolla, CA). Benjamin-Hochberg adjustment was used to control the false discovery rate.

Fig.Ā 1
figure 1

Raw data processing workflow. TPA total peak area (sum of transition peaks), nTPA normalized TPA (ratio native/labeled standard), TPA healthy cut-off determined by mean +3Ɨ SD

Results

LOD/LOQ determination in a pooled exosome matrix

Exosomes isolated from the serum of 16Ā TB suspects and 20Ā TB culture-confirmed patients were pooled to create a matrix for the determination of the LOD and LOQ for the 41 peptides in our two MRM assays. Most peptides displayed sub-nanomolar (nM) LODs (Additional file 2). The rationale for creating a pooled matrix for the purposes of establishing LOD and LOQ was to maximize the likelihood of selecting the most intense transitions for query as we have previously identified several transitions which are highly influenced by the matrix (data not shown), and also to reduce bias created by a single representative sample. Individual transitions were subject to elimination from final analyses if LOD was greater than 1Ā nM, with a maximum removal of two transitions per peptide.

Identification and geographic diversity of Mtb peptides in serum exosomes

Out of 40 patients with active TB, we were able to identify 35 (87.5%) of them by the identification of at least one Mtb peptide monitored in MRM-1 in their serum exosomes. MRM-2 was less successful, allowing for the identification of only 23 (57.5%) of the TB positive patients. Overall, if both assay 1 and 2 results are combined, we see at least one peptide present in 37 out of the 40 active TB subjects (92.5%). The gain of 2 new positive IDs were based on peptide TAVEQAAAELGDTGR from the Mtb protein MrsA (Rv3441c) from MRM-2. Three patients out of the 40 (7.5%) did not have any bacterial peptide identified by either MRM assay, and they varied by smear status, HIV status, and region of origin (Fig.Ā 2). Sixteen out of 17 (94%) patients missed by sputum smear microscopy were identified in our MRM assays by at least one peptide (Fig.Ā 2; Additional file 4). There was no statistical difference (pĀ =Ā 0.14) between the numbers of peptides identified in HIV positive patients (3.2Ā Ā±Ā 1.7) when compared to those without co-infection (2.4Ā Ā±Ā 1.6). When stratified by geographical location, patients from South Africa displayed the highest numbers of peptides identified per patient (3.6Ā Ā±Ā 1.1), followed by Vietnam (2.9Ā Ā±Ā 1.5) and lastly, Peru and Bangladesh (2.1Ā Ā±Ā 1.9 and 2.1Ā Ā±Ā 1.8, respectively) (Additional file 5). The difference between the South African subgroup was significant when compared to the Peru (pĀ =Ā 0.040) and Bangladesh (pĀ =Ā 0.036) subgroups.

Fig.Ā 2
figure 2

Stratification of patient samples by HIV status, sputum smear microscopy status, and geography. Each column indicates a single patient (nĀ =Ā 40). Black circles indicate the presence of each protein (rows) determined by a normalized TPA of one or more peptides above the healthy threshold. Geography is indicated by: B Bangladesh, P Peru, S South Africa, V Vietnam

The frequency by which each peptide/protein was identified in our cohort of 40 samples is summarized in TableĀ 4. The top ranking peptide for assay 1 was: GSLVEGGIGGTEAR (from the Mtb protein Cfp2); this peptide was identified in 24 out of 40 (60%) samples. Interestingly, this peptide was identified in 90% of all of the South African samples, but in only 50% of the other three locations (Fig.Ā 2). The difference between the normalized TPA in TB patients and healthy individuals is statistically different by t test for the GSLVEGGIGGTEAR peptide with a p valueĀ <Ā 0.001 (Fig.Ā 3). In addition, the peptides LYASAEATDSK, FLSAATSSTPR, and EALALALDQER from Mtb proteins Mpt32, Mpt64, and BrfB, respectively, showed a significant (pĀ <Ā 0.05) difference in signal between those with active disease and the healthy cohort (Fig.Ā 3/Additional file 6). When these 4 significant peptides were stratified by the HIV status of each patient, the normalized TPA values between the HIV positive and HIV negative groups were not significant. Similarly, when stratified by smear status, the differences in signal were not statistically different between groups. One or more of these 4 peptides was represented in 70% of the samples screened. In addition to these four peptides we also identified 16 additional peptides that were detected in at least 1 or up to 7 individual active TB patients and thus also represent candidates for TB disease biomarkers (Additional file 7).

TableĀ 4 Percentage of samples positively identified with each Mtb peptide (left) and protein (right)
Fig.Ā 3
figure 3

Comparison of normalized total peak areas from the TB patients and healthy individuals. The solid lines indicate the mean nTPA with SD. The dotted line represents the threshold created by the average normalized value of the healthy plus 3 times the standard deviation; TB patients were designated as positive for a peptide if the TPA was above this threshold. Of the 41 peptides screened, four were statistically more abundant in those with active disease (p values detailed in Additional file 6)

Eight peptides in MRM-1 and 9 peptides in MRM-2 showed no discrimination between healthy individuals and those with active tuberculosis (Additional file 8, Panel 1 and 2, respectively). This includes: VQFQGGGPHAVYLLDGLR (Ag85c), FLEGLTLR (Ag85c), AGANLFELENFVAR (BrfB), SLADPNVSFANK (Cfp2), VVADLTPQNQALLNAR (GlcB), TVSLPVGADEDDIK (HspX), SLENYIAQTR (Mpt64), and GVTEETTTGVLR (SahH) from MRM-1 and TVGDVVAYIQK (AcpM), AADMWGPSSDPAWER (Ag85b), ITQDLLDR (DnaK), LAAAWGGSGSEAYQGVQQK (Esat-6), FLLDQAITSAGR (GarA), GGYFPVAPNDQYVDLR (GlnA1), YVLEELR (MrsA), VIQGFMIQGGDPTGTGR (PpiA), and HTIFGEVIDAESQR (PpiA) from MRM-2. As none of these peptides showed a signal above the background level determined by the healthy serum exosomes, these 17 peptides will be excluded as markers of active tuberculosis in future studies. Several of the eliminated peptides correlate with poor LOD values. Since multiple peptides are targeted per protein in our assays, the exclusion of these 17 peptides did not eliminate any of the Mtb proteins in our original candidate list. Lastly, three peptides could not be fully evaluated due to the contribution of the matrix effect, as a baseline value in the healthy patient cohort could not be determined (Additional file 9); the suppression of the 3 standard peptides in all 20 healthy exosome samples supports our concern of matrix variation and the use of a pooled matrix for the LOD calculations.

Discussion

The burden of mycobacterial proteins circulating in human serum during active TB infection is currently unknown, but based on prior proteomic and molecular-based investigations, it is believed to be of low abundance [20]. In this study, we have shown that intracellular mycobacteria contribute to the protein content of exosomes [3, 18] and purification of these microvesicles from the serum can serve as an effective method of depletion of abundant host proteins (e.g. albumin and immunoglobulins) while concentrating the bacterial products of interest. Recent studies have shown that proteins from other mycobacterial species, such as M. avium [21], as well as other intracellular pathogens, such as Helicobacter pylori [22] can be readily detected in exosomes. In the current study, we sought to confirm and build upon our previous discovery studies in which Mtb peptides for 33 unique Mtb proteins were detected in human serum or patients with active and/or latent TB [11]. Our previous study utilized 17 unrefined targeted MRM assays, here we optimized the assays to include all selected peptides in two different methods. In this study, we included isotopically labeled standards for the first time, to confirm peptide identification based on retention time and fragmentation pattern while also providing a normalizing tool allowing direct comparison of peak areas of a given peptide from sample to sample. Ideally, we would have included all of the significant peptides determined in our original publication [11], however this was not feasible, due to poor performance of several of the isotope-labeled peptide standards during the development of the MRM assays. Specifically, the Mpt64 protein, dictated the selection of alternative peptides; of the two replacement Mpt64 peptides, one remained statistically significant (Fig.Ā 3), while the second was nondiscriminatory (Additional file 8). Both the addition of the isotope-labeled standards, as well as iRT, retention time standard peptides, provide extra QC tools that allow us to monitor run-to-run variance. Even though RT variability was well within the 2Ā min monitoring window (Additional file 10), the inclusion of iRTs and labeled standards allowed us to monitor run-to-run variation and to determine the need for re-injection in a few cases were RT shifted across the board due potentially to matrix effects. Monitoring of RT variability also allowed us to see a pattern in which the late eluting peptides seemed to have a more stable RT. Although we did not follow up on the potential reasons for this trend, it may be an important factor to consider during development of future MRMs.

In this study, we also optimized data analysis to increase the stringency of which a given peptide signal is considered positive in any given sample. Matching retention times to those of the labeled standard, as well as the requirement to only count peptides with at least three detected transitions were important parameters in our algorithm. However, for many peptides, the most restrictive filter was the inclusion of the threshold determined by the normalized TPA in the healthy donor samples which represents high noise due to the complexity of the exosome matrix and which could be mistaken for a positive result. Overall, the optimized MRM method including the addition of labeled standards and iRTS and the optimized data analysis algorithm, allowed us to reduce our original candidate biomarker list from 41 to 20 highly significant peptides in the context of active TB disease.

As noted in our method section, the MRM methods presented in this study were optimized with a window of 2 or 4Ā min for scanning of each peptide. Retention times were optimized using the labeled peptides spiked into a mixture of serum exosomes (our pooled matrix). While for the majority of peptides, the retention time was stable within that window from sample to sample, there are limitations to this type of scheduling, and a few peptides, such as FALNAANAR, AGANLFELENFVAR, VQFQGGGPHAVYLLDGLR and PGLPVEYLQVPSPSMGR, were not routinely detected within our sample window as a result. This is likely because the retention time for these peptides was very dependent on the matrix background and often resulted in either increased or decreased retention time outside of the scanning window. These peptides were also very susceptible to sample-to-sample variation in the total peak area measured for the labeled standard. While scheduling the scanning window for each peptide allowed us to include more peptides in a single method without a reduction in sensitivity, some peptides, such as the ones mentioned above, may benefit from an unscheduled method in which all channels (one for each peptide) are continuously scanned throughout the length of the chromatography. The variation observed here demonstrates the effect that matrix background can have in peptide detection, which was further evidenced by the ion suppression we observed for three labeled peptide standards (TTGDPPFPGQPPPVANDTR, PGLPVEYLQVPSPSMGR, and SVFDDGLAFDGSSIR) in the healthy controls.

Importantly, the variation in retention time and ion suppression described above demonstrates the importance of having isotopically labeled peptides spiked in each sample such that the correct peak is identified (i.e. in cases of retention time drift) and to correctly normalize the area under the peak of the native peptide (i.e. ion suppression dependent on matrix background). The inclusion of isotopically labeled peptides is thus a major improvement from our previous study, strengthening the data set by significantly increasing the confidence for the detection of bona-fide mycobacterial peptides. For future studies, we believe that the methods can be improved even further by optimizing the method without any scheduling and by increasing the number of healthy controls so that more accurate mean TPA and more tolerant cut-offs can be applied as outliers appear to have a significant impact on the threshold.

Blood-based assays are highly desirable for TB diagnostics, as they are independent of an individualā€™s ability to produce sputum (i.e. children), as well as remove the occupational hazard of transmission during sputum collection. Our study emphasize the benefit of blood based biomarkers by highlighting the prevalence of Mtb proteins in serum exosomes in those failed to be identified by the first-line diagnostic screening by sputum microscopy. Downstream translation of our assay to a point-of-care platform is in line with the end goal of creating an assay that detects active infection faster than sputum culture. Of equal importance, the identification of mycobacterial proteins in serum exosomes is not affected by HIV co-infection; serum biomarkers may have an added benefit in detecting disseminated disease often displayed in HIV positive patients.

Inclusion of samples from a variety of geographic locations was another strength of this study, as it provided us with a diverse set of samples from persons with disease caused by different Mtb strains [23]. Our original work focused on a clinical cohort from Uganda [11], at this time we did not fully appreciate the potential effect of geography on our results which impacted the selection of protein candidates for this current study. Thus, it is possible that we may have missed proteins or peptides that may perform better in other geographical locations, including those analyzed in the present study. As the mycobacterial contribution to the exosome is reflective of the bacterial secretome, we sought to identify biomarkers common to Mtb complex, independent of lineage or strain. However, we found that Cfp2 peptides performed best in the South Africa cohort (90%) and Mpt32 peptides underperformed as a biomarker in Bangladesh (40%, compared to 60ā€“70% in other regions). Cfp2 is a secreted protein for which differences in abundance and secretion among Mtb strains had been previously reported by our group [24], supporting the idea that some of the differences we observed among geographical locations may be associated with specific circulating Mtb strains. However, other factors such as those related to the host (i.e. matrix effects favoring ionization of specific peptides), as well as the heterogeneity of the infection, may also account for the differences seen here, but additional studies are needed.

Conclusions

Through improved sample preparation, enhanced MRM assays and optimized data analysis work flow, in this study, we were able to find the proverbial Mtb peptide ā€œneedleā€ in a serum ā€œhaystackā€. Whereas the methodology as it stands is not suitable for the purpose of TB diagnosis in the field setting, it provides confirmation of the presence of mycobacterial protein candidates in serum exosomes. Additional verification and validation studies with larger and more geographically diverse sample sets, as well as translation to a point-of-care assay are future goals.

Abbreviations

MRM:

multiple reaction monitoring

MS:

mass spectrometry

TB:

tuberculosis

Mtb:

Mycobacterium tuberculosis

FIND:

Foundation for Innovative novel Diagnostics

ACN:

acetonitrile

FA:

formic acid

iRT:

indexed retention time

TPA:

total peak area

QC:

quality control

LOD:

limit of detection

LOQ:

limit of quantification

FDR:

false discovery rate

References

  1. Tuberculosis Fact Sheet [http://www.who.int/mediacentre/factsheets/fs104/en/].

  2. Steingart KR, Ng V, Henry M, Hopewell PC, Ramsay A, Cunningham J, Urbanczik R, Perkins MD, Aziz MA, Pai M. Sputum processing methods to improve the sensitivity of smear microscopy for tuberculosis: a systematic review. Lancet Infect Dis. 2006;6(10):664ā€“74.

    ArticleĀ  PubMedĀ  Google ScholarĀ 

  3. Giri PK, Kruh NA, Dobos KM, Schorey JS. Proteomic analysis identifies highly antigenic proteins in exosomes from M. tuberculosis-infected and culture filtrate protein-treated macrophages. Proteomics. 2010;10(17):3190ā€“202.

    ArticleĀ  CASĀ  PubMedĀ  PubMed CentralĀ  Google ScholarĀ 

  4. Kruh-Garcia NA, Wolfe LM, Dobos KM. Deciphering the role of exosomes in tuberculosis. Tuberculosis (Edinb). 2015;95(1):26ā€“30.

    ArticleĀ  CASĀ  Google ScholarĀ 

  5. Kruh-Garcia NA, Schorey SJ, Dobos KM. Exosomes: new tuberculosis biomarkersā€”prospects from the bench to the clinic. In: Cardona P-J, editor. Understanding tuberculosis: global experiences and innovative approaches to the diagnosis. Paris: InTech; 2012.

    Google ScholarĀ 

  6. Admyre C, Telemo E, Almqvist N, Lotvall J, Lahesmaa R, Scheynius A, Gabrielsson S. Exosomesā€”nanovesicles with possible roles in allergic inflammation. Allergy. 2008;63(4):404ā€“8.

    ArticleĀ  CASĀ  PubMedĀ  Google ScholarĀ 

  7. Hessels D, Schalken JA. Urinary biomarkers for prostate cancer: a review. Asian J Androl. 2013;15(3):333ā€“9.

    ArticleĀ  CASĀ  PubMedĀ  PubMed CentralĀ  Google ScholarĀ 

  8. Hoorn EJ, Pisitkun T, Zietse R, Gross P, Frokiaer J, Wang NS, Gonzales PA, Star RA, Knepper MA. Prospects for urinary proteomics: exosomes as a source of urinary biomarkers. Nephrology (Carlton). 2005;10(3):283ā€“90.

    ArticleĀ  CASĀ  Google ScholarĀ 

  9. Mathivanan S, Ji H, Simpson RJ. Exosomes: extracellular organelles important in intercellular communication. J Proteomics. 2010;73(10):1907ā€“20.

    ArticleĀ  CASĀ  PubMedĀ  Google ScholarĀ 

  10. Simpson RJ, Jensen SS, Lim JW. Proteomic profiling of exosomes: current perspectives. Proteomics. 2008;8(19):4083ā€“99.

    ArticleĀ  CASĀ  PubMedĀ  Google ScholarĀ 

  11. Kruh-Garcia NA, Wolfe LM, Chaisson LH, Worodria WO, Nahid P, Schorey JS, Davis JL, Dobos KM. Detection of Mycobacterium tuberculosis peptides in the exosomes of patients with active and latent M. tuberculosis infection using MRM-MS. PLoS ONE. 2014;9(7):e103811.

    ArticleĀ  PubMedĀ  PubMed CentralĀ  Google ScholarĀ 

  12. Monkongdee P, McCarthy KD, Cain KP, Tasaneeyapan T, Nguyen HD, Nguyen TN, Nguyen TB, Teeratakulpisarn N, Udomsantisuk N, Heilig C, et al. Yield of acid-fast smear and mycobacterial culture for tuberculosis diagnosis in people with human immunodeficiency virus. Am J Respir Crit Care Med. 2009;180(9):903ā€“8.

    ArticleĀ  PubMedĀ  Google ScholarĀ 

  13. Tornheim JA, Dooley KE: Tuberculosis associated with HIV infection. Microbiol Spectr 2017, 5(1).

  14. Fisk TL, Hon HM, Lennox JL, Fordham von Reyn C, Horsburgh CR Jr. Detection of latent tuberculosis among HIV-infected patients after initiation of highly active antiretroviral therapy. AIDS. 2003;17(7):1102ā€“4.

    ArticleĀ  PubMedĀ  Google ScholarĀ 

  15. Girardi E, Palmieri F, Zaccarelli M, Tozzi V, Trotta MP, Selva C, Narciso P, Petrosillo N, Antinori A, Ippolito G. High incidence of tuberculin skin test conversion among HIV-infected individuals who have a favourable immunological response to highly active antiretroviral therapy. AIDS. 2002;16(14):1976ā€“9.

    ArticleĀ  PubMedĀ  Google ScholarĀ 

  16. Escher C, Reiter L, MacLean B, Ossola R, Herzog F, Chilton J, MacCoss MJ, Rinner O. Using iRT, a normalized retention time for more targeted measurement of peptides. Proteomics. 2012;12(8):1111ā€“21.

    ArticleĀ  CASĀ  PubMedĀ  PubMed CentralĀ  Google ScholarĀ 

  17. MacLean B, Tomazela DM, Shulman N, Chambers M, Finney GL, Frewen B, Kern R, Tabb DL, Liebler DC, MacCoss MJ. Skyline: an open source document editor for creating and analyzing targeted proteomics experiments. Bioinformatics. 2010;26(7):966ā€“8.

    ArticleĀ  CASĀ  PubMedĀ  PubMed CentralĀ  Google ScholarĀ 

  18. Kruh-Garcia NA, Murray M, Prucha JG, Dobos KM. Antigen 85 variation across lineages of Mycobacterium tuberculosis-implications for vaccine and biomarker success. J Proteom. 2014;97:141ā€“50.

    ArticleĀ  CASĀ  Google ScholarĀ 

  19. Sharma V, Eckels J, Taylor GK, Shulman NJ, Stergachis AB, Joyner SA, Yan P, Whiteaker JR, Halusa GN, Schilling B, et al. Panorama: a targeted proteomics knowledge base. J Proteom Res. 2014;13(9):4205ā€“10.

    ArticleĀ  CASĀ  Google ScholarĀ 

  20. Liu C, Zhao Z, Fan J, Lyon CJ, Wu H-J, Nedelkov D, Zelazny AM, Olivier KN, Cazares LH, Holland SM, et al. Quantification of circulating Mycobacterium tuberculosis antigen peptides allows rapid diagnosis of active disease and treatment monitoring. Proc Natl Acad Sci. 2017;114(15):3969ā€“74.

    ArticleĀ  CASĀ  PubMedĀ  PubMed CentralĀ  Google ScholarĀ 

  21. Wang JJ, Chen C, Xie PF, Pan Y, Tan YH, Tang LJ. Proteomic analysis and immune properties of exosomes released by macrophages infected with Mycobacterium avium. Microbes Infect. 2014;16(4):283ā€“91.

    ArticleĀ  CASĀ  PubMedĀ  Google ScholarĀ 

  22. Shimoda A, Ueda K, Nishiumi S, Murata-Kamiya N, Mukai SA, Sawada S, Azuma T, Hatakeyama M, Akiyoshi K. Exosomes as nanocarriers for systemic delivery of the Helicobacter pylori virulence factor CagA. Sci Rep. 2016;6:18346.

    ArticleĀ  CASĀ  PubMedĀ  PubMed CentralĀ  Google ScholarĀ 

  23. Gagneux S, Small PM. Global phylogeography of Mycobacterium tuberculosis and implications for tuberculosis product development. Lancet Infect Dis. 2007;7(5):328ā€“37.

    ArticleĀ  PubMedĀ  Google ScholarĀ 

  24. Mehaffy C, Hess A, Prenni JE, Mathema B, Kreiswirth B, Dobos KM. Descriptive proteomic analysis shows protein variability between closely related clinical isolates of Mycobacterium tuberculosis. Proteomics. 2010;10(10):1966ā€“84.

    ArticleĀ  CASĀ  PubMedĀ  PubMed CentralĀ  Google ScholarĀ 

Download references

Authorā€™s contributions

CM: MRM method development, data analysis, manuscript writing. KMD: Study design, data analysis, manuscript writing. PN: Results interpretation and manuscript writing. NAK-G: Study design, data analysis, manuscript writing. All authors read and approved the final manuscript.

Acknowledgements

We thank the Foundation for Innovative Novel Diagnostics for provision of samples, Lisa Wolfe and Jessica Prenni from the Proteomics and Metabolomics Facility (CSU) for guidance with assay development and instrumentation, Kaitlyn Nagel (CSU) for sample processing, Eric Vittinghoff (UCSF) for assistance with data clean-up method development, and Ann Hess (CSU) for assistance with the statistical analysis. Finally, we would like to acknowledege support from the Colorado State University Libraries Open Access Research and Scholarship Fund.

Competing interests

The authors declare that they have no competing interests.

Availability of data and materials

The final Skyline data files have been deposited to the Panorama Repository (https://panoramaweb.org/data/labkey/files/CSU-Kruh-GarciaLab) [19].

Declarations

This work was determined to be not human subjectā€™s research upon initial review from Colorado State Universityā€™s Institutional Review Board, as all human samples were obtained as de-identified samples from biorepositories. Specific to the tuberculosis patient serum samples; all were acquired from the not-for-profit foundation, Foundation for Innovative New Diagnostics (FIND; Geneva, Switzerland). Subjects cannot be identified, directly or through identifiers linked to the subjects. As stated on their webpage ā€œAll specimens are collected at participating qualified clinics, under a protocol that is approved by an Institutional Review Board (IRB) and under Informed Consent (IC) from patients meeting the inclusion criteria for the collection of samples as indicated in the study protocol.ā€ http://www.finddx.org/specimen-banks/. Similarly, Bioreclamation IVT samples were collected under IRB approved protocols.

Funding

This research was supported by funding from The Bill and Melinda Gates Foundation (OPP1039688), K. Dobos (PI) and the National Institute of Health (R01A1104589), P. Nahid (PI). The Bill and Melinda Gates Foundation did not have any role in study design, sample collection, analysis, and interpretation of data or writing of the manuscript.

Publisherā€™s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Nicole A. Kruh-Garcia.

Additional files

12014_2017_9156_MOESM1_ESM.xlsx

Additional file 1. Table of native and isotope-labeled peptide multiple reaction monitoring parameters including parent and transition m/z, retention time and collision energy.

12014_2017_9156_MOESM2_ESM.xlsx

Additional file 2. Limit of Detection/Quantitation (LOD/LOQ) summary for transitions in each peptide in mixed exosome matrix. Values are in fmol/Ī¼l.

12014_2017_9156_MOESM3_ESM.xlsx

Additional file 3. Table of normalized TPA cut-off values determined by the healthy donor group. False discovery rates (FDR) are calculated as the false positives (healthy controls with nTPA above the cut-off)/total number of healthy controls.

12014_2017_9156_MOESM4_ESM.tif

Additional file 4. Graph detailing the patient breakdown by smear microscopy and HIV status and the number of peptides identified in each class.

Additional file 5. Number of peptides identified per patient by geography.

12014_2017_9156_MOESM6_ESM.xlsx

Additional file 6. Statistical analysis summary. Includes results from two-tailed t-test analysis, before and after Benjamin-Hochberg adjustment.

12014_2017_9156_MOESM7_ESM.tif

Additional file 7. Graphs depicting the sixteen peptides in which several active TB patients displayed TPAs above the healthy threshold; the overall nTPA mean (solids line) between the two groups was not statistically significant.

12014_2017_9156_MOESM8_ESM.doc

Additional file 8. Summary of seventeen peptides which failed to discriminate between active TB and healthy controls; panel 1 and 2 are the peptides from MRM assay 1 and 2, respectively.

12014_2017_9156_MOESM9_ESM.tif

Additional file 9. Three peptides for which no cut-off threshold was determined due to indeterminate results in all 20 healthy samples.

12014_2017_9156_MOESM10_ESM.tif

Additional file 10. A. Peptides retention time (line) and their coefficient of variation (bars) for all samples (darker) and final samples (lighter). Mtb peptides: orange, iRTs: black. B. Boxplot representing retention time deviation from the mean for the 20 Mtb peptides and 9 iRTs. All: includes all samples and replicates including truncated peaks. Final: includes only samples/replicates included in final analysis.

Rights and permissions

Open Access This 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.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Mehaffy, C., Dobos, K.M., Nahid, P. et al. Second generation multiple reaction monitoring assays for enhanced detection of ultra-low abundance Mycobacterium tuberculosis peptides in human serum. Clin Proteom 14, 21 (2017). https://doi.org/10.1186/s12014-017-9156-y

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12014-017-9156-y

Keywords