Validation of high rates of nucleotide substitution in geminiviruses: phylogenetic evidence from East African cassava mosaic viruses

Cassava brown streak virus evolves with a nucleotide-substitution rate that is typical for the family Potyviridae

The ipomoviruses (family Potyviridae) that cause cassava brown streak disease (cassava brown streak virus [CBSV] and Uganda cassava brown streak virus [UCBSV]) are damaging plant pathogens that affect the sustainability of cassava production in East and Central Africa. However, little is known about the rate at which the viruses evolve and when they emerged in Africa - which inform how easily these viruses can host shift and resist RNAi approaches for control. We present here the rates of evolution determined from the coat protein gene (CP) of CBSV (Temporal signal in a UCBSV dataset was not sufficient for comparable analysis). Our BEAST analysis estimated the CBSV CP evolves at a mean rate of 1.43 × 10-3 nucleotide substitutions per site per year, with the most recent common ancestor of sampled CBSV isolates existing in 1944 (95% HPD, between years 1922 - 1963). We compared the published measured and estimated rates of evolution of CPs from ten families of plant viruses and showed that CBSV is an average-evolving potyvirid, but that members of Potyviridae evolve more quickly than members of Virgaviridae and the single representatives of Betaflexiviridae, Bunyaviridae, Caulimoviridae and Closteroviridae.

Taxonomic Classification of Geminiviruses Based on Pairwise Sequence Comparisons

Geminiviridae is the largest and one of the most diverse families of plant viruses, comprising 14 genera demarcated based on host range, type of insect vector, and phylogenetic relationships. The use of unbiased, whole-genome multiple displacement amplification techniques coupled with high-throughput sequencing has greatly expanded our knowledge of geminivirus diversity over the last two decades. As a result, a large number of new species have been described in recent years. Species demarcation criteria in the family are entirely based on sequence comparisons, but the specific cutoff values vary for each genus. The purpose of this chapter is to provide a step-by-step pipeline to classify new species in the family Geminiviridae.

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.

SEGS-1 a cassava genomic sequence increases the severity of African cassava mosaic virus infection in Arabidopsis thaliana

Cassava is a major crop in Sub-Saharan Africa, where it is grown primarily by smallholder farmers. Cassava production is constrained by Cassava mosaic disease (CMD), which is caused by a complex of cassava mosaic begomoviruses (CMBs). A previous study showed that SEGS-1 (sequences enhancing geminivirus symptoms), which occurs in the cassava genome and as episomes during viral infection, enhances CMD symptoms and breaks resistance in cassava. We report here that SEGS-1 also increases viral disease severity in Arabidopsis thaliana plants that are co-inoculated with African cassava mosaic virus (ACMV) and SEGS-1 sequences. Viral disease was also enhanced in Arabidopsis plants carrying a SEGS-1 transgene when inoculated with ACMV alone. Unlike cassava, no SEGS-1 episomal DNA was detected in the transgenic Arabidopsis plants during ACMV infection. Studies using Nicotiana tabacum suspension cells showed that co-transfection of SEGS-1 sequences with an ACMV replicon increases viral DNA accumulation in the absence of viral movement. Together, these results demonstrated that SEGS-1 can function in a heterologous host to increase disease severity. Moreover, SEGS-1 is active in a host genomic context, indicating that SEGS-1 episomes are not required for disease enhancement.

Amesuviridae: a new family of plant-infecting viruses in the phylum Cressdnaviricota, realm Monodnaviria

The phylum Cressdnaviricota comprises viruses with single-stranded, circular DNA genomes that encode an HUH-type endonuclease (known as Rep). The phylum includes two classes, eight orders, and 11 families. Here, we report the creation of a twelfth family in the order Mulpavirales, class Arfiviricetes of the phylum Cressdnaviricota. The family Amesuviridae comprises viruses that infect plants and is divided into two genera: Temfrudevirus, including the species Temfrudevirus temperatum (with temperate fruit decay-associated virus as a member), and Yermavirus, including the species Yermavirus ilicis (with yerba mate-associated circular DNA virus as a member). Both viruses encode Rep proteins with HUH endonuclease and SH3 superfamily helicase domains. Phylogenetic analysis indicates that the replicative module of amesuviruses constitutes a well-supported monophyletic clade related to Rep proteins from viruses in the order Mulpavirales. Furthermore, both viruses encode a single capsid protein (CP) related to geminivirus CPs. Phylogenetic incongruence between the replicative and structural modules of amesuviruses suggests a chimeric origin resulting from remote recombination events between ancestral mulpavirales and geminivirids. The creation of the family Amesuviridae has been ratified by the International Committee on Taxonomy of Viruses (ICTV).

A shift of vector specificity acquired by a begomovirus through natural homologous recombination

Recombination is common in plant viruses such as geminiviruses, but the ecological and pathogenic consequences have been explored only in a few cases. Here, we found that a new begomovirus, tomato yellow leaf curl Shuangbai virus (TYLCSbV), probably originated from the recombination of Ageratum yellow vein China virus (AYVCNV) and tobacco curl shoot virus (TbCSV). Agrobacterium-mediated inoculation showed that TYLCSbV and AYVCNV have similar levels of infectivity on tomato and tobacco plants. However, the two viruses exhibit contrasting specificities for vector transmission, that is, TYLCSbV was efficiently transmitted by the whitefly Bemisia tabaci Mediterranean (MED) rather than by the whitefly B. tabaci Middle East-Asia Minor 1 (MEAM1), whereas AYVCNV was more efficiently transmitted by MEAM1. We also showed that the transmission efficiencies of TYLCSbV and AYVCNV are positively correlated with the accumulation of the viruses in whitefly whole bodies and organs/tissues. The key coat protein amino acids that determine their accumulation are between positions 147 and 256. Moreover, field surveys suggest that MED has displaced MEAM1 in some regions where TYLCSbV was collected. Viral competition assays indicated that TYLCSbV outcompeted AYVCNV when transmitted by MED, while the outcome was the opposite when transmitted by MEAM1. Our findings suggest that recombination has resulted in a shift of vector specificity that could provide TYLCSbV with a potential selective transmission advantage, and the population shift of whitefly cryptic species could have influenced virus evolution towards an extended trajectory of transmission.

Sentinel plot surveillance of cotton leaf curl disease in Pakistan- a case study at the cultivated cotton-wild host plant interface

A sentinel plot case study was carried out to identify and map the distribution of begomovirus-betasatellite complexes in sentinel plots and commercial cotton fields over a four-year period using molecular and high-throughput DNA 'discovery' sequencing approaches. Samples were collected from 15 study sites in the two major cotton-producing areas of Pakistan. Whitefly- and leafhopper-transmitted geminiviruses were detected in previously unreported host plant species and locations. The most prevalent begomovirus was cotton leaf curl Kokhran virus-Burewala (CLCuKoV-Bu). Unexpectedly, a recently recognized recombinant, cotton leaf curl Multan virus-Rajasthan (CLCuMuV-Ra) was prevalent in five of 15 sites. cotton leaf curl Alabad virus (CLCuAlV) and cotton leaf curl Kokhran virus-Kokhran, 'core' members of CLCuD-begomoviruses that co-occurred with CLCuMuV in the 'Multan' epidemic were detected in one of 15 sentinel plots. Also identified were chickpea chlorotic dwarf virus and 'non-core' CLCuD-begomoviruses, okra enation leaf curl virus, squash leaf curl virus, and tomato leaf curl New Delhi virus. Cotton leaf curl Multan betasatellite (CLCuMuB) was the most prevalent CLCuD-betasatellite, and less commonly, two 'non-core' betasatellites. Recombination analysis revealed previously uncharacterized recombinants among helper virus-betasatellite complexes consisting of CLCuKoV, CLCuMuV, CLCuAlV and CLCuMuB. Population analyses provided early evidence for CLCuMuV-Ra expansion and displacement of CLCuKoV-Bu in India and Pakistan from 2012-2017. Identification of 'core' and non-core CLCuD-species/strains in cotton and other potential reservoirs, and presence of the now predominant CLCuMuV-Ra strain are indicative of ongoing diversification. Investigating the phylodynamics of geminivirus emergence in cotton-vegetable cropping systems offers an opportunity to understand the driving forces underlying disease outbreaks and reconcile viral evolution with epidemiological relationships that also capture pathogen population shifts.

An Insight into Emerging Begomoviruses and their Satellite Complex causing Papaya Leaf Curl Disease

Papaya leaf curl disease (PaLCD) was primarily detected in India and causes major economic damage to agriculture crops grown globally, seriously threatening food security. Begomoviruses are communicated by the vector Bemisia tabaci, and their transmission efficiency and persistence in the vector are the highest, exhibiting the widest host range due to adaptation and evolution. Symptoms induced during PaLCD include leaf curl, leaf yellowing, interveinal chlorosis, and reduced fruit quality and yield. Consequently, plants have evolved several multi-layered defense mechanisms to resist Begomovirus infection and distribution. Subsequently, Begomovirus genomes organise circular ssDNA of size ~2.5-2.7 kb of overlapping viral transcripts and carry six-seven ORFs encoding multifunctional proteins, which are precisely evolved by the viruses to maintain the genome-constraint and develop complex but integrated interactions with a variety of host components to expand and facilitate successful infection cycles, i.e., suppression of host defense strategies. Geographical distribution is continuing to increase due to the advent and evolution of new Begomoviruses, and sweep to new regions is a future scenario. This review summarizes the current information on the biological functions of papaya-infecting Begomoviruses and their encoded proteins in transmission through vectors and modulating host-mediated responses, which may improve our understanding of how to challenge these significant plant viruses by revealing new information on the development of antiviral approaches against Begomoviruses associated with PaLCD.

Cassava begomovirus species diversity changes during plant vegetative cycles

Cassava is a root crop important for global food security and the third biggest source of calories on the African continent. Cassava production is threatened by Cassava mosaic disease (CMD), which is caused by a complex of single-stranded DNA viruses (family: Geminiviridae, genus: Begomovirus) that are transmitted by the sweet potato whitefly (Bemisia tabaci). Understanding the dynamics of different cassava mosaic begomovirus (CMB) species through time is important for contextualizing disease trends. Cassava plants with CMD symptoms were sampled in Lake Victoria and coastal regions of Kenya before transfer to a greenhouse setting and regular propagation. The field-collected and greenhouse samples were sequenced using Illumina short-read sequencing and analyzed on the Galaxy platform. In the field-collected samples, African cassava mosaic virus (ACMV), East African cassava mosaic virus (EACMV), East African cassava mosaic Kenya virus (EACMKV), and East African cassava mosaic virus-Uganda variant (EACMV-Ug) were detected in samples from the Lake Victoria region, while EACMV and East African mosaic Zanzibar virus (EACMZV) were found in the coastal region. Many of the field-collected samples had mixed infections of EACMV and another begomovirus. After 3 years of regrowth in the greenhouse, only EACMV-like viruses were detected in all samples. The results suggest that in these samples, EACMV becomes the dominant virus through vegetative propagation in a greenhouse. This differed from whitefly transmission results. Cassava plants were inoculated with ACMV and another EACMV-like virus, East African cassava mosaic Cameroon virus (EACMCV). Only ACMV was transmitted by whiteflies from these plants to recipient plants, as indicated by sequencing reads and copy number data. These results suggest that whitefly transmission and vegetative transmission lead to different outcomes for ACMV and EACMV-like viruses.

Molecular Characterization of Three Apple Geminivirus Isolates in Crabapples Detected in Inner Mongolia, China

Apple geminivirus 1 (AGV) in the genus Maldovirus of the family Geminiviridae was first identified infecting apple trees in the year 2015 in China. In this work, we characterized three isolates of the AGV in the Chinese pearleaf crabapple (Malus asiatica) in Inner Mongolia Autonomous Region. The viruses were detected by Illumina sequencing and its existence was confirmed by reverse transcription-polymerase chain reaction (RT-PCR) amplification of an AGV fragment. Between the three AGV isolates and the initially characterized AGV isolate PL2015, the nucleotide sequence identities of the complete genome ranged from 91.2 to 91.7%, of the coat protein gene (V1) ranged from 95.4% to 97.3%, and of the replicase gene (C1) ranged from 87.3% to 88.0%. Phylogenetic analysis indicated that the three isolates formed a monophyletic group together with the AGV, separated from the current genera in the family Geminiviridae. This is the first description of the AGV infecting crabapples.

Begomovirus populations in single plants are complex and may include both well-adapted and poorly-adapted viruses

Begomoviruses (single-stranded DNA plant viruses transmitted by whiteflies) are economically important pathogens causing epidemics worldwide. Tomato-infecting begomoviruses emerged in Brazil in the 1990's following the introduction of Bemisia tabaci Middle East-Asia Minor 1. It is believed that these viruses evolved from indigenous viruses infecting non-cultivated hosts. However, tomato-infecting viruses are rarely found in non-cultivated hosts, and vice-versa. It is possible that viral populations in a given host are composed primarily of viruses which are well adapted to this host, but also include a small proportion of poorly adapted viruses. Following transfer to a new host, the composition of the viral population would shift rapidly, with the viruses which are better adapted to the new host becoming predominant. To test this hypothesis, we collected tomato and Sida plants growing next to each other at two locations in 2014 and 2018. Total DNA was extracted from tomato and Sida samples from each location and year and used as a template for high-throughput sequencing. Reads were mapped following a highly stringent set of criteria. For the 2014 samples, >98% of the Sida reads mapped to Sida micrantha mosaic virus (SiMMV), but 0.1% of the reads mapped to tomato severe rugose virus (ToSRV). Conversely, >99% of the tomato reads mapped to ToSRV, with 0.18% mapping to SiMMV. For the 2018 samples, 41% of the Sida reads mapped to three Sida-adapted viruses and 0.1% of the reads mapped to ToSRV, while 99.9% of the tomato reads mapped to ToSRV. These results are consistent with the hypothesis that viral populations in a single plant are composed primarily of the virus that is better adapted to the host but also include a small proportion of viruses that are poorly adapted.

