Screening metagenomic data for viruses using the e-probe diagnostic nucleic acid assay

A Step Towards Validation of High-Throughput Sequencing for the Identification of Plant Pathogenic Oomycetes

The advancement in high-throughput sequencing (HTS) technology allows the detection of pathogens without the need for isolation or template amplification. Plant regulatory agencies worldwide are adopting HTS as a prescreening tool for plant pathogens in imported plant germplasm. The technique is a multipronged process and, often, the bioinformatic analysis complicates detection. Previously, we developed E-probe diagnostic nucleic acid analysis (EDNA), a bioinformatic tool that detects pathogens in HTS data. EDNA uses custom databases of signature nucleic acid sequences (e-probes) to reduce computational effort and subjectivity when determining pathogen presence in a sample. E-probes of Pythium ultimum and Phytophthora ramorum were previously validated only using simulated HTS data. However, HTS samples generated from infected hosts or pure culture may vary in pathogen concentration, sequencing bias, and data quality, suggesting that each pathosystem requires further validation. Here, we used metagenomic and genomic HTS data generated from infected hosts and pure culture, respectively, to further validate and curate e-probes of Pythium ultimum and Phytophthora ramorum. E-probe length was found to be a determinant of diagnostic sensitivity and specificity; 80-nucleotide e-probes increased the diagnostic specificity to 100%. Curating e-probes to increase specificity affected diagnostic sensitivity only for 80-nucleotide Pythium ultimum e-probes. Comparing e-probes with alternative databases and bioinformatic tools in their speed and ability to find Pythium ultimum and Phytophthora ramorum demonstrated that, although pathogen sequence reads were detected by other methods, they were less specific and slower when compared with e-probes.

An In Silico Detection of a Citrus Viroid from Raw High-Throughput Sequencing Data

E-probe Diagnostic for Nucleic Acid Analysis (EDNA) is a user-friendly bioinformatic tool that has been adapted for the detection and identification of citrus exocortis viroid (CEVd). Here, we describe the procedures for RNA extraction from citrus tissues, library and sequencing preparation, and the utilization of EDNA Mi-Finder online platform on raw high-throughput sequencing (HTS) data.

Towards the validation of high-throughput sequencing (HTS) for routine plant virus diagnostics: measurement of variation linked to HTS detection of citrus viruses and viroids

BACKGROUND

High-throughput sequencing (HTS) has been applied successfully for virus and viroid discovery in many agricultural crops leading to the current drive to apply this technology in routine pathogen detection. The validation of HTS-based pathogen detection is therefore paramount.

METHODS

Plant infections were established by graft inoculating a suite of viruses and viroids from established sources for further study. Four plants (one healthy plant and three infected) were sampled in triplicate and total RNA was extracted using two different methods (CTAB extraction protocol and the Zymo Research Quick-RNA Plant Miniprep Kit) and sent for Illumina HTS. One replicate sample of each plant for each RNA extraction method was also sent for HTS on an Ion Torrent platform. The data were evaluated for biological and technical variation focussing on RNA extraction method, platform used and bioinformatic analysis.

RESULTS

The study evaluated the influence of different HTS protocols on the sensitivity, specificity and repeatability of HTS as a detection tool. Both extraction methods and sequencing platforms resulted in significant differences between the data sets. Using a de novo assembly approach, complemented with read mapping, the Illumina data allowed a greater proportion of the expected pathogen scaffolds to be inferred, and an accurate virome profile was constructed. The complete virome profile was also constructed using the Ion Torrent data but analyses showed that more sequencing depth is required to be comparative to the Illumina protocol and produce consistent results. The CTAB extraction protocol lowered the proportion of viroid sequences recovered with HTS, and the Zymo Research kit resulted in more variation in the read counts obtained per pathogen sequence. The expression profiles of reference genes were also investigated to assess the suitability of these genes as internal controls to allow for the comparison between samples across different protocols.

CONCLUSIONS

This study highlights the need to measure the level of variation that can arise from the different variables of an HTS protocol, from sample preparation to data analysis. HTS is more comprehensive than any assay previously used, but with the necessary validations and standard operating procedures, the implementation of HTS as part of routine pathogen screening practices is possible.

High Throughput Sequencing For Plant Virus Detection and Discovery

Over the last decade, virologists have discovered an unprecedented number of viruses using high throughput sequencing (HTS), which led to the advancement of our knowledge on the diversity of viruses in nature, particularly unraveling the virome of many agricultural crops. However, these new virus discoveries have often widened the gaps in our understanding of virus biology; the forefront of which is the actual role of a new virus in disease, if any. Yet, when used critically in etiological studies, HTS is a powerful tool to establish disease causality between the virus and its host. Conversely, with globalization, movement of plant material is increasingly more common and often a point of dispute between countries. HTS could potentially resolve these issues given its capacity to detect and discover. Although many pipelines are available for plant virus discovery, all share a common backbone. A description of the process of plant virus detection and discovery from HTS data are presented, providing a summary of the different pipelines available for scientists' utility in their research.

