Evidence of pervasive biologically functional secondary structures within the genomes of eukaryotic single-stranded DNA viruses

The enigmatic roles of Anelloviridae and Redondoviridae in humans

Anelloviridae and Redondoviridae are virus families with small, circular, single-stranded DNA genomes that are common components of the human virome. Despite their small genome size of less than 5000 bases, they are remarkably successful - anelloviruses colonize over 90% of adult humans, while the recently discovered redondoviruses have been found at up to 80% prevalence in some populations. Anelloviruses are present in blood and many organs, while redondoviruses are found mainly in the ororespiratory tract. Despite their high prevalence, little is known about their biology or pathogenic potential. In this review, we discuss anelloviruses and redondoviruses and explore their enigmatic roles in human health and disease.

Primordial Capsid and Spooled ssDNA Genome Structures Unravel Ancestral Events of Eukaryotic Viruses

Marine algae viruses are important for controlling microorganism communities in the marine ecosystem and played fundamental roles during the early events of viral evolution. Here, we have focused on one major group of marine algae viruses, the single-stranded DNA (ssDNA) viruses from the Bacilladnaviridae family. We present the capsid structure of the bacilladnavirus Chaetoceros tenuissimus DNA virus type II (CtenDNAV-II), determined at 2.4-Å resolution. A structure-based phylogenetic analysis supported the previous theory that bacilladnaviruses have acquired their capsid protein via horizontal gene transfer from a ssRNA virus. The capsid protein contains the widespread virus jelly-roll fold but has additional unique features; a third β-sheet and a long C-terminal tail. Furthermore, a low-resolution reconstruction of the CtenDNAV-II genome revealed a partially spooled structure, an arrangement previously only described for dsRNA and dsDNA viruses. Together, these results exemplify the importance of genetic recombination for the emergence and evolution of ssDNA viruses and provide important insights into the underlying mechanisms that dictate genome organization. IMPORTANCE Single-stranded DNA (ssDNA) viruses are an extremely widespread group of viruses that infect diverse hosts from all three domains of life, consequently having great economic, medical, and ecological importance. In particular, bacilladnaviruses are highly abundant in marine sediments and greatly influence the dynamic appearance and disappearance of certain algae species. Despite the importance of ssDNA viruses and the last couple of years' advancements in cryo-electron microscopy, structural information on the genomes of ssDNA viruses remains limited. This paper describes two important achievements: (i) the first atomic structure of a bacilladnavirus capsid, which revealed that the capsid protein gene presumably was acquired from a ssRNA virus in early evolutionary events; and (ii) the structural organization of a ssDNA genome, which retains a spooled arrangement that previously only been observed for double-stranded viruses.

An extended APOBEC3A mutation signature in cancer

APOBEC mutagenesis, a major driver of cancer evolution, is known for targeting TpC sites in DNA. Recently, we showed that APOBEC3A (A3A) targets DNA hairpin loops. Here, we show that DNA secondary structure is in fact an orthogonal influence on A3A substrate optimality and, surprisingly, can override the TpC sequence preference. VpC (non-TpC) sites in optimal hairpins can outperform TpC sites as mutational hotspots. This expanded understanding of APOBEC mutagenesis illuminates the genomic Twin Paradox, a puzzling pattern of closely spaced mutation hotspots in cancer genomes, in which one is a canonical TpC site but the other is a VpC site, and double mutants are seen only in trans, suggesting a two-hit driver event. Our results clarify this paradox, revealing that both hotspots in these twins are optimal A3A substrates. Our findings reshape the notion of a mutation signature, highlighting the additive roles played by DNA sequence and DNA structure.

Characterization of two newly emerged torque teno sus virus isolates from a large-scale pig farm in China, in 2018

Torque teno sus virus (TTSuV) infection is common in China's pig herd. Although of uncertain pathogenicity, TTSuVs have been reported as a worsening factor of other porcine diseases, including porcine circovirus associated disease (PCVAD), porcine respiratory diseases complex (PRDC) or porcine dermatitis and nephropathy syndrome (PDNS). To better understand the genetic diversity in TTSuVs, the complete genomes of two newly emerged isolates, referred to as HeN1-A9 and HeN1-A11, collected from pig samples at a large-scale pig farm in China, were analyzed. Phylogenetic relationships of TTSuV sequences separated TTSuV1 and TTSuVk2a groups and divided TTSuV1 into two major subtypes, including TTSuV1a and TTSuV1b; HeN1-A9 and HeN1-A11 strains classified into the TTSuV1a subtype. Recombination analysis demonstrated HeN1-A9 and HeN1-A11 were generated via recombination in the overlapping ORF1/ORF3 region of TTSuV1a genome, which we report for the first time. Furthermore, we found that HeN1-A9 could be replicated in cultured MARC-145 cells for 18 passages. Our findings may be useful for elucidating the characteristics and epidemic status of TTSuVs in China.