Phylostems: a new graphical tool to investigate temporal signal of heterochronous sequences datasets

Motivation

Molecular tip-dating of phylogenetic trees is a growing discipline that uses DNA sequences sampled at different points in time to co-estimate the timing of evolutionary events with rates of molecular evolution. Importantly, such inferences should only be performed on datasets displaying sufficient temporal signal, a feature important to test prior to any tip-dating inference. For this purpose, the most popular method considered to-date has been the 'root-to-tip regression' which consist in fitting a linear regression of the number of substitutions accumulated from the root to the tips of a phylogenetic tree as a function of sampling times. The main limitation of the regression method, in its current implementation, relies in the fact that the temporal signal can only be tested at the whole-tree scale (i.e. its root).

Results

To overcome this limitation we introduce Phylostems, a new graphical user-friendly tool developed to investigate temporal signal within every clade of a phylogenetic tree. We provide a 'how to' guide by running Phylostems on an empirical dataset and supply guidance for results interpretation.

Availability and implementation

Phylostems is freely available at https://pvbmt-apps.cirad.fr/apps/phylostems.

Analysis of the genome of grapevine red blotch virus and related grabloviruses indicates diversification prior to the arrival of Vitis vinifera in North America

In this study 163 complete whole-genome sequences of the emerging pathogen grapevine red blotch virus (GRBV; genus Grablovirus, family Geminiviridae) were used to reconstruct phylogenies using Bayesian analyses on time-tipped (heterochronous) data. Using different combinations of priors, Bayes factors identified heterochronous datasets (3×200 million chains) generated from strict clock and exponential tree priors as being the most robust. Substitution rates of 3.2×10^-5 subsitutions per site per year (95% HPD 4.3-2.1×10^-5) across the whole of the GRBV genome were estimated, suggesting ancestral GRBV diverged from ancestral wild Vitis latent virus 1 around 9 000 years ago, well before the first documented arrival of Vitis vinifera in North America. Whole-genome analysis of GRBV isolates in a single infected field-grown grapevine across 12 years identified 12 single nucleotide polymorphisms none of which were fixed substitutions: an observation not discordant with the in silico estimate. The substitution rate estimated here is lower than those estimated for other geminiviruses and is the first for a woody-host-infecting geminivirus.

Identification and characterization of a new geminivirus from soybean plants and determination of V2 as a pathogenicity factor and silencing suppressor

BACKGROUND

Soybean is one of the four major crops in China. The occurrence of viruses in soybean causes significant economic losses.

RESULTS

In this study, the soybean leaves from stay-green plants showing crinkle were collected for metatranscriptomic sequencing. A novel geminivirus, tentatively named soybean geminivirus A (SGVA), was identified in soybean stay-green plants. Sequence analysis of the full-length SGVA genome revealed a genome of 2762 nucleotides that contain six open reading frames. Phylogenetic analyses revealed that SGVA was located adjacent to the clade of begomoviruses in both the full genome-based and C1-based phylogenetic tree, while in the CP-based phylogenetic tree, SGVA was located adjacent to the clade of becurtoviruses. SGVA was proposed as a new recombinant geminivirus. Agroinfectious clone of SGVA was constructed. Typical systemic symptoms of curly leaves were observed at 11 dpi in Nicotiana benthamiana plants and severe dwarfism was observed after 3 weeks post inoculation. Expression of the SGVA encoded V2 and C1 proteins through a potato virus X (PVX) vector caused severe symptoms in N. benthamiana. The V2 protein inhibited local RNA silencing in co-infiltration assays in GFP transgenic 16C N. benthamiana plants. Further study revealed mild symptoms in N. benthamiana plants inoculated with SGVA-ZZ _V2-STOP and SGVA-ZZ _V2-3738AA mutants. Both the relative viral DNA and CP protein accumulation levels significantly decreased when compared with SGVA-inoculated plants.

CONCLUSIONS

This work identified a new geminivirus in soybean stay-green plants and determined V2 as a pathogenicity factor and silencing suppressor.

Phylodynamics and Coat Protein Analysis of Babaco Mosaic Virus in Ecuador

Babaco is a fast-growing herbaceous shrub with great commercial potential because of the organoleptic properties of its fruit. Babaco mosaic virus (BabMV) is a potexvirus in the family Alphaflexiviridae affecting babaco in all the provinces that produce this crop in Ecuador. BabMV was recently described but it has been affecting babaco for decades and, since many potexviruses are serologically indistinguishable, it may have been previously misidentified as papaya mosaic virus. Based on the coat protein (CP) gene, we aimed to study the distribution and epidemiological patterns of BabMV in babaco and chamburo over the years and to model its three-dimensional structure. Sequences of the CP were obtained from thirty-six isolates from plants collected in the main babaco-producing provinces of Ecuador between 2016 and 2021. The evolution rate of BabMV was estimated at 1.21 × 10−3 nucleotide substitutions site−1 year−1 and a time of origin of the most recent common ancestor around 1958.80. From molecular dynamics simulations, compared to other proteins of BabMV—RDRP, TGB1, and Alkb domain—the CP exhibited a higher flexibility with the C and N terminals as the most flexible regions. The reconstructed viral distribution provides dispersion patterns which have implications for control approaches of BabMV.

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.

Chilli leaf curl disease populations in India are highly recombinant, and rapidly segregated

Capsicum annuum, a valuable spice and vegetable crop belonging to the Solanaceae family, is extensively grown across the Indian subcontinent. Chilli production is restricted by a begomoviral infection named as chilli leaf curl disease (ChiLCD) mainly in tropical and subtropical regions which leads to considerable economic losses, thus affecting chilli cultivation. Here, we studied the genetic diversity with structural evaluation of chilli leaf curl disease and satellite molecules infecting Chilli in India. We retrieved 121 reference sequences of ChiLCD including DNA-A, DNA-B, beta-satellite and alpha-satellites from GenBank reported from India. The population diversity and genetic variation were estimated through various parameters which decipher the four major groups of phylogenetic divergence for DNA-A and five groups of beta-satellite showing percentage similarity with isolates within and across India. Further, transitional and transversional bias for ORFs were observed highest in C4 and REn genes, respectively, and for DNA-A and DNA-B, these values were 1.07 and 1.22, respectively. The recombination breakpoints for DNA-A were estimated 49 majorly in V1, C1,C2 and C4 genome region and highest 22 breakpoints were determined for Rep (AC1) of ORFs, similarly 9 events for beta-satellite were found less around βC1ORF. Moreover, the evolution and genetic variability were also contributed through parameters such as nucleotide substitution which were found within the range of RNA viruses for DNA-A, DNA-B, for all 6 ORFs (relaxed clock) and beta-satellite. Additionally, total numbers of mutations (η) for DNA-A, DNA-B, alpha-satellites and beta-satellites were 2505, 419, 807 and 1288 detected, respectively, while it was found 987 highest for Rep gene among all ORFs. Further, neutrality tests determine the dominant nature of population expansion and purifying selection for all the genes of begomovirus associated with ChiLCD and satellite molecules supporting conserved nature of gene. The combined Tajima's D and Fu and Li'S D* negative values in tests indicated that population are under purified selection and an excess of low-frequency polymorphism. Our analysis indicates the potential contribution of genetic mutations and recombination of ChiLCD which leads to rapid adaptation and evolution of begomovirus and its satellite molecules accelerating its host range and diversity within and across the Indian subcontinent.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s13205-022-03139-w.

