Screening of Microbial Volatile Organic Compounds for Detection of Disease in Cattle: Development of Lab-scale Method

The primary hurdle for diagnosis of some diseases is the long incubation required to culture and confirm the presence of bacteria. The concept of using microbial VOCs as “signature markers” could provide a faster and noninvasive diagnosis. Finding biomarkers is challenging due to the specificity required in complex matrices. The objectives of this study were to (1) build/test a lab-scale platform for screening of microbial VOCs and (2) apply it to Mycobacterium avium paratuberculosis; the vaccine strain of M. bovis Bacillus Calmette-Guérin; and M. kansasii to demonstrate detection times greater those typically required for culture. SPME-GC-MS was used for sampling, sample preparation, and analyses. For objective (1), a testing platform was built for headspace sampling of bacterial cultures grown in standard culture flasks via a biosecure closed-loop circulating airflow system. For (2), results show that the suites of VOCs produced by Mycobacteria ssp. change over time and that individual strains produce different VOCs. The developed method was successful in discriminating between strains using a pooled multi-group analysis, and in timepoint-specific multi- and pair-wise comparisons. The developed testing platform can be useful for minimally invasive and biosecure collection of biomarkers associated with human, wildlife and livestock diseases for development of diagnostic point-of-care and field surveillance.

diagnosis in humans and animals. These VOC emissions may be collected from breath, blood, skin, urine, feces and other bodily secretions. Studies have identified potential VOC biomarkers related to multiple diseases such as cholera, cancer, diabetes, uremia, schizophrenia, asthma, liver disease, chronic lung disease, Pseudomoniasis, tuberculosis, and others [13][14][15][16] . Several methods of collecting and analyzing VOCs for potential diagnosis of M. tuberculosis and M. bovis have been described. Closed loop stripping analysis (CLSA) -gas chromatographymass spectrometry (GC-MS) was used to detect VOCs from multiple strains of M. tuberculosis from cultures 17 , and select ion flow tube (SIFT)-MS has been used to measure VOCs present in the headspace of M. bovis BCG cultures 18 and to measure VOCs in the breath of children with cystic fibrosis 14 , NH 3 levels in human breath for Helicobacter pylori screening, and detect acetonitrile levels in smokers' breath 19 .
Electronic 'e-nose' technology has been used for detection of VOCs present in sputum samples collected from humans infected with M. tuberculosis 14 , and VOCs from the headspace of cultures of M. bovis BCG and M. smegmatis 18 . 'E-nose' has also been reported capable of differentiating between M. tuberculosis, three other bacteria, and a control 19 , monitoring smokers' habits by measuring breath CO, detecting H. pylori presence in association with chronic gastritis, and detecting N 2 O produced by respiratory inflammation 16 . Thermal desorption (TD)-GC-MS has been used to determine VOCs present in the headspace of M. bovis BCG cultures 18 , and cattle breath 20 ; and solid-phase microextration (SPME)-GC-MS was used to determine biomarkers M. tuberculosis and M. bovis in cultures and M. tuberculosis in human breath 21,22 . Other methods that have been used to collect VOCs associated with diseases that could be utilized in the future for TB detection include ion mobility spectrometry (IMS) 14,23 , proton transfer reaction (PTR)-MS [14][15][16] , and laser spectroscopy 24 . Select ion flow tube-MS,  'e-noses' , IMS, PTR-MS have the advantage of being fast and potentially mobile. The downside of these methods includes decreased sensitivity, and the inability to chemically identify or profile VOCs. GC-MS-based methods may be slower, more expensive, and the instrumentation is not typically mobile; however, they have the advantage of being able to reproducibly identify and profile known and unknown microbial VOCs at low concentration ranges 14,16 .
Solid-phase microextration is an attractive technology for collecting microbial VOCs due to its simplicity, ease of use, and ability to sample and pre-concentrate a wide range of potential target compounds. Dynamic headspace extraction of VOCs can increase mass transfer to the SPME fiber compared to static headspace extraction, thus, reducing sampling times while improving mass loading of the fiber 25,26 . In contrast, both CLSA and TD are capable of pre-concentrating VOCs but require extra equipment and are more labor intensive than SPME. Syhre et al. 21 collected VOCs from seven mycobacterial and 16 other respiratory pathogen cultures using three different SPME fiber types; 100 µm polydimethylsiloxane (PDMS), 2 cm 50/30 µm divinylbenzene (DVB)/Carboxen/ PDMS, and 70 µm Carbowax/DVB. The SPME fiber coated with 2 cm 50/30 µm DVB/Carboxen/PDMS was found to recover higher concentrations of all target VOCs. Other studies have utilized the 50/30 µm DVB/Carboxen/ PDMS-type and Carboxen/PDMS-type SPME fiber coatings for microbial VOC collection in human and cattle breath samples 22,27,28 , in bovine nasals excretions 29 , and in bovine fecal excretions from cattle vaccinated with M. bovis challenge 30 . SPME has been also used for the in-vivo and in-vitro collection of rumen gases 31,32 , VOCs emitted from wildlife marking fluids 33,34 , decaying carcasses 35 .
The objectives of this research were to (1) design, build, and test a lab-scale dynamic VOC sampling platform specifically capable of simultaneous biosecure SPME collection of headspace VOCs emitted from controlled bacterial cultures and a media control; (2) apply this screening method to M. avium paratuberculosis (MAP); the vaccine strain of M. bovis Bacillus Calmette-Guérin (BCG); and M. kansasii cultures to demonstrate a proof-of-concept detection method that is faster than standard culture methods.
Our first working hypothesis was that by sampling the recirculating culture headspace air, we would be able to detect trace-levels of microbial VOCs early in the incubation process with minimal background interference. Our second hypothesis was that SPME sampling of microbial VOCs followed by GC-MS analysis would be suitable to detect differences in microbial VOCs emitted by different cultured strains of mycobacteria.
If successful, this proof-of-concept identification of VOC biomarkers would allow differentiation between the microbial agents prior to the eight weeks often required for cultured mycobacterial strain identification. The knowledge gained from this work could be directly applied for diagnosis of hTB and bTB, for detection and identification of VOC biomarkers produced by culture of other pathogenic bacteria, and as a reference library for pathogen-produced VOCs present in other samples such as breath, feces, urine, blood, and other biofluids.

