Skip to main content

Identification and genetic analysis of Kadipiro virus isolated in Shandong province, China

Abstract

Background

Kadipiro virus (KDV) belongs to the Reoviridae family, which consists of segmented, non-enveloped, double-stranded RNA viruses. It has previously been isolated from Culex, Anopheles, Armigeres and Aedes mosquitoes in Indonesia and China. Here, we describe the isolation and characterization of SDKL1625 from Anopheles sinensis mosquitoes in Shandong province, China.

Methods

In this study, we isolated Kadipiro virus in Aedes albopictus C6/36 cell culture and the complete genome sequencing was made by next generation sequencing.

Results

We isolated and characterized a Kadipiro virus from Anopheles sinensis mosquitoes in 2016 in Shandong province, China. Nucleotide and amino acid homology analysis of SDKL1625 showed higher levels of sequence identity with QTM27331 (Odonata, China, 2016) than with JKT-7075 (Culex fuscocephalus, Indonesia, 1981). The SDKL1625 has 86–97% amino acid identity with the JKT-7075, 88–99% amino acid identity with the QTM27331. Among the 12 fragments, VP1, VP2, VP4, VP6, VP7, VP9 and VP12 showed high amino acid identity (> 90%) and VP5 showed the lowest identity (86% and 88%).

Conclusions

This is the first identification of KDV from mosquito in China. Virus morphology and genome organization were also determined, which will further enrich our understanding of the molecular biological characteristics of KDV and seadornaviruses.

Background

Kadipiro virus (KDV) is a species of the new 12-segmented RNA virus genus Seadornavirus, which was assigned to the Reoviridae family in the eighth and ninth report of the International Committee on the Taxonomy of Viruses [1, 2]. The Seadornavirus genus includes three species, Banna virus (BAV), Kadipiro virus (KDV) and Liao ning virus (LNV). New species have been isolated in recent years, including the Balaton [3] and Mangshi viruses [4].

BAV is the prototype species of the genus Seadornavirus and was first isolated in 1987 from patients with encephalitis in Xishuangbanna prefecture, Yunnan province, China [5]. To date, BAV has been isolated from mosquito species, humans with encephalitis, pigs, cattle and ticks in Southeast Asia, particularly China and Indonesia [6, 7]. BAV is an emerging pathogen that causes human viral encephalitis, exhibits rich genetic diversity and great differentiation potential, and does not exhibit a species barrier [8]. The discovery of LNV, a new virus strain obtained from mosquitoes in China, enriched the variety of the Seadornavirus genus. LNV was able to replicate not only in mosquito cell lines, but also in mammalian cell lines, and thus may be able to infect mammals [9]. Phylogenetic analysis show LNV might be an emerging virus that evolved rapidly and is widely distributed in the northern part of China [10, 11].

KDV was first isolated from Culex fuscocephalus in Java, Indonesia in 1981 by replication in the C6/36 cell line, grouped into the genus Coltivirus within the family Reoviridae as the enough genetic information could not be obtained to permit exact classification [7, 12]. Since then, KDV has also been isolated from Culex tritaeniorhynchus, Anopheles sinensis and Armigeres subalbatus in northwestern Yunnan province, China [13]. To date, only two full genome KDV sequences have been recorded in GenBank [14, 15]. Recently, a viral metagenomics study showed that a plasma sample of a febrile adult, who enrolled in a study of acute human immunodeficiency virus type 1 (HIV-1) infection in coastal Kenya, exhibited translated protein matches to KDV virus. This discovery indicated that the tropism of KDV includes humans [16].

Here, we have isolated a strain of KDV from Anopheles sinensis mosquitoes in Kenli county, Shandong province, China. We describe the isolation and characterization, as well as genomic analyses, of this virus and its phylogenetic relationship with members of Seadornavirus, which will further enrich our understanding of the molecular biological characteristics of KDV and seadornaviruses.

Methods

Mosquito collection

