Next Article in Journal
Genetic Diversity and Distribution of Virulence-Associated Genes in Y. enterocolitica and Y. enterocolitica-Like Isolates from Humans and Animals in Poland
Next Article in Special Issue
Evidence of Antimicrobial Resistance and Presence of Pathogenicity Genes in Yersinia enterocolitica Isolate from Wild Boars
Previous Article in Journal
Plant Health and Sound Vibration: Analyzing Implications of the Microbiome in Grape Wine Leaves
Previous Article in Special Issue
Molecular Typing of Neisseria gonorrhoeae Clinical Isolates in Russia, 2018–2019: A Link Between penA Alleles and NG-MAST Types
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Development of a Multiplex Real-Time PCR Assay for Predicting Macrolide and Tetracycline Resistance Associated with Bacterial Pathogens of Bovine Respiratory Disease

1
Department of Statistics, University of Nebraska-Lincoln, Lincoln, NE 68583, USA
2
School of Veterinary Medicine and Biomedical Sciences, University of Nebraska-Lincoln, Lincoln, NE 68583, USA
3
U.S. Meat Animal Research Center, Agricultural Research Service, United States Department of Agriculture, Clay Center, NE 68933, USA
4
Department of Food Science and Technology, University of Nebraska-Lincoln, Lincoln, NE 68588, USA
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Submission received: 21 December 2020 / Revised: 5 January 2021 / Accepted: 12 January 2021 / Published: 13 January 2021
(This article belongs to the Special Issue Detection and Characterization of Drug-Resistant Organisms)

Abstract

:
Antimicrobial resistance (AMR) in bovine respiratory disease (BRD) is an emerging concern that may threaten both animal and public health. Rapid and accurate detection of AMR is essential for prudent drug therapy selection during BRD outbreaks. This study aimed to develop a multiplex quantitative real-time polymerase chain reaction assay (qPCR) to provide culture-independent information regarding the phenotypic AMR status of BRD cases and an alternative to the gold-standard, culture-dependent test. Bovine clinical samples (297 lung and 111 nasal) collected in Nebraska were subjected to qPCR quantification of macrolide (MAC) and tetracycline (TET) resistance genes and gold-standard determinations of AMR of BRD pathogens. Receiver operating characteristic curve analysis was used to classify AMR based on the qPCR results. For lung tissues, the qPCR method showed good agreement with the gold-standard test for both MACs and TETs, with a sensitivity of 67–81% and a specificity higher than 80%. For nasal swabs, qPCR results passed validation criteria only for TET resistance detection, with a sensitivity of 88%, a specificity of 80% and moderate agreement. The culture-independent assay developed here provides the potential for more rapid AMR characterization of BRD cases directly from clinical samples at equivalent accuracy and higher time efficiency compared with the gold-standard, culture-based test.

1. Introduction

Bovine respiratory disease (BRD) is one of the most common and costly cattle diseases and affects 97% of feedlots, or 16% of cattle, with economic losses estimated to exceed 1 billion USD/year [1,2]. BRD is a disease complex with multiple contributing factors, including environment, viruses, bacteria, and the host [3]. Disease onset is commonly initiated by viral infections, which may suppress host defense mechanisms, allowing opportunistic bacterial pathogens to replicate and colonize deeper in the lung. Environmental factors such as crowding, poor ventilation, weather, and weaning increase stress and also reduce host immunity [4]. The coexistence of multiple viral and bacterial pathogens often increases the severity of BRD outbreaks, leading to higher mortality and morbidity rates and associated economic burdens [5,6].
Among pathogenic bacteria causing BRD in feedlot cattle and neonatal calves, Mannheimia haemolytica is the most frequently isolated, followed by Histophilus somni and Pasteurella multocida [1]. The multifactorial nature of BRD is often a challenge for disease management. Administration of antibiotic drugs for treatment and/or metaphylaxis (control and prevention) purposes may potentially select for resistance among bacterial pathogen populations [7]. Recent studies have shown that the frequency of antimicrobial resistance (AMR) in BRD bacterial pathogens in feedlots is increasing, particularly resistance to macrolides (MACs) and tetracyclines (TETs) [8], which are among the most commonly used antibiotics on beef cattle feedlots in the United States [9,10,11,12]. The development of AMR could compromise the effectiveness of antibiotics for the treatment and control of BRD, potentially leading to higher mortality and lower productivity among cattle populations. Increasing AMR in BRD pathogens may also threaten public health, as the resistance genes found in BRD pathogens are contained within integrative conjugative elements (ICEs) that could potentially be transferred from BRD pathogens to zoonotic bacteria [10]. Hence, tools for rapidly and accurately detecting potential AMR are critical components of outbreak monitoring by veterinarians and producers. Such tools could be used as adjuncts to existing culture-independent BRD pathogen detection methods that use similar technology and workflows [13]. In addition, the use of these tools would enable rapid evaluation of the potential efficacy of antimicrobial interventions and ultimately benefit the beef industry and support antimicrobial stewardship by enabling informed selection of optimal drugs for BRD therapies.
Classic culture-based methods for determining AMR commonly involve isolating and identifying organisms on solid media and assessing the growth inhibition of isolated bacterial strains under a series of concentrations of the target antibiotic to determine the minimum inhibitory concentration (MIC) [14]. As described by the Clinical and Laboratory Standards Institute (CLSI), isolates are classified as susceptible (S), intermediate (I), or resistant (R) to an antibiotic based on the relationships between MIC measurements and “breakpoints”, which are usually determined by taking into consideration the clinical outcomes of infections when the antibiotic is used [15]. The culture-based approach is widely used and is still considered the “gold-standard” test by national and international surveillance programs for monitoring AMR [14]. However, this method requires the growth, isolation, and identification of the target pathogens in pure culture, which can be time consuming and challenging due to sample contamination or overgrowth with environmental organisms [13].
In recent decades, molecular methods for detecting and quantifying AMR genes have shown promise as potential alternatives to the gold-standard AMR detection [16]. Quantitative real-time polymerase chain reaction (qPCR), which uses hydrolysis probes to generate a fluorescence signal, enables real-time assessment of DNA amplification and quantification of the genetic materials in the original sample. The qPCR approach also permits multiplexed detection of several gene targets in a single reaction. A recent study demonstrated that this approach is extremely useful for detecting opportunistic bacterial pathogens in bovine clinical samples and has advantages over culture-based approaches, especially when multiple pathogens coexist in the presence of normal flora [13]. These advantages make qPCR a powerful tool for the detection of genes that confer AMR using existing sample types and workflows. The detection of AMR genes in clinical samples (lungs and respiratory swabs) would provide veterinarians and clinicians with information on the presence or absence of AMR genes within hours to guide therapy selection. Although the detection of AMR genes in a sample does not necessarily indicate that the resistance genes are carried by BRD pathogens, such detection can be used as an indicator for the potential risks of antimicrobial-resistant infections in BRD outbreaks. Although the use of qPCR for AMR detection in tandem with culture-independent pathogen detection assays has not yet been evaluated, such a combination would provide information on the co-presence of BRD pathogens and BRD-associated resistance genes [13].
Unlike the gold-standard method, which generates a phenotypic classification of AMR, qPCR assays generate cycle threshold (Ct) values, which are continuous values representing the number of cycles at which the fluorescence signal exceeds the threshold value (i.e., the background signal level) and are used to estimate the number of gene copies in the original sample. To accurately predict the AMR classification index using a molecular testing method, an extensive epidemiological evaluation of the method based on a comparison with gold-standard methods must be conducted. Hence, the objectives of the present study were to (1) develop a qPCR method for the quantification of MAC and TET resistance genes in BRD clinical samples; (2) determine the optimal cutoff to support translation of qPCR results to a phenotypic classification of AMR; and (3) evaluate the validity of AMR phenotypic classification by the qPCR assay.

2. Results

2.1. Characteristics of the Multiplex qPCR Assay

