- 1Co-Innovation Center for the Sustainable Forestry in Southern China, Nanjing Forestry University, Nanjing, Jiangsu, China
- 2Jiangsu Provincial Key Construction Laboratory of Special Biomass Resource Utilization, Nanjing Xiaozhuang University, Nanjing, China
Introduction: Phytophthora sojae is among the most devastating pathogens of soybean (Glycine max) and severely impacts soybean production in several countries. The resulting disease can be difficult to diagnose and other Phytophthora species can also infect soybean. Accurate diagnosis is important for management of the disease caused by P. sojae.
Methods: In this study, recombinase polymerase amplification (RPA) in combination with the CRISPR/Cas12a system were used for detection of P. sojae. The assay was highly specific to P. sojae.
Results: The test results were positive for 29 isolates of P. sojae, but negative for 64 isolates of 29 Phytophthora species, 7 Phytopythium and Pythium species, 32 fungal species, and 2 Bursaphelenchus species. The method was highly sensitive, detecting as little as 10 pg.µL−1 of P. sojae genomic DNA at 37°C in 20 min. The test results were visible under UV light and readout coming from fluorophores. In addition, P. sojae was detected from natural inoculated hypocotyls of soybean seedlings using this novel assay. The rapidity and accuracy of the method were verified using 30 soybean rhizosphere samples.
Discussion: In conclusion, the RPA-CRISPR/Cas12a detection assay developed here is sensitive, efficient, and convenient, and has potential for further development as a kit for monitoring root rot of soybean in the field.
Introduction
Phytophthora sojae is among the most destructive pathogens of soybean (Glycine max). An estimated annual worldwide loss of US$1–2 billion is caused by P. sojae, which was first identified as a soybean pathogen in the 1950s in Indiana, Ohio, and North Carolina (Schmitthenner, 1985; Tyler, 2007; Dorrance, 2018). The pathogen has gradually spread to many soybean-producing regions of the world, including Brazil, Canada, China, and Argentina, where it has become a major impediment to soybean production (Wrather and Koenning, 2006).
China is the third-largest producer of soybeans in the world, following the United States and Brazil (Zhang et al., 2010; Roese and Pereira Goulart, 2013). P. sojae was first detected in Heilongjiang Province of China in 1989 (Su and Shen, 1993; Cui et al., 2010). After assessing its potential risk to economic and agricultural security, the Ministry of Agriculture of the People’s Republic of China listed P. sojae as a quarantine pathogen in 2007. The disease has caused severe soybean losses in China and has shown a trend for gradual expansion in the country (Chen et al., 2004).
Rapid detection of P. sojae is a crucial step towards effective control of soybean root rot and seedling blight. Traditionally, methods to identify P. sojae involve direct isolation from diseased plant tissues and baiting from soil on semi-selective media (Davison, 1998). Subsequent identification based on pathogen morphology and DNA sequences are usually time-consuming and requires trained personnel. The development of a PCR-based assay for rapid and sensitive detection of P. sojae may facilitate efficient pathogen identification and lead to effective disease management. A variety of molecular detection methods for phytopathogens have been developed, including PCR (Wang et al., 2006; Bienapfl et al., 2011; Xiong et al., 2019), real-time fluorescent PCR, loop-mediated isothermal amplification (LAMP), and recombinase polymerase amplification–lateral flow dipstick (RPA-LFD) assays (Cullen et al., 2001; Cullen et al., 2002; Lees et al., 2002; Dai et al., 2012; Dai et al., 2015; Rojas et al., 2017; Dai et al., 2019). However, these detection procedures rely on thermal-cycling instruments and other specialized equipment for sample treatment and result output, which limits the application of these methods to diagnose diseases in the field and by gel electrophoresis (Cullen et al., 2001; Cullen et al., 2002; Lees et al., 2002; Filion et al., 2003; Zhu et al., 2021).
In recent years, several molecular detection technologies based on isothermal amplification reactions have been developed and applied to pathogen detection, such as nuclear acid sequence-based amplification (NASBA), rolling circle amplification (RCA), helicase-dependent amplification (HAD), LAMP, and recombinase polymerase amplification (RPA) (Polstra et al., 2002; Andresen et al., 2009; Clancy et al., 2015; Magriñá Lobato and O’Sullivan, 2018; Huang et al., 2019; Zhang et al., 2020; Xu Z. C. et al., 2022). In 2018, a new DNA detection method [DNA endonuclease targeted CRISPR trans reporter (DETECTR)] by combining CRISPR/Cas with isothermal amplification was developed (Chen et al., 2018). DETECTR uses RPA to amplify targeted double-stranded DNA (dsDNA) at a constant temperature. When Cas12a specifically recognises and cleaves the amplified product under the control of a programmable specific CRISPR RNA(crRNA), non-specific trans-cleavage activity is activated, which cleaves the single-stranded DNA (ssDNA) reporter in the system, generating a fluorescence signal. DETECTR is used for the rapid detection of viruses, and the sensitivity can reach the attomole level (Barrangou and Marraffini, 2014; Le et al., 2021; Wang et al., 2021; Yue et al., 2021).
Because the adjustable properties of CRISPR-Cas effectors, such as simplicity of design, easy operation, collateral cleavage activity, and high biocompatibility, CRISPR-Cas systems have been widely used in genome editing and transcriptional regulation (Kadam et al., 2023). The first biological evidence that CRISPR-Cas systems play a role in adaptive immunity was reported in 2007 when S. thermophilus CRISPR loci were shown to acquire novel spacers derived from the invasive phage DNA (Barrangou et al., 2007). In recent years, technologies based on CRISPR/Cas12a and Cas13a collateral cleavage activity have been used to detect COVID-19, nCOV-2019, African swine fever virus (ASFV), Zika virus (ZIKV), Dengue virus (DENV), and other important viruses (Myhrvold et al., 2018; Bai et al., 2019; Broughton et al., 2020). In addition, this novel technique can be used for rapid detection of Salmonella, Brucella, Helicobacter, and other bacteria (Dai et al., 2022; Liu et al., 2022; Xu J. H. et al., 2022). To date, the method has been used to detect not only bacteria, viruses, and mycoplasmas, but also phytopathogenic fungi causing diseases such as citus black star disease, citrus scab, and wheat blast (Shi, 2021; Shin et al., 2021; Li et al., 2022).
In this study, a novel RPA-CRISPR/Cas12a assay was developed for detection of P. sojae, targeting the Ypt1 gene, comprising 10 min for crude DNA extraction, 5 min for the RPA reaction, and 15 min for the CRISPR/Cas12a reaction. The results can be observed as green fluorescence under a blue LED transilluminator at a wavelength of 470 nm or detected using a multifunctional microplate reader (fluorescence excitation wavelength [λex] 485 nm, fluorescence emission wavelength [λem] 520 nm). The specificity of the assay was evaluated by testing against P. sansomeana, P. melonis, P. vignae, and other oomycete and fungal species. The analytical sensitivity and feasibility of this method were confirmed using artificially inoculated samples.
Materials and methods
Maintenance of isolates and DNA extraction
Twenty-nine isolates of Phytophthora sojae were recovered from roots, stems, and rhizosphere soil of diseased soybean plants in the Jiangsu, Anhui, and Heilongjiang provinces in China from 2002 to 2022 (Table 1). Sixty-four isolates of 29 Phytophthora species, 7 isolates of five species of other oomycetes, 32 isolates of 31 fungal species, and 2 isolates of Bursaphelenchus species used in this study were obtained from a collection maintained at the Department of Plant Pathology of Nanjing Forestry University (Nanjing, China) (Table 1). Phytophthora sojae and other oomycete isolates were grown on 10% clarified V8 juice agar in a 70-mm Petri dish and maintained in the dark at 20–25°C. After culture for 3–5 days, the mycelium was scraped from the medium surface. The fungal isolates were grown on potato glucose agar in 90-mm Petri dishes and maintained at 20°C in the dark. Genomic DNA (gDNA) was extracted from all isolates using the DNAsecure Plant Kit (Tiangen Biotech, Beijing, China) and quantified using a Nanodrop ND-1000 spectrophotometer (NanoDrop Technologies, Wilmington, DE, USA). All gDNA samples were stored at-20°C until use.
Table 1 Information and Crisp-cas12a deteciton results of Phytophthora and other oomycete and fungal isolates used in this study.
Design of RPA primers, crRNA, and ssDNA reporter
The Ypt1 gene was selected as the target for the design of gene-specific RPA primers (Dai et al., 2019). The RPA primers were designed using Primer Premier 6.0 (Premier Biosoft, Palo Alto, CA, USA) and based on the recommendations in the DNA Amplification Kit manual which were listed in Supplementary Figure S1. The CHOPCHOP web tool (http://chopchop.cbu.uib.no/) was used to design the CRISPR RNA (crRNA) ssDNA reporters (Zhao et al., 2022). The crRNA sequence did not overlap with that of the RPA primers and target conserved regions of the RPA amplicon (crRNA: UAAUUUCUACUAAGUGUAGAUCGAUCCAGUUGCAGUUGCUGACAAUA) (Figure S1). The 5′-end of the ssDNA reporter was tagged with 6-FAM and the 3′-end was tagged with the BHQ-1 quencher (5′ 6-FAM-TTATT-BHQ-1 3′) (Chen et al., 2018; Li et al., 2018). The crRNA and ssDNA reporter were synthesized by GenScript (Nanjing, China) and stored at -80°C until use.
RPA-CRISPR/Cas12a assay
The 30-min assay included a two-step approach (5 min for the RPA reaction and 15 min for the CRISPR/Cas12a assay). First, the pair of RPA primers (Ypt1RPA-F/Ypt1RPA-R) were used to amplify the Ypt1 gene of P. sojae in a 5-min RPA step. Second, the CRISPR/Cas12a system was used to detect and visualize the amplified products within 15 min.
The detection process involved in the RPA-CRISPR/Cas12a analysis is shown in Figure 1. The RPA assay was performed in a 50-µL reaction mixture using the Test Strip Kit (LeShang Ltd., WuXi, China) in accordance with the manufacturer’s instructions. Each reaction mixture initially contained 2 µL of each forward and reverse primer (Ypt1RPA-F/Ypt1RPA-R, 10 µM), 25 µL rehydration buffer, 2 µL gDNA (10 ng.µL−1), and 16 µL double-distilled H2O (ddH2O), amounting to 47 µL in volume. After centrifugation of this mixture at 4000 rpm for 5 s, 3 μL of activator (supplied with the kit) was added to the lid of the reaction unit. The reaction unit was closed tightly, centrifuged at 4,000 rpm for 5 s, and then shaken manually for 3 s so that the mixture and activator were well mixed. Reactions were conducted at 37°C. After 4 min, the reactions were shaken manually and centrifuged at 4000 rpm for 5 s. The reaction tube was incubated at 37°C for 20 min. The RPA products were then analyzed using the CRISPR/Cas12a system.
Figure 1 Schematic diagram of the RPA-CRISPR/Cas12a assay for detection of Phytophthora sojae. (A) Recombinase polymerase amplification. (B) Cas12a protein can combine with each amplicon and target-specific crRNA to form a complex with indiscriminate ssDNA cleavage activity. The FAM-labeled ssDNA reporter is cleaved and produces visible green fluorescence under excitation at a wavelength of 470 nm. (C) Positive result: visible green fluorescence. Negative result: no green fluorescence visible.
The ssDNA reporter labeled with 6-FAM and BHQ-1 was added to the reaction tube of the CRISPR/Cas12a system. The Cas12a/crRNA recognized the target protospacer-adjacent motif sequence, the cleavage activity of Cas12a was generated, and the ssDNA reporter was cleaved, resulting in fluorescence. Various concentrations of crRNA (40 nM, 300 nM, 0.5 µM, 1 µM, 2 µM, 5 µM, and 10 µM) and ssDNA reporter (40 nM, 300 nM, 1 µM, 2 µM, 5 µM, and 10 µM) were used to screen the optimal concentration combination (Table S1). The optimal incubation time of the RPA reaction and the cleavage time of Cas12a were tested separately with eight time points (5, 10, 15, 20, 25, 30, 35, and 40 min). The CRISPR/Cas12a reaction was performed in a 50-μL reaction mixture, comprising 38 μL ddH2O, 5 μL reaction buffer, 3 μL crRNA, 1 μL Cas12a, 1 μL ssDNA reporter, and 2 μL RPA product. The reaction mixture was centrifuged at 4000 rpm for 5 s, then incubated at 37°C. The results of the RPA-CRISPR/Cas12a assay can be detected by two methods. Strong fluorescence signals were detected by a multifunctional microplate reader (λex 485 nm, λem 520 nm) or visible green fluorescence was detected under a blue LED transilluminator at 470 nm wavelength, whereas no fluorescent signal or visible green fluorescence was detected in the negative controls. All RPA-CRISPR/Cas12a reactions were repeated at least three times. The STDEVP function was used to analyze the three results (number 1, number 2, and number 3) obtained by repeating CRISPR/Cas12a analysis three times to calculate the standard deviation. Statistical analysis was performed using GraphPad Prism 8 software (GraphPad Software Inc., San Diego, CA, USA). The experimental group and control group were compared by performing the Student’s t-test for a difference analysis by calculating P value. P <0.05 (*) was considered statistically significan.
Specificity and sensitivity of the RPA-CRISPR/Cas12a assay
Specificity of the RPA-CRISPR/Cas12a assay was evaluated against all isolates (22 isolates of P.sojae, 64 isolates of 29 other Phytophthora species, 7 isolates of other oomycetes, 32 isolates of fungal species, and 2 isolates of Bursaphelenchus species) listed in Table 1. Purified gDNA samples (100 ng) were used as templates and a positive control (100 ng P. sojae isolate) and double-distilled H2O as a no template control (NTC) were included in each set of reactions. To determine sensitivity, 10-fold dilutions of P. sojae gDNA (isolate P6497) ranging from 100 to 0.0001 ng.µL−1 were used as DNA templates in the RPA-CRISPR/Cas12a assay. A negative control (NC; ddH2O) was included in each set of reactions. All RPA-CRISPR/Cas12a reactions were repeated at least three times. The STDEVP function was used to analyze the three results (number 1, number 2, and number 3) obtained by repeating CRISPR/Cas12a analysis three times to calculate the standard deviation. Statistical analysis was performed using GraphPad Prism 8 software (GraphPad Software Inc., San Diego, CA, USA). The experimental group and control group were compared by performing the Student’s t-test for a difference analysis by calculating P value. P <0.05 (*) was considered statistically significan.
Detection of P. sojae in naturally infected soybean plants using the RPA-CRISPR/Cas12a assay
After inoculation of P. sojae, soybean plants seedlings had severe wilting with discoloration. Using a five-point sampling method, five soybean plants were collected at the booting stage at each sampling point from a field infected with soybean root rot. For sample processing, each of the ground tissues of 25 soybean plants were cut into 2-cm-long segments, and 250 g fresh samples were taken from each plant. The samples were ground for 1 min with a grinder (XB-M101, Lingrui, China). A sample (1 g) was placed in a 1.5-mL centrifuge tube containing 1 mL PEG-OH extraction buffer (1 g NaOH was dissolved in 900 mL ddH2O, 60 mL PEG200 was added, and the volume adjusted to 1 L with ddH2O, and sterilized with a bacterial filter) and incubated at room temperature (25°C) for 8 min. The supernatant was collected as the template for the RPA-CRISPR/Cas12a assay. Purified gDNA (10 ng.µL−1) extracted from healthy soybean plants and ddH2O were included in each repeat as a NC and NTC, respectively. The STDEVP function was used to analyze the three results (number 1, number 2, and number 3) obtained by repeating CRISPR/Cas12a analysis three times to calculate the standard deviation. Statistical analysis was performed using GraphPad Prism 8 software (GraphPad Software Inc., San Diego, CA, USA). The experimental group and control group were compared by performing the Student’s t-test for a difference analysis by calculating P value. P <0.05 (*) was considered statistically significan.
Results
Optimization of the RPA-CRISPR/Cas12a assay for detection of P. sojae
Various concentrations of crRNA and ssDNA reporter for the RPA-CRISPR/Cas12a assay were used to determine the optimal concentration (Table S1). The visible green fluorescence and fluorescence intensity results indicated that 10 μM crRNA and 10 μM ssDNA reporter were the optimal concentrations for the RPA-CRISPR/Cas12a assay (Figures 2A, B). Subsequently, 10 μM crRNA and 10 μM ssDNA reporter were used to test for the optimal RPA reaction time (assessed at the time points 5, 10, 15, 20, 25, 30, 35, and 40 min). Strong green fluorescence was observed with the blue LED transilluminator after 5 min and the fluorescence intensity increased from 5 min (Figures 3A, B). Therefore, the optimal RPA reaction time was 5 min. Assessment of the optimal Cas12a cleavage time (at 5, 10, 15, 20, 25, 30, 35, and 40 min) using the product of the 5 min RPA reaction revealed that the optimal cleavage time was 15 min (Figures 3C, D). In conclusion, green fluorescence was clearly visible after 20 min (with 5 min for the RPA reaction and 15 min for Cas12a cleavage) (Figure 3).
Figure 2 Optimization of the crRNA and ssDNA reporter concentrations for the RPA-CRISPR/Cas12a assay. The crRNA and ssDNA reporter concentrations were set to the following values: 1: 40 nM crRNA, 40 nM ssDNA reporter; 2: 300 nM crRNA, 40 nM ssDNA reporter; 3: 500 nM crRNA, 300 nM ssDNA reporter; 4: 1 µM crRNA, 1 µM ssDNA reporter; 5: 1 µM crRNA, 2 µM ssDNA reporter; 6: 2 µM crRNA, 1 µM ssDNA reporter; 7: 5 µM crRNA, 1 µM ssDNA reporter; 8: 1 µM crRNA, 5 µM ssDNA reporter; 9: 10 µM crRNA, 5 µM ssDNA reporter; and 10: 10 µM crRNA, 10 µM ssDNA reporter. NC, Negative control (double-distilled H2O). (A) Fluorescence detection using a multifunctional microplate reader (λex: 485 nm, λem: 520 nm). (B) Visible green fluorescence detection under a blue LED transilluminator at a wavelength of 470 nm. By measuring the fluorescence value, one-way ANOVA was performed on the fluorescence value with fluorescence and the fluorescence value of the negative control to obtain the result of P <0.0001. This indicates that the difference between fluorescence and non fluorescence is very significant, represented by "****".
Figure 3 Optimization of the recombinase polymerase amplification (RPA) reaction time and Cas12a cleavage time for the RPA-CRISPR/Cas12a assay. (A, B) RPA reaction times: 1, 3, 5, 7, 9, 11, 13, and 15: 5, 10, 15, 20, 25, 30, 35, and 40 min; 2, 4, 6, 8, 10, 12, 14, and 16: NC (negative control, double-distilled H2O). (C, D) Cas12a cleavage times: 1, 3, 5, 7, 9, 11, 13, and 15: 5, 10, 15, 20, 25, 30, 35, and 40 min; 2, 4, 6, 8, 10, 12, 14, and 16: NC (negative control, double-distilled H2O). (A, C) Fluorescence detection using a multifunctional microplate reader (λex: 485 nm, λem: 520 nm). (B, D) Visible green fluorescence detection under a blue LED transilluminator at a wavelength of 470 nm.
Specificity of the RPA-CRISPR/Cas12a assay for rapid detection of P. sojae
A PCR amplification product of approximately 239 bp was amplified in the RPA reaction from the gDNA of P. sojae with the primers Ypt1RPA-F and Ypt1RPA-R. No PCR amplicons were detected in the reactions with gDNA of P. sansomeana, P. melonis, P. vignae, P. cinnamomi, P. cryptogea, P. citrophthora, Alternaria alternata, Botrytis cinerea, Colletotrichum truncatum, Endothia parasitica, Fusarium oxysporium, F. solani, and the NC (Figure 4). In the specificity assessment for the RPA-CRISPR/Cas12a assay, a multifunctional microplate reader detected a strong fluorescence signal for gDNA of P. sojae, whereas no fluorescence signal was observed for gDNAs of other oomycetes and nematodes (Figures 5B, D). Green fluorescence from gDNA of P. sojae was clearly observed under a blue LED transilluminator at a wavelength of 470 nm, whereas no green fluorescence was detected for the gDNAs of other oomycete isolates, fungal isolates, and nematodes (Figures 5A, C). The specificity tests were replicated three times and each replicate yielded identical results.
Figure 4 Primers designed based on the Ypt1 gene sequence were screened for specificity using conventional PCR assays, and detected by 1.5% agarose gel electrophoresis. PCR amplicons (239 bp) were only detected in the P. sojae sample, indicating specificity for detection of P. sojae DNA. Marker DL2000 (Takara Shuzo, Shiga, Japan). NC, Negative control (double-distilled H2O).
Figure 5 Evaluation of the specificity of the RPA-CRISPR/Cas12a assay. (A, B) Evaluation using genomic DNA isolated from 1: Phytophthora sojae, 2: P. melonis, 3: P. cambivora, 4: P. cinnamomi, 5: P. citricola, 6: P. cactorum, 7: P. pini, and 8: negative control (NC; no template). (C, D) Evaluation using genomic DNA from 1: P. sojae, 2: Alternaria alternata, 3: Botrytis cinerea, 4: Colletotrichum truncatum, 5: Endothia parasitica, 6: Fusarium oxysporium, 7: F. solani, and 8: negative control (NC; no template).
Sensitivity of the RPA-CRISPR/Cas12a assay for rapid detection of P. sojae
To evaluate the sensitivity of the RPA-CRISPR/Cas12a assay for detection of P. sojae, 2 µL of P. sojae gDNA of various concentrations (10 ng.µL−1, 10 ng.µL−1, 1 ng.µL−1, 100 pg.µL−1, 10 pg.µL−1, 1 pg.µL−1, 100 fg.µL−1 and 10 fg.µL−1) were used. The gDNA concentrations of 10 ng.µL−1, 1 ng.µL−1, 100 pg.µL−1, and 10 pg.µL−1 resulted in visible green fluorescence, whereas the remaining gDNA concentrations and the NC did not produce visible fluorescence. All results were consistent among three replicates. These results indicated that the minimum detectable gDNA concentration was 10 pg.μL−1 gDNA (Figure 6). The RPA-CRISPR/Cas12a assay was repeated three times for each concentration of gDNA template under identical experimental conditions as described above.
Figure 6 Evaluation of the sensitivity of the RPA-CRISPR/Cas12a assay for detection of Phytophthora sojae. The pathogen was detected at a minimum genomic DNA concentration of 10 pg µL−1 using (A) a multifunctional microplate reader (λex: 485 nm, λem: 520 nm) or (B) a blue LED transilluminator at a wavelength of 470 nm. 1–8: 10 ng µL−1, 1 ng µL−1, 100 pg µL−1, 10 pg µL−1, 1 pg µL−1, 100 fg µL−1, 10 fg µL−1 and NC (negative control, double-distilled H2O). By measuring the fluorescence value, one-way ANOVA was performed on the fluorescence value with fluorescence and the fluorescence value of the negative control to obtain the result of P <0.0001. This indicates that the difference between fluorescence and non fluorescence is very significant, represented by "****".
Detection of P. sojae in naturally infected soybean plants using the RPA-CRISPR/Cas12a assay
Based on combination of the simplified PEG-NaOH method and the RPA-
CRISPR/Cas12a assay, we detected P. sojae in naturally infected soybean samples. Using the RPA-CRISPR/Cas12a assay, green fluorescence was detected based on DNA extracted from the positive control and diseased samples, whereas no fluorescence was detected based on the DNA obtained from the NIS (non-inoculated samples) or the NC. All samples generated visible green fluorescence that was detected by a multifunctional microplate reader and yielded identical results (Figures 7A, B). These results demonstrated that the RPA-CRISPR/Cas12a assay combined with the simplified PEG-NaOH method could be used for effective detection of P. sojae in field-collected samples.
Figure 7 Detection of Phytophthora sojae in naturally infected soybean root samples using the RPA-CRISPR/Cas12a assay combined with a simple DNA extraction method. (A) Strong fluorescence signals were detected by a multifunctional microplate reader (λex: 485 nm, λem: 520 nm). 1: Positive control, PC; 2-4: naturally infected soybean root samples; 5: healthy soybean root sample; NC: negative control. (B) Visible green fluorescence was detected under a blue LED transilluminator at a wavelength of 470 nm. 1-4: Naturally infected soybean root samples; 5: healthy soybean root sample; NC: negative control. By measuring the fluorescence value, one-way ANOVA was performed on the fluorescence value with fluorescence and the fluorescence value of the negative control to obtain the result of P <0.0001. This indicates that the difference between fluorescence and non fluorescence is very significant, represented by "****".
Discussion
P. sojae is a destructive soilborne pathogen of soybean, present in most soybean-growing regions of the world, that can be difficult to diagnose and causes severe economic losses. The pathogen thrives in wet conditions and in compacted or heavy clay soils. Motile, water-borne zoospores were released from sporangia and were attracted to soybean root exudates (Morris and Ward, 1992). This pathogen can infect seeds, seedlings, and plants at all stages of growth when the soil conditions favor pathogen reproduction and infection. Symptoms are usually apparent 1–2 weeks after heavy rains and are most common where soils are poorly drained. Soybean root rot is among the most serious diseases in China. Accurate and rapid detection of P. sojae in plants and the soil is a critical step towards effective prevention and management of soybean root and crown rot and seedling damping-off.
Molecular detection technology is useful to rapidly monitor the incidence of P. sojae. P. sojae is easily mistaken for P. megasperma. Indeed, P. sojae and P. medicaginis were formerly known as P. megasperma f. sp. glycinea and P. megasperma f. sp. medicaginis, respectively. Molecular analysis of the mitochondrial region and isozyme analysis provided evidence that these were distinct species (Forster et al., 1989). Given the unique symptoms of soybean root rot in infected tissues, we aimed to establish a rapid detection technology to assist in early diagnosis of soybean root rot and the formulation of disease control strategies (Wang and Fang, 2005).
This is the first report of a RPA-CRISPR/Cas12a assay for rapid detection of P. sojae. The total amplification time was limited to 20 min (5 min for the RPA reaction and 15 min for the CRISPR/Cas12a assay) and the test results were visible under UV light and readout coming from fluorophores. The crRNA and ssDNA concentrations strongly influence the results of the RPA-CRISPR/Cas12a assay. Accordingly, different combinations of crRNA and ssDNA reporter concentrations were screened and the fluorescence intensity peaked with crRNA and ssDNA reporter concentrations of 10 μM and 10 μM, respectively. As demonstrated in the specificity evaluation, the present novel RPA-CRISPR/Cas12a system specifically detected gDNA of P. sojae, whereas no positive reactions to gDNA of 29 other Phytophthora species, including P. medicaginis, P. megasperma, P. sansomeana, P. melonis, and P. vignae (Table 1). This indicates that the method meets the requirements for specific detection of P. sojae. In addition, the lowest detectable gDNA concentration of P. sojae was 10 pg.µL−1 using the RPA-CRISPR/Cas12a assay. The feasibility of using this assay was investigated using soybean samples from a field infected with soybean root rot. The results demonstrated that the assay showed high accuracy, indicating its potential for early diagnosis of P. sojae.
This method has several advantages compared with conventional methods. First, the RPA-CRISPR/Cas12a reactions can be performed at a constant and relatively low temperature (37°C). Human body temperature, USB-powered holding boxes or thermostatic heaters can provide the required temperature conditions without the need for complex temperature control equipment, thereby eliminating a requirement for expensive specialized instruments, such as thermocyclers. Second, the entire diagnostic process can be completed within approximately 20 min. In contrast, a PCR-based detection method requires at least 2.5 h, including 90 min for the PCR reaction and 30 min for gel electrophoresis (Shin et al., 2021). Third, Some substances in the RPA reaction interfere with the antibodies on the test paper may causes non-specific binding and false positive signals to occur, which are not sufficiently diluted. Binding RPA to CRISPR-Cas12a allows us to detect the target twice, once for the recognition of the RPA primer and the second time for the recognition of the RPA amplification product by CRISPR-Cas12a when the RPA reaction is performed, which effectively avoids the problem of false positives during RPA amplification. Forth, it is not necessary to use high quality and purified DNA for amplification template which is required for conventional PCR-based diagnosis (Shin et al., 2021). That maked it an ideal choice for establishing a rapid on-site detection technology platform. A novel RPA-LFD assay was previously developed for accurate, simple, and rapid detection of P. sojae (Dai et al., 2019). The sensitivity of the RPA-CRISPR/Cas12a assay reported here is comparable, if not higher, than most previously developed detection methods. The detection limit for the RPA-CRISPR/Cas12a assay was 10 pg in a 50-μL reaction. Although the sensitivity was limited to 10 pg.μL−1, which is same with the RPA-LFD assay, the RPA-CRISPR/Cas12a assay sensitivity was sufficient to reliably detect P. sojae. The sensitivity is higher than that of a LAMP assay (Dai et al., 2012). Compared with the RPA-LFD assay, the present RPA-CRISPR/Cas12a detection system exhibits a number of promising advantages, including multiplexed detection capability, low requirement for sophisticated instruments and specialist personnel, and enhanced testing accuracy (Dai et al., 2019). Compared with the RPA-LFD method, the RPA-CRISPR/Cas12a technology enables isothermal label-free detection of target genes by designing RPA primers and performing isothermal amplification of the products obtained as amplification reporters to generate the CRISPR/Cas12a cleavage products has. Thus, the method offers the opportunity for development of an accurate, user-friendly, inexpensive platform for point-of-care testing application of CRISPR-based diagnostics.
This research aimed to develop a RPA-CRISPR/Cas12a assay for detection of P. sojae. The assay is easy to perform and portable, and does not require expensive equipment. Although these diverse properties of the CRISPR/Cas12a system provide potential for the development of versatile tools for pathogen detection (Kim et al., 2016a; Kim et al., 2016b; Li et al., 2018), there remain challenges to overcome, including few currently identified orthologs of Cas12a, limited genomic targeting coverage, and relatively low editing efficiency (Kleinstiver et al., 2016; Kim et al., 2017; Tu et al., 2017). In addition, different crRNA scaffolds affect the activities of the Cas12a–crRNA complex (Li et al., 2017). These scaffolds will help to elucidate the exact mechanisms of the reactions and enable improvements in the future. Meanwhile, the lyophilization application of CRISPR reaction needs to be further explored, and the reagents required for CRISPR reaction need to be lyophilized before use, which makes the application of CRISPR assay limited. The lyophilization application of the whole CRISPR reaction system has been reported (Nguyen et al., 2021; Rybnicky et al., 2022), and the lyophilization application of CRISPR reaction can be further explored in the future, so that the reagents can be stored for a long time at room temperature and really put into large-scale application. In addition, consecutive base changes and deletions were not investigated in this study, and there is still much room to explore whether RPA-CRISPR/Cas12a can recognize single-base differences.
The proposed CRISPR/Cas12a assay can be used for detection and identification of P. sojae in the field and laboratory. The RPA-CRISPR/Cas12a can complete the rapid detection of P. sojae at 37°C for 30 min.The accuracy of RPA-CRISPR/Cas12a was verified by collecting 30 soybean rhizosphere samples and has potentially significant applications to detect P. sojae. In-field diagnosis of soybean root rot would help to create an effective and accurate control strategy for farmers.
Conclusions
The present novel RPA-CRISPR/Cas12a-based assay exhibits high specificity and sensitivity, and a short detection time. The test results were visible under UV light and readout coming from fluorophores (Table S3) without the need for expensive equipment and facilitates the early detection of P. sojae. The assay is sensitive, efficient, and convenient. Practitioners could consider improvements to this assay to increase the sensitivity and expand detection to other pathogens.
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.
Author contributions
YG conceptualized and designed the research, analyzed the data, interpreted the results, performed the experiments and wrote the manuscript. HX participated discussed the experimental design. TD and TL revised the manuscript and directed the project. All authors contributed to the article and approved the submitted version.
Funding
This work was supported by the Jiangsu University Natural Science Research Major Project (21KJA220003), and the Priority Academic Program Development of Jiangsu Higher Education Institutions.
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/fcimb.2023.1208837/full#supplementary-material
References
Andresen, D., Von Nickisch-Rosenegk, M., Bier, F. F. (2009). Helicase-dependent amplification: use in OnChip amplification and potential for point-of-care diagnostics. Expert Rev. Mol. Diagn. 9, 645–650. doi: 10.1586/erm.09.46
Bai, J., Lin, H., Li, H., Zhou, Y., Liu, J., Zhong, G., et al. (2019). Cas12a-based on-site and rapid nucleic acid detection of African swine fever. Front. Microbiol. 10. doi: 10.3389/fmicb.2019.02830
Barrangou, R., Fremaux, C., Deveau, H., Richards, M., Boyaval, P., Moineau, S., et al. (2007). CRISPR provides acquired resistance against viruses in prokaryotes. Science. 315, 1709–1712. doi: 10.1126/science.1138140
Barrangou, R., Marraffini, L. A. (2014). CRISPR-cas systems: prokaryotes upgrade to adaptive immunity. Mol. Cell. 54, 234–244. doi: 10.1016/j.molcel.2014.03.011
Bienapfl, J. C., Malvick, D. K., Percich, J. A. (2011). Specific molecular detection of Phytophthora sojae using conventional and real-time PCR. Fungal Biol. 115, 733–740. doi: 10.1016/j.funbio.2011.05.007
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
Chen, J. S., Ma, E. B., Harrington, L. B., Costa, M. D., Tian, X. R., Palefsky, J. M., et al. (2018). CRISPR-Cas12a target binding unleashes indiscriminate single-stranded DNase activity. Science 360, 436–439. doi: 10.1126/science.aar6245
Chen, Q., Weng, Q., Wang, Y., Zheng, X. (2004). Identification and sequencing of ribosomal DNA-ITS of Phytophthora sojae fujian. Acta Phytopathol. Sin. 34, 112–116. doi: 10.13926/j.cnki.apps.2004.02.003
Clancy, E., Higgins, O., Forrest, M. S., Boo, T. W., Cormican, M., Barry, T., et al. (2015). Development of a rapid recombinase polymerase amplification assay for the detection of Streptococcus pneumoniae in whole blood. BMC Infect. Dis. 15, 481. doi: 10.1186/s12879-015-1212-5
Cui, L., Yin, W., Tang, Q., Dong, S., Zheng, X., Zhang, Z., et al. (2010). Distribution, pathotypes, and metalaxyl sensitivity of Phytophthora sojae from heilongjiang and fujian provinces in China. Plant Dis. 94, 881–884. doi: 10.1094/PDIS-94-7-0881
Cullen, D. W., Lees, A. K., Toth, I. K., Duncan, J. M. (2001). Conventional PCR and real-time quantitative PCR detection of Helminthosporium solani in soil and on potato tubers. Euro. J. Plant Pathol. 107, 387–398. doi: 10.1023/A:1011247826231
Cullen, D. W., Lees, A. K., Toth, I. K., Duncan, J. M. (2002). Detection of Colletotrichum coccodes from soil and potato tubers by conventional and quantitative real-time PCR. Plant Pathol. 51, 281–292. doi: 10.1046/j.1365-3059.2002.00690.x
Dai, T. T., Lu, C. C., Lu, J., Dong, S., Ye, W., Wang, Y., et al. (2012). Development of a loop-mediated isothermal amplification assay for detection of Phytophthora sojae. FEMS Microbiol. Lett. 334, 27–34. doi: 10.1111/j.1574-6968.2012.02619.x
Dai, B., Xiang, A., Qu, D., Chen, G., Wang, L., Wang, W., et al. (2022). Rapid and sensitive assay of Helicobacter pylori with one-tube RPA-CRISPR/Cas12 by portable array detector for visible analysis of thermostatic nucleic acid amplification. Front. Microbiol. 13. doi: 10.3389/FMICB.2022.858247
Dai, M., Feng, M., Liu, D., Cao, W., Liao, M. (2015). Development and application of SYBR Green I real-time PCR assay for the separate detection of subgroup J Avian leukosis virus and multiplex detection of avian leukosis virus subgroups A and B. Virol. J. 12, 52. doi: 10.1186/s12985-015-0291-7
Dai, T. T., Yang, X., Hu, T., Jiao, B. B., Xu, Y., Zheng, X. B., et al. (2019). Comparative evaluation of a novel recombinase polymerase amplification-lateral flow dipstick (rpa-lfd) assay, lamp, conventional pcr, and leaf-disc baiting methods for detection of Phytophthora sojae. Front. Microbiol. 10. doi: 10.3389/fmicb.2019.01884
Davison, E. M. (1998). “Book review: Phytophthora diseases worldwide,” in Plant pathology, vol. 47 . Eds. Erwin, D. C., Ribeiro, O. K., 224–225. doi: 10.1046/j.1365-3059.1998.0179a.x
Dorrance, A. E. (2018). Management of Phytophthora sojae of soybean: a review and future perspectives. Can. J. Plant Pathol. 40, 210–219. doi: 10.1080/07060661.2018.1445127
Filion, M., St-Arnaud, M., Jabaji-Hare, S. H. (2003). Direct quantification of fungal DNA from soil substrate using real-time PCR. J. Microbiol. Meth. 53, 67–76. doi: 10.1016/S0167-7012(02)00225-7
Forster, H., Kinscherf, T. G., Leong, S. (1989). F1-ATPase α subunit in the mitochondrial genome of the Phytiaceae. Nucleic Acids Res. 17, 7109. doi: 10.1093/nar/17.17.7109
Huang, C., Huang, P. T., Yao, J. Y., Li, Z. W., Weng, L. B., Guo, X. G. (2019). Pooled analysis of nuclear acid sequence-based amplification for rapid diagnosis of Mycoplasma pneumoniae infection. J. Clin. Lab. Anal. 33, e22879. doi: 10.1002/jcla.22879
Kadam, U. S., Cho, Y., Park, T. Y., Hong, J. C. (2023). Aptamer-based crispr-cas powered diagnostics of diverse biomarkers and small molecule targets. Appl. Biol. Chem. 66, 13. doi: 10.1186/s13765-023-00771-9
Kim, Y., Cheong, S. A., Lee, J. G., Lee, S. W., Lee, M. S., Baek, I. J., et al. (2016a). Generation of knockout mice by Cpf1-mediated gene targeting. Nat. Biotechnol. 34, 808–810. doi: 10.1038/nbt.3614
Kim, D., Kim, J., Hur, J. K., Been, K. W., Yoon, S. H., Kim, J. S. (2016b). Genome-wide analysis reveals specificities of Cpf1 endonucleases in human cells. Nat. Biotechnol. 34, 863–868. doi: 10.1038/nbt.3609
Kim, H. K., Song, M., Lee, J., Menon, A. V., Jung, S., Kang, Y. M., et al. (2017). In vivo high-throughput profiling of CRISPR-Cpf1 activity. Nat. Meth. 14, 153–159. doi: 10.1038/nmeth.4104
Kleinstiver, B. P., Tsai, S. Q., Prew, M. S., Nguyen, N. T., Welch, M. M., Lopez, J. M., et al. (2016). Genome-wide specificities of CRISPR-cas Cpf1 nucleases in human cells. Nat. Biotechnol. 34, 869–874. doi: 10.1038/nbt.3620
Le, C., Liu, Y., López-Orozco, J., Joyce, M. A., Le, X. C., Tyrrell, D. L. (2021). CRISPR technique incorporated with single-cell RNA sequencing for studying hepatitis b infection. Anal. Chem. 93, 10756–10761. doi: 10.1021/acs.analchem.1c02227
Lees, A. K., Cullen, D. W., Sullivan, L., Nicolson, M. J. (2002). Development of conventional and quantitative real-time PCR assays for the detection and identification of Rhizoctonia solani AG-3 in potato and soil. Plant Pathol. 51, 293–302. doi: 10.1046/j.1365-3059.2002.00712.x
Li, S. Y., Cheng, Q. X., Wang, J. M., Li, X. Y., Zhang, Z. L., Gao, S., et al. (2018). CRISPR-Cas12a-assisted nucleic acid detection. Cell Discovery 4, 20. doi: 10.1038/s41421-018-0028-z
Li, F., Xiao, J., Yang, H., Yao, Y., Li, J., Zheng, H., et al. (2022). Development of a rapid and efficient RPA-CRISPR/Cas12a assay for Mycoplasma pneumoniae detection. Front. Microbiol. 13. doi: 10.3389/FMICB.2022.858806
Li, B., Zhao, W., Luo, X., Zhang, X., Li, C., Zeng, C., et al. (2017). Engineering CRISPR–Cpf1 crRNAs and mRNAs to maximize genome editing efficiency. Nat. BioMed. Eng. 1, 66. doi: 10.1038/s41551-017-0066
Liu, L., Zhao, G., Li, X., Xu, Z., Lei, H., Shen, X. (2022). Development of rapid and easy detection of Salmonella in food matrics using RPA-CRISPR/Cas12a method. LWT 162, 113443. doi: 10.1016/J.LWT.2022.113443
Magriñá Lobato, I., O’Sullivan, C. K. (2018). Recombinase polymerase amplification: basics, applications and recent advances. Trends Anal. Chem. 98, 19–35. doi: 10.1016/j.trac.2017.10.015
Morris, P. F., Ward, E. W. B. (1992). Chemoattraction of zoospores of the soybean pathogen, Phytophthora sojae, by isoflavones. Physiol. Mol. Plant Pathol. 40, 17–22. doi: 10.1016/0885-5765(92)90067-6
Myhrvold, C., Freije, C. A., Gootenberg, J. S., Abudayyeh, O. O., Metsky, H. C., Durbin, A. F., et al. (2018). Field-deployable viral diagnostics using CRISPR-Cas13. Science. 360, 444–448. doi: 10.1126/science.aas8836
Nguyen, P. Q., Soenksen, L. R., Donghia, N. M., Angenent-Mari, N. M., Puig, H. D., Huang, A., et al. (2021). Wearable materials with embedded synthetic biology sensors for biomolecule detection. Nat. Biotechnol. 39, 1366–1374. doi: 10.1186/s13765-023-00771-9
Polstra, A., Goudsmit, J., Cornelissen, M. (2002). Development of real-time NASBA assays with molecular beacon detection to quantify mRNA coding for HHV-8 lytic and latent genes. BMC Infect. Dis. 2, 18. doi: 10.1186/1471-2334-2-18
Roese, A. D., Pereira Goulart, A. C. (2013). First report of root and stem rot of soybean caused by Phytophthora sojae in mato grosso do sul state, Brazil. Summa Phytopathol. 39, 293. doi: 10.1590/S0100-54052013000400012
Rojas, J. A., Miles, T. D., Coffey, M. D., Martin, F. N., Chilvers, M. I. (2017). Development and application of qPCR and RPA genus- and species-specific detection of Phytophthora sojae and P. sansomeana root rot pathogens of soybean. Plant Dis. 101, 1171–1181. doi: 10.1094/PDIS-09-16-1225-RE
Rybnicky, G. A., Dixon, R. A., Kuhn, R. M., Karim, A. S., Jewett, M. C. (2022). Development of a freeze-dried crispr-cas12 sensor for detecting Wolbachia in the secondary science classroom. ACS SYNTH Biol. 11, 835–842. doi: 10.1021/acssynbio.1c00503
Schmitthenner, A. F. (1985). Problems and progress in control of phytophthora root rot of soybean. Plant Dis. 69, 362–368. doi: 10.1094/PD-69-362
Shi, J. (2021). CIRSPR: a unfulfilled journey from gene editing in ‘Blind box’ to ‘Precision targeting’ genome editing. J. For. Res. Univ. (Nat. Sci. Ed.) 45, 12–14. doi: 10.12302/j.issn.1000-2006.202111027
Shin, K., Kwon, S. H., Lee, S. C., Moon, Y. E. (2021). Sensitive and rapid detection of citrus scab using an RPA-CRISPR/Cas12a system combined with a lateral flow assay. Plants. 10, 2132. doi: 10.3390/PLANTS10102132
Su, Y., Shen, C. (1993). The discovery and biological characteristic studies of Phytophthora megasperma f. sp. glycinea on soyabean in China. Acta Phytopathologica Sinica. 23, 341–347. doi: 10.13926/j.cnki.apps.1993.04.016
Tu, M., Lin, L., Cheng, Y., He, X. B., Sun, H. H., Xie, H. H., et al. (2017). A ‘new lease of life’: FnCpf1 possesses DNA cleavage activity for genome editing in human cells. Nucleic Acids Res. 45, 11295–11304. doi: 10.1093/nar/gkx783
Tyler, B. M. (2007). Phytophthora sojae: root rot pathogen of soybean and model oomycete. Mol. Plant Pathol. 8, 1–8. doi: 10.1111/j.1364-3703.2006.00373.x
Wang, W. Y., Fang, T. J. (2005). Development of multiplex and quantitative PCR assay to detect genetically modified roundup ready soybean in foods. J. Food Drug Anal. 13, 132–138. doi: 10.38212/2224-6614.2537
Wang, Y., Zhang, Y., Chen, J., Wang, M., Zhang, T., Luo, W., et al. (2021). Detection of SARS-CoV-2 and its mutated variants via CRISPR-Cas13-Based transcription amplification. Analytical Chem. 93, 3393–3402. doi: 10.1021/acs.analchem.0c04303
Wang, Y. C., Zhang, W. L., Zheng, X. B. (2006). Rapid and sensitive detection of Phytophthora sojae in soil and infected soybeans by species-specific polymerase chain reaction assays. Phytopathology. 96, 1315–1321. doi: 10.1094/PHYTO-96-1315
Wrather, J. A., Koenning, S. R. (2006). Estimates of disease effects on soybean yields in the United States 2003 to 2005. J. Nematology. 38, 173–180.
Xiong, Q., Xu, J., Zheng, X. Y., Zhu, Y., Zhang, C., Wang, X. L., et al. (2019). Development of seven novel specific SCAR markers for rapid identification of Phytophthora sojae: the cause of root- and stem-rot disease of soybean. Eur. J. Plant Pathol. 153, 517–531. doi: 10.1007/s10658-018-1579-4
Xu, Z. C., Chen, D. J., Li, T., Yan, J. Y., Zhu, J., He, T., et al. (2022). Microfluidic space coding for multiplexed nucleic acid detection via CRISPR-Cas12a and recombinase polymerase amplification. Nat. Commun. 13, 6480. doi: 10.1038/s41467-022-34086-y
Xu, J. H., Ma, J. F., Li, Y. W., Kang, L., Yuan, B., Li, S. Q., et al. (2022). A general RPA-CRISPR/Cas12a sensing platform for Brucella spp. detection in blood and milk samples. Sensors Actuat: B-Chem. 364, 131864. doi: 10.1016/J.SNB.2022.131864
Yue, H., Huang, M., Tian, T., Xiong, E., Zhou, X. (2021). Advances in clustered, regularly interspaced short palindromic repeats (CRISPR)-based diagnostic assays assisted by Micro/Nanotechnologies. ACS Nano. 15, 7848–7859. doi: 10.1021/acsnano.1c02372
Zhang, S. Z., Xu, P. F., Wu, J. J., Xue, A. G., Zhang, J. X., Li, W. B., et al. (2010). Races of Phytophthora sojae and their virulences on soybean cultivars in heilongjiang, China. Plant Dis. 94, 87–91. doi: 10.1094/PDIS-94-1-0087
Zhang, Y. M., Zhang, Y., Xie, K. B. (2020). Evaluation of CRISPR/Cas 12a-based DNA detection for fast pathogen diagnosis and GMO test in rice. Mol. Breeding. 40, 495–4504. doi: 10.1007/s11032-019-1092-2
Zhao, G., Wang, J., Yao, C. Y., Xie, P. C., Li, X. M., Xu, Z. L., et al. (2022). Alkaline lysis-recombinase polymerase amplification combined with CRISPR/Cas12a assay for the ultrafast visual identification of pork in meat products. Food Chem. 383, 132318. doi: 10.1016/j.foodchem.2022.132318
Keywords: Phytophthora sojae, CRISPR/Cas12a, recombinase polymerase amplification, disease diagnosis, detection
Citation: Guo Y, Xia H, Dai T and Liu T (2023) RPA-CRISPR/Cas12a mediated isothermal amplification for visual detection of Phytophthora sojae. Front. Cell. Infect. Microbiol. 13:1208837. doi: 10.3389/fcimb.2023.1208837
Received: 19 April 2023; Accepted: 03 May 2023;
Published: 26 May 2023.
Edited by:
Hokyoung Son, Seoul National University, Republic of KoreaReviewed by:
Ulhas S. Kadam, Gyeongsang National University, Republic of KoreaJing Si, Beijing Forestry University, China
Copyright © 2023 Guo, Xia, Dai and Liu. 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: Tingting Dai, 13770647123@163.com; Tingli Liu, tlliu2100@163.com