Could Phylogenetic Analysis Be Used for Feline Leukemia Virus (FeLV) Classification?

Abstract

The surface envelope (SU) protein determines the cell tropism and consequently the pathogenesis of the feline leukemia virus (FeLV) in felids. Recombination of exogenous FeLV (exFeLV) with endogenous retroviruses (enFeLV) allows the emergence of more pathogenic variants. Currently, phenotypic testing through interference assays is the only method to distinguish among subgroups—namely, FeLV-A, -B, -C, -E, and -T. This study proposes a new method for FeLV classification based on molecular analysis of the SU gene. A total of 404 publicly available SU sequences were used to reconstruct a maximum likelihood tree. However, only 63 of these sequences had available information about phenotypic tests or subgroup assignments. Two major clusters were observed: (a) clade FeLV-A, which includes FeLV-A, FeLV-C, FeLV-E, and FeLV-T sequences, and (b) clade enFeLV, which includes FeLV-B and enFeLV strains. We found that FeLV-B, FeLV-C, FeLV-E, and FeLV-T SU sequences share similarities to FeLV-A viruses and most likely arose independently through mutation or recombination from this strain. FeLV-B and FeLV-C arose from recombination between FeLV-A and enFeLV viruses, whereas FeLV-T is a monophyletic subgroup that has probably originated from FeLV-A through combined events of deletions and insertions. Unfortunately, this study could not identify polymorphisms that are specifically linked to the FeLV-E subgroup. We propose that phylogenetic and recombination analysis together can explain the current phenotypic classification of FeLV viruses.

1. Introduction

Feline leukemia virus (FeLV) is an exogenous gammaretrovirus horizontally transmitted among domestic and wild felids and is known to cause a range of diseases, including lymphadenopathy, anemia, bone marrow suppression, immune suppression, lymphoma, leukemia, and ultimately death [1,2,3,4]. FeLV is shed in saliva, nasal secretions, urine, feces, and milk of infected cats and is mainly transmitted horizontally through saliva, blood, and other body fluids by close contact between cats. Transmission can also take place from an infected mother cat to her kittens, either before birth or while nursing [5]. FeLV has the highest case fatality rate in domestic cats out of all major feline viruses in spite of remaining one of the few retroviral diseases with an available effective vaccine [6,7]. A clear reduction in incidence was observed in countries where the vaccine is efficiently established. Nonetheless, FeLV remains one of the most common viral agents in countries where vaccination is not a common practice [8].

enFeLV is the endogenous counterpart of exogenous FeLV and represents generally nonfunctional elements, naturally present in the genome of the genus Felis with different numbers of copies [4,9,10]. These endogenous elements invaded the feline genome prior to the speciation of cats and are transmitted from parent to offspring as integral components of chromosomes [11]. Despite not forming infectious particles or inducing disease in the host, the expression of enFeLV transcripts has been observed in many tissue types of healthy or FeLV-positive cats [12,13]. In addition, enFeLV seems to act on the biology of exogenous retroviral infection either increasing or decreasing susceptibility to FeLV infection [14,15,16]. FeLV and enFeLV are approximately 86% similar at the nucleotide level and can readily recombine to generate new FeLV variants [17].

The entry of FeLV into the host cell is mediated by the binding of the envelope glycoprotein (SU or gp70) to host receptors, which vary depending on the FeLV subgroup. Currently, FeLV classification is primarily based on phenotypic tests to evaluate receptor usage, such as viral interference assays (IAs). Six different groups have been described so far: FeLV-A, FeLV-B, FeLV-C, FeLV-D, FeLV-E, and FeLV-T [17,18]. Variations in SU protein affect the affinity to the cell receptor and can ultimately influence the clinical manifestation [19]. FeLV-A is the most commonly described subgroup of FeLV and has been reported to be less pathogenic than other FeLV subgroups [6]. FeLV-B arises through recombination events between the envelope (env) gene of FeLV-A and enFeLV transcripts [20]. This process occurs during retroviral transcriptase-directed DNA synthesis following copackaging of endogenous and exogenous viral genomes. Furthermore, this might allow for a change in cellular tropism and a shift in disease progression [21]. The emergence of FeLV-B following FeLV-A infection is considered to result in both higher morbidity and mortality. While FeLV-A is commonly transmitted horizontally, similar processes have only rarely been observed with FeLV-B [4,22,23]. FeLV-C is a less common subgroup that arises from de novo mutations in the env gene of FeLV-A and has been associated with the development of aplastic anemia [6,24]. FeLV-D is the result of recombination between exogenous viruses and domestic cat endogenous retrovirus (ERV-DC) that are divergent from enFeLV [25]. FeLV-E is a recently described subgroup that was isolated from natural thymic lymphoma in cats [26]. Ultimately, FeLV-T is a T-cytopathic virus, isolated from a cat with FeLV feline acquired immune deficiency syndrome (FeLV-FAIDS) that has emerged via the mutation of FeLV-A strains [26]. Although these recombination events and mutation in env typically have led to increased pathogenicity, the FeLV variants are associated with decreased ability for horizontal transmission requiring FeLV-A as a helper virus [1,20,25].

Viral interference assays have determined the virus biological activity and coinfection profiles of FeLV viruses [6]. Even with the advancement of sequencing technologies, these classical techniques continue to be the main method to identify historic and novel subgroups [17]. Data linking the current phenotypic classification of FeLV samples n to their phylogenetic history are lacking. Only a few studies have aimed to characterize FeLV subgroups using a phylogenetic approach [3,6,17]. The objective of this study was to investigate the phylogenetic relationship underlying the FeLV subgroups and explore more deeply the recombination process driving FeLV evolution.

2. Materials and Methods

2.1. Dataset

Partial and complete FeLV SU sequences were selected from the GenBank (http://www.ncbi.nlm.nih.gov/nucleotide/ accessed date on 21 March 2021). To ensure the selection of high-quality data, sequences were chosen only when meeting the following criteria: (a) one sequence per host, (b) minimum length of 1050 base pairs (bp), and (c) absence of premature stop codons. Nucleotide alignments were generated using the MUSCLE software [27] and visualized in AliView [28].

2.2. Recombination Analysis

Recombinant sequences were detected using seven different methods (RDP, GENECONV, Chimaera, MaxChi, Bootscan, SiScan, and 3Seq) in RDP4 v4.85 [29]. Default parameters were used and recombination events were only accepted if three or more of these methods detected breakpoints with significant p-values (p < 0.05). Sequences were also imported into SplitsTree4 [30] and tested for recombination using the PHI test [31]. PHI test calculates the pairwise homoplasy index (PHI) as the mean of the refined incompatibility scores obtained for nearby nucleotide sites across the sequences.

Recombinant sequences were submitted to bootscanning analysis in Simplot 3.5.1 [32] to detect putative recombination among subgroups. The Kimura 2-parameter model was used as a distance model on a sliding window of 200 nucleotides (nt) by increments of 20 nt. Recombinant samples identified by RDP4 were used as query sequences in Simplot. For the bootscan, non-recombinant sequences were used as references.

2.3. Phylogenetic Reconstruction

A maximum likelihood (ML) tree incorporating the best-fitted nucleotide substitution model (GTR + F + R5) was reconstructed in the IQtree web server [33]. The robustness of the resulting tree was evaluated by rapid bootstrapping (UFB) with 1000 pseudoreplicate datasets. The phylogenetic tree was visualized in FigTree v1.4.4 [34] and drawn using R [35].

3. Results

3.1. Dataset Curation

Four hundred and fifteen FeLV SU sequences met our criteria for data quality and were downloaded from the GenBank (1368 base pairs, Supplementary Table S1). However, only 63 of these sequences had information available on phenotypic tests or subgroup assignments (39 FeLV-A, 5 FeLV-B, 1 FeLV-C, 7 FeLV-D, 1 FeLV-E, 2 FeLV-T, and 4 enFeLV). When aligning FeLV-D samples with other FeLV/enFeLV sequences, vast differences between sequences were observed. Therefore, the phylogenetic analysis of this subgroup could not be rendered possible, resulting ultimately in the exclusion of the seven FeLV-D sequences from the following analyses. Even though the exclusion does not preclude phylogenetic analyses, the reduced number of sequences may interfere with the extent of conclusions regarding FeLV evolution.

3.2. Different Patterns of Recombination Were Identified

The results for the PHI test support that FeLV evolution might be modeled by recombination (p > 0.01). To detect putative recombinant sequences and minimize the potentially disruptive impact of recombination on phylogenetic interpretation [29,36], the dataset was analyzed in RDP4. Overall, 55 FeLV SU sequences displayed at least 1 of 24 different recombination patterns and were subsequently treated as recombinant sequences. These sequences comprise one FeLV-A, four FeLV-B, and a single FeLV-C sequence (Figure 1 and Figure 2, and Supplementary Table S1). After excluding those 55 recombinant sequences from the alignment, no significant recombination signal was detected using the PHI test (p = 0.8244), therefore suggesting that all recombinant sequences were effectively identified. Simplot analysis allowed for the inference of the most likely origin of each recombinant fragment (Figure 2).

3.3. FeLV SU Gene Can Be Divided into Two Large Clades

Phylogenetic reconstruction was performed using all FeLV samples included in this study (Figure 3, Supplementary Table S1). The incorporation of recombinant sequences in this phylogenetic analysis aims to identify the impact of recombination on FeLV subgroups. In addition, it helps clarify the evolutionary relationship between these samples. The phylogenetic tree displays the subdivision of FeLV SU sequences into two well-supported (UFB > 90) main clades (Figure 3): (a) clade FeLV-A, which includes all FeLV-A sequences intermingled with FeLV-C, FeLV-E, and FeLV-T; (b) clade enFeLV, which encompasses all FeLV-B and enFeLV sequences included in this study. The split of the dataset into two different clades is corroborated by the analysis of the SU gene structure (Figure 2 and Figure 3).

3.4. The SU Envelope Gene Has Five Variable Regions; the Pattern of the Variable Regions Is Conserved Intraclade

While most of the sequences inside clade FeLV-A present insertions in VRA and VRC regions, sequences in clade enFeLV have a 30 bp and 15 bp deletion in this region, respectively.

A similar pattern was observed in the VRB region: The 69 bp deletion is exclusively shared by viruses of the FeLV-A clade but not by viruses in the enFeLV clade. Additional variable regions (here labeled as VRD and VRE) were identified inside the proline-rich region (PRR) and in the C-domain region. VRD was exclusively observed in sequences grouping within the enFeLV clade, while VRE insertion was exhibited by 13 sequences within clade FeLV-A (Figure 3 and Supplementary Table S1).

3.5. The FeLV-A Clade

In total, 315 sequences clustered in clade FeLV-A, including 39 FeLV-A sequences (39/39 publicly available), 2 FeLV-T (2/2 publicly available), 1 FeLV-C (1/1 publicly available), 1 FeLV-E (1/1 publicly available), and 272 sequences without information about subgroup assignment. Seven recombinant sequences with unique patterns of recombination could be identified within this clade (Figure 1 and Supplementary Table S1): one FeLV-A (AB635503, recombination pattern 20), one FeLV-C (M14331, recombination pattern 18), and other five sequences without subgroup assignment (recombination patterns 8, 8 + 14, 10, 18, and 24). Four of these samples presented patterns of recombination that seem to exclusively include fragments derived from the FeLV-A clade, suggesting that sequences joining clade FeLV-A can recombine with themselves. Conversely, samples with recombination patterns 8 (AB635618), 8 + 14 (AB635637), and 10 (AY662460) have the C-domain derived from endogenous viruses (clade enFeLV), despite presenting most of their SU derived from clade FeLV-A.

3.6. The enFeLV Clade

Clade enFeLV contains 89 FeLV SU sequences, including 4 enFeLV (4/4 publicly available), 5 FeLV-B (5/5 publicly available), and 80 sequences without information on subgroup assignment (Figure 3). Among these samples, 4 FeLV-B and other 44 samples without subgroup assignment presented evidence for recombination, altogether displaying 19 different patterns of recombination. Despite being recombinants, the majority of sequences (96%) preserved the RBD region derived from endogenous strains (absence of VRA and VRC but presence of VRB). In addition, most of these samples (78%) exhibited an enFeLV profile within the PRR region (Supplementary Table S1). Together, these results suggest that recombination appears to be an important factor for the evolution of FeLV-B and enFeLV samples. This is supported by the fact that 54% of the samples within this clade are recombinant sequences. Moreover, specific signatures inside RBD and PRR regions might be critical for these viruses, with almost 80% of them presenting conserved motifs at the gene level.

3.7. Geographical Distribution of FeLV-A Clade

A second tree was reconstructed using only non-recombinant FeLV sequences (Table S1). In addition to the enFeLV basal clade, this phylogeny suggests that clade FeLV-A is divided into two well-defined supported groups. Curiously, these clades group sequences from different geographical origins: While the superior clade includes only sequences from Japan, the second clade includes sequences from Brazil, the United States, and the United Kingdom. On the other hand, the non-recombinant clade enFeLV mostly grouped sequences from Japan and China (Figure 4, Supplementary Table S1). Unfortunately, FeLV sequences are not systematically collected on a global scale, therefore preventing further analysis on the relationship between phylogenetic classification and geographical distribution.

4. Discussion

In the last 40 years, viral interference assays have been the only method used to distinguish and define FeLV subgroups [17,36]. Analysis of the superinfection interference properties of FeLV variants has led to the identification of four interference types among FeLV subgroups. Eventually, these subtypes were associated with specific clinical phenotypes: FeLV-A uses the feline thiamine transport protein 1 (feThTr1) [37] and has been linked to macrocytic anemia, immunosuppression, and lymphoma [2]; FeLV-B utilizes the inorganic phosphate-sodium symporter fePit1 (SLC20A1), and the closely related fePit2 (SLC20A2) [38,39,40] and is tumorigenic [2]; FeLV-C uses a heme export protein called FLVCR [41,42] and has been associated with the development of aplastic anemia [43]; finally, FeLV-T requires Pit1 and a soluble cofactor expressed by endogenous FeLV-related sequences (FeLIX) [18,44] and has been connected to fatal immunodeficiency [45]. FeLV-E was only recently described and its receptor remains unidentified [26]. In spite of being the only described method to date, phenotypic tests are laborious and time consuming. This, in turn, creates a technical barrier and restricts extensive FeLV classification, ultimately interfering with molecular epidemiology or evolutionary studies. With the advancement of next-generation sequencing, genetic tests have become cheaper, faster, more reliable, and, finally, also more precise than phenotypic characterization. However, despite the clear advantages of sequencing technologies, only a few studies thus far have tried to correlate phenotypic classification to phylogenetic reconstruction for FeLV viruses [3,17,46,47]. Here, we used molecular data to characterize FeLV subgroups and outline a proposition for a new classification system based on phylogenetic analysis.

Despite being classified into different subgroups, FeLV-A, FeLV-C, FeLV-E, and FeLV-T SU sequences grouped within the same clade in the phylogenetic tree (clade FeLV-A, Figure 3). All sequences within this clade presented a similar gene organization, sharing an insertion of 30 bp in VRA and 15 bp in VRC, in addition to a 69 bp deletion in VRB and a 27 bp deletion in VRD (Figure 3). FeLV-A sequences are the most abundant strains within this cluster and were found scattered through the whole clade. Therefore, it is reasonable to assume that the most recent common ancestor of clade FeLV-A presented a gene organization related to FeLV-A and used feThTr1 as its host receptor. The inclusion of other FeLV subgroups with different host receptors, but similar gene organization within this clade may suggest that FeLV-C, FeLV-E, and FeLV-T arose in separate monophyletic events from FeLV-A viruses. This process most likely includes mutation and/or recombination. In addition, it seems that the evolution of sequences inside clade FeLV-A is somehow related to its respective geographical origin. However, insufficient data regarding sampling locations of the majority of the sequences analyzed here prevent further conclusions.

Apparently, FeLV-C has arisen by means of recombination (Figure 1). Our analysis suggests that FeLV-C is the result of a recombination event that included mostly the FeLV-A SU gene and only a small fragment of endogenous strains (Figure 1). The recombination breakpoint in the C domain might be the explanation for the fact that FeLV-C has a different receptor (FLVCR) than the one used by most samples grouped within the clade FeLV-A (feThTr1) [37,41]. A previous study investigating the role of the C-domain region of FeLV-C strain proposes that this region is critical for efficient SU binding. The region may function as a second receptor-binding domain, which, in addition to RBD, interacts with the host receptor to initiate virus infection [48]. The phenotypic results for FeLV-C seem to corroborate the molecular analyses performed in this study, suggesting that recombination may be the main evolutionary event defining the receptor change in this variant.

FeLV-T, in contrast, may have arisen from the accumulation of mutations in FeLV-A virions, allowing the change of host receptor from FETHTR1 to Pit1. Both FeLV-T samples included in this study were grouped in a monophyletic clade and, unlike all sequences inside clade FeLV-A, exhibited an 18 bp deletion in VRA and an 18 bp insertion in the VRE region. Even though the deletion of 18 bp in VRA is smaller, it still lies in the same position as the 30 bp deletion presented by FeLV-B viruses, which could explain the use of the same receptor (Pit1) for these two strains. Intriguingly, Shojima et al. (2006) in an ex vivo model using different cell lines proposed that the host ranges of FeLV-T and FeLV-B were not identical. Furthermore, they suggested a different Pit1 usage at the post-binding level for the two strains [49]. Another 11 sequences with no information on subgroups also presented the 18 bp deletion in VRE. Only one of these sequences (M87886) has the 18 bp insertion in the VRA and is grouped in the same monophyletic clade as FeLV-T. This unclassified sequence most likely represents another strain of FeLV-T. Together, these results support the idea that polymorphisms in VRA are notoriously more important to defining the host receptor for strains grouped within clade FeLV-A than polymorphisms in VRE.

Unfortunately, this study was not able to find genetic polymorphisms or evolutionary mechanisms in the FeLV-E strain that could explain its receptor usage. We used entropy analysis to compare FeLV-E with all FeLV-A sequences, concluding that no amino acid site could be specifically associated with the FeLV-A strain (p > 0.05). In addition, no recombination breakpoint was detected in this sequence. The fact that only one sequence of this subgroup has been sampled so far may have prevented any conclusions about the relationship between phenotypic testing and phylogenetic clustering.

Clade enFeLV encompassed FeLV-B, enFeLV, and most of the recombinant sequences identified in this study (87%). The similarity of the gene structure seems to be the main factor driving the clustering of these sequences. In general, these samples presented a 30 bp deletion in the VRA (100% of the sequences in the clade), a 15 bp deletion in VRC (95%), and 2 insertions of 13 and 14 bp in the VRD (88%). Despite sharing a similar structure, enFeLV and FeLV-B sequences were found segregated within the enFeLV clade in our phylogenetic analysis. Notably, FeLV-B and all the recombinant sequences without subgroup assignment were basal to a monophyletic clade that includes, among other sequences, all four enFeLV samples.

Corroborating previous studies, our analysis shows that FeLV-B arose several times through recombination between FeLV-A and endogenous strains most likely after copackaging of expressed transcripts into a single virion [50,51]. These results suggest that FeLV-B strains must include the RBD and PRR regions derived from endogenous viruses alongside a short and variable length fragment derived from exogenous strains at the 3′ end of the SU gene (Figure 1 and Figure 2). Pandey et al. (1991), examining exogenous/endogenous recombination during in vitro infections, have revealed that replication efficiency and cellular tropism depend on the length and region of the enFeLV sequence incorporated into the FeLV-B recombinant [52]. In addition, this study shows that amino acid changes inside VRA and VRB are the main factors for the ability of FeLV-B to bind receptor Pit1 and/or Pit2. Sample JF957363 was identified as the only non-recombinant FeLV-B sample in our analysis. However, when subjecting this sample to more in-depth analysis, we found that this isolate carries the enFeLV SU, whereas long terminal repeats (LTRs) and Gag were derived from exogenous viruses. This accounts for the relationship with the enFeLV sequences in the phylogenetic tree [53].

Frequent recombination generated FeLV strains combining endogenous and exogenous viruses. Only seven recombinant sequences were grouped within clade FeLV-A of which all maintained the RBD and PRR regions from FeLV-A strains. The C domain is, however, variable and can apparently be either derived from exogenous or endogenous viruses (Figure 1 and Supplementary Table S1). Once the RBD is completely exogenous, these recombinants probably display the same interference pattern as the exogenous FeLV-A viruses. On the other hand, clade enFeLV clustered a total of 48 recombinant sequences. Out of these sequences, the majority present the RBD and the PRR regions derived from endogenous strains. Following the pattern observed for recombinant sequences within clade FeLV-A, the C-domain region in these sequences was found to be derived from exogenous or endogenous viruses. It is interesting to highlight that the farther the sequence is from the root of the enFeLV clade in the phylogenetic tree, the greater is the size of the fragment derived from exogenous viruses in these recombinant sequences (Figure 2).

Here, we used both genetic and evolutionary approaches to explain the current phenotypic classification of FeLV viruses. We found that FeLV-B, FeLV-C, FeLV-E, and FeLV-T SU gene sequences share similarities to FeLV-A viruses. They most likely arose independently through mutation or recombination from the FeLV-A strain. FeLV-B is a classical product of recombination between FeLV-A and enFeLV strains that presents the RBD and PRR regions derived from endogenous viruses. Despite also arising from recombination events, FeLV-C incorporated both RBD and PRR regions from FeLV-A strains. FeLV-T is a monophyletic subgroup that probably originated from FeLV-A through combined events of deletions in the RBD region and insertions in the C domain of the SU gene. Unfortunately, this study could not identify polymorphisms that are specifically linked to the FeLV-E subgroup and could explain its phenotypic results for receptor usage. To conclude, we propose that the use of phylogenetic analysis can be efficiently used to classify FeLV subgroups.

Author Contributions

Conceptualization, methodology, investigation, validation, formal analysis, investigation, data curation, writing—original draft preparation, visualization, L.C.-O.; methodology, formal analysis, investigation, C.T.; writing—review and editing, funding acquisition, P.M.R.; investigation, writing—review and editing, supervision, funding acquisition, A.C.F.; conceptualization, methodology, investigation, validation, formal analysis, investigation, data curation, writing—original draft preparation, writing—review and editing, visualization, supervision, D.M.J. All authors have read and agreed to the published version of the manuscript.

Funding

This study was financed in part by the Coordenação de Aperfeiçoamento de Pessoal de Nível Superior, Brasil (CAPES), and by Conselho Nacional de Desenvolvimento Científico e Tecnológico, CNPq project 429113/2018-8. L.C.O. is a CAPES PhD fellow. A.C.F. is a PQ2 research fellow of CNPq.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Figures

Figure 1. Scheme of the recombinant patterns obtained in RPD4. The number of sequences representing each recombinant pattern is indicated by colored circles. Sequences belonging to clade enFeLV are indicated in yellow and sequences belonging to clade FeLV-A in green. The recombination segments are shown for each pattern. Colored circles with letters A, B or C represent FeLV-A, FeLV-B or FeLV-C sequences identified by phenotypic characterization. RBD: receptor binding site, PRR: proline-rich region, VRA, VRC, VRB, VRD, VRE: variable region A, C, B, D and E respectively.

Figure 2. Bootscan of some recombinant sequences. Recombinant sequences were submitted to bootscanning analysis in Simplot 3.5.1 [32]. The Kimura 2-parameter model was used as a distance model on a sliding window of 200 nucleotides (nt) by increments of 20 nt. For the bootscan analysis, recombinant samples identified by RDP4 were used as query sequences in Simplot and, non-recombinant sequences were used as references. A phylogeny was reconstructed by using the best-fitted nucleotide substitution model (GTR + F + R5) in the IQtree web server with some recombinant sequences (AB635493, AB635503, AB635637, AB635638, AB635673, AB635698, AB635749, M14331, MF681671, and V01172) and non-recombinant sequences according to the RDP4 analyses. Branches of some non-recombinant sequences are colored according to the clades found in the Figure 4.

Figure 3. Maximum likelihood tree of 404 FeLV SU env sequences and alignment. Branches colored in red represent recombinant sequences detected by RDP4 software. Colored circles at the tip of the branches represent sequences with a previously defined subgroup based on interference assays. Red circles denote sequences without information on subgroup assignment: yellow: enFeLV; green: FeLV-A; blue: FeLV-A; aquamarine: FeLV-B; blue: FeLV-C; purple FeLV-E; pink: FeLV-T. RBD: receptor binding site, PRR: proline-rich region, VRA, VRC, VRB, VRD, VRE: variable region A, C, B, D and E respectively.

Figure 4. Phylogenetic reconstruction of the non-recombinant FeLV SU sequences and geographical distribution of the clade FeLV-A. Maximum likelihood tree reconstructed by incorporating the best-fitted nucleotide substitution model (GTR+F+R5) in IQtree web server. The robustness of the resulting tree was evaluated by rapid bootstrapping (UFB) with 1000 pseudoreplicate datasets. The phylogenetic tree was visualized in FigTree v1.4.4 [34] and drawn using R [35]. Three main clades are presented. The yellow clade corresponds to non-recombinant sequences from clade enFeLV and the blue and purple to the non-recombinant sequences from clade A (Figure 1).

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/v14020249/s1. Table S1: Description of the FeLV SU sequences used in this study.

References

1. Levy, L.S. Advances in understanding molecular determinants in FeLV pathology. Vet. Immunol. Immunopathol.; 2008; 123, pp. 14-22. [DOI: https://dx.doi.org/10.1016/j.vetimm.2008.01.008] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/18289704]

2. Hartmann, K. Clinical aspects of feline retroviruses: A review. Viruses; 2012; 4, pp. 2684-2710. [DOI: https://dx.doi.org/10.3390/v4112684] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/23202500]

3. Watanabe, S.; Kawamura, M.; Odahara, Y.; Anai, Y.; Ochi, H.; Nakagawa, S.; Endo, Y.; Tsujimoto, H.; Nishigaki, K. Phylogenetic and Structural Diversity in the Feline Leukemia Virus Env Gene. PLoS ONE; 2013; 8, e61009. [DOI: https://dx.doi.org/10.1371/journal.pone.0061009]

4. Powers, J.A.; Chiu, E.S.; Kraberger, S.J.; Roelke-Parker, M.; Lowery, I.; Erbeck, K.; Troyer, R.; Carver, S.; VandeWoude, S. Feline Leukemia Virus (FeLV) Disease Outcomes in a Domestic Cat Breeding Colony: Relationship to Endogenous FeLV and Other Chronic Viral Infections. J. Virol.; 2018; 92, e00649-18. [DOI: https://dx.doi.org/10.1128/JVI.00649-18]

5. Vobis, M.; D’Haese, J.; Mehlhorn, H.; Mencke, N. Evidence of horizontal transmission of feline leukemia virus by the cat flea (Ctenocephalides felis). Parasitol. Res.; 2003; 91, pp. 467-470. [DOI: https://dx.doi.org/10.1007/s00436-003-0949-8]

6. Willett, B.J.; Hosie, M.J. Feline leukaemia virus: Half a century since its discovery. Vet. J.; 2013; 195, pp. 16-23. [DOI: https://dx.doi.org/10.1016/j.tvjl.2012.07.004]

7. Addie, D.D.; Toth, S.; Reid, S.; Jarrett, O.; Dennis, J.M.; Callanan, J.J. Long-term impact on a closed household of pet cats of natural infection with feline coronavirus, feline leukaemia virus and feline immunodeficiency virus. Vet. Rec.; 2000; 146, pp. 419-424. [DOI: https://dx.doi.org/10.1136/vr.146.15.419]

8. Patel, M.; Carritt, K.; Lane, J.; Jayappa, H.; Stahl, M.; Bourgeois, M. Comparative efficacy of feline leukemia virus (FeLV) inactivated whole-virus vaccine and canarypox virus-vectored vaccine during virulent FeLV challenge and immunosuppression. Clin. Vaccine Immunol.; 2015; 22, pp. 798-805. [DOI: https://dx.doi.org/10.1128/CVI.00034-15]

9. Coffin, J.M. Evolution of Retroviruses: Fossils in Our Dna. Proc. Am. Philos. Soc.; 2004; 148, pp. 264-280.

10. Polani, S.; Roca, A.L.; Rosensteel, B.B.; Kolokotronis, S.O.; Bar-Gal, G.K. Evolutionary dynamics of endogenous feline leukemia virus proliferation among species of the domestic cat lineage. Virology; 2010; 405, pp. 397-407. [DOI: https://dx.doi.org/10.1016/j.virol.2010.06.010]

11. Koshy, R.; Gallo, R.C.; Wong-Staal, F. Characterization of the endogenous feline leukemia virus-related DNA sequences in cats and attempts to identify exogenous viral sequences in tissues of virus-negative leukemic animals. Virology; 1980; 103, pp. 434-445. [DOI: https://dx.doi.org/10.1016/0042-6822(80)90202-0]

12. Busch, M.P.; Devi, B.G.; Soe, L.H.; Roy-Burman, P.; Perbal, B.; Marcel, A. Baluda Characterization of the expression of cellular retrovirus genes and oncogenes in feline cells. Hematol. Oncol.; 1983; 1, pp. 61-75. [DOI: https://dx.doi.org/10.1002/hon.2900010108] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/6329934]

13. Tandon, R.; Cattori, V.; Willi, B.; Lutz, H.; Hofmann-Lehmann, R. Quantification of endogenous and exogenous feline leukemia virus sequences by real-time PCR assays. Vet. Immunol. Immunopathol.; 2008; 123, pp. 129-133. [DOI: https://dx.doi.org/10.1016/j.vetimm.2008.01.027] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/18295344]

14. Rasmussen, H.B. Interactions between exogenous and endogenous retroviruses. J. Biomed. Sci.; 1997; 4, pp. 1-8. [DOI: https://dx.doi.org/10.1007/BF02255587] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/11725127]

15. Tandon, R.; Cattori, V.; Willi, B.; Meli, M.L.; Gomes-Keller, M.A.; Lutz, H.; Hofmann-Lehmann, R. Copy number polymorphism of endogenous feline leukemia virus-like sequences. Mol. Cell. Probes; 2007; 21, pp. 257-266. [DOI: https://dx.doi.org/10.1016/j.mcp.2007.01.003] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/17329079]

16. Tandon, R.; Cattori, V.; Pepin, A.C.; Riond, B.; Meli, M.L.; McDonald, M.; Doherr, M.G.; Lutz, H.; Hofmann-Lehmann, R. Association between endogenous feline leukemia virus loads and exogenous feline leukemia virus infection in domestic cats. Virus Res.; 2008; 135, pp. 136-143. [DOI: https://dx.doi.org/10.1016/j.virusres.2008.02.016] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/18407367]

17. Chiu, E.S.; Hoover, E.A.; VandeWoude, S. A Retrospective Examination of Feline Leukemia Subgroup Characterization: Viral Interference Assays to Deep Sequencing. Viruses; 2018; 10, 29. [DOI: https://dx.doi.org/10.3390/v10010029] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/29320424]

18. Lauring, A.S.; Anderson, M.M.; Overbaugh, J. Specificity in receptor usage by T-cell-tropic feline leukemia viruses: Implications for the in vivo tropism of immunodeficiency-inducing variants. J. Virol.; 2001; 75, pp. 8888-8898. [DOI: https://dx.doi.org/10.1128/JVI.75.19.8888-8898.2001]

19. Ramsey, I.K.; Spibey, N.; Jarrett, O. The receptor binding site of feline leukemia virus surface glycoprotein is distinct from the site involved in virus neutralization. J. Virol.; 1998; 72, pp. 3268-3277. [DOI: https://dx.doi.org/10.1128/JVI.72.4.3268-3277.1998]

20. Erbeck, K.; Gagne, R.B.; Kraberger, S.; Chiu, E.S.; Roelke-Parker, M.; VandeWoude, S. Feline Leukemia Virus (FeLV) Endogenous and Exogenous Recombination Events Result in Multiple FeLV-B Subtypes during Natural Infection. J. Virol.; 2021; 95, pp. 1-10. [DOI: https://dx.doi.org/10.1128/JVI.00353-21]

21. Odelberg, S.J.; Weiss, R.B.; Hata, A.; White, R. Template-switching during DNA synthesis by Thermus aquaticus DNA polymerase I. Nucleic Acids Res.; 1995; 23, pp. 2049-2057. [DOI: https://dx.doi.org/10.1093/nar/23.11.2049] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/7596836]

22. Jarrett, O.; Hardy, W.D.; Golder, M.C.; Hay, D. The frequency of occurrence of feline leukaemia virus subgroups in cats. Int. J. Cancer; 1978; 21, pp. 334-337. [DOI: https://dx.doi.org/10.1002/ijc.2910210314] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/204584]

23. Sarma, P.S.; Log, T. Subgroup classification of feline leukemia and sarcoma viruses by viral interference and neutralization tests. Virology; 1973; 54, pp. 160-169. [DOI: https://dx.doi.org/10.1016/0042-6822(73)90125-6]

24. Jarrett, O.; Golder, M.C.; Toth, S.; Onions, D.E.; Stewart, M.F. Interaction between feline leukaemia virus subgroups in the pathogenesis of erythroid hypoplasia. Int. J. Cancer; 1984; 34, pp. 283-288. [DOI: https://dx.doi.org/10.1002/ijc.2910340222]

25. Anai, Y.; Ochi, H.; Watanabe, S.; Nakagawa, S.; Kawamura, M.; Gojobori, T.; Nishigaki, K. Infectious endogenous retroviruses in cats and emergence of recombinant viruses. J. Virol.; 2012; 86, pp. 8634-8644. [DOI: https://dx.doi.org/10.1128/JVI.00280-12]

26. Miyake, A.; Watanabe, S.; Hiratsuka, T.; Ito, J.; Ngo, M.H.; Makundi, I.; Kawasaki, J.; Endo, Y.; Tsujimoto, H.; Nishigakia, K. Novel Feline Leukemia Virus Interference Group Based on the env Gene. J. Virol.; 2016; 90, pp. 4832-4837. [DOI: https://dx.doi.org/10.1128/JVI.03229-15]

27. Edgar, R.C. MUSCLE: Multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res.; 2004; 32, pp. 1792-1797. [DOI: https://dx.doi.org/10.1093/nar/gkh340]

28. Larsson, A. AliView: A fast and lightweight alignment viewer and editor for large datasets. Bioinformatics; 2014; 30, pp. 3276-3278. [DOI: https://dx.doi.org/10.1093/bioinformatics/btu531]

29. Martin, D.P.; Murrell, B.; Golden, M.; Khoosal, A.; Muhire, B. RDP4: Detection and analysis of recombination patterns in virus genomes. Virus Evol.; 2015; 1, pp. 1-5. [DOI: https://dx.doi.org/10.1093/ve/vev003]

30. Huson, D.H.; Bryant, D. Application of phylogenetic networks in evolutionary studies. Mol. Biol. Evol.; 2006; 23, pp. 254-267. [DOI: https://dx.doi.org/10.1093/molbev/msj030]

31. Bruen, T.C.; Philippe, H.; Bryant, D. A simple and robust statistical test for detecting the presence of recombination. Genetics; 2006; 172, pp. 2665-2681. [DOI: https://dx.doi.org/10.1534/genetics.105.048975] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/16489234]

32. Lole, K.S.; Bollinger, R.C.; Paranjape, R.S.; Gadkari, D.; Kulkarni, S.S.; Novak, N.G.; Ingersoll, R.; Sheppard, H.W.; Ray, S.C. Full-length human immunodeficiency virus type 1 genomes from subtype C-infected seroconverters in India, with evidence of intersubtype recombination. J. Virol.; 1999; 73, pp. 152-160. [DOI: https://dx.doi.org/10.1128/JVI.73.1.152-160.1999] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/9847317]

33. Trifinopoulos, J.; Nguyen, L.T.; von Haeseler, A.; Minh, B.Q. W-IQ-TREE: A fast online phylogenetic tool for maximum likelihood analysis. Nucleic Acids Res.; 2016; 44, pp. W232-W235. [DOI: https://dx.doi.org/10.1093/nar/gkw256]

34. Rambaut, A. FigTree v1.4: Tree Figure Drawing Tool 2009.

35. R Core Team. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2017.

36. Sarma, P.S.; Log, T. Viral interference in feline leukemia-sarcoma complex. Virology; 1971; 44, pp. 352-358. [DOI: https://dx.doi.org/10.1016/0042-6822(71)90266-2]

37. Mendoza, R.; Anderson, M.M.; Overbaugh, J. A Putative Thiamine Transport Protein Is a Receptor for Feline Leukemia Virus Subgroup A. J. Virol.; 2006; 80, pp. 3378-3385. [DOI: https://dx.doi.org/10.1128/JVI.80.7.3378-3385.2006] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/16537605]

38. Rudra-Ganguly, N.; Ghosh, A.K.; Roy-Burman, P. Retrovirus receptor PiT-1 of the Felis catus. Biochim. Biophys. Acta-Gene Struct. Expr.; 1998; 1443, pp. 407-413. [DOI: https://dx.doi.org/10.1016/S0167-4781(98)00241-3]

39. Takeuchi, Y.; Vile, R.G.; Simpson, G.; O’Hara, B.; Collins, M.K.; Weiss, R.A. Feline leukemia virus subgroup B uses the same cell surface receptor as gibbon ape leukemia virus. J. Virol.; 1992; 66, pp. 1219-1222. [DOI: https://dx.doi.org/10.1128/jvi.66.2.1219-1222.1992]

40. Anderson, M.M.; Lauring, A.S.; Robertson, S.; Dirks, C.; Overbaugh, J. Feline Pit2 Functions as a Receptor for Subgroup B Feline Leukemia Viruses. J. Virol.; 2001; 75, pp. 10563-10572. [DOI: https://dx.doi.org/10.1128/JVI.75.22.10563-10572.2001]

41. Tailor, C.S.; Willett, B.J.; Kabat, D. A Putative Cell Surface Receptor for Anemia-Inducing Feline Leukemia Virus Subgroup C Is a Member of a Transporter Superfamily. J. Virol.; 1999; 73, pp. 6500-6505. [DOI: https://dx.doi.org/10.1128/JVI.73.8.6500-6505.1999]

42. Quigley, J.G.; Yang, Z.; Worthington, M.T.; Phillips, J.D.; Sabo, K.M.; Sabath, D.E.; Berg, C.L.; Sassa, S.; Wood, B.L.; Abkowitz, J.L. Identification of a human heme exporter that is essential for erythropoiesis. Cell; 2004; 118, pp. 757-766. [DOI: https://dx.doi.org/10.1016/j.cell.2004.08.014]

43. Riedel, N.; Hoover, E.A.; Gasper, P.W.; Nicolson, M.O.; Mullins, J.I. Molecular analysis and pathogenesis of the feline aplastic anemia retrovirus, feline leukemia virus C-Sarma. J. Virol.; 1986; 60, pp. 242-250. [DOI: https://dx.doi.org/10.1128/jvi.60.1.242-250.1986] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/3018287]

44. Anderson, M.M.; Lauring, A.S.; Burns, C.C.; Overbaugh, J. Identification of a Cellular Cofactor Required for Infection by Feline Leukemia Virus. Science; 2000; 287, pp. 1828-1830. [DOI: https://dx.doi.org/10.1126/science.287.5459.1828] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/10710311]

45. Hoover, E.A.; Mullins, J.I.; Quackenbush, S.L.; Gasper, P.W. Experimental transmission and pathogenesis of immunodeficiency syndrome in cats. Blood; 1987; 70, pp. 1880-1892. [DOI: https://dx.doi.org/10.1182/blood.V70.6.1880.1880] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/2823940]

46. Brown, M.A.; Cunningham, M.W.; Roca, A.L.; Troyer, J.L.; Johnson, W.E.; O’Brien, S.J. Genetic characterization of feline leukemia virus from Florida panthers. Emerg. Infect. Dis.; 2008; 14, 252. [DOI: https://dx.doi.org/10.3201/eid1402.070981] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/18258118]

47. Ortega, C.; Valencia, A.C.; Duque-Valencia, J.; Ruiz-Saenz, J. Prevalence and Genomic Diversity of Feline Leukemia Virus in Privately Owned and Shelter Cats in Aburrá Valley, Colombia. Viruses; 2020; 12, 464. [DOI: https://dx.doi.org/10.3390/v12040464]

48. Rey, M.A.; Prasad, R.; Tailor, C.S. The C domain in the surface envelope glycoprotein of subgroup C feline leukemia virus is a second receptor-binding domain. Virology; 2008; 370, pp. 273-284. [DOI: https://dx.doi.org/10.1016/j.virol.2007.09.011]

49. Shojima, T.; Nakata, R.; Miyazawa, T. Host cell range of T-lymphotropic feline leukemia virus in vitro. Biochem. Biophys. Res. Commun.; 2006; 345, pp. 1466-1470. [DOI: https://dx.doi.org/10.1016/j.bbrc.2006.05.039]

50. Stewart, M.A.; Warnock, M.; Wheeler, A.; Wilkie, N.; Mullins, J.I.; Onions, D.E.; Neil, J.C. Nucleotide sequences of a feline leukemia virus subgroup A envelope gene and long terminal repeat and evidence for the recombinational origin of subgroup B viruses. J. Virol.; 1986; 58, pp. 825-834. [DOI: https://dx.doi.org/10.1128/jvi.58.3.825-834.1986]

51. Roy-Burman, P. Endogenousenv elements: Partners in generation of pathogenic feline leukemia viruses. Virus Genes; 1995; 11, pp. 147-161. [DOI: https://dx.doi.org/10.1007/BF01728655]

52. Pandey, R.; Ghosh, A.K.; Kumar, D.V.; Bachman, B.A.; Shibata, D.; Roy-Burman, P. Recombination between feline leukemia virus subgroup B or C and endogenous env elements alters the in vitro biological activities of the viruses. J. Virol.; 1991; 65, pp. 6495-6508. [DOI: https://dx.doi.org/10.1128/jvi.65.12.6495-6508.1991]

53. Stewart, H.; Jarrett, O.; Hosie, M.J.; Willett, B.J. Are endogenous feline leukemia viruses really endogenous?. Vet. Immunol. Immunopathol.; 2011; 143, pp. 325-331. [DOI: https://dx.doi.org/10.1016/j.vetimm.2011.06.011] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/21715022]