Next Article in Journal
Mosquito Vector Competence for Japanese Encephalitis Virus
Next Article in Special Issue
Discovery of Known and Novel Viruses in Wild and Cultivated Blueberry in Florida through Viral Metagenomic Approaches
Previous Article in Journal
Emergence of a Novel Dengue Virus 3 (DENV-3) Genotype-I Coincident with Increased DENV-3 Cases in Yangon, Myanmar between 2017 and 2019
Previous Article in Special Issue
Novel Ampeloviruses Infecting Cassava in Central Africa and the South-West Indian Ocean Islands
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Viruses
Volume 13
Issue 6
10.3390/v13061153
viruses-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Diversity and Distribution of Viruses Infecting Wild and Domesticated Phaseolus spp. in the Mesoamerican Center of Domestication

by
Elizabeth Chiquito-Almanza
1,
Juan Caballero-Pérez
2,
Jorge A. Acosta-Gallegos
3,
Victor Montero-Tavera
1,
Luis Antonio Mariscal-Amaro
4 and
José Luis Anaya-López
1,*
1
Biotechnology Department, National Institute for Forestry Agriculture and Livestock Research (INIFAP), Celaya, Guanajuato 38110, Mexico
2
Faculty of Chemistry, Autonomous University of Querétaro, Santiago de Querétaro 76017, Mexico
3
Bean Breeding Program, National Institute for Forestry Agriculture and Livestock Research (INIFAP), Celaya, Guanajuato 38110, Mexico
4
Forestry and Plant Protection Program, National Institute for Forestry Agriculture and Livestock Research (INIFAP), Celaya, Guanajuato 38110, Mexico
*
Author to whom correspondence should be addressed.
Viruses 2021, 13(6), 1153; https://doi.org/10.3390/v13061153
Submission received: 15 April 2021 / Revised: 26 May 2021 / Accepted: 11 June 2021 / Published: 16 June 2021
(This article belongs to the Special Issue Plant Virus Surveillance and Metagenomics)
Browse Figures
Review Reports Versions Notes

Abstract

:
Viruses are an important disease source for beans. In order to evaluate the impact of virus disease on Phaseolus biodiversity, we determined the identity and distribution of viruses infecting wild and domesticated Phaseolus spp. in the Mesoamerican Center of Domestication (MCD) and the western state of Nayarit, Mexico. We used small RNA sequencing and assembly to identify complete or near-complete sequences of forty-seven genomes belonging to nine viral species of five genera, as well as partial sequences of two putative new endornaviruses and five badnavirus- and pararetrovirus-like sequences. The prevalence of viruses in domesticated beans was significantly higher than in wild beans (97% vs. 19%; p < 0.001), and all samples from domesticated beans were positive for at least one virus species. In contrast, no viruses were detected in 80–83% of the samples from wild beans. The Bean common mosaic virus and Bean common mosaic necrosis virus were the most prevalent viruses in wild and domesticated beans. Nevertheless, Cowpea mild mottle virus, transmitted by the whitefly Bemisia tabaci, has the potential to emerge as an important pathogen because it is both seed-borne and a non-persistently transmitted virus. Our results provide insights into the distribution of viruses in cultivated and wild Phaseolus spp. and will be useful for the identification of emerging viruses and the development of strategies for bean viral disease management in a center of diversity.

1. Introduction

Virus diseases are a significant threat to agriculture production. The emergence of indigenous and newly introduced viruses is a potential threat to biodiversity in centers of origin and plant domestication, reviewed in [1]; therefore, studies on the dynamics of pathogens and their hosts at sites where wild and domesticated hosts coexist are extremely valuable. Understanding the diversity and evolution of the prevalent viruses in wild and domesticated Phaseolus spp. is needed to quantify the impact of viruses on bean cultivation and the effect of human management on population heterogeneity, biodiversity loss and disease risk. Common bean (Phaseolus vulgaris L.) is likely the most susceptible plant species in the Leguminosae to virus infection [2] and one of the most widely cultivated pulses in the world [3].
There are two major genetic pools of common bean, the Mesoamerican and Andean, reviewed in [4]. Genomic and phenotypic data support the Mesoamerican origin of the common bean (P. vulgaris L.) [5,6], with the most likely center of origin in Mexico [6], and one of the Mesoamerican centers of domestication (MCD) is located in the Lerma-Santiago basin in the states of Jalisco and Guanajuato, Mexico [7]. In Mexico, there have not been comprehensive studies on the distribution and identification of viral species that infect domesticated and wild beans. This information is required for strategic breeding and use of cultivars, for informed decision-making on imposing quarantine, and for developing diagnostic tools and predicting the emergence of new viral strains [8].
The analysis of RNA sequences from plant samples through next-generation sequencing (NGS) is an effective method for virus detection, reviewed in [9]. Among NGS approaches, small RNA sequencing and assembly (sRSA) offers a feasible alternative for the simultaneous detection of multiple known and novel DNA and RNA viruses without prior knowledge or the need for antibodies or virus-specific primers or probes and is 10-times more sensitive than reverse-transcription real-time PCR (RT-qPCR) [10,11]. The sRSA is based on RNA interference (RNAi) of the plant antiviral defense system, which targets double-stranded RNA (dsRNA). Since DNA viruses also induce the production of small interfering RNA (siRNA) in plants, the small RNA (sRNA) from infected plant tissues can be extracted and sequenced, and the sequences of virus-derived from small interfering RNAs (vsiRNA) can then be reassembled in silico into partial or full viral genomes of both DNA and RNA viruses, which allows their identification by comparison with viral sequences available in databases [10].
In aims to expand the knowledge on the diversity, distribution, and dissemination of viruses infecting domesticated and wild beans in the MCD and surrounding areas, symptomatic and asymptomatic samples from field-grown common bean and wild bean populations were collected in Guanajuato, Jalisco, and Nayarit, Mexico, and were analyzed through the sRSA approach using the Illumina platform. The data obtained revealed the presence of nine known viral species, two putative novel alphaendornaviruses, and five virus-like sequences of putative badnaviruses and pararetroviruses, and significant differences in the prevalence of viruses between domesticated and wild Phaseolus spp. Our results provide insights into the distribution of viruses in domesticated and wild Phaseolus spp., as well as information that could be useful for the identification of emerging viruses and the development of strategies for common bean viral disease management.

2. Materials and Methods

2.1. Plant Material

Foliar tissues from 627 plants of wild (Phaseolus spp.) and domesticated bean (P. vulgaris L.) with and without symptoms of viral infection were collected in the MCD and the western State of Nayarit. The foliar tissues were collected and preserved in silica gel with a moisture indicator, as reported in Chiquito-Almanza et al. [12]. Domesticated common bean was collected from irrigated and rain-fed commercial fields in Guanajuato and Jalisco during the spring-summer growing season of 2014 and in Nayarit, Mexico, during the fall-winter growing seasons of 2013–2014 and 2014–2015. Based on their distribution and symptom diversity, including asymptomatic plants, 42 samples of wild bean and 30 samples of domesticated bean were selected for this study. The sample IDs, species, cultivar, growth stage, type of bean, collection date, and coordinates were uploaded to the National Center for Biotechnology Information (NCBI) under the BioProject number PRJNA362733; sample accession numbers are indicated in Table S1.

2.2. Library Construction and Small RNA Sequencing

Total RNA was purified from 0.4 g of dried leaf tissue using TRIzol reagent following the manufacturer’s instructions (Invitrogen, Carlsbad, CA, USA); RNA quantity was determined using a NanoDrop 8000™ analyzer (Thermo Fisher Scientific, Waltham, MA, USA), and RNA quality was checked by agarose gel electrophoresis. Total RNA (3 μg) was shipped to Macrogen Inc. (Seoul, Korea), where sRNA libraries were constructed with the TruSeq® Small RNA Library Preparation Kit (Illumina, San Diego, CA, USA) and sequenced on an Illumina HiSeq 2000 system in 50-bp single-end mode. Seventy-two samples were used to construct 53 sRNA libraries, 26 of which were from domesticated bean species (25 from P. vulgaris and one from P. acutifolius), while 27 were from the wild bean species P. acutifolius, P. coccineus, P. leptostachyus, P. lunatus, P. microcarpus, and P. vulgaris. All libraries were constructed from a single sample except for the composite libraries CJ16-25Lc and CJ25-56Lc from domesticated beans and SG5-136Lc, SG7-245Lc, SJ9-34Lc, SJ19-24678Lc, and SN41-2345Lc from wild beans, which were made of the equimolar mixtures of total RNA from three to six plants collected from the same spot as the single sample in the population (Table S1).

2.3. Small RNA Bioinformatic Analysis