In silico analyses showed 100% identity of the primer and probe loci to strains of M. haemolytica and P. multocida (Table S1). For the H. somni strains, 100% identity to ICEtetR gene targets was observed for sequences in GenBank; the remaining three targets were not found in H. somni strains available at the time of analysis. These resistance genes are not unique to BRD pathogens and are found on the chromosomes and/or plasmids of non-BRD-causing microbes (Table S2). The sequence identity of the primers and probes with non-BRD pathogen bacterial species does not indicate off-target binding to other genes, as it is expected that these AMR gene sequences are found in various other bacterial species.
Preliminary evaluation of the assay using the reference strains indicated that the multiplex qPCR assay detected all of the targets in the M. haemolytica reference strain and had sensitivity for all four targets of <3.2 CFU/mL, with Ct values ranging from 37.34 (ICEtetR) to 35.70 (msrE) at the lowest detectable concentration (Table 1). The assay also correctly classified all M. haemolytica reference strains that had been determined to contain the target genes (Table S3). Due to the lack of sequences for non-M. haemolytica BRD pathogens in GenBank or other databases, a panel of H. somni (35 strains) and P. multocida (5 strains) from BRD cases that had been subjected to MIC testing was also evaluated (Table S4). The detection of MAC resistance genes by the assay corresponded with phenotypic MAC and TET resistance in 26 and 28 of 35 H. somni strains, respectively. For P. multocida, genotypic results that agreed with the phenotypes for both TET and MAC resistance were obtained for all five strains. To validate assay specificity, we also tested ATCC reference strains of Actinobacillus pleuropneumoniae, Mannheimia granulomatis, H. somni, and B. trehalosi that do not exhibit phenotypic resistance. No qPCR targets were detected in any of these reference strains. The susceptibility of two bovine respiratory strains of B. trehalosi was also evaluated; as expected, the tetracycline-resistant (TETr) and macrolide-resistant (MACr) strain was positive for all four targets, whereas the susceptible strain was negative for all targets.

2.2. Phenotypic Antimicrobial Resistance Based on the Gold-Standard Test

BRD pathogen detection. Along with the sample distribution, Table 2 shows the prevalence of BRD pathogens by sample type. In addition to the individual BRD pathogens examined in this study, a new pathogen category, “positive for at least one BRD pathogen”, was evaluated. Overall, of the submitted samples for BRD diagnosis, at least one BRD pathogen was detected in 97.8%. In general, the prevalence of samples with at least one BRD pathogen was higher for lung samples (99.7%) than nasal samples (97.3%). Consistent with previous reports, M. haemolytica was the most frequently isolated BRD pathogen, comprising 63.7% of tested samples, followed by P. multocida (42.5%) and H. somni (29.6%). However, the pathogen distribution varied by sample type (Figure 1). Specifically, among lung samples, M. haemolytica was the most prevalent pathogen (64.3%), followed by P. multocida and H. somni at similar levels (32.3% and 31.3%). By contrast, both M. haemolytica and P. multocida were relatively highly prevalent in nasal samples (64.9% and 72.1%), whereas H. somni was not (26.1%).
Phenotypic AMR data manipulation. The distributions of AMR classifications determined by the MIC test by BRD pathogen and sample type are presented in Table 3. BRD pathogens isolated from clinical samples were classified as S, I, or R based on the comparison of the MIC test results with the breakpoints listed in Table 4. Strains classified as intermediate were much less prevalent, resulting in an inadequate number of samples in this category for statistical analysis. This was the major reason for combining R and I into the newly defined category R+I, resulting in a binary classification of AMR status based on the gold-standard MIC test, i.e., R+I and S. AMR status was recorded for individual pathogens as well as for the category of “positive for at least one pathogen”. In addition, resistance to two drugs in the MAC class was evaluated, and “resistance” to tilmicosin, tulathromycin, or MACs as a class (either tilmicosin or tulathromycin) was recorded.
Resistant BRD pathogen isolation. Based on the newly defined category R+I, the prevalence of TETr BRD pathogens was higher than that of MACr pathogens for all pathogens and sample types. TET resistance prevalence ranged from 19.2% to 52.7%, while MAC resistance prevalence was between 8.8% and 39.3%. The only exception was M. haemolytica isolated from nasal samples, in which the prevalence of TET resistance (6.9%) was slightly lower than that of MAC resistance (8.3%). In addition, the occurrence of antimicrobial-resistant BRD pathogens was always higher in lung samples than in nasal samples. Specifically, 44.6% of lung samples and 23.1% of nasal samples possessed at least one TETr pathogen, while 36.9% of lung samples and 19.4% of nasal samples possessed at least one MACr pathogen. The same trends were observed for all resistant BRD pathogens.

2.3. Optimal Ct Cutoff Value Determination

Real-time PCR data manipulation. Three gene targets conferring resistance to MACs, i.e., msrE, mphE, and erm42, were detected using the qPCR method, and the average observed Ct value for the three targets was used. The use of the average value was deemed reasonable because of the high similarity of the distributions of the Ct values for the three genes (Figure 2). Although the distribution of the Ct values for mphE was slightly shifted from those for erm42 and msrE, the percentile distributions were similar, so the average value for the three genes is a good representation of the inherent distribution of any of the three genes. Hence, all qPCR results for MAC drugs reported in this study are based on the average measurement of the three gene targets.
Optimal Ct cutoff value determination. As shown in Table 5, for the lung samples, the optimal Ct cutoff value for TET was 36.06, with more than 80% sensitivity (Se) and specificity (Sp) and a good level of agreement with the MIC test (κ = 0.64). The prevalence of lung samples possessing at least one BRD pathogen resistant to TET was approximately 45%, leading to positive predictive value (PPV) and negative predictive value (NPV) estimates of 0.79 and 0.84; these values indicate that 79% of tested samples classified as resistant actually possessed at least one BRD pathogen resistant to TET and 84% of samples classified as susceptible were truly susceptible. Compared with TET resistance, the optimal Ct cutoff values for MAC resistance were lower, with lower Se but higher Sp of approximately 90%. The qPCR assay showed a good level of agreement with phenotypic MAC resistance, with κ ranging from 0.61 to 0.64. Similar to TETr strain detection, both the NPV and PPV of the qPCR assay for MAC resistance were approximately 80%, ensuring relatively low false-negative and false-positive rates. Figure 3 shows the receiver operating characteristic (ROC) plots for all of the antimicrobials for the lung and nasal samples, respectively. In general, the performance of the qPCR approach for detecting resistance was lower for nasal samples than for lung samples. Se and Sp were both lower for nasal samples, with ranges of 50%–88% and 79%–83%, respectively; and the concordance with the MIC tests was fair to moderate, as indicated by κ values ranging from 0.30 to 0.58.

2.4. Validation of the Computational Approach

Determination of the required sample size for calculating the optimal cutoff. The optimal Ct cutoff value was considered valid if the total sample size and the number of samples positive for resistant BRD pathogens were both greater than or equal to the minimum requirements established by considering the importance of both PPV and NPV. Estimates of the sample size necessary for the validity evaluation are listed Table S5. Clearly, the optimal Ct obtained for the lung samples for both TET and MACs satisfied the minimum requirement. For the nasal samples, only the optimal Ct values obtained for oxytetracycline and tilmicosin satisfied the requirements for both sample size and number of positive samples, whereas the optimal Ct values for tulathromycin and the combined MAC group did not satisfy the second requirement, i.e., the number of resistant samples included for determining the optimal Ct value was insufficient. If only NPV is of importance, then the optimal Ct obtained for all nasal samples is valid. However, increasing PPV or minimizing the false-positive rate helps reduce the overall misclassification rate and thus cannot be ignored completely.
Cross-validation. To validate the consistency of the optimal Ct cutoff value estimation, 5-fold and 10-fold cross-validation (CV) were conducted with a focus on the category of “resistance to at least one BRD pathogen”. Because more consistent estimates in terms of the optimal Ct cutoff and the κ value were obtained, 5-fold CV was used to optimize the optimal Ct cutoff value for the qPCR method. Table S6 shows the results of the 5-fold CV along with a comparison of the average level of agreement between the training and test sets.
The optimal Ct values for the lung samples for both TETs and MACs were consistent in terms of required diagnostic accuracy. Overall, for the lung samples, the optimal Ct obtained from the overall data fell within the 95% average optimal Ct obtained using 5-fold CV. In addition, the average κ evaluated based on the training and test sets indicated a good level of agreement with the results of the MIC tests for both groups of antimicrobials. However, the optimal Ct obtained from the nasal samples was sufficiently validated only for the TET group. For the MAC group, a disparity in the average κ was observed between the training and test datasets.
Table 6 summarizes the optimal Ct values and the level of agreement based on the minimum requirements for total sample size and resistant sample size and 5-fold CV. The optimal Ct for predicting the phenotypic AMR classification based on qPCR results was sufficiently validated for the lung samples containing MACr and/or TETr BRD pathogens and for the nasal samples containing TETr BRD pathogens. The optimal Ct values for the lung samples were 35.66 and 33.12 for TET and MAC resistance, respectively, with a good level of agreement between the gold-standard and the qPCR approach. The optimal Ct value obtained for the nasal samples was 33.27 for TET resistance, with a moderate level of agreement between the gold-standard and the qPCR approach.

3. Discussion

