Evolutionary Rates and Phylogeographical Analysis of Odontoglossum Ringspot Virus Based on the 166 Coat Protein Gene Sequences

First report of Odontoglossum ringspot virus in Vanilla (Orchidaceae) in Madagascar

Vanilla (Vanilla planifolia, Orchidaceae) is Madagascar's leading agricultural export resource which provides 80% of world's consumption. During a phytosanitary survey conducted from November 2019 to March 2021 in the main vanilla production regions of Madagascar, 250 plots were indexed for cymbidium mosaic virus (CymMV, Potexvirus genus) and odontoglossum ringspot virus (ORSV, Tobamovirus genus) the two most prevalent viruses of cultivated orchids worldwide (Zettler et al., 1990). For each plot, bulk samples (ten leaves taken at random) were assayed using Immunostrips (AGDIA, ISK 13301). A quarter of the plots (63/250) tested positive for CymMV. The highest prevalence of CymMV was observed in the SAVA region (57 out of 153 plots = 37,2%) where the virus has been reported since 1997 (Grisoni et al., 2010). Six plots located in the district of Mahanoro (Atsinanana) tested positive for ORSV. A few plants in these plots showed chlorotic often annular spots on their leaves. They were individually tested positive for ORSV, and negative for CymMV and potyviruses (Immunostrips AGDIA ISK 27200), the other two viruses reported so far in vanilla in Madagascar. To confirm the diagnosis of ORSV, leaf samples from five of the six infected plots were analysed by Tube Capture-RT-PCR (Grisoni et al., 2017) using two pairs of primers flanking the ORSV coat protein (CP) gene: OrCP1 (GGTCGGTAATGGTGTTAG) / OrCP2 (TGCATTATCGTATGCTCC), and CPOR-F(ATGTCTTACACTATTACAGACC) / CPOR-R(TTAGGAAGAGGTCCAAGTAAG). The five samples gave amplicons of the expected size (820 nt and 476 nt, respectively) and were sequenced with Sanger technology (Macrogen, The Netherlands). The ORSV-CP sequences of the Mahanoro isolates showed very close similarity to 198 ORSV-CP sequences from GenBank (95.8% to 99.6% nucleotide and 94.5 to 100% amino-acid identities), and less than 75.4% nucleotide (80.1% amino-acid) identities with Bell pepper mosaic virus (DQ355023), the tobamovirus closest to ORSV. The five ORSV-CP sequences from vanilla were deposited in GenBank under accessions numbers OM847399 to OM847403. These data confirmed that ORSV infects vanilla vines in Madagascar. To our knowledge, this is the first report of this virus in Madagascar and of its ability to infect symptomatically V. planifolia. The five ORSV isolates from vanilla had more than 98.7 % nucleotide identities of CP gene and clustered into a monophyletic group in maximum likelihood phylogenetic tree, suggesting a single origin of these isolates. To further investigate the origin of ORSV in Madagascar, we made use of RNA sequences isolated at different points in time to infer the timing of evolutionary events (Rieux et al., 2016). We estimated the CP gene substitution rate to 4.8E-4 subst/site/year [95%HPD 2.1E-4 - 8.7E-4] which is close to the estimate of He et al. (2019) based on a slightly different sequences set (1.25E-3 subst/site/year). We dated the initial contamination of vanilla plts by ORSV between 2004 and 2013. Both ORSV and CymMV have deleterious effects on many ornamental orchids, and the pathogenicity of CymMV is exacerbated when co-infecting with ORSV (Lee et al., 2021). Therefore, ORSV represents a new threat to the Malagasy vanilla crop, especially in regions where CymMV is already rife. Given the economic importance of vanilla cultivation in the country, the implementation of prophylactic measures aimed at preventing the spread of ORSV, in particular through the sanitary control of cuttings, should be a priority for the vanilla industry.

Phylodynamics and Codon Usage Pattern Analysis of Broad Bean Wilt Virus 2

Broad bean wilt virus 2 (BBWV-2), which belongs to the genus Fabavirus of the family Secoviridae, is an important pathogen that causes damage to broad bean, pepper, yam, spinach and other economically important ornamental and horticultural crops worldwide. Previously, only limited reports have shown the genetic variation of BBWV2. Meanwhile, the detailed evolutionary changes, synonymous codon usage bias and host adaptation of this virus are largely unclear. Here, we performed comprehensive analyses of the phylodynamics, reassortment, composition bias and codon usage pattern of BBWV2 using forty-two complete genome sequences of BBWV-2 isolates together with two other full-length RNA1 sequences and six full-length RNA2 sequences. Both recombination and reassortment had a significant influence on the genomic evolution of BBWV2. Through phylogenetic analysis we detected three and four lineages based on the ORF1 and ORF2 nonrecombinant sequences, respectively. The evolutionary rates of the two BBWV2 ORF coding sequences were 8.895 × 10⁻⁴ and 4.560 × 10⁻⁴ subs/site/year, respectively. We found a relatively conserved and stable genomic composition with a lower codon usage choice in the two BBWV2 protein coding sequences. ENC-plot and neutrality plot analyses showed that natural selection is the key factor shaping the codon usage pattern of BBWV2. Strong correlations between BBWV2 and broad bean and pepper were observed from similarity index (SiD), codon adaptation index (CAI) and relative codon deoptimization index (RCDI) analyses. Our study is the first to evaluate the phylodynamics, codon usage patterns and adaptive evolution of a fabavirus, and our results may be useful for the understanding of the origin of this virus.