Vector-Enabled Metagenomic (VEM) Surveys Using Whiteflies (Aleyrodidae) Reveal Novel Begomovirus Species in the New and Old Worlds

Phylogeographic analysis of Begomovirus coat and replication-associated proteins

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

From Metagenomics to Ecogenomics: NGS-Based Approaches for Discovery of New Circular DNA Single-Stranded Viral Species

Viruses comprise the most abundant genetic material in the biosphere; however, global viral genomic population (virome) has been largely underestimated. Recently, high-throughput sequencing (HTS) has provided a powerful tool for the detection of known viruses and the discovery of novel viral species from environmental and individual samples using metagenomics and ecogenomics approaches, respectively. Viruses with circular DNA single-stranded (ssDNA) genomes belonging to the begomovirus genera (family Geminiviridae) constitute the largest group of emerging plant viruses worldwide. The knowledge of begomoviruses viromes is mostly restricted to crop plant systems; nevertheless, it has been described that noncultivated plants specifically at the interface between wild and cultivated plants are important reservoirs leading to viral evolution and the emergence of new diseases. Here we present a protocol that allows the identification and isolation of known and novel begomoviruses species infecting cultivated and noncultivated plant species. The method consists of circular viral molecules enrichment by rolling circle amplification (RCA) from begomovirus-positive total plant DNA, followed by NGS-based metagenomic sequencing. Subsequently, metagenomic reads are processed for taxonomic classification using Viromescan software and a customized Geminiviridae family database, and begomovirus-related reads are used for contigs assembly and annotation using Spades software and Blastn algorithm, respectively. Then, the obtained begomovirus-related signatures are used as templates for specific primers design and implemented for PCR-based ecogenomic identification of individual samples harboring the corresponding viral species. Lastly, full-length begomovirus genomes are obtained by RCA-based amplification from total plant DNA of selected individual samples, cloning, and viral molecular identity corroborated by Sanger sequencing. Conclusively, the identification and isolation of a novel monopartite begomovirus species native to the New World (NW) named Gallium leaf deformation virus (GLDV) is shown.

Translating virome analyses to support biosecurity, on-farm management, and crop breeding

Virome analysis via high-throughput sequencing (HTS) allows rapid and massive virus identification and diagnoses, expanding our focus from individual samples to the ecological distribution of viruses in agroecological landscapes. Decreases in sequencing costs combined with technological advances, such as automation and robotics, allow for efficient processing and analysis of numerous samples in plant disease clinics, tissue culture laboratories, and breeding programs. There are many opportunities for translating virome analysis to support plant health. For example, virome analysis can be employed in the development of biosecurity strategies and policies, including the implementation of virome risk assessments to support regulation and reduce the movement of infected plant material. A challenge is to identify which new viruses discovered through HTS require regulation and which can be allowed to move in germplasm and trade. On-farm management strategies can incorporate information from high-throughput surveillance, monitoring for new and known viruses across scales, to rapidly identify important agricultural viruses and understand their abundance and spread. Virome indexing programs can be used to generate clean germplasm and seed, crucial for the maintenance of seed system production and health, particularly in vegetatively propagated crops such as roots, tubers, and bananas. Virome analysis in breeding programs can provide insight into virus expression levels by generating relative abundance data, aiding in breeding cultivars resistant, or at least tolerant, to viruses. The integration of network analysis and machine learning techniques can facilitate designing and implementing management strategies, using novel forms of information to provide a scalable, replicable, and practical approach to developing management strategies for viromes. In the long run, these management strategies will be designed by generating sequence databases and building on the foundation of pre-existing knowledge about virus taxonomy, distribution, and host range. In conclusion, virome analysis will support the early adoption and implementation of integrated control strategies, impacting global markets, reducing the risk of introducing novel viruses, and limiting virus spread. The effective translation of virome analysis depends on capacity building to make benefits available globally.

Novel Virus Identification through Metagenomics: A Systematic Review

Metagenomic Next Generation Sequencing (mNGS) allows the evaluation of complex microbial communities, avoiding isolation and cultivation of each microbial species, and does not require prior knowledge of the microbial sequences present in the sample. Applications of mNGS include virome characterization, new virus discovery and full-length viral genome reconstruction, either from virus preparations enriched in culture or directly from clinical and environmental specimens. Here, we systematically reviewed studies that describe novel virus identification through mNGS from samples of different origin (plant, animal and environment). Without imposing time limits to the search, 379 publications were identified that met the search parameters. Sample types, geographical origin, enrichment and nucleic acid extraction methods, sequencing platforms, bioinformatic analytical steps and identified viral families were described. The review highlights mNGS as a feasible method for novel virus discovery from samples of different origins, describes which kind of heterogeneous experimental and analytical protocols are currently used and provides useful information such as the different commercial kits used for the purification of nucleic acids and bioinformatics analytical pipelines.

An experimental strategy for preparing circular ssDNA virus genomes for next-generation sequencing

The ability of begomoviruses to evolve rapidly threatens many crops and underscores the importance of detecting these viruses quickly and to understand their genome diversity. This study presents an improved protocol for the enhanced amplification and enrichment of begomovirus DNA for use in next generation sequencing of the viral genomes. An enhanced rolling circle amplification (RCA) method using EquiPhi29 polymerase was combined with size selection to generate a cost-effective, short-read sequencing method. This improved short-read sequencing produced at least 50 % of the reads mapping to the target viral reference genomes, African cassava mosaic virus and East African cassava mosaic virus. This study provided other insights into common misconceptions about RCA and lessons that could be learned from the sequencing of single-stranded DNA virus genomes. This protocol can be used to examine the viral DNA as it moves from host to vector, thus producing valuable information for viral DNA population studies, and would likely work well with other circular Rep-encoding ssDNA viruses (CRESS) DNA viruses.

A Bipartite Geminivirus with a Highly Divergent Genomic Organization Identified in Olive Trees May Represent a Novel Evolutionary Direction in the Family Geminiviridae

Olea europaea Geminivirus (OEGV) was recently identified in olive in Italy through HTS. In this work, we used HTS to show the presence of an OEGV isolate in Portuguese olive trees and suggest the evolution direction of OEGV. The bipartite genome (DNA-A and DNA-B) of the OEGV-PT is similar to Old World begomoviruses in length, but it lacks a pre-coat protein (AV2), which is a typical feature of New World begomoviruses (NW). DNA-A genome organization is closer to NW, containing four ORFs; three in complementary-sense AC1/Rep, AC2/TrAP, AC3/REn and one in virion-sense AV1/CP, but no AC4, typical of begomoviruses. DNA-B comprises two ORFs; MP in virion sense with higher similarity to the tyrosine phosphorylation site of NW, but in opposite sense to begomoviruses; BC1, with no known conserved domains in the complementary sense and no NSP typical of bipartite begomoviruses. Our results show that OEGV presents the longest common region among the begomoviruses, and the TATA box and four replication-associated iterons in a completely new arrangement. We propose two new putative conserved regions for the geminiviruses CP. Lastly, we highlight unique features that may represent a new evolutionary direction for geminiviruses and suggest that OEGV-PT evolution may have occurred from an ancient OW monopartite Begomovirus that lost V2 and C4, gaining functions on cell-to-cell movement by acquiring a DNA-B component.

Characterisation of the Viral Community Associated with the Alfalfa Weevil (Hypera postica) and Its Host Plant, Alfalfa (Medicago sativa)

Advances in viral metagenomics have paved the way of virus discovery by making the exploration of viruses in any ecosystem possible. Applied to agroecosystems, such an approach opens new possibilities to explore how viruses circulate between insects and plants, which may help to optimise their management. It could also lead to identifying novel entomopathogenic viral resources potentially suitable for biocontrol strategies. We sampled the larvae of a natural population of alfalfa weevils (Hypera postica), a major herbivorous pest feeding on legumes, and its host plant alfalfa (Medicago sativa). Insect and plant samples were collected from a crop field and an adjacent meadow. We characterised the diversity and abundance of viruses associated with weevils and alfalfa, and described nine putative new virus species, including four associated with alfalfa and five with weevils. In addition, we found that trophic accumulation may result in a higher diversity of plant viruses in phytophagous pests compared to host plants.

