Next Article in Journal
Lack of Middle East Respiratory Syndrome Coronavirus Transmission in Rabbits
Next Article in Special Issue
Characterization of Signal Sequences Determining the Nuclear/Nucleolar Import and Nuclear Export of the Caprine Arthritis-Encephalitis Virus Rev Protein
Previous Article in Journal
Characterization of a New Member of Alphacoronavirus with Unique Genomic Features in Rhinolophus Bats
Previous Article in Special Issue
Comparative Analysis of Different Serological and Molecular Tests for the Detection of Small Ruminant Lentiviruses (SRLVs) in Belgian Sheep and Goats
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Characterization of EIAV env Quasispecies during Long-Term Passage In Vitro: Gradual Loss of Pathogenicity

1
State Key Laboratory of Veterinary Biotechnology, Harbin Veterinary Research Institute of Chinese Academy of Agricultural Sciences, Harbin 150069, China
2
Key Laboratory of Special Animal Epidemic Disease of Ministry of Agriculture, Institute of Special Animal and Plant Sciences, Chinese Academy of Agricultural Sciences, Changchun 130112, China
3
Department of Microbiology, Harbin Medical University, Harbin 150069, China
*
Authors to whom correspondence should be addressed.
These authors contribute equally to the work.
Submission received: 15 February 2019 / Revised: 8 April 2019 / Accepted: 17 April 2019 / Published: 24 April 2019
(This article belongs to the Special Issue Nonprimate Lentivirus)

Abstract

:
As the only widely used live lentiviral vaccine, the equine infectious anima virus (EIAV) attenuated vaccine was developed by in vitro passaging of a virulent strain for 121 generations. In our previous study, we observed that the attenuated vaccine was gradually selected under increased environmental pressure at the population level (termed a quasispecies). To further elucidate the potential correlation between viral quasispecies evolution and pathogenesis, a systematic study was performed by sequencing env using several methods. Some key mutations were identified within Env, and we observed that increased percentages of these mutations were accompanied by an increased passage number and attenuated virulence. Phylogenetic analysis revealed that env mutations related to the loss of virulence might have occurred evolutionarily. Among these mutations, deletion of amino acid 236 in the V4 region of Env resulted in the loss of one N-glycosylation site that was crucial for virulence. Notably, the 236-deleted sequence represented a “vaccine-specific” mutation that was also found in wild EIAVLN40 strains based on single genome amplification (SGA) analysis. Therefore, our results suggest that the EIAV attenuated vaccine may originate from a branch of quasispecies of EIAVLN40. Generally, the presented results may increase our understanding of the attenuation mechanism of the EIAV vaccine and provide more information about the evolution of other lentiviruses.

1. Introduction

As the simplest member of the lentivirus family, equine infectious anemia virus (EIAV) shares similar features with other lentiviruses, including its genomic structure, life cycle, cell tropism and antigen evolution [1,2]. The uniqueness of EIAV lies in the fact that most infected horses become asymptomatic carriers after recurrent febrile episodes, which are associated with a high viral load and anemia, and thus the horses achieve natural immunologic control [3]. Notably, an attenuated vaccine strain has been developed that successfully controls EIAV prevalence in China. The vaccine was generated by serial passage of a virulent EIAV strain (EIAVLN40) for 16 generations in vivo, followed by passaging for as many as 121 passages in vitro [4,5]. Due to these peculiarities, the EIAV attenuation system may provide an ideal model for elucidating the correlates of genomic evolution and immunologic control.
For RNA viruses, mutation swarms are generated rapidly through error-prone viral replication due to inaccuracy of the RNA polymerase [6]. This mutant swarm, which is also termed a viral quasispecies, appears as a whole population with a mutation-selection balance and is significantly involved in viral pathogenesis [7,8]. Viral quasispecies are also essential for viral survival because they yield beneficial viral phenotypes under in vivo selection. Interestingly, the mutation spectra are higher after long-term passage in vitro, which causes a change in virulence [9,10,11]. Hence, the relationship between RNA viral quasispecies evolution and viral pathogenesis is an intriguing problem.
As discussed previously, multi-position mutations located in the LTR and in env appeared during the EIAV attenuation process in vitro, especially in the viral env gene, which developed eight hyper variable regions located in the V3 and V4 regions [1,12,13,14]. Further analysis showed that these variations were related to viral pathogenesis during the EIAV attenuation process. A virulence-correlated parallel in env variation was also observed during the EIAV attenuation process [4]. In another study, we observed that the LTR showed a similar pattern at the population level [15]. However, the LTR is a noncoding region which could not reflect the predominant antigen gene of the EIAV evolutionary pattern. Hence, env was chosen as the target of evolutionary selection in this study to address this important question. The aim of this study was to characterize env quasispecies evolution and further investigate the related mechanism between virulence attenuation and env quasispecies. Our results will be of great interest for understanding the evolutionary mechanism of the EIAV attenuated vaccine and defining vaccine development strategies for other lentiviruses.

2. Materials and Methods

2.1. Study Subjects