The method described in this study enables AMR gene detection in clinical samples in a culture-independent manner and thus can rapidly provide clinicians, veterinarians, and cattle producers preliminary information on the presence of potential AMR in BRD cases. Previous examinations of AMR in BRD pathogens have relied on whole-genome sequencing (WGS) approaches to mine genomic sequences for resistance genes. Although this approach is much more comprehensive, it is less practical for use with clinical samples and in situations where a clinician may need rapid results for drug decisions. Compared with isolates, working with clinical samples from livestock environments poses a significant challenge due to the large variations in sample types and quality, environment, and collection requirements. Working with clinical samples requires a targeted and quantitative approach such as qPCR, which has been shown to increase the detection of pathogens compared with culture-dependent methods such as WGS [13].
In developing this assay for use in a rapid format with clinical samples, our goal was to select robust gene targets with spatial and temporal consistency in BRD pathogen strain genomes and that confer resistance to highly relevant drug classes. Consequently, the performance of the method for both TETs and MACs was evaluated using ATCC- and clinical sample-sourced strains that have been subjected to MIC testing and/or WGS. Overall, the developed multiplex qPCR method performed well with high specificity for detecting the target AMR genes. In the initial validation, evaluation of the phenotypic–genotypic agreement using less comprehensive methods showed high levels of agreement for M. haemolytica, P. multocida and B. trehalosi. By contrast, for H. somni, disagreement between the results of the multiplex qPCR assay and phenotypic AMR classification using MIC was observed for 9 of 35 strains tested. Potential explanations for this discrepancy include alternative mechanisms of phenotypic resistance or inactive resistance genes in these strains.
The results indicate that this method may be particularly useful for assessing the presence of TET or MAC resistance in BRD pathogens in general or in M. haemolytica, one of the most dominant bacterial pathogens. TETs and MACs were chosen as the target drug classes because they are often used as first-line treatments for BRD. More than 70% of drugs used for the treatment or prevention of BRD in the United States are MACs (including tulathromycin, tilmicosin, tildipirosin, and gamithromycin), followed by TETs (mainly oxytetracycline, ~9%) [9,17,18]. Determinants conferring resistance to MACs or TETs in the population of BRD pathogens are frequently carried on ICE elements, revealing a potential challenge in the use of these important drugs for BRD disease treatment and/or prevention [19,20]. The ability to rapidly identify AMR to these drugs would be critically useful to veterinarians in the field, who could rapidly adjust their treatment or prevention therapies in response to the detection of such resistance [19].
One of the largest challenges in developing our assay was ensuring the concordance between phenotypic AMR status and the detection of genes conferring resistance in complex clinical samples. Investigations of concordance between genotypic and phenotypic resistance based on WGS or selected genes are increasing [21,22], and like the present study, most previous works have found substantial variation in the level of concordance depending on the combination of bacterial species, resistance gene(s) and resistance type. To construct our assay, we selected genes conferring resistance to TETs and MACs that are most commonly found in BRD pathogens and observed good agreement with corresponding phenotypic resistance for some combinations. However, discrepancies between genotypic and phenotypic resistance are not surprising, given the challenges of working with clinical samples outlined above and the large number of potential resistance genes. In one of the most comprehensive databases, ResFinder, more than 140 and 180 genes or gene variants encoding TET or MAC resistance, respectively, have been reported to date [23]. A more comprehensive gene panel may increase the predictive capability of our proposed qPCR assay of resistance in BRD bacterial pathogens. Our in silico analysis showed that the gene targets are potentially carried by other non-BRD pathogen bacteria, which might cause false-positive results if these bacteria were present in sufficient numbers in clinical samples.
The performance of the PCR-based rapid detection method was most optimal for predicting collective phenotypic resistance among all coexisting BRD pathogens in a single sample. In the MIC test, AMR classification is pathogen dependent, whereas qPCR detects genes in a culture-independent manner regardless of bacterial source. Therefore, a collective phenotypic resistance status for a sample, defined as possessing at least one resistant bacterial pathogen, enables a more rational comparison between the results of MIC testing and qPCR. For predicting phenotypic resistance in individual pathogens, the assay showed the most promise for M. haemolytica in lung tissues for both TETs and MACs. Lower agreement was observed for nasal swabs or for H. somni isolates, likely due to challenges in the use of the gold-standard method with nasal swab samples and the isolation of H. somni from all sample types. Nasal swab samples are frequently contaminated with environmental bacteria, which makes recovery of BRD pathogens challenging. Additionally, because H. somni is sensitive to oxygen and a fastidious organism, recovery is low, and false-negative testing results are common. Overall, the results indicate that the newly developed PCR-based rapid detection method is valid for determining the potential likelihood of resistance to TETs and MACs in at least one BRD pathogen or M. haemolytica, the most prevalent BRD pathogen, isolated from bovine lung tissues.
The ROC approach used in our study is in essence a classification method. The gold-standard, culture-based MIC approach gives a dichotomous outcome of resistance or susceptibility to antimicrobials, while the qPCR assay uses a Ct cutoff to transpose the continuous values produced into a finding of resistance or susceptibility. If the purpose of the qPCR assay was to rule out the presence of resistance in a clinical sample for BRD diagnosis based on a susceptible result, a higher Ct cutoff could be used to generate a “rule-out” test with a high Se or low false negatives in detecting resistance. Conversely, a lower Ct cutoff would be more likely to identify infections with the resistance of interest; such a “rule-in” test would have a higher Sp and lower false positives, resulting in higher confidence in the presence of resistance in a sample when a resistant result is observed. In the present study, the ROC approach was used to determine the optimal Ct cutoff by prioritizing a balance between Se and Sp. However, because the qPCR approach provides continuous information, it allows for flexibility in assay design for multiple purposes (such as a priority for “rule-in” or “rule-out” purpose) in real-world scenarios by varying the Ct cutoff.
Although the properties of diagnostic tests are generally expressed by Se and Sp, PPV and NPV are also considered important indicators of a test’s usefulness. Unlike Se and Sp, PPV and NPV depend heavily on prevalence, which in the present study was the prevalence of resistant BRD pathogens determined using the gold-standard method. A low prevalence of AMR would lead to a higher false-positive rate and imply lower PPV, whereas a high prevalence of AMR would lead to a higher false-negative rate and imply lower NPV. In the present study, AMR was more likely to be detected in lung samples than in nasal samples (Table 5), which is largely attributable to the fact that the prevalence of AMR in the samples in this study was lower among nasal samples (Table 3). As a result, the nasal samples had a lower prevalence and PPV, which is the primary reason why the optimal Ct cutoff values obtained for nasal samples were not valid. By contrast, the prevalence of AMR was considerably higher among the lung samples, providing a more balanced outcome in favor of the estimation of Se, Sp, PPV and NPV. Compared with the lung samples, the estimates of Se and Sp for the nasal samples were similar, but PPV was much lower, resulting in invalid estimates of the optimal Ct cutoff values. These results suggest that the most effective application of the qPCR assay for phenotypic AMR classification is situations in which the presence of AMR is not uncommon, such as during BRD outbreaks.
Although this novel application for predicting AMR phenotypes based on genotypic data could potentially be translated into a rapid assay for use in a variety of platforms, including pen side or in veterinary clinics, several limitations remain to be addressed. First, the data on phenotypic resistance were recategorized into binary outcomes due to the insufficient number of samples with intermediate resistance, and the ROC approach is less optimal for discriminating among more than two classifiers. Second, the classification was conducted for specific pathogen–sample–resistance type combinations in a series of stratified analyses but did not include predictors other than the Ct value, such as historical or clinical information like antimicrobial treatments, which might greatly influence the development of resistance. Third, to estimate the optimal Ct value, the ROC approach employed area under the curve (AUC), which is a widely used methodology but may be further improved because its use of a threshold of 50% for classification ignores the actual probability. To potentially address these limitations and increase the accuracy of diagnosis of AMR among BRD infections using PCR-based rapid detection methods, alternative classification methods, such as the H measure approach, warrant investigation [24].
In conclusion, the rapid multiplex real time-based PCR detection assay reported here, which can estimate levels of AMR genes potentially associated with multiple species of BRD pathogens, holds promise for the rapid detection of these genes in complex samples. By combining this gene detection assay with a pathogen-specific assay that can quantify pathogen abundance in samples, veterinarians could rapidly assess the risk of BRD caused by a potentially resistant pathogen, thereby enabling judicious use of antimicrobials in cattle production systems. Further expansion of the number of gene targets may enable more robust AMR assessments of complex samples.

4. Materials and Methods

4.1. Sample Collection and Distribution

