Skip to main content

Analysis of reassortant and intragenic recombination in Cypovirus

Abstract

Cypoviruses (CPVs) are RNA viruses with segmented double-stranded genome and major pathogens of various insects, including economic insects like silkworms and pest insects for agricultural crops and forests. Genome reassortment and recombination are common phenomenon for viruses as a mechanism to expand host range and increase virulence. In the present study, we analyzed the reassortant and recombination events for CPVs. The results showed that two genome segments (S1 and S4) of BmCPV1-YN shared higher nucleotide identity with the corresponding segment of BmCPV1-I while others were all more closely to BmCPV1-SZ, suggesting BmCPV1-YN was originated from reassortant events between BmCPV1-I and BmCPV1-SZ. Recombination analyses revealed that S6 of BmCPV1-YN was a recombinant segment derived from BmCPV1-I and BmCPV1-SZ, and S10 of DpCPV1 was a recombinant segment emerged from BmCPV1-I and LdCPV1. Our findings provide the evidence for the fact that CPVs could undergo reassortant and recombinant events and enrich the knowledge about etiology and molecular epidemiology of CPVs.

Main text

Cypoviruses (CPVs) are important viral pathogens that infect the midgut epithelial cells of insects, cause larval developmental retardation and decreased fecundity, and have been isolated from more than 250 insect species of Lepidoptera, Hymenoptera, Diptera and Coleoptera [1, 2]. CPVs are segmented double-stranded RNA (dsRNA) viruses, belonging to the genus Cypovirus in the subfamily Spinareoviridae and family Reoviridae, which have so far been classified into 22 elecropherotypes based on their electrophoretic migration patterns [3, 4]. There are usually 10 dsRNA segments in the genome, but some CPVs consist of 9, 11, 12 and even 16 segments [4,5,6]. For each segment, a high identity was revealed among viruses classified in the same electropherotype, however, no nucleotide sequence identity was found among CPVs of distinct electropherotypes [5]. Reassortant and recombination have been confirmed in various RNA viruses and several DNA viruses as a mechanism to adapt to changing environment, expand host range and increase virulence, resulting in the genetic diversity of viruses [7,8,9]. Reassortant and intragenic recombination events have also been reported in most viruses of Reoviridae, such as bluetongue virus (BTV), African horse sickness virus (AHSV), rotaviruses (RVs) and so on [10,11,12]. It is essential for two or more viruses to infect the same host cells, so that their genetic material might be exchanged for reassortment. Interestingly, it has been confirmed that insects could be co-infected with several CPVs [13], but there was so far no analysis or evidence of reassortant and recombination for CPVs. In the present study, we first analyzed the majority of genomic sequences of CPVs deposited in National Center for Biotechnology Information (NCBI) and then identified their reassortant and intragenic recombination events.

About 200 sequences of different segment of CPVs were acquired from GenBank database, and 15 CPV isolates have been determined for their complete genome sequences (Table S1). To better understand the relationship of CPV strains, firstly, we conducted the sequence comparative analysis by using ClustalW in Lasergene software (DNASTAR Inc., Madison, USA). Higher nucleotide identity was found among CPV strains of the same electropherotype and isolated from the same host species, for example, each genomic segment of CPV1s isolated from silkworm, Bombyx mori (i.e. BmCPV1-I, BmCPV1-SZ and BmCPV1-YN) share 88.6–99.8% nucleotide identity with each other. The lower nucleotide identity was found among CPVs isolated from different host species (each segment of BmCPV1-I and DpCPV1 share 78.7–90.1% nucleotide identity) although they are classified in the same electropherotype. No or few nucleotide sequence identities was found among CPVs of different electropherotypes just as described before [5]. Besides, the protein encoded by the same segment may be different among different CPV electropherotypes, for example, S1 of BmCPV1 encodes major capsid protein whereas S1 of HaCPV14 encodes RNA-dependent RNA polymerase (RdRP). Therefore, it is more meaningful to analyze the reassortant and recombination events among CPV isolates within the same electropherotype.