Inferring the presence of aflatoxin-producing Aspergillus flavus strains using RNA sequencing and electronic probes as a transcriptomic screening tool

E-probe Diagnostic for Nucleic acid Analysis (EDNA) is a bioinformatic tool originally developed to detect plant pathogens in metagenomic databases. However, enhancements made to EDNA increased its capacity to conduct hypothesis directed detection of specific gene targets present in transcriptomic databases. To target specific pathogenicity factors used by the pathogen to infect its host or other targets of interest, e-probes need to be developed for transcripts related to that function. In this study, EDNA transcriptomics (EDNAtran) was developed to detect the expression of genes related to aflatoxin production at the transcriptomic level. E-probes were designed from genes up-regulated during A. flavus aflatoxin production. EDNAtran detected gene transcripts related to aflatoxin production in a transcriptomic database from corn, where aflatoxin was produced. The results were significantly different from e-probes being used in the transcriptomic database where aflatoxin was not produced (atoxigenic AF36 strain and toxigenic AF70 in Potato Dextrose Broth).

Recent Advances on the Multiplex Molecular Detection of Plant Viruses and Viroids

Plant viruses are still one of the main contributors to economic losses in agriculture. It has been estimated that plant viruses can cause as much as 50 billion euros loss worldwide, per year. This situation may be worsened by recent climate change events and the associated changes in disease epidemiology. Reliable and early detection methods are still one of the main and most effective actions to develop control strategies for plant viral diseases. During the last years, considerable progress has been made to develop tools with high specificity and low detection limits for use in the detection of these plant pathogens. Time and cost reductions have been some of the main objectives pursued during the last few years as these increase their feasibility for routine use. Among other strategies, these objectives can be achieved by the simultaneous detection and (or) identification of several viruses in a single assay. Nucleic acid-based detection techniques are especially suitable for this purpose. Polyvalent detection has allowed the detection of multiple plant viruses at the genus level. Multiplexing RT polymerase chain reaction (PCR) has been optimized for the simultaneous detection of more than 10 plant viruses/viroids. In this short review, we provide an update on the progress made during the last decade on techniques such as multiplex PCR, polyvalent PCR, non-isotopic molecular hybridization techniques, real-time PCR, and array technologies to allow simultaneous detection of multiple plant viruses. Also, the potential and benefits of the powerful new technique of deep sequencing/next-generation sequencing are described.

Nanopore Sequencing as a Surveillance Tool for Plant Pathogens in Plant and Insect Tissues

Plant pathogens are constantly emerging and spreading into new areas and there are often limited postdiagnosis treatment options for infection, making surveillance key to their control. Here we present results from a study testing the efficacy of a portable nanopore-based massively parallel sequencing (MPS) technology for use in the detection of diverse plant pathogens in selected samples. The Oxford MinION device was coupled with whole transcriptome amplification (WTA) to sequence the metatranscriptome of plant and insect tissues infected with either Candidatus Liberibacter asiaticus or plum pox virus. Results showed that this methodology is useful for detecting unsuspected viral and bacterial pathogens in plant and insect tissues. The percentage of generated reads assigned to plum pox virus was 95% from infected tissue and 3% from the viruliferous insect, Myzus persicae. Diaphorina citri sequencing led to 22% of the reads mapping as Ca. L. asiaticus. Plum pox virus and Ca. L. asiaticus were detected in both tissue and insect samples near the beginning of each sequencing run, demonstrating the capability of this methodology to obtain results rapidly. This approach also proved the capability of this system to determine the major components of the insect vector's microbiome and the specific strain of small-genome, high-titer pathogens.

Combining Analysis of DNA in a Crude Virion Extraction with the Analysis of RNA from Infected Leaves to Discover New Virus Genomes

This metagenome approach is used to identify plant viruses with circular DNA genomes and their transcripts. Often plant DNA viruses that occur in low titers in their host or cannot be mechanically inoculated to another host are difficult to propagate to achieve a greater titer of infectious material. Infected leaves are ground in a mild buffer with optimal pH and ionic composition recommended for purifying most bacilliform Para retroviruses. Urea is used to break up inclusion bodies that trap virions and to dissolve cellular components. Differential centrifugation provides further separation of virions from plant contaminants. Then proteinase K treatment removes the capsids. Then the viral DNA is concentrated and used for next-generation sequencing (NGS). The NGS data are used to assemble contigs which are submitted to NCBI-BLASTn to identify a subset of virus sequences in the generated dataset. In a parallel pipeline, RNA is isolated from infected leaves using a standard column-based RNA extraction method. Then ribosome depletion is carried out to enrich for a subset of mRNA and virus transcripts. Assembled sequences derived from RNA sequencing (RNA-seq) were submitted to NCBI-BLASTn to identify a subset of virus sequences in this dataset. In our study, we identified two related full-length badnavirus genomes in the two datasets. This method is preferred to another common approach which extracts the aggregate population of small RNA sequences to reconstitute plant virus genomic sequences. This latter metagenomic pipeline recovers virus related sequences that are retro-transcribing elements inserted into the plant genome. This is coupled to biochemical or molecular assays to further discern the actively infectious agents. The approach documented in this study, recovers sequences representative of replicating viruses that likely indicate active virus infection.

