- 1Department for Disinfection and Infection Control, Chinese PLA Center for Disease Control and Prevention, Beijing, China
- 2School of Public Health, China Medical University, Shenyang, Liaoning, China
- 3College of Public Health, Zhengzhou University, Zhengzhou, Henan, China
Mucormycosis, an invasive fungal disease with severe consequences, poses a significant threat to immunocompromised individuals. However, the timely and accurate identification of Mucorales infection continues to present difficulties. In this study, novel detection techniques utilizing recombinase polymerase amplification (RPA) and quantitative real-time polymerase chain reaction (qPCR) were developed, specifically targeting the mitochondrial rnl gene, in order to address this challenge. The specificity of the RPA and qPCR assay was assessed by adding genomic DNAs extracted from 14 non-targeted strains, as well as human and mouse blood. No false-positive results were observed. Additionally, genomic DNAs from 13 species in five genera of order Mucorales were tested and yielded positive results in both methods. To further evaluate the sensitivity of the assays, DNAs from Rhizopus oryzae, Mucor racemosus, Absidia glauca, Rhizomucor miehei, and Cunninghamella bertholletiae were utilized, with concentrations ranging from 1 ng/μL to 1 fg/μL. The limit of detection (LoD) for the RPA assay was determined to be 1 pg., with the exception of Rhizomucor miehei which had a LoD of 1 ng. The LoD for the qPCR assay varied between 10 fg and 1 pg., depending on the specific species being tested. Sensitivity analysis conducted on simulated clinical samples revealed that the LoD for RPA and qPCR assays were capable of detecting DNA extracted from 103 and 101 colony forming units (CFU) conidia in 200 μL of blood and serum, respectively. Consequently, the real-time RPA and qPCR assays developed in this study exhibited favorable sensitivity and specificity for the diagnosis of mucormycosis.
Introduction
Mucormycosis is a perilous and extensively damaging invasive fungal disease caused by filamentous fungi belonging to the order Mucorales, primarily affecting individuals with compromised immune systems (Millon et al., 2019; Guegan et al., 2020). Additional predisposing risk factors include diabetes mellitus, recipients of hematopoietic stem cell transplants, severe traumatic injuries, burns, and hematological malignancies. The majority of mucormycosis cases arise from the inhalation of fungal sporangiospores or direct inoculation through wounds, subsequently leading to the growth of angioinvasive hyphae (Skiada et al., 2018). The order Mucorales encompasses a vast array of species, totaling more than 250 across 55 genera (Badali et al., 2021). Among these, the primary pathogens responsible for human mucormycosis are Rhizopus spp., Mucor spp., Lichtheimia spp. (formerly known as Absidia). Additionally, Rhizomucor spp., Cunninghamella spp., Apophysomyces spp., and Saksenae spp. have also been identified as causative agents (Badali et al., 2021). It is worth noting that the etiology of mucormycosis exhibits regional variation, with Rhizopus spp., Mucor spp., and Lichtheimia spp. being the most prevalent in European patients. In India, Rhizopus species are the predominant pathogens, including Apophysomyces elegans, Apophysomyces variabilis and Rhizopus homothallicus, which have also emerged as significant contributors (Skiada et al., 2018).
The incidence of mucormycosis has experienced a substantial increase in recent decades due to the expanding number of susceptible populations (Ambrosioni et al., 2010; Saegeman et al., 2010; Guinea et al., 2017). Based on calculations from the Leading International Fungal Education portal, the annual prevalence of mucormycosis worldwide, excluding India, is estimated to be approximately 10,000 cases (Chander et al., 2018). However, in India, the prevalence is approximately 70 times higher than the global data, with an incidence rate of approximately 0.14 per 1,000 population. The Corona Virus Disease 2019 (COVID-19) pandemic has resulted in a notable increase in the incidence of mucormycosis. Atul Patel’s research revealed a 2.1-fold rise in mucormycosis cases in India between September and December 2020 compared to the corresponding period in 2019 (Patel et al., 2021), with COVID-19 being identified as the primary contributing factor. A comprehensive analysis of 80 cases of COVID-19-associated mucormycosis (CAM) across 18 countries, demonstrated a mortality rate of 49%, with surviving patients experiencing significant long-term impairments, such as vision loss in 46% (Hoenigl et al., 2022) of cases. Furthermore, a recent study conducted in France reported a remarkably high mortality rate of 88% within a 12-week timeframe, underscoring the potential lethality of mucormycosis as a complication of COVID-19 (Danion et al., 2022).
The most prevalent clinical presentations of mucormycosis include rhino-orbito-cerebral, pulmonary, cutaneous, gastrointestinal, and disseminated forms. Typically, symptoms manifest abruptly, leading to rapid disease progression and a high mortality rate within a short period of time (Steinbrink and Miceli, 2021). Consequently, early detection and subsequent administration of antifungal therapy are imperative for a favorable prognosis. Presently, the diagnosis of mucormycosis primarily relies on mycological culture and histopathological examination. Under direct microscopy, the clinical specimens of Mucorales are characterized with the hyphae of a variable width (ranging from 6 to 25 μM), non-septate or pauci-septate and an irregular, ribbon-like branching, and wide-angle (90°) bifurcations (Skiada et al., 2018). The potential impact of friable hyphae damage during tissue manipulation on the low sensitivity and high rate of false negative outcomes has been identified as a contributing factor (Bougioukas et al., 2022). Additionally, the suboptimal sensitivity of culture, with only approximately 50% of specimens yielding positive results even in cases of positive microscopy, may lead to delays or missed diagnoses of infections (Skiada et al., 2020). Furthermore, clinical and radiological features are unable to differentiate between various causes of invasive fungal diseases (IFD), including Aspergillosis, though recent studies have indicated that the presence of a reverse halo sign may suggest a Mucorales infection (Baldin et al., 2018). The identification of mucormycosis is not aided by the detection of galactomannan or beta-glucan (Guegan et al., 2020). Given that a delay of 6 days in diagnosis has led to a significant increase in 30-day mortality rates from 35 to 66%, it is crucial to establish prompt and accurate methods for early detection (Werthman-Ehrenreich, 2021).
Molecular methods commonly used in research include conventional polymerase chain reaction (PCR), sequencing of the amplicons, restriction fragment length polymorphism (RFLP), quantitative real-time PCR (qPCR), and high-resolution melting (HRM) analysis of PCR products. The primers used in these methods target specific regions such as the internal transcribed spacer (ITS1 and ITS2) region (Bernal-Martínez et al., 2013), the 18S rRNA gene (Scherer et al., 2018), the 28S rRNA gene (Springer et al., 2016), the high-affinity iron permease (FTR1) gene (Nyilasi et al., 2008), the mitochondrial gene rnl (Caramalho et al., 2019), the cytochrome b gene (Hata et al., 2008), and the Mucorales-specific spore coating encoding protein (cotH) gene (Baldin et al., 2018). Numerous studies have utilized formalin-fixed, paraffin-embedded (FFPE) or fresh tissue samples, serum, blood and BALF samples (Millon et al., 2013; Bellanger et al., 2018; Scherer et al., 2018) to investigate this matter, with varying levels of sensitivity (62 to 100%) and specificity (89 to 100%). A detailed overview of molecular tools published Supplementary Table S1 provides further insight into this topic (Bialek et al., 2005; Hsiao et al., 2005; Nagao et al., 2005; Machouart et al., 2006; Kasai et al., 2008; Massire et al., 2013; Lengerova et al., 2014; Alanio et al., 2015; Salehi et al., 2016; Springer et al., 2016; Yang et al., 2016; Ino et al., 2017; Ashraf et al., 2022; Bigot et al., 2022).
Presently, there exist a minimum of three commercially available qPCR methods that have been documented for the purpose of detecting Mucorales. One such method, the MucorGenius kit, has been subjected to testing and has demonstrated commendable levels of sensitivity (ranging from 75 to 90%) and specificity (100%) when applied to clinical samples (Dannaoui, 2022). Additionally, Dolatabadi et al. (2014) have devised a rolling circle amplification (RCA) technique that relies on the ITS region for the molecular identification of six distinct species of Mucorales. This technique exhibits a rapid turnaround time of 2 h and is visualized through the utilization of agarose gel electrophoresis (Dolatabadi et al., 2014). When compared to conventional detection methods, these techniques have significantly decreased the time required and are likely to enhance outcomes by facilitating timely intervention with appropriate treatment. However, the clinical application of these techniques is restricted due to their reliance on a thermal cycler. Additionally, certain studies lack comprehensive data on the limit of detection (LoD), sensitivity, specificity, and cross reactivity. The limited number of clinically validated cases necessitates further investigation for a comprehensive evaluation of these established methods.
Recombinase polymerase amplification (RPA) is an isothermal nucleic acid amplification technique that functions within a consistent temperature range of 25–40°C, completing the amplification of the target sequence in a time frame of 10–20 min (Ghosh et al., 2018). The initiation of amplification occurs through the pairing of primers with homologous sequences on the template by a recombinase. Subsequently, the reaction is stabilized by the single strand binding protein, which binds to the displaced DNA strand and prevents primer dissociation. Finally, the extension step is carried out by the strand displacing polymerase (Wang et al., 2021). The primary methods utilized for detecting RPA products include real-time fluorescence, lateral flow strips (LFS), and agarose gel electrophoresis. The majority of necessary reagents are readily available in a dehydrated pellet form. Recent literature has demonstrated the successful application of RPA-based techniques in the identification of various fungal species, such as Candida albicans (Wang et al., 2021), Candida tropicalis (Wang et al., 2022), Candida glabrata (Wang et al., 2022), Cryptococcus neoformans/C. gattii (Ma et al., 2019), Aspergillus fumigatus (Li et al., 2022), and Fusarium graminearum (Lei et al., 2022; Liang et al., 2022; Wang et al., 2022), thereby fulfilling the demand for rapid, specific, and highly sensitive detection within a timeframe of 20–30 min (Supplementary Table S2). However, to date, there have been no published studies investigating the use of RPA for identifying fungi belonging to the order Mucorales. In this study, two detection approaches were devised for the identification of order Mucorales using probe-based real-time RPA and qPCR in order to establish a technical reference for the prompt clinical detection of mucormycosis.
Results
Specificity validation
The specificity of the RPA and qPCR assay, which targeted the mitochondrial rnl gene, was evaluated using 13 species from various genera within the order Mucorales, as well as 14 non-targeted strains (Table 1) and samples of human and mouse blood. All 13 Mucorales strains exhibited positive signals (Figures 1A,C), while the 14 non-Mucorales strains and the human and mouse blood samples did not yield amplification in either the RPA or qPCR assays (Figures 1B,D). These findings suggest that there were no instances of false-positive amplifications with non-Mucorales species in either assay, indicating that the primer and probe combination, Mt. RPA F/R/P and Mt. qPCR F/R/P, demonstrated good interspecies specificity.
Figure 1. Analytical specificity of Mucorales in real-time RPA and real-time PCR assay. Genomic DNAs from 13 species in five genera of order Mucorales and 16 non-Mucorales stains, were used as template for real-time RPA reaction mixtures (A,B) and real-time PCR reaction mixtures (C,D).
Sensitivity assay
The sensitivity of the real-time RPA and qPCR assay was assessed by conducting tests on 10-fold serial dilutions of genomic DNAs from R. oryzae, M. racemosus, A. glauca, Rh. miehei and C. bertholletiae. The real-time RPA assay demonstrated a detection limit of 1 pg./reaction for the mitochondrial rnl gene, with the exception of Rh. miehei which had a detection limit of 1 ng/reaction (Figure 2A). The detection limits of qPCR assay were 10 fg/reaction for M. racemosus and A. glauca, 100 fg/reaction for R. oryzae and C. bertholletiae, and 1 pg./reaction for Rh. miehei (Figure 2B). The presence of A. fumigatus or C. albicans genomic DNA did not impact the detection limits of either assay in any of the aforementioned strains (Supplementary Tables S3A,B).
Figure 2. Linearity and dynamic range. (A) Real-time RPA: 10−3-100 ng/μL genomic DNAs were introduced. (B) Real-time PCR: 10−5-100 ng/μL genomic DNAs were introduced.
Detection of simulated clinical samples
The simulated clinical samples comprised of a ten-fold dilution series ranging from 1 × 106 to 1 × 101 conidia of R. oryzae, M. racemosus, A. glauca, Rh. miehei, and C. bertholletiae, mixed with 200 μL of blood or serum. In terms of the simulated serum samples, the limits of detection (LoDs) of RPA were determined to be 103 CFU conidia/reaction for M. racemosus and A. glauca, 104 CFU for C. bertholletiae, 105 CFU for R. oryzae, and 106 CFU for Rh. miehei (Table 2A). Conversely, the LoDs of qPCR were 101 CFU conidia/reaction for R. oryzae, M. racemosus, A. glauca, and C. bertholletiae, and 103 CFU for Rh. miehei (Table 2B). Regarding simulated blood samples, the LoDs of RPA were determined to be 104 CFU conidia/reaction for M. racemosus, 105 CFU for R. oryzae, A. glauca and C. bertholletiae, and 106 CFU in Rh. miehei (Table 2A), while the LoDs of qPCR were found to be 102 CFU conidia/reaction for R. oryzae, M. racemosus, A. glauca, and C. bertholletiae, and 103 CFU for Rh. miehei (Table 2B).
Table 2. LoDs (CFU) of simulated clinical samples tested by (A) real-time RPA and (B) real-time PCR①.
The analysis revealed that the RPA assay exhibited lower sensitivity compared to the qPCR assay in terms of the LoDs in simulated clinical samples (Supplementary Table S4). Overall, both the RPA and qPCR assays demonstrated lower LoDs for these five species in simulated serum samples as opposed to blood samples.
Discussion
The conventional diagnostic techniques for mucormycosis, such as pathology and culture, are predominantly invasive and may not be appropriate for specific patient populations due to their reduced sensitivity. To address this limitation, molecular assays have been proposed as supplementary tools to conventional diagnostic procedures for the detection and identification of mucormycosis (Skiada et al., 2018). In this regard, Springer et al. have developed a probe-based Mucorales-specific real-time PCR assay (Muc18S), which successfully detected Mucorales DNA in all patients with probable or proven invasive mucormycosis, achieving a detection rate of 100%. The utilization of serum samples for detection allowed for an earlier diagnosis, approximately 21 days prior to the use of tissue samples (Springer et al., 2016). Millon et al. (2022) reported a sensitivity and specificity of 85.2 and 89.8%, respectively, for serum detection. Additionally, the first positive result in serum was observed 4 days before the first mycological or histological positive specimen, and 1 day prior to the initial imaging procedure (Millon et al., 2022).
Caramalho et al. (2019) conducted a study demonstrating the potential of rnl, a representative gene of Mitochondrial (mtDNA), as a promising marker for diagnostic applications. This is due to the protective nature of mitochondrial DNA against degradation, its stability across multiple generations, and its higher copy numbers compared to nuclear DNA (Caramalho et al., 2019). While nuclear ribosomal genes including ITS, 18S and 28S, appear to be less suitable as pan-Mucorales markers. In particular, The ITS region shows intraspecific variability within 3.24% of the Mucoromycotina (Thiyagaraja et al., 2022). Moreover, heterogeneity in the sequences of ITS 1 and ITS 2 has been observed in R. microspores (Baldin et al., 2018).
The objective of our research was to develop a real-time RPA and qPCR assay utilizing probes, with the aim of rapidly detecting Mucorales by targeting the mitochondrial rnl gene. The results of this assay were obtained within 20 min and 40 min, respectively. The selection of R. oryzae, M. racemosus, A. glauca, Rh. miehei and C. bertholletiae for this study was based on their clinical significance. To determine the specificity of the assay, we utilized genomic DNA from 13 species within 5 genera of Mucorales, 14 non-Mucorales strains, as well as human and mouse blood samples. The real-time RPA and qPCR assays demonstrated high specificity, as no cross-reactivity was observed with other organisms. Our findings indicate that the real-time RPA assay has a detection limit of 1 pg. genomic DNA/reaction or 103 CFU conidia/reaction in simulated clinical specimens. Similarly, the LoD of qPCR assay was down to 10 fg genomic DNA/reaction or 10 CFU conidia/reaction in simulated clinical specimens. The sensitivity of RPA was found to be lower than that of qPCR; however, the RPA reaction time was half as long. Importantly, both methods exhibited resistance to interference from other fungal genomic DNA (A. fumigatus and C. albicans), and their sensitivity remained unaffected.
As evidenced by the data presented in Table 2 and Supplementary Table S4, the LoDs for the simulated clinical samples was found to be lower than that of the control groups (saline groups), suggesting the influence of serum and blood. Furthermore, when comparing the sensitivity of simulated serum samples to that of simulated blood samples, it was observed that the latter exhibited a less effective sensitivity, potentially attributable to the complex composition of blood and the loss incurred during DNA extraction. Notably, in both genomic DNAs and simulated clinical specimens, the analytical sensitivity of Rh. miehei was significantly lower compared to the other four species.
In comparison to traditional culture methods, molecular-based assays have demonstrated the ability to provide earlier diagnosis, making them a valuable adjunct to conventional diagnostic procedures (Skiada et al., 2018). Springer et al. (2016) have described two independent probe-based real-time PCR tests for the detection of 12 Mucorales strains, with the LoDs ranging from 3 to 64 fg. Millon et al. (2013) have developed three qPCR assays for the detection of Mucor/Rhizopus, Lichtheimia, and Rhizomucor, with LoDs ranging from 3.7 to 15 fg/10 μL. Additionally, Dolatabadi et al. (2014) have developed a rolling circle amplification (RCA) method that has shown positive results for six of the most virulent species, with an LoD as low as 3.2 × 105 copies of amplicons. In comparison to other molecular-based assays utilized for the identification of Mucorales, the qPCR assay developed in this study demonstrated an LoD as low as 10 fg, while the RPA assay exhibited an LoD as low as 1 pg. Both the RPA and qPCR methods devised in this research possess a brief turnaround time, with approximately 20 and 40 min required, respectively. In contrast, other molecular-based assays necessitated 2–24 h to obtain results (Supplementary Table S1).
The primary drawback of these two assays is their inability to identify the specific genus of clinically relevant Mucorales strains. Accurate determination of Mucorales genera is crucial for guiding comprehensive diagnosis and selecting appropriate antifungal medication, which directly impacts clinical outcomes. While amphotericin B, posaconazole and isavuconazole antifungal therapies are generally effective against Mucorales strains in vitro, certain genera within this order have been linked to suboptimal prognosis in terms of clinical response, potentially due to reduced susceptibility (Badali et al., 2021). For instance, studies have demonstrated that the minimal inhibitory concentration of amphotericin B against Cunninghamella can be as high as 40% (Almyroudis et al., 2007), resulting in a higher mortality rate compared to diseases caused by other genera (71 versus 44%) (Cornely et al., 2019; Jeong et al., 2019). Arendrup et al. (2015) discovered that triazole antifungals are less effective against M. circinelloides compared to other clinically relevant Mucorales species, particularly isavuconazole.
Additionally, when compared to the qPCR assay developed in this study and other molecular detection methods, the real-time RPA assay is relatively more expensive. Despite being conducted in a category II biosafety cabinet, the risk of cross contamination between samples remains significant due to the presence of numerous amplicon (Ghosh et al., 2018). Therefore, it is imperative to implement rigorous quality control measures and contamination precautions. Through further optimization, specifically enhancing analytical sensitivity and conducting comprehensive evaluations of specificity, RPA has the potential to emerge as a promising alternative to qPCR or other isothermal methods for the rapid detection of pathogenic fungi. The absence of reliance on rapid temperature ramping, as observed in PCR, renders RPA more suitable for integration into microfluidic lab-on-chip devices and amenable to field use. Consequently, RPA exhibits promising potential for rapid detection (Wang et al., 2021).
The utilization of our real-time RPA and qPCR assay in clinical samples has yet to be established, and this will be a focal point of our future research endeavors. Recent investigations have documented instances of combined aspergillosis and mucormycosis cases (Loubet et al., 2022), particularly as secondary complications in critically ill COVID-19 patients (Johnson et al., 2021; Paul et al., 2022). Consequently, the screening for specific infections or co-infections, such as aspergillosis, presents an intriguing challenge (Muthu et al., 2022; Ramani et al., 2022).
In this study, our real-time RPA and qPCR system serves as the experimental basis for the swift identification of order Mucorales in clinical specimens, thereby enhancing the diagnosis and treatment of Mucorales infection and exhibiting promising potential for practical implementation.
Materials and methods
Mucorales strains and growth conditions
The Mucorales strains utilized in this research were obtained from the BeNa Culture Collection, China Center of Industrial Culture Collection, and the laboratory of the Department of Disinfection and Infection Control, Chinese PLA Center for Disease Control and Prevention (as indicated in Table 1). The fungi were cultivated on potato dextrose agar at a temperature of 28°C and allowed to incubated for 2–5 days.
Sporangiospores were harvested by washing the cultures with phosphate-buffered saline (PBS) and then separated from hyphal components using a 40 μm nylon cell strainer. Following centrifugation at 3,700 × g for 5 min, the sporangiospores were washed and resuspended in PBS for subsequent specificity testing. Additionally, 5 × 106 CFU sporangiospores of Rhizopus oryzae, Mucor racemosus, Absidia glauca, Rhizomucor miehei, and Cunninghamella bertholletiae (hereinafter abbreviated as R. oryzae, M. racemosus, A. glauca, Rh. miehei, and C. bertholletiae, respectively) were inoculated into 150 mL of potato dextrose medium and incubated at 28°C with agitation. Once the majority of sporangiospores had swollen, the germlings were separated from the medium by filtration through sterile gauze and then ground with liquid nitrogen for use in analytical sensitivity experiments.
Human and mouse blood
This study was carried out in accordance with the recommendations of Management of Laboratory Animal, Laboratory Animal Welfare and Ethics Committee, Academy of Military Medical Sciences, China, and the protocol was approved by this committee (IACUC-13-2016-002). Mice blood was acquired from C57BL/6 mice, which were obtained from Animal Center of Academy of Military Medical Sciences. Human simulated blood (Product No. A7930) was purchased from Beijing Solarbio Science & Technology Co., Ltd.
DNA extraction methods
The extraction of DNA from fungal cultures was conducted using the Biospin Fungus Genomic DNA Extraction Kit (Bioer Technology, Hangzhou, China), and simulated clinical samples were processed using the DNeasy® Blood & Tissue Kit (Qiagen, Madrid, Spain) according to the manufacturer’s instructions. Elution was performed using 50 ml of elution buffer. The quantification of DNA was carried out using the DeNovix Spectrophotometer (DeNovix Inc., America). The identification of ITS sequences was accomplished using the pairwise sequence alignment tool of BLAST from NCBI.
Primer and probe design
The primers and probes utilized in this study were designed based on a consensus sequence obtained from the mitochondrial rnl (encoding large-subunit-ribosomal-RNA) gene of R. oryzae, M. circinelloides and A. glauca sequences (GenBank accession numbers AY863212.1, KR809877.1, and KU196782.1). The software MEGA6 (version 6.06) was employed for this purpose. The design of primers and probes for real-time RPA and qPCR was performed using Primer Premier 5 (version 5.0.0). To ensure specificity, these primers (listed in Table 3) were subjected to BLAST analysis to confirm no cross-reactivity with human and non-targeted fungi. All primers were synthesized by Tsingke Biotechnology (Beijing, China).
RPA and real-time PCR assay methods
The RPA reaction was conducted in a 50 μL volume using the Nucleic Acid Amplification Kit (Qitian, Jiangsu, China). The reaction mixture consisted of 5 μL Buffer VI, 2.1 μL (10 μM) of each primer, 0.6 μL of the probe, and 2.5 μL of magnesium acetate. All reagents, except for the DNA template and magnesium acetate were combined in a master mix. The master mix was briefly vortexed for a few seconds and then distributed into individual reaction tubes containing a freeze-dried pellet. To each tube, 1 μL of genomic DNA or 5 μL of DNA extracted from simulated clinical samples was added. Subsequently, magnesium acetate was pipetted onto the tube lid and centrifuged into the bottom of the tube. Then, the tubes were immediately inserted into Light Cycler 96 instrument (Roche Diagnostics, Mannheim, Germany). The reaction was conducted at 40°C for 20 min.
A qPCR reaction with a volume of 25 μL was carried out in the Light Cycler 96, consisting of 12.5 μL 2 × FsatFire qPCR PreMix (TIANGEN, Beijing, China), 10 μM each primer, 10 μM probe, RNase-free water and either 2 μL genomic DNA or 5 μL DNA elution extracted obtained from simulated clinical samples. The cycling conditions encompassed an initial preincubation step at 95°C for 1 min, followed by an amplification program consisting of 40 cycles: denaturation at 95°C for 5 s and annealing/extension at 58°C for 15 s. In each experimental trial, both negative and positive controls were incorporated. The threshold for determining positivity was established as a Ct value of <38.
Analytical sensitivity and specificity evaluation
The analytical sensitivity of real-time RPA and qPCR techniques was evaluated by employing 10-fold dilutions of DNAs from various fungal species, including R. oryzae, M. racemosus, A. glauca, Rh. miehei, and C. bertholletiae, ranging from 1 fg/μL to 1 ng/μL. Each assay was conducted a minimum of three times and the average CT value was calculated to ensure the reliability of the threshold cycle value outcomes.
The specificity and potential cross-reactivity were assessed by utilizing DNA samples extracted from 13 species within the Mucorales order, encompassing 5 genera, as well as 14 non-targeted strains, and samples of human and mouse blood. The strains of order Mucorales are Rhizopus spp. (Rhizopus oryzae, Rhizopus microspores), Mucor spp. (Mucor mucedo, Mucor circinelloides, Mucor racemosus, Mucor pusillus, Mucor fragilis), Lichtheimia corymbifera, Absidia glauca, Rhizomucor spp. (Rhizomucor miehei, Rhizomucor pusillus), Cunninghamella bertholletiae, and Actinomucor elegans. The non-Mucorales stains were Aspergillus fumigatus, Aspergillus terreus, Trichoderma aureoviride, Penicillium chrysogenum, Candida albicans, Candida glabrata, Rhodotorula mucilaginosa, Cryptococcus neoformans, Alternaria alternata, Paraconiothyrium brasiliense, Fusarium oxysporum, Cladosporium cladosporioides, Acrophialophora levis, Escherichia coli, and mammalian whole blood from mice and human.
In order to assess the potential impact of contamination from other strains on the sensitivity of detection, we introduced clinically significant pathogenic fungi of A. fumigatus or C. albicans (10 ng genomic DNA) into 10-fold dilutions of R. oryzae, M. racemosus, A. glauca, Rh. miehei and C.bertholletiae genomic DNA (ranging from 1 ng/μL to10 fg/μL).
Preparation of simulated serum and blood samples
To simulate clinical patient samples, a 10-fold dilution series of 1 × 106 to 1 × 101 sporangiospores of R. oryzae, M. racemosus, A. glauca, Rh. miehei, and C. bertholletiae were introduced into 200 μL of simulated blood and serum. DNA from saline samples containing the aforementioned sporangiospores series served as the control group. All procedures were conducted within a class II biosafety cabinet to mitigate the risk of contamination arising from environmental fungal spores.
Statistical analysis
Statistical analysis was conducted using GraphPad Prism 5 software (GraphPad Software, La Jolla, CA). The data in Table 2 and Supplementary Tables S3, S4 was represented as the mean ± standard deviation of the mean from at least three independent experiments, while the data in Figure 2 was represented as the mean.
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
Ethical approval was not required for the studies on humans in accordance with the local legislation and institutional requirements because only commercially available established cell lines were used. The animal study was approved by Laboratory Animal Welfare and Ethics Committee, Academy of Military Medical Sciences, China. The study was conducted in accordance with the local legislation and institutional requirements.
Author contributions
RX: Investigation, Methodology, Project administration, Validation, Writing – review & editing, Formal analysis, Data curation, Software, Writing – original draft. DL: Formal analysis, Investigation, Methodology, Project administration, Validation, Writing – review & editing, Conceptualization, Resources, Supervision. JZ: Methodology, Project administration, Resources, Writing – review & editing. HZ: Data curation, Investigation, Software, Writing – review & editing. HC: Investigation, Methodology, Resources, Supervision, Writing – review & editing. YJ: Conceptualization, Data curation, Project administration, Software, Writing – review & editing. FC: Conceptualization, Data curation, Formal analysis, Methodology, Project administration, Resources, Supervision, Writing – review & editing. LH: Conceptualization, Funding acquisition, Investigation, Project administration, Resources, Supervision, Writing – review & editing.
Funding
The author(s) declare financial support was received for the research, authorship, and/or publication of this article. This work was supported by the Special Key Project of Biosafety Technologies (2022YFC2604000) for the National Major Research & Development Program of China, the National Natural Science Foundation of China under Grant (nos. 81971914, 82172293). The funders had no role in study design, data collection and analysis, decision to publish or preparation of manuscript.
Acknowledgments
The authors express their gratitude to the Special Key Project of Biosafety Technologies for the National Major Research & Development Program of China and the National Natural Science Foundation of China for their generous support.
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.1273073/full#supplementary-material
Abbreviations
RPA, recombinase polymerase amplification; qPCR, quantitative real-time polymerase chain reaction; LoD, limit of detection; CFU, colony forming unit; HSCT, hematopoietic stem cell transplant; COVID-19, corona virus disease 2019; CAM, COVID-19-associated mucormycosis; IFD, invasive fungal disease; PCR, polymerase chain reaction; HRM, high-resolution melting; ITS, internal transcribed spacer; FTR, high-affinity iron permease; cotH, spore coating encoding protein; BALF, bronchoalveolar lavage fluid; FFPE, formalin-fixed, paraffin-embedded tissue; RCA, rolling circle amplification; ESI-MS, electrospray-ionization mass spectrometry; LFS, lateral flow strips; PBS, phosphate-buffered saline; Ct value, threshold cycle value.
References
Alanio, A., Garcia-Hermoso, D., Mercier-Delarue, S., Lanternier, F., Gits-Muselli, M., Menotti, J., et al. (2015). Molecular identification of Mucorales in human tissues: contribution of PCR electrospray-ionization mass spectrometry. Clin. Microbiol. Infect. 21, 594.e1–594.e5. doi: 10.1016/j.cmi.2015.01.017
Almyroudis, N. G., Sutton, D. A., Fothergill, A. W., Rinaldi, M. G., and Kusne, S. (2007). In vitro susceptibilities of 217 clinical isolates of zygomycetes to conventional and new antifungal agents. Antimicrob. Agents Chemother. 51, 2587–2590. doi: 10.1128/AAC.00452-07
Ambrosioni, J., Bouchuiguir-Wafa, K., and Garbino, J. (2010). Emerging invasive zygomycosis in a tertiary care center: epidemiology and associated risk factors. Int. J. Infect. Dis. 14, e100–e103. doi: 10.1016/j.ijid.2009.11.024
Arendrup, M. C., Jensen, R. H., and Meletiadis, J. (2015). In vitro activity of Isavuconazole and comparators against clinical isolates of the Mucorales order. Antimicrob. Agents Chemother. 59, 7735–7742. doi: 10.1128/AAC.01919-15
Ashraf, M. J., Kord, M., Morovati, H., Ansari, S., Shekarkhar, G., Badali, H., et al. (2022). Evaluating a semi-nested PCR to support histopathology reports of fungal rhinosinusitis in formalin-fixed paraffin-embedded tissue samples. J. Clin. Lab. Anal. 36:e24209. doi: 10.1002/jcla.24209
Badali, H., Cañete-Gibas, C., McCarthy, D., Patterson, H., Sanders, C., David, M. P., et al. (2021). Epidemiology and antifungal susceptibilities of mucoralean Fungi in clinical samples from the United States. J. Clin. Microbiol. 59:e0123021. doi: 10.1128/JCM.01230-21
Baldin, C., Soliman, S. S. M., Jeon, H. H., Alkhazraji, S., Gebremariam, T., Gu, Y., et al. (2018). PCR-based approach targeting Mucorales-specific gene family for diagnosis of Mucormycosis. J. Clin. Microbiol. 56:e00746-18. doi: 10.1128/JCM.00746-18
Bellanger, A.-P., Berceanu, A., Rocchi, S., Valot, B., Fontan, J., Chauchet, A., et al. (2018). Development of a quantitative PCR detecting Cunninghamella bertholletiae to help in diagnosing this rare and aggressive mucormycosis. Bone Marrow Transplant. 53, 1180–1183. doi: 10.1038/s41409-018-0194-5
Bernal-Martínez, L., Buitrago, M. J., Castelli, M. V., Rodriguez-Tudela, J. L., and Cuenca-Estrella, M. (2013). Development of a single tube multiplex real-time PCR to detect the most clinically relevant Mucormycetes species. Clin. Microbiol. Infect. 19, E1–E7. doi: 10.1111/j.1469-0691.2012.03976.x
Bialek, R., Konrad, F., Kern, J., Aepinus, C., Cecenas, L., Gonzalez, G. M., et al. (2005). PCR based identification and discrimination of agents of mucormycosis and aspergillosis in paraffin wax embedded tissue. J. Clin. Pathol. 58, 1180–1184. doi: 10.1136/jcp.2004.024703
Bigot, J., Godmer, A., Prudenté, L., Angebault, C., Brissot, E., Bige, N., et al. (2022). Diagnosis of mucormycosis using an intercalating dye-based quantitative PCR. Med. Mycol. 60:myac015. doi: 10.1093/mmy/myac015
Bougioukas, L., Noyes, C. D., Peterson, K., Smith, L. M., and Hale, A. J. (2022). Successful treatment of rhino-cerebral mucormycosis with dual antifungal therapy and minimal surgical debridement. IDCases. 30:e01615. doi: 10.1016/j.idcr.2022.e01615
Caramalho, R., Madl, L., Rosam, K., Rambach, G., Speth, C., Pallua, J., et al. (2019). Evaluation of a novel mitochondrial Pan- marker for the detection, identification, quantification, and growth stage determination of Mucormycetes. J Fungi (Basel). 5:98. doi: 10.3390/jof5040098
Chander, J., Kaur, M., Singla, N., Punia, R. P. S., Singhal, S. K., Attri, A. K., et al. (2018). Mucormycosis: Battle with the deadly enemy over a five-year period in India. J Fungi (Basel). 4:46. doi: 10.3390/jof4020046
Cornely, O. A., Alastruey-Izquierdo, A., Arenz, D., Chen, S. C. A., Dannaoui, E., Hochhegger, B., et al. (2019). Global guideline for the diagnosis and management of mucormycosis: an initiative of the European Confederation of Medical Mycology in cooperation with the mycoses study group education and research consortium. Lancet Infect. Dis. 19, e405–e421. doi: 10.1016/S1473-3099(19)30312-3
Danion, F., Letscher-Bru, V., Guitard, J., Sitbon, K., Dellière, S., Angoulvant, A., et al. (2022). Coronavirus disease 2019-associated Mucormycosis in France: a rare but deadly complication. Open Forum Infect. Dis. 9:ofab566. doi: 10.1093/ofid/ofab566
Dannaoui, E. (2022). Recent developments in the diagnosis of Mucormycosis. J Fungi (Basel). 8:457. doi: 10.3390/jof8050457
Dolatabadi, S., Najafzadeh, M. J., and de Hoog, G. S. (2014). Rapid screening for human-pathogenic Mucorales using rolling circle amplification. Mycoses 57, 67–72. doi: 10.1111/myc.12245
Ghosh, D. K., Kokane, S. B., Kokane, A. D., Warghane, A. J., Motghare, M. R., Bhose, S., et al. (2018). Development of a recombinase polymerase based isothermal amplification combined with lateral flow assay (HLB-RPA-LFA) for rapid detection of "Candidatus Liberibacter asiaticus". PLoS One 13:e0208530. doi: 10.1371/journal.pone.0208530
Guegan, H., Iriart, X., Bougnoux, M.-E., Berry, A., Robert-Gangneux, F., and Gangneux, J.-P. (2020). Evaluation of MucorGenius® mucorales PCR assay for the diagnosis of pulmonary mucormycosis. J. Infect. 81, 311–317. doi: 10.1016/j.jinf.2020.05.051
Guinea, J., Escribano, P., Vena, A., Muñoz, P., Martínez-Jiménez, M. D. C., Padilla, B., et al. (2017). Increasing incidence of mucormycosis in a large Spanish hospital from 2007 to 2015: epidemiology and microbiological characterization of the isolates. PLoS One 12:e0179136. doi: 10.1371/journal.pone.0179136
Hata, D. J., Buckwalter, S. P., Pritt, B. S., Roberts, G. D., and Wengenack, N. L. (2008). Real-time PCR method for detection of zygomycetes. J. Clin. Microbiol. 46, 2353–2358. doi: 10.1128/JCM.02331-07
Hoenigl, M., Seidel, D., Carvalho, A., Rudramurthy, S. M., Arastehfar, A., Gangneux, J.-P., et al. (2022). The emergence of COVID-19 associated mucormycosis: a review of cases from 18 countries. Lancet Microbe. 3, e543–e552. doi: 10.1016/S2666-5247(21)00237-8
Hsiao, C. R., Huang, L., Bouchara, J.-P., Barton, R., Li, H. C., and Chang, T. C. (2005). Identification of medically important molds by an oligonucleotide array. J. Clin. Microbiol. 43, 3760–3768. doi: 10.1128/JCM.43.8.3760-3768.2005
Ino, K., Nakase, K., Nakamura, A., Nakamori, Y., Sugawara, Y., Miyazaki, K., et al. (2017). Management of pulmonary mucormycosis based on a polymerase chain reaction (PCR) diagnosis in patients with hematologic malignancies: a report of four cases. Intern. Med. 56, 707–711. doi: 10.2169/internalmedicine.56.7647
Jeong, W., Keighley, C., Wolfe, R., Lee, W. L., Slavin, M. A., Kong, D. C. M., et al. (2019). The epidemiology and clinical manifestations of mucormycosis: a systematic review and meta-analysis of case reports. Clin. Microbiol. Infect. 25, 26–34. doi: 10.1016/j.cmi.2018.07.011
Johnson, A. K., Ghazarian, Z., Cendrowski, K. D., and Persichino, J. G. (2021). Pulmonary aspergillosis and mucormycosis in a patient with COVID-19. Med Mycol Case Rep. 32, 64–67. doi: 10.1016/j.mmcr.2021.03.006
Kasai, M., Harrington, S. M., Francesconi, A., Petraitis, V., Petraitiene, R., Beveridge, M. G., et al. (2008). Detection of a molecular biomarker for zygomycetes by quantitative PCR assays of plasma, bronchoalveolar lavage, and lung tissue in a rabbit model of experimental pulmonary zygomycosis. J. Clin. Microbiol. 46, 3690–3702. doi: 10.1128/JCM.00917-08
Lei, R., Li, Y., Li, L., Wang, J., Cui, Z., Ju, R., et al. (2022). A CRISPR/Cas12a-based portable platform for rapid detection of Leptosphaeria maculans in Brassica crops. Front. Plant Sci. 13:976510. doi: 10.3389/fpls.2022.976510
Lengerova, M., Racil, Z., Hrncirova, K., Kocmanova, I., Volfova, P., Ricna, D., et al. (2014). Rapid detection and identification of mucormycetes in bronchoalveolar lavage samples from immunocompromised patients with pulmonary infiltrates by use of high-resolution melt analysis. J. Clin. Microbiol. 52, 2824–2828. doi: 10.1128/JCM.00637-14
Li, Z., Wang, M., Xu, T., Zhan, Y., Chen, F., Lin, Y., et al. (2022). Development and clinical implications of a novel CRISPR-based diagnostic test for pulmonary Aspergillus fumigatus infection. J. Microbiol. Immunol. Infect. 55, 749–756. doi: 10.1016/j.jmii.2021.11.008
Liang, X., Zhang, X., Haseeb, H. A., Tang, T., Shan, J., Yin, B., et al. (2022). Development and evaluation of a novel visual and rapid detection assay for toxigenic Fusarium graminearum in maize based on recombinase polymerase amplification and lateral flow analysis. Int. J. Food Microbiol. 372:109682. doi: 10.1016/j.ijfoodmicro.2022.109682
Loubet, D., Sarton, B., Lelièvre, L., Grouteau, G., Iriart, X., Chauvin, P., et al. (2022). Fatal mucormycosis and aspergillosis coinfection associated with haemophagocytic lymphohistiocytosis: a case report and literature review. J Mycol Med. 33:101325. doi: 10.1016/j.mycmed.2022.101325
Ma, Q., Yao, J., Yuan, S., Liu, H., Wei, N., Zhang, J., et al. (2019). Development of a lateral flow recombinase polymerase amplification assay for rapid and visual detection of Cryptococcus neoformans/C. gattii in cerebral spinal fluid. BMC Infect. Dis. 19:108. doi: 10.1186/s12879-019-3744-6
Machouart, M., Larché, J., Burton, K., Collomb, J., Maurer, P., Cintrat, A., et al. (2006). Genetic identification of the main opportunistic Mucorales by PCR-restriction fragment length polymorphism. J. Clin. Microbiol. 44, 805–810. doi: 10.1128/JCM.44.3.805-810.2006
Massire, C., Buelow, D. R., Zhang, S. X., Lovari, R., Matthews, H. E., Toleno, D. M., et al. (2013). PCR followed by electrospray ionization mass spectrometry for broad-range identification of fungal pathogens. J. Clin. Microbiol. 51, 959–966. doi: 10.1128/JCM.02621-12
Millon, L., Caillot, D., Berceanu, A., Bretagne, S., Lanternier, F., Morio, F., et al. (2022). Evaluation of serum Mucorales polymerase chain reaction (PCR) for the diagnosis of Mucormycoses: the MODIMUCOR prospective trial. Clin. Infect. Dis. 75, 777–785. doi: 10.1093/cid/ciab1066
Millon, L., Larosa, F., Lepiller, Q., Legrand, F., Rocchi, S., Daguindau, E., et al. (2013). Quantitative polymerase chain reaction detection of circulating DNA in serum for early diagnosis of mucormycosis in immunocompromised patients. Clin. Infect. Dis. 56, e95–e101. doi: 10.1093/cid/cit094
Millon, L., Scherer, E., Rocchi, S., and Bellanger, A.-P. (2019). Molecular Strategies to Diagnose Mucormycosis. J Fungi (Basel). 5:24. doi: 10.3390/jof5010024
Muthu, V., Agarwal, R., Patel, A., Kathirvel, S., Abraham, O. C., Aggarwal, A. N., et al. (2022). Definition, diagnosis, and management of COVID-19-associated pulmonary mucormycosis: Delphi consensus statement from the fungal infection study forum and academy of pulmonary sciences, India. Lancet Infect. Dis. 22, e240–e253. doi: 10.1016/S1473-3099(22)00124-4
Nagao, K., Ota, T., Tanikawa, A., Takae, Y., Mori, T., Udagawa, S.-i., et al. (2005). Genetic identification and detection of human pathogenic Rhizopus species, a major mucormycosis agent, by multiplex PCR based on internal transcribed spacer region of rRNA gene. J. Dermatol. Sci. 39, 23–31. doi: 10.1016/j.jdermsci.2005.01.010
Nyilasi, I., Papp, T., Csernetics, A., Krizsán, K., Nagy, E., and Vágvölgyi, C. (2008). High-affinity iron permease (FTR1) gene sequence-based molecular identification of clinically important zygomycetes. Clin. Microbiol. Infect. 14, 393–397. doi: 10.1111/j.1469-0691.2007.01932.x
Patel, A., Agarwal, R., Rudramurthy, S. M., Shevkani, M., Xess, I., Sharma, R., et al. (2021). Multicenter epidemiologic study of coronavirus disease-associated Mucormycosis, India. Emerg. Infect. Dis. 27, 2349–2359. doi: 10.3201/eid2709.210934
Paul, M., Sasidharan, J., Taneja, J., Chatterjee, K., Abbas, S. Z., Chowdhury, V., et al. (2022). Invasive mucormycosis and aspergillosis coinfection associated with post-COVID-19 pneumonia in a tertiary care hospital. Med Mycol J. 63, 59–64. doi: 10.3314/mmj.21-00019
Ramani, P., Krishnan, R. P., Pandiar, D., Benitha, J. G., Ramalingam, K., and Gheena, S. (2022). Chronic invasive aspergillosis with fulminant Mucormycosis sparing palate in a post-COVID-19 patient - a case report. Ann Maxillofac Surg. 12, 102–105. doi: 10.4103/ams.ams_296_21
Saegeman, V., Maertens, J., Meersseman, W., Spriet, I., Verbeken, E., and Lagrou, K. (2010). Increasing incidence of mucormycosis in university hospital, Belgium. Emerg. Infect. Dis. 16, 1456–1458. doi: 10.3201/eid1609.100276
Salehi, E., Hedayati, M. T., Zoll, J., Rafati, H., Ghasemi, M., Doroudinia, A., et al. (2016). Discrimination of aspergillosis, Mucormycosis, Fusariosis, and Scedosporiosis in formalin-fixed paraffin-embedded tissue specimens by use of multiple real-time quantitative PCR assays. J. Clin. Microbiol. 54, 2798–2803. doi: 10.1128/JCM.01185-16
Scherer, E., Iriart, X., Bellanger, A. P., Dupont, D., Guitard, J., Gabriel, F., et al. (2018). Quantitative PCR (qPCR) detection of Mucorales DNA in bronchoalveolar lavage fluid to diagnose pulmonary Mucormycosis. J. Clin. Microbiol. 56:e00289-18. doi: 10.1128/JCM.00289-18
Skiada, A., Lass-Floerl, C., Klimko, N., Ibrahim, A., Roilides, E., and Petrikkos, G. (2018). Challenges in the diagnosis and treatment of mucormycosis. Med. Mycol. 56, S93–S101. doi: 10.1093/mmy/myx101
Skiada, A., Pavleas, I., and Drogari-Apiranthitou, M. (2020). Epidemiology and diagnosis of Mucormycosis: an update. J Fungi (Basel). 6:265. doi: 10.3390/jof6040265
Springer, J., Goldenberger, D., Schmidt, F., Weisser, M., Wehrle-Wieland, E., Einsele, H., et al. (2016). Development and application of two independent real-time PCR assays to detect clinically relevant Mucorales species. J. Med. Microbiol. 65, 227–234. doi: 10.1099/jmm.0.000218
Springer, J., Lackner, M., Ensinger, C., Risslegger, B., Morton, C. O., Nachbaur, D., et al. (2016). Clinical evaluation of a Mucorales-specific real-time PCR assay in tissue and serum samples. J. Med. Microbiol. 65, 1414–1421. doi: 10.1099/jmm.0.000375
Steinbrink, J. M., and Miceli, M. H. (2021). Mucormycosis. Infect. Dis. Clin. N. Am. 35, 435–452. doi: 10.1016/j.idc.2021.03.009
Thiyagaraja, V., Ertz, D., Lücking, R., Wanasinghe, D. N., Aptroot, A., Cáceres, M. E. S., et al. (2022). Taxonomic and phylogenetic reassessment of Pyrgidium (Mycocaliciales) and investigation of ascospore morphology. J Fungi (Basel). 8:966. doi: 10.3390/jof8090966
Wang, L., Fu, J., Cai, G., Cheng, X., Zhang, D., Shi, S., et al. (2022). Rapid and visual RPA-Cas12a fluorescence assay for accurate detection of dermatophytes in cats and dogs. Biosensors (Basel). 12:636. doi: 10.3390/bios12080636
Wang, F., Ge, D., Wang, L., Li, N., Chen, H., Zhang, Z., et al. (2021). Rapid and sensitive recombinase polymerase amplification combined with lateral flow strips for detecting Candida albicans. Anal. Biochem. 633:114428. doi: 10.1016/j.ab.2021.114428
Wang, K., Huo, L., Li, Y., Zhu, L., Wang, Y., and Wang, L. (2022). Establishment of a rapid diagnosis method for Candida glabrata based on the ITS2 gene using recombinase polymerase amplification combined with lateral flow strips. Front. Cell. Infect. Microbiol. 12::953302, doi: 10.3389/fcimb.2022.953302
Wang, L., Xu, A., Zhou, P., Zhao, M., Xu, C., Wang, Y., et al. (2022). Rapid detection of in clinical samples from different sources using RPA-LFS. Front. Cell. Infect. Microbiol. 12:898186. doi: 10.3389/fcimb.2022.898186
Werthman-Ehrenreich, A. (2021). Mucormycosis with orbital compartment syndrome in a patient with COVID-19. Am. J. Emerg. Med. 42, 264.e5–264.e8. doi: 10.1016/j.ajem.2020.09.032
Keywords: mucormycosis, order Mucorales, RPA, real-time PCR, COVID-19, diagnosis
Citation: Xu R, Li D, Zhao J, Zhong H, Chen H, Jia Y, Chen F and Han L (2023) Rapid detection of Mucorales based on recombinase polymerase amplification and real-time PCR. Front. Microbiol. 14:1273073. doi: 10.3389/fmicb.2023.1273073
Edited by:
Ana Cristina Sampaio, Universidade de Trás-os-Montes e Alto Douro, PortugalReviewed by:
Mahdi Abastabar, Mazandaran University of Medical Sciences, IranHua Zhou, Zhejiang University, China
Copyright © 2023 Xu, Li, Zhao, Zhong, Chen, Jia, Chen and Han. 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: Li Han, hanlicdc@163.com; Fangyan Chen, chenfangyan1986@163.com
†These authors have contributed equally to this work