Next Article in Journal
Swine Influenza Virus Induces RIPK1/DRP1-Mediated Interleukin-1 Beta Production
Previous Article in Journal
Neutralizing Epitopes and Residues Mediating the Potential Antigenic Drift of the Hemagglutinin-Esterase Protein of Influenza C Virus
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Viruses
Volume 10
Issue 8
10.3390/v10080418
viruses-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Xenotropic Mouse Gammaretroviruses Isolated from Pre-Leukemic Tissues Include a Recombinant

by
Devinka Bamunusinghe
†,
Matthew Skorski
,
Alicia Buckler-White
and
Christine A. Kozak
*
Laboratory of Molecular Microbiology, National Institute of Allergy and Infectious Diseases, Bethesda, MD 20892-0460, USA
*
Author to whom correspondence should be addressed.
†
Current address: Qiagen, Germantown, MD 20874, USA.
Viruses 2018, 10(8), 418; https://doi.org/10.3390/v10080418
Submission received: 15 June 2018 / Revised: 28 July 2018 / Accepted: 3 August 2018 / Published: 9 August 2018
(This article belongs to the Section Animal Viruses)
Browse Figures
Versions Notes

Abstract

:
Naturally-occurring lymphomagenesis is induced by mouse leukemia viruses (MLVs) carried as endogenous retroviruses (ERVs). Replicating the ecotropic MLVs recombines with polytropic (P-ERVs) and xenotropic ERVs (X-ERVs) to generate pathogenic viruses with an altered host range. While most recovered nonecotropic recombinants have a polytropic host range, the X-MLVs are also present in the pre-leukemic tissues. We analyzed two such isolates from the AKR mice to identify their ERV progenitors and to look for evidence of recombination. AKR40 resembles the active X-ERV Bxv1, while AKR6 has a Bxv1-like backbone with substitutions that alter the long terminal repeat (LTR) enhancer and the envelope (env). AKR6 has a modified xenotropic host range, and its Env residue changes all lie outside of the domain that governs the receptor choice. The AKR6 segment spanning the two substitutions, but not the entire AKR6 env-LTR, exists as an ERV, termed Xmv67, in AKR, but not in the C57BL/6 mice. This suggests that AKR6 is the product of one, not two, recombination events. Xmv67 originated in the Asian mice. These data indicate that the recombinant X-MLVs that can be generated during lymphomagenesis, describe a novel X-ERV subtype found in the AKR genome, but not in the C57BL/6 reference genome, and identify residues in the envelope C-terminus that may influence the host range.

1. Introduction

Laboratory mice carry mouse leukemia viruses (MLVs) of three host range subgroups (reviewed in [1]). The MLVs with an ecotropic host range (E-MLVs) infect only the rodent cells, while the various xenotropic and polytropic MLVs (X-, P-MLVs, and collectively, X/P-MLVs) infect different subsets of the mouse taxa and other mammalian species [2,3]. The E-MLVs use the CAT-1 receptor [4], and all X/P-MLVs use the XPR1 receptor [5,6,7]. There are functionally distinct XPR1 variants in the genus Mus, but the majority of mouse strains carry Xpr1n, which allows for the infection of P-MLVs but not X-MLVs [8].
All three of the host range subgroups of the MLVs exist as infectious viruses and as endogenous retroviruses (ERVs), which are germline DNA copies acquired during past virus infections [9,10]. The specific ERVs can be present or absent in the mouse strains, and can be shared by the strains sharing a common ancestry. Many of the E-MLV ERVs (E-ERVs, termed Emvs) [9] and some of the X-ERVs (Xmvs) are capable of producing an infectious virus [1]. There are two subclasses of P-ERVs, the polytropic murine viruses (Pmvs) and modified polytropic murine viruses (Mpmvs) [11]. None of the infectious P-MLVs derive directly from Pmvs or Mpmvs, although these P-ERVs can contribute to the generation of intersubgroup recombinant viruses that have the distinctive P-MLV host range [12,13,14,15].
Recombinant P-MLVs are routinely generated de novo in laboratory mouse strains carrying replicating E-MLVs, and the appearance of these recombinants is linked to virus-induced lymphomas and leukemias [16,17,18]. High virus, high-leukemic mouse strains, like AKR, HRS, and C58, carry active Emvs, and the E-MLVs they generate can recombine with the products of nonecotropic ERVs to produce P-MLVs with an increased virulence and altered host range [18,19].
Our recent analyses of nonecotropic MLVs from leukemic and pre-leukemic tissues focused on the E-MLV derived P-MLVs [15], but the P-ERV substitutions have also been found in the lymphomagenic viruses from the mice infected with amphotropic MLVs [20], and there are multiple examples of intersubtype recombination involving Xmvs. For example, the melanoma associated MelARV MLV is a recombinant of the Emv2 E-MLV and Xmv45/XmvIV1 [21], the XMRV X-MLV is a recombinant of two nonecotropic ERV sequences [22], and some P-MLV recombinants can incorporate Xmv- and Pmv-derived sequences [15]. The analyses of the other infectious and endogenous MLVs of the inbred and wild mice suggest that some are recombinants that include X-MLV-related sequences [23,24].
Studies on the virological events that lead to naturally occurring leukemogenesis have largely been conducted in AKR mice [16,17,18]. Here, we sequenced two X-MLVs isolated from a pre-leukemic AKR mouse thymus, to identify their ERV progenitors and to determine if they had been subjected to recombination. One of these viruses, AKR6, is a recombinant derived from the Bxv1 Xmv, with substitutions derived from a distinctive ERV found in AKR, but not in the C57BL/6 (B6) reference genome. We demonstrated that this ERV subtype was acquired by Asian mice, and identified multiple envelope (Env) replacement mutations that may be responsible for the altered AKR6 xenotropic host range.

2. Materials and Methods

2.1. Viruses and Sequencing

