Integrated Bacterial Identification and Antimicrobial ... - ACS Publications

Oct 13, 2017 - The Sidney Kimmel Comprehensive Cancer Center at Johns Hopkins, Baltimore, Maryland 21287, United States. •S Supporting Information. ...
0 downloads 13 Views 1MB Size
Subscriber access provided by Gothenburg University Library

Article

Integrated Bacterial Identification and Antimicrobial Susceptibility Testing using PCR and High-Resolution Melt Pornpat Athamanolap, Kuangwen Hsieh, Liben Chen, Samuel Yang, and Tza-Huei Wang Anal. Chem., Just Accepted Manuscript • DOI: 10.1021/acs.analchem.7b02809 • Publication Date (Web): 13 Oct 2017 Downloaded from http://pubs.acs.org on October 15, 2017

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a free service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are accessible to all readers and citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

Analytical Chemistry is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 15

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Analytical Chemistry

Integrated Bacterial Identification and Antimicrobial Susceptibility Testing using PCR and HighResolution Melt Pornpat Athamanolap1, Kuangwen Hsieh2, Liben Chen2, Samuel Yang3, Tza-Huei Wang1-2, 4* 1

Department of Biomedical Engineering, Johns Hopkins School of Medicine, Baltimore, MD 21205, USA, Department of Mechanical Engineering, Johns Hopkins University, Baltimore, MD 21218, USA, 3 Department of Emergency Medicine, Stanford University, Stanford, CA 94305, USA 4 Johns Hopkins Institute for NanoBioTechnology, Johns Hopkins University, Baltimore, MD 21218, USA and The Sidney Kimmel Comprehensive Cancer Center at Johns Hopkins, Baltimore, MD 21287, USA *Corresponding author: [email protected], Phone: (+1)410-516-7086, Fax: (+1)410-516-7254 2

Abstract Accurate and timely diagnostics are critical for managing bacterial infections. The current gold standard, culture-based diagnostics can provide clinicians with comprehensive diagnostic information including bacterial identity and antimicrobial susceptibility, but they often require several days of turnaround time, which leads to compromised clinical outcome and promotes the spread of antibiotic resistance. Nucleic acid amplification tests such as PCR have significantly accelerated the detection of specific bacteria, but generally lack the capacities for broad-based bacterial identification or antimicrobial susceptibility testing. Here we report an integrated assay based on PCR and high resolution melt (HRM) for rapid diagnosis for bacterial infections. In our assay, we measure bacterial growth in the presence or absence of certain antibiotics with real-time quantitative PCR or digital PCR to determine antimicrobial susceptibility. In addition, we use HRM and a machine learning algorithm to identify bacterial species based on melt-curve profiles of the 16S rRNA gene in an automated fashion. As a demonstration, we correctly identified the bacterial species and their antimicrobial susceptibility profiles for multiple unknown samples in blinded tests within ~6.5 hours.

Introduction Bacterial infections can result in many serious or life-threatening complications such as sepsis and urinary tract infections (UTI). Successful disease management hinges on timely treatment with proper antibiotics, which requires rapid diagnostics that can provide both bacterial identification (ID) and antimicrobial susceptibility testing (AST) profiles1–3. Unfortunately, the current gold standard for bacteria ID and AST remain heavily reliant on conventional culture method. This traditional method begins with Gram staining, followed by culturing for ID; culturing-based AST is performed either in parallel with or after the ID process. This culturing workflow involves tedious processes and requires at least 48 hours to produce final results4. Consequently, although bacterial culture can provide accurate ID and comprehensive AST, the lengthy turnaround time still leads to empirical use of broad-spectrum antibiotics, which has fueled the spread of antibiotics resistance5,6. Therefore, faster diagnostics that are capable of acquiring both ID and AST remain highly desirable7. Toward developing faster diagnostic methods, significant efforts have focused on detecting specific bacterial genes and proteins with mass spectrometry8,9, Raman spectroscopy10,11, immunoassays12 and nucleic acid amplification tests (NAATs)13–28. Among these, NAATs represent the most mature approach and offer high sensitivity and significantly shortened turnaround time for ID, exemplified by FDA-approved

ACS Paragon Plus Environment

Analytical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

NAAT platforms such as Cephid’s GeneXpert® for detecting a specific pathogen (e.g., MRSA)29,30 and bioMérieux’s BioFire FilmArray® for detecting a panel of pathogens (e.g., respiratory panel) via nested multiplex PCR and amplicon melting analysis for specificity confirmation31. Notably, these NAATs either use specific primers or probes for detecting a specific bacterial strain in each test, or use multiple primer pairs to detect several bacterial species, which render the assay inevitably more complex to design, perform, and optimize32. Consequently, expanding the panel of pathogens of interests for broader-scale identification would be challenging. As an alternative in enhancing the ID capability of NAATs, we and others have been developing PCR assays that can achieve broad-based bacteria detection33–42 by targeting the 16S ribosomal RNA (rRNA) gene that is present in all prokaryotic cells. The 16S rRNA gene has conserved regions that are consistent among bacteria and variable regions that are specific to the bacteria’s genus and species43–45. As such, a pair of “universal” primers can be designed to hybridize to conserved regions that flank hypervariable regions, allowing the hypervariable regions to be PCR amplified. PCR products can then be analyzed via various techniques such as sequencing, microarrays, mass spectrometry, and high-resolution melt (HRM) to identify the bacteria. Among these techniques, HRM can be performed directly after PCR without opening reaction tubes and transferring PCR products, and thus presents a simpler, faster, more costeffective, and less contamination-prone option. We have therefore employed HRM for identifying bacteria by differentiating their unique, species-specific melting profile42,46,47. Moreover, we have developed a nested one-versus-one support vector machine learning algorithm (ovoSVM)48,49 to differentiate the unique melting profiles of 37 clinically relevant bacterial species in an automated fashion. Nevertheless, although the PCR-HRM method offers an effective means to broad-based bacteria identification, it does not provide AST information and thus must be integrated with an AST method. Genotypic approaches offer some capacity for determining the antimicrobial susceptibility profile of the causative bacteria by detecting the presence of resistant genes or mutations, and this genetic information can be readily detected by PCR. However, only a few genes or mutations have been found to correlate with known mechanisms of resistance50. Moreover, because genotypic approaches rely on detecting known resistant genes and mutations, they would not be able to detect new variants of antibiotic resistant mechanisms51,52 that continue to emerge. As such, phenotypic approaches in which bacteria are grown directly with antibiotics remain the most broadly applicable way for obtaining AST profiles. Provided that our PCR-HRM assay and more broadly genetic detection assays all require extracting and amplifying DNA from lysed cells, performing growth-based phenotypic AST after these assays is impossible. This has necessitated the development of a new workflow that can seamlessly integrate these two inherently incompatible assays. Herein, we report a simple yet comprehensive bacterial diagnostic assay by effectively coupling our broad-based bacterial ID with growth-based phenotypic AST. This is achieved by simply incubating bacteria in the presence and absence of an antibiotic before performing real-time quantitative PCR (qPCR) and HRM. In doing so, bacteria of interest are still identified via HRM, which examines the profiles of melt curves; the antimicrobial susceptibility profiles are determined by qPCR, which quantifies the difference in bacterial DNA quantity from the two samples to assess the status of growth. Compared to earlier studies that coupled growth-based phenotypic AST with genetic amplification methods for bacteria detection53,54, our current work offers true capacity for bacterial species identification from a broader panel of bacteria. As a demonstration, we use our assay to correctly identify the bacteria species and determine the antimicrobial susceptibility of multiple unknown bacterial samples in blinded tests. Finally, as an