Soybean Thrips (Thysanoptera: Thripidae) Harbor Highly Diverse Populations of Arthropod, Fungal and Plant Viruses

Soybean thrips (Neohydatothrips variabilis) are one of the most efficient vectors of soybean vein necrosis virus, which can cause severe necrotic symptoms in sensitive soybean plants. To determine which other viruses are associated with soybean thrips, the metatranscriptome of soybean thrips, collected by the Midwest Suction Trap Network during 2018, was analyzed. Contigs assembled from the data revealed a remarkable diversity of virus-like sequences. Of the 181 virus-like sequences identified, 155 were novel and associated primarily with taxa of arthropod-infecting viruses, but sequences similar to plant and fungus-infecting viruses were also identified. The novel viruses were predicted to have positive-sense RNA, negative-stranded RNA, double-stranded RNA, and single-stranded DNA genomes. The assembled sequences included 100 contigs that represented at least 95% coverage of a virus genome or genome segment. Sequences represented 12 previously described arthropod viruses including eight viruses reported from Hubei Province in China, and 12 plant virus sequences of which six have been previously described. The presence of diverse populations of plant viruses within soybean thrips suggests they feed on and acquire viruses from multiple host plant species that could be transmitted to soybean. Assessment of the virome of soybean thrips provides, for the first time, information on the diversity of viruses present in thrips.

Reconstruction and Characterization of Full-Length Begomovirus and Alphasatellite Genomes Infecting Pepper through Metagenomics

In northwestern Argentina (NWA), pepper crops are threatened by the emergence of begomoviruses due to the spread of its vector, Bemisia tabaci (Gennadius). The genus Begomovirus includes pathogens that can have a monopartite or bipartite genome and are occasionally associated with sub-viral particles called satellites. This study characterized the diversity of begomovirus and alphasatellite species infecting pepper in NWA using a metagenomic approach. Using RCA-NGS (rolling circle amplification-next generation sequencing), 19 full-length begomovirus genomes (DNA-A and DNA-B) and one alphasatellite were assembled. This ecogenomic approach revealed six begomoviruses in single infections: soybean blistering mosaic virus (SbBMV), tomato yellow spot virus (ToYSV), tomato yellow vein streak virus (ToYVSV), tomato dwarf leaf virus (ToDfLV), sida golden mosaic Brazil virus (SiGMBRV), and a new proposed species, named pepper blistering leaf virus (PepBLV). SbBMV was the most frequently detected species, followed by ToYSV. Moreover, a new alphasatellite associated with ToYSV, named tomato yellow spot alphasatellite 2 (ToYSA-2), was reported for the first time in Argentina. For the Americas, this was the first report of an alphasatellite found in a crop (pepper) and in a weed (Leonurus japonicus). We also detected intra-species and inter-species recombination.

High-Throughput Sequencing Reveals Differential Begomovirus Species Diversity in Non-Cultivated Plants in Northern-Pacific Mexico

Plant DNA viruses of the genus Begomovirus have been documented as the most genetically diverse in the family Geminiviridae and present a serious threat for global horticultural production, especially considering climate change. It is important to characterize naturally existing begomoviruses, since viral genetic diversity in non-cultivated plants could lead to future disease epidemics in crops. In this study, high-throughput sequencing (HTS) was employed to determine viral diversity of samples collected in a survey performed during 2012–2016 in seven states of Northern-Pacific Mexico, areas of diverse climatic conditions where different vegetable crops are subject to intensive farming. In total, 132 plant species, belonging to 34 families, were identified and sampled in the natural ecosystems surrounding cultivated areas (agro-ecological interface). HTS analysis and subsequent de novo assembly revealed a number of geminivirus-related DNA signatures with 80 to 100% DNA similarity with begomoviral sequences present in the genome databank. The analysis revealed DNA signatures corresponding to 52 crop-infecting and 35 non-cultivated-infecting geminiviruses that, interestingly, were present in different plant species. Such an analysis deepens our knowledge of geminiviral diversity and could help detecting emerging viruses affecting crops in different agro-climatic regions.

