Nonstructural p26 proteins encoded by the 3’-proximal genes of velariviruses and criniviruses are orthologs

The 3’-most genes in RNA-2 of the Crinivirus genus members (family Closteroviridae) code for non-structural p26 proteins that share amino acid sequence similarity [Stewart LR, Hwang MS, Falk BW (2009) Virus Res 145:293-299]. In this study, sensitive bioinformatic tools have been used to identify the homologous p26 proteins encoded by the 3’ genes in monopartite genomes of the members of Velarivirus, another Closteroviridae genus, and mint vein banding-associated virus, an unassigned member of the family. The p26 proteins showed similarity in their predicted secondary structures, but an amino acid sequence alignment showed no strictly conserved positions, thus indicating a high plasticity of these non-structural proteins. The implications of the sequence analysis for possible functions of the crinivirus and velarivirus p26 proteins are discussed. Electronic supplementary material The online version of this article (10.1007/s00705-019-04491-8) contains supplementary material, which is available to authorized users.

The family Closteroviridae includes about 50 filamentous plant viruses with large positive-sense RNA genomes that show numerous traces of recombination events, such as gene duplication and gene capture [1][2][3]. Closteroviruses contain up to 12 genes, most of which are arranged in two conserved modules: the replicative module, which encodes proteins responsible for RNA synthesis and membrane modification, and the five-gene block, which encodes proteins involved in particle formation and cell-to-cell movement [2][3][4] (Fig. 1). In addition, closteroviruses carry variable accessory genes in the 3' part of their genome (Fig. 1). Some products of the 3' genes are conserved in some members of the Closteroviridae [3,5], whereas others have no apparent homologs and are species-specific.
The members of the genus Crinivirus have bipartite genomes with RNA-1 and RNA-2 bearing the replicase module and the five-gene block, respectively, plus some additional 3' ORFs [6,7] (Fig. 1). The 3' genes in RNA-2 of criniviruses code for proteins with a molecular weight of about 26 kDa (p26) that have certain similarity in amino acid sequence and predicted secondary structure [3,5] (Fig. 1). In this study, our purpose was to identify the p26-related proteins in members of the other genera and unassigned virus species of the family Closteroviridae.
Initial BLAST and PSI-BLAST searches did not reveal any putative p26 protein homologs outside the genus Crinivirus when the crinivirus p26 sequences were used as a query (Supplementary Table S1). At the next step, we used HMMER 2.0, a sensitive tool for detecting remote protein homologs [10]. The hidden Markov model for the crinivirus p26 proteins was constructed and used for directed search for the 3' ORF products in members of other Closteroviridae genera (Velarivirus, Closterovirus, and Ampelovirus). Low probability values (indicating a statistically significant similarity) were obtained for the 27-to 29-kDa proteins encoded Handling Editor: Jesús Navas-Castillo.

Electronic supplementary material
The online version of this article (https ://doi.org/10.1007/s0070 5-019-04491 -8) contains supplementary material, which is available to authorized users. * A. A. Agranovsky aaa@genebee.msu.su by the 3'-most ORFs of velariviruses (below, also referred to as p26 proteins) ( Table 1). Figure 2 shows a sequence alignment of the p26 proteins of criniviruses and velariviruses. Although some positions in the alignment are occupied by similar amino acid residues, none of them is strictly conserved (Fig. 2). Analysis of secondary structure suggests that the p26 proteins of velariviruses and criniviruses are alpha-helical with a few beta-strands (Fig. 2), which corroborates the previous data for the crinivirus proteins [5].
Six alpha-helices in the p26 proteins have a similar location (Fig. 2), suggesting that these proteins share a common three-dimensional structure.
We also performed an additional HMM database search using the combined multiple alignment of crinivirus and velarivirus p26 proteins (Fig. 2) as a query for HMMER2.0. The 24-kDa protein of mint vein banding-associated virus (MVBaV), an unassigned member of the family Closteroviridae [13], was detected as a possible remote homolog (Fig. 3). A significant probability value (0.009) was obtained, supporting the relatedness of the MVBaV p24 to the p26 of criniviruses and velariviruses (Fig. 3). Attempts to include the MVBaV p24 into the multiple alignment using T-Coffee [9] were not successful due to the lack of detectable similarity in the C-terminal regions (Fig. 3). Additional HHPred database searches did not reveal any putative homologs of crinivirus/velarivirus p26 proteins among the available HMM profiles [11].
The data presented here indicate that the p26 genes are not Crinivirus taxon-specific as has been thought but are conserved across the monopartite genomes of Velarivirus members and mint vein banding-associated virus. Taking into account the 3'-proximal location of the p26 genes   ( Fig. 1), similarity of the predicted secondary structures, and statistically significant similarity of the amino acid sequences (Table 1, Fig. 2), it is likely that the p26 proteins of criniviruses and velariviruses are orthologs that may perform the same or similar function(s). It should be noted that members of the genera Crinivirus and Velarivirus are markedly different from each other in their biological properties and the genome structure. Criniviruses have divided genomes, are transmitted by whiteflies, and infect herbaceous hosts, whereas velariviruses possess monopartite genomes, have no known vectors, and infect woody hosts [2] (Fig. 1). The absence of conserved amino acid positions in the p26 alignment suggests high plasticity of these non-structural proteins, which may indicate their involvement in the response of virus systems to rapidly changing environmental conditions. On the other hand, the p26 protein of lettuce infectious yellows virus (LIYV), the type member of the genus Crinivirus, induces specific ultrastructures in the infected cells -conical plasmalemma deposits over plasmadesmata -that are thought to be associated with the vascular transport of the virus [7,14,15]. In support of this, a knockout LIYV p26 mutant proved to be unable to spread systemically in a Nicotiana benthamiana host [15]. Although the Closteroviridae members other than LIYV do not induce plasmalemma deposits, the involvement of