The Apis Mellifera Filamentous Virus Genome
Total Page:16
File Type:pdf, Size:1020Kb
View metadata, citation and similar papers at core.ac.uk brought to you by CORE provided by Bern Open Repository and Information System (BORIS) Viruses 2015, 7, 3798-3815; doi:10.3390/v7072798 OPEN ACCESS viruses ISSN 1999-4915 www.mdpi.com/journal/viruses Article The Apis mellifera Filamentous Virus Genome Laurent Gauthier 1;:;*, Scott Cornman 2;:;;, Ulrike Hartmann 1, François Cousserans 1, Jay D. Evans 2, Joachim R. de Miranda 3 and Peter Neumann 1;4 1 Agroscope, Swiss Bee Research Centre, Schwarzenburgstrasse 161, CH-3003 Bern, Switzerland; E-Mails: [email protected] (U.H.); [email protected] (F.C.); [email protected] (P.N.) 2 Bee Research Laboratory, Beltsville, MD 20705, USA; E-Mails: [email protected] (S.C.); [email protected] (J.D.E.) 3 Department of Ecology, Swedish University of Agricultural Sciences, Uppsala 750 07, Sweden; E-Mail: [email protected] 4 Institute of Bee Health, Vetsuisse Faculty, University of Bern, CH-3001 Bern, Switzerland : These authors contributed equally to this work. ; Current address: Fort Collins Science Center, USGS, Fort Collins, CO 80526, USA. * Author to whom correspondence should be addressed; E-Mail: [email protected]; Tel.: +41-584-622-155; Fax: +41-584-638-227. Academic Editors: Elke Genersch and Sebastian Gisder Received: 31 March 2015 / Accepted: 2 July 2015 / Published: 9 July 2015 | downloaded: 13.3.2017 Abstract: A complete reference genome of the Apis mellifera Filamentous virus (AmFV) was determined using Illumina Hiseq sequencing. The AmFV genome is a double stranded DNA molecule of approximately 498,500 nucleotides with a GC content of 50.8%. It encompasses 247 non-overlapping open reading frames (ORFs), equally distributed on both strands, which cover 65% of the genome. While most of the ORFs lacked threshold sequence alignments to reference protein databases, twenty-eight were found to display significant homologies with proteins present in other large double stranded DNA viruses. Remarkably, 13 ORFs had strong similarity with typical baculovirus domains such as PIFs (per os infectivity factor genes: pif-1, pif-2, pif-3 and p74) and BRO (Baculovirus Repeated Open http://boris.unibe.ch/82558/ Reading Frame). The putative AmFV DNA polymerase is of type B, but is only distantly related to those of the baculoviruses. The ORFs encoding proteins involved in nucleotide source: metabolism had the highest percent identity to viral proteins in GenBank. Other notable Viruses 2015, 7 3799 features include the presence of several collagen-like, chitin-binding, kinesin and pacifastin domains. Due to the large size of the AmFV genome and the inconsistent affiliation with other large double stranded DNA virus families infecting invertebrates, AmFV may belong to a new virus family. Keywords: Apis mellifera; filamentous virus; genome annotation 1. Introduction Apis mellifera filamentous virus (AmFV) was originally described as a rickettsia disease of honeybees (Apis mellifera)[1], but was later characterized as a large enveloped DNA virus [2]. It derives its name from the long (3150 ˆ 40 nm) filamentous nucleoprotein that folds in three superimposed figure-8 loops (Figure1) into a 450 nm ˆ 170 nm rod-shaped virion [2,3]. The nucleoprotein contains a central core of double stranded DNA (dsDNA) wrapped by two major nucleoproteins while the tri-laminate virion envelope contains lipids, two major proteins and several minor proteins [4]. One major diagnostic feature of AmFV infection is that the hemolymph of severely infected adult honeybees becomes milky-white due to cellular degradation and the large number of virions present, which is also a primary feature of insect ascoviruses infecting Lepidopteran hosts [5]. A. mellifera filamentous virus epizooties were identified in various parts of the world through electron microscopic examination of bee hemolymph samples [6]. AmFV was shown to multiply in different bee tissues, principally in the fat body and ovarian tissues of adult bees [7]. Highly infected bees display signs of weakness and are usually gathered at the hive entrance. Although it may occasionally induce such overt, colony-level symptoms, AmFV is usually considered a weakly pathogenic virus of honeybees, with low impact on host lifespan [7,8]. Experimental infections can be established by feeding adult bees AmFV particles, but only when co-infected with Nosema apis spores [9]. Historically, there has been a strong association between Nosema sp. and AmFV infections [7]. AmFV has occasionally been associated with winter honeybee colony mortality [3], but its annual peak incidence is mostly in spring, coinciding with the peak incidence of N. apis [9]. Filamentous viruses with similar morphology and tissue distribution have also been observed in a number of parasitic wasps, where they may be associated with super-parasitism behavior [10]. Here we present the complete nucleotide sequence and annotation of the AmFV genome, as well as insights into the distribution of this virus in honeybee colonies. Although some coding sequences are homologous with other large dsDNA viruses of invertebrates, such as baculoviruses, the size of AmFV genome, the proportion of putatively non-coding sequences and the large number of open reading frames (ORFs) lacking database matches all suggest a very distinct virus that merits classification in a new family of insect viruses. 2. Materials and Methods 2.1. AmFV DNA Isolation, Sequencing, and Assembly Viral particles were purified from six infected honeybee workers, collected from a single colony in the spring of 2009 in Bern, Switzerland, that were displaying the characteristic “milky hemolymph” Viruses 2015, 7 3 2. Materials and Methods 2.1. AmFV DNA Isolation, Sequencing, and Assembly Viruses 2015, 7 3800 Viral particles were purified from six infected honeybee workers, collected from a single colony in the spring of 2009 in Bern, Switzerland, that were displaying the characteristic “milky hemolymph” AmFV clinical signs. After homogenization in 1.5 mLmL of 50 mM Tris pH 7.5 buffer, viralviral particlesparticles werewere r separated byby centrifugationcentrifugation onon aa 10%–60%10%–60% NicodenzNicodenz® (Sigma-Aldrich,(Sigma-Aldrich, St Louis,Louis, MI,MI, USA)USA) gradient.gradient. Two visiblevisible bandsbands werewere extractedextracted fromfrom thethe gradientgradient andand analyzedanalyzed byby transmission electron microscopy (TEM).(TEM). The presence presence of of AmFV AmFV enveloped enveloped virions virions (Figure (Figure 1)1 )was was confirmed confirmed in in the the lower lower band band while while the theupper upper band band contained contained mostly mostly degraded degraded nucleoproteins. nucleoproteins. Total DNA Total was DNA purified was purified from the from lower the band lower by bandphenol by-chloroform phenol-chloroform extraction extraction and precipitation and precipitation with ethanol. with The ethanol. DNA Thewas resuspended DNA was resuspended in 5 mM Tris in 5pH mM 7.5 Trisand pHquantified 7.5 and using quantified the Qubit using® system the Qubit (Lifer systemTechnologies (Life Technologies,, Carlsbad, CA, Carlsbad, USA) One CA, hundred USA) Onenanograms hundred of nanograms DNA were of sent DNA to wereFasteris sent Life to Fasteris Science Life Co. Science(Geneva, Co. Switzerland) (Geneva, Switzerland) where a paired where-end a paired-endlibrary was libraryprepared was using prepared the TruSeq using® the Nano TruSeq DNAr Nano LT kit DNA (Illumina LT kit, San (Illumina, Diego, SanCA, Diego, USA) CA,with USA)bead- withbased bead-based size selection size selectionof 550 bp of fragments. 550 bp fragments. The libr Theary was library sequenced was sequenced on an onIllumina an Illumina HiSeq HiSeq 2500 2500Analyzer Analyzer,, producing producing 100-bp 100-bp paired pairedreads that reads were that quality were filtered quality and filtered adapter and screened adapter screenedby the sequence by the sequenceprovider and provider assembled and assembled de novo atde a novorangeat of a kmers range ofusing kmers the usingVELVET the VELVET program program[11]. [11]. Alignments from the scaffold datasetdataset werewere producedproduced usingusing thethe SeqMan (Lasergene 7.0;7.0; DnastarDnastar Inc.,Inc., Madison, WI, USA)USA) software. Gaps and ambiguous regions were amplifiedamplified by PCR using a high fidelityfidelity polymerase (HIFI,(HIFI, Kapa-Biosystems,Kapa-Biosystems, Wilmington,Wilmington, MA,MA, USA)USA) andand appropriateappropriate primers.primers. The resultingresulting amplicons werewere sequencedsequenced onon bothboth strandsstrands usingusing thethe SangerSanger sequencingsequencing method.method. JellyfishJellyfish [[12]12] waswas usedused to to evaluate evaluate the the heterogeneity heterogeneity of of the the isolated isolated sample sample by by plotting plotting the the coverage coverage of allof kmersall kmers of length of length 59 occurring 59 occurring in the rawin th reads.e raw Batchesreads. Batches of kmers of were kmers assembled were assembled with CLC Genomicswith CLC v.Genomics 7.1 (Qiagen, v. 7.1 Valencia, (Qiagen CA,, Valencia, USA) using CA, USA a kmer) using of 59 a and kmer searched of 59 and by BLASTX searched againstby BLASTX the UniRef90 against databasethe UniRef90 [13]. database [13]. Figure 1.1. MorphologyMorphology of of AmFV AmFV virions virions:: ( (AA,,BB)) Electron Electron micrographs micrographs of “milky” beebee hemolymph containing characteristic AmFV nucleoproteins (np) and enveloped virionsvirions (v).(v). Approximate