Characterization of a novel citrus tristeza virus genotype within three cross-protecting source GFMS12 sub-isolates in South Africa by means of Illumina sequencing

Population dynamics of Citrus tristeza virus (CTV) in single aphid-transmitted sub-isolates of the South African GFMS12 isolate

Grapefruit trees in South Africa have been cross protected against severe stem pitting genotypes of Citrus tristeza virus (CTV) since the 1920s using a mild strain initially called 'Nartia' but later referred to as grapefruit mild strain 12 (GFMS12). In the current study, the GFMS12 isolate was used as the source for single aphid transmissions (SAT) using Toxoptera citricida, commonly called the brown citrus aphid (BrCA). The BrCA-transmitted CTV sub-isolates were analyzed by the heteroduplex mobility assay (HMA), serological assays, genetic marker analysis (GMA), and selected sub-isolates were biologically indexed. Reverse transcription PCR of genomic regions was conducted using universal primers followed by cloning the PCR products, HMA and sequence analysis; nine genotypes of CTV were identified in the complex of GFMS12, including both severe and mild genotypes. A single BrCA transmitted up to six CTV genotypes simultaneously in one sub-isolate. The HMA was found to be a rapid, reliable tool for the identification of genotypes and can be useful in the development of CTV management strategies and budwood certification programs.

High throughput sequencing from Angolan citrus accessions discloses the presence of emerging CTV strains

BACKGROUND

Citrus industry is worldwide dramatically affected by outbreaks of Citrus tristeza virus (CTV). Controls should be applied to nurseries, which could act as diversity hotspots for CTV. Early detection and characterization of dangerous or emerging strains of this virus greatly help to prevent outbreaks of disease. This is particularly relevant in those growing regions where no dedicated certification programs are currently in use.

METHODS

Double-stranded RNA extracted from Citrus spp. samples, collected in two locations in Angola, were pooled and submitted to a random-primed RNA-seq. This technique was performed to acquire a higher amount of data in the survey, before the amplification and sequencing of genes from single plants. To confirm the CTV infection in individual plants, as suggested by RNA-seq information from the pooled samples, the analysis was integrated with multiple molecular marker amplification (MMM) for the main known CTV strains (T30, T36, VT and T3).

RESULTS

From the analysis of HTS data, several assembled contigs were identified as CTV and classified according to their similarity to the established strains. By the MMM amplification, only five individual accessions out of the eleven pooled samples, resulted to be infected by CTV. Amplified coat protein genes from the five positive sources were cloned and sequenced and submitted to phylogenetic analysis, while a near-complete CTV genome was also reconstructed by the fusion of three overlapping contigs.

CONCLUSION

Phylogenetic analysis of the ORF1b and CP genes, retrieved by de novo assembly and RT-PCR, respectively, revealed the presence of a wide array of CTV strains in the surveyed citrus-growing spots in Angola. Importantly, molecular variants among those identified from HTS showed high similarity with known severe strains as well as to recently described and emerging strains in other citrus-growing regions, such as S1 (California) or New Clade (Uruguay).

Citrus Tristeza Virus Isolates of the Same Genotype Differ in Stem Pitting Severity in Grapefruit

Two isolates of the T68 genotype of citrus tristeza virus (CTV) were derived from a common source, GFMS12, by single aphid transmission. These isolates, named GFMS12-8 and GFMS12-1.3, induced stem pitting with differing severity in 'Duncan' grapefruit (Citrus × paradisi [Macfad.]). Full-genome sequencing of these isolates showed only minor nucleotide sequence differences totaling 45 polymorphisms. Numerous nucleotide changes, in relatively close proximity, were detected in the p33 open reading frame (ORF) and the leader protease domains of ORF1a. This is the first report of full-genome characterization of CTV isolates of a single genotype, derived from the same source, but showing differences in pathogenicity. The results demonstrate the development of intragenotype heterogeneity known to occur with single-stranded RNA viruses. Identification of genetic variability between isolates showing different pathogenicity will enable interrogation of specific genome regions for potential stem pitting determinants.

A New Era for Mild Strain Cross-Protection

Societal and environmental pressures demand high-quality and resilient cropping plants and plant-based foods grown with the use of low or no synthetic chemical inputs. Mild strain cross-protection (MSCP), the pre-immunization of a plant using a mild strain of a virus to protect against subsequent infection by a severe strain of the virus, fits with future-proofing of production systems. New examples of MSCP use have occurred recently. New technologies are converging to support the discovery and mechanism(s) of action of MSCP strains thereby accelerating the popularity of their use.

Analysis of Genotype Composition of Citrus tristeza virus Populations Using Illumina Miseq Technology