Bovine lung tissue and nasal swab samples were collected over a one-year period in 2018 from submissions, mainly for BRD diagnostic testing, to the University of Nebraska-Lincoln Veterinary Diagnostic Center and combined with an archived sample collection that had been stored at −80 °C from 2012 to 2017. Although detailed clinical information was not available for many submissions, most of these samples were likely collected from non-healthy, clinically ill animals since they were submitted to a veterinary diagnostic lab for BRD diagnosis. A total of 416 bovine clinical samples were collected, primarily consisting of lung tissues (n = 297) and nasal swabs (n = 111). Omitted samples that were excluded from further data analysis were either tissues less related to BRD diagnosis (n = 4, i.e., skin and liver samples) or lacked records about tissue type (n = 4). The samples were analyzed for the existence of target BRD bacterial pathogens. If the pathogens were detected, phenotypic AMR characteristics were subsequently determined using the gold-standard culture-based method. In addition, target genes conferring AMR in the samples were quantified using the culture-independent qPCR rapid detection method developed in this present study.

4.2. Reference Strains

ATCC-sourced strains, a collection of field isolates from diagnostic cases, and off-target controls were used for a preliminary evaluation of the target specificity of the qPCR assay. M. haemolytica strains that had been previously subjected to WGS and with established presence or absence of the AMR targets and ICE elements were used to capture the known genomic diversity of this pathogen (Table S3 [19,25]). H. somni, P. multocida, and B. trehalosi strains isolated from clinical cases that had been subjected to MIC testing were also used for evaluation purposes along with ATCC reference strains (Table S4). Bibersteinia trehalosi was included in the reference set because it is occasionally isolated from BRD cases and has shown AMR [26]. This collection of reference strains was employed only for evaluating assay performance and specificity and was excluded from further data analysis of the clinical samples.

4.3. Molecular-Based Rapid Detection Assay

Nucleic acid extraction. For initial validation, nucleic acids were extracted from reference strains in pure subculture, as described previously [13]. Briefly, a single purified colony was resuspended in nuclease-free water, boiled at 100 °C for 10 min and clarified by centrifugation at 15,700 RCF for 2 min. Swabs, which were submitted to the laboratory in various liquid transport media from referring veterinarians, were vortexed vigorously, and 100 µL of the transport medium was combined with lysis solution and extracted using a MagMax Total Nucleic isolation kit (AM1840) on a Kingfisher flex instrument (Thermo Scientific) following the manufacturer’s instructions. Lung tissue samples were added to a filter Whirl-Pak (Nasco, Fort Atkinson, WI) with 1–5 mL of sterile phosphate-buffered saline (PBS) and placed in a stomacher for 30–60 s. Following stomaching, 2 mL of suspension was removed from the bag, centrifuged at 15,700 RCF for 2 min and extracted using the MagMax procedure and kit, as described above.
qPCR assay. A multiplexed hydrolysis nucleic acid probe assay was designed based on whole-genome sequences of M. haemolytica strains possessing ICE elements containing AMR genes. Four gene targets carried on the ICEs were selected, including three that confer resistance to MACs (msrE, mphE, and erm42) and one that is the repressor gene for TET class drug resistance (ICEtetR). The primer and probe sequences were designed by using PrimerQuest software (IDT, Iowa City) using genomic sequences from a reference strain of M. haemolytica [19] (Table 7). In silico analyses of the primer and probe sequences were conducted to evaluate the use of these targets in addition to M. haemolytica strains with the ICEs and four AMR genes and in strains of P. multocida and H. somni.
Following in silico analysis, primers and probes were purchased from Integrated DNA Technologies (IDT) and used in a single 4-plex reaction in a rotary-based real-time PCR instrument (Qiagen Rotorgene Q). Protocols were optimized for primer and probe concentrations and annealing temperatures. The PCR master mix had a total volume of 25 µL comprising 12.5 μL of 2× Quantifast Multiplex PCR Master Mix (Qiagen), 1 μL of each primer probe mix (4 μL total) containing F (10 μM), R (10 μM), and P (10 μM), 6.5 μL of nuclease-free water and 2 µL of template DNA. The thermocycling conditions were as follows: 95 °C for 5 min, followed by 45 cycles of 95 °C for 15 s, and 59 °C for 40 s. The Ct threshold was set at a fixed value of 0.1 for all detection channels following dynamic tube normalization.
Following optimization, assay performance was tested using reference strains, and limit of detection analysis was conducted using a reference M. haemolytica strain known through WGS to have all four genes (Reference strain 1621 in Table S3) [19]. The assay was further validated using a panel of M. haemolytica strains (Table S3) previously characterized by WGS and representing all known major genotypes and subtypes of M. haemolytica and different combinations of AMR genes [19]. The assay gene targets have also been reported in strains of H. somni and P. multocida isolated from cattle, and therefore strains of H. somni and P. multocida with phenotypic resistance patterns similar to those of the M. haemolytica strain collection were also included in this validation (Table S4). In addition, analyses of non-target closely related strains, including B. trehalosi, and other ATCC reference strains were performed (Table S4).

4.4. Preliminary Assay Validation and Analytical Sensitivity

Assay performance was evaluated for all four targets using DNA extracted from M. haemolytica reference strain 1621. Briefly, the strain was grown in pure culture on tryptic soy agar with 5% sheep blood (BAP) (Remel, Lenexa, KS). Colonies were picked and cultured in 10 mL of brain heart infusion broth in a 100 mL flask shaken at 200 RPM for 12 h. Serial dilutions (1:10) of this culture were plated on BAP to enumerate colonies and estimate colony forming units/mL, and the remaining portion of each dilution was used for nucleic acid extraction, as described above, to evaluate assay performance and limit of detection (Table 1). Following the determinations of limit of detection and dynamic range, the target specificity of the assay was evaluated using a diverse collection of reference strains and controls (Table S4).

4.5. Pathogen Isolation and Determination of Phenotypic Antimicrobial Resistance Characteristics

BRD pathogen isolation and confirmation. All clinical samples were processed for bacterial isolation and identification of the BRD pathogens M. haemolytica, P. multocida, and H. somni, as described previously [13]. Briefly, samples were documented, processed, and tested by trained personnel following validated and approved standard operating procedures in an American Association of Veterinary Laboratory Diagnosticians (AAVLD)-accredited diagnostic laboratory. For bacterial culture, diseased portions of fresh lung tissues were excised using flame-sterilized scissors and dipped in alcohol. After flame sterilizing the exterior, the sample was bisected, and the cut surface was used to inoculate culture media. Swabs were used directly to inoculate plates if submitted in liquid medium or otherwise were resuspended in approximately 1.5 mL of PBS. Colonies were isolated by streaking on culture media, including BAP, chocolate agar, and MacConkey agar (Remel, Lenexa, KS), incubated in 5% CO2 and examined at 24 and 48 h after inoculation. Colonies with morphologies consistent with M. haemolytica, P. multocida, or H. somni were identified by MALDI-TOF MS (Bruker Biotyper) testing using the manufacturer’s validated procedures for definitive identification of these organisms.
MIC determination. For each confirmed BRD pathogen strain, the presence of phenotypic AMR was evaluated using oxytetracycline to represent the TET class and tilmicosin and tulathromycin to represent the MAC class. As the gold-standard approach, AMR to TETs and MACs was evaluated by MIC testing using the broth microdilution method according to CLSI guidelines [27]. By comparing the MIC values with the corresponding breakpoints for bovine respiratory pathogens (Thermo Sensititre using the BOPO6F panel), which are listed in Table 4, BRD pathogen isolates were classified as S, I or R to each drug [28]. In this study, intermediate and resistant were both re-defined as “resistant” (R+I) at the strain level. A sample was classified as containing “resistant” BRD pathogens if at least one strain of antimicrobial-resistant BRD pathogens evaluated in this study was detected.

4.6. Predicting Phenotypic Antimicrobial Resistance

Statistical analysis was performed to optimize and validate the capability of the developed qPCR assay to classify samples as either S or R+I to a given antimicrobial. First, optimal Ct cutoff values were determined using ROC analysis to maximize the assay accuracy (both Se and Sp as the primary accuracy indices) compared to the gold-standard test. Optimal cutoff values were determined for clinical samples containing M. haemolytica, P. multocida or H. somni resistant to OTC, tilmicosin or tulathromycin, as well as for samples classified as “possessing at least one resistant BRD pathogen”. Second, the diagnostic accuracy and validity of the developed qPCR assay were assessed by determining the sufficiency of the sample size collected in this study for detecting MAC and TET resistance in samples possessing BRD pathogens and by using the cross-validation technique. All statistical analyses were performed in statistical software R version 3.5.1 [29].

4.6.1. Optimal Ct Cutoff Value Determination

