Samples
(i) Clinical samples Analytical clinical samples included conventional blood culture samples and blood samples that had been treated using a pre-analytical treatment protocol (detailed under “Pre-treatment protocol” in the Materials and methods section), and both of these were collected according to the Chinese guidelines of Specimen collection and transport in clinical microbiology (WS/T 640-2018). Versa TREK Redox bottles and the Versa TREK automated blood culture system (Thermo Fisher Scientific, MA, USA) were used for blood culture. Blood culture samples were collected consecutively during January 2019 to December 2020 at Huadong hospital based on the routine identification results. Inclusion criteria were as follows: 1) all patients had a positive blood culture for at least one of the target pathogens, 2) clinical manifestations of bloodstream infection, and 3) inpatients with complete clinical microbiological data for analysis. Repeat blood cultures within 14 days of the initial blood culture and yielding an identical pathogen were excluded. In total, 426 clinical positive blood culture samples were collected in the current study. Moreover, 100 clinical negative blood culture samples were also used to test the performance of the new methods.Meanwhile, duplicate sets of blood samples were collected from 214 enrolled patients who volunteered to provide double samples, and one set of these samples was subjected to conventional testing and the other one was subjected to the pre-analytical treatment protocol and identification using the BSI-HMGS.
(ii) Simulated samples Positive whole blood samples were simulated to evaluate the efficiency of the pre-analytical treatment protocol and BSI-HMGS. The simulation workflow was as follows: 1) The bacterial/fungal reference strains were incubated on Columbia blood agar plates (Haemophilus influenzae on chocolate agar plate) and cultured at 35°C in an atmosphere of 5% CO2 for 24 h. Then the colonies were picked from the solid mediums and diluted with 0.45% NaCl to reach turbidity of 0.5 McFarland turbidity unit (MCF) using a DensiCHEK turbidimeter (BioMerieux, Lyon, France), and these suspensions were serially diluted to the desired concentrations and incubated on appropriate solid culture mediums until colony-forming units (CFU) were visible (the bacterial suspension was totally diluted about 1.0 × 105 fold and the fungal suspension was totally diluted about 1.0 × 103 fold). The final countable colonies were shown in Supplemental Figure S1 of Additional file 1. Accurate CFU/mL of each 0.5 MCF suspension was determined by plate count. 2) Then diluted each 0.5 MCF suspension with 0.45% NaCl and spiked in sterile defibrinated sheep blood (Yuanye Bio-Technology, Shanghai, China) to yield load 10 CFU/mL of simulated whole blood sample [21]. 3) Finally, injected 10 mL the spiked blood samples to Versa TREK Redox bottles and loaded the bottles in the Versa TREK blood culture system for further pre-analytical treatment. The concentration–response curves, including the information of time to positivises, from the Versa TREK system for these bottles were shown as Supplemental Figure S2 in Additional file 1, which was consistent with that reported in previous studies [22]. Five clinical isolates from different patients were used as simulation bacteria/fungi for each pathogen. In total, 720 simulated samples (30 for each organism) were obtained in this study.
(iii) Reference strains Twenty-four reference strains relative to the pathogens in the panel were used to establish and optimize the BSI-HMGS. Including: Klebsiella pneumoniae (ATCC700603), Burkholderia cepacia (ATCC25416), Proteus mirabilis (ATCC12453), Moraxella catarrhalis (ATCC25238), Serratia marcescens (ATCC8100), H. influenzae (ATCC49766), Pseudomonas aeruginosa (ATCC27853), Enterobacter cloacae (ATCC13047), Escherichia coli (ATCC25922), Acinetobacter baumannii (ATCC19606), Stenotrophomonas maltophilia (ATCC17666), Salmonella typhimurium (ATCC14028), Staphylococcus hominis (ATCC27844), Staphylococcus aureus (ATCC29213), Streptococcus salivarius(ATCC13419), Streptococcus pneumoniae (ATCC49619), Streptococcus pyogenes (ATCC12344), Streptococcus agalactiae (ATCC12388), Enterococcus faecalis (ATCC29212), Enterococcus faecium (ATCC35667), Candida albicans (ATCC14053), Candida tropicalis (ATCC66029), Candida parapsilosis (ATCC22019) and Candida glabrata (ATCC15126).
Conventional tests
Positive clinical blood culture samples in the Versa TREK system were transferred to Columbia blood agar, Chocolate agar, and Sabouraud’s agar. Negative samples were incubated for at least 5 days in the automated blood culture system, and all incubation and sub-cultivation procedures were performed following the Clinical and Laboratory Standards Institute (CLSI, M47-A) guidelines. The pathogens cultured in the medium were identified by Matrix-Assisted Laser Desorption Ionization time-of-flight mass spectrometry (MALDI-TOF MS, BioMerieux), and these results were used for comparison. The samples were also detected using BSI-HMGS. If the BSI-HMGS result was not consistent with the corresponding result in MALDI-TOF MS, sequencing (Shenggong Bioengineering Technology, Shanghai, China) was performed for verification and this result was determined as the final reference result. If BSI-HMGS yielded the same detection result as MALDI-TOF MS, this result was directly determined as the reference result.
Pre-treatment protocol
To accelerate BSI pathogen identification using PCR-based methods, a protocol for the enrichment of pathogens from whole blood sample and depleting PCR inhibitors prior to extraction of microbial DNA was developed. The details are as follows:
(i) Incubation To increase pathogen load in whole blood samples from patients with suspected BSI, one set of the blood samples from the enrolled patients was collected into Versa TREK Redox bottles and incubated in the Versa TREK automated blood culture system for 5 h, as were the simulated whole blood samples.
(ii) Concentration The amount of blood sample utilized for nucleic acid extraction can affect the detection sensitivity of PCR-based methods. Therefore, we concentrated the pathogen from 5 mL blood culture samples for the subsequent BSI-HMGS assay. First, 5 mL of suspension from step (i) was added to a Serum separation tube (SST) (Becton, Dickinson and Company, NJ, USA). Second, the SST was centrifuged at 900 g for 10 min and the supernatant was discarded; the pathogen cells would be concentrated on the surface of the separation gel. Third, 1 mL of 0.45% NaCl was added to the SST, and shaken intermittently using a vortex mixer for 15 s to resuspend the pathogen cells. Finally, the suspension above the separation gel was transferred to a 2 mL centrifugation tube.
(iii) Inhibitor removal The existence of PCR inhibitors in blood culture bottles, including 3 dominant vendors Versa TREK Redox, BACTEC Vials, and BacT/AlerT, has been proven in the following studies (detailed under “Confirmation of the presence of PCR inhibitors in the blood culture system” in the Materials and methods section), even after nucleic acid extraction. Thus, the 1 mL suspension obtained in step (ii) was purified to remove any inhibitors. First, 1 mL of selective lysis buffer (Fosun Pharmaceutical, Shanghai, China), mainly containing 10% sodium dodecyl sulphate (SDS) and 1 μL glycogen (Yuanye Bio-Technology, 20 mg/mL) as a carrier for the precipitation of DNA, was added to the suspension, agitated for 15 s, and centrifuged at 15850 g for 1 min. Then, the supernatant was discarded and the sediment was washed using 1 mL wash buffer (Fosun Pharmaceutical), mainly containing Tris buffer (PH 8.0), in a manner similar to that used for the lysis buffer. The effect of the 10% SDS treatment on the target organisms was further evaluated to confirm its availability: each organism suspension was treated with the lysis buffer and 0.45% NaCl in the same way above, respectively, then the sediment was examined by smear microscopy (See Supplemental Figure S3, Additional file 1).
Nucleic acid extraction
Total nucleic acids were extracted from a 300 μL suspension of blood culture samples or from the pre-treated sample indicated above, using a Nucleic Acid Extraction or Purification Kit (HEALTH Gene Technologies, Ningbo, China) according to the manufacturer’s instructions. This kit was based on the principle of guanidine isolation and silica dioxide adsorption using a magnetic bead-based method. Meanwhile, a modified fragment of the kanamycin resistance gene (Kanr) was inserted into the pcDNA3.1 vector to generate a fusion plasmid that served as internal control (IC) for the detection system, and the IC plasmid (3 μL, 1.0 × 105 copies/μL) was added to every clinical sample prior to extraction as a positive control throughout the whole experiment. The extracts were eluted with 60 μL of DNase/Rnase-free H2O (ddH2O). The concentration of each extract was determined using a Thermo Nanodrop 2000 spectrophotometer (Thermo Fisher Scientific). The extracts were stored at -80°C until further analysis. Another nucleic acid extraction kit (Tiangen Biotech, Beijing, China) based on an absorbent column principle was also tested to verify its DNA purification efficiency.
BSI-HMGS primer design
Twenty-five pairs of primers targeting one IC gene and twenty-four conserved genes of specific pathogens, including 12 gram-negative bacterial species/complex (K. pneumoniae, B. cepacia, P. mirabilis, M. catarrhalis, S. marcescens, H. influenzae, P. aeruginosa, Enterobacter cloacae complex, E. coli, A. baumannii, S. maltophilia and S. enterica), 8 gram-positive bacterial genus/species (Staphylococcus, S. aureus, Streptococcus, S. pneumoniae, S. pyogenes, S. agalactiae, E. faecalis, and E. faecium), and 4 fungi (C. albicans, C. tropicalis, C. parapsilosis, and C. glabrata), were designed using Premier 5.0 software (Premier Biosoft International, CA, USA) (Table 1). To achieve sufficient coverage of pathogens and to ensure that the primer pairs work appropriately, two pairs of primers were targeted at genus level and one was targeted at the bacterial complex, including primer pairs for Staphylococcus and Streptococcus generaand the Enterobacter cloacae complex, respectively. The primer pair for Staphylococcus targeted a specific gene present in S. aureus, S. lugdunensis and eight kinds of predominant coagulase-negative staphylococci including S. capitis, S. epidermidis, S. hominis, S. haemolyticus, S. intermedius, S. saprophyticus, S. simulans, and S. warneri. Thus, identification of S. aureus by BSI-HMGS would show two amplified peaks, one at the position for S. aureus and the other at the position for Staphylococcus. The primer pair for Streptococcus targeted a specific gene present in eight Streptococcus species including S. agalactiae, S. pneumoniae, S. pyogenes, S. milleri, S. mitis, S. mutans, S. salivarius, and S. sanguinis. Thus, identification of S. agalactiae, S. pneumoniae, and S. pyogenes by BSI-HMGS would also show two amplified peaks each, one at the position for Streptococcusand the other at their specific position. In addition, the E. cloacae complex in the BSI-HMGS panel includes seven species: E. arbornii, E. cloacae, E. hokalii, E. kodoi, E. lukalii, E. mori, and E. nimipressuralis, which would show a single peak at the position for E. cloacae complex when detected by BSI-HMGS. Overall, 42 bacteria and fungi could be detected with these 25 pairs of primers. All amplification product sizes ranged from 100–400 bp, with at least 3 base differences between each product. The peaks of primers dimers would show between 0–100 nt. The detailed primer information is listed in Table 1.
Multiplex PCR and fragment analysis for the BSI-HMGS assay
Multiplex PCR was performed in a 10 μL volume containing 2 μL PCR Buffer (HEALTH Gene Technologies), 0.5 μL uracil-N-glycosylase (UNG), 0.4 μL Phoenix™ Hot Start Taq DNA Polymerase (Qiagen, Dusseldorf, Germany), 1 μL primer pool, 1.1 μL ddH2O, and 5 μL template. The primer pool comprised 25 mixed pairs of primers at different proportions to achieve optimum sensitivity for all targets (Table 1). The mixture was subjected to the following amplification conditions: 42℃ for 5 min; 94℃ for 8 min; 35 cycles of 94℃ for 30 s, 60℃ for 30 s, and 70℃ for 1 min; and 70℃ for 1 min. The products were stored at 4°C until further use. The ddH2O was simultaneously detected as a negative control throughout the whole study. The multiplex PCR takes approximately 1 h and 24 min for 1-96 samples in one run.
The PCR products were then prepared for fragment analysis using the Genetic Analyzer 3500xL (Applied Biosystems, CA, USA). For each reaction, 1 μL PCR product was added to 9 μL highly deionized formamide containing 4% GeneScan 500 LIZ (Applied Biosystems) as internal standard size marker, and analysed using high-resolution capillary gel electrophoresis. The data obtained was further analysed using the Gene Mapper ID-Xsoftware v1.6 (Applied Biosystems). Finally, detection was considered positive when the peak height was greater than 500 relative fluorescence units (rfu). The fragment analysis takes approximately 40 min for 1-24 samples in one run.
Plasmids construction
Plasmids of all the target genes were used to establish and optimize the BSI-HMGS. The process for construction and transformation of the 25 plasmids was as follows: 1) Ligation: 2 μL PCR products, 0.5 μL PMD-18T vector (Tiangen Biotech) and 2.5 μL solution Ⅰ were mixed and kept under 16℃ for at least 45 min. 2) Transformation: 5 μL above product mixed with 1mL competent cells (Tiangen Biotech) were processed as follows: incubation on ice-bath for 30 min; 42℃ thermal shock for 90 s; ice-bath for 2 min; 600 μL LB broth without ampicillin resistance was added and then placed in shaker under 37℃ for 45 min. 3) Inoculation: 200 μL transformed cells were inoculated under 37℃ for 12 h. 4) Identification: the target clone was identified by sequencing. DNA copy number was calculated by the following formula: [(6.02 × 1023 copy number/mol) × plasmid concentration (g/mL)/L × (MW g/mol) = copies/mL (MW: average molecular weight). Purification of plasmids was performed according to the manufacturer’s protocol.
Confirmation of PCR inhibitors in the blood culture system
To determine whether the blood culture matrix components in the blood culture bottles contained PCR inhibitors, PCR-based detection was conducted to test the mixtures of the matrix components and specific bacterial suspensions. Six kinds of matrix components from three predominant vendors (aerobic/anaerobic) including Versa TREK Redox (Thermo Fisher Scientific), BACTEC Vials (Becton, Dickinson and Company, NJ, USA), and BacT/AlerT (BioMerieux), were tested. The matrix component in each bottle (100 μL) was mixed with specific bacterial suspension (200 μL, 1.0 × 104 CFU/mL). Then the mixtures were extracted, amplified, and identified using common multiplex PCR Buffer and Taq DNA Polymerase (Roche, Basel, Switzerland), respectively. The mixture of bacterial suspension and matrix components from Versa TREK Redox, which was used for the blood culture of clinical samples in the current study, was subjected to the pre-analytical inhibitor removal protocol and then identified by the BSI-HMGS to verify its ability to counter possible inhibition. The mixture of bacterial suspension and ddH2O was also detected as a positive control. Two bacteria, P. mirabilis and S. aureus, were used for this validation. Three replicates were conducted for each test.
Performance evaluation of BSI-HMGS and the pre-treatment protocol
The 24 reference strains in the panel as well as five pathogens detected in BSIs but not included in the BSI-HMGS panel, were tested to verify that the BSI-HMGS is specific to target pathogens. Quantified bacterial/fungal suspensions were tested to evaluated the limit of detection (LOD) of BSI-HMGS. Then, serial tenfold dilutions of plasmids mixed with equal amounts of templates were used to test the ability of BSI-HMGS to detect all pathogens simultaneously. Different two-plasmid mixtures (including gram-negative bacterial mixed with gram-positive bacterial plasmids, gram-negative bacterial with gram-negative bacterial plasmids, gram-positive bacterial with gram-positive bacterial plasmids, gram-negative bacterial with fungal plasmids, and gram-positive bacterial with fungal plasmids) were examined using BSI-HMGS to verify its capability of detecting polymicrobial infection. The concentrations of the two plasmids in each mixture had a tenfold difference. The performance of BSI-HMGS on blood culture samples was evaluated by comparison with the reference results. Clinical whole blood samples from the enrolled patients and simulated whole blood samples were used to test the performance of BSI-HMGS combined with the pre-treatment protocol.
Statistics
The level of agreement between the assays was determined by calculating the overall percent agreement (OPA), positive percent agreement (PPA) and negative percent agreement (NPA) using the following formulas: OPA = (TP + TN)/(TP + FP + TN + FN); PPA = TP/(TP + FN); NPA = TN/(TN + FP) (where TP: true positive; FN: false negative; TN: true negative; FP: false positive, and all these true or false were determined based on the reference results). All statistical analyses were performed using Stata/SE 14.0 (Stata Corp College Station, TX, USA) for Mac.