Evolutionary dynamics of bipartite begomoviruses revealed by complete genome analysis

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.

Molecular characterization of a tetra segmented ssDNA virus infecting Botrytis cinerea worldwide

BACKGROUND

Family Genomoviridae was recently established, and only a few mycoviruses have been described and characterized, and almost all of them (Sclerotinia sclerotiorum hypovirulence-associated DNA virus 1, Fusarium graminearum gemyptripvirus 1 and Botrytis cinerea gemydayirivirus 1) induced hypovirulence in their host. Botrytis cinerea ssDNA virus 1 (BcssDV1), a tetrasegmented single-stranded DNA virus infecting the fungus Botrytis cinerea, has been molecularly characterized in this work.

METHODS

BcssDV1 was detected in Spanish and Italian B. cinerea field isolates obtained from grapevine. BcssDV1 variants genomes were molecularly characterized via NGS and Sanger sequencing. Nucleotide and amino acid sequences were used for diversity and phylogenetic analysis. Prediction of protein tertiary structures and putative associated functions were performed by AlphaFold2 and DALI.

RESULTS

BcssDV1 is a tetrasegmented single-stranded DNA virus. The mycovirus was composed by four genomic segments of approximately 1.7 Kb each, which are DNA-A, DNA-B, and DNA-C and DNA-D, that coded, respectively, for the rolling-circle replication initiation protein (Rep), capsid protein (CP) and two hypothetical proteins. BcssDV1 was present in several Italian and Spanish regions with high incidence and low variability among the different viral variants. DNA-A and DNA-D were found to be the more conserved genomic segments among variants, while DNA-B and DNA-C segments were shown to be the most variable ones. Tertiary structures of the proteins encoded by each segment suggested specific functions associated with each of them.

CONCLUSIONS

This study presented the first complete sequencing and characterization of a tetrasegmented ssDNA mycovirus, its incidence in Spain and Italy, its presence in other countries and its high conservation among regions.

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.

Constant ratio between the genomic components of bipartite begomoviruses during infection and transmission

The genomic components of multipartite viruses are encapsidated in separate virus particles, and the frequencies of genomic components represent one of the key genetic features. Many begomoviruses of economic significance are bipartite, and the details of the association between their genomic components remain largely unexplored. We first analyzed the temporal dynamics of the quantities of DNA-A and DNA-B and the B/A ratio of the squash leaf curl China virus (SLCCNV) in plants and found that while the quantities of DNA-A and DNA-B varied significantly during infection, the B/A ratio remained constant. We then found that changes in the B/A ratio in agrobacteria inoculum may significantly alter the B/A ratio in plants at 6 days post inoculation, but the differences disappeared shortly thereafter. We next showed that while the quantities of DNA-A and DNA-B among plants infected by agrobacteria, sap transmission and whitefly-mediated transmission differed significantly, the B/A ratios were similar. Further analysis of gene expression revealed that the ratio of the expression of genes encoded by DNA-A and DNA-B varied significantly during infection. Finally, we monitored the temporal dynamics of the quantities of DNA-A and DNA-B and the B/A ratio of another bipartite begomovirus, and a constant B/A ratio was similarly observed. Our findings highlight the maintenance of a constant ratio between the two genomic components of bipartite begomoviruses during infection and transmission, and provide new insights into the biology of begomoviruses.

Genome segment ratios change during whitefly transmission of two bipartite cassava mosaic begomoviruses

Cassava mosaic disease is caused by a complex of whitefly-transmitted begomoviruses, which often occur in co-infections. These viruses have bipartite genomes consisting of DNA-A and DNA-B that are encapsidated into separate virions. Individual viruses exist in plants and whitefly vectors as populations comprising both genome segments, which can occur at different frequencies. Both segments are required for infection, and must be transmitted for virus spread to occur. Cassava plants infected with African cassava mosaic virus (ACMV) and/or East African cassava mosaic Cameroon virus (EACMCV), in which the ratios of DNA-A:DNA-B titers differed between plants, were used to examine how titers of the segments in a plant relate to their respective probabilities of acquisition by whiteflies and to the titers of each segment acquired and subsequently transmitted by whiteflies. The probabilities of acquiring each segment of ACMV did not reflect their relative titers in the source plant but they did for EACMCV. However, for both viruses, DNA-A:DNA-B ratios acquired by whiteflies differed from those in the source plant and the ratios transmitted by the whitefly did not differ from one - the ratio at which the highest probability of transmitting both segments is expected.