Quality check for raw sRNA reads was performed using the sRNA_clean.pl script (https://github.com/kentnf/VirusDetect/tree/master/tools/sRNA_clean, accessed on 16 October 2020) to trim the 3’ adaptor sequences and remove the sRNA reads without adaptor sequences and those of low quality or shorter than 15 nt after trimming. In order to determinate the size distribution of sRNA, trimmed sRNA sequences between 15–40 nt were aligned to the common bean genome (NCBI: GCA_000499845.1) using the Bowtie v1.1.1 program [13]; the average number of the sRNA mapped and unmapped to host were plotted. The de novo assembled of reads into contigs and identification of viruses were conducted according to Zheng et al. [14] using the online VirusDetect pipeline v1.6 with default parameters (http://virusdetect.feilab.net/cgi-bin/virusdetect/index.cgi, accessed on 18 October 2020); the common bean genome and plant virus reference database in the VirusDetect software was used as a reference to subtract host sRNA and for reference-guided assemblies through aligning sRNA sequences, and other parameters applied were coverage = 0.1 and depth = 5. SnapGene V4.1.9. (GSL Biotech LLC, Chicago, IL, USA) was used for the annotation and reconstruction of complete and near-complete genomes of viruses from the pooled sRNA reads using the complete genomes of the viruses with the highest bit score as a reference (Table S1). The resulting viral genomes were compared to the de novo assembled contigs to ensure consistency. The percentage of coverage of the reconstructed genome upon the viral genomes used as reference was also determined. The assembled genomes were considered complete if they cover 100% of the complete genome used as reference. The distribution patterns in domesticated and wild beans from MCD and Nayarit were analyzed through a heatmap that was built using the function heatmap2 from the package gplots [15] on R 4.0.3 [16]. Code is available in Supplementary Materials.

2.4. Validation of Viral Sequences by RT-PCR and Sequencing

The viruses identified by sRSA in each library were confirmed by PCR or RT-PCR in every single sample that integrated both single and composite libraries. The simultaneous extraction of both RNA and DNA viruses was performed using the cetyl trimethyl ammonium bromide (CTAB) method, according to Chiquito-Almanza et al. [12]. The absence of PCR inhibitors in total nucleic acids was tested by PCR amplification of the 26S rRNA reference gene, according to Montero-Tavera et al. [17]. Reverse transcription (RT) was conducted using 1 µg of total nucleic acids, random hexamers, and SuperScript III Reverse Transcriptase (Invitrogen, Carlsbad, CA, USA) according to the manufacturer’s instructions. The resulting cDNA was used as a template for the PCR detection of DNA and RNA viruses. Specific primers for each viral species were designed with GEMI v1.5.0 software [18] using multiple alignments of assembled contigs and complete viral genome sequences available in NCBI GenBank; the sequences were aligned using the algorithm ClustalW in BIOEDIT v7.2.5 [19]. The designed primers were analyzed with OligoAnalyzer [20], and their specificity was confirmed in silico using Primer-BLAST [21]. PCR amplifications were performed in a volume of 50 µL according to Chiquito-Almanza et al. [12]. The PCR program included one denaturation cycle at 94 °C for 3 min, followed by 35 cycles at 94 °C for 45 s, the optimal Tm of each primer pair (Table S2) for 1 min, and 72 °C for 2 min, with a final extension step at 72 °C for 7 min. Amplicons were analyzed by 1% agarose gel electrophoresis and sequenced directly using the dideoxy terminator method. The nucleotide sequences of the PCR products were compared with BLASTn [22] to the non-redundant database to validate their sequence identity.

2.5. Estimation of Virus Prevalence

Virus prevalence in each region was calculated considering individual plants as the percentage of the number of individual samples with at least one virus divided by the number of total individual samples. The Mann–Whitney U-test was used to determine the significance of the differences between wild and domesticated beans concerning viral prevalence in the MCD and Nayarit using the number of samples of individual plants that contained at least one virus from the total of individual samples tested as input data.

2.6. Phylogenetic Analysis

The nucleotide sequences of the complete and near-complete reconstructed genomes were used to establish phylogenetic relationships with other members of their respective viral families that displayed best bit scores in BLAST searches. The phylogenetic analysis was performed using the complete genome nucleotide sequences of the isolates of Bean common mosaic virus (BCMV), Bean common mosaic necrosis virus (BCMNV), Peanut mottle virus (PeMoV), Phaseolus vulgaris alphaendornavirus 1 (PvEV-1), Phaseolus vulgaris alphaendornavirus 2 (PvEV-2), Cowpea mild mottle virus (CPMMV), the DNA A components of Bean golden yellow mosaic virus (BGYMV), and Bean latent virus (BLV); the amino acid sequences of the RNA-dependent RNA polymerase (RdRp) of the Bean-associated cytorhabdovirus (BaCV) isolate CN3; the putative novel alphaendornaviruses Phaseolus leptostachyus alphaendornavirus (PlepEV) and Phaseolus lunatus alphaendornavirus (PluEV). The RdRp domains of PlepEV and PluEV were predicted using the NCBI Conserved Domain-Search tool (https://www.ncbi.nlm.nih.gov/Structure/cdd/wrpsb.cgi, accessed on 7 December 2020). The MEGA X package [23] was used to perform the alignment of nucleotide and amino acid sequences for the selection of the best-fit substitution model for each set of data and the construction of phylogenetic trees. The nucleotide and amino acid sequences were aligned using the ClustalW algorithm with the default settings, and the phylogenetic trees were constructed using the maximum likelihood method with 1000 bootstraps. Nodes with less than 70% bootstrap values were collapsed with TreeGraph 2 ver. 2.14.0-771b [24].

3. Results

3.1. Viral Species Present in the Domesticated and Wild Bean

The sequencing of the sRNA libraries generated between 3.4 and 52.3 million raw reads per library. Trimming adapters and low-quality reads resulted in 2 to 49.2 million clean reads. In order to determine the size distribution of sRNA, trimmed reads between 15–40 nt were aligned to the common bean genome. Between 34% and 89% of the sRNA in each library were from the host; the unmapped sRNA reads with highly enriched vsiRNAs were between 11% and 66% of the total sRNA counts in each library (Table S3). The most abundant class of unmapped to host reads was 21–24 nt in length with a peak of 21 followed by 24 nt (Figure 1).
Forty-seven complete or nearly complete genomes of nine known viral species, two sequences of 3581 bp and 2878 bp of two putative novel alphaendornaviruses, and five badnavirus- and pararetrovirus-like sequences were assembled (Table S1). The viral species were included in five genera, where the subgroup of the (+) ssRNA viruses included the carlavirus, CPMMV, and the potyviruses BCMV, BCMNV and PeMoV. The subgroup of dsRNA viruses included the alphaendornaviruses PvEV-1 and PvEV-2, and the putative novel viruses PlepEV and PluEV; the subgroup of ssDNA included the begomoviruses BGYMV and BLV, and the (−) ssRNA subgroup included the cytorhabdovirus BaCV (Table 1). The geographic distribution of the identified viruses is shown in Figure 2. With the exception of five partial genomes of PvEV-1 and four of PvEV-2 from libraries CG3Ls, CG9Ls, CG13Ls, SG1Ls, CJ17-4Ls that were not submitted to the GenBank because their sequences showed low quality; the GenBank accession numbers of the assembled genomes and amplified fragments are indicated in Table S1. The phylogenetic tree of all the viruses identified is shown in Figures S1–S4.
Among the viruses identified, BaCV and BLV correspond to species recently described; thus, the characterization of their genomes was performed. The isolate of BLV identified in sample CN30 from Nayarit is a novel bipartite begomovirus whose complete genome was reconstructed and characterized by Martínez-Marrero et al. [25]. Meanwhile, the genome of the cytorhabdovirus identified in sample CN3, hereinafter referred to as BaCV isolate CN3 (BaCV-CN3; GenBank accession MT792847), consisted of 13,427 nt. It presents ORFs of the five canonical rhabdovirus structural proteins: a nucleoprotein (N) ORF1 (451 amino acids, aa), a phosphoprotein (P) ORF2 (445 aa), a matrix protein (M) ORF5 (214 aa), a glycoprotein (G) ORF6 (519 aa), and an RNA-dependent RNA polymerase (L) ORF7 (2113 aa), as well as a putative movement protein (P3) ORF3 (197 aa), and the hypothetical small protein P4 ORF4 (67 aa) (Figure S5A). Except for the 3′/N region, the P3/P4 region, which does not have the conserved consensus sequences, and the region L/5′ that had 12 of the 17 conserved nucleotides, all the intergenic regions have the highly conserved sequence 3′-TAAGAAAAACYGGGAKC-5′ (Figure S5B). BaCV-CN3 is closely related to the cytorhabdovirus BaCV strain BaCV-BR-GO ([26], GenBank accession MK202584], BaCV-LUZ ([27], GenBank accession MT811775), and Papaya cytorhabdovirus (PpVE) from Ecuador ([28], GenBank accession MH282832] (Figure S5C).
Except for the sample CG32Ls of P. acutifolius, all the samples of domesticated beans were from P. vulgaris. The number of virus genera found in the MCD was lower than in the western state of Nayarit (Figure 3A). The viruses identified in the MCD were members of the potyvirus (BCMV, BCMNV, and PeMoV) and alphaendornavirus genera (PvEV-1, PvEV-2, and PlepEV), which are seed-borne except for the putative PlepEV, for which seed-borne transmission is unknown. In Nayarit, a higher diversity of viral species was identified; in addition to BCMV and BCMNV, we identified the carlavirus CPMMV, the cytorhabdovirus BaCV, the putative alphaendornavirus PluEV, and the begomoviruses BGYMV and BLV (Table 1 and Figure 3B).

3.2. Mixed Infections and Prevalence of Viruses in Domesticated and Wild Bean

In order to determine the presence of mixed infections, single and multiple infections were recorded. In samples from the MCD, there were 12 domesticated beans with single infections (10 with BCMV and two with PvEV-1) and eight with double infections (five with PvEV1 + PvEV-2 and three with BCMV + BCMNV). Conversely, among the wild bean samples from this region, 24 samples harbored no viruses, and five showed single infection (two with BCMV and one each with PvEV-1, PlepEV, and BCMNV) (Figure 3C). In the samples of wild bean from the MCD, BCMV was found infecting P. vulgaris (samples SJ9 and SG27-3) and P. leptostachyus (sample SJ8-5), and the latter sample was also infected with PeMoV; BCMNV was found infecting P. leptostachyus (sample SJ12), while PvEV-1 and PlepEV infected two samples, one of P. vulgaris (sample SG1-1) and one of P. leptostachyus (sample SG29) (Table S1).
In the samples from domesticated beans from Nayarit, we detected single infections (2 with BCMV and one with BCMNV), double infections (one each with BCMNV + CPMMV, BCMNV + BaCV, BCMNV + BGYMV, and BCMNV + BCMV) and triple infections (one each with BCMNV + BCMV + BGYMV and BCMNV + BCMV + BLV). In the wild beans, no virus was detected in 10 samples (Figure 3C), and single infections of BCMNV and PluEV were found in two samples. The viruses identified in wild beans from Nayarit were BCMNV, infecting P. vulgaris (sample SN44-5), and a putative PluEV, infecting P. lunatus (SN35) (Table S1).
These results show that BCMV and BCMNV were the most frequent viruses identified in domesticated and wild beans (Figure 3B and Table S1), and single and double infections are frequent in the MCD in comparison with Nayarit, where single, double, and triple infections occur (Figure 3C). Furthermore, in all the samples of domesticated beans, at least one virus species was identified, while in wild beans, no virus was found in 34 samples (Figure 3B,C). In order to quantify these differences, the prevalence of viruses in domesticated and wild beans in each region was calculated using Mann–Whitney U-tests to compare the differences. Because all the viruses identified by sRSA were confirmed by PCR or RT-PCR in every single sample that integrated both single and composite libraries, their prevalence was calculated considering individual plants.
Twenty-one and 30 of the samples used for the construction of the 30 and 42 sRNA libraries of domesticated and wild beans, respectively, were from the MCD. The prevalence of viruses in the domesticated bean was significantly higher than in wild bean (p < 0.001). In the MCD, 20 of the 21 samples of domesticate bean (95.2%) were infected by at least one viral species, while in wild beans, only six of 30 samples (20%) were infected; in Nayarit, these percentages were 100% (nine of the nine samples) and 16.7% (two of the 12 samples), respectively (Figure 3D).

3.3. Viral Distribution Patterns in Domesticated and Wild Bean

In order to determine the distribution patterns of the viruses in domesticated and wild beans, a heatmap analysis was performed, using the coverage % of the assembled viral genome sequences upon the viral genome as reference (Figure 4). With the exception of the partial sequences of the putative novel alphaendornaviruses PluEV and PlepEV for which only a fragment was assembled, the % coverage of the assembled genomes ranged 83–100% (Figure 4, Table S1).
The dendrogram on the left side of Figure 4 shows that grouping was based on the viruses that were found in each library, which is in accord with the domestication level (domesticated and wild). The viruses were grouped into three clusters: clusters 1 and 2 formed compact groups, while cluster 3 is integrated by two subgroups, clusters 3a and 3b (Figure 4). Clusters 1, 2, and 3a include the viruses mainly identified in libraries constructed from plant samples of domesticated beans, while in cluster 3b, with the exception of the putative alphaendornaviruses PluEV and PlepEV where no viruses were detected, these exclusively include samples of wild beans from both, the MCD and Nayarit.
Cluster 1 includes the alphaendornaviruses PvEV-1 and PvEV-2; both species were found exclusively in samples from the MCD, mainly in domesticated beans from Guanajuato (libraries CG3Ls, CG9Ls, CG13Ls, cg18Ls, CG17Ls, and CG1Ls) and Jalisco (library CJ14-4Ls); meanwhile, PvEV-1 was also identified in one sample of wild bean from Guanajuato (library SG1Ls).
Cluster 2 encompasses the higher diversity of viral species. Among the six viruses in this cluster, BCMNV is present in all the libraries and includes mainly samples of domesticated beans from Nayarit, and two samples of wild beans, one from Nayarit (SN44-5Ls) and one from Jalisco (SJ12Ls). This result confirms the high prevalence of BCMNV in common beans grown in Nayarit and the presence of this virus in wild and domesticated bean from Nayarit and Jalisco (library CJ25-56Lc in cluster 3a).
Cluster 3a includes all libraries infected with BCMV, except those from libraries CN61Ls and CN30Ls from cluster 2. In addition, PeMov was identified in mixed infection along with BCMV in library SJ8-5Ls. Although most of the libraries in this cluster are from samples of domesticated beans, it also includes samples of wild beans (SJ9-34Lc and SG27-3Ls). Samples in this cluster were collected from Nayarit, Guanajuato and Jalisco. This result indicates the wide distribution of BCMV in the common bean that is grown in this region and confirms the presence of this virus in wild beans in the MCD.

3.4. Identification of Novel Viral Species and Viral-Like Sequences

Several assembled contigs without detectable nucleotide homology but with amino acid identity to reference virus protein sequences were identified in 11 of the libraries (Table S4). In order to determine if these sequences corresponded to novel viral species, the assembled contigs were used as references to design specific primers. Through standard DNA amplification techniques and sequencing, the partial nucleotide sequences of two putative novel alphaendornaviruses and five virus-like sequences were reconstructed. Below we describe the characteristics of the reconstructed sequences of both putative novel alphaendornaviruses and the five virus-like sequences.

3.4.1. Identification of Two Putative Alphaendornaviruses

The partial nucleotide sequences of two putative novel alphaendornaviruses were reconstructed from the wild asymptomatic samples SG29 from P. leptostachyus collected in Guanajuato and SN35 from P. lunatus collected in Nayarit. For identification purposes, the names Phaseolus leptostachyus alphaendornavirus (PlepEV) and Phaseolus lunatus alphaendornavirus (PluEV) were assigned. The partial nucleotide sequence of PlepEV consisted of 3581 nt (GenBank accession MT792848) and included one partial ORF of 3540 nt that was translated into a fragment of a putative polyprotein of 1179 aa, sharing 45.9% and 43.9% aa identity with the carboxy-terminal polyprotein of the unclassified alphaendornavirus Geranium carolinianum endornavirus (GcEV; unpublished; GenBank accession MH577297) and the Phaseolus vulgaris endornavirus 3 isolate LA (PvEV-3; [29]; GenBank accession MG242064), respectively. The query coverage of both aa sequences was 96%. Furthermore, the partial nucleotide sequence of PluEV consisted of 2878 nt (GenBank accession MT792849), which includes one partial ORF of 2874 nt that was translated into a putative polyprotein fragment of 958 aa, sharing 65.1% aa identity (98% query cover) with the carboxy-terminal polyprotein of the PvEV-3 isolate LA and 46% aa identity (97% query cover) with the carboxy-terminal polyprotein of GcEV.
The aa identity between the predicted polyproteins of PlepEV and PluEV was 44.9%. The nucleotide sequences of the fragments amplified by PCR with specific primers (Table S2) for PlepEV (1076 bp) and PluEV (1090 bp) exhibited 100% identity with the respective reconstructed sequences. The predicted amino acid sequences of both reconstructed sequences of PlepEV and PluEV included the full RdRp domain. The phylogenetic analysis of this domain with different alpha- and betaendornaviruses show that PlepEV and PluEV clustered together with PvEV-3 and GcEV in the clade corresponding to the alphaendornavirus genus (Figure 5).

3.4.2. Identification of Virus-Like Sequences

Five virus-like sequences with amino acid identity to the viral reference sequences were reconstructed (Table S5). For identification purposes, the sequences were named based on the sample from where they were obtained and the genus of the virus to which they match aa identity, as Phaseolus vulgaris badnavirus-like sequence CG17 (PvBV-CG17), Phaseolus coccineus badnavirus-like sequence SJ10-3 (PcBV-SJ10-3), Phaseolus coccineus pararetrovirus-like sequence SJ10-3 (PcPRV-SJ10-3), Phaseolus lunatus badnavirus-like sequence SN35 (PluBV-SN35), and Phaseolus lunatus pararetrovirus-like sequence SN35 (PluPRV-SN35).
The sequence of PvBV-CG17 consisted of 1629 bp and was reconstructed from a sample from common bean cv. Azufrado Higuera with symptoms of local chlorosis lesions; in this sample, PvEV-1 (GenBank accession MG640416) was also identified. The sequence of PcBV-SJ10-3, consisting of 2604 bp, was reconstructed from a wild P. coccineus sample with symptoms of local chlorosis lesions, while the sequence PluBV-SN35, consisting of 1378 bp, was reconstructed from a symptomless wild P. lunatus sample; in both samples, SJ10-3 and SN35, the pararetrovirus-like sequences PcPRV-SJ10-3 (2588 bp) and PluPRV-SN35 (1440 bp) were also identified; these translated into 862 aa and 480 aa proteins, respectively, with 49% and 32% aa identity to Citrus endogenous pararetrovirus polyprotein (GenBank accession KF800044).
Both the PvBV-CG17 and PluBV-SN35 sequences exhibited one partial ORF, which was translated into putative proteins of 543 and 456 aa, respectively, with 39.0% and 35.6% aa identity to the putative P3 polyprotein (ORF3) of the unclassified Badnavirus-like pelargonium vein banding virus (PVBV; GenBank accessionGQ428155). On the other hand, the sequence of PcBV-SJ10-3 exhibited two partial ORFs, one was translated into a putative protein of 58 aa with 41.3% aa identity to the ORF2 protein of the badnavirus cacao swollen shoot virus (CSSV; GenBank accession MF642734), and the other was translated into a putative protein of 808 aa with 43.4% aa identity to the P3 putative polyprotein of PVBV. Although the reconstructed sequences of PvBV-CG17, PcBV-SJ10-3, and PluBV-SN35 shared low nucleotide similarity with the available sequences of badnavirus in the GenBank, suggesting that they could correspond to a novel badnavirus species, the reconstructed fragment did not include the RT/RNAse H domain required for the demarcation of the species. Additionally, since badnaviruses and pararetroviruses had been reported as integrated fragmented sequences, the circularity of the genome needs to be confirmed.

4. Discussion

4.1. Identification and Distribution of Viruses in Domesticated and Wild Bean

In addition to the MCD, samples of wild and domesticated Phaseolus spp. were collected in the western state of Nayarit since the occurrence of virus damaging in the common bean crop is particularly high in this area [12,30], and it is located near the MCD. The collection of samples at a site was based on the presence of viral symptoms, and symptomless plants were also sampled. However, in most of the samples of wild Phaseolus spp. with presumptive viral symptoms, no viruses were identified; the symptoms observed in those samples could be due to nutrient deficiencies or the effect of other biotic or abiotic stresses that may show similarities to mosaic or mottling [31]. Nine species were identified in the MCD and Nayarit, with a higher diversity of viral species in Nayarit, which includes six species (BCMV, BCMNV, BGYMV, BLV, CPMMV, and BaCV). BCMNV was mainly found in Nayarit, and BCMV was widely distributed in both the MCD and Nayarit. PeMoV was identified in one sample from Jalisco; PvEV1 and PvEV-2 were found exclusively in Guanajuato (Table 1, Figure 2 and Figure 4).
The number of virus genera found in the MCD was lower than in the western state of Nayarit (Figure 3A); a possible explanation for the higher incidence of viruses in Nayarit could be that as an important bean producing area during the fall-winter cycle, Nayarit was one of the chosen locations for the establishment of international common bean trials in Mexico, which facilitated the introduction of the seed-borne viruses [32,33,34]. Additionally, the abundant presence of the whitefly B. tabaci represents a constraint to vegetable production [35]. Moreover, the detection of a high proportion of mixed infections in Nayarit suggests that the climatic conditions and cultural practices, such as the exchange of seeds between growers or the use of the grain of susceptible cultivars as seed, favor the dissemination of viruses that is further enhanced by the presence of the vector.

4.2. Prevalence of Viruses in Domesticated and Wild Bean

The evidence indicates that the incidence of viral diseases in the wild populations sampled is low in contrast to domesticated bean, in which almost all samples showed single, double, or triple infections (Figure 3B,C). The higher prevalence of viruses in the domesticated bean (Figure 3D) could be the result of anthropogenic influence. Cultivation is a major cause of environmental heterogeneity [36], and monoculture provides suitable conditions for the fast emergence and spread of plant viruses and their vectors [1], which is a situation that does not occur in wild bean populations that are in distinct ecological conditions. Studies of virus infection in chiltepin populations (Capsicum annuum var. glabriusculum), a wild pepper plant under different anthropogenic influences, including cultivation and incipient domestication, support the fact that ecological conditions are a valuable predictor of virus prevalence and that human management is associated with a reduction of habitat biodiversity and genetic diversity of host species; in addition, as host plant density is increased, the viral disease risk is enhanced due to plant-to-plant proximity. Furthermore, virus infection decreases the survival and fecundity of wild chiltepin populations, hence negatively affecting fitness [36]. It is also likely that the variation of traits, such as resistance to viruses and other quantitative traits found in wild beans, are not present in the domesticated common bean due to genetic diversity bottlenecks during bean domestication [37,38]. Nonetheless, the rate of the discovery of putative novel viruses and virus-like sequences in wild beans was high in comparison to the domesticated bean, in which almost all the species identified were already known. These results agree with those of Bernardo et al. [39], who found that prevalence and identified family-level virus diversity were greatest in cultivated areas, and unknown virus species were found primarily from uncultivated plants. Ultimately, the number of studies on wild bean populations, as compared to the domesticated form, is low.

4.3. Characteristics of Known Viruses Identified

The viruses BCMV, BCMNV, BGYMV, and CPMMV have been documented infecting common bean and wild bean in Guanajuato, Jalisco, and Nayarit, Mexico [12,30,40]. In contrast, the viruses PvEV-1, PvEV-2, PeMoV, and BaCV have not been reported infecting Phaseolus spp. in Mexico before.
Predicting epidemics of plant virus diseases and the causes that facilitate pathogen emergence is complex and depends on multiple factors, reviewed in [1,41,42]. Besides these factors, the spread from the main reservoir into a new environment, the establishment of infection, and effective mechanisms of transmission between hosts are steps required for emergence to occur [43]. In this context, the virulence and the mechanism of transmission of a virus, particularly its capacity for seed-borne and insect transmission, are some of the features intrinsic to the virus that could contribute to define its potential to become an emerging pathogen. BGYMV, BaCV, and CPMMV are transmitted by B. tabaci. BGYMV is an important non-seed-transmitted virus in areas with a high incidence of B. tabaci; the disease incidence and the extent of losses vary depending on the populations of the whitefly, the susceptibility of the bean cultivar, and environmental conditions, particularly rainfall, which affects both B. tabaci populations and cultural practices [44]. The vector that spreads BLV and the affectation that may cause to the common bean production are unknown; under experimental conditions, this begomovirus does not cause symptoms in common bean [25].
Based on the high nucleotide sequence identity between the cytorhabdovirus BaCV-BR-GO and PpVE, it was proposed that the virus species would be named Papaya cytorhabdovirus with strains PpVE infecting papayas and BaCV infecting beans [45]; thus, BaCV-CN3 is likely a strain of BaCV since it infects common bean and is closely related to the strains BaCV-BR-GO, BaCV-LUZ, as shown in the phylogenetic tree based on the RdRp amino acid sequences (Figure S1C). The vector of PpVE is unknown, and although BaCV is the first rhabdovirus reported to be transmitted to common beans by whitefly B. tabaci MEAM1 [27], the transmissibility of BaCV-CN3 by this vector was not determined. The symptoms induced by BaCV-CN3 could neither be confirmed because the sample from which it was identified was co-infected with BCMNV (GenBank accession MK069983). Interestingly, BaCV and PpVE have also been identified in mixed infections. Pinheiro-Lima et al. [27] found that 85 of the 91 common bean samples that tested positive to BaCV were co-infected with CPMMV, BGMV, or with CPMMV and BGMV, and although six samples were “single infected,” it is also possible that these samples were infected with another virus distinct to the ones tested; similarly, all the plants that tested positive to PpVE were also positive for Papaya ringspot virus (PRSV), and no differences in leaf symptoms were observed between PRSV single-infected plants and plants co-infected with both viruses [28]. Therefore, it is important to investigate if this group of viruses requires a helper virus, are seed-borne transmitted, and cause important economic losses.
The presence of endornaviruses PvEV-1 and PvEV-2 in almost all the domesticated and wild samples from Guanajuato (Figure 4) could be related to the domestication process of the common bean, which likely occurred in the Lerma-Santiago basin in the states of Jalisco and Guanajuato, Mexico [7]. Endornaviruses are persistent dsRNA viruses that are transmitted vertically [46] and develop a highly symbiotic relationship with the host plant [47]. PvEV-1 and PvEV-2 are transmitted to the progeny of P. vulgaris plants at rates close to 100% [48] and have been detected in many genotypes of Mesoamerican origin but rarely in genotypes of the Andean pool, suggesting that these two endornaviruses likely originated from infected wild species that were selected during common bean domestication and introduced into some cultivars during the breeding process [49].
The main host of PeMoV is peanuts (Arachis hypogea), but it also infects susceptible genotypes of P. vulgaris, causing various reactions that range from mosaic, banding, leaf deformation and blistering, to systemic necrosis in genotypes that have only the dominant I gene; these reactions are similar to those caused by necrotic strains of BCMV and BCMNV. Although PeMoV is mechanically and non-persistently transmitted by aphids, transmission to common bean seeds is low (1%) [50].
Among the viruses identified in this study, BCMV and BCMNV likely present the highest phytosanitary risk due to their high incidence in the region and their efficient mechanisms of seed and vector transmission [51]. BCMV and BCMNV are considered one of the major production constraints of common bean in Latin America, the Caribbean, and East, West, and Southern Africa and the most worldwide spread viruses [52,53]; the percentage of seed transmission of these viruses is as high as 50% [51] and hence are readily spread within and between countries [50,54]. BCMV can remain viable in bean seeds for more than 30 years [54], and under field conditions, the infected seedlings originated from BCMV- or BCMNV-infected seeds act as primary foci, and subsequently, the virus spreads through several aphid species in a non-persistent manner [55]. Yield losses of 53–83% in beans have been attributed to BCMV [56], and even symptomless BCMV infections can induce yield losses higher than 50% in susceptible common bean genotypes [2]. The high incidence of BCMV and BCMNV in wild Phaseolus spp. from the MCD and the western state of Nayarit suggest that the transmission of these viruses through insect vectors is highly efficient. Further studies are needed to know how and from where these viruses came into the wild bean populations of the MCD and Nayarit and if transmission occurred from cultivated into wild bean populations. The other virus that has the potential to emerge as an important pathogen is CPMMV, which is emerging as a soybean pathogen in Brazil, and it is re-emerging in cultivated common beans [57], particularly in genetically modified common bean cultivars resistant to BGMV [58]. CPMMV is the only Carlavirus that is not transmitted by aphids, and it is nonpersistently transmitted by B. tabaci (MEAM1), with a period of access to acquisition of 10 min and an inoculation period of 5 min [59]. Although the seed transmission of CPMMV has been reported, the different rates of seed transmission in cowpea, soybean, and common bean suggest that this characteristic depends on the isolate, the strain, and the cultivar of the host, which can make the control of the disease difficult and facilitate its spread due to the existence of asymptomatic infections [50,58]. Furthermore, symptoms are rather variable, such that some infections of CPMMV have been confounded as new and different viral infections [58,60]. In soybean, CPMMV causes stem necrosis disease, but distinct symptoms have been observed in the same cultivar within the same cultivation field [58]. In common beans, CPMMV causes angular mosaic disease, including symptoms of systemic mosaic and leaf distortion, although the infection may be asymptomatic depending on the cultivar [40,50]. In Mexico, this virus has been identified in infections of common bean (P. vulgaris) and wild bean (P. vulgaris, P. leptostachyus, and P. coccineus), and some RT-PCR-positive samples did not present symptoms [40].
CPMMV is classified in Annex 1/A1 in the European Union and in the United States as a quarantine pest, where it is considered to represent a potential risk to soybeans, beans, cowpea, and some species of the Solanaceae family [61,62]. Mexico does not regulate CPMMV, but its occurrence in beans could have negative consequences for sales abroad. To determine the phytosanitary risk of CPMMV, it is necessary to understand its viral genomic diversity, host range, capacity for seed-borne transmission, distribution, prevalence, and incidence, especially in areas with high whitefly vector pressure and where common beans are grown from seed or grain produced in situ [40].

4.4. Putative Novel Viruses Identified

The assembled sequences of PlepEV and PluEV share amino acid identity with the RdRp domains of PvEV-3 and GcEV, suggesting that they could be novel alphaendornaviruses; however, their complete genome sequences will be required for their classification. The plant samples in which PlepEV and PluEV were identified were symptomless, and like other endornaviruses that have been reported in economically important crops, their effect on the host has not been fully determined. For instance, PvEV-1 and PvEV-2 were not associated with visible pathogenic effects on P. vulgaris cv Black Turtle Soup [63]; also, the interaction of these viruses with other common bean viruses and pathogens should be tested for possible interactions.
The badnavirus- and pararetrovirus-like sequences reconstructed in this study were obtained by sequencing since multiple single-nucleotide variations (SNVs) in the assembled contigs made scaffolding difficult. Some badnaviruses occur as integrated, complete, fragmented, and/or rearranged genomic sequences of the original circular viral DNA in their host plant genomes, in which case they are referred to as endogenous pararetroviruses or endogenous viral elements [64]. Their integration takes place by illegitimate recombination into host genomes, and their presence is not necessarily associated with infection. The identification of putative badnaviruses or pararetroviruses highlights the need for further specific studies to determine their biology and phytosanitary risks.

5. Conclusions

Higher viral diversity was found in Nayarit in comparison with the MCD. BCMV and BCMNV were the most prevalent viruses in wild and domesticated beans; these viruses, along with CPMMV, have the potential to emerge as an important pathogen because they are seed-borne and non-persistently transmitted viruses. Integrated disease management approaches using virus-resistant common bean cultivars and efficient management of insect vectors and weeds can reduce the dispersion and the production losses associated with these viruses. More studies are needed to confirm the putative viruses identified in this study and determine the effect of human management on population heterogeneity, biodiversity loss, and disease risk that these viruses represent for the domesticated common bean and wild bean populations.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/v13061153/s1, Figure S1: Phylogeny of potyvirus species identified in domesticated and wild bean in the Mesoamerican center of domestication of common bean, and Nayarit, Mexico, Figure S2: Phylogeny of begomovirus species identified in domesticated and wild bean in the Mesoamerican center of domestication of common bean, and Nayarit, Mexico, Figure S3: Phylogeny of alphaendornavirus species identified in domesticated and wild bean in the Mesoamerican center of domestication of common bean, and Nayarit, Mexico, Figure S4: Phylogeny of isolates of CPMMV identified in domesticated and wild bean in the Mesoamerican center of domestication of common bean, and Nayarit, Mexico, Figure S5: Genome of organization and phylogeny of bean-associated cytorhabdovirus isolate CN3 (BaCV-CN3), Table S1: Summary of samples, libraries and viruses identified in domesticated and wild bean from the Mesoamerican center of domestication of common bean and Nayarit, Mexico, Table S2: Primers used for the confirmation of the viruses identified by sRSA, Table S3: Number of raw and clean reads obtained for high throughput sequencing of the sRNA libraries of domesticated and wild bean from the Mesoamerican center of domestication common bean and Nayarit, Mexico, Table S4: Contigs identified through Blastx with amino acid identity to reference viruses assembled from small RNA libraries from samples of domesticated and wild bean plants collected in the Mesoamerican center of domestication of common bean and Nayarit, Mexico, Table S5: Nucleotide sequences of the virus-like identified in small RNA libraries of domesticated and wild bean from the Mesoamerican center of domestication of common bean and Nayarit, Mexico.

Author Contributions

Conceptualization, J.A.A.-G., E.C.-A. and J.L.A.-L.; software and formal analysis, J.C.-P.; investigation, validation, and data curation, E.C.-A.; project administration, supervision, and methodology, J.L.A.-L.; formal analysis, writing—original draft preparation, J.A.A.-G., J.C.-P., J.L.A.-L., E.C.-A. and V.M.-T.; writing—review and editing, J.L.A.-L., L.A.M.-A. and V.M.-T. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by INIFAP, grant number 1551434856, and by CONACYT project No. PDCPN-2013/214641.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The sRNA reads were uploaded to the NCBI Sequence Read Archive (SRA) under accession numbers SRX3217971-SRX3218007. The sequences of the viral genomes reconstructed, as well as the PCR-amplified fragments for the validation of the viruses detected, were uploaded to the GenBank, and the accession numbers are indicated in Table S1.

Acknowledgments

The authors acknowledge Wilmer J. Cuellar for the training on the sRSA tech.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. Jones, R.A. Plant virus emergence and evolution: Origins, new encounter scenarios, factors driving emergence, effects of changing world conditions, and prospects for control. Virus Res. 2009, 141, 113–130. [Google Scholar] [CrossRef]
  2. Morales, J.F. Common beans. In Natural Resistance Mechanisms of Plants to Viruses, 1st ed.; Loebenstein, G., Carr, J.P., Eds.; Springer: Dordrecht, The Netherlands, 2006; pp. 367–382. [Google Scholar]
  3. FAOSTAT (Food and Agriculture Organization of the United Nations). Available online: http://www.fao.org/faostat/es/#data/QC (accessed on 16 October 2020).
  4. Bellucci, E.; Bitocchi, E.; Rau, D.; Rodriguez, M.; Biagetti, E.; Giardini, A.; Attene, G.; Nanni, L.; Papa, R. Genomics of origin, domestication and evolution of Phaseolus vulgaris. In Genomics of Plant Genetic Resources, 1st ed.; Tuberosa, R., Graner, A., Frison, E., Eds.; Springer: Dordrecht, The Netherlands, 2014; Volume 1, pp. 483–507. [Google Scholar]
  5. Rendón-Anaya, M.; Montero-Vargas, J.M.; Saburido-Álvarez, S.; Vlasova, A.; Capella-Gutierrez, S.; Ordaz-Ortiz, J.J.; Aguilar, O.M.; Vianello-Brondani, R.P.; Santalla, M.; Delaye, L.; et al. Genomic history of the origin and domestication of common bean unveils its closest sister species. Genome Biol. 2017, 18, 60. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Bitocchi, E.; Nanni, L.; Bellucci, E.; Rossi, M.; Giardini, A.; Zeuli, P.S.; Logozzo, G.; Stougaard, J.; McClean, P.; Attene, G.; et al. Mesoamerican origin of the common bean (Phaseolus vulgaris L.) is revealed by sequence data. Proc. Natl. Acad. Sci. USA 2012, 109, E788–E796. [Google Scholar] [CrossRef] [Green Version]
  7. Kwak, M.; Kami, J.A.; Gepts, P. The putative Mesoamerican domestication center of Phaseolus vulgaris is located in the Lerma–Santiago Basin of Mexico. Crop Sci. 2009, 49, 554–563. [Google Scholar] [CrossRef] [Green Version]
  8. Mwaipopo, B.; Nchimbi-Msolla, S.; Njau, P.J.; Mark, D.; Mbanzibwa, D.R. Comprehensive surveys of Bean common mosaic virus and Bean common mosaic necrosis virus and molecular evidence for occurrence of other Phaseolus vulgaris viruses in Tanzania. Plant Dis. 2018, 102, 2361–2370. [Google Scholar] [CrossRef] [Green Version]
  9. Jones, S.; Baizan-Edge, A.; MacFarlane, S.; Torrance, L. Viral diagnostics in plants using next generation sequencing: Computational analysis in practice. Front. Plant Sci. 2017, 8, 1770. [Google Scholar] [CrossRef] [PubMed]
  10. Kreuze, J. siRNA deep sequencing and assembly: Piecing together viral infections. In Detection and Diagnostics of Plant Pathogens. Plant Pathology in the 21st Century, 1st ed.; Gullino, M.L., Bonants, P.J.M., Eds.; Springer: Dordrecht, The Netherlands, 2014; pp. 21–38. [Google Scholar]
  11. Santala, J.; Valkonen, J. Sensitivity of Small RNA-based detection of plant viruses. Front. Microbiol. 2018, 9, 939. [Google Scholar] [CrossRef] [Green Version]
  12. Chiquito-Almanza, E.; Acosta-Gallegos, J.A.; García-Álvarez, N.C.; Garrido-Ramírez, E.R.; Montero-Tavera, V.; Guevara-Olvera, L.; Anaya-López, J.L. Simultaneous detection of both RNA and DNA viruses infecting dry bean and occurrence of mixed infections by BGYMV, BCMV and BCMNV in the Central-west region of Mexico. Viruses 2017, 9, 63. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Langmead, B.; Trapnell, C.; Pop, M.; Salzberg, S.L. Ultrafast and memory-efficient alignment of short DNA sequences to the human genome. Genome Biol. 2009, 10, R25. [Google Scholar] [CrossRef] [Green Version]
  14. Zheng, Y.; Gao, S.; Padmanabhan, C.; Li, R.; Galvez, M.; Gutierrez, D.; Fuentes, S.; Ling, K.-S.; Kreuze, J.; Fei, Z. VirusDetect: An automated pipeline for efficient virus discovery using deep sequencing of small RNAs. Virology 2017, 500, 130–138. [Google Scholar] [CrossRef]
  15. Warnes, G.R.; Bolker, B.; Bonebakker, L.; Gentleman, R.; Huber, W.; Liaw, T.A.; Lumley, T.; Maechler, M.; Magnusson, A.; Moeller, S.; et al. Gplots: Various R Programming Tools for Plotting Data. R Package Version 3.1.1. 2020. Available online: https://CRAN.R-project.org/package=gplots (accessed on 17 May 2021).
  16. R Core Team. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2020; Available online: https://www.R-project.org/ (accessed on 17 May 2021).
  17. Montero-Tavera, V.; Escobedo-Landín, M.A.; Acosta-Gallegos, J.A.; Anaya-López, J.L.; Ruiz-Nieto, J.E. 26S: Novel reference gene from leaves and roots of common bean for biotic stress expression studies based on PCR. Legume Res. 2017, 40, 429–433. [Google Scholar] [CrossRef]
  18. Sobhy, H.; Colson, P. Gemi: PCR primers prediction from multiple alignments. Comp. Funct. Genom. 2012, 2012, 783138. [Google Scholar] [CrossRef] [PubMed]
  19. Hall, T.A. BioEdit: A user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT. Nucleic Acids Symp. Ser. 1999, 41, 95–98. [Google Scholar] [CrossRef]
  20. Owczarzy, R.; Tataurov, A.V.; Wu, Y.; Manthey, J.A.; McQuisten, K.A.; Almabrazi, H.G.; Pedersen, K.F.; Lin, Y.; Garretson, J.; McEntaggart, N.O.; et al. IDT SciTools: A suite for analysis and design of nucleic acid oligomers. Nucleic Acids Res. 2008, 36, W163–W169. [Google Scholar] [CrossRef]
  21. Ye, J.; Coulouris, G.; Zaretskaya, I.; Cutcutache, I.; Rozen, S.; Madden, T. Primer-BLAST: A tool to design target-specific primers for polymerase chain reaction. BMC Bioinform. 2012, 13, 134. [Google Scholar] [CrossRef] [Green Version]
  22. Altschul, S.F.; Gish, W.; Miller, W.; Myers, E.W.; Lipman, D.J. Basic local alignment search tool. J. Mol. Biol. 1990, 215, 403–410. [Google Scholar] [CrossRef]
  23. Kumar, S.; Stecher, G.; Li, M.; Knyaz, C.; Tamura, K. MEGA X: Molecular Evolutionary Genetics Analysis across computing platforms. Mol. Biol. Evol. 2018, 35, 1547–1549. [Google Scholar] [CrossRef]
  24. Stöver, B.C.; Müller, K.F. TreeGraph 2: Combining and visualizing evidence from different phylogenetic analyses. BMC Bioinform. 2010, 11, 7. [Google Scholar] [CrossRef] [Green Version]
  25. Martínez-Marrero, N.; Avalos-Calleros, J.A.; Chiquito-Almanza, E.; Acosta-Gallegos, J.A.; Ambriz-Granados, S.; Anaya-López, J.L.; Argüello-Astorga, G.R. A new begomovirus isolated from a potyvirus-infected bean plant causes asymptomatic infections in bean and N. benthamiana. Arch. Virol. 2020, 165, 1659–1665. [Google Scholar] [CrossRef]
  26. Alves-Freitas, D.M.; Pinheiro-Lima, B.; Faria, J.C.; Lacorte, C.; Ribeiro, S.G.; Melo, F.L. Double-stranded RNA high-throughput sequencing reveals a new cytorhabdovirus in a bean golden mosaic virus-resistant common bean transgenic line. Viruses 2019, 11, 90. [Google Scholar] [CrossRef] [Green Version]
  27. Pinheiro-Lima, B.; Pereira-Carvalho, R.C.; Alves-Freitas, D.M.T.; Kitajima, E.W.; Vidal, A.H.; Lacorte, C.; Godinho, M.T.; Fontenele, R.S.; Faria, J.C.; Abreu, E.F.M.; et al. Transmission of the Bean-associated cytorhabdovirus by the whitefly Bemisia tabaci MEAM1. Viruses 2020, 12, 1028. [Google Scholar] [CrossRef]
  28. Medina-Salguero, A.X.; Cornejo-Franco, J.F.; Grinstead, S.; Mollov, D.; Mowery, J.D.; Flores, F.; Quito-Avila, D.F. Sequencing, genome analysis and prevalence of a cytorhabdovirus discovered in Carica papaya. PLoS ONE 2019, 14, e0215798. [Google Scholar] [CrossRef]
  29. Okada, R.; Alcalá-Briseño, R.I.; Escalante, C.; Sabanadzovic, S.; Valverde, R.A. Genomic sequence of a novel endornavirus from Phaseolus vulgaris and occurrence in mixed infections with two other endornaviruses. Virus Res. 2018, 257, 63–67. [Google Scholar] [CrossRef]
  30. Chiquito-Almanza, E.; Acosta-Gallegos, J.A.; García-Álvarez, N.C.; Cuellar, W.; Martínez-Martínez, T.O.; Anaya-López, J.L. Detection of virus damaging the dry bean crop in Nayarit, Mexico. J. Chem. Biol. Physical. Sci. 2014, 4, 48–55. [Google Scholar]
  31. Narayanasamy, P. Viral and viroid pathogens. In Microbial Plants Pathogens-Detection and Disease Diagnosis; Springer: Dordrecht, The Netherlands, 2011; Volume 3, pp. 12–13. [Google Scholar]
  32. Voysest, O. Resultados del Primer Vivero Internacional de Rendimiento y Adaptación de Frijol (Phaseolus vulgaris L.) IBYAN 1976. In Frijol Arbustivo; Centro Internacional de Agricultura Tropical (CIAT): Cali, Colombia, 1979; pp. 157–159. [Google Scholar]
  33. Voysest, O. Vivero Internacional de Rendimiento y Adaptación de Frijol (Phaseolus vulgaris L.) IBYAN 1977. In Frijol de Grano Negro, Frijol de Grano de Diversos Colores; Centro Internacional de Agricultura Tropical (CIAT): Cali, Colombia, 1980; pp. 38–41. [Google Scholar]
  34. Voysest, O. Vivero Internacional de Rendimiento y Adaptación de Frijol (Phaseolus vulgaris L.) IBYAN 1978. In Frijol Arbustivo con Grano Negro, con Grano de Diversos Colores; Centro Internacional de Agricultura Tropical (CIAT): Cali, Colombia, 1981; pp. 30–47. [Google Scholar]
  35. Ortiz-Catón, M.; Median-Torres, R.; Valdivia-Bernal, R.; Ortiz-Catón, A.; Alvarado-Casillas, S.; Rodríguez-Blanco, J.R. Mosquitas blancas plaga primaria de hortalizas en Nayarit. Rev. Fuente 2010, 2, 31–40. [Google Scholar]
  36. Fraile, A.; McLeish, M.J.; Pagán, I.; González-Jara, P.; Piñero, D.; García-Arenal, F. Environmental heterogeneity and the evolution of plant-virus interactions: Viruses in wild pepper populations. Virus Res. 2017, 241, 68–76. [Google Scholar] [CrossRef] [PubMed]
  37. Sonnante, G.; Stockton, T.; Nodari, R.O.; Becerra Velásquez, V.L.; Gepts, P. Evolution of genetic diversity during the domestication of common bean (Phaseolus vulgaris L.). Theor. Appl. Genet. 1994, 89, 629–635. [Google Scholar] [CrossRef] [PubMed]
  38. Schmutz, J.; McClean, P.E.; Mamidi, S.; Wu, G.A.; Cannon, S.B.; Grimwood, J.; Jenkins, J.; Shu, S.; Song, Q.; Chavarro, C.; et al. A reference genome for common bean and genome-wide analysis of dual domestications. Nat. Genet. 2014, 46, 707–713. [Google Scholar] [CrossRef] [Green Version]
  39. Bernardo, P.; Charles-Dominique, T.; Barakat, M.; Ortet, P.; Fernandez, E.; Filloux, D.; Hartnady, P.; Rebelo, T.A.; Cousins, S.R.; Mesleard, F.; et al. Geometagenomics illuminates the impact of agriculture on the distribution and prevalence of plant viruses at the ecosystem scale. ISME J. 2018, 12, 173–184. [Google Scholar] [CrossRef] [Green Version]
  40. Chiquito-Almanza, E.; Caballero-Pérez, J.; Guevara-Olvera, L.; Acosta-García, G.; Pérez-Pérez, M.C.I.; Acosta-Gallegos, J.A.; Anaya-López, J.L. First report of Cowpea mild mottle virus infecting cultivated and wild Phaseolus in the Central-Western region of Mexico. Plant Dis. 2018, 102, 1047. [Google Scholar] [CrossRef]
  41. Jones, R.A.; Salam, M.U.; Maling, T.J.; Diggle, A.J.; Thackray, D.J. Principles of predicting plant virus disease epidemics. Annu. Rev. Phytopathol. 2010, 48, 179–203. [Google Scholar] [CrossRef] [PubMed]
  42. Rojas, M.R.; Gilbertson, R.L. Emerging plant viruses: A diversity of mechanisms and opportunities. In Plant Virus Evolution, 1st ed.; Roossinck, M.J., Ed.; Springer: Berlin/Heidelberg, Germany, 2008; pp. 27–51. [Google Scholar]
  43. Elena, S.F.; Fraile, A.; García-Arenal, F. Evolution and emergence of plant viruses. Adv. Virus Res. 2014, 88, 161–191. [Google Scholar] [CrossRef] [Green Version]
  44. Morales, F.J.; Anderson, P.K. The emergence and dissemination of whitefly-transmitted geminiviruses in Latin America. Arch. Virol. 2001, 146, 415–441. [Google Scholar] [CrossRef]
  45. Freitas-Astua, J.; Dietzgen, R.G.; Walker, P.J.; Blasdell, K.R.; Breyta, R.; Fooks, A.R. Create Twelve New Species in the Genus Cytorhabdovirus, Family Rhabdoviridae. Approved ICTV Proposal. 2019. Available online: https://ictv.global/ICTV/proposals/2019.030M.zip (accessed on 15 October 2020).
  46. Valverde, R.A.; Khalifa, M.E.; Okada, R.; Fukuhara, T.; Sabanadzovic, S. ICTV virus taxonomy profile: Endornaviridae. J. Gen. Virol. 2019, 100, 1204–1205. [Google Scholar] [CrossRef]
  47. Fukuhara, T. Endornaviruses: Persistent dsRNA viruses with symbiotic properties in diverse eukaryotes. Virus Genes 2019, 55, 165–173. [Google Scholar] [CrossRef] [PubMed]
  48. Okada, R.; Yong, C.K.; Valverde, R.A.; Sabanadzovic, S.; Aoki, N.; Hotate, S.; Kiyota, E.; Moriyama, H.; Fukuhara, T. Molecular characterization of two evolutionally distinct endornaviruses coinfecting common bean (Phaseolus vulgaris). J. Gen. Virol. 2013, 94, 220–229. [Google Scholar] [CrossRef] [Green Version]
  49. Khankhum, S.; Valverde, R.A. Physiological traits of endornavirus-infected and endornavirus-free common bean (Phaseolus vulgaris) cv Black Turtle Soup. Arch. Virol. 2018, 163, 1051–1056. [Google Scholar] [CrossRef]
  50. Morales, F.L.; Castaño, M. Enfermedades Virales del Frijol Común en América Latina; Centro Internacional de Agricultura Tropical (CIAT): Cali, Colombia, 2008; pp. 9–44. [Google Scholar]
  51. Morales, F.J.; Castaño, M. Seed transmission characteristics of selected bean common mosaic virus strains in differential bean cultivars. Plant Dis. 1987, 71, 51–53. [Google Scholar] [CrossRef]
  52. Boland, G.J.; Melzer, M.S.; Hopkin, A.; Higgins, V.; Nassuth, A. Climate change and plant diseases in Ontario. Can. J. Plant Pathol. 2004, 26, 335–350. [Google Scholar] [CrossRef]
  53. Beebe, S.; Ramirez, J.; Jarvis, A.; Rao, E.M.; Mosquera, G.; Bueno, J.M.; Blair, M.W. Genetic improvement of common beans and the challenges of climate change. In Crop Adaptation to Climate Change; Yadav, S.S., Redden, R., Hatfield, J.L., Lotze-Campen, H., Hall, A.J.W., Eds.; Wiley-Blackwell: Oxford, UK, 2011; pp. 356–369. [Google Scholar]
  54. Sastry, K.S. Seed-Borne Plant Virus Diseases; Springer: New Delhi, India, 2013; p. 88. [Google Scholar]
  55. Morales, F.J.; Bos, L. Bean Common Mosaic Virus. Descriptions of Plant Viruses, 337. Association of Applied Biologists, UK. 1998. Available online: http://www.dpvweb.net/dpv/showadpv.php?dpvno=337 (accessed on 28 May 2020).
  56. Sastry, K.S. Plant Virus and Viroid Diseases in the Tropics: Introduction of Plant Viruses and Sub-Viral Agents, Classification, Assessment of Loss, Transmission and Diagnosis; Springer: Dordrecht, The Netherlands, 2013; p. 110. [Google Scholar]
  57. Lamas, N.S.; Matos, V.O.R.L.; Alves-Freitas, D.M.T.; Melo, F.L.; Costa, A.F.; Faria, J.C.; Ribeiro, S.G. Occurrence of Cowpea mild mottle virus in common bean and associated weeds in Northeastern Brazil. Plant Dis. 2017, 101, 1828. [Google Scholar] [CrossRef]
  58. Zanardo, L.G.; Carvalho, C.M. Cowpea mild mottle virus (Carlavirus, Betaflexiviridae): A review. Trop. Plant Pathol. 2017, 42, 417–430. [Google Scholar] [CrossRef]
  59. Marubayashi, J.M.; Yuki, V.A.; Wutke, E.B. Transmission of Cowpea mild mottle virus by Bemisia tabaci biotype B from plants of beans and soybean. Summa Phytopathol. 2010, 36, 158–160. [Google Scholar] [CrossRef] [Green Version]
  60. Zanardo, L.G.; Silva, F.N.; Bicalho, A.A.C.; Urquiza, G.P.C.; Lima, A.T.M.; Almeida, A.M.R.; Zerbini, F.M.; Carvalho, C.M. Molecular and biological characterization of Cowpea mild mottle virus isolates infecting soybean in Brazil and evidence of recombination. Plant Pathol. 2014, 63, 456–465. [Google Scholar] [CrossRef]
  61. European and Mediterranean Plant Protection Organization (EPPO). Cowpea Mild Motte Virus. Categorization. Available online: https://gd.eppo.int/taxon/CPMMV0/categorization (accessed on 10 August 2020).
  62. Brown, J.K.; Verle-Rodrigues, J.C. Recovery Plan. Cowpea Mild Mottle Virus. 2014. Available online: https://www.ars.usda.gov/ARSUserFiles/00000000/opmp/CPMMV%20Recovery%20Plan.pdf (accessed on 10 August 2020).
  63. Khankhum, S.; Valverde, R.A.; Pastor-Corrales, M.; Osorno, J.M.; Sabanadzovic, S. Two endornaviruses show differential infection patterns between gene pools of Phaseolus vulgaris. Arch. Virol. 2015, 160, 1131–1137. [Google Scholar] [CrossRef]
  64. Staginnus, C.; Richert-Poggeler, K.R. Endogenous pararetroviruses: Two-faced travelers in the plant genome. Trends Plant Sci. 2006, 11, 485–491. [Google Scholar] [CrossRef] [PubMed]
Viruses 13 01153 g001 550
Figure 1. Small RNA size distribution of reads mapped and unmapped to host. Average number of reads of mapped and unmapped to host are shown in orange and blue lines, respectively.
Figure 1. Small RNA size distribution of reads mapped and unmapped to host. Average number of reads of mapped and unmapped to host are shown in orange and blue lines, respectively.
Viruses 13 01153 g001
Viruses 13 01153 g002 550
Figure 2. Virus distribution infecting Phaseolus spp. in the Mesoamerican center of domestication of common bean (Jalisco and Guanajuato States) and Nayarit, Mexico. The colors indicate the corresponding viral species. Circles and triangles correspond to the viruses identified in domesticated and wild Phaseolus spp., respectively. Geographic coordinates are indicated in Table S1.
Figure 2. Virus distribution infecting Phaseolus spp. in the Mesoamerican center of domestication of common bean (Jalisco and Guanajuato States) and Nayarit, Mexico. The colors indicate the corresponding viral species. Circles and triangles correspond to the viruses identified in domesticated and wild Phaseolus spp., respectively. Geographic coordinates are indicated in Table S1.
Viruses 13 01153 g002
Viruses 13 01153 g003 550
Figure 3. Virus detection, mixed infections, and prevalence of viruses in domesticated and wild Phaseolus spp. samples from the Mesoamerican center of domestication of common bean (Jalisco and Guanajuato States) and Nayarit, Mexico. (A) Genus and number of viral species detected in domesticated and wild bean in the MCD and Nayarit. (B) The classification of viral species detected by their mechanism of transmission, number of samples, and host (domesticated and wild bean). (C) A summary of samples organized by the number of detected viruses per sample and region. (D) Virus prevalence associated with domesticated (Dom.) and wild bean (Wild) in the MCD and Nayarit. Significant differences in virus prevalence are indicated by ** = p-value < 0.01.
Figure 3. Virus detection, mixed infections, and prevalence of viruses in domesticated and wild Phaseolus spp. samples from the Mesoamerican center of domestication of common bean (Jalisco and Guanajuato States) and Nayarit, Mexico. (A) Genus and number of viral species detected in domesticated and wild bean in the MCD and Nayarit. (B) The classification of viral species detected by their mechanism of transmission, number of samples, and host (domesticated and wild bean). (C) A summary of samples organized by the number of detected viruses per sample and region. (D) Virus prevalence associated with domesticated (Dom.) and wild bean (Wild) in the MCD and Nayarit. Significant differences in virus prevalence are indicated by ** = p-value < 0.01.
Viruses 13 01153 g003
Viruses 13 01153 g004 550
Figure 4. Viral distribution patterns in domesticated and wild bean samples from the Mesoamerican center of domestication of common beans (Guanajuato and Jalisco States), and Nayarit, Mexico. On the left side, the infection profiles are clustered and colored according to the region of provenance. Coverage % of reconstructed genome upon genomes used as reference is shown. DB: domesticated bean; WB: wild bean.
Figure 4. Viral distribution patterns in domesticated and wild bean samples from the Mesoamerican center of domestication of common beans (Guanajuato and Jalisco States), and Nayarit, Mexico. On the left side, the infection profiles are clustered and colored according to the region of provenance. Coverage % of reconstructed genome upon genomes used as reference is shown. DB: domesticated bean; WB: wild bean.
Viruses 13 01153 g004
Viruses 13 01153 g005 550
Figure 5. Phylogenetic tree of putative new alphaendornavirus PlepEV and PluEV: (A) Alphaendornavirus clade and (B) Betaendornavirus clade. Phylogenetic trees were constructed using an alignment of the amino acid sequences of the RdRp domains with the Maximum-Likelihood (LG + G + I model) with 1000 bootstraps. Grapevine leafroll-associated virus 1 (genus Ampelovirus, family Closteroviridae) was included as an outgroup. Nodes with less than 70% bootstrap values are collapsed. PlepEV and PluEV are underlined.
Figure 5. Phylogenetic tree of putative new alphaendornavirus PlepEV and PluEV: (A) Alphaendornavirus clade and (B) Betaendornavirus clade. Phylogenetic trees were constructed using an alignment of the amino acid sequences of the RdRp domains with the Maximum-Likelihood (LG + G + I model) with 1000 bootstraps. Grapevine leafroll-associated virus 1 (genus Ampelovirus, family Closteroviridae) was included as an outgroup. Nodes with less than 70% bootstrap values are collapsed. PlepEV and PluEV are underlined.
Viruses 13 01153 g005
Table
Table 1. Distribution, numbers of assembled genomes, and transmission of viruses infecting Phaseolus spp. in the MCD and Nayarit, Mexico.
Table 1. Distribution, numbers of assembled genomes, and transmission of viruses infecting Phaseolus spp. in the MCD and Nayarit, Mexico.
Virus NameAcronymGenusFamilyGenome *MCDNay.Genomes
Assembled §		Transmission
Gto.Jal.SeedVectors
Phaseolus vulgaris alphaendornavirus 1PvEV-1AlphaendornavirusEndornaviridaedsRNA++ 8+U
Phaseolus vulgaris alphaendornavirus 2PvEV-2AlphaendornavirusEndornaviridaedsRNA++ 5+U
Bean common mosaic virusBCMVPotyvirusPotyviridaessRNA (+)+++18+Aphids
Bean common mosaic necrosis virusBCMNVPotyvirusPotyviridaessRNA (+) ++10+Aphids
Peanut mottle virusPeMoVPotyvirusPotyviridaessRNA (+) + 1+Aphids
Cowpea mild mottle virusCPMMVCarlavirusBetaflexiviridaessRNA (+) +1+B. tabaci
Bean-associated cytorhabdovirusBaCVCytorhabdovirusRhabdoviridaessRNA (-) +1UB. tabaci
Bean golden yellow mosaic virus DNA ABGYMV DNA ABegomovirusGeminiviridaessDNA +2−B. tabaci
Bean golden yellow mosaic virus DNA BBGYMV DNA BBegomovirusGeminiviridaessDNA +2−B. tabaci
Bean latent virus DNA ABLV DNA ABegomovirusGeminiviridaessDNA +1UU
Bean latent virus DNA BBLV DNA BBegomovirusGeminiviridaessDNA +1UU
Phaseolus leptostachyus alphaendornavirusPlepEVPutative Alphaendornavirus EndornaviridaedsRNA+ UU
Phaseolus lunatus alphaendornavirusPluEVPutative Alphaendornavirus EndornaviridaedsRNA + UU
Badnavirus-like sequencePvBV-CG17Putative Badnavirus ¥CaulimoviridaedsDNA+ UU
Badnavirus-like sequencePcBV-SJ10-3Putative Badnavirus ¥CaulimoviridaedsDNA + UU
Badnavirus-like sequencePluBV-SN35Putative Badnavirus ¥CaulimoviridaedsDNA + UU
Pararetrovirus-like sequencePcPRV-SJ10-3Not assigned ¥CaulimoviridaedsDNA + UU
Pararetrovirus-like sequencePluPRV-SN35Not assigned ¥CaulimoviridaedsDNA + UU
: Putative 1. ¥: Putative 2. * Baltimore class d (W) Domestic. §: Number of complete or near-complete genomes assembled. MCD: Mesoamerican Center of Domestication. Gto.: Guanajuato; Jal.: Jalisco; Nay.: Nayarit. U: Unknown.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Chiquito-Almanza, E.; Caballero-Pérez, J.; Acosta-Gallegos, J.A.; Montero-Tavera, V.; Mariscal-Amaro, L.A.; Anaya-López, J.L. Diversity and Distribution of Viruses Infecting Wild and Domesticated Phaseolus spp. in the Mesoamerican Center of Domestication. Viruses 2021, 13, 1153. https://doi.org/10.3390/v13061153

AMA Style

Chiquito-Almanza E, Caballero-Pérez J, Acosta-Gallegos JA, Montero-Tavera V, Mariscal-Amaro LA, Anaya-López JL. Diversity and Distribution of Viruses Infecting Wild and Domesticated Phaseolus spp. in the Mesoamerican Center of Domestication. Viruses. 2021; 13(6):1153. https://doi.org/10.3390/v13061153

Chicago/Turabian Style

Chiquito-Almanza, Elizabeth, Juan Caballero-Pérez, Jorge A. Acosta-Gallegos, Victor Montero-Tavera, Luis Antonio Mariscal-Amaro, and José Luis Anaya-López. 2021. "Diversity and Distribution of Viruses Infecting Wild and Domesticated Phaseolus spp. in the Mesoamerican Center of Domestication" Viruses 13, no. 6: 1153. https://doi.org/10.3390/v13061153

APA Style

Chiquito-Almanza, E., Caballero-Pérez, J., Acosta-Gallegos, J. A., Montero-Tavera, V., Mariscal-Amaro, L. A., & Anaya-López, J. L. (2021). Diversity and Distribution of Viruses Infecting Wild and Domesticated Phaseolus spp. in the Mesoamerican Center of Domestication. Viruses, 13(6), 1153. https://doi.org/10.3390/v13061153

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop