Email updates

Keep up to date with the latest news and content from BMC Veterinary Research and BioMed Central.

Open Access Research article

Study of compartmentalization in the visna clinical form of small ruminant lentivirus infection in sheep

Hugo Ramírez15, Ramsés Reina1, Luigi Bertolotti2, Amaia Cenoz1, Mirna-Margarita Hernández1, Beatriz San Román1, Idoia Glaria1, Ximena de Andrés1, Helena Crespo1, Paula Jáuregui1, Julio Benavides3, Laura Polledo4, Valentín Pérez4, Juan F García-Marín4, Sergio Rosati2, Beatriz Amorena1 and Damián de Andrés16*

  • * Corresponding author: Damián de Andrés ancad@unavarra.es

  • † Equal contributors

Author Affiliations

1 Instituto de Agrobiotecnología, CSIC-UPNA-Gobierno de Navarra, 31192 Mutilva, Navarra, Spain

2 Dipartimento di Produzioni Animali, Epidemiologia, Ecologia, Facoltá di Medicina Veterinaria, Universitá degli Studi di Torino, Grugliasco (TO), Italy

3 Instituto de Ganadería de Montaña (CSIC-ULE), León, Spain

4 Facultad de Veterinaria, Universidad de León, León, Spain

5 Laboratorio de Virología, Genética y Biología Molecular, FESC-UNAM, Veterinary C-4, 54700 Cuautitlán Izcalli, Estado de México, Mexico

6 Instituto de Agrobiotecnología, CSIC-UPNA-Gobierno de Navarra, Ctra Mutilva s/n, 31192 Mutilva, Navarra, Spain

For all author emails, please log on.

BMC Veterinary Research 2012, 8:8  doi:10.1186/1746-6148-8-8

The electronic version of this article is the complete one and can be found online at: http://www.biomedcentral.com/1746-6148/8/8


Received:5 July 2011
Accepted:26 January 2012
Published:26 January 2012

© 2012 Ramirez et al; BioMed Central Ltd.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

Background

A central nervous system (CNS) disease outbreak caused by small ruminant lentiviruses (SRLV) has triggered interest in Spain due to the rapid onset of clinical signs and relevant production losses. In a previous study on this outbreak, the role of LTR in tropism was unclear and env encoded sequences, likely involved in tropism, were not investigated. This study aimed to analyze heterogeneity of SRLV Env regions - TM amino terminal and SU V4, C4 and V5 segments - in order to assess virus compartmentalization in CNS.

Results

Eight Visna (neurologically) affected sheep of the outbreak were used. Of the 350 clones obtained after PCR amplification, 142 corresponded to CNS samples (spinal cord and choroid plexus) and the remaining to mammary gland, blood cells, bronchoalveolar lavage cells and/or lung. The diversity of the env sequences from CNS was 11.1-16.1% between animals and 0.35-11.6% within each animal, except in one animal presenting two sequence types (30% diversity) in the CNS (one grouping with those of the outbreak), indicative of CNS virus sequence heterogeneity. Outbreak sequences were of genotype A, clustering per animal and compartmentalizing in the animal tissues. No CNS specific signature patterns were found.

Conclusions

Bayesian approach inferences suggested that proviruses from broncoalveolar lavage cells and peripheral blood mononuclear cells represented the common ancestors (infecting viruses) in the animal and that neuroinvasion in the outbreak involved microevolution after initial infection with an A-type strain. This study demonstrates virus compartmentalization in the CNS and other body tissues in sheep presenting the neurological form of SRLV infection.

Keywords:
Compartmentalization; Visna; Small ruminant lentivirus; Spinal cord; Choroid plexus; Sheep

Background

Small ruminant lentiviruses (SRLVs) cause arthritis, mastitis, interstitial pneumonia (Maedi) and leukoencephalitis (Visna) in sheep and goats [1]. The clinical form Visna in sheep was first described in Iceland during the Visna Maedi Virus (VMV) epidemic in the first half of the last century [2]. From then, this neurologic form has been reported only sporadically [3-5]. However, in the last few years, the neurological form of SRLV disease has been diagnosed in numerous sheep of North-Western Spain, causing production losses [6]. In this outbreak, clinical signs usually appear at 1-2 years of age, but have been detected in animals as young as 4 months [7]. However, in other geographic areas of the world, Visna appears in animals over 2 years of age [2,8]. Also, in the animals affected by Visna within the outbreak, lesions frequently involve the white matter of the spinal cord [9,10] whereas in reported experimental Visna cases, lesions usually appear mainly in the periventricular areas of the brain [5].

Central nervous system (CNS) infection may occur in the early phases of sheep SRLV infection, as it occurs in human CNS with human immunodeficiency virus infections (HIV; [11]). In HIV infections, viral quasispecies change in their adaptation to certain cell types, due to differences in selective pressure that may be exerted, for example, by the immune system [12]. Consequently, the flow of genes among viral subpopulations is significantly restricted, leading to genetically distinct subpopulations and development of compartmentalization [13-15]. Hence, the genetic heterogeneity between subpopulations of the individual could be explained by both an independent microevolution and the presence of related but phylogenetically distinct infecting virus genotypes [15,16]. Viral populations distributed in compartments may have different phenotypic characteristics, such as cell tropism [17], drug resistance [18-21] and pathogenicity [22]. Further studies on viral diversification, compartmentalization and adaptation of the virus to the brain are needed [13]. The CNS provides a unique environment for the replication of lentiviruses. The blood-brain barrier presents specialized capillary endothelial cells, linked to each other by highly selective narrow bridges that may restrict viral traffic. During viral infection, a single strain appears to be able to cross the blood-brain barrier and establish the infection in the brain, which is otherwise isolated from the peripheral blood and the immune response [14,23]. Viral genome structure may also contain key features determining the fate of infection.

In SRLV infection, evidence of the relationship between SRLV-derived clinical manifestation and strain genetic features has been described [24], but there is no consensus on the specific viral genetic region that determines cell or tissue tropism [25]. Tissue tropism may be dependent on the virus promoter sequence of both VMV and HIV [26,27]. However, sequence divergence found in viral promoter does not ensure changes in transcriptional activities and/or biological characteristics. LTR alone does not exert this function and other viral genes could be involved [25].