ACS Paragon Plus Environment

Page 2 of 15

Page 3 of 15

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Analytical Chemistry

additional benefit, we can also implement an absolute quantitative digital PCR in our workflow to further enhance the diagnosis of polymicrobial infections.

Experimental Section Bacterial Strains, Antibiotic, and Storage A total of six bacteria strains were purchased from the American Type Culture Collection (ATCC; Manassas, VA): Escherichia coli (E. coli) ATCC 25922, Staphylococcus aureus (S. aureus) ATCC 29213, Enterococcus faecalis (E. faecalis) ATCC 29212, Proteus mirabilis (P. mirabilis) ATCC 12453, E. coli BAA2471, and S. aureus BAA-44. Of these, E. coli BAA-2471 and S. aureus BAA-44 are multi-drug resistant. The six bacteria strains were separately plated on tryptic soy agar plates (TSA; BD Diagnostics, Sparks, MD) at 37 °C overnight. An isolated colony from each plate was grown at 37 °C to log phase in tryptic soy broth (TSB; BD Diagnostics, Sparks, MD). The bacteria were then mixed with 20% glycerol (v/v; Sigma-Aldrich, St. Louis, MO), aliquoted, and frozen at -80 °C until use. After the aliquoted stocks were completely frozen, one aliquot for each strain was thawed and counted via plating in TSA, which accounts the potential effect of adding glycerol and freezing bacteria and thus provides a more accurate estimate of the stock concentrations. Gentamicin (10 mg/mL) was purchased from Quality Biological Inc. (Gaithersburg, MD), stored at 4 °C, and used without further purification. Integrated ID-AST Assay Samples containing either gentamicin or no gentamicin with bacteria in culture broth were freshly prepared before each experiment. Gentamicin was diluted to desired concentration with molecular grade water (Quality Biological Inc., Gaithersburg, MD). Aliquots of bacteria stock were thawed in a 37 °C water bath for ~5 minutes and were subsequently diluted with Mueller-Hinton II broth (MHII; BD Diagnostics, Sparks, MD) to 5 × 106 CFU/mL. Finally, 5 µL of diluted bacteria and 5 µL of gentamicin or water were added to 40 µL MHII broth, achieving final concentrations of 5 × 105 CFU/mL, which adheres with the CLSI standard cell concentration for AST. Subsequently, to initiate the integrated ID-AST assay, all samples were briefly cultured in an incubator at 37˚C for 2 to 4 hours. Cultured bacteria samples were then immediately aliquoted into a PCR mix to perform qPCR and HRM. The optimal conditions of PCR mixture contained in 25 l consisted of 2 l cultured bacteria sample, 1X Gold Buffer (Thermo Fisher Scientific, Waltham, MA), 1X Evagreen dye (Biotium, Freemont, CA), 3.5 mM MgCl2 (Thermo Fisher Scientific, Waltham, MA), 200 M of each deoxynucleotide triphosphate (dNTP) (Thermo Fisher Scientific, Waltham, MA), 0.3 M forward primer 5’-GYGGCGNACGGGTGAGTAA-3’ (Integrated DNA Technologies, Coralville, IA), 0.3 M reverse primer 5’AGCTGACGACANCCATGCA-3’ (Integrated DNA Technologies, Coralville, IA), 0.025 U/l of Amplitaq Gold LD (Thermo Fisher Scientific, Waltham, MA), and Ultra Pure PCR water (Quality Biological Inc., Gaithersburg, MD). qPCR and HRM were performed using a CFX96 Touch Real-time PCR Detection System (Bio-Rad, Hercules, CA) and analyzed using CFX Manager software. The qPCR cycling conditions were as follows: 10 min at 95˚C, followed by 50 cycles of (95˚C for 15 s, 65˚C for 15 s and a 60-s extension step at 72˚C). HRM was performed immediately after amplification with an extra extension step at 72˚C for 7 mins by temperature range of 75˚C - 95˚C at 0.2˚C increments and a 2-s hold time before each measurement.

ACS Paragon Plus Environment