Contribution of historical herbarium small RNAs to the reconstruction of a cassava mosaic geminivirus evolutionary history

Emerging viral diseases of plants are recognised as a growing threat to global food security. However, little is known about the evolutionary processes and ecological factors underlying the emergence and success of viruses that have caused past epidemics. With technological advances in the field of ancient genomics, it is now possible to sequence historical genomes to provide a better understanding of viral plant disease emergence and pathogen evolutionary history. In this context, herbarium specimens represent a valuable source of dated and preserved material. We report here the first historical genome of a crop pathogen DNA virus, a 90-year-old African cassava mosaic virus (ACMV), reconstructed from small RNA sequences bearing hallmarks of small interfering RNAs. Relative to tip-calibrated dating inferences using only modern data, those performed with the historical genome yielded both molecular evolution rate estimates that were significantly lower, and lineage divergence times that were significantly older. Crucially, divergence times estimated without the historical genome appeared in discordance with both historical disease reports and the existence of the historical genome itself. In conclusion, our study reports an updated time-frame for the history and evolution of ACMV and illustrates how the study of crop viral diseases could benefit from natural history collections.

A New Type of Satellite Associated with Cassava Mosaic Begomoviruses

Cassava mosaic disease (CMD), which is caused by single-stranded DNA begomoviruses, severely limits cassava production across Africa. A previous study showed that CMD symptom severity and viral DNA accumulation increase in cassava in the presence of a DNA sequence designated SEGS-2 (sequence enhancing geminivirus symptoms). We report here that when SEGS-2 is coinoculated with African cassava mosaic virus (ACMV) onto Arabidopsis thaliana, viral symptoms increase. Transgenic Arabidopsis with an integrated copy of SEGS-2 inoculated with ACMV also display increased symptom severity and viral DNA levels. Moreover, SEGS-2 enables Cabbage leaf curl virus (CaLCuV) to infect a geminivirus-resistant Arabidopsis thaliana accession. Although SEGS-2 is related to cassava genomic sequences, an earlier study showed that it occurs as episomes and is packaged into virions in CMD-infected cassava and viruliferous whiteflies. We identified SEGS-2 episomes in SEGS-2 transgenic Arabidopsis. The episomes occur as both double-stranded and single-stranded DNA, with the single-stranded form packaged into virions. In addition, SEGS-2 episomes replicate in tobacco protoplasts in the presence, but not the absence, of ACMV DNA-A. SEGS-2 episomes contain a SEGS-2 derived promoter and an open reading frame with the potential to encode a 75-amino acid protein. An ATG mutation at the beginning of the SEGS-2 coding region does not enhance ACMV infection in A. thaliana. Together, the results established that SEGS-2 is a new type of begomovirus satellite that enhances viral disease through the action of an SEGS-2-encoded protein that may also be encoded by the cassava genome. IMPORTANCE Cassava is an important root crop in the developing world and a food and income crop for more than 300 million African farmers. Cassava is rising in global importance and trade as the demands for biofuels and commercial starch increase. More than half of the world's cassava is produced in Africa, where it is primarily grown by smallholder farmers, many of whom are from the poorest villages. Although cassava can grow under high temperature, drought, and poor soil conditions, its production is severely limited by viral diseases. Cassava mosaic disease (CMD) is one of the most important viral diseases of cassava and can cause up to 100% yield losses. We provide evidence that SEGS-2, which was originally isolated from cassava crops displaying severe and atypical CMD symptoms in Tanzanian fields, is a novel begomovirus satellite that can compromise the development of durable CMD resistance.

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.

Kinship networks of seed exchange shape spatial patterns of plant virus diversity

By structuring farmers’ informal networks of seed exchange, kinship systems play a key role in the dynamics of crop genetic diversity in smallholder farming systems. However, because many crop diseases are propagated through infected germplasm, local seed systems can also facilitate the dissemination of seedborne pathogens. Here, we investigate how the interplay of kinship systems and local networks of germplasm exchange influences the metapopulation dynamics of viruses responsible for the cassava mosaic disease (CMD), a major threat to food security in Africa. Combining anthropological, genetic and plant epidemiological data, we analyzed the genetic structure of local populations of the African cassava mosaic virus (ACMV), one of the main causal agents of CMD. Results reveal contrasted patterns of viral diversity in patrilineal and matrilineal communities, consistent with local modes of seed exchange. Our results demonstrate that plant virus ecosystems have also a cultural component and that social factors that shape regional seed exchange networks influence the genetic structure of plant virus populations.

This study combines ethnobotanical and epidemiological data to understand how social networks of seed exchange influence the genetic structure of the African cassava mosaic virus in Gabon. Results reveal contrasted patterns of viral diversity in patrilineal and matrilineal communities, consistent with cultural differences in modes of seed exchange.

Population diversity of cassava mosaic begomoviruses increases over the course of serial vegetative propagation

Cassava mosaic disease (CMD) represents a serious threat to cassava, a major root crop for more than 300 million Africans. CMD is caused by single-stranded DNA begomoviruses that evolve rapidly, making it challenging to develop durable disease resistance. In addition to the evolutionary forces of mutation, recombination and reassortment, factors such as climate, agriculture practices and the presence of DNA satellites may impact viral diversity. To gain insight into the factors that alter and shape viral diversity in planta, we used high-throughput sequencing to characterize the accumulation of nucleotide diversity after inoculation of infectious clones corresponding to African cassava mosaic virus (ACMV) and East African cassava mosaic Cameroon virus (EACMCV) in the susceptible cassava landrace Kibandameno. We found that vegetative propagation had a significant effect on viral nucleotide diversity, while temperature and a satellite DNA did not have measurable impacts in our study. EACMCV diversity increased linearly with the number of vegetative propagation passages, while ACMV diversity increased for a time and then decreased in later passages. We observed a substitution bias toward C→T and G→A for mutations in the viral genomes consistent with field isolates. Non-coding regions excluding the promoter regions of genes showed the highest levels of nucleotide diversity for each genome component. Changes in the 5' intergenic region of DNA-A resembled the sequence of the cognate DNA-B sequence. The majority of nucleotide changes in coding regions were non-synonymous, most with predicted deleterious effects on protein structure, indicative of relaxed selection pressure over six vegetative passages. Overall, these results underscore the importance of knowing how cropping practices affect viral evolution and disease progression.

Evolutionary dynamics of bipartite begomoviruses revealed by complete genome analysis

