Kai Sun1†
Yi Liu1†Xin Zhou2Chuanlin Yin1Pengjun Zhang1Qianqian Yang1Lingfeng Mao3
Xuping Shentu1*
Xiaoping Yu1*- 1Zhejiang Provincial Key Laboratory of Biometrology and Inspection and Quarantine, College of Life Sciences, China Jiliang University, Hangzhou, China
- 2Ausper Biopharma, Hangzhou, China
- 3Hangzhou Baiyi Technology Co., Ltd., Hangzhou, China
Plant viruses threaten crop yield and quality; thus, efficient and accurate pathogen diagnostics are critical for crop disease management and control. Recent advances in sequencing technology have revolutionized plant virus research. Metagenomics sequencing technology, represented by next-generation sequencing (NGS), has greatly enhanced the development of virus diagnostics research because of its high sensitivity, high throughput and non-sequence dependence. However, NGS-based virus identification protocols are limited by their high cost, labor intensiveness, and bulky equipment. In recent years, Oxford Nanopore Technologies and advances in third-generation sequencing technology have enabled direct, real-time sequencing of long DNA or RNA reads. Oxford Nanopore Technologies exhibit versatility in plant virus detection through their portable sequencers and flexible data analyses, thus are wildly used in plant virus surveillance, identification of new viruses, viral genome assembly, and evolution research. In this review, we discuss the applications of nanopore sequencing in plant virus diagnostics, as well as their limitations.
Introduction
Because plant viruses cause devastating diseases and enormous economic losses worldwide in agricultural systems every year, they are major threats to sustainable and productive agriculture (Scholthof et al., 2011). Virus diseases are characterized by various symptoms, including ringspots, mosaic pattern development, leaf yellowing and distortion, plus impaired growth. However, diagnosis based on symptoms is unreliable because it lacks specificity. For example, some symptoms are quite subtle and can be easily confused with nutrient deficiencies and herbicide injury. Additionally, mixed infections often cause nonspecific and severe symptoms that further complicate virus identification. The expansion of international trade in agricultural products has accelerated the emergence and spread of plant viruses (Elena et al., 2014; Gilbertson et al., 2015). Therefore, rapid and effective detection of plant viruses is key to preventing and controlling infections.
Traditionally, plant viruses have been detected using morphological, PCR, and immunological analyses. However, these methods have several drawbacks, including long detection times, complicated preprocessing steps, and low detection efficiencies. Furthermore, these tests often require prior knowledge of viral genomic information or immunological characteristics that are not available for novel viruses.
Currently, metagenomic sequencing, represented by next-generation sequencing (NGS), has been widely used for detection and identification of plant viruses (Hadidi et al., 2016; Roossinck, 2017). NGS employs sequencing by synthesis, which determines the DNA sequences by capturing markers on newly-synthesized nucleotides (Goldberg et al., 2015). These methods do not require previous knowledge of viral sequences and can be used to sequence millions or billions of nucleotides in parallel. This permits detection of all viruses, including those not previously described (Mehetre et al., 2021). In theory, NGS can detect all viruses in a single assay and performance is only limited by the completeness of reference databases against which the sequences are compared. The resulting sequence information can also provide insight into virus population structure, ecology and evolution, as well as reveal virus variants that may contribute to various disease etiologies (Maree et al., 2018). However, NGS technology is not applicable for field detection of viruses, and although it can generate large datasets, the quality of sequencing needs improvement because the error rate of conventional NGS is about 1% and the length of sequence reads is generally short (35–700 bp)when compared with Sanger sequencing platforms (Fox et al., 2014; Kumar et al., 2019). Moreover, NGS sample preparation and actual sequencing are time-consuming and NGS equipment is expensive and often bulky (Gu et al., 2019), which limits widespread adoption for on-site pathogen detection or field deployment.
Nanopore sequencing, developed by Oxford Nanopore Technologies (ONT), is a third-generation sequencing technology that overcomes NGS technology shortcomings, by combining genetic engineering technology with computer based analyses (Lin et al., 2021). ONT nanopore sequencing has several advantages, including single-molecule sequencing, long sequencing read lengths, rapid sequencing speeds, and real-time monitoring of sequencing data (Laver et al., 2015; Deamer et al., 2016). Table 1 shows the performance of four high-throughput sequencing platforms. The most widely accepted nanopore sequencing platform is ONT MinION, a highly portable instrument (10 cm × 2 cm × 3.3 cm, 90 g), powered via a computer USB port (Figure 1). Unlike other sequencing technologies, sequencing and real-time data analysis can be completed on a personal computer. Coupled with its portability, nanopore technology has unprecedented applicability for plant virus detection. Here, we review the technical principles of nanopore sequencing and introduce applications for on-site, rapid, reliable, and sensitive plant virus detections.
Figure 1. MinION sequencing device. MinION is a highly portable instrument (10 cm × 2 cm × 3.3 cm, about 90 g; B). It can be connected to an ONTs “MinIT” decoder (C) or to any computer with a USB port (A). The entire nanopore sequencing process and real-time data analysis can be done on a PC. Nanopore sequencing is achieved by adding a sequencing library to the flowCell sequencing chip (D). When DNA or RNA molecules pass through the nanopore, there is a shift in the nanopore current, which is measured by a sensor. The data is further transmitted and decoded into base signals that determine the nucleic acid sequence.
Table 1. Comparison of the performance of four high-throughput sequencing platforms.
Principle and advantages of nanopore sequencing technology
ONT nanopore sequencing is a promising, highly sensitive platform for identification of single strand DNA or RNA molecules. Unlike previous sequencing technologies, ONT nanopore sequencing identifies nucleotides based on electrical signals rather than optical signals (Deamer et al., 2016). This technology relies on a protein that forms a channel (nanopore) on a membrane. Upon application of an electrophoretic force, negatively charged biomolecules such as DNA or RNA pass through the nanopores and temporarily affecting the strength of the current flowing through the nanopore. The signal is then recorded and further analyzed to determine the nucleic acid sequence. The conversion of electrical signals into base sequences is called basecalling. It is the last step in the analysis of entire DNA/RNA sequences and uses deep learning to eliminate sequencing noise and signal errors. This step is crucial for improving sequencing accuracy and requires efficient algorithms and large datasets for computational training (Lin et al., 2021).
Nanopore sequencing allows direct DNA and RNA sequencing that does not require DNA synthesis reactions because nucleotides are recorded directly from the current signals. Whole-genome sequencing and methylation calling using Nanopore sequencing has been carried out with different species (Xu and Seki, 2020). For example, Kim et al. (2020) used direct RNA sequencing to obtain the COVID-19 genome sequence and detected at least 41 RNA modification sites on viral transcripts.
Nanopores can generate very long read lengths because the entire DNA or RNA fragment is analyzed. Thus, the read length directly correlates with the length of the DNA or RNA fragment. Users can choose appropriate nucleic acid preparation methods based on their experimental goals. Data released at the 2020 Nanopore scientific society conference showed that the longest sequence read lengths from nanopore sequencing have reached 4.1 Mb can be obtained from a single sequence read. Ultra-long read lengths are associated with higher quality and precision (Petersen et al., 2019). Nanopore sequencing permits real-time sequencing. Thus, the user can obtain data and experimental quality analysis as soon as sequencing begins, monitor the status of the sequencing, and stop the process as soon as sufficient data is available. Real-time sequence alignments can also be performed using the EPI2ME workflow, which enables rapid identification of pathogen-related sequences in the samples. Additionally, ONT tools1 are available to the user. These include Medaka, Tombo, Pomoxis, and Nanopolishare, which provide sequence correction, identification of modified nucleotides, genome assembly, and error correction in the genome assembly (Lu et al., 2016; Payne et al., 2019). Real-time Nanopore sequencing has a good application prospect in plant virus detection. Because early identification of viruses allows growers to take quick actions and effective sanitary measures. Real-time and on-site sequencing reduces overall costs and gives crop protection officers and farmers in rural communities’ information that is critical for sustainable crop production and management of pests and diseases (Boykin et al., 2019).
ONT provide a variety of sequencing platforms to meet various sequencing requirements, including MinION, a portable sequencer from ONT weighing about 100 g. The high throughput sequencers, GridION (which uses five MinION chips at a time) or PromethION (with 48 chips containing 3,000 nanopore channels each) allows for customization of chip number. Flongle, a one-off sequencing platform enables direct, real-time DNA or RNA sequencing and costs about $90, compared with the $900 MinION sequencing chip, and can generate about 1 GB of data in 24 h (Petersen et al., 2019). Thus, Flongle is a cost-effective and rapid option for smaller sequencing experiments.
The portability and real-time sequencing advantages of MinION provide an ideal platform for field sequencing because the method can provide results within 24 h of receiving samples with sequencing taking 15 to 60 min. In 2015, (Quick et al., 2016). shipped a nanopore sequencer in standard airline luggage to Guinea and conducted real-time monitoring of the West African Ebola outbreak. In a different trial, Boykin et al. (2019) detected cassava mosaic virus infections through metagenomics sequencing of cassava leaf, stem, and tuber samples using MinION at a cassava farm in the south of the Sahara desert. The whole process (collection of infected plant samples, nucleic acid database construction and sequencing analyses) was completed within 3 h after the researchers arrived at the farm.
Methods for plant virus detection
In general, nanopore sequencing involves wet and dry lab steps. Wet lab work includes nucleic acid extraction and sequencing library preparation, which may take several hours. Dry lab work includes bioinformatics analysis, which can be completed within hours or days depending on the amount of data to be processed.
The wet lab process
Viral nucleic acids, including genomic DNA or RNA, transcripts, and replicative intermediates can be noticed by a nanopore sequencer. For plant DNA viruses, several library preparation strategies have been established by ONT. For example, DNA from plant extracts can be directly sequenced using a ligation sequencing kit (SQK-LSK109) or amplified using PCR to obtain sufficient amplicons using a rapid barcoding kit (SQK-RBK004). Plant DNA viruses, both the ssDNA geminiviruses and the reverse-transcribing pararetroviruses, have circular genomes. Researchers often use roll-circle amplification (RCA) to enrich virus-derived DNA fragments and then break the circles into linear DNA strands before sequencing on a nanopore sequencing platform (Figure 2).
Figure 2. Overview of Nanopore sequencing and library preparation for plant virus detection. (A) Workflow for plant virus detection from sample collection to virus identification. (B) Direct RNA sequencing using the ONT-SQK-RNA001 library method with poly (A) RNA as a template. A reverse transcription step with oligo-dT primers was used to circumvent secondary structure of the RNA. A sequencing adapter was ligated to the mRNA using T4 DNA ligase. Since only the RNA strand is motor-ligated, only the RNA molecule is sequenced. (C) cDNA-PCR sequencing using the ONT-SQK-PCS108 library preparation method. Poly (A) RNA was used as a template for first strand cDNA synthesis with oligo-dT30VN primers (VNP). When first strand cDNA synthesis reaches the end of the RNA molecule, few non-templated Cs were added to the end of the cDNA by the reverse transcriptase. Then, a strand-switching primer (SSP) present in the reaction binds to the non-templated Cs, followed by enrichment. (D) DNA sequencing: the ONT-SQK-LSK108 library preparation method. The Barcoding Kit such as ONT-EXP-NBD103 can be used to tag the native DNA or amplicons DNA (after PCR or RCA reactions).
However, nanopore methods that detect DNA cannot detect RNA viruses directly. Thus, RNA is directly sequenced or first converted into cDNA before being used to detect plant viruses. Current ONT library preparation kits include a direct cDNA sequencing kit (SQK-DCS109), a PCR-cDNA sequencing kit, with (SQK-PCB109) and without (SQK-PCS109) barcoding, and a direct RNA sequencing (SQK-RNA002) kit (Liefting et al., 2021). These kits are designed for polyadenylated [poly(A)] RNA samples (Figure 2). Hence, plant viruses that lack poly(A) tails, cannot be sequenced using these strategies. However, multiple strategies have been developed for sequencing RNA that lack poly (A) tails. In one strategy, a poly (A)-tailing reaction is carried out on the total RNA using E. coli poly (A) polymerase and then using the resulting poly (A) RNA as input for the ONT direct cDNA sequencing kit (SQK-DCS109). In another strategy, double-stranded (ds) cDNA is synthesized using random hexamers and used as input in the end-prep step of the direct cDNA sequencing kit (SQK-DCS109). These methods have proven to be efficient in viral metagenomics and enable the use of nanopore sequencing to detect all types of plant viruses (Palanga et al., 2016; Claverie et al., 2018).
The dry-lab process
A number of bioinformatics tools have been used for virus sequence analysis and diagnosis. The virus detection workflow needs to: (1) parse input electropherogram files to obtain base sequences, (2) conduct quality control measures on raw data files, including trimming of poor quality reads and adaptor sequences, and (3) assemble and map reads in order to identify known and novel viruses (Figure 3). Recent software for the analysis of nanopore data are shown on Table 2.
Table 2. Bioinformatics tools for the identification of plant viruses using nanopore sequencing.
Figure 3. Outline of potential stages in the workflow for nanopore-seq analysis for plant virus detection. Representative bioinformatics tools are shown in the workflow.
Examples of plant virus detection using nanopore sequencing
Although Nanopore sequencing technology is relatively new, several studies have used the technology to detect plant viruses (Table 3).
Table 3. Summery of viruses detected in plants using nanopore sequencing technology.
Peach (Prunus persica)
Through whole transcriptome analysis with the Oxford MinION technology, plum pox virus (PPV) was detected in P. persica (peach; Bronzato Badial et al., 2018). Briefly, total RNA was extracted from plant samples and amplified using a TransPlex WTA kit (Sigma Aldrich, St. Louis, MO). A cDNA library was then prepared and loaded into the MinION flow cell and sequenced for 24 h and the raw electrical data was processed with Metrichor (Oxford Nanopore). FASTA sequences and read event data, including strand translocation times, were extracted and analyzed by the R package, poRe (Watson et al., 2015). The GraphMap Aligner was used for nucleotide base alignment (Sović et al., 2016) and the resulting contigs were compared with the complete set of high quality, annotated, viral genomes on NCBI RefSeq.2 This method permitted the simultaneous detection of PPV in peach plant after a 38 s run, demonstrating the capability of this methodology to obtain results rapidly. Results showed that this methodology is useful for detecting unsuspected viral pathogens in plant.
Cassava (Manihot manifera)
Two studies Nanopore sequencing have investigated virus infections in cassava. Boykin et al. (2019) combined PDQeX DNA purification technology with the MinION and MinIT mobile sequencing devices to generate an effective point-of-need field diagnostic system with cassava samples from Tanzania, Uganda, and Kenya. Barcoded DNA libraries were generated with the rapid barcoding kit, SQK-RBK004 DNA sequencing was conducted with a MinION connected to MinIT, a portable live base calling server. Reads across all of the sequencing runs (mean read length: 35–948 bp) were compared with a cassava mosaic disease reference database using BLASTn. The study detected several viruses in real time, including African cassava mosaic virus (ACMV) and the East African cassava mosaic virus (EACMV). This case study was designed to show the possibility to go from sample to diagnosis, in a regional setting, on farm in 3 h versus the normal 6 months with conventional methods. The results of this research show that it is possible to use a range of battery powered devices to achieve DNA extraction, long read sequencing and analysis all under a tree on the farm while the farmers wait for results (Boykin et al., 2019). In the other study, nanopore sequencing detected Sri Lankan cassava mosaic virus (SLCMV) in cassava samples from Thailand (Leiva et al., 2020). In this case, total DNA was extracted from an infected cassava plant with cetyltrimethylammonium bromide and Phi29 DNA polymerase was used for RCA. The DNA library was prepared with a ligation sequencing kit (SQK-LSK109) and then sequenced in a FLO-MIN106D (R9.5) flow cell. The assembly was assembled from the raw reads via de novo assembly (using Canu v1.8) and reference assembly (using Minimap2 Li, 2018) and Pilon (Walker et al., 2014) with distinct SLCMV genome sequences. The resulting contigs contained 72,800 reads for DNA-A, with an average coverage of 15,000 X, and 70,681 reads for DNA-B, with an average coverage of 6,000 X. This study provides a well-arranged sequencing process for obtaining complete geminivirus genomes using nanopore technology.
Yam (Dioscoriaalata)
Multiple virus infections have been identified in yams (Dioscoriaalata) with a nanopore cDNA-seq analysis (Filloux et al., 2018) in which a MinION sequencing library was prepared with the SQK-PCS108 cDNA-PCR kit (Oxford Nanopore Technologies Ltd). Briefly, cDNA was first generated from total RNA from an infected yam by revers transcription (strand-switching method) and then amplified by PCR. The PCR product was ligated with adapters, loaded onto a R9.4 Flow Cell (FLO-MIN106 R9.4) and sequenced for 48 h with MinION. The MinION reads were assembled de novo using Canu v1.6 (Koren et al., 2017). and compared against the GenBank database with DIAMOND (Buchfink et al., 2015). The study also compared the use of small RNA Illumina sequencing with nanopore sequencing for virus detection. The two methods detected three viruses, Dioscorea bacilliform virus (DBV), yam mild mosaic virus (YMMV) and yam chlorotic necrosis virus (YCNV). However, MinION sequencing failed to detect ampelovirus-related sequences, which were detected by Illumina. The study determined that this failure resulted from the fact that the Nanopore cDNA sequencing protocol used specifically targeted poly (A) sequences and hence missed viruses that lack 3′ poly (A) tails. Furthermore, this research shows that the consensus sequence obtained either by de novo assembly or after mapping the MinION reads on the virus genomic sequence was >99.8% identical with the Sanger-derived reference sequence (Filloux et al., 2018). These degrees of sequencing accuracy demonstrate that the Nanopore sequencing approach can be used to both reliably detect and accurately sequence nearly full-length virus genomes with positive-sense single-strand polyadenylated RNA.
Cowpea (Vigna unguiculata)
Nanopore sequencing technology was used to identify a novel bipartite begomovirus in cowpea plants which has been designated cowpea bright yellow mosaic virus (CoBYMV; Naito et al., 2019). In this case, total DNA was extracted from cowpea samples exhibiting bright golden mosaic symptoms and subjected to RCA using Phi29 DNA polymerase. The DNA library was prepared with the ligation kit 1D and sequenced with a Spot-on flowcell (FLO-MIN106) on the MinION Mk1B. The raw FAST5 reads were uploaded to the online server for base calling through the Metrichor EPI2ME platform. The base-called MinION reads were converted to FASTQ and assembled using Canu v1.7 (Koren et al., 2017). The resulting contigs were then compared against the viral protein RefSeq database using Blastx in Geneious version 9.1.2 (Kearse et al., 2012). The assembled contigs of the putative complete DNA-A and DNA-B sequences were confirmed by PCR and Sanger sequencing. The MinION derived consensus sequences had 100% nucleotide sequence identity with the corresponding Sanger sequences, except for one nucleotide in a T-rich region of DNA B. This study also demonstrated that a portable nanopore sequencing device is a rapid and accurate alternative tool for the characterization of novel plant virus.
Wheat (Triticum aestivum)
Wheat streak mosaic virus (WSMV), triticum mosaic virus (TriMV), and high plains wheat mosaic virus (HPWoMV) together comprise the wheat streak mosaic complex (WSM), which causes losses of up to 5% of U.S (Fellers et al., 2019). wheat production. Use of resistant wheat cultivars consisting primarily of introgressed Wsm1 and Wsm2 resistance genes is the most economical and effective method for disease control. However, in 2015, the ‘Clara CL’ resistant variety harboring Wsm2 resistance developed severe mosaic symptoms and stunted growth in Hamilton County, KS (Fellers et al., 2019). To investigate virus variants that could circumvent the Wsm2 resistance, infected wheat tissue was sequenced by nanopore sequencing with the 1D cDNA sequencing kit (SQK-LSK108) in two cDNA reactions, one of which used XhoI-oligo-d(T)20 for poly (A) viruses and the other used random 6-mer primers for poly (A) minus RNA viruses. A DNA library was prepared from 500 ng of amplified cDNA and the samples were loaded onto a MinION 107 v9.5 Flow Cell. Raw data base-calling and adapter trimming was carried out by MinKNOW, with albacore v1.7.3. and porechop v0.2.3 (Wick, R, University of Melbourne). Clean data was then aligned to a cereal virus reference file with CLC Genomics Workbench v11 (Qiagen). Non-wheat virus assemblies were compared against the NCBI’s virus genome database by BLAST analyses, which revealed the presence of mixed infections of WSMV, TriMV and BYDV-PAV. The study also revealed that susceptibility resulted from a nucleotide sequence difference in the WSMV field isolates at position 6,833 that resulted in threonine (T) substitutions for valine (V) or methionine (M). One or more of these differences may explain why resistant wheat cultivars exhibited WSMV symptoms. These results demonstrate that ONT can more accurately identify causal virus agents and has sufficient resolution to provide evidence of causal variants.
Tomato (Solanum lycopersicum), cucumber (Cucumis sativus) and butternut squash (Cucurbita moschata)
To evaluate Nanopore sequencing for diagnosis of plant diseases, a recent study described DNA or RNA sequencing of symptomatic plant tissues infected with known viruses (tomato with tomato yellow leaf curl virus or tomato brown rugose fruit virus, watermelon samples with watermelon chlorotic stunt virus or an unknown virus that elicited yellowish leaves and mottling, and butternut squash with mosaic leaf symptoms; Chalupowicz et al., 2019). A DNA library was prepared from nucleic acid extracted from these plants with the ligation sequencing kit 1D (SQK-LSK 108). Direct RNA library preparations were carried out with the Oxford Nanopore direct RNA sequencing protocol. The sequences were run on the MinION flow cell (FLO-Min106 version R9.4, Oxford Nanopore Technologies) and assembled on the MinION sequencer (MK 1B version, Oxford Nanopore Technologies). MinION basecalling was conducted for 6 h in real-time on a local computer with MinKNOW 2.0 version 18.03.1 software (Oxford Nanopore Technologies Ltd.). The basecalled reads were then uploaded by use of EPI2ME desktop agent software (version 2.52.1202033). The “What’s in my Pot?” (WIMP) workflow was used to assign taxonomic classifications of the MinION reads based on the NCBI RefSeq virus database.3 DNA viruses were identified within the DNA sequencing samples, and the CGMMV and ZYMV RNA viruses were detected during RNA sequencing of unknown samples. These findings highlight the ability of the Nanopore platform to detect viruses through direct DNA or RNA sequencing.
Summary and future prospective
Nanopore sequencing technology allows direct sequencing of DNA or RNA samples and provides fast and real-time dynamic monitoring of sequencing data in the field. However, compared with plant bacterial and fungal pathogens, the detection of plant viruses using nanopore sequencing faces greater challenges. As we all known that the genomic meterials of fungi and bacteria are DNA. In the detection of plant bacteria and fungi, DNA direct sequencing or marker gene DNA amplification sequencing method is mainly used. The most commonly used marker genes in metataxonomic are the 16S rRNA for bacteria and archaea, and the 18S rRNA for fungi, while there are no marker genes for viruses (Ciuffreda et al., 2021). In addition, viruses have more diverse genome types classified into DNA viruses and RNA viruses. Furthermore, RNA genomes can be divided into those carrying Poly (A) tails and those without Poly(A) tails. There is currently no general library construction method to achieve sequencing of all types of plant viruses using nanopore sequencing technology. So, researchers must choose a suitable library construction strategy according to the characteristics of target viruses.
To improve field detection efficiency, pre-screened databases that only contain data of the suspected disease and host genome information are often used (Boykin et al., 2019). Consequently, existing or new field viruses may not be identified depending on the status of the databases until subsequent data analysis when the scientist has returned to the lab or is within range of a good internet connection capable of uploading large amounts of data to the cloud. Nonetheless, the current nanopore sequencing technology is the most economical and convenient sequencing platform for detection of plant viruses in the field. However, the use of nanopore sequencing to detect human pathogens is at a more advanced stage compared with plant viruses. Therefore, we have discussed the latest application of nanopore sequencing technology in clinical diagnosis of pathogens with a view towards providing information to improve the detection of plant viruses.
The combination of nanopore sequencing technology and isothermal amplification technology has been shown to shorten detection times and improve detection efficiency. Boykin et al. (2019) developed a simple and accurate molecular diagnosis method combining Loop-mediated isothermal amplification (LAMP) and nanopore sequencing for detection of malaria (Imai et al., 2017). This protocol included 18S rRNA specific LAMP primers for human plasmodium and tested blood samples from 63 malaria patients. The LAMP product was sequenced using MinION, which revealed that the sequence obtained was consistent with that of the reference plasmid sequence and nested PCR results. Bi et al. (2021) developed a method for multiplex isothermal amplification-based sequencing and real-time analysis of multiple viral genomes, referred to as the nanopore sequencing of isothermal rapid viral amplification for near real-time analysis. The design of the method provided the ability to simultaneously detect SARS-CoV-2, influenza A, human adenovirus, and human coronavirus, as well as to identify mutations in up to 96 samples in real time. These methods can avoid the false positives caused by isothermal amplification, and also increase the content of targeted nucleic acids in the sequencing library.
The application of bioinformatics tools in field applications has continuously improved sequencing accuracy, read length, and throughput. For example, a new software named UNCALLED developed by Kovaka et al. (2021) has improved the sequencing of target nucleic acid amounts by manipulating the nanopore current. The software has the potential to quickly match the nanopore current signal flow to a preset sequence reference database. This allows matching targeted nucleic acids to pass through the nanopore to obtain the sequences of full-length molecules. In cases where the sequence does not match the reference sequence, the software will reverse the voltage across the nanopore and physically eject the nucleic acid molecules to create room for the next nucleic acid. In the study, UNCALLED was employed to eliminate the sequence information of known bacterial genomes in metagenomic populations, in order to enrich the nucleic acids of other species by nearly five folds. A total of 148 human genes associated with hereditary cancer were enriched by nearly 30 folds, and the numbers of single-nucleotide polymorphisms, insertions, and deletions, structural variations, and methylation variants in the analyses were more than twice the number of those detected by the 50X coverage in short read-lengths of the whole genome sequencing.
High-throughput sequencing can be used to characterize virus evolution and to identify adaptations that affect transmission or pathogenicity. Coupled with a high replication rate, plant RNA viruses can form a virus population containing a group of genetically related but different haplotypes. Many attempts have been made to reconstruct viral haplotypes using NGS reads. However, the short length of NGS reads cannot cover distant single-nucleotide variants, making it difficult to reconstruct complete or near-complete haplotypes (Cai and Sun, 2022). Long-read ONT data can be used to study sample-specific genomic details, including structural variation (SV) and haplotypes. Given that single long reads can encompass multiple variants, including both single-nucleotide variants (SNVs) and SVs, it is possible to perform phasing of haplotype-resolved analyses with appropriate bioinformatics software, such as LongShot for SNV detection and WhatsHap for haplotyping/phasing (Wang et al., 2021).
Although single-nucleotide variations (SNVs) are considered to be key drivers of virus adaptation, RNA recombination events that delete or insert nucleic acid sequences have also been implicated in this process. Jaworski et al. (2021) designed an approach called ‘Tiled-ClickSeq’ to provide complete genome coverage, including the 5′UTR, at high resolution and virus specificity on both Illumina and Nanopore platforms. Using the developed method, these platforms simultaneously analyzed multiple SARS-CoV-2 isolates and clinical samples and characterized minority variants, sub-genomic mRNAs, structural variants, and D-RNAs. This indicates that the Tiled-ClickSeq is a convenient and robust platform for SARS-CoV-2 genomics because it captures a full range of RNA species in a single and simple assay. Using the nanopore long-read sequencing platform, Westergren and Colleagues revealed that the complex human adenovirus type 2 (Ad2) transcriptome has a flexible splicing machinery that generates several mRNAs from the early and late transcription units (Westergren Jakobsson et al., 2020). More than 900 alternatively spliced mRNAs generated from the Ad2 transcriptome were identified among which 850 were novel mRNAs. The nanopore sequencing provides access to genomic regions that cannot be accessed by traditional sequencing methods. Thus, application of the nanopore technology for detection and analyses of plant virus can reveal structural variations and can promote early diagnosis, treatment, and monitoring of plant viruses.
Nanopore sequencing technology provides a very convenient method for investigating the role of epigenetics in the regulation of plant virus activities. Epigenetic covalent nucleotide modifications including 5-methylcytosine (5-mC), N4-methylcytosine (4-mC), and N6-methyladenine (6-mA), which do not alter the primary DNA/RNA sequence, participate in the regulation of virus replication. Long-read nanopore sequencing enables direct observation of modified nucleotides by assessing deviating current signals as has been revealed by findings from animal viruses studies. Goldsmith et al. (2021) proposed a nanopore sequencing method for detection of 5mCpG on the HBV genome, which does not rely on bisulfite conversion or PCR. They found that assessment of 5mCpG levels in HBV by bisulfite-quantitative methyl-specific qPCR and nanopore sequencing were highly correlated in their ~2000 times coverage of the viral genome. Zhang et al. (2021) conducted methylated RNA immunoprecipitation sequencing and nanopore direct RNA sequencing analyses to show that SARS-CoV-2 RNA contained m6A modifications. Moreover, the results indicated that SARS-CoV-2 infection not only increased the expression of methyl transferase-like 3 (METTL3) but also altered the distribution of the protein, and these modifications of METTL3 expression affected virus replication. Collectively, the findings described above demonstrate that nanopore technology can effectively detect nucleotide modifications during viral infection.
Currently, the biggest obstacle to application of nanopore sequencing technology is its comparatively lower read accuracy when compared to short read sequencers. Nanopore sequencing platforms have sequencing accuracy rates between 97 and 99%, depending on the flow cell chemistry and the mode of basecalling step. In recent years, ONT has improved the read quality of nanopore sequencing by changing the chemical reagents used in library preparation kits and flow cells and developing improved algorithms for basecalling (Javaran et al., 2021). For example, higher raw read accuracy (98.3%) was achieved by a new basecaller software, Bonito (Silvestre-Ryan and Holmes, 2021). Additional improvements include sequencing kits based on a new chemistry, Q20+, which has been tested in nanopore flow cells and will be released in the near future. This technology uses a refined motor protein (E8.1), which increases the raw read accuracy to 99.3%.4 Bioinformatic analysis has significant impact to the accuracy of Nanopore sequencing. Different computational pipelines of the same nanopore data may lead to different results. Normally, MinION pipeline contains primer trimming, alignment, variant calling and consensus generation (Yang et al., 2020). For example, there are multiple updates to the Guppy algorithm in recent years. At the London Calling meeting of 2020, the median single read accuracy for Guppy 3.6.0 when sequencing a mixture of reference microbial genomes or the human genome was reported at 96.5%. While the current version of Guppy (6.1.1) released in March 2022 increased the sequencing accuracy to 99.2%.5 If the current raw reading accuracy is as claimed by Oxford Nanopore Technologies, we think the remaining errors do not have a great impact on virus identification and this technology is still an interesting option for virus detection from a variety of biological samples.
Another major hindrance to the application of nanopore sequencing is the lack of a bioinformatics system that is fast, accurate, and easy to use. Although many bioinformatics tools have been developed for nanopore sequencing, their adaptability, accuracy, robustness, and efficiency are far from being satisfactory. In typical use cases, bioinformatics software requires users who are proficient in using commands and codes to perform analysis operations in the Linux environment. A webserver is the best option for non-expert users, but this requires strong computer hardware at the remote site and is associated with difficulties in uploading large files.
Taken together, these results indicate that the nanopore sequencing technologies have high diagnostic efficiency for plant viruses and can enrich pathogenic resources in plant disease databases. In the future, these technologies are expected to provide accurate, rapid, and on-site diagnoses for numerous phytosanitary requirements.
Author contributions
YL and KS: writing—original draft preparation. CY and PZ: writing—review and editing. XY: project administration. KS and XS: funding acquisition. All authors contributed to the article and approved the submitted version.
Funding
This research was supported by Key Research Program of Zhejiang Province of China (No. 2022C02047); National Natural Science Foundation of China (No. 32102165); Zhejiang Provincial Natural Science Foundation (No. LQ22C140003); and Zhejiang Province Natural Science Foundation of China (No. LGC21C140001).
Acknowledgments
We thank Andrew O. Jackson for his critical reading of the manuscript.
Conflict of interest
XZ was employed by the company Ausper Biopharma. LM was employed by the company Hangzhou Baiyi Technology Co., Ltd.
The remaining 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.
Footnotes
1. ^https://github.com/nanoporetech
2. ^ftp://ftp.ncbi.nlm.nih.gov/refseq/release/viral
3. ^ftp://ftp.ncbi.nlm.nih.gov/genomes/refseq/
4. ^https://nanoporetech.com/q20plus-chemistry
5. ^https://community.nanoporetech.com/posts/guppy-v6-1-1-release
References
Amoia, S. S., Minafra, A., Nicoloso, V., Loconsole, G., and Chiumenti, M. (2022). A new jasmine virus C isolate identified by Nanopore sequencing is associated to yellow mosaic symptoms of Jasminum officinale in Italy. Plants 11:309. doi: 10.3390/plants11030309
Ben Chehida, S., Filloux, D., Fernandez, E., Moubset, O., Hoareau, M., Julian, C., et al. (2021). Nanopore sequencing is a credible alternative to recover complete genomes of Geminiviruses. Microorganisms. 9:903. doi: 10.3390/microorganisms9050903
Bi, C., Ramos-Mandujano, G., Tian, Y., Hala, S., Xu, J., Mfarrej, S., et al. (2021). Simultaneous detection and mutation surveillance of SARS-CoV-2 and multiple respiratory viruses by rapid field-deployable sequencing. Med 2, 689–700.e4. doi: 10.1016/j.medj.2021.03.015
Boykin, L. M., Sseruwagi, P., Alicai, T., Ateka, E., Mohammed, I. U., Stanton, J. L., et al. (2019). Tree lab: portable genomics for early detection of plant viruses and pests in sub-Saharan Africa. Genes 10:632. doi: 10.3390/genes10090632
Bronzato Badial, A., Sherman, D., Stone, A., Gopakumar, A., Wilson, V., Schneider, W., et al. (2018). Nanopore sequencing as a surveillance tool for plant pathogens in plant and insect tissues. Plant Dis. 102, 1648–1652. doi: 10.1094/PDIS-04-17-0488-RE
Buchfink, B., Xie, C., and Huson, D. H. (2015). Fast and sensitive protein alignment using DIAMOND. Nat. Methods 12, 59–60. doi: 10.1038/nmeth.3176
Cai, D., and Sun, Y. (2022). Reconstructing viral haplotypes using long reads. Bioinformatics 38, 2127–2134. doi: 10.1093/bioinformatics/btac089
Camacho, C., Coulouris, G., Avagyan, V., Ma, N., Papadopoulos, J., Bealer, K., et al. (2009). BLAST+: architecture and applications. BMC 10:421. doi: 10.1186/1471-2105-10-421
Chalupowicz, L., Dombrovsky, A., Gaba, V., Luria, N., Reuven, M., Beerman, A., et al. (2019). Diagnosis of plant diseases using the Nanopore sequencing platform. Plant Pathol. 68, 229–238. doi: 10.1111/ppa.12957
Ciuffreda, L., Rodríguez-Pérez, H., and Flores, C. (2021). Nanopore sequencing and its application to the study of microbial communities. Comput. Struct. Biotechnol. J. 19, 1497–1511. doi: 10.1016/j.csbj.2021.02.020
Claverie, S., Bernardo, P., Kraberger, S., Hartnady, P., Lefeuvre, P., Lett, J. M., et al. (2018). From spatial Metagenomics to molecular characterization of plant viruses: A Geminivirus case study. Adv. Virus Res. 101, 55–83. doi: 10.1016/bs.aivir.2018.02.003
De Coster, W., D'Hert, S., Schultz, D. T., Cruts, M., and Van Broeckhoven, C. (2018). NanoPack: visualizing and processing long-read sequencing data. Bioinformatics 34, 2666–2669. doi: 10.1093/bioinformatics/bty149
Deamer, D., Akeson, M., and Branton, D. (2016). Three decades of nanopore sequencing. Nat. Biotechnol. 34, 518–524. doi: 10.1038/nbt.3423
Della Bartola, M., Byrne, S., and Mullins, E. (2020). Characterization of potato virus Y isolates and assessment of Nanopore sequencing to detect and genotype potato viruses. Viruses 12:478. doi: 10.3390/v12040478
Elena, S. F., Fraile, A., and García-Arenal, F. (2014). Evolution and emergence of plant viruses. Adv. Virus Res. 88, 161–191. doi: 10.1016/B978-0-12-800098-4.00003-9
Fellers, J. P., Webb, C., Fellers, M. C., Shoup Rupp, J., and De Wolf, E. (2019). Wheat virus identification Within infected tissue using Nanopore sequencing technology. Plant Dis. 103, 2199–2203. doi: 10.1094/PDIS-09-18-1700-RE
Filloux, D., Fernandez, E., Loire, E., Claude, L., Galzi, S., Candresse, T., et al. (2018). Nanopore-based detection and characterization of yam viruses. Sci. Rep. 8:17879. doi: 10.1038/s41598-018-36042-7
Fox, E. J., Reid-Bayliss, K. S., Emond, M. J., and Loeb, L. A. (2014). Accuracy of next generation sequencing platforms. Next Gener. Seq. Appl. 1:1000106. doi: 10.4172/jngsa.1000106
Gilbertson, R. L., Batuman, O., Webster, C. G., and Adkins, S. (2015). Role of the insect Supervectors Bemisia tabaci and Frankliniella occidentalis in the emergence and global spread of plant viruses. Annu Rev Virol. 2, 67–93. doi: 10.1146/annurev-virology-031413-085410
Goldberg, B., Sichtig, H., Geyer, C., Ledeboer, N., and Weinstock, G. M. (2015). Making the leap from research laboratory to clinic: challenges and opportunities for next-generation sequencing in infectious disease diagnostics. MBio 6, e01888–e01815. doi: 10.1128/mBio.01888-15
Goldsmith, C., Cohen, D., Dubois, A., Martinez, M. G., Petitjean, K., Corlu, A., et al. (2021). Cas 9-targeted nanopore sequencing reveals epigenetic heterogeneity after de novo assembly of native full-length hepatitis B virus genomes. Microb Genom. 7:507. doi: 10.1099/mgen.0.000507
Grubaugh, N. D., Gangavarapu, K., Quick, J., Matteson, N. L., De Jesus, J. G., Main, B. J., et al. (2019). An amplicon-based sequencing framework for accurately measuring intrahost virus diversity using primal Seq and iVar. Genome Biol. 20:8. doi: 10.1186/s13059-018-1618-7
Gu, W., Miller, S., and Chiu, C. Y. (2019). Clinical metagenomic next-generation sequencing for pathogen detection. Annu. Rev. Pathol. 14, 319–338. doi: 10.1146/annurev-pathmechdis-012418-012751
Hadidi, A., Flores, R., Candresse, T., and Barba, M. (2016). Next-generation sequencing and genome editing in plant virology. Front. Microbiol. 7:1325. doi: 10.3389/fmicb.2016.01325
Imai, K., Tarumoto, N., Misawa, K., Runtuwene, L. R., Sakai, J., Hayashida, K., et al. (2017). A novel diagnostic method for malaria using loop-mediated isothermal amplification (LAMP) and MinION™ nanopore sequencer. BMC Infect. Dis. 17:621. doi: 10.1186/s12879-017-2718-9
Javaran, V. J., Moffett, P., Lemoyne, P., Xu, D., Adkar-Purushothama, C. R., and Fall, M. L. (2021). Grapevine virology in the third-generation sequencing era: From virus detection to viral Epitranscriptomics. Plan. Theory 10:2355. doi: 10.3390/plants10112355
Jaworski, E., Langsjoen, R. M., Mitchell, B., Judy, B., Newman, P., Plante, J. A., et al. (2021). Tiled-ClickSeq for targeted sequencing of complete coronavirus genomes with simultaneous capture of RNA recombination and minority variants. Elife 10:828. doi: 10.1101/2021.03.10.434828
Kearse, M., Moir, R., Wilson, A., Stones-Havas, S., Cheung, M., Sturrock, S., et al. (2012). Geneious basic: an integrated and extendable desktop software platform for the organization and analysis of sequence data. Bioinformatics 28, 1647–1649. doi: 10.1093/bioinformatics/bts199
Kim, D., Lee, J. Y., Yang, J. S., Kim, J. W., Kim, V. N., and Chang, H. (2020). The architecture of SARS-CoV-2 Transcriptome. Cell 181, 914–921.e10. doi: 10.1016/j.cell.2020.04.011
Koren, S., Walenz, B. P., Berlin, K., Miller, J. R., Bergman, N. H., and Phillippy, A. M. (2017). Canu: scalable and accurate long-read assembly via adaptive k-mer weighting and repeat separation. Genome Res. 27, 722–736. doi: 10.1101/gr.215087.116
Kovaka, S., Fan, Y., Ni, B., Timp, W., and Schatz, M. C. (2021). Targeted nanopore sequencing by real-time mapping of raw electrical signal with UNCALLED. Nat. Biotechnol. 39, 431–441. doi: 10.1038/s41587-020-0731-9
Kumar, K. R., Cowley, M. J., and Davis, R. L. (2019). Next-generation sequencing and emerging technologies. Semin. Thromb. Hemost. 45, 661–673. doi: 10.1055/s-0039-1688446
Laver, T., Harrison, J., O'Neill, P. A., Moore, K., Farbos, A., Paszkiewicz, K., et al. (2015). Assessing the performance of the Oxford Nanopore technologies MinION. Biomol Detect Quantif. 3, 1–8. doi: 10.1016/j.bdq.2015.02.001
Lee, S., Nguyen, L. T., Hayes, B. J., and Ross, E. (2021). Prowler: a novel trimming algorithm for Oxford Nanopore sequence data. Bioinformatics 37, 3936–3937. doi: 10.1093/bioinformatics/btab630
Leiva, A. M., Jimenez, J., Sandoval, H., Perez, S., and Cuellar, W. J. (2022). Complete genome sequence of a novel secovirid infecting cassava in the Americas. Arch. Virol. 167, 665–668. doi: 10.1007/s00705-021-05325-2
Leiva, A. M., Siriwan, W., Lopez-Alvarez, D., Barrantes, I., Hemniam, N., Saokham, K., et al. (2020). Nanopore-based complete genome sequence of a Sri Lankan cassava mosaic virus (Geminivirus) strain from Thailand. Microbiol. Resour. Announc. 9:19. doi: 10.1128/MRA.01274-19
Li, H. (2018). Minimap2: pairwise alignment for nucleotide sequences. Bioinformatics 34, 3094–3100. doi: 10.1093/bioinformatics/bty191
Li, H., Handsaker, B., Wysoker, A., Fennell, T., Ruan, J., Homer, N., et al. (2009). The sequence alignment/map format and SAMtools. Bioinformatics 25, 2078–2079. doi: 10.1093/bioinformatics/btp352
Liefting, L. W., Waite, D. W., and Thompson, J. R. (2021). Application of Oxford Nanopore technology to plant virus detection. Viruses 13:424. doi: 10.3390/v13081424
Lin, B., Hui, J., and Mao, H. (2021). Nanopore technology and its applications in gene sequencing. Biosensors 11:214. doi: 10.3390/bios11070214
Lu, H., Giordano, F., and Ning, Z. (2016). Oxford Nanopore MinION sequencing and genome assembly. Genom. Proteom. Bioinform. 14, 265–279. doi: 10.1016/j.gpb.2016.05.004
Maree, H. J., Fox, A., Al Rwahnih, M., Boonham, N., and Candresse, T. (2018). Application of HTS for routine plant virus diagnostics: state of the art and challenges. Front. Plant Sci. 9:1082. doi: 10.3389/fpls.2018.01082
Martins, T. P., Souza, T. A., da Silva, P. S., Nakasu, E. Y. T., Melo, F. L., Inoue-Nagata, A. K., et al. (2021). Nanopore sequencing of tomato mottle leaf distortion virus, a new bipartite begomovirus infecting tomato in Brazil. Arch. Virol. 166, 3217–3220. doi: 10.1007/s00705-021-05220-w
Mehetre, G. T., Leo, V. V., Singh, G., Sorokan, A., Maksimov, I., Yadav, M. K., et al. (2021). Current developments and challenges in plant viral diagnostics: A systematic review. Viruses 13:412. doi: 10.3390/v13030412
Naito, F. Y. B., Melo, F. L., Fonseca, M. E. N., Santos, C. A. F., Chanes, C. R., Ribeiro, B. M., et al. (2019). Nanopore sequencing of a novel bipartite New World begomovirus infecting cowpea. Arch. Virol. 164, 1907–1910. doi: 10.1007/s00705-019-04254-5
Palanga, E., Filloux, D., Martin, D. P., Fernandez, E., Gargani, D., Ferdinand, R., et al. (2016). Metagenomic-based screening and molecular characterization of cowpea-infecting viruses in Burkina Faso. PLoS One 11:e0165188. doi: 10.1371/journal.pone.0165188
Payne, A., Holmes, N., Rakyan, V., and Loose, M. (2019). Bulk Vis: a graphical viewer for Oxford nanopore bulk FAST5 files. Bioinformatics 35, 2193–2198. doi: 10.1093/bioinformatics/bty841
Petersen, L. M., Martin, I. W., Moschetti, W. E., Kershaw, C. M., and Tsongalis, G. J. (2019). Third-generation sequencing in the clinical laboratory: exploring the advantages and challenges of Nanopore sequencing. J. Clin. Microbiol. 58:19. doi: 10.1128/JCM.01315-19
Quick, J., Loman, N. J., Duraffour, S., Simpson, J. T., Severi, E., Cowley, L., et al. (2016). Real-time, portable genome sequencing for Ebola surveillance. Nature 530, 228–232. doi: 10.1038/nature16996
Roossinck, M. J. (2017). Deep sequencing for discovery and evolutionary analysis of plant viruses. Virus Res. 239, 82–86. doi: 10.1016/j.virusres.2016.11.019
Scholthof, K. B., Adkins, S., Czosnek, H., Palukaitis, P., Jacquot, E., Hohn, T., et al. (2011). Top 10 plant viruses in molecular plant pathology. Mol. Plant Pathol. 12, 938–954. doi: 10.1111/j.1364-3703.2011.00752.x
Silvestre-Ryan, J., and Holmes, I. (2021). Pair consensus decoding improves accuracy of neural network basecallers for nanopore sequencing. Genome Biol. 22:38. doi: 10.1186/s13059-020-02255-1
Sović, I., Šikić, M., Wilm, A., Fenlon, S. N., Chen, S., and Nagarajan, N. (2016). Fast and sensitive mapping of nanopore sequencing reads with GraphMap. Nat. Commun. 7:11307. doi: 10.1038/ncomms11307
Walker, B. J., Abeel, T., Shea, T., Priest, M., Abouelliel, A., Sakthikumar, S., et al. (2014). Pilon: an integrated tool for comprehensive microbial variant detection and genome assembly improvement. PLoS One 9:e112963. doi: 10.1371/journal.pone.0112963
Wang, Y., Zhao, Y., Bollas, A., Wang, Y., and Au, K. F. (2021). Nanopore sequencing technology, bioinformatics and applications. Nat. Biotechnol. 39, 1348–1365. doi: 10.1038/s41587-021-01108-x
Watson, M., Thomson, M., Risse, J., Talbot, R., Santoyo-Lopez, J., Gharbi, K., et al. (2015). poRe: an R package for the visualization and analysis of nanopore sequencing data. Bioinformatics 31, 114–115. doi: 10.1093/bioinformatics/btu590
Westergren Jakobsson, A., Segerman, B., Wallerman, O., Lind, S. B., Zhao, H., Rubin, C. J., et al. (2020). The human adenovirus type 2 Transcriptome: An amazing complexity of alternatively spliced mRNAs. J. Virol. 95:20. doi: 10.1128/JVI.01869-20
Xu, L., and Seki, M. (2020). Recent advances in the detection of base modifications using the Nanopore sequencer. J. Hum. Genet. 65, 25–33. doi: 10.1038/s10038-019-0679-0
Yang, M., Cousineau, A., Liu, X., Luo, Y., Sun, D., Li, S., et al. (2020). Direct Metatranscriptome RNA-seq and multiplex RT-PCR amplicon sequencing on Nanopore MinION-promising strategies for multiplex identification of viable pathogens in food. Front. Microbiol. 11:514. doi: 10.3389/fmicb.2020.00514
Keywords: nanopore sequencing, NGS, virus detection, plant pathogens, virus quarantine
Citation: Sun K, Liu Y, Zhou X, Yin C, Zhang P, Yang Q, Mao L, Shentu X and Yu X (2022) Nanopore sequencing technology and its application in plant virus diagnostics. Front. Microbiol. 13:939666. doi: 10.3389/fmicb.2022.939666
Edited by:
Mengji Cao, Southwest University, ChinaReviewed by:
Yanni Sun, City University of Hong Kong, Hong Kong SAR, ChinaPiero Attilio Bianco, University of Milan, Italy
Copyright © 2022 Sun, Liu, Zhou, Yin, Zhang, Yang, Mao, Shentu and Yu. 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: Xuping Shentu, c3R4cEBjamx1LmVkdS5jbg==; Xiaoping Yu, eXhwQGNqbHUuZWR1LmNu
†These authors have contributed equally to this work