Evolutionary Analyses of Base-Pairing Interactions in DNA and RNA Secondary Structures

Pairs of nucleotides within functional nucleic acid secondary structures often display evidence of coevolution that is consistent with the maintenance of base-pairing. Here, we introduce a sequence evolution model, MESSI (Modeling the Evolution of Secondary Structure Interactions), that infers coevolution associated with base-paired sites in DNA or RNA sequence alignments. MESSI can estimate coevolution while accounting for an unknown secondary structure. MESSI can also use graphics processing unit parallelism to increase computational speed. We used MESSI to infer coevolution associated with GC, AU (AT in DNA), GU (GT in DNA) pairs in noncoding RNA alignments, and in single-stranded RNA and DNA virus alignments. Estimates of GU pair coevolution were found to be higher at base-paired sites in single-stranded RNA viruses and noncoding RNAs than estimates of GT pair coevolution in single-stranded DNA viruses. A potential biophysical explanation is that GT pairs do not stabilize DNA secondary structures to the same extent that GU pairs do in RNA. Additionally, MESSI estimates the degrees of coevolution at individual base-paired sites in an alignment. These estimates were computed for a SHAPE-MaP-determined HIV-1 NL4-3 RNA secondary structure. We found that estimates of coevolution were more strongly correlated with experimentally determined SHAPE-MaP pairing scores than three nonevolutionary measures of base-pairing covariation. To assist researchers in prioritizing substructures with potential functionality, MESSI automatically ranks substructures by degrees of coevolution at base-paired sites within them. Such a ranking was created for an HIV-1 subtype B alignment, revealing an excess of top-ranking substructures that have been previously identified as having structure-related functional importance, among several uncharacterized top-ranking substructures.

Redondoviridae, a Family of Small, Circular DNA Viruses of the Human Oro-Respiratory Tract Associated with Periodontitis and Critical Illness

The global virome is largely uncharacterized but is now being unveiled by metagenomic DNA sequencing. Exploring the human respiratory virome, in particular, can provide insights into oro-respiratory diseases. Here, we use metagenomics to identify a family of small circular DNA viruses-named Redondoviridae-associated with human diseases. We first identified two redondovirus genomes from bronchoalveolar lavage samples from human lung donors. We then queried thousands of metagenomic samples and recovered 17 additional complete redondovirus genomes. Detections were exclusively in human samples and mostly from respiratory tract and oro-pharyngeal sites, where Redondoviridae was the second most prevalent eukaryotic DNA virus family. Redondovirus sequences were associated with periodontal disease, and abundances decreased with treatment. Some critically ill patients in a medical intensive care unit were found to harbor high levels of redondoviruses in respiratory samples. These results suggest that redondoviruses colonize human oro-respiratory sites and can bloom in several human disorders.

Argonaute bypasses cellular obstacles without hindrance during target search

Argonaute (Ago) proteins are key players in both gene regulation (eukaryotes) and host defense (prokaryotes). Acting on single-stranded nucleic-acid substrates, Ago relies on base pairing between a small nucleic-acid guide and its complementary target sequences for specificity. To efficiently scan nucleic-acid chains for targets, Ago diffuses laterally along the substrate and must bypass secondary structures as well as protein barriers. Using single-molecule FRET in conjunction with kinetic modelling, we reveal that target scanning is mediated through loose protein-nucleic acid interactions, allowing Ago to slide short distances over secondary structures, as well as to bypass protein barriers via intersegmental transfer. Our combined single-molecule experiment and kinetic modelling approach may serve as a platform to dissect search processes and study the effect of sequence on search kinetics for other nucleic acid-guided proteins.

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.

Recombinant Goose Circoviruses Circulating in Domesticated and Wild Geese in Poland

