Evidence That Lokern Virus
Total Page:16
File Type:pdf, Size:1020Kb
Tangudu et al. Virology Journal (2018) 15:122 https://doi.org/10.1186/s12985-018-1031-6 RESEARCH Open Access Evidence that Lokern virus (family Peribunyaviridae) is a reassortant that acquired its small and large genome segments from Main Drain virus and its medium genome segment from an undiscovered virus Chandra S. Tangudu1, Jermilia Charles1 and Bradley J. Blitvich1,2* Abstract Background: Lokern virus (LOKV) is a poorly characterized arthropod-borne virus belonging to the genus Orthobunyavirus (family Peribunyaviridae). All viruses in this genus have tripartite, single-stranded, negative-sense RNA genomes, and the three RNA segments are designated as small, (S), medium (M) and large (L). A 559 nt. region of the M RNA segment of LOKV has been sequenced and there are no sequence data available for its S or L RNA segments. The purpose of this study was to sequence the genome of LOKV. Methods: The genome of LOKV was fully sequenced by unbiased high-throughput sequencing, 5′ and 3′ rapid amplification of cDNA ends, reverse transcription-polymerase chain reaction and Sanger sequencing. Results: The S and L RNA segments of LOKV consist of 952 and 6864 nt. respectively and both have 99.0% nucleotide identity with the corresponding regions of Main Drain virus (MDV). In contrast, the 4450-nt. M RNA segment has only 59.0% nucleotide identity with the corresponding region of MDV and no more than 72.7% nucleotide identity with all other M RNA segment sequences in the Genbank database. Phylogenetic data support these findings. Conclusions: This study provides evidence that LOKV is a natural reassortant that acquired its S and L RNA segments from MDV and its M RNA segment from an undiscovered, and possibly extinct, virus. The availability of complete genome sequence data facilitates the accurate detection, identification and diagnosis of viruses and viral infections, and this is especially true for viruses with segmented genomes because it can be difficult or even impossible to differentiate between reassortants and their precursors when incomplete sequence data are available. Keywords: Lokern virus, Main drain virus, Orthobunyavirus, Peribunyaviridae, Reassortant, Genome sequence, High- throughput sequencing * Correspondence: [email protected] 1Department of Veterinary Microbiology and Preventive Medicine, College of Veterinary Medicine, Iowa State University, Ames, IA, USA 22116 Veterinary Medicine, Iowa State University, Ames, Iowa 50011, USA © The Author(s). 2018 Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated. Tangudu et al. Virology Journal (2018) 15:122 Page 2 of 8 Background has been sequenced (Genbank Accession no. AY593736.1) Lokern virus (LOKV) is a poorly characterized virus in the and no sequence data are available for the S or L RNA genus Orthobunyavirus (family Peribunyaviridae). The segments. The objective of this study is to fully sequence genus has previously been separated into serogroups, with the genome of LOKV and to determine its genetic and LOKV classified in the Bunyamwera (BUN) serogroup [1]. phylogenetic relationships with other BUN serogroup Other viruses in this serogroup include Batai virus viruses. (BATV), Bunyamwera virus (BUNV), Cache Valley virus (CVV), Maguari virus (MAGV), Ngari virus (NRIV), Methods Northway virus (NORV), and Tensaw virus (TENV). Virus LOKV was originally isolated from Culex tarsalis mosqui- LOKV (isolate FMS 4332) was obtained from the World toes in California in 1962 [2]. Additional isolations have Reference Center for Emerging Viruses and Arboviruses been made from arthropods of various species as well as at the University of Texas Medical Branch in Galveston, rabbits and hares in the western half of the United States, TX. The virus had been passaged 13 times in suckling with most isolates recovered from Culicoides spp. midges mouse brains prior to receipt and underwent an additional [2–5]. Serologic data indicate that the geographic range of passage in African Green Monkey kidney (Vero) cells at LOKV also includes the eastern half of the United States Iowa State University. [6]. LOKV is not a recognized pathogen of humans or ver- tebrate animals but its ability to cause disease has not High-throughput sequencing been widely investigated. Total RNA was extracted from subconfluent monolayers All viruses in the genus Orthobunyavirus possess a tri- of LOKV-infected Vero cells using Trizol Reagent (Ther- partite, single-stranded, negative-sense RNA genome [7]. moFisher, Carlsbad, CA), fragmented using RNase III The three RNA segments are designated as small, (S), (New England BioLabs, Ipswich, MA) and assessed for medium (M) and large (L) and are approximately 1.0, quality using an Agilent 2100 Bioanalyzer (Agilent, Santa 4.5 and 6.9 kb in length, respectively. Complementary Clara, CA). Libraries were constructed using the Ion Total and highly conserved nucleotide sequences are located RNA-Seq Kit v2 (ThermoFisher) and assessed for quality at the 5′ and 3′ termini of each RNA segment. Base and analyzed at the Genomic Technologies Facility at pairing of these sequences results in the formation of Iowa State University using an Ion Proton Sequencer panhandle structures assumed to be required for the ini- (ThermoFisher). Ion-Torrent reads were mapped to the tiation of genome transcription and replication [8, 9]. genome of Chlorocebus sabaeus using Bowtie 2 in order The S RNA segment codes the nucleocapsid (N) protein to identify host cell sequences [17]. Unmapped reads were and a nonstructural protein NSS in overlapping reading analyzed using the sortMeRNA program to remove frames [10]. The M RNA segment codes for a polypro- rRNA-related reads [18]. All remaining reads with Phred tein that is post-translationally processed to generate the values ≥ 33 were subjected to de novo SPAdes (ver 3.5.0) N- and C-terminal glycoproteins, Gn and Gc, and a non- assembly [19]. Contigs were mapped to a reference ortho- structural protein NSM [11]. The L RNA segment codes bunyavirus genome (MDV for S and L RNA segment for the L protein which functions as the RNA-dependent sequences; CVV for M RNA segment sequences) using RNA polymerase [12]. LASTZ [20]. Alignment files were manually verified on The exchange of RNA segments (genome reassortment) Tablet [21]. between closely-related orthobunyaviruses can occur when a vector or host is co-infected with two or more vi- Reverse transcription-polymerase chain reaction and sanger ruses [13, 14]. Reassortment is considered to be a major sequencing means of evolution for these viruses because the resulting Reverse transcription-polymerase chain reaction (RT- progeny can possess new and advantageous genetic and PCR) and Sanger sequencing were performed to verify phenotypic traits. Reassortment sometimes generates vi- LOKV sequences generated by high-throughput sequen- ruses that are more pathogenic than their donor viruses. cing and to close gaps between contigs. Complementary One example is NRIV which acquired its S and L RNA DNAs were generated using Superscript III reverse tran- segments from BUNV and its M RNA segment from scriptase (ThermoFisher) and PCRs were performed using BATV [15, 16]. NRIV has been associated with large out- high-fidelity Taq polymerase (ThermoFisher) in accord- breaks of hemorrhagic fever in humans in East Africa ance to the manufacturer’s instructions. Primers were de- while the most serious disease manifestations associated signed from the sequences generated by high-throughput with its precursor viruses are low-grade fever and myalgia. sequencing and are available upon request. RT-PCR prod- The genome of LOKV has not been fully sequenced. A ucts were purified using the purelink gel extraction kit 559 nt. region of the M RNA segment that spans part of (ThermoFisher) and sequenced using a 3730 × 1 DNA se- the 5′ untranslated region (UTR) and Gn coding region quencer (Applied Biosystems, Foster City, CA). Tangudu et al. Virology Journal (2018) 15:122 Page 3 of 8 5′ and 3′ rapid amplification of cDNA ends Table 1 Genetic relatedness of the small, medium and large The 5′ and 3′ ends of each RNA segment were identi- genome segments of Lokern virus and its closest known relatives fied using 5′ and 3′ rapid amplification of cDNA ends, Virus Lokern virus respectively. Briefly, 20 pmol of pre-adenylated DNA % nucleotide identity % amino acid identity ′ a adaptor (5 -rApp/TGGAATTCTCGGGTGCCAAGGT/ Small Medium Large N NSS ML ′ μ ddC-3 ) was ligated to 1 g of viral genomic and Abbey Lake virus 70.9 59.8 71.4 73.4 64.0 55.0 76.5 anti-genomic RNAs using 200 U of T4 RNA ligase 2 (New Batai virus 79.5 68.7 73.1 88.4 90.1 71.5 80.7 England BioLabs) and incubated for 1 h at 25 °C. Ligated b–– –– products were recovered by ethanol precipitation and Birao virus 78.3 87.6 85.0 complementary cDNAs were created using Superscript III Bozo virus 76.3 ––79.0 78.2 –– reverse transcriptase and an adapter-specific primer. PCRs Bunyamwera virus 80.0 63.4 73.5 89.7 86.2 62.4 81.1 were performed using adapter- and gene-specific primers, Cache Valley virus 83.5 72.5 76.1 89.3 94.1 78.4 85.5 and RT-PCR products were purified and subjected to Cholul virus 83.3 58.6 – 89.3 94.1 54.4 – Sanger sequencing. Fort Sherman virus 81.8 72.7 76.3 91.0 93.1 77.9 86.0 – – Nucleotide and amino acid sequence alignments Germiston virus 71.1 59.3 73.8 67.0 54.7 Nucleotide and predicted amino acid sequences of LOKV Ilesha virus 79.6 67.4 73.4 86.7 85.2 67.0 80.2 were aligned with other sequences in the Genbank data- Main Drain virus 99.0 59.0 99.0 99.6 100.0 53.9 99.0 base by application of BLASTn and BLASTp, respectively.