Next Article in Journal
Pandemic Risk Assessment for Swine Influenza A Virus in Comparative In Vitro and In Vivo Models
Next Article in Special Issue
The Development of a Sensitive Droplet Digital Polymerase Chain Reaction Test for Quantitative Detection of Goose Astrovirus
Previous Article in Journal
No Remdesivir Resistance Observed in the Phase 3 Severe and Moderate COVID-19 SIMPLE Trials
Previous Article in Special Issue
Overview of Modern Commercial Kits for Laboratory Diagnosis of African Swine Fever and Swine Influenza A Viruses
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Brief Report

Genetic Characterization of Small Ruminant Lentiviruses Isolated from Dairy Sheep in Greece

by
Aphrodite I. Kalogianni
1,*,
Ilias Bouzalas
2,
Sofia Marka
3,
Maria-Eleftheria Zografaki
3,
Sofia Mavrikou
3 and
Athanasios I. Gelasakis
1
1
Laboratory of Anatomy and Physiology of Farm Animals, Department of Animal Science, School of Animal Biosciences, Agricultural University of Athens (AUA), Iera Odos 75 Str., 11855 Athens, Greece
2
Veterinary Research Institute, Hellenic Agricultural Organization-DEMETER, Campus of Thermi, 57001 Thessaloniki, Greece
3
Laboratory of Cell Technology, Department of Biotechnology, School of Applied Biology and Biotechnology, Agricultural University of Athens (AUA), EU-CONEXUS European University, 11855 Athens, Greece
*
Author to whom correspondence should be addressed.
Viruses 2024, 16(4), 547; https://doi.org/10.3390/v16040547
Submission received: 12 February 2024 / Revised: 15 March 2024 / Accepted: 29 March 2024 / Published: 31 March 2024
(This article belongs to the Special Issue Viral Diseases of Livestock and Diagnostics)

Abstract

:
The high genetic heterogeneity of small ruminant lentiviruses (SRLV) renders the genetic characterization of the circulating strains crucial for the epidemiological investigation and the designation of effective diagnostic tools. In Greece, research data regarding the genetic diversity of the circulating SRLV strains is scarce, hindering the implementation of efficient surveillance and control programs. The objective of the study was to genetically characterize SRLV strains isolated from intensive dairy sheep farms in Greece and evaluate the variability of the immunodominant regions of the capsid protein. For this reason, a total of 12 SRLV-infected animals from four intensive dairy sheep farms with purebred Chios and Lacaune ewes were used for the amplification and sequencing of an 800 bp gag-pol fragment. The phylogenetic analyses revealed a breed-related circulation of strains; Chios ewes were infected with strains belonging exclusively to a separate group of genotype A, whereas strains belonging to subtype B2 were isolated from Lacaune ewes. Immunodominant epitopes of capsid protein were quite conserved among the strains of the same genotype, except for the Major Homology Region which showed some unique mutations with potential effects on viral evolution. The present study contributes to the extension of the current knowledge regarding the genetic diversity of SRLV strains circulating in sheep in Greece. However, broader genetic characterization studies are warranted for the exploration of possible recombinant events and the more comprehensive classification of the circulating strains.

1. Introduction

The Maedi-visna virus (MVV) and caprine arthritis-encephalitis virus (CAEV) belong to small ruminant lentiviruses (SRLV), a subgroup of viruses of the family Retroviridae which cause persistent and insidious infections in sheep and goats [1]. These infections can lead to the manifestation of clinical signs such as pneumonia, mastitis, arthritis, and encephalitis [2], compromising health, welfare, and productivity of the affected animals. The global spreading of SRLV infections renders their control imperative for the limitation of the monetary losses in the small ruminant industry. Considering the lack of effective treatment and vaccines against SRLV infections, the eradication of SRLV infections is based on the early diagnosis and the implementation of general management practices and biosecurity measures to prevent the spread of the virus between and within herds [3].
Nevertheless, the development of a “gold standard” assay for the early diagnosis of SRLV infections is a challenging endeavor considering the high genetic variability of the virus [4]. The genetic characterization of the circulating SRLV strains is an important prerequisite for the implementation of efficient surveillance and eradication programs, not only for the “mapping” of the spatial spreading of various strains, but also for the development of sensitive and specific diagnostic tools for application on a case-specific basis. For this reason, many countries have attempted the molecular characterization of the circulating SRLV strains; the phylogenetic analyses, genotyping, and classification of the isolates have identified five genotypes (A–E) and many subtypes thereof in sheep and goats [5,6,7,8,9,10]. Among them, genotypes A and B have a global spread and have been initially isolated from sheep and goats, respectively, although cross-species transmission has been previously reported in many countries [11,12,13]. On the other hand, genotypes C, D, and E have been identified in specific regions and species, and their spreading is considered more restricted compared to genotypes A and B [7,9,14,15].
Phylogenetic analyses contribute significantly to the investigation of the SRLV genetic diversity, the tracing of transmission patterns, and the assessment of the genetic distance and relationships among different viral isolates. The commonly targeted genomic regions for phylogenetic analyses in SRLV are the three major coding regions of the SRLV genome, namely, gag, pol, and env genes, as well as the non-coding LTR regions [8]. Among these genomic regions, the gag gene, which encodes structural proteins including the matrix (p16MA), capsid (p25CA), and nucleocapsid (p14NC) proteins, is highly preferred due to its conserved function. The viral capsid protein contains three immunodominant epitopes, namely, epitope 2, including the double glycine “GG” motif; epitope 3; and the Major Homology Region (MHR), which determine the genetic variability between the SRLV strains [16]. Considering that the viral capsid protein is commonly used as an antigen for the detection of SRLV-specific antibodies [17,18,19], the high divergence in linear epitopes may adversely affect the diagnostic performance of serological diagnostic tools [7].
Although SRLV infections have been previously reported in Greece [20,21,22,23], relevant research data regarding circulating strains and their genetic heterogeneity are scarce [24,25,26]. Hence, the objective of this study was the molecular characterization of SRLV strains isolated from dairy sheep in Greece based on the partial characterization of the gag-pol genomic region, with the genetic analyses presented including phylogenetic analyses and heterogeneity of linear epitopes of capsid protein between the different strains.

2. Materials and Methods

2.1. Animals, Blood Samples, and DNA Extraction

A total of seven ewes that were constantly found seropositive and 13 ewes that were constantly found seronegative with a commercial indirect ELISA test (ELISA, CAEV/MVV Total Ab Test, IDEXX, Westbrook, ME, USA) during a two-year epidemiological study conducted in four intensive dairy sheep farms in Greece [23,27] were used for the purpose of the present study. From each animal, a whole blood sample (~9 mL) was collected from the jugular vein in an EDTA-anticoagulated tube and was further processed for the molecular analysis. The protocol of this study was approved by the Animal Research Ethics Committee of the Agricultural University of Athens following the national animal welfare regulations (protocol no 9).
Leukocyte pellets were isolated from whole blood samples. A total of 2 mL blood was mixed with 13 mL of ACK lysis buffer (8.02 g NH4Cl, 0.84 g NaHCO3, and 0.37 g EDTA per liter, pH 7.2–7.4) in a FalconTM 15 mL conical centrifuge tube. After 15 min of incubation at room temperature, the mixtures were centrifuged at 450× g for 15 min. The supernatant was discarded, and the leukocyte pellets were diluted in 1 mL ACK lysis buffer and transferred in 2 mL microcentrifuge tubes. The mixtures were centrifuged at 450× g for 4 min and the supernatant was discarded. The leukocyte pellets were resuspended in 1 mL PBS (phosphate buffered saline, pH = 7.4), and were centrifuged at 650× g for 4 min. The leukocyte pellets were used for the genomic DNA extraction with a commercial kit (PureLink® Genomic DNA Kit, Life technologies Corp., Carlsbad, CA, USA), according to the manufacturer’s instructions. The procedure of DNA extraction was based on the selective binding of DNA to a silica-based membrane in the presence of chaotropic salts (ethanol and guanidinium hydrochloride). DNA concentration and purity were measured in a spectrophotometer (Q5000, Quawell Technology, Inc., San Jose, CA, USA) considering A260/A280 and A260/A230 ratios. Samples with A260/280 and A260/230 ratios of ~1.8 and 2.0–2.2, respectively, were considered as pure for DNA and acceptable for PCR and sequencing. DNA samples were stored at −20 °C until they were ready to be assayed.

2.2. PCR, Sequencing, and Phylogenetic Analyses

A nested PCR protocol, previously described by Grego et al., 2007 [7], was used after a slight modification. In brief, the first round of PCR amplified a 1.3 kb DNA fragment of the gag-pol region. The reaction was carried out in 50 μL final volume which contained a 5 μL DNA sample (~1 μg), 400 nM of each primer (F1: TGGTGARKCTAGMTAGAGACATG; R1: CATAGGRGGHGCGGACGGCASCA), 25 μL of OneTaq 2× Master Mix (New England, Biolabs Inc., Hitchin, UK), and 18 μL of DNase-free water. The amplification protocol consisted of initial denaturation at 94 °C for 5 min, followed by 40 cycles of denaturation at 95 °C for 30 s, annealing at 55 °C for 30 s, and extension at 72 °C for 1.5 min, terminating with a final extension at 72 °C for 5 min. The second round of PCR was carried out with 4 μL from PCR amplicons from the first PCR and a different primer set (F2: CAAACWGTRGCAATGCAGCATGG; R2: GCGGACGGCASCACACG) for the amplification of a ~800 bp gag-pol gene fragment (988–1802 positions in genome of reference strain K1514 with GenBank accession number M60609). The amplification protocol of the second PCR consisted of initial denaturation at 94 °C for 5 min, followed by 40 cycles of denaturation at 95 °C for 20 s, annealing at 60 °C for 30 s, and extension at 72 °C for 1.5 min, terminating with a final extension at 72 °C for 5 min. All amplifications were carried out in the same thermal cycler (Labcycler, Sensoquest GmbH, Göttingen, Germany) in duplicate, and a negative control reaction with DNase-free water instead of DNA template was used each time to determine any possible nucleic acid contamination. All PCR products were mixed with DNA Gel Loading Dye (Invitrogen, Vilnius, Lithuania) and were electrophoresed in 1.5% agarose gel in 1× RNase-free TAE buffer (Invitrogen, Vilnius, Lithuania), stained with SYBRsafe DNA gel staining (Invitrogen, Carlsbad, CA, USA), and visualized under ultraviolet (UV) light. For the evaluation of the amplicon length, a 100 bp ladder was used (GeneRuler 100 bp DNA Ladder, Thermoscientific, Vilnius, Lithuania).
PCR bands of desired length were gel-extracted (PureLink™ Quick Gel Extraction Kit, Invitrogen, Vilnius, Lithuania) and sequenced in both directions with Sanger dideoxy sequencing method on an ABI PRISM 3730xl Genetic Analyzer in an external laboratory (Cemia SA, Larissa, Greece) using the primers F2 and R2.
All the phylogenetic analyses were performed using the MEGA11 software [28]. The nucleotide sequences were trimmed according to the obtained chromatographs and the consensus sequences generated were submitted to GenBank with accession numbers OR283217–OR283228. After this, a similarity search was conducted with the NCBI (National Center of Biotechnology Information) database using BLAST (Basic Local Alignment Search Tool), and the sequences most homologous to them and all the representative sequences from SRLV strains of genotypes A, B, C, and E from different geographical regions were aligned by CLUSTAL W [29]. After the alignment, pairwise genetic distances were calculated with the p-distance model applying the gamma distribution parameter between the nucleotide sequences. A phylogenetic tree was constructed using the Maximum Likelihood (ML) method [30] and the General Time Reversiblemodel with gamma distribution (G) and invariant sites (I) with 100 bootstrap replicates [31] based on the lowest Bayesian Information Criterion (BIC) [32].
Also, nucleotide sequences were translated into amino acid sequences and aligned using the MEGA v.11 software, for the calculation of their pairwise genetic distance (p-distance with gamma distribution parameter). Also, the immunodominant regions of gag protein (epitope 2, epitope 3, and MHR) were comparatively assessed.

3. Results

3.1. Pairwise Nucletide and Amino Acid Sequence Comparisons

A total of 12 animals were found positive with the nested PCR protocol, and their PCR products were sequenced. In Table 1, information about the breed, age, and serological status of these animals, the accession number of the respective sequences, the farms, and their geographical location are summarized. All PCR negative animals (8/20) were constantly seronegative during the study.
The mean pairwise genetic distance of partial nucleotide gag-pol sequences obtained in this study was 15.0% and varied from 1.2% between the s-8m and s-98m sequences to 29.2% between the s-40sx and s-42fl sequences (Table 2). Also, the intra-farm pairwise comparisons were 19.0% for farm A (s-16sx, s-56sx, and s-40sx sequences), 5.0% for farm B (s-11s, s-14s, s-35s, and s-83s sequences), 7.0% for farm C (s-10fl, s-21fl, and s-42fl sequences), and 1.2% for farm D (s-8m and s-98m, sequences), with an average distance of 8.1% (inter-farm mean genetic distance).
The overall mean pairwise genetic distance of partial gag-pol sequences from the present study, the SRLV representative reference strains, and the strains isolated from nearby countries (Italy, Turkey, Jordan, and Lebanon) was 24.0%. As presented in Table 3, gag-pol sequences from the present study presented the highest genetic distance with the 1GA strain (genotype C) (23.5–32.2%), and the Roccaverano and Seui strains (genotype E) (35.5–41.4%). The strains s-56sx, s-11s, s-14s, s-35s, s-83s, s-10fl, s-21fl, and s-42fl were the most similar to the Greek strains SRLV-Greece-S1 and SRLV-Greece-S2 (genotype A) (genetic distance 9.1–12.6%), whereas the strains s-16sx, s-40sx, s-8m, and s-98m were most related to Italian strains SRLV001, SRLV042, lt-Pi1, and lt-007 of genotype B2 (genetic distance 6.2–8.9%).

3.2. Phylogenetic Analysis

The gag-pol sequences isolated in the present study were clustered in two groups in the phylogenetic tree (Figure 1). The first group included s-8m, s-98m, s-16sx, and s-40sx sequences and was strongly associated with the Spanish (FJ195346) and the Italian (MG554402, EU010126, AY265456, and MH374288) strains of genotype B2. The other group included s-10fl, s-21fl, s-42fl, s-11s, s-14s, s-35s, s-83s, and s-56sx sequences and was clustered with other Greek (AY530289 and AY530290) strains of genotype A.

3.3. Comparative Assessment of immunodominant Linear Epitopes of Capsid Protein

In Figure 2a–c, the alignment of the amino acid sequences is presented, and epitope 2, including the double glycine “GG” motif; epitope 3; and the MHR are comparatively assessed. Strains from the present study belonging to genotype A, according to the afore-mentioned phylogenetic tree, demonstrated high conservation regarding epitopes 2 and 3, whereas more alterations were found in MHR. Also, all strains belonging to genotype A presented the asparagine-valine “NV” motif, whereas the strains of genotype B presented the typical double glycine “GG” motif. Also, in the s-42fl strain, glutamic acid (E) was replaced by aspartic acid (D) in epitope 2, and arginine (R) and serine (S) were replaced by glycine (G) and tyrosine (Y), respectively, in epitope 3. In the MHR, i) asparagine (N) was replaced by serine (S) in the s-42fl strain and by lysine (K) in the s-10fl and s-21fl strains; ii) aspartic acid (D) was replaced by asparagine (N) in the s-10fl and s-21fl strains, and by glutamic acid (E) in the s-35s and s-42fl strains; iii) glutamine (Q) was replaced by histidine (H) in s-11s strains; and iv) lysine (K) was replaced by arginine (R) in s-35s strains. Strains from the present study belonging to genotype B2 presented identical amino acid residues in all epitopes, except for the s-40sx strain in epitope 2, where 3 amino acids were substituted: tryptophan (W) with glycine (G), and the next two consecutive arginines (R) with two glutamic acids (E). Also, strains from the present study belonging to genotype A demonstrated high similarity with a few discrepancies with the reference strains of the same genotype and with the strains that presented low genetic distance in the phylogenetic tree, while strains from the present study belonging to genotype B2 demonstrated absolute identity with the other reference strains belonging to the same genotype, except for the s-40sx strain.

4. Discussion

This is the first time that an extended gag-pol fragment (~800 bp), including all the immunodominant regions of capsid protein, from SRLV strains isolated in Greece is analyzed for their genetic characterization. Although numerous sequencing and phylogenetic analyses have been performed for the genotyping and classification of SRLV strains worldwide, based on different regions of the viral genome (mainly gag, pol, and env genes) [5,6,7,8,9,13,33,34,35,36,37,38,39], phylogenetic data regarding circulating SRLV strains in Greece are limited and derived from previous studies based on partial pol (~390 bp) [25,26] and gag (574 bp) [24] sequences, which in the latter case did not include all the immunodominant regions of capsid protein in contrast to our study.
In the present study, all the seropositive animals also tested PCR positive (infected seropositive animals), while five seronegative animals were found PCR positive (infected seronegative animals), indicating the significance of the combination of serological and molecular testing for the detection of SRLV-infected animals [25,40,41,42,43]. Infected seronegative animals were found in both breeds and on all farms. This underpins the possibility of late seroconversion rather than poor diagnostic performance of the serological test due to the circulation strains with antigenic variation. Indeed, the studied antigenic epitopes did not present any alteration that could be linked with the putative misdiagnosed animals. This is also supported by the age distribution of the infected seropositive animals compared to the infected seronegative animals (2.0 to 7.0 vs. 1.5 to 3.0 years, respectively), which indicates the underdiagnosis of the younger animals by using ELISA.
The phylogenetic analysis in the present study was based on the gag-pol region of the viral genome, which is considered relatively conserved (more conserved than the env and less conserved than the pol regions) and suitable for designation of PCR primers, sequencing, and meaningful comparisons between different sequences [18,44]. At the same time, this region encodes for significant viral structural (gag gene) and functional (pol gene) proteins and is under selective pressure presenting stable genetic divergence for phylogenetic analyses and genetic characterization of the isolates. Also, in the present study the size of gag-pol gene fragment was ~800 bp, making it suitable for both sequencing and phylogenetic analysis, as it provides enough sequence information while being manageable in terms of sequencing efforts and computational analysis. The phylogenetic analyses in the present study revealed the circulation of strains belonging to genotypes A and B (subtype B2) in the studied farms. This is consistent with previous studies that have also reported the existence of strains belonging to genotypes A [24,25,45] and B [25] in Greece. Specifically, phylogenetic analyses of partial gag and LTR sequences identified strains more similar to MMV-like strains (A genotype) [24,45], whereas strains from both genotypes A and B were identified with phylogenetic analyses based on partial pol sequences [25].
Gag-pol sequences from our study were clustered into two discrete groups. The sequences from the first group were grouped within the same clade as the strains previously isolated in Greece [24], although the latter were shorter sequences including a 574-nucleotide fragment of the gag gene. Also, the sequences from the present study were closely related to the strains isolated from the nearby countries of the Middle East (Lebanon and Jordan) (subtypes A22–A24) (see Figure 1). The genetic distance between the sequences from the present study and the other Greek strains was lower than 15.0%, whereas it was greater than 15.0% in the case of sequences from subtypes A22–A24. Considering that a 15% genetic distance is considered the limit for the categorization of variants in the same subtype [8], it could be suggested that this indicates the identification of a new subtype into the genotype A, including the Greek strains. On the other hand, the other group of gag-pol sequences was clearly clustered with B2 strains isolated in Italy and Spain [7,46,47,48], presenting very low values of nucleotide distance (<10%) as derived from the nucleotide pairwise comparisons. The genetic association of SRLV strains of genotypes A and B with these specific regions is consistent with previous phylogenetic studies indicating the Middle East and Central Europe as the ancestral regions of the A and B genotypes, respectively [49,50]. Nevertheless, further investigation and bioinformatics analyses of whole-genome sequence data are imperative for the exploration of the origin of SRLV strains isolated in the present study.
Our results also confirmed the infection of sheep with strains of genotype B, which were initially isolated from goats [50]. This finding is in agreement with previous studies, which reported cross-species transmission of A and B genotype strains between sheep and goats [11,12,13]. In particular, subtype B2, which was found in our study, has been reported in sheep in several studies in the past [7,8,46,47]. Moreover, the results obtained from the constructed phylogenetic tree indicated the circulation of strains of both A and B2 genotypes only in the single mixed farm (farm A), where both Chios and Lacaune sheep were reared. On the other hand, in Chios farms (farms B and C), all the sequences were associated with strains of genotype A, whereas the Lacaune farm (farm D) sequences were associated with strains of subtype B2. This is in accordance with the results from mean pairwise genetic distances between the sequences; all sequences originating from the same farm presented low genetic heterogeneity < 15.0%, except for the gag-pol sequences in farm A (19.0%). Also, the identification of strains of genotypes A and B2 in Chios and Lacaune sheep, respectively, suggests a breed-related circulation of strains that may be associated with the origin of animals. Considering that the Chios breed is an indigenous Greek breed with a limited number of breeding stocks across the country, and Lacaune breeding stocks have been massively imported during the last decades from France, the circulation of strains with high genetic similarity among the sheep of the same breed is possibly related to their common origin. This is also confirmed by the fact that in farm A, where Lacaune and Chios sheep were housed together, genotype A was identified only in Chios sheep and genotype B2 only in Lacaune sheep, implying that the initial breeding animals were infected when introduced into the farm, and then their lambs were vertically infected with the same viral strain. However, the possibility of breed-related infectivity among the different SRLV genotypes cannot be excluded. These assumptions need further investigation which must include the molecular characterization of strains isolated from a large number of animals from this farm to explore the possibility of a recombinant event characterizing more viral genomic regions.
Immunodominant epitopes of capsid protein have been previously characterized in many studies for the assessment of their heterogeneity [6,7,9,24,34,51,52,53,54]. Considering that immunodominant linear epitopes of capsid protein are strongly associated with the type-specific host immune response [46,53,55], their role in serological diagnosis is crucial. The investigation of genetic heterogeneity in these epitopes among the different SRLV strains could reveal diagnostic drawbacks of the existent serological diagnostic tools and contribute to the development of highly efficient ones [56,57]. Also, serotyping diagnostic tools discriminating the different SRLV genotypes have been developed using subunits of the viral capsid protein [58,59]. These tests could be used for the implementation of epidemiological studies in countries where the genetic characterization of the circulating SRLV strains has not been performed. Based on this, the exploration of immunodominant linear epitopes of capsid protein in SRLV strains in our country could be exploited for the development of highly sensitive serological tests and serotyping tools instead of the costly and time-consuming genetic analyses. Additionally, it is the first time that a comparative analysis of all the immunodominant linear epitopes of capsid protein is performed in SRLV sequences isolated in Greece. A previous study in our country analyzed shorter amino acid sequences, including only epitope 2 and the MHR [24], without, however, comparing the specific regions of epitopes, instead comparing the total amino acid sequences. Our findings confirmed the presence of the asparagine-valine “NV” and the double glycine “GG” motifs clustered in genotypes A and B, respectively, as previously described in similar studies [34,39,49,54,60,61]. Also, epitopes 2 and 3 were found highly conserved when the strains of genotype A from the present study were compared to the reference strains of the same clade, with the exception of one isolate (s-42fl) in which three new amino acid replacements were identified [one conservative substitution in epitope 2 (D/E) and two non-conservative substitutions in epitope 3 (G/R and Y/S)]. Regarding MHR, in our study, medium heterogeneity was presented among the sequences belonging to genotype A, whereas three new replacements were also observed (H/Q, N/D, and K/N). This finding is in accordance with previous studies, which reported moderate variability in the MHR of isolated strains [39,51,60]. Even though genotype A strains in our study demonstrated high similarity with the strains previously isolated in Greece, the new unique replacements that were identified could lead to a significant divergence related either to a poor diagnostic performance of serological diagnostic tools based on specific epitopes or to alternative infectivity. On the other hand, strains belonging to subtype B2 demonstrated absolute similarity in all epitopes with the reference strains of this subtype, apart from one isolate that demonstrated a specific triple amino acid substitution in epitope 2 (GEE/WRR).
Although it is the first time that the genetic characterization of SRLV isolates from our county included both phylogenetic analyses and assessment of heterogeneity of linear epitopes of capsid protein, the number of the studied animals does not allow the generalization of the results. However, our findings facilitate the extraction of useful conclusions regarding the epidemiology of SRLV infections, as the farms used in the study (i) are representative of the intensive farming system, (ii) are located in different counties of Greece with well-developed dairy sheep farming sectors, (iii) are rearing purebred sheep of the most commonly used dairy sheep breeds in the country (one indigenous and one imported breed). Genetic characterization of different and more extended viral genomic regions or even whole-genome sequencing of isolated strains are warranted for the confirmation of the results from the present study and the comprehensive classification of strains. Also, further studies including both sheep and goats of various breeds and regions of Greece could contribute to (i) the extension of the current knowledge regarding the genetic variability of SRLV strains circulating in different breeds reared under various farming systems in the country, (ii) the elucidation of viral evolution, (iii) the exploration of possible recombination events, and (iv) the designation of effective surveillance and eradication programs against SRLV infections.

5. Conclusions

The present study constitutes the first partial genetic characterization of SRLV strains circulating in dairy sheep in Greece based on an extended gag-pol fragment including all the immunodominant regions of the capsid protein. The phylogenetic analysis revealed the circulation of strains belonging to genotype A, clustered in a distinct clade with the other Greek isolates, and strains grouped clearly with strains of subtype B2. Chios sheep were found to be infected with A genotype, whereas Lacaune sheep were infected with B2 strains, implying a breed-associated circulation of strains and a specific transmission pattern possibly related to the available breeding stocks. Even though amino acid sequences of immunodominant epitopes of capsid protein were sufficiently conserved, especially in B2 strains, remarkable mutations were found in the MHR of A strains, potentially affecting their virulence or challenging serological diagnosis.

Author Contributions

Conceptualization, A.I.K. and A.I.G.; methodology, A.I.K., I.B. and A.I.G.; software, A.I.K.; validation, A.I.K., I.B. and A.I.G.; formal analysis, A.I.K.; investigation, A.I.K. and A.I.G.; data curation, A.I.K., S.M. (Sofia Mavrikou), S.M. (Sofia Marka), and M.-E.Z.; writing—original draft preparation, A.I.K.; writing—review and editing, I.B., S.M. (Sofia Mavrikou), S.M. (Sofia Marka), M.-E.Z. and A.I.G.; visualization, A.I.K.; supervision, A.I.G.; project administration, A.I.K. and A.I.G. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

This study was conducted in accordance with the national animal welfare regulations, i.e., the Presidential Decree 56/13 “Bringing Greek legislation into line with Directive 2010/63/EC of the European Parliament and of the Council of 22 September 2010 (L 276/33/20.10.2010) regarding the protection of animals used for experimental and other scientific purposes”. The applied diagnostic veterinary procedures are not within the context of relevant EU legislation for animal experimentations (Directive 86/609/EC) and may be performed to diagnose animal diseases and improve animal welfare. The protocol of this study was approved by the Animal Research Ethics Committee of the Agricultural University of Athens (protocol code 9 and date of approval 9 January 2022).

Informed Consent Statement

Informed consent was obtained from all farm owners involved in this study.

Data Availability Statement

The data presented in this study are available on request from the corresponding author. The data are not publicly available due to privacy restrictions.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Pépin, M.; Vitu, C.; Russo, P.; Mornex, J.F.; Peterhans, E. Maedi-Visna Virus Infection in Sheep: A Review. Vet. Res. 1998, 29, 341–367. [Google Scholar] [PubMed]
  2. 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]
  3. Peterhans, E.; Greenland, T.; Badiola, J.; Harkiss, G.; Bertoni, G.; Amorena, B.; Eliaszewicz, M.; Juste, R.A.; Kraßnig, 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]
  4. Kalogianni, A.I.; Stavropoulos, I.; Chaintoutis, S.C.; Bossis, I.; Gelasakis, A.I. Serological, Molecular and Culture-based Diagnosis of Lentiviral Infections in Small Ruminants. Viruses 2021, 13, 1711. [Google Scholar] [CrossRef] [PubMed]
  5. L’Homme, Y.; Leboeuf, A.; Arsenault, J.; Fras, M. Identification and Characterization of an Emerging Small Ruminant Lentivirus Circulating Recombinant Form (CRF). Virology 2015, 475, 159–171. [Google Scholar] [CrossRef] [PubMed]
  6. 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]
  7. 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]
  8. 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]
  9. 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]
  10. Kalogianni, A.I.; Bossis, I.; Ekateriniadou, L.V.; Gelasakis, A.I. Etiology, Epizootiology and Control of Maedi-Visna in Dairy Sheep: A Review. Animals 2020, 10, 616. [Google Scholar] [CrossRef]
  11. Leroux, C.; Cruz, J.C.M.; Mornex, J.-F. SRLVs: A Genetic Continuum of Lentiviral Species in Sheep and Goats with Cumulative Evidence of Cross Species Transmission. Curr. HIV Res. 2010, 8, 94–100. [Google Scholar] [CrossRef] [PubMed]
  12. Da Cruz, J.C.M.; Singh, D.K.; Lamara, A.; Chebloune, Y. Small Ruminant Lentiviruses (SRLVs) Break the Species Barrier to Acquire New Host Range. Viruses 2013, 5, 1867–1884. [Google Scholar] [CrossRef] [PubMed]
  13. Michiels, R.; Adjadj, N.R.; De Regge, N. Phylogenetic Analysis of Belgian Small Ruminant Lentiviruses Supports Cross Species Virus Transmission and Identifies New Subtype B5 Strains. Pathogens 2020, 9, 183. [Google Scholar] [CrossRef]
  14. Gjerset, B.; Rimstad, E.; Teige, J.; Soetaert, K.; Jonassen, C.M. Impact of Natural Sheep-Goat Transmission on Detection and Control of Small Ruminant Lentivirus Group C Infections. Vet. Microbiol. 2009, 135, 231–238. [Google Scholar] [CrossRef] [PubMed]
  15. Gjerset, B.; Storset, A.K.; Rimstad, E. Genetic Diversity of Small-Ruminant Lentiviruses: Characterization of Norwegian Isolates of Caprine Arthritis Encephalitis Virus. J. Gen. Virol. 2006, 87, 573–580. [Google Scholar] [CrossRef]
  16. Rosati, S.; Mannelli, A.; Merlo, T.; Ponti, N. Characterization of the Immunodominant Cross-Reacting Epitope of Visna Maedi Virus and Caprine Arthritis-Encephalitis Virus Capsid Antigen. Virus Res. 1999, 61, 177–183. [Google Scholar] [CrossRef]
  17. Gomez-Lucia, E.; Barquero, N.; Domenech, A. Maedi-Visna Virus: Current Perspectives. Vet. Med. Res. Rep. 2018, 9, 11–21. [Google Scholar] [CrossRef]
  18. De Andrés, D.; Klein, D.; Watt, N.J.; Berriatua, E.; Torsteinsdottir, S.; Blacklaws, B.A.; Harkiss, G.D. Diagnostic Tests for Small Ruminant Lentiviruses. Vet. Microbiol. 2005, 107, 49–62. [Google Scholar] [CrossRef]
  19. Herrmann-Hoesing, L.M. Diagnostic Assays Used to Control Small Ruminant Lentiviruses. J. Vet. Diagn. Investig. 2010, 22, 843–855. [Google Scholar] [CrossRef]
  20. Karanikolaou, K.; Angelopoulou, K.; Papanastasopoulou, M.; Koumpati-Artopiou, M.; Papadopoulos, O.; Koptopoulos, G. Detection of Small Ruminant Lentiviruses by PCR and Serology Tests in Field Samples of Animals from Greece. Small Rumin. Res. 2005, 58, 181–187. [Google Scholar] [CrossRef]
  21. Brellou, G.D.; Angelopoulou, K.; Poutahidis, T.; Vlemmas, I. Detection of Maedi-Visna Virus in the Liver and Heart of Naturally Infected Sheep. J. Comp. Pathol. 2007, 136, 27–35. [Google Scholar] [CrossRef] [PubMed]
  22. Angelopoulou, K.; Brellou, G.D.; Vlemmas, I. Detection of Maedi-Visna Virus in the Kidneys of Naturally Infected Sheep. J. Comp. Pathol. 2006, 134, 329–335. [Google Scholar] [CrossRef] [PubMed]
  23. Kalogianni, A.I.; Bouzalas, I.; Bossis, I.; Gelasakis, A.I. Seroepidemiology of Maedi-Visna in Intensively Reared Dairy. Animals 2023, 13, 2273. [Google Scholar] [CrossRef] [PubMed]
  24. Angelopoulou, K.; Karanikolaou, K.; Papanastasopoulou, M.; Koumpati-Artopiou, M.; Vlemmas, I.; Papadopoulos, O.; Koptopoulos, G. First Partial Characterisation of Small Ruminant Lentiviruses from Greece. Vet. Microbiol. 2005, 109, 1–9. [Google Scholar] [CrossRef] [PubMed]
  25. Chassalevris, T.; Chaintoutis, S.C.; Apostolidi, E.D.; Giadinis, N.D.; Vlemmas, I.; Brellou, G.D.; Dovas, C.I. A Highly Sensitive Semi-Nested Real-Time PCR Utilizing Oligospermine-Conjugated Degenerate Primers for the Detection of Diverse Strains of Small Ruminant Lentiviruses. Mol. Cell. Probes 2020, 51, 101528. [Google Scholar] [CrossRef] [PubMed]
  26. Eltahir, Y.M.; Dovas, C.I.; Papanastassopoulou, M.; Koumbati, M.; Giadinis, N.; Verghese-Nikolakaki, S.; Koptopoulos, G. Development of a Semi-Nested PCR Using Degenerate Primers for the Generic Detection of Small Ruminant Lentivirus Proviral DNA. J. Virol. Methods 2006, 135, 240–246. [Google Scholar] [CrossRef] [PubMed]
  27. Kalogianni, A.I.; Bouzalas, I.; Bossis, I.; Gelasakis, A.I. A Longitudinal Cohort Study of Risk Factors Associated with Small Ruminant Lentivirus Seropositivity in Intensively Reared Dairy Ewes in Greece. Pathogens 2023, 12, 1200. [Google Scholar] [CrossRef]
  28. Tamura, K.; Stecher, G.; Kumar, S. MEGA11: Molecular Evolutionary Genetics Analysis Version 11. Mol. Biol. Evol. 2021, 38, 3022–3027. [Google Scholar] [CrossRef]
  29. Thompson, J.D.; Higgins, D.G.; Gibson, T.J. CLUSTAL W: Improving the Sensitivity of Progressive Multiple Sequence Alignment through Sequence Weighting, Position-Specific Gap Penalties and Weight Matrix Choice. Nucleic Acids Res. 1994, 22, 4673–4680. [Google Scholar] [CrossRef]
  30. Huelsenbeck, J.P.; Crandall, K.A. Phylogeny Estimation and Hypothesis Testing Using Maximum Likelihood. Annu. Rev. Ecol. Syst. 1997, 28, 437–466. [Google Scholar] [CrossRef]
  31. Felsenstein, J. Confidence Limits on Phylogenies: An Approach Using the Bootstrap. Evolution 1985, 39, 783. [Google Scholar] [CrossRef] [PubMed]
  32. Nei, M.; Kumar, S. Molecular Evolution and Phylogenetics; Oxford University Press: New York, NY, USA, 2000; ISBN 9780199881222. [Google Scholar]
  33. Olech, M.; Rachid, A.; Croisé, 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]
  34. Olech, M.; Hodor, D.; Toma, C.; Negoescu, A.; Taulescu, M. First Molecular Characterization of Small Ruminant Lentiviruses Detected in Romania. Animals 2023, 13, 3718. [Google Scholar] [CrossRef] [PubMed]
  35. 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]
  36. Ramírez, H.; Glaria, I.; de Andrés, X.; Martínez, H.; Hernández, M.; Reina, R.; Iráizoz, E.; Crespo, H.; Berriatua, E.; Vázquez, J.; et al. Recombinant Small Ruminant Lentivirus Subtype B1 in Goats and Sheep of Imported Breeds in Mexico. Vet. J. 2011, 190, 169–172. [Google Scholar] [CrossRef] [PubMed]
  37. Mendiola, W.P.S.; Tórtora, J.L.; Martínez, H.A.; García, M.M.; Cuevas-Romero, S.; Cerriteño, J.L.; Ramírez, H. Genotyping Based on the LTR Region of Small Ruminant Lentiviruses from Naturally Infected Sheep and Goats from Mexico. BioMed Res. Int. 2019, 2019, 4279573. [Google Scholar] [CrossRef]
  38. Martín-Sánchez, A.M.; Ciro-Gómez, G.; Sayas, E.; Vilella-Esplá, J.; Ben-Abda, J.; Pérez-Alvarez, J.Á. Date Palm By-Products as a New Ingredient for the Meat Industry: Application to Pork Liver Pâté. MESC 2013, 93, 880–887. [Google Scholar] [CrossRef] [PubMed]
  39. Olech, M.; Kuźmak, J. Molecular Characterization of Small Ruminant Lentiviruses in Polish Mixed Flocks Supports Evidence of Cross Species Transmission, Dual Infection, a Recombination Event, and Reveals the Existence of New Subtypes within Group A. Viruses 2021, 13, 2529. [Google Scholar] [CrossRef]
  40. Barquero, N.; Arjona, A.; Domenech, A.; Toural, C.; De Las Heras, A.; Fernández-Garayzabal, J.F.; Ruiz-Santa Quiteria, J.A.; Gomez-Lucia, E. Paper: Diagnostic Performance of PCR and ELISA on Blood and Milk Samples and Serological Survey for Small Ruminant Lentiviruses in Central Spain. Vet. Rec. 2011, 168, 20. [Google Scholar] [CrossRef]
  41. De Regge, N.; Cay, B. Development, Validation and Evaluation of Added Diagnostic Value of a q(RT)-PCR for the Detection of Genotype A Strains of Small Ruminant Lentiviruses. J. Virol. Methods 2013, 194, 250–257. [Google Scholar] [CrossRef]
  42. Barquero, N.; Domenech, A.; Arjona, A.; Fernández-Garayzabal, J.F.; Ruiz-Santa-Quiteria, J.A.; Gomez-Lucia, E. Comparison of Two PCR and One ELISA Techniques for the Detection of Small Ruminant Lentiviruses (SRLVs) in Milk of Sheep and Goats. Res. Vet. Sci. 2013, 94, 817–819. [Google Scholar] [CrossRef] [PubMed]
  43. Dolfini, T.; Conrad, L.F.; Flores, I.V.C.; Ravazzolo, A.P. Comparison of Primer Pairs: Greater Degeneracy Improves Small Ruminant Lentivirus (SRLVs) Detection by Seminested PCR. Small Rumin. Res. 2015, 123, 189–192. [Google Scholar] [CrossRef]
  44. Zanoni, R.G.; Nauta, I.M.; Kuhnert, P.; Pauli, U.; Pohl, B.; Peterhans, E. Genomic Heterogeneity of Small Ruminant Lentiviruses Detected by PCR. Vet. Microbiol. 1992, 33, 341–351. [Google Scholar] [CrossRef] [PubMed]
  45. 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]
  46. Grego, E.; Bertolotti, L.; Carrozza, M.L.; Profiti, M.; Mazzei, M.; Tolari, F.; Rosati, S. Genetic and Antigenic Characterization of the Matrix Protein of Two Genetically Distinct Ovine Lentiviruses. Vet. Microbiol. 2005, 106, 179–185. [Google Scholar] [CrossRef]
  47. 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]
  48. Germain, K.; Valas, S. Distribution and Heterogeneity of Small Ruminant Lentivirus Envelope Subtypes in Naturally Infected French Sheep. Virus Res. 2006, 120, 156–162. [Google Scholar] [CrossRef]
  49. Molaee, V.; Bazzucchi, M.; De Mia, G.M.; Otarod, V.; Abdollahi, D.; Rosati, S.; Lühken, G. Phylogenetic Analysis of Small Ruminant Lentiviruses in Germany and Iran Suggests Their Expansion with Domestic Sheep. Sci. Rep. 2020, 10, 2243. [Google Scholar] [CrossRef]
  50. Carrozza, M.L.; Niewiadomska, A.M.; Mazzei, M.; Abi-Said, M.R.; Hué, S.; Hughes, J.; Gatseva, A.; Gifford, R.J. Emergence and Pandemic Spread of Small Ruminant Lentiviruses. Virus Evol. 2023, 9, vead005. [Google Scholar] [CrossRef]
  51. Zhao, L.; Zhang, L.; Shi, X.; Duan, X.; Li, H.; Liu, S. Next-Generation Sequencing for the Genetic Characterization of Maedi/Visna Virus Isolated from the Northwest of China. J. Vet. Sci. 2021, 22, e66. [Google Scholar] [CrossRef]
  52. Laamanen, I.; Jakava-Viljanen, M.; Sihvonen, L. Genetic Characterization of Maedi-Visna Virus (MVV) Detected in Finland. Vet. Microbiol. 2007, 122, 357–365. [Google Scholar] [CrossRef]
  53. Grego, E.; Profiti, M.; Giammarioli, M.; Giannino, L.; Rutili, D.; Woodall, C.; Rosati, S. Genetic Heterogeneity of Small Ruminant Lentiviruses Involves Immunodominant Epitope of Capsid Antigen and Affects Sensitivity of Single-Strain-Based Immunoassay. Clin. Diagn. Lab. Immunol. 2002, 9, 828–832. [Google Scholar] [CrossRef]
  54. L’Homme, Y.; Ouardani, M.; Lévesque, V.; Bertoni, G.; Simard, C.; Pisoni, G. Molecular Characterization and Phylogenetic Analysis of Small Ruminant Lentiviruses Isolated from Canadian Sheep and Goats. Virol. J. 2011, 8, 271. [Google Scholar] [CrossRef] [PubMed]
  55. Lacerenza, D.; Giammarioli, M.; Grego, E.; Marini, C.; Profiti, M.; Rutili, D.; Rosati, S. Antibody Response in Sheep Experimentally Infected with Different Small Ruminant Lentivirus Genotypes. Vet. Immunol. Immunopathol. 2006, 112, 264–271. [Google Scholar] [CrossRef] [PubMed]
  56. Reina, R.; Grego, E.; Profiti, M.; Glaria, I.; Robino, P.; Quasso, A.; Amorena, B.; Rosati, S. Development of Specific Diagnostic Test for Small Ruminant Lentivirus Genotype E. Vet. Microbiol. 2009, 138, 251–257. [Google Scholar] [CrossRef] [PubMed]
  57. Olech, M. The Genetic Variability of Small-Ruminant Lentiviruses and Its Impact on Tropism, the Development of Diagnostic Tests and Vaccines and the Effectiveness of Control Programmes. J. Vet. Res. 2023, 67, 479–502. [Google Scholar] [CrossRef] [PubMed]
  58. Nogarol, C.; Bertolotti, L.; Klevar, S.; Profiti, M.; Gjerset, B.; Rosati, S. Serological Characterization of Small Ruminant Lentiviruses: A Complete Tool for Serotyping Lentivirus Infection in Goat. Small Rumin. Res. 2019, 176, 42–46. [Google Scholar] [CrossRef]
  59. Acevedo Jiménez, G.E.; Tórtora Pérez, J.L.; Rodríguez Murillo, C.; Arellano Reynoso, B.; Ramírez Álvarez, H. Serotyping versus Genotyping in Infected Sheep and Goats with Small Ruminant Lentiviruses. Vet. Microbiol. 2021, 252, 108931. [Google Scholar] [CrossRef] [PubMed]
  60. Arcangeli, C.; Torricelli, M.; Sebastiani, C.; Lucarelli, D.; Ciullo, M.; Passamonti, F.; Giammarioli, M.; Biagetti, M. Genetic Characterization of Small Ruminant Lentiviruses (SRLVs) Circulating in Naturally Infected Sheep in Central Italy. Viruses 2022, 14, 686. [Google Scholar] [CrossRef]
  61. Gayo, E.; Cuteri, V.; Preziuso, S. Genetic Characterization and Phylogenetic Analysis of Small Ruminant Lentiviruses Detected in Spanish Assaf Sheep with Different Mammary Lesions. Viruses 2018, 10, 315. [Google Scholar] [CrossRef]
Figure 1. Phylogenetic tree indicates the relationship of gag nucleotide sequences (~800 bp) of the present study with the available database (Genbank) SRLV strains originating from different geographical areas. This unrooted tree was inferred in MEGA11 by using the Maximum Likelihood method and General Time Reversible. Bootstrap values are based on 100 repetitions and are shown at the nodes. The tree with the highest log likelihood (−16884.54) is shown. The percentage of trees in which the associated taxa clustered together is shown above the branches. A discrete Gamma distribution was used to model evolutionary rate differences among sites [5 categories (+G, parameter = 0.6938)]. The rate variation model allowed for some sites to be evolutionarily invariable ([+I], 17.37% sites). This analysis involved 65 nucleotide sequences. Sequences from the present study are shown with a red rhombus symbol (GenBank: OR283217 to OR283228). Database derived sequences are denoted with their strain name, GenBank accession number and genotype in parenthesis, the abbreviation of the country (CH: China; BR: Brazil; CA: Canada; GR: Greece; IC: Iceland; IT: Italy; J: Jordan; LB: Lebanon; ME: Mexico; NW: Norway; PL: Poland; PR: Portugal; SA: South-Africa; SC: Scotland; SP: Spain; SW; Switzerland; TR: Turkey; USA: United States of America), and the host (s: sheep; g: goat).
Figure 1. Phylogenetic tree indicates the relationship of gag nucleotide sequences (~800 bp) of the present study with the available database (Genbank) SRLV strains originating from different geographical areas. This unrooted tree was inferred in MEGA11 by using the Maximum Likelihood method and General Time Reversible. Bootstrap values are based on 100 repetitions and are shown at the nodes. The tree with the highest log likelihood (−16884.54) is shown. The percentage of trees in which the associated taxa clustered together is shown above the branches. A discrete Gamma distribution was used to model evolutionary rate differences among sites [5 categories (+G, parameter = 0.6938)]. The rate variation model allowed for some sites to be evolutionarily invariable ([+I], 17.37% sites). This analysis involved 65 nucleotide sequences. Sequences from the present study are shown with a red rhombus symbol (GenBank: OR283217 to OR283228). Database derived sequences are denoted with their strain name, GenBank accession number and genotype in parenthesis, the abbreviation of the country (CH: China; BR: Brazil; CA: Canada; GR: Greece; IC: Iceland; IT: Italy; J: Jordan; LB: Lebanon; ME: Mexico; NW: Norway; PL: Poland; PR: Portugal; SA: South-Africa; SC: Scotland; SP: Spain; SW; Switzerland; TR: Turkey; USA: United States of America), and the host (s: sheep; g: goat).
Viruses 16 00547 g001
Figure 2. Amino acid sequence alignment of gag-pol sequences from the present study (OR283217–OR283228) and reference SRLV strains from GenBank belonging to genotypes A, B, C, and E. Immunodominant regions are boxed: epitope 2 and “GG” motif (a), Major Homology Region (b), and epitope 3 (c). Dots represent identical amino acid residues with reference strain K1514, dashes indicate deletions, and missing information sites are denoted by a question mark.
Figure 2. Amino acid sequence alignment of gag-pol sequences from the present study (OR283217–OR283228) and reference SRLV strains from GenBank belonging to genotypes A, B, C, and E. Immunodominant regions are boxed: epitope 2 and “GG” motif (a), Major Homology Region (b), and epitope 3 (c). Dots represent identical amino acid residues with reference strain K1514, dashes indicate deletions, and missing information sites are denoted by a question mark.
Viruses 16 00547 g002aViruses 16 00547 g002b
Table 1. Information of the PCR positive samples of the study.
Table 1. Information of the PCR positive samples of the study.
Sample ID BreedAge
(Years)
Serological StatusFarmLocationAccession NumberGenotype
s-16sxLacaune3.5+ACentral GreeceOR283228 B
s-40sxLacaune2.0AOR283225 B
s-56sxChios7.0+AOR283223 A
s-11sChios2.5+BEpirusOR283218 A
s-14sChios7.0+BOR283221 A
s-35sChios1.5BOR283220 A
s-83sChios7.0+BOR283224 A
s-10flChios4.0+CPeloponneseOR283217 A
s-21flChios3.0COR283219 A
s-42flChios3.0COR283222 A
s-8mLacaune3.0DWestern GreeceOR283226 B
s-98mLacaune3.0+DOR283227 B
+: ELISA positive animal; −: ELISA negative animal.
Table 2. Nucleotide (plain numbers) and amino acid (bold numbers) distance values (%) between partial gag-pol sequences obtained in this study.
Table 2. Nucleotide (plain numbers) and amino acid (bold numbers) distance values (%) between partial gag-pol sequences obtained in this study.
s-16sxs-40sxs-56sxs-11ss-14ss-35ss-83ss-10fls-21fls-42fls-8m
s-40sx3.9----------
2.2----------
s-56sx26.827.3---------
13.915.2---------
s-11s24.224.96.8--------
14.615.83.8--------
s-14s24.725.25.44.7-------
13.014.22.52.0-------
s-35s24.725.05.75.24.8------
13.815.03.43.22.0------
s-83s23.024.34.74.74.73.6-----
12.613.83.43.22.02.4-----
s-10fl25.426.16.55.14.14.75.1----
14.215.43.82.82.42.43.6----
s-21fl25.426.16.55.14.14.75.1----
14.215.43.82.82.42.43.6----
s-42fl28.929.211.911.610.511.211.411.111.1--
15.917.15.56.14.95.36.56.16.1--
s-8m4.95.625.223.024.423.922.624.324.326.5-
2.42.814.314.613.814.613.415.015.015.9-
s-98m5.66.325.623.224.423.822.324.424.426.41.2
2.42.814.314.313.514.313.514.714.715.90.0
Table 3. Nucleotide (plain numbers) and amino acid (bold numbers) distance values (%) between partial gag-pol sequences obtained in this study and representative SRLV strains.
Table 3. Nucleotide (plain numbers) and amino acid (bold numbers) distance values (%) between partial gag-pol sequences obtained in this study and representative SRLV strains.
s-16sxs-40sxs-56sxs-11ss-14ss-35ss-83ss-10fls-21fls-42fls-8ms-98m
EV1 (A1/SC)27.929.324.221.222.421.721.922.822.824.726.525.8
18.719.912.213.312.913.311.214.114.115.519.118.4
K1514 (A1/USA)23.324.822.019.321.219.619.820.520.523.824.223.8
14.215.48.08.46.47.66.88.88.88.514.614.7
SA-OMVV (A1/SA)24.125.721.319.619.819.120.620.020.021.523.923.9
18.319.58.411.610.011.610.412.012.012.218.718.0
KV1772 (A1/IC)23.425.022.219.521.319.820.020.720.723.624.324.0
13.815.07.68.06.07.26.48.48.48.114.214.3
OvLv 85/34 (A2/USA)20.522.719.018.718.517.717.218.218.220.120.921.4
15.416.67.68.46.47.66.88.88.89.315.014.7
SRLV007 (A8/IT)27.328.920.718.619.818.818.119.619.621.026.226.2
17.418.610.512.811.612.812.012.012.013.818.217.6
SRLV009 (A18/IT)26.829.622.720.421.120.420.921.321.323.726.324.8
19.020.211.814.413.214.013.614.814.815.919.418.8
SRLV038 (A19/IT)25.026.920.118.418.717.719.117.617.619.524.524.4
17.418.67.110.09.210.09.610.010.011.017.817.1
It-561 (A19/IT)27.427.520.719.820.220.820.619.519.520.324.925.3
18.219.49.711.510.711.511.111.511.511.519.019.0
sh248 (A22/LB)28.328.724.223.521.921.522.421.721.724.027.428.5
15.817.011.311.610.011.610.011.211.212.617.017.1
sh173 (A22/J)30.030.426.224.823.422.824.923.723.725.128.628.1
15.817.010.510.89.611.211.211.211.211.016.616.3
SRLV-Greece-S1 (A/GR)22.424.212.610.011.510.09.610.410.411.123.122.1
10.612.22.21.70.62.21.12.82.83.311.111.1
SRLV-Greece-S2 (A/GR)22.925.512.610.011.010.49.110.410.411.323.422.4
11.713.33.93.32.22.82.83.93.93.912.212.2
CAEV-Co (B1/USA)15.115.027.225.223.926.124.824.924.926.013.213.1
8.49.214.315.313.714.114.114.514.515.57.26.5
GANSU (B1/CH)15.115.426.924.123.724.424.223.923.927.914.114.5
10.010.817.316.915.314.515.716.116.117.69.69.3
Ov496 (B2/SP)9.99.626.824.224.624.924.225.225.228.58.79.5
2.83.215.216.114.515.314.115.715.716.72.02.0
SRLV042 (B2/IT)7.48.528.124.925.725.424.525.925.928.36.26.9
10.110.517.419.419.020.218.519.819.821.79.38.5
It-007 (B2/IT)7.78.628.125.225.125.424.626.426.429.46.47.2
10.410.816.919.318.519.718.119.319.321.29.68.9
It-Pi1 (B2/IT)8.38.927.526.025.826.125.226.126.129.36.37.1
15.716.120.322.922.123.321.722.922.924.914.914.2
SRLV001 (B2/IT)8.88.627.625.026.526.625.327.327.329.26.97.2
10.811.217.719.719.320.518.920.120.121.210.09.3
Volterra (B3/IT)23.926.226.426.527.726.326.726.226.228.822.823.3
14.615.417.920.319.520.719.120.320.321.514.613.9
TR-DM (B3/TR)25.225.727.928.627.426.928.528.528.531.623.923.9
14.114.916.518.517.718.918.118.518.520.413.713.0
1GA (C/NW)25.925.732.032.230.431.330.431.531.531.623.524.2
20.621.923.325.023.823.823.824.224.225.420.620.0
Roccaverano (E1/IT)35.537.841.439.839.938.737.938.438.440.136.437.3
27.228.526.728.627.828.227.428.228.228.327.627.0
Seui (E2/IT)38.340.440.837.436.737.136.136.436.440.338.639.3
27.228.527.529.428.229.427.829.029.029.928.027.5
SC: Scotland; USA: United States of America; SA: South-Africa; IC: Iceland; IT: Italy; LB: Lebanon; J: Jordan; GR: Greece; CH: China; SP: Spain; TR: Turkey; NW: Norway.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Kalogianni, A.I.; Bouzalas, I.; Marka, S.; Zografaki, M.-E.; Mavrikou, S.; Gelasakis, A.I. Genetic Characterization of Small Ruminant Lentiviruses Isolated from Dairy Sheep in Greece. Viruses 2024, 16, 547. https://doi.org/10.3390/v16040547

AMA Style

Kalogianni AI, Bouzalas I, Marka S, Zografaki M-E, Mavrikou S, Gelasakis AI. Genetic Characterization of Small Ruminant Lentiviruses Isolated from Dairy Sheep in Greece. Viruses. 2024; 16(4):547. https://doi.org/10.3390/v16040547

Chicago/Turabian Style

Kalogianni, Aphrodite I., Ilias Bouzalas, Sofia Marka, Maria-Eleftheria Zografaki, Sofia Mavrikou, and Athanasios I. Gelasakis. 2024. "Genetic Characterization of Small Ruminant Lentiviruses Isolated from Dairy Sheep in Greece" Viruses 16, no. 4: 547. https://doi.org/10.3390/v16040547

APA Style

Kalogianni, A. I., Bouzalas, I., Marka, S., Zografaki, M. -E., Mavrikou, S., & Gelasakis, A. I. (2024). Genetic Characterization of Small Ruminant Lentiviruses Isolated from Dairy Sheep in Greece. Viruses, 16(4), 547. https://doi.org/10.3390/v16040547

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