Several key evolutionary events marked the evolution of geminiviruses, culminating with the emergence of divided (bipartite) genomes represented by viruses classified in the genus Begomovirus. This genus represents the most abundant group of multipartite viruses, contributing significantly to the observed abundance of multipartite species in the virosphere. Although aspects related to virus-host interactions and evolutionary dynamics have been extensively studied, the bipartite nature of these viruses has been little explored in evolutionary studies. Here, we performed a parallel evolutionary analysis of the DNA-A and DNA-B segments of New World begomoviruses. A total of 239 full-length DNA-B sequences obtained in this study, combined with 292 DNA-A and 76 DNA-B sequences retrieved from GenBank, were analysed. The results indicate that the DNA-A and DNA-B respond differentially to evolutionary processes, with the DNA-B being more permissive to variation and more prone to recombination than the DNA-A. Although a clear geographic segregation was observed for both segments, differences in the genetic structure between DNA-A and DNA-B were also observed, with cognate segments belonging to distinct genetic clusters. DNA-B coding regions evolve under the same selection pressures than DNA-A coding regions. Together, our results indicate an interplay between reassortment and recombination acting at different levels across distinct subpopulations and segments.

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.

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.

Metagenomic analyses and genetic diversity of Tomato leaf curl Arusha virus affecting tomato plants in Kenya

BACKGROUND

Tomato production is threatened worldwide by the occurrence of begomoviruses which are associated with tomato leaf curl diseases. There is little information on the molecular properties of tomato begomoviruses in Kenya, hence we investigated the population and genetic diversity of begomoviruses associated with tomato leaf curl in Kenya.

METHODS

Tomato leaf samples with virus-like symptoms were obtained from farmers' field across the country in 2018 and Illumina sequencing undertaken to determine the genetic diversity of associated begomoviruses. Additionally, the occurrence of selection pressure and recombinant isolates within the population were also evaluated.

RESULTS

Twelve complete begomovirus genomes were obtained from our samples with an average coverage of 99.9%. The sequences showed 95.7-99.7% identity among each other and 95.9-98.9% similarities with a Tomato leaf curl virus Arusha virus (ToLCArV) isolate from Tanzania. Analysis of amino acid sequences showed the highest identities in the regions coding for the coat protein gene (98.5-100%) within the isolates, and 97.1-100% identity with the C4 gene of ToLCArV. Phylogenetic algorithms clustered all Kenyan isolates in the same clades with ToLCArV, thus confirming the isolates to be a variant of the virus. There was no evidence of recombination within our isolates. Estimation of selection pressure within the virus population revealed the occurrence of negative or purifying selection in five out of the six coding regions of the sequences.

CONCLUSIONS

The begomovirus associated with tomato leaf curl diseases of tomato in Kenya is a variant of ToLCArV, possibly originating from Tanzania. There is low genetic diversity within the virus population and this information is useful in the development of appropriate management strategies for the disease in the country.

Intra-host evolution of the ssDNA virus tomato severe rugose virus (ToSRV)

To evaluate and quantify the evolutionary dynamics of the bipartite begomovirus tomato severe rugose virus (ToSRV) in a cultivated and a non-cultivated host, plants of tomato and Nicandra physaloides were biolistically inoculated with an infectious clone and systemically infected leaves were sampled at 30, 75 and 120 days after inoculation. Total DNA was extracted and sequenced in the Illumina HiSeq 2000 platform. The datasets were trimmed with the quality score limit set to 0.01, and the assembly was performed using the infectious clone sequence as reference. SNPs were filtered using a minimum p-value of 0.001 and the sum frequencies were used to calculate the deviation from the original clone sequence. Nucleotide substitution rates were calculated for the two DNA components in both hosts: 1.73 × 10^-3 and 3.07 × 10^-4 sub/site/year for the DNA-A and DNA-B, respectively, in N. physaloides, and 8.05 × 10^-4 and 7.02 × 10^-5 sub/site/year the for DNA-A and DNA-B, respectively, in tomato. These values are in the same range of those estimated for viruses with single-stranded RNA genomes and for other begomoviruses. Strikingly, the number of substitutions decreased over time, suggesting the presence of bottlenecks during systemic infection. Determination of Shannon's entropy indicated different patterns of variation in the DNA-A and the DNA-B, suggesting distinct evolutionary forces acting upon each component.

Molecular Characterization and Genomic Function of Grapevine Geminivirus A

A new grapevine geminivirus A (GGVA) isolate (named as GGVA-17YM1) and its associated defective genome (GGVA-D) were identified from a grapevine sample collected in Yuanmou, Yunnan Province, using sRNA high throughput sequencing and traditional Sanger sequencing. To explore the pathogenicity of GGVA and GGVA-D, infectious clones of GGVA-17YM1 and GGVA-D-17YM1 were constructed. Infection assays indicated that Nicotiana benthamiana plants inoculated with GGVA alone or a combination of GGVA and GGVA-D exhibited upward curled apical leaves and dwarfism. Southern blotting and quantitative real-time polymerase chain reaction analysis revealed that GGVA-D increased the accumulation level of GGVA DNA. Transient expression using a PVX-derived recombinant vector indicated that C2 and C4 encoded by GGVA are involved in symptom induction in N. benthamiana. Furthermore, the V2 protein inhibited local RNA silencing in co-infiltration assays in GFP transgenic N. benthamiana plants. Subsequently, full-length genome sequencing resulted in the identification of 11 different isolates of GGVA and 9 associated defective DNA molecules. Phylogenetic analysis based on whole genome sequences showed that all GGVA isolates, including our sequences, clustered into two distinct branches with no geographical grouping. Analyses of molecular variation indicated single nucleotide polymorphisms (SNPs) with more transitions (55.97%) than transversions (44.03%). Furthermore, the main variants for ORF C1, C3, or V1 were synonymous mutations, and non-synonymous mutations for ORF C2, C4, and V2. Genetic selection analysis indicated that negative selection acted on four ORFs (V1, C1, C2, and C3), while V2 and C4 were under positive selection. Our results contribute to the characterization of the genetic diversity of GGVA and provide insights into its pathogenicity.

Insertion and deletion evolution reflects antibiotics selection pressure in a Mycobacterium tuberculosis outbreak

In genome evolution, genetic variants are the source of diversity, which natural selection acts upon. Treatment of human tuberculosis (TB) induces a strong selection pressure for the emergence of antibiotic resistance-conferring variants in the infecting Mycobacterium tuberculosis (MTB) strains. MTB evolution in response to treatment has been intensively studied and mainly attributed to point substitutions. However, the frequency and contribution of insertions and deletions (indels) to MTB genome evolution remains poorly understood. Here, we analyzed a multi-drug resistant MTB outbreak for the presence of high-quality indels and substitutions. We find that indels are significantly enriched in genes conferring antibiotic resistance. Furthermore, we show that indels are inherited during the outbreak and follow a molecular clock with an evolutionary rate of 5.37e-9 indels/site/year, which is 23 times lower than the substitution rate. Inherited indels may co-occur with substitutions in genes along related biological pathways; examples are iron storage and resistance to second-line antibiotics. This suggests that epistatic interactions between indels and substitutions affect antibiotic resistance and compensatory evolution in MTB.