Results and Discussion
Reducing interfering background Vocs in lab-scale testing platform. To minimize the VOC background and improve the likelihood of finding unique microbial VOCs for Objective 2 the prototype sampling platform was exposed to an initial 21 h bake-out at 50 °C, which only decreased the background slightly. After bake-out of the Neoprene tubing, system background emissions were reduced by 75% (Fig. 1). Background emissions present in the microbial growth media were also determined as reproducible and necessary background populated by numerous VOCs with relatively low peak areas.
Multi-group analysis of single ion data. A total of 123 fragmentation ions were identified by XCMS Online, with 33 meeting the criteria for statistical significance (α = 0.05; fold change ≥1.5). Visualization of the statistically significant ion feature characteristics are depicted by XCMS Online as a cloud plot (Fig. 2) 36 . Briefly, all sample chromatograms were aligned and overlaid onto the x-axis. Significant ion features are identified as circles, with ions with greatest m/z ratios located furthest from the x-axis. Circle size is proportional to the degree of fold change (larger circle = greater fold change), while color intensity corresponds to the statistical significance (p-value) of the fold change as calculated by a Welch t-test with unequal variances (darker color = lower p-value).
Identification of culture specific peak area data. Statistically significant ion intensities were matched by retention time to 77 total ion chromatographic (TIC) peaks. For each culture, trends in the changes of potential VOCs TIC intensities across the three time-points in each of the replicates were graphed. This allowed identification of ten VOC compounds with consistent, repeatable changes across time in each of the three replicates (3 studies) (Figs 3, 4 and 5).
There has been little research exploring the metabolic pathways of M. bovis BCG, MAP, and M. kansasii, therefore, published literature exploring the metabolome of M. tuberculosis complex and environmental strains of mycobacteria were considered sources for comparison (Table 3).

conclusions
This study demonstrates a proof-of-concept for the detection and use of microbial VOCs as a means to discriminate between mycobacterial cultures associated with one to three week post-culture inoculation, a time span preceding the time required to currently identify pathogenic M. tuberculosis and M. bovis in diagnostic cultures. To accomplish this task (Objective 1), a sampling system was designed, built, and tested for controlled collection of microbial volatiles using closed-loop headspace airflow over microbial cultures, sampling and sample preparation with SPME, analyte separation and identification via GC-MS, and tentative compound identification using www.nature.com/scientificreports www.nature.com/scientificreports/ novel metabolomics databases and the NIST W8N08 library. The capability of this system to produce results more efficiently than some currently utilized diagnostic modalities such as culture exemplifies its potential use as a diagnostic tool. The lab-scale testing platform concept can be useful for minimally invasive and biosecure collection of marker volatiles associated with human, wildlife and production animal diseases for development of diagnostic non-invasive point-of-care tools, field surveillance technologies and strategies.
Objective 2 provided a comprehensive assessment of VOCs collected from the headspace of three different mycobacterial cultures and one control media sample. Discrimination between mycobacterial cultures was successful one, two, and three weeks post-culture inoculation, a time span preceding the time required to currently identify pathogenic M. tuberculosis and M. bovis in diagnostic cultures. Unique VOCs representing potential biomarkers were identified in two mycobacterial cultures (e.g., BCG, M. kansasii). No unique plausible biomarkers were identified for the MAP cultures; however, discrimination from the two other mycobacterial cultures was possible when the VOC profiles of all the cultures were examined in context.
From a diagnostic perspective, detection of VOCs produced by pathogenic mycobacteria at early states of culture growth could improve disease diagnosis and treatment, especially in developing countries where access to sophisticated laboratory diagnostics is limited. The capability to differentiate between human and zoonotic mycobacteria under such circumstances could improve the capability of physicians to more accurately diagnose tuberculosis patients, to differentiate between hTB and zoonotic bTB, and to appropriately dispense medication targeted toward the etiological disease agent. Lab-scale testing platform. The lab-scale testing platform was designed, built and tested using several general lab use components. The reusable testing platform consisted of four Omegaflex FPU100 peristaltic pumps (Omega Engineering Inc., Stamford, CT, USA) to circulate air (average flow rate of 162 mL min −1 ) through each modified culture bottle in closed loops with each separate loop containing an inline ~13 mL, glass sampling bulb (28526-U, Supelco, Bellefonte, PA, USA) for SPME extraction (Fig. 6). Each closed loop was constructed from PTFE Teflon tubing from the pumps to the culture bottles, to the sampling bulbs, back to the pumps; and Neoprene tubing sections in the peristaltic pumps. The inlet and outlet fittings on each culture bottle included 0.22 micron filters to prevent introduction of mycobacteria into the closed loop system and rendering the system to be self-contained and biosecure. The platform's compact size enabled the entire system to be placed in an incubator for optimal temperature for growing mycobacteria. The materials used to constuct the platform allowed relatively easy decontamination after and between-the-trials use.

