Improvements to pairwise sequence comparison (PASC): a genome-based web tool for virus classification

Two Novel Bacteriophage Species Against Hybrid Intestinal Pathogenic Escherichia coli/Extraintestinal Pathogenic Escherichia coli Strains

Background

Colibacillosis caused by Escherichia coli is one of the main problems in the swine industry. In addition, the emergence of antimicrobial resistance and the combination of virulence genes among pathotypes have led to the emergence of more virulent pathogenic E. coli strains. Phage therapy has become a promising approach to address these issues.

Materials and Methods

Virulence genes for intestinal pathogenic E. coli (IPEC) and extraintestinal pathogenic E. coli (ExPEC) were investigated in pathogenic E. coli isolated from pigs. In addition, two potential phages, vB_EcoM-RPN187 and vB_EcoM-RPN226, isolated in our previous study, were further characterized in this study.

Results

Both phages were lytic and were highly effective at 20-37°C. Interestingly, they infected the hybrid IPEC/ExPEC strains. vB_EcoM-RPN187 and vB_EcoM-RPN226 possess 167 kbp of linear double-stranded DNA without virulence or antibiotic resistance genes and may be classified as new phage species in the genera Mosigvirus and Tequatrovirus, respectively.

Conclusion

Both phages could be promising candidates for phage therapy against pathogenic E. coli.

Diving deep into fish bornaviruses: Uncovering hidden diversity and transcriptional strategies through comprehensive data mining

Recently, we discovered two novel orthobornaviruses in colubrid and viperid snakes using an in silico data-mining approach. Here, we present the results of a screening of more than 100,000 nucleic acid sequence datasets of fish samples from the Sequence Read Archive (SRA) for potential bornaviral sequences. We discovered the potentially complete genomes of seven bornavirids in datasets from osteichthyans and chondrichthyans. Four of these are likely to represent novel species within the genus Cultervirus, and we propose that one genome represents a novel genus within the family of Bornaviridae. Specifically, we identified sequences of Wǔhàn sharpbelly bornavirus in sequence data from the widely used grass carp liver and kidney cell lines L8824 and CIK, respectively. A complete genome of Murray-Darling carp bornavirus was identified in sequence data from a goldfish (Carassius auratus). The newly discovered little skate bornavirus, identified in the little skate (Leucoraja erinacea) dataset, contained a novel and unusual genomic architecture (N-Vp1-Vp2-X-P-G-M-L), as compared to other bornavirids. Its genome is thought to encode two additional open reading frames (tentatively named Vp1 and Vp2), which appear to represent ancient duplications of the gene encoding the viral glycoprotein (G). The datasets also provided insights into the possible transcriptional gradients of these bornavirids and revealed previously unknown splicing mechanisms.

Virus taxonomy and the role of the International Committee on Taxonomy of Viruses (ICTV)

The taxonomy of viruses is developed and overseen by the International Committee on Taxonomy of Viruses (ICTV), which scrutinizes, approves and ratifies taxonomic proposals, and maintains a list of virus taxa with approved names (https://ictv.global). The ICTV has approximately 180 members who vote by simple majority. Taxon-specific Study Groups established by the ICTV have a combined membership of over 600 scientists from the wider virology community; they provide comprehensive expertise across the range of known viruses and are major contributors to the creation and evaluation of taxonomic proposals. Proposals can be submitted by anyone and will be considered by the ICTV irrespective of Study Group support. Thus, virus taxonomy is developed from within the virology community and realized by a democratic decision-making process. The ICTV upholds the distinction between a virus or replicating genetic element as a physical entity and the taxon category to which it is assigned. This is reflected by the nomenclature of the virus species taxon, which is now mandated by the ICTV to be in a binomial format (genus + species epithet) and is typographically distinct from the names of viruses. Classification of viruses below the rank of species (such as, genotypes or strains) is not within the remit of the ICTV. This article, authored by the ICTV Executive Committee, explains the principles of virus taxonomy and the organization, function, processes and resources of the ICTV, with the aim of encouraging greater understanding and interaction among the wider virology community.

Genomic Analysis of a New Freshwater Cyanophage Lbo240-yong1 Suggests a New Taxonomic Family of Bacteriophages

A worldwide ecological issue, cyanobacterial blooms in marine and freshwater have caused enormous losses in both the economy and the environment. Virulent cyanophages-specifically, infecting and lysing cyanobacteria-are key ecological factors involved in limiting the overall extent of the population development of cyanobacteria. Over the past three decades, reports have mainly focused on marine Prochlorococcus and Synechococcus cyanophages, while information on freshwater cyanophages remained largely unknown. In this study, a novel freshwater cyanophage, named Lbo240-yong1, was isolated via the double-layer agar plate method using Leptolyngbya boryana FACHB-240 as a host. Transmission electron microscopy observation illustrated the icosahedral head (50 ± 5 nm in diameter) and short tail (20 ± 5 nm in length) of Lbo240-yong1. Experimental infection against 37 cyanobacterial strains revealed that host-strain-specific Lbo240-yong1 could only lyse FACHB-240. The complete genome of Lbo240-yong1 is a double-stranded DNA of 39,740 bp with a G+C content of 51.99%, and it harbors 44 predicted open reading frames (ORFs). A Lbo240-yong1 ORF shared the highest identity with a gene of a filamentous cyanobacterium, hinting at a gene exchange between the cyanophage and cyanobacteria. A BLASTn search illustrated that Lbo240-yong1 had the highest sequence similarity with the Phormidium cyanophage Pf-WMP4 (89.67% identity, 84% query coverage). In the proteomic tree based on genome-wide sequence similarities, Lbo240-yong1, three Phormidium cyanophages (Pf-WMP4, Pf-WMP3, and PP), one Anabaena phage (A-4L), and one unclassified Arthronema cyanophage (Aa-TR020) formed a monophyletic group that was more deeply diverging than several other families. Pf-WMP4 is the only member of the independent genus Wumpquatrovirus that belongs to the Caudovircetes class. Pf-WMP3 and PP formed the independent genus Wumptrevirus. Anabaena phage A-4L is the only member of the independent Kozyakovvirus genus. The six cyanopodoviruses share similar gene arrangements. Eight core genes were found in them. We propose, here, to set up a new taxonomic family comprising the six freshwater cyanopodoviruses infecting filamentous cyanobacteria. This study enriched the field's knowledge of freshwater cyanophages.

Subspecific rodent taxa as the relevant host taxonomic level for mammarenavirus host specificity

Mastomys natalensis-borne mammarenaviruses appear specific to subspecific M. natalensis taxa rather than to the whole species. Yet mammarenaviruses carried by M. natalensis are known to spill over and jump hosts in northern sub-Saharan Africa. Phylogeographic studies increasingly show that, like M. natalensis, small mammals in sub-Saharan Africa are often genetically structured into several subspecific taxa. Other mammarenaviruses may thus also form virus-subspecific host taxon associations. To investigate this, and if mammarenaviruses carried by M. natalensis in southern Africa are less prone to spill-over, we screened 1225 non-M. natalensis samples from Tanzania where many small mammal taxa meet. We found mammarenavirus RNA in 6 samples. Genetic/genomic characterisation confirmed they were not spill-over from M. natalensis. We detected host jumps among rodent tribe members and an association between mammarenaviruses and subspecific taxa of Mus minutoides and Grammomys surdaster, indicating host genetic structure may be crucial to understand virus distribution and host specificity.

Four principles to establish a universal virus taxonomy

A universal taxonomy of viruses is essential for a comprehensive view of the virus world and for communicating the complicated evolutionary relationships among viruses. However, there are major differences in the conceptualisation and approaches to virus classification and nomenclature among virologists, clinicians, agronomists, and other interested parties. Here, we provide recommendations to guide the construction of a coherent and comprehensive virus taxonomy, based on expert scientific consensus. Firstly, assignments of viruses should be congruent with the best attainable reconstruction of their evolutionary histories, i.e., taxa should be monophyletic. This fundamental principle for classification of viruses is currently included in the International Committee on Taxonomy of Viruses (ICTV) code only for the rank of species. Secondly, phenotypic and ecological properties of viruses may inform, but not override, evolutionary relatedness in the placement of ranks. Thirdly, alternative classifications that consider phenotypic attributes, such as being vector-borne (e.g., "arboviruses"), infecting a certain type of host (e.g., "mycoviruses," "bacteriophages") or displaying specific pathogenicity (e.g., "human immunodeficiency viruses"), may serve important clinical and regulatory purposes but often create polyphyletic categories that do not reflect evolutionary relationships. Nevertheless, such classifications ought to be maintained if they serve the needs of specific communities or play a practical clinical or regulatory role. However, they should not be considered or called taxonomies. Finally, while an evolution-based framework enables viruses discovered by metagenomics to be incorporated into the ICTV taxonomy, there are essential requirements for quality control of the sequence data used for these assignments. Combined, these four principles will enable future development and expansion of virus taxonomy as the true evolutionary diversity of viruses becomes apparent.

Large-scale investigation of zoonotic viruses in the era of high-throughput sequencing

Zoonotic diseases considerably impact public health and socioeconomics. RNA viruses reportedly caused approximately 94% of zoonotic diseases documented from 1990 to 2010, emphasizing the importance of investigating RNA viruses in animals. Furthermore, it has been estimated that hundreds of thousands of animal viruses capable of infecting humans are yet to be discovered, warning against the inadequacy of our understanding of viral diversity. High-throughput sequencing (HTS) has enabled the identification of viral infections with relatively little bias. Viral searches using both symptomatic and asymptomatic animal samples by HTS have revealed hidden viral infections. This review introduces the history of viral searches using HTS, current analytical limitations, and future potentials. We primarily summarize recent research on large-scale investigations on viral infections reusing HTS data from public databases. Furthermore, considering the accumulation of uncultivated viruses, we discuss current studies and challenges for connecting viral sequences to their phenotypes using various approaches: performing data analysis, developing predictive modeling, or implementing high-throughput platforms of virological experiments. We believe that this article provides a future direction in large-scale investigations of potential zoonotic viruses using the HTS technology.

Novel Aeromonas Phage Ahy-Yong1 and Its Protective Effects against Aeromonas hydrophila in Brocade Carp (Cyprinus aka Koi)

Aeromonas hydrophila is a zoonotic pathogen and an important fish pathogen. A new lytic phage, Ahy-yong1, against multi-antibiotic-resistant pathogen A. hydrophila was isolated, identified, and tentatively used in therapy. Ahy-yong1 possesses a head of approximately 66 nm in diameter and a short tail of approximately 26 nm in length and 32 nm in width. Its complete dsDNA genome is 43,374 bp with a G + C content of 59.4%, containing 52 predicted opening reading frames (ORFs). Taxonomic analysis indicated Ahy-yong1 as a new species of the Ahphunavirus genus of the Autographiviridae family of the Caudoviricetes class. Ahy-yong1 was active only against its indicator host strain among the 35 strains tested. It is stable at 30-40 °C and at pH 2-12. Aeromonas phage Ahy-yong1 revealed an effective biofilm removal capacity and an obvious protective effect in brocade carp (Cyprinus aka Koi). The average cumulative mortality for the brocade carp in the blank groups intraperitoneally injected with PBS was 1.7% ± 2.4%;for the control groups treated with A. hydrophila (10⁸ CFU/fish) via intraperitoneal injection, it was 100.00%;and for the test group I, successively treated with A. hydrophila (10⁸ CFU/fish) and Aeromonas phage Ahy-yong1 (10⁷ PFU/fish) via intraperitoneal injection witha time interval of 2 hours, it was only 43.4% ± 4.7%. Furthermore, the cumulative mortality of the test group II, successively treated with Aeromonas phage Ahy-yong1 (10⁷ PFU/fish) and A. hydrophila (10⁸ CFU/fish), was only 20.0% ± 8.2%, and that of the test group III, simultaneously treated with Aeromonas phage Ahy-yong1 (10⁷ PFU/fish) and A. hydrophila (10⁸ CFU/fish), was only 30.0% ± 8.2%. The results demonstrated that phage Ahy-yong1 was very effective in the therapies against A. hydrophila A18, prophylaxis was more effective than rescue, and earlier treatment was better for the reduction of mortality. This study enriches knowledge about Aeromonas phages.

Boid Inclusion Body Disease Is Also a Disease of Wild Boa Constrictors

Reptarenaviruses cause boid inclusion body disease (BIBD), a potentially fatal disease, occurring in captive constrictor snakes boas and pythons worldwide. Classical BIBD, characterized by the formation of pathognomonic cytoplasmic inclusion bodies (IBs), occurs mainly in boas, whereas in pythons, for example, reptarenavirus infection most often manifests as central nervous system signs with limited IB formation. The natural hosts of reptarenaviruses are unknown, although free-ranging/wild constrictor snakes are among the suspects. Here, we report BIBD with reptarenavirus infection in indigenous captive and wild boid snakes in Costa Rica using histology, immunohistology, transmission electron microscopy, and next-generation sequencing (NGS). The snakes studied represented diagnostic postmortem cases of captive and wild-caught snakes since 1989. The results from NGS on archival paraffin blocks confirm that reptarenaviruses were already present in wild boa constrictors in Costa Rica in the 1980s. Continuous sequences that were de novo assembled from the low-quality RNA obtained from paraffin-embedded tissue allowed the identification of a distinct pair of reptarenavirus S and L segments in all studied animals; in most cases, reference assembly could recover almost complete segments. Sampling of three prospective cases in 2018 allowed an examination of fresh blood or tissues and resulted in the identification of additional reptarenavirus segments and hartmanivirus coinfection. Our results show that BIBD is not only a disease of captive snakes but also occurs in indigenous wild constrictor snakes in Costa Rica, suggesting boa constrictors to play a role in natural reptarenavirus circulation. IMPORTANCE The literature describes cases of boid inclusion body disease (BIBD) in captive snakes since the 1970s, and in the 2010s, others and ourselves identified reptarenaviruses as the causative agent. BIBD affects captive snakes globally, but the origin and the natural host of reptarenaviruses remain unknown. In this report, we show BIBD and reptarenavirus infections in two native Costa Rican constrictor snake species, and by studying archival samples, we show that both the viruses and the disease have been present in free-ranging/wild snakes in Costa Rica at least since the 1980s. The diagnosis of BIBD in wild boa constrictors suggests that this species plays a role in the circulation of reptarenaviruses. Additional sample collection and analysis would help to clarify this role further and the possibility of, e.g., vector transmission from an arthropod host.

A lipopolysaccharide-dependent phage infects a pseudomonad phytopathogen and can evolve to evade phage resistance

Bacterial pathogens are major causes of crop diseases, leading to significant production losses. For instance, kiwifruit canker, caused by the phytopathogen Pseudomonas syringae pv. actinidiae (Psa), has posed a global challenge to kiwifruit production. Treatment with copper and antibiotics, whilst initially effective, is leading to the rise of bacterial resistance, requiring new biocontrol approaches. Previously, we isolated a group of closely related Psa phages with biocontrol potential, which represent environmentally sustainable antimicrobials. However, their deployment as antimicrobials requires further insight into their properties and infection strategy. Here, we provide an in-depth examination of the genome of ΦPsa374-like phages and show that they use lipopolysaccharides (LPS) as their main receptor. Through proteomics and cryo-electron microscopy of ΦPsa374, we revealed the structural proteome and that this phage possess a T = 9 capsid triangulation, unusual for myoviruses. Furthermore, we show that ΦPsa374 phage resistance arises in planta through mutations in a glycosyltransferase involved in LPS synthesis. Lastly, through in vitro evolution experiments we showed that phage resistance is overcome by mutations in a tail fibre and structural protein of unknown function in ΦPsa374. This study provides new insight into the properties of ΦPsa374-like phages that informs their use as antimicrobials against Psa.

A Novel Freshwater Cyanophage Mae-Yong1326-1 Infecting Bloom-Forming Cyanobacterium Microcystis aeruginosa

Microcystis aeruginosa is a major harmful cyanobacterium causing water bloom worldwide. Cyanophage has been proposed as a promising tool for cyanobacterial bloom. In this study, M. aeruginosa FACHB-1326 was used as an indicator host to isolate cyanophage from Lake Taihu. The isolated Microcystis cyanophage Mae-Yong1326-1 has an elliptical head of about 47 nm in diameter and a slender flexible tail of about 340 nm in length. Mae-Yong1326-1 could lyse cyanobacterial strains across three orders (Chroococcales, Nostocales, and Oscillatoriales) in the host range experiments. Mae-Yong1326-1 was stable in stability tests, maintaining high titers at 0–40 °C and at a wide pH range of 3–12. Mae-Yong 1326-1 has a burst size of 329 PFU/cell, which is much larger than the reported Microcystis cyanophages so far. The complete genome of Mae-Yong1326-1 is a double-stranded DNA of 48, 822 bp, with a G + C content of 71.80% and long direct terminal repeats (DTR) of 366 bp, containing 57 predicted ORFs. No Mae-Yong1326-1 ORF was found to be associated with virulence factor or antibiotic resistance. PASC scanning illustrated that the highest nucleotide sequence similarity between Mae-Yong1326-1 and all known phages in databases was only 17.75%, less than 70% (the threshold to define a genus), which indicates that Mae-Yong1326-1 belongs to an unknown new genus. In the proteomic tree based on genome-wide sequence similarities, Mae-Yong1326-1 distantly clusters with three unclassified Microcystis cyanophages (MinS1, Mwe-Yong1112-1, and Mwes-Yong2). These four Microcystis cyanophages form a monophyletic clade, which separates at a node from the other clade formed by two independent families (Zierdtviridae and Orlajensenviridae) of Caudoviricetes class. We propose to establish a new family to harbor the Microcystis cyanophages Mae-Yong1326-1, MinS1, Mwe-Yong1112-1, and Mwes-Yong2. This study enriched the understanding of freshwater cyanophages.

Lack of robust evidence for a Wolbachia infection in Anopheles gambiae from Burkina Faso

The endosymbiont Wolbachia can have major effects on the reproductive fitness, and vectorial capacity of host insects and may provide new avenues to control mosquito-borne pathogens. Anopheles gambiae s.l is the major vector of malaria in Africa but the use of Wolbachia in this species has been limited by challenges in establishing stable transinfected lines and uncertainty around native infections. High frequencies of infection of Wolbachia have been previously reported in An. gambiae collected from the Valle du Kou region of Burkina Faso in 2011 and 2014. Here, we re-evaluated the occurrence of Wolbachia in natural samples, collected from Valle du Kou over a 12-year time span, and in addition, expanded sampling to other sites in Burkina Faso. Our results showed that, in contrast to earlier reports, Wolbachia is present at an extremely low prevalence in natural population of An. gambiae. From 5341 samples analysed, only 29 were positive for Wolbachia by nested PCR representing 0.54% of prevalence. No positive samples were found with regular PCR. Phylogenetic analysis of 16S rRNA gene amplicons clustered across supergroup B, with some having similarity to sequences previously found in Anopheles from Burkina Faso. However, we cannot discount the possibility that the amplicon positive samples we detected were due to environmental contamination or were false positives. Regardless, the lack of a prominent native infection in An. gambiae s.l. is encouraging for applications utilizing Wolbachia transinfected mosquitoes for malaria control.

Co-Infection of Mammarenaviruses in a Wild Mouse, Tanzania

Mammarenaviruses are bi-segmented RNA viruses. They encompass viruses responsible for several severe diseases in humans. While performing a de novo assembly of a new virus found in a wild single-striped grass mouse in Tanzania, we found a single S but two divergent L segments. Natural co-infections, common within reptarenaviruses in captivity, were never reported for mammarenaviruses and never in a wild sample. This finding can have implications for virus evolution as co-infection could trigger viral recombination/reassortment in natural reservoirs.

Genome sequence of the novel freshwater Microcystis cyanophage Mwe-Yong1112-1

The freshwater cyanophage Mwe-Yong1112-1 was isolated using Microcystis wesenbergii as a host and found to have an icosahedral head, about 45 nm in diameter, and a flexible tail, approximately 133 nm in length and 4.5 nm in width. The complete genome of the cyanophage is 39,679 bp in length with a G+C content of 66.6%. Mwe-Yong1112-1 shared the highest pairwise average nucleotide identity (ANI) value of 67.7% (below the ≥95% boundary to define a species) and the highest nucleotide sequence similarity of 17.48% (below the >70% boundary to define a genus) with the most closely related phage. In a proteomic tree, Mwe-Yong1112-1 and three unclassified phages formed a monophyletic clade between the families Saparoviridae and Pyrstoviridae, but Mwe-Yong1112-1 occupied a separate branch from the other three phages, suggesting that it represents a new evolutionary lineage. This study enriches the available information about freshwater cyanophages.

Classification of new morbillivirus and jeilongvirus sequences from bats sampled in Brazil and Malaysia

As part of a broad One Health surveillance effort to detect novel viruses in wildlife and people, we report several paramyxovirus sequences sampled primarily from bats during 2013 and 2014 in Brazil and Malaysia, including seven from which we recovered full-length genomes. Of these, six represent the first full-length paramyxovirid genomes sequenced from the Americas, including two that are the first full-length bat morbillivirus genome sequences published to date. Our findings add to the vast number of viral sequences in public repositories, which have been increasing considerably in recent years due to the rising accessibility of metagenomics. Taxonomic classification of these sequences in the absence of phenotypic data has been a significant challenge, particularly in the subfamily Orthoparamyxovirinae, where the rate of discovery of novel sequences has been substantial. Using pairwise amino acid sequence classification (PAASC), we propose that five of these sequences belong to members of the genus Jeilongvirus and two belong to members of the genus Morbillivirus. We also highlight inconsistencies in the classification of Tupaia virus and Mòjiāng virus using the same demarcation criteria and suggest reclassification of these viruses into new genera. Importantly, this study underscores the critical importance of sequence length in PAASC analysis as well as the importance of biological characteristics such as genome organization in the taxonomic classification of viral sequences.

Genomic analysis of a novel active prophage of Hafnia paralvei

Little is known about the prophages in Hafniaceae bacteria. A novel Hafnia phage, yong2, was induced from Hafnia paralvei by treatment with mitomycin C. The phage has an elliptical head with dimensions of approximately 45 × 38 nm and a long noncontractile tail of approximately 157 × 4 nm. The complete genome of Hafnia phage yong2 is a 39,546-bp double-stranded DNA with a G+C content of 49.9%, containing 59 open reading frames (ORFs) and having at least one fixed terminus (GGGGCAGCGACA). In phylogenetic analysis, Hafnia phage yong2 clustered with four predicted Hafnia prophages and one predicted Enterobacteriaceae prophage. These prophages and members of the family Drexlerviridae together formed two distinct subclades nested within a clade, suggesting the existence of a novel class of prophages with conserved sequences and a unique evolutionary status not yet studied before in Hafniaceae and Enterobacteriaceae bacteria.

A World of Viruses Nested within Parasites: Unraveling Viral Diversity within Parasitic Flatworms (Platyhelminthes)

Because parasites have an inextricable relationship with their host, they have the potential to serve as viral reservoirs or facilitate virus host shifts. And yet, little is known about viruses infecting parasitic hosts except for blood-feeding arthropods that are well-known vectors of zoonotic viruses. Herein, we uncovered viruses of flatworms (phylum Platyhelminthes, group Neodermata) that specialize in parasitizing vertebrates and their ancestral free-living relatives. We discovered 115 novel viral sequences, including 1 in Macrostomorpha, 5 in Polycladida, 44 in Tricladida, 1 in Monogenea, 15 in Cestoda, and 49 in Trematoda, through data mining. The majority of newly identified viruses constitute novel families or genera. Phylogenetic analyses show that the virome of flatworms changed dramatically during the transition of neodermatans to a parasitic lifestyle. Most Neodermata viruses seem to codiversify with their host, with the exception of rhabdoviruses, which may switch hosts more often, based on phylogenetic relationships. Neodermata rhabdoviruses also have a position ancestral to vertebrate-associated rhabdo viruses, including lyssaviruses, suggesting that vertebrate-associated rhabdoviruses emerged from a flatworm rhabdovirus in a parasitized host. This study reveals an extensive diversity of viruses in Platyhelminthes and highlights the need to evaluate the role of viral infection in flatworm-associated diseases.

IMPORTANCE Little is known about the diversity of parasite-associated viruses and how these viruses may impact parasite fitness, parasite-host interactions, and virus evolution. The discovery of over a hundred viruses associated with a range of free-living and parasitic flatworms, including parasites of economic and clinical relevance, allowed us to compare the viromes of flatworms with contrasting lifestyles. The results suggest that flatworms acquired novel viruses after their transition to a parasitic lifestyle and highlight the possibility that they acquired viruses from their hosts and vice versa. An interesting example is the discovery of flatworm rhabdoviruses that have a position ancestral to rabies viruses and other vertebrate-associated rhabdoviruses, demonstrating that flatworm-associated viruses have emerged in a vertebrate host at least once in history. Therefore, parasitic flatworms may play a role in virus diversity and emergence. The roles that parasite-infecting viruses play in parasite-associated diseases remain to be investigated.

Jingchuvirales: a New Taxonomical Framework for a Rapidly Expanding Order of Unusual Monjiviricete Viruses Broadly Distributed among Arthropod Subphyla

Technical advances in metagenomics and metatranscriptomics have dramatically accelerated virus discovery in recent years. "Chuviruses" were first described in 2015 as obscure negative-sense RNA viruses of diverse arthropods. Although "chuviruses" first appeared to be members of the negarnaviricot order Mononegavirales in phylogenetic analyses using RNA-directed RNA polymerase sequences, further characterization revealed unusual gene orders in genomes that are nonsegmented, segmented, and/or possibly circular. Consequently, a separate order, Jingchuvirales, was established to include a monospecific family, Chuviridae. Recently, it has become apparent that jingchuvirals are broadly distributed and are therefore likely of ecological and economic importance. Here, we describe recent and ongoing efforts to create the necessary taxonomic framework to accommodate the expected flood of novel viruses belonging to the order.

First Characterization of a Hafnia Phage Reveals Extraordinarily Large Burst Size and Unusual Plaque Polymorphism

A unique lytic phage infecting Hafnia paralvei was isolated and identified. Hafnia phage Ca belongs to the family Autographiviridae, possessing an icosahedral head with a diameter of 55 nm and a short non-contractile tail. Unusually, the burst size of Hafnia phage Ca of 10,292 ± 1,097 plaque-forming units (PFUs)/cell is much larger than other dsDNA phages reported before. Compared to the genome of the related phage, Hafnia phage Ca genome contains extra genes including DNA mimic ocr, dGTP triphosphohydrolase inhibitor, endonuclease, endonuclease VII, and HNH homing endonuclease gene. Extraordinarily, the phage developed different sizes of plaques when a single plaque was picked out and inoculated on a double-layer Luria broth agar plate with its host. Furthermore, varied packaging tightness for the tails of Hafnia phage Ca was observed (tail length: 4.35–45.92 nm). Most of the tails appeared to be like a cone with appendages, some were dot-like, bun-like, table tennis racket handle-like, and ponytail-like. Although the complete genome of Hafnia phage Ca is 40,286 bp, an incomplete genome with a deletion of a 397-bp fragment, containing one ORF predicted as HNH homing endonuclease gene (HEG), was also found by high throughput sequencing. Most of the genome of the virus particles in large plaques is complete (>98%), while most of the genome of the virus particles in small plaques is incomplete (>98%), and the abundance of both of them in medium-sized plaques is similar (complete, 40%; incomplete, 60%). In an experiment to see if the phage could be protective to brocade carps intramuscularly injected with H. paralvei LY-23 and phage Ca, the protection rate of Hafnia phage Ca to brocade carp (Cyprinus aka Koi) against H. paralvei was 33.38% (0.01 < p < 0.05). This study highlights some new insights into the peculiar biological and genomic characteristics of phage.

Isolation, Characterization and Draft Genome Analysis of Bacteriophages Infecting Acidovorax citrulli

Bacterial fruit blotch and seedling blight, caused by Acidovorax citrulli, is one of the most destructive diseases of melon and watermelon in many countries. Pathogen-free seed and cultural practices are major pillars of the disease control. However, use of bacteriophages as natural biocontrol agents might also contribute to the disease management. Therefore, we isolated 12 bacteriophages specific to A. citrulli, from phyllosphere and rhizosphere of diseased watermelon plants. The phage strains were characterized based on their host range, plaque and virion morphology, thermal inactivation point, adsorption rate, one step growth curve, restriction fragment length polymorphism (RFLP), and genomic analysis. Transmission electron microscopy of three phage strains indicated that they belong to the order Caudovirales, family Siphoviridae. All phages lysed 30 out of 32 tested A. citrulli strains isolated in Serbia, and did not lyse other less related bacterial species. They produced clear plaques, 2 mm in diameter, on bacterial lawns of different A. citrulli strains after 24 h of incubation. The thermal inactivation point was 66 or 67°C. They were stable at pH 5–9, but were sensitive to chloroform and inactivated in either 5 or 10 min exposure to ultraviolet (UV) light. RFLP analysis using EcoRI, BsmI and BamHI enzymes did not show genetic differences among the tested phages. Adsorption rate and one step growth curve were determined for the Acidovorax phage ACF1. Draft genome sequence of the ACF1 phage was 59.377 bp in size, with guanine-cytosine (GC) content 64.5%, including 89 open reading frames. This phage shared a very high genomic identity with Acidovorax phage ACPWH, isolated in South Korea. Evaluation of systemic nature of ACF1 strain showed that it can be absorbed by roots and translocated to upper parts of watermelon plants where it survived up to 10 days.

VirStrain: a strain identification tool for RNA viruses

Viruses change constantly during replication, leading to high intra-species diversity. Although many changes are neutral or deleterious, some can confer on the virus different biological properties such as better adaptability. In addition, viral genotypes often have associated metadata, such as host residence, which can help with inferring viral transmission during pandemics. Thus, subspecies analysis can provide important insights into virus characterization. Here, we present VirStrain, a tool taking short reads as input with viral strain composition as output. We rigorously test VirStrain on multiple simulated and real virus sequencing datasets. VirStrain outperforms the state-of-the-art tools in both sensitivity and accuracy.

A Novel Freshwater Cyanophage, Mae-Yong924-1, Reveals a New Family

Cyanobacterial blooms are a worldwide ecological issue. Cyanophages are aquatic viruses specifically infecting cyanobacteria. Little is known about freshwater cyanophages. In this study, a freshwater cyanophage, Mae-Yong924-1, was isolated by the double-layer agar plate method using Microcystis aeruginosa FACHB-924 as an indicator host. Mae-Yong924-1 has several unusual characteristics: a unique shape, cross-taxonomic order infectivity and a very unique genome sequence. Mae-Yong924-1 contains a nearly spherical head of about 100 nm in diameter. The tail or tail-like structure (approximately 40 nm in length) is like the tassel of a round Chinese lantern. It could lyse six diverse cyanobacteria strains across three orders including Chroococcales, Nostocales and Oscillatoriales. The genome of the cyanophage is 40,325 bp in length, with a G + C content of 48.32%, and 59 predicted open reading frames (ORFs), only 12 (20%) of which were functionally annotated. Both BLASTn and BLASTx scanning resulted in “No significant similarity found”, i.e., the Mae-Yong924-1 genome shared extremely low homology with sequences in NCBI databases. Mae-Yong924-1 formed a root node alone and monopolized a root branch in the proteomic tree based on genome-wide sequence similarities. The results suggest that Mae-Yong924-1 may reveal a new unknown family apparently distinct from other viruses.

Type I CRISPR-Cas provides robust immunity but incomplete attenuation of phage-induced cellular stress

During infection, phages manipulate bacteria to redirect metabolism towards viral proliferation. To counteract phages, some bacteria employ CRISPR-Cas systems that provide adaptive immunity. While CRISPR-Cas mechanisms have been studied extensively, their effects on both the phage and the host during phage infection remains poorly understood. Here, we analysed the infection of Serratia by a siphovirus (JS26) and the transcriptomic response with, or without type I-E or I-F CRISPR-Cas immunity. In non-immune Serratia, phage infection altered bacterial metabolism by upregulating anaerobic respiration and amino acid biosynthesis genes, while flagella production was suppressed. Furthermore, phage proliferation required a late-expressed viral Cas4 homologue, which did not influence CRISPR adaptation. While type I-E and I-F immunity provided robust defence against phage infection, phage development still impacted the bacterial host. Moreover, DNA repair and SOS response pathways were upregulated during type I immunity. We also discovered that the type I-F system is controlled by a positive autoregulatory feedback loop that is activated upon phage targeting during type I-F immunity, leading to a controlled anti-phage response. Overall, our results provide new insight into phage-host dynamics and the impact of CRISPR immunity within the infected cell.

Bioinformatics of virus taxonomy: foundations and tools for developing sequence-based hierarchical classification

The genome sequence is the only characteristic readily obtainable for all known viruses, underlying the growing role of comparative genomics in organizing knowledge about viruses in a systematic evolution-aware way, known as virus taxonomy. Overseen by the International Committee on Taxonomy of Viruses (ICTV), development of virus taxonomy involves taxa demarcation at 15 ranks of a hierarchical classification, often in host-specific manner. Outside the ICTV remit, researchers assess fitting numerous unclassified viruses into the established taxa. They employ different metrics of virus clustering, basing on conserved domain(s), separation of viruses in rooted phylogenetic trees and pair-wise distance space. Computational approaches differ further in respect to methodology, number of ranks considered, sensitivity to uneven virus sampling, and visualization of results. Advancing and using computational tools will be critical for improving taxa demarcation across the virosphere and resolving rank origins in research that may also inform experimental virology.

A k-mer based approach for classifying viruses without taxonomy identifies viral associations in human autism and plant microbiomes

Viruses are an underrepresented taxa in the study and identification of microbiome constituents; however, they play an essential role in health, microbiome regulation, and transfer of genetic material. Only a few thousand viruses have been isolated, sequenced, and assigned a taxonomy, which limits the ability to identify and quantify viruses in the microbiome. Additionally, the vast diversity of viruses represents a challenge for classification, not only in constructing a viral taxonomy, but also in identifying similarities between a virus' genotype and its phenotype. However, the diversity of viral sequences can be leveraged to classify their sequences in metagenomic and metatranscriptomic samples, even if they do not have a taxonomy. To identify and quantify viruses in transcriptomic and genomic samples, we developed a dynamic programming algorithm for creating a classification tree out of 715,672 metagenome viruses. To create the classification tree, we clustered proportional similarity scores generated from the k-mer profiles of each of the metagenome viruses to create a database of metagenomic viruses. The resulting Kraken2 database of the metagenomic viruses can be found here: https://www.osti.gov/biblio/1615774 and is compatible with Kraken2. We then integrated the viral classification database with databases created with genomes from NCBI for use with ParaKraken (a parallelized version of Kraken provided in Supplemental Zip 1), a metagenomic/transcriptomic classifier. To illustrate the breadth of our utility for classifying metagenome viruses, we analyzed data from a plant metagenome study identifying genotypic and compartment specific differences between two Populus genotypes in three different compartments. We also identified a significant increase in abundance of eight viral sequences in post mortem brains in a human metatranscriptome study comparing Autism Spectrum Disorder patients and controls. We also show the potential accuracy for classifying viruses by utilizing both the JGI and NCBI viral databases to identify the uniqueness of viral sequences. Finally, we validate the accuracy of viral classification with NCBI databases containing viruses with taxonomy to identify pathogenic viruses in known COVID-19 and cassava brown streak virus infection samples. Our method represents the compulsory first step in better understanding the role of viruses in the microbiome by allowing for a more complete identification of sequences without taxonomy. Better classification of viruses will improve identifying associations between viruses and their hosts as well as viruses and other microbiome members. Despite the lack of taxonomy, this database of metagenomic viruses can be used with any tool that utilizes a taxonomy, such as Kraken, for accurate classification of viruses.

Epidemiological assessment of cassava mosaic disease in Burkina Faso

Surveys were conducted in 2016 and 2017 across the main cassava-growing regions of Burkina Faso to assess the status of cassava mosaic disease (CMD) and to determine the virus strains causing the disease, using field observation and phylogenetic analysis. CMD incidence varied between regions and across years but was lowest in Hauts-Bassins (6.0%, 2016 and 5.4%, 2017) and highest in Centre-Sud (18.5%, 2016) and in Boucle du Mouhoun (51.7%, 2017). The lowest CMD severity was found in Est region (2.0) for both years and the highest in Sud-Ouest region (3.3, 2016) and Centre-Sud region (2.8, 2017). The CMD infection was primarily associated with contaminated cuttings in all regions except in Hauts-Bassins, where whitefly-borne infection was higher than cuttings-borne infection in 2016. PCR screening of 687 samples coupled with sequence analysis revealed the presence of African cassava mosaic-like (ACMV-like) viruses and East African cassava mosaic-like (EACMV-like) viruses as single infections at 79.5% and 1.1%, respectively. Co-infections of ACMV-like and EACMV-like viruses were detected in 19.4% of the tested samples. In addition, 86.7% of the samples positive for EACMV-like virus were found to be positive for East African cassava mosaic Cameroon virus (EACMCMV). Phylogenetic analysis revealed the segregation of cassava mosaic geminiviruses (CMGs) from Burkina Faso into three clades specific to ACMV, African cassava mosaic Burkina Faso virus (ACMBFV), and EACMCMV, confirming the presence of these viruses. The results of this study show that EACMCMV occurrence may be more prevalent in Burkina Faso than previously thought.

A novel cripavirus of an ectoparasitoid wasp increases pupal duration and fecundity of the wasp's Drosophila melanogaster host

We identified a 9332-nucleotide-long novel picornaviral genome sequence in the transcriptome of an agriculturally important parasitoid wasp (Pachycrepoideus vindemmiae (Rondani, 1875)). The genome of the novel virus, Rondani's wasp virus 1 (RoWV-1), contains two long open reading frames encoding a nonstructural and a structural protein, respectively, and is 3'-polyadenylated. Phylogenetic analyses firmly place RoWV-1 into the dicistrovirid genus Cripavirus. We detected RoWV-1 in various tissues and life stages of the parasitoid wasp, with the highest virus load measured in the larval digestive tract. We demonstrate that RoWV-1 is transmitted horizontally from infected to uninfected wasps but not vertically to wasp offspring. Comparison of several important biological parameters between the infected and uninfected wasps indicates that RoWV-1 does not have obvious detrimental effects on wasps. We further demonstrate that RoWV-1 also infects Drosophila melanogaster (Meigen, 1830), the hosts of the pupal ectoparasitoid wasps, and thereby increases its pupal developmental duration and fecundity, but decreases the eclosion rate. Together, these results suggest that RoWV-1 may have a potential benefit to the wasp by increasing not only the number of potential wasp hosts but also the developmental time of the hosts to ensure proper development of wasp offspring.

Phylogenomic Characterization of Lopma Virus and Praja Virus, Two Novel Rodent-Borne Arteriviruses

Recent years have witnessed the discovery of several new viruses belonging to the family Arteriviridae, expanding the known diversity and host range of this group of complex RNA viruses. Although the pathological relevance of these new viruses is not always clear, several well-studied members of the family Arteriviridae are known to be important animal pathogens. Here, we report the complete genome sequences of four new arterivirus variants, belonging to two putative novel species. These new arteriviruses were discovered in African rodents and were given the names Lopma virus and Praja virus. Their genomes follow the characteristic genome organization of all known arteriviruses, even though they are only distantly related to currently known rodent-borne arteriviruses. Phylogenetic analysis shows that Lopma virus clusters in the subfamily Variarterivirinae, while Praja virus clusters near members of the subfamily Heroarterivirinae: the yet undescribed forest pouched giant rat arterivirus and hedgehog arterivirus 1. A co-divergence analysis of rodent-borne arteriviruses confirms that they share similar phylogenetic patterns with their hosts, with only very few cases of host shifting events throughout their evolutionary history. Overall, the genomes described here and their unique clustering with other arteriviruses further illustrate the existence of multiple rodent-borne arterivirus lineages, expanding our knowledge of the evolutionary origin of these viruses.

Diverse Bacteriophages Infecting the Bacterial Striped Catfish Pathogen Edwardsiella ictaluri

Bacteriophages infecting Edwardsiella ictaluri have been less investigated, although the host bacterium is one of the most important fish pathogens causing enteric septicemia of catfish (ESC). We present here two distinctly novel bacteriophages vB_EiM_PVN06 and vB_EiA_PVN09 infecting Edwardsiella ictaluri E1, with their geographical origins from the Mekong Delta, Vietnam. Bacteriophage vB_EiM_PVN06 native to a mud sample reveals complete differences of biological properties with the phage vB_EiA_PVN09 originated from a viscus of a healthy catfish (Pangasianodon hypophthalmus) cultured in the same area. Morphological analyses combined with genomic data indicate that phage vB_EiM_PVN06 is classified to Myoviridae family and shares high similarity with E. ictaluri phage PEi21 genome, while vB_EiA_PVN09 is a member of Teseptimavirus genus, Autographiviridae family, and mostly closes to phage vB_EcoP_IME390. The vB_EiA_PVN09 is a T7-like bacteriophage, which has been firstly found infecting to E. ictaluri, and host range analysis also evidences for the cross-infection of this phage to Escherichia coli K12 and Escherichia coli DH5α. Together, our research highlights the diversity of bacteriophages infecting the pathogen E. ictaluri and suggests further explorations of lytic phages in environmental niches, to be exploited in feasible strategies of phage therapy in ESC disease control.

Identification of a Novel Quinvirus in the Family Betaflexiviridae That Infects Winter Wheat

Yellow mosaic disease in winter wheat is usually attributed to the infection by bymoviruses or furoviruses; however, there is still limited information on whether other viral agents are also associated with this disease. To investigate the wheat viromes associated with yellow mosaic disease, we carried out de novo RNA sequencing (RNA-seq) analyses of symptomatic and asymptomatic wheat-leaf samples obtained from a field in Hokkaido, Japan, in 2018 and 2019. The analyses revealed the infection by a novel betaflexivirus, which tentatively named wheat virus Q (WVQ), together with wheat yellow mosaic virus (WYMV, a bymovirus) and northern cereal mosaic virus (a cytorhabdovirus). Basic local alignment search tool (BLAST) analyses showed that the WVQ strains (of which there are at least three) were related to the members of the genus Foveavirus in the subfamily Quinvirinae (family Betaflexiviridae). In the phylogenetic tree, they form a clade distant from that of the foveaviruses, suggesting that WVQ is a member of a novel genus in the Quinvirinae. Laboratory tests confirmed that WVQ, like WYMV, is potentially transmitted through the soil to wheat plants. WVQ was also found to infect rye plants grown in the same field. Moreover, WVQ-derived small interfering RNAs accumulated in the infected wheat plants, indicating that WVQ infection induces antiviral RNA silencing responses. Given its common coexistence with WYMV, the impact of WVQ infection on yellow mosaic disease in the field warrants detailed investigation.

Taxonomy classification using genomic footprint of mitochondrial sequences

BACKGROUND

Advancement in the sequencing technology yields a huge number of genomes of a multitude of organisms in our planet. One of the fundamental tasks for processing and analyzing these sequences is to organize them in the existing taxonomic orders. <P> Method: Recently we proposed a novel approach, GenFooT, of taxonomy classification using the concept of genomic footprint (GFP). The technique is further refined and enhanced in this work leading to improved accuracies in the task of taxonomic classification on various benchmark datasets. GenFooT maps a genome sequence in a 2D coordinate space and extracts features from that representation. It uses two hyper-parameters, namely block size and number of fragments of genomic sequence while computing the feature. In this work, we propose an analysis for choosing values of those parameters adaptively from the sequences. The enhanced version of GenFooT is named GenFooT2. <P> Results and Conclusion: We have experimented GenFooT2 on ten different biological datasets of genomic sequences of various organisms belonging to different taxonomy ranks. Our experimental results indicate more than 3% improved classification performance of the proposed features with Logistic regression classifier than the GenFooT. We also performed the statistical test to compare the performance of GenFooT2 with the state-of-the-art methods including our previous method GenFooT. The experimental results as well as the statistical test exhibit that the performance of the proposed GenFooT2 is significantly better.

Understanding the Genetic Diversity of Picobirnavirus: A Classification Update Based on Phylogenetic and Pairwise Sequence Comparison Approaches

Picobirnaviruses (PBVs) are small, double stranded RNA viruses with an ability to infect a myriad of hosts and possessing a high degree of genetic diversity. PBVs are currently classified into two genogroups based upon classification of a 200 nt sequence of RdRp. We demonstrate here that this phylogenetic marker is saturated, affected by homoplasy, and has high phylogenetic noise, resulting in 34% unsolved topologies. By contrast, full-length RdRp sequences provide reliable topologies that allow ancestralism of members to be correctly inferred. MAFFT alignment and maximum likelihood trees were established as the optimal methods to determine phylogenetic relationships, providing complete resolution of PBV RdRp and capsid taxa, each into three monophyletic groupings. Pairwise distance calculations revealed these lineages represent three species. For RdRp, the application of cutoffs determined by theoretical taxonomic distributions indicates that there are five genotypes in species 1, eight genotypes in species 2, and three genotypes in species 3. Capsids were also divided into three species, but sequences did not segregate into statistically supported subdivisions, indicating that diversity is lower than RdRp. We thus propose the adoption of a new nomenclature to indicate the species of each segment (e.g., PBV-C1R2).

Identification of a novel filovirus in a common lancehead (Bothrops atrox (Linnaeus, 1758))

I performed metaviromic analysis of publicly available RNA-seq data from reptiles to understand the diversity of filoviruses (family Filoviridae). I identified a coding-complete sequence of a filovirus from the common lancehead (Bothrops atrox (Linnaeus, 1758)), tentatively named Tapajós virus (TAPV). Although the genome organization of TAPV is similar to mammalian filoviruses, our phylogenetic analysis showed that TAPV forms a cluster with a fish filovirus. However, TAPV is still distantly related to all the known filoviruses, suggesting that TAPV can be assigned as a species of a novel genus in Filoviridae. To our knowledge, this is the first report identifying a filovirus in reptiles, and thus contributes to a deeper understanding of the diversity and evolution of filoviruses.

Pantoea Bacteriophage vB_PagS_MED16-A Siphovirus Containing a 2'-Deoxy-7-amido-7-deazaguanosine-Modified DNA

A novel siphovirus, vB_PagS_MED16 (MED16) was isolated in Lithuania using Pantoea agglomerans strain BSL for the phage propagation. The double-stranded DNA genome of MED16 (46,103 bp) contains 73 predicted open reading frames (ORFs) encoding proteins, but no tRNA. Our comparative sequence analysis revealed that 26 of these ORFs code for unique proteins that have no reliable identity when compared to database entries. Based on phylogenetic analysis, MED16 represents a new genus with siphovirus morphology. In total, 35 MED16 ORFs were given a putative functional annotation, including those coding for the proteins responsible for virion morphogenesis, phage-host interactions, and DNA metabolism. In addition, a gene encoding a preQ0 DNA deoxyribosyltransferase (DpdA) is present in the genome of MED16 and the LC-MS/MS analysis indicates 2'-deoxy-7-amido-7-deazaguanosine (dADG)-modified phage DNA, which, to our knowledge, has never been experimentally validated in genomes of Pantoea phages. Thus, the data presented in this study provide new information on Pantoea-infecting viruses and offer novel insights into the diversity of DNA modifications in bacteriophages.

Multiple Mammarenaviruses Circulating in Angolan Rodents

Rodents are a speciose group of mammals with strong zoonotic potential. Some parts of Africa are still underexplored for the occurrence of rodent-borne pathogens, despite this high potential. Angola is at the convergence of three major biogeographical regions of sub-Saharan Africa, each harbouring a specific rodent community. This rodent-rich area is, therefore, strategic for studying the diversity and evolution of rodent-borne viruses. In this study we examined 290 small mammals, almost all rodents, for the presence of mammarenavirus and hantavirus RNA. While no hantavirus was detected, we found three rodent species positive for distinct mammarenaviruses with a particularly high prevalence in Namaqua rock rats (Micaelamys namaquensis). We characterised four complete virus genomes, which showed typical mammarenavirus organisation. Phylogenetic and genetic distance analyses revealed: (i) the presence of a significantly divergent strain of Luna virus in Angolan representatives of the ubiquitous Natal multimammate mouse (Mastomys natalensis), (ii) a novel Okahandja-related virus associated with the Angolan lineage of Micaelamys namaquensis for which we propose the name Bitu virus (BITV) and (iii) the occurrence of a novel Mobala-like mammarenavirus in the grey-bellied pygmy mouse (Mus triton) for which we propose the name Kwanza virus (KWAV). This high virus diversity in a limited host sample size and in a relatively small geographical area supports the idea that Angola is a hotspot for mammarenavirus diversity.

Genome Sequence of Ruloma Virus, a Novel Paramyxovirus Clustering Basally to Members of the Genus Jeilongvirus

We report here the complete genome sequence of ruloma virus, a novel paramyxovirus detected in a Machangu’s brush-furred rat from Tanzania. Ruloma virus has the longest orthoparamyxovirus genome reported to date and forms a sister clade to all currently known members of the genus Jeilongvirus.

We report here the complete genome sequence of ruloma virus, a novel paramyxovirus detected in a Machangu’s brush-furred rat from Tanzania. Ruloma virus has the longest orthoparamyxovirus genome reported to date and forms a sister clade to all currently known members of the genus Jeilongvirus.

A Primer on the Analysis of High-Throughput Sequencing Data for Detection of Plant Viruses

High-throughput sequencing (HTS) technologies have become indispensable tools assisting plant virus diagnostics and research thanks to their ability to detect any plant virus in a sample without prior knowledge. As HTS technologies are heavily relying on bioinformatics analysis of the huge amount of generated sequences, it is of utmost importance that researchers can rely on efficient and reliable bioinformatic tools and can understand the principles, advantages, and disadvantages of the tools used. Here, we present a critical overview of the steps involved in HTS as employed for plant virus detection and virome characterization. We start from sample preparation and nucleic acid extraction as appropriate to the chosen HTS strategy, which is followed by basic data analysis requirements, an extensive overview of the in-depth data processing options, and taxonomic classification of viral sequences detected. By presenting the bioinformatic tools and a detailed overview of the consecutive steps that can be used to implement a well-structured HTS data analysis in an easy and accessible way, this paper is targeted at both beginners and expert scientists engaging in HTS plant virome projects.

Molecular and Phylogenetic Characterization of Novel Papillomaviruses Isolated from Oral and Anogenital Neoplasms of Japanese Macaques (Macaca fuscata)

Papillomaviruses (PVs) are a diverse group of host species-specific DNA viruses, etiologically linked with various benign and malignant neoplasms of cutaneous and mucosal epithelia. Here, we describe the detection and characterization of the first two PVs naturally infecting Japanese macaques (Macaca fuscata), including the determination of their etiological association(s) with the development of original neoplasms. The molecular and phylogenetic analyses were performed on complete genome sequences of Macaca fuscata PV types 1 (MfuPV1) and 2 (MfuPV2), which were completely sequenced in samples of a malignant oral tumor and benign anogenital neoplasm of Japanese macaques, respectively. Subsequently, two type-specific quantitative real-time PCRs were developed to estimate viral loads of MfuPV1 and MfuPV2 and to evaluate their etiological roles. The in silico molecular analyses revealed that both viral genomes encode characteristic PV proteins with conserved functional domains and have a non-coding genomic region with regulatory sequences to regulate and complete the viral life cycle. However, additional experimental evidence is needed to finally confirm the presence and biological functionality of the molecular features of both novel PVs. While MfuPV1, together with PVs identified in other macaques, is classified into the Alphapapillomavirus (Alpha-PV) species 12, MfuPV2 is most likely a representative of the novel viral species within the Alpha-PV genus. Their relatively high viral loads suggest that both PVs are etiologically linked with the development of the original neoplasms.

Characterization of a novel roseophage and the morphological and transcriptional response of the sponge symbiont Ruegeria AU67 to infection

Sponges have recently been recognized to contain complex communities of bacteriophages; however, little is known about how they interact with their bacterial hosts. Here, we isolated a novel phage, called Ruegeria phage Tedan, and characterized its impact on the bacterial sponge symbiont Ruegeria AU67 on a morphological and molecular level. Phage Tedan was structurally, genomically and phylogenetically characterized to be affiliated with the genus Xiamenvirus of the family Siphoviridae. Through microscopic observations and transcriptomic analysis, we show that phage Tedan upon infection induces a process leading to metabolic and morphological changes in its host. These changes would render Ruegeria AU67 better adapted to inhabit the sponge holobiont due to an improved utilization of ecologically relevant energy and carbon sources as well as a potential impediment of phagocytosis by the sponge through cellular enlargement. An increased survival or better growth of the bacterium in the sponge environment will likely benefit the phage reproduction. Our results point towards the possibility that phages from host-associated environments require, and have thus evolved, different strategies to interact with their host when compared to those phages from free-living or planktonic environments.

Pantoea Bacteriophage vB_PagS_AAS23: A Singleton of the Genus Sauletekiovirus

A cold-adapted siphovirus, vB_PagS_AAS23 (AAS23) was isolated in Lithuania using the Pantoea agglomerans strain AUR for the phage propagation. The double-stranded DNA genome of AAS23 (51,170 bp) contains 92 probable protein encoding genes, and no genes for tRNA. A comparative sequence analysis revealed that 25 of all AAS23 open reading frames (ORFs) code for unique proteins that have no reliable identity to database entries. Based on the phylogenetic analysis, AAS23 has no close relationship to other viruses publicly available to date and represents a single species of the genus Sauletekiovirus within the family Drexlerviridae. The phage is able to form plaques in bacterial lawns even at 4 °C and demonstrates a depolymerase activity. Thus, the data presented in this study not only provides the information on Pantoea-infecting bacteriophages, but also offers novel insights into the diversity of cold-adapted viruses and their potential to be used as biocontrol agents.

A Mixed Infection of Helenium Virus S With Two Distinct Isolates of Butterbur Mosaic Virus, One of Which Has a Major Deletion in an Essential Gene

Multiple carlaviruses infect various ornamental plants, often having limited host ranges and causing minor symptoms, yet often reducing yield or quality. In this study we have identified a mixed infection of butterbur mosaic virus (ButMV) and helenium virus S (HelVS) from a plant of veronica (Veronica sp.) showing foliar mosaic and distortion. Carlavirus-like particles were observed by transmission electron microscopy (TEM), and RNA from partially purified virions was amplified by random RT-PCR, yielding clones of 439–1,385 bp. Two partially overlapping clones including coat protein (CP) sequence, and two of four partial replicase clones, were closely related to ButMV-J (AB517596), previously reported only from butterbur (Petasites japonicus) in Japan. Two other partial replicase clones showed lower identity to multiple carlaviruses. Generic primers which amplify the 3′-terminal region of multiple carlaviruses yielded clones of three distinct sequences: (1) with 98% nt identity to HelVS; (2) ButMV-A, showing 82% nt identity to ButMV-J; and (3) ButMV-B, with 78% nt identity to each of ButMV-J and ButMV-A. Further amplification of upstream fragments revealed that ButMV-B had an internal deletion in TGB1, confirmed using isolate-specific primers. Near-complete genomes of both ButMV-A and ButMV-B were obtained from next-generation sequencing (NGS), confirming the deletion within ButMV-B, which is presumably maintained through complementation by ButMV-A. HelVS was previously reported only from Helenium hybrids and Impatiens holstii. A near-complete HelVS genome was obtained for the first time by NGS from the same sample. Additional Veronica hybrids infected with HelVS were identified by TEM and RT-PCR, including cv. ‘Sunny Border Blue’ which was also subjected to NGS. This resulted in assembly of an 8,615 nt near-complete HelVS genome, with high identity to that from the mixed infection. The predicted CP sequence has 96% amino acid (aa) identity to HelVS from helenium (Q00556). Other ORFs show a maximum of 54% (TGB3) to 68% (NABP) aa identity to the equivalent ORFs of other carlaviruses. These results demonstrate for the first time maintenance by complementation of a carlavirus isolate with a major deletion in an essential gene, and confirm that HelVS is a distinct species in the genus Carlavirus.

VIRIDIC-A Novel Tool to Calculate the Intergenomic Similarities of Prokaryote-Infecting Viruses

Nucleotide-based intergenomic similarities are useful to understand how viruses are related with each other and to classify them. Here we have developed VIRIDIC, which implements the traditional algorithm used by the International Committee on Taxonomy of Viruses (ICTV), Bacterial and Archaeal Viruses Subcommittee, to calculate virus intergenomic similarities. When compared with other software, VIRIDIC gave the best agreement with the traditional algorithm, which is based on the percent identity between two genomes determined by BLASTN. Furthermore, VIRIDIC proved best at estimating the relatedness between more distantly-related phages, relatedness that other tools can significantly overestimate. In addition to the intergenomic similarities, VIRIDIC also calculates three indicators of the alignment ability to capture the relatedness between viruses: the aligned fractions for each genome in a pair and the length ratio between the two genomes. The main output of VIRIDIC is a heatmap integrating the intergenomic similarity values with information regarding the genome lengths and the aligned genome fraction. Additionally, VIRIDIC can group viruses into clusters, based on user-defined intergenomic similarity thresholds. The sensitivity of VIRIDIC is given by the BLASTN. Thus, it is able to capture relationships between viruses having in common even short genomic regions, with as low as 65% similarity. Below this similarity level, protein-based analyses should be used, as they are the best suited to capture distant relationships. VIRIDIC is available at viridic.icbm.de, both as a web-service and a stand-alone tool. It allows fast analysis of large phage genome datasets, especially in the stand-alone version, which can be run on the user's own servers and can be integrated in bioinformatics pipelines. VIRIDIC was developed having viruses of Bacteria and Archaea in mind; however, it could potentially be used for eukaryotic viruses as well, as long as they are monopartite.

Genomic analysis of Ranavirus and exploring alternative genes for phylogenetics

Ranaviruses can infect both captive and wild cold-blooded vertebrates, leading to significant economic and environmental losses. With the cases of ranavirus infection increasing, many ranavirus genomic sequences were published, but little is known about ranavirus taxonomy on a whole-genome level. In this study, 44 ranaviruses core genes were identified in 32 ranaviruses genome sequences by using PanX. The neighbour-joining phylogenetic trees (NJ-tree) based on 44 ranaviruses core genes and 24 iridoviridae core genes and composition vector phylogenetic tree (CV-Tree) based on whole genome were constructed. The three of phylogenetic trees showed that 32 ranavirus isolates can be divided into 4 different subgroups including SGIV-like, EHNV-like, FV3-like and CMTV-like, and subgroups taxonomic position of three phylogenetic trees were consistent. However, the phylogenetic position of ToRV could not be determined if it belongs to FV3-like or CMTV-like group. Subsequently, we carried out dot plot analysis and confirmed that ToRV should belong to CMTV-like group. Based on dot plot analysis and phylogenetic trees, the taxonomic classification of ranaviruses was confirmed. Finally, four genes which are suitable for the construction of phylogenetic tree were selected from ranavirus core genes by recombination analysis, substitution saturation analysis and single-gene phylogenetic analysis. Phylogenetic tree based on concatenated sequences of the four selected genes showed that the classification of subgroups was identical with three of the phylogenetic trees. Conclusion: Our results confirmed taxonomic identification of ranaviruses; the four selected genes used in phylogenetic analysis will make taxonomic identification more convenient and accurate.

Metatranscriptomic Identification of Diverse and Divergent RNA Viruses in Green and Chlorarachniophyte Algae Cultures

Our knowledge of the diversity and evolution of the virosphere will likely increase dramatically with the study of microbial eukaryotes, including the microalgae within which few RNA viruses have been documented. By combining total RNA sequencing with sequence and structural-based homology detection, we identified 18 novel RNA viruses in cultured samples from two major groups of microbial algae: the chlorophytes and the chlorarachniophytes. Most of the RNA viruses identified in the green algae class Ulvophyceae were related to the Tombusviridae and Amalgaviridae viral families commonly associated with land plants. This suggests that the evolutionary history of these viruses extends to divergence events between algae and land plants. Seven Ostreobium sp-associated viruses exhibited sequence similarity to the mitoviruses most commonly found in fungi, compatible with horizontal virus transfer between algae and fungi. We also document, for the first time, RNA viruses associated with chlorarachniophytes, including the first negative-sense (bunya-like) RNA virus in microalgae, as well as a distant homolog of the plant virus Virgaviridae, potentially signifying viral inheritance from the secondary chloroplast endosymbiosis that marked the origin of the chlorarachniophytes. More broadly, these data suggest that the scarcity of RNA viruses in algae results from limited investigation rather than their absence.

Metatranscriptomic Identification of Diverse and Divergent RNA Viruses in Green and Chlorarachniophyte Algae Cultures

Our knowledge of the diversity and evolution of the virosphere will likely increase dramatically with the study of microbial eukaryotes, including the microalgae within which few RNA viruses have been documented. By combining total RNA sequencing with sequence and structural-based homology detection, we identified 18 novel RNA viruses in cultured samples from two major groups of microbial algae: the chlorophytes and the chlorarachniophytes. Most of the RNA viruses identified in the green algae class Ulvophyceae were related to the Tombusviridae and Amalgaviridae viral families commonly associated with land plants. This suggests that the evolutionary history of these viruses extends to divergence events between algae and land plants. Seven Ostreobium sp-associated viruses exhibited sequence similarity to the mitoviruses most commonly found in fungi, compatible with horizontal virus transfer between algae and fungi. We also document, for the first time, RNA viruses associated with chlorarachniophytes, including the first negative-sense (bunya-like) RNA virus in microalgae, as well as a distant homolog of the plant virus Virgaviridae, potentially signifying viral inheritance from the secondary chloroplast endosymbiosis that marked the origin of the chlorarachniophytes. More broadly, these data suggest that the scarcity of RNA viruses in algae results from limited investigation rather than their absence.

Characterization of virus species associated with sweetpotato virus diseases in Burkina Faso

Sweetpotato (Ipomoea batatas) production in sub-Saharan Africa is severely affected by viral diseases caused by several interacting viruses, including sweet potato feathery mottle virus (SPFMV), sweet potato chlorotic stunt virus (SPCSV), and sweet potato leaf curl virus (SPLCV). However, the aetiology of viral symptoms on sweetpotato is rarely established in most countries in Africa. Here, we aimed to investigate and characterize the incidence of sweetpotato viruses in Burkina Faso. We performed a countrywide survey in 18 districts of Burkina Faso and collected 600 plants, with and without symptoms, from 80 fields. Viral strains were identified using nitrocellulose membrane-ELISA, PCR, and reverse transcription-PCR. Three scions from each of 50 selected plants with symptoms were grafted to healthy Ipomoea setosa and then serological and molecular tests were performed on the 150 recorded samples. Three viruses were detected: 24% of samples were positive for SPFMV, 18% for SPLCV, and 2% for SPCSV. Across all diagnostic tests, 40% of all plant samples were virus-negative. Coinfections were found in 16% of samples. Partial sequences were obtained, including 13 that matched SPFMV, one that matched SPLCV, and one that matched SPCSV. All identified SPFMV isolates belonged to either phylogroup B or A-II. The SPCSV-positive isolates had 98% gene sequence homology with SPCSV-West Africa for the coat protein. Begomovirus-positive isolates clustered with SPLCV-United States. This first study of sweetpotato viral diseases in Burkina Faso indicates widespread occurrence and suggests a need for further epidemiological investigations, breeding programmes focused on virus-resistant varieties, and improved farming practices to control disease spread.

Genome Analysis of Melon Necrotic Spot Virus Incursions and Seed Interceptions in Australia

Melon necrotic spot virus (MNSV) was detected in field-grown Cucumis melo (rockmelon) and Citrullus lanatus (watermelon) plants in the Sunraysia district of New South Wales and Victoria, Australia, in 2012, 2013, and 2016, and in two watermelon seed lots tested at the Australian border in 2016. High-throughput sequencing was used to generate near full-length genomes of six isolates detected during the incursions and seed testing. Phylogenetic analysis of the genomes suggests that there have been at least two incursions of MNSV into Australia and none of the field isolates were the same as the isolates detected in seeds. The analysis indicated that one watermelon field sample (L10), the Victorian rockmelon field sample, and two seed interception samples may have European origins. The results showed that two isolates (L8 and L9) from watermelon were divergent from the type MNSV strain (MNSV-GA, D12536.2) and had 99% nucleotide identity to two MNSV isolates from human stool collected in the United States (KY124135.1, KY124136.1). These isolates also had high nucleotide pairwise identity (96%) to a partial sequence from a Spanish MNSV isolate (KT962848.1). The analysis supported the identification of three previously described MNSV genotype groups: EU-LA, Japan melon, and Japan watermelon. To account for the greater diversity of hosts and geographic regions of the MNSV isolates used in this study, it is suggested that the genotype groups EU-LA, Japan melon, and Japan watermelon be renamed to groups I, II, and III, respectively. The divergent isolates L8 and L9 from this study and the stool isolates from the United States formed a fourth genotype group, group IV. Soil collected from the site of the Victorian rockmelon MNSV outbreak was found to contain viable MNSV and the virus vector, a chytrid fungus, Olpidium bornovanus (Sahtiyanci) Karling, 18 months after the initial MNSV detection. This is a first report of O. bornovanus from soil sampled from an MNSV-contaminated site in Australia.

Complete genome analysis of Pantoea agglomerans-infecting bacteriophage vB_PagM_AAM22

A novel myovirus, vB_PagM_AAM22 (AAM22), was isolated in Lithuania using Pantoea agglomerans as the host for phage propagation. The 49,744-bp genome of AAM22 has a G + C content of 48.4% and contains 96 probable protein-encoding genes and no genes for tRNA. In total, 34 ORFs were given a putative functional annotation, including genes associated with virion morphogenesis, DNA metabolism, and phage-host interactions. Based on comparative phylogenetic analysis, AAM22 cannot be assigned to any genus currently recognized by the ICTV and is a potential candidate to form a new genus within the family Myoviridae.