Widely conserved recombination patterns among single-stranded DNA viruses

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

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

Bacilladnaviridae: refined taxonomy and new insights into the biology and evolution of diatom-infecting DNA viruses

Bacilladnaviruses are single-stranded DNA viruses that infect diatoms that, so far, have been primarily identified in marine organisms and environments. Using a viral metagenomics approach, we discovered 13 novel bacilladnaviruses originating from samples of mud-flat snail (Amphibola crenata; n=3 genomes) and benthic sediments (n=10 genomes) collected from the Avon-Heathcote Estuary in New Zealand. Comparative genomics and phylogenetic analysis of the new bacilladnavirus sequences in the context of the previously classified members of the family helped refine and further expand the Bacilladnaviridae taxonomy. Here, based on the replication-associated protein phylogeny and pairwise identities, we established 4 new genera - Aberdnavirus, Keisodnavirus, Puahadnavirus and Seawadnavirus - and 13 new species within the family. Comparison of the bacilladnavirus capsid protein sequences suggests that the positively charged N-terminal region (R-arm) is required for encapsidation of the larger genomes, whereas the smaller bacilladnavirus genomes can be packaged in the absence of the R-arm subdomain. Furthermore, analysis of the bacilladnavirus genomes revealed that members of three genera encode a highly derived variant of a phospholipase A1, which is predicted to be involved in the lysis of the infected diatoms and/or facilitates the entry of the virions into the host cells. Collectively, our results allow refining of the taxonomy of bacilladnaviruses and provide new insights into the biology and evolution of this understudied group of diatom viruses.

Diverse Anelloviruses Identified in Leporids from Arizona (USA)

The communities of viruses studied in rabbits and hares (family Leporidae) have largely been those with clinical significance. Consequently, less is known broadly about other leporid viruses. Anelloviruses (family Anelloviridae) are likely commensal members of the single-stranded DNA virome in mammals. Here, we employ a viral metagenomic approach to identify DNA viruses of leporids and the ticks feeding on them in Arizona, USA. We characterize five complete anellovirus genomes from four leporids belonging to the black-tailed jackrabbit (Lepus californicus, n = 3) and the desert cottontail (Sylvilagus audubonii, n = 1). All five anellovirus genomes share > 69% orf1 gene pairwise nucleotide identity with those found in Iberian hares and thus belong to the species Aleptorquevirus lepor1. Accordingly, we expand the known host range of this anellovirus species to include Iberian hares in Europe and black-tailed jackrabbit and desert cottontail in the USA. We also sequenced the complete mitochondrial genomes of the four leporid hosts (Sylvilagus audubonii, n = 1; Lepus californicus, n = 3) and two ticks (Dermacentor parumapertus, n = 2) found feeding on two black-tailed jackrabbits. These results expand the diversity of anelloviruses in leporids while giving insight into the host genetics of leporids and ticks in Arizona, USA.

Novel lineage of anelloviruses with large genomes identified in dolphins

Anellovirus infections are ubiquitous in mammals but lack any clear disease association, suggesting a commensal virus-host relationship. Although anelloviruses have been identified in numerous mammalian hosts, their presence in members of the family Delphinidae has yet to be reported. Here, using a metagenomic approach, we characterize complete anellovirus genomes (n = 69) from four Delphinidae host species: short-finned pilot whale (Globicephala macrorhynchus, n = 19), killer whale (Orcinus orca, n = 9), false killer whale (Pseudorca crassidens, n = 6), and pantropical spotted dolphin (Steno attenuatus, n = 1). Sequence comparison of the open reading frame 1 (ORF1) encoding the capsid protein, the only conserved gene shared by all anelloviruses, shows that the Delphinidae anelloviruses form a novel genus-level clade that encompasses 22 unique species-level groupings. We provide evidence that different Delphinidae species can be co-infected by multiple anelloviruses belonging to distinct species groupings. Notably, the ORF1 protein of the Delphinidae anelloviruses is considerably larger than those encoded by all previously described anelloviruses from other hosts (spanning 14 vertebrate orders and including 27 families). Comprehensive analysis of the ORF1 sequences and predicted protein structures showed that the increased size of these proteins results from divergent elaborations within the capsid-distal P2 subdomain and elongation of the C-terminal domain of ORF1. Comparative structural and phylogenetic analyses suggest that acquisition of the P2 subdomain and its diversification occurred convergently in the anelloviruses associated with primate and Delphinidae hosts. Collectively, our results further the appreciation of diversity and evolution of the ubiquitous and enigmatic viruses in the family Anelloviridae.

IMPORTANCE

Anelloviruses are ubiquitous in mammals, but their infection has not yet been linked to any disease, suggesting a commensal virus-host relationship. Here, we describe the first anelloviruses associated with diverse species of dolphins. The dolphinid anelloviruses represent a new genus (tentatively named "Qoptorquevirus") and encode open reading frame 1 (ORF1) (capsid) proteins that are considerably larger than those encoded by previously described anelloviruses from other hosts. Comprehensive analysis of the ORF1 sequences and predicted protein structures revealed the underlying structural basis for such an extravagant ORF1 size and suggested that ORF1 size increased convergently in the anelloviruses associated with primate and Delphinidae hosts, respectively. Collectively, our results provide insights into the diversity and evolution of Anelloviridae. Further exploration of the anellovirus diversity, especially in the host species that have not yet been sampled, is expected to further clarify their evolutionary trajectory and explain the unusual virus-host commensal relationship.

Expanding the genomic diversity of human anelloviruses

Anelloviruses are a group of small, circular, single-stranded DNA viruses that are found ubiquitously across mammalian hosts. Here, we explored a large number of publicly available human microbiome datasets and retrieved a total of 829 anellovirus genomes, substantially expanding the known diversity of these viruses. The majority of new genomes fall within the three major human anellovirus genera: Alphatorquevirus, Betatorquevirus, and Gammatorquevirus, while we also present new genomes of the under-sampled Hetorquevirus, Memtorquevirus, and Samektorquevirus genera. We performed recombination analysis and show evidence of extensive recombination across all human anelloviruses. Interestingly, more than 95% of the detected events are between members of the same genus and only 15 inter-genus recombination events were detected. The breakpoints of recombination cluster in hotspots at the ends and outside of the ORF1 gene, while a recombination coldspot was detected within the gene. Our analysis suggests that anellovirus evolution is governed by homologous recombination; however, events between distant viruses or ones producing chimaeric ORF1s likely lead to nonviable recombinants. The large number of genomes further allowed us to examine how essential genomic features vary across anelloviruses. These include functional domains in the ORF1 protein and the nucleotide motif of the replication loop region, required for the viruses' rolling-circle replication. A subset of the genomes assembled in both this and previous studies are completely lacking these essential elements, opening up the possibility that anellovirus intracellular populations contain nonstandard viral genomes. However, low-read depth of the metagenomically assembled contigs may partly explain the lack of some features. Overall, our study highlights key features of anellovirus genomics and evolution, a largely understudied group of viruses whose potential in virus-based therapeutics is recently being explored.

Phylogeographic analysis of Begomovirus coat and replication-associated proteins

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

Circular Single-Stranded DNA: Discovery, Biological Effects, and Applications

The field of nucleic acid therapeutics has witnessed a significant surge in recent times, as evidenced by the increasing number of approved genetic drugs. However, current platform technologies containing plasmids, lipid nanoparticle-mRNAs, and adeno-associated virus vectors encounter various limitations and challenges. Thus, we are devoted to finding a novel nucleic acid vector and have directed our efforts toward investigating circular single-stranded DNA (CssDNA), an ancient form of nucleic acid. CssDNAs are ubiquitous, but generally ignored. Accumulating evidence suggests that CssDNAs possess exceptional properties as nucleic acid vectors, exhibiting great potential for clinical applications in genetic disorders, gene editing, and immune cell therapy. Here, we comprehensively review the discovery and biological effects of CssDNAs as well as their applications in the field of biomedical research for the first time. Undoubtedly, as an ancient form of DNA, CssDNA holds immense potential and promises novel insights for biomedical research.

Investigating pigeon circovirus infection in a pigeon farm: molecular detection, phylogenetic analysis and complete genome analysis

BACKGROUND

Pigeon circovirus infections in pigeons (Columba livia domestica) have been reported worldwide. Pigeons should be PiCV-free when utilized as qualified experimental animals. However, pigeons can be freely purchased as experimental animals without any clear guidelines to follow. Herein, we investigated the status quo of PiCV infections on a pigeon farm in Beijing, China, which provides pigeons for experimental use.

RESULTS

PiCV infection was verified in at least three types of tissues in all forty pigeons tested. A total of 29 full-length genomes were obtained and deposited in GenBank. The whole genome sequence comparison among the 29 identified PiCV strains revealed nucleotide homologies of 85.8-100%, and these sequences exhibited nucleotide homologies of 82.7-98.9% as compared with those of the reference sequences. The cap gene displayed genetic diversity, with a wide range of amino acid homologies ranging from 64.5% to 100%. Phylogenetic analysis of the 29 full-genome sequences revealed that the PiCV strains in this study could be further divided into four clades: A (17.2%), B (10.4%), C (37.9%) and D (34.5%). Thirteen recombination events were also detected in 18 out of the 29 PiCV genomes obtained in this study. Phylogenetic research using the rep and cap genes verified the recombination events, which occurred between clades A/F, A/B, C/D, and B/D among the 18 PiCV strains studied.

CONCLUSIONS

In conclusion, PiCV infection, which is highly genetically varied, is extremely widespread on pigeon farms in Beijing. These findings indicate that if pigeons are to be used as experimental animals, it is necessary to evaluate the impact of PiCV infection on the results.

Detection of circovirus in free-ranging brown rats (Rattus norvegicus)

Accidentally found, two poisoned brown rats from Hungary were surveyed for presence of circoviral DNA, using specific nested primers, designed against the rep gene of the virus. Both specimens were positive. The whole genomes were amplified using inverse PCR based on the Rep sequence parts and sequenced by the primer walking method. Genomic analyses revealed that these novel rat viruses, together with tawny owl-associated circovirus reported by Italian researchers in 2022, are sequence variations of the same virus from genus Circovirus. In phylogenetic reconstructions, these circovirus strains detected from brown rats clustered closest to circoviruses derived from faeces samples of various predatory mammals. Molecular data as well as the phylogenetic analyses of the complete derived replication-associated protein and the capsid protein, as well as the prey preference of the host species of the recently described tawny owl-associated virus suggest that brown rat could be the evolutionary adapted host of the viruses described in this paper (brown rat circovirus types 1 and 2) and the previously reported tawny owl-associated virus. Possible pathogenic and zoonotic role of these viruses need further studies.

Unusual outbreaks of curly top disease in processing tomato fields in northern California in 2021 and 2022 were caused by a rare strain of beet curly top virus and facilitated by extreme weather events

In the western United States, curly top disease (CTD) is caused by beet curly top virus (BCTV). In California, CTD causes economic loss to processing tomato production in central and southern areas but, historically, not in the north. Here, we document unusual CTD outbreaks in processing tomato fields in the northern production area in 2021 and 2022, and show that these were caused by the rare spinach curly top strain (BCTV-SpCT). These outbreaks were associated with proximity of fields to foothills and unusually hot, dry, and windy spring weather conditions, possibly by altering migrations of the beet leafhopper (BLH) vector from locations with BCTV-SpCT reservoirs. Support for this hypothesis came from the failure to observe CTD outbreaks and BLH migrations in 2023, when spring weather conditions were cool and wet. Our results show the climate-induced emergence of a rare plant virus strain to cause an economically important disease in a new crop and location.

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.

Prevalence, genomic characteristics, and pathogenicity of fowl adenovirus 2 in Southern China

In recent years, the occurrence of fowl adenovirus 2 (FAdV-2) has been on the rise in China, posing a significant threat to the poultry industry. This study aimed to investigate the epidemiology, phylogenetic relationship, genomic characteristics, and pathogenicity of FAdV-2. The epidemiological analysis revealed the detection of multiple FAdV serotypes, including FAdV-1, FAdV-2, FAdV-3, FAdV-4, FAdV-8a, FAdV-8b, and FAdV-11 serotypes. Among them, FAdV-2 exhibited the highest proportion, accounting for 21.05% (8/38). The complete genomes of these 8 FAdV-2 strains were sequenced. Genetic evolution analysis indicated that these FAdV-2 strains formed a separate branch within the FAdV-D group, sharing 94.60 to 97.90% nucleotide similarity with the reference FAdV-2 and FAdV-11 strains. Notably, the recombination analysis revealed that 5 out of the 8 FAdV-2 strains, exhibited recombination events between FAdV-2 and FAdV-11. The recombination regions involved Hexon, Fiber, ORF19 genes and 3' end. Furthermore, pathogenicity experiments demonstrated that recombinant FAdV-2 XX strain is capable of inducing mortality rate of 66.70% and causing more severe hepatitis hydropericardium syndrome (HHS) in 6-wk-old specific-pathogen-free chickens. These findings contribute to our understanding of the prevalence, genomic characteristics, and the pathogenicity of FAdV-2, providing foundations for FAdV-2 vaccine development.

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.