In Silico Identification of Cassava Genome-Encoded MicroRNAs with Predicted Potential for Targeting the ICMV-Kerala Begomoviral Pathogen of Cassava

Cassava mosaic disease (CMD) is caused by several divergent species belonging to the genus Begomovirus (Geminiviridae) transmitted by the whitefly Bemisia tabaci cryptic species group. In India and other parts of Asia, the Indian cassava mosaic virus-Kerala (ICMV-Ker) is an emergent begomovirus of cassava causing damage that results in reduced yield loss and tuber quality. Double-stranded RNA-mediated interference (RNAi) is an evolutionary conserved mechanism in eukaryotes and highly effective, innate defense system to inhibit plant viral replication and/or translation. The objective of this study was to identify and characterize cassava genome-encoded microRNAs (mes-miRNA) that are predicted to target ICMV-Ker ssDNA-encoded mRNAs, based on four in silico algorithms: miRanda, RNA22, Tapirhybrid, and psRNA. The goal is to deploy the predicted miRNAs to trigger RNAi and develop cassava plants with resistance to ICMV-Ker. Experimentally validated mature cassava miRNA sequences (n = 175) were downloaded from the miRBase biological database and aligned with the ICMV-Ker genome. The miRNAs were evaluated for base-pairing with the cassava miRNA seed regions and to complementary binding sites within target viral mRNAs. Among the 175 locus-derived mes-miRNAs evaluated, one cassava miRNA homolog, mes-miR1446a, was identified to have a predicted miRNA target binding site, at position 2053 of the ICMV-Ker genome. To predict whether the cassava miRNA might bind predicted ICMV-Ker mRNA target(s) that could disrupt viral infection of cassava plants, a cassava locus-derived miRNA-mRNA regulatory network was constructed using Circos software. The in silico-predicted cassava locus-derived mes-miRNA-mRNA network corroborated interactions between cassava mature miRNAs and the ICMV-Ker genome that warrant in vivo analysis, which could lead to the development of ICMV-Ker resistant cassava plants.

Genetic Diversity and Geographic Distribution of Cucurbit-Infecting Begomoviruses in the Philippines

Cucurbits are important economic crops worldwide. However, the cucurbit leaf curl disease (CuLCD), caused by whitefly-transmitted begomoviruses constrains their production. In Southeast Asia, three major begomoviruses, Tomato leaf curl New Delhi virus (ToLCNDV), Squash leaf curl China virus (SLCCNV) and Squash leaf curl Philippines virus (SLCuPV) are associated with CuLCD. SLCuPV and SLCCNV were identified in Luzon, the Philippines. Here, the genetic diversity and geographic distribution of CuLCD-associated begomoviruses in the Philippines were studied based on 103 begomovirus detected out of 249 cucurbit samples collected from 60 locations throughout the country in 2018 and 2019. The presence of SLCCNV and SLCuPV throughout the Philippines were confirmed by begomovirus PCR detection and viral DNA sequence analysis. SLCuPV was determined as a predominant CuLCD-associated begomovirus and grouped into two strains. Interestingly, SLCCNV was detected in pumpkin and bottle gourd without associated viral DNA-B and mixed-infected with SLCuPV. Furthermore, the pathogenicity of selected isolates of SLCCNV and SLCuPV was confirmed. The results provide virus genetic diversity associated with CuLCD for further disease management, especially in developing the disease-resistant cultivars in the Philippines as well as Southeast Asia.

Limited role of recombination in the global diversification of begomovirus DNA-B proteins

Approximately half of the characterized begomoviruses have bipartite genomes, but the second genomic segment, the DNA-B, is understudied relative to the DNA-A, which is homologous to the entire genome of monopartite begomoviruses. We examined the evolutionary history of the two proteins encoded by the DNA-B, the genes of which make up ∼60% of the DNA-B segment, from all bipartite begomovirus species. Our dataset of 131 movement protein (MP) and nuclear shuttle protein (NSP) sequences confirmed the deep split between Old World (OW) and New World (NW) species, and showed strong support for deep, congruent branches among the OW sequences of the MP and NSP. NW sequences were much less diverse and had poor phylogenetic resolution; over half of nodes in both the NSP and MP NW clades were supported by <50% bootstrap support. This poor resolution hampered our ability to detect incongruent phylogenies between the MP and NSP datasets, and we found no statistical evidence for recombination within our MP and NSP datasets. Finally, we quantified the sequence diversity between the NW and OW proteins, showing that the NW MP has particularly low diversity, suggesting it has been subject to different evolutionary pressures than the NW NSP.

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.