The AKR6 and AKR40 X-MLVs were obtained from J. Hartley (National Institute of Allergy and Infectious Diseases, Bethesda, MD, USA) and had been isolated from the thymus tissue of two- and six-month old AKR mice, respectively [25]. The ferret MA139 cells, obtained from J. Hartley, were infected with AKR6, and were passaged several times. The viral genome was sequenced from the DNA of the infected ferret cells, using primers developed to sequence a set of P-MLVs [15].
The DNAs from various inbred mouse strains and wild-derived or wild-trapped mice were obtained from The Jackson Laboratory (Bar Harbor, ME, USA), the RIKEN BioResource Center through the National BioResource Project of the MEXT/AMED (Ibaraki, Japan), M. Potter (NCI, Bethesda, MD, USA), S. Rasheed (University of Southern California, Los Angeles, CA, USA), R. Abe (Naval Med. Res. Inst., Bethesda, MD, USA), S. Chattopadhyay, and H. Morse III (NIAID, Bethesda, MD, USA), or were isolated from the mice obtained from The Jackson Laboratory or from mice maintained in our laboratory (Table S1). To identify the AKR mouse ERVs related to AKR6, we used forward and reverse primers designed to amplify the AKR6 env and/or the long terminal repeat (LTR), as follows: F1, GAGCAGCATGGAAGGTCC; F2, CGGAATCCTCTATTGGAC; F3, CACTGCCCTAGTGACCACACG; R1, TGCAAACAGCAAAAGGCTTTATTGG; R2, AGTCTAACCCTCTTCGGT. Four primers were designed from AKR6; and F2 was based on the F1–R2 PCR product.
The AKR6 and AKR40 sequences were compared with previously described ERVs and MLVs. The highlighter plots were constructed using appropriate ERVs and MLVs (www.hiv.lanl.gov/content/sequence/HIGHLIGHT/highlighter_top.html) [26].
The sequences of AKR6 and AKR40 were deposited in GenBank under accession nos. MH450109-10. The 5′ 1.3 kb end of AKR6 env was previously sequenced (DQ199948) [27]. Three partial AKR6-related ERVs amplified from the AKR mouse DNA were deposited under MH443069-71.

2.2. Phylogenetic Trees

Phylogenetic trees were constructed in MEGA7 [28], using the maximum likelihood method based on the Kimura 2-parameter model with 500 bootstrap replicates [29]. Four trees were based on the two segments of AKR6 that diverge from Bxv1, the 5′ 357 bp of env, and the reverse transcriptase (RT) domain of pol. The alignments were generated by MUSCLE [30] and were corrected manually. The initial trees for the heuristic search were obtained automatically by applying neighbor-join and BioNJ algorithms to a matrix of pairwise distances estimated using the Maximum Composite Likelihood (MCL) approach, and then selecting the topology with a superior log likelihood value. A discrete Gamma distribution was used to model the evolutionary rate differences among the sites. All of the positions with less than a 95% site coverage were eliminated. That is, fewer than 5% alignment gaps, missing data, and ambiguous bases were allowed at any position. The trees included up to 14 Xmvs and representative Pmvs and Mpmvs from the B6 reference genome [10,31], as well as the following: PreXMRV-1 (FR871849), PreXMRV-2 (FR871850), NZB-9-1 (K02730, EU035300), CasE#1 (KU324802), Cz524 (KU324804), Kyushu (KU324806), CAST-X (KU324803), and a set of Y chromosome-linked ERVs that, together with the chromosome 5 Xmv45/XmvIV1, comprise the distinct XmvIV subset of Xmvs [21,23,32].

2.3. Inbred Strain Distribution and Subspecies Origins of AKR6-Related ERVs

We used the BLAST/BLAT tool [33] on the Ensembl website (https://useast.ensembl.org/Mus_musculus/Tools/Blast?db=core) [34] to search for the AKR6 related ERV sequences in the B6 mouse reference genome, and in the draft sequences of the AKR and 10 other mouse genomes generated by the Mouse Genomes Group at the Wellcome Trust Sanger Institute (https://www.sanger.ac.uk/science/data/mouse-genomes-project).
The DNA from 44 laboratory strains and from 43 wild-derived and wild-caught mice of the M. musculus subspecies were screened for the AKR6-related X-ERVs, with primers F2 and R1, and for a specific AKR6-related ERV found in AKR mice, Xmv67, using a primer from the viral LTR [35] and the flanking sequence primer, GGCTCTTGTCTGGCTTTAACC.
The Mouse Phylogeny Viewer (MPV) (http://msub.csbio.unc.edu) [36] was used to determine the subspecies origins of the AKR6-related ERV found in AKR strain of mice. Because the NCBI37/mm9 B6 assembly was used to create the MPV database, we first determined the B6 coordinates for this ERV using cellular sequences flanking the ERV in the AKR genome to do a BLAT search of NCBI37/mm9 with the UCSC Genome browser (http://genome.ucsc.edu/).

3. Results and Discussion

3.1. Sequence Analysis of Two X-MLVs

The AKR6 and AKR40 MLV isolates both show a xenotropic host range, in that they can infect mink cells, but fail to infect laboratory mouse cells [37]. An analysis of their full-length genomes shows that, unlike the P-MLVs isolated from these tissues, neither of the X-MLVs incorporate the E-MLV-derived sequences. Highlighter plots show the genome-wide distribution of the nucleotide differences between the two sequenced X-MLVs and their closest homolog, the Bxv1 Xmv (Figure 1a). Bxv1 is a full-length non-defective ERV that is normally quiescent, but can be induced to produce an infectious virus [38]; Bxv1 was first identified in B6, but is also present in the AKR genome [10,39]. The AKR40 genome is 99.9% identical to Bxv1. While AKR6 shows a comparably high homology to Bxv1 throughout most of its genome, there are two segments of 129 bp and 1.5 kb that show a reduced identity to Bxv1, of 94.3% and 94.5%, respectively (Figure 1b); the intervening 191 bp are 100% identical to Bxv1. The 129 bp replacement is in the LTR, and the 1.5 kb replacement is in env.
Four phylogenetic trees were constructed to compare the AKR6 to the previously described X/P-MLVs and ERVs, based on RTpol (Figure 2a), the AKR6 1.5 kb and 129 bp replacements (Figure 2b,c), as well as the env segment upstream of the AKR6 replacement (Figure 2d). While the B6 polytropic Pmvs and Mpmvs cluster together in all four trees, the B6 Xmvs show much more sequence variation, as previously shown for their full length genomes [31]. In all four trees, AKR40 clusters with Bxv1. AKR6 also clusters with Bxv1 in the RTpol and 5’env trees (Figure 2a,d), but the AKR6 env substitution groups with the XmvIV subset of the Xmv ERVs [23,32] (Figure 2b,c). BLAST searches confirm that the best sequence matches for the AKR6 env replacement in the reference genome are among the XmvIV ERVs, with ~98% identity. The XmvIV ERVs do not produce an infectious virus, but related segments have been found in other recombinant viruses from inbred mice, including the melanoma-associated MLV, MelARV, and several P-MLVs as well as wild mouse MLVs [15,21,23].
The 129 bp LTR substitution shows closest identity to two solo LTRs in the B6 reference genome (Chr15:74933718-74934248; 17:35000213-35000742) with a 99.4% identity to both (Figure 2c). The high sequence identity is, however, restricted to this 129 bp segment, as the identity at the 3′ and 5′ ends of these LTRs is reduced to 96.7% and 94.3%, respectively (Figure 1b).
Both replacements in the AKR6 genome alter the functionally important viral segments. The smaller 129 bp replacement spans the U3 LTR enhancer segment containing the four transcription factor binding sites, which form the retroviral enhancer framework, namely, LVb, core, NF1, and GRE [40]. The generation of virulent P-MLVs is associated with alterations in two enhancer elements, core and NF1, and these changes in the P-MLVs can result from mutation or recombination [15]. For AKR6, however, the enhancer replacement introduces the more benign variant of the core site (TGTGGTCAA vs. TGTGGTCGA) [41], without changing the other sites.
The larger 1.5 kb AKR6 replacement in env begins in the middle of the receptor binding domain of env (RBDenv), extending from the end of the VRC variable region to the beginning of the R peptide in the cytoplasmic tail of the transmembrane subunit of env (TMenv) (Figure 1b). Our previous analysis of 16 P-MLVs showed that these viruses, derived from the recombination of E-MLVs and P-ERVs, showed 16 different recombination patterns that collectively covered the whole viral genome. All 16 shared two and only two common replacements, both in env, at the 5′ end of the surface (SU), and the 3′ cytoplasmic tail R peptide (Figure 1b) [15]. These two segments show a substantial sequence divergence between the E-MLV backbone of these P-MLVs and their P-ERV-derived substitutions. The replicating viruses carry either both P-MLV segments or both E-MLV types, but never a combination, suggesting a sequence-based constraint. Compared to these P-MLVs, AKR6 shows a reciprocal replacement pattern in that the two env segments that are replaced in all of the P-MLVs are retained by AKR6 from its Bxv1 progenitor, while the intervening sequence is replaced (Figure 1b). The R peptide has been implicated in the regulation of the Env ectodomain [42], so the fact that these two env segments are either coordinately retained (in AKR6) or replaced (in P-MLVs), suggests restrictions on the pairing of some variants of these two env segments. The maintenance of both Bxv1-like env segments in the AKR6 recombinant suggests that the comparable segments in its ERV progenitor either produce nonfunctional proteins in either or both of these segments, or that these segments encode variants that are not compatible with either of the Bxv1 env segments.

3.2. ERV Progenitors of the AKR6 Replacements

We used in silico searches and PCR to look for possible progenitors of AKR6 in the AKR mouse genome (Figure 1b). Firstly, we screened the available AKR genome draft sequence using BLAST/BLAT. The only good match for the LTR enhancer replacement, at 17:34443814-34444343, is a solo LTR that is orthologous to one also found in the B6 genome, as described above. Like its B6 ortholog, this AKR LTR shows 100% homology to the enhancer region of the AKR6 LTR, but has only 96.7% and 94.3% identity to the 5′ and 3′ ends of this LTR.
We also identified a single sequence in the AKR draft genome related to the env replacement (11:98208845-98210232). This partial ERV has no ortholog in B6 and was not identified among a set of AKR Xmvs [10]. Therefore, we gave it the next consecutive Xmv number, Xmv67. In this AKR assembly, Xmv67 is flanked by large sequencing gaps and contains two internal gaps, both of which are associated with duplications (Figure 1b). The 3′ 721 bp of this sequence is identical to the corresponding AKR6 env segment, but the sequence at its 5′ end is more divergent, with only a 93.7% identity to AKR6. This ERV segment encodes an env ORF, but is distinct from the other X-ERVs (Figure 2d).
Next, we coupled AKR6- and Xmv-specific primers to amplify the extended segments containing the AKR6-related env and/or LTR replacement segment from the AKR genome. We consistently identified the same overlapping amplicons (Figure 1b). F1–R2 contains the 5′ end of the env replacement segment, and F3-R1 contains much of the AKR6 TMenv, the LTR substitution, and the intervening sequence, which is shared with Bxv1. The 5′ end of F1–R2, and the 3′ end of F3–R1, are not found in the AKR6 viral genome (Figure 1b), nor are they closely related to the B6 Xmvs (shown for 5′ env, Figure 2d). A primer designed from the divergent segment at the 5′ end of F1–R2 was used to generate an extended amplicon, F2–R1, produced by AKR but not B6. The segment from 357-2190 of F2–R1 that is identical to AKR6, includes the two segments, in env and LTR, that distinguish AKR6 from Bxv1, as well as the intervening Bxv1-like segment. The presence of all three of these segments in this single ERV-derived PCR product indicates that the AKR6 substitution spanning the env and LTR exists as an ERV in AKR, but is not found in the B6 reference genome. Thus, AKR6 is a recombinant likely derived from one, not two, recombination events, and furthermore, the sequence differences between AKR6 and Xmv67 at the 5′ and 3′ ends, show that AKR6 does not exist as an intact ERV in AKR.

3.3. Acquisition of AKR6-Related ERVs in M. musculus Subspecies

The classical inbred strains of the laboratory mouse are mosaics of 3 M. musculus subspecies: domesticus, musculus, and castaneus [43]. Xmvs and the single autosomal copy of the XmvIV subgroup can be traced to their Asian M. m. musculus and M. m. castaneus progenitors [23,35]. To determine if the Xmv67 AKR6 progenitor has a similar wild mouse origin or is a recent acquisition, we designed a primer from the flanking nonviral sequences downstream of Xmv67, and typed a set of mouse DNA for a cell-virus junction product specific to Xmv67. The strain distribution of this ERV (Figure 3) does not correspond to that of any of the previously described ERVs identified in AKR [10]. The inbred strains were also typed for the AKR6-related F2–R1 internal viral fragment (Figure 3); this sequence is found in ten strains that lack the Xmv67 insertion, indicating that F2–R1 identifies a distinctive subset of Xmvs found in some inbred strains.
Of the 44 strains typed for Xmv67, 39 are included in the MPV database, which used single nucleotide polymorphisms (SNPs) to infer the subspecific origins of DNA segments along each chromosome for 100 inbred strains [44]. Figure 4a shows a 10 Mb segment containing the Xmv67 insertion site. This site lies in a genomic segment of AKR that originated in M. m. musculus, and this site is also clearly M. m. musculus-derived in seven of the eight other Xmv67-positive strains, but in none of the negative strains. This ERV is also detected in three of the four wild-derived Japanese mice that have M. m. musculus-derived sequences at this locus (Figure 4a).
Figure 4a shows two discrepancies, the inbred strain, YBR, which places its Xmv67 locus in a M. m. domesticus-derived segment, and the wild-derived strain, MOLF, which lacks Xmv67 but is M. m. musculus in this chromosomal segment. Because small segments of intersubspecies recombination can be missed, we examined the SNP polymorphisms around this integration site, and showed that YBR but not MOLF shares the SNPs flanking the ERV insertion site that are exclusive to the Xmv67-positive mice (Figure 4b).
The Asian mouse origin of Xmv67 was confirmed by PCR-typing of the DNA from the M. musculus subspecies (Figure 4c). Xmv67 was identified in 10 of the 13 M. m. molossinus, the Japanese mouse that is a natural hybrid of M. m. musculus and M. m. castaneus, in 1 of the 16 M. m. musculus; and in none of the 7 M. m. castaneus and 8 M. m. domesticus. These data indicate that Xmv67, like all other X-ERVs, originated in Asian mice, but unlike the other Xmvs, some of which are found in castaneus as well as molossinus, Xmv67 may derive from the musculus ancestor of molossinus.

3.4. Env Replacement Mutations Linked to the Modified AKR6 Host Range

The AKR6 X-MLV was initially described as a representative X-MLV in host range [37], but this virus shows differences in entry-related phenotypes compared with other X-MLVs (Figure 5a). Firstly, AKR6 and one other X/P-MLV, CasE#1, are poorly infectious for an XPR1 receptor with the double mutation K500E/F584I [35]. Secondly, the AKR6, CasE#1, and XMRV MLVs uniquely fail to infect the Chinese hamster cells in which the glycosylation is inhibited [45,46]. These previous studies did not attempt to identify the Env residues responsible for these host range differences.
AKR6 is identical to Bxv1 in the N-terminus of Env containing VRA, which determines the receptor choice [47]. The AKR6 substitution begins in VRC and spans the rest of SUenv including the variant segments VRB and the proline rich region (PRR). A comparison of the env sequences of the relevant X-MLVs identified shared replacements at 18 sites for the viruses that differentially infect XPR1 K500E/F584I, 12 of which also distinguished viruses that vary in sensitivity to glycosylation inhibition in hamster cells (Figure 5b). Four of these sites map to VRB and three adjacent sites are in PRR; previous studies on MLVs and other gammaretroviruses implicate both VRB and PRR in receptor interactions [47,48]. The remaining replacement mutations are scattered, but one site implicated in the mutant receptor usage, T217A, is at one of two adjacent sites in the RBD C-terminus, previously shown to permit some X/P-MLVs to infect human cells [49]. This analysis identifies multiple sites in the C-terminal end of the RBD that may be responsible for the altered receptor properties of AKR6 and other MLVs.

4. Conclusions

We show that the recombinant MLVs generated during pathogenesis can include X-MLVs, and that intrasubtype recombination can modify their phenotype, in this case viral tropism. The substitution in AKR6 that extends from the VRC of RBDenv through the LTR enhancer, derives from an Xmv ERV variant found in AKR, but not B6 mice. This ERV substitution is most closely related to, but not identical to, the XmvIV subgroup of Xmvs and, like all other Xmvs, is traceable to the ERVs present in Japanese mice, known contributors to the fancy mouse colonies that gave rise to the classical laboratory strains.
The inbred strains differ in their MLV ERV complements, largely due to the different insertions of the same set of MLV subtypes [10,31]. The AKR6 env-LTR ERV progenitor is a novel X-ERV subtype present in at least one copy in AKR, but most other laboratory strains lack this ERV. The only closely related ERVs retained in the B6 reference genome are represented by two solo LTRs. The 5’ end of the AKR Xmv67 Env is distinctive, as are segments of the AKR6-related LTRs, and that these sequences are not found in replicating viruses suggests they may not be functional.
The presence of the Xmv-derived viruses in the AKR pre-leukemic thymus is not surprising. While the Bxv1 Xmv is normally quiescent, its expression can be induced by immunological stimulation [50], and the virus spread of X/P-MLVs is not strictly governed by the receptor compatibility. The greater abundance of recombinant P-MLVs in preleukemic tissues may be due to the fact that these viruses have three alternative methods of spread, they can use the AKR XPR1 receptor, they are packaged preferentially into E-MLV virions that rely on the CAT-1 receptor for entry [51], and they can also use the CAT-1 receptor if soluble E-MLV SU-Env is present (transactivation) [52]. X-MLVs, on the other hand, cannot use the AKR XPR1 receptor and are not efficiently transactivated by the E-MLV SU-Env [52], but the X-MLV genomes can be co-packaged with other MLV subtypes, resulting in P-MLVs that incorporate X-MLV sequences [15] or, as shown there, can also produce recombinant X-MLVs like AKR6. The routine generation of recombinant X-MLVs and P-MLVs that include Xmv-derived segments indicates that these recombinants are a common event that may have biological consequences during virus-induced leukemogenesis.

Supplementary Materials

The following are available online at https://www.mdpi.com/1999-4915/10/8/418/s1, Table S1: Wild-derived or wild-caught M. musculus mice typed for Xmv67.

Author Contributions

Conceptualization, C.A.K.; methodology, D.B., M.S., and A.B.-W.; formal analysis, D.B., M.S., A.B.-W., and C.A.K.; resources, C.A.K.; writing (original draft preparation, review, and editing), C.A.K.; supervision, C.A.K.; funding acquisition, C.A.K.

Funding

This research was funded by the Intramural Research Program of the National Institute of Allergy and Infectious Diseases, under grant number AI000300-36 to C.A.K.

Acknowledgments

We thank Janet Hartley (NIAID, Bethesda, MD, USA) for providing the virus isolates.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; and in the decision to publish the results.

References

  1. Kozak, C.A. Origins of the endogenous and infectious laboratory mouse gammaretroviruses. Viruses 2015, 7, 1–26. [Google Scholar] [CrossRef] [PubMed]
  2. Oie, H.K.; Russell, E.K.; Dotson, J.H.; Rhoads, J.M.; Gazdar, A.F. Host range properties of murine xenotropic and ecotropic type-C viruses. J. Natl. Cancer Inst. 1976, 56, 423–426. [Google Scholar] [CrossRef] [PubMed]
  3. Yan, Y.; Liu, Q.; Wollenberg, K.; Martin, C.; Buckler-White, A.; Kozak, C.A. Evolution of functional and sequence variants of the mammalian XPR1 receptor for mouse xenotropic gammaretroviruses and the human-derived XMRV. J. Virol. 2010, 84, 11970–11980. [Google Scholar] [CrossRef] [PubMed]
  4. Albritton, L.M.; Tseng, L.; Scadden, D.; Cunningham, J.M. A putative murine ecotropic retrovirus receptor gene encodes a multiple membrane-spanning protein and confers susceptibility to virus infection. Cell 1989, 57, 659–666. [Google Scholar] [CrossRef]
  5. Tailor, C.S.; Nouri, A.; Lee, C.G.; Kozak, C.; Kabat, D. Cloning and characterization of a cell surface receptor for xenotropic and polytropic murine leukemia viruses. Proc. Natl. Acad. Sci. USA 1999, 96, 927–932. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Battini, J.L.; Rasko, J.E.; Miller, A.D. A human cell-surface receptor for xenotropic and polytropic murine leukemia viruses: Possible role in G protein-coupled signal transduction. Proc. Natl. Acad. Sci. USA 1999, 96, 1385–1390. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Yang, Y.-L.; Guo, L.; Xu, S.; Holland, C.A.; Kitamura, T.; Hunter, K.; Cunningham, J.M. Receptors for polytropic and xenotropic mouse leukaemia viruses encoded by a single gene at Rmc1. Nat. Genet. 1999, 21, 216–219. [Google Scholar] [CrossRef] [PubMed]
  8. Baliji, S.; Liu, Q.; Kozak, C.A. Common inbred strains of the laboratory mouse that are susceptible to infection by mouse xenotropic gammaretroviruses and the human-derived retrovirus XMRV. J. Virol. 2010, 84, 12841–12849. [Google Scholar] [CrossRef] [PubMed]
  9. Jenkins, N.A.; Copeland, N.G.; Taylor, B.A.; Lee, B.K. Organization, distribution, and stability of endogenous ecotropic murine leukemia virus DNA sequences in chromosomes of Mus musculus. J. Virol. 1982, 43, 26–36. [Google Scholar] [PubMed]
  10. Frankel, W.N.; Stoye, J.P.; Taylor, B.A.; Coffin, J.M. A linkage map of endogenous murine leukemia proviruses. Genetics 1990, 124, 221–236. [Google Scholar] [PubMed]
  11. Stoye, J.P.; Coffin, J.M. The four classes of endogenous murine leukemia virus: Structural relationships and potential for recombination. J. Virol. 1987, 61, 2659–2669. [Google Scholar] [PubMed]
  12. Lung, M.L.; Hartley, J.W.; Rowe, W.P.; Hopkins, N.H. Large RNase T1-resistant oligonucleotides encoding p15E and the U3 region of the long terminal repeat distinguish two biological classes of mink cell focus-forming type C viruses of inbred mice. J. Virol. 1983, 45, 275–290. [Google Scholar] [PubMed]
  13. Chattopadhyay, S.K.; Cloyd, M.W.; Linemeyer, D.L.; Lander, M.R.; Rands, E.; Lowy, D.R. Cellular origin and role of mink cell focus-forming viruses in murine thymic lymphomas. Nature 1982, 295, 25–31. [Google Scholar] [CrossRef] [PubMed]
  14. Khan, A.S.; Rowe, W.P.; Martin, M.A. Cloning of endogenous murine leukemia virus-related sequences from chromosomal DNA of BALB/c and AKR/J mice: Identification of an env progenitor of AKR-247 mink cell focus-forming proviral DNA. J. Virol. 1982, 44, 625–636. [Google Scholar] [PubMed]
  15. Bamunusinghe, D.; Liu, Q.; Plishka, R.; Dolan, M.A.; Skorski, M.; Oler, A.J.; Yedavalli, V.R.K.; Buckler-White, A.; Hartley, J.W.; Kozak, C.A. Recombinant origins of pathogenic and nonpathogenic mouse gammaretroviruses with polytropic host range. J. Virol. 2017. [Google Scholar] [CrossRef] [PubMed]
  16. Hartley, J.W.; Wolford, N.K.; Old, L.J.; Rowe, W.P. New class of murine leukemia-virus associated with development of spontaneous lymphomas. Proc. Natl. Acad. Sci. USA 1977, 74, 789–792. [Google Scholar] [CrossRef] [PubMed]
  17. Herr, W.; Gilbert, W. Somatically acquired recombinant murine leukemia proviruses in thymic leukemias of AKR/J mice. J. Virol. 1983, 46, 70–82. [Google Scholar] [PubMed]
  18. Stoye, J.P.; Moroni, C.; Coffin, J.M. Virological events leading to spontaneous AKR thymomas. J. Virol. 1991, 65, 1273–1285. [Google Scholar] [PubMed]
  19. Cloyd, M.W.; Hartley, J.W.; Rowe, W.P. Lymphomagenicity of recombinant mink cell focus-inducing murine leukemia viruses. J. Exp. Med. 1980, 151, 542–552. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  20. Rasheed, S.; Pal, B.K.; Gardner, M.B. Characterization of a highly oncogenic murine leukemia virus from wild mice. Int. J. Cancer 1982, 29, 345–350. [Google Scholar] [CrossRef] [PubMed]
  21. Pothlichet, J.; Mangeney, M.; Heidmann, T. Mobility and integration sites of a murine C57BL/6 melanoma endogenous retrovirus involved in tumor progression in vivo. Int. J. Cancer 2006, 119, 1869–1877. [Google Scholar] [CrossRef] [PubMed]
  22. Paprotka, T.; Delviks-Frankenberry, K.A.; Cingoz, O.; Martinez, A.; Kung, H.J.; Tepper, C.G.; Hu, W.S.; Fivash, M.J., Jr.; Coffin, J.M.; Pathak, V.K. Recombinant origin of the retrovirus XMRV. Science 2011, 333, 97–101. [Google Scholar] [CrossRef] [PubMed]
  23. Bamunusinghe, D.; Naghashfar, Z.; Buckler-White, A.; Plishka, R.; Baliji, S.; Liu, Q.; Kassner, J.; Oler, A.; Hartley, J.; Kozak, C.A. Sequence diversity, intersubgroup relationships, and origins of the mouse leukemia gammaretroviruses of laboratory and wild mice. J. Virol. 2016, 90, 4186–4198. [Google Scholar] [CrossRef] [PubMed]
  24. Tomonaga, K.; Coffin, J.M. Structure and distribution of endogenous nonecotropic murine leukemia viruses in wild mice. J. Virol. 1998, 72, 8289–8300. [Google Scholar] [PubMed]
  25. Cloyd, M.W.; Hartley, J.W.; Rowe, W.P. Cell-surface antigens associated with recombinant mink cell focus-inducing murine leukemia viruses. J. Exp. Med. 1979, 149, 702–712. [Google Scholar] [CrossRef] [PubMed]
  26. Keele, B.F.; Giorgi, E.E.; Salazar-Gonzalez, J.F.; Decker, J.M.; Pham, K.T.; Salazar, M.G.; Sun, C.; Grayson, T.; Wang, S.; Li, H.; et al. Identification and characterization of transmitted and early founder virus envelopes in primary HIV-1 infection. Proc. Natl. Acad. Sci. USA 2008, 105, 7552–7557. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Van Hoeven, N.S.; Miller, A.D. Use of different but overlapping determinants in a retrovirus receptor accounts for non-reciprocal interference between xenotropic and polytropic murine leukemia viruses. Retrovirology 2005, 2, 76. [Google Scholar] [CrossRef] [PubMed]
  28. 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] [PubMed]
  29. Kimura, M. A simple method for estimating evolutionary rates of base substitutions through comparative studies of nucleotide sequences. J. Mol. Evol. 1980, 16, 111–120. [Google Scholar] [CrossRef] [PubMed]
  30. Edgar, R.C. MUSCLE: Multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 2004, 32, 1792–1797. [Google Scholar] [CrossRef] [PubMed]
  31. Jern, P.; Stoye, J.P.; Coffin, J.M. Role of APOBEC3 in genetic diversity among endogenous murine leukemia viruses. PLoS Genet. 2007, 3, e183. [Google Scholar] [CrossRef] [PubMed]
  32. Lamont, C.; Culp, P.; Talbott, R.L.; Phillips, T.R.; Trauger, R.J.; Frankel, W.N.; Wilson, M.C.; Coffin, J.M.; Elder, J.H. Characterization of endogenous and recombinant proviral elements of a highly tumorigenic AKR cell line. J. Virol. 1991, 65, 4619–4628. [Google Scholar] [PubMed]
  33. Kent, W.J. BLAT—The BLAST-like alignment tool. Genome Res. 2002, 12, 656–664. [Google Scholar] [CrossRef] [PubMed]
  34. Zerbino, D.R.; Achuthan, P.; Akanni, W.; Amode, M.R.; Barrell, D.; Bhai, J.; Billis, K.; Cummins, C.; Gall, A.; Giron, C.G.; et al. Ensembl 2018. Nucleic Acids Res. 2018, 46, D754–D761. [Google Scholar] [CrossRef] [PubMed]
  35. Bamunusinghe, D.; Liu, Q.; Lu, X.; Oler, A.; Kozak, C.A. Endogenous gammaretrovirus acquisition in Mus musculus subspecies carrying functional variants of the XPR1 virus receptor. J. Virol. 2013, 87, 9845–9855. [Google Scholar] [CrossRef] [PubMed]
  36. Wang, J.R.; de Villena, F.P.; McMillan, L. Comparative analysis and visualization of multiple collinear genomes. BMC Bioinform. 2012, 13 (Suppl. 3), S13. [Google Scholar]
  37. Cloyd, M.W.; Thompson, M.M.; Hartley, J.W. Host range of mink cell focus-inducing viruses. Virology 1985, 140, 239–248. [Google Scholar] [CrossRef]
  38. Kozak, C.; Rowe, W.P. Genetic mapping of xenotropic leukemia virus-inducing loci in two mouse strains. Science 1978, 199, 1448–1449. [Google Scholar] [CrossRef] [PubMed]
  39. Kozak, C.A.; Rowe, W.P. Genetic mapping of xenotropic murine leukemia virus-inducing loci in five mouse strains. J. Exp. Med. 1980, 152, 219–228. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  40. Golemis, E.A.; Speck, N.A.; Hopkins, N. Alignment of U3 region sequences of mammalian type C viruses: Identification of highly conserved motifs and implications for enhancer design. J. Virol. 1990, 64, 534–542. [Google Scholar] [PubMed]
  41. Morrison, H.L.; Soni, B.; Lenz, J. Long terminal repeat enhancer core sequences in proviruses adjacent to c-myc in T-cell lymphomas induced by a murine retrovirus. J. Virol. 1995, 69, 446–455. [Google Scholar] [PubMed]
  42. Spies, C.P.; Ritter, G.D., Jr.; Mulligan, M.J.; Compans, R.W. Truncation of the cytoplasmic domain of the simian immunodeficiency virus envelope glycoprotein alters the conformation of the external domain. J. Virol. 1994, 68, 585–591. [Google Scholar] [PubMed]
  43. Yang, H.; Bell, T.A.; Churchill, G.A.; Pardo-Manuel de Villena, F. On the subspecific origin of the laboratory mouse. Nat. Genet. 2007, 39, 1100–1107. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Yang, H.; Wang, J.R.; Didion, J.P.; Buus, R.J.; Bell, T.A.; Welsh, C.E.; Bonhomme, F.; Yu, A.H.; Nachman, M.W.; Pialek, J.; et al. Subspecific origin and haplotype diversity in the laboratory mouse. Nat. Genet. 2011, 43, 648–655. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Miller, D.G.; Miller, A.D. Tunicamycin treatment of CHO cells abrogates multiple blocks to retrovirus infection, one of which is due to a secreted inhibitor. J. Virol. 1992, 66, 78–84. [Google Scholar] [PubMed]
  46. Yan, Y.; Liu, Q.; Kozak, C.A. Six host range variants of the xenotropic/polytropic gammaretroviruses define determinants for entry in the XPR1 cell surface receptor. Retrovirology 2009, 6, 87. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  47. Battini, J.L.; Heard, J.M.; Danos, O. Receptor choice determinants in the envelope glycoproteins of amphotropic, xenotropic, and polytropic murine leukemia viruses. J. Virol. 1992, 66, 1468–1475. [Google Scholar] [PubMed]
  48. Argaw, T.; Wilson, C.A. Detailed mapping of determinants within the porcine endogenous retrovirus envelope surface unit identifies critical residues for human cell infection within the proline-rich region. J. Virol. 2012, 86, 9096–9104. [Google Scholar] [CrossRef] [PubMed]
  49. Bahrami, S.; Duch, M.; Pedersen, F.S. Change of tropism of SL3-2 murine leukemia virus, using random mutational libraries. J. Virol. 2004, 78, 9343–9351. [Google Scholar] [CrossRef] [PubMed]
  50. Sherr, C.J.; Lieber, M.M.; Todaro, G.J. Mixed splenocyte cultures and graft versus host reactions selectively induce an S-tropic murine type C virus. Cell 1974, 1, 55–58. [Google Scholar] [CrossRef]
  51. Rosenke, K.; Lavignon, M.; Malik, F.; Kolokithas, A.; Hendrick, D.; Virtaneva, K.; Peterson, K.; Evans, L.H. Profound amplification of pathogenic murine polytropic retrovirus release from coinfected cells. J. Virol. 2012, 86, 7241–7248. [Google Scholar] [CrossRef] [PubMed]
  52. Barnett, A.L.; Davey, R.A.; Cunningham, J.M. Modular organization of the Friend murine leukemia virus envelope protein underlies the mechanism of infection. Proc. Natl. Acad. Sci. USA 2001, 98, 4113–4118. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Viruses 10 00418 g001 550
Figure 1. Comparative analysis of AKR6 and AKR40 X-MLVs. (a) Highlighter plots use vertical lines to identify nucleotide mismatches relative to the Bxv1 Xmv. Brackets indicate two substituted segments in AKR6 relative to Bxv1. (b) Sequence relationships between AKR6 and related endogenous retroviruses (ERVs) and PCR-amplified segments from AKR mice. The env-U3-R segment of the viral genome diagramed at the top notes the locations of the three variable domains (VRA, VRB, and VRC), the proline rich region (PRR) and the receptor binding domain (RBD) within the surface (SU) domain of Env, the cytoplasmic tail (CT) and R-peptide (R) in the transmembrane (TM) domain, and the enhancer framework (EF) within U3. Six black lines show segmental sequence comparisons to AKR6 for Bxv1 and 5 AKR-derived AKR6-related sequences, two from the AKR genome representing a solo long terminal repeat (LTR) and a partial env with two duplications associated with sequencing gaps (red triangles), and three PCR products from AKR amplified using the indicated primers (in red). The orange line positions the two replacements common to the P-MLVs. The numbers above each line show % identity to AKR6, and the numbers below the line identify the start and stop positions of each segment relative to AKR6. Identity > 99% is shown with a thick line.
Figure 1. Comparative analysis of AKR6 and AKR40 X-MLVs. (a) Highlighter plots use vertical lines to identify nucleotide mismatches relative to the Bxv1 Xmv. Brackets indicate two substituted segments in AKR6 relative to Bxv1. (b) Sequence relationships between AKR6 and related endogenous retroviruses (ERVs) and PCR-amplified segments from AKR mice. The env-U3-R segment of the viral genome diagramed at the top notes the locations of the three variable domains (VRA, VRB, and VRC), the proline rich region (PRR) and the receptor binding domain (RBD) within the surface (SU) domain of Env, the cytoplasmic tail (CT) and R-peptide (R) in the transmembrane (TM) domain, and the enhancer framework (EF) within U3. Six black lines show segmental sequence comparisons to AKR6 for Bxv1 and 5 AKR-derived AKR6-related sequences, two from the AKR genome representing a solo long terminal repeat (LTR) and a partial env with two duplications associated with sequencing gaps (red triangles), and three PCR products from AKR amplified using the indicated primers (in red). The orange line positions the two replacements common to the P-MLVs. The numbers above each line show % identity to AKR6, and the numbers below the line identify the start and stop positions of each segment relative to AKR6. Identity > 99% is shown with a thick line.
Viruses 10 00418 g001
Viruses 10 00418 g002 550
Figure 2. Phylogenetic trees of four domains of the MLV genome. (a) RTpol, (b,c) the env and U3-LTR enhancer substitutions in AKR6, (d) the 5’ VRA-containing end of env. The percentage of trees in which the associated taxa clustered together is shown next to the branches. The trees are drawn to scale, with the branch lengths measured in the number of substitutions per site. Red arrows identify AKR40 and AKR6 MLVs, and green arrows identify AKR6-related PCR products or mined sequences from the AKR/J genome. All of the Xmvs, Pmvs, and Mpmvs are from the B6 reference genome.
Figure 2. Phylogenetic trees of four domains of the MLV genome. (a) RTpol, (b,c) the env and U3-LTR enhancer substitutions in AKR6, (d) the 5’ VRA-containing end of env. The percentage of trees in which the associated taxa clustered together is shown next to the branches. The trees are drawn to scale, with the branch lengths measured in the number of substitutions per site. Red arrows identify AKR40 and AKR6 MLVs, and green arrows identify AKR6-related PCR products or mined sequences from the AKR/J genome. All of the Xmvs, Pmvs, and Mpmvs are from the B6 reference genome.
Viruses 10 00418 g002
Viruses 10 00418 g003 550
Figure 3. Distribution of Xmv67 and F2–R1 in common inbred mouse strains.
Figure 3. Distribution of Xmv67 and F2–R1 in common inbred mouse strains.
Viruses 10 00418 g003
Viruses 10 00418 g004 550
Figure 4. ERVs related to the AKR6 substitutions map to chromosomal regions that originated in Asian mice. (a) Horizontal tracks represent a 10-Mb segment of chromosome 11 in selected strains and subspecies typed for Xmv67. Chromosomal regions originating from M. m. domesticus are in blue and M. m. musculus in red. PCR typing is shown on the right. The position of the AKR ERV insertion site at Chr11:97 (based on NCBI37/mm9) is marked by an inverted triangle. (b) Polymorphic SNPs surrounding the Xmv67 insertion site are shown for a subset of the tested mice. Yellow highlights the AKR-specific SNPs. (c) Distribution of Xmv67 in subspecies of M. musculus.
Figure 4. ERVs related to the AKR6 substitutions map to chromosomal regions that originated in Asian mice. (a) Horizontal tracks represent a 10-Mb segment of chromosome 11 in selected strains and subspecies typed for Xmv67. Chromosomal regions originating from M. m. domesticus are in blue and M. m. musculus in red. PCR typing is shown on the right. The position of the AKR ERV insertion site at Chr11:97 (based on NCBI37/mm9) is marked by an inverted triangle. (b) Polymorphic SNPs surrounding the Xmv67 insertion site are shown for a subset of the tested mice. Yellow highlights the AKR-specific SNPs. (c) Distribution of Xmv67 in subspecies of M. musculus.
Viruses 10 00418 g004
Viruses 10 00418 g005 550
Figure 5. Multiple replacement mutations alter the C-terminal end of the AKR6 SU Env. (a) AKR6 host range differs from other X-MLVs [35,45,46]. (b) Alignment of the amino acid sequences in the replaced segment of the AKR6 Env compared to X-MLVs with shared or different host range properties. Red lines indicate the variable VRB and PRR domains. Red dots mark sites that distinguish viruses able to use a mutant XPR1 receptor and blue dots mark viruses that differ in their sensitivity to glycosylation inhibition.
Figure 5. Multiple replacement mutations alter the C-terminal end of the AKR6 SU Env. (a) AKR6 host range differs from other X-MLVs [35,45,46]. (b) Alignment of the amino acid sequences in the replaced segment of the AKR6 Env compared to X-MLVs with shared or different host range properties. Red lines indicate the variable VRB and PRR domains. Red dots mark sites that distinguish viruses able to use a mutant XPR1 receptor and blue dots mark viruses that differ in their sensitivity to glycosylation inhibition.
Viruses 10 00418 g005

Share and Cite

MDPI and ACS Style

Bamunusinghe, D.; Skorski, M.; Buckler-White, A.; Kozak, C.A. Xenotropic Mouse Gammaretroviruses Isolated from Pre-Leukemic Tissues Include a Recombinant. Viruses 2018, 10, 418. https://doi.org/10.3390/v10080418

AMA Style

Bamunusinghe D, Skorski M, Buckler-White A, Kozak CA. Xenotropic Mouse Gammaretroviruses Isolated from Pre-Leukemic Tissues Include a Recombinant. Viruses. 2018; 10(8):418. https://doi.org/10.3390/v10080418

Chicago/Turabian Style

Bamunusinghe, Devinka, Matthew Skorski, Alicia Buckler-White, and Christine A. Kozak. 2018. "Xenotropic Mouse Gammaretroviruses Isolated from Pre-Leukemic Tissues Include a Recombinant" Viruses 10, no. 8: 418. https://doi.org/10.3390/v10080418

APA Style

Bamunusinghe, D., Skorski, M., Buckler-White, A., & Kozak, C. A. (2018). Xenotropic Mouse Gammaretroviruses Isolated from Pre-Leukemic Tissues Include a Recombinant. Viruses, 10(8), 418. https://doi.org/10.3390/v10080418

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