Different HIV and feline immunodeficiency virus (FIV) studies, based on a hypervariable region of the viral surface protein (SU), have revealed tissue compartmentalization of the virus [14]. In HIV and FIV Env, the V3 region is a determinant of cell tropism and replication efficiency, since it is thought to be related to the adsorption and fusion of the virus to the cell [23,28-30]. Furthermore, in chronically infected individuals, particular amino acids of the V3 region modulate neurotropism and neurovirulence [14]. In SRLV, five variable (V1 to V5) and four conserved (C1 to C4) regions have been identified in the SU protein [31] and the function of V4 hypervariable region has been found analogous to that of the V3 region of HIV [32,33].

In this work involving Env SU regions (V4, C4 and V5), and the amino terminal sequence of the Env TM protein, we first determined the existence of viral sequence diversity in the CNS cells [choroid plexus cells (CPx), spinal cord cells (SC)], broncho-alveolar lavage macrophages (BAL), lung tissue (L), mammary gland tissue (MG), and peripheral blood mononuclear cells (PBMC) and analyzed the phylogenetic relationships of SRLV within animals belonging to an outbreak of clinical Visna [9]. Phylogeography approaches were also applied to identify the common ancestor sequences of the viral quasispecies. Finally, the existence of compartmentalization and/or presumptive positive selection of viral genomes was assessed, with special emphasis on CNS comparisons in order to better understand genetic diversity involved neurotropism and neurovirulence.

Results

Clinical signs and lesions

The animals included in this study presented symptoms compatible with the nervous form of VMV, such as hindleg ataxia or even paralysis, sternal or lateral recumbency and pedaling. Histopathological examination indicated the presence of severe lesions in the CNS in all the animals. In three of these animals the CNS lesions were located only in the SC, which were characterized by non-suppurative perivascular cuffs and demyelination, when affecting the neuroparenchyma. Seven of the animals showed also mild interstitial pneumonia, denoted by lymphoid infiltration in the lung alveolar walls. Only two of the animals had lesions in the MG, characteristic of a mild non-purulent interstitial mastitis (Table 1). No other pathogens were identified in the brain or other organs in the animals under study.

Table 1. Diagnosis and lesions. Results on PCR (LTR region), ELISA (Elitest) and histopathology obtained from eight sheep with Visna clinical symptoms belonging to the Castilla-y-León outbreak

ELISA and PCR

All the animals had anti-SRLV antibodies as assessed by Elitest [34], with OD ratios ranging from 3 to 18 (Table 1). In spite of the availability of PBMC, BAL, MG (except for the male), L, CPx and SC samples from all the animals for PCR testing, only some samples allowed proviral amplification (Table 1). PCR amplification detected proviral LTR DNA mostly in tissues of the CNS (SC and CPx), whereas samples from L or MG yielded few positives and inconsistent results, respectively. LTR SRLV sequences obtained in this study had a classical genotype A organization and characteristics.

The env PCR amplification was successful in all the animals, SC being always amplified (data not shown). BAL and CPx yielded PCR positive results in 4 of the 8 animals. The number of clones was also variable among samples and animals.

Sequence and phylogeny studies

A representative sample of the env amplicons was sequenced and genotyped, revealing that, like in LTR, env SRLV sequences obtained in this study had a classical genotype A characteristics, closer to VMV than to caprine arthritis encephalitis virus (CAEV) prototypes (not shown). Within each animal, the mean diversity of env sequences from CNS ranged from 0.35 to 11.6% (Table 2). The analysis of nucleotide diversity showed a quite uniform distribution of distances, except for the animal No. 333 having two distant sequence sub-sets (333a and 333b) in different tissues, including CNS. One of them grouped with those of the outbreak but not the other (30% diversity between both sequence sub-sets). Diversity observed among sequences from the CNS in different animals ranged between 11.1 and 16.1%.

Table 2. Nucleotide and amino acid diversity, ω statistic and ancestor tissue. Nucleotide diversity, amino acid diversity, ω statistic (ratio between average number of non-synonymous substitutions per non-synonymous site, dN; and average number of synonymous substitutions per synonymous site, dS), and ancestor tissue for SRLV clones obtained from eight sheep of the Visna outbreak

Phylogenetic analysis of SRLV sequences within animals is shown in Figure 1 and globally in Additional file 1: 'Phylogenetic tree including 350 partial env sequences from 8 animals'. With the exception of animal 333, sequences belonging to the same animal clearly clustered together. In addition, within the animal, sequences showed a strong tissue compartmentalization, confirming the presence of different viral sub-populations in a single host. This was compatible with either a clonal evolution after initial infection with a single strain or a polyclonal origin of infecting viral strains. Sequence distribution in tissues differed between animals. In addition, MG and L sequences frequently formed clearly defined clades, while PBMC sequences were the most variable and were usually located in different clades within an animal (Figure 1). Individual phylogram study revealed the existence of sequences in CNS, L, and/or MG that could derive from BAL and/or PBMC. In animals 166 and 368, sequences from CNS (Cpx and SC) clustered close to each other. In both animals, these sequences, when present, were also close to BAL, like SC sequences from animal 336. In fact, in sheep 166, both CNS and BAL sequences belonging to the late clade appeared to descend from BAL, L and/or PBMC (early clade). In sheep 333, there were differences in sequence distribution between the two sequence sub-sets. In one of them BAL and PBMC sequences were close to those from CNS (SC and CPx) and MG; and the other revealed the proximity of PBMC sequences to CNS (CPx), MG and L. In animal 223 CPx sequences belonged to a cluster different from that of SC sequences of the CNS, confirming the existence of different viral sequence clusters within the CNS, and PBMC sequences were close to the tree root. Also, sequences from a particular tissue appeared to have evolved within the tissue; this was the case of MG sequences in animals 336 and 698 or, less markedly SC sequences in animal 697, among others. Thus, once in the target organ, the virus seemed to have evolved, allowing the appearance of new clones within the same tissue.

thumbnailFigure 1. Phylogenetic relationships among sequences belonging to the same animal. Animal identification number is reported for each subtree. Taxa names include the tissue from which the sequence was obtained and the clone number. Tissues are coded as follow: broncho-alveolar macrophages, BAL; choroid plexus, CPx; spinal cord, SC; peripheral blood mononuclear cells, PBMC; mammary gland, MG; lung, L. Posterior probabilities of each node are reported above branches. Sequences from animal 698 with stop codons are reported in grey.

Additional file 1. Phylogenetic tree including 350 partial env sequences from 8 animals. Taxa names include the animal number, the tissue from which the sequence was obtained and the clone number. Tissues are identified by a three-letter code as follows: broncho-alveolar macrophages BAL, choroid plexus CPx, spinal cord SCO, peripheral blood mononuclear cells PBM, mammary gland MGL, lung LUN. Animal clusters are depicted in different colors. Sequences presenting stop codons within the 698 animal clade are reported in gray. Posterior probability supporting tree's nodes are reported.

Format: PDF Size: 1.7MB Download file

This file can be viewed with: Adobe Acrobat ReaderOpen Data

The common ancestor state was then reconstructed, using Bayesian methods recently developed in phylogeography [35] (Table 2). Evolutionary parameters were estimated under the Bayesian Skyline Plot uncorrelated lognormal relaxed clock model. The convergence of chains was reached (ESS > 200 in all analyses) and Maximum Clade Credibility trees were obtained for each animal. The results showed that viral evolution within animals followed the usual infection paths. Specifically, in animals 166, 336 and 368, where sequences from BAL were present, bronchoalveolar macrophages were identified as the most likely tissue infected by the viral common ancestor. In animals 223 and 333, in which BAL samples were unavailable, PBMC had the highest probability to be the tissue infected by the viral ancestor. These data strongly suggest that the infection starts from the lungs and later reaches other tissues through the blood. An interesting aspect of these results is that CNS, lung and mammary gland are terminal clades in this evolutionary model, associated to BAL and PBMC taxa (Table 2).

In this evolutionary context, the average number of non-synonymous substitutions per non-synonymous site (dN) and the number of synonymous substitutions per synonymous site (dS) were determined for proviral sequences obtained from each animal. The dN rates were lower than the dS rates in all the animals, indicating the presence of within-host purifying selection (against amino acid sequence variability), but in favour of synonimous substitutions leading to nucleotide heterogeneity and compartmentalization (Table 2).

Compartmentalization regarding "segregation" of sequences between different animals and between tissues of a particular animal was assessed (Tables 3 and 4) according to parsimony score (PS), association index (AI) and monophyletic clade (MC) values. Results on PS and AI showed weak phylogeny-trait association and significant monophyletic clade (MC) value, altogether indicative of tissue compartmentalization in six of the studied animals. In particular, CNS sequences belonged to a different cluster compared to other tissues, within the animal. This was observed in the six animals that had at least two tissues (one of them being CNS) from which the number of sequences obtained was sufficient (≥ 28 clones per animal) to carry out the study, demonstrating and confirming the existence of CNS tissue compartmentalization. Lack of significance in PS, AI and MC analyses was observed only rarely, when a low number of sequences was available (animals 292 and 697; Table 4).

Table 3. SRLV sequences segregation. SRLV sequence compartmentalization (segregation) between the eight Visna-affected animals, using Bayesian MCMC approach for determination of Association Index (AI), Parsimony Score (PS) and Monophyletic Clade (MC) values

Table 4. Tissue compartmentalization of SRLV sequences. SRLV sequence compartmentalization between tissues within the animal (n = 7)*, using Bayesian MCMC approach for determination of Association Index (AI), Parsimony Score (PS) and Monophyletic Clade (MC) values

Overall, these results indicate strong sequence segregation between animals and the existence of tissue compartmentalization within the animal.

A signature pattern analysis was performed in search of conserved amino acid motifs specific to a particular tissue, but we were unable to find this pattern in the genetic region under analysis. Also analysis of the potential Asn-X-Ser/Thr glycosylation site was carried out [36]. This site was present in sequences from all the animals (except No. 223) and was equally distributed in different tissues, discarding glycosylation pattern as a possible explanation for compartmentalization.

Discussion

This study provides evidence of VMV compartmentalization within the CNS and other tissues and confirms its existence in the mammary gland as previously described between blood and colostrum [32] amongst seropositive naturally infected small ruminants. The sheep had Visna-like severe lesions in the CNS, with non-suppurative meningo-encephalitis [10]. According to the specific feature of this outbreak [6,7], most animals presented lesions in the spinal cord and 50 percent in the brain. Most of them also presented lung lesions. VMV sequence distribution differed between animals, as revealed by phylogenetic analysis. The different compartmentalization patterns may be explained by differences in host genetic susceptibility, anatomical features, the genome of the infecting viruses, as well as high rates of mutation and replication, generating a quasispecies [37], and variants which may escape from surveillance [38]. Compartmentalization may involve genetic drift and founder effects [23], and/or being a consequence of different selective pressures [12].

Respiratory secretions containing infected alveolar macrophages can be responsible for SRLV transmission, especially in intensive rearing systems such as those applied in the Assaf breed of the animals under study [39]. Consequently, blood, lymph, and then target tissues, such as L, MG, CPx and SC of CNS, become infected [40]. Following purifying selection (as observed in this study) the virus evolves and compartmentalizes, producing a quasispecies. Thereafter, the virus may go back through lymphoid and blood cells to other tissues. In line with this scheme and using a novel methodology previously applied to investigate geographical origin of sequences [35], we found as expected, that intra-host ancestor sequences were not located in the CNS or other tissues, but in BAL or PBMC (when available). PBMC sequences were present in different clades, displaying the broadest variation per tissue, which is compatible with migration of virus from blood to tissues and vice versa.

However, particular ancestor/primary sequences could not be always identified. This may be due to the unavailability of PCR amplification in some samples. Thus, clones analyzed in the ancestor tissues (BAL and PBMC) may have not included the complete spectrum of ancestor viral sequences. Besides, all samples were obtained at a single time point (necropsy) even though the time points of infection by a particular virus may have not been simultaneous at different body sites. Also, there may have been re-infection events from inside or outside the body (like sheep 333, apparently with two main infections). Finally, body sites other than BAL (such as nasal and conjunctival mucosa) may have been the initial source of infection.

The virus under study belonged to the genotype A, also involved in other Visna infections [41]. Like in CNS studies of HIV [42] and in contrast with findings on FIV [14] infections, we were unable to find any CNS-specific signature pattern. In these lentiviruses, V3-V5 regions of SU encoding neutralizing epitopes may mutate, leading to viral resistance to existing antibodies also affecting cell tropism [18,30], including neurotropism and neurovirulence [14,43]. Amongst the SRLV sequences under study (TM and V4-V5), V4 - claimed to have a function analogous to that of HIV-V3 region [32,33] - was the most variable. Although variable and constant regions within the SRLV SU protein may trigger the production of neutralizing antibodies [31], the main sites/epitopes involved in this process have not been identified. Additional factors or genetic regions, such as those encoding viral proteins (Vif) that interact with host molecules (APOBEC; [44]), might be involved in neuroinvasion.

Of the 16 MG clones of animal 698, 12 showed a non-synonymous transition C to T at the same position, inserting a stop codon along the env region encoding the SU protein (Figure 1). The finding of stop codons along with codifying sequences suggests the emergence of quasispecies, the presence of a high rate of viral evolution and thus the existence of defective interfering particles/heteroclites. However, numerically realistic model studies suggest that these particles are unlikely to survive or influence viral dynamics [45].

Work on HIV and simian immunodeficiency virus (SIV) has shown that viral evolution in the brain is peculiar in that this immunologically privileged site allows a low level of viral replication [13] before disease onset [46], but it may be considered a "sanctuary" where the virus can evolve freely, avoiding the effects of adaptive immunity and/or therapy [12]. In line with this, the least conserved Env amino acid sequence was found in SC and the most conserved in L (data not shown). The diagnostic LTR-PCR used in this study only yielded consistent amplifications when using CNS as DNA source, since all animals were seropositive and yielded PCR-specific amplicons in CNS (Table 1). This strongly suggests that CNS is a preferred site of infection in the Visna affected animals involved. SRLV replication occurred in brain, since viral antigens p25 and gp130 were detected in CNS by immunohistochemistry and two fully replicative viral isolates, able to replicate in SCP cells, were obtained from SC and cerebrospinal fluid from two of the diseased sheep used (Nos. 697 and 698, respectively) [47]. This is compatible with studies on HIV infections, where viral replication in the CNS occurs in diseased individuals, leading then to divergent trends in different compartments [28].

This study shows that different SRLVs can coexist in the CNS from a Visna diseased sheep. This was clearly observed in animal 333, infected by two viruses (333a and 333b sequences), one closer to the remaining outbreak sequences than the other. Under the hypothesis that the strain's genetic makeup (in the SU genetic regions analyzed) and induced pathology are associated, the set of sequences most similar to other sequences in the outbreak would be causing the nervous clinical signs. The CNS 333a sub-set sequences were essentially from CPx (n = 7) and to a lower extent from SC (n = 2), whereas 333b ones were similarly represented in CPx (n = 9) and SC (n = 11). This would be consistent with a recent infection in the case of 333a, in which the virus from infected CPx cells (blood-CSF barrier) reached only a few SC sites. In the case of 333b sequences, the infection appeared to have been longer established, as it extended to various SC sites after CPx initial infection. Whether different viruses infected different cells of the same tissue (such as choroid plexus) or different viruses coinfected the same cell type remains unknown. In any case, this study demonstrates in VMV infections that, like in those by SIV, the physical continuity of the brain tissue does not necessarily result in phylogenetic proximity [13]. Disruption of the blood-brain or the blood-CSF (CPx) barriers may have taken place, as proposed in FIV and HIV infections [43,48]. In HIV infections, CPx may allow a bidirectional productive lentiviral infection of CNS and peripheral organs, having often a mixture of viral sequences [43,46]. Similarly in this study, VMV-infected cells may have crossed the CPx barrier to colonize CSF and SC, and gone back from there to blood through blood-brain or blood-CSF (CPx) barriers.

Conclusions

The results obtained in this work provide evidence of type A SRLV compartmentalization in the CNS and other tissues. These viruses can be found in tissues that allow horizontal transmission among the animals of the outbreak. Ancestor tissue estimation strongly suggests that infection appears to start in alveolar macrophages (BAL) or blood (PBMC) when BAL was unavailable. Likely, PBMC became infected after intake of virus or infected cells/particles through respiratory and/or mammary secretions, including colostrum/milk. Moreover, the identification of sequences from mammary gland (MG) as terminal tips in the trees supports this hypothesis, indicating the two main compartments from which the virus can be transmitted.

Through this route, circulating A viruses cause neuropathological disease in the Assaf animals under study. Subsequently and according to phylogeny observations, viruses may either reach different tissues including the CNS, lung, and mammary gland and compartmentalize or be transmitted from ancestor tissues (BAL/PBMC) to new hosts.

Methods

Animals and diagnostic tests

This study involved seven sheep of the Assaf breed - six females (223, 292, 333, 336, 697 and 698) and one male (166) - and a Milchschaf breed ewe (368), ranging in age from 2 to 5 years. Sheep were from different dairy farms (except for 333 and 336, which were from the same flock) located in the Autonomous Community of Castilla-y-León (North-Western Spain), where several cases of Visna Maedi nervous form had been previously reported [9]. All sheep showed clinical nervous symptoms and post-mortem analysis was consistent with Visna lesions (Table 1).

Animal handling, euthanasia and experimental procedures were carried out in compliance with the current European and national (RD 1201/2005) regulations, with the approval of the Comité de Ética y Experimentación Animal of the University of León and authorization of the Gobierno de Castilla y León. Euthanasia was done by intravenous injection of barbiturate overdose followed by exsanguination. Complete necropsy analysis was performed in the animals. Animals were submitted to necropsy and used for clinical and pathology studies at the University of León following the farmers' request and consent.

For serological and molecular diagnosis of SRLV infection, animals were analyzed serologically using a commercially available ELISA (Elitest-MVV®; HYPHEN Biomed, France); and by PCR on viral promoter LTR (about 300 bp), using primers and a procedure previously described [49] for proviral DNA amplification.

Samples

Whole blood from each animal was collected to isolate plasma for ELISA and peripheral blood mononuclear cells (PBMC) for PCR analysis. These were isolated from EDTA-blood samples by Ficoll-Hypaque gradient centrifugation (d = 1.077; Lymphoprep Axis-Shield, Oslo). One cubic centimetre samples of spinal cord (SC), choroid plexus (CPx), mammary gland (MG), lung (L), and bronchoalveolar lavage (BAL) were also collected post-mortem and embedded in RNAlater (Qiagen) until use. Samples from brain, SC, L and MG were also obtained and fixed in buffered formalin 10% for histopathological studies. Genomic DNA was extracted from PBMC, BAL and tissue samples with QIAamp® DNA Blood Mini Kit (Qiagen), following the manufacturer's instructions. Tissue samples were previously lysed in buffer (100 mM Tris-HCl pH 8.5; 5 mM EDTA; 400 mM NaCl; 0.2% SDS) containing Proteinase K (50 μg/ml per 107 PBMC, Sigma) at 56°C for 1 h.

PCRs and cloning for compartmentalization

To avoid template resampling [50], DNA was diluted in fivefold series and amplified by PCR in triplicate (limiting-dilution PCRs). Diluted DNA was used in a seminested PCR that was carried out with primers previously described [51] slightly modified for this study. Briefly, the reaction mix of the first round of seminested PCR consisted in: Buffer Certamp 1 × (Biotools), 2 mM MgCl2 (Biotools), 225 μM of each dNTP (Biotools), 600 nM of each primer, 0.04 U/μl of Certamp enzyme (Biotools) and sample DNA diluted 5-fold serially in a final volume of 50 μl. PCR conditions were: initial denaturation step for 5 min at 95°C followed by 45 cycles of 95°C for 60 s, annealing at 47°C for 60 s and extension at 72°C for 60 s, followed by a final extension step of 10 min at 72°C. Primers were: #563 FW 5'-GAYATGRYRGARCAYATGACA-3' (7272-7292 nt); #564 RV 5'-GCYAYRTGCTGYACCATGGCATA-3' (8089-8067 nt). For the second round, 5 μl of the first PCR were transferred to the new reaction mix. Reaction conditions for the second round of PCR were identical to those described for the first round, with the exception that the #567 FW primer (#567 FW 5'-GGGNACNARIACRAAYTGGAC-3'; 7481-7501 nt) was used and the annealing temperature was increased to 49°C. A 607 bp amplicon was obtained encompassing the C-terminal part of surface protein (SU) including V4, C4 and V5 regions and the N-terminal part of transmembrane (TM) protein encoded by the env gene (7481-8089 nt of CAEV CO, M33677). Amplified products were chilled at 4°C. A negative control reaction without template DNA was included for each primer pair to control potential contaminations. In order to minimize the risk of contamination, mixing of reagents, sample preparation, and analysis of the amplified products were performed in three separate areas.

PCR amplicons were resolved on 1.5% agarose gels, and bands of the expected sizes were excised and purified with QIAquick Gel Extraction Kit (Qiagen), and subsequently cloned into the vector pGEMT-easy® (Promega) following the manufacturer's instructions. Fifty to 100 transformed colonies (Escherichia coli XL1-Blue) were picked per sample and grown overnight at 37°C in LB medium (1% Bacto Tryptone, 0.5% yeast extract, 1% NaCl) with 100 μg/ml ampicillin. After selection of positive clones upon screening by restriction enzyme digestion (EcoRI), plasmid DNA was purified from bacterial culture using the Quantum Prep® Plasmid Miniprep kit (Bio-Rad) according to the manufacturer's instructions. Purified plasmid DNA was sequenced using an ABI PRISM 310 Genetic Analyzer (Applied Biosystems).

Sequence analysis

The total number of sequences obtained in the group of 8 animals was 350: 54 from CPx, 88 from SC, 39 from L, 46 from PBMC, 69 from BAL and 54 from MG.

Nucleotide sequences were analyzed with Chromas 2.23 (Technelysium, Helensvale, Australia) and edited using BioEdit software [52]. The corresponding amino acid sequences were obtained according to the Expasy server. Multiple alignments for comparison and phylogenetic studies were made with the ClustalW [53] respecting the reading frame. The nucleotide distance matrix was calculated with the PAUP* software [54]. The presence of natural selection was evaluated as the estimation of ω, the ratio between the number of non-synonymous substitutions per non-synonymous site (dN) and the number of synonymous substitutions per synonymous site (dS). The ω estimation was carried out using Single-Likelihood Ancestor Counting (SLAC, [55]) method implemented in Datamonkey webserver [56].

A phylogenetic tree was created evaluating the best model of molecular evolution (ModelTest software; [57]) and using Bayesian heuristic approaches (BEAST (16) and MrBayes software; [58]).

To explore the overall compartmentalization in the sequence data we employed a Bayesian Markov chain Monte Carlo (MCMC) approach (program BaTS; [59]). This analysis was based on the trees output from the MrBayes BEAST analysis described above, with 10% removed as burn-in and employing 1000 replications. From these trees we computed the significance of the parsimony score (PS) and association index (AI) statistics of the strength of geographical clustering by phylogeny [59]. In addition, we computed the monophyletic clade (MC) statistic that measures the strength of clustering within individual geographic regions. Thus, common ancestors were not inferred from the tree topology directly, since tree clusters were based only on minimal differences between the sequences compared (independently of the pathways used in diversification) whereas the ancestor determination applied was a complex process, determining the most parsimonious pathway to maximally relate the set of sequences under study while minimizing the phyletic diversity. Importantly, the common ancestor approach accounts for uncertainty in the underlying phylogeny by using a large number of plausible trees. In order to understand correctly the results, we have considered that low AI and PS values represent strong phylogeny-trait association and MC is positively correlated with the strength of the phylogeny-trait association [59]. Thus, high AIs and PSs and low MCs show sequences' compartmentalization. To better understand the clonal origin of collected sequences, Bayesian MCMC approaches described by Lemey and coauthors [35] and developed for phylogeographical applications were applied. Specifically, according to Lemey and collaborators, the state of each sequence was represented by its geographical origin and the analyses lead to infer about the state of each internal node: this approach allowed to follow the viral evolution through the space of assigned states. We used BEAST package software [60] to infer about the state of internal nodes within the phylogeny of each animal. In this work, we assigned as state the animal tissue from which each sequence was obtained. Parameters were estimated using the Bayesian method implemented in the BEAST software package and the convergence of chains was assessed by calculating the Effective Sample Size (ESS) of the runs. All parameter estimates showed significant ESS (> 200). Maximum Clade Credibility (MCC) phylogeny and information about common ancestors' state were obtained using Tree Annotator.

A signature pattern analysis was carried out to search for possible amino acid positions that would provide a conserved pattern within the SNC Env sequences relative to the PBMC or BAL Env sequences. The analysis was calculated with the VESPA software [61] that determines the frequency of an amino acid at a specific position and then determined whether there was a distinct pattern for one set of sequences.

Nucleotide sequences accession number

The SRLV-env sequence data generated in this study were deposited into GenBank database and are available under accession numbers: BAL 166 (JF288187-JF288213), Cpx 166 (JF288214-JF288223), L 166 (JF288224-JF288235), PBMC 166 (JF288236-JF288255), SC 166 (JF288256-JF288276), Cpx 223 (JF288277-JF288289), L 223 (JF288290-JF288293), MG 223 (JF288294-JF288299), PBMC 223 (JF288300-JF288304), SC 223 (JF288305-JF288308), SC 292 (JF288309-JF288315), BAL 333 (JF28831-JF288329), Cpx 333 (JF288330-JF288345), L 333 (JF288346-JF288353), MG 333 (JF288354-JF288359), PBMC 333 (JF288360-JF288376), SC 333 (JF288377-JF288389), BAL 336 (JF288390-JF288403), L 336 (JF288404-JF288415), MG 336 (JF288416-JF288440), SC 336 (JF288441-JF288453), BAL 368 (JF288454-JF288467), Cpx 368 (JF288468-JF288482), L 368 (JF288483-JF288485), MG 368 (JF288486), SC 368 (JF288487-JF288496), SC 697 (JF288497-JF288508), MG 698 (JF288509-JF288524), PBMC 698 (JF288525-JF288528), SC 698 (JF288529-JF288536).

Competing interests

The authors declare that they have no personal, financial or non financial competing interests.

Authors' contributions

HR, AC, M-MH, BS-R, IG, XA, HC and PJ carried out the immunoassays, molecular genetic studies and participated in the sequence alignment. LB performed the statistical analysis and helped to draft the manuscript. JB, LP, VP and J-F G-M carried out the histopathological study. HR, RA, BA and DA conceived of the study, and participated in its design and coordination and helped to draft the manuscript. All authors read and approved the final manuscript.

Acknowledgements

Authors, gratefully acknowledge Dr. Tony L. Goldberg for help in sequence analyses. This study was supported by grants from Spanish CICYT (AGL2007-66874-C04, AGL2010-22341-C04) and Government of Navarra Department of Agriculture (IIQ010449.RI1 and IIQ14064.RI1). H. Ramirez was supported by PASPA program, UNAM and Public University of Navarra. L. Bertolotti acknowledges financial support by Ricerca Sanitaria Finalizzata 2009-Regione Piemonte.

References

  1. de Andres D, Klein D, Watt NJ, Berriatua E, Torsteinsdottir S, Blacklaws BA, Harkiss GD: Diagnostic tests for small ruminant lentiviruses.

    Vet Microbiol 2005, 107:49-62. PubMed Abstract | Publisher Full Text OpenURL

  2. Sigurdsson B, Palsson PA: Visna of sheep; a slow, demyelinating infection.

    Br J Exp Pathol 1958, 39:519-528. PubMed Abstract | PubMed Central Full Text OpenURL

  3. Dawson M: Pathogenesis of maedi-visna.

    Vet Rec 1987, 120:451-454. PubMed Abstract | Publisher Full Text OpenURL

  4. Watt NJ, Roy DJ, McConnell I, King TJ: A case of visna in the United Kingdom.

    Vet Rec 1990, 126:600-601. PubMed Abstract OpenURL

  5. Zink MC, Johnson LK: Pathobiology of lentivirus infections of sheep and goats.

    Virus Res 1994, 32:139-154. PubMed Abstract | Publisher Full Text OpenURL

  6. Benavides J, Gomez N, Gelmetti D, Ferreras MC, Garcia-Pariente C, Fuertes M, Garcia-Marin JF, Perez V: Diagnosis of the nervous form of Maedi-Visna infection with a high frequency in sheep in Castilla y Leon, Spain.

    Vet Rec 2006, 158:230-235. PubMed Abstract | Publisher Full Text OpenURL

  7. Benavides J, Garcia-Pariente C, Ferreras MC, Fuertes M, Garcia-Marin JF, Perez V: Diagnosis of clinical cases of the nervous form of Maedi-Visna in 4- and 6-month-old lambs.

    Vet J 2007, 174:655-658. PubMed Abstract | Publisher Full Text OpenURL

  8. Pritchard GC, Done SH, Dawson M: Multiple cases of maedi and visna in a flock in East Anglia.

    Vet Rec 1995, 137:443. PubMed Abstract OpenURL

  9. Benavides J, Fuertes M, Garcia-Pariente C, Ferreras MC, Garcia Marin JF, Perez V: Natural cases of visna in sheep with myelitis as the sole lesion in the central nervous system.

    J Comp Pathol 2006, 134:219-230. PubMed Abstract | Publisher Full Text OpenURL

  10. Benavides J, Garcia-Pariente C, Fuertes M, Ferreras MC, Garcia-Marin JF, Juste RA, Perez V: Maedi-visna: the meningoencephalitis in naturally occurring cases.

    J Comp Pathol 2009, 140:1-11. PubMed Abstract | Publisher Full Text OpenURL

  11. McArthur JC, Brew BJ, Nath A: Neurological complications of HIV infection.

    Lancet Neurol 2005, 4:543-555. PubMed Abstract | Publisher Full Text OpenURL

  12. Caragounis EC, Gisslen M, Lindh M, Nordborg C, Westergren S, Hagberg L, Svennerholm B: Comparison of HIV-1 pol and env sequences of blood, CSF, brain and spleen isolates collected ante-mortem and post-mortem.

    Acta Neurol Scand 2008, 117:108-116. PubMed Abstract | Publisher Full Text OpenURL

  13. Chen MF, Westmoreland S, Ryzhova EV, Martin-Garcia J, Soldan SS, Lackner A, Gonzalez-Scarano F: Simian immunodeficiency virus envelope compartmentalizes in brain regions independent of neuropathology.

    J Neurovirol 2006, 12:73-89. PubMed Abstract | Publisher Full Text OpenURL

  14. Liu P, Hudson LC, Tompkins MB, Vahlenkamp TW, Meeker RB: Compartmentalization and evolution of feline immunodeficiency virus between the central nervous system and periphery following intracerebroventricular or systemic inoculation.

    J Neurovirol 2006, 12:307-321. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  15. Zarate S, Pond SL, Shapshak P, Frost SD: Comparative study of methods for detecting sequence compartmentalization in human immunodeficiency virus type 1.

    J Virol 2007, 81:6643-6651. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  16. Philpott S, Burger H, Tsoukas C, Foley B, Anastos K, Kitchen C, Weiser B: Human immunodeficiency virus type 1 genomic RNA sequences in the female genital tract and blood: compartmentalization and intrapatient recombination.

    J Virol 2005, 79:353-363. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  17. Thompson KA, Churchill MJ, Gorry PR, Sterjovski J, Oelrichs RB, Wesselingh SL, McLean CA: Astrocyte specific viral strains in HIV dementia.

    Ann Neurol 2004, 56:873-877. PubMed Abstract | Publisher Full Text OpenURL

  18. Si-Mohamed A, Kazatchkine MD, Heard I, Goujon C, Prazuck T, Aymard G, Cessot G, Kuo YH, Bernard MC, Diquet B, et al.: Selection of drug-resistant variants in the female genital tract of human immunodeficiency virus type 1-infected women receiving antiretroviral therapy.

    J Infect Dis 2000, 182:112-122. PubMed Abstract | Publisher Full Text OpenURL

  19. Strain MC, Letendre S, Pillai SK, Russell T, Ignacio CC, Gunthard HF, Good B, Smith DM, Wolinsky SM, Furtado M, et al.: Genetic composition of human immunodeficiency virus type 1 in cerebrospinal fluid and blood without treatment and during failing antiretroviral therapy.

    J Virol 2005, 79:1772-1788. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  20. Tashima KT, Flanigan TP, Kurpewski J, Melanson SM, Skolnik PR: Discordant Human Immunodeficiency Virus Type 1 drug resistance mutations, including K103N, observed in cerebrospinal fluid and plasma.

    Clin Infect Dis 2002, 35:82-83. PubMed Abstract | Publisher Full Text OpenURL

  21. Venturi G, Catucci M, Romano L, Corsi P, Leoncini F, Valensin PE, Zazzi M: Antiretroviral resistance mutations in human immunodeficiency virus type 1 reverse transcriptase and protease from paired cerebrospinal fluid and plasma samples.

    J Infect Dis 2000, 181:740-745. PubMed Abstract | Publisher Full Text OpenURL

  22. Donaldson YK, Bell JE, Holmes EC, Hughes ES, Brown HK, Simmonds P: In vivo distribution and cytopathology of variants of human immunodeficiency virus type 1 showing restricted sequence variability in the V3 loop.

    J Virol 1994, 68:5991-6005. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  23. Korber BT, Kunstman KJ, Patterson BK, Furtado M, McEvilly MM, Levy R, Wolinsky SM: Genetic differences between blood- and brain-derived viral sequences from human immunodeficiency virus type 1-infected patients: evidence of conserved elements in the V3 region of the envelope protein of brain-derived sequences.

    J Virol 1994, 68:7467-7481. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  24. Angelopoulou K, Brellou GD, Greenland T, Vlemmas I: A novel deletion in the LTR region of a Greek small ruminant lentivirus may be associated with low pathogenicity.

    Virus Res 2006, 118:178-184. PubMed Abstract | Publisher Full Text OpenURL

  25. Murphy B, McElliott V, Vapniarsky N, Oliver A, Rowe J: Tissue tropism and promoter sequence variation in caprine arthritis encephalitis virus infected goats.

    Virus Res 2010, 151:177-184. PubMed Abstract | Publisher Full Text OpenURL

  26. Agnarsdottir G, Thorsteinsdottir H, Oskarsson T, Matthiasdottir S, Haflidadottir BS, Andresson OS, Andresdottir V: The long terminal repeat is a determinant of cell tropism of maedi-visna virus.

    J Gen Virol 2000, 81:1901-1905. PubMed Abstract | Publisher Full Text OpenURL

  27. Murphy B, Jasmer DP, White SN, Knowles D: Localization of a TNF-activated transcription site and interactions with the gamma activated site within the CAEV U3 70 base pair repeat.

    Virology 2007, 364:196-207. PubMed Abstract | Publisher Full Text OpenURL

  28. Abbate I, Cappiello G, Longo R, Ursitti A, Spano A, Calcaterra S, Dianzani F, Antinori A, Capobianchi MR: Cell membrane proteins and quasispecies compartmentalization of CSF and plasma HIV-1 from aids patients with neurological disorders.

    Infect Genet Evol 2005, 5:247-253. PubMed Abstract | Publisher Full Text OpenURL

  29. Del Prete GQ, Leslie GJ, Haggarty B, Jordan AP, Romano J, Hoxie JA: Distinct molecular pathways to X4 tropism for a V3-truncated human immunodeficiency virus type 1 lead to differential coreceptor interactions and sensitivity to a CXCR4 antagonist.

    J Virol 2010, 84:8777-8789. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  30. Motokawa K, Hohdatsu T, Imori A, Arai S, Koyama H: Mutations in feline immunodeficiency (FIV) virus envelope gene V3-V5 regions in FIV-infected cats.

    Vet Microbiol 2005, 106:33-40. PubMed Abstract | Publisher Full Text OpenURL

  31. Valas S, Benoit C, Baudry C, Perrin G, Mamoun RZ: Variability and immunogenicity of caprine arthritis-encephalitis virus surface glycoprotein.

    J Virol 2000, 74:6178-6185. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  32. Pisoni G, Moroni P, Turin L, Bertoni G: Compartmentalization of small ruminant lentivirus between blood and colostrum in infected goats.

    Virology 2007, 369:119-130. PubMed Abstract | Publisher Full Text OpenURL

  33. Skraban R, Matthiasdottir S, Torsteinsdottir S, Agnarsdottir G, Gudmundsson B, Georgsson G, Meloen RH, Andresson OS, Staskus KA, Thormar H, Andresdottir V: Naturally occurring mutations within 39 amino acids in the envelope glycoprotein of maedi-visna virus alter the neutralization phenotype.

    J Virol 1999, 73:8064-8072. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  34. Saman E, Van Eynde G, Lujan L, Extramiana B, Harkiss G, Tolari F, Gonzalez L, Amorena B, Watt N, Badiola J: A new sensitive serological assay for detection of lentivirus infections in small ruminants.

    Clin Diagn Lab Immunol 1999, 6:734-740. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  35. Lemey P, Rambaut A, Drummond AJ, Suchard MA: Bayesian phylogeography finds its roots.

    PLoS Comput Biol 2009, 5:e1000520. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  36. Pisoni G, Bertoni G, Manarolla G, Vogt HR, Scaccabarozzi L, Locatelli C, Moroni P: Genetic analysis of small ruminant lentiviruses following lactogenic transmission.

    Virology 2010, 407:91-99. PubMed Abstract | Publisher Full Text OpenURL

  37. Domingo E, Gomez J: Quasispecies and its impact on viral hepatitis.

    Virus Res 2007, 127:131-150. PubMed Abstract | Publisher Full Text OpenURL

  38. Martin-Garcia J, Cao W, Varela-Rohena A, Plassmeyer ML, Gonzalez-Scarano F: HIV-1 tropism for the central nervous system: brain-derived envelope glycoproteins with lower CD4 dependence and reduced sensitivity to a fusion inhibitor.

    Virology 2006, 346:169-179. PubMed Abstract | Publisher Full Text OpenURL

  39. Leginagoikoa I, Daltabuit-Test M, Alvarez V, Arranz J, Juste RA, Amorena B, de Andres D, Lujan LL, Badiola JJ, Berriatua E: Horizontal Maedi-Visna virus (MVV) infection in adult dairy-sheep raised under varying MVV-infection pressures investigated by ELISA and PCR.

    Res Vet Sci 2006, 80:235-241. PubMed Abstract | Publisher Full Text OpenURL

  40. Narayan O: Role of macrophages in the immunopathogenesis of visna-maedi of sheep.

    Prog Brain Res 1983, 59:233-235. PubMed Abstract | Publisher Full Text OpenURL

  41. Shah C, Boni J, Huder JB, Vogt HR, Muhlherr J, Zanoni R, Miserez R, Lutz H, Schupbach J: Phylogenetic analysis and reclassification of caprine and ovine lentiviruses based on 104 new isolates: evidence for regular sheep-to-goat transmission and worldwide propagation through livestock trade.

    Virology 2004, 319:12-26. PubMed Abstract | Publisher Full Text OpenURL

  42. Ohagen A, Devitt A, Kunstman KJ, Gorry PR, Rose PP, Korber B, Taylor J, Levy R, Murphy RL, Wolinsky SM, Gabuzda D: Genetic and functional analysis of full-length human immunodeficiency virus type 1 env genes derived from brain and blood of patients with AIDS.

    J Virol 2003, 77:12336-12345. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  43. Liu P, Hudson LC, Tompkins MB, Vahlenkamp TW, Colby B, Rundle C, Meeker RB: Cerebrospinal fluid is an efficient route for establishing brain infection with feline immunodeficiency virus and transfering infectious virus to the periphery.

    J Neurovirol 2006, 12:294-306. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  44. Larue RS, Lengyel J, Jonsson SR, Andresdottir V, Harris RS: Lentiviral Vif degrades the APOBEC3Z3/APOBEC3H protein of its mammalian host and is capable of cross-species activity.

    J Virol 2010, 84:8193-8201. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  45. Nelson GW, Perelson AS: Modeling defective interfering virus therapy for AIDS: conditions for DIV survival.

    Math Biosci 1995, 125:127-153. PubMed Abstract | Publisher Full Text OpenURL

  46. Burkala EJ, He J, West JT, Wood C, Petito CK: Compartmentalization of HIV-1 in the central nervous system: role of the choroid plexus.

    AIDS 2005, 19:675-684. PubMed Abstract | Publisher Full Text OpenURL

  47. Glaria I, Reina R, Ramirez H, de Andres X, Crespo H, Jauregui P, Salazar E, Lujan L, Perez MM, Benavides J, et al.: Visna/Maedi virus genetic characterization and serological diagnosis of infection in sheep from a neurological outbreak.

    Vet Microbiol 2011.

    doi:10.1016/j.vetmic.2011.08.027

    OpenURL

  48. Fletcher NF, Meeker RB, Hudson LC, Callanan JJ: The neuropathogenesis of feline immunodeficiency virus infection: Barriers to overcome.

    Vet J 2010, 188:260-269. PubMed Abstract | Publisher Full Text OpenURL

  49. Reina R, Mora MI, Glaria I, Garcia I, Solano C, Lujan L, Badiola JJ, Contreras A, Berriatua E, Juste R, et al.: Molecular characterization and phylogenetic study of Maedi Visna and Caprine Arthritis Encephalitis viral sequences in sheep and goats from Spain.

    Virus Res 2006, 121:189-198. PubMed Abstract | Publisher Full Text OpenURL

  50. Liu SL, Rodrigo AG, Shankarappa R, Learn GH, Hsu L, Davidov O, Zhao LP, Mullins JI: HIV quasispecies and resampling.

    Science 1996, 273:415-416. PubMed Abstract | Publisher Full Text OpenURL

  51. Mordasini F, Vogt HR, Zahno ML, Maeschli A, Nenci C, Zanoni R, Peterhans E, Bertoni G: Analysis of the antibody response to an immunodominant epitope of the envelope glycoprotein of a lentivirus and its diagnostic potential.

    J Clin Microbiol 2006, 44:981-991. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  52. Hall T: BioEdit: a user-friendly biological sequence alignment editor ansd analysis progveam for Windows 95/98/NT.

    Nucl Acids Symp Ser 1999, 41:95-98. OpenURL

  53. Thompson JD, Gibson TJ, Plewniak F, Jeanmougin F, Higgins DG: The CLUSTAL_X windows interface: flexible strategies for multiple sequence alignment aided by quality analysis tools.

    Nucleic Acids Res 1997, 25:4876-4882. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  54. Swofford D: Phylogenetic Analysis Using Parsimony (*and Other Methods) 4.0. Sunderland: Sinauer Associates; 2003. OpenURL

  55. Pond SL, Frost SD: Not so different after all: a comparison of methods for detecting amino acid sites under selection.

    Mol Biol Evol 2005, 22:1208-1222. PubMed Abstract | Publisher Full Text OpenURL

  56. Pond SL, Frost SD: Datamonkey: rapid detection of selective pressure on individual sites of codon alignments.

    Bioinformatics 2005, 21:2531-2533. PubMed Abstract | Publisher Full Text OpenURL

  57. Posada D, Crandall KA: MODELTEST: testing the model of DNA substitution.

    Bioinformatics 1998, 14:817-818. PubMed Abstract | Publisher Full Text OpenURL

  58. Ronquist F, Huelsenbeck JP: MrBayes 3: bayesian phylogenetic inference under mixed models.

    Bioinformatics 2003, 19:1572-1574. PubMed Abstract | Publisher Full Text OpenURL

  59. Parker J, Rambaut A, Pybus OG: Correlating viral phenotypes with phylogeny: accounting for phylogenetic uncertainty.

    Infect Genet Evol 2008, 8:239-246. PubMed Abstract | Publisher Full Text OpenURL

  60. Drummond AJ, Rambaut A: BEAST: bayesian evolutionary analysis by sampling trees.

    BMC Evol Biol 2007, 7:214. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  61. Korber B, Myers G: Signature pattern analysis: a method for assessing viral sequence relatedness.

    AIDS Res Hum Retroviruses 1992, 8:1549-1560. PubMed Abstract | Publisher Full Text OpenURL