Roles of two distinct alphasatellites modulating geminivirus pathogenesis

BACKGROUND

Alphasatellites are small coding DNA satellites frequently associated with a begomovirus/betasatellite complex, where they are known to modulate virulence and symptom development. Two distinct alphasatellites, namely, Cotton leaf curl Multan alphasatellite (CLCuMuA), and Gossypium darwinii symptomless alphasatellite (GDarSLA) associated with Cotton leaf curl Multan virus-India (CLCuMuV-IN) and Ludwigia leaf distortion betasatellite (LuLDB) were found to be associated with yellow mosaic disease of hollyhock (Alcea rosea) plants. In this study, we show that alphasatellites CLCuMuA and GDarSLA attenuate and delay symptom development in Nicotiana benthamiana. The presence of either alphasatellites reduce the accumulation of the helper virus CLCuMuV-IN. However, the levels of the associated betasatellite, LuLDB, remains unchanged. These results suggest that the alphasatellites could contribute to the host defence and understanding their role in disease development is important for developing resistance strategies.

METHODS

Tandem repeat constructs of two distinct alphasatellites, namely, CLCuMuA and GDarSLA associated with CLCuMuV-IN and LuLDB were generated. N. benthamiana plants were co-agroinoculated with CLCuMuV and its associated alphasatellites and betasatellite molecules and samples were collected at 7, 14 and 21 days post inoculation (dpi). The viral DNA molecules were quantified in N. benthamiana plants by qPCR. The sequences were analysed using the MEGA-X tool, and a phylogenetic tree was generated. Genetic diversity among the CLCuMuA and GDarSLA was analysed using the DnaSP tool.

RESULTS

We observed a reduction in symptom severity and accumulation of helper virus in the presence of two alphasatellites isolated from naturally infected hollyhock plants. However, no reduction in the accumulation of betasatellite was observed. The phylogenetic and genetic variability study revealed the evolutionary dynamics of these distinct alphasatellites , which could explain the role of hollyhock-associated alphasatellites in plants.

CONCLUSIONS

This study provides evidence that alphasatellites have a role in symptom modulation and suppress helper virus replication without any discernible effect on the replication of the associated betasatellite.

Small but mighty: Functional landscape of the versatile geminivirus-encoded C4 protein

The fast-paced evolution of viruses enables them to quickly adapt to the organisms they infect by constantly exploring the potential functional landscape of the proteins encoded in their genomes. Geminiviruses, DNA viruses infecting plants and causing devastating crop diseases worldwide, produce a limited number of multifunctional proteins that mediate the manipulation of the cellular environment to the virus’ advantage. Among the proteins produced by the members of this family, C4, the smallest one described to date, is emerging as a powerful viral effector with unexpected versatility. C4 is the only geminiviral protein consistently subjected to positive selection and displays a number of dynamic subcellular localizations, interacting partners, and functions, which can vary between viral species. In this review, we aim to summarize our current knowledge on this remarkable viral protein, encompassing the different aspects of its multilayered diversity, and discuss what it can teach us about geminivirus evolution, invasion requirements, and virulence strategies.

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.

Capsicum-infecting begomoviruses as global pathogens: host-virus interplay, pathogenesis, and management

Whitefly-transmitted begomoviruses are among the major threats to the cultivation of Capsicum spp. (Family: Solanaceae) worldwide. Capsicum-infecting begomoviruses (CIBs) have a broad host range and are commonly found in mixed infections, which, in turn, fuels the emergence of better-adapted species through intraspecies and interspecies recombination. Virus-encoded proteins hijack host factors to breach the well-coordinated antiviral response of plants. Epigenetic modifications of histones associated with viral minichromosomes play a critical role in this molecular arms race. Moreover, the association of DNA satellites further enhances the virulence of CIBs as the subviral agents aid the helper viruses to circumvent plant antiviral defense and facilitate expansion of their host range and disease development. The objective of this review is to provide a comprehensive overview on various aspects of CIBs such as their emergence, epidemiology, mechanism of pathogenesis, and the management protocols being employed for combating them.

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.