All samples were stored at the Harbin Veterinary Research Institute, Chinese Academy of Agricultural Sciences (CAAS). These samples included a virulent strain (EIAVLN40) that caused an infection in horses with 100% mortality at a dose of 1 × 105 TCID50; this strain was initially isolated from an EIA positive horse and passed for 16 generations in horses, resulting in three representative strains (EIAVDLV34, EIAVDLV62, and EIAVDLV92) with 100%, 100%, and 9.1% mortality that stemmed from EIAV passage in vitro for 34, 62, and 92 generations, respectively (Figure S1) [16]. Additionally, we included a vaccine strain (EIAVDLV121) that provided 85% protection against EIAVLN40 challenge and a full-length infectious molecular clone (pLGFD3-8) constructed from a vaccine strain and then passaged in fetal donkey dermal (FDD) cells for 3 generations. Animal experiments showed that this infectious clone was avirulent [17].

2.2. Viral RNA Extraction and cDNA Synthesis

Viral RNA was extracted from 140 μL of the virus samples using the Viral RNA Mini Kit (QIAamp, Dusseldorf, Germany) according to the manufacturer’s protocol. Reverse transcription of RNA to single-stranded cDNA was performed using the SuperScriptTM IV Reverse Transcriptase System (Life Technologies, Carlsbad, USA). First, 1 μg of RNA, dNTPs (0.5 mM each), and 0.5 μM of the env NR primer (5’-CAGCTACAATGGCAGCTATTATAGCAG-3’; nucleotides (nt) 6702 to 6676 of the EIAV sequence) were incubated for 5 min at 65 °C to denature the RNA secondary structure. First-strand cDNA synthesis was carried out in 20 μL reaction mixtures with 5× SSIV buffer, 5 mM DTT, 2 U/μL of an RNase inhibitor (RNaseOUT) (Life Technologies, Carlsbad, USA), and 10 U/μL of SuperScriptTM IV (Life Technologies, Carlsbad, USA). The reaction mixture was incubated at 55°C for 15 min and heated to 80°C for 10 min to complete the reverse transcription reaction, followed by RNase H (Life Technologies, Carlsbad, USA) digestion at 37 °C for 20 min. The synthesized cDNA was used immediately for PCR or stored at −80 °C.

2.3. Bulk PCR

The full-length hypervariable region of the env cassette was amplified by nested PCR from the viral cDNA. The specific method was as follows. First, 0.6 μL of bulk cDNA was used for the first-round PCR in a 20 μL volume. The PCR was performed using the KOD FX (Toyobo, Osaka, Japan) system, which included 2× PCR buffer, 0.4 mM dNTPs, and 0.1 μM of the env NF (5’-GAAGGCCATCAGGGAGGGAAG-3’; (nt) 5240 to 5260) and env NR primers. The cycling conditions were as follows: 94 °C for 2 min, followed by 35 cycles of 94 °C for 30 s, 58 °C for 30 s, and 68 °C for 1 min 30 s, and a final extension of 68 °C for 2 min. The second-round PCR was performed using 0.6 μL of the first-round PCR product and the V3-1 (5’-ACAAACATATACAGGACATCT-3’; (nt) 5808 to 5828) and V4-1 (5’-CTCCAATATTCCAAGAAATAC-3’; (nt) 6218 to 6198) primers under the same conditions used for the first-round PCR. The final PCR products were analyzed by 1% agarose gel electrophoresis. The recycled products were ligated into pMD18-T (TaKaRa, Dalian, China). The constructed plasmids were transformed into E. coli DH5α (TaKaRa, Dalian, China) and spread onto a plate. A total of 98 positive clones identified by PCR were sequenced by COME (Jilin, China).

2.4. Probe Design and Real-Time PCR

Two specific probes and primers were synthesized by General Biosystems (Anhui, China). The first probe was the 6-carboxy-fluorescein (FAM) probe (probe-FAM: 5’-TGCAGCAAAGTAAAAACACTTGGAT-3’; (nt) 6010 to 6034), which targeted the env-Δ236D-phenotype sequence, and the second was the hexachloride fluorescein (HEX) probe (probe-HEX: 5’-TGCAGCAAAGCGATAATAACACTTG-3’; (nt) 6010 to 6034), which targeted the env-236D-phenotype sequence. The env-Δ236D-phenotype and env-236D-phenotype sequences were named depending on the presence of the 236D deletion in the EIAV Env. The DNA fragments of env-Δ236D-phenotype and env-236D-phenotype were amplified with primer F (5’-ATAGGAGGTAGACTAAATGGTTCAGG-3’; (nt) 5634-5668) and primer R (5’-AAACAAAAAGAATGGAGGTTGGACA-3’; (nt) 6143-6119) and ligated into the pMD-18T (TaKaRa, Dalian, China) vector. Then, the two plasmids were serially diluted to 102 to 109 copies/μL as standards to construct real-time PCR standard curves. The experiments were performed using 1 μM of each primer, 10 μL of premix Ex Taq (TaKaRa, Dalian, China), 0.25 μM of the HEX-labeled or FAM-labeled probe, 2 μL of the cDNA samples and RNase-free water to a final volume of 20 μL. The samples were amplified using TaqMan-based real-time PCR on the ABI Prism 7500 sequence detection system (Applied Biosystems GmbH, Thermo Fisher Scientific, Wilmington, USA) with one cycle at 95 °C (5 min), followed by 40 cycles at 95 °C (15 s), 58 °C (30 s), and 72 °C (20 s). The calculated efficiencies for all primers were determined by dilution experiments and ranged from 97% to 98%; thus, the target sequences were amplified with similar efficiencies. All samples were run with at least three duplicates.

2.5. SGA (Single Genome Amplification)

According to the Poisson distribution [18], the DNA dilution yielding positive PCR products in no more than 30% of the wells will contain one amplicon resulting from single molecule amplification in over 80% of these wells. The viral cDNA was serially diluted to yield a single copy by nested PCR as described for the bulk PCR. All products derived from cDNA dilutions conforming to the above description were sequenced with the V3-1 and V4-1 primers by GENEWIZ (Suzhou, China).

2.6. Data Processing and Analysis

The sequences generated in this study have been submitted to GenBank (GenBank acc. no. MK268242 to MK268339 for bulk PCR and MK278920 to MK279296 for SGA). A preliminary sequence analysis was conducted using the SeqMan and MegAlign programs from the Lasergene DNA & Protein analysis software (version 7.0, DNASTAR Inc., Madison, WI, USA), the BioEdit Sequence Alignment Editor (version 7.2.5), and ClustalX (version 2.1). A maximum likelihood tree was constructed using MEGA 7.0 [19]. The amino acid alignment was automatically generated by MargFreq (version 1.02).

3. Results

3.1. General Characteristics of EIAV env Variability during the Course of EIAV Attenuation Using Bulk PCR

To verify the characteristics of EIAV env variation, some env genes derived from EIAV strains during the critical stages of the attenuation process were analyzed using bulk PCR [20,21]. A phylogenetic tree was constructed based on 98 env sequences amplified from EIAV strains with differing virulence (EIAVLN40, EIAVDLV34, EIAVDLV62, EIAVDLV92, and EIAVDLV121) (the morbidity of all strains showed in Figure S1). The general features of env variability are shown in Figure 1. Obviously, individual clusters accompanied by a decrease in virulence were generated. In other words, the evolutionary pattern of the env gene was consistent with the loss of virulence. Notably, the sequences derived from EIAVDLV62 displayed a separated distribution between EIAVLN40 and EIAVDLV121. This result suggested that EIAVDLV62 was an important viral evolutionary state that was associated with EIAV pathogenicity. The details of the sequences that also displayed several mutations in each cluster are presented in Table 1. The isolates from each generation displayed the same mutation types (insertions and deletions) but had different mutation frequencies and sites. We observed that the proportion of strains carrying this mutation gradually increased with the serial passage number in vitro (EIAVDLV34, EIAVDLV62, EIAVDLV92, and EIAVDLV121 = 6.25%, 18.5%, 42.1%, and 100%, respectively). We also noted that most mutations were random except for a successful amino acid deletion (236D) located in the V4 domain of Env compared to the reference strain sequence (EIAVLN40), as shown in Figure 2. This mutation was also stable in horses inoculated with EIAVDLV121 [22]. Thus, we speculated that this site might represent a vaccine-specific mutation.

3.2. Precise Characteristics of the EIAV env Mutation Distribution during the Attenuation Process

Once the vaccine-specific mutation was verified, we wanted to investigate whether this mutation could be detected during serial passage using a higher sensitivity method. These data will further benefit our understanding of the relationship between env quasispecies and virulence. Firstly, we adopted double-probe real-time PCR to accurately determine the proportion of strains with the env-Δ236D-phenotype and env-236D-phenotype sequences (which refer to virus sequences without and with the 236D mutation) in different EIAV generations [23,24]. An env-236D-phenotype-specific probe (HEX) and an env-Δ236D-phenotype-specific probe (FAM) were designed based on the V4 domain sequences (Figure S2). To avoid the influence of the initial viral copy number on the detection sensitivity, samples containing different initial viral copy numbers were analyzed in parallel (102 copies/μL, 103 copies/μL and 104 copies/μL). As presented in Figure S3, the env-Δ236D-phenotype sequence was detected with low copy numbers from all virulent strains (EIAVDLV34 and EIAVDLV62) except for EIAVLN40, and the proportions of strains with the env-Δ236D-phenotype gradually increased with the passage number (from 2.47% to 67.5% at the 103 copies/μL concentration), whereas the proportion of strains with the env-236D-phenotype decreased concomitantly.
Although a gradual increase in the percentage of vaccine-specific sequences was observed using double-probe real-time PCR, the distribution of this mutation in EIAVLN40 was still unclear. The SGA method can avoid the recombination events induced by the polymerase and better detect minor nucleotide changes at the same position; thus, this technique has been performed recently for viral quasispecies analysis [25,26,27]. Thus, we were interested in using SGA to resequence env to obtain more information about env quasispecies during long-term passage. Samples obtained separately from EIAVLN40 (n = 88), EIAVDLV34 (n = 98), EIAVDLV121 (n = 95) and infectious clone pLGFD3-8 (n = 100) were amplified by SGA and aligned. As expected, the 236D mutation was found in the EIAV strains derived from each generation. The proportion of sequences with the 236D mutation best matched the results derived from the bulk PCR (17.0%, 30.6%, 77.9%, and 100%, respectively) (Table 2). Importantly, we detected 236D (15/72) in the EIAVLN40 reservoir, which was not observed using bulk PCR. Additionally, the SGA method exhibited a relatively higher mutation detection rate than bulk PCR. Hence, more variants located at the same position in the EIAV population were observed; for instance, multiple substitutions (E, V, K or Q) instead of a single substitution (E) were detected at position 229 in EIAVLN40 using the SGA method. In particular, in addition to position 236, we found that three other positions, which were separately located at S235R, N237K, and N246K, presented gradual and successful mutation accumulation with the increasing passage number, as presented in Table 2. Further analysis showed that the deletion of amino acid 237D resulted in the loss of a N-glycosylation site (237NNTW240) and that the N237K/N246K mutations would also be involved in N-glycosylation site formation (Figure 2A,B). We also noted that these key mutations (235, 236, 237, and 247) were absent or present in combination; for instance, the frequency of the combined mutations at sites 236 and 237 was 85.8% (109/127). The proportions of mutation that concurred at 235, 236, and 237 was 57.5% (73/127). Meanwhile, the sequences with mutations at 235, 236, 237, and 246 showed a lower frequency than the other sequences, as shown in Figure 2C. Moreover, the proportion of these combined mutations at sites 235–236, 235–236–237 and 235–236–237–246 gradually increased as the virulence of the strains decreased, especially in 235–236 (shown in Figure 2D). Hence, we speculated the existence of a link between the sites 236 and 237. According to these results, we speculated that long-term passage in vitro not only had a significant effect on the env quasispecies compositions but also might contribute to viral attenuation.

3.3. The Relationship between EIAV Attenuation In Vitro and env Gene Quasispecies

To better analyze the relationship between virulence and env quasispecies, we constructed a maximum likelihood tree using the MEGA 7.0 program based on the SGA results. As presented in Figure 3, two obvious major groups of env sequences have evolved, as represented by the env-236D-phenotype (indicated with solid circles) and env-Δ236D-phenotype (indicated with hollow circles). Then, serial passaging events continued to lead to accumulation of the env-Δ236D-phenotype in the whole population, which made the variant distribution even more complex (i.e., the variant distribution of EIAVDLV34). Over the course of 121 generations in vitro, the cluster distribution of related sequences taken from EIAVDLV121 showed a pattern similar to that of EIAVLN40 except for reversal of the predominant phenotype from env-236D to env-Δ236D at the individual population level. We observed that two typical env phenotypes displayed a gradual transition pattern. These results demonstrated that a predominant phenotype transition occurred in the natural virus population under selection, which was consistent with the LTR investigation results [15].

4. Discussion

In contrast to the use of artificial mutagenesis to induce RNA viral attenuation, the EIAV live-attenuated vaccine was developed based on natural evolution through consecutive passaging in vitro. There is growing evidence that analyzing virulence determinants from an evolutionary perspective rather than nucleotide mutations in specific virulence genes is very important, particularly for viruses undergoing consecutive passages in vitro [28,29]. Although a quasispecies is well accepted as a key factor for viral virulence [30], the relationship between these factors remains an open question. In other words, not all viruses will lose their virulence even through serial passage in vitro, such as the CDV Rockborn strain generated in canine macrophage cells [31].
In our previous study, we observed that the EIAV attenuated vaccine also existed as a complex quasispecies [15]. Furthermore, our results suggested that the attenuated vaccine might be a result driven by selection of the dominant EIAV variant at the whole population level. However, whether a similar evolutionary pattern existed in the predominant virulence gene of EIAV was not clear. Here, we used the EIAV attenuation system to investigate the relationship between virulence and EIAV env gene variation. In the present study, we used three different methods to answer this question. Regardless of the method used, some stable mutations located in the V3 and V4 regions of Env were verified. Compared with those of the other methods, SGA showed the highest detection sensitivity for minor nucleotide changes; for instance, the mutations located at positions 235, 236, 237, 246, 248, and 250 could be verified only using the SGA method. In other words, more variants were detected at the same positions. Several key mutations, including S235R, N237K, and N246K, attracted our attention. Not only did the percentages of these mutations parallel with the increased passage number but these mutations (237NNTW240 and N237K/N246K) could also result in loss of an N-glycosylation site or formation of an adjacent N-glycosylation site. We confirmed that this position in combination with three other N-glycosylation site mutations in Env resulted in reversion of virulence and influenced the neutralizing activity [32]. Further analysis showed that inter-linkages existed among these key mutations. Based on these results, we speculated that the observed pattern of Env variation was consistent with the loss of viral virulence.
In particular, we noted that the env-Δ236D-phenotype, which could present as a vaccine-specific mutation, was present in its parental virus reservoir (EIAVLN40) at a lower level. Hence, our observations indicate that the EIAV vaccine strain (env-Δ236D-phenotype) may preexist in the EIAV population as a minority population; indeed, another study from our laboratory detected vaccine-specific LTR sequences in horses infected with virulent strains [15]. According to these results, we proposed a model describing env evolution during the attenuation process, as shown in Figure 4. In this model, the env-236D-phenotype is the predominant hypervirulent EIAV quasispecies that causes typical EIA symptoms in infected horses and induces 100% mortality. The env-Δ236D-phenotype is an avirulent quasispecies that may be a beneficial phenotype for new culture environments (equine macrophages), but the underlying reason is complicated (possibly due to its high fitness or harboring a selectable trait) [6]. To adapt to new culture conditions, the avirulent phenotype selectively accumulated [8,33]. The shape of the phylogenetic tree further reflected this phenomenon in which the env-Δ236D-phenotype gradually became predominant within the whole population but progressively lost its virulence. Correspondingly, the env-Δ236D-phenotype became dominant with ongoing passages, and thus the strain lost its virulence in animals.
In summary, our study characterized EIAV env quasispecies and observed a virulence-correlated parallel in env variation during EIAV passage in vitro. However, this strain in the context of the ensemble of the predominant population progressively lost its virulence in vivo. Hence, our data support the view of quasispecies theory and provide a role for the env quasispecies in EIAV pathogenesis, which partly explains the attenuation-related mechanism of EIAV.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/1999-4915/11/4/380/s1, Figure S1: The morbidity of varied EIAV strains. Figure S2: WebLogo presentation of variability in EIAV env. Figure S3: Quantification of the env-Δ236D-phenotype and env-236D-phenotype sequences located in the V4 region of env by double-probe real-time PCR.

