Comparison of microbial molecular diagnosis efficiency within unstable template metagenomic DNA samples between qRT-PCR and chip-based digital PCR platforms
Article information
Abstract
Accurate and efficient microbial diagnosis is crucial for effective molecular diagnostics, especially in the field of human healthcare. The gold standard equipment widely employed for detecting specific microorganisms in molecular diagnosis is quantitative real-time polymerase chain reaction (qRT-PCR). However, its limitations in low metagenomic DNA yield samples necessitate exploring alternative approaches. Digital PCR, by quantifying the number of copies of the target sequence, provides absolute quantification results for the bacterial strain. In this study, we compared the diagnostic efficiency of qRT-PCR and digital PCR in detecting a particular bacterial strain (Staphylococcus aureus), focusing on skin-derived DNA samples. Experimentally, specific primer for S. aureus were designed at transcription elongation factor (greA) gene and the target amplicon were cloned and sequenced to validate efficiency of specificity to the greA gene of S. aureus. To quantify the absolute amount of microorganisms present on the skin, the variable region 5 (V5) of the 16S rRNA gene was used, and primers for S. aureus identification were used to relative their amount in the subject’s skin. The findings demonstrate the absolute convenience and efficiency of digital PCR in microbial diagnostics. We suggest that the high sensitivity and precise quantification provided by digital PCR could be a promising tool for detecting specific microorganisms, especially in skin-derived DNA samples with low metagenomic DNA yields, and that further research and implementation is needed to improve medical practice and diagnosis.
Introduction
The development of the digital polymerase chain reaction (PCR) method, third-generation PCR equipment, has recently changed the research trend of the existing microbial molecular diagnosis field. Especially, a new digital PCR equipment ‘LOAA (Lab On An Array) digital real-time PCR analyzer (Optolane, Seongnam, Korea)’ has 1,000 times more fluorescence detection sensitivity than quantitative real-time polymerase chain reaction (qRT-PCR), second-generation PCR equipment, and is a high-performance PCR equipment that can detect each PCR amplification reaction in more than 20,000 nano-sized PCR reaction wells present on a single semiconductor chip-based micro electro mechanical system (MEMS) [1,2]. Above all things, unlike qRT-PCR, LOAA digital PCR is very useful in the field of specific pathogens or cancer diagnosis because it has independent ‘absolute gene quantitative analysis systems’ within equipment without a separate standard curve quantification analysis [3,4]. In the case of LOAA digital PCR, about 20,000 PCR reaction wells are distributed on a semiconductor chip, and PCR amplification and fluorescence detection proceeds from 0 or 1 template DNA molecule included per well. This characteristic valid a ‘Poisson distribution’ principle that can calculate positive and negative reactions on PCR amplification detected in each well in detail, enabling absolute quantities in the sample for a specific gene to be identified without a separate standard curve analysis. The formula for absolute quantification method within a sample of a specific gene based on the 'Poisson distribution' principle is as follows:
However, although the digital PCR generates such high-efficiency molecular diagnostic outputs, the cost burden of one-time semiconductor-based experimental consumables is high, so if the quality and quality of template DNA are unstable during the experiment, consumer preference may decrease according to the experimental results. For example, metagenomic DNA (mDNA) samples isolated from skin applied in the human skin healthcare field are difficult to derive clear research results because their quantity and quality are often unstable due to various factors (e.g., Whether to wash, use skin care products, UV exposure, and individual living environment) affecting the growth environment of microorganisms on the human skin surface [5,6]. For this reason, the results of publishing research cases related to molecular microbial diagnosis using digital PCR in these unstable template DNA conditions, isolated from human skin, are insufficient, and for digital PCR to be optimized in various research fields in the future, it is necessary to identify these issues.
In this study, we aimed the comparing and verifying the quantitative efficiency of Staphylococcus aureus related to atopic disease occurrence within mDNA samples collected from lesions and non-lesions of seven atopic patients using qRT-PCR and digital PCR. We compared the specific gene detection and quantitative efficiency of the two PCR equipment under unstable template DNA conditions to evaluate the applicability of digital PCR to research fields such as skin microbiome studies that are highly affected by template DNA quality. In addition, by comparing the relative frequency difference of S. aureus between lesions and non-lesions skin sites in atopic patients, we demonstrated the reliability of the comparison results of quantitative efficiency for particular microorganisms between the two equipment [7]. Ultimately, we suggest through this study that digital PCR has high utilization value for various human healthcare industries related to molecular microbial diagnosis.
Methods and Results
Particular bacterial species-specific primer and probe sets design method
In this study, we tried to verify the detection and quantification efficiency of digital PCR equipment for specific bacterial species within the unstable template mDNA sample. Therefore, for this verification, we selected S. aureus, which is expected to exist on skin sites samples of atopic patients, as a specific bacterial species considering previous studies showing that S. aureus affects the atopic disease [8-10]. First, we designed a S. aureus–specific primer and probe set to detect and quantify a particular bacterial species present within mDNA samples extracted from various microbial pools collected from the surface of human skin, and to compare the detection and quantification efficiency between qRT-PCR and digital PCR. The experimental verification process for the primer and probe set design we conducted is as follows.
Bacterial gene selection for targeting particular species
The greA gene encodes a transcription elongation factor that affects bacterial gene transcription by regulating gene promoters, thereby regulating the environmental adaptation of bacteria [11]. In addition, the greA gene is recognized as a bacterial housekeeping gene, which, like the bacterial 16S ribosomal RNA, is an evolutionarily conserved transcription factor widely distributed in prokaryotes [12].
Retrieval of coding sequence region base information from reference database
We obtained sequence information (FASTA format) of the greA gene coding sequence region (CDS) for 22 different S. aureus strains (at the strain level) annotated in the National Center for Biotechnology Information (NCBI) reference database (Supplementary Table 1). This process is essential to improve primer and probe binding accuracy and specificity, as the exact strain information of S. aureus present within the human skin-derived mDNA samples applied in this study is unclear.
Selection of target-specific sequence regions for PCR reaction
To select primer sequences that could detect all 22 different S. aureus strains, multiple sequence alignment method (MSA) of each CDS information was performed using BioEdit 7.2.5v software (Supplementary Fig. 1A). Two consistent regions identified through the MSA method were selected as forward and reverse primer sequences (Table 1).
In silico test for pre-validating primer binding specificity
We used the Oligo calc (http://biotools.nubic.northwetern.edu/) and ‘Oligo Analysis (http://www.operon.com/tools/oligo-analysis-tool.aspx)’ open web tools to pre-simulate the suitability of the selected primer pairs for experimental application (including Tm values, GC%, and probability of primer dimer formation). Next, the NCBI nucleotide Basic Local Alignment Search Tool (BLAST) tool was used to confirm species-specificity for each primer sequence and amplification region included in each primer pair (Supplementary Fig. 1B).
Specific probe design
Finally, we designed a specific probe sequence region within between each selected primer sequence (Table 1, Supplementary Fig. 1A). The fluorescent reporter dyes and quencher applied in the probe design were 6-FAM (6-carboxy fluorescein) and SFCQ1 (SFC probe, Cheongju, Korea), which were attached to the 5' and 3' end regions of the selected probe sequences, respectively.
Comprehensive molecular genetical validation for primer/probe-specificity
To demonstrate the species-specificity of the pre-designed S. aureus–specific primers and probe through experimental validation, we set up a comparison group (Table 2). Information about the comparison group is as follows—positive control-1 (PC1): genomic DNA (gDNA) of a single strain of S. aureus ATCC6538; for gDNA extraction of PC, HiGene Genomic DNA Prep Kit For microorganisms (BIOFACT, Daejeon, Korea) was used, and all experimental procedures were performed according to the official protocol guide provided in the kit; negative control (NC) : ‘Siga-Microbial community DNA mix MBD0026’ contains genomic DNA from 10 bacterial species (Akkermansia muciniphila, Bacillus subtilis, Burkholderia pyrrocinia, Escherichia coli, Enterococcus faecalis, Pseudomonas aeruginosa, Proteus mirabilis, Proteus vulgaris, Porphyromonas gingivalis, and Salmonella enterica) are included in uniform proportions within one sample tube (Sigma-Aldrich, St. Louis, MO, USA); positive control-2 (PC2): genomic DNA of the S. aureus ATCC6538 was added to the NC; this control was set up to confirm that the specific-primers within the different microbial gDNA pools specifically bind to S. aureus. Therefore, we validated the experimental suitability (binding sensitivity) of S. aureus–specific primers and probes for application in qRT-PCR and digital PCR through a comprehensive molecular genetic experimental process. The experimental procedure was as follows.
General PCR validation for confirming primer binding specificity to S. aureus
First, we performed a general PCR validation to confirm that our pre-designed primers specifically bind to the CDS region of the greA gene on the S. aureus genome (Fig. 1A). PCR verification confirmed the presence of a DNA amplicon band approximately 146 bp long in PC1 and PC2, and no amplicon DNA was identified in the NC. This shows that the primers we designed bind specifically to S. aureus and do not bind to the NC (T100 Thermal Cycler, Bio-Rad, Hercules, CA, USA). The running condition of the PCR is as follows: pre-denaturation 95.0°C, 5 min, denaturation 95.0°C, 30 s, annealing 59°C, 40 s, elongation 72.0°C, 15 s, final-extension 72.0°C, 30 s, and total PCR cycle was 30. Additionally, we confirmed that specific primers are bound normally to S. aureus genomic DNA within a complex microbial DNA pool through PCR reaction results of PC2 and could confirm the potential applicable for specific detection of S. aureus even within skin-derived mDNA samples to be applied in this study.
Optimizing the primer and probe experimental conditions for qRT-PCR and digital PCR
Next, to set the optimal experimental conditions of primers and probes for application to qRT-PCR and digital PCR, preferentially, we performed the qRT-PCR validation using PC1 template DNA sample with S. aureus–specific primer, and the primer concentration was set the 25 pmol, probe concentration was set 20 pmol, and 10 pmol conditions, respectively. The qRT-PCR (CFX Opus 96 Real-Time PCE System, Bio-Rad; BioFACT 2× Multi-Star Real-Time PCR Master Mix For Probe, UDG system) experimental conditions were as follows (Supplementary Table 2); pre-denaturation 95.0°C 10 min, denaturation 95.0°C 10 s, annealing/Flour detection 61°C 10 s, elongation 72.0°C 15 s, and total PCR cycle was 45). As a result, the resolution of detection fluorescence value (RFU; relative fluorescence units) for the S. aureus greA gene detected on the probe concentration condition of 20 pmol about PC1 and PC2 samples were found to be higher compared to the 10 pmol concentration condition, confirming that this experimental condition was the most optimal (Fig. 1B).
Standard curve analysis about positive control (S. aureus) using qRT-PCR
We performed a standard curve analysis using the Ct value reflected about the positive control to contrast the results of qRT-PCR–based detection and quantification of S. aureus within mDNA samples applied in this study (Fig. 2, Supplementary Table 3). We performed qRT-PCR in triplicate to confirm the reliability of the standard curve analysis results and performed a standard curve analysis using 10-fold serial-diluted template DNA (PC1) samples (101, 100, 10-1, and 10-2 diluted samples, respectively). Standard curve analysis for PC1 based on qRT-PCR showed that the average Ct values for each dilution factor were 10 ng, 24.67; 1 ng, 28.90; 0.1 ng: 33.20; 0.01 ng, 38.24; and the R2 value of trend line was calculated on the standard curve graph was approximately 0.99.
Primer binding specificity cross-validation through Sanger sequencing validation
Finally, we performed a Sanger sequencing (ABI 3500 Genetic Analyzer, Thermo Fisher Scientific, Waltham, MA, USA) and NCBI nucleotide BLAST test to verify that the primer designed in this study correctly targeted the S. aureus greA gene CDS region and used it for PCR amplification (Fig. 1C). In this verification step, since the PCR amplification region to be identified is too short and unsuitable for conducting the Sanger sequencing, we performed the sequencing process by expanding the gene reading region through the TA cloning method (TOPcloner TA Kit, Enzynomics, Daejeon, Korea [13]). Competent cells required in the transformation process for gene cloning were performed using DH5α Chemically Competent E. coli (Enzynomics), and white colonies (potentially successful transformation and greA gene ligation) identified in 37°C incubation were selectively harvested and applied for sequencing process. The primer pair applied for Sanger sequencing to read the greA gene amplicon region was the M13 region-specific primer set contained in the TA cloning plasmid vector sequence. As a result, we confirmed from the base-pair sequence reading data that both forward and reverse primer sequence information were included. And then, to verify that the generated sequencing data contained CDS regions (PCR amplicon region) of the S. aureus greA gene, we input the sequencing data into the NCBI nucleotide BLAST search engine and compared it with annotated bacterial classification information within the NCBI reference database. As a result of the BLAST test, we confirmed that only bacterial identification information for S. aureus was matched from the NCBI reference database. We could also confirm that the BLAST test detected various S. aureus strains (at the strain level) rather than a single bacterial strain. This result showed that the S. aureus–specific primer, considered up to the strain level, has broad detection efficiency for various S. aureus strains within any type of clinical sample.
Detection and quantification for S. aureus using qRT-PCR method
We used human skin-derived mDNA samples as template DNA samples to confirm the efficiency of detection and quantification of specific microbial DNA within low molecule density template DNA using qRT-PCR. For this process, a total of 14 skin-derived clinical samples were obtained from seven patients with atopic dermatitis who visited the medical center (Kyung Hee University College of Medicine, Seoul, Republic of Korea), divided into lesion (case group) and non-lesion (control group) skin sites. Diagnosing atopic diseases of all participants was performed according to a dermatologist’s examination. Approval for the study protocol, informed consent forms, and related supporting documents was granted by the institutional review boards at the Korean Skin Research Center in South Korea (IRB No. HBABN01-220509-HRBR-E0113-01). All mDNA was extracted using the QIAamp PowerFecal Pro DNA kit (Qiagen, Hilden, Germany), and all experimental procedures were performed according to the official protocol guide provided with the kit. The average concentration values of extracted mDNA for each comparison group were checked with the control group 2.05 ng/μL and test group 4.10 ng/μL, respectively (Supplementary Table 4), by using a Thermo Scientific NanoDrop One/Onec Microvolume UV-Vis Spectrophotometer (Thermo Fisher Scientific). Next, we used qRT-PCR to determine the expected microbial frequency of S. aureus present within the 14 mDNA samples with a comparative analysis between groups to confirm the detection and quantification efficiency. The experimental procedure was as follows.
Normalization of bacterial quantification within each skin microbiome sample
Before the detection and quantitative analysis of S. aureus using qRT-PCR, we performed a qRT-PCR-based standard curve analysis using a primer pair targeting the V5 hyper-variable region included on the bacterial 16S ribosomal RNA (forward primer: 5′-GGATTAGATACCCTGGTA-3′, reverse primer: 5′-CCGTCAATTCMTTTRAGTTT-3′) to standardize the amount of potential bacterial DNA distributed within each 14 mDNA samples to equal condition (Table 3, Fig. 3). In this process, we normalized the concentration of all mDNAs to 10ng, and then we're going to confirm the uniformity of concentration and reliability of the experiment through standard curve graph analysis. We performed standard curve analysis about expected bacterial DNA quantity by diluting each normalized mDNA through a 10-fold serial dilution method (101, 100, 10-1, 10-2, respectively). The qRT-PCR experimental condition of standard curve analysis for 14 mDNA samples with bacterial 16S V5 primer pair was as follows: pre-denaturation 95.0°C 10 min, denaturation 95.0°C, 10 s, annealing/Flour detection 60°C, 15 s, elongation 72.0°C, 15 s, and total PCR cycle was 45. As a result, we confirmed that the average Ct values reflecting the potential bacterial DNA concentration for each dilution factor identified within the control group were 10 ng, 25.24; 1 ng, 28.67; 0.1 ng, 31.82; and the r2 value about the trend line calculated in the standard curve graph was 0.94. In the test group, we confirmed that the average Ct values for each diluted DNA sample were 10 ng, 25.72; 1 ng, 28.70; 0.1 ng, 30.81; and the r2 value was 0.9166. In the case of the ‘Test-7’ sample, considering that the concentration and quality of mDNA were low and unstable at 0.82 ng (A260/280: 2.37, A260/230: 0.01), we judged that it was difficult to quantify bacterial DNA by qRT-PCR because there were few potential bacterial communities present in this sample. Therefore, through this standard curve analysis, we could confirm that the potential bacterial DNA quantities within the 13 mDNA samples, except for the ‘Test-7’ sample, were equally normalized and suitable for detection and quantification analysis for S. aureus.
Evaluation of absolute and relative quantification effect for S. aureus via qRT-PCR
We evaluated the absolute and relative quantification effects of S. aureus in the skin microbiome sample through qRT-PCR using S. aureus–specific primer and probe set. Before the experiment, we confirmed that the distilled water used for making the qRT-PCR mixture solution was free of microbial contamination by checking the Ct value: N/A result reflected in the non-template DNA control (NTC). The qRT-PCR amplification conditions applied in this experiment were as follows: template DNA (mDNA) 10 ng; probe concentration 20 pmol, GreA primer concentration 25 pmol (this concentration condition was based on previous experimental results showing the most optimal running condition). As a result, we were not able to detect and quantify S. aureus within both comparison groups, which was not the expected result for this experiment (Table 4). Bacterial DNA normalization results confirmed by standard curve analysis demonstrated that the bacterial DNA in the skin microbiome sample was normalized to approximately the same amount, but due to the unstable quality of the mDNA (low concentration and low purity) and the low frequency of potential S. aureus expected to be present in it. Therefore, we judged that absolute and relative quantification of S. aureus within mDNA samples using the qRT-PCR method could not be meaningfully performed.
Quantification of S. aureus within each skin microbiome sample using LOAA digital PCR
We evaluated the applicability of the digital PCR platform to the microbial molecular diagnosis field by validating the detection and quantification effect for a proportion of S. aureus, which is expected to be present in each skin microbiome sample set in this study, using LOAA digital PCR equipment (Table 5, Fig. 4). The experimental group was set in the same as the process of detecting and quantifying S. aureus through the qRT-PCR method. Before conducting the digital PCR reaction for each comparison group, we confirmed that the number of expected copies of the greA gene within S. aureus genomic DNA in the PC1 and NC groups was 69,109 copies/µL and 0.00 copies/µL, respectively. Next, the average copies number of S. aureus greA gene identified in the control and test groups was 6.58 copies/µL and 36.28 copies/µL, respectively. However, in the case of the ‘Tes-7’ sample identified in the qRT-PCR results, it was confirmed that a low number of gene copies was detected, unlike the other mDNA samples in the test group, similar to the ‘16S V5 standard curve analysis’ result (Table 3), which was evaluated to have little bacterial genomic DNA density within this template DNA sample. Additionally, we could confirm that the potential bacterial frequency of S. aureus within the mDNA samples isolated from lesion skin sites of atopic patients was relatively high compared to the non-lesion skin sites. Considering previous studies showing that the dominant rate of S. aureus identified in atopic patients' lesions sites is higher than that of non-lesion sites, we could confirm the reliability of the S. aureus quantitative analysis results between each comparison group derived through LOAA digital PCR [14].
Discussion
The digital PCR platform is spotlighted as third-generation PCR equipment that can absolutely quantify the copy number of a particular gene to be identified within the sample without calculating a separate standard curve [15]. However, due to the high-cost burden of consumables of digital PCR, consumers’ careful consideration of the quality of DNA samples or primer and probe suitability is required. Therefore, the present study evaluated the application possibility of microbial molecular diagnosis of digital PCR platform by confirming the effect of detection and quantification for specific microorganisms within mDNA samples in unstable conditions. First, we selected ‘Staphylococcus aureus’ as a particular microbe, and by designing the S. aureus–specific primer and probe set, we successfully evaluated its high experimental accuracy for targeted-species specificity through comprehensive molecular genetic experimental validation. Additionally, since we designed the primer pair considering the strain level of S. aureus, we suggest that the designed primer can detect a broad range of strain-level single strains present in various types of clinical samples. As a result of comparing the detection and quantification efficiency of S. aureus between qRT-PCR and digital PCR, unlike qRT-PCR, which could not confirm the relative proportion of S. aureus in any mDNA sample, we confirmed that the copy number of S. aureus greA gene was calculated in digital PCR. S. aureus inhabits human skin sites at a high proportion, with an average frequency of 20%–30% (approximately 50%–60% for atopic patients) [16]. However, we judged that detecting S. aureus, even with qRT-PCR with high fluorescence sensitivity, was challenging because the quality and quantity of the mDNA samples applied in this study were very low. Despite these conditions, we estimated that the semiconductor-based LOAA digital PCR equipment was capable of S. aureus greA gene amplification not detected by qRT-PCR because specific gene amplification is possible for each bacterial genomic DNA within each over 20,000 nano-size PCR well. Additionally, the applicability of digital PCR to microbial-related diagnosis and clinical research for particular diseases was confirmed by validating the relative frequency difference of S. aureus between atopic patients' lesions and non-lesions sites.
In summary, we evaluated in this study that digital PCR has a simpler experimental process spent quantifying and detecting specific genes compared to qRT-PCR and is significantly less time-consuming (qRT-PCR takes about 5 h, digital PCR takes about 1 h). Additionally, we verified that digital PCR has excellent performance on equipment in that it enables its independent absolute quantification application. Above all things, we could confirm that gene amplification, detection, and quantification are possible using digital PCR, even within DNA samples with unstable quality and quantity, such as skin-derived mDNA. However, digital PCR has high prices for semiconductor-based PCR response detection chips, limitations of not having more than three multiple fluorescent channel functionality support in equipment, and experimental inefficiency of ‘only one sample per equipment operation.’ However, considering the digital PCR's operating excellence in deriving high-quality output data, we suggest that digital PCR has high potential application value in the microbiome-based human healthcare field if it is used for cross-validation for quantification and detection for specific microorganisms to be identified.
Notes
Authors’ Contribution
Conceptualization: DGL, DKY,KH. Data curation: DK, JJ, MK, JJ. Formal analysis: DK, Junhyeon Jeon, MK, YMH, Jinuk Jeong. Funding acquisition: KH. Methodology: DK, Junhyeon Jeon, MK, YMH, Jinuk Jeong. Writing - original draft: DK, Junhyeon Jeon, MK, Jinuk Jeong. Writing - review & editing: DGL, DKY, KH.
Conflicts of Interest
No potential conflict of interest relevant to this article was reported.
Availability
Data analysis related to qRT-PCR and digital PCR experiments was performed through the official computational applications of Bio-Rad (CFX-96) and Optolane (LOAA digital PCR) companies, and is available at https://www.bio-rad.com/https://optolane.com/ko.
Acknowledgements
This research was supported by Basic Science Research Capacity Enhancement Project through Korea Basic Science Institute (National Research Facilities and Equipment Center) grant funded by the Ministry of Education (Grant No. 2019R1A6C1010033). This study has been done with the participation of undergraduate students working at total-omics analysis research institute of Dankook University. The research institute has been supported by the VIP system as a part of Support Program for University Development 2023 of Dankook University. Following are results of a study on the "Leaders in INdustry-university Cooperation 3.0" Project, supported by the Ministry of Education and National Research Foundation of Korea.
Supplementary Materials
Supplementary data can be found with this article online at http://www.genominfo.org.