- 1Department of Clinical Laboratory, China-Japan Friendship Hospital, Beijing, China
- 2Laboratory of Clinical Microbiology and Infectious Diseases, Department of Pulmonary and Critical Care Medicine, China-Japan Friendship Hospital, Beijing, China
- 3National Center for Clinical Laboratories, Beijing Hospital, National Center of Gerontology, Beijing, China
- 4Department of Blood Transfusion, China-Japan Friendship Hospital, Beijing, China
Introduction: The severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) and influenza viruses can cause respiratory illnesses with similar clinical symptoms, making their differential diagnoses challenging. Additionally, in critically ill SARS-CoV-2–infected patients, co-infections with other respiratory pathogens can lead to severe cytokine storm and serious complications. Therefore, a method for simultaneous detection of SARS-CoV-2 and influenza A and B viruses will be clinically beneficial.
Methods: We designed an assay to detect five gene targets simultaneously via asymmetric PCR-mediated melting curve analysis in a single tube. We used specific probes that hybridize to corresponding single-stranded amplicons at low temperature and dissociate at high temperature, creating different detection peaks representing the targets. The entire reaction was conducted in a closed tube, which minimizes the risk of contamination. The limit of detection, specificity, precision, and accuracy were determined.
Results: The assay exhibited a limit of detection of <20 copies/μL for SARS-CoV-2 and influenza A and <30 copies/μL for influenza B, with high reliability as demonstrated by a coefficient of variation for melting temperature of <1.16% across three virus concentrations. The performance of our developed assay and the pre-determined assay showed excellent agreement for clinical samples, with kappa coefficients ranging from 0.98 (for influenza A) to 1.00 (for SARS-CoV-2 and influenza B). No false-positive, and no cross-reactivity was observed with six common non-influenza respiratory viruses.
Conclusion: The newly developed assay offers a straightforward, cost-effective and nucleic acid contamination-free approach for simultaneous detection of the SARS-CoV-2, influenza A, and influenza B viruses. The method offers high analytical sensitivity, reliability, specificity, and accuracy. Its use will streamline testing for co-infections, increase testing throughput, and improve laboratory efficacy.
1 Introduction
Infection with severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) causes coronavirus disease 2019 (COVID-19), which has become a sustained global pandemic (Lai et al., 2020; Sharma et al., 2020). As of May 4, 2023, the COVID-19 pandemic has affected nearly every country, with more than 752 million cases and 6.9 million deaths reported worldwide (World Health Organization, 2023). Individuals infected with SARS-CoV-2 can exhibit a range of clinical symptoms, from an asymptomatic course to severe acute respiratory distress syndrome, with most individuals showing mild symptoms or remaining asymptomatic (Gao et al., 2021). Although symptomatic individuals are a common source of transmission, undiagnosed asymptomatic cases can transmit the virus, thus increasing the prevalence of COVID-19 infection (Gao et al., 2021). The clinical symptoms of SARS-CoV-2 are very similar to those caused by the influenza A and B viruses, including fever, cough, dyspnea and fatigue (Mayuramart et al., 2021; Havasi et al., 2022). However, the prevention and treatment strategies for different viral infections are distinct. Several clinical studies have reported that co-infection with SARS-CoV-2 and one or more influenza viruses or other pathogens, particularly among severe COVID-19 cases, can exacerbate patients’ condition (Nowak et al., 2020; Wu et al., 2020; Huang et al., 2021). Thus, rapid detection and accurate identification of SARS-CoV-2 and other influenza viruses are critical for providing precise treatment, controlling the infection source and implementing effective prevention measures. Accordingly, an accurate and convenient assay that can simultaneously detect SARS-CoV-2 and influenza A and B viruses will be of great clinical benefit.
The feasibility of a detection system to be applied in clinical laboratories is influenced by various factors, including instrument cost, laboratory space, complexity, testing capability, handling time, and result time (Burd, 2010). An efficient detection method offers enhanced detection capability and reduces the burden on the detection system, which helps free up staff resources and reduce turnaround times. Generally, virus isolation based on cell culture is the “gold standard” for diagnosing viral diseases, but it is time-consuming and complex (Chen et al., 2020). At present, numerous nucleic acid-based molecular detection methods have been developed for the simultaneous detection of SARS-CoV-2 and influenza viruses. Reverse transcription-polymerase chain reaction (RT-PCR) is a widely used method that provides high sensitivity, specificity, and good repeatability (Jokela et al., 2010; Tang et al., 2020; Zhang et al., 2020). To date, several commercial detection kits based on RT-PCR have been developed and used, such as the CDC’s Influenza SARS-CoV-2 (Flu SC2) Multiplex Assay (Shu et al., 2021), the Roche Cobas® SARS-CoV-2 and Influenza A/B Assay (Kosai et al., 2022), and Cepheid’s Xpert Xpress SARS-CoV-2/Flu/RSV (Leung et al., 2021). The use of multiple detection methods based on RT-PCR mainly depends on two detection modes. The first mode is uses different fluorophore-labeled probes that target various genes, and the results are determined according to the colors of fluorescent probes. However, this mode has the disadvantage of crosstalk between fluorescence spectra, which can reduce the accuracy of the results, and this method requires labeling with fluorescent dye probes of different colors and reacting in detection instruments with multiple fluorescent channels, which are associated with higher costs. The second mode involves dividing a sample into several reaction units in which different target genes will be amplified independently. However, this strategy requires a greater amount of viral RNA input and increases the workload of laboratory staff, thus reducing the efficiency of detection.
Except for RT-PCR, reverse-transcription loop-mediated isothermal amplification (RT-LAMP) is another molecular technique to diagnose COVID-19 that grabbed the attention of many scientists due to its simplicity and high specificity (Aoki et al., 2021; de Oliveira Coelho et al., 2021). It is a rapid, and cost-effective technique that diagnoses viral diseases by targeting and amplifying six regions in the cDNA, using 4–6 primers (Notomi et al., 2015). However, the high number of primers increases the possibility of forming primer dimers, leading to false-positive results. So, the main limitation lies in primer design, a time-consuming and complex process requiring high expertise (Alhamid et al., 2022). Droplet digital PCR (ddPCR) was previously developed to reduce false-negatives results limitation reported in RT-PCR because it can accurately detect the virus in samples with low viral load (Hindson et al., 2011). Nevertheless, ddPCR requires expensive instrumentations and is time-consuming (Falzone et al., 2020; Behera et al., 2021). Clustered regularly interspaced short palindromic repeat (CRISPR)-Cas system has been developed to detect influenza or SARS-CoV-2 infections (Broughton et al., 2020). CRISPR-based technology offers advantages due to its simplicity, speed, and cost-effectiveness. It can provide results within 30 min up to 1 h and does not require expensive instruments (Broughton et al., 2020; Aman et al., 2021). However, the possible “off-target” phenomenon can affect the judgment of the test results. The non-specific collateral cleavage of Cas12 and Cas13 systems may influence other target pathogens, which can make it challenging to develop multiplex detection (Zhang et al., 2022). Whole-genome sequencing methods have been utilized to determine the origin of SARS-CoV-2 (Ren et al., 2020) and monitor the emergence of mutations and new variants (Harilal et al., 2020; John et al., 2021). However, it requires complex equipment and highly trained personnel. The high cost and long assay time make it unsuitable for large-scale testing. Moreover, sequencing errors may occur due to a large number of reads or low viral loads in clinical samples (Slatko et al., 2018). Therefore, a new method with high accuracy, efficiency and cost-effectiveness is needed for simultaneous detection of SARS-CoV-2 and influenza A/B viruses.
Multiplex asymmetric PCR based detection has evolved a rapid and convenient tool in wide variety of bioanalyses (Zhu et al., 2007; Yu et al., 2022; Kong et al., 2023). It uses many pairs of unequal concentrations of the primers for amplification, generating single stranded amplicons. This facilitates the simultaneous amplification of many targets of interest in one reaction, thus enhancing assay throughput and allowing more efficient use of each sample. Previous study has reported the development of multiplex polymorphism detection in a single tube for screening folate metabolism genes (Yu et al., 2022). Additionally, it has been utilized to identify seven immune-escape RBD mutations of Omicron (Kong et al., 2023), as well as detecting drug resistance genes in Enterobacteriaceae (Zhu et al., 2007). In the current study, we aimed to employ multiplex reverse transcription (RT) asymmetric PCR combined with melting curve analysis (MCA) to simultaneously detect SARS-CoV-2, influenza A/B viruses in a single reaction, which could be applied in poorly equipped primary hospitals and laboratories, as well as in scenarios demanding urgent diagnoses.
2 Materials and methods
2.1 Primer and probe design
The one-pot RT-asymmetric PCR-combined MCA assay was designed to target the ORF1ab gene and nucleocapsid (N) gene for the detection of SARS-CoV-2. Detection of influenza A and B was achieved by targeting the matrix (M) and hemagglutinin (HA) genes, respectively. The human RPP30 gene was selected as an internal reference and can serve as an indicator of the specimen quality. Primer and probe sequences targeted highly conserved regions of the influenza A, B and SARS-CoV-2 genome and were based on the published literature (van Elden et al., 2001; Duchamp et al., 2010; Dunn et al., 2014; Wei et al., 2020). The SARS-CoV-2 ORF1ab gene and N gene probe, with a carboxy-X-rhodamine (ROX) reporter/QSY quencher, and probes for the influenza A M gene, influenza B HA gene, and human RPP30 gene probe, with a carboxyfluorescein (FAM) reporter/QSY quencher, were purchased from Sangon (Sangon Biotech, Shanghai, China). Five pairs of primers and probes were used to simultaneously detect the three viruses in a single-tube. The optimal sequences of primers and probes are listed in Table 1.
Table 1. Primers and probes for the detection of SARS-CoV-2, influenza A, influenza B in the RT- asymmetric PCR-combined MCA assay.
2.2 Simulated single-stranded DNA (ssDNA) synthesis and viral RNA extraction
Simulated ssDNAs were designed to simulate the single-stranded amplicons that are generated through asymmetric PCR and were used to simulate the MCA process of the single-stranded amplicon and probe and to evaluate the melting temperature (Tm) value of the detection peak. The simulated ssDNAs were purchased from Sangon (Sangon Biotech) and are listed in Supplementary Table 1.
RNA templates were extracted from MS2-based virus-like particles (MS2-VLPs) corresponding to different viruses as well as the collected clinical samples. The specific virus sequences packaged in MS2-VLPs were quantified by digital PCR, and the results were regarded as the required positive standards. Viral RNA was extracted and purified with the Gene Rotex96 nucleic acid extraction and purification system (Tianlong Science and Technology Co., Ltd., Xi’an, China) using an RNA extraction kit (Tianlong Science and Technology Co., Ltd.) in accordance with the manufacturer’s instructions. The concentration and purity of the extracted RNA samples were assessed using a NanoDrop™ 2000 spectrophotometer (Thermo Fisher Scientific, Waltham, MA, United States). The qualified RNA samples were subsequently frozen and stored at −80°C.
2.3 One-pot RT-asymmetric PCR-combined MCA assay
The assay was performed on the SLAN-96P real-time PCR system (Shanghai Hongshi Medical Technology Co., Ltd). The Abstart One-step RT-PCR Mix (Sangon Biotech, Shanghai, China) was used to reverse transcription and PCR amplification. Each 50-μL PCR mixture contained 2 μL Abstart Taq® Polymerase with dNTP (0.125 U/μL) supplied in 5 × One-step RT-PCR Buffer (with 3.3 mM Mg2+), 5 × Solution I buffer, with 0.01–0.2 μM limiting primers, 0.08–1.6 μM excess primers, 0.005–0.2 μM probes, and 15 μL RNA template. The concentration details for the primers and probes are presented in Supplementary Table 2. The reactions were incubated at 42°C for 30 min and 95°C for 5 min followed by 50 cycles of 94°C for 30 s, 60°C for 60 s, and 72°C for 30 s, and then incubated at 72°C for 10 min, 95°C for 2 min, and 30°C for 2 min. The reaction mixture was subjected to MCA from 30 to 90°C and finally incubated at 40°C for 20 s. To avoid contamination of the PCR products, the whole reaction is performed in a closed PCR amplification tube that is never opened during assay performance.
2.4 Limits of detection (LoDs)
MS2-VLPs have been widely used as quality control materials for detecting pathogenic RNA viruses (Sun et al., 2013; Wang et al., 2015a; Zhang et al., 2016) due to their similarity in structure to natural viruses, stability, RNase resistance, non-infectivity and inability to replicate by itself both in vivo and in vitro (Sun et al., 2013; Wang et al., 2015a,b). Previous studies showed that MS2-VLPs were used to evaluate the LOD and precision of 7 commercial real-time PCR kits for Zaire ebolavirus (Wang et al., 2015b). Recent studies showed that SARS-CoV-2-MS2-VLPs was used to assess the LOD of rRT-PCR assays when detecting SARS-CoV-2 variants (Chen et al., 2022), as well as be used to evaluate the external quality assessment for molecular detection of SARS-CoV-2 in clinical laboratories (Wang et al., 2021). In our study, the analytical sensitivity was determined quantitatively using MS2-VLPs that contained the SARS-CoV-2 ORF1ab gene and N gene, the influenza A M gene, and the influenza B HA gene, separately. For each type of MS2-VLPs, the concentration was calibrated by digital PCR to 1.0 × 103 copies/μL. The three types of MS2-VLPs were mixed equally and serially diluted 10-fold with TE buffer to 500, 250, 125, 62.5, 31.3, 15.6, and 7.8 copies/μL. Twenty replicates of each dilution were tested and the lower LoD was defined as the concentration in copies/μL of the lowest dilution that could be detected with 95% probability and determined by probit analysis.
2.5 Specificity
The specificity of each primer and probe oligonucleotide sequence was evaluated by conducting a BLAST analysis1 using the sequences of SARS-CoV-2, influenza A, and influenza B against the nr/nt database of the National Center for Biotechnology Information and the Global Initiative on Sharing All Influenza Database2 to ensure that the primers and probes accurately matched the target gene sequence. Supplementary Figure 1 showed supporting information to better present the blastn analysis.
To determine the specificity of the developed assay, a species-specific sample panel was created that included clinical samples containing other respiratory viruses from patients who had clinical symptoms overlapping with those of influenza infection and COVID-19 (n = 58). The clinical samples were obtained from the Department of Respiratory and Critical Care Medicine at China-Japan Friendship Hospital (Beijing, China) and this study was approved by the Institutional Review Board of the China-Japan Friendship Hospital (No. 2022-KY-248), China. The sample panel contained human rhinovirus (n = 25), human parainfluenza virus (n = 12), adenovirus (n = 4), respiratory syncytial virus (n = 6), and human bocavirus (n = 1), all of which were positively detected via their respective RT-PCR assays. Approximately, 1 ng/μL of each viral RNA (15 μL) was used in the assay.
2.6 Precision
The precision of the developed assay was assessed in terms of both intra-assay precision (repeatability) and inter-assay precision (reproducibility). The variability of the developed assay was evaluated by detecting three concentrations (500, 125, and 31.3 copies/μL) of equally mixed positive MS2-VLPs. To observe intra-assay variability, each concentration was analyzed five times in one reaction. To observe inter-assay variability, each concentration was analyzed five times in independent reactions performed by different users on separate days. The variability of the Tm value was evaluated using variable analysis.
2.7 Accuracy
To evaluate the accuracy of the developed assay, a total of 345 clinical samples were collected and tested, including samples containing influenza A, influenza B, SARS-CoV-2, and co-infection viruses. The pre-determined testing for SARS-CoV-2 was conducted by using the Novel Coronavirus (2019-nCoV) Nucleic Acid Detection Kit [BioGerm Medical Technology Co., Ltd. (Shanghai, China)], while the influenza viruses were using the Influenza A Virus and Influenza B Virus Detection Kit [Coyote Bioscience Co., Ltd. (Beijing, China)]. The clinical samples were obtained from the Department of Respiratory and Critical Care Medicine at China-Japan Friendship Hospital (Beijing, China) and this study was approved by the Institutional Review Board of the China-Japan Friendship Hospital (No. 2022-KY-248). The positive clinical samples selected for this study were identified through pre-determined testing, with Ct values spanning the range of positivity (as indicated in Table 2). Detailed Ct values information for the positive samples were listed in Supplementary Table 3. For influenza A, a total of 111 samples were collected, including 40 positives and 71 negatives, and the specimen types included nasopharyngeal swabs (n = 73), oropharyngeal swabs (n = 38). For influenza B, 105 previously tested samples (29 positives and 76 negatives) were tested, and the specimen types included nasopharyngeal swabs (n = 78) and oropharyngeal swabs (n = 27). For SARS-CoV-2, a total of 129 previously diagnosed samples (42 positives and 87 negatives) were tested, including nasopharyngeal swabs (n = 102), and oropharyngeal swabs (n = 27). In addition, two samples with co-infection of influenza A and SARS-CoV-2 and one sample with co-infection of influenza B and SARS-CoV-2 were tested to assess the ability of the assay to detect co-infection, and the specimen types included nasopharyngeal swab (n = 1) and oropharyngeal swabs (n = 2). The detailed clinical samples information was presented in Supplementary Table 4.
2.8 Statistical analysis
Basic statistical values, including mean, standard deviation, and coefficient of variation for the mean Tm (°C) value, were calculated using Excel (Microsoft Corp., Redmond, WA, USA). The LoD for each virus was calculated using probit regression analysis, which determines the concentration at which the target is successfully detected in 95% of replicates (Miller et al., 2019). Overall percent agreement, positive percent agreement (PPA), negative percent agreement (NPA), and Cohen’s kappa coefficient with associated 95% confidence intervals (95% CI) were calculated using VassarStats online.3 Cohen’s kappa values were interpreted according to Landis and Koch (1977).
3 Results
3.1 Principle of the one-pot RT- asymmetric PCR-combined MCA assay
We developed a rapid and cost-efficient method for simultaneous detection of the SARS-CoV-2 and influenza A/B viruses. The RNA templates were extracted from both MS2-VLPs and clinical samples and then added into the prefabricated reaction solution for one-pot detection. As shown in Figure 1, the assay consisted of two steps. First, the extracted viral RNA was initially reverse transcribed into cDNA. Due to the insufficient amount of the generated cDNA, we then used cDNA as the template to amplify and generate enough amount of single stranded amplicons through asymmetric PCR. The asymmetric PCR uses unequal concentrations of the primers for amplification. In the early stage of amplification, both primers are available the double-stranded amplicons can be generated exponentially through normal PCR. When the limiting primer is depleted in the reaction mixture, only the still available excess primer is able to continue the amplification of target cDNA and thereby produces single stranded amplicons. Therefore, we generated sufficient amount of distinct detectable viral single-stranded amplicons in one reaction by increasing different viral pairs of primers. Subsequently, the designed fluorescent-labeled probes specifically bind to the corresponding viral single-stranded amplicons to produce detectable peaks by MCA. The feasibility of the method was verified using simulated ssDNAs, and the assay results can be directly visualized as a specific single peak by MCA (shown in Figure 2A). This assay can detect five targeted genes using probes labeled with two different fluorophores, allowing for the simultaneous detection of three viruses, which has many advantages, including short amplification time, low cost, easy accessibility, and high-throughput detection capability.
Figure 2. Validation performance for one-pot RT-asymmetric PCR-combined MCA assay. (A) The feasibility analysis confirmed by using simulated ssDNAs. (B) The determined LoD for each target using dilutions of SARS-CoV-2, influenza A and B virus MS2-VLPs. (C) The LoD [95% CI: (Lower bound, upper bound)] of the assay for the virus targeted genes. (D–L) The intra-assay precision and inter-assay precision assessed by three concentrations (500, 125, and 31.3 copies/μL) of equally mixed positive MS2-VLPs. (D–F) 500 copies/μL of equally mixed positive MS2-VLPs was analyzed five times in independent reactions performed by different users on three days, respectively. (G–I) 125 copies/μL of equally mixed positive MS2-VLPs was analyzed five times in independent reactions performed by different users on three days, respectively. (J–L) 31.3 copies/μL of equally mixed positive MS2-VLPs was analyzed five times in independent reactions performed by different users on three days, respectively.
3.2 Turnaround time
Considering the differences in operational time caused by various sample preparation and RNA extraction methods, our evaluated turnaround time merely refers to the RNA detection time. In a single experiment, it is feasible to simultaneously detect 96 samples, yielding results in roughly 3 h. Since only the extracted viral RNA needs to be added, the hands-on time is reduced, which in turn reduces exposure risk, particularly given the specificity of COVID-19.
3.3 LoDs
The LoD for each target was determined using probit analysis with dilutions of SARS-CoV-2 and influenza A and B virus MS2-VLPs (Table 3 and Figure 2B). The calculated LoD values were 18.37 copies/μL (95% confidence interval [CI]: 15.395–26.758 copies/μL) for the SARS-CoV-2 ORF1ab gene, 14.96 copies/μL (95% CI: 12.862–19.805 copies/μL) for the SARS-CoV-2 N gene, 15.60 copies/μL (95% CI: 12.516–30.316 copies/μL) for the influenza A M gene, and 26.70 copies/μL (95% CI: 21.304–42.783 copies/μL) for the influenza B HA gene. The SARS-CoV-2 N gene showed the lowest LoD and was the most sensitive among the four targets. Additionally, the analytical sensitivity was slightly and nominally lower for influenza B than for influenza A and SARS-CoV-2 (Figure 2C). These results indicate that the assay has high sensitivity, which supports its evaluation as a potential diagnostic tool.
3.4 Specificity
The specificity of the developed assay was verified in two ways. First, silico analysis and blastn analysis were performed, and the results showed that the primers and probes could match the target gene sequences correctly, with no evidence of non-target matches observed. Second, the assay’s specificity was evaluated by measuring cross-reactivity against 58 human respiratory clinical samples known to contain several different diagnostic respiratory viruses. Table 4 showed that the assay returned negative results for the other respiratory viruses and the blank control, with no amplification of the corresponding nucleic acids. These results showed that a specificity of 100% was achieved, with no cross-reactivity with any of the viruses tested.
3.5 Precision
All four targeted viral genes and the internal control gene were detected in the tested samples. We conducted a quantitative analysis of the coefficient of variation of the Tm value for each target by calculating the mean and standard deviation Tm values. The percent coefficient of variation (%CV) was calculated as a parameter representing intra-assay variability. Testing mixtures included 500 copies/μL (Figures 2D–F), 125 copies/μL (Figures 2G–I), and 31.3 copies/μL (Figures 2J–L) viral concentrations of the influenza A M gene, influenza B HA gene, and SARS-CoV-2 ORF1ab and N gene. The obtained intra-assay %CV values ranged from 0.12 to 1.16%. The inter-assay variability among three runs ranged from 0.34 to 0.99% for all targets, indicating that the assay achieved reliable detection across different viral loads. Supplementary Table 5 showed the detailed results of intra-assay precision and inter-assay precision.
3.6 Assay performance with clinical samples
To assess the accuracy of the developed multiplex assay, a total of 345 clinical samples were tested. The overall percent agreement, as well as positive and negative percent agreement for influenza A were estimated at 99.1% (95% CI: 94.4–99.9%), 97.5% (95% CI: 85.3–99.9%), and 100% (95% CI: 93.6–100%), respectively. Only one specimen showed discordant results between our developed assay and the pre-determined assays. This specimen exhibited a weakly positive result for influenza A in the pre-determined assay (with Ct value of Influenza A gene:34.676, Ct value of internal gene:34.191 and Ct value less than 35 regarded as positive), whereas our assay result was negative. This discrepancy may be attributed to the very low presence of viral RNA in the oropharyngeal swab sample and potential RNA degradation. For SARS-CoV-2 and influenza B, our developed assay yielded results identical to the pre-determined assay. Cohen’s kappa statistic ranged from 0.98 (influenza A) to 1.00 (SARS-CoV-2 and influenza B), indicating nearly perfect agreement between the two assays. Detailed concordance metrics between two assays were shown in Table 5. The results for clinical positive samples were presented in Figures 3A–C. Moreover, this method accurately detected multiple viruses within samples, including cases of co-infection with two different viruses, such as two cases of co-infection with influenza A and SARS-CoV-2 (Figures 3D, E) and one case with influenza B and SARS-CoV-2 (Figure 3F). Additionally, we also observed no false positive results among the 234 samples that were predetermined to be negative for SARS-CoV-2, influenza A, and influenza B.
Figure 3. Accuracy analysis of one-pot RT-asymmetric PCR-combined MCA assay. (A–C) Detection peak of Influenza A, Influenza B, and SARS-CoV-2 positive clinical samples, respectively. (D, E) Detection peak of co-infection of Influenza A and SARS-CoV-2 positive clinical samples. (F) Detection peak of co-infection of Influenza B and SARS-CoV-2 positive clinical samples.
4 Discussion
From 2019 to 2023, SARS-CoV-2 has caused a global pandemic resulting in millions of deaths. Similar to other respiratory viruses, the influenza viruses cause annual epidemics, which also result in millions of infections and deaths worldwide. The symptoms of influenza vary from mild to severe, including extra-pulmonary complications such as viral myocarditis and encephalitis. Patients infected with SARS-CoV-2 or influenza viruses demonstrate similar clinical symptoms, which complicates the differential diagnosis of COVID-19 from respiratory illnesses caused by influenza viruses (Ma et al., 2020). Nucleic acid tests are necessary for a definitive diagnosis of infection with these viruses, but sequential testing for different pathogens can delay the application of mitigation strategies. In critically ill COVID-19 patients, co-infections with other respiratory pathogens pose a concern, as recent studies suggest they can lead to serious complications such as acute respiratory distress syndrome, fulminant myocarditis, acute kidney injury, and multiple organ failure (Kim et al., 2020; Ma et al., 2020). However, antiviral agents that target influenza A and B viruses can be effective at lowering patients’ risks of developing pneumonia, requiring hospitalization, and dying while hospitalized, making the simultaneous detection of SARS-CoV-2 and influenza A/B viruses as co-infecting pathogens critically important for informing the clinical management of patients.
The multiplex real-time reverse transcription PCR is the most widely used method for simultaneous detection of influenza A, influenza B and SARS-CoV-2 (Chung et al., 2021; Ni et al., 2021; Norz et al., 2021; Pabbaraju et al., 2021). This approach relies on probes labeled with different fluorophores to identify various virus gene targets. However, due to the requirement for distinct labeled fluorescein to emit fluorescence under specific excitation wavelengths, this method needs to be tested on the multi-channel fluorescence PCR instrument, which is associated with higher hardware costs. Considering the crosstalk factor between fluorescence spectra, there should be no more than four types of fluorescein-labeled probes in a single reaction (Xu et al., 2021). Therefore, when multiple targets are detected, the samples may be divided into several reaction units in which different target genes are amplified and detected independently (Mancini et al., 2021). This may also lead to an increase in the number of reagents and labor costs. In comparison with this method, our developed assay uses two types of fluorescein-labeled probes, which can be implemented on the dual-channel fluorescence PCR instrument, reducing device costs. Additionally, it can simultaneously detect five gene targets in a single reaction, thus reducing reagent and labor costs.
In the present study, our established assay was based on asymmetric PCR-mediated melting curve method in a one-pot detection system with an uncomplicated experimental procedure. The extracted viral RNA was first reverse transcribed into cDNA as template, and different pairs of viral primers in unequal concentrations were used to amplify and generate sufficient amount of different viral single-stranded amplicons through asymmetric PCR. Then, specific probes hybridize to their corresponding single-stranded amplicons at low temperature and dissociate as the temperature is increased during the MCA, resulting in different detection peaks representing the corresponding targets. In comparison to other fluorescent PCR methods, this approach utilizes two types of fluorophore-labeled probes and can be performed in a two-channel fluorescent PCR system, reducing detection costs. Moreover, we demonstrated that the method can detect five different gene targets in a single tube, substantially enhancing the sample detection throughput. The assay was validated and demonstrated to be highly accurate, sensitive and reliable. The LoDs for SARS-CoV-2 and influenza A were <50 copies per reaction, and that for influenza B was <80 copies per reaction, which is equivalent to LoD for the current commercial RT-PCR kit. No cross-reactivity was observed between gene targets, even at high viral loads, and no cross-reactivity was observed with six other common non-influenza respiratory viruses. The intra-assay and inter-assay precision with three viral concentrations for four targeted genes ranged from 0.12 to 1.16%, suggesting similar reliability with the current RT-PCR method. The assay demonstrated remarkable consistency with the pre-determined testing, accurately detecting positive samples covering different viral load levels. In addition, the entire reaction is conducted in a closed-tube, reducing the risk of contamination. Overall, the developed assay is simple, low cost, nucleic acid contamination-free, and can simultaneously detect and differentiate five gene targets in a two-channel fluorescent PCR system, which will streamline testing for coinfections, increase testing throughput, and improve the laboratory turnaround time.
Due to the strict management regulations on clinical samples of SARS-CoV-2 during the pandemic period, the number of clinical samples were limited, especially the co-infection clinical samples, which is a limitation of our study. Besides, this method incorporates a melting curve analysis, which entails an extra step compared to the traditional RT-PCR approach, thereby rendering it slightly more time-consuming. Therefore, we will further optimize the cycling conditions for reverse transcription and multiple asymmetric PCR amplification to achieve more advantage in terms of detection time. In addition, the broad temperature interval between detection peaks is another limitation, which resulted in only three gene targets being detected using the same fluorophore-labeled probe. To overcome this limitation, we plan to refine the probe and PCR reaction conditions, with the goal of reducing the temperature interval between each detection peak. We also aim to use fewer fluorescently labeled probes to detect more gene targets in the future.
In conclusion, the method presented in this study offers several benefits for the simultaneous detection of SARS-CoV-2 and influenza A/B viruses. These include: (a) comparable analytical sensitivity (with a LoD < 50 copies per reaction for SARS-CoV-2 and influenza A, and <80 copies per reaction for influenza B); (b) satisfactory reliability (with a coefficient of variation of Tm values between intra-assay and inter-assay precision for three virus concentrations of <1.16%); (c) high consistent performance in clinical specimen validation experiments; (d) prevention of PCR product contamination (as the entire PCR reaction is conducted in a closed tube); and (e) simple operation, cost-effectiveness and high testing throughput (with the ability to detect and differentiate three respiratory viruses and the human internal control target in a single reaction using two fluorophore-labeled probes, and the ability to conduct high-throughput sample testing in a 96- or 384-well PCR analyzer). Therefore, this study presents a significant and clinically validated assay that can be implemented for the detection of SARS-CoV-2 and influenza A/B viruses. We anticipate that this assay will prove beneficial during upcoming influenza seasons when influenza may co-circulate with SARS-CoV-2, as rapid detection of co-infections can provide valuable time for local health authorities to contain transmission and enable clinicians to provide appropriate treatments.
Data availability statement
The original contributions presented in the study are included in the article/Supplementary material, further inquiries can be directed to the corresponding authors.
Ethics statement
The studies involving humans were approved by the Institutional Review Board of the China-Japan Friendship Hospital. The studies were conducted in accordance with the local legislation and institutional requirements. The participants provided their written informed consent to participate in this study.
Author contributions
All authors listed have made a substantial, direct, and intellectual contribution to the work, and approved it for publication.
Funding
This study was supported by the National High Level Hospital Clinical Research Funding, Elite Medical Professionals project of China-Japan Friendship Hospital (No. ZRJY2021-GG03).
Acknowledgments
We thank National Center for Clinical Laboratories to provide us different types of MS2-VLPs as research support. We also thank Medjaden Inc., for its assistance in the preparation of this manuscript.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2023.1249085/full#supplementary-material
Supplementary Figure 1 | The supporting information presents the blastn analysis.
Footnotes
- ^ http://blast.ncbi.nlm.nih.gov/Blast.cgi
- ^ https://www.gisaid.org
- ^ http://vassarstats.net/index.html
References
Alhamid, G., Tombuloglu, H., Rabaan, A. A., and Al-Suhaimi, E. (2022). SARS-CoV-2 detection methods: a comprehensive review. Saudi J. Biol. Sci. 29:103465. doi: 10.1016/j.sjbs.2022.103465
Aman, R., Marsic, T., Sivakrishna Rao, G., Mahas, A., Ali, Z., Alsanea, M., et al. (2021). iSCAN-V2: a one-pot RT-RPA-CRISPR/Cas12b assay for point-of-care SARS-CoV-2 detection. Front. Bioeng. Biotechnol. 9:800104. doi: 10.3389/fbioe.2021.800104
Aoki, M. N., de Oliveira Coelho, B., Goes, L. G. B., Minoprio, P., Durigon, E. L., Morello, L. G., et al. (2021). Colorimetric RT-LAMP SARS-CoV-2 diagnostic sensitivity relies on color interpretation and viral load. Sci. Rep. 11:9026. doi: 10.1038/s41598-021-88506-y
Behera, B. C., Mishra, R. R., and Thatoi, H. (2021). Recent biotechnological tools for diagnosis of corona virus disease: a review. Biotechnol. Prog. 37:e3078. doi: 10.1002/btpr.3078
Broughton, J. P., Deng, X., Yu, G., Fasching, C. L., Servellita, V., Singh, J., et al. (2020). CRISPR-Cas12-based detection of SARS-CoV-2. Nat. Biotechnol. 38, 870–874. doi: 10.1038/s41587-020-0513-4
Burd, E. M. (2010). Validation of laboratory-developed molecular assays for infectious diseases. Clin. Microbiol. Rev. 23, 550–576. doi: 10.1128/CMR.00074-09
Chen, Q., He, Z., Mao, F., Pei, H., Cao, H., and Liu, X. (2020). Diagnostic technologies for COVID-19: a review. RSC Adv. 10, 35257–35264. doi: 10.1039/d0ra06445a
Chen, Y., Han, Y., Yang, J., Ma, Y., Li, J., and Zhang, R. (2022). Impact of SARS-CoV-2 variants on the analytical sensitivity of rRT-PCR assays. J. Clin. Microbiol. 60:e0237421. doi: 10.1128/jcm.02374-21
Chung, H. Y., Jian, M. J., Chang, C. K., Lin, J. C., Yeh, K. M., Chen, C. W., et al. (2021). Novel dual multiplex real-time RT-PCR assays for the rapid detection of SARS-CoV-2, influenza A/B, and respiratory syncytial virus using the BD MAX open system. Emerg. Microbes Infect. 10, 161–166. doi: 10.1080/22221751.2021.1873073
de Oliveira Coelho, B., Sanchuki, H. B. S., Zanette, D. L., Nardin, J. M., Morales, H. M. P., Fornazari, B., et al. (2021). Essential properties and pitfalls of colorimetric reverse transcription loop-mediated isothermal amplification as a point-of-care test for SARS-CoV-2 diagnosis. Mol. Med. 27:30. doi: 10.1186/s10020-021-00289-0
Duchamp, M. B., Casalegno, J. S., Gillet, Y., Frobert, E., Bernard, E., Escuret, V., et al. (2010). Pandemic A(H1N1)2009 influenza virus detection by real time RT-PCR: is viral quantification useful? Clin. Microbiol. Infect. 16, 317–321. doi: 10.1111/j.1469-0691.2010.03169.x
Dunn, J., Obuekwe, J., Baun, T., Rogers, J., Patel, T., and Snow, L. (2014). Prompt detection of influenza A and B viruses using the BD Veritor System Flu A+B, Quidel(R) Sofia(R) Influenza A+B FIA, and Alere BinaxNOW(R) Influenza A&B compared to real-time reverse transcription-polymerase chain reaction (RT-PCR). Diagn. Microbiol. Infect. Dis. 79, 10–13. doi: 10.1016/j.diagmicrobio.2014.01.018
Falzone, L., Musso, N., Gattuso, G., Bongiorno, D., Palermo, C. I., Scalia, G., et al. (2020). Sensitivity assessment of droplet digital PCR for SARS-CoV-2 detection. Int. J. Mol. Med. 46, 957–964. doi: 10.3892/ijmm.2020.4673
Gao, Z., Xu, Y., Sun, C., Wang, X., Guo, Y., Qiu, S., et al. (2021). A systematic review of asymptomatic infections with COVID-19. J. Microbiol. Immunol. Infect. 54, 12–16. doi: 10.1016/j.jmii.2020.05.001
Harilal, D., Ramaswamy, S., Loney, T., Suwaidi, H. A., Khansaheb, H., Alkhaja, A., et al. (2020). SARS-CoV-2 whole genome amplification and sequencing for effective population-based surveillance and control of viral transmission. Clin. Chem. 66, 1450–1458. doi: 10.1093/clinchem/hvaa187
Havasi, A., Visan, S., Cainap, C., Cainap, S. S., Mihaila, A. A., and Pop, L. A. (2022). Influenza A, influenza B, and SARS-CoV-2 similarities and differences–a focus on diagnosis. Front. Microbiol. 13:908525. doi: 10.3389/fmicb.2022.908525
Hindson, B. J., Ness, K. D., Masquelier, D. A., Belgrader, P., Heredia, N. J., Makarewicz, A. J., et al. (2011). High-throughput droplet digital PCR system for absolute quantitation of DNA copy number. Anal. Chem. 83, 8604–8610. doi: 10.1021/ac202028g
Huang, B. R., Lin, Y. L., Wan, C. K., Wu, J. T., Hsu, C. Y., Chiu, M. H., et al. (2021). Co-infection of influenza B virus and SARS-CoV-2: a case report from Taiwan. J. Microbiol. Immunol. Infect. 54, 336–338. doi: 10.1016/j.jmii.2020.06.011
John, G., Sahajpal, N. S., Mondal, A. K., Ananth, S., Williams, C., Chaubey, A., et al. (2021). Next-generation sequencing (NGS) in COVID-19: a tool for SARS-CoV-2 diagnosis, monitoring new strains and phylodynamic modeling in molecular epidemiology. Curr. Issues Mol. Biol. 43, 845–867. doi: 10.3390/cimb43020061
Jokela, P., Piiparinen, H., Luiro, K., and Lappalainen, M. (2010). Detection of human metapneumovirus and respiratory syncytial virus by duplex real-time RT-PCR assay in comparison with direct fluorescent assay. Clin. Microbiol. Infect. 16, 1568–1573. doi: 10.1111/j.1469-0691.2010.03191.x
Kim, D., Quinn, J., Pinsky, B., Shah, N. H., and Brown, I. (2020). Rates of co-infection between SARS-CoV-2 and other respiratory pathogens. JAMA 323, 2085–2086. doi: 10.1001/jama.2020.6266
Kong, X., Gao, P., Jiang, Y., Lu, L., Zhao, M., Liu, Y., et al. (2023). Discrimination of SARS-CoV-2 omicron variant and its lineages by rapid detection of immune-escape mutations in spike protein RBD using asymmetric PCR-based melting curve analysis. Virol. J. 20:192. doi: 10.1186/s12985-023-02137-5
Kosai, K., Kaku, N., Horie, M., Kodama, H., Akamatsu, N., Narita, Y., et al. (2022). Clinical evaluation of a fully automated and high-throughput molecular testing system for detection of influenza virus. Virol. J. 19:188. doi: 10.1186/s12985-022-01916-w
Lai, C. C., Shih, T. P., Ko, W. C., Tang, H. J., and Hsueh, P. R. (2020). Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) and coronavirus disease-2019 (COVID-19): the epidemic and the challenges. Int. J. Antimicrob. Agents 55:105924. doi: 10.1016/j.ijantimicag.2020.105924
Landis, J. R., and Koch, G. G. (1977). The measurement of observer agreement for categorical data. Biometrics 33, 159–174.
Leung, E. C., Chow, V. C., Lee, M. K., Tang, K. P., Li, D. K., and Lai, R. W. (2021). Evaluation of the Xpert Xpress SARS-CoV-2/Flu/RSV assay for simultaneous detection of SARS-CoV-2, influenza A and B viruses, and respiratory syncytial virus in nasopharyngeal specimens. J. Clin. Microbiol. 59:e02965-20. doi: 10.1128/JCM.02965-20
Ma, S., Lai, X., Chen, Z., Tu, S., and Qin, K. (2020). Clinical characteristics of critically ill patients co-infected with SARS-CoV-2 and the influenza virus in Wuhan, China. Int. J. Infect. Dis. 96, 683–687. doi: 10.1016/j.ijid.2020.05.068
Mancini, F., Barbanti, F., Scaturro, M., Fontana, S., Di Martino, A., Marsili, G., et al. (2021). Multiplex real-time reverse-transcription polymerase chain reaction assays for diagnostic testing of severe acute respiratory syndrome coronavirus 2 and seasonal influenza viruses: a challenge of the phase 3 pandemic setting. J. Infect. Dis. 223, 765–774. doi: 10.1093/infdis/jiaa658
Mayuramart, O., Nimsamer, P., Rattanaburi, S., Chantaravisoot, N., Khongnomnan, K., Chansaenroj, J., et al. (2021). Detection of severe acute respiratory syndrome coronavirus 2 and influenza viruses based on CRISPR-Cas12a. Exp. Biol. Med. (Maywood) 246, 400–405. doi: 10.1177/1535370220963793
Miller, S., Naccache, S. N., Samayoa, E., Messacar, K., Arevalo, S., Federman, S., et al. (2019). Laboratory validation of a clinical metagenomic sequencing assay for pathogen detection in cerebrospinal fluid. Genome Res. 29, 831–842. doi: 10.1101/gr.238170.118
Ni, M., Xu, H., Luo, J., Liu, W., and Zhou, D. (2021). Simultaneous detection and differentiation of SARS-CoV-2, influenza A virus and influenza B virus by one-step quadruplex real-time RT-PCR in patients with clinical manifestations. Int. J. Infect. Dis. 103, 517–524. doi: 10.1016/j.ijid.2020.12.027
Norz, D., Hoffmann, A., Aepfelbacher, M., Pfefferle, S., and Lutgehetmann, M. (2021). Clinical evaluation of a fully automated, laboratory-developed multiplex RT-PCR assay integrating dual-target SARS-CoV-2 and influenza A/B detection on a high-throughput platform. J. Med. Microbiol. 70:001295. doi: 10.1099/jmm.0.001295
Notomi, T., Mori, Y., Tomita, N., and Kanda, H. (2015). Loop-mediated isothermal amplification (LAMP): principle, features, and future prospects. J. Microbiol. 53, 1–5. doi: 10.1007/s12275-015-4656-9
Nowak, M. D., Sordillo, E. M., Gitman, M. R., and Paniz Mondolfi, A. E. (2020). Coinfection in SARS-CoV-2 infected patients: where are influenza virus and rhinovirus/enterovirus? J. Med. Virol. 92, 1699–1700. doi: 10.1002/jmv.25953
Pabbaraju, K., Wong, A. A., Ma, R., Zelyas, N., and Tipples, G. A. (2021). Development and validation of a multiplex reverse transcriptase-PCR assay for simultaneous testing of influenza A, influenza B and SARS-CoV-2. J. Virol. Methods 293:114151. doi: 10.1016/j.jviromet.2021.114151
Ren, L. L., Wang, Y. M., Wu, Z. Q., Xiang, Z. C., Guo, L., Xu, T., et al. (2020). Identification of a novel coronavirus causing severe pneumonia in human: a descriptive study. Chin. Med. J. (Engl.) 133, 1015–1024. doi: 10.1097/CM9.0000000000000722
Sharma, A., Tiwari, S., Deb, M. K., and Marty, J. L. (2020). Severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2): a global pandemic and treatment strategies. Int. J. Antimicrob. Agents 56:106054. doi: 10.1016/j.ijantimicag.2020.106054
Shu, B., Kirby, M. K., Davis, W. G., Warnes, C., Liddell, J., Liu, J., et al. (2021). Multiplex real-time reverse transcription PCR for influenza A virus, influenza B virus, and severe acute respiratory syndrome coronavirus 2. Emerg. Infect. Dis. 27, 1821–1830. doi: 10.3201/eid2707.210462
Slatko, B. E., Gardner, A. F., and Ausubel, F. M. (2018). Overview of next-generation sequencing technologies. Curr. Protoc. Mol. Biol. 122:e59. doi: 10.1002/cpmb.59
Sun, Y., Jia, T., Sun, Y., Han, Y., Wang, L., Zhang, R., et al. (2013). External quality assessment for Avian Influenza A (H7N9) Virus detection using armored RNA. J. Clin. Microbiol. 51, 4055–4059. doi: 10.1128/JCM.02018-13
Tang, Y. W., Schmitz, J. E., Persing, D. H., and Stratton, C. W. (2020). Laboratory diagnosis of COVID-19: current issues and challenges. J. Clin. Microbiol. 58, e512–e520. doi: 10.1128/JCM.00512-20
van Elden, L. J., Nijhuis, M., Schipper, P., Schuurman, R., and van Loon, A. M. (2001). Simultaneous detection of influenza viruses A and B using real-time quantitative PCR. J. Clin. Microbiol. 39, 196–200. doi: 10.1128/JCM.39.1.196-200.2001
Wang, G., Sun, Y., Zhang, K., Jia, T., Hao, M., Zhang, D., et al. (2015a). External quality assessment of molecular detection of ebola virus in China. PLoS One 10:e0132659. doi: 10.1371/journal.pone.0132659
Wang, G., Zhang, R., Han, Y., Chang, L., Zhang, D., Zhang, L., et al. (2015b). The evaluation of 7 commercial real-time PCR kits for Zaire ebolavirus using virus-like particle-encapsulated EBOV RNA. Diagn. Microbiol. Infect. Dis. 83, 355–358. doi: 10.1016/j.diagmicrobio.2015.07.025
Wang, Z., Chen, Y., Yang, J., Han, Y., Shi, J., Zhan, S., et al. (2021). External quality assessment for molecular detection of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) in clinical laboratories. J. Mol. Diagn. 23, 19–28. doi: 10.1016/j.jmoldx.2020.10.008
Wei, Q., Wang, Y., Ma, J., Han, J., Jiang, M., Zhao, L., et al. (2020). Description of the first strain of 2019-nCoV, C-Tan-nCoV Wuhan strain–national pathogen resource center, China, 2020. China CDC Wkly. 2, 81–82.
World Health Organization (2023). WHO coronavirus (COVID-19) dashboard. Available online at: https://covid19.who.int/ (access May 4, 2023).
Wu, X., Cai, Y., Huang, X., Yu, X., Zhao, L., Wang, F., et al. (2020). Co-infection with SARS-CoV-2 and influenza A virus in patient with Pneumonia, China. Emerg. Infect. Dis. 26, 1324–1326. doi: 10.3201/eid2606.200299
Xu, Y., Chen, B., He, M., and Hu, B. (2021). A homogeneous nucleic acid assay for simultaneous detection of SARS-CoV-2 and influenza A (H3N2) by single-particle inductively coupled plasma mass spectrometry. Anal. Chim. Acta 1186:339134. doi: 10.1016/j.aca.2021.339134
Yu, W., Yao, J., and Zhang, Z. (2022). Simultaneous detection of three genotypes of gene methylene tetrahydrofolate reductase and methionine synthase reductase based on multiplex asymmetric real-time PCR-HRM biosensing. Anal. Chem. 94, 13052–13060. doi: 10.1021/acs.analchem.2c02096
Zhang, L., Hao, M., Zhang, K., Zhang, R., Lin, G., Jia, T., et al. (2016). External quality assessment for the molecular detection of MERS-CoV in China. J. Clin. Virol. 75, 5–9. doi: 10.1016/j.jcv.2015.12.001
Zhang, N., Wang, L., Deng, X., Liang, R., Su, M., He, C., et al. (2020). Recent advances in the detection of respiratory virus infection in humans. J. Med. Virol. 92, 408–417. doi: 10.1002/jmv.25674
Zhang, Y., Huang, Z., Zhu, J., Li, C., Fang, Z., Chen, K., et al. (2022). An updated review of SARS-CoV-2 detection methods in the context of a novel coronavirus pandemic. Bioeng. Transl. Med. 8:e10356. doi: 10.1002/btm2.10356
Keywords: SARS-CoV-2, influenza viruses, simultaneous detection, differentiation, melting curve analysis
Citation: Gao P, Fan Y, Kong X, Zhang R, Chen L, Jiang Y, Liu Y, Zhao M, Deng G, Cao Y and Ma L (2023) A novel assay based on DNA melting temperature for multiplexed identification of SARS-CoV-2 and influenza A/B viruses. Front. Microbiol. 14:1249085. doi: 10.3389/fmicb.2023.1249085
Received: 28 June 2023; Accepted: 20 November 2023;
Published: 19 December 2023.
Edited by:
Maurizio Sanguinetti, Catholic University of the Sacred Heart, ItalyReviewed by:
Alan Clavelland Ochioni, Federal University of Rio de Janeiro, BrazilSuman Chakraborty, Indian Institute of Technology Kharagpur, India
Copyright © 2023 Gao, Fan, Kong, Zhang, Chen, Jiang, Liu, Zhao, Deng, Cao and Ma. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Liang Ma, bGlhbmdtYTMyMUAxNjMuY29t; Yongtong Cao, Y2FveW9uZ3Rvbmc5MkBzaW5hLmNvbQ==