Isolation and full-length genome analysis of mosquito-borne Manzanilla virus from Yunnan Province, China

There have been four strains on Manzanilla virus (MANV) identified to date. Here, we identify a novel MANV strain (DHL10M107) isolated from Culex tritaeniorhynchus Giles mosquitoes from Ruili city, Dehong prefecture, Yunnan Province, in the People’s Republic of China. The DHL10M107 L, M and S genes were sequenced at the nucleotide and deduced amino acid levels. The L, M and S gene sequences of DHL10M107 clustered with the MANV strains VN04-2108, TRVL3587, SA An 4165, and AV 782. DHL10M107 was most closely related to VN04-2108. Nucleotide homology ranged between 96 and 99% between DHL10M107 and VN04-2108. In terms of amino acid homology, all of the amino acid differences were in the L (96.3% homologous) and M (97.7% homologous) fragments. DHL10M107 is likely a MANV isolated from mosquitos in the Yunnan Province. This is the first reported isolation of MANV in mainland China.


Findings
Manzanilla virus (MANV) belongs to the Simbu serogroup of the genus Orthobunyavirus of family Bunyaviridae [1][2][3]. MANV is a single-stranded negative sense RNA virus that contains three RNA fragments: large (L), medium (M) and small (S). Four strains of MNV have been isolated from to date. Anderson et al. were the first group to isolate MANV (TRVL3587) [1] from the blood of a Howler monkey (Alouatta seniculus insularis) in Trinidad in 1954. More recently, in 2004 Bryant et al. isolated a MANV strain identified as the Cat Que virus (VN04-2108) from mosquitoes (Culex sp.) in Vietnam [2]. In 2014, Ladner et al. [3] reclassified two viruses as MANV strains, the Ingwavuma virus (SA An 4165) isolated from a South African spectacled weaver (Hyphanturgus ocularis) in 1959 [4], and the Mermet virus (AV 782) isolated from a North American purple martin (Progne subis) in the United States in 1964 [5].
In August 2010, mosquitoes were collected using the Kongfu Xiaoshuai trap (Wuhan Jixing Environmental Protection, Scientific, and Technological LLC) from cattle barns in the suburb of Ruili city, Dehong prefecture, Yunnan Province, People's Republic of China. A total of 425 mosquitoes representing 10 species in four genera were collected. Of these, 318 (74.82%) were Culex tritaeniorhynchus Giles and the remaining 107 were from the other nine species (25.18%). The mosquito samples were sorted by species into 15 pools and ground together. The supernatant was used to inoculate Baby hamster kidney (BHK-21) cells to isolate any viruses present as previously described [15]. One of the pools containing C. tritaeniorhynchus (DHL10M107) caused significant cytopathic effects (CPE) in BHK-21 cells, characterized by the cells shrinking, rounding, and eventually floating.
The viral RNA was extracted from the culture supernatant of the DHL10M107 isolate, using the QIAamp viral RNA mini kit (Qiagen, Valencia, CA, USA) per the manufacturer's instructions. cDNA was prepared with Ready-To-Go You-Prime First-Strand Beads™ (American Amersham Pharmacia Biotech, Piscataway, NJ, USA), and amplified with MANV specific primers (see Additional file 1: Table s1). The reactions contained 1 µL cDNA template, 10 µL 2 × GCI Buffer, 0.5 µL 2.5 mM dNTP, 0.4 µL each of the upstream and downstream primers (10 mM), 0.4 µL ExTaq enzyme (5 U/µL), and 7.3 µL deionized water. The reaction conditions were: pre-denaturing at 95°C for 3 min, followed by 32 cycles of 95°C for 30 s, 58°C for 40 s, and 72°C for 1-2 min, and extension at 72°C for 10 min. The amplification was assessed using 2 µL of the product for 1% agarose gel electrophoresis. The whole genome sequence of the DHL10M107 strain was obtained (Beijing Liuhe Genomics Technology Co., Ltd. Shenzhen Branch). The L, M and S sequences of DHL10M107 strain were submitted to GenBank under accession nos. KP016012-KP016014. The sequence fragments were spliced, edited and corrected using SeqMan (DNAstar software package) and the phylogenetic analysis and nucleotide sequence analysis were performed using ClustalX (version 1.8), MEGA5.0, DNAStar alignment package, and MegAlign software.
Phylogenetic analysis of the viral L and M gene sequences indicated that DHL10M107 clustered with MANV (VN04-2108) and an Oya virus (OYV) strain (SC0806) isolated from mosquitos in Sichuan province, China [17], although it was more closely related to VN04-2108 (Figure 1a, b). However, in terms of the S fragments, while DHL10M107 was in the same clade as VN04-2108 and SC0806, it formed a secondary clade, and the S fragments of VN04-2108 and SC0806 were more closely related (Figure 1c). Regardless of the genome fragment, DHL10M107 was closely related to three MANV strains (TRVL3587, SA An 4165, and AV 782; Figure 1). The OYV reported by Kono et al. [18] had a partial sequence for the S fragment in GenBank (accession no. AB075611). Similar to Figure 1c, the phylogenetic tree constructed using this sequence indicated AB075611, DHL10M107, VN04-2108, and SC0806 were in the same clade and had similar evolutionary relationships (Figure 1d). We then compared the homology of the open reading frame (ORF) and amino acid sequences in DHL10M107 to five other Orthobunyavirus strains (Table 1). There was minimal nucleotide and amino acid homology between DHL10M107 and Tahyna virus, BATV, and Oropouche virus (Table 1). In terms of nucleotide homology, the similarity of the DHL10M107 L fragment to VN04-2108 and SC0806 was 96.3 and 92.6%, respectively. The M fragment (VN04-2108: 97.7% and SC0806: 97.2%) and S fragment (VN04-2108: 98.7% and SC0806: 98.9%) had similarly high levels of homology. Amino acid homology was between 99 and 100% for DHL10M107, VN04-2108, and SC0806 ( Table 1). All of the amino acid differences were in the L and M fragments. In the L fragment, DHL10M107 differed from VN04-2108 by six amino acids and SC0806 by 11 amino acids. In the M fragment the difference was four (VN04-2108) and five amino acids (SC0806; Table 2).
In conclusion, the high degree of similarity in the whole genome sequence between DHL10M107 and the four MANV strains in GenBank, but not other viruses in the Simbu serogroup, confirmed it was a MANV strain. This is the first report of MANV in mainland China. Currently, the records pertaining to MANV in the International Committee on Taxonomy of Virus (ICTV) [19] excludes Oya virus. Here, we identified two virus strains named OYV in GenBank. One was a partial OYV S fragment sequence (AB075611) isolated from the lungs of pigs thought to have Nipah virus infection in Malaysia [8]. The second strain was the entire sequence of OYV SC0806 (JX983192, JX983193 and JX983194) isolated from mosquitos in Sichuan province, China [7].
Phylogenetic analysis indicated that both OYV sequences clustered with the MANV strains. They were most closely related to VN04-2108 and DHL10M107, suggesting that OYV (SC0806 and AB075611) was also a member of the MANV complex.