Arbuscular Mycorrhizal Symbiosis Primes Tolerance to Cucumber Mosaic Virus in Tomato

Tomato plants can establish symbiotic interactions with arbuscular mycorrhizal fungi (AMF) able to promote plant nutrition and prime systemic plant defenses against pathogens attack; the mechanism involved is known as mycorrhiza-induced resistance (MIR). However, studies on the effect of AMF on viral infection, still limited and not conclusive, indicate that AMF colonization may have a detrimental effect on plant defenses against viruses, so that the term "mycorrhiza-induced susceptibility" (MIS) has been proposed for these cases. To expand the case studies to a not yet tested viral family, that is, Bromoviridae, we investigated the effect of the colonization by the AMF Funneliformis mosseae on cucumber mosaic virus (CMV) infection in tomato by phenotypic, physiological, biochemical, and transcriptional analyses. Our results showed that the establishment of a functional AM symbiosis is able to limit symptoms development. Physiological and transcriptomic data highlighted that AMF mitigates the drastic downregulation of photosynthesis-related genes and the reduction of photosynthetic CO2 assimilation rate caused by CMV infection. In parallel, an increase of salicylic acid level and a modulation of reactive oxygen species (ROS)-related genes, toward a limitation of ROS accumulation, was specifically observed in CMV-infected mycorrhizal plants. Overall, our data indicate that the AM symbiosis influences the development of CMV infection in tomato plants and exerts a priming effect able to enhance tolerance to viral infection.

Evolutionary timescale and geographical movement of cucumber mosaic virus, with focus on Iranian strains

Cucumber mosaic virus (CMV) is a geographically widespread plant virus with a very broad host range. The virus has been detected in diverse crops all over Iran. In this study, we estimated the timescale of the evolution of CMV subgroup I and the geographical movement of the virus with a focus on Iranian strains. Analyses using the MP and CP genes and their concatenation revealed that the CMV population within subgroup I had a single ancestor dating back to about 450-550 years ago. The Iranian strains formed three clusters in a maximum-clade-credibility phylogenetic tree. It was found that the most recent common ancestor of the Iranian strains within each cluster dates back to less than 100 years ago. Our results also suggest that both short- and long-distance migration of Iranian CMV strains has occurred in the last 100 years.

A Bipartite Begomovirus Infecting Boerhavia erecta (Family Nyctaginaceae) in the Dominican Republic Represents a Distinct Phylogenetic Lineage and has a High Degree of Host Specificity

Boerhavia erecta plants in and around agricultural fields in the Azua Valley of the southeastern Dominican Republic often show striking golden mosaic symptoms. Leaf samples from B. erecta plants showing these symptoms were collected in 2012 and 2013, and PCR tests with degenerate primers revealed begomovirus DNA-A and DNA-B components. The complete sequences of the DNA-A and DNA-B components of four isolates show a high degree of sequence identity (>96%) and a genome organization typical of New World (NW) bipartite begomoviruses. Sequence comparisons and phylogenetic analyses revealed that these isolates composed a new phylogenetic lineage of NW bipartite begomoviruses. The most closely related begomovirus is Merremia mosaic virus, a weed-infecting species from Puerto Rico. Because DNA-A sequence identities are well below the 91% threshold, these isolates represent a new begomovirus species, for which the name Boerhavia golden mosaic virus (BoGMV) is proposed. Infectious cloned BoGMV DNA-A and DNA-B components induced golden mosaic symptoms in agroinoculated B. erecta plants, thereby fulfilling Koch's postulates for this disease. Agroinoculation and mechanical transmission experiments revealed that BoGMV has an unusually narrow host range, limited to members of the family Nyctaginaceae and not including the permissive host Nicotiana benthamiana. The inability of BoGMV to infect N. benthamiana was due to a deficiency in cell-to-cell movement but not to a unique amino acid residue in the movement protein.

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.

Disease management at the wildlife-livestock interface: Using whole-genome sequencing to study the role of elk in Mycobacterium bovis transmission in Michigan, USA

The role of wildlife in the persistence and spread of livestock diseases is difficult to quantify and control. These difficulties are exacerbated when several wildlife species are potentially involved. Bovine tuberculosis (bTB), caused by Mycobacterium bovis, has experienced an ecological shift in Michigan, with spillover from cattle leading to an endemically infected white-tailed deer (deer) population. It has potentially substantial implications for the health and well-being of both wildlife and livestock and incurs a significant economic cost to industry and government. Deer are known to act as a reservoir of infection, with evidence of M. bovis transmission to sympatric elk and cattle populations. However, the role of elk in the circulation of M. bovis is uncertain; they are few in number, but range further than deer, so may enable long distance spread. Combining Whole Genome Sequences (WGS) for M. bovis isolates from exceptionally well-observed populations of elk, deer and cattle with spatiotemporal locations, we use spatial and Bayesian phylogenetic analyses to show strong spatiotemporal admixture of M. bovis isolates. Clustering of bTB in elk and cattle suggests either intraspecies transmission within the two populations, or exposure to a common source. However, there is no support for significant pathogen transfer amongst elk and cattle, and our data are in accordance with existing evidence that interspecies transmission in Michigan is likely only maintained by deer. This study demonstrates the value of whole genome population studies of M. bovis transmission at the wildlife-livestock interface, providing insights into bTB management in an endemic system.

Resistance of advanced cassava breeding clones to infection by major viruses in Uganda

Cassava brown streak disease (CBSD) and cassava mosaic disease (CMD) are two viral diseases that cause severe yield losses in cassava of up to 100%, thereby persistently threatening food and income security in sub-Saharan Africa. For effective management of these diseases, there is a critical need to develop and deploy varieties with dual resistance to CBSD and CMD. In this study, we determined the response of advanced breeding lines to field infection by cassava brown streak viruses (CBSVs) and cassava mosaic begomoviruses (CMBs). This aim helped in identifying superior clones for downstream breeding. In total, 220 cassava clones, three in uniform yield trials (UYTs) and 217 in a crossing block trial (CBT), were evaluated for virus and disease resistance. Field data were collected on disease incidence and severity. To detect and quantify CBSVs, 448 and 128 leaf samples from CBSD symptomatic and symptomless plants were analyzed by reverse transcription PCR and real-time quantitative PCR, respectively. In addition, 93 leaf samples from CMD symptomatic plants in the CBT were analyzed by conventional PCR using CMB species-specific primers. In the CBT, 124 (57%) cassava clones did not express CMD symptoms. Of the affected plants, 44 (55%) had single African cassava mosaic virus infection. Single Cassava brown streak virus (CBSV) infections were more prevalent (81.6%) in CBT clones than single Ugandan cassava brown streak virus (UCBSV) infection (3.2%). Of the three advanced clones in the UYT, NAROCASS 1 and NAROCASS 2 had significantly lower (P < 0.05) CBSD severity, incidence, and CBSV load than MH04/0300. In the UYT, only 22% of samples tested had CBSVs, and all showed a negative result for CMBs. The low disease incidence, severity, and viral load associated with NAROCASS 1 and NAROCASS 2 is evidence of their tolerance to both CBSD and CMD. Therefore, these two cassava clones should be utilized in CBSD and CMD management in Uganda, including their utilization as progenitors in further virus resistance breeding.

Geographic Distribution and Genetic Diversity of Rice Stripe Mosaic Virus in Southern China

Rice stripe mosaic virus (RSMV) transmitted by the leafhopper Recilia dorsalis is a tentative new species in the genus Cytorhabdovirus identified recently in South China. To explore its geographic distribution and genetic diversity, field investigation and viral whole-genome sequencing were conducted in this study. The results indicated that RSMV was present in the rice samples collected across southern China. Twelve representative samples from different geographical regions were selected for viral whole-genome sequencing and the viral genome variation was analyzed in combination with a previously reported RSMV isolate. Identity analysis showed that the genome sequences of 13 RSMV isolates were highly conserved with nucleotide identities over 99.4%. There was a strong negative selection pressure during the evolution of RSMV with more transitions (72.08%) than transversions (27.92%) found between the RSMV isolates. Among the seven genes encoded by RSMV, the P gene was the most variable, followed by N, M, L, and G; the P3 and P6 amino acid sequences were not found to be mutated and no mutations were found in the non-coding region. A phylogenetic tree based on the RSMV whole-genome nucleotide sequences revealed that all RSMV isolates clustered in two groups based on geographical origin. Notably, the L proteins of the Guangxi and Hainan isolates had five and one specific amino acid sites, respectively, suggesting that the L gene has undergone environmental adaptive variation during the dispersal of RSMV.

Eukaryotic Circular Rep-Encoding Single-Stranded DNA (CRESS DNA) Viruses: Ubiquitous Viruses With Small Genomes and a Diverse Host Range

While single-stranded DNA (ssDNA) was once thought to be a relatively rare genomic architecture for viruses, modern metagenomics sequencing has revealed circular ssDNA viruses in most environments and in association with diverse hosts. In particular, circular ssDNA viruses encoding a homologous replication-associated protein (Rep) have been identified in the majority of eukaryotic supergroups, generating interest in the ecological effects and evolutionary history of circular Rep-encoding ssDNA viruses (CRESS DNA) viruses. This review surveys the explosion of sequence diversity and expansion of eukaryotic CRESS DNA taxonomic groups over the last decade, highlights similarities between the well-studied geminiviruses and circoviruses with newly identified groups known only through their genome sequences, discusses the ecology and evolution of eukaryotic CRESS DNA viruses, and speculates on future research horizons.

An unusually high substitution rate in transplant-associated BK polyomavirus in vivo is further concentrated in HLA-C-bound viral peptides

Infection with human BK polyomavirus, a small double-stranded DNA virus, potentially results in severe complications in immunocompromised patients. Here, we describe the in vivo variability and evolution of the BK polyomavirus by deep sequencing. Our data reveal the highest genomic evolutionary rate described in double-stranded DNA viruses, i.e., 10⁻³–10⁻⁵ substitutions per nucleotide site per year. High mutation rates in viruses allow their escape from immune surveillance and adaptation to new hosts. By combining mutational landscapes across viral genomes with in silico prediction of viral peptides, we demonstrate the presence of significantly more coding substitutions within predicted cognate HLA-C-bound viral peptides than outside. This finding suggests a role for HLA-C in antiviral immunity, perhaps through the action of killer cell immunoglobulin-like receptors. The present study provides a comprehensive view of viral evolution and immune escape in a DNA virus.

Little is known about the mechanisms of evolution and viral immune escape in double-stranded DNA (dsDNA) viruses. Here, we study the evolution of BK polyomavirus and observe the highest genomic evolutionary rate described so far for a dsDNA virus, in the range of RNA viruses, which usually evolve rapidly. Furthermore, the prediction of viral peptides to determine immune escape suggests a specific role of HLA-C in antiviral immunity. These findings are helpful for future advances in antiviral therapies and provide a step forward in our understanding of in vivo viral evolution in humans.

Virus discovery in all three major lineages of terrestrial arthropods highlights the diversity of single-stranded DNA viruses associated with invertebrates

Viruses encoding a replication-associated protein (Rep) within a covalently closed, single-stranded (ss)DNA genome are among the smallest viruses known to infect eukaryotic organisms, including economically valuable agricultural crops and livestock. Although circular Rep-encoding ssDNA (CRESS DNA) viruses are a widespread group for which our knowledge is rapidly expanding, biased sampling toward vertebrates and land plants has limited our understanding of their diversity and evolution. Here, we screened terrestrial arthropods for CRESS DNA viruses and report the identification of 44 viral genomes and replicons associated with specimens representing all three major terrestrial arthropod lineages, namely Euchelicerata (spiders), Hexapoda (insects), and Myriapoda (millipedes). We identified virus genomes belonging to three established CRESS DNA viral families (Circoviridae, Genomoviridae, and Smacoviridae); however, over half of the arthropod-associated viral genomes are only distantly related to currently classified CRESS DNA viral sequences. Although members of viral and satellite families known to infect plants (Geminiviridae, Nanoviridae, Alphasatellitidae) were not identified in this study, these plant-infecting CRESS DNA viruses and replicons are transmitted by hemipterans. Therefore, members from six out of the seven established CRESS DNA viral families circulate among arthropods. Furthermore, a phylogenetic analysis of Reps, including endogenous viral sequences, reported to date from a wide array of organisms revealed that most of the known CRESS DNA viral diversity circulates among invertebrates. Our results highlight the vast and unexplored diversity of CRESS DNA viruses among invertebrates and parallel findings from RNA viral discovery efforts in undersampled taxa.

Whitefly-transmitted viruses threatening cassava production in Africa

Emerging plant viruses are one of the greatest problems facing crop production worldwide, and have severe consequences in the developing world where subsistence farming is a major source of food production, and knowledge and resources for management are limited. In Africa, evolution of two viral disease complexes, cassava mosaic begomoviruses (CMBs) (Geminiviridae) and cassava brown streak viruses (CBSVs) (Potyviridae), have resulted in severe pandemics that continue to spread and threaten cassava production. Identification of genetically diverse and rapidly evolving CMBs and CBSVs, extensive genetic variation in the vector, Bemisia tabaci (Hemiptera: Aleyrodidae), and numerous secondary endosymbiont profiles that influence vector phenotypes suggest that complex local and regional vector-virus-plant-environment interactions may be driving the evolution and epidemiology of these viruses.

The diversity, evolution and epidemiology of plant viruses: A phylogenetic view

During the past four decades, the scientific community has seen an exponential advance in the number, sophistication, and quality of molecular techniques and bioinformatics tools for the genetic characterization of plant virus populations. Predating these advances, the field of Phylogenetics has significantly contributed to understand important aspects of plant virus evolution. This review aims at summarizing the impact of Phylogenetics in the current knowledge on three major aspects of plant virus evolution that have benefited from the development of phylogenetic inference: (1) The identification and classification of plant virus diversity. (2) The mechanisms and forces shaping the evolution of plant virus populations. (3) The understanding of the interaction between plant virus evolution, epidemiology and ecology. The work discussed here highlights the important role of phylogenetic approaches in the study of the dynamics of plant virus populations.