Circoviruses are circular single-stranded DNA (ssDNA) viruses that infect a variety of animals, both domestic and wild. Circovirus infection in birds is associated with immunosuppression and this in turn predisposes the infected animals to secondary infections that can lead to mortality. Farmed geese (Anser anser) in many parts of the world are infected with circoviruses. The majority of the current genomic information for goose circoviruses (GoCVs) (n = 40) are from birds sampled in China and Taiwan, and only two genome sequences are available from Europe (Germany and Poland). In this study, we sampled 23 wild and 19 domestic geese from the Gopło Lake area in Poland. We determined the genomes of GoCV from 21 geese; 14 domestic Greylag geese (Anser anser), three wild Greylag geese (A. anser), three bean geese (A. fabalis), and one white fronted goose (A. albifrons). These genomes share 83–95% nucleotide pairwise identities with previously identified GoCV genomes, most are recombinants with exchanged fragment sizes up to 50% of the genome. Higher diversity levels can be seen within the genomes from domestic geese compared with those from wild geese. In the GoCV capsid protein (cp) and replication associated protein (rep) gene sequences we found that episodic positive selection appears to largely mirror those of beak and feather disease virus and pigeon circovirus. Analysis of the secondary structure of the ssDNA genome revealed a conserved stem-loop structure with the G-C rich stem having a high degree of negative selection on these nucleotides.

DNA secondary structure formation by DNA shuffling of the conserved domains of the Cry protein of Bacillus thuringiensis

Background

The Cry toxins, or δ-endotoxins, are a diverse group of proteins produced by Bacillus thuringiensis. While DNA secondary structures are biologically relevant, it is unknown if such structures are formed in regions encoding conserved domains of Cry toxins under shuffling conditions. We analyzed 5 holotypes that encode Cry toxins and that grouped into 4 clusters according to their phylogenetic closeness. The mean number of DNA secondary structures that formed and the mean Gibbs free energy \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \left(\overline{\varDelta G}\right) $$\end{document}ΔG¯ were determined by an in silico analysis using different experimental DNA shuffling scenarios. In terms of spontaneity, shuffling efficiency was directly proportional to the formation of secondary structures but inversely proportional to ∆G.

Results

The results showed a shared thermodynamic pattern for each cluster and relationships among sequences that are phylogenetically close at the protein level. The regions of the cry11Aa, Ba and Bb genes that encode domain I showed more spontaneity and thus a greater tendency to form secondary structures (<∆G). In the region of domain III; this tendency was lower (>∆G) in the cry11Ba and Bb genes. Proteins that are phylogenetically closer to Cry11Ba and Cry11Bb, such as Cry2Aa and Cry18Aa, maintained the same thermodynamic pattern. More distant proteins, such as Cry1Aa, Cry1Ab, Cry30Aa and Cry30Ca, featured different thermodynamic patterns in their DNA.

Conclusion

These results suggest the presence of thermodynamic variations associated to the formation of secondary structures and an evolutionary relationship with regions that encode highly conserved domains in Cry proteins. The findings of this study may have a role in the in silico design of cry gene assembly by DNA shuffling techniques.

The blood DNA virome in 8,000 humans

The characterization of the blood virome is important for the safety of blood-derived transfusion products, and for the identification of emerging pathogens. We explored non-human sequence data from whole-genome sequencing of blood from 8,240 individuals, none of whom were ascertained for any infectious disease. Viral sequences were extracted from the pool of sequence reads that did not map to the human reference genome. Analyses sifted through close to 1 Petabyte of sequence data and performed 0.5 trillion similarity searches. With a lower bound for identification of 2 viral genomes/100,000 cells, we mapped sequences to 94 different viruses, including sequences from 19 human DNA viruses, proviruses and RNA viruses (herpesviruses, anelloviruses, papillomaviruses, three polyomaviruses, adenovirus, HIV, HTLV, hepatitis B, hepatitis C, parvovirus B19, and influenza virus) in 42% of the study participants. Of possible relevance to transfusion medicine, we identified Merkel cell polyomavirus in 49 individuals, papillomavirus in blood of 13 individuals, parvovirus B19 in 6 individuals, and the presence of herpesvirus 8 in 3 individuals. The presence of DNA sequences from two RNA viruses was unexpected: Hepatitis C virus is revealing of an integration event, while the influenza virus sequence resulted from immunization with a DNA vaccine. Age, sex and ancestry contributed significantly to the prevalence of infection. The remaining 75 viruses mostly reflect extensive contamination of commercial reagents and from the environment. These technical problems represent a major challenge for the identification of novel human pathogens. Increasing availability of human whole-genome sequences will contribute substantial amounts of data on the composition of the normal and pathogenic human blood virome. Distinguishing contaminants from real human viruses is challenging.