Author Contributions

Conceived and designed the experiments: Y.L. and X.W. Performed the experiments: C.L., X.-F.W., Y.W., and J.C. Analyzed the data: C.L. Wrote the paper: Z.Z., Y.L. and X.W.

Funding

This research was funded by National Natural Science Foundation of China (31672533, 31302066, and 31672578).

Acknowledgments

We thank Han GZ for technical assistance on evolutionary analysis.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Yuan, Q.; Liu, C.; Liang, Z.; Chen, X.; Diao, D.; Kong, X. The comparison of genetic variation in the envelope protein between various immunodeficiency viruses and equine infectious anemia virus. Virol. Sin. 2012, 27, 241–247. [Google Scholar] [CrossRef]
  2. Cook, R.F.; Leroux, C.; Issel, C.J. Equine infectious anemia and equine infectious anemia virus in 2013: A review. Vet. Microbiol. 2013, 167, 181–204. [Google Scholar] [CrossRef]
  3. Ma, J.; Shi, N.; Jiang, C.G.; Lin, Y.Z.; Wang, X.F.; Wang, S.; Lv, X.L.; Zhao, L.P.; Shao, Y.M.; Kong, X.G.; et al. A proviral derivative from a reference attenuated eiav vaccine strain failed to elicit protective immunity. Virology 2011, 410, 96–106. [Google Scholar] [CrossRef] [PubMed]
  4. Wang, X.; Wang, S.; Lin, Y.; Jiang, C.; Ma, J.; Zhao, L.; Lv, X.; Wang, F.; Shen, R.; Kong, X.; et al. Genomic comparison between attenuated chinese equine infectious anemia virus vaccine strains and their parental virulent strains. Arch. Virol. 2011, 156, 353–357. [Google Scholar] [CrossRef]
  5. Wang, X.F.; Lin, Y.Z.; Li, Q.; Liu, Q.; Zhao, W.W.; Du, C.; Chen, J.; Wang, X.; Zhou, J.H. Genetic evolution during the development of an attenuated EIAV vaccine. Retrovirology 2016, 13, 9. [Google Scholar] [CrossRef]
  6. Bingham, R.J.; Dykeman, E.C.; Twarock, R. RNA virus evolution via a quasispecies-based model reveals a drug target with a high barrier to resistance. Viruses 2017, 9, 347. [Google Scholar] [CrossRef] [PubMed]
  7. Colizzi, E.S.; Hogeweg, P. Evolution of functional diversification within quasispecies. Genome Biol. Evol. 2014, 6, 1990–2007. [Google Scholar] [CrossRef]
  8. Park, J.M.; Niestemski, L.R.; Deem, M.W. Quasispecies theory for evolution of modularity. Phys. Rev. E Stat. Nonlin. Soft Matter Phys. 2015, 91, 012714. [Google Scholar] [CrossRef]
  9. Domingo, E.; Sheldon, J.; Perales, C. Viral quasispecies evolution. Microbiol. Mol. Biol. Rev. 2012, 76, 159–216. [Google Scholar] [CrossRef]
  10. Andino, R.; Domingo, E. Viral quasispecies. Virology 2015, 479–480, 46–51. [Google Scholar] [CrossRef] [PubMed]
  11. Domingo, E.; Martin, V.; Perales, C.; Grande-Perez, A.; Garcia-Arriaza, J.; Arias, A. Viruses as quasispecies: Biological implications. Curr. Top. Microbiol. 2006, 299, 51–82. [Google Scholar]
  12. Craigo, J.K.; Ezzelarab, C.; Cook, S.J.; Chong, L.; Horohov, D.; Issel, C.J.; Montelaro, R.C. Envelope determinants of equine lentiviral vaccine protection. PLoS ONE 2013, 8, e66093. [Google Scholar] [CrossRef]
  13. Du, J.S.; Wang, X.F.; Ma, J.; Wang, J.X.; Qin, Y.Y.; Zhu, C.H.; Liu, F.; Shao, Y.M.; Zhou, J.H.; Qiao, W.T.; et al. Structural and biochemical insights into the v/i505t mutation found in the eiav gp45 vaccine strain. Retrovirology 2014, 11. [Google Scholar] [CrossRef]
  14. Duan, L.; Du, J.; Wang, X.; Zhou, J.; Wang, X.; Liu, X. Structural and functional characterization of EIAV gp45 fusion peptide proximal region and asparagine-rich layer. Virology 2016, 491, 64–72. [Google Scholar] [CrossRef]
  15. Wang, X.F.; Liu, Q.; Wang, Y.H.; Wang, S.; Chen, J.; Lin, Y.Z.; Ma, J.; Zhou, J.H.; Wang, X. Characterization of equine infectious anemia virus long terminal repeat quasispecies in vitro and in vivo. J. Virol. 2018, 92. [Google Scholar] [CrossRef]
  16. Lin, Y.Z.; Shen, R.X.; Zhu, Z.Y.; Deng, X.L.; Cao, X.Z.; Wang, X.F.; Ma, J.; Jiang, C.G.; Zhao, L.P.; Lv, X.L.; et al. An attenuated eiav vaccine strain induces significantly different immune responses from its pathogenic parental strain although with similar in vivo replication pattern. Antiviral Res. 2011, 92, 292–304. [Google Scholar] [CrossRef]
  17. Jiang, C.G.; Gao, X.; Ma, J.; Lin, Y.Z.; Wang, X.F.; Zhao, L.P.; Hua, Y.P.; Liu, D.; Zhou, J.H. C-terminal truncation of the transmembrane protein of an attenuated lentiviral vaccine alters its in vitro but not in vivo replication and weakens its potential pathogenicity. Virus Res. 2011, 158, 235–245. [Google Scholar] [CrossRef]
  18. Hayat, M.J.; Higgins, M. Understanding poisson regression. J. Nurs. Educ. 2014, 53, 207–215. [Google Scholar] [CrossRef]
  19. Kumar, S.; Stecher, G.; Tamura, K. Mega7: Molecular evolutionary genetics analysis version 7.0 for bigger datasets. Mol. Biol. Evol. 2016, 33, 1870–1874. [Google Scholar] [CrossRef]
  20. Craigo, J.K.; Montelaro, R.C. Eiav envelope diversity: Shaping viral persistence and encumbering vaccine efficacy. Curr. HIV Res. 2010, 8, 81–86. [Google Scholar] [CrossRef]
  21. Etemad, B.; Ghulam-Smith, M.; Gonzalez, O.; White, L.F.; Sagar, M. Single genome amplification and standard bulk pcr yield hiv-1 envelope products with similar genotypic and phenotypic characteristics. J. Virol. Methods 2015, 214, 46–53. [Google Scholar] [CrossRef] [PubMed]
  22. Wang, X.F.; Wang, S.; Liu, Q.; Lin, Y.Z.; Du, C.; Tang, Y.D.; Na, L.; Wang, X.; Zhou, J.H. A unique evolution of the s2 gene of equine infectious anemia virus in hosts correlated with particular infection statuses. Viruses 2014, 6, 4265–4279. [Google Scholar] [CrossRef]
  23. Ntziora, F.; Paraskevis, D.; Haida, C.; Magiorkinis, E.; Manesis, E.; Papatheodoridis, G.; Manolakopoulos, S.; Beloukas, A.; Chryssoy, S.; Magiorkinis, G.; et al. Quantitative detection of the m204v hepatitis B virus minor variants by amplification refractory mutation system real-time PCR combined with molecular beacon technology. J. Clin. Microbiol. 2009, 47, 2544–2550. [Google Scholar] [CrossRef]
  24. Chen, Q.; Belmonte, I.; Buti, M.; Nieto, L.; Garcia-Cehic, D.; Gregori, J.; Perales, C.; Ordeig, L.; Llorens, M.; Soria, M.E.; et al. New real-time-pcr method to identify single point mutations in hepatitis C virus. World J. Gastroenterol. 2016, 22, 9604–9612. [Google Scholar] [CrossRef] [PubMed]
  25. Salazar-Gonzalez, J.F.; Bailes, E.; Pham, K.T.; Salazar, M.G.; Guffey, M.B.; Keele, B.F.; Derdeyn, C.A.; Farmer, P.; Hunter, E.; Allen, S.; et al. Deciphering human immunodeficiency virus type 1 transmission and early envelope diversification by single-genome amplification and sequencing. J. Virol. 2008, 82, 3952–3970. [Google Scholar] [CrossRef]
  26. Geng, Q.M.; Li, H.P.; Bao, Z.Y.; Liu, Y.J.; Zhuang, D.M.; Li, L.; Liu, S.Y.; Li, J.Y. Indinavir resistance evolution in one human immunodeficiency virus type 1 infected patient revealed by single-genome amplification. Virol. Sin. 2010, 25, 316–328. [Google Scholar] [CrossRef]
  27. Guinoiseau, T.; Moreau, A.; Hohnadel, G.; Ngo-Giang-Huong, N.; Brulard, C.; Vourc’h, P.; Goudeau, A.; Gaudy-Graffin, C. Deep sequencing is an appropriate tool for the selection of unique hepatitis C virus (HCV) variants after single genomic amplification. PLoS ONE 2017, 12, e0174852. [Google Scholar] [CrossRef]
  28. Presloid, J.B.; Novella, I.S. Rna viruses and rnai: Quasispecies implications for viral escape. Viruses 2015, 7, 3226–3240. [Google Scholar] [CrossRef] [PubMed]
  29. Liu, F.; Wu, X.; Li, L.; Zou, Y.; Liu, S.; Wang, Z. Evolutionary characteristics of morbilliviruses during serial passages in vitro: Gradual attenuation of virus virulence. Comp. Immunol. Microbiol. Infect. Dis. 2016, 47, 7–18. [Google Scholar] [CrossRef] [PubMed]
  30. Vignuzzi, M.; Stone, J.K.; Arnold, J.J.; Cameron, C.E.; Andino, R. Quasispecies diversity determines pathogenesis through cooperative interactions in a viral population. Nature 2006, 439, 344–348. [Google Scholar] [CrossRef]
  31. Martella, V.; Blixenkrone-Moller, M.; Elia, G.; Lucente, M.S.; Cirone, F.; Decaro, N.; Nielsen, L.; Banyai, K.; Carmichael, L.E.; Buonavoglia, C. Lights and shades on an historical vaccine canine distemper virus, the rockborn strain. Vaccine 2011, 29, 1222–1227. [Google Scholar] [CrossRef]
  32. Han, X.; Zhang, P.; Yu, W.; Xiang, W.; Li, X. Amino acid mutations in the env gp90 protein that modify n-linked glycosylation of the chinese eiav vaccine strain enhance resistance to neutralizing antibodies. Virus Genes 2016, 52, 814–822. [Google Scholar] [CrossRef] [PubMed]
  33. Glebova, O.; Knyazev, S.; Melnyk, A.; Artyomenko, A.; Khudyakov, Y.; Zelikovsky, A.; Skums, P. Inference of genetic relatedness between viral quasispecies from sequencing data. BMC Genomics 2017, 18, 918. [Google Scholar] [CrossRef]