Specific Nucleotides in the Common Region of the Begomovirus Tomato Rugose Mosaic Virus (ToRMV) Are Responsible for the Negative Interference over Tomato Severe Rugose Virus (ToSRV) in Mixed Infection

Mixed infection between two or more begomoviruses is commonly found in tomato fields and can affect disease outcomes by increasing symptom severity and viral accumulation compared with single infection. Viruses that affect tomato include tomato severe rugose virus (ToSRV) and tomato rugose mosaic virus (ToRMV). Previous work showed that in mixed infection, ToRMV negatively affects the infectivity and accumulation of ToSRV. ToSRV and ToRMV share a high degree of sequence identity, including cis-elements in the common region (CR) and their specific recognition sites (iteron-related domain, IRD) within the Rep gene. Here, we investigated if divergent sites in the CR and IRD are involved in the interaction between these two begomoviruses. ToSRV clones were constructed containing the same nucleotides as ToRMV in the CR (ToSRV-A(ToR:CR)), IRD (ToSRV-A(ToR:IRD)) and in both regions (ToSRV-A(ToR:CR+IRD)). When plants were co-inoculated with ToRMV and ToSRV-A(ToR:IRD), the infectivity and accumulation of ToSRV were negatively affected. In mixed inoculation of ToRMV with ToSRV-A(ToR:CR), high infectivity of both viruses and high DNA accumulation of ToSRV-A(ToR:CR) were observed. A decrease in viral accumulation was observed in plants inoculated with ToSRV-A(ToR:CR+IRD). These results indicate that differences in the CR, but not the IRD, are responsible for the negative interference of ToRMV on ToSRV.

The phylogenetic and phylogeographic landscape of the beak and feather disease virus, 1996-2022

The beak and feather disease virus (BFDV), causative agent of Psittacine beak and feather disease (PBFD), is a highly fatal and widespread virus that infects both the wild and captive Psittaciformes around the world. The BFDV genome is a ssDNA of approximately 2 kb in size, making it among the smallest known pathogenic viruses. Though, the virus is placed in Circoviridae family of the Circovirus genus, there is no classification system on clade and sub-clade level according to the International Committee on Taxonomy of Viruses and the strains are grouped on the bases of geographic locations. Thus, we provide the latest and robust phylogenetic classification of BFDVs in this study based on full-length genomic sequences, grouping all the available 454 strains detected during 1996-2022 into two distinct clades, e.g., GI and GII. The GI clade is further divided into six sub-clades (GI a-f), while GII into two sub-clades (GII a and b). In addition, the phylogeographic network identified high variability among the BFDV strains, showing several branches, where all the branches are connected to four strains, e.g., BFDV-ZA-PGM-70A(GenBank ID: HM748921.1, 2008-South Africa), BFDV-ZA-PGM-81A(GenBank ID: JX221009.1, 2008-South Africa), BFDV14(GenBank ID: GU015021.1, 2010-Thailand) and BFDV-isolate-9IT11(GenBank ID: KF723390.1, 2014-Italy). Furthermore, we identified 27 recombination events in the rep (replication-associated protein) and cap (capsid protein) coding regions using the complete genomes of BFDVs. Similarly, the amino acids variability analysis indicated that both the rep and cap regions are highly variable with values exceeding the variability coefficient estimation limit of 1.00, speculating the possible amino acids drift with the emergence of new strains. The findings provided in this study may offer the latest phylogenetic, phylogeographic and evolutionary landscape of the BFDVs.

Non-structural genes of novel lemur adenoviruses reveal codivergence of virus and host

Adenoviruses (AdVs) are important human and animal pathogens and are frequently used as vectors for gene therapy and vaccine delivery. Surprisingly, there are only scant data regarding primate AdV origin and evolution, especially in the most basal primate hosts. We detect and sequence AdVs from faeces of two Madagascan lemur species. Complete genome sequence analyses define a new AdV species with a particularly large gene encoding a protein of unknown function in the early gene region 3. Unexpectedly, the new AdV species is not most similar to human or other simian AdVs but to bat adenovirus C. Genome characterisation shows signals of virus-host codivergence in non-structural genes, which show lower diversity than structural genes. Outside a lemur species mixing zone, recombination less frequently separates structural genes, as in human adenovirus C. The evolutionary history of lemur AdVs likely involves both a host switch and codivergence with the lemur hosts.

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.

The Role of Extensive Recombination in the Evolution of Geminiviruses

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

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

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

Epidemiology, genetic diversity, and association of canine circovirus infection in dogs with respiratory disease

Although canine circovirus (CanineCV)-associated with gastroenteritis has been well documented, the virus is also detectable in the respiratory discharge of dogs with respiratory disease. In this study, an epidemiological approach was used to explore the association between the presence of CanineCV and respiratory symptoms in dogs. Respiratory swabs were collected from 76 healthy dogs and 114 dogs with respiratory illness and tested for CanineCV using conventional PCR (cPCR). Furthermore, lung tissues collected from 15 necropsied dogs showing pneumonia were tested using the real-time PCR (qPCR) and in situ hybridization (ISH) technique. A total of 8.95% (17/190) of dogs were CanineCV positive, with a significant association (p = 0.013) in dogs with respiratory signs. Four necropsied dogs were qPCR positive with the CanineCV-DNA labeling localized in tracheobronchial lymphoid cells (3/4), pulmonary parenchyma, capillary endothelia, and mononuclear cells harboring in alveoli (2/4). Full-length genome sequences of seven CanineCV strains were analyzed, indicating that the detected CanineCV genome clustered in the CanineCV-4 genotype. Genetic recombination was also evident in the replicase (Rep) gene. Although the role of CanineCV primarily affecting lung lesions could not be determined from this study, the presence of CanineCV DNA in pulmonary-associated cells indicated the potential association of the virus with canine respiratory disease; thus, linking causality must be examined in further studies.

Conserved recombination patterns across coronavirus subgenera

