Complex recombination patterns arising during geminivirus coinfections preserve and demarcate biologically important intra-genome interaction networks

Virome analysis of field-collected chilli samples reveals diverse viruses

BACKGROUND

Chilli (Capsicum annuum L.), an important spice crop, is susceptible to diverse viral infections. Traditional detection methods including PCR and its variants had difficulty in identifying the complete spectrum of viruses, especially in mixed infections. High-throughput sequencing (HTS) has emerged as a successful tool for comprehensive virome analyses, enabling the identification of the known and novel viruses in the infected samples. Using HTS, we investigated virome analyses to identify known and novel viruses in chilli.

METHODS

In 2021-22, 19 leaf samples were collected from chili plants in farmer fields in Karnataka, India, showing symptoms such as leaf curling, vein banding, mosaic, mottling, filiform, leathery, dull-colored, and bunchy leaves. Total RNA was extracted, pooled at equimolar concentrations, and subjected to virome profiling. rRNA-depleted RNA was used to prepare mRNA and sRNA libraries, which were sequenced on the Illumina NovaSeq 6000 platform. Bioinformatics tools were used to analyze the sequencing data and identify plant viruses.

RESULTS

Viral disease incidences varied from 26.6 to 47.5% in the farmer fields surveyed. Virome analyses revealed complete/ near-complete genomes of six different viruses: chilli leaf curl virus (ChiLCV), cucumber mosaic virus (CMV), groundnut bud necrosis orthotospovirus (GBNV), pepper cryptic virus-2 (PCV-2), pepper vein yellows virus (PeVYV) and bell pepper alphaendornavirus (BPEV). The viral copy number of ChiLCV was found to be the highest (45.36%) and had the least mutational frequency (SNPs) and was also associated with five satellites. Recombination breakpoints were observed in ChiLCV (coat protein and AC4 regions), CMV RNA2 (2a protein) and PeVYV (P0, P3 and P5 proteins), indicating their origins from intra- and interspecific recombination events. Identified viruses in the pooled RNA sample were confirmed by PCR. Further, novel loop-mediated isothermal amplification (LAMP) diagnostic assays were developed for diagnosing the identified viruses for future use. Among the six viruses identified in chilli, PeVYV and BPEV are the first reports from India.

CONCLUSIONS

This study presents the first virome profiling of chili using HTS and identified known and previously unreported viruses in farmer fields of Karnataka, India. Understanding viral diversity provides insights for developing diagnostic tools and effective management strategies.

Macroptilium yellow vein virus (Begomovirus macroptilivenae) associated with Macroptilium plants is geographically structured into divergent subpopulations

The genus Begomovirus (family Geminiviridae) includes plant viruses that have a negative impact on the production of several economically important crops worldwide. Begomoviruses infect a wide range of non-cultivated plants, which may serve as alternative hosts or inoculum sources for epidemics in crops. In this study, we evaluated the genetic variability and population structure of Macroptilium yellow vein virus (MacYVV; Begomovirus macroptilivenae), associated with non-cultivated Macroptilium spp. in the state of Alagoas (AL), northeastern Brazil. Complete DNA-A sequences (n = 29) were obtained from symptomatic plants collected in different municipalities of AL during 2019 and 2020. Pairwise sequence comparisons indicated a low degree of genetic variability, with 93.9-100% identity to 18 previously characterized sequences. To better understand the population structuring of MacYVV in AL, the complete DNA-A data set (n = 47) was analyzed. Although the DAPC and STRUCTURE approaches yielded slightly divergent numbers of putative MacYVV subpopulations, the results are consistent with the phylogenetic relationships among MacYVV isolates, where five well-supported subgroups correlated with geographic location were observed. The low variability of the MacYVV population could be related to the small number of recombination events detected and be reflective of a long-term virus-host relationship in which the virus is very well adapted to the host.

Novel genetic variants of banana streak MY virus and banana streak IM virus naturally infecting banana in Northeast India

Divergent banana streak viruses (BSV) were characterized from banana plants exhibiting diverse symptoms in the Northeast region (NER) of India. Using rolling circle amplification (RCA), the complete genome sequences of seven episomal banana streak MY virus (BSMYV) isolates, including two novel variants, and two new banana streak IM virus (BSIMV) isolates were characterized. The novel BSMYV genetic variants were associated with conspicuous necrosis on newly emerged leaves, peduncle distortion, pseudostem internal necrosis, in addition to common streak symptoms. For complete genome nucleotide sequences, BSMYV-IN4 and IN5 shared 77-79% identity with other BSMYVs, while BSMYV-IN7 and IN8 exhibited identities of 77-97%. This study reports for the first time, the complete genomes of two banana streak IM virus (BSIMV-IN1 and -IN2) infecting triploid banana hybrids exhibiting leaf distortion, stunted rosette-like growth, and necrosis, sharing 87% sequence identity with reference BSIMV genome (GenBank accession no. HQ593112). Phylogenetic inference based on complete genomes revealed the distinct and congruent placement of BSMYV-IN4 and IN5 within the BSMYV cluster. Pairwise sequence comparisons of the conserved RT/RNase H nucleotide (nt) sequences revealed that the BSMYV-IN7 and IN4 isolates showed 85% and 97% identity to BSMYV (AY805074), respectively, which shared highest nt identity with BSMYV-IN6, IN9, and IN10, at 100%. The RT/RNase H nt sequences of BSIMV-IN1 and IN2 had 98% identity with the BSIMV (HQ593112), but were characterized as novel variants of BSIMV based on complete genomes. An analysis of relative synonymous codon usage (RSCU) pattern in the ORFIII polyprotein of BSMYV and BSIMV isolates revealed AGA and AGG (arginine) as the most frequently overrepresented codons (>1.5), evolutionary conserved in the genome of both species. A total of 14 recombination events were detected among the 36 BSV genomes, with recombination breakpoints mainly located in the ORFI, III, and IGR genomic regions. A novel phylogenetic cluster, comprised of BSMYV-IN4 and IN5 within the clade I was probably derived from heterologous recombination between parents resembling banana streak VN virus (BSVNV; AY750155) and banana streak GF virus (BSGFV; KJ013507) isolates. The present study conclusively reports the infection of genetically and symptomatically distinct variants of BSMYV and BSIMV infecting banana hybrids in NER India.

Phylogeographic analysis of Begomovirus coat and replication-associated proteins

Begomoviruses are globally distributed plant pathogens that significantly limit crop production. These viruses are traditionally described according to phylogeographic distribution and categorized into two groups: begomoviruses from the Africa, Asia, Europe and Oceania (AAEO) region and begomoviruses from the Americas. Monopartite begomoviruses are more common in the AAEO region, while bipartite viruses predominate in the Americas, where the begomoviruses lack the V2/AV2 gene involved in inter-cellular movement and RNA silencing suppression found in AAEO begomoviruses. While these features are generally accepted as lineage-defining, the number of known species has doubled due to sequence-based discovery since 2010. To re-evaluate the geographic groupings after the rapid expansion of the genus, we conducted phylogenetic analyses for begomovirus species representatives of the two longest and most conserved begomovirus proteins: the coat and replication-associated proteins. Both proteins still largely support the broad AAEO and Americas begomovirus groupings, except for sweet potato-infecting begomoviruses that form an independent, well-supported clade for their coat protein regardless of the region they were isolated from. Our analyses do not support more fine-scaled phylogeographic groupings. Monopartite and bipartite genome organizations are broadly interchanged throughout the phylogenies, and the absence of the V2/AV2 gene is highly reflective of the split between Americas and AAEO begomoviruses. We observe significant evidence of recombination within the Americas and within the AAEO region but rarely between the regions. We speculate that increased globalization of agricultural trade, the invasion of polyphagous whitefly vector biotypes and recombination will blur begomovirus phylogeographic delineations in the future.

Development, Design, and Application of Efficient siRNAs Against Cotton Leaf Curl Virus-Betasatellite Complex to Mediate Resistance Against Cotton Leaf Curl Disease

Cotton leaf curl disease (CLCuD), caused by the Cotton leaf curl virus, is one of the most irrepressible diseases in cotton due to high recombination in the virus. RNA interference (RNAi) is widely used as a biotechnological approach for sequence-specific gene silencing guided by small interfering RNAs (siRNAs) to generate resistance against viruses. The success of RNAi depends upon the fact that the target site of the designed siRNA must be conserved even if the genome undergoes recombination. Thus, the present study designs the most efficient siRNA against the conserved sites of the Cotton leaf curl Multan virus (CLCuMuV) and the Cotton leaf curl Multan betasatellite (CLCuMB). From an initial prediction of 9 and 7 siRNAs against CLCuMuV and CLCuMB, respectively, the final selection was made for 2 and 1 siRNA based on parameters such as no off-targets, good GC content, high validity score, and targeting coding region. The target sites of siRNA were observed to lie in the AC3 and an overlapping region of AC2-AC1 of CLCuMuV and βC1 of CLCuMB; all target sites showed a highly conserved nature in recombination analysis. Docking the designed siRNAs with the Argonaute-2 protein of Gossypium hirsutum showed stable binding. Finally, BLASTn of siRNA-target positions in genomes of other BGVs indicated the suitability of designed siRNAs against a broad range of BGVs. The designed siRNAs of the present study could help gain complete control over the virus, though experimental validation is highly required to suggest predicted siRNAs for CLCuD resistance.

Supplementary Information

The online version contains supplementary material available at 10.1007/s12088-024-01191-z.

Reassortments in single-stranded DNA multipartite viruses: Confronting expectations based on molecular constraints with field observations

Single-stranded DNA multipartite viruses, which mostly consist of members of the genus Begomovirus, family Geminiviridae, and all members of the family Nanoviridae, partly resolve the cost of genomic integrity maintenance through two remarkable capacities. They are able to systemically infect a host even when their genomic segments are not together in the same host cell, and these segments can be separately transmitted by insect vectors from host to host. These capacities potentially allow such viruses to reassort at a much larger spatial scale, since reassortants could arise from parental genotypes that do not co-infect the same cell or even the same host. To assess the limitations affecting reassortment and their implications in genome integrity maintenance, the objective of this review is to identify putative molecular constraints influencing reassorted segments throughout the infection cycle and to confront expectations based on these constraints with empirical observations. Trans-replication of the reassorted segments emerges as the major constraint, while encapsidation, viral movement, and transmission compatibilities appear more permissive. Confronting the available molecular data and the resulting predictions on reassortments to field population surveys reveals notable discrepancies, particularly a surprising rarity of interspecific natural reassortments within the Nanoviridae family. These apparent discrepancies unveil important knowledge gaps in the biology of ssDNA multipartite viruses and call for further investigation on the role of reassortment in their biology.

Genetic Diversity, Evolutionary Dynamics, and Ongoing Spread of Pedilanthus Leaf Curl Virus

Pedilanthus leaf curl virus (PeLCV) is a monopartite begomovirus (family Geminiviridae) discovered just a few decades ago. Since then, it has become a widely encountered virus, with reports from ca. 25 plant species across Pakistan and India, indicative of its notable evolutionary success. Viruses mutate at such a swift rate that their ecological and evolutionary behaviors are inextricably linked, and all of these behaviors are imprinted on their genomes as genetic diversity. So, all these imprints can be mapped by computational methods. This study was designed to map the sequence variation dynamics, genetic heterogeneity, regional diversity, phylogeny, and recombination events imprinted on the PeLCV genome. Phylogenetic and network analysis grouped the full-length genome sequences of 52 PeLCV isolates into 7 major clades, displaying some regional delineation but lacking host-specific demarcation. The progenitor of PeLCV was found to have originated in Multan, Pakistan, in 1977, from where it spread concurrently to India and various regions of Pakistan. A high proportion of recombination events, distributed unevenly throughout the genome and involving both inter- and intraspecies recombinants, were inferred. The findings of this study highlight that the PeLCV population is expanding under a high degree of genetic diversity (π = 0.073%), a high rate of mean nucleotide substitution (1.54 × 10-3), demographic selection, and a high rate of recombination. This sets PeLCV apart as a distinctive begomovirus among other begomoviruses. These factors could further exacerbate the PeLCV divergence and adaptation to new hosts. The insights of this study that pinpoint the emergence of PeLCV are outlined.

Geminiviral C4/AC4 proteins: An emerging component of the viral arsenal against plant defence

Virus infection triggers a plethora of defence reactions in plants to incapacitate the intruder. Viruses, in turn, have added additional functions to their genes so that they acquire capabilities to neutralize the above defence reactions. In plant-infecting viruses, the family Geminiviridae comprises members, majority of whom encode 6-8 genes in their small single-stranded DNA genomes. Of the above genes, one which shows the most variability in its amino acid sequence is the C4/AC4. Recent studies have uncovered evidence, which point towards a wide repertoire of functions performed by C4/AC4 revealing its role as a major player in suppressing plant defence. This review summarizes the various plant defence mechanisms against viruses and highlights how C4/AC4 has evolved to counter most of them.

The Role of Extensive Recombination in the Evolution of Geminiviruses

Mutation, recombination and pseudo-recombination are the major forces driving the evolution of viruses by the generation of variants upon which natural selection, genetic drift and gene flow can act to shape the genetic structure of viral populations. Recombination between related virus genomes co-infecting the same cell usually occurs via template swapping during the replication process and produces a chimeric genome. The family Geminiviridae shows the highest evolutionary success among plant virus families, and the common presence of recombination signatures in their genomes reveals a key role in their evolution. This review describes the general characteristics of members of the family Geminiviridae and associated DNA satellites, as well as the extensive occurrence of recombination at all taxonomic levels, from strain to family. The review also presents an overview of the recombination patterns observed in nature that provide some clues regarding the mechanisms involved in the generation and emergence of recombinant genomes. Moreover, the results of experimental evolution studies that support some of the conclusions obtained in descriptive or in silico works are summarized. Finally, the review uses a number of case studies to illustrate those recombination events with evolutionary and pathological implications as well as recombination events in which DNA satellites are involved.

A Trait-Based Approach to Predicting Viral Host-Range Evolvability

Predicting the evolution of virus host range has proven to be extremely difficult, in part because of the sheer diversity of viruses, each with unique biology and ecological interactions. We have not solved this problem, but to make the problem more tractable, we narrowed our focus to three traits intrinsic to all viruses that may play a role in host-range evolvability: mutation rate, recombination rate, and phenotypic heterogeneity. Although each trait should increase evolvability, they cannot do so unbounded because fitness trade-offs limit the ability of all three traits to maximize evolvability. By examining these constraints, we can begin to identify groups of viruses with suites of traits that make them especially concerning, as well as ecological and environmental conditions that might push evolution toward accelerating host-range expansion.

How To Be a Successful Monopartite Begomovirus in a Bipartite-Dominated World: Emergence and Spread of Tomato Mottle Leaf Curl Virus in Brazil

Begomoviruses are members of the family Geminiviridae, a large and diverse group of plant viruses characterized by a small circular single-stranded DNA genome encapsidated in twinned quasi-icosahedral virions. Cultivated tomato (Solanum lycopersicum L.) is particularly susceptible and is infected by >100 bipartite and monopartite begomoviruses worldwide. In Brazil, 25 tomato-infecting begomoviruses have been described, most of which are bipartite. Tomato mottle leaf curl virus (ToMoLCV) is one of the most important of these and was first described in the late 1990s but has not been fully characterized. Here, we show that ToMoLCV is a monopartite begomovirus with a genomic DNA similar in size and genome organization to those of DNA-A components of New World (NW) begomoviruses. Tomato plants agroinoculated with the cloned ToMoLCV genomic DNA developed typical tomato mottle leaf curl disease symptoms, thereby fulfilling Koch's postulates and confirming the monopartite nature of the ToMoLCV genome. We further show that ToMoLCV is transmitted by whiteflies, but not mechanically. Phylogenetic analyses placed ToMoLCV in a distinct and strongly supported clade with other begomoviruses from northeastern Brazil, designated the ToMoLCV lineage. Genetic analyses of the complete sequences of 87 ToMoLCV isolates revealed substantial genetic diversity, including five strain groups and seven subpopulations, consistent with a long evolutionary history. Phylogeographic models generated with partial or complete sequences predicted that the ToMoLCV emerged in northeastern Brazil >700 years ago, diversifying locally and then spreading widely in the country. Thus, ToMoLCV emerged well before the introduction of MEAM1 whiteflies, suggesting that the evolution of NW monopartite begomoviruses was facilitated by local whitefly populations and the highly susceptible tomato host. IMPORTANCE Worldwide, diseases of tomato caused by whitefly-transmitted geminiviruses (begomoviruses) cause substantial economic losses and a reliance on insecticides for management. Here, we describe the molecular and biological properties of tomato mottle leaf curl virus (ToMoLCV) from Brazil and establish that it is a NW monopartite begomovirus indigenous to northeastern Brazil. This answered a long-standing question regarding the genome of this virus, and it is part of an emerging group of these viruses in Latin America. This appears to be driven by widespread planting of the highly susceptible tomato and by local and exotic whiteflies. Our extensive phylogenetic studies placed ToMoLCV in a distinct strongly supported clade with other begomoviruses from northeastern Brazil and revealed new insights into the origin of Brazilian begomoviruses. The novel phylogeographic analysis indicated that ToMoLCV has had a long evolutionary history, emerging in northeastern Brazil >700 years ago. Finally, the tools used here (agroinoculation system and ToMoLCV-specific PCR test) and information on the biology of the virus (host range and whitefly transmission) will be useful in developing and implementing integrated pest management (IPM) programs targeting ToMoLCV.

Evolutionary Dynamics of Begomoviruses and Its Satellites Infecting Papaya in India

The genus Begomovirus represents a group of multipartite viruses that significantly damage many agricultural crops, including papaya, and influence overall production. Papaya leaf curl disease (PaLCD) caused by the complex begomovirus species has several important implications and substantial losses in papaya production in many developing countries, including India. The increase in the number of begomovirus species poses a continuous threat to the overall production of papaya. Here, we attempted to map the genomic variation, mutation, evolution rate, and recombination to know the disease complexity and successful adaptation of PaLCD in India. For this, we retrieved 44 DNA-A and 26 betasatellite sequences from GenBank reported from India. An uneven distribution of evolutionary divergence has been observed using the maximum-likelihood algorithm across the branch length. Although there were phylogenetic differences, we found high rates of nucleotide substitution mutation in both viral and sub-viral genome datasets. We demonstrated frequent recombination of begomovirus species, with a maximum in intra-species recombinants. Furthermore, our results showed a high degree of genetic variability, demographic selection, and mean substitution rate acting on the population, supporting the emergence of a diverse and purifying selection of viruses and associated betasatellites. Moreover, variation in the genetic composition of all begomovirus datasets revealed a predominance of nucleotide diversity principally driven by mutation, which might further accelerate the advent of new strains and species and their adaption to various hosts with unique pathogenicity. Therefore, the finding of genetic variation and selection emphases on factors that contribute to the universal spread and evolution of Begomovirus and this unanticipated diversity may also provide guidelines toward future evolutionary trend analyses and the development of wide-ranging disease control strategies for begomoviruses associated with PaLCD.

Role of Diversity and Recombination in the Emergence of Chilli Leaf Curl Virus

Chilli leaf curl virus (ChiLCV), (Genus Begomovirus, family Geminiviridae) and associated satellites pose a serious threat to chilli production, worldwide. This study highlights the factors accountable for genetic diversity, recombination, and evolution of ChiLCV, and associated chilli leaf curl alphasatellite (ChiLCA) and chilli leaf curl betasatellite (ChiLCB). Phylogenetic analysis of complete genome (DNA-A) sequences of 132 ChiLCV isolates from five countries downloaded from NCBI database clustered into three major clades and showed high population diversity. The dN/dS ratio and Tajima D value of all viral DNA-A and associated betasatellite showed selective control on evolutionary relationships. Negative values of neutrality tests indicated purified selection and an excess of low-frequency polymorphism. Nucleotide diversity (π) for C4 and Rep genes was higher than other genes of ChiLCV with an average value of π = 18.37 × 10⁻² and π = 17.52 × 10⁻² respectively. A high number of mutations were detected in TrAP and Rep genes, while ChiLCB has a greater number of mutations than ChiLCA. In addition, significant recombination breakpoints were detected in all regions of ChiLCV genome, ChiLCB and, ChiLCA. Our findings indicate that ChiLCV has the potential for rapid evolution and adaptation to a range of geographic conditions and could be adopted to infect a wide range of crops, including diverse chilli cultivars.

Banana bunchy top virus genetic diversity in Pakistan and association of diversity with recombination in its genomes

Banana Bunchy top virus (BBTV) is a multipartite circular single strand DNA virus that belongs to genus Babuvirus and family Nanoviridae. It causes significant crop losses worldwide and also in Pakistan. BBTV is present in Pakistan since 1988 however, till now only few (about twenty only) sequence of genomic components have been reported from the country. To have insights into current genetic diversity in Pakistan fifty-seven genomic components including five complete genomes (comprises of DNA-R, -U3, -S, -M, -C and -N components) were sequenced in this study. The genetic diversity analysis of populations from Pakistan showed that DNA-R is highly conserved followed by DNA-N, whereas DNA-U3 is highly diverse with the most diverse Common Region Stem-loop (CR-SL) in BBTV genome, a functional region, which previously been reported to have undergone recombination in Pakistani population. A Maximum Likelihood (ML) phylogenetic analysis of entire genomes of isolates by using sequence of all the components concatenated together with the reported genomes around the world revealed deeper insights about the origin of the disease in Pakistan. A comparison of the genetic diversity of Pakistani and entire BBTV populations around the world indicates that there exists a correlation between genetic diversity and recombination. Population genetics analysis indicated that the degree of selection pressure differs depending on the area and genomic component. A detailed analysis of recombination across various components and functional regions suggested that recombination is closely associated with the functional parts of BBTV genome showing high genetic diversity. Both genetic diversity and recombination analyses suggest that the CR-SL is a recombination hotspot in all BBTV genomes and among the six components DNA-U3 is the only recombined component that has extensively undergone inter and intragenomic recombination. Diversity analysis of recombinant regions results on average one and half fold increase and, in some cases up to four-fold increase due to recombination. These results suggest that recombination is significantly contributing to the genetic diversity of BBTV populations around the world.

Whitefly-Mediated Transmission and Subsequent Acquisition of Highly Similar and Naturally Occurring Tomato Yellow Leaf Curl Virus Variants

Begomoviruses are whitefly-transmitted viruses that infect many agricultural crops. Numerous reports exist on individual host plants harboring two or more begomoviruses. Mixed infection allows recombination events to occur among begomoviruses. However, very few studies have examined mixed infection of different isolates/variants/strains of a Begomovirus species in hosts. In this study, the frequency of mixed infection of tomato yellow leaf curl virus (TYLCV) variants in field-grown tomato was evaluated. At least 60% of symptomatic field samples were infected with more than one TYLCV variant. These variants differed by a few nucleotides and amino acids, resembling a quasispecies. Subsequently, in the greenhouse, single and mixed infection of two TYLCV variants (variant #2 and variant #4) that shared 99.5% nucleotide identity and differed by a few amino acids was examined. Plant-virus variant-whitefly interactions including transmission of one and/or two variants, variants' concentrations, competition between variants in inoculated tomato plants, and whitefly acquisition of one and/or two variants were assessed. Whiteflies transmitted both variants to tomato plants at similar frequencies; however, the accumulation of variant #4 was greater than that of variant #2 in tomato plants. Despite differences in variants' accumulation in inoculated tomato plants, whiteflies acquired variant #2 and variant #4 at similar frequencies. Also, whiteflies acquired greater amounts of TYLCV from singly infected plants than from mixed-infected plants. These results demonstrated that even highly similar TYLCV variants could differentially influence component (whitefly-variant-plant) interactions.

Gene Overlapping as a Modulator of Begomovirus Evolution

In RNA viruses, which have high mutation—and fast evolutionary— rates, gene overlapping (i.e., genomic regions that encode more than one protein) is a major factor controlling mutational load and therefore the virus evolvability. Although DNA viruses use host high-fidelity polymerases for their replication, and therefore should have lower mutation rates, it has been shown that some of them have evolutionary rates comparable to those of RNA viruses. Notably, these viruses have large proportions of their genes with at least one overlapping instance. Hence, gene overlapping could be a modulator of virus evolution beyond the RNA world. To test this hypothesis, we use the genus Begomovirus of plant viruses as a model. Through comparative genomic approaches, we show that terminal gene overlapping decreases the rate of virus evolution, which is associated with lower frequency of both synonymous and nonsynonymous mutations. In contrast, terminal overlapping has little effect on the pace of virus evolution. Overall, our analyses support a role for gene overlapping in the evolution of begomoviruses and provide novel information on the factors that shape their genetic diversity.

Interspecies Recombination Has Driven the Macroevolution of Cassava Mosaic Begomoviruses

Begomoviruses (family Geminiviridae, genus Begomovirus) significantly hamper crop production and threaten food security around the world. The frequent emergence of new begomovirus genotypes is facilitated by high mutation frequencies and the propensity to recombine and reassort. Homologous recombination has been especially implicated in the emergence of novel cassava mosaic begomovirus (CMB) genotypes, which cause cassava mosaic disease (CMD). Cassava (Manihot esculenta) is a staple food crop throughout Africa and an important industrial crop in Asia, two continents where production is severely constrained by CMD. The CMD species complex is comprised of 11 bipartite begomovirus species with ample distribution throughout Africa and the Indian subcontinent. While recombination is regarded as a frequent occurrence for CMBs, a revised, systematic assessment of recombination and its impact on CMB phylogeny is currently lacking. We assembled data sets of all publicly available, full-length DNA-A (n = 880) and DNA-B (n = 369) nucleotide sequences from the 11 recognized CMB species. Phylogenetic networks and complementary recombination detection methods revealed extensive recombination among the CMB sequences. Six out of the 11 species descended from unique interspecies recombination events. Estimates of recombination and mutation rates revealed that all species experience mutation more frequently than recombination, but measures of population divergence indicate that recombination is largely responsible for the genetic differences between species. Our results support that recombination has significantly impacted the CMB phylogeny and has driven speciation in the CMD species complex. IMPORTANCE Cassava mosaic disease (CMD) is a significant threat to cassava production throughout Africa and Asia. CMD is caused by a complex comprised of 11 recognized virus species exhibiting accelerated rates of evolution, driven by high frequencies of mutation and genetic exchange. Here, we present a systematic analysis of the contribution of genetic exchange to cassava mosaic virus species-level diversity. Most of these species emerged as a result of genetic exchange. This is the first study to report the significant impact of genetic exchange on speciation in a group of viruses.

Molecular phylogenetics and evolutionary analysis of a highly recombinant begomovirus, Cotton leaf curl Multan virus, and associated satellites

Cotton leaf curl Multan virus (CLCuMuV) and its associated satellites are a major part of the cotton leaf curl disease (CLCuD) caused by the begomovirus species complex. Despite the implementation of potential disease management strategies, the incessant resurgence of resistance-breaking variants of CLCuMuV imposes a continuous threat to cotton production. Here, we present a focused effort to map the geographical prevalence, genomic diversity, and molecular evolutionary endpoints that enhance disease complexity by facilitating the successful adaptation of CLCuMuV populations to the diversified ecosystems. Our results demonstrate that CLCuMuV populations are predominantly distributed in China, while the majority of alphasatellites and betasatellites exist in Pakistan. We demonstrate that together with frequent recombination, an uneven genetic variation mainly drives CLCuMuV and its satellite's virulence and evolvability. However, the pattern and distribution of recombination breakpoints greatly vary among viral and satellite sequences. The CLCuMuV, Cotton leaf curl Multan alphasatellite, and Cotton leaf curl Multan betasatellite populations arising from distinct regions exhibit high mutation rates. Although evolutionarily linked, these populations are independently evolving under strong purifying selection. These findings will facilitate to comprehensively understand the standing genetic variability and evolutionary patterns existing among CLCuMuV populations across major cotton-producing regions of the world.

Revealing the Complexity of Sweepovirus-Deltasatellite-Plant Host Interactions: Expanded Natural and Experimental Helper Virus Range and Effect Dependence on Virus-Host Combination

Sweepoviruses are begomoviruses (genus Begomovirus, family Geminiviridae) with ssDNA genomes infecting sweet potato and other species of the family Convolvulaceae. Deltasatellites (genus Deltasatellite, family Tolecusatellitidae) are small-size non-coding DNA satellites associated with begomoviruses. In this study, the genetic diversity of deltasatellites associated with sweepoviruses infecting Ipomoea indica plants was analyzed by further sampling the populations where the deltasatellite sweet potato leaf curl deltasatellite 1 (SPLCD1) was initially found, expanding the search to other geographical areas in southern continental Spain and the Canary Islands. The sweepoviruses present in the samples coinfected with deltasatellites were also fully characterized by sequencing in order to define the range of viruses that could act as helper viruses in nature. Additionally, experiments were performed to assess the ability of a number of geminivirids (the monopartite tomato leaf deformation virus and the bipartite NW begomovirus Sida golden yellow vein virus, the bipartite OW begomovirus tomato leaf curl New Delhi virus, and the curtovirus beet curly top virus) to transreplicate SPLCD1 in their natural plant hosts or the experimental host Nicotiana benthamiana. The results show that SPLCD1 can be transreplicated by all the geminivirids assayed in N. benthamiana and by tomato leaf curl New Delhi virus in zucchini. The presence of SPLCD1 did not affect the symptomatology caused by the helper viruses, and its effect on viral DNA accumulation depended on the helper virus-host plant combination.

Weed-infecting viruses in a tropical agroecosystem present different threats to crops and evolutionary histories

In the Caribbean Basin, malvaceous weeds commonly show striking golden/yellow mosaic symptoms. Leaf samples from Malachra sp. and Abutilon sp. plants with these symptoms were collected in Hispaniola from 2014 to 2020. PCR tests with degenerate primers revealed that all samples were infected with a bipartite begomovirus, and sequence analyses showed that Malachra sp. plants were infected with tobacco leaf curl Cuba virus (TbLCuCV), whereas the Abutilon sp. plants were infected with a new bipartite begomovirus, tentatively named Abutilon golden yellow mosaic virus (AbGYMV). Phylogenetic analyses showed that TbLCuCV and AbGYMV are distinct but closely related species, which are most closely related to bipartite begomoviruses infecting weeds in the Caribbean Basin. Infectious cloned DNA-A and DNA-B components were used to fulfilled Koch's postulates for these diseases of Malachra sp. and Abutilon sp. In host range studies, TbLCuCV also induced severe symptoms in Nicotiana benthamiana, tobacco and common bean plants; whereas AbGYMV induced few or no symptoms in plants of these species. Pseudorecombinants generated with the infectious clones of these viruses were highly infectious and induced severe symptoms in N. benthamiana and Malachra sp., and both viruses coinfected Malachra sp., and possibly facilitating virus evolution via recombination and pseudorecombination. Together, our results suggest that TbLCuCV primarily infects Malachra sp. in the Caribbean Basin, and occasionally spills over to infect and cause disease in crops; whereas AbGYMV is well-adapted to an Abutilon sp. in the Dominican Republic and has not been reported infecting crops.

New World Cactaceae Plants Harbor Diverse Geminiviruses

The family Cactaceae comprises a diverse group of typically succulent plants that are native to the American continent but have been introduced to nearly all other continents, predominantly for ornamental purposes. Despite their economic, cultural, and ecological importance, very little research has been conducted on the viral community that infects them. We previously identified a highly divergent geminivirus that is the first known to infect cacti. Recent research efforts in non-cultivated and asymptomatic plants have shown that the diversity of this viral family has been under-sampled. As a consequence, little is known about the effects and interactions of geminiviruses in many plants, such as cacti. With the objective to expand knowledge on the diversity of geminiviruses infecting cacti, we used previously acquired high-throughput sequencing results to search for viral sequences using BLASTx against a viral RefSeq protein database. We identified two additional sequences with similarity to geminiviruses, for which we designed abutting primers and recovered full-length genomes. From 42 cacti and five scale insects, we derived 42 complete genome sequences of a novel geminivirus species that we have tentatively named Opuntia virus 2 (OpV2) and 32 genomes of an Opuntia-infecting becurtovirus (which is a new strain of the spinach curly top Arizona virus species). Interspecies recombination analysis of the OpV2 group revealed several recombinant regions, in some cases spanning half of the genome. Phylogenetic analysis demonstrated that OpV2 is a novel geminivirus more closely related to viruses of the genus Curtovirus, which was further supported by the detection of three recombination events between curtoviruses and OpV2. Both OpV2 and Opuntia becurtoviruses were identified in mixed infections, which also included the previously characterized Opuntia virus 1. Viral quantification of the co-infected cactus plants compared with single infections did not show any clear trend in viral dynamics that might be associated with the mixed infections. Using experimental Rhizobium-mediated inoculations, we found that the initial accumulation of OpV2 is facilitated by co-infection with OpV1. This study shows that the diversity of geminiviruses that infect cacti is under-sampled and that cacti harbor diverse geminiviruses. The detection of the Opuntia becurtoviruses suggests spill-over events between viruses of cultivated species and native vegetation. The threat this poses to cacti needs to be further investigated.

Population Genomics of Plant Viruses: The Ecology and Evolution of Virus Emergence

The genomics era has revolutionized studies of adaptive evolution by monitoring large numbers of loci throughout the genomes of many individuals. Ideally, the investigation of emergence in plant viruses requires examining the population dynamics of both virus and host, their interactions with each other, with other organisms and the abiotic environment. Genetic mechanisms that affect demographic processes are now being studied with high-throughput technologies, traditional genetics methods, and new computational tools for big-data. In this review, we discuss the utility of these approaches to monitor and detect changes in virus populations within cells and individuals, and over wider areas across species and communities of ecosystems. The advent of genomics in virology has fostered a multidisciplinary approach to tackling disease risk. The ability to make sense of the information now generated in this integrated setting is by far the most substantial obstacle to the ultimate goal of plant virology to minimize the threats to food security posed by disease. To achieve this goal, it is imperative to understand and forecast how populations respond to future changes in complex natural systems.

Fitness advantage of inter-species TYLCV recombinants induced by beneficial intra-genomic interactions rather than by specific mutations

Tomato yellow leaf curl virus (TYLCV) and its related viruses are prone to recombination. It was reported that random homologous recombination between 20% diverging TYLCV related species is rarely deleterious and may be associated with a fitness advantage. Indeed, TYLCV-IS76, a recombinant between the 20% divergent TYLCV and tomato yellow leaf curl Sardinia virus (TYLCSV), exhibited a higher fitness than that of parental viruses. As this typical fitness advantage was observed with TYLCV-IS76 representatives of different pedigrees, it was thought that it is induced by beneficial intra-genomic interactions rather than by specific mutations. This hypothesis was further supported with TYLCV-IS141, a TYLCV recombinant with a short TYLCSV inherited fragment of around 141 nts, slightly longer than that of TYLCV-IS76. Indeed, the typical fitness advantage was detected irrespective of the position of the recombination breakpoint (loci 76 or 141) and the sequences of the TYLCV and TYLCSV inherited fragments.

Deciphering how plant pathogenic bacteria disperse and meet: Molecular epidemiology of Xanthomonas citri pv. citri at microgeographic scales in a tropical area of Asiatic citrus canker endemicity

Although some plant pathogenic bacteria represent a significant threat to agriculture, the determinants of their ecological success and evolutionary potential are still poorly understood. Refining our understanding of bacterial strain circulation at small spatial scales and the biological significance and evolutionary consequences of co-infections are key questions. The study of bacterial population biology can be challenging, because it requires high-resolution markers that can be genotyped with a high throughput. Here, we overcame this difficulty for Xanthomonas citri pv. citri, a genetically monomorphic bacterium causing Asiatic citrus canker (ACC). Using a genotyping method that did not require cultivating the bacterium or purifying DNA, we deciphered the pathogen's spatial genetic structure at several microgeographic scales, down to single lesion, in a situation of ACC endemicity. In a grove where copper was recurrently applied for ACC management, copper-susceptible and copper-resistant X. citri pv. citri coexisted and the bacterial population structured as three genetic clusters, suggesting a polyclonal contamination. The range of spatial dependency, estimated for the two largest clusters, was four times greater for the cluster predominantly composed of copper-resistant bacteria. Consistently, the evenness value calculated for this cluster was indicative of increased transmission. Linkage disequilibrium was high even at a tree scale, probably due to a combination of clonality and admixture. Approximately 1% of samples exhibited within-lesion multilocus polymorphism, explained at least in part by polyclonal infections. Canker lesions, which are of major biological significance as an inoculum source, may also represent a preferred niche for horizontal gene transfer. This study points out the potential of genotyping data for estimating the range of spatial dependency of plant bacterial pathogens, an important parameter for guiding disease management strategies.

Life on the Edge: Geminiviruses at the Interface Between Crops and Wild Plant Hosts

Viruses constitute the largest group of emerging pathogens, and geminiviruses (plant viruses with circular, single-stranded DNA genomes) are the major group of emerging plant viruses. With their high potential for genetic variation due to mutation and recombination, their efficient spread by vectors, and their wide host range as a group, including both wild and cultivated hosts, geminiviruses are attractive models for the study of the evolutionary and ecological factors driving virus emergence. Studies on the epidemiological features of geminivirus diseases have traditionally focused primarily on crop plants. Nevertheless, knowledge of geminivirus infection in wild plants, and especially at the interface between wild and cultivated plants, is necessary to provide a complete view of their ecology, evolution, and emergence. In this review, we address the most relevant aspects of geminivirus variability and evolution in wild and crop plants and geminiviruses' potential to emerge in crops.

Tomato Yellow Leaf Curl Sardinia Virus, a Begomovirus Species Evolving by Mutation and Recombination: A Challenge for Virus Control

The tomato leaf curl disease (TYLCD) is associated with infections of several species of begomoviruses (genus Begomovirus, family Geminiviridae) and causes severe damage to tomatoes throughout tropical and sub-tropical regions of the world. Among others, the Tomato yellow leaf curl Sardinia virus (TYLCSV) species causes damage in the Mediterranean Basin since early outbreaks occurred. Nevertheless, scarce information is available about the diversity of TYLCSV. Here, we study this aspect based on the sequence information accessible in databases. Isolates of two taxonomically differentiated TYLCSV strains can be found in natural epidemics. Their evolution is mostly associated with mutation combined with selection and random genetic drift and also with inter-species recombination which is frequent in begomoviruses. Moreover, a novel putative inter-strain recombinant is reported. Although no significantly new biological behaviour was observed for this latter recombinant, its occurrence supports that as shown for other related begomoviruses, recombination continues to play a central role in the evolution of TYLCD-associated viruses and the dynamism of their populations. The confrontation of resistant tomatoes with isolates of different TYLCD-associated viruses including the novel recombinant demonstrates the existence of a variable virus x plant genotype interaction. This has already been observed for other TYLCD-associated viruses and is a challenge for the control of their impact on tomato production.

Notes on recombination and reassortment in multipartite/segmented viruses

Besides evolving through nucleotide substitution, viruses frequently also evolve by genetic recombination which can occur when related viral variants co-infect the same cells. Viruses with segmented or multipartite genomes can additionally evolve via the reassortment of genomic components. Various computational techniques are now available for identifying and characterizing recombination and reassortment. While these techniques have revealed both that all well studied segmented and multipartite virus species show some capacity for reassortment, and that recombination is common in many multipartite species, they have indicated that recombination is either rare or does not occur in species with segmented genomes. Reassortment and recombination can make it very difficult to study segmented/multipartite viruses using metagenomics-based approaches. Notable challenges include, both the accurate identification and assignment of genomic components to individual genomes, and the differentiation between natural 'real' recombination events and artifactual 'fake' recombination events arising from the inaccurate de novo assembly of genome component sequences determined using short read sequencing.

Pervasive Chimerism in the Replication-Associated Proteins of Uncultured Single-Stranded DNA Viruses

Numerous metagenomic studies have uncovered a remarkable diversity of circular replication-associated protein (Rep)-encoding single-stranded (CRESS) DNA viruses, the majority of which are uncultured and unclassified. Unlike capsid proteins, the Reps show significant similarity across different groups of CRESS DNA viruses and have conserved domain organization with the N-terminal nuclease and the C-terminal helicase domain. Consequently, Rep is widely used as a marker for identification, classification and assessment of the diversity of CRESS DNA viruses. However, it has been shown that in certain viruses the Rep nuclease and helicase domains display incongruent evolutionary histories. Here, we systematically evaluated the co-evolutionary patterns of the two Rep domains across classified and unclassified CRESS DNA viruses. Our analysis indicates that the Reps encoded by members of the families Bacilladnaviridae, Circoviridae, Geminiviridae, Genomoviridae, Nanoviridae and Smacoviridae display largely congruent evolutionary patterns in the two domains. By contrast, among the unclassified CRESS DNA viruses, 71% appear to have chimeric Reps. Such massive chimerism suggests that unclassified CRESS DNA viruses represent a dynamic population in which exchange of gene fragments encoding the nuclease and helicase domains is extremely common. Furthermore, purging of the chimeric sequences uncovered six monophyletic Rep groups that may represent new families of CRESS DNA viruses.

Molecular evolution and invasion pattern of Cryphonectria hypovirus 1 in Europe: Mutation rate, and selection pressure differ between genome domains

Understanding virus evolution is a fundamental goal of virology, evolutionary biology, and disease epidemiology. We provide a detailed analysis of evolution and origin of Cryphonectria hypovirus 1 (CHV1) populations in Europe, based on the complete genome sequence of all European subtypes. Phylogenetic analyses divided European strains into two closely related clades. Strains of the subtype I belong to the first, while strains of the subtypes F1, D and E belong to the second clade suggesting that the subtypes F1, D and E are more closely related than previously thought. Strains of the subtype F2 appeared to be recombinant; subtypes F1/D/E contributed a larger fraction of sequence while subtype I contributed a smaller fraction. The p29 was the most variable domain, while the replication-associated large ORF B protein was the most conserved domain within the CHV1. Low sequence similarity, predominant negative selection and frequent recombination characterise the evolution of CHV1.

Geminiviruses and Plant Hosts: A Closer Examination of the Molecular Arms Race

Geminiviruses are plant-infecting viruses characterized by a single-stranded DNA (ssDNA) genome. Geminivirus-derived proteins are multifunctional and effective regulators in modulating the host cellular processes resulting in successful infection. Virus-host interactions result in changes in host gene expression patterns, reprogram plant signaling controls, disrupt central cellular metabolic pathways, impair plant's defense system, and effectively evade RNA silencing response leading to host susceptibility. This review summarizes what is known about the cellular processes in the continuing tug of war between geminiviruses and their plant hosts at the molecular level. In addition, implications for engineered resistance to geminivirus infection in the context of a greater understanding of the molecular processes are also discussed. Finally, the prospect of employing geminivirus-based vectors in plant genome engineering and the emergence of powerful genome editing tools to confer geminivirus resistance are highlighted to complete the perspective on geminivirus-plant molecular interactions.

The diversification of begomovirus populations is predominantly driven by mutational dynamics

Begomoviruses (single-stranded DNA plant viruses) are responsible for serious agricultural threats. Begomovirus populations exhibit a high degree of within-host genetic variation and evolve as quickly as RNA viruses. Although the recombination-prone nature of begomoviruses has been extensively demonstrated, the relative contribution of recombination and mutation to the genetic variation of begomovirus populations has not been assessed. We estimated the genetic variability of begomovirus datasets from around the world. An uneven distribution of genetic variation across the length of the cp and rep genes due to recombination was evident from our analyses. To estimate the relative contributions of recombination and mutation to the genetic variability of begomoviruses, we mapped all substitutions over maximum likelihood trees and counted the number of substitutions on branches which were associated with recombination (η_r) and mutation (η_μ). In addition, we also estimated the per generation relative rates of both evolutionary mechanisms (r/μ) to express how frequently begomovirus genomes are affected by recombination relative to mutation. We observed that the composition of genetic variation in all begomovirus datasets was dominated by mutation. Additionally, the low correlation between the estimates indicated that the relative contributions of recombination and mutation are not necessarily a function of their relative rates. Our results show that, although a considerable fraction of the genetic variation levels could be assigned to recombination, it was always lower than that due to mutation, indicating that the diversification of begomovirus populations is predominantly driven by mutational dynamics.

iTRAQ analysis of the tobacco leaf proteome reveals that RNA-directed DNA methylation (RdDM) has important roles in defense against geminivirus-betasatellite infection

Geminiviruses have caused serious losses in crop production. To investigate the mechanisms underlying host defenses against geminiviruses, an isobaric tags for relative and absolute quantification (iTRAQ)-based quantitative proteomic approach was used to explore the expression profiles of proteins in Nicotiana benthamiana (N. benthamiana) leaves in response to tomato yellow leaf curl China virus (TYLCCNV) with its betasatellite (TYLCCNB) at an early phase. In total, 4155 proteins were identified and 272 proteins were changed differentially in response to TYLCCNV/TYLCCNB infection. Bioinformatics analysis indicated that S-adenosyl-l-methionine cycle II was the most significantly up-regulated biochemical process during TYLCCNV/TYLCCNB infection. The mRNA levels of three proteins in S-adenosyl-l-methionine cycle II were further analyzed by qPCR, each was found significantly up-regulated in TYLCCNV/TYLCCNB-infected N. benthamiana. This result suggested a strong promotion of the biosynthesis of available methyl groups during geminivirus infections. We further tested the potential role of RdDM in N. benthamiana by virus-induced gene silencing (VIGS) and found that a disruption in RdDM resulted in more severe infectious symptoms and higher accumulation of viral DNA after TYLCCNV/TYLCCNB infection. Although the precise functions of these proteins still need to be determined, our proteomic results enhance the understanding of plant antiviral mechanisms.

BIOLOGICAL SIGNIFICANCE

One of the major limitations to crop growth in the worldwide is the prevalence of geminiviruses. They are able to infect food and cash crops and cause serious crop failures and economic losses worldwide, especially in Africa and Asia. Tomato yellow leaf curl China virus (TYLCCNV), which causes severe viral diseases in China, is a monopartite geminivirus associated with the betasatellite (TYLCCNB). However, the mechanisms underlying the TYLCCNV/TYLCCNB defense in plants are still not fully understood at the molecular level. In this study, the combined proteomic, bioinformatic and VIGS analyses revealed that TYLCCNV/TYLCCNB invasion caused complex proteomic alterations in the leaves of N. benthamiana involving the processes of stress and defense, energy production, photosynthesis, protein homeostasis, metabolism, cell structure, signal transduction, transcription, transportation, and cell growth/division. Promotion of available methyl groups via the S-adenosyl-l-methionine cycle II pathway in N. benthamiana appeared crucial for antiviral responses. These findings enhance our understanding in the proteomic aspects of host antiviral defenses against geminiviruses, and also demonstrate that the combination of proteomics with bioinformatics and VIGS analysis is an effective approach to investigate systemic plant responses to geminiviruses and to shed light on plant-virus interactions.

Evolutionary ecology of virus emergence

The cross-species transmission of viruses into new host populations, termed virus emergence, is a significant issue in public health, agriculture, wildlife management, and related fields. Virus emergence requires overlap between host populations, alterations in virus genetics to permit infection of new hosts, and adaptation to novel hosts such that between-host transmission is sustainable, all of which are the purview of the fields of ecology and evolution. A firm understanding of the ecology of viruses and how they evolve is required for understanding how and why viruses emerge. In this paper, I address the evolutionary mechanisms of virus emergence and how they relate to virus ecology. I argue that, while virus acquisition of the ability to infect new hosts is not difficult, limited evolutionary trajectories to sustained virus between-host transmission and the combined effects of mutational meltdown, bottlenecking, demographic stochasticity, density dependence, and genetic erosion in ecological sinks limit most emergence events to dead-end spillover infections. Despite the relative rarity of pandemic emerging viruses, the potential of viruses to search evolutionary space and find means to spread epidemically and the consequences of pandemic viruses that do emerge necessitate sustained attention to virus research, surveillance, prophylaxis, and treatment.

Grapevine red blotch-associated virus is Present in Free-Living Vitis spp. Proximal to Cultivated Grapevines

Red blotch is an emerging disease of grapevine associated with grapevine red blotch-associated virus (GRBaV). The virus spreads with infected planting stocks but no vector of epidemiological significance has been conclusively identified. A vineyard block of red-blotch-affected Vitis vinifera 'Cabernet franc' clone 214 was observed in California, with a clustering of infected, symptomatic vines focused along one edge of the field proximal to a riparian habitat with free-living Vitis spp. No genetic heterogeneity was observed in a 587-nucleotide region of the GRBaV genome in a population of 44 Cabernet franc clone 214 isolates. By contrast, genetic differences were observed in isolates from other cultivars and clones growing in adjacent blocks. GRBaV was confirmed infecting four free-living vines, two of which were shown to be V. californica × V. vinifera hybrids. The genomes of three free-living GRBaV vine isolates and seven from V. vinifera cultivars were compared; free-living vine isolates were shown to be more similar to each other and a 'Merlot' isolate than to the other cultivated vine isolates. The finding that GRBaV is present in free-living Vitis spp. indicates the virus can be spread by natural (nonhuman-mediated) means, and we hypothesize that in-field spread of GRBaV is occurring.

Mapping genetic determinants of viral traits with FST and quantitative trait locus (QTL) approaches

The genetic determinism of viral traits can generally be dissected using either forward or reverse genetics because the clonal reproduction of viruses does not require the use of approaches based on laboratory crosses. Nevertheless, we hypothesized that recombinant viruses could be analyzed as sexually reproducing organisms, using either a quantitative trait loci (QTL) approach or a locus-by-locus fixation index (F_ST). Locus-by-locus F_ST analysis, and four different regressions and interval mapping algorithms of QTL analysis were applied to a phenotypic and genotypic dataset previously obtained from 47 artificial recombinant genomes generated between two begomovirus species. Both approaches assigned the determinant of within-host accumulation—previously identified using standard virology approaches—to a region including the 5׳ end of the replication-associated protein (Rep) gene and the upstream intergenic region. This study provides a proof of principle that QTL and population genetics tools can be extended to characterize the genetic determinants of viral traits.

•

F_ST and QTL approaches were used to map the genetic determinants of viral traits.

•

Genetic determinants were detected using 47 begomovirus recombinant genomes.

•

Genetic determinants of begomovirus infectivity and accumulation were identified.

•

Proof of principles that F_ST and QTL can be used in virology.

Wild waterfowl migration and domestic duck density shape the epidemiology of highly pathogenic H5N8 influenza in the Republic of Korea

•

Phylogeographic analyses of H5N8, including 49 new sequences from South Korea.

•

H5N8 movement was mostly among areas dense in wild and domestic ducks.

•

New viral introductions to South Korea occurred at time of wild bird migration.

•

H5N8 epidemiology is shaped by wild waterfowl migration and domestic duck density.

•

H5N8 may have entered Europe at least twice, and Asia at least three times.

Highly pathogenic avian influenza (HPAI) viruses threaten human and animal health yet their emergence is poorly understood, partly because sampling of the HPAI Asian-origin H5N1 lineage immediately after its identification in 1996 was comparatively sparse. The discovery of a novel H5N8 virus in 2013 provides a new opportunity to investigate HPAI emergence in greater detail. Here we investigate the origin and transmission of H5N8 in the Republic of Korea, the second country to report the new strain. We reconstruct viral spread using phylogeographic methods and interpret the results in the context of ecological data on poultry density, overwintering wild bird numbers, and bird migration patterns. Our results indicate that wild waterfowl migration and domestic duck density were important to H5N8 epidemiology. Specifically, we infer that H5N8 entered the Republic of Korea via Jeonbuk province, then spread rapidly among western provinces where densities of overwintering waterfowl and domestic ducks are higher, yet rarely persisted in eastern regions. The common ancestor of H5N8 in the Republic of Korea was estimated to have arrived during the peak of inward migration of overwintering birds. Recent virus isolations likely represent re-introductions via bird migration from an as-yet unsampled reservoir. Based on the limited data from outside the Republic of Korea, our data suggest that H5N8 may have entered Europe at least twice, and Asia at least three times from this reservoir, most likely carried by wild migrating birds.

Global-scale computational analysis of genomic sequences reveals the recombination pattern and coevolution dynamics of cereal-infecting geminiviruses

Genetic diversity and recombination patterns were evaluated for 229 isolates of Wheat dwarf virus (WDV), which are important cereal-infecting geminiviruses. Recombination hot spots were concentrated at the boundary of the genes encoding for the replication protein (Rep), the coat protein (cp) and the movement protein (mp), as well as inside Rep and cp and in the short intergenic regions (SIR). Phylogenomic analyses confirmed that the global population of WDV clustered into two groups according to their specific host: wheat and barley, and the crucial regions for the division of two groups were mp and the large intergenic regions (LIR). The computationally inferred pattern of coevolution between amino acid residues and the predicted 3D structure for the viral proteins provided further differences among the strains or species at the genome and protein level. Pervasive interaction between Rep and Rep A proteins in WDV-wheat-specific group reflected their important and complex function in the replication and transcription of WDV. Furthermore, significant predicted interactions between CP and Rep and CP and Rep A proteins in the WDV-wheat-specific group are thought to be crucial for successful encapsidation and movement of the virus during infection.

Grapevine red blotch-associated virus Is Widespread in the United States

Grapevine red blotch disease has been recognized since 2008 as affecting North American grape production. The presence of the newly described Grapevine red blotch-associated virus (GRBaV) is highly correlated with the disease. To more effectively detect and monitor the presence of the virus, a sample processing strategy and multiplex polymerase chain reaction assay were developed. A total of 42 of 113 vine samples collected in or received from seven of the United States were shown to harbor the virus, demonstrating the virus is widely distributed across North America. Phylogenetic analyses of a viral replication-associated protein (Rep) gene fragment from the 42 isolates of GRBaV demonstrated distinct clades of the virus (1 and 2), with clade 1 showing the greatest variability. The full-length genome of six virus isolates was sequenced, and phylogenetic analyses of 14 whole genomes recapitulated results seen for the Rep gene. A comparison of GRBaV genomes revealed evidence of recombination underlying some of the variation seen among GRBaV genomes within clade 1. Phylogenetic analyses of coat and replicase-associated protein sequences among single-stranded DNA viruses showed GRBaV to group within the family Geminiviridae. This grouping is distinct from members of the families Nanoviridae and Circoviridae, with limited significant affinities to both recognized genera and novel plant-infecting, gemini-like viruses.

Contrasting genetic structure between two begomoviruses infecting the same leguminous hosts

Begomoviruses are whitefly-transmitted, ssDNA plant viruses and are among the most damaging pathogens causing epidemics in economically important crops worldwide. Wild/non-cultivated plants play a crucial epidemiological role, acting as begomovirus reservoirs and as ‘mixing vessels' where recombination can occur. Previous work suggests a higher degree of genetic variability in begomovirus populations from non-cultivated hosts compared with cultivated hosts. To assess this supposed host effect on the genetic variability of begomovirus populations, cultivated (common bean, Phaseolus vulgaris, and lima bean, Phaseolus lunatus) and non-cultivated (Macroptilium lathyroides) legume hosts were sampled from two regions of Brazil. A total of 212 full-length DNA-A genome segments were sequenced from samples collected between 2005 and 2012, and populations of the begomoviruses Bean golden mosaic virus (BGMV) and Macroptilium yellow spot virus (MaYSV) were obtained. We found, for each begomovirus species, similar genetic variation between populations infecting cultivated and non-cultivated hosts, indicating that the presumed genetic variability of the host did not a priori affect viral variability. We observed a higher degree of genetic variation in isolates from MaYSV populations than BGMV populations, which was explained by numerous recombination events in MaYSV. MaYSV and BGMV showed distinct distributions of genetic variation, with the BGMV population (but not MaYSV) being structured by both host and geography.

Tomato mottle wrinkle virus, a recombinant begomovirus infecting tomato in Argentina

Begomoviruses seriously threaten tomato production in South America. Here, we present the molecular characterization of a novel tomato-infecting begomovirus isolated in Argentina and demonstrate its infectivity. After cloning and sequencing the complete genome of this new virus, pairwise genetic distance and phylogenetic analyses revealed that it is a novel virus that is closely related to other begomoviruses found in Argentina, Brazil and Bolivia. We have proposed naming the virus tomato mottle wrinkle virus (ToMoWrV), based on symptoms produced upon its biolistic inoculation into tomato plants. Recombination analysis revealed that ToMoWrV is a recombinant, with parental sequences likely belonging to the South American begomoviruses soybean blistering mosaic virus (SoBlMV) and tomato yellow vein streak virus (ToYVSV).

SnRK1 phosphorylation of AL2 delays Cabbage leaf curl virus infection in Arabidopsis

Geminivirus AL2/C2 proteins play key roles in establishing infection and causing disease in their plant hosts. They are involved in viral gene expression, counter host defenses by suppressing transcriptional gene silencing, and interfere with the host signaling involved in pathogen resistance. We report here that begomovirus and curtovirus AL2/C2 proteins interact strongly with host geminivirus Rep-interacting kinases (GRIKs), which are upstream activating kinases of the protein kinase SnRK1, a global regulator of energy and nutrient levels in plants. We used an in vitro kinase system to show that GRIK-activated SnRK1 phosphorylates recombinant AL2/C2 proteins from several begomoviruses and to map the SnRK1 phosphorylation site to serine-109 in the AL2 proteins of two New World begomoviruses: Cabbage Leaf Curl Virus (CaLCuV) and Tomato mottle virus. A CaLCuV AL2 S109D phosphomimic mutation did not alter viral DNA levels in protoplast replication assays. In contrast, the phosphomimic mutant was delayed for symptom development and viral DNA accumulation during infection of Arabidopsis thaliana, demonstrating that SnRK1 contributes to host defenses against CaLCuV. Our observation that serine-109 is not conserved in all AL2/C2 proteins that are SnRK1 substrates in vitro suggested that phosphorylation of viral proteins by plant kinases contributes to the evolution of geminivirus-host interactions.

IMPORTANCE

Geminiviruses are single-stranded DNA viruses that cause serious diseases in many crops. Dicot-infecting geminiviruses carry genes that encode multifunctional AL2/C2 proteins that are essential for infection. However, it is not clear how AL2/C2 proteins are regulated. Here, we show that the host protein kinase SnRK1, a central regulator of energy balance and nutrient metabolism in plants, phosphorylates serine-109 in AL2 proteins of three subgroups of New World begomoviruses, resulting in a delay in viral DNA accumulation and symptom appearance. Our results support SnRK1's antiviral role and reveal a novel mechanism underlying this function. Phylogenetic analysis suggested that AL2 S109 evolved as begomoviruses migrated from the Old World to the New World and may have provided a selective advantage as begomoviruses adapted to a different environment and different plant hosts. This study provides new insights into the interaction of viral pathogens with their plant hosts at the level of viral protein modification by the host.

SnRK1 phosphorylation of AL2 delays Cabbage leaf curl virus infection in Arabidopsis

Geminivirus AL2/C2 proteins play key roles in establishing infection and causing disease in their plant hosts. They are involved in viral gene expression, counter host defenses by suppressing transcriptional gene silencing, and interfere with the host signaling involved in pathogen resistance. We report here that begomovirus and curtovirus AL2/C2 proteins interact strongly with host geminivirus Rep-interacting kinases (GRIKs), which are upstream activating kinases of the protein kinase SnRK1, a global regulator of energy and nutrient levels in plants. We used an in vitro kinase system to show that GRIK-activated SnRK1 phosphorylates recombinant AL2/C2 proteins from several begomoviruses and to map the SnRK1 phosphorylation site to serine-109 in the AL2 proteins of two New World begomoviruses: Cabbage Leaf Curl Virus (CaLCuV) and Tomato mottle virus. A CaLCuV AL2 S109D phosphomimic mutation did not alter viral DNA levels in protoplast replication assays. In contrast, the phosphomimic mutant was delayed for symptom development and viral DNA accumulation during infection of Arabidopsis thaliana, demonstrating that SnRK1 contributes to host defenses against CaLCuV. Our observation that serine-109 is not conserved in all AL2/C2 proteins that are SnRK1 substrates in vitro suggested that phosphorylation of viral proteins by plant kinases contributes to the evolution of geminivirus-host interactions.

IMPORTANCE

Geminiviruses are single-stranded DNA viruses that cause serious diseases in many crops. Dicot-infecting geminiviruses carry genes that encode multifunctional AL2/C2 proteins that are essential for infection. However, it is not clear how AL2/C2 proteins are regulated. Here, we show that the host protein kinase SnRK1, a central regulator of energy balance and nutrient metabolism in plants, phosphorylates serine-109 in AL2 proteins of three subgroups of New World begomoviruses, resulting in a delay in viral DNA accumulation and symptom appearance. Our results support SnRK1's antiviral role and reveal a novel mechanism underlying this function. Phylogenetic analysis suggested that AL2 S109 evolved as begomoviruses migrated from the Old World to the New World and may have provided a selective advantage as begomoviruses adapted to a different environment and different plant hosts. This study provides new insights into the interaction of viral pathogens with their plant hosts at the level of viral protein modification by the host.

Extensive recombination-induced disruption of genetic interactions is highly deleterious but can be partially reversed by small numbers of secondary recombination events

Although homologous recombination can potentially provide viruses with vastly more evolutionary options than are available through mutation alone, there are considerable limits on the adaptive potential of this important evolutionary process. Primary among these is the disruption of favorable coevolved genetic interactions that can occur following the transfer of foreign genetic material into a genome. Although the fitness costs of such disruptions can be severe, in some cases they can be rapidly recouped by either compensatory mutations or secondary recombination events. Here, we used a maize streak virus (MSV) experimental model to explore both the extremes of recombination-induced genetic disruption and the capacity of secondary recombination to adaptively reverse almost lethal recombination events. Starting with two naturally occurring parental viruses, we synthesized two of the most extreme conceivable MSV chimeras, each effectively carrying 182 recombination breakpoints and containing thorough reciprocal mixtures of parental polymorphisms. Although both chimeras were severely defective and apparently noninfectious, neither had individual movement-, encapsidation-, or replication-associated genome regions that were on their own "lethally recombinant." Surprisingly, mixed inoculations of the chimeras yielded symptomatic infections with viruses with secondary recombination events. These recombinants had only 2 to 6 breakpoints, had predominantly inherited the least defective of the chimeric parental genome fragments, and were obviously far more fit than their synthetic parents. It is clearly evident, therefore, that even when recombinationally disrupted virus genomes have extremely low fitness and there are no easily accessible routes to full recovery, small numbers of secondary recombination events can still yield tremendous fitness gains. Importance: Recombination between viruses can generate strains with enhanced pathological properties but also runs the risk of producing hybrid genomes with decreased fitness due to the disruption of favorable genetic interactions. Using two synthetic maize streak virus genome chimeras containing alternating genome segments derived from two natural viral strains, we examined both the fitness costs of extreme degrees of recombination (both chimeras had 182 recombination breakpoints) and the capacity of secondary recombination events to recoup these costs. After the severely defective chimeras were introduced together into a suitable host, viruses with between 1 and 3 secondary recombination events arose, which had greatly increased replication and infective capacities. This indicates that even in extreme cases where recombination-induced genetic disruptions are almost lethal, and 91 consecutive secondary recombination events would be required to reconstitute either one of the parental viruses, moderate degrees of fitness recovery can be achieved through relatively small numbers of secondary recombination events.

Patterns of recombination in HIV-1M are influenced by selection disfavouring the survival of recombinants with disrupted genomic RNA and protein structures

Genetic recombination is a major contributor to the ongoing diversification of HIV. It is clearly apparent that across the HIV-genome there are defined recombination hot and cold spots which tend to co-localise both with genomic secondary structures and with either inter-gene boundaries or intra-gene domain boundaries. There is also good evidence that most recombination breakpoints that are detectable within the genes of natural HIV recombinants are likely to be minimally disruptive of intra-protein amino acid contacts and that these breakpoints should therefore have little impact on protein folding. Here we further investigate the impact on patterns of genetic recombination in HIV of selection favouring the maintenance of functional RNA and protein structures. We confirm that chimaeric Gag p24, reverse transcriptase, integrase, gp120 and Nef proteins that are expressed by natural HIV-1 recombinants have significantly lower degrees of predicted folding disruption than randomly generated recombinants. Similarly, we use a novel single-stranded RNA folding disruption test to show that there is significant, albeit weak, evidence that natural HIV recombinants tend to have genomic secondary structures that more closely resemble parental structures than do randomly generated recombinants. These results are consistent with the hypothesis that natural selection has acted both in the short term to purge recombinants with disrupted RNA and protein folds, and in the longer term to modify the genome architecture of HIV to ensure that recombination prone sites correspond with those where recombination will be minimally deleterious.

Identifying geographic hot spots of reassortment in a multipartite plant virus

Reassortment between different species or strains plays a key role in the evolution of multipartite plant viruses and can have important epidemiological implications. Identifying geographic locations where reassortant lineages are most likely to emerge could be a valuable strategy for informing disease management and surveillance efforts. We developed a predictive framework to identify potential geographic hot spots of reassortment based upon spatially explicit analyses of genome constellation diversity. To demonstrate the utility of this approach, we examined spatial variation in the potential for reassortment among Cardamom bushy dwarf virus (CBDV; Nanoviridae, Babuvirus) isolates in Northeast India. Using sequence data corresponding to six discrete genome components for 163 CBDV isolates, a quantitative measure of genome constellation diversity was obtained for locations across the sampling region. Two key areas were identified where viruses with highly distinct genome constellations cocirculate, and these locations were designated as possible geographic hot spots of reassortment, where novel reassortant lineages could emerge. Our study demonstrates that the potential for reassortment can be spatially dependent in multipartite plant viruses and highlights the use of evolutionary analyses to identify locations which could be actively managed to facilitate the prevention of outbreaks involving novel reassortant strains.

A novel begomovirus isolated from sida contains putative cis- and trans-acting replication specificity determinants that have evolved independently in several geographical lineages

A novel begomovirus isolated from a Sida rhombifolia plant collected in Sinaloa, Mexico, was characterized. The genomic components of sida mosaic Sinaloa virus (SiMSinV) shared highest sequence identity with DNA-A and DNA-B components of chino del tomate virus (CdTV), suggesting a vertical evolutionary relationship between these viruses. However, recombination analysis indicated that a short segment of SiMSinV DNA-A encompassing the plus-strand replication origin and the 5´-proximal 43 codons of the Rep gene was derived from tomato mottle Taino virus (ToMoTV). Accordingly, the putative cis- and trans-acting replication specificity determinants of SiMSinV were identical to those of ToMoTV but differed from those of CdTV. Modeling of the SiMSinV and CdTV Rep proteins revealed significant differences in the region comprising the small β1/β5 sheet element, where five putative DNA-binding specificity determinants (SPDs) of Rep (i.e., amino acid residues 5, 8, 10, 69 and 71) were previously identified. Computer-assisted searches of public databases led to identification of 33 begomoviruses from three continents encoding proteins with SPDs identical to those of the Rep encoded by SiMSinV. Sequence analysis of the replication origins demonstrated that all 33 begomoviruses harbor potential Rep-binding sites identical to those of SiMSinV. These data support the hypothesis that the Rep β1/β5 sheet region determines specificity of this protein for DNA replication origin sequences.