Recent research describing the strain-specific mechanisms underlying experimental CTV superinfection exclusion has far-reaching implications for the manner in which cross-protecting sources should be selected for. The strain composition of both cross-protecting sources and field populations needs to be sufficiently characterized to improve control of severe stem-pitting and decline isolates. Many of the biological, serological, and molecular techniques used in previous studies yield very limited information about the strain composition of populations and the relative titer of their components. In this chapter we describe a protocol for the characterization of CTV populations, based on the use of the next-generation sequencing Illumina MiSeq platform of p33 gene amplicons.

Bioinformatic Tools and Genome Analysis of Citrus tristeza virus

High-throughput sequencing (HTS) is a powerful tool employed by plant virologists for the detection of viruses, the characterization of virus genomes and the study of host-pathogen interactions. Virus detection has been an important application of this technology, which has resulted in the discovery of novel viruses or viral strains as well as for the detection of known viruses in a plant sample. Here we describe the entire process that needs to be considered for the genome analysis of Citrus tristeza virus (CTV) by HTS, including the experimental design, sample preparation, nucleic acid purification, HTS library construction, and bioinformatic analysis.

Analyses of 3' half genome of citrus tristeza virus reveal existence of distinct virus genotypes in citrus growing regions of India

Citrus tristeza virus (CTV, genus Closterovirus) is one of the most serious pathogens responsible for huge loss of citrus trees worldwide. Four Indian CTV isolates, Kat1 (C. reticulata/Central India), D1 (C. sinensis/North India), B5 (Citrus limettoides/South India) and G28 (C. lemon/Northeast India) collected from different regions of India were characterized based on sequencing of 3' half genome (~ 8.4 kb) comprising 10 open reading frames (ORFs2-11) and 3' UTR and the sequences were submitted to NCBI database as Acc. No KJ914662, HQ912022, HQ912023 and KJ914661, respectively. The present and previously reported Indian isolates Kpg3 and B165 were analyzed and compared with other Asian and international CTV isolates. The Indian CTV isolates had 92-99% nt identities among them. The phylogenetic analysis generated overall ten genogroups/lineages. Of them, all the Asian isolates fell into seven genogroups, whereas the Indian isolates into four. Indian isolates Kat1, D1 and Kpg3 grouped together, termed "Kpg3Gr", along with Florida severe isolate T3. The Indian isolates B5, and G28 were found to be two distinct and separate lineages, indicating that these isolates are two new CTV entities. Based on phylogenetic analysis, Kpg3Gr was identified as "Indian VT" subtype which is distinct from the Asian and the Western VT subtype within diversified VT genotype. The recombination detecting-program, RDP4 detected Indian isolates Kat1, B5, B165 and G28 as recombinants, where G28 as strong recombinant. The present study determined the occurrence of at least four CTV genotypes, B5 (distinct), B165 (T68 type) G28 (distinct) and Kpg3Gr in citrus growing regions of India.

Recent Advances on Detection and Characterization of Fruit Tree Viruses Using High-Throughput Sequencing Technologies

Perennial crops, such as fruit trees, are infected by many viruses, which are transmitted through vegetative propagation and grafting of infected plant material. Some of these pathogens cause severe crop losses and often reduce the productive life of the orchards. Detection and characterization of these agents in fruit trees is challenging, however, during the last years, the wide application of high-throughput sequencing (HTS) technologies has significantly facilitated this task. In this review, we present recent advances in the discovery, detection, and characterization of fruit tree viruses and virus-like agents accomplished by HTS approaches. A high number of new viruses have been described in the last 5 years, some of them exhibiting novel genomic features that have led to the proposal of the creation of new genera, and the revision of the current virus taxonomy status. Interestingly, several of the newly identified viruses belong to virus genera previously unknown to infect fruit tree species (e.g., Fabavirus, Luteovirus) a fact that challenges our perspective of plant viruses in general. Finally, applied methodologies, including the use of different molecules as templates, as well as advantages and disadvantages and future directions of HTS in fruit tree virology are discussed.

Molecular and biological characterization of a novel mild strain of citrus tristeza virus in California

Strain differentiating marker profiles of citrus tristeza virus (CTV) isolates from California have shown the presence of multiple genotypes. To better define the genetic diversity involved, full-length genome sequences from four California CTV isolates were determined by small-interfering RNA sequencing. Phylogenetic analysis and nucleotide sequence comparisons differentiated these isolates into the genotypes VT (CA-VT-AT39), T30 (CA-T30-AT4), and a new strain called S1 (CA-S1-L and CA-S1-L65). S1 isolates had three common recombination events within portions of genes from VT, T36 and RB strains and were transmissible by Aphis gossypii. Virus indexing showed that CA-VT-AT39 could be classified as a severe strain, whereas CA-T30-AT4, CA-S1-L and CA-S1-L65 were mild. CA-VT-AT39, CA-S1-L, and CA-S1-L65 reacted with monoclonal antibody MCA13, whereas CA-T30-AT4 did not. RT-PCR and RT-qPCR detection assays for the S1 strain were developed and used to screen MCA13-reactive isolates in a CTV collection from central California collected from 1968 to 2011. Forty-two isolates were found to contain the S1 strain, alone or in combinations with other genotypes. BLAST and phylogenetic analysis of the S1 p25 gene region with other extant CTV sequences from the NCBI database suggested that putative S1-like isolates might occur elsewhere (e.g., China, South Korea, Turkey, Bosnia and Croatia). This information is important for CTV evolution, detection of specific strains, and cross-protection.

