- 1College of Veterinary Medicine, Guangdong Laboratory for Lingnan Modern Agriculture, South China Agricultural University, Guangzhou, China
- 2Key Laboratory of Animal Vaccine Development, Ministry of Agriculture and Rural Affairs, Guangzhou, China
- 3Guangdong Provincial Key Laboratory of Zoonosis Prevention and Control, Guangzhou, China
- 4Guangzhou Collaborative Innovation Center on Science-Tech of Ecology and Landscape, Guangzhou Zoo, Guangzhou, China
Duck Tembusu virus (DTMUV) is an emerging pathogen that poses a serious threat to the duck industry in China. Currently, polymerase chain reaction (PCR), quantitative PCR (qPCR) and reverse transcription loop-mediated isothermal amplification (RT-LAMP) are commonly used for DTMUV detection. However, these methods require complex steps and special equipment and easily cause false-positive results. Therefore, we urgently need to establish a simple, sensitive and specific method for the clinical field detection of DTMUV. In this study, we developed an RT-LAMP-based CRISPR-Cas12a assay targeting the C gene to detect DTMUV with a limited detection of 3 copies/μL. This assay was specific for DTMUV without cross-reaction with other common avian viruses and only required some simple pieces of equipment, such as a thermostat water bath and blue/UV light transilluminator. Furthermore, this assay showed 100% positive predictive agreement (PPA) and negative predictive agreement (NPA) relative to SYBR Green qPCR for DTMUV detection in 32 cloacal swabs and 22 tissue samples, supporting its application for clinical field detection.
1 Introduction
Tembusu virus (TMUV) is a single-stranded, positive-sense RNA virus, that belongs to the genus Flavivirus in the family Flaviviridae (Zhang et al., 2017). The genome of TMUV is formed by a 5′-terminal noncoding region (NCR), a single open reading frame (ORF) encoding polyprotein and a 3′-terminal noncoding region (NCR). The polyprotein is cleaved into three structural proteins (capsid C, membrane precursor prM, and envelope E) and seven nonstructural proteins (NS1, NS2A, NS2B, NS3, NS4A, NS4B, and NS5) by viral and host proteases (Hamel et al., 2021; Xu et al., 2023). The virus was first isolated from Culex tritaeniorhynchus mosquitoes in Malaysia in 1955 and then infected a wide range of birds, including ducks, chickens, geese, pigeons and sparrows, due to its cross-species transmission (Yu et al., 2018). Among these hosts, ducks are the most susceptible animals to TMUV infection, and duck-origin TMUV is abbreviated as DTMUV (Tang et al., 2012; Zhu et al., 2012; Thontiravong et al., 2015; Yan et al., 2017; Tunterak et al., 2018). Domestic ducks infected with DTMUV exhibit diarrhea, anorexia, encephalitis, paralysis, and a significant decline in egg production, resulting in a mortality rate ranging from 5% to 30% (Yan P. X. et al., 2011; Liang et al., 2016; Li et al., 2020). Since the first report of DTMUV in China in 2010, it has infected more than 130 million domestic ducks and presents a grave danger to the duck industry (Lei et al., 2017). Herein, the diagnosis and prevention of DTMUV infection are essential for the duck industry. In addition to clinical diagnosis methods, the molecular diagnosis methods, including PCR-based assays, qPCR-based assays and isothermal amplification-based assays are commonly used for DTMUV infection. Currently, several PCR-based assays, targeting the E, NS3 and NS5 genes have been established for DTMUV detection (Zang et al., 2012; Li, 2013; Wang et al., 2017). In addition, some qPCR-based assays targeting the E, NS1 and 3′-NCR genes were also exploited to detect DTMUV (Yan et al., 2012; Yun et al., 2012; Chen G. Q. et al., 2018; Li et al., 2023). However, these PCR-based and qPCR-based assays are mainly applied to laboratory testing of DTMUV. To develop clinical field detection methods, several RT-LAMP assays targeting the E or NS5 gene were established for DTMUV detection (Jiang et al., 2012; Wu et al., 2017; Wang et al., 2022). The assays effectively amplify target genes with four to six primers within 1 h, and their products are visualized by the naked eye via turbidity or dye indicators (Yan et al., 2020; Hongjaisee et al., 2021; Soroka et al., 2021). However, these assays easily cause false-positive results owing to their nonspecific amplification (Soroka et al., 2021). Therefore, we urgently need to establish a sensitive and specific molecular diagnostic method for the clinical field detection of DTMUV.
Recently, the Cas12a (Cpf1), a member of the clustered regularly interspaced short palindromic repeats (CRISPR)-associated system (Cas), was used for the molecular diagnosis of pathogenic microorganisms (Chertow, 2018; Li et al., 2018b; Liang et al., 2019). The Cas12a recognizes the target DNA under the guidance of CRISPR-derived RNA (crRNA) complementary to the target DNA sequence with T nucleotide-rich PAM (5′-TTTV-3′) and some suboptimal PAMs (5′-TCTV-3′, 5′-TTV-3′, 5′-CCCV-3′ and 5′-TCCV-3′) (Dronina et al., 2022). Subsequently, the Cas12a performs the collateral cleavage on nonspecific single-stranded DNA (ssDNA) (e.g., the fluorescent ssDNA reporter HEX-N12-BHQ1) when the Cas12a, crRNA, and target DNA form a ternary complex (Gootenberg et al., 2018; Li et al., 2018b). Based on this collateral cleavage activity, a number of Cas12a-based molecular diagnostic methods coupled with isothermal amplification techniques containing RPA, RCA, RAA and LAMP have been developed to detect pathogens, including RNA viruses, DNA viruses and bacteria (Li et al., 2018b; Xu et al., 2020; Mayuramart et al., 2021; Chen et al., 2022; Qin et al., 2022; de Sousa et al., 2023; Fan et al., 2023; Mao et al., 2023; Wu et al., 2023; Zhao et al., 2023). For instance, an RT-RAA-based CRISPR-Cas12a assay was established to detect SARS-CoV-2 and influenza viruses, and the assay was effective for screening of these viruses (Mayuramart et al., 2021). Another RT-RAA-based CRISPR-Cas12a assay was also developed to detect SARS-CoV-2 on automated centrifugal microfluidics, and this assay achieved ultrasensitive SARS-CoV-2 detection using two crRNAs to simultaneously match two sites of the target DNA (Chen et al., 2022). In addition, a LAMP-based CRISPR-Cas12a assay was exploited to detect African swine fever virus (ASFV), and this assay performed high sensitivity and specificity for the diagnosis of ASFV (Qin et al., 2022). Furthermore, an RCA-based CRISPR-Cas12a assay was exploited for methicillin-resistant Staphylococcus aureus detection, and this assay was applied to point-of-care testing (Xu et al., 2020). Altogether, these assay-based CRISPR-Cas12a methods have higher specificity and sensitivity than conventional PCR-based assays, qPCR-based assays and isothermal amplification-based assays, and are suitable for clinical field detection.
In this study, we established an RT-LAMP-based CRISPR-Cas12a assay to detect DTMUV infection. In addition, we evaluated its performance in terms of specificity and sensitivity, and validated its clinical field detection with cloacal swabs and tissue samples taken from DTMUV-infected ducks.
2 Materials and methods
2.1 Viruses and clinical samples
Duck Tembusu virus (DTMUV, S132), H5N6 highly pathogenic avian influenza virus (H5N6 AIV, S230), H7N9 highly pathogenic avian influenza virus (H7N9 AIV, Q26), H9N2 avian influenza virus (H9N2 AIV, S448), infectious bronchitis virus (IBV, M41), and Newcastle disease virus (NDV, S391) were isolated by our laboratory. The 54 clinical samples, including 32 cloacal swabs and 22 tissue samples, were randomly taken from the ducklings infected with DTMUV in the animal infection experiment.
2.2 RNA extraction and cDNA synthesis
Total RNA was extracted from these viruses and clinical samples using a Fast Total RNA Extraction Kit (Fastagen). Then, these total RNAs were utilized to synthesize cDNAs with PrimeScript™ II 1st Strand cDNA Synthesis Kit (Takara). Finally, the obtained RNAs and cDNAs were preserved at a temperature of −80°C for further experimental investigation.
2.3 SYBR Green qPCR assay
The SYBR Green qPCR assay previously established by our laboratory was used for DTMUV detection in the 54 clinical samples mentioned above. Each 20 μL SYBR Green qPCR consisted of 10 μL 2× GoTaq® qPCR Master Mix (Promega), 10 μM NS5-F (forward primer), 10 μM NS5-R (reverse primer) and 5 μL cDNA (Table 1). Then, the reaction was conducted on a Bio-Rad CFX96 system (Bio-Rad) using the following thermal profile: 95°C for 2 min and 40 cycles of 95°C for 15 s and 58°C for 30 s. A sample with an average Ct value ≤35 was positive, and a sample with an average Ct value >35 or undetectable DTMUV RNA level was negative.
In addition, the limited detection of this SYBR Green qPCR assay was determined. Briefly, the copy numbers of standard plasmid DNA containing the C gene of DTMUV previously constructed by our laboratory were calculated via the NEBioCalculator with the formula: DNA copy numbers = (M × 6.022 × 1023)/(N × 109 × 660), in which M represented the amount of DNA in nanograms, N was the number of base pairs in a single standard plasmid DNA. Then, this DNA was diluted in a tenfold series (105 to 100 copies/μL) based on its copy numbers. Finally, 1 μL of tenfold serial dilution was subjected to the SYBR Green qPCR assay with three repeats.
2.4 Establishment and optimization of the RT-LAMP assay
Thirty whole genome sequences of DTMUV isolated in China were downloaded from GenBank and aligned by MegAlign (LaserGene 7.0 software, DNAStar) to find conserved regions. Based on the conserved region, the RT-LAMP primers, including FIP, BIP, F3, B3, Loop-F and Loop-B, were designed by PrimerExplorer V5 and synthesized by GENEWIZ. According to the recommendation of the NEB company, we initially established an RT-LAMP reaction containing 2.5 μL 10× isothermal amplification buffer (NEB), 1.6 μM each of FIP and BIP primers, 0.2 μM each of F3 and B3 primers, 0.4 μM each of Loop-F and Loop-B primers, 6 mM MgSO4 (NEB), 1.4 mM dNTP (Promega), 7.5 units WarmStart® RTx Reverse Transcriptase (NEB), 8 units Bst 2.0 WarmStart® DNA Polymerase (NEB) and 2 μL RNA. The total volume of this reaction was adjusted to 25 μL with sterile deionized water. Then, the reaction was performed at 65°C for 50 min.
The RT-LAMP assay was optimized via a variable-controlling approach with different concentrations of loop primers, MgSO4 and dNTPs, different ratios of outer primers (F3 and B3 primers) to inner primers (FIP and BIP primers), different reaction times, and different reaction temperatures. Two microlitre of RT-LAMP product was utilized for 1% agarose gel electrophoresis, and the grayscale value of each agarose gel electrophoretic band was measured by ImageJ software (National Institutes of Health).
2.5 Evaluation of specificity and sensitivity for the RT-LAMP assay
To evaluate the specificity of this RT-LAMP assay, H5N6 AIV, H7N9 AIV, H9N2 AIV, IBV, and NDV were utilized to determine whether these viruses were able to be falsely detected. Briefly, 1 μL RNA of each virus was subjected to this assay with three repeats, and 2 μL of RT-LAMP product was utilized for 1% agarose gel electrophoresis.
To evaluate the sensitivity of this assay, the copy numbers of DTMUV RNA were firstly calculated by NEBioCalculator with the formula: RNA copy numbers = (M × 6.022 × 1023)/(N × 109 × 340), in which M represented the amount of RNA in nanograms, N was the length of the RNA in base. Then, the DTMUV RNA was diluted to 109, 108, 107, 106, 105, 104, 103, 102, 10, 9, 8, 7, 6, 5, 4, 3 and 2 copies/μL, and 1 copy/μL based on its copy numbers. Finally, 1 μL of tenfold serial dilution was subjected to this RT-LAMP assay with three repeats, and 2 μL of RT-LAMP product was utilized for 1% agarose gel electrophoresis.
2.6 Preparation of Cas12a crRNAs
CRISPR-derived RNAs (crRNAs) complementary to the target sequences located on the C gene of DTMUV with T nucleotide-rich PAM (5′-TTTV-3′) and suboptimal PAM (5′-TTV-3′) were designed by CRISPR-DT software. The crRNAs were chemically synthesized by GenScript (Nanjing, China).
2.7 CRISPR-Cas12a-mediated cleavage assay
The 20 μL of CRISPR-Cas12a-mediated cleavage reaction contained 100 nM Cas12a (NEB), 125 nM crRNA, 500 nM HEX-N12-BHQ1 (HuicH), 40 units RNasin® Ribonuclease Inhibitors (Promega), 2 μL NEBuffer™ r2.1 (NEB), and 2 μL RT-LAMP product. The reaction was performed at 37°C for 30 min, and the result was obtained via a blue/UV light transilluminator (Thermo Fisher). In addition, the product of this assay was transferred into a 96-well plate containing 80 μL sterile-deionized water and then monitored for fluorescence intensity by a microplate reader (Promega) at wavelengths from 533 to 559 nm.
2.8 Optimization of the cleavage time of CRISPR-Cas12a
To optimize the cleavage time of CRISPR-Cas12a, the cleavage reactions were performed at 37°C, and fluorescence results were read at 5, 10, 15, 20, 25 and 30 min, respectively. The reactions corresponding to each cleavage time were conducted with three repeats. The results were obtained using a blue/UV light transilluminator. In addition, the products of these reactions were transferred into a 96-well plate containing 80 μL sterile-deionized water and then monitored for fluorescence intensity by a microplate reader at wavelengths from 533 to 559 nm.
2.9 Evaluation of specificity and sensitivity for the RT-LAMP-based CRISPR-Cas12a assay
To evaluate the specificity of the RT-LAMP-based CRISPR-Cas12a assay, H5N6 AIV, H7N9 AIV, H9N2 AIV, IBV and NDV were used to determine whether these viruses could be falsely detected. Briefly, 1 μL RNA of each virus was subjected to this assay with three repeats. The results were obtained using a blue/UV light transilluminator.
To evaluate the sensitivity of this RT-LAMP-based CRISPR-Cas12a assay, the copy number of DTMUV RNA was calculated by the same formula mentioned above. Then, the DTMUV RNA was diluted to 105, 104, 103, 102, 10, 9, 8, 7, 6, 5, 4, 3 and 2 copies/μL, and 1 copy/μL based on its copy numbers. Finally, 1 μL of tenfold serial dilution was subjected to this assay with three repeats. The results were obtained using the blue/UV light transilluminator.
2.10 RT-LAMP-based CRISPR-Cas12a assay for clinical sample detection
The 54 clinical samples taken from DTMUV-infected ducklings, including 22 cloacal swabs and 32 tissue samples, were used to detect DTMUV by RT-LAMP-based CRISPR-Cas12a assay. This assay was run in parallel with the SYBR Green qPCR assay for these clinical samples. The results were obtained with the blue/UV light transilluminator. Then, the results were compared with those obtained with the SYBR Green qPCR assay as mentioned above.
2.11 Statistical analysis
All data were analyzed with GraphPad Prism 5.0 (La Jolla, CA, United States). Significant differences were evaluated using a two-way analysis of variance (ANOVA). p-values <0.05 were considered statistically significant. Statistical significance is indicated as follows: *p < 0.05, **p < 0.01, and ***p < 0.001, and no significant difference is represented as “ns.”
2.12 Ethics statement
All animal experiments were approved via the Laboratory Animal Ethics Committee of South China Agricultural University. Highly pathogenic avian influenza viruses were conducted in the BSL-3 facility.
3 Results
3.1 Workflow of the RT-LAMP-based CRISPR-Cas12a assay
In this study, cloacal swabs or tissue samples were used for rapid extraction of RNA within 20 min. Then, the target DNA was amplified at 63°C for 60 min by the RT-LAMP assay using the extracted RNA and cleaved at 37°C for 20 min by CRISPR Cas12a under the guidance of specific crRNA. Under this condition, the collateral cleavage activity of Cas12a for ssDNA without selectivity was activated, and then the fluorescent ssDNA reporter HEX-N12-BHQ1 was cleaved by it. The resulting fluorescence was detected by a simple blue light transilluminator, UV light transilluminator or microplate reader (Figure 1).
3.2 Establishment of the RT-LAMP assay
Thirty whole genome sequences of DTMUV isolated in China were downloaded from NCBI GenBank and aligned to find their conservative regions. As a result, the nucleotide sequences of the C genes for these DTMUV were relatively conserved (Supplementary Figure S1B). Based on their conserved regions, a set of specific RT-LAMP primers, including F3, B3, FIP, BIP, Loop-F and Loop-B, were designed (Supplementary Figures S1A,B; Table 1). The RT-LAMP reaction system recommended by the NEB company was used as the standard for optimizing each component in this system. The results showed that the optimal concentrations of loop primers, dNTPs and MgSO4 were 0.6 μM, 1.4 mM and 8 mM, respectively (Figures 2B–D). In addition, the optimal ratio of outer primers to inner primers was 1:8 in this RT-LAMP reaction system (Figure 2A). Meanwhile, the optimal reaction time and reaction temperature for this system are 60 min and 63°C, respectively (Figures 2E,F).
Figure 2. Optimization of the RT-LAMP assay. (A) Optimization of the ratio of outer primers to inner primers. Lanes 1–5 represent 1:1, 1:2, 1:4, 1:8, and 1:16, respectively. (B) Optimization of the concentration of loop primers. Lanes 1–5 represent 0.1, 0.2, 0.4, 0.6, and 0.8 μM, respectively. (C) Optimization of the concentration of dNTPs. Lanes 1–5 represent 0.8, 1, 1.2, 1.4, and 1.6 mM, respectively. (D) Optimization of the concentration of MgSO4. Lanes 1–5 represent 2, 4, 6, 8, and 10 mM, respectively. (E) Optimization of reaction time. Lanes 1–5 represent 30, 40, 50, 60, and 70 min, respectively. (F) Optimization of reaction temperature. Lanes 1–5 represent 59, 61, 63, 65, and 67°C, respectively. In each figure, Lane 6 represents the negative control and the red asterisk indicates the optimal conditions for this assay. Each reaction was performed with three repeats, and one representative readout is shown in the figures. The grayscale value of each agarose gel electrophoretic band was measured by ImageJ software (National Institutes of Health).
3.3 Specificity and sensitivity of the RT-LAMP assay
To evaluate the specificity of the RT-LAMP assay, the RNAs of H5N6 AIV, H7N9 AIV, H9N2 AIV, IBV, and NDV were used for this assay with specific primers of DTMUV. The results showed that the RT-LAMP reaction with DTMUV RNA generated the expected product, and no product was observed for the RT-LAMP reaction with other viral RNAs, which indicated that this RT-LAMP assay was specific for DTMUV (Figure 3A). To evaluate the sensitivity of this assay, the DTMUV RNA was serially diluted from 109 copies/μL to 1 copy/μL. We found that its limit of detection was 7 copies/μL (Figure 3B).
Figure 3. Specificity and sensitivity of the RT-LAMP assay. (A) Specificity of the RT-LAMP assay. The plus sign (+) and minus sign (−) represent positive samples and negative controls, respectively. Each reaction was performed with three repeats, and one representative readout is shown in the figure. (B) Sensitivity of the RT-LAMP assay. Serial dilutions of DTMUV RNA from 109 copies/μL to 1 copy/μL were used in this assay. Each dilution was performed with three repeats, and one representative readout is shown in the figure.
3.4 CRISPR-Cas12a-mediated cleavage assay
In this study, three crRNAs, crRNA-1, crRNA-2 and crRNA-3, located on the target DNA were designed and tested for the Cas12a-mediated cleavage assay (Table 1). Among them, crRNA-1 performed better than crRNA-2 and crRNA-3 because of its lower fluorescence background for the negative control in this assay (Figures 4A–D). crRNA-3 was not specific enough for the target DNA due to its high fluorescence background for the negative control (Figures 4C,D).
Figure 4. CRISPR-Cas12a-mediated cleavage assay. (A–C) crRNA-1, crRNA-2 and crRNA-3 targeting the C gene of DTMUV were used for the CRISPR-Cas12a-mediated cleavage assay. Each crRNA was performed with three repeats in this assay. (D) The fluorescence intensity value of each reaction was measured by a microplate reader. All data were analyzed with GraphPad Prism 5.0 (La Jolla, CA, United States). Significant differences were evaluated using a two-way analysis of variance (ANOVA). A p-value <0.05 was considered statistically significant. Statistical significance is indicated in the figure as follows: **p < 0.01 and ***p < 0.001.
3.5 Optimization of the cleavage time of the CRISPR-Cas12a-mediated assay
To find the optimal cleavage time of this assay, the fluorescence intensity values corresponding to different time points were measured. We found that the average fluorescence intensity values at 5, 10, 15 and 20 min were 104.57, 104.72, 104.88 and 104.98, respectively, which indicated that the cleavage efficiency of CRISPR-Cas12a significantly increased with the extension of cleavage time within 20 min (p < 0.01) (Figures 5A–D,G). However, the average fluorescence intensities at 20, 25 and 30 min showed no significant change (Figures 5D–G). In summary, the best cleavage time of this assay was 20 min.
Figure 5. Optimization for cleavage time of the CRISPR-Cas12a-mediated assay. (A–F) Cleavage results of the CRISPR-Cas12a-mediated assay at different times. The reaction was performed with three repeats at each cleavage time. (G) The fluorescence intensity value of each reaction was measured by a microplate reader. All data were analyzed with GraphPad Prism 5.0 (La Jolla, CA, United States). Significant differences were evaluated using a two-way analysis of variance (ANOVA). A p-value <0.05 was considered to indicate a significant difference. Statistical significance is indicated in the figure as follows: **p < 0.01, ***p < 0.001, and no significant difference is represented as “ns.”
3.6 Specificity and sensitivity of the RT-LAMP-based CRISPR-Cas12a assay
To evaluate the specificity of the RT-LAMP-based CRISPR-Cas12a assay, the RNAs of H5N6 AIV, H7N9 AIV, H9N2 AIV, IBV and NDV were used for this assay with specific primers and crRNA of DTMUV. We found that only the CRISPR-Cas12a reaction with the RT-LAMP product corresponding to DTMUV RNA emitted fluorescence, which indicated that this RT-LAMP-based CRISPR-Cas12a assay was specific for DTMUV (Figure 6A). To evaluate the sensitivity of this assay, the DTMUV RNA was serially diluted from 105 copies/μL to 1 copy/μL. We found that its limited detection reached 3 copies/μL for DTMUV RNA (Figure 6B). Furthermore, this assay exhibited higher sensitivity compared to SYBR Green qPCR, whose limited detection was 100 copies/μL for standard plasmid DNA of DTMUV (Figure 6C).
Figure 6. Specificity and sensitivity of the RT-LAMP-based CRISPR-Cas12a assay. (A) Specificity of the RT-LAMP-based CRISPR-Cas12a assay. The plus sign (+) and minus sign (−) represent positive samples and negative controls, respectively. Each reaction was performed with three repeats, and one representative readout is shown in the figure. (B) Sensitivity of the RT-LAMP-based CRISPR-Cas12a assay. Serial dilutions of DTMUV RNA from 105 copies/μL to 1 copy/μL were used in this assay. Each reaction was performed with three repeats, and one representative readout is shown in the figure. (C) Sensitivity of SYBR Green qPCR. Tenfold serial dilutions of DTMUV RNA from 105 copies/μL to 1 copy/μL were used in this assay. Each reaction was performed with three repeats. A sample with an average Ct value ≤35 was positive, and a sample with an average Ct value >35 or undetectable DTMUV RNA level was negative.
3.7 Evaluation of the RT-LAMP-based CRISPR-Cas12a assay on clinical samples
To further evaluate the RT-LAMP-based CRISPR-Cas12a assay, we tested 54 clinical samples (32 cloacal swabs and 22 tissue samples) taken from DTMUV-infected ducks. All clinical samples were first tested using SYBR Green qPCR previously established by our laboratory. Of these samples, 23 of 32 cloacal swabs (Sample numbers 2–4, 6–7, 9–10, 13–16, 18–26, 28–29, and 32) were positive for DTMUV infection by SYBR Green qPCR, but the other 9 cloacal swabs (Sample numbers 1, 5, 8, 11–12, 17, 27, 30, and 31) were negative (Figure 7A). Then, these samples were used to evaluate the RT–LAMP-based CRISPR-Cas12a assay. Our results showed that 23 of 32 cloacal swabs (Sample numbers 2–4, 6–7, 9–10, 13–16, 18–26, 28–29, and 32) were positive for DTMUV infection by the RT-LAMP-based CRISPR-Cas12a assay, but the remaining 9 cloacal swabs were negative (Figure 7B). Thus, the PPA and NPA of the RT-LAMP-based CRISPR-Cas12a assay compared to the SYBR Green qPCR were all 100% for DTMUV detection in 32 total cloacal swabs. Thereafter, we also tested 22 tissue samples by the two methods. As a result, 17 of 22 tissue samples (Sample numbers 1–3, 6, and 10–22) were positive for DTMUV infection by the RT-LAMP-based CRISPR-Cas12a assay and SYBR Green qPCR, but 5 of 22 tissue samples (Sample numbers 4–5 and 7–9) were negative, which indicated that there were 100%PPA and 100%NPA for detection of the 22 tissue samples using these two methods (Figures 8A,B). In conclusion, there was 100% concordance between the RT-LAMP-based CRISPR-Cas12a assay and SYBR Green qPCR for DTMUV detection in cloacal swabs and tissue samples.
Figure 7. Evaluation of the RT-LAMP-based CRISPR-Cas12a assay with cloacal swabs. (A) SYBR Green qPCR assay on cloacal swabs. Numbers 1–32 represent cloacal swab identification numbers. Numbers 33 and 34 represent the positive sample and negative control, respectively. A sample with an average Ct value ≤35 was positive, and a sample with an average Ct value >35 or undetectable DTMUV RNA level was negative. (B) RT-LAMP-based CRISPR-Cas12a assay on cloacal swabs. Lanes 1–32 represent cloacal swabs. Lanes 33 and 34 represent the positive sample and negative control, respectively. In each figure, the red asterisk indicates the negative samples identified by the two assays.
Figure 8. Evaluation of the RT-LAMP-based CRISPR-Cas12a assay on tissue samples. (A) SYBR Green qPCR assay on tissue samples. Numbers 1–22 represent tissue sample identification numbers. Numbers 23 and 24 represent positive samples and negative controls, respectively. A sample with an average Ct value ≤35 was positive, and a sample with an average Ct value >35 or undetectable DTMUV RNA level was negative. (B) RT-LAMP-based CRISPR-Cas12a assay on tissue samples. Lanes 1–22 represent tissue samples. Lanes 23 and 24 represent the positive sample and negative control, respectively. In each figure, the red asterisk indicates the negative samples identified by the two assays.
4 Discussion
Molecular diagnostic methods, which rapidly detect DTMUV, facilitate the early detection of virus infection to avoid further outbreaks. Indeed, the PCR-based assays are more commonly used for virus detection among various diagnostic methods. Since the first outbreak of DTMUV in China in 2010, several PCR-based assays targeting the E, NS3 and NS5 genes have been developed to detect DTMUV, with limits of detection reaching 10 copies/μL, 1.82TCID50 and 0.001EID50/μL for viral RNA, respectively (Zang et al., 2012; Li, 2013; Wanzhe et al., 2016; Wang et al., 2017; Ninvilai et al., 2020). However, these PCR-based assays are not convenient for detecting a large number of clinical samples due to their complex steps and requirements for special equipment. Compared with them, the qPCR assays achieve high-throughput pathogen detection and provide a real-time readout via the qPCR instrument. Currently, some qPCR assays, such as one-step RT-qPCR assay targeting the 3′ terminal noncoding region (NCR), TaqMan-based qPCR assays targeting the E or NS1 gene and SYBR Green qPCR assay targeting the E gene, have been exploited to detect DTMUV RNA or cDNA, and these assays have limits of detection from 10 to 2.67 × 102 copies/μL (Yan L. P. et al., 2011; Yan et al., 2012; Yun et al., 2012; Wan et al., 2016; Chen G. Q. et al., 2018; Zhang et al., 2020; Li et al., 2023). Nonetheless, both PCR-based assays and qPCR assays take at least 2 to 4 h and require special equipment, which suggests that they are not convenient for the field detection of DTMUV. In addition, several RT-LAMP assays targeting the E or NS5 gene were also established to detect DTMUV owing to their rapidity, visualization and convenience, and their limits of detection were from 10 to 4 × 102 copies/μL for DTMUV RNA (Jiang et al., 2012; Wu et al., 2017; Wang et al., 2022). These assays are performed in a thermostatic water bath (61°C to 65°C), and their amplification can be completed within 50 min. However, the RT-LAMP assay easily causes false-positive results owning to cross-contamination and nonspecific amplification. Overall, these conventional detection methods, including PCR-based assays, qPCR assays and RT-LAMP assays, are difficult to widely use in clinical field detection due to their complex steps, requirements for special equipment or false-positive results.
In recent years, some studies have proven that CRISPR-Cas12a performs collateral cleavage on nontargeted ssDNAs when the Cas12a, crRNA and target DNA form a ternary complex (Li et al., 2018a). Under this condition, some isothermal amplification methods coupled with CRISPR-Cas12a had been established for pathogens detection (Chen J. S. et al., 2018; Li et al., 2018b; Tao et al., 2020; Hu et al., 2022). For example, an RPA-based CRISPR-Cas12a assay was developed for human papilloma virus (HPV) detection with a fluorescence readout, and its limited detection was 10 pM for HPV DNA (Chen J. S. et al., 2018). Within 1 h, this assay accurately identified the presence or absence of HPV in 25 patient samples. Meanwhile, another RT-LAMP-based CRISPR-Cas12a assay for hepatitis C virus (HCV) detection was also established, whose detectable concentration can reach 10 ng/μL for HCV RNA, and this assay can detect HCV within 1 h 30 min (Kham-Kjing et al., 2022). Therefore, these novel detection assays do not require special equipment and are more sensitive and specific than conventional PCR, qPCR and RT-LAMP assays. In this study, we developed an RT-LAMP-based CRISPR-Cas12a assay targeting the C gene for DTMUV detection. This assay shows several advantages over the existing methods based on conventional PCR, qPCR, and RT-LAMP. First, this assay only required some simple equipment, such as a water bath for isothermal amplification and CRISPR-Cas12a-mediated cleavage, and a blue or a UV light transilluminator for readout. Second, our RT-LAMP-based CRISPR-Cas12a assay has higher specificity than RT-LAMP solely because the Cas12a-crRNA complex specifically recognizes the target gene and discriminates nucleotide mismatches to avoid false-positive results (Chen J. S. et al., 2018). Third, our assay detects DTMUV RNA with a limited detection of 3 copies/μL, which is more sensitive than reported qPCR assays of DTMUV, including the SYBR Green qPCR assay (100 copies/μL) in our study (Yan L. P. et al., 2011; Yan et al., 2012; Yun et al., 2012; Wan et al., 2016; Chen G. Q. et al., 2018; Zhang et al., 2020; Li et al., 2023). Fourth, our assay can be completed within 2 h from RNA extraction to readout, which suggests that our assay is more suitable for rapid detection of DTMUV relative to conventional PCR and qPCR.
In addition, we tested 32 cloacal swabs and 22 tissue samples to further evaluate our RT-LAMP-based CRISPR-Cas12a assay. This assay showed 100% PPA and 100% NPA to SYBR Green qPCR for DTMUV detection in these clinical samples. Therefore, our assay can be applied for clinical field detection.
In summary, we successfully developed an RT-LAMP-based CRISPR-Cas12a assay to detect DTMUV RNA. This is the first report of DTMUV detection based on the combination of RT-LAMP with CRISPR-Cas12a. This novel method is suitable for clinical field detection owning to its visualization, high sensitivity and specificity.
Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary material.
Ethics statement
The animal study was approved by the Laboratory Animal Ethics Committee of South China Agricultural University. The study was conducted in accordance with the local legislation and institutional requirements.
Author contributions
YD: Data curation, Methodology, Software, Validation, Writing – original draft, Writing – review & editing, Investigation, Formal analysis. ZH: Writing – review & editing, Investigation. XL: Writing – review & editing, Investigation. MT: Writing – review & editing, Investigation. WL: Writing – review & editing, Investigation. SF: Writing – review & editing, Investigation. LZ: Writing – review & editing, Investigation. JZ: Writing – review & editing, Software. SY: Writing – review & editing, Investigation. FS: Writing – review & editing, Resources. PJ: Writing – original draft, Writing – review & editing, Conceptualization, Data curation, Formal analysis, Funding acquisition, Investigation, Methodology, Project administration, Resources, Supervision, Validation, Visualization.
Funding
The author(s) declare financial support was received for the research, authorship, and/or publication of this article. This work was funded by grant from the Laboratory of Lingnan Modern Agriculture Project (NT2021007).
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.
The author(s) declared that they were an editorial board member of Frontiers, at the time of submission. This had no impact on the peer review process and the final decision.
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.1301653/full#supplementary-material
References
Chen, G. Q., Li, Y. W., Wan, K. M., Chen, L., Qian, W., Zhu, T., et al. (2018). Establishment of a SYBR Green qPCR detection method for duck Tembusu virus. Anim. Husb. Vet. Med. 50, 78–81.
Chen, J. S., Ma, E. B., Harrington, L. B., Da Costa, M., 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, Y., Zong, N., Ye, F. D., Mei, Y. X., Qu, J. X., and Jiang, X. Y. (2022). Dual-CRISPR/Cas12a-assisted RT-RAA for ultrasensitive SARS-CoV-2 detection on automated centrifugal microfluidics. Anal. Chem. 94, 9603–9609. doi: 10.1021/acs.analchem.2c00638
Chertow, D. S. (2018). Next-generation diagnostics with CRISPR. Science 360, 381–382. doi: 10.1126/science.aat4982
de Sousa, K. A. F., Nonaka, C. K. V., Khouri, R., Rocha, C. A. G., Regis-Silva, C. G., and Souza, B. S. D. (2023). Rapid detection of SARS-CoV-2 based on the LAMP assay associated with the CRISPR-Cas12a system. Diagnostics 13:2233. doi: 10.3390/diagnostics13132233
Dronina, J., Samukaite-Bubniene, U., and Ramanavicius, A. (2022). Towards application of CRISPR-Cas12a in the design of modern viral DNA detection tools (review). J. Nanobiotechnology 20:41. doi: 10.1186/s12951-022-01246-7
Fan, Z. Y., Mei, Y. X., Xing, J. W., Chen, T., Hu, D., Liu, H., et al. (2023). Loop-mediated isothermal amplification (LAMP)/Cas12a assay for detection of Ralstonia solanacearum in tomato. Front. Bioeng. Biotechnol. 11:1188176. doi: 10.3389/fbioe.2023.1188176
Gootenberg, J. S., Abudayyeh, O. O., Kellner, M. J., Joung, J., Collins, J. J., and Zhang, F. (2018). Multiplexed and portable nucleic acid detection platform with Cas13, Cas12a, and Csm6. Science 360, 439–444. doi: 10.1126/science.aaq0179
Hamel, R., Phanitchat, T., Wichit, S., Vargas, R. E. M., Jaroenpool, J., Diagne, C. T., et al. (2021). New insights into the biology of the emerging Tembusu virus. Pathogens 10:1010. doi: 10.3390/pathogens10081010
Hongjaisee, S., Doungjinda, N., Khamduang, W., Carraway, T. S., Wipasa, J., Debes, J. D., et al. (2021). Rapid visual detection of hepatitis C virus using a reverse transcription loop-mediated isothermal amplification assay. Int. J. Infect. Dis. 102, 440–445. doi: 10.1016/j.ijid.2020.10.082
Hu, M. L., Qiu, Z. Q., Bi, Z. R., Tian, T., Jiang, Y. Z., and Zhou, X. M. (2022). Photocontrolled crRNA activation enables robust CRISPR-Cas12a diagnostics. Proc. Natl. Acad. Sci. U. S. A. 119:e2202034119. doi: 10.1073/pnas.2202034119
Jiang, T., Liu, J., Deng, Y. Q., Su, J. L., Xu, L. J., Liu, Z. H., et al. (2012). Development of RT-LAMP and real-time RT-PCR assays for the rapid detection of the new duck Tembusu-like BYD virus. Arch. Virol. 157, 2273–2280. doi: 10.1007/s00705-012-1431-7
Kham-Kjing, N., Ngo-Giang-Huong, N., Tragoolpua, K., Khamduang, W., and Hongjaisee, S. (2022). Highly specific and rapid detection of hepatitis C virus using RT-LAMP-coupled CRISPR-Cas12 assay. Diagnostics 12:1524. doi: 10.3390/diagnostics12071524
Lei, W. W., Guo, X. F., Fu, S. H., Feng, Y., Tao, X. Y., Gao, X. Y., et al. (2017). The genetic characteristics and evolution of Tembusu virus. Vet. Microbiol. 201, 32–41. doi: 10.1016/j.vetmic.2017.01.003
Li, G. P. (2013). Establishment and application of RT-PCR diagnosis method for duck Tambusu virus. China Anim. Husb. Vet. Med. 40, 69–71.
Li, S. Y., Cheng, Q. X., Liu, J. K., Nie, X. Q., Zhao, G. P., and Wang, J. (2018a). CRISPR-Cas12a has both cis- and trans-cleavage activities on single-stranded DNA. Cell Res. 28, 491–493. doi: 10.1038/s41422-018-0022-x
Li, S. Y., Cheng, Q. X., Wang, J. M., Li, X. Y., Zhang, Z. L., Gao, S., et al. (2018b). CRISPR-Cas12a-assisted nucleic acid detection. Cell Discov. 4:20. doi: 10.1038/s41421-018-0028-z
Li, H. Q., Wan, C. H., Wang, Z. Z., Tan, J., Tan, M. F., Zeng, Y. B., et al. (2023). Rapid diagnosis of duck Tembusu virus and goose astrovirus with TaqMan-based duplex real-time PCR. Front. Microbiol. 14:1146241. doi: 10.3389/fmicb.2023.1146241
Li, L. L., Zhang, Y., Dong, J. W., Zhang, J. Q., Zhang, C. H., Sun, M. H., et al. (2020). The truncated E protein of DTMUV provide protection in young ducks. Vet. Microbiol. 240:108508. doi: 10.1016/j.vetmic.2019.108508
Liang, M. D., Li, Z. L., Wang, W. S., Liu, J. K., Liu, L. S., Zhu, G. L., et al. (2019). A CRISPR-Cas12a-derived biosensing platform for the highly sensitive detection of diverse small molecules. Nat. Commun. 10:3672. doi: 10.1038/s41467-019-11648-1
Liang, T., Liu, X. X., Cui, S. L., Qu, S. H., Wang, D., Liu, N., et al. (2016). Generation of a reliable full-length cDNA of infectious Tembusu virus using a PCR-based protocol. Virus Res. 213, 255–259. doi: 10.1016/j.virusres.2015.12.022
Mao, G. B., Luo, X., Ye, S. L., Wang, X., He, J., Kong, J. L., et al. (2023). Fluorescence and colorimetric analysis of African swine fever virus based on the RPA-assisted CRISPR/Cas12a strategy. Anal. Chem. 95, 8063–8069. doi: 10.1021/acs.analchem.3c01033
Mayuramart, O., Nimsamer, P., Rattanaburi, S., Chantaravisoot, N., Khongnomnan, K., Chansaenroj, J., et al. (2021). Detection of severe acute respiratory syndrome coronavirus 2 and influenza viruses based on CRISPR-Cas12a. Exp. Biol. Med. 246, 400–405. doi: 10.1177/1535370220963793
Ninvilai, P., Tunterak, W., Prakairungnamthip, D., Oraveerakul, K., and Thontiravong, A. (2020). Development and validation of a universal one-step RT-PCR assay for broad detection of duck Tembusu virus. Avian Dis. 64, 294–299. doi: 10.1637/aviandiseases-D-19-00199
Qin, C., Liu, J. J., Zhu, W. Q., Zeng, M. C., Xu, K., Ding, J. M., et al. (2022). One-pot visual detection of African swine fever virus using CRISPR-Cas12a. Front. Vet. Sci. 9:962438. doi: 10.3389/fvets.2022.962438
Soroka, M., Wasowicz, B., and Rymaszewska, A. (2021). Loop-mediated isothermal amplification (LAMP): the better sibling of PCR? Cells 10:1931. doi: 10.3390/cells10081931
Tang, Y., Diao, Y., Gao, X., Yu, C., Chen, L., and Zhang, D. (2012). Analysis of the complete genome of Tembusu virus, a Flavivirus isolated from ducks in China. Transbound. Emerg. Dis. 59, 336–343. doi: 10.1111/j.1865-1682.2011.01275.x
Tao, D. G., Liu, J. J., Nie, X. W., Xu, B. R., Tran-Thi, T. N., Niu, L. L., et al. (2020). Application of CRISPR-Cas12a enhanced fluorescence assay coupled with nucleic acid amplification for the sensitive detection of African swine fever virus. ACS Synth. Biol. 9, 2339–2350. doi: 10.1021/acssynbio.0c00057
Thontiravong, A., Ninvilai, P., Tunterak, W., Nonthabenjawan, N., Chaiyavong, S., Angkabkingkaew, K., et al. (2015). Tembusu-related Flavivirus in ducks, Thailand. Emerg. Infect. Dis. 21, 2164–2167. doi: 10.3201/eid2112.150600
Tunterak, W., Prakairungnamthip, D., Ninvilai, P., Bunyapisitsopa, S., Oraveerakul, K., Sasipreeyajan, J., et al. (2018). Serological evidence of duck Tembusu virus infection in free-grazing ducks, Thailand. Transbound. Emerg. Dis. 65, 1943–1950. doi: 10.1111/tbed.12975
Wan, C. H., Chen, C. T., Fu, Q. L., Cheng, L. F., Fu, G. H., Chen, H. M., et al. (2016). Development of a TaqMan-based real-time PCR method for avian Tembusu virus. J. Northwest A&F Univ. 44, 49–54. doi: 10.13207/j.cnki.jnwafu.2016.08.008
Wang, Q. P., Su, K. H., Yuan, Z. J., Chang, Z. S., Zhang, W., Zhang, Y. F., et al. (2022). Establishment and preliminary application of a rapid real-time fluorescent RT-LAMP for duck Tembusu virus. China Anim. Health Inspect. 2022, 39, 103–109. doi: 10.3969/j.issn.1005-944X.2022.05.020
Wang, Y. J., Zhu, S. Y., Hong, W. M., Wang, A. P., and Zuo, W. Y. (2017). A multiplex PCR for detection of six viruses in ducks. J. Virol. Methods 248, 172–176. doi: 10.1016/j.jviromet.2017.07.004
Wanzhe, Y., Jianuan, L., Peng, L., Jiguo, S., Ligong, C., and Juxiang, L. (2016). Development of a nano-particle-assisted PCR assay for detection of duck Tembusu virus. Lett. Appl. Microbiol. 62, 63–67. doi: 10.1111/lam.12509
Wu, Z., Xia, W. L., Cai, S. D., Guo, C., Yuan, W., Wang, Y., et al. (2017). Establishment and application of visual RT-LAMP for rapid detection of duck Tembusu virus. China Anim. Husb. Vet. Med. 44, 3334–3339. doi: 10.16431/j.cnki.1671-7236.2017.11.032
Wu, S. Y., Yuan, J., Xu, A., Wang, L., Li, Y. B., Lin, J. H., et al. (2023). A lab-on-a-tube biosensor combining recombinase-aided amplification and CRISPR-Cas12a with rotated magnetic extraction for Salmonella detection. Micromachines 14:830. doi: 10.3390/mi14040830
Xu, L. Q., Dai, Q. Q., Shi, Z. Y., Liu, X. T., Gao, L., Wang, Z. Z., et al. (2020). Accurate MRSA identification through dual-functional aptamer and CRISPR-Cas12a assisted rolling circle amplification. J. Microbiol. Methods 173:105917. doi: 10.1016/j.mimet.2020.105917
Xu, B. F., Liu, X. P., Yan, D. W., Teng, Q. Y., Yuan, C. X., Zhang, Z. F., et al. (2023). Generation and characterization of chimeric Tembusu viruses containing pre-membrane and envelope genes of Japanese encephalitis virus. Front. Microbiol. 14:1140141. doi: 10.3389/fmicb.2023.1140141
Yan, C., Cui, J., Huang, L., Du, B., Chen, L., Xue, G., et al. (2020). Rapid and visual detection of 2019 novel coronavirus (SARS-CoV-2) by a reverse transcription loop-mediated isothermal amplification assay. Clin. Microbiol. Infect. 26, 773–779. doi: 10.1016/j.cmi.2020.04.001
Yan, L. P., Peng, S., Yan, P. X., Zhou, J. W., Teng, Q. Y., Li, G. X., et al. (2012). Comparison of real-time reverse transcription loop-mediated isothermal amplification and real-time reverse transcription polymerase chain reaction for duck Tembusu virus. J. Virol. Methods 182, 50–55. doi: 10.1016/j.jviromet.2012.03.007
Yan, Z., Shen, H., Wang, Z., Lin, W., Xie, Q., Bi, Y., et al. (2017). Isolation and characterization of a novel Tembusu virus circulating in Muscovy ducks in South China. Transbound. Emerg. Dis. 64, E15–E17. doi: 10.1111/tbed.12525
Yan, L. P., Yan, P. X., Zhou, J. W., Teng, Q. Y., and Li, Z. J. (2011). Establishing a TaqMan-based real-time PCR assay for the rapid detection and quantification of the newly emerged duck Tembusu virus. Virol. J. 8:464. doi: 10.1186/1743-422X-8-464
Yan, P. X., Zhao, Y. S., Zhang, X., Xu, D. W., Dai, X. G., Teng, Q. Y., et al. (2011). An infectious disease of ducks caused by a newly emerged Tembusu virus strain in mainland China. Virology 417, 1–8. doi: 10.1016/j.virol.2011.06.003
Yu, G. L., Lin, Y., Tang, Y., and Diao, Y. X. (2018). Evolution of Tembusu virus in ducks, chickens, geese, sparrows, and mosquitoes in northern China. Viruses 10:485. doi: 10.3390/v10090485
Yun, T., Ni, Z., Hua, J. G., Ye, W. C., Chen, L., Zhang, S., et al. (2012). Development of a one-step real-time RT-PCR assay using a minor-groove-binding probe for the detection of duck Tembusu virus. J. Virol. Methods 181, 148–154. doi: 10.1016/j.jviromet.2012.01.019
Zang, S., Yun, T., Ye, W. C., Deng, X. H., Cui, Y. S., Hua, J. G., et al. (2012). Development and application of one step RT-PCR technique for detection of duck Tembusu virus. Acta Agric. Zhejiangensis 21, 37–40.
Zhang, W., Chen, S., Mahalingam, S., Wang, M. S., and Cheng, A. C. (2017). An updated review of avian-origin Tembusu virus: a newly emerging avian Flavivirus. J. Gen. Virol. 98, 2413–2420. doi: 10.1099/jgv.0.000908
Zhang, X. Y., Yao, M., Tang, Z. H., Xu, D. N., Luo, Y., Gao, Y. F., et al. (2020). Development and application of a triplex real-time PCR assay for simultaneous detection of avian influenza virus, Newcastle disease virus, and duck Tembusu virus. BMC Vet. Res. 16:203. doi: 10.1186/s12917-020-02399-z
Zhao, F. R., Wang, P., Wang, H. X., Liu, S. R., Sohail, M., Zhang, X., et al. (2023). CRISPR/Cas12a-mediated ultrasensitive and on-site monkeypox viral testing. Anal. Methods 15, 2105–2113. doi: 10.1039/d2ay01998a
Keywords: duck, Tembusu virus, RT-LAMP, CRISPR-Cas12a, clinical field detection
Citation: Ding Y, Huang Z, Li X, Tang M, Li W, Feng S, Zhao L, Zhang J, Yuan S, Shan F and Jiao P (2023) Development of a reverse transcription loop-mediated isothermal amplification based clustered regularly interspaced short palindromic repeats Cas12a assay for duck Tembusu virus. Front. Microbiol. 14:1301653. doi: 10.3389/fmicb.2023.1301653
Edited by:
Chengming Wang, Auburn University, United StatesReviewed by:
Lanjing Wei, University of Kansas, United StatesLiu Sidang, Shandong Agricultural University, China
Copyright © 2023 Ding, Huang, Li, Tang, Li, Feng, Zhao, Zhang, Yuan, Shan and Jiao. 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: Fen Shan, c2hhbmZlbl9nenpvb0BzaW5hLmNvbQ==; Peirong Jiao, cHJqaWFvQHNjYXUuZWR1LmNu