Novel sequencing technologies offer insight into the virome in human samples. Here, we identify the viral DNA sequences in blood of over 8,000 individuals undergoing whole genome sequencing. This approach serves to identify 94 viruses; however, many are shown to reflect widespread DNA contamination of commercial reagents or of environmental origin. While this represents a significant limitation to reliably identify novel viruses infecting humans, we could confidently detect sequences and quantify abundance of 19 human viruses in 42% of individuals. Ancestry, sex, and age were important determinants of viral prevalence. This large study calls attention on the challenge of interpreting next generation sequencing data for the identification of novel viruses. However, it serves to categorize the abundance of human DNA viruses using an unbiased technique.

Comparative analysis of common regions found in babuviruses and alphasatellite molecules

Viruses in the genus Babuvirus have multi-component ssDNA genomes and often associate with alphasatellite molecules containing two common motifs, a common-region stem-loop (CR-SL) involved in initiation of rolling-circle replication and a common-region major (CR-M) motif involved in secondary-strand synthesis. We compared known babuvirus genome components and alphasatellite CR-SL and CR-M sequences, defining five divergent CR-SL sequence classes. We identified iterated sequence elements in babuvirus genome components that have particularly conserved sequences and spatial arrangements between known babuviruses.

The Strange Lifestyle of Multipartite Viruses

Multipartite viruses have one of the most puzzling genetic organizations found in living organisms. These viruses have several genome segments, each containing only a part of the genetic information, and each individually encapsidated into a separate virus particle. While countless studies on molecular and cellular mechanisms of the infection cycle of multipartite viruses are available, just as for other virus types, very seldom is their lifestyle questioned at the viral system level. Moreover, the rare available “system” studies are purely theoretical, and their predictions on the putative benefit/cost balance of this peculiar genetic organization have not received experimental support. In light of ongoing progresses in general virology, we here challenge the current hypotheses explaining the evolutionary success of multipartite viruses and emphasize their shortcomings. We also discuss alternative ideas and research avenues to be explored in the future in order to solve the long-standing mystery of how viral systems composed of interdependent but physically separated information units can actually be functional.

Molecular characterization and prevalence of two capulaviruses: Alfalfa leaf curl virus from France and Euphorbia caput-medusae latent virus from South Africa

Little is known about the prevalence, diversity, evolutionary processes, genomic structures and population dynamics of viruses in the divergent geminivirus lineage known as the capulaviruses. We determined and analyzed full genome sequences of 13 Euphorbia caput-medusae latent virus (EcmLV) and 26 Alfalfa leaf curl virus (ALCV) isolates, and partial genome sequences of 23 EcmLV and 37 ALCV isolates. While EcmLV was asymptomatic in uncultivated southern African Euphorbia caput-medusae, severe alfalfa disease symptoms were associated with ALCV in southern France. The prevalence of both viruses exceeded 10% in their respective hosts. Besides using patterns of detectable negative selection to identify ORFs that are probably functionally expressed, we show that ALCV and EcmLV both display evidence of inter-species recombination and biologically functional genomic secondary structures. Finally, we show that whereas the EcmLV populations likely experience restricted geographical dispersion, ALCV is probably freely moving across the French Mediterranean region.

Impact of the MRN Complex on Adeno-Associated Virus Integration and Replication during Coinfection with Herpes Simplex Virus 1