Molecular characterization and phylogenetic analysis of a Squash leaf curl virus isolate from Baja California Sur, Mexico

Background

The begomovirus, squash leaf curl virus (SLCuV) is one of the causal agents of squash leaf curl (SLC) disease, which is among the most destructive diseases of cucurbit crops in tropical, subtropical, and semiarid regions worldwide. This disease was originally reported in the American continent with subsequent spread to the Mediterranean basin. Up to now, SLCuV has only been detected by PCR in Mexico. This study provides the first complete sequence of a Mexican SLCuV isolate from Baja California Sur (BCS). In addition, the genome of the virus was characterized, establishing its phylogenetic relationship with other SLCuV isolates.

Methods

The full genome (DNA-A and DNA-B) was amplified by rolling circle amplification, cloned and sequenced and the open reading frames (ORF) were annotated. Virus identification was performed according to the International Committee on Taxonomy of Viruses (ICTV) criteria for begomovirus species demarcation. To infer evolutionary relationship with other SLCuV isolates, phylogenetic and recombination analyses were performed.

Results

The SLCuV-[MX-BCS-La Paz-16] genome (DNA-A and DNA-B) had 99% identity with SLCuV reference genomes. The phylogenetic analysis showed that SLCuV-[MX-BCS-La Paz-16] is closely related to SLCuV isolates from the Middle East (Egypt, Israel, Palestine and Lebanon). No evidence of interspecific recombination was determined and iterons were 100% identical in all isolates in the SLCuV clade.

Conclusions

SLCuV-[MX-BCS-La Paz-16] showed low genetic variability in its genome, which could be due to a local adaptation process (isolate environment), suggesting that SLCuV isolates from the Middle East could have derived from the southwestern United States of America (USA) and northwestern Mexico.

Genome wide molecular evolution analysis of begomoviruses reveals unique diversification pattern in coat protein gene of Old World and New World viruses

Begomoviruses (Family-Geminiviridae) are plant infecting single stranded DNA viruses known to evolve very fast. Here, we have analysed the DNA-A sequences of 302 begomoviruses reported as 'type isolates' from different countries following the list of International Committee on Taxonomy of Viruses till 2017. Phylogenetic analysis was performed which revealed two major evolutionarily distinct groups namely Old World (OW) and New World (NW) viruses. Our work present evidence that cp gene has varied degree of diversification among the viruses reported from NW and OW. The NW viruses are more conserved in their cp gene sequences than that of OW viruses irrespective of host plant families. Further analysis reveals that cp gene differs in its recombination pattern among OW and NW viruses whereas rep gene is highly recombination prone in both OW and NW viruses. The sequence conservation in cp gene in NW viruses is a result of meagre recombination and subsequent low substitution rate in comparison to OW viruses. Our results demonstrated that the cp gene in NW viruses is less likely to possess nuclear localisation sequences than OW cp gene. Further we present evidence that the NW-cp is under the influence of strong purifying selection. We propose that the precoat protein (pcp) gene present exclusively in the 5' of cp gene in OW viruses is highly diversified and strong positive selection working on pcp gene might be attributing largely to the diversity of OW-cp gene.

Tomato Twisted Leaf Virus: A Novel Indigenous New World Monopartite Begomovirus Infecting Tomato in Venezuela