Data preparation for ROC analysis. To determine the optimal Ct, ROC curves were plotted by comparing the phenotypic AMR classification determined by the gold-standard test with the continuous Ct value generated by the qPCR assay using the ‘pROC’ package in statistical software R [30]. ROC curve analysis is a common method for evaluating diagnostic tests with a binary classifier [31], but this method does not work optimally if there is no detection or anomaly observed in either the qPCR assay or MIC test. In this study, no detection by qPCR means that the target gene is either present in very small amounts or completely absent. Hence, samples recorded as “no gene copies detected” were imputed with a Ct value of 40, the highest amplification cycle observed in this study. For the gold-standard test, it is not easy to reinvestigate missing results. Hence, another four samples were further excluded from data analysis due to a lack of MIC records, resulting in a total of 404 samples (296 lung and 108 nasal) for inclusion in the following data analysis.
ROC analysis to determine the optimal Ct cutoff. In this study, ROC curves were generated by plotting the true-positive rate (i.e., Se) against the false-positive rate (i.e., 1–Sp) to determine the diagnostic equivalency between the qPCR assay and the gold-standard test over different cutoff values. Table 8 shows a schematic representation of the comparison of the outcomes of the qPCR assay and the culture-based gold-standard test and the calculation of Se and Sp. The optimal cutoff translating the continuous Ct value into a binary classification of R+I or S was determined by optimizing the AUC [32]. AUC can take values between 0.5 and 1, with values closer to 1 indicating high performance of the qPCR approach in predicting the resistance classification with minimum classification error compared with the gold-standard test.
Performance of qPCR. As a post-analysis statistic, the determined optimal Ct cutoff value was used to calculate Cohen’s Kappa (κ) in order to measure the agreement between the two tests [33]. κ ranges from −1 to 1, with higher values implying greater concordance between the two tests. The following interpretations of κ were used: κ < 0.2, poor agreement; 0.2 < κ ≤ 0.4, fair agreement; 0.4 < κ ≤ 0.6, moderate agreement; 0.6 < κ ≤ 0.8, good agreement; and κ > 0.8, very good agreement [34]. κ was calculated using Equation (1):
κ = A c t u a l   a g r e e m e n t   b e y o n d   c h a n c e   P o t e n t i a l   a g r e e m e n t   b e y o n d   c h a n c e   = p 0 p e 1 p e , w h e r e   p 0 =   p r o p o r t i o n   o f   t r u e   a g r e e m e n t = T P + T N T P + F N + F P + T N   a n d p e = p r o p o r t i o n   o f   r a n d o m   a g r e e m e n t = T P F N T P F P + F P T N F N T N ( T P + F N + F P + T N ) 2 , T P F N = T P + F N ,   T P F P = T P + F P ,   F P T N = F P + T N   a n d   F N T N = F N + T N , T P ,   T N ,   F N   a n d   F P   a r e   r e p o r t e d   i n   Table   8 .
The PPV and NPV were also calculated to measure the true-positive and true-negative results of the multiplex qPCR assay, as shown in Table 8. In this study, PPV is the proportion of samples showing a test result of “resistant” based on the multiplex qPCR assay that were actually identified as resistant using the gold-standard test, while NPV is the proportion of samples showing a test result of “susceptible” in the qPCR assay that were also identified as susceptible using the gold-standard test.

4.6.2. Diagnostic Accuracy Evaluation

Developing a multiplex qPCR assay with relatively high diagnostic accuracy for classifying AMR status in clinical samples is contingent on a balance of samples with positive and negative responses to the outcomes of interest, i.e., R+I and S [35]. To determine whether this balance was reached, the minimum required sample size and the required proportion of resistant samples were estimated and compared with the sample sizes and the proportions of resistant samples included in this study to validate the optimal Ct cutoff values determined using the ROC approach. Once the sufficiency of the sample size and required number of resistant samples were determined to be adequate, cross-validation (CV) was conducted to evaluate the computational approach used in ROC analysis to determine the optimal Ct value. Cross-validation is one of the most common resampling methods to evaluate the performance of a diagnostic test [36]. Here, the 5-fold CV technique was used to validate the optimal Ct value of the multiplex qPCR assay for classifying AMR equivalent to the gold standard. Further details of the accuracy evaluation techniques are provided in the Supplementary Materials, Text.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/2076-0817/10/1/64/s1, Table S1: M. haemolytica in silico primer and probes analysis; Table S2: Primer binding other species; Table S3: Isolates used in MALDI study (M. haemolytica WGS isolates); Table S4: Target and non-target isolates used to validate assay specificity, with respective threshold cycle (Ct) values obtained from multiplex qPCR and MIC values for antimicrobial resistance determination; Table S5: Assessment of sufficiency of sample size for determining the optimal cycle threshold value; Table S6: Comparison of 5-fold cross-validation with the optimal threshold cycle (Ct) cutoff value obtained using ROC curves on the overall data; Text: Elaboration of accuracy evaluation techniques used in the present study.

Author Contributions