Adeno-associated virus (AAV) is a helper-dependent parvovirus that requires coinfection with adenovirus (AdV) or herpes simplex virus 1 (HSV-1) to replicate. In the absence of the helper virus, AAV can persist in an episomal or integrated form. Previous studies have analyzed the DNA damage response (DDR) induced upon AAV replication to understand how it controls AAV replication. In particular, it was shown that the Mre11-Rad50-Nbs1 (MRN) complex, a major player of the DDR induced by double-stranded DNA breaks and stalled replication forks, could negatively regulate AdV and AAV replication during coinfection. In contrast, MRN favors HSV-1 replication and is recruited to AAV replication compartments that are induced in the presence of HSV-1. In this study, we examined the role of MRN during AAV replication induced by HSV-1. Our results indicated that knockdown of MRN significantly reduced AAV DNA replication after coinfection with wild-type (wt) HSV-1 or HSV-1 with the polymerase deleted. This effect was specific to wt AAV, since it did not occur with recombinant AAV vectors. Positive regulation of AAV replication by MRN was dependent on its DNA tethering activity but did not require its nuclease activities. Importantly, knockdown of MRN also negatively regulated AAV integration within the human AAVS1 site, both in the presence and in the absence of HSV-1. Altogether, this work identifies a new function of MRN during integration of the AAV genome and demonstrates that this DNA repair complex positively regulates AAV replication in the presence of HSV-1.

IMPORTANCE

Viral DNA genomes trigger a DNA damage response (DDR), which can be either detrimental or beneficial for virus replication. Adeno-associated virus (AAV) is a defective parvovirus that requires the help of an unrelated virus such as adenovirus (AdV) or herpes simplex virus 1 (HSV-1) for productive replication. Previous studies have demonstrated that the cellular Mre11-Rad50-Nbs1 (MRN) complex, a sensor and regulator of the DDR, negatively regulates AAV replication during coinfection with AdV, which counteracts this effect by inactivating the complex. Here, we demonstrate that MRN positively regulates AAV replication during coinfection with HSV-1. Importantly, our study also indicates that MRN also favors integration of AAV genomes within the human AAVS1 site. Altogether, this work indicates that MRN differentially regulates AAV replication depending on the helper virus which is present and identifies a new function of this DNA repair complex during AAV integration.

Functionally conserved architecture of hepatitis C virus RNA genomes

Hepatitis C virus (HCV) infects over 170 million people worldwide and is a leading cause of liver disease and cancer. The virus has a 9,650-nt, single-stranded, messenger-sense RNA genome that is infectious as an independent entity. The RNA genome has evolved in response to complex selection pressures, including the need to maintain structures that facilitate replication and to avoid clearance by cell-intrinsic immune processes. Here we used high-throughput, single-nucleotide resolution information to generate and functionally test data-driven structural models for three diverse HCV RNA genomes. We identified, de novo, multiple regions of conserved RNA structure, including all previously characterized cis-acting regulatory elements and also multiple novel structures required for optimal viral fitness. Well-defined RNA structures in the central regions of HCV genomes appear to facilitate persistent infection by masking the genome from RNase L and double-stranded RNA-induced innate immune sensors. This work shows how structure-first comparative analysis of entire genomes of a pathogenic RNA virus enables comprehensive and concise identification of regulatory elements and emphasizes the extensive interrelationships among RNA genome structure, viral biology, and innate immune responses.

SDT: a virus classification tool based on pairwise sequence alignment and identity calculation

The perpetually increasing rate at which viral full-genome sequences are being determined is creating a pressing demand for computational tools that will aid the objective classification of these genome sequences. Taxonomic classification approaches that are based on pairwise genetic identity measures are potentially highly automatable and are progressively gaining favour with the International Committee on Taxonomy of Viruses (ICTV). There are, however, various issues with the calculation of such measures that could potentially undermine the accuracy and consistency with which they can be applied to virus classification. Firstly, pairwise sequence identities computed based on multiple sequence alignments rather than on multiple independent pairwise alignments can lead to the deflation of identity scores with increasing dataset sizes. Also, when gap-characters need to be introduced during sequence alignments to account for insertions and deletions, methodological variations in the way that these characters are introduced and handled during pairwise genetic identity calculations can cause high degrees of inconsistency in the way that different methods classify the same sets of sequences. Here we present Sequence Demarcation Tool (SDT), a free user-friendly computer program that aims to provide a robust and highly reproducible means of objectively using pairwise genetic identity calculations to classify any set of nucleotide or amino acid sequences. SDT can produce publication quality pairwise identity plots and colour-coded distance matrices to further aid the classification of sequences according to ICTV approved taxonomic demarcation criteria. Besides a graphical interface version of the program for Windows computers, command-line versions of the program are available for a variety of different operating systems (including a parallel version for cluster computing platforms).

The influence of secondary structure, selection and recombination on rubella virus nucleotide substitution rate estimates

BACKGROUND

Annually, rubella virus (RV) still causes severe congenital defects in around 100 000 children globally. An attempt to eradicate RV is currently underway and analytical tools to monitor the global decline of the last remaining RV lineages will be useful for assessing the effectiveness of this endeavour. RV evolves rapidly enough that much of this information might be inferable from RV genomic sequence data.

METHODS

Using BEASTv1.8.0, we analysed publically available RV sequence data to estimate genome-wide and gene-specific nucleotide substitution rates to test whether current estimates of RV substitution rates are representative of the entire RV genome. We specifically accounted for possible confounders of nucleotide substitution rate estimates, such as temporally biased sampling, sporadic recombination, and natural selection favouring either increased or decreased genetic diversity (estimated by the PARRIS and FUBAR methods), at nucleotide sites within the genomic secondary structures (predicted by the NASP method).

RESULTS

We determine that RV nucleotide substitution rates range from 1.19 × 10(-3) substitutions/site/year in the E1 region to 7.52 × 10(-4) substitutions/site/year in the P150 region. We find that differences between substitution rate estimates in different RV genome regions are largely attributable to temporal sampling biases such that datasets containing higher proportions of recently sampled sequences, will tend to have inflated estimates of mean substitution rates. Although there exists little evidence of positive selection or natural genetic recombination in RV, we show that RV genomes possess pervasive biologically functional nucleic acid secondary structure and that purifying selection acting to maintain this structure contributes substantially to variations in estimated nucleotide substitution rates across RV genomes.

CONCLUSION

Both temporal sampling biases and purifying selection favouring the conservation of RV nucleic acid secondary structures have an appreciable impact on substitution rate estimates but do not preclude the use of RV sequence data to date ancestral sequences. The combination of uniformly high substitution rates across the RV genome and strong temporal structure within the available sequence data, suggests that such data should be suitable for tracking the demographic, epidemiological and movement dynamics of this virus during eradication attempts.

The RXL motif of the African cassava mosaic virus Rep protein is necessary for rereplication of yeast DNA and viral infection in plants

Geminiviruses, single-stranded DNA plant viruses, encode a replication-initiator protein (Rep) that is indispensable for virus replication. A potential cyclin interaction motif (RXL) in the sequence of African cassava mosaic virus Rep may be an alternative link to cell cycle controls to the known interaction with plant homologs of retinoblastoma protein (pRBR). Mutation of this motif abrogated rereplication in fission yeast induced by expression of wildtype Rep suggesting that Rep interacts via its RXL motif with one or several yeast proteins. The RXL motif is essential for viral infection of Nicotiana benthamiana plants, since mutation of this motif in infectious clones prevented any symptomatic infection. The cell-cycle link (Clink) protein of a nanovirus (faba bean necrotic yellows virus) was investigated that activates the cell cycle by binding via its LXCXE motif to pRBR. Expression of wildtype Clink and a Clink mutant deficient in pRBR-binding did not trigger rereplication in fission yeast.

Appearances can be deceptive: revealing a hidden viral infection with deep sequencing in a plant quarantine context

Comprehensive inventories of plant viral diversity are essential for effective quarantine and sanitation efforts. The safety of regulated plant material exchanges presently relies heavily on techniques such as PCR or nucleic acid hybridisation, which are only suited to the detection and characterisation of specific, well characterised pathogens. Here, we demonstrate the utility of sequence-independent next generation sequencing (NGS) of both virus-derived small interfering RNAs (siRNAs) and virion-associated nucleic acids (VANA) for the detailed identification and characterisation of viruses infecting two quarantined sugarcane plants. Both plants originated from Egypt and were known to be infected with Sugarcane streak Egypt Virus (SSEV; Genus Mastrevirus, Family Geminiviridae), but were revealed by the NGS approaches to also be infected by a second highly divergent mastrevirus, here named Sugarcane white streak Virus (SWSV). This novel virus had escaped detection by all routine quarantine detection assays and was found to also be present in sugarcane plants originating from Sudan. Complete SWSV genomes were cloned and sequenced from six plants and all were found to share >91% genome-wide identity. With the exception of two SWSV variants, which potentially express unusually large RepA proteins, the SWSV isolates display genome characteristics very typical to those of all other previously described mastreviruses. An analysis of virus-derived siRNAs for SWSV and SSEV showed them to be strongly influenced by secondary structures within both genomic single stranded DNA and mRNA transcripts. In addition, the distribution of siRNA size frequencies indicates that these mastreviruses are likely subject to both transcriptional and post-transcriptional gene silencing. Our study stresses the potential advantages of NGS-based virus metagenomic screening in a plant quarantine setting and indicates that such techniques could dramatically reduce the numbers of non-intercepted virus pathogens passing through plant quarantine stations.