Begomoviruses are one of the major groups of plant viruses with an important economic impact on crop production in tropical and subtropical regions. The global spread of its polyphagous vector, the whitefly Bemisia tabaci, has contributed to the emergence and diversification of species within this genus. In this study, we found a putative novel begomovirus infecting tomato plants in Venezuela without a cognate DNA-B component. This begomovirus was genetically characterized and compared with related species. Furthermore, its infectivity was demonstrated by agroinoculation of infectious clones in tomato (Solanum lycopersicum) and Nicotiana benthamiana plants. The name Tomato twisted leaf virus (ToTLV) is proposed. ToTLV showed the typical genome organization of the DNA-A component of New World bipartite begomoviruses. However, the single DNA component of ToTLV was able to develop systemic infection in tomato and N. benthamiana plants, suggesting a monopartite nature of its genome. Interestingly, an additional open reading frame ORF was observed in ToTLV encompassing the intergenic region and the coat protein gene, which is not present in other closely related begomoviruses. A putative transcript from this region was amplified by strand-specific reverse transcription-PCR. Along with recent studies, our results showed that the diversity of monopartite begomoviruses from the New World is greater than previously thought.

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

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

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

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

From Spatial Metagenomics to Molecular Characterization of Plant Viruses: A Geminivirus Case Study

The number of plant viruses that are known likely remains only a vanishingly small fraction of all extant plant virus species. Consequently, the distribution and population dynamics of plant viruses within even the best-studied ecosystems have only ever been studied for small groups of virus species. Even for the best studied of these groups very little is known about virus diversity at spatial scales ranging from an individual host, through individual local host populations to global host populations. To date, metagenomics studies that have assessed the collective or metagenomes of viruses at the ecosystem scale have revealed many previously unrecognized viral species. More recently, novel georeferenced metagenomics approaches have been devised that can precisely link individual sequence reads to both the plant hosts from which they were obtained, and the spatial arrangements of these hosts. Besides illuminating the diversity and the distribution of plant viruses at the ecosystem scale, application of these "geometagenomics" approaches has enabled the direct testing of hypotheses relating to the impacts of host diversity, host spatial variations, and environmental conditions on plant virus diversity and prevalence. To exemplify how such top-down approaches can provide a far deeper understanding of host-virus associations, we provide a case-study focusing on geminiviruses within two complex ecosystems containing both cultivated and uncultivated areas. Geminiviruses are a highly relevant model for studying the evolutionary and ecological aspects of viral emergence because the family Geminiviridae includes many of the most important crop pathogens that have emerged over the past century. In addition to revealing unprecedented degrees of geminivirus diversity within the analyzed ecosystems, the geometagenomics-based approach enabled the focused in-depth analysis of the complex evolutionary dynamics of some of the highly divergent geminivirus species that were discovered.

The first transcriptomes from field-collected individual whiteflies ( Bemisia tabaci, Hemiptera: Aleyrodidae): a case study of the endosymbiont composition

Background: Bemisia tabaci species ( B. tabaci), or whiteflies, are the world's most devastating insect pests. They cause billions of dollars (US) of damage each year, and are leaving farmers in the developing world food insecure. Currently, all publically available transcriptome data for B. tabaci are generated from pooled samples, which can lead to high heterozygosity and skewed representation of the genetic diversity. The ability to extract enough RNA from a single whitefly has remained elusive due to their small size and technological limitations. Methods: In this study, we optimised a single whitefly RNA extraction procedure, and sequenced the transcriptome of four individual adult Sub-Saharan Africa 1 (SSA1) B. tabaci. Transcriptome sequencing resulted in 39-42 million raw reads. De novo assembly of trimmed reads yielded between 65,000-162,000 Contigs across B. tabaci transcriptomes. Results: Bayesian phylogenetic analysis of mitochondrion cytochrome I oxidase (mtCOI) grouped the four whiteflies within the SSA1 clade. BLASTn searches on the four transcriptomes identified five endosymbionts; the primary endosymbiont Portiera aleyrodidarum and four secondary endosymbionts: Arsenophonus, Wolbachia, Rickettsia, and Cardinium spp. that were predominant across all four SSA1 B. tabaci samples with prevalence levels of between 54.1 to 75%. Amino acid alignments of the NusG gene of P. aleyrodidarum for the SSA1 B. tabaci transcriptomes of samples WF2 and WF2b revealed an eleven amino acid residue deletion that was absent in samples WF1 and WF2a. Comparison of the protein structure of the NusG protein from P. aleyrodidarum in SSA1 with known NusG structures showed the deletion resulted in a shorter D loop. Conclusions: The use of field-collected specimens means time and money will be saved in future studies using single whitefly transcriptomes in monitoring vector and viral interactions. Our method is applicable to any small organism where RNA quantity has limited transcriptome studies.