Recombination contributes to the genetic diversity found in coronaviruses and is known to be a prominent mechanism whereby they evolve. It is apparent, both from controlled experiments and in genome sequences sampled from nature, that patterns of recombination in coronaviruses are non-random and that this is likely attributable to a combination of sequence features that favour the occurrence of recombination break points at specific genomic sites, and selection disfavouring the survival of recombinants within which favourable intra-genome interactions have been disrupted. Here we leverage available whole-genome sequence data for six coronavirus subgenera to identify specific patterns of recombination that are conserved between multiple subgenera and then identify the likely factors that underlie these conserved patterns. Specifically, we confirm the non-randomness of recombination break points across all six tested coronavirus subgenera, locate conserved recombination hot- and cold-spots, and determine that the locations of transcriptional regulatory sequences are likely major determinants of conserved recombination break-point hotspot locations. We find that while the locations of recombination break points are not uniformly associated with degrees of nucleotide sequence conservation, they display significant tendencies in multiple coronavirus subgenera to occur in low guanine-cytosine content genome regions, in non-coding regions, at the edges of genes, and at sites within the Spike gene that are predicted to be minimally disruptive of Spike protein folding. While it is apparent that sequence features such as transcriptional regulatory sequences are likely major determinants of where the template-switching events that yield recombination break points most commonly occur, it is evident that selection against misfolded recombinant proteins also strongly impacts observable recombination break-point distributions in coronavirus genomes sampled from nature.

Molecular detection and phylogenetic analysis of pigeon circovirus from racing pigeons in Northern China

Background

Pigeon circovirus (PiCV) infections in pigeons (Columba livia) have been reported worldwide. Currently, pigeon racing is becoming increasingly popular and considered to be a national sport in China, and even, the greatest competitions of racing pigeons are taking place in China. However, there are still no epidemiologic data regarding PiCV infections among racing pigeons in China. The purpose of our study was to provide information of prevalence, genetic variation and evolution of PiCV from racing pigeons in China.

Results

To trace the prevalence, genetic variation and evolution of PiCV in sick and healthy racing pigeons, 622 samples were collected from 11 provinces or municipalities in China from 2016 to 2019. The results showed that the positive rate of PiCV was 19.3% (120/622) at the sample level and 59.0% (23/39) at the club level, thus suggesting that the virus was prevalent in Chinese racing pigeons. A sequence analysis revealed that the cap genes of the PiCV strains identified in our study displayed a high genetic diversity and shared nucleotide homologies of 71.9%–100% and amino acid homologies of 71.7%–100%. 28 and 36 unique amino acid substitutions were observed in the Cap and Rep proteins derived from our PiCV strains, respectively. A cladogram representation of PiCV strains phylogeny based on 90 cap gene sequences showed that the strains in this study could be further divided into seven clades (A, B, C, E, G, H, and I) and some of them were closely related to worldwide strains from different types of pigeons. A large number of recombination events (31 events) were also detected in the PiCV genomes from Chinese racing pigeons.

Conclusions

These findings indicate that PiCV strains circulating in China exhibit a high genetic diversity and also contribute to information of prevalence, genetic variation and evolution of PiCV from racing pigeons in China.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08425-8.

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

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

WGS- versus ORF5-Based Typing of PRRSV: A Belgian Case Study

Porcine reproductive and respiratory syndrome virus (PRRSV) is the causative agent of one of the most widespread and economically devastating diseases in the swine industry. Typing circulating PRRSV strains by means of sequencing is crucial for developing adequate control strategies. Most genetic studies only target the highly variable open reading frame (ORF) 5, for which an extensive database is available. In this study, we performed whole-genome sequencing (WGS) on a collection of 124 PRRSV-1 positive serum samples that were collected over a 5-year period (2015–2019) in Belgium. Our results show that (nearly) complete PRRSV genomes can be obtained directly from serum samples with a high success rate. Analysis of the coding regions confirmed the exceptionally high genetic diversity, even among Belgian PRRSV-1 strains. To gain more insight into the added value of WGS, we performed phylogenetic cluster analyses on separate ORF datasets as well as on a single, concatenated dataset (CDS) containing all ORFs. A comparison between the CDS and ORF clustering schemes revealed numerous discrepancies. To explain these differences, we performed a large-scale recombination analysis, which allowed us to identify a large number of potential recombination events that were scattered across the genome. As PRRSV does not contain typical recombination hot-spots, typing PRRSV strains based on a single ORF is not recommended. Although the typing accuracy can be improved by including multiple regions, our results show that the full genetic diversity among PRRSV strains can only be captured by analysing (nearly) complete genomes. Finally, we also identified several vaccine-derived recombinant strains, which once more raises the question of the safety of these vaccines.

Redondovirus Diversity and Evolution on Global, Individual, and Molecular Scales

Redondoviridae is a newly established family of circular Rep-encoding single-stranded (CRESS) DNA viruses found in the human ororespiratory tract. Redondoviruses were previously found in ∼15% of respiratory specimens from U.S. urban subjects; levels were elevated in individuals with periodontitis or critical illness. Here, we report higher redondovirus prevalence in saliva samples: four rural African populations showed 61 to 82% prevalence, and an urban U.S. population showed 32% prevalence. Longitudinal, limiting-dilution single-genome sequencing revealed diverse strains of both redondovirus species (Brisavirus and Vientovirus) in single individuals, persistence over time, and evidence of intergenomic recombination. Computational analysis of viral genomes identified a recombination hot spot associated with a conserved potential DNA stem-loop structure. To assess the possible role of this site in recombination, we carried out in vitro studies which showed that this potential stem-loop was cleaved by the virus-encoded Rep protein. In addition, in reconstructed reactions, a Rep-DNA covalent intermediate was shown to mediate DNA strand transfer at this site. Thus, redondoviruses are highly prevalent in humans, found in individuals on multiple continents, heterogeneous even within individuals and encode a Rep protein implicated in facilitating recombination. IMPORTANCE Redondoviridae is a recently established family of DNA viruses predominantly found in the human respiratory tract and associated with multiple clinical conditions. In this study, we found high redondovirus prevalence in saliva from urban North American individuals and nonindustrialized African populations in Botswana, Cameroon, Ethiopia, and Tanzania. Individuals on both continents harbored both known redondovirus species. Global prevalence of both species suggests that redondoviruses have long been associated with humans but have remained undetected until recently due to their divergent genomes. By sequencing single redondovirus genomes in longitudinally sampled humans, we found that redondoviruses persisted over time within subjects and likely evolve by recombination. The Rep protein encoded by redondoviruses catalyzes multiple reactions in vitro, consistent with a role in mediating DNA replication and recombination. In summary, we identify high redondovirus prevalence in humans across multiple continents, longitudinal heterogeneity and persistence, and potential mechanisms of redondovirus evolution by recombination.

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.