Figure 1. The evolutionary history was inferred by applying MEGA 7.0 program based on the general time reversible model. Whole genome sequences were obtained by bulk PCR. Bootstrap values above 80% are indicated. EIAVLN40, EIAVDLV34, EIAVDLV62, EIAVDLV92, and EIAVDLV121 are presented with black circles, black triangles, gray circles, blue circles and red circles, respectively.
Figure 1. The evolutionary history was inferred by applying MEGA 7.0 program based on the general time reversible model. Whole genome sequences were obtained by bulk PCR. Bootstrap values above 80% are indicated. EIAVLN40, EIAVDLV34, EIAVDLV62, EIAVDLV92, and EIAVDLV121 are presented with black circles, black triangles, gray circles, blue circles and red circles, respectively.
Viruses 11 00380 g001
Figure 2. Distribution of the consensus mutations located in the V3 and V4 regions and schematic structural representations. (A) Red letters are defined as critical consensus mutation sites displayed within the Env region during the EIAV attenuation process. (B) Schematic figure of the EIAVDLV121 V3 and V4 regions. N-Glycosylation sites are shown as branched lines. (C) Comparison of the proportions of several combined mutations at 235, 236, 237, and 246 for the whole mutants. (D) The percentages of these combined mutations at the key stages of vaccine development (EIAVLN40, EIAVDLV34, and EIAVDLV121).
Figure 2. Distribution of the consensus mutations located in the V3 and V4 regions and schematic structural representations. (A) Red letters are defined as critical consensus mutation sites displayed within the Env region during the EIAV attenuation process. (B) Schematic figure of the EIAVDLV121 V3 and V4 regions. N-Glycosylation sites are shown as branched lines. (C) Comparison of the proportions of several combined mutations at 235, 236, 237, and 246 for the whole mutants. (D) The percentages of these combined mutations at the key stages of vaccine development (EIAVLN40, EIAVDLV34, and EIAVDLV121).
Viruses 11 00380 g002
Figure 3. Construction of a maximum likelihood tree using the MEGA 7.0 program based on the Tamura-Nei model. Bootstrap values above 80% are indicated. EIAVLN40, EIAVDLV34, and EIAVDLV121 are shown in black, green and red, respectively. Solid and hollow circles represent the env-236D-phenotype and env-Δ236D-phenotype sequences, respectively. All samples were collected from three virus strains (EIAVLN40, EIAVDLV34, and EIAVDLV121). A total of 381 sequences were analyzed.
Figure 3. Construction of a maximum likelihood tree using the MEGA 7.0 program based on the Tamura-Nei model. Bootstrap values above 80% are indicated. EIAVLN40, EIAVDLV34, and EIAVDLV121 are shown in black, green and red, respectively. Solid and hollow circles represent the env-236D-phenotype and env-Δ236D-phenotype sequences, respectively. All samples were collected from three virus strains (EIAVLN40, EIAVDLV34, and EIAVDLV121). A total of 381 sequences were analyzed.
Viruses 11 00380 g003
Figure 4. A schematic evolutionary model comparing the EIAV virulent strain to the vaccine strain during long-term passage in vitro. The virulent virus strains are presented with black circles, and the avirulent virus strains are presented with red circles. In this model, the env-Δ236D-phenotype sequences preexisted in EIAVLN40 as a minority population. With serial passages in vitro, sequences containing the env-Δ236D-phenotype gradually became the predominant EIAV quasispecies and generated high replicative fidelity under the new environmental selection conditions. Virulence was simultaneously attenuated as a result, as discussed in the text.
Figure 4. A schematic evolutionary model comparing the EIAV virulent strain to the vaccine strain during long-term passage in vitro. The virulent virus strains are presented with black circles, and the avirulent virus strains are presented with red circles. In this model, the env-Δ236D-phenotype sequences preexisted in EIAVLN40 as a minority population. With serial passages in vitro, sequences containing the env-Δ236D-phenotype gradually became the predominant EIAV quasispecies and generated high replicative fidelity under the new environmental selection conditions. Virulence was simultaneously attenuated as a result, as discussed in the text.
Viruses 11 00380 g004
Table 1. Amino acid differences in Env V4 among different Equine Infectious Anima Virus (EIAV) strains obtained using bulk PCR.
Table 1. Amino acid differences in Env V4 among different Equine Infectious Anima Virus (EIAV) strains obtained using bulk PCR.
22222222222222222222222222222222Mutation Rate
(236D/-)
Morbidity
(Horses)
22222222233333333334444444444555
12345678901234567890123456789012
EIAVLN40N29A29I29T29I29L29V29P29E29V29E29M28Q29Q28S29D29N28N29T29W29I29P29K29R29C27N29E29T29W29A29R29V290.00%100%
I1 R1 S1 Y2
EIAVDLV34S9A16I15T15I16L16V16P16E16V10E16M16Q16Q14S16D14N16N16T16W16I16P16K16G8C16N14E13T16W16A16R12V166.25%100%
N7 V1A1 I6 K1 N1 R8 K2K3 K4
R1 -1
EIAVDLV62N9A16I16T16I16L16V16P16E16V9E16M16Q16Q14S16D9N16N16T16W16I16P16K16R11C16N13E13T16W16A16R12V1618.5%100%
S7 I7 R2 N4 G5 K3K2 K4
-3 G1
EIAVDLV92N11A19I19T16I19L19V19P19E19V18E19M19Q19Q17S19D3N13N18T16W19I18P19K19R17C19N14E18T18W19A19R12V1942.1%9.1%
S8 A3 I1 R1 N8K5D1A3 V1 G2 K5K1A1 K7
K1 -8R1
EIAVDLV121N11A17I18T17I18L18V18P18E17V17E17M18Q18Q18R15-18K17N18T18W18I17P18K18R17C18K9E10T18W18A18K13V17100%0
S7T1 A1 K1I1K1 S3 S1 V1 G1 N8K8 R5I1
E1
The additional amino acid D, which consistently appeared in virulent EIAV strains, is marked with dark grey, and N-linked glycosylation sites are marked with light grey.
Table 2. Amino acid differences in Env V4 among different EIAV strains detected using single genome amplification (SGA).
Table 2. Amino acid differences in Env V4 among different EIAV strains detected using single genome amplification (SGA).
22222222222222222222222222222222Mutation Rate
(236D/-)
22222222233333333334444444444555
12345678901234567890123456789012
EIAVLN40N86A88I88T88I88L88V88P88E85V87E88M88Q87Q88S80D72N71N88T88W88I88P88K87R80C88N86E87T87W88A88R81V8817.0%
S2 V1I1 H1 R8-15K14 N1G8 K2K1I1 K7
K1 N1H1
Q1
EIAVDLV34N51A97I98T98I98L97V98P97E94V60E96M98Q98Q82S86D63N75N98T98W95I95P97K96R57C97N85E79T96W94A94R75V9430.6%
S47T1 V1 L1K3I38K2 K8R10-30K22 S1S2L1Q1G41W1K10K17A1C3C2K20I4
D1 R8K1N5Y1 C2C1 P1 S3G2I1R1S1T3
N1 T1
EIAVDLV121N76A95I95T92I95L95V94P95E92I10E93M94Q94Q91R54-74K69N95T95W95I95P95K95R68C95N74E83T94W95A93K55V9577.9%
S18 A3 I1 K1V84K1V1H1R3S40D19N23 G25 K21K11A1 T2R40
K1 D1L1V1 K1N1N2S2 K2 G1
Q1 E1
pLGFD3-8N100A100I100T100I100L100V100P100E99V97E99M100Q98Q100S87-100K100N100T97W98I99P100K100G88C99K99K99T100W99A100R97V100100%
Q1I3K1 H1 R13 I3C2V1 R12G1N1E1 C1 K2
R1 M1
The additional amino acid D, which consistently appeared in virulent EIAV strains, is marked with dark grey, and N-linked glycosylation sites are marked with light grey.

Share and Cite

MDPI and ACS Style

Liu, C.; Wang, X.-F.; Wang, Y.; Chen, J.; Zhong, Z.; Lin, Y.; Wang, X. Characterization of EIAV env Quasispecies during Long-Term Passage In Vitro: Gradual Loss of Pathogenicity. Viruses 2019, 11, 380. https://0-doi-org.brum.beds.ac.uk/10.3390/v11040380

AMA Style

Liu C, Wang X-F, Wang Y, Chen J, Zhong Z, Lin Y, Wang X. Characterization of EIAV env Quasispecies during Long-Term Passage In Vitro: Gradual Loss of Pathogenicity. Viruses. 2019; 11(4):380. https://0-doi-org.brum.beds.ac.uk/10.3390/v11040380

Chicago/Turabian Style

Liu, Cong, Xue-Feng Wang, Yan Wang, Jie Chen, Zhaohua Zhong, Yuezhi Lin, and Xiaojun Wang. 2019. "Characterization of EIAV env Quasispecies during Long-Term Passage In Vitro: Gradual Loss of Pathogenicity" Viruses 11, no. 4: 380. https://0-doi-org.brum.beds.ac.uk/10.3390/v11040380

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop