Novel virus discovery and genome reconstruction from field RNA samples reveals highly divergent viruses in dipteran hosts.

We investigated whether small RNA (sRNA) sequenced from field-collected mosquitoes and chironomids (Diptera) can be used as a proxy signature of viral prevalence within a range of species and viral groups, using sRNAs sequenced from wild-caught specimens, to inform total RNA deep sequencing of samples of particular interest. Using this strategy, we sequenced from adult Anopheles maculipennis s.l. mosquitoes the apparently nearly complete genome of one previously undescribed virus related to chronic bee paralysis virus, and, from a pool of Ochlerotatus caspius and Oc. detritus mosquitoes, a nearly complete entomobirnavirus genome. We also reconstructed long sequences (1503-6557 nt) related to at least nine other viruses. Crucially, several of the sequences detected were reconstructed from host organisms highly divergent from those in which related viruses have been previously isolated or discovered. It is clear that viral transmission and maintenance cycles in nature are likely to be significantly more complex and taxonomically diverse than previously expected

Novel virus discovery and genome reconstruction from field RNA samples reveals highly divergent viruses in dipteran hosts.

We investigated whether small RNA (sRNA) sequenced from field-collected mosquitoes and chironomids (Diptera) can be used as a proxy signature of viral prevalence within a range of species and viral groups, using sRNAs sequenced from wild-caught specimens, to inform total RNA deep sequencing of samples of particular interest. Using this strategy, we sequenced from adult Anopheles maculipennis s.l. mosquitoes the apparently nearly complete genome of one previously undescribed virus related to chronic bee paralysis virus, and, from a pool of Ochlerotatus caspius and Oc. detritus mosquitoes, a nearly complete entomobirnavirus genome. We also reconstructed long sequences (1503-6557 nt) related to at least nine other viruses. Crucially, several of the sequences detected were reconstructed from host organisms highly divergent from those in which related viruses have been previously isolated or discovered. It is clear that viral transmission and maintenance cycles in nature are likely to be significantly more complex and taxonomically diverse than previously expected

Codon usage statistics for the RdRp ORFs in the narnavirus-like contigs KF298275 and KF298276.

<p>In-frame forward read-direction stop codons (red) are necessarily absent. Reverse complements of in-frame but reverse read-direction stop codons are highlighted in orange; a single UUA codon corresponds to the UAA stop codon of the >1000-codon reverse-read-direction ORF in KF298276.</p

Common nucleotide polymorphisms observed in AACV and CAZV assemblies.

<p>Only nucleotide variations with a frequency in the mapped reads of at least 20% of the most common nucleotide present at that position are shown. The nucleotide present in the chosen reference assembly is shown first.</p

Neighbour-joining (NJ) phylogenetic tree for the orbivirus-like KF298266 and KF298273 and related sequences.

<p>Nodes with <80% bootstrap support have been collapsed.</p

Analysis of CBPV and AACV sequences.

<p>(A) Map of the CBPV RNA1 genome segment. The region covered by contig KF298264 (AACV) is indicated by the orange bar. The read coverage density is indicated in green. (B) Analysis of variability at synonymous sites in an alignment of the currently available full-length CBPV sequences (EU122229 and EU122231) and AACV (KF298264). Shown are the degree of variability at synonymous sites in a 75-codon sliding window, relative to the average in the ORF1-ORF3 frameshift fusion (obs/exp), and the corresponding statistical significance (<i>p</i>-value). (C) Positions of stop codons in the three forward reading frames in the three sequences (KF298264 - top row of triangles in each panel; CBPV EU122231 and EU122229 - bottom two rows of triangles in each panel). (D, E and F) Corresponding figures for RNA2 (KF298265 - top row of triangles in each panel; CBPV EU122232 and EU122230 - bottom two rows of triangles in each panel). Note that AACV RNA2 lacks a homolog of the CBPV ORF1, and has a shorter 3' UTR than CBPV, as indicated by gaps in the orange bar.</p

Alignments of entomobirnavirus 5' and 3' UTR sequences, indicating that the CAZV sequences are nearly complete.

<p>Initiation and termination codons are highlighted in green and red respectively. Note that both entomobirnavirus segments have upstream AUG codons, although it is not known whether these are utilised. ESV - JN589003/JN589002, CYV - JQ659254/JQ659255, MXV - JX403941/JX403942, CAZV - KF298271/KF298272.</p

Bayesian maximum likelihood (ML) phylogenetic tree for the flavivirus-like KF298267 and related sequences.

<p>The tree is midpoint-rooted and, for clarity, only posterior probability values <u>></u>80% are shown.</p