The Begomovirus Species Melon Chlorotic Leaf Curl Virus is Composed of Two Highly Divergent Strains that Differ in Their Genetic and Biological Properties

Since the early 1990s, squash production in Costa Rica has been affected by a whitefly-transmitted disease characterized by stunting and yellow mottling of leaves. The squash yellow mottle disease (SYMoD) was shown to be associated with a bipartite begomovirus, originally named squash yellow mild mottle virus (SYMMoV). It was subsequently established that SYMMoV is a strain of melon chlorotic leaf curl virus (MCLCuV), a bipartite begomovirus that causes a chlorotic leaf curl disease of melons in Guatemala. In the present study, the complete sequences of the DNA-A and DNA-B components of a new isolate of the strain MCLCuV-Costa Rica (MCLCuV-CR) were determined. Comparisons of full-length DNA-A sequences revealed 97% identity with a previously characterized isolate of MCLCuV-CR and identities of 90 to 91% with those of isolates of the strain MCLCuV-Guatemala (MCLCuV-GT), which is below or at the current begomovirus species demarcation threshold of 91%. A more extensive analysis of the MCLCuV-CR and -GT sequences revealed substantial divergence in both components and different histories of recombination for the DNA-A components. The cloned full-length DNA-A and DNA-B components of this new MCLCuV-CR isolate were infectious and induced SYMoD in a range of squashes and in pumpkin, thereby fulfilling Koch's postulates for this disease. However, in contrast to MCLCuV-GT, MCLCuV-CR induced mild symptoms in watermelon and no symptoms in melon and cucumber. Taken together, our results indicate that MCLCuV-CR and -GT have substantially diverged, genetically and biologically, and have evolved to cause distinct diseases of different cucurbit crops. Taxonomically, these viruses are at the strain/species boundary, but retain the designation as strains of Melon chlorotic leaf curl virus under current International Committee on Taxonomy guidelines.

High molecular diversity and divergent subpopulations of the begomovirus cnidoscolus mosaic leaf deformation virus associated with Cnidoscolus urens

Begomoviruses have circular, single-stranded DNA genomes encapsidated into twinned quasi-icosahedral particles and are transmitted by whiteflies of the Bemisia tabaci sibling group. Begomoviruses infect cultivated and non-cultivated plants, causing great losses in economically important crops worldwide. To better understand the genetic diversity of begomoviruses infecting the non-cultivated host Cnidoscolus urens, leaf samples exhibiting virus-like symptoms were collected in different localities in the state of Alagoas, Brazil, during 2015 and 2016. Forty-two complete DNA-A sequences were cloned and sequenced by the Sanger method. Based on nucleotide sequence comparisons, the 42 new isolates were identified as the bipartite begomovirus cnidoscolus mosaic leaf deformation virus (CnMLDV). The CnMLDV isolates were clustered in two phylogenetic groups (clusters I and II) corresponding to their sampling areas, and the high value of Wright's F fixation index observed for the DNA-A sequences suggests population structuring. At least seven independent intraspecies recombination events were predicted among CnMLDV isolates, with recombination breakpoints located in the common region (CR) and in the CP and Rep genes. Also, a high per site nucleotide diversity (π) was observed for CnMLDV isolates, with CP being significantly more variable than Rep. Despite the high genetic variability, strong negative or purifying selection was identified as the main selective force acting upon CP and Rep.

Global genome analysis reveals a vast and dynamic anellovirus landscape within the human virome

Anelloviruses are a ubiquitous component of healthy human viromes and remain highly prevalent after being acquired early in life. The full extent of "anellome" diversity and its evolutionary dynamics remain unexplored. We employed in-depth sequencing of blood-transfusion donor(s)-recipient pairs coupled with public genomic resources for a large-scale assembly of anellovirus genomes and used the data to characterize global and personal anellovirus diversity through time. The breadth of the anellome is much greater than previously appreciated, and individuals harbor unique anellomes and transmit lineages that can persist for several months within a diverse milieu of endemic host lineages. Anellovirus sequence diversity is shaped by extensive recombination at all levels of divergence, hindering traditional phylogenetic analyses. Our findings illuminate the transmission dynamics and vast diversity of anelloviruses and set the foundation for future studies to characterize their biology.

Illuminating the Virosphere Through Global Metagenomics

Viruses are the most abundant biological entity on Earth, infect cellular organisms from all domains of life, and are central players in the global biosphere. Over the last century, the discovery and characterization of viruses have progressed steadily alongside much of modern biology. In terms of outright numbers of novel viruses discovered, however, the last few years have been by far the most transformative for the field. Advances in methods for identifying viral sequences in genomic and metagenomic datasets, coupled to the exponential growth of environmental sequencing, have greatly expanded the catalog of known viruses and fueled the tremendous growth of viral sequence databases. Development and implementation of new standards, along with careful study of the newly discovered viruses, have transformed and will continue to transform our understanding of microbial evolution, ecology, and biogeochemical cycles, leading to new biotechnological innovations across many diverse fields, including environmental, agricultural, and biomedical sciences.

Discovery and Characterization of a Novel Bipartite Botrexvirus From the Phytopathogenic Fungus Botryosphaeria dothidea

In this study, we describe a novel positive, single-stranded (+ss) RNA mycovirus, named Botryosphaeria dothidea botrexvirus 1 (BdBV1), from a phytopathogenic fungus Botryosphaeria dothidea showing abnormal morphology and attenuated virulence. BdBV1 is phylogenetically related to Botrytis virus X (BotVX) and is the second potential member of the proposed genus Botrexvirus in the family Alphaflexiviridae. However, it differs from the monopartite BotVX in that BdBV1 possesses a bipartite genome comprised of two ssRNA segments (RNA1 and RNA2 with lengths of 5,035 and 1,063 nt, respectively). BdBV1 RNA1 and RNA2 encode putative RNA-dependent RNA polymerase (RdRp) and coat protein (CP) genes, which share significant identity with corresponding genes in both fungal and plant viruses. Moreover, open reading frames (ORFs) 2-4 of BdBV1 RNA1 shared no detectable identity with any known viral proteins. Immunosorbent electron microscopy (ISEM) analysis using an antibody against the virus CP generated in vitro revealed that BdBV1 is encapsidated in filamentous particles. A comparison of the biological effects of BdBV1 infection on symptoms and growth in isogenic lines of virus-free and virus-infected B. dothidea revealed that BdBV1 is probably involved in reduced growth and virulence of the host fungus. This study describes and characterizes a novel bipartite botrexvirus, which is closely related to uni- and multi-partite fungal and plant viruses and contributes useful information to a better understanding of virus evolution.

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.