Each genomic segment of five strains in CPV1, three strains in CPV5 and three strains in CPV14 were compared with each other to analyze the reassortant events. Based on the results of nucleotide identity by using ClustalW in Lasergene software (DNASTAR Inc., Madison, USA), the genomic sequences of 8 segments (S2, S3, S5-S10) in BmCPV1-YN share the highest identity (98.3–99.8%) with the corresponding segments of BmCPV1-SZ, whereas S1 and S4 of BmCPV1-YN show the highest identity (98.3, 98.2%) with BmCPV1-I (Table 1), suggesting BmCPV1-YN might undergo reassortant event between BmCPV1-I and BmCPV1-SZ. To further confirm the result, the phylogenetic tree was constructed based on each segment of CPV1 genomes by the MEGA 6 software using the maximum likelihood method with a General Time Reversible (GTR) model of nucleotide substitution and bootstrap tests of 1000 replicates [14]. As shown in Fig. 1, BmCPV1-I, BmCPV1-SZ and BmCPV1-YN are clustered into a larger branch, 8 segments (S2, S3, S5-S10) of BmCPV1-YN are closely related to BmCPV1-SZ, while S1 and S4 of BmCPV1-YN are closer to BmCPV1-I. These findings clearly indicate that BmCPV1-YN is a reassortant virus emerged from segment exchange between BmCPV1-I and BmCPV1-SZ, with 8 genome segments (S2, S3, S5-S10) from BmCPV1-SZ and 2 segments (S1 and S4) from BmCPV1-I. Complete genome sequences of three strains in electropherotype 5, HaCPV5-C, OpCPV5 and TpCPV5, were also analyzed. Each segment of HaCPV5-C genome exhibits higher identity (82.9–99.1%) with OpCPV5, except for S10. The sequence of S10 of HaCPV5-C shows highest identity (99.8%) with EsCPV5 (GenBank accession no. J04338), but the relationship between these two viruses is still unclear due to the unavailability of other segment sequence of EsCPV5. Compared to the high nucleotide identity between HaCPV5-C and OpCPV5, TpCPV5 shows lower identity (40–80%) with either HaCPV5-C or OpCPV5, and TpCPV5 is clustered into a separate branch based on the phylogenetic tree (data not shown), which indicates that no reassortant events among these three strains. According to the similar analysis, no reassortment is identified among HaCPV14, LdCPV14 and TaCPV14. The intragenic recombination of individual segment within the same electropherotype was further analyzed by using RDP v4.16 and confirmed by using the software SimPlot v3.5.1 [15] and the blast results in NCBI. Two recombinant segments, S6 of BmCPV1-YN and S10 of DpCPV1, were identified (Fig. 2). S6 of BmCPV1-YN is a recombinant segment derived from the recombination event between BmCPV1-I and BmCPV1-SZ, and S10 of DpCPV1 was emerged from the recombination event between BmCPV1-I and LdCPV1.

Table 1 Nucleotide identities (%) between BmCPV1-YN and BmCPV1-I, BmCPV1-SZ
Fig. 1
figure1

Phylogenetic analyses of CPV1. The phylogenetic trees were constructed based on the nucleotide sequences of each segment of CPV1. Reliability of the tree was assessed by bootstrap analysis of 1000 replications. Scale bar indicates nucleotide substitutions per site

Fig. 2
figure2

Intragenic recombination analysis. A S6 of BmCPV1-YN; B S10 of DpCPV1. Comparisons of genetic similarity between recombinant segment and parental segment were made using SimPlot. The S6 of BmCPV1-YN (a) and S10 of DpCPV1 (b) were chosen as query sequence

Reassortant and recombination are very important mechanism for virus evolution to escape host immune surveillance, adapt to changing environment, expand host range and so on, which has been analyzed and reported in most viruses of Reoviridae. The recombination events of BTV was documented in 2010 [10] and Ngoveni gave the evidence of intragenic recombination in AHSV in 2019 [12]. A study by Chen firstly reported on the occurrence of reassortant and recombinant virus between human and pig RVA strains [16]. As dsRNA viruses belonging to Reoviridae, some CPVs can infect several kinds of insects, while some insects can also be infected by different CPVs. Therefore, the hypothesis of reassortant and recombination events among CPVs was proposed. To our best knowledge, nucleotide similarity among CPVs of different electropherotypes is lower, therefore only the reassortant and recombination events in CPVs within the same electropherotype were analyzed in the present study. Our results show that BmCPV1-YN is a reassortant virus between BmCPV1-I and BmCPV1-SZ, and segment 6 of BmCPV1-YN is produced by recombination events between BmCPV1-I and BmCPV1-SZ (Fig. 3a). So far, only three CPV strains were isolated from the silkworm, Bombyx mori. As in other viruses reported in the similar papers, more reassortant and intragenic recombination events might be also possible in BmCPVs if more virus strains are identified and more sequences are deposited to NCBI in the future. It is also possible that BmCPV1-YN might be a resultant virus derived from reassortant events of some other intermediate cypoviruses, but we could not compare BmCPVs with CPVs isolated from other insect because it is not clear whether those CPVs with available genome sequence are infectious to Bombyx mori. The evidence of interspecies reassortant and recombination events of rotavirus have been reported [16]. Here the segment 10 of DpCPV1 is identified as a recombinant sequence derived from BmCPV1-I and LdCPV1 (Fig. 3b). It has been reported that carboxy terminus of BmCPV polyhedron is very important for its nuclear localization and modification of crystallization pattern [17], so the recombination events within the S10 terminus of DpCPV1 may be also associated with the function of polyhedron. Although more than 200 sequences were analyzed in the present study, the rate of reassortant and recombination events among CPVs are relatively low, which may be due to few available CPVs sequences within same electropherotype, lower degree of genetic compatibility or the poor competitiveness of nascent virus. Besides, for influenza A viruses and rotaviruses, the reassortant has contributed to viral pathogenic and zoonotic potentials, so the pathogenicity of reassortant and recombinant CPVs and ongoing epidemiological surveillance are required to be explored in the future.

