|
|
||||||||
Agence Française de Securite Sanitaire des Aliments (AFSSA), Laboratoire d'Etudes et de Recherches Caprines, Niort, France
Correspondence
Stephen Valas
s.valas{at}niort.afssa.fr
| ABSTRACT |
|---|
|
|
|---|
The GenBank/EMBL/DDBJ accession numbers for the sequences determined in this work are EU262438–EU262464 (gag) and EU262465–EU262481 (env).
| INTRODUCTION |
|---|
|
|
|---|
The development of an efficacious diagnosis of SRLV infection remains one of the biggest challenges in the control and eradication programmes that have been implemented in several European countries. Major obstacles to diagnosis development include delayed seroconversion and long-term-seronegative-infected animals, low viral loads impairing PCR tests and SRLV variability. As epidemiological studies have revealed the co-circulation of highly divergent genetic subtypes in restricted geographical regions (Germain & Valas, 2006
; Gjerset et al., 2007
; Grego et al., 2002
, 2007
; Leroux et al., 1997
; Reina et al., 2006
; Shah et al., 2004a
), up-to-date information on the prevalent SRLV subtypes circulating in each country is of crucial importance to improve the efficiency of the diagnosis tests. In addition, the spread of multiple subtypes in a specific geographical region, combined with the ability of SRLV strains to cross the species barrier, may give rise to dual infection and recombination. However, only one case of dual infection and viral recombination has been reported recently in naturally infected goats (Pisoni et al., 2007
), and whether interspecies transmission contributes to the evolution of genetic and biological properties of SRLV isolates is currently unknown.
Current tools for the characterization of SRLV diversity include sequencing of short genome segments, typically in the gag/pol and env regions, and, more rarely, full-genome sequencing. However, DNA sequencing methodologies are unsuitable for genetic analyses performed on large numbers of samples, limiting the detection of viral variants that are different from the local predominant subtypes. Large-scale, non-sequencing methods are thus required. Recently, the heteroduplex mobility assay (HMA) methodology has been adapted as a high-throughput genotyping tool for the study of SRLV genetic diversity (Germain & Valas, 2006
). This method was originally designed for the V1–V2 region of the env gene to discriminate among SRLV groups. In the present study, the HMA was validated as an accurate approach to estimate the genetic divergence among intra-group SRLV env sequences. Moreover, the HMA was adapted for the SRLV gag region spanning the nucleocapsid-coding sequence. French SRLV field isolates among sheep and goats from mixed flocks were investigated using HMA analysis performed in parallel on both the env and gag genes. The results of gag/env HMA subtyping were confirmed by genetic and phylogenetic analyses.
| METHODS |
|---|
|
|
|---|
Nested PCR.
A 394 bp fragment spanning the V1–V2 variable region of the env gene and a 238 bp fragment spanning the nucleocapsid-coding region of the gag gene were amplified by nested PCR. The gag and env fragments corresponded to nt 1605–1842 and nt 6352–6745 in the Cork genome, respectively (Saltarelli et al., 1990
). PCR amplification of the gag fragment was performed using the primer pairs P14/P15 in the first round and P40/P41 in the second round (Shah et al., 2004a
). The PCR primer pairs specific for the env gene were A5/A31 in the first round and A51/B31 in the second round (Germain & Valas, 2006
). First-round amplifications were performed using 500 ng genomic DNA in a 50 µl reaction containing 2 mM MgCl2, 0.2 mM each dNTP, 300 nM each primer and 2.5 U Taq DNA polymerase (iTaq; Bio-Rad). The PCR cycling conditions were 3 min at 95 °C, followed by 35 cycles of 40 s at 94 °C, 60 s at 57 °C and 60 s at 72 °C, with a final extension of 4 min at 72 °C. Second-round PCRs were performed with 2 µl of the first-round reaction used as template. The cycling conditions were the same for amplification of the gag fragment, but the hybridization step for amplification of the env fragment was carried out for 50 s at 60 °C. All PCRs were carried out in a Bio-Rad Icycler thermocycler.
HMA.
Heteroduplex formation reactions were carried out as described previously (Germain & Valas, 2006
). The relative mobilities of the heteroduplexes were calculated by dividing the migration distance of the heteroduplex bands (or the midpoint between two bands) by the migration distance of the corresponding homoduplex band. The same experimental conditions were carried out for gag and env HMA. A panel of 33 reference SRLV isolates was used in the env HMA analysis, comprising the prototype strains of subtypes A1 (Islandic strain K1514), B1 (American strain Cork) and B2 (Switzerland strain 5720), three French caprine lentiviruses (strains 680, 032 and 786) representative of subtype B1, and 27 French ovine lentiviruses belonging to either subtype B1 (six strains) or subtype B2 (21 strains), which have recently been characterized in our laboratory (Germain & Valas, 2006
).
DNA sequence analysis.
Specific PCR products were purified from agarose gels using a QIAquick Gel Extraction kit (Qiagen). Sequencing of both strands using a standard ABI BigDye terminator reaction (MilleGen) was performed on uncloned PCR fragments or on PCR products cloned into the linear pDrive vector using a cloning kit (Qiagen). Nucleotide sequences were aligned for subsequent phylogenetic analysis using CLUSTAL W version 1.8 (Thompson et al., 1994
). Manual rearrangements of the alignments, including gap exclusion and length adjustment, were applied to achieve optimal results. Pairwise genetic distances were calculated using MEGA version 3.1 with the Tamura–Nei substitution model. Phylogenetic construction was carried out using the neighbour-joining method (Saitou & Nei, 1987
) implemented in MEGA with the Tamura–Nei gamma distance (Tamura & Nei, 1993
). The statistical confidence of the topologies was assessed using 1000 bootstrap replicates (Felsenstein, 1985
).
Subtype reference strains.
The GenBank accession numbers of the previously published SRLV sequences are as follows: Cork (M33677
[GenBank]
), 680 (AJ400718
[GenBank]
), 032 (AJ400720
[GenBank]
), 786 (AJ400721
[GenBank]
), 4668 (AY445885
[GenBank]
), S93 (AF338226
[GenBank]
), 85/34 (U64439), WLC1 (AY362038
[GenBank]
), SA-OMVV (M31646
[GenBank]
), K1514 (M60609
[GenBank]
), EV-1 (S51392), 5560 (AY454175
[GenBank]
), 5561 (AY454176
[GenBank]
), 5720 (AY454218
[GenBank]
and AY842774
[GenBank]
), A-S76 (AY842749
[GenBank]
), A-S450 (AY842748
[GenBank]
), A-S2013 (AY842750
[GenBank]
), B-S6164 (AY842757
[GenBank]
), B-3056 (DQ149845
[GenBank]
), C-S192 (AY842753
[GenBank]
), C-S269 (AY842754
[GenBank]
), C-S3034 (AY842755
[GenBank]
), D-S536 (AY842738
[GenBank]
), D-S9103 (AY842740
[GenBank]
), D-S553 (AY842739
[GenBank]
), E-S46 (AY842763
[GenBank]
), E-S26 (AY842761
[GenBank]
) and E-S218 (AY842762
[GenBank]
), 5692 (AY454208
[GenBank]
) and CA1GA (AF322109
[GenBank]
).
| RESULTS |
|---|
|
|
|---|
Evaluation of the gag HMA
A 238 bp fragment spanning the coding region of the p14 nucleocapsid protein was selected for HMA analysis. The choice of this gag region was guided by (i) the availability of sequences representative of the different SRLV subtypes, (ii) the performance of the PCR primers allowing detection of all subtypes (Shah et al., 2004a
) and (iii) the variation in size of the PCR products due to group- and subtype-specific deletions. Indeed, alignment of reference sequences of SRLV subtypes revealed that variable domains with small deletions and insertions surrounded the highly conserved central domain of the p14 protein containing the zinc finger motifs (Fig. 1a
). Subtype reference sequences of CAEV-like viruses differed from those of VMV-like viruses by deletions in both terminal domains, except for subtype B1 lacking the deletion in the C-terminal domain. These genetic features between subtypes were confirmed by analysis of all available SRLV sequences retrieved from databases, showing that all CAEV-like viruses (35/35) contained a deletion in the N-terminal domain of the p14 protein, compared with only 1/20 VMV-like virus. Given that subtype B2 was represented by only one gag sequence, the p14-coding sequence of nine French ovine viruses that have previously been classified into subtype B2 on the basis of V1–V2 env sequences (Germain & Valas, 2006
) was determined. As shown in Fig. 1(b)
, eight French isolates contained a deletion in the N-terminal domain, whilst the remaining strain, E-S218, showed no deletion. However, phylogenetic analysis revealed that strain E-S218 belonged to group A on the basis of the gag sequence (see Fig. 4b
), indicating a possible recombinant origin between SRLV groups A and B. In addition, all viruses classified in subtype B2 (9/9) contained a deletion in the C-terminal domain of the p14 protein, whilst 88 % (29/33) of subtype B1 viruses showed no deletion within this domain.
|
|
|
|
A positive gag HMA PCR product was obtained for 73/75 (97.3 %) of the analysed blood samples. PCR products were tested individually by HMA using as the reference standard the Cork CAEV strain representative of subtype B1. As described above, the samples producing HMA patterns with one homoduplex band and two subsets of heteroduplexes with fast or moderate mobilities were assigned to subtype B1, whilst those forming HMA patterns with one homoduplex band and slow-migrating heteroduplexes were assigned to group A. The samples producing HMA patterns with two homoduplex bands were classified into subtype B2. The SRLV gag subtyping results are presented in Table 1
. Complete correlation was observed between the gag and env HMA subtyping for 67/73 samples (91.8 %). Five goat samples, four in flock B and one in flock D, gave discordant results by the gag/env HMA. These samples were classified into subtypes B1 and B2 by the env and gag HMA, respectively. Of 13 sheep samples analysed in flock E, 12 were classified into subtype B2 by the gag/env HMA, whereas the remaining sample was assigned to subtype B2 by the env HMA and group A by the gag HMA.
|
| DISCUSSION |
|---|
|
|
|---|
The env HMA has two main limitations for the genetic subtyping of SRLV. Firstly, the env gene exhibits the highest variability within the lentiviral genome, impeding the PCR amplification of some SRLV strains. Secondly, it cannot differentiate alone between pure subtypes and recombinant forms. In order to circumvent these problems, we developed a gag HMA based on a 238 bp fragment spanning the p14 nucleocapsid-coding region. The PCR primers used for the amplification of this gag fragment have been used previously for successful amplification of SRLV sequences representative of all subtypes (Shah et al., 2004a
). Moreover, sequence alignment revealed both group- and subtype-specific features in the p14-coding region that could be used as markers in the HMA. Indeed, all of the CAEV-like strains showed a short amino acid deletion in the N-terminal domain of the protein compared with all of the VMV-like strains, except for the prototype VMV K1514 strain. In addition, most of the reference sequences of subtype B1 differed from those of subtype B2, which contained an additional amino acid deletion in the C-terminal domain of the protein. This distinct feature between subtypes B1 and B2 was confirmed by sequence determination of the corresponding fragments of eight additional pure subtype B2 strains. These deletions, combined with base substitutions, allowed the gag HMA to assign subtypes B1 and B2 correctly.
The specificity of the gag/env HMA was verified subsequently using SRLV field isolates originating from seropositive sheep and goats from mixed flocks in France. The gag HMA subtype results correlated fully with sequence and phylogenetic analyses. Concordant results between the env HMA subtyping and phylogenetic analysis was observed for 30/31 (96.7 %) samples, one sample being classified into subtype B1 by the env HMA and into subtype B2 by sequencing. Overall, the HMA subtyping results revealed a predominance of subtypes B1 and B2 in the French mixed flocks investigated in this study. However, the relative distribution of these subtypes differs significantly in sheep and goats. Subtype B1 was found predominantly in goats, whereas subtype B2 was prevalent in sheep. Similar observations have been described for SRLV epidemiology in other European countries, including Italy and Switzerland (Grego et al., 2002
; Shah et al., 2004a
), suggesting that some SRLV subtypes may be better adapted to spread in goats and others in sheep. The HMA analysis also provided direct evidence for interspecies transmission of subtypes B1 and B2 in three out of four flocks for which PCR products were obtained from both sheep and goat samples, indicating that natural interspecies transmissions occur regularly. This was confirmed by DNA analysis showing that sheep and goat samples producing no heteroduplex when tested pairwise by HMA were phylogenetically closely related and displayed extremely low V1–V2 sequence divergences, ranging from 1.2 to 1.8 %. For comparison, these genetic distances are of the same extent as the viral quasispecies diversity within individual hosts (Herrmann et al., 2004
; Leroux et al., 1997
). To our knowledge, this is the first report of env sequences representative of subtype B2 in goats, confirming that no SRLV subtype is particularly prone to crossing the species barrier.
Comparison of gag and env HMA subtyping results revealed a different classification of viral sequences from the same sample in five goats from flock B and one sheep from flock E. The goat samples were classified as subtype B1 by the env HMA and subtype B2 by the gag HMA. The sheep sample was classified as subtype B2 by env HMA and group A by gag HMA. Phylogenetic analysis performed on both gag and env sequences from the sheep sample and three out of five goat samples confirmed the subtype assignment inferred by HMA. This is the first time, to our knowledge, that sequences from different genes and classified into different SRLV subtypes have been found within the same host in sheep and goats. Thus, these data indicate the relevance of subtyping two different gene fragments. The different phylogenetic classification of gag and env sequences obtained from the same sample strongly suggests the occurrence of inter-subtype recombination. Dual infection with different subtypes rather than infection with recombinant forms is unlikely, as both gag and env sequences were concomitantly amplified from all samples from these two flocks, excluding the possibility of PCR amplification failure due to genetic variability of viruses. In addition, all samples gave no heteroduplex in the gag/env HMA when tested with themselves, indicating the presence of only one major variant in each sample (data not shown). Overall, these results suggest that dual infection rarely occurs, although such events can lead to inter-subtype recombination in sheep and goats. This is consistent with a recent study showing that only two goats from a commercial flock of about 300 goats that were seropositive for SRLV reacted strongly to both group A- and group B-specific antigens, one of them being infected with a recombinant virus carrying a mosaic A/B env sequence (Pisoni et al., 2007
). The situation is quite different for other lentiviruses, particularly human and feline immunodeficiency viruses, for which dual infections and recombinant forms have been detected at an appreciable frequency (Bachmann et al., 1997
; Quinones-Mateu & Arts, 1999
; Robertson et al., 1995
; Troyer et al., 2004
). However, the frequency of dual infection in sheep and goats may be underestimated, as the viral load in peripheral blood is usually low and is subjected to temporal fluctuations (Klevjer-Anderson et al., 1984
), limiting the simultaneous detection of divergent SRLV strains by PCR from blood samples.
In addition to the high incidence of genetic diversity and evolution of SRLV, recombination between SRLV groups A and B may yield viruses with distinct biological properties. Indeed, prototypic VMV and CAEV strains differ in their cytopathic phenotype and have been classified as lytic and non-lytic viruses, respectively (Querat et al., 1984
). These SRLV strains also differ in their in vitro species tropism. In particular, the Icelandic VMV strain K1514 can use receptors for entry and/or cell-to-cell fusion in human cells (Bruett & Clements, 2001
; Hötzel & Cheevers, 2001
), whilst the lack of functional receptors is the only barrier that prevents CAEV replication in these cells (Mselli-Lakhal et al., 2000
), suggesting that recombinant virus between highly divergent SRLV subtypes may achieve the complete virus life cycle in cells other than those of small ruminants. Importantly, SRLV field isolates belonging to subtypes including the prototypic CAEV and VMV strains are still circulating in European countries (Barros et al., 2004
; Grego et al., 2002
; Leroux et al., 1997
; Rolland et al., 2002
; Shah et al., 2004a
). Up-to-date information on circulating SRLV isolates may thus be of crucial importance. In this regard, HMAs may offer an important contribution to improve the molecular biology-based surveillance of SRLV, and the accurate subtype prediction for both the gag and env region using this assay may be extremely helpful to estimate the frequency of inter-subtype genetic recombination.
| ACKNOWLEDGEMENTS |
|---|
| REFERENCES |
|---|
|
|
|---|
Barros, S. C., Ramos, F., Duarte, M., Fagulha, T., Cruz, B. & Fevereiro, M. (2004). Genomic characterization of a slow/low maedi visna virus. Virus Genes 29, 199–210.[CrossRef][Medline]
Bruett, L. & Clements, J. E. (2001). Functional murine leukaemia virus vectors pseudotyped with the visna virus envelope show expanded visna virus cell tropism. J Virol 75, 11464–11473.
Felsenstein, J. (1985). Confidence limits on phylogenies: an approach using the bootstrap. Evolution 39, 783–791.[CrossRef]
Germain, K. & Valas, S. (2006). Distribution and heterogeneity of small ruminant lentivirus envelope subtypes in naturally infected French sheep. Virus Res 120, 156–162.[CrossRef][Medline]
Gjerset, B., Storset, A. K. & Rimstad, E. (2006). Genetic diversity of small-ruminant lentiviruses: characterization of Norwegian isolates of caprine arthritis encephalitis virus. J Gen Virol 87, 573–580.
Gjerset, B., Jonassen, C. M. & Rimstad, E. (2007). Natural transmission and comparative analysis of small ruminant lentiviruses in the Norwegian sheep and goat populations. Virus Res 125, 153–161.[CrossRef][Medline]
Grego, E., Profiti, M., Giammarioli, M., Giannino, L., Rutili, D., Woodall, C. & Rosati, S. (2002). Genetic heterogeneity of small ruminant lentiviruses involves immunodominant epitope of capsid antigen and affects sensitivity of single-strain-based immunoassay. Clin Diagn Lab Immunol 9, 828–832.[CrossRef][Medline]
Grego, E., Bertolotti, L., Quasso, A., Profiti, M., Lacerenza, D., Muz, D. & Rosati, S. (2007). Genetic characterization of small ruminant lentiviruses in Italian mixed flocks: evidence for a novel genotype circulating in a local goat population. J Gen Virol 88, 3423–3427.
Herrmann, L. M., Hötzel, I., Cheevers, W. P., On Top, K. P., Lewis, G. S. & Knowles, D. P. (2004). Seven new progressive pneumonia virus (OPPV) field isolates from Dubois Idaho sheep comprise part of OPPV clade II based on surface envelope glycoprotein (SU) sequences. Virus Res 102, 215–220.[CrossRef][Medline]
Hötzel, I. & Cheevers, W. P. (2001). Host range of small ruminant lentivirus cytopathic variants determined with a selectable caprine arthritis-encephalitis virus pseudotype system. J Virol 75, 7384–7391.
Klevjer-Anderson, P., Adams, D. S., Anderson, L. W., Banks, K. L. & McGuire, T. C. (1984). A sequential study of virus expression in retrovirus-induced arthritis of goats. J Gen Virol 65, 1519–1525.
Leroux, C., Chastang, J., Greenland, T. & Mornex, J. F. (1997). Genomic heterogeneity of small ruminant lentiviruses: existence of heterogeneous populations in sheep and of the same lentiviral genotypes in sheep and goats. Arch Virol 142, 1125–1137.[CrossRef][Medline]
Mselli-Lakhal, L., Favier, C., Leung, K., Guiguen, F., Grezel, D., Miossec, P., Mornex, J. F., Narayan, O., Querat, G. & Chebloune, Y. (2000). Lack of functional receptors is the only barrier that prevents caprine arthritis-encephalitis virus from infecting human cells. J Virol 74, 8343–8348.
Pisoni, G., Quasso, A. & Moroni, P. (2005). Phylogenetic analysis of small-ruminant lentivirus subtype B1 in mixed flocks: evidence for natural transmission from goats to sheep. Virology 339, 147–152.[CrossRef][Medline]
Pisoni, G., Bertoni, G., Puricelli, M., Maccalli, M. & Moroni, P. (2007). Demonstration of coinfection with and recombination by caprine arthritis-encephalitis virus and maedi-visna virus in naturally infected goats. J Virol 81, 4948–4955.
Querat, G., Barban, V., Sauze, N., Filippi, P., Vigne, R., Russo, P. & Vitu, C. (1984). Highly lytic and persistent lentiviruses naturally present in sheep with progressive pneumonia are genetically distinct. J Virol 52, 672–679.
Quinones-Mateu, M. & Arts, E. J. (1999). Recombination in HIV-1: update and implications. AIDS Rev 1, 89–100.
Reina, R., Mora, M. I., Glaria, I., García, I., Solano, C., Luján, L., Badiola, J. J., Contreras, A., Berriatua, E. & other authors (2006). Molecular characterization and phylogenetic study of maedi visna and caprine arthritis encephalitis viral sequences in sheep and goats from Spain. Virus Res 121, 189–198.[CrossRef][Medline]
Robertson, D. L., Sharp, M., McCutchan, F. E. & Hahn, B. H. (1995). Recombination in HIV-1. Nature 374, 124–126.[Medline]
Rolland, M., Mooney, J., Valas, S., Perrin, G. & Mamoun, R. Z. (2002). Characterization of an Irish caprine lentivirus strain: SRLV phylogeny revisited. Virus Res 85, 29–39.[CrossRef][Medline]
Saitou, N. & Nei, M. (1987). The neighbor-joining method: a new method for reconstructing phylogenetic trees. Mol Biol Evol 4, 406–425.[Abstract]
Saltarelli, M., Querat, G., Konings, D. A., Vigne, R. & Clements, J. E. (1990). Nucleotide sequence and transcriptional analysis of molecular clones of CAEV which generate infectious virus. Virology 179, 347–364.[CrossRef][Medline]
Shah, C., Böni, J., Huder, J. B., Vogt, H.-R., Mühlherr, J., Zanoni, R., Miserez, R., Lutz, H. & Schüpbach, J. (2004a). 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 319, 12–26.[CrossRef][Medline]
Shah, C., Huder, J. B., Böni, J., Schönmann, M., Mühlherr, J., Lutz, H. & Schüpbach, J. (2004b). Direct evidence for natural transmission of small-ruminant lentiviruses of subtype A4 from goats to sheep and vice versa. J Virol 78, 7518–7522.
Tamura, K. & Nei, M. (1993). Estimation of the number of nucleotide substitutions in the control region of mitochondrial DNA in humans and chimpanzees. Mol Biol Evol 10, 512–526.[Abstract]
Thompson, J. D., Higgins, D. G. & Gibson, T. J. (1994). CLUSTAL W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. Nucleic Acids Res 22, 4673–4680.
Troyer, J. L., Pecon-Slattery, J., Roelke, M. E., Black, L., Packer, C. & O'Brien, S. J. (2004). Patterns of feline immunodeficiency virus multiple infection and genome divergence in a free-ranging population of African lions. J Virol 78, 3777–3791.
Zanoni, R. G. (1998). Phylogenetic analysis of small ruminant lentiviruses. J Gen Virol 79, 1951–1961.[Abstract]
Received 29 January 2008;
accepted 11 April 2008.
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| INT J SYST EVOL MICROBIOL | MICROBIOLOGY | J GEN VIROL |
| J MED MICROBIOL | ALL SGM JOURNALS | |