Characterization of Citrus tristeza virus Single-Variant Sources in Grapefruit in Greenhouse and Field Trials

Citrus tristeza virus (CTV) is endemic to southern Africa and the stem pitting syndrome that it causes was a limiting factor in grapefruit production prior to the introduction of cross-protection in the Citrus Improvement Scheme. This disease mitigation strategy, using various field-derived CTV sources, has significantly extended the productive lifespan of grapefruit orchards in South Africa. CTV commonly occurs as a population of various strains, masking the phenotypic effect of individual strains. Likewise, current South African CTV cross-protection sources are strain mixtures, obscuring an understanding of which strains are influencing cross-protection. The severity of various CTV strains has mostly been assessed on sensitive indicator hosts, but their effect on commercial varieties has seldom been investigated. Single-variant CTV isolates were used to investigate the phenotypic expression of CTV strains in commercial grapefruit varieties as well as CTV indicator hosts. They were biologically characterized for their ability to cause stem pitting and their rate of translocation and titer in the different hosts, monitored by enzyme-linked immunosorbent assay. Complete genome sequences for three CTV strain variants were generated. Isolates of CTV strains VT, T68, RB, and HA16-5 did not induce severe stem pitting in four grapefruit hosts in a glasshouse trial. Viral titers of the strains differed in the grapefruit hosts, but the RB isolate reached a higher titer in the grapefruit hosts compared with the VT, T68, and HA16-5 isolates. Additionally, horticultural assessment of two grapefruit varieties inoculated with the RB isolate in two field trials demonstrated that mild stem pitting did not negatively influence the horticultural performance of the grapefruit trees over an eight-year assessment period. 'Star Ruby' trees containing the CTV source GFMS35 showed less stem pitting than trees inoculated with the RB isolate, but had smaller canopy volumes and lower yields than trees containing the RB isolate. This suggests that the influence of CTV sources on tree performance is not limited to the effect of stem pitting.

Comparison of multiple viral population characterization methods on a candidate cross-protection Citrus tristeza virus (CTV) source

Citrus tristeza virus (CTV) is the most economically important virus found on citrus and influences production worldwide. The 3' half of the RNA genome is generally conserved amongst sources, whereas the 5' portion is more divergent, allowing for the classification of the virus into a number of genotypes based on sequence diversity. The acknowledged genotypes of CTV are continually being expanded, and thus far include T36, T30, T3, VT, B165, HA16-5, T68 and RB. The genotype composition of the CTV populations of a potential cross protection source in Mexican lime was studied whilst comparing different techniques of viral population characterization. Cloning and sequencing of an ORF 1a fragment, genotype specific RT-PCRs and Illumina sequencing of the p33 gene as well as RNA template enrichment through immuno-capture was done. Primers used in the cloning and sequencing proved to be biased towards detection of the VT genotype. RT-PCR and Illumina sequencing using the two different templates provided relatively comparable results, even though the immuno-captured enriched template provided less than expected CTV specific data, while the RT-PCRs and p33 sequencing cannot be used to make inferences about the rest of the genome; which may vary due to recombination. The source was found to contain multiple genotypes, including RB and VT. When choosing a characterization method, the features of the virus under study should be considered. It was found that Illumina sequencing offers an opportunity to gain a large amount of information regarding the entire viral genome, but challenges encountered are discussed.

Identification of viruses and viroids by next-generation sequencing and homology-dependent and homology-independent algorithms

A fast, accurate, and full indexing of viruses and viroids in a sample for the inspection and quarantine services and disease management is desirable but was unrealistic until recently. This article reviews the rapid and exciting recent progress in the use of next-generation sequencing (NGS) technologies for the identification of viruses and viroids in plants. A total of four viroids/viroid-like RNAs and 49 new plant RNA and DNA viruses from 18 known or unassigned virus families have been identified from plants since 2009. A comparison of enrichment strategies reveals that full indexing of RNA and DNA viruses as well as viroids in a plant sample at single-nucleotide resolution is made possible by one NGS run of total small RNAs, followed by data mining with homology-dependent and homology-independent computational algorithms. Major challenges in the application of NGS technologies to pathogen discovery are discussed.