Differential Shape of Geminivirus Mutant Spectra Across Cultivated and Wild Hosts With Invariant Viral Consensus Sequences

Geminiviruses (family Geminiviridae) possess single-stranded circular DNA genomes that are replicated by cellular polymerases in plant host cell nuclei. In their hosts, geminivirus populations behave as ensembles of mutant and recombinant genomes, known as viral quasispecies. This favors the emergence of new geminiviruses with altered host range, facilitating new or more severe diseases or overcoming resistance traits. In warm and temperate areas several whitefly-transmitted geminiviruses of the genus Begomovirus cause the tomato yellow leaf curl disease (TYLCD) with significant economic consequences. TYLCD is frequently controlled in commercial tomatoes by using the dominant Ty-1 resistance gene. Over a 45 day period we have studied the diversification of three begomoviruses causing TYLCD: tomato yellow leaf curl virus (TYLCV), tomato yellow leaf curl Sardinia virus (TYLCSV) and tomato yellow leaf curl Malaga virus (TYLCMaV, a natural recombinant between TYLCV and TYLCSV). Viral quasispecies resulting from inoculation of geminivirus infectious clones were examined in plants of susceptible tomato (ty-1/ty-1), heterozygous resistant tomato (Ty-1/ty-1), common bean, and the wild reservoir Solanum nigrum. Differences in virus fitness across hosts were observed while viral consensus sequences remained invariant. However, the complexity and heterogeneity of the quasispecies were high, especially in common bean and the wild host. Interestingly, the presence or absence of the Ty-1 allele in tomato did not lead to differences in begomovirus mutant spectra. However, the fitness decrease of TYLCSV and TYLCV in tomato at 45 dpi might be related to an increase in CP (Coat protein) mutation frequency. In Solanum nigrum the recombinant TYLCMaV, which showed lower fitness than TYLCSV, at 45 dpi actively explored Rep (Replication associated protein) ORF but not the overlapping C4. Our results underline the importance of begomovirus mutant spectra during infections. This is especially relevant in the wild reservoir of the viruses, which has the potential to maintain highly diverse mutant spectra without modifying their consensus sequences.

Evolutionary Determinants of Host and Vector Manipulation by Plant Viruses

Plant viruses possess adaptations for facilitating acquisition, retention, and inoculation by vectors. Until recently, it was hypothesized that these adaptations are limited to virus proteins that enable virions to bind to vector mouthparts or invade their internal tissues. However, increasing evidence suggests that viruses can also manipulate host plant phenotypes and vector behaviors in ways that enhance their own transmission. Manipulation of vector-host interactions occurs through virus effects on host cues that mediate vector orientation, feeding, and dispersal behaviors, and thereby, the probability of virus transmission. Effects on host phenotypes vary by pathosystem but show a remarkable degree of convergence among unrelated viruses whose transmission is favored by the same vector behaviors. Convergence based on transmission mechanism, rather than phylogeny, supports the hypothesis that virus effects are adaptive and not just by-products of infection. Based on this, it has been proposed that viruses manipulate hosts through multifunctional proteins that facilitate exploitation of host resources and elicitation of specific changes in host phenotypes. But this proposition is rarely discussed in the context of the numerous constraints on virus evolution imposed by molecular and environmental factors, which figure prominently in research on virus-host interactions not dealing with host manipulation. To explore the implications of this oversight, we synthesized available literature to identify patterns in virus effects among pathogens with shared transmission mechanisms and discussed the results of this synthesis in the context of molecular and environmental constraints on virus evolution, limitations of existing studies, and prospects for future research.

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 variation of tomato yellow leaf curl virus in the insect vector Bemisia tabaci

Insect vectors play significant roles in geminivirus spread and evolution in nature. To date little is known about the population dynamics of begomoviruses in their insect vector Bemisia tabaci. In this study we analyzed the genetic variation of tomato yellow leaf curl virus (TYLCV) in its host plant, Solanum lycopersicum, in its transmission vector B. tabaci raised on TYLCV-infected S. lycopersicum plants, and in B. tabaci after being transferred from S. lycopersicum to Gossypium hirsutum. We found that the levels of variability of TYLCV remained stable in S. lycopersicum plants, but increased significantly in both invasive and indigenous species of B. tabaci. We also presented evidence that the elevated mutation frequencies in TYLCV populations from vector whiteflies were caused mainly by mutations that occurred at several distinct sites within the TYLCV genome. Simultaneous introduction of mutations in the hot spots did not affect the ability of TYLCV to be transmitted by B. tabaci, but reduced its pathogenicity in both S. lycopersicum and Nicotiana benthamiana. Our findings provide new information on population variability of TYLCV in its insect vector, extending the knowledge of the influence of insect vector on plant virus population dynamics.

The Evolutionary History and Spatiotemporal Dynamics of the NC Lineage of Citrus Tristeza Virus

Citrus tristeza virus (CTV) is a major pathogen affecting citrus trees worldwide. However, few studies have focused on CTV's evolutionary history and geographic behavior. CTV is locally dispersed by an aphid vector and long distance dispersion due to transportation of contaminated material. With the aim to delve deeper into the CTV-NC (New Clade) genotype evolution, we estimated an evolution rate of 1.19 × 10^-3 subs/site/year and the most common recent ancestor in 1977. Furthermore, the place of origin of the genotype was in the United States, and a great expansion of the population was observed in Uruguay. This expansion phase could be a consequence of the increment in the number of naïve citrus trees in Uruguayan orchards encompassing citrus industry growth in the past years.

Evolutionary dynamics and genomic features of the Elizabethkingia anophelis 2015 to 2016 Wisconsin outbreak strain

An atypically large outbreak of Elizabethkingia anophelis infections occurred in Wisconsin. Here we show that it was caused by a single strain with thirteen characteristic genomic regions. Strikingly, the outbreak isolates show an accelerated evolutionary rate and an atypical mutational spectrum. Six phylogenetic sub-clusters with distinctive temporal and geographic dynamics are revealed, and their last common ancestor existed approximately one year before the first recognized human infection. Unlike other E. anophelis, the outbreak strain had a disrupted DNA repair mutY gene caused by insertion of an integrative and conjugative element. This genomic change probably contributed to the high evolutionary rate of the outbreak strain and may have increased its adaptability, as many mutations in protein-coding genes occurred during the outbreak. This unique discovery of an outbreak caused by a naturally occurring mutator bacterial pathogen provides a dramatic example of the potential impact of pathogen evolutionary dynamics on infectious disease epidemiology.

Elizabethkingia anophelis is an emerging pathogen of high antimicrobial resistance. Perrin and colleagues sequenced isolates of a 2015/2016 E. anophelis outbreak in Wisconsin and found substantial genetic diversity, accelerated evolutionary rate and a disruptive mutation in the DNA repair gene mutY.