Conceptualization, J.D.L. and B.W.; formal analysis, E.D. and J.D.L.; funding acquisition, J.D.L. and B.W.; investigation, C.A.D.; methodology, E.D., J.D.L., J.C. and B.W.; project administration, J.D.L. and B.W.; resources, E.L.W. and M.L.C.; software, E.D.; visualization, E.D. and C.A.D.; writing—original draft, E.D., J.D.L. and B.W.; writing—review and editing, E.L.W., M.L.C. and J.C. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the Nebraska Beef Council (#109924). Partial support for J.D.L. and C.A.D. was provided by the Nebraska Experiment Station with funds from the Animal Health and Disease Research (section 1433) and Hatch capacity funding programs (accessions# 1002196, 1007070 and 1017646) and Multistate Research Project of Antimicrobial Resistance (NC 1206, accession# 1014035) through the USDA National Institute of Food and Agriculture. M.L.C. and E.L.W. received support from the United States Department of Agriculture (USDA)/Agricultural Research Service (ARS) in support of the project “Genomic Intervention Strategies to Prevent and/or Treat Respiratory Diseases of Ruminants” (project# 3040-32000-034-00-D).

Institutional Review Board Statement

This use of animal derived materials and samples was reviewed and approved by the Institutional Animal Care and Use Committee of the University of Nebraska-Lincoln, and AAALAC accredited institution. (ID # 1970, approval 26 May 2020).

Informed Consent Statement

Not applicable.

Data Availability Statement

All data are available in the main text or as electronic supplementary.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results. The use of product and company names is necessary to accurately report the methods and results; however, the United States Department of Agriculture (USDA) neither guarantees nor warrants the standard of the products, and the use of names by the USDA implies no approval of the product to the exclusion of others that may also be suitable. The USDA is an equal opportunity provider and employer.

References

  1. Griffin, D.; Chengappa, M.M.; Kuszak, J.; McVey, D.S. Bacterial pathogens of the bovine respiratory disease complex. Vet. Clin. North Am. Food Anim. Pr. 2010, 26, 381–394. [Google Scholar] [CrossRef] [PubMed]
  2. Amat, S. Bovine respiratory disease in feedlot cattle: Antimicrobial resistance in bovine respiratory bacterial pathogens and alternative antimicrobial approaches. Bact. Cattle Dis. 2019. [Google Scholar] [CrossRef] [Green Version]
  3. Faber, R.; Hartwig, N.; Busby, D.; BreDahl, R. The Costs and Predictive Factors of Bovine Respiratory Disease in Standardized Steer Tests. 2000. Available online: https://lib.dr.iastate.edu/beefreports_1999/24 (accessed on 23 September 2020).
  4. Campbell, J. Bovine Respiratory Disease Complex. Available online: https://www.merckvetmanual.com/respiratory-system/respiratory-diseases-of-cattle/bovine-respiratory-disease-complex (accessed on 30 November 2020).
  5. D’Amours, G.H.; Ward, T.I.; Mulvey, M.R.; Read, R.R.; Morck, D.W. Genetic diversity and tetracycline resistance genes of Histophilus somni. Vet. Microbiol. 2011, 150, 362–372. [Google Scholar] [CrossRef] [PubMed]
  6. Arcangioli, M.A.; Duet, A.; Meyer, G.; Dernburg, A.; Bezille, P.; Poumarat, F.; Le Grand, D. The role of Mycoplasma bovis in bovine respiratory disease outbreaks in veal calf feedlots. Vet. J. 2008, 177, 89–93. [Google Scholar] [CrossRef]
  7. Anholt, R.M.; Klima, C.; Allan, N.; Matheson-Bird, H.; Schatz, C.; Ajitkumar, P.; Otto, S.J.G.; Peters, D.; Schmid, K.; Olson, M.; et al. Antimicrobial susceptibility of bacteria that cause bovine respiratory disease complex in Alberta, Canada. Front. Vet. Sci. 2017, 4, 207. [Google Scholar] [CrossRef] [Green Version]
  8. Portis, E.; Lindeman, C.; Johansen, L.; Stoltman, G. A ten-year (2000–2009) study of antimicrobial susceptibility of bacteria that cause bovine respiratory disease complex—Mannheimia haemolytica, Pasteurella multocida, and Histophilus somni—in the United States and Canada. J. Vet. Diagn. Investig. 2012, 24, 932–944. [Google Scholar] [CrossRef] [Green Version]
  9. USDA. Feedlot 2011 Part IV: Health and Health Management on U.S. Feedlots with a Capacity of 1000 or More Head. Available online: https://www.aphis.usda.gov/animal_health/nahms/feedlot/downloads/feedlot2011/Feed11_dr_PartIV_1.pdf (accessed on 28 August 2020).
  10. Klima, C.L.; Zaheer, R.; Cook, S.R.; Booker, C.W.; Hendrick, S.; Alexander, T.W.; McAllister, T.A. Pathogens of bovine respiratory disease in North American feedlots conferring multidrug resistance via integrative conjugative elements. J. Clin. Microbiol. 2014, 52, 438–448. [Google Scholar] [CrossRef] [Green Version]
  11. Noyes, N.R.; Benedict, K.M.; Gow, S.P.; Booker, C.W.; Hannon, S.J.; McAllister, T.A.; Morley, P.S. Mannheimia haemolytica in feedlot cattle: Prevalence of recovery and associations with antimicrobial use, resistance, and health outcomes. J. Vet. Intern. Med. 2015, 29, 705–713. [Google Scholar] [CrossRef]
  12. Timsit, E.; Hallewell, J.; Booker, C.; Tison, N.; Amat, S.; Alexander, T.W. Prevalence and antimicrobial susceptibility of Mannheimia haemolytica, Pasteurella multocida, and Histophilus somni isolated from the lower respiratory tract of healthy feedlot cattle and those diagnosed with bovine respiratory disease. Vet. Microbiol. 2017, 208, 118–125. [Google Scholar] [CrossRef]
  13. Loy, J.D.; Leger, L.; Workman, A.M.; Clawson, M.L.; Bulut, E.; Wang, B. Development of a multiplex real-time PCR assay using two thermocycling platforms for detection of major bacterial pathogens associated with bovine respiratory disease complex from clinical samples. J. Vet. Diagn. Investig. 2018, 30, 837–847. [Google Scholar] [CrossRef] [Green Version]
  14. McLain, J.E.; Cytryn, E.; Durso, L.M.; Young, S. Culture-based methods for detection of antibiotic resistance in agroecosystems: Advantages, challenges, and gaps in knowledge. J. Environ. Qual. 2016, 45, 432–440. [Google Scholar] [CrossRef] [PubMed]
  15. DeDonder, K.D.; Apley, M.D. A literature review of antimicrobial resistance in pathogens associated with bovine respiratory disease. Anim. Health Res. Rev. 2015, 16, 125–134. [Google Scholar] [CrossRef] [PubMed]
  16. Boolchandani, M.; D’Souza, A.W.; Dantas, G. Sequencing-based methods and resources to study antimicrobial resistance. Nat. Rev. Genet. 2019, 20, 356–370. [Google Scholar] [CrossRef] [PubMed]
  17. Dennis, E.J. Antimicrobial Used for Metaphylaxis and Implications for Product Diversification in the Animal Health Sector. NebGuide G2319 Animal Agriculture/Beef. 2019. Available online: https://extensionpublications.unl.edu/assets/pdf/g2319.pdf (accessed on 28 August 2020).
  18. USDA NAHMS. Antimicrobial Use and Stewardship on U.S. Feedlots, 2017. 2019. Available online: https://www.aphis.usda.gov/animal_health/nahms/amr/downloads/amu-feedlots_1.pdf (accessed on 28 August 2020).
  19. Clawson, M.L.; Murray, R.W.; Sweeney, M.T.; Apley, M.D.; DeDonder, K.D.; Capik, S.F.; Larson, R.L.; Lubbers, B.V.; White, B.J.; Kalbfleisch, T.S.; et al. Genomic signatures of Mannheimia haemolytica that associate with the lungs of cattle with respiratory disease, an integrative conjugative element, and antibiotic resistance genes. BMC Genom. 2016, 17, 982. [Google Scholar] [CrossRef] [Green Version]
  20. Livermore, D.M.; Winstanley, T.G.; Shannon, K.P. Interpretative reading: Recognizing the unusual and inferring resistance mechanisms from resistance phenotypes. J. Antimicrob. Chemother. 2001, 48 (Suppl. S1), 87–102. [Google Scholar] [CrossRef] [Green Version]
  21. Owen, J.R.; Noyes, N.; Young, A.E.; Prince, D.J.; Blanchard, P.C.; Lehenbauer, T.W.; Aly, S.S.; Davis, J.H.; O’Rourke, S.M.; Abdo, Z.; et al. Whole-genome sequencing and concordance between antimicrobial susceptibility genotypes and phenotypes of bacterial isolates associated with bovine respiratory disease. G3 Genes Genomes Genet. 2017, 7, 3059. [Google Scholar] [CrossRef] [Green Version]
  22. Ruegg, P.L.; Oliveira, L.; Jin, W.; Okwumabua, O. Phenotypic antimicrobial susceptibility and occurrence of selected resistance genes in gram-positive mastitis pathogens isolated from Wisconsin dairy cows. J. Dairy Sci. 2015, 98, 4521–4534. [Google Scholar] [CrossRef]
  23. Bortolaia, V.; Kaas, R.S.; Ruppe, E.; Roberts, M.C.; Schwarz, S.; Cattoir, V.; Philippon, A.; Allesoe, R.L.; Rebelo, A.R.; Florensa, A.F.; et al. ResFinder 4.0 for predictions of phenotypes from genotypes. J. Antimicrob. Chemother. 2020, 75, 3491–3500. [Google Scholar] [CrossRef]
  24. Hand, D.J. Measuring classifier performance: A coherent alternative to the area under the ROC curve. Mach. Learn. 2009, 77, 103–123. [Google Scholar] [CrossRef] [Green Version]
  25. Loy, J.D.; Clawson, M.L. Rapid typing of Mannheimia haemolytica major genotypes 1 and 2 using MALDI-TOF mass spectrometry. J. Microbiol. Methods 2017, 136, 30–33. [Google Scholar] [CrossRef] [Green Version]
  26. Narayanan, S.; Bates, H.; Confer, A.; Couger, B.; Ramachandran, A. Whole-genome sequence of multidrug-resistant Bibersteinia trehalosi strain OADDL-BT1. Microbiol. Resour. Announc. 2019, 8, e01690-18. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. CLSI. VET01-A4 Performance Standards for Antimicrobial Disk and Dilution Susceptibility Test for Bacteria Isolated from Animals: Approved Standard-Fourth Edition; Watts, J.E.A., Ed.; CLSA: Wayne, PA, USA, 2013; Volume 33. [Google Scholar]
  28. CLSI. VET01-S2 Performance Standards for Antimicrobial Disk and Dilution Susceptibility Tests for Bacteria Isolated from Animals; Second Information Supplement; Watts, J.E.A., Ed.; CLSA: Wayne, PA, USA, 2013; Volume 33. [Google Scholar]
  29. Team, R.C. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2019. [Google Scholar]
  30. Robin, X.; Turck, N.; Hainard, A.; Tiberti, N.; Lisacek, F.; Sanchez, J.C.; Muller, M. pROC: An open-source package for R and S+ to analyze and compare ROC curves. BMC Bioinform. 2011, 12, 77. [Google Scholar] [CrossRef] [PubMed]
  31. Zweig, M.H.; Campbell, G. Receiver-operating characteristic (ROC) plots: A fundamental evaluation tool in clinical medicine. Clin. Chem. 1993, 39, 561–577. [Google Scholar] [CrossRef]
  32. Melo, F. Area under the ROC Curve. In Encyclopedia of Systems Biology; Dubitzky, W.W.O., Cho, K.H., Yokota, H., Eds.; Springer: New York, NY, USA, 2013. [Google Scholar]
  33. Dohoo, I.R.; Martin, W.; Stryhn, H.E. Veterinary Epidemiologic Research, 2nd ed.; VER Inc.: Charlottetown, PEI, Canada, 2009; pp. 97–98. [Google Scholar]
  34. Mabmud, S.M. Cohen’s Kappa. In Encyclopedia of Research Design; Salkind, N.J., Ed.; SAGE Publications, Inc.: Thousand Oaks, CA, USA, 2010. [Google Scholar] [CrossRef]
  35. Steinberg, D.M.; Fine, J.; Chappell, R. Sample size for positive and negative predictive value in diagnostic research using case-control designs. Biostatistics 2009, 10, 94–105. [Google Scholar] [CrossRef] [Green Version]
  36. Riley, R.D.; Ahmed, I.; Debray, T.P.A.; Willis, B.H.; Noordzij, J.P.; Higgins, J.P.; Deeks, J.J. Summarising and validating test accuracy results across multiple studies for use in clinical practice. Stat. Med. 2015, 34, 2081–2103. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Prevalence of major bovine respiratory disease pathogens tested in this study by sample type.
Figure 1. Prevalence of major bovine respiratory disease pathogens tested in this study by sample type.
Pathogens 10 00064 g001
Figure 2. Consistent distributions of cycle threshold values quantified by the multiplex qPCR assay among three genes conferring resistance to macrolides.
Figure 2. Consistent distributions of cycle threshold values quantified by the multiplex qPCR assay among three genes conferring resistance to macrolides.
Pathogens 10 00064 g002
Figure 3. Receiver operating characteristic curves for determining optimal cycle threshold (Ct) cutoff values for lung samples (a) and nasal samples (b). Determined optimal Ct cutoff values with corresponding sensitivity and specificity were provided for predicting the presence of bovine respiratory disease pathogens resistant to tetracycline and macrolide drugs, respectively. Combined macrolide refers to resistance to either tilmicosin or tulathromycin.
Figure 3. Receiver operating characteristic curves for determining optimal cycle threshold (Ct) cutoff values for lung samples (a) and nasal samples (b). Determined optimal Ct cutoff values with corresponding sensitivity and specificity were provided for predicting the presence of bovine respiratory disease pathogens resistant to tetracycline and macrolide drugs, respectively. Combined macrolide refers to resistance to either tilmicosin or tulathromycin.
Pathogens 10 00064 g003aPathogens 10 00064 g003b
Table 1. Sensitivity and limit of detection in colony forming units per reaction (CFU/rxn) using a Mannheimia haemolytica reference strain (1621) confirmed to have targets by whole-genome sequencing.
Table 1. Sensitivity and limit of detection in colony forming units per reaction (CFU/rxn) using a Mannheimia haemolytica reference strain (1621) confirmed to have targets by whole-genome sequencing.
TargetCFU/rxnCt 1CD 2RE 3
ICETetR320026.78
ICETetR32030.23
ICETetR3233.53
ICETetR3.237.340.9820.95
erm42320026.32
erm4232029.85
erm423232.95
erm423.236.170.9861.00
mph E320025.61
mph E32029.05
mph E3232.30
mph E3.236.560.9850.88
msr E320025.44
msr E32028.82
msrE3232.04
msr E3.235.700.9850.98
1 Ct = threshold cycle, 2 CD = correlation of determination, R2, 3 RE = reaction efficiency (10^(−1/slope)−1).
Table 2. Sample distribution by clinical sample type and pathogen.
Table 2. Sample distribution by clinical sample type and pathogen.
Sample Type 1Total Sample Size 2Occurrence of BRD Pathogens
No. of Positive Samples (Prevalence, 95% CI)
M. haemolyticaP. multocidaH. somniPositive for at Least One BRD Pathogen
Lung sample297191
(64.3%,
58.7–69.5%)
96
(32.3%,
27.3–37.8%)
93
(31.3%,
26.3–36.8%)
296
(99.7%,
98.1–100.0%)
Nasal sample11172
(64.9%,
55.6–73.1%)
80
(72.1%,
63.1–79.6%)
29
(26.1%,
18.9–35.0%)
108
(97.3%,
92.4–99.1%)
Others—skin/liver41
(25.0%,
4.6–69.9%)
0
(0.0%,
0.0%–49.0%)
1
(25.0%,
4.6–69.9%)
2
(50%,
15.0–85.0%)
Missing41
(25.0%,
4.6–69.9%)
1
(25.0%,
4.6–69.9%)
1
(25.0%,
4.6–69.9%)
1
(25.0%,
4.6–69.9%)
Total416265
(63.7%,
59.0–68.2%)
177
(42.5%,
37.9–47.3%)
123
(29.6%,
25.4–34.1%)
407
(97.8%,
95.9–98.8%)
1 Lung samples refer to lung tissues and abdominal fluid clinical samples; nasal samples are collected as nasal swabs; other samples are either skin or liver tissues; missing samples are those lacking records about tissue types. 2 The distribution of the samples by the pathogens is not exclusive, as more than one pathogen can be isolated from a single sample.
Table 3. Summary of minimum inhibitory concentration test results by sample type, pathogen, and antimicrobial.
Table 3. Summary of minimum inhibitory concentration test results by sample type, pathogen, and antimicrobial.
PathogenSampleClassAntibioticsNo. of Total SamplesNo. (Percentage) of Samples with 1
RIR+IS
M. haemolyticaLungTetracyclineOxytetracycline 19178 (40.8%)3 (1.6%)81 (42.4%)110 (57.6%)
MacrolideTilmicosin19169 (36.1%)5 (2.6%)74 (38.7%)117 (61.3%)
MacrolideTulathromycin19164 (33.5%)3 (1.6%)67 (35.1%)124 (64.9%)
MacrolideTilmicosin or tulathromycin19175 (39.3%)0 (0%)75 (39.3%)116 (60.7%)
NasalTetracyclineOxytetracycline 723 (4.2%)2 (2.7%)5 (6.9%)67 (93.1%)
MacrolideTilmicosin722 (2.8%)0 (0%)2 (2.8%)70 (97.2%)
MacrolideTulathromycin726 (8.3%)0 (0%)6 (8.3%)66 (91.7%)
MacrolideTilmicosin or tulathromycin726 (8.3%)0 (0%)6 (8.3%)66 (91.7%)
P. multocidaLungTetracyclineOxytetracycline 9628 (29.2%)2 (2.1%)30 (31.3%)66 (68.7%)
MacrolideTilmicosin9515 (15.8%)1 (1.0%)16 (16.8%)79 (83.2%)
MacrolideTulathromycin969 (9.4%)1 (1.1%)10 (10.5%)86 (89.5%)
MacrolideTilmicosin or tulathromycin9515 (15.8%)1 (1.0%)16 (16.8%)79 (83.2%)
NasalTetracyclineOxytetracycline 7815 (19.2%)0 (0%)15 (19.2%)63 (80.8%)
MacrolideTilmicosin806 (7.5%)0 (0%)6 (7.5%)74 (92.5%)
MacrolideTulathromycin781 (1.3%)3 (3.8%)4 (5.1%)74 (94.9%)
MacrolideTilmicosin or tulathromycin807 (8.8%)0 (0%)7 (8.8%)73 (91.2%)
H. somniLungTetracyclineOxytetracycline 9341 (44.1%)8 (8.6%)49 (52.7%)44 (47.3%)
MacrolideTilmicosin9320 (21.5%)1 (1.1%)21 (22.6%)72 (77.4%)
MacrolideTulathromycin9320 (21.5%)8 (8.6%)28 (30.1%)65 (69.9%)
MacrolideTilmicosin or tulathromycin9326 (28.0%)6 (6.4%)32 (34.4%)61 (65.6%)
NasalTetracyclineOxytetracycline 2913 (44.8%)0 (0%)13 (44.8%)16 (55.2%)
MacrolideTilmicosin292 (6.9%)1 (3.4%)3 (10.3%)26 (89.7%)
MacrolideTulathromycin285 (17.9%)2 (7.1%)7 (25.0%)21 (75.0%)
MacrolideTilmicosin or tulathromycin286 (21.4%)3 (10.7%)9 (32.1%)19 (67.9%)
At least one BRD pathogenLungTetracyclineOxytetracycline 296124 (41.9%)8 (2.7%)132 (44.6%)164 (55.4%)
MacrolideTilmicosin29595 (32.2%)6 (2.0%)101 (34.2%)194 (65.8%)
MacrolideTulathromycin29687 (29.4%)8 (2.7%)95 (32.1%)197 (67.9%)
MacrolideTilmicosin or tulathromycin295105 (35.6%)4 (1.3%)109 (36.9%)182 (63.1%)
NasalTetracyclineOxytetracycline 10823 (21.3%)2 (1.8%)25 (23.1%)83 (76.9%)
MacrolideTilmicosin10811 (10.2%)0 (0.0%)11 (10.2%)97 (89.8%)
MacrolideTulathromycin10811 (10.2%)5 (4.6%)16 (14.8%)92 (85.2%)
MacrolideTilmicosin or tulathromycin10819 (17.6%)2 (1.8%)21 (19.4%)87 (80.6%)
1 S, R, I denotes susceptible, resistant, and intermediate resistant to the drug classified based on MIC test and CLSI breakpoints, and the categories of R and I are combined into a new category of “resistant”, or R+I.
Table 4. Phenotypic classification of antimicrobial resistance based on minimum inhibitory concentration measurements by antibiotics tested in this study.
Table 4. Phenotypic classification of antimicrobial resistance based on minimum inhibitory concentration measurements by antibiotics tested in this study.
ClassAntibioticsAntimicrobial Resistance Classification
Susceptible (S)Intermediate (I)Resistant (R)“Resistant” (R+I)
TetracyclineOxytetracycline≤2>2 and ≤8>8>2
MacrolideTilmicosin≤8>8 and ≤32>32>8
MacrolideTulathromycin≤16>16 and ≤32>32>16
Table 5. Optimal cycle threshold (Ct) cutoff value by sample type and antibiotics.
Table 5. Optimal cycle threshold (Ct) cutoff value by sample type and antibiotics.
BRD PathogenSampleClassAntibioticsNo. of Total SamplesNo. of Samples with 1Optimal Cycle Threshold (Ct)Se 2
(%)
Sp 2
(%)
AUC 2 (%)Kappa (κ)Prevalence
(%)
PPV 2
(%)
NPV 2
(%)
R+IS
M. haemolyticaLungTetracyclineOxytetracycline 1918111031.0077.7895.4589.990.7542.4192.6585.37
MacrolideTilmicosin1917411733.0474.3293.1686.140.6938.7487.3085.16
MacrolideTulathromycin1916712432.8979.1092.7488.940.7335.0885.4889.15
MacrolideTilmicosin or tulathromycin1917511633.0473.3393.1085.450.6839.2787.3084.38
NasalTetracyclineOxytetracycline 7256732.26100.0079.1092.540.346.9421.8398.35
MacrolideTilmicosin7227021.42100.00100.00100.001.002.7852.2599.17
MacrolideTulathromycin7266630.7350.0086.3663.890.258.3325.0095.00
MacrolideTilmicosin or tulathromycin7266630.7350.0086.3663.890.258.3325.0095.00
P. multocidaLungTetracyclineOxytetracycline 96306636.1083.3386.3686.520.6731.2573.5391.94
MacrolideTilmicosin95167932.9156.2589.8772.230.4516.8452.9491.03
MacrolideTulathromycin96108632.9180.0089.5384.940.5310.4247.0697.47
MacrolideTilmicosin or tulathromycin95167932.9156.2589.8772.230.4516.8452.9491.03
NasalTetracyclineOxytetracycline 78156329.3566.6792.0688.040.5919.2366.6792.06
MacrolideTilmicosin8067436.4066.6762.1659.010.107.5012.5095.83
MacrolideTulathromycin7847431.4750.0078.3858.780.115.1311.1196.67
MacrolideTilmicosin or tulathromycin8077332.2242.8676.7156.360.118.7515.0093.33
H. somniLungTetracyclineOxytetracycline 93494436.2881.6372.7375.720.5552.6976.9278.05
MacrolideTilmicosin93217233.0861.9079.1767.290.3722.5846.4387.69
MacrolideTulathromycin93286531.6760.7187.6975.360.5030.1168.0083.82
MacrolideTilmicosin or tulathromycin93326133.0862.5086.8976.950.5134.4171.4381.54
NasalTetracyclineOxytetracycline 29131632.8592.3156.2562.500.4744.8363.1690.00
MacrolideTilmicosin2932630.8866.6765.3852.560.1510.3418.1894.44
MacrolideTulathromycin2872127.19100.0038.1068.03-0.3625.0010.5467.85
MacrolideTilmicosin or tulathromycin2891926.83100.0036.8466.67−0.3932.1412.4559.49
At least one BRD pathogenLungTetracyclineOxytetracycline 29613216436.0681.0682.9385.290.6444.5979.2684.47
MacrolideTilmicosin29510119433.0869.3189.6981.250.6134.2477.7884.88
MacrolideTulathromycin2969520132.8972.6390.5584.290.6432.0978.4187.50
MacrolideTilmicosin or tulathromycin29510918633.0867.8991.4081.940.6236.9582.2282.93
NasalTetracyclineOxytetracycline 108258332.8188.0079.5290.070.5623.1556.4195.65
MacrolideTilmicosin108119731.8263.6482.4771.420.3010.1929.1795.24
MacrolideTulathromycin108169231.4743.7583.7057.540.2414.8131.8289.53
MacrolideTilmicosin or tulathromycin108218731.8242.8682.7657.850.2419.4437.5085.71
1 S, R, I denotes susceptible, resistant, and intermediate resistant to the drug classified based on MIC test and CLSI breakpoints, and the categories of R and I are combined into a new category of “resistant”, or R+I. 2 Se = sensitivity; Sp = specificity; AUC = area under the curve; PPV = positive predictive value; NPV = negative predictive value.
Table 6. Optimal cycle threshold (Ct) cutoff value for lung and nasal samples.
Table 6. Optimal cycle threshold (Ct) cutoff value for lung and nasal samples.
ClassAntibioticsLung SampleNasal Sample
Optimal Cycle Threshold (Ct)Kappa (κ)Optimal Cycle Threshold (Ct)Kappa (κ)
TetracyclineOxytetracycline 35.660.6133.270.49
MacrolideTilmicosin33.120.61--
MacrolideTulathromycin32.640.63--
MacrolideTilmicosin or tulathromycin33.120.62--
Table 7. Oligonucleotide sequences for primers and probes used in the assay.
Table 7. Oligonucleotide sequences for primers and probes used in the assay.
TargetPrimer/ProbeSequence (5′-3′)SizeReference
tetRICEtetR-F TTTGGCTTTCTTGATGCTCTTG 71This paper
ICEtetR-R GTGATGCTGGGTTTAGTCTATCT
ICEtetR-P (CY5/TAO-IAB-RQ) CGCAATAGAGCTTAATGCATACACGGC
erm42erm42-FGCCGTTAATGCTATTGAGTTCG105This paper
erm42-RCGGCTTCAATAATAGACACATTTGA
erm42-P (FAM/ZEN-IAB-FQ)AGTGTATTGGCTGATAAGTTGAGCCATGA
msrEmsrE-FGGGTGGTTACTCGGATTACTTG88This paper
msrE-R CTCCCGTTCCTTCATCATCAG
msrE-P (Texas Red/IAB-RQ) AGCGACAACACCAAGCCGTAGAAT
mphEmphE-F TTGGAAACCCGCTACAGAAA 113This paper
mphE-R GCTCCATCCTTTGAAGCTAGT
mphE-P (JOE/ZEN-IAB-FQ) TGATGTTCTATGGGCAGATTTCACCCA
Table 8. Schematic representation of the outcomes of the multiplex qPCR assay detecting antimicrobial resistance compared to the gold standard.
Table 8. Schematic representation of the outcomes of the multiplex qPCR assay detecting antimicrobial resistance compared to the gold standard.
Multiplex qPCR AssayCulture-Based Gold-Standard Test
“R” 1S 2
“R” 1True Positive (TP)False Positive (FP)PPV 3 = TP/(TP + FP)
S 2False Negative (FN)True Negative (TN)NPV 3 = T N/(TN + FN)
Se 3 = TP/(TP + FN)Sp 3 = TN/(TN + FP)
1 “R” indicated by the gold-standard test refers to samples containing pathogens resistant (R) or intermediate resistant (I) to antibiotics tested (“resistant” or R+I); hence classification given qPCR results based on the comparison to the gold-standard test also refers to the redefined “resistant” in this study. 2 S indicated by both the gold-standard test and qPCR assay refers to samples containing pathogens susceptible to antibiotics tested. 3 Se = sensitivity; Sp = specificity; PPV = positive predictive value; and NPV = negative predictive value.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Dutta, E.; Loy, J.D.; Deal, C.A.; Wynn, E.L.; Clawson, M.L.; Clarke, J.; Wang, B. Development of a Multiplex Real-Time PCR Assay for Predicting Macrolide and Tetracycline Resistance Associated with Bacterial Pathogens of Bovine Respiratory Disease. Pathogens 2021, 10, 64. https://0-doi-org.brum.beds.ac.uk/10.3390/pathogens10010064

AMA Style

Dutta E, Loy JD, Deal CA, Wynn EL, Clawson ML, Clarke J, Wang B. Development of a Multiplex Real-Time PCR Assay for Predicting Macrolide and Tetracycline Resistance Associated with Bacterial Pathogens of Bovine Respiratory Disease. Pathogens. 2021; 10(1):64. https://0-doi-org.brum.beds.ac.uk/10.3390/pathogens10010064

Chicago/Turabian Style

Dutta, Enakshy, John Dustin Loy, Caitlyn A. Deal, Emily L. Wynn, Michael L. Clawson, Jennifer Clarke, and Bing Wang. 2021. "Development of a Multiplex Real-Time PCR Assay for Predicting Macrolide and Tetracycline Resistance Associated with Bacterial Pathogens of Bovine Respiratory Disease" Pathogens 10, no. 1: 64. https://0-doi-org.brum.beds.ac.uk/10.3390/pathogens10010064

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop