Next Article in Journal
Small RNA NGS Revealed the Presence of Cherry Virus A and Little Cherry Virus 1 on Apricots in Hungary
Next Article in Special Issue
Enrichment Preferences of FIV-Infected and Uninfected Laboratory-Housed Cats
Previous Article in Journal
Antiviral Effects of Clinically-Relevant Interferon-α and Ribavirin Regimens against Dengue Virus in the Hollow Fiber Infection Model (HFIM)
Previous Article in Special Issue
Non-Primate Lentiviral Vectors and Their Applications in Gene Therapy for Ocular Disorders
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Viruses
Volume 10
Issue 6
10.3390/v10060315
viruses-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

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

Genetic Characterization and Phylogenetic Analysis of Small Ruminant Lentiviruses Detected in Spanish Assaf Sheep with Different Mammary Lesions

by
Elena Gayo
1,
Vincenzo Cuteri
2,
Laura Polledo
3,
Giacomo Rossi
2,
Juan F. García Marín
1 and
Silvia Preziuso
2,*
1
Pathological Anatomy Section, Animal Health Department, School of Veterinary Medicine, University of Leon, via Profesor Pedro Carmenes s/n Campus de Vegazana, 24071 León, Spain
2
School of Biosciences and Veterinary Medicine, University of Camerino, Via Circonvallazione 93/95, 62024 Matelica (MC), Italy
3
Micros Veterinaria, INDEGSAL, via Profesor Pedro Carmenes s/n Campus de Vegazana, 24071 León, Spain
*
Author to whom correspondence should be addressed.
Viruses 2018, 10(6), 315; https://doi.org/10.3390/v10060315
Submission received: 7 April 2018 / Revised: 30 May 2018 / Accepted: 7 June 2018 / Published: 9 June 2018
(This article belongs to the Special Issue Nonprimate Lentivirus)
Browse Figures
Versions Notes

Abstract

:
Small Ruminant Lentiviruses (SRLVs) are widespread in many countries and cause economically relevant, slow, and persistent diseases in sheep and goats. Monitoring the genetic diversity of SRLVs is useful to improve the diagnostic tools used in the eradication programs. In this study, SRLVs detected in Spanish Assaf sheep with different grades of lymphoproliferative mastitis were sequenced. Genetic characterization showed that most samples belonged to type A and were closer to Spanish SRLV isolates previously classified as A2/A3. Four samples belonged to subtype B2 and showed higher homology with Italian B2 strains than with Spanish B2 isolates. Amino acid sequences of immuno-dominant epitopes in the gag region were very conserved while more alterations were found in the LTR sequences. No significant correlations were found between grades of mastitis and alterations in the sequences although samples with similar histological features were phylogenetically closer to each other. Broader genetic characterization surveys in samples with different grades of SRLV-lesions are required for evaluating potential correlations between SRLV sequences and the severity of diseases.

1. Introduction

Small Ruminant Lentiviruses (SRLVs) include Visna-maedi virus (VMV) and Caprine arthritis encephalitis virus (CAEV), which cause slow inflammatory diseases in sheep and goats named, respectively, Visna-Maedi (VM) and Caprine Arthritis Encephalitis (CAE). Viruses are transmitted mainly via the respiratory route and by colostrum intake [1,2] and cause persistent infections with a long incubation period. According to the current nomenclature based on gag–pol and pol sequences, SRLV can be subdivided into genotypes A–E with subtypes present in A, B, and E [3,4,5]. Typically, VMV was believed to infect specifically sheep and was included in genotype A while CAEV was considered goat-specific and was included in genotype B. However, several investigations showed that cross-infection may occur ([6] and other studies reviewed in [5]). Infected animals can develop neurological, pulmonary, arthritic, and/or mammary diseases that affect considerable animal welfare and production. Different patterns of inflammatory mononuclear cell accumulation are observed usually in the central nervous system, lung, joint, and/or udder and the predominant clinical manifestation depends on the severity and extension of the lesions reached in the affected organs [7]. Mammary lesions often consist of mononuclear cell infiltration with scattered hyperplastic lymphoid follicles [7,8]. Moderate lesions are described when diffused infiltration of lymphocytes within lobules with distortion of acini are observed and mild lesions are reported when occasional aggregates of lymphocytes in inter-acinar stroma are present [9]. Recently, minimal lesions consisting of a few small clusters to multifocal small groups of inflammatory cells with the presence of a few small lymphoid aggregates and/or one small lymphoid follicle have been reported in SRLV infected sheep [10].
One of the most important productive impacts of SRLV disease is due to the premature removal of diseased animals because of low milk production and quality with consequent economic losses in the milk-related and lamb/kids-related industry [11,12,13]. Due to the significant economic impact of diseases, VM and CAE are included in the OIE List and specific control and eradication programs are carried out in many countries [14,15]. Availability of sensitive and specific diagnostic tests is of great importance for a correct discrimination between infected and non-infected animals. Due to the high rate of genetic diversity, new genotypes and subtypes might escape the diagnostic detection with the possible consequence of invalidating any eradication program in place [1]. Therefore, genotype and subtype surveys of the circulating SRLVs should be encouraged. Many studies describe phylogenetic analysis of SRLVs found during epidemiological surveys or in outbreaks of diseases, but only a few of them describe the histopathological lesions observed in target organs [16,17,18]. For example, mild mammary lesions with a multifocal method to diffuse mononuclear inflammatory interstitial infiltrates have been observed in sheep with arthritis and infected by B2 SRLV [17]. Different histopathological scores have been reported in mammary glands but not in lungs, synovial membranes of joints, or the chorioid plexus of five goats infected by A4 SRLV [16]. A 13–14 nucleotide deletion in the R region of the LTR has been observed in sheep with a decreased pathology in the lung but not in the udder even though SRLV subgenotypes were not known [18]. To our knowledge, correlations among histopathological grading of mammary lesions and SRLV genotypes and subtypes are yet to be investigated.
The aim of this work was to carry out genetic characterization and phylogenetic analysis of SRLV detected in Spanish sheep showing different histopathological grades of mastitis.

2. Material and Methods

2.1. Samples

A total of 35 udder samples were collected randomly at the slaughterhouse in the region of Castilla y Leon, Northwestern Spain, from Assaf sheep (1–4 years of age) belonging to seropositive flocks between March 2017 and May 2017. Nineteen samples were collected at the slaughterhouse named M (samples M1–M19) and 16 were collected at the slaughterhouse named Q (samples Q1–Q16), which was about 45 Km far from the slaughterhouse M. A first aliquot of each sample was stored at −20 °C and DNA was obtained from 25 mg of each sample by using the Genomic DNA isolation Kit (Norgen Biotek Corp., Thorold, ON, Canada) and following the manufacturer’s instructions when eluting the DNA in 100 µL final volume. A second aliquot of samples was fixed in 10% neutral buffered formalin for 48 h at room temperature and embedded in paraffin wax (FFPE) for histopathology and immunohistochemistry (IHC). In addition, 2 FFPE mammary samples (N16-426 and N-17-44) of sheep (Assaf breed, 3 years old) with histological mastitis referable to SRLV disease were available for this study. DNA from these latter samples was obtained from 4 slides 10 µm thick of each sample by Recover All Total Nucleic Acid Isolation (Ambion, Waltham, MA, USA) following the instructions and eluting DNA with 60 µL of elution solution warmed up at 95 °C.

2.2. Histopathology and Immunohistochemistry

Slides 4 µm thick were obtained from FFPE samples and were stained with haematoxylin and eosin (HE) for histopathology. Grading of histopathological lesions of mammary glands was carried out by three independent pathologists, which was previously described [10]. Briefly, “no lesions” was defined when no inflammatory cells were observed and “minimal lesions” (+) consisted of a few small clusters.
Multifocal small groups of inflammatory cells with the presence of a few small lymphoid aggregates and/or one small lymphoid follicle known as ‘moderate lesions’ (++) were characterized by the multifocal method to diffuse interstitial non-supportive inflammation and/or the presence of two to 15 lymphoid aggregates/follicles and ‘severe lesions’ (+++) consisted of a marked diffuse interstitial mastitis and/or the presence of >15 lymphoid aggregates/follicles. When differences between the severity of interstitial inflammatory infiltrates and the presence of lymphoid follicles were observed in the same organ of an animal, the most severe lesion was considered as the score for the lesion in that target organ [10].
Serial sections (4 µm) were used for IHC, which was previously reported [19]. A monoclonal antibody to the SRLV core protein p28 (VMRD Inc., Pullman, WA, USA) diluted 1:1000 was used. A technique based on an avidin-biotin-peroxidase complex (VECTASTAIN® ELITE® ABC Kits, Vector Laboratories, Burlingame, CA, USA) with diaminobenzidine as the chromogen (DAB Peroxidase substrate kit—Vector Laboratories, Burlingame, CA, USA) was used to stain the antigen.

2.3. PCR

A nested PCR was used to amplify about 800 bp long sequences of SRLV gag-pol genes, which was reported previously [4]. Primers GAG-F1 and POL-R1 were used in the first PCR. The product obtained was used as a template in a second PCR with primers GAG-F2 and POL-R2. The PCR mixture included 50 µL 2× Taq PCR Master Mix, 500 nM each primer, 4 µL DNA, and PCR grade water up to 100 µL final volume. PCR conditions were 94 °C for 5 min, 45 cycles of 94 °C for 1 min, 55 °C for 1 min, 72 °C for 2 min, and a final extension of 72 °C for 10 min. The second PCR was carried out with the same conditions but 5 µL of the first PCR products were used as the template and the annealing temperature was 60 °C instead of 55 °C [4]. LTR sequences (203 bp long) were amplified by nested PCR with primers described elsewhere [20]. The PCR reaction mix was described above, but 2 µL of DNA (first PCR) or 2 µL of the first PCR products (second PCR) were used as the template. PCR conditions were 94 °C for 5 min, 35 cycles of 94 °C for 30 s, 55 °C (or 50 °C in the second PCR) for 30 s, 72 °C for 40 s, and a final extension of 72 °C for 7 min. PCR products were visualized in 1.5% agarose gel and positive samples were submitted to an external laboratory for sequencing (BMR Genomics, Padova, Italy). Both the sense and antisense strands were sequenced by performing two independent reactions for each PCR product. Nucleotide gag-pol sequences were deposited in GenBank (Accession numbers MH179145—MH179153 and MH179156—MH179159).

2.4. Sequence Analysis

Nucleotide sequences were manually checked and edited with the program BioEdit. A preliminary analysis by BLASTn was carried out to detect regions of similarity with sequences included in databases. Sequences of strains considered to be prototypes of different genotypes and SRLV sequences highly similar to those found in the samples were included in the study (Figure 1 and Figure 3). Sequences were aligned by MUSCLE [21] and phylogenetic trees were inferred with the program MEGA 7.0.21 [22]. The best-fitting nucleotide substitution models were estimated and the General Time Reversible model [23] with a gamma distribution with invariant sites (gag sequences) or a Kimura 2-parameter (LTR sequences) model [24] with gamma-distributed rates among sites were used with bootstrap values based on 1000 repetitions. Phylogeny was estimated by both the neighbor-joining algorithm (NJ) and the maximum likelihood (ML) method. Correlations among sequence alterations and histological features were evaluated by using the Fisher’s exact test. Pairwise distances between sequences of samples and sequences of reference strains belonging to different genotypes were calculated by MEGA 7.0.21 with the p-distance model [22].

3. Results

Histological examination and grading of mammary lesions (see Figure 1) resulted in five samples with severe lesions, 13 samples with moderate lesions, 11 samples with minimal lesions, and five samples without lesions (see Table 1). Grading of M5 sample was not possible due to a concomitant purulent mastitis.
IHC results were used to distinguish SRLV infected from uninfected sheep. Three out of the five samples without lesions (M13, M18 and Q2) were negative by IHC and by both gag-pol and LTR PCR. Therefore, they were considered negative (see Table 1). The remaining 32 samples of groups M and Q were positive by IHC (see Figure 2). Sixteen out of the 32 M and Q samples and the 2 N samples were positive by gag-pol PCR, but good-quality sequences were obtained only from 15 samples. LTR PCR products were obtained in all but two IHC-positive samples.

3.1. Analysis of Gag Sequences

Genotyping was carried out by phylogenetic analysis of partial gag gene sequences, according to the taxonomic classification proposed by Shah et al. [3]. All sequences were different from each other and nine samples were type A and four samples were type B. In particular, samples M5, M12, M15, M19, Q7, Q8, Q10, N16-426M, and N17-44M clustered within genotype A were more closely related with strains 292, 160, 166, and 697, which were previously detected in the same Spanish region (see Figure 3). Only strain 697 had been fully sequenced. Since similar values located with this isolate intermingled between A2 and A3, the isolate 697 had been assigned to the A2/A3 subtype [25]. Samples M1, M3, M17, and Q1 resulted of genotype B and subtype B2 (see Figure 3). Additionally, phylogenic and BLAST analysis showed that they were more related to B2 viruses detected in Italy than in Spain (Ov496). These results were confirmed by the pairwise distances comparison (see Table 2). Samples M12, M15, M19, Q7, Q8, Q10, N16-426, and N17-44 were more closely related to the A2/A3 Spanish strain HQ848062.1 (0.105–0.142). Moreover, samples M1, M3, M17, and Q1 were more closely related to the B2 strains FJ195346.1 and EU010126.1. In particular, they were more closely related to the Italian strain EU010126.1 (0.064–0.081) than to the Spanish strain FJ195346.1 (0.092–0.102).
Nucleotide sequences were translated into amino acid sequences and the results of the alignment and comparison with the most representative sequences are reported in Figure 4a,b. The set of primers used in this study amplifies a partial sequence of the gag gene codifying for the majority of the capsid protein (CA). Comparing amino acid sequence alterations of the partial gag protein obtained, the “GG” motifs of the four type B sequences were glycine-glycine (GG) like type B reference SRLVs while those of the 11 type A samples were asparagine-valine (NV) like other type A reference SRLVs (see Figure 4a). In type A samples, sequences of epitopes 2 and 3 of reference isolates and of most samples were conserved (see Figure 4a). Only arginine (R) replaced lysine (K) in samples N16-426 and N17-44 and serine (S) replaced asparagine (N) in sample Q7. Type B samples had highly conserved epitope 3 sequences since there are only two alterations in M1 (isoleucine (I) instead of serine (S) and lysine (K) instead of glutamic acid (E). More alterations were found in epitope 2 where three out of four type B samples showed one alteration in comparison with type B reference isolates (see Figure 4a,b).
In the Major Homology Region (MHR), which is usually a highly conserved sequence in the gag gene of all retroviruses, some alterations were present in type A samples. In particular, all but one type A samples showed one or two alterations comparing to the A2/A3 reference strain 697 (see Figure 4b). Samples M5, M12, M15, Q10, and N16-426 had isoleucine (I) instead of valine (V) at the fourth position as type A1 reference strain SA-OMVV. This latter had also a serine (S) instead of asparagine (N) at the ninth position. Sample Q7 had not only this alteration, but also glutamic acid (E) instead of aspartic acid (D) at the 14th position, which shows the same alterations found in the A2/A3 Spanish strain 160. Sample N17-44 had isoleucine (I) instead of serine (S) at the 11th position unlike the other samples and reference strains. Unusual alterations were found in sample M19 where asparagine (N) replaced serine (S) at the 11th position and lysine (K) substituted glutamic acid (E) at the 21st position. Type B samples had highly conserved MHR sequences and showed the same amino-acidic sequences of B2 reference isolates even though some single alterations were present in the nucleotide sequences (see Figure S1). A significant correlation among sequence alterations and severity of mastitis was not found (p > 0.05).

3.2. Analysis of LTR Sequences

The alignment and phylogenetic analysis of LTR nucleotide sequences showed that most samples were closer to the reference Spanish A2/A3 strain 697 (Figure 5 and Figure 6a,b). Comparing samples with the 697 reference strain, Q10 and Q11 showed a 23 nt deletion (9133–9160 nt) in the R region, which appeared longer than the 13 nt deletion present in the reference strains 697, EV1, and in other Spanish strains (C3, 160, 292). Sample Q12 had similar deletions than 697 while other samples showed 2–8 nt deletions in the same tract. No significant differences were found among sequences of samples with a different grade of mastitis (p > 0.05). Similarly, as previously described, the TATAbox, the polyadenylation signal, and the AML (vis), which is a site possessing the consensus sequence for the AML/PEBP2/CBF transcriptional factor family [26], were conserved and only a substitution G with A was present in three samples at position 9065 of AML (vis).

4. Discussion

This study describes for the first time gag and LTR sequences of SRLVs detected in Spanish Assaf sheep with different grades of histopathological mastitis and their phylogenetic relationships in the context of known SRLV sequences.
Although initially genotypes B viruses were thought to infect only goats, it is not unusual to find reports about infections by type B viruses in sheep and by type A viruses in goats [27,28,29,30]. Even in this study, both genotypes A and B SRLVs have been found in ovine samples. It was not known if the examined sheep had contact with goats, but infection with A genotypes in goats was not reared in contact with sheep, which was reported [31].
The gag-pol phylogenetic tree and the pairwise genetic distances comparison revealed that most sequences of the samples were closer to the Spanish A2/A3 isolate 697 while four sequences belonged to subtype B2. Isolate 697 has been previously detected in sheep with neurological diseases from the Spanish region of Castilla y León in Spain and has been classified as A2/A3 because differences between A2 and A3 are often not large enough to separate the two groups [3,25]. Partial sequences of viruses classified as A2/A3 have been detected further seven sheep with neurological signs in Spain [25]. These findings suggest that, in Northern Spain, subtype A2/A3 SRLVs is genetically related to SRLVs, which caused nervous diseases. However, in our cases, neurological signs were not reported. Samples M1, M3, M17, and Q1 were closer to Italian than to Spanish B2 isolates and mammary lesions were found from moderate to severe in three out of four samples while, in sample Q1, lesions were not found. B2 SRLV has been detected in Spain for the first time in SRLV-seropositive adult sheep of the Rasa Aragonesa breed, which shows clinical signs of arthritis [17,32]. Mammary histological lesions were present in 10 out of 13 animals with arthritis, which suggests that udders can be involved even if clinical signs might remain unrecognized until the losses of milk production are significant. B2 viruses have been detected in Italian small ruminants during epidemiological surveys, but data about clinical signs or histological lesions are not reported [4,33].
In addition, good-quality sequences about 800 bp long were obtained from FFPE samples. Fixation in formalin and embedding in paraffin at high temperatures is thought to degrade DNA. Fragmentation of DNA molecules can interfere with their amplification by PCR and with consequent sequencing. In our case, good-quality DNA has been extracted and amplified by PCR from archival FFPE samples, which suggests that this method could be attempted for studying FFPE samples as well as for retrospective investigations.
Analysis of the genetic sequences is important not only for evaluating the spread of SRLV types and subtypes but also for monitoring antigenic variability. Actually, remarkable antigenic variation might be responsible for the misdiagnosis of highly divergent genotypes [34]. The gag gene encodes nucleocapsid, capsid, and matrix proteins. Indirect diagnostic assays usually use the capsid protein as the antigen, which helps monitor immuno-dominant epitopes of gag-encoded structural proteins. This is useful for detecting antigenic variability in the field and forevaluating and improving the sensitivity of indirect diagnostic tools. Alterations in the amino acid sequences of immuno-dominant epitope regions suggest altered antigenicity, which may affect the sensitivity of serological tests such as ELISA and AGID. The gag-pol set of primers used in this study allowed sequencing only of epitopes 2 and 3. Amino acid sequences of epitopes 2 and 3 of type A2/A3 viruses were quite conserved and limited alterations only in three and one samples, respectively. Epitope 2 of B2 isolates had more alterations, which shows single amino acid alterations in three out of four sequences. In addition, more variability was found in the MHR of A2/A3 viruses, which show all but two samples and at least one alteration in comparison with the reference A2/A3 strain 697. The MHR is usually conserved in many retroviruses and is essential for viral assembly [35]. Mutations in the MHR sequence of HIV-1 cause capsid assembly that reduces infectivity [36]. While some studies have been carried out on MHR of human retroviruses, the consequences of MHR mutations on infectivity of SRLVs should be better investigated. MHR of B2 viruses and GG motif of both A2/A3 and B2 viruses in the gag amino acid sequences, AML (vis) motif, TATA-box, and poly-A of both A2/A3 and B2 viruses in the LTR nucleotide sequences were highly conserved, which was previously reported in strains from different geographic areas [31,37,38].
Most LTR sequences showed higher homology with A2/A3 Spanish SRLV isolates. Samples Q10 and Q11 showed a 23 nt deletion in the R region, which appeared longer than the 13 nt deletion observed in type A2/3 reference isolate 697. A 13 nt deletion in this region has been found in sequences of clinically affected sheep and a correlation among this deletion and the appearance of clinical signs has been suggested [32]. On the contrary, a similar deletion has been found in SRLVs infecting asymptomatic sheep and the lungs of animals infected with viruses carrying the deletion were significantly less affected than sheep infected with viruses without deletion [18]. In the present study, significant correlations among deletions in the R region of the LTR and severity of mammary lesions were not found (p > 0.05). Samples with deletions were from sheep with more severe mammary lesions and sheep with moderate to severe mastitis did not show this deletion (see Figure 6b).
Although the severity of mammary lesions was not significantly related to the viral genotype, SRLV sequences from samples with similar grades of lesions (no/minimal and moderate/severe) were most closely related to each other (see Figure 3 and Figure 5 and Table 2). Considering the high economic impact of SRLV diseases, some countries aim to eradicate the diseases by identifying and prematurely culling infected animals. Selecting animals on the basis of serological results could determine the selection of SRLV variants with significant alterations in the antigen sequences. Permanent and extensive surveys should be encouraged in different countries to evaluate the antigenic variability of SRLV and to monitor the sensitivity and specificity of diagnostic tests in detecting these variants. In particular, seronegative animals should be investigated for infections by new viral genotypes not detected by traditional serological tests. Histological screening of different target organs at the slaughterhouse could be a useful tool for selecting samples with lesions, which suggests an SRLV disease in seronegative animals.

5. Conclusions

In conclusion, this is the first study investigating the association between the SRLV sequence analysis and histopathological grading of mammary lesions in sheep. Circulation of SRLVs of types A2/A3 and B2 in Spanish Assaf sheep was confirmed and new viral variants have not been found, but moderate alterations were present in some immuno-dominant epitopes and in the MHR tract.
No significant correlation was found among histological features and alterations in the sequences. Although some sequences obtained from samples with similar grades of mammary lesions appeared closer to each other, more extensive and interdisciplinary studies are required for establishing the existence of viral clusters with a higher or lower pathogenicity for specific target organs.

Supplementary Materials

The following are available online at https://www.mdpi.com/1999-4915/10/6/315/s1, Figure S1: Nucleotide alignment of MHR sequences in the gag gene of the samples and of some reference strains.

Author Contributions

S.P. and E.G. conceived and designed the experiments. E.G. performed many experiments. S.P., G.R., L.P., J.F.G.M., and V.C. performed some experiments. E.G. and S.P. analyzed the data and wrote the paper. All authors read and approved the paper.

Acknowledgments

This work has been partially supported by the Spanish government (LE361A12-1 project and FPU13/01081 grant) and by the University of Camerino (FAR-Preziuso).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Peterhans, E.; Greenland, T.; Badiola, J.; Harkiss, G.; Bertoni, G.; Amorena, B.; Eliaszewicz, M.; Juste, R.A.; Krassnig, R.; Lafont, J.P.; et al. Routes of transmission and consequences of small ruminant lentiviruses (SRLVs) infection and eradication schemes. Vet. Res. 2004, 35, 257–274. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Preziuso, S.; Renzoni, G.; Allen, T.E.; Taccini, E.; Rossi, G.; DeMartini, J.C.; Braca, G. Colostral transmission of maedi visna virus: Sites of viral entry in lambs born from experimentally infected ewes. Vet. Microbiol. 2004, 104, 157–164. [Google Scholar] [CrossRef] [PubMed]
  3. Shah, C.; Böni, J.; Huder, J.B.; Vogt, H.R.; Mühlherr, J.; Zanoni, R.; Miserez, R.; Lutz, H.; Schüpbach, 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. [Google Scholar] [CrossRef] [PubMed]
  4. Grego, E.; Bertolotti, L.; Quasso, A.; Profiti, M.; Lacerenza, D.; Muz, D.; Rosati, S. Genetic characterization of small ruminant lentivirus in Italian mixed flocks: Evidence for a novel genotype circulating in a local goat population. J. Gen. Virol. 2007, 88, 3423–3427. [Google Scholar] [CrossRef] [PubMed]
  5. Ramírez, H.; Reina, R.; Amorena, B.; de Andrés, D.; Martínez, H.A. Small ruminant lentiviruses: Genetic variability; tropism and diagnosis. Viruses 2013, 5, 1175–1207. [Google Scholar] [CrossRef] [PubMed]
  6. Leroux, C.; Chastang, J.; Greenland, T.; Mornex, J.F. Genomic heterogeneity of small ruminant lentiviruses: Existence of heterogeneous populations in sheep and of the same lentiviral genotypes in sheep and goats. Arch. Virol. 1997, 142, 1125–1137. [Google Scholar] [CrossRef] [PubMed]
  7. Minguijón, E.; Reina, R.; Pérez, M.; Polledo, L.; Villoria, M.; Ramírez, H.; Leginagoikoa, I.; Badiola, J.J.; García-Marín, J.F.; de Andrés, D.; et al. Small ruminant lentivirus infections and diseases. Vet. Microbiol. 2015, 181, 75–89. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  8. Lujan, L.; García Marin, J.F.; Fernandez de Luco, D.; Vargas, A.; Badiola, J.J. Pathological changes in the lungs and mammary glands of sheep and their relationship with maedi-visna infection. Vet. Rec. 1991, 129, 51–54. [Google Scholar] [CrossRef] [PubMed]
  9. Benavides, J.; Fuertes, M.; García-Pariente, C.; Otaola, J.; Delgado, L.; Giraldez, J.; García Marín, J.F.; Carmen Ferreras, M.; Pérez, V. Impact of maedi-visna in intensively managed dairy sheep. Vet. J. 2013, 197, 607–612. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  10. Gayo, E.; Polledo, L.; Balseiro, A.; Martínez, C.P.; García Iglesias, M.J.; Preziuso, S.; Rossi, G.; García Marín, J.F. Inflammatory Lesion Patterns in Target Organs of Visna/Maedi in Sheep and their Significance in the Pathogenesis and Diagnosis of the Infection. J. Comp. Pathol. 2018, 159, 49–56. [Google Scholar] [CrossRef] [PubMed]
  11. Anderson, B.C.; Bulgin, M.S.; Adams, S.; Duelke, B. Firm udder in periparturient ewes with lymphocytic accumulations; retrovirus infection; and milk unavailable at the teat. J. Am. Vet. Med. Assoc. 1985, 186, 391–393. [Google Scholar] [PubMed]
  12. Ryan, D.P.; Greenwood, P.L.; Nicholls, P.J. Effect of caprine arthritis-encephalitis virus infection on milk cell count and N-acetyl-beta-glucosaminidase activity in dairy goats. J. Dairy Res. 1993, 60, 299–306. [Google Scholar] [CrossRef] [PubMed]
  13. Martinez-Navalon, B.; Peris, C.; Gomez, E.A.; Peris, B.; Roche, M.L.; Caballero, C.; Goyena, E.; Berriatua, E. Quantitative estimation of the impact of caprine arthritis encephalitis virus infection on milk production by dairy goats. Vet. J. 2013, 197, 311–317. [Google Scholar] [CrossRef] [PubMed]
  14. Sihvonen, L.; Nuotio, L.; Rikula, U.; Hirvelä-Koski, V.; Kokkonen, U. Preventing the spread of maedi-visna in sheep through a voluntary control programme in Finland. Prev. Vet. Med. 2000, 47, 213–220. [Google Scholar] [CrossRef]
  15. Tavella, A.; Bettini, A.; Ceol, M.; Zambotto, P.; Stifter, E.; Kusstatscher, N.; Lombardi, R.; Nardeli, S.; Beato, M.S.; Capello, K.; et al. Achievements of an eradication programme against caprine arthritis encephalitis virus in South Tyrol; Italy. Vet. Rec. 2018, 182, 51. [Google Scholar] [CrossRef] [PubMed]
  16. Deubelbeiss, M.; Blatti-Cardinaux, L.; Zahno, M.L.; Zanoni, R.; Vogt, H.R.; Posthaus, H.; Bertoni, G. Characterization of small ruminant lentivirus A4 subtype isolates and assessment of their pathogenic potential in naturally infected goats. Virol. J. 2014, 11, 65. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. Pérez, M.; Biescas, E.; Reina, R.; Glaria, I.; Marín, B.; Marquina, A.; Salazar, E.; Álvarez, N.; de Andrés, D.; Fantova, E.; et al. Small ruminant lentivirus-induced arthritis: Clinicopathologic findings in sheep infected by a highly replicative SRLV B2 genotype. Vet. Pathol. 2015, 52, 132–139. [Google Scholar] [CrossRef] [PubMed]
  18. Angelopoulou, K.; Poutahidis, T.; Brellou, G.D.; Greenland, T.; Vlemmas, I. A deletion in the R region of long terminal repeats in small ruminant lentiviruses is associated with decreased pathology in the lung. Vet. J. 2008, 175, 346–355. [Google Scholar] [CrossRef] [PubMed]
  19. Preziuso, S.; Taccini, E.; Rossi, G.; Renzoni, G.; Braca, G. Experimental Maedi Visna Virus Infection in sheep: A morphological; immunohistochemical and PCR study after three years of infection. Eur. J. Histochem. 2003, 47, 373–378. [Google Scholar] [CrossRef] [PubMed]
  20. Ryan, S.; Tiley, L.; McConnell, I.; Blacklaws, B. Infection of dendritic cells by the Maedi-Visna lentivirus. J. Virol. 2000, 74, 10096–10103. [Google Scholar] [CrossRef] [PubMed]
  21. Edgar, R.C. MUSCLE: Multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 2004, 32, 1792–1797. [Google Scholar] [CrossRef] [PubMed]
  22. Kumar, S.; Stecher, G.; Tamura, K. MEGA7: Molecular Evolutionary Genetics Analysis Version 7.0 for Bigger Datasets. Mol. Biol. Evol. 2016, 33, 1870–1874. [Google Scholar] [CrossRef] [PubMed]
  23. Nei, M.; Kumar, S. Molecular Evolution and Phylogenetics; Oxford University Press: New York, NY, USA, 2000. [Google Scholar]
  24. Kimura, M. A simple method for estimating evolutionary rates of base substitutions through comparative studies of nucleotide sequences. J. Mol. Evol. 1980, 16, 111–120. [Google Scholar] [CrossRef] [PubMed]
  25. Glaria, I.; Reina, R.; Ramírez, H.; de Andrés, X.; Crespo, H.; Jauregui, P.; Salazar, E.; Luján, L.; Pérez, M.M.; Benavides, J.; et al. Visna/Maedi virus genetic characterization and serological diagnosis of infection in sheep from a neurological outbreak. Vet. Microbiol. 2012, 155, 137–146. [Google Scholar] [CrossRef] [PubMed]
  26. Sutton, K.A.; Lin, C.T.; Harkiss, G.D.; McConnell, I.; Sargan, D.R. Regulation of the long terminal repeat in visna virus by a transcription factor related to the AML/PEBP2/CBF superfamily. Virology 1997, 229, 240–250. [Google Scholar] [CrossRef] [PubMed]
  27. Kuhar, U.; Barlic-Maganja, D.; Grom, J. Phylogenetic analysis of small ruminant lentiviruses detected in Slovenia. Vet. Microbiol. 2012, 162, 201–206. [Google Scholar] [CrossRef] [PubMed]
  28. Olech, M.; Rachid, A.; Croise, B.; Kuźmak, J.; Valas, S. Genetic and antigenic characterization of small ruminant lentiviruses circulating in Poland. Virus Res. 2012, 163, 528–536. [Google Scholar] [CrossRef] [PubMed]
  29. Fras, M.; Leboeuf, A.; Labrie, F.M.; Laurin, M.A.; Singh Sohal, J.; L’Homme, Y. Phylogenetic analysis of small ruminant lentiviruses in mixed flocks: Multiple evidence of dual infection and natural transmission of types A2 and B1 between sheep and goats. Infect. Genet. Evol. 2013, 19, 97–104. [Google Scholar] [CrossRef] [PubMed]
  30. Santry, L.A.; de Jong, J.; Gold, A.C.; Walsh, S.R.; Menzies, P.I.; Wootton, S.K. Genetic characterization of small ruminant lentiviruses circulating in naturally infected sheep and goats in Ontario, Canada. Virus Res. 2013, 175, 30–44. [Google Scholar] [CrossRef] [PubMed]
  31. Reina, R.; Mora, M.I.; Glaria, I.; García, I.; Solano, C.; Luján, L.; Badiola, J.J.; 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. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Glaria, I.; Reina, R.; Crespo, H.; de Andrés, X.; Ramírez, H.; Biescas, E.; Pérez, M.M.; Badiola, J.; Luján, L.; Amorena, B.; et al. Phylogenetic analysis of SRLV sequences from an arthritic sheep outbreak demonstrates the introduction of CAEV-like viruses among Spanish sheep. Vet. Microbiol. 2009, 138, 156–162. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Giammarioli, M.; Bazzucchi, M.; Puggioni, G.; Brajon, G.; Dei Giudici, S.; Taccori, F.; Feliziani, F.; De Mia, G.M. Phylogenetic analysis of small ruminant lentivirus (SRLV) in Italian flocks reveals the existence of novel genetic subtypes. Virus Genes 2011, 43, 380–384. [Google Scholar] [CrossRef] [PubMed]
  34. Cardinaux, L.; Zahno, M.L.; Deubelbeiss, M.; Zanoni, R.; Vogt, H.R.; Bertoni, G. Virological and phylogenetic characterization of attenuated small ruminant lentivirus isolates eluding efficient serological detection. Vet. Microbiol. 2013, 162, 572–581. [Google Scholar] [CrossRef] [PubMed]
  35. Purdy, J.G.; Flanagan, J.M.; Ropson, I.J.; Rennoll-Bankert, K.E.; Craven, R.C. Critical role of conserved hydrophobic residues within the major homology region in mature retroviral capsid assembly. J. Virol. 2008, 82, 5951–5961. [Google Scholar] [CrossRef] [PubMed]
  36. Chu, H.H.; Chang, Y.F.; Wang, C.T. Mutations in the alpha-helix directly C-terminal to the major homology region of human immunodeficiency virus type 1 capsid protein disrupt Gag multimerization and markedly impair virus particle production. J. Biomed. Sci. 2006, 13, 645–656. [Google Scholar] [CrossRef] [PubMed]
  37. Angelopoulou, K.; Brellou, G.D.; 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. [Google Scholar] [CrossRef] [PubMed]
  38. Kokawa, S.; Oba, M.; Hirata, T.; Tamaki, S.; Omura, M.; Tsuchiaka, S.; Nagai, M.; Omatsu, T.; Mizutani, T. Molecular characteristics and prevalence of small ruminant lentiviruses in goats in Japan. Arch. Virol. 2017, 162, 3007–3015. [Google Scholar] [CrossRef] [PubMed]
Viruses 10 00315 g001 550
Figure 1. SRLV lesions in mammary gland of sheep. (A) Minimal lesion (+) with small focal lymphocyte aggregates within the mammary interstitium. 40×. (B) Black arrow indicates moderate (++) focal inflammatory lesion surrounded by minimal lymphocytic infiltrates. 40×. (C) Large lymphoid follicle in a severe lesion (+++). Hematoxylin and eosin (HE) staining. 20×.
Figure 1. SRLV lesions in mammary gland of sheep. (A) Minimal lesion (+) with small focal lymphocyte aggregates within the mammary interstitium. 40×. (B) Black arrow indicates moderate (++) focal inflammatory lesion surrounded by minimal lymphocytic infiltrates. 40×. (C) Large lymphoid follicle in a severe lesion (+++). Hematoxylin and eosin (HE) staining. 20×.
Viruses 10 00315 g001
Viruses 10 00315 g002 550
Figure 2. IHC against p28 of CAEV/VMV in mammary gland. (A) Scattered positive cells within minimal inflammatory lesions (+). 20×. (B) Positive macrophage-like cell in moderate lesion (++). 40×.
Figure 2. IHC against p28 of CAEV/VMV in mammary gland. (A) Scattered positive cells within minimal inflammatory lesions (+). 20×. (B) Positive macrophage-like cell in moderate lesion (++). 40×.
Viruses 10 00315 g002
Viruses 10 00315 g003 550
Figure 3. Phylogenetic analysis of the SRLV partial gag-pol region. Sequences of different SRLV genotypes and subtypes available in GenBank were used as reference isolates. Reference sequences are indicated with their accession number and country of origin (CAN: Canada; CHI: China; ENG: England; ICE: Iceland; ITA: Italy; MEX: Mexico; NOR: Norway; POR: Portugal; SOA: South Africa; SPA: Spain; SWI: Switzerland; USA: the U.S.A.). Samples are indicated with their codes and are labeled on the basis of the score of the mammary lesions observed (▲ severe, ■ moderate, ▼ minimal, ● no lesions, ♦ not available). The evolutionary history was inferred by using the Maximum Likelihood method based on the General Time Reversible model with a gamma distribution with invariant sites and with bootstrap values based on 1000 repetitions. The tree is unrooted.
Figure 3. Phylogenetic analysis of the SRLV partial gag-pol region. Sequences of different SRLV genotypes and subtypes available in GenBank were used as reference isolates. Reference sequences are indicated with their accession number and country of origin (CAN: Canada; CHI: China; ENG: England; ICE: Iceland; ITA: Italy; MEX: Mexico; NOR: Norway; POR: Portugal; SOA: South Africa; SPA: Spain; SWI: Switzerland; USA: the U.S.A.). Samples are indicated with their codes and are labeled on the basis of the score of the mammary lesions observed (▲ severe, ■ moderate, ▼ minimal, ● no lesions, ♦ not available). The evolutionary history was inferred by using the Maximum Likelihood method based on the General Time Reversible model with a gamma distribution with invariant sites and with bootstrap values based on 1000 repetitions. The tree is unrooted.
Viruses 10 00315 g003
Viruses 10 00315 g004a 550Viruses 10 00315 g004b 550
Figure 4. Alignment (MUSCLE) of deduced amino acid sequences of partial gag-p25 of the SRLV sequences obtained and of some SRLV reference strains; (a) positions from 209 to 278, (b) posiions from 279 to 368. Two immuno-dominant epitopes of this capsid protein, the GG motif, and the Major Homology Region (MHR) are within squares. The score of mammary lesion of each sample is reported. Legend: (·) homology, (−) deletion, (+++) severe lesions, (++) moderate lesions, (+) minimal lesions, (no) no lesions, (n. c.) not classified, CAN: Canada; CHI: China; ENG: England; ICE: Iceland; ITA: Italy; MEX: Mexico; NOR: Norway; POR: Portugal; SOA: South Africa; SPA: Spain; SWI: Switzerland; USA.
Figure 4. Alignment (MUSCLE) of deduced amino acid sequences of partial gag-p25 of the SRLV sequences obtained and of some SRLV reference strains; (a) positions from 209 to 278, (b) posiions from 279 to 368. Two immuno-dominant epitopes of this capsid protein, the GG motif, and the Major Homology Region (MHR) are within squares. The score of mammary lesion of each sample is reported. Legend: (·) homology, (−) deletion, (+++) severe lesions, (++) moderate lesions, (+) minimal lesions, (no) no lesions, (n. c.) not classified, CAN: Canada; CHI: China; ENG: England; ICE: Iceland; ITA: Italy; MEX: Mexico; NOR: Norway; POR: Portugal; SOA: South Africa; SPA: Spain; SWI: Switzerland; USA.
Viruses 10 00315 g004aViruses 10 00315 g004b
Viruses 10 00315 g005 550
Figure 5. Phylogenetic analysis of the SRLV partial LTR. Sequences of different SRLV genotypes and subtypes available in GenBank were used as reference isolates. Reference sequences are indicated with their accession number and country of origin (CAN: Canada; CHI: China; ENG: England; ICE: Iceland; ITA: Italy; MEX: Mexico; NOR: Norway; POR: Portugal; SOA: South Africa; SPA: Spain; SWI: Switzerland; USA). Samples are indicated with their codes and are labeled on the basis of the score of the mammary lesions observed (▲ severe, ■ moderate, ▼ minimal, ● no lesions, ♦ not available). The phylogenetic analysis was performed with a maximum likelihood (ML) method using the Kimura 2-parameter model with a gamma distribution and with bootstrap values based on 1000 repetitions. Sequences are not deposited in GenBank because they do not reach the minimum length required.
Figure 5. Phylogenetic analysis of the SRLV partial LTR. Sequences of different SRLV genotypes and subtypes available in GenBank were used as reference isolates. Reference sequences are indicated with their accession number and country of origin (CAN: Canada; CHI: China; ENG: England; ICE: Iceland; ITA: Italy; MEX: Mexico; NOR: Norway; POR: Portugal; SOA: South Africa; SPA: Spain; SWI: Switzerland; USA). Samples are indicated with their codes and are labeled on the basis of the score of the mammary lesions observed (▲ severe, ■ moderate, ▼ minimal, ● no lesions, ♦ not available). The phylogenetic analysis was performed with a maximum likelihood (ML) method using the Kimura 2-parameter model with a gamma distribution and with bootstrap values based on 1000 repetitions. Sequences are not deposited in GenBank because they do not reach the minimum length required.
Viruses 10 00315 g005
Viruses 10 00315 g006 550
Figure 6. Alignment (MUSCLE) of nucleotide sequences of the LTR region of the SRLV sequences obtained and of some SRLV reference strains. LTR sequence of the isolate Icelandic 1514 (M10609.1) was used as a reference; (a) positions from 9057 to 9126, (b) posiions from 9127 to 9206. Sequence of reference strains EV1 (S51392.1) was also used. Sequences previously found in Spain in goats (C3—DQ084355.1) and in sheep (697—HQ864615.1, 160—HQ864610.1, 292—HQ864611.1 and 368—HQ864614.1) and showing high homology with our samples were also included. Sequences of the AML (vis), the TATA box element, and the polyadenylation signal are within black squares. Sequences of samples with similar scores of mammary lesions are within colored squares. Legend: (·) homology, (−) deletion, (blue squares) severe lesions, (red squares) moderate lesions, (light blue squares) minimal lesions, and (green squares) no lesions.
Figure 6. Alignment (MUSCLE) of nucleotide sequences of the LTR region of the SRLV sequences obtained and of some SRLV reference strains. LTR sequence of the isolate Icelandic 1514 (M10609.1) was used as a reference; (a) positions from 9057 to 9126, (b) posiions from 9127 to 9206. Sequence of reference strains EV1 (S51392.1) was also used. Sequences previously found in Spain in goats (C3—DQ084355.1) and in sheep (697—HQ864615.1, 160—HQ864610.1, 292—HQ864611.1 and 368—HQ864614.1) and showing high homology with our samples were also included. Sequences of the AML (vis), the TATA box element, and the polyadenylation signal are within black squares. Sequences of samples with similar scores of mammary lesions are within colored squares. Legend: (·) homology, (−) deletion, (blue squares) severe lesions, (red squares) moderate lesions, (light blue squares) minimal lesions, and (green squares) no lesions.
Viruses 10 00315 g006
Table
Table 1. List of ovine mammary gland samples collected for this study. Samples are classified on the basis of the grade of mastitis observed by histopathology. Lesions in sample M5 were not classified due to concomitant purulent mastitis was present. “+”: positive result, “−”: negative result.
Table 1. List of ovine mammary gland samples collected for this study. Samples are classified on the basis of the grade of mastitis observed by histopathology. Lesions in sample M5 were not classified due to concomitant purulent mastitis was present. “+”: positive result, “−”: negative result.
SampleGrade of MastitisIHCLTR PCRGag-pol PCRGenotype
M1moderate+++B2
M2minimal++−−
M3severe+++B2
M4moderate++−−
M5not classified+++A2/A3
M6moderate++−−
M7moderate+++−
M8moderate+++−
M9severe+++−
M10minimal++−−
M11moderate+++−
M12minimal+++A2/A3
M13no−−−−
M14moderate++−−
M15minimal+++A2/A3
M16minimal++−−
M17moderate+++B2
M18no−−−−
M19moderate+++A2/A3
Q1no+++B2
Q2no−−−−
Q3minimal+−−−
Q4minimal++−−
Q5minimal++−−
Q6severe++−−
Q7no+++A2/A3
Q8severe+++A2/A3
Q9moderate++−−
Q10severe+++A2/A3
Q11moderate++−−
Q12moderate+++−
Q13minimal+−−−
Q14moderate++−−
Q15minimal++−−
Q16minimal++−−
N16-426severe+++A2/A3
N17-44moderate+++A2/A3
Table
Table 2. Pairwise nucleotidic genetic distances (p-distance model) of the partial gag-pol region of some SRLV reference strains and SRLV strains sequenced in this study.
Table 2. Pairwise nucleotidic genetic distances (p-distance model) of the partial gag-pol region of some SRLV reference strains and SRLV strains sequenced in this study.
SampleGenotypeM5M12M15M19Q7Q8Q10N16-426N17-44M1M3M17Q1
M12 0.112-
M15 0.0130.013-
M19 0.0130.1220.104-
Q7 0.1300.1330.1290.109-
Q8 0.1000.1150.0970.1220.122-
Q10 0.1090.1190.1050.1290.1270.105-
N16-426 0.1270.1380.1290.1150.1370.1270.135-
N17-44 0.1520.1550.1520.1380.1630.1570.1520.094-
M1 0.2080.2130.2080.2130.2160.2090.2160.2170.216-
M3 0.2160.2220.2110.2210.2310.2110.2210.2240.2170.040-
M17 0.2140.2210.2090.2260.2320.2170.2170.2310.2270.0480.041-
Q1 0.2090.2160.2030.2170.2220.2210.2190.2160.2160.0810.0860.076-
M10608.1A10.1710.1600.1680.1370.1570.1660.1450.1750.1750.2140.2080.2290.224
S51392A10.1730.1650.1680.1730.1810.1710.1630.1760.1680.2420.2470.2450.244
AY101611.1A20.1550.1350.1530.1500.1470.1530.1520.1650.1650.1860.1880.1910.201
HQ848062.1A2/A30.1170.1350.1200.1050.1220.1270.1350.1420.1240.2060.2090.2190.213
AY454176.1A30.1300.1380.1300.1300.1370.1330.1450.1450.1520.1780.1910.1930.198
AY454161.1A40.1380.1660.1350.1530.1850.1430.1500.1730.1710.2110.2160.2140.214
AY454175.1A50.1530.1450.1550.1480.1710.1600.1470.1760.1650.2090.2030.2080.206
AY454208.1A70.1760.1570.1700.1450.1580.1580.1660.1610.1400.1980.2090.2010.209
FR694908A90.1660.1680.1700.1660.1780.1800.1710.1810.1610.2130.2240.2240.231
FR693825A110.1830.1710.1800.1520.1570.1760.1830.1890.1780.2130.2170.2310.232
M33677B10.2040.2140.2010.2110.2240.2110.2110.2220.2110.1200.1190.1200.129
FJ195346.1B20.2260.2310.2210.2140.2320.2260.2370.2210.2170.0920.0920.1020.100
EU010126.1B20.2080.2220.2060.2260.2320.2190.2270.2240.2160.0640.0680.0680.081
JF502417.1B30.2130.2210.2140.2090.2170.2190.2320.2130.2240.1850.1860.1960.178
AF322109.1C0.2620.2520.2570.2500.2520.2450.2550.2440.2360.2130.2140.2140.216
EU293537.2E10.2920.2950.2920.2970.3010.2930.3000.2930.2900.2930.2850.2900.303
GQ381130.1E20.2900.3010.2930.2900.3130.2930.2980.2880.2900.2900.2900.2970.293

Share and Cite

MDPI and ACS Style

Gayo, E.; Cuteri, V.; Polledo, L.; Rossi, G.; García Marín, J.F.; Preziuso, S. Genetic Characterization and Phylogenetic Analysis of Small Ruminant Lentiviruses Detected in Spanish Assaf Sheep with Different Mammary Lesions. Viruses 2018, 10, 315. https://doi.org/10.3390/v10060315

AMA Style

Gayo E, Cuteri V, Polledo L, Rossi G, García Marín JF, Preziuso S. Genetic Characterization and Phylogenetic Analysis of Small Ruminant Lentiviruses Detected in Spanish Assaf Sheep with Different Mammary Lesions. Viruses. 2018; 10(6):315. https://doi.org/10.3390/v10060315

Chicago/Turabian Style

Gayo, Elena, Vincenzo Cuteri, Laura Polledo, Giacomo Rossi, Juan F. García Marín, and Silvia Preziuso. 2018. "Genetic Characterization and Phylogenetic Analysis of Small Ruminant Lentiviruses Detected in Spanish Assaf Sheep with Different Mammary Lesions" Viruses 10, no. 6: 315. https://doi.org/10.3390/v10060315

APA Style

Gayo, E., Cuteri, V., Polledo, L., Rossi, G., García Marín, J. F., & Preziuso, S. (2018). Genetic Characterization and Phylogenetic Analysis of Small Ruminant Lentiviruses Detected in Spanish Assaf Sheep with Different Mammary Lesions. Viruses, 10(6), 315. https://doi.org/10.3390/v10060315

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

Article Metrics

Back to TopTop