Evolution of eukaryotic single-stranded DNA viruses of the Bidnaviridae family from genes of four other groups of widely different viruses

Single-stranded (ss)DNA viruses are extremely widespread, infect diverse hosts from all three domains of life and include important pathogens. Most ssDNA viruses possess small genomes that replicate by the rolling-circle-like mechanism initiated by a distinct virus-encoded endonuclease. However, viruses of the family Bidnaviridae, instead of the endonuclease, encode a protein-primed type B DNA polymerase (PolB) and hence break this pattern. We investigated the provenance of all bidnavirus genes and uncover an unexpected turbulent evolutionary history of these unique viruses. Our analysis strongly suggests that bidnaviruses evolved from a parvovirus ancestor from which they inherit a jelly-roll capsid protein and a superfamily 3 helicase. The radiation of bidnaviruses from parvoviruses was probably triggered by integration of the ancestral parvovirus genome into a large virus-derived DNA transposon of the Polinton (polintovirus) family resulting in the acquisition of the polintovirus PolB gene along with terminal inverted repeats. Bidnavirus genes for a receptor-binding protein and a potential novel antiviral defense modulator are derived from dsRNA viruses (Reoviridae) and dsDNA viruses (Baculoviridae), respectively. The unusual evolutionary history of bidnaviruses emphasizes the key role of horizontal gene transfer, sometimes between viruses with completely different genomes but occupying the same niche, in the emergence of new viral types.

Pigeon circoviruses display patterns of recombination, genomic secondary structure and selection similar to those of beak and feather disease viruses

Pigeon circovirus (PiCV) has a ~2 kb genome circular ssDNA genome. All but one of the known PiCV isolates have been found infecting pigeons in various parts of the world. In this study, we screened 324 swab and tissue samples from Polish pigeons and recovered 30 complete genomes, 16 of which came from birds displaying no obvious pathology. Together with 17 other publicly available PiCV complete genomes sampled throughout the Northern Hemisphere and Australia, we find that PiCV displays a similar degree of genetic diversity to that of the related psittacine-infecting circovirus species, beak and feather disease virus (BFDV). We show that, as is the case with its pathology and epidemiology, PiCV also displays patterns of recombination, genomic secondary structure and natural selection that are generally very similar to those of BFDV. It is likely that breeding facilities play a significant role in the emergence of new recombinant PiCV variants and given that ~50 % of the domestic pigeon population is infected subclinically, all pigeon breeding stocks should be screened routinely for this virus.

SDT: a virus classification tool based on pairwise sequence alignment and identity calculation

The perpetually increasing rate at which viral full-genome sequences are being determined is creating a pressing demand for computational tools that will aid the objective classification of these genome sequences. Taxonomic classification approaches that are based on pairwise genetic identity measures are potentially highly automatable and are progressively gaining favour with the International Committee on Taxonomy of Viruses (ICTV). There are, however, various issues with the calculation of such measures that could potentially undermine the accuracy and consistency with which they can be applied to virus classification. Firstly, pairwise sequence identities computed based on multiple sequence alignments rather than on multiple independent pairwise alignments can lead to the deflation of identity scores with increasing dataset sizes. Also, when gap-characters need to be introduced during sequence alignments to account for insertions and deletions, methodological variations in the way that these characters are introduced and handled during pairwise genetic identity calculations can cause high degrees of inconsistency in the way that different methods classify the same sets of sequences. Here we present Sequence Demarcation Tool (SDT), a free user-friendly computer program that aims to provide a robust and highly reproducible means of objectively using pairwise genetic identity calculations to classify any set of nucleotide or amino acid sequences. SDT can produce publication quality pairwise identity plots and colour-coded distance matrices to further aid the classification of sequences according to ICTV approved taxonomic demarcation criteria. Besides a graphical interface version of the program for Windows computers, command-line versions of the program are available for a variety of different operating systems (including a parallel version for cluster computing platforms).