How Do Pathogens Evolve Novel Virulence Activities?

This article is part of the Top 10 Unanswered Questions in MPMI invited review series.We consider the state of knowledge on pathogen evolution of novel virulence activities, broadly defined as anything that increases pathogen fitness with the consequence of causing disease in either the qualitative or quantitative senses, including adaptation of pathogens to host immunity and physiology, host species, genotypes, or tissues, or the environment. The evolution of novel virulence activities as an adaptive trait is based on the selection exerted by hosts on variants that have been generated de novo or arrived from elsewhere. In addition, the biotic and abiotic environment a pathogen experiences beyond the host may influence pathogen virulence activities. We consider host-pathogen evolution, host range expansion, and external factors that can mediate pathogen evolution. We then discuss the mechanisms by which pathogens generate and recombine the genetic variation that leads to novel virulence activities, including DNA point mutation, transposable element activity, gene duplication and neofunctionalization, and genetic exchange. In summary, if there is an (epi)genetic mechanism that can create variation in the genome, it will be used by pathogens to evolve virulence factors. Our knowledge of virulence evolution has been biased by pathogen evolution in response to major gene resistance, leaving other virulence activities underexplored. Understanding the key driving forces that give rise to novel virulence activities and the integration of evolutionary concepts and methods with mechanistic research on plant-microbe interactions can help inform crop protection.[Formula: see text] Copyright © 2021 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

Global Patterns of Recombination across Human Viruses

Viral recombination is a major evolutionary mechanism driving adaptation processes, such as the ability of host-switching. Understanding global patterns of recombination could help to identify underlying mechanisms and to evaluate the potential risks of rapid adaptation. Conventional approaches (e.g., those based on linkage disequilibrium) are computationally demanding or even intractable when sequence alignments include hundreds of sequences, common in viral data sets. We present a comprehensive analysis of recombination across 30 genomic alignments from viruses infecting humans. In order to scale the analysis and avoid the computational limitations of conventional approaches, we apply newly developed topological data analysis methods able to infer recombination rates for large data sets. We show that viruses, such as ZEBOV and MARV, consistently displayed low levels of recombination, whereas high levels of recombination were observed in Sarbecoviruses, HBV, HEV, Rhinovirus A, and HIV. We observe that recombination is more common in positive single-stranded RNA viruses than in negatively single-stranded RNA ones. Interestingly, the comparison across multiple viruses suggests an inverse correlation between genome length and recombination rate. Positional analyses of recombination breakpoints along viral genomes, combined with our approach, detected at least 39 nonuniform patterns of recombination (i.e., cold or hotspots) in 18 viral groups. Among these, noteworthy hotspots are found in MERS-CoV and Sarbecoviruses (at spike, Nucleocapsid and ORF8). In summary, we have developed a fast pipeline to measure recombination that, combined with other approaches, has allowed us to find both common and lineage-specific patterns of recombination among viruses with potential relevance in viral adaptation.

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

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

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

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

New World Cactaceae Plants Harbor Diverse Geminiviruses

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

Redondoviridae and periodontitis: a case–control study and identification of five novel redondoviruses from periodontal tissues

Redondoviridae is a family of DNA viruses recently identified in the human oro-respiratory tract. However, the characteristics of this new virus family are not yet fully understood. The aim of the present study was to investigate the relationship between redondoviruses and chronic periodontitis. In addition, the complete circular genome, phylogenetic relationship, and biological characteristics of novel redondoviruses were analyzed. The gingival tissues of healthy individuals (n = 120) and periodontitis patients (n = 120) were analyzed using nested polymerase chain reaction assays. The prevalence of redondovirus infection in the periodontitis group was 71.67%. Logistic regression analysis revealed an association between redondoviruses and chronic periodontitis after controlling the confounding factors (odds ratio = 2.53). Five novel redondoviruses, named ‘human periodontal circular-like virus (HPeCV)’, were identified in patients with periodontitis and detailed genetic analysis of the viruses was performed. The 3,035–3,056 bp genome contained a capsid protein, a replication-associated protein, an open reading frame 3 protein, and a stem-loop structure. Phylogenetic analysis demonstrated that HPeCV-1, HPeCV-10, and HPeCV-25 formed a cluster. Recombination may be common in the genomes of HPeCVs. Potential antigenic epitopes in the capsid protein, which may be involved in the host immune response, were predicted. In conclusion, periodontitis patients had a significantly higher prevalence of redondoviruses than healthy controls. Genetic characterization enhanced the current understanding of the genetic diversity and pathogenicity of redondoviruses as well as their association with periodontitis in humans. The data presented in this article will expand the current understanding of the epidemiology, genetic diversity, and pathogenicity of redondoviruses.

Deformed wing virus type a and b in managed honeybee colonies of Argentina

Apis mellifera is infected by more than 24 virus species worldwide, mainly positive-sense, single-stranded RNA viruses of the Dicistroviridae and Iflaviridae families. Among the viruses that infect honey bees, Deformed wing virus is the most prevalent and is present as three master variants DWV-A, B, and C. Given that the ectoparasitic mite Varroa destructor vectors these virus variants, recombination events between them are expected, and variants and their recombinants can co-exist in mites and honeybees at the same time. In this study, we detect, through RT-qPCR, the presence of DWV-A and B in the same samples of adult bees from colonies of Argentina. Total RNA was extracted from pools of ten adult bees from 45 apiaries distributed across the main beekeeping Provinces of Argentina (Buenos Aires, Santa Fe, Córdoba, Santiago del Estero, Río Negro, and Mendoza); then RT-qPCR reactions were performed to detect DWV-A and B, with specific primer pairs. After the amplifications, PCR products (204 and 660 bp amplicons for DWV-B, and ~250 bp for DWV-A) were purified and sequenced to verify that they corresponded to reported sequences, analyzing them using the Blast software. Of the 45 samples analyzed by RT-qPCR, over 90% were infected with DWV-A and 47% were also positive for DWV-B, where it was found in high prevalence specifically in colonies of A. mellifera of the Buenos Aires Province. Future studies will determine the impact of this type of the virus and its ability to recombine with the other DWV types in the apiaries of our country.

Analysis of the beak and feather disease viral genome indicates evidence of multiple introduction events into Saudi Arabia

Psittacine beak and feather disease (PBFD), caused by beak and feather disease virus (BFDV) is a highly contagious disease in wild and captive psittacine populations and has an almost global presence. However, the BFDV infection in Saudi Arabia remains largely unknown. In the present study, we report the full genome sequence of BFDV strains from Saudi Arabia and its genetic diversity. The complete genome sequences were analyzed for 14 BFDV-infected birds representing 6 psittacine species. The complete genome sequence of BFDV strains was compared with 201 previously reported sequences to evaluate their diversity and possible recombination events, if any. Our analysis revealed that newly sequenced BFDV genomes from Saudi Arabia belonged to six different strains. Phylogenetic analysis suggested that the isolated BFDV genomes were highly recombinant with a high degree of diversity. It is evident from the study that psittacine species in Saudi Arabia are at risk from the spread of BFDV. As per the CITES trade database, about 190,000 parrots have been imported to Saudi Arabia since 1975 over a thousand instances. Presumably, during any of these trade events or unregulated trade of birds has predisposed the introduction of BFDV to Saudi Arabia. Understanding the epidemiology of BFDV is necessitated to address the threat posed by the virus to the psittacine population of Saudi Arabia.

Recombination and Positive Selection Differentially Shaped the Diversity of Betacoronavirus Subgenera

The Betacoronavirus genus of mammal-infecting viruses includes three subgenera (Sarbecovirus, Embecovirus, and Merbecovirus), in which most known human coronaviruses, including SARS-CoV-2, cluster. Coronaviruses are prone to host shifts, with recombination and positive selection possibly contributing to their high zoonotic potential. We analyzed the role of these two forces in the evolution of viruses belonging to the Betacoronavirus genus. The results showed that recombination has been pervasive during sarbecovirus evolution, and it is more widespread in this subgenus compared to the other two. In both sarbecoviruses and merbecoviruses, recombination hotspots are clearly observed. Conversely, positive selection was a less prominent force in sarbecoviruses compared to embecoviruses and merbecoviruses and targeted distinct genomic regions in the three subgenera, with S being the major target in sarbecoviruses alone. Overall, the results herein indicate that Betacoronavirus subgenera evolved along different trajectories, which might recapitulate their host preferences or reflect the origins of the presently available coronavirus sequences.

Bovine Coronavirus: Variability, Evolution, and Dispersal Patterns of a No Longer Neglected Betacoronavirus

Bovine coronavirus (BoCV) is an important pathogen of cattle, causing severe enteric disease and playing a role in the bovine respiratory disease complex. Similar to other coronaviruses, a remarkable variability characterizes both its genome and biology. Despite their potential relevance, different aspects of the evolution of BoCV remain elusive. The present study reconstructs the history and evolution of BoCV using a phylodynamic approach based on complete genome and spike protein sequences. The results demonstrate high mutation and recombination rates affecting different parts of the viral genome. In the spike gene, this variability undergoes significant selective pressures—particularly episodic pressure—located mainly on the protein surface, suggesting an immune-induced selective pressure. The occurrence of compensatory mutations was also identified. On the contrary, no strong evidence in favor of host and/or tissue tropism affecting viral evolution has been proven. The well-known plasticity is thus ascribable to the innate broad viral tropism rather than mid- or long-term adaptation. The evaluation of the geographic spreading pattern clearly evidenced two clusters: a European cluster and an American–Asian cluster. While a relatively dense and quick migration network was identified in the former, the latter was dominated by the primary role of the United States (US) as a viral exportation source. Since the viral spreading pattern strongly mirrored the cattle trade, the need for more intense monitoring and preventive measures cannot be underestimated as well as the need to enforce the vaccination of young animals before international trade, to reduce not only the clinical impact but also the transferal and mixing of BoCV strains.

Insights into the genetic diversity, recombination, and systemic infections with evidence of intracellular maturation of hepadnavirus in cats

Hepatitis B virus (HBV) is a human pathogen of global concern, while a high diversity of viruses related to HBV have been discovered in other animals during the last decade. Recently, the novel mammalian hepadnavirus, tentatively named domestic cat hepadnavirus (DCH), was detected in an immunocompromised cat. Herein, a collection of 209 cat sera and 15 hepato-diseased cats were screened for DCH using PCR, resulting in 12.4% and 20% positivity in the tested sera and necropsied cats, respectively. Among the DCH-positive sera, a significantly high level of co-detection with retroviral infection was found, with the highest proportion being co-detection with feline immunodeficiency virus (FIV). Full-length genome characterization of DCH revealed the genetic diversity between the nine Thai DCH sequences obtained, and that they phylogenetically formed three distinct monophyletic clades. A putative DCH recombinant strain was found, suggesting a possible role of recombination in DCH evolution. Additionally, quantitative PCR was used to determine the viral copy number in various organs of the DCH-moribund cats, while the pathological findings were compared to the viral localization in hepatocytes, adjacent to areas of hepatic fibrosis, by immunohistochemical (IHC) and western blot analysis. In addition to the liver, positive-DCH immunoreactivity was found in various other organs, including kidneys, lung, heart, intestine, brain, and lymph nodes, providing evidence of systemic infection. Ultrastructure of infected cells revealed electron-dense particles in the nucleus and cytoplasm of hepatocytes, bronchial epithelial cells, and fibroblasts. We propose the intracellular development mechanism of this virus. Although the definitive roles of pathogenicity of DCH remains undetermined, a contributory role of the virus associated with systemic diseases is possible.

A Novel Divergent Geminivirus Identified in Asymptomatic New World Cactaceae Plants

Cactaceae comprise a diverse and iconic group of flowering plants which are almost exclusively indigenous to the New World. The wide variety of growth forms found amongst the cacti have led to the trafficking of many species throughout the world as ornamentals. Despite the evolution and physiological properties of these plants having been extensively studied, little research has focused on cactus-associated viral communities. While only single-stranded RNA viruses had ever been reported in cacti, here we report the discovery of cactus-infecting single-stranded DNA viruses. These viruses all apparently belong to a single divergent species of the family Geminiviridae and have been tentatively named Opuntia virus 1 (OpV1). A total of 79 apparently complete OpV1 genomes were recovered from 31 different cactus plants (belonging to 20 different cactus species from both the Cactoideae and Opuntioideae clades) and from nine cactus-feeding cochineal insects (Dactylopius sp.) sampled in the USA and Mexico. These 79 OpV1 genomes all share > 78.4% nucleotide identity with one another and < 64.9% identity with previously characterized geminiviruses. Collectively, the OpV1 genomes display evidence of frequent recombination, with some genomes displaying up to five recombinant regions. In one case, recombinant regions span ~40% of the genome. We demonstrate that an infectious clone of an OpV1 genome can replicate in Nicotiana benthamiana and Opuntia microdasys. In addition to expanding the inventory of viruses that are known to infect cacti, the OpV1 group is so distantly related to other known geminiviruses that it likely represents a new geminivirus genus. It remains to be determined whether, like its cactus hosts, its geographical distribution spans the globe.

Free to Circulate: An Update on the Epidemiological Dynamics of Porcine Circovirus 2 (PCV-2) in Italy Reveals the Role of Local Spreading, Wild Populations, and Foreign Countries

Porcine circovirus 2 (PCV-2) is one of the most impactful and widespread pathogens of the modern swine industry. Unlike other DNA viruses, PCV-2 is featured by a remarkable genetic variability, which has led to the emergence and recognition of different genotypes, some of which (PCV-2a, 2b, and 2d) have alternated over time. Currently, PCV-2d is considered the most prevalent genotype, and some evidence of differential virulence and vaccine efficacy have been reported. Despite the potential practical relevance, the data on PCV-2 epidemiology in Italy are quite outdated and do not quantify the actual circulation of this genotype in Italy. In the present study, 82 complete ORF2 sequences were obtained from domestic pigs and wild boars sampled in Northern Italy in the period 2013–2018 and merged with those previously obtained from Italy and other countries. A combination of phylogenetic, haplotype network, and phylodynamic analyses were used to genotype the collected strains and evaluate the temporal trend and the spatial and host spread dynamics. A rising number of PCV-2d detections was observed in domestic pigs, particularly since 2013, reaching a detection frequency comparable to PCV-2b. A similar picture was observed in wild boars, although a lower sequence number was available. Overall, the present study demonstrates the extreme complexity of PCV-2 molecular epidemiology in Italy, the significant spread across different regions, the recurrent introduction from foreign countries, and the frequent occurrence of recombination events. Although a higher viral flux occurred from domestic to wild populations than vice versa, wild boars seem to maintain PCV-2 infection and spread it over relatively long distances.

Genotyping Porcine Circovirus 3 (PCV-3) Nowadays: Does It Make Sense?

The discovery of a globally distributed porcine circovirus (Porcine circovirus 3; PCV-3) has led to intense research activity and the production of a large amount of molecular data. Different research groups have proposed several, not always concordant, genotypes for this virus. While such categories could aid an easier interpretation of PCV-3 molecular epidemiology, any classification, to be useful in practical settings, must be univocal and of help in the understanding of underlying biological features and epidemiology. Based on these premises, the possibility of defining PCV-3 genotypes was evaluated on the broadest available dataset of PCV-3 complete genome (n = 357) and open reading frame 2 (ORF2, n = 653) sequences. Genetic distance and phylogenetic clustering were selected as the main objective criteria. Additional factors, including the number of within-cluster sequences, host and geographic clustering, concordance between different genomic regions, and analysis method were also taken in account to generate a classification that could be effectively applied in research and diagnostic settings. A maximum within-genotype genetic distance of 3% at the complete genome and 6% at the ORF2 levels, bootstrap support higher than 90%, and concordance between analysis methods allowed us to clearly define two clades which could be potentially defined as genotypes. Further subdivision was not suggested due to the absence of a meaningful association between PCV-3 and its biological/epidemiological features. Nevertheless, since one of the clades included two strains only, thus far we formally propose the definition of only one PCV-3 genotype (PCV-3a). The established criteria will allow us to automatically recognize other genotypes when more strain sequences are characterized.

Evolution and phylogeography of Culex pipiens densovirus

Viruses of the Parvoviridae family infect a wide range of animals including vertebrates and invertebrates. So far, our understanding of parvovirus diversity is biased towards medically or economically important viruses mainly infecting vertebrate hosts, while invertebrate infecting parvoviruses—namely densoviruses—have been largely neglected. Here, we investigated the prevalence and the evolution of the only mosquito-infecting ambidensovirus, Culex pipiens densovirus (CpDV), from laboratory mosquito lines and natural populations collected worldwide. CpDV diversity generally grouped in two clades, here named CpDV-1 and -2. The incongruence of the different gene trees for some samples suggested the possibility of recombination events between strains from different clades. We further investigated the role of selection on the evolution of CpDV genome and detected many individual sites under purifying selection both in non-structural and structural genes. However, some sites in structural genes were under diversifying selection, especially during the divergence of CpDV-1 and -2 clades. These substitutions between CpDV-1 and -2 clades were mostly located in the capsid protein encoding region and might cause changes in host specificity or pathogenicity of CpDV strains from the two clades. However, additional functional and experimental studies are necessary to fully understand the protein conformations and the resulting phenotype of these substitutions between clades of CpDV.

Fowl Adenovirus (FAdV) Recombination with Intertypic Crossovers in Genomes of FAdV-D and FAdV-E, Displaying Hybrid Serological Phenotypes

After analyzing 27 new genomes from fowl adenovirus (FAdV) field isolates and so-far unsequenced prototypes, we report the first evidence for recombination in FAdVs. Recombination was confined to species FAdV-D and FAdV-E, accommodating the largest number of, and the intraspecies-wise most differentiated, types. The majority of detected events occurred in FAdV-E, involving segments with parental origin of all constitutive types. Together with the diversity of breakpoints, this suggests widespread recombination in this species. With possible constraints through species-specific genes and diversification patterns, the recombinogenic potential of FAdVs attains particular interest for inclusion body hepatitis (IBH), an important disease in chickens, caused by types from the recombination-prone species. Autonomously evolving, recombinant segments were associated with major sites under positive selection, among them the capsid protein hexon and fiber genes, the right-terminal ORFs 19, 25, and the ORF20/20A family. The observed mosaicism in genes indicated as targets of adaptive pressures points toward an immune evasion strategy. Intertypic hexon/fiber-recombinants demonstrated hybrid neutralization profiles, retrospectively explaining reported controversies on reference strains B3-A, T8-A, and X11-A. Furthermore, cross-neutralization supported sequence-based evidence for interdomain recombination in fiber and contributed to a tentatively new type. Overall, our findings challenge the purported uniformity of types responsible for IBH, urging more complete identification strategies for FAdVs. Finally, important consequences arise for in vivo studies investigating cross-protection against IBH.