Materials and
Microbial volatile organic compounds. A 2 cm 50/30 µm DVB/Carboxen/PDMS fiber (57348-U, Supelco, Bellefonte, PA, USA) SPME fiber was used to extract and pre-concentrate VOCs from circulating gas. Samples were collected by dynamic headspace extraction with SPME for 1 h at 37 °C using the lab-scale testing platform. After the sampling period, the SPME fiber with extracted VOCs was inserted into the 250 °C GC injector for 2 min for thermal desorption, sample introduction, analyte separation and analysis. Corporation, Round Rock, TX, USA) used for analysis was equipped with two columns connected in series. The non-polar pre-column was 30 m × 0.53 mm i.d.; film thickness, 0.50 µm with 5% phenyl polysilphenylene siloxane stationary phase (SGE BPX-5) and operated with constant pressure mode at 13.5 psi (0.92 atm). The polar analytical column was 30 m × 0.53 mm bonded polyethylene glycol (PEG) embedded in a synthetic glass (SGE SolGel-Wax) at a film thickness of 0.50 µm. System automation and data acquisition software were MultiTraxTM V. 10.1 (Microanalytics, Volatile Analysis Corporation, Round Rock, TX, USA) and ChemStation ™ (Agilent Technologies, Santa Clara, CA, USA). The GC run parameters were as follows: injector, 250 °C; column, 40 °C initial, 3 min hold, 7 °C min −1 ramp to 240 °C final, 8.43 min hold; carrier gas, UHP-grade helium (99.999%). The GC was operated in a constant pressure mode where the mid-point pressure, i.e., pressure between pre-column and analytical column, was always at 5.7 psi (0.39 atm) and the heart-cut sweep pressure was 5.0 psi. This type of columns configuration (in series with different polarity) does not lend itself to a classic retention index (RI) approach for tentative compound identification with n-alkanes. The MS full scan range was 34 to 150 m z −1 . Spectra were collected at 2 scans s −1 using full scan. The quadrupole MS was set to electron ionization (EI) mode with ionization energy of 70 eV. MS tuning was performed using the default autotune setting using perfluorotributylamine (PFTBA) daily.
Reducing interfering background Vocs in lab-scale testing platform. After construction, the platform's background VOCs were initially baked-out for 21 h at 50 °C (the maximum temperature that some plastic components in the platform could withstand) The Neoprene tubing in the peristaltic pumps was removed from the platform and baked-out separately at 110 °C for 18 h to remove any VOCs the might be released into the closed platform loop.

Data analysis.
Cultures were identified in raw data as culture 1-4 and sample collection time-point (Week 1, Week 2, Week 3) until the data analysis was complete to prevent bias. Total ion chromatograms (TICs) from all cultures at all time-points were analyzed using the multi-group comparison feature in XCMS Online to identify peak ion abundances that differed between the cultures at each time-point 37 . Statistically significant ion intensities were matched to GC column retention time chromatographic peaks using Agilent Mass Hunter software (Agilent Technologies, Santa Clara, CA, USA). Peak areas were determined using the TICs. Peak area data for each culture at each replicate and time-point were evaluated to determine the suite of VOCs best suited to provide optimal discrimination in multi-group and pair-wise comparisons. Peaks were tentatively identified using AMDIS deconvolution software 38 , the National Institute of Standards and Technology (NIST) W8N08 database 39 , the Kyoto Encyclopedia of Genes and Genomes (KEGG) database [40][41][42] , and the Human Metabolome Database (HMDB) [43][44][45][46] . Tentative compound identifications were made based on ≥65% match with compounds present in these libraries. Tentative metabolic sources for each compound were explored using KEGG, HMDB, and review of peer-reviewed literature 47 .