Mosquito samples were collected on September 7, 2016 in residential regions of the village of Dongsui (river estuary town) (37°38′31.94”N, 118°51′41.69”E; altitude, 2 m) in Kenli county, Dongying, in northeast Shandong province, China.

Mosquitoes were collected from a pigsty at night using ultraviolet light traps (12 V; 300 mA; Wuhan Lucky Star Environmental Protection Tech Co., Ltd., Hubei, China). Traps were set from 7:00 pm to 7:00 am, which spanned the period from sunset to sunrise. The captured mosquitoes were frozen at − 20 °C for at least 40 min and identified based on morphology. The adult mosquitoes captured were identified according to morphological keys in Lu Baolin et al. [17, 18]. The mosquito species were selected on a chilled plate and only female mosquitoes were collected. Approximately 50–100 mosquitoes of the same species were pooled in cryogenic vials and stored in liquid nitrogen until further use [19].

Virus isolation

Mosquitoes were removed from the liquid nitrogen and immediately homogenized and centrifuged as previously reported [20]. The supernatants were inoculated in a monolayer of C6/36 or BHK-21 cells in 24-well plates and incubated at 28 °C and 37 °C, respectively. Cells were observed for cytopathic effects (CPEs) every 8 h for 6–7 days after incubation for 24 h. A specimen was regarded as a positive isolate if it caused CPEs in three successive cell passages. Infected cell supernatants were stored at − 80 °C for further analysis.

Virus electron microscopy

For negative staining, cell culture supernatants were absorbed on Formvar and carbon-coated grids for 1 min and stained with 1% (w/v) phosphotungstic acid (pH 6.8) for 1 min. Grids were air dried and observed under a Tecnai 12 transmission electron microscope (FEI, Eindhoven, Netherlands).

For ultrathin sections, cell pellets were fixed in a solution of 2% formaldehyde and 2.5% glutaraldehyde, post-fixed in 1% osmium tetroxide, dehydrated in an ethanol gradient, embedded in epoxy resin, and polymerized at 60 °C for 24 h. Ultrathin sections (80 nm) were obtained from the resin blocks, mounted on copper grids, and stained with uranyl acetate and lead citrate. Finally, the ultrathin sections were observed using transmission electron microscopy.

Polyacrylamide gel electrophoresis (PAGE)

RNA was subject to PAGE on a discontinuous system with 5% acrylamide (acrylamide/bisacrylamide 29:1; Bio-Rad Laboratories, Hercules, CA, USA) stacking gel (with Tris-HCl buffer, pH 6.8) and 8% acrylamide resolving gel (with Tris-HCl buffer, pH 8.8). Following electrophoresis, the gel was stained with a Protein Silver Stain Kit (Leagene, Beijing, China) according to the manufacturer’s instructions and observed under white light. The KDV YN0557 was selected as positive control [13].

Extraction of nucleic acids and RT-PCR assays for detection of virus RNA

Viral RNA was extracted from specimens with positive CPEs using the QIAamp Viral RNA Mini Kit (Qiagen, Valencia, CA, USA). Complementary DNA (cDNA) synthesis was performed using Ready-to-Go RT-PCR Beads (Amersham Biosciences Co., Piscataway, NJ, USA). Viral cDNA was subjected to PCR amplification with primers targeting the 12th segment of KDV [13] (KDV-12-F: 5’-GACGCTTTGAGATTATCTCGAC-3′, KDV-12-R: 5’-GCTCAATCGCATTCTCACC-3′) and GoTaqGreen Master Mix. PCR products were sequenced using BigDye Terminator chemistry (Applied Biosystems, Foster City, CA, USA).

Deep sequencing

The library was prepared using NEBNext Ultra RNA Library Prep Kit for Illumina (New England Biolabs, Hitchin, UK). The library was quantified on Qubit 3.0 fluorometer and Agilent Bioanalyzer 2100. Paired-end sequencing (2 × 250 base pairs) was performed on an Illumina MiSeq v2 platform. For analysis, reads were filtered based on their length and mean quality values (Q30). KDV sequence reads were mapped against the published sequence. Mapping was performed using the Burrows-Wheeler aligner-MEM program [21]. The complete nucleotide sequences of the viruses reported in this study have been submitted to the GenBank database under accession numbers MG590140-MG590151.

Phylogenetic analysis

Phylogenetic analysis was conducted based on the protein sequence deduced from the conserved RNA-dependent RNA polymerase gene of selected BAV, LNV and KDV. A multiple sequence alignment matrix was created using ClustalW software (http://clustalw.ddbj.nig.ac.jp/) with default settings. The aligned matrix data were confirmed manually, and the amino acid sequences were analyzed using the maximum-likelihood method. The statistical significance of the resulting tree was evaluated using a bootstrap test with 1000 replicates.

Results

Mosquito collection

On September 7, 2016, a total of 3000 mosquitoes representing two species, Culex tritaeniorhynchus (2300) and Anopheles sinensis (700) were collected in Kenli county, Shandong province, China.

Virus isolation and identification

The SDKL1625 of KDV was isolated from one (Anopheles sinensis) of 33 mosquito pools captured in Kenli county, Shandong province, China in 2016. The virus was shown to cause CPEs in C6/36 cells beginning at 72 h after inoculation in first-blind passage. The main characteristic of the CPEs was cell shedding, along with cell shrinking, or floating with eventual detachment from the growth surface (Fig. 1a). No CPEs were observed in BHK-21 cells during three blind passages.

Fig. 1
figure 1

Virus isolation and identification. a Cytopathic effects (CPEs) caused by Kadipiro virus (KDV; SDKL1625) in C6/36 cells 72 h after inoculation. The non-infected control and infected cells are shown. b Electrophoretic migration patterns of the dsRNA of Kadipiro virus as determined by polyacrylamide gel electrophoresis. Lane 1: Positive control; lane 2: SDKL1625. c Electron micrograph of KDV shows virus particles with a diameter of approximately 70 nm. Arrows indicate KDV particles in infected cell culture supernatant; ultrathin sections show large, electron-dense viral inclusion body (arrow) within the cytoplasm composed of a large number of KDVs in infected C6/36 cells

RT-PCR assays were performed to detect the presence of the virus. Reactions with the KDV primers were positive with a product length of 370 bp. The YN0557 of KDV was isolated from mosquitoes in Yunnan province, China in 2009, and the RNA of this virus was used as a positive control in agarose gel electrophoresis [13]. After sequence analysis using the National Center for Biotechnology Information BLAST program (https://blast.ncbi.nlm.nih.gov/Blast.cgi), we preliminarily identified SDKL1625 as KDV.

The double-stranded RNA of the SDKL1625, which was purified and extracted from infected cell cultures, was analyzed by PAGE. PAGE analysis revealed that SDKL1625 exhibited the same pattern of bands as YN0557. The banding pattern consisted of 12 segments, forming a “6–5-1” pattern, which is an important characteristic of KDV (Fig. 1b).

Virus electron microscopy

The KDV virus particles from the supernatant of infected C6/36 cells appeared spherical with a diameter of approximately 70 nm and the virion surface was covered with projections, presenting a clear double-layered capsid structure. Ultrathin sections showed electron-dense virus particles within the cytoplasm of KDV-infected C6/36 cells aggregated as viral inclusion bodies, which are thought to be the main site of replication, similar to the Banna virus (Fig. 1c) [22].

Genome characterization

The complete genome sequences of SDKL1625 were characterized by next generation sequencing. A maximum-likelihood phylogenetic tree was inferred from the aligned amino acid sequences of segment 1 of the SDKL1625 together with two KDVs and 15 other representatives of the genus Seadornavirus (Fig. 2). In the phylogenetic tree, the SDKL1625 sequences clustered with the KDV sequences. Nucleotide and amino acid sequences from genome segments 1–12 of SDKL1625 were compared with homologous KDV sequences of JKT-7075 and QTM27331 (Table 1). Nucleotide sequence identities between SDKL1625 and JKT-7075 ranged from 81% (VP1 and VP3) to 91% (VP7 and VP9), and ranged from 86% (VP5) to 99% (VP1, VP3, VP6, VP7 and VP9) between SDKL1625 and QTM27331. The amino acid identities of SDKL1625 with the JKT-7075 and QTM27331 ranged from 86% (VP5) to 97% (VP2) and from 88% (VP5) to 99% (VP1, VP3, VP6, VP7 and VP9), respectively.

Fig. 2
figure 2

Molecular phylogenetic analyzed by the maximum-likelihood method. A maximum likelihood phylogenetic tree of Banna virus (BAV), KDV and Liao ning virus (LNV) was inferred based on the amino acid sequences of the first segments. BAV, LNV and KDV are labeled in red, green and blue, respectively. The newly discovered KDV strain is labeled with a solid circle. Numbers associated with branches indicate the percentage of 1000 bootstrap replicates that support the existence of these branches. Branches with < 60% bootstrap support have been collapsed

Table 1 Sequence Identities Between discovered KDVs

Discussion

KDV is a member of the genus Seadornavirus, family Reoviridae. There are three species of viruses in this genus: BAV, KDV and LNV. Each of these viruses has been isolated from Aedes, Anopheles and Culex mosquito populations, but only BAV has been shown to cause infection in humans [9, 23]. KDV has been isolated from Culex, Anopheles, Armigeres and Aedes mosquitoes in Indonesia, and from Culex, Anopheles and Armigeres mosquitoes in China [13]. This range includes the tropics and subtropics. Here, we isolated the SDKL1625 from Anopheles sinensis mosquitoes in Shandong province in eastern China. This is the first time KDV has been isolated from the northern temperate zone.

The genome of SDKL1625 shows variation compared with the previously discovered KDV isolates QTM27331 and JKT-7075. Phylogenetic analysis showed that SDKL1625 was more closely related to QTM27331 (Odonata, China, 2013) than to JKT-7075 (Culex fuscocephalus, Indonesia, 1981). SDKL1625 was isolated from Anopheles sinensis mosquitoes in China, which have similar origin but a greater genetic distance from JKT-7075. Indeed, SDKL1625 and QTM27331 collected in China have high degrees of sequence identity. Although QTM27331 was isolated from Odonata, we don’t consider Odonata as the natural host of KDV just because mosquitoes belong to the diet of the Odonata. Previous KDV isolates from Yunnan province were cultured with cell lines and characterized morphologically, however, the full genome was not sequenced [13]. The sequence data in this study represent the first complete genome sequence of KDV isolated in China. A recent study of a plasma sample from a febrile adult indicates that the tropism of KDV includes humans [16]. Thus, it will be important to conduct further investigation of KDV, and in particular, the association between KDV and human diseases.

Mosquito-borne viruses, such as BAV and LNV, were originally discovered in one location, but later, more extensive surveys provided a more accurate representation of the distributions of these viruses. BAV was first isolated in Yunnan province in 1987 [5], and isolates including three genera and 10 species have now been obtained from mosquitoes in various provinces of China (Gansu, Shanxi, Liaoning, Yunnan and Beijing) [8]. LNV was first isolated in Jilin province but additional isolates were later obtained from several regions in northern China [10]. KDV has now been isolated in both northern and southwestern China, suggesting that the distribution of KDV in China was originally underestimated.

This study described the isolation and characterization, as well as genomic analyses, of KDV SDKL1625 and its phylogenetic relationship with members of the genus Seadornavirus. These discoveries have enriched our knowledge of the distribution of KDV in China, and prompted further arbovirus investigation. Serologic screening should be conducted to determine the presence, prevalence and distribution of KDV.

Conclusions

We describe the isolation, genome sequencing and characterization of the Kadipiro virus during arbovirus investigation in Shandong province. The genome of SDKL1625 revealed that the isolated virus was genetically close to strain isolated from Odonata in China. The sequence data in this study represent the first complete genome sequence of KDV isolated in China, which will further enrich our understanding of the molecular biological characteristics of KDV and seadornaviruses as well as the distribution of KDV in China.

Abbreviations

BAV:

Banna virus

CPE:

cytopathic effect

KDV:

Kadipiro virus

LNV:

Liao ning virus

PAGE:

Polyacrylamide gel electrophoresis

References

  1. The Double Stranded RNA Viruses. In Virus Taxonomy. San Diego: Academic Press; 2005:441–605.

  2. Attoui HMP, Becnel J, Belaganahalli S, Bergoin M, Brussaard CPCJ, Ciarlet M, del Vas M, Dermody TS, Dormitzer PR, Duncan R, Fang Q, Graham R, Guglielmi KM, HB HRM, Makkay A, Marzachı` C, MR MJ, Mohd Jaafar F, Mori H, Noordeloos AA, Omura T, RS PJT, Maan M, Stoltz D, Suzuki N, Upadhyaya NM, ZH WC. Family - Reoviridae. San Diego: Elsevier; 2012.

    Google Scholar 

  3. Reuter G, Boros Á, Delwart E, et al. Novel seadornavirus (family Reoviridae) related to Banna virus in Europe. Arch Virol. 2013;158:2163–7.

    Article  CAS  PubMed  Google Scholar 

  4. Wang JL, Li HC, He YW, et al. Isolation and genetic characterization of Mangshi virus: a newly discovered Seadornavirus of the Reoviridae family found in Yunnan Province, China. PLoS One. 2015;10(12):e0143601.

    Article  PubMed  PubMed Central  Google Scholar 

  5. Xu PT, Wang YM, Zuo JM, et al. Isolation of orbiviruses from sera of patients with unknown fever and encephalitis collected from Yunnan Province. Bing du xue bao. 1990;6(1):27–33. In China

    Google Scholar 

  6. Attoui H, Mohd JF, de Micco P, et al. Coltiviruses and seadornaviruses in North America, Europe, and Asia. Emerg Infect Dis. 2005;11(11):1673–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  7. Attoui H, Charrel RN, Billoir F, et al. Comparative sequence analysis of American, European and Asian isolates of viruses in the genus Coltivirus. J Gen Virol. 1998;79(10):2481–9.

    Article  CAS  PubMed  Google Scholar 

  8. Liu H, Gao XY, Fu SH, et al. Molecular evolution of emerging Banna virus. Infect Genet Evol. 2016;45:250–5.

    Article  CAS  PubMed  Google Scholar 

  9. Attoui H, Mohd JF, Belhouchet M, et al. Liao ning virus, a new Chinese seadornavirus that replicates in transformed and embryonic mammalian cells. J Gen Virol. 2006;87(1):199–208.

    Article  CAS  PubMed  Google Scholar 

  10. Lu Z, Liu H, Fu SH, et al. Liao ning virus in China. Virol J. 2011;8:282.

    Article  PubMed  PubMed Central  Google Scholar 

  11. Lv XJ, Mohd JF, Sun XH, et al. Isolates of Liao ning virus from wild-caught mosquitoes in the Xinjiang province of China in 2005. PLoS One. 2012;7(5):e37732.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  12. Attoui H, Billoir F, Biagini P, et al. Complete sequence determination and genetic analysis of Banna virus and Kadipiro virus: proposal for assignment to a new genus (Seadornavirus) within the family Reoviridae. J Gen Virol. 2000;81(6):1507–15.

    Article  CAS  PubMed  Google Scholar 

  13. Sun XH, Meng WS, Fu SH, et al. The first report of Kadipiro virus isolation in China. Bing du xue bao. 2009;25(3):173–7. In China

    CAS  PubMed  Google Scholar 

  14. Shi M, Lin XD, Tian JH, et al. Redefining the invertebrate RNA virosphere. Nature. 2016;540:539–43.

    Article  CAS  Google Scholar 

  15. Brown SE, Gorman BM, Tesh RB, et al. Coltiviruses isolated from mosquitoes collected in Indonesia. Virology. 1993;196(1):363–7.

    Article  CAS  PubMed  Google Scholar 

  16. Ngoi CN, Siqueira J, Li LL, et al. The plasma virome of febrile adult Kenyans shows frequent parvovirus B19 infections and a novel arbovirus (Kadipiro virus). J Gen Virol. 2016;97(12):3359–67.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  17. Lu BL. Insecta Vol.8 Diptera: Culicidae I. In: Editorial Committee of Fauna Sinica, Academia Sinica. Fauna Sinica. Beijing: Science; 1997. p. 432–3.

    Google Scholar 

  18. Lu BL. Insecta Vol.9 Diptera: Culicidae II. In: Editorial Committee of Fauna Sinica, Academia Sinica. Fauna Sinica. Beijing: Science; 1997. p. 31–3.

    Google Scholar 

  19. Li WJ, Wang JL, Li MH, et al. Mosquitoes and mosquito-borne arboviruses in the Qinghai-Tibet plateau--focused on the Qinghai area, China. Am J Trop Med Hyg. 2010;82(4):705–11.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  20. Wang HY, Wang HY, Fu SH, et al. Isolation and identification of a distinct strain of Culex Flavivirus from mosquitoes collected in mainland China. Virol J. 2012;9:73.

    Article  Google Scholar 

  21. Li H, Durbin R. Fast and accurate long-read alignment with burrows-wheeler transform. Bioinformatics. 2010;26(5):589–95.

    Article  PubMed  PubMed Central  Google Scholar 

  22. Jaafar FM, Attoui H, Mertens PPC, et al. Structural organization of an encephalitic human isolate of Banna virus (genus Seadornavirus, family Reoviridae). J Gen Virol. 2005;86(4):1147–57.

    Article  CAS  Google Scholar 

  23. Liu H, Li MH, Zhai YG, et al. Banna Virus, China, 1987-2007. Emerg Infect Dis. 2010;16(3):514–7.

    Article  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgments

We thank Han Yujie from Kenli County Center for Disease Control and Prevention for assistance with the mosquito collection.

Funding

This work was supported by grants from the National Key Research and Development Program (2016YFD0500401); the Special National Project on Research and Development of Key Biosafety Technologies (2016YFC1201904); the Development Grant of State Key Laboratory for Infectious Disease Prevention and Control (2015SKLID505, 2014SKLID103); the Zhi Fei research foundation programme on disease control and prevention technology from Shandong preventive medicine association (LYH2017–05). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Availability of data and materials

Sequence data obtained in this study is available in the GenBank (accession numbers MG590140-MG590151).

Author information

Authors and Affiliations

Authors

Contributions

LAG, LXJ, LGF and TZX collected the mosquitoes sample. ZWJ, FSH, SN, HY and LWW performed virus isolation and RT-PCR. ZWJ and WQY performed polyacrylamide gel electrophoresis experiment and analysis. SJD performed the virus electron microscopy. LF did the Bioinformatics analysis. ZWJ and LF were major contributors in writing the manuscript. WHY, ZL, LGD and XAQ conceived the study and drafted the manuscript. All authors read and approved the final manuscript.

Corresponding authors

Correspondence to Li Zhao or Huanyu Wang.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

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.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Zhang, W., Li, F., Liu, A. et al. Identification and genetic analysis of Kadipiro virus isolated in Shandong province, China. Virol J 15, 64 (2018). https://doi.org/10.1186/s12985-018-0966-y

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12985-018-0966-y

Keywords