The first transcriptomes from field-collected individual whiteflies ( Bemisia tabaci, Hemiptera: Aleyrodidae): a case study of the endosymbiont composition

Background: Bemisia tabaci species ( B. tabaci), or whiteflies, are the world’s most devastating insect pests. They cause billions of dollars (US) of damage each year, and are leaving farmers in the developing world food insecure. Currently, all publically available transcriptome data for B. tabaci are generated from pooled samples, which can lead to high heterozygosity and skewed representation of the genetic diversity. The ability to extract enough RNA from a single whitefly has remained elusive due to their small size and technological limitations.

Methods: In this study, we optimised a single whitefly RNA extraction procedure, and sequenced the transcriptome of four individual adult Sub-Saharan Africa 1 (SSA1) B. tabaci. Transcriptome sequencing resulted in 39-42 million raw reads. De novo assembly of trimmed reads yielded between 65,000-162,000 Contigs across B. tabaci transcriptomes.

Results: Bayesian phylogenetic analysis of mitochondrion cytochrome I oxidase (mtCOI) grouped the four whiteflies within the SSA1 clade. BLASTn searches on the four transcriptomes identified five endosymbionts; the primary endosymbiont Portiera aleyrodidarum and four secondary endosymbionts: Arsenophonus, Wolbachia, Rickettsia, and Cardinium spp. that were predominant across all four SSA1 B. tabaci samples with prevalence levels of between 54.1 to 75%. Amino acid alignments of the NusG gene of P. aleyrodidarum for the SSA1 B. tabaci transcriptomes of samples WF2 and WF2b revealed an eleven amino acid residue deletion that was absent in samples WF1 and WF2a. Comparison of the protein structure of the NusG protein from P. aleyrodidarum in SSA1 with known NusG structures showed the deletion resulted in a shorter D loop.

Conclusions: The use of field-collected specimens means time and money will be saved in future studies using single whitefly transcriptomes in monitoring vector and viral interactions. Our method is applicable to any small organism where RNA quantity has limited transcriptome studies.

The first transcriptomes from field-collected individual whiteflies (Bemisia tabaci, Hemiptera: Aleyrodidae)

Background: Bemisia tabaci species (B. tabaci), or whiteflies, are the world’s most devastating insect pests. They cause billions of dollars (US) of damage each year, and are leaving farmers in the developing world food insecure. Currently, all publically available transcriptome data for B. tabaci are generated from pooled samples, which can lead to high heterozygosity and skewed representation of the genetic diversity. The ability to extract enough RNA from a single whitefly has remained elusive due to their small size and technological limitations. Methods: In this study, we optimised a single whitefly RNA extraction procedure, and sequenced the transcriptome of four individual adult Sub-Saharan Africa 1 (SSA1) B. tabaci. Transcriptome sequencing resulted in 39-42 million raw reads. De novo assembly of trimmed reads yielded between 65,000-162,000 Contigs across B. tabaci transcriptomes. Results: Bayesian phylogenetic analysis of mitochondrion cytochrome I oxidase (mtCOI) grouped the four whiteflies within the SSA1 clade. BLASTn searches on the four transcriptomes identified five endosymbionts; the primary endosymbiont Portiera aleyrodidarum and four secondary endosymbionts: Arsenophonus, Wolbachia, Rickettsia, and Cardinium spp. that were predominant across all four SSA1 B. tabaci samples with prevalence levels of between 54.1 to 75%. Amino acid alignments of the NusG gene of P. aleyrodidarum for the SSA1 B. tabaci transcriptomes of samples WF2 and WF2b revealed an eleven amino acid residue deletion that was absent in samples WF1 and WF2a. Comparison of the protein structure of the NusG protein from P. aleyrodidarum in SSA1 with known NusG structures showed the deletion resulted in a shorter D loop. Conclusions: The use of field-collected specimens means time and money will be saved in future studies using single whitefly transcriptomes in monitoring vector and viral interactions. Our method is applicable to any small organism where RNA quantity has limited transcriptome studies.

Plasma virome of cattle from forest region revealed diverse small circular ssDNA viral genomes

Background

Free-range cattle are common in the Northeast China area, which have close contact with farmers and may carry virus threatening to cattle and farmers.

Methods

Using viral metagenomics we analyzed the virome in plasma samples collected from 80 cattle from the forested region of Northeast China.

Results

The virome of cattle plasma is composed of the viruses belonging to the families including Parvoviridae, Papillomaviridae, Picobirnaviridae, and divergent viral genomes showing sequence similarity to circular Rep-encoding single stranded (CRESS) DNA viruses. Five such CRESS-DNA genomes were full characterized, with Rep sequences related to circovirus and gemycircularvirus. Three bovine parvoviruses belonging to two different genera were also characterized.

Conclusion

The virome in plasma samples of cattle from the forested region of Northeast China was revealed, which further characterized the diversity of viruses in cattle plasma.

Electronic supplementary material

The online version of this article (10.1186/s12985-018-0923-9) contains supplementary material, which is available to authorized users.

Discovery of the first maize-infecting mastrevirus in the Americas using a vector-enabled metagenomics approach

The genus Mastrevirus (family Geminiviridae) is composed of single-stranded DNA viruses that infect mono- and dicotyledonous plants and are transmitted by leafhoppers. In South America, there have been only two previous reports of mastreviruses, both identified in sweet potatoes (from Peru and Uruguay). As part of a general viral surveillance program, we used a vector-enabled metagenomics (VEM) approach and sampled leafhoppers (Dalbulus maidis) in Itumbiara (State of Goiás), Brazil. High-throughput sequencing of viral DNA purified from the leafhopper sample revealed mastrevirus-like contigs. Using a set of abutting primers, a 2746-nt circular genome was recovered. The circular genome has a typical mastrevirus genome organization and shares <63% pairwise identity with other mastrevirus isolates from around the world. Therefore, the new mastrevirus was tentatively named "maize striate mosaic virus". Seventeen maize leaf samples were collected in the same field as the leafhoppers, and ten samples were found to be positive for this mastrevirus. Furthermore, the ten genomes recovered from the maize samples share >99% pairwise identity with the one from the leafhopper. This is the first report of a maize-infecting mastrevirus in the Americas, the first identified in a non-vegetatively propagated mastrevirus host in South America, and the first mastrevirus to be identified in Brazil.

Discovery of two small circular ssDNA viruses associated with the whitefly Bemisia tabaci

The complete genome sequences of two novel small circular DNA viruses isolated from sweet-potato whiteflies collected in central-West (AdDF) and Southeast (AdO) regions of Brazil were determined by Next Generation Sequencing (NGS), and confirmed by cloning and Sanger sequencing. The genomes are 2,199 and 2,211 nt-long, respectively, encoding a putative coat protein (CP) and a replication-associated protein (Rep) and showing a genomic organization typical of viruses from the family Genomoviridae. Phylogenetic analysis with deduced amino acid sequences of Rep indicates that the virus from AdO is closely related to other members of the genus Gemycircularvirus, while the virus from AdDF is distantly related to other genomovirus. It was thus classified in a putative new genus, for which the name "Gemybolavirus" is proposed. These new genomoviruses are tentatively named "Bemisia associated gemybolavirus AdDF", and "Bemisia associated gemycircularvirus AdO".

Genome Sequence of Cauliflower Mosaic Virus Identified in Earwigs (Doru luteipes) through a Metagenomic Approach

Here we report the first complete genome sequence of a cauliflower mosaic virus from Brazil, obtained from the gut content of the predator earwig (Doru luteipes). This virus has a genome of 8,030 nucleotides (nt) and shares 97% genome-wide identity with an isolate from Argentina.

Genome Sequences of Beet curly top Iran virus, Oat dwarf virus, Turnip curly top virus, and Wheat dwarf virus Identified in Leafhoppers

Implementation of a vector-enabled metagenomics approach resulted in the identification of various geminiviruses. We identified the genome sequences of Beet curly top Iran virus, Turnip curly top viruses, Oat dwarf viruses, the first from Iran, and Wheat dwarf virus from leafhoppers feeding on beet, parsley, pumpkin, and turnip plants.

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.

Diverse circular replication-associated protein encoding viruses circulating in invertebrates within a lake ecosystem

Over the last five years next-generation sequencing has become a cost effective and efficient method for identifying known and unknown microorganisms. Access to this technique has dramatically changed the field of virology, enabling a wide range of environmental viral metagenome studies to be undertaken of organisms and environmental samples from polar to tropical regions. These studies have led to the discovery of hundreds of highly divergent single stranded DNA (ssDNA) virus-like sequences encoding replication-associated proteins. Yet, few studies have explored how viruses might be shared in an ecosystem through feeding relationships. Here we identify 169 circular molecules (160 CRESS DNA molecules, nine circular molecules) recovered from a New Zealand freshwater lake, that we have tentatively classified into 51 putatively novel species and five previously described species (DflaCV-3, -5, -6, -8, -10). The CRESS DNA viruses identified in this study were recovered from molluscs (Echyridella menzeisii, Musculium novaezelandiae, Potamopyrgus antipodarum and Physella acuta) and insect larvae (Procordulia grayi, Xanthocnemis zealandica, and Chironomus zealandicus) collected from Lake Sarah, as well as from the lake water and benthic sediments. Extensive diversity was observed across most CRESS DNA molecules recovered. The putative capsid protein of one viral species was found to be most similar to those of members of the Tombusviridae family, thus expanding the number of known RNA-DNA hybrid viruses in nature. We noted a strong association between the CRESS DNA viruses and circular molecules identified in the water and browser organisms (C. zealandicus, P. antipodarum and P. acuta), and between water sediments and undefended prey species (C. zealandicus). However, we were unable to find any significant correlation of viral assemblages to the potential feeding relationships of the host aquatic invertebrates.

Begomovirus-Associated Satellite DNA Diversity Captured Through Vector-Enabled Metagenomic (VEM) Surveys Using Whiteflies (Aleyrodidae)

Monopartite begomoviruses (Geminiviridae), which are whitefly-transmitted single-stranded DNA viruses known for causing devastating crop diseases, are often associated with satellite DNAs. Since begomovirus acquisition or exchange of satellite DNAs may lead to adaptation to new plant hosts and emergence of new disease complexes, it is important to investigate the diversity and distribution of these molecules. This study reports begomovirus-associated satellite DNAs identified during a vector-enabled metagenomic (VEM) survey of begomoviruses using whiteflies collected in various locations (California (USA), Guatemala, Israel, Puerto Rico, and Spain). Protein-encoding satellite DNAs, including alphasatellites and betasatellites, were identified in Israel, Puerto Rico, and Guatemala. Novel alphasatellites were detected in samples from Guatemala and Puerto Rico, resulting in the description of a phylogenetic clade (DNA-3-type alphasatellites) dominated by New World sequences. In addition, a diversity of small (~640-750 nucleotides) satellite DNAs similar to satellites associated with begomoviruses infecting Ipomoea spp. were detected in Puerto Rico and Spain. A third class of satellite molecules, named gammasatellites, is proposed to encompass the increasing number of reported small (<1 kilobase), non-coding begomovirus-associated satellite DNAs. This VEM-based survey indicates that, although recently recovered begomovirus genomes are variations of known genetic themes, satellite DNAs hold unexplored genetic diversity.