Analytical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Machine Learning Algorithm and Melt Curve Database New experimentally generated melt curves were exported from Bio-Rad CFX software and imported into an in-house Matlab program, where the melt curves were first translated with a translation function and then compared with the melt curve database created in our previous study49 to perform bacterial identification. Specifically, the nested ovoSVM algorithm was used to classify bacterial species against the stored melt curves database, which contains 4 replicates of 37 organisms. Each unknown melt curve was input into multiple SVM classifiers that distinguished between all combinations of two organisms out of 37 organisms. By using all negative derivative RFUs at each temperature point as classification features, each SVM classifier output which species the input curve was closer to. The number of times that were classified to each species were counted. Finally, an unknown curve was identified to an organism that had the highest score when comparing to all other organisms. Antimicrobial Susceptibility Testing through qPCR Real-time amplification results were obtained via Bio-Rad CFX software. Cq values were determined using Regression mode with Baseline Subtracted Curve Fit. Cq was calculated based on the differences between Cqs from sample incubated with and without antibiotic [Cq_drug – Cq_no drug]. Polymicrobial Digital PCR-HRM Experiment S. aureus and E. faecalis were mixed together at ~106 CFU/mL and ~2x104 CFU/mL respectively to mimic a polymicrobial infection. The experiment was performed in triplicate by first splitting each mixed sample into two reactions to incubate with and without 32 µg/mL of gentamicin (Quality Biological Inc., Gaithersburg, MD) for two hours in MHII broth. Bacterial DNA was then extracted using QuickExtract Bacterial DNA Extraction Kit (Epicentre, Madison, WI) according to manufacturer’s protocol. In order to analyze individual melt curves that were derived from single bacterial species, we incorporated our previously developed technique, universal digital HRM (U-dHRM)46. To achieve either zero or only one bacterial DNA in each well, the extracted DNA was diluted 60000-fold before mixing with PCR mastermix containing similar condition as mentioned previously with the total final volume of 10 l for each reaction for the whole 96-well microtiter plate. Subsequent qPCR and HRM with similar cycling condition as mentioned above was performed using a CFX96 Touch Real-time PCR Detection System (Bio-Rad, Hercules, CA). Output melt curves generated by CFX software were entered into our in-house implemented Matlab program for bacterial identification with DovoSVM. The amount of DNA molecule of each species was determined by counting the number of positive reactions that were identified as the target species; S. aureus or E. faecalis. Antimicrobial susceptibility could then be defined based on different amount of detected DNA from broth with and without gentamicin.

Results and Discussion Assay Overview Our approach integrates qPCR and HRM with phenotypic AST to achieve bacterial ID/AST in one assay (Fig. 1). Samples containing unknown bacteria are first divided into two portions, where one portion is incubated in culture media with the antibiotic condition of interest and the other is incubated in culture media without any antibiotics as the control. Due to the high sensitivity of PCR, both samples only need

ACS Paragon Plus Environment

Page 4 of 15

Page 5 of 15

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Analytical Chemistry

Figure 1. Overview of our integrated bacterial ID/AST assay. The procedure starts with growth of a bacteria sample in MHII broth with and without antibiotics (e.g., gentamicin). Both cultivated samples are then added into PCR mixtures to perform broad-based qPCR-HRM on a real-time thermocycler. After qPCR-HRM, the species of the bacteria in the sample is identified by matching the newly generated melt curve to our species-specific melt curve database of 37 common bacterial species via our machine learning algorithm. AST is determined by measuring the differences of Cqs between reactions with and without antibiotics (∆Cq) obtained from qPCR results. That is, bacteria that do not grow in the presence of the antibiotic would have less bacterial DNA than the no-antibiotic control, resulting in a measurable ∆Cq.

to be briefly cultivated (e.g., 2 to 4 hours) before they are aliquoted into qPCR mixture containing EvaGreen, which to facilitates HRM after qPCR. qPCR is performed using a single pair of universal primers that amplify a ~1000-bp hypervariable region of the 16S rRNA gene to enable broad-based detection, as demonstrated previously49. After qPCR and HRM, bacterial species are identified by matching their distinct, species-specific melt curves to the melt curves of 37 clinically relevant bacteria in our previously established database using our ovoSVM algorithm as classifiers. Subsequently, the growth of bacteria is determined by quantifying the amount of bacterial DNA with qPCR based on the quantification cycle (Cq) value55. Antimicrobial susceptibility is then determined by comparing the difference in Cq (Cq) between the samples cultivated with and without the antibiotic. Susceptible strains do not grow in the presence of the antibiotic, resulting in less DNA than the no drug control and thus a noticeably large Cq. In contrast, resistant strains are expected to have similar growth profile regardless of the presence of antibiotic and therefore show minimal change in Cq between the two samples. Broad-based bacterial identification based on HRM using machine learning algorithm In our assay, accurate bacteria ID is predicated on the exact shape of HRM profiles (i.e., melt curves), and the shape must be independent from the instrument used to generate the melt curves. As the first step toward accurate ID, we have developed a translation function (Fig. 2, step 1) that allows us to match the melt curves of the same bacterial species generated from different instruments (i.e., Bio-Rad CFX platform and LightScanner). To do so, the translation function stretches and normalizes the melt curve from one instrument and fine-tunes these parameters until the translated melt curve matches that from the other instrument. In this work, we first established the initial translation parameters through matching the melt curves of E. coli, E. faecalis, and S. aureus that were generated from a Bio-Rad instrument to the melt curves of the same three species previously generated from a LightScanner instrument and stored in our database. We subsequently validated and finalized these translation parameters, as they allowed us to translate the Bio-Rad melt curve of P. mirabilis and fully match the stored LightScanner melt curve of P. mirabilis (Fig. S-1). Our work here demonstrates that if we were to use PCR instruments from other manufacturers and generate new melt curves in the future, these new melt curves would still be matched to existing melt curves of the same species in our database with similar translation functions with new,

ACS Paragon Plus Environment

Analytical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

fine-tuned parameters. For identifying a bacteria species, our identification workflow starts by obtaining its melt curve through qPCR-HRM (Fig. 2). This melt curve is then pre-processed and normalized in the translation function and then transferred into the ovoSVM classifiers (Fig. 2, step 2) to compare with our previously built melt curve database, which contains 16S melt curves (~1000 bp) of 37 clinically relevant bacteria species49 (Examples: E. coli – green, S. aureus – orange, P. mirabilis – yellow and E. faecalis – blue in Fig. 2). The output of the classifiers is alikeness scores indicating how close the unknown melt curve is to the melt curve of each species. The alikeness score is subsequently ranked and the unknown melt curve is finally identified as species with the highest alikeness score (Fig. 2, step 3).