Fig. 3
figure3

Reassortant and intragenic recombination model. a BmCPV1-YN; b DpCPV1. The individual segment is shown as “rectangle” with different colors and the length of “rectangle” indicates the size of segments

Availability of data and materials

All sequences obtained in this study are available in GenBank.

Abbreviations

AHSV:

African horse sickness virus

Bm:

Bombyx mori

BTV:

bluetongue virus

CPV:

cypovirus

Dp:

Dendrolimus punctatus

dsRNA:

double-stranded RNA

Es:

Euxoa scandens

Ha:

Heliothis armigera

Ld:

Lymantria dispar

NCBI:

National Center for Biotechnology Information

Op:

Orgyia pseudotsugata

RDP:

Recombination Detection Program

RdRP:

RNA-dependent RNA polymerase

RVs:

rotaviruses

Ta:

Thyrinteina arnobia

Tp:

Thaumetopoea pityocampa

References

  1. 1.

    Belloncik S. Cytoplasmic polyhedrosis viruses-Reoviridae. Adv Virus Res. 1989;37:173–209.

    Article  CAS  Google Scholar 

  2. 2.

    Magnoler A. Effects of a cytoplasmic polyhedrosis on larval and postlarval stages of the gypsy moth, Porthetria dispar. J Invertebr Pathol. 1974;23:263–74.

    Article  CAS  Google Scholar 

  3. 3.

    Hill CL, Booth TF, Prasad BV, Grimes JM, Mertens PP, Sutton GC, et al. The structure of a cypovirus and the functional organization of dsRNA viruses. Nat Struct Biol. 1999;6:565–8.

    Article  CAS  Google Scholar 

  4. 4.

    Zhou Y, Qin T, Xiao Y, Qin F, Lei C, Sun X. Genomic and biological characterization of a new Cypovirus isolated from Dendrolimus punctatus. PLoS One. 2014;9:e113201.

    Article  CAS  Google Scholar 

  5. 5.

    Li Y, Zhang J, Li Y, Tan L, Chen W, Luo H, et al. Phylogenetic analysis of Heliothis armigera cytoplasmic polyhedrosis virus type 14 and a series of dwarf segments found in the genome. J Gen Virol. 2007;88:991–7.

    Article  CAS  Google Scholar 

  6. 6.

    Sinha-Datta U, Chavali VRM, Ghosh AK. Molecular cloning and characterization of Antheraea mylitta cytoplasmic polyhedrosis virus polyhedrin gene and its variant forms. Biochem Biophys Res Commun. 2005;332:710–8.

    Article  CAS  Google Scholar 

  7. 7.

    Chen Y, Trovao NS, Wang G, Zhao W, He P, Zhou H, et al. Emergence and evolution of novel Reassortant influenza a viruses in canines in southern China. MBio. 2018;9:3.

    Google Scholar 

  8. 8.

    Franzo G, Cortey M, Segalés J, Hughes J, Drigo M. Phylodynamic analysis of porcine circovirus type 2 reveals global waves of emerging genotypes and the circulation of recombinant forms. Mol Phylogenet Evol. 2016;100:269–80.

    Article  Google Scholar 

  9. 9.

    Zhang Z, Zhou L, Ge X, Guo X, Han J, Yang H. Evolutionary analysis of six isolates of porcine reproductive and respiratory syndrome virus from a single pig farm: MLV-evolved and recombinant viruses. Infect Genet Evol. 2018;66:111–9.

    Article  Google Scholar 

  10. 10.

    He CQ, Ding NZ, He M, Li SN, Wang XM, He HB, et al. Intragenic recombination as a mechanism of genetic diversity in bluetongue virus. J Virol. 2010;84:11487–95.

    Article  CAS  Google Scholar 

  11. 11.

    Esona MD, Roy S, Rungsrisuriyachai K, Sanchez J, Vasquez L, Gomez V, et al. Characterization of a triple-recombinant, reassortant rotavirus strain from the Dominican Republic. J Gen Virol. 2017;98:134–42.

    Article  CAS  Google Scholar 

  12. 12.

    Ngoveni HG, van Schalkwyk A, Koekemoer JJO. Evidence of intragenic recombination in African horse sickness virus. Viruses. 2019;11:654.

    Article  CAS  Google Scholar 

  13. 13.

    Li Y, Tan L, Li Y, Chen W, Zhang J, Hu Y. Identification and genome characterization of Heliothis armigera cypovirus types 5 and 14 and Heliothis assulta cypovirus type 14. J Gen Virol. 2006;87:387–94.

    Article  CAS  Google Scholar 

  14. 14.

    Tamura K, Stecher G, Peterson D, Filipski A, Kumar S. MEGA6: molecular evolutionary genetics analysis version 6.0. Mol Biol Evol. 2013;30:2725–9.

    Article  CAS  Google Scholar 

  15. 15.

    Lole KS, Bollinger RC, Paranjape RS, Gadkari D, Kulkarni SS, Novak NG, et al. Full-length human immunodeficiency virus type 1 genomes from subtype C-infected seroconverters in India, with evidence of intersubtype recombination. J Virol. 1999;73:152–60.

    Article  CAS  Google Scholar 

  16. 16.

    Chen D, Zhou L, Tian Y, Wu X, Feng L, Zhang X, et al. Genetic characterization of a novel G9P[23] rotavirus a strain identified in southwestern China with evidence of a reassortment event between human and porcine strains. Arch Virol. 2019;164:1229–32.

    Article  CAS  Google Scholar 

  17. 17.

    Cao G, Meng X, Xue R, Zhu Y, Zhang X, Pan Z, et al. Characterization of the complete genome segments from BmCPV-SZ, a novel Bombyx mori cypovirus 1 isolate. Can J Microbiol. 2012;58:872–83.

    Article  CAS  Google Scholar 