In Silico Probe-Based Detection of Citrus Viruses in NGS Data

The conservation of plant biosecurity relies on the rapid identification of pathogenic organisms, including viruses. With next-generation sequencing (NGS), it is possible to identify multiple viruses within a metagenomic sample. In this study, we explored the use of electronic probes (e-probes) for the simultaneous detection of 11 recognized citrus viruses. E-probes were designed and screened against raw sequencing data to minimize the bioinformatic processing time required. The e-probes were able to accurately detect their cognate viruses in simulated datasets, without any false negatives or positives. The efficiency of the e-probe-based approach was validated with NGS datasets generated from different RNA preparations: double-stranded RNA (dsRNA) from 'Mexican' lime infected with different Citrus tristeza virus (CTV) genotypes, dsRNA from field samples, and small RNA and total RNA from grapefruit infected with the CTV T3 genotype. A set of probes was made available that is able to accurately detect CTV in sequence data regardless of the input dataset or the genotype that plants are infected with.

Detection, discrimination and discovery of a new Tobacco streak virus strain

Soybean plants that exhibited symptoms of virus infection were sampled from different counties of Oklahoma. These plants were tested serologically for 15 viruses known to infect soybean plants. Fifty-seven samples that exhibited typical virus-like symptoms did not test positive for any of the 15 viruses used in a dot-immunobinding assay (DIBA). Four samples were pooled and used for next generation sequencing using the 454-Roche protocol. Sequence and phylogenetic analysis of the sequences obtained revealed infection with a distinct strain of Tobacco streak virus (TSV). TSV was one of the 15 viruses initially tested for using DIBA and had tested negative. TSV belongs to the genus Ilarvirus and has been reported as a causal agent of bud blight in soybean crops in Brazil and the United States. Out of 10 reported primer pairs for TSV reverse transcription-polymerase chain reaction (RT-PCR), only two had the potential, based on sequence similarity, to amplify part of the genome of the distinct strain of TSV found in Oklahoma and only one was actually able to amplify the region. In this study, a new primer pair, specific to all known TSV and capable of amplifying the Oklahoma strain (TSV-OK), was designed from a highly conserved region of coat protein (CP) sequences and end-point PCR and quantitative RT-PCR detection methods were developed and their sensitivity assayed. This is the first report of specific primers designed from this highly conserved region in the CP of TSV for detection of TSV. Twenty-three of the 57 DIBA soybean samples that initially tested negative were retested with the new specific end-point PCR method and found positive for TSV infection.

Validation of high throughput sequencing and microbial forensics applications

High throughput sequencing (HTS) generates large amounts of high quality sequence data for microbial genomics. The value of HTS for microbial forensics is the speed at which evidence can be collected and the power to characterize microbial-related evidence to solve biocrimes and bioterrorist events. As HTS technologies continue to improve, they provide increasingly powerful sets of tools to support the entire field of microbial forensics. Accurate, credible results allow analysis and interpretation, significantly influencing the course and/or focus of an investigation, and can impact the response of the government to an attack having individual, political, economic or military consequences. Interpretation of the results of microbial forensic analyses relies on understanding the performance and limitations of HTS methods, including analytical processes, assays and data interpretation. The utility of HTS must be defined carefully within established operating conditions and tolerances. Validation is essential in the development and implementation of microbial forensics methods used for formulating investigative leads attribution. HTS strategies vary, requiring guiding principles for HTS system validation. Three initial aspects of HTS, irrespective of chemistry, instrumentation or software are: 1) sample preparation, 2) sequencing, and 3) data analysis. Criteria that should be considered for HTS validation for microbial forensics are presented here. Validation should be defined in terms of specific application and the criteria described here comprise a foundation for investigators to establish, validate and implement HTS as a tool in microbial forensics, enhancing public safety and national security.

Metagenomic search strategies for interactions among plants and multiple microbes

Plants harbor multiple microbes. Metagenomics can facilitate understanding of the significance, for the plant, of the microbes, and of the interactions among them. However, current approaches to metagenomic analysis of plants are computationally time consuming. Efforts to speed the discovery process include improvement of computational speed, condensing the sequencing reads into smaller datasets before BLAST searches, simplifying the target database of BLAST searches, and flipping the roles of metagenomic and reference datasets. The latter is exemplified by the e-probe diagnostic nucleic acid analysis approach originally devised for improving analysis during plant quarantine.