Figure 2. Broad-based bacteria identification based on HRM analysis. A melt curve of an unknown bacterial species generated from a qPCR-HRM machine (Bio-Rad) is first translated into the same format as previously stored HRM database (LightScanner) with our translation function. The translated unknown melt curve is then used as the input in our machine learning algorithm, one-versus-one support vector machine (ovoSVM) and compared with the melt curves of the 37 clinically relevant bacterial species in our HRM database to perform bacterial species ID. The ovoSVM outputs alikeness score between the translated unknown melt curve and melt curves from each bacterial species. After ranking the score from all species, the unknown sample is identified as species with the highest score.

Our method is capable of identifying bacterial species even from different strains, which will be useful in clinically relevant settings where samples could be any strains. This is because the 16S primers were designed to target long amplicon covering multiple variable regions that are sufficient to differentiate bacterial species. As a result, melt curves that are generated from different strains but the same species will exhibit similar melt profiles, as demonstrated by our previous studies42,49. As an additional verification in this work, we tested two strains of E. coli (ATCC 25922 and BAA-2471), and two strains of S. aureus (ATCC 29213 and BAA-44), via our PCR-HRM assay and our ovoSVM algorithm, which were not performed in those previous studies. The melt curves from both strains of E. coli closely overlap each other, both having a single melting temperature peak at ~89.8˚C (Fig. S-2A and S-2B) and matching correctly with the E. coli melt curve in our database. Likewise, both strains of S. aureus also have the same melt curves with a small rise before the main melting temperature peak at ~89˚C (Fig. S-2C and S-2D) and are also correctly

ACS Paragon Plus Environment

Page 6 of 15

Page 7 of 15

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Analytical Chemistry

identified as S. aureus by our algorithm. These results clearly support that our method can accurately identify the bacteria species despite different strains. Antimicrobial Susceptibility Testing

Figure 3. Growth-based antimicrobial susceptibility testing (AST) through real-time quantitative PCR (qPCR). Real-time PCR curves of (A) gentamicin-susceptible E. coli, (B) gentamicin-resistant E. coli, (C) gentamicin-susceptible S. aureus, and (D) gentamicin-resistant S. aureus (MRSA) after incubating with and without gentamicin for two hours. In each case, the ∆Cq is calculated from the difference between the Cq of the sample without gentamicin (red) and the Cq of the sample with gentamicin (black). Significantly large ∆Cqs of 4.3 and 2.4 were obtained from susceptible E. coli and susceptible S. aureus, while small ∆Cqs of -0.2 and 0.8 were measured from gentamicin-resistant E. coli and MRSA respectively. Statistically significant ∆Cqs (p < 0.05) are observed between gentamicin-susceptible E. coli and gentamicinresistant E. coli (E), as well as between gentamicin-susceptible S. aureus and gentamicin-resistant S. aureus (F). These results provide further support that antibiotic-susceptible and resistant strains can be effectively differentiated by our method. . (G) Reference ∆Cqs for six different bacterial strains, including gentamicin-susceptible strains (green) of E. coli, S. aureus, E. faecalis, and P. mirabilis, as well as gentamicinresistant strains (orange) of E. coli and S. aureus (MRSA) are obtained after four hours of cultivation in culture broth with and without gentamicin. These reference ∆Cqs define the ∆Cq range for determining the susceptibility or resistance to gentamicin for future bacterial samples under the same cultivation conditions.

Bacterial strains that belong to the same species may exhibit different susceptibility to an antibiotic. This information can be effectively acquired by our integrated growth-based phenotypic test with qPCR analysis. In our method, we incubate bacteria in the presence and absence of an antibiotic and then we quantify the copy number of the 16S rRNA gene, which is highly correlated with bacterial growth (Fig. S3). As a specific example, we incubated gentamicin-susceptible and gentamicin-resistant strains of E. coli with no antibiotics or with 1 µg/mL gentamicin for two hours and then quantified bacterial growth via qPCR. For the gentamicin-susceptible E. coli, we measured a Cq of 27.86 ± 0.26 for the no-gentamicin control (Fig. 3A, red) and a significantly higher Cq of 32.19 ± 0.21 for the gentamicin-treated sample (Fig. 3A, black). The ∆Cq of 4.33 indicates that gentamicin effectively inhibited E. coli growth. In contrast, the measured Cq from gentamicin-resistant strain of E. coli for the no-gentamicin control was 25.10 ± 0.56 (Fig. 3C, red), which is similar to the Cq of 24.94 ± 0.66 measured for the gentamicin-treated sample (Fig. 3C, black). The small Cq of -0.16 indicates that gentamicin did not affect the growth of resistant strain of E. coli. The difference in Cqs (p = 0.0012) represents an antimicrobial susceptibility profile that is indeed statistically distinct for the two strains of E. coli (Fig. 3E). Likewise, our method correctly identified gentamicin susceptibility profile of S. aureus. For the gentamicin-susceptible S. aureus strain, we measured a Cq of 29.40 ± 0.17 from the no-gentamicin control

ACS Paragon Plus Environment

Analytical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