Download references

Acknowledgements

This work was supported by the National Natural Science Foundation of China (Grant No. 31572463) and the available sequences deposited by researchers from all over the world.

Funding

This study was financially supported by the National Natural Science Foundation of China (Grant No. 31572463).

Author information

Affiliations

Authors

Contributions

Zhendong Zhang and Xijie Guo interpreted the data and drafted the work. Xijie Guo coordinated the study. Ning Li and Chengxiang Hou confirmed the results of reassortant and recombination events analysis. Kun Gao and Xudong Tang contributed to collect the sequence of cypoviruses deposited into GenBank. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Xijie Guo.

Ethics declarations

Ethics approval and consent to participate

“Not applicable” in this section.

Consent for publication

“Not applicable” in this section.

Competing interests

All authors have declared no conflict of interest.

Additional information

Publisher’s Note

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

Supplementary information

Additional file 1: Table S1.

GenBank accession number of each segment of 15 CPV strains. The GenBank accession number of each segment of 15 CPV strains were acquired from GenBank database. Note: BmCPV1-I, Bombyx mori CPV 1 strain I, Japan; BmCPV1-SZ, Bombyx mori CPV 1 isolate Suzhou, China; BmCPV1-YN, Bombyx mori CPV 1 isolate Yunnan, China; DpCPV1, Dendrolimus punctatus CPV 1; LdCPV1, Lymantria dispar CPV 1; IiCPV2, Inachis io CPV 2; HaCPV5-C, Heliothis armigera CPV 5 isolate China; OpCPV5, Orgyia pseudotsugata CPV 5; TpCPV5, Thaumetopoea pityocampa CPV 5; HaCPV14, Heliothis armigera CPV 14; LdCPV14, Lymantria dispar CPV 14; TaCPV14, Thyrinteina arnobia CPV 14; TnCPV15, Trichoplusia ni CPV 15 (segment 11: NC_002566); DpCPV22, Dendrolimus punctatus CPV 22 (segment 11–16: KJ191114, KJ191115, KJ191116, KJ191117, KJ191118, KJ191119); DnCPV-NC, Daphnis nerii CPV isolate Nanchang, China.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. 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 in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Zhang, Z., Li, N., Hou, C. et al. Analysis of reassortant and intragenic recombination in Cypovirus. Virol J 17, 48 (2020). https://doi.org/10.1186/s12985-020-01321-1

Download citation

Keywords

  • Cypovirus
  • Reassortment
  • Recombination