(Fig. 3B, red) and a Cq of 31.80 ± 0.60 for the gentamicin-treated sample (Fig. 3B, black), or a noticeable Cq of 2.40. In contrast to the susceptible strain, MRSA, a resistant S. aureus strain, exhibited a smaller Cq of 0.81, which was calculated from Cqs of 30.23 ± 0.16 and 31.04 ± 0.43 from the no-gentamicin control and the gentamicin-treated sample, respectively (Fig. 3D). The difference in Cqs between the two strains of S. aureus is statistically significant (p = 0. 0241), indicating that the gentamicin susceptibility profiles for these two strains are indeed distinct (Fig. 3F). We further expanded the susceptibility profiles to gentamicin to a total of six strains of bacteria, including four gentamicin-susceptible strains (E. coli, P. mirabilis, E. faecalis and S. aureus) and two gentamicin-resistant strains (E. coli and S. aureus). We characterized AST profiles to gentamicin of all the six strains by calculating their Cqs using our method with four hours of incubation time (Fig. 3G). Here, we lengthened the time of cultivation to four hours to ensure that our assay condition is consistent and generally applicable to multiple species and that the growth of slow-growing bacteria such as E. faecalis can be detectable via our qPCR protocol. Of note, we could have shortened the cultivation to as brief as one hour (Fig. S-3) if we were targeting only fast-growing bacteria such as E. coli. For the gentamicinsusceptible E. coli, P. mirabilis, E. faecalis and S. aureus, we measured relatively large Cqs of 8.73 ± 0.14, 7.28 ± 0.48, 2.41 ± 0.29, and 6.96 ± 0.51, respectively (n = 3). On the other hand, for the gentamicinresistant E. coli and S. aureus, we measured small Cqs of -0.06 ± 0.27 and 0.71 ± 0.33, respectively (n = 3). These AST profiles were established in order to be used as reference Cq ranges in determining the susceptibility or resistance of future bacterial samples. In addition, we note that the high sensitivity of qPCR also allows us to measure bacterial growth under a wide range of antibiotic concentrations and determine minimum inhibitory concentrations (MICs). For example, we determined the MIC for E. coli against gentamicin (a commonly used intravenous antibiotic for bacterial infections) as 0.25 µg/mL, which is consistent with the standard MIC and the traditional broth microdilution results (Fig. S-4). Blinded Tests To validate our assay and demonstrate its usefulness in a clinically relevant scenario, we performed blinded tests where we identified bacteria species and antimicrobial susceptibility without prior knowledge of the samples. We randomly picked three out of the six strains that we tested to be our unknown samples in performing blinded tests. In this blinded test scenario, the protocol simply started with splitting each sample into two portions to briefly incubate bacterial cells in MHII broth in the presence and absence of gentamicin followed by real-time PCR amplification and HRM analysis in triplicate. In result analysis, bacterial species identification was first determined by our translation function and HRM-based identification against the stored HRM database of 37 bacteria species using ovoSVM classifiers. Three samples were identified as E. faecalis, S. aureus and E. coli respectively. Afterwards, Cq of each testing sample was calculated from qPCR result. We obtained Cqs of 2.93, 0.20 and 8.79 for E. faecalis, S. aureus and E. coli respectively. Each sample was determined to be susceptible or resistant if the Cq was within ±2SD of the reference ΔCq of its own species (Fig. 3G). After comparing the Cq of each sample with our reference Cq ranges, final results of the three samples were E. faecalis - susceptible, S. aureus - resistant (MRSA), and E. coli - susceptible, respectively (Fig. 4). We correctly identified the species and the antimicrobial susceptibility of 3 out of 3 target bacteria within 6.5 hours, thus demonstrating our assay’s capacity in rapid and precise bacterial identification and antimicrobial susceptibility testing.

ACS Paragon Plus Environment

Page 8 of 15

Page 9 of 15

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Analytical Chemistry

Figure 4. Validation of our integrated workflow via blinded tests. Three samples with unknown bacteria species and unknown gentamicin susceptibility are all correctly diagnosed by our assay. After briefly incubating the samples with and without gentamicin followed by performing qPCR-HRM, bacterial species is identified based on HRM profile using ovoSVM and AST profile is determined based on comparing ∆Cq to reference values. Sample 1 is identified as E. faecalis (blue) and susceptible to gentamicin, Sample 2 is identified as S. aureus (orange) and resistant to gentamicin, and Sample 3 is identified as E. coli (green) and susceptible to gentamicin, respectively.

Digital PCR for Polymicrobial Infections The concept of briefly cultivating bacteria and then quantifying bacterial DNA can be further expanded by using digital PCR (dPCR) and digital HRM (dHRM). In fact, by enumerating individual molecule of bacterial DNA, dPCR and dHRM provide additional advantages. First, whereas multiple bacteria presented polymicrobial infections may generate a composite melt curve that is challenging to resolve via HRM, they are possibly decoupled and enumerated with digital analysis. By individually assessing a single molecule of bacterial DNA, each sub-population is detected and quantified, which also allows for distinguishing between infectious pathogen and contaminant by the bacterial load56. Second, counting individual molecules improves detection sensitivity, especially when analyzing a sample with only a few pathogens that may not be able to reliably detect the growth with Cq values from qPCR. Due to the better sensitivity, this digital concept can also reduce the incubation time in order to detect bacterial growth. The workflow of this digital assay starts with an analogous process with the previous described method, which begins by splitting the heterogeneous sample and incubating with and without an antibiotic. We then perform cell lysis and DNA extraction for both portions in parallel. The post-DNA extracted samples are then partitioned into separate 96-well plates to perform dPCR-dHRM. In result analysis, reactions containing bacterial DNA show a strong fluorescent signal and are considered as positive reactions, which subsequently species identified using HRM based classification. Based on the ID result, the total number of DNA molecules is counted for each species and AST profile of each bacterial species is then determined by comparing the number of DNA molecules between the bacteria grown in broth with and without an antibiotic as illustrated in Figure S-5.

ACS Paragon Plus Environment

Analytical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 10 of 15

Figure 5. Demonstration of using digital PCR-HRM (dPCR-HRM) for bacterial ID/AST from mixed bacterial samples. Samples containing gentamicin-susceptible S. aureus and gentamicin-susceptible E. faecalis were incubated in the presence and absence of gentamicin for 2 hours. Microtiter plates were used to perform dPCR-dHRM with extracted DNA from the mixed samples. By using our DovoSVM algorithm, 27 and 9 melt curves were respectively identified as S. aureus (orange) and E. faecalis (blue) from the no-gentamicin sample as shown in the color-coded melt curves (A) and in the 96-well plate (B). For the gentamicin-treated sample, 5 and 4 melt curves were identified as S. aureus (blue) and E. faecalis (orange), respectively, as shown in the color-coded melt curves (C) and in the 96-well plate (D). (E) Average absolute counts of each bacterial species under each condition are calculated from triplicate experiment. Based on the differences in absolute counts of the two bacterial species between growing with and without gentamicin, our digital bacterial ID/AST approach correctly identified that the mixed samples contained susceptible S. aureus and susceptible E. faecalis as expected.

To demonstrate this digital approach, we utilized a simple limiting dilution of heterogeneous samples and coupled with “quasi-digital” assay using microtiter plates to analyze individual bacterial DNA molecule. Here, we mixed gentamicin-susceptible S. aureus with gentamicin-susceptible E. faecalis as the sample. By using digital melt curve classification for bacterial species identification (see below), we identified on average 33  7.9 positive reactions as S. aureus and 10  1.5 positive reactions as E. faecalis from triplicated 96-well plates of no-gentamicin control (a representative HRM result and plate layout are shown in Fig. 5A and 5B; S. aureus – orange and E. faecalis – blue). For the gentamicin-treated sample, we identified on average 7  1.5 positive reactions as S. aureus and 3  0.6 positive reactions as E. faecalis from triplicated 96-well plates (a representative HRM result and plate layout are shown in Fig. 5C and 5D; S. aureus – orange and E. faecalis – blue). Based on these average absolute counts of both bacterial species (Fig. 5E), there were approximately 5-fold and 3-fold decreases for S. aureus and E. faecalis total counts in gentamicin-treated sample, indicating that gentamicin had effectively inhibited the growth of both bacterial species. These results demonstrate that digital assay allows for enumerating multiple species presented in a sample, which also enables AST measurement for each individual organism by comparing the total number of molecules of each bacterial species between drug and no drug control. There are some important considerations that are essential in analyzing dPCR-dHRM. First, due to the limited dynamic range of a 96- well plate, we optimized the colony-forming unit (CFU) number of each species so that both of them could be detected at single molecule level and provided positive counts that were indicatives of their AST profiles. Second, because of only a single molecule presented in each PCR reaction, the resulting melt curves exhibit higher variations in shapes. A new training database was generated based on single molecule-derived melt curves. Digital melt curves in this database were categorized into three groups; S. aureus, E. faecalis and unspecific group, which was collected from water control runs and Sanger sequencing-verified contaminants (Fig. S-6A). The ovoSVM algorithm has also been amended for digital melt curve classification (DovoSVM) by incorporating screening criteria and identifying if the curve matches to one of the three categories or no match as described in Figure S-6B

ACS Paragon Plus Environment

Page 11 of 15

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Analytical Chemistry

(result shown in Fig. S-6C). These additional modules help us to achieve accurate digital melt curve-based classification with a small training data set.

Conclusions In summary, we have introduced an integrated, PCR-based workflow that achieves bacterial identification and phenotypic AST in a single test. Rapid broad-based 16S PCR allows for measuring bacterial growth after brief exposure to antibiotics and that correlates with susceptibility. By integrating with HRM, any presented bacteria can be detected and species-identified with machine learning algorithm. Our integrated workflow was validated with a blinded test in which we correctly identified the bacteria species and gentamicin susceptibility profiles of 3 unknown bacterial samples. Importantly, our assay could be completed in less than 6.5 hours, representing significant time saving. Finally, as an added benefit, we demonstrate that ID and AST of multiple bacteria can be achieved by implementing dPCR and dHRM in our assay. While we have demonstrated and highlighted the feasibility of the proposed approach, the method should be further enhanced in several areas. First, sample preparation process can be incorporated to retrieve viable bacterial cells from various types of sample. Second, both the HRM database and the AST references should be expanded to cover more bacterial species and different antibiotics, which would make our assay applicable to various infections. Factors that are possibly affect bacterial growth including bacterial species, bacterial load, bacterial viability, culture broth, antibiotic type, and antibiotic dosage can be further studied and profiled for better understanding and improving assay design. Lastly, the dynamic range of digital experiment can be increased by performing in a higher throughput format such as microfluidic digital PCR57. This dPCR also provides superior detection sensitivity and precise quantification capability, which results in shorter incubation time comparing to traditional qPCR. We envision that this work has potential to be clinically applied in guiding antibiotic treatments, which ultimately improves global healthcare.

Supporting Information The Supporting Information is available free of charge on the ACS Publications website. Additional experimental results and data analysis including melt curves comparison between different platforms, similar melt curves from different strains of the same bacterial species, bacterial growth measurement via qPCR, MIC determination by using our proposed method and more details on digital bacterial ID/AST workflow (PDF)

Acknowledgments This work has been supported by National Institutes of Health (R01AI117032), National Science Foundation (1159771 and 1033744) and the Royal Thai Government Scholar program. We thank Siwat Jakaratanopas for helping in the table of content graphic preparation.

ACS Paragon Plus Environment

Analytical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

References (1)

Angus, D. C.; Linde-Zwirble, W. T.; Lidicker, J.; Clermont, G.; Carcillo, J.; Pinsky, M. R. Crit. Care Med. 2001, 29, 1303– 1310.

(2)

Cohen, J.; Vincent, J. L.; Adhikari, N. K. J.; Machado, F. R.; Angus, D. C.; Calandra, T.; Jaton, K.; Giulieri, S.; Delaloye, J.; Opal, S.; Tracey, K.; van der Poll, T.; Pelfrene, E. Lancet Infect. Dis. 2015, 581–614.

(3)

Davenport, M.; Mach, K. E.; Shortliffe, L. M. D.; Banaei, N.; Wang, T.-H.; Liao, J. C. Nat. Rev. Urol. 2017, 14, 296–310.

(4)

Dellinger, R. P.; Levy, M. M.; Rhodes, A.; Annane, D.; Gerlach, H.; Opal, S. M.; Sevransky, J. E.; Sprung, C. L.; Douglas, I. S.; Jaeschke, R.; Osborn, T. M.; Nunnally, M. E.; Townsend, S. R.; Reinhart, K.; Kleinpell, R. M.; Angus, D. C.; Deutschman, C. S.; Machado, F. R.; Rubenfeld, G. D.; Webb, S. A.; Beale, R. J.; Vincent, J.-L.; Moreno, R. Crit. Care Med. 2013, 41, 580–637.

(5)

Frieden, T. Centers Dis. Control Prev. 2013, 114.

(6)

Costelloe, C.; Metcalfe, C.; Lovering, A.; Mant, D.; Hay, A. D. Bmj 2010, 340, c2096.

(7)

House, W. The White House, Washington, DC 2015.

(8)

Carbonnelle, E.; Mesquita, C.; Bille, E.; Day, N.; Dauphin, B.; Beretti, J. L.; Ferroni, A.; Gutmann, L.; Nassif, X. Clin. Biochem. 2011, 104–109.

(9)

Stevenson, L. G.; Drake, S. K.; Murray, P. R. J. Clin. Microbiol. 2010, 48, 444–447.

(10)

Jarvis, R. M.; Goodacre, R. Anal. Chem. 2004, 76, 40–47.

(11)

Liu, C.-Y.; Han, Y.-Y.; Shih, P.-H.; Lian, W.-N.; Wang, H.-H.; Lin, C.-H.; Hsueh, P.-R.; Wang, J.-K.; Wang, Y.-L. Sci. Rep. 2016, 6, 23375.

(12)

Peruski, A. H.; Peruski, L. F. J. Clin. Diagn. Lab. Immunol. 2003, 10, 506–513.

(13)

Thomas, L. C.; Gidding, H. F.; Ginn, A. N.; Olma, T.; Iredell, J. J. Microbiol. Methods 2007, 68, 296–302.

(14)

Wellinghausen, N.; Nöckler, K.; Sigge, A.; Bartel, M.; Essig, A.; Poppert, S. J. Clin. Microbiol. 2006, 44, 1828–1830.

(15)

Westh, H.; Lisby, G.; Breysse, F.; Böddinghaus, B.; Chomarat, M.; Gant, V.; Goglio, A.; Raglio, A.; Schuster, H.; Stuber, F.; Wissing, H.; Hoeft, A. Clin. Microbiol. Infect. 2009, 15, 544–551.

(16)

Wojewoda, C. M.; Sercia, L.; Navas, M.; Tuohy, M.; Wilson, D.; Hall, G. S.; Procop, G. W.; Richter, S. S. J. Clin. Microbiol. 2013, 51, 2072–2076.

(17)

Peter, H.; Berggrav, K.; Thomas, P.; Pfeifer, Y.; Witte, W.; Templeton, K.; Bachmann, T. T. J. Clin. Microbiol. 2012, 50, 3990–3998.

(18)

Peters, R. P. H.; Van Agtmael, M. A.; Gierveld, S.; Danner, S. A.; Groeneveld, A. B. J.; Vandenbroucke-Grauls, C. M. J. E.; Savelkoul, P. H. M. J. Clin. Microbiol. 2007, 45, 3641–3646.

(19)

Sango, A.; McCarter, Y. S.; Johnson, D.; Ferreira, J.; Guzman, N.; Jankowski, C. A. J. Clin. Microbiol. 2013, 51, 4008–4011.

(20)

Shang, S.; Chen, G.; Wu, Y.; Du, L.; Zhao, Z. Pediatr. Res. 2005, 58, 143–148.

(21)

Jordan, J. A.; Jones-Laughner, J.; Durso, M. B. J. Clin. Microbiol. 2009, 47, 368–372.

(22)

Beal, S. G.; Ciurca, J.; Smith, G.; John, J.; Lee, F.; Doern, C. D.; Gander, R. M. J. Clin. Microbiol. 2013, 51, 3988–3992.

(23)

Blaschke, A. J.; Heyrend, C.; Byington, C. L.; Fisher, M. A.; Barker, E.; Garrone, N. F.; Thatcher, S. A.; Pavia, A. T.; Barney, T.; Alger, G. D.; Daly, J. A.; Ririe, K. M.; Ota, I.; Poritz, M. A. Diagn. Microbiol. Infect. Dis. 2012, 74, 349–355.

(24)

Boedicker, J. Q.; Li, L.; Kline, T. R.; Ismagilov, R. F. Lab Chip 2008, 8, 1265–1272.

(25)

Brolund, A.; Tegmark Wisell, K.; Edquist, P. J.; Elfström, L.; Walder, M.; Giske, C. G. J. Microbiol. Methods 2010, 82, 229– 233.

ACS Paragon Plus Environment

Page 12 of 15

Page 13 of 15

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Analytical Chemistry

(26)

Buchan, B. W.; Ginocchio, C. C.; Manii, R.; Cavagnolo, R.; Pancholi, P.; Swyers, L.; Thomson, R. B.; Anderson, C.; Kaul, K.; Ledeboer, N. A. PLoS Med. 2013, 10.

(27)

Forrest, G. N.; Mankes, K.; Jabra-Rizk, M. A.; Weekes, E.; Johnson, J. K.; Lincalis, D. P.; Venezia, R. A. J. Clin. Microbiol. 2006, 44, 3381–3383.

(28)

Grif, K.; Fille, M.; Würzner, R.; Weiss, G.; Lorenz, I.; Gruber, G.; Eschertzhuber, S.; Nachbaur, D.; Lass-Flörl, C.; Orth, D. Wien. Klin. Wochenschr. 2012, 124, 266–270.

(29)

Raja, S.; Ching, J.; Xi, L.; Hughes, S. J.; Chang, R.; Wong, W.; McMillan, W.; Gooding, W. E.; McCarty, K. S.; Chestney, M.; Luketich, J. D.; Godfrey, T. E. Clin. Chem. 2005, 51, 882–890.

(30)

Ulrich, M. P.; Christensen, D. R.; Coyne, S. R.; Craw, P. D.; Henchal, E. A.; Sakai, S. H.; Swenson, D.; Tholath, J.; Tsai, J.; Weir, A. F.; Norwood, D. A. J. Appl. Microbiol. 2006, 100, 1011–1016.

(31)

Poritz, M. A.; Blaschke, A. J.; Byington, C. L.; Meyers, L.; Nilsson, K.; Jones, D. E.; Thatcher, S. A.; Robbins, T.; Lingenfelter, B.; Amiott, E.; Herbener, A.; Daly, J.; Dobrowolski, S. F.; Teng, D. H. F.; Ririe, K. M. PLoS One 2011, 6.

(32)

Edwards, M. C.; Gibbs, R. A. Genome Res. 1994, 3.

(33)

Barghouthi, S. A. Indian J. Microbiol. 2011, 51, 430–444.

(34)

Clifford, R. J.; Milillo, M.; Prestwood, J.; Quintero, R.; Zurawski, D. V; Kwak, Y. I.; Waterman, P. E.; Lesho, E. P.; Mc Gann, P. PLoS One 2012, 7, e48558.

(35)

Jeng, K.; Yang, S.; Won, H.; Gaydos, C. A.; Hsieh, Y. H.; Kecojevic, A.; Carroll, K. C.; Hardick, J.; Rothmana, R. E. J. Clin. Microbiol. 2012, 50, 1122–1124.

(36)

Masek, B. J.; Hardick, J.; Won, H.; Yang, S.; Hsieh, Y. H.; Rothman, R. E.; Gaydos, C. A. J. Mol. Diagnostics 2014, 16, 261– 266.

(37)

Perng, C. L.; Chen, H. Y.; Chiueh, T. S.; Wang, W. Y.; Huang, C. T.; Sun, J. R. J. Med. Microbiol. 2012, 61, 944–951.

(38)

Rothman, R.; Ramachandran, P.; Yang, S.; Hardick, A.; Won, H.; Kecojevic, A.; Quianzon, C.; Hsieh, Y. H.; Gaydos, C. Acad. Emerg. Med. 2010, 17, 741–747.

(39)

Sauer, P.; Gallo, J.; Kesselová, M.; Kolár, M.; Koukalová, D. Biomed. Pap. Med. Fac. Univ. Palacky. Olomouc. Czech. Repub. 2005, 149, 285–288.

(40)

Srinivasan, R.; Karaoz, U.; Volegova, M.; MacKichan, J.; Kato-Maeda, M.; Miller, S.; Nadarajan, R.; Brodie, E. L.; Lynch, S. V. PLoS One 2015, 10.

(41)

Won, H.; Rothman, R.; Ramachandran, P.; Hsieh, Y. H.; Kecojevic, A.; Carroll, K. C.; Aird, D.; Gaydos, C.; Yang, S. J. Clin. Microbiol. 2010, 48, 3410–3413.

(42)

Yang, S.; Ramachandran, P.; Rothman, R.; Hsieh, Y. H.; Hardick, A.; Won, H.; Kecojevic, A.; Jackman, J.; Gaydos, C. J. Clin. Microbiol. 2009, 47, 2252–2255.

(43)

Yarza, P.; Yilmaz, P.; Pruesse, E.; Glöckner, F. O.; Ludwig, W.; Schleifer, K.-H.; Whitman, W. B.; Euzéby, J.; Amann, R.; Rosselló-Móra, R. Nat. Rev. Microbiol. 2014, 12, 635–645.

(44)

Fettweis, J. M.; Serrano, M. G.; Sheth, N. U.; Mayer, C. M.; Glascock, A. L.; Brooks, J. P.; Jefferson, K. K.; Buck, G. A. BMC Genomics 2012, 13, S17.

(45)

Clarridge, J. E. Clin. Microbiol. Rev. 2004, 17, 840–862.

(46)

Fraley, S. I.; Hardick, J.; Masek, B. J.; Athamanolap, P.; Rothman, R. E.; Gaydos, C. A.; Carroll, K. C.; Wakefield, T.; Wang, T. H.; Yang, S. Nucleic Acids Res. 2013, 41, e175.

(47)

Andini, N.; Wang, B.; Athamanolap, P.; Hardick, J.; Masek, B. J.; Thair, S.; Hu, A.; Avornu, G.; Peterson, S.; Cogill, S.; others. Sci. Rep. 2017, 7.

(48)

Athamanolap, P.; Parekh, V.; Fraley, S. I.; Agarwal, V.; Shin, D. J.; Jacobs, M. A.; Wang, T. H.; Yang, S. PLoS One 2014, 9, e109094.

ACS Paragon Plus Environment

Analytical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

(49)

Fraley, S. I.; Athamanolap, P.; Masek, B. J.; Hardick, J.; Carroll, K. C.; Hsieh, Y.-H.; Rothman, R. E.; Gaydos, C. A.; Wang, T.-H.; Yang, S. Sci. Rep. 2016, 6, 19218.

(50)

Walsh, C. Nature 2000, 406, 775–781.

(51)

Cirz, R. T.; Chin, J. K.; Andes, D. R.; De Crécy-Lagard, V.; Craig, W. A.; Romesberg, F. E. PLoS Biol. 2005, 3, 1024–1033.

(52)

Davies, J.; Davies, D. Microbiol. Mol. Biol. Rev. 2010, 74, 417–433.

(53)

Waldeisen, J. R.; Wang, T.; Mitra, D.; Lee, L. P. PLoS One 2011, 6.

(54)

Mezger, A.; Gullberg, E.; Göransson, J.; Zorzet, A.; Herthnek, D.; Tano, E.; Nilsson, M.; Andersson, D. I. J. Clin. Microbiol. 2015, 53, 425–432.

(55)

Bustin, S. A.; Benes, V.; Garson, J. A.; Hellemans, J.; Huggett, J.; Kubista, M.; Mueller, R.; Nolan, T.; Pfaffl, M. W.; Shipley, G. L.; Vandesompele, J.; Wittwer, C. T. Clin. Chem. 2009, 55, 611–622.

(56)

Park, S. Y.; Park, K. H.; Bang, K. M.; Chong, Y. P.; Kim, S. H.; Lee, S. O.; Choi, S. H.; Jeong, J. Y.; Woo, J. H.; Kim, Y. S. J. Infect. 2012, 65, 119–127.

(57)

Ottesen, E. A.; Hong, J. W.; Quake, S. R.; Leadbetter, J. R. Science 2006, 314, 1464–1467.

ACS Paragon Plus Environment

Page 14 of 15

Page 15 of 15

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Analytical Chemistry

For Table of Contents Only

ACS Paragon Plus Environment