Next Article in Journal
Phage Therapy Potentiates Second-Line Antibiotic Treatment against Pneumonic Plague
Previous Article in Journal
Effective Reduction of SARS-CoV-2 RNA Levels Using a Tailor-Made Oligonucleotide-Based RNA Inhibitor
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Genetic Characterization of Small Ruminant Lentiviruses (SRLVs) Circulating in Naturally Infected Sheep in Central Italy

by
Chiara Arcangeli
1,†,
Martina Torricelli
1,†,
Carla Sebastiani
1,*,
Daniele Lucarelli
2,
Marcella Ciullo
1,
Fabrizio Passamonti
3,
Monica Giammarioli
1 and
Massimo Biagetti
1
1
Istituto Zooprofilattico Sperimentale dell’Umbria e delle Marche “Togo Rosati”, 06126 Perugia, Italy
2
TranslaTUM Center for Translational Cancer Research, Technical University of Munich, 81675 Munich, Germany
3
Dipartimento di Medicina Veterinaria, Università degli Studi di Perugia, 06126 Perugia, Italy
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Viruses 2022, 14(4), 686; https://doi.org/10.3390/v14040686
Submission received: 28 February 2022 / Revised: 22 March 2022 / Accepted: 24 March 2022 / Published: 25 March 2022
(This article belongs to the Topic Veterinary Infectious Diseases)

Abstract

:
Small ruminant lentiviruses (SRLVs) represent a very heterogeneous group of ss-RNA viruses that infect sheep and goats worldwide. They cause important, deleterious effects on animal production and limit the animal trade. SRLVs show a high genetic variability due to high mutation rate and frequent recombination events. Indeed, five genotypes (A–E) and several subtypes have been detected. The aim of this work was to genetically characterize SRLVs circulating in central Italy. On this basis, a phylogenetic study on the gag-pol genetic region of 133 sheep, collected from 19 naturally infected flocks, was conducted. In addition, to evaluate the frequency of mutation and the selective pressure on this region, a WebLogo 3 analysis was performed, and the dN/dS ratio was computed. The results showed that 26 samples out of 133 were clustered in genotype A and 106 samples belonged to genotype B, as follows: A9 (n = 8), A11 (n = 10), A24 (n = 7), B1 (n = 2), B2 (n = 59), and B3 (n = 45). No recombination events were found. Mutations were localized mainly in the VR-2 region, and the dN/dS ratio of 0.028 indicated the existence of purifying selection. Since the genetic diversity of SRLVs could make serological identification difficult, it is important to perform molecular characterization to ensure a more reliable diagnosis, to maintain flock health status, and for the application of local and national control programs.

1. Introduction

Small ruminant lentiviruses (SRLVs), which include Visna-maedi virus (VMV) and caprine arthritis encephalitis virus (CAEV), belong to the Retroviridae family and Lentivirus genus [1]. VMV and CAEV, initially considered specific to sheep and goats, respectively, indiscriminately infected both species, representing a “genetic continuum” due to the genetic closeness of the two viruses [2,3].
SRLVs stably infect the monocyte/macrophage lineage, modulating cellular responses, differentiation pathways, and cytokine secretion in order to ensure prolonged viral replication [4].
VMV and CAEV are not associated with immunodeficiency differently from other lentiviruses, but they cause a multisystem and chronic disease that affects many organs, such as the lungs, mammary glands, joints, and the nervous system [5,6]. The clinical patterns of infection may include interstitial pneumonia, dyspnoea, indurative mastitis (hard udder), arthritis, encephalitis, lymphadenopathy, and chronic weight loss [7]. Historically, the lactogenic route was considered the primary mode of virus transmission; however, housing type, animal density, and the length of time healthy and infected animals are raised in close contact, all greatly promote the spread of SRLVs. The intrauterine and seminal routes could represent potential sources of viral diffusion, which, however, requires further confirmatory investigations [8].
These viruses cause huge economic losses, affecting both production and animal welfare. Epizootic data suggest an increased seroprevalence, especially in Europe (40.9%), where sheep and goat farming are of great significance [2,9]. The profit losses are represented mostly by reductions in birth and growth weight, milk production, as well as by premature culling and trade restrictions [3]. Furthermore, the subclinical course of the disease makes the identification of infected animals difficult. Thus, early, and accurate diagnosis of SRLVs is required in order to adopt adequate control programs to reduce the prevalence of infection [9] and control the disease. The World Organisation for Animal Health (OIE) has included SRLVs in the list of notifiable terrestrial and aquatic animal diseases and has recognized the enzyme-linked immunosorbent assay (ELISA) as the method of choice for SRLV diagnosis [10]. Anyway, SRLVs’ mutation rate and antigenic heterogeneity determine the occurrence of new genetic subtypes that can escape diagnosis using current monovalent serological tests. Therefore, PCR combined with ELISA might allow a more reliable diagnosis, especially when used as a confirmatory test to detect the infection before seroconversion and for viral characterization based on specific genetic regions [7].
Similar to other retroviruses, the SRLV genome is characterised by two identical positive-sense single-stranded RNA subunits (8.4–9.2 kb) [11], consisting of three structural genes (gag, pol, and env) as well as three accessory genes (vif, vpr-like, and rev) flanked by non-coding long terminal repeat sequences (LTRs). The last ones provide the signals required for transcription, integration, and polyadenylation of viral RNA [12]. The gag gene encodes for three structural proteins, the matrix (p16MA), the capsid (p25CA) and the nucleocapsid (p14NC) proteins; the pol gene codes for protease (PR), retrotranscriptase (RT), dUTPase (DU) and integrase (IN) enzymes; the env gene encodes for the transmembrane (gp46TM) and surface (gp135SU) glycoproteins [13,14,15]. Among them, gag and pol genes exhibit less genetic variability than the env gene [8,16]. The high genetic heterogeneity originates from reverse transcriptase low fidelity, from a lack of 3′ exonuclease proofreading capability, and from the absence of repair mechanisms during replication [17].
Cross-species infections, genetic drift, genetic recombination in host cells during co-infections, and the formation of quasi-species can determine the introduction of new subtypes [18], leading to different pathogenicity, disease progression, and susceptibility to bacterial infections [3,19,20]. In particular, recombination represents a huge evolutionary advantage because, by eliminating deleterious mutations and assembling beneficial fragments, it ensures adaptability, integration, and persistence of the virus in the host [21]. All these aspects contribute to making a vaccine or specific treatment development a major challenge.
Likewise, genetic variability is the key feature of the SRLVs, so its deep knowledge allows to perform a more correct and accurate serological and molecular diagnosis, in addition to a better comprehension of host-virus interactions.
With regard to phylogenetic classification, based on gag and pol genes, Shah [22] classified SRLVs into five genetic groups (A–E) and into several subtypes.
Indeed, in accordance with human immunodeficiency virus (HIV) classification criteria [23], differences in the nucleotide sequences range from 25% to 37% and from 15% to 27% among genotypes and among subtypes, respectively.
Genotype A consists of MVV-like strains and represents the most common and heterogeneous group, containing 24 distinct subtypes (A1–A24) [24,25,26].
Genotype B, represented by the subtypes B1–B3, includes CAEV-like strains and exhibits lower genetic variability compared to genotype A [26,27]. It has to be noted that subtype B4, previously identified by Santry [3], was later reclassified as a recombinant strain [28]. Subtype B5 was identified by analysing only the pol region, but it belongs to subtype B1 considering the gag-pol overlapping region [25].
On the other hand, genotypes C, D, and E (subtypes E1–E2) are restricted to specific geographical areas. Genotype C has been identified only in Norwegian sheep and goats [29], genotype D was found only in Swiss and Spanish sheep and genotype E is closely associated with goats of northern Italy (strain E1), Sardinia, and Umbria regions (strain E2) [30,31,32]. Genotype D has been initially identified by analysing only the pol gene [22], but a subsequent additional phylogenetic analysis including also the gag region, on the same isolates, allowed to assign these strains to genotype A [11,26].
Many epidemiological and phylogenetic investigations conducted in different goat and sheep populations worldwide have shown that the following subtypes have a host preference: subtypes A2, A15, and A16 have been found only in sheep, while subtypes A7, A8, A10, A14, A17, E1, and E2 are present only in goats [25].
Regarding Italy, SRLVs are widespread, in particular in certain regions where sheep and goat breeding are of considerable importance [33,34]. Therefore, viral characterization is of great significance to allow the correct identification of flocks. A recent genetic and epidemiological study identified at least three genotypes (A, B, and E) and 14 subtypes (A3, A5, A8, A9, A11, A19, A20, A23, A24, B1, B2, B3, E1, and E2) circulating in Italy [24].
The aim of this work was to carry out a genetic characterization analysis of SRLVs circulating in central Italy, evaluating mutation frequencies and selection pressure on the SRLV gag region. This phylogenetic analysis was useful to characterize viral strains in a geographical area where sheep breeding represents a considerable economic resource.

2. Materials and Methods

2.1. SRLV Samples Collection and DNA Extraction

For the genetic characterization of SRLV strains circulating in the studied population, 206 sheep that resulted SRLV seropositive in a previous investigation [35] were analysed by gag-pol PCR. These samples were collected in the period 2019–2020 from 19 naturally infected sheep flocks located in central Italy. Genomic DNA was extracted from blood clots as described by Arcangeli [35].
Aliquots of the blood samples used for the analysis were taken during obligatory routine animal sanitary controls by authorised veterinarians; therefore, no ethical approval was required.

2.2. SRLV Proviral Amplification and Bioinformatic Analysis

The phylogenetic study was based on the partial sequence of the gag-pol genetic region as previously described by Shah [22]. A gag-pol fragment of 800 bp was amplified using the nested-PCR protocol described by Grego [31]. PCR products were detected by 1% agarose gel electrophoresis, purified from gel bands using the QIAquick Gel Extraction Kit (Qiagen, Hilden, Germany) according to manufacturer’s instructions, and then used as templates in sequencing reactions with BrilliantDyeTM Terminator Cycle Sequencing Kit v3.1 (NimaGen BV, Nijmegen, The Netherlands). For each sample, sequencing reactions were performed in three replicates in both sense and antisense strands and subsequently run in an ABI 3500 Genetic Analyzer (Applied Biosystem, Foster City, CA, USA).
Sequences dataset was analysed and nucleotide sequences were aligned to 154 published SRLV reference strains retrieved from GenBank (http://www.ncbi.nlm.nih.gov/genbank/, accessed on 10 January 2022) using BioEdit v7.2.5 software by ClustalW algorithm [36].
The phylogenetic tree was inferred using maximum likelihood (ML) method and Bayesian (BI) inference to improve robustness of the analysis. ML analysis were performed in MEGAX [37], using the General Time Reversible (GTR) statistical model [38] with gamma distribution + I (G + I) [39]. Cluster robustness was evaluated by performing 10,000 bootstrap replications, and branches with bootstrap values greater than 70% were clustered. BI analysis was evaluated with two runs consisting of Markov chains using BEAST v.1.8.4 with the GTR + G + I substitution model, selected as the best-fit nucleotide substitution model. A consensus tree was created using TreeAnnotator v.1.8.4 and the trees were displayed and edited using FigTree v.1.4.0 [40]. MEGAX was used to calculate pairwise distances between samples and reference strain sequences applying the p-distance model [37].
SRLV sequences were analysed by SplitsTree4 [41] software applying the Phi test of SplitsTree v.4 to assess the presence of recombination events.
The sequences detected in this study were submitted to the GenBank database and are available under the accession numbers from OK325451 to OK325583. All data generated and analysed are reported in Supplementary Table S1.

2.3. Weblogo Diagrams

In order to determine the amino acid profile and sequence variability of SRLVs, multiple alignment of the partial gag protein sequences of the analysed samples was performed by WebLogo v.3 software (http://weblogo.threeplusone.com/, accessed on 10 January 2022). Particularly, this analysis was carried out on the protein sequences deduced from the gag genetic fragments (nucleotide: 601–1236; numbering according to reference strain ItPi1—AY265456) evaluated in the phylogenetic investigation. The derived graphical representation consists of stacks of symbols, each one representing an amino acid of the protein chain, where the height of the symbols reflects the relative frequency of the amino acid residue at that specific position. In the graphical output, only one letter is shown when the residue is invariable, whereas different letters, corresponding to the most common amino acid substitutions, are present when the residue is variable [42,43].

2.4. Estimation of Nonsynonymous and Synonymous Substitution Rates

In order to estimate the selection pressure on SRLV gag gene, the non-synonymous (dN) and synonymous (dS) substitution ratio was calculated using SNAP (Synonymous Non-synonymous Analysis Program) software v2.1.1 implemented in the Los Alamos National Laboratory HIV-sequence database, Los Alamos, NM, USA (https://www.hiv.lanl.gov/content/sequence/SNAP/SNAP.html, accessed on 10 January 2022) [44,45]. The result derived from the dN/dS ratio (ω) allowed us to evaluate the following type of viral selection that occurred in the population: ω < 1 indicates a purifying selection, ω > 1 a positive selection while ω = 1 neutrality.

3. Results

3.1. Phylogenetic Analysis

A total of 133 SRLV partial gag-pol sequences were aligned with each other and with reference strains of genotypes A–E. Phylogenetic analysis was performed using the ML and the BI methods and resulted in the same classification of strains (data not shown).
Our results revealed that 26 samples belonged to genotype A and 106 samples belonged to genotype B (Figure 1). Specifically, we detected the following subtypes: A9 (n = 8), A11 (n = 10), A24 (n = 7), B1 (n = 2), B2 (n = 59), and B3 (n = 45). Interestingly, only a sample belonging to genotype A did not cluster into any of the already known subtypes, showing a mean genetic distance ranging from 0.129–0.237 among subtypes and thus labelled as unassigned.
Sequences belonging to both genotypes and specifically to their subtypes showed high nucleotide and amino acid sequence identities with the reference strains (Table S2). Samples closely related to subtypes A9, A11, and A24 showed percentage ranges of nucleotide identity with reference strains of 87.6–94.6%, 87.6–92.3%, and 82.5–98.4%, respectively. For B1, B2, and B3 subtypes, ranges of 82.8–90.5%, 83.8–95.2%, and 80.5–92.6% were observed, respectively.
Furthermore, the same samples were compared to their reference strains in order to evaluate the amino acid identity. We found percentages ≥ 90.0% for the subtypes belonging to genotype A and ≥ 83.1% for the subtypes belonging to genotype B.
In addition, nucleotide and amino acid sequence identity among strains of the same subtype with each other resulted also significant. In particular, the identity percentage ranged from 89.0% to 98.8% for the A9 subtype, from 87.4% to 98.4% for the A11 subtype, from 87.6% to 98.8% for the A24 subtype, from 84.1% to 99.8% for the B2 subtype, and from 79.8% to 99.6% for the B3 subtype. Finally, the two sequences clustered within the B1 subtype presented 98.8% of their identity.
Similarly, regarding the amino acid sequence identity among the samples with each other, the percentage was ≥93.8% for the subtypes belonging to genotype A and ≥85.3% for the subtypes belonging to genotype B.
These results were confirmed by the pairwise nucleotide distance comparison (Table S3). No sequences clustered within the C–E genotype strains and no recombination events were detected.

3.2. Comparative Analysis of Immunodominant Regions

In order to evaluate the level of conservation of the immunodominant regions, the deduced amino acid sequences of the gag fragment were aligned with the corresponding sequences of reference strains belonging to known subtypes of genotypes A and B (Figure 2). In particular, immunodominant regions in SRLV sequences include epitopes 2 and 3, major homology region (MHR), double glycine (GG) motif, and variability region-2 (VR-2).
Despite the moderate nucleotide heterogeneity, the amino acid sequences were quite conserved, and most of the found nucleotide mutations were synonymous.
Samples clustering in genotype A showed the presence of an asparagine-valine (NV) motif, and likewise, sequences clustering in genotype B showed the presence of the typical GG-motif [26].
Within genotype A, the epitopes 2 and 3 of the references and samples presented a high grade of conservation. In epitope 2, a lysine (K) was replaced by an arginine (R) only in sample 27358.13_2019, while in epitope 3 a threonine (T) was replaced by a serine (S) in samples 24109.15_2019, 24109.18_2019, 27358.71_2019 and a glutamic acid (E) was replaced by an aspartic acid (D) in the sample 24109.12_2019 (Figure 2, Table S1).
With regard to genotype B, more alterations were found in epitopes 2 and 3, where 11 and 25 samples showed at least one substitution, respectively. In the MHR region, which is usually highly conserved in retroviruses, the subtypes A9, A11, A24, and B2 were moderately variable, while in the remaining subtypes B1 and B3, the variability was more evident (Figure 2). Concerning the VR-2 region, the higher variability was observed for the A9 and B3 subtypes.
From Figure 2, it is possible to observe that sheep belonging to the same flock are characterised by the same subtype, except for flocks 1 and 10, which present different genotypes and subtypes, respectively. This could probably depend on trade activities, but these supporting data are not available.
Furthermore, the relative frequency of the amino acid substitutions can be deduced from the height of the single letter corresponding to the specific amino acid (Figure 3).
In order to determine the selective pressure on the SRLV gag gene, the ratio between synonymous and non-synonymous substitutions (dN/dS) was calculated. The value of 0.028 obtained from the analysis suggested the existence of a purifying selection (dN/dS < 1).

4. Discussion

SRLVs are widely spread in sheep and goats worldwide and are endemic in most European countries, causing multisystem diseases that affect animal production and welfare [46]. No particular attention is yet paid to SRLV infection and only sporadic control programs are applied in countries with high breeding density [47,48,49,50].
SRLVs are also widespread in the Italian sheep/goat population as well, and several subtypes have been identified and molecularly characterised by now. In a recent report [24], authors identified the following subtypes mainly circulating in Italy: A3, A5, A8, A9, A11, A19, A20, A23, A24, B1, B2, B3, E1, and E2. Some subtypes (i.e., A3 and A5), initially found in other European countries (Spain, Switzerland, Germany, Poland, Slovenia, and Turkey), have also been detected in the Italian sheep population, suggesting that animal movements represent a major source of viral circulation.
In Italy, there are no mandatory eradication programmes yet, but some northern and southern regions have implemented voluntary programmes to receive the SRLV-free status due to their high sheep and goat breeding activity [33,34,51].
In this paper, SRLV strains detected in 19 naturally infected sheep flocks in central Italy were characterized. The presence of some of the most common circulating subtypes in Italy, in particular the A9, A11, A24, B1, B2 and B3 strains, was confirmed. The new A24 subtype identified in sheep in Umbria, Lazio, and Marche regions by Bazzucchi [24] has also been found in seven animals analysed in this study. This evidence suggests that the A24 strain is currently circulating in the regions of central Italy, probably due to the intense commercial activity among these regions. In addition, the A11 subtype, which represents the most frequent genotype A strain, followed by the A9 subtype, was found. However, genotype B was the most represented in the studied population and specifically the B1, B2 and B3 subtypes were detected. A sample belonging to genotype A showed moderate genetic variability compared to the other samples of the same genotype. Although it resulted genetically close to subtype A24, the genetic distance value did not allow the clustering to any of the known subtypes. This finding could suggest the existence of a new subtype, but further investigations are necessary to confirm this hypothesis. Therefore, the sample was classified as unassigned.
From the results obtained in this study, it is evident that amino acid sequences are quite conserved, and the nucleotide mutations present in the samples are predominantly synonymous. In addition, higher levels of homology were found in sequences derived from animals from the same farms or coming from different farms but with a history of commercial trades and activities.
In order to estimate the selection pressure on SRLV gag gene sequences object of this study, the dN/dS ratio was calculated. The obtained value indicated the existence of purifying selection, that is, selection against new deleterious mutations. This phenomenon is essential for the virus to preserve its biological functions [52].
In our sheep population originating from non-mixed herds, the B1 subtype did not represent the most frequent one, while subtypes B2 and B3 were predominantly found, also according to Bazzucchi [24]. These authors observed a preferential tropism of the B1 strain for goats differently from the B2 and B3, which are known to be predominant in sheep [24].
Analysis of the deduced amino acid sequences revealed a high degree of gag protein conservation, especially of epitopes 2 and 3 (Figure 2). Conservation of these immunodominant regions is important for maintaining cross-reactivity in serological tests based on gag antigen [53].
However, MHR, which is usually quite conserved in all retroviruses [19,54], in our samples exhibited a certain grade of variability in both genotypes (Figure 2). HIV studies conducted in this region revealed that some non-synonymous substitutions impair the capsid assembly, reducing the infectivity of HIV-1 [55]. On the other hand, it is not clear whether the several substitutions detected in the MHR region in this study could exert the same effects. Little is known about mutations in this region that may lead to reduced infectivity; therefore, further studies are still necessary to support this hypothesis.
Insertions at amino acid positions 172–173 described by Molaee [26] were not detected in this study, confirming that this profile is specific to the A21 and A22 subtypes, found in German and Iranian flocks, and to other species of lentiviruses (i.e., BIV, SIV, FIV, EIAV, and HIV) [26].
The great variability of the subtypes circulating in Italy could depend on many conditions, including different managerial practises of breeding, environmental factors, and the presence in some territories of large wild ruminant populations [24,56]. The slow progression typical of lentiviruses causes an underestimation of infection prevalence, promoting virus spread and disease development. Since VM and CAE result in significant economic damage as a result of trade restrictions and impaired milk and meat production, it is essential to contain the infection [46].
For this purpose, several approaches could be applied, among which the most important is certainly serological screening for the identification of infected individuals. To this end, it would be appropriate to sensitise farmers, who still underestimate this serious problem, to make a proper assessment of their barns in terms of production data and analysis of the animal health status. Therefore, serological screening together with a careful evaluation of the flock may allow increased protection against infection.
In addition, for the application of local and national control programs, molecular characterization of viral isolates represents another tool essential for the development of more effective diagnostic tests.

5. Conclusions

The SRLVs are among the main causative agents that negatively affect the production and welfare of sheep and goats, and consequently, the profitability of affected farms. Therefore, accurate screening surveys for the early detection of infected animals together with more strict farm management are essential in order to contain the spread of infection. Furthermore, knowledge of SRLV genetic variability and of circulating genotypes and subtypes is important for epidemiological studies, for monitoring the effectiveness of control programs, for the development of new and more performant diagnostic assays, and in order to better clarify the evolution of these diseases.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/v14040686/s1, Table S1: SRLV sequences dataset analysed in this study; Table S2: (a) Nucleotide identity matrix, (b) Amino acid identity matrix; Table S3: Pairwise distances. References [57,58,59,60,61,62,63,64,65,66,67] are cited in the supplementary materials.

Author Contributions

Conceptualization, M.B.; formal analysis, C.A., D.L. and M.G.; investigation, C.A., M.T., C.S., D.L. and M.C.; resources, M.B.; data curation, C.A. and C.S.; writing—original draft preparation, C.A. and M.B.; writing—review and editing, C.A., M.T., C.S., F.P. and M.B.; supervision, M.B.; project administration, M.B. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Aliquots of the blood samples used for the analysis were taken during obligatory routine animal sanitary controls by authorised veterinarians; therefore, no ethical approval was required.

Informed Consent Statement

Not applicable.

Data Availability Statement

Accession numbers: the partial SRLV sequences generated in this study have been deposited in the NCBI GenBank database www.ncbi.nlm.nih.gov/genbank/, accessed on 28 September 2021 (Supplementary Table S1).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Krupovic, M.; Blomberg, J.; Coffin, J.M.; Dasgupta, I.; Fan, H.; Geering, A.D.; Gifford, R.; Harrach, B.; Hull, R.; Johnson, W.; et al. Ortervirales: New Virus Order Unifying Five Families of Reverse-Transcribing Viruses. J. Virol. 2018, 92, e00515-18. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. 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] [PubMed] [Green Version]
  3. 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]
  4. De Pablo-Maiso, L.; Echeverría, I.; Rius-Rocabert, S.; Luján, L.; Garcin, D.; de Andrés, D.; Nistal-Villán, E.; Reina, R. Sendai Virus, a Strong Inducer of Anti-Lentiviral State in Ovine Cells. Vaccines 2020, 8, 206. [Google Scholar] [CrossRef] [PubMed]
  5. Blacklaws, B.A. Small ruminant lentiviruses: Immunopathogenesis of visna-maedi and caprine arthritis and encephalitis virus. Comp. Immunol. Microbiol. Infect. Dis. 2012, 35, 259–269. [Google Scholar] [CrossRef]
  6. 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]
  7. 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. 2020, 252, 108931. [Google Scholar] [CrossRef]
  8. Furtado Araújo, J.; Andrioli, A.; Pinheiro, R.R.; Sider, L.H.; de Sousa, A.L.M.; de Azevedo, D.A.A.; Peixoto, R.M.; Lima, A.M.C.; Damasceno, E.M.; Souza, S.C.R.; et al. Vertical transmissibility of small ruminant lentivirus. PLoS ONE 2020, 15, e0239916. [Google Scholar] [CrossRef]
  9. De Miguel, R.; Arrieta, M.; Rodríguez-Largo, A.; Echeverría, I.; Resendiz, R.; Pérez, E.; Ruiz, H.; Pérez, M.; de Andrés, D.; Reina, R.; et al. Worldwide Prevalence of Small Ruminant Lentiviruses in Sheep: A Systematic Review and Meta-Analysis. Animals 2021, 11, 784. [Google Scholar] [CrossRef]
  10. Nowicka, D.; Czopowicz, M.; Mickiewicz, M.; Szaluś-Jordanow, O.; Witkowski, L.; Bagnicka, E.; Kaba, J. Diagnostic performance of ID screen MVV-CAEV Indirect Screening ELISA in identifying small ruminant lentiviruses-infected goats. Pol. J. Vet. Sci. 2014, 17, 501–506. [Google Scholar] [CrossRef]
  11. Ramírez, H.; Reina, R.; Amorena, B.; Andrés, D.D.; Martínez, H.A. Small Ruminant Lentiviruses: Genetic Variability, Tropism and Diagnosis. Viruses 2013, 5, 1175–1207. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. 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]
  13. Gomez-Lucia, E.; Barquero, N.; Domenech, A. Maedi-Visna virus: Current perspectives. Vet. Med. Res. Rep. 2018, 9, 11–21. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Querat, G.; Audoly, G.; Sonigo, P.; Vigne, R. Nucleotide sequence analysis of SA-OMVV, a visna-related ovine lentivirus: Phylogenetic history of lentiviruses. Virology 1990, 175, 434–447. [Google Scholar] [CrossRef]
  15. Vigne, R.; Filippi, P.; Quérat, G.; Sauze, N.; Vitu, C.; Russo, P.; Delori, P. Precursor Polypeptides to Structural Proteins of Visna Virus. J. Virol. 1982, 42, 1046–1056. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Olech, M.; Valas, S.; Kuźmak, J. Epidemiological survey in single-species flocks from Poland reveals expanded genetic and antigenic diversity of small ruminant lentiviruses. PLoS ONE 2018, 13, e0193892. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. Wainhobson, S. Running the gamut of retroviral variation. Trends Microbiol. 1996, 4, 135–141. [Google Scholar] [CrossRef]
  18. Olech, M.; Kuźmak, J. Molecular Characterization of Small Ruminant Lentiviruses of Subtype A5 Detected in Naturally Infected but Clinically Healthy Goats of Carpathian Breed. Pathogens 2020, 9, 992. [Google Scholar] [CrossRef]
  19. Gayo, E.; Cuteri, V.; Polledo, L.; Rossi, G.; García Marín, J.; 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] [Green Version]
  20. Urbanska, D.; Puchała, R.; Jarczak, J.; Czopowicz, M.; Kaba, J.; Horbańczuk, K.; Bagnicka, E. Does Small Ruminant Lentivirus Infection in Goats Predispose to Bacterial Infection of the Mammary Gland? A Preliminary Study. Animals 2021, 11, 1851. [Google Scholar] [CrossRef]
  21. Murphy, F.A. Veterinary Virology; Academic Press: San Diego, CA, USA, 1999; ISBN 9780080552033. [Google Scholar]
  22. 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] [Green Version]
  23. Peeters, M.; Esu-Williams, E.; Vergne, L.; Montavon, C.; Mulanga-Kabeya, C.; Harry, T.; Ibironke, A.; Lesage, D.; Patrel, D.; Delaporte, E. Predominance of subtype A and G HIV type 1 in Nigeria, with geographical differences in their distribution. AIDS Res. Hum. Retrovir. 2000, 16, 315–325. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Bazzucchi, M.; Pierini, I.; Gobbi, P.; Pirani, S.; Torresi, C.; Iscaro, C.; Feliziani, F.; Giammarioli, M. Genomic Epidemiology and Heterogeneity of SRLV in Italy from 1998 to 2019. Viruses 2021, 13, 2338. [Google Scholar] [CrossRef] [PubMed]
  25. 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] [Green Version]
  26. 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] [Green Version]
  27. Bertolotti, L.; Mazzei, M.; Puggioni, G.; Carrozza, M.L.; Dei Giudici, S.; Muz, D.; Juganaru, M.; Patta, C.; Tolari, F.; Rosati, S. Characterization of new small ruminant lentivirus subtype B3 suggests animal trade within the Mediterranean Basin. J. Gen. Virol. 2011, 92, 1923–1929. [Google Scholar] [CrossRef]
  28. 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] [Green Version]
  29. 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]
  30. 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]
  31. 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]
  32. Reina, R.; Bertolotti, L.; Dei Giudici, S.; Puggioni, G.; Ponti, N.; Profiti, M.; Patta, C.; Rosati, S. Small ruminant lentivirus genotype E is widespread in Sarda goat. Vet. Microbiol. 2010, 144, 24–31. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Nardelli, S.; Bettini, A.; Capello, K.; Bertoni, G.; Tavella, A. Eradication of caprine arthritis encephalitis virus in the goat population of South Tyrol, Italy: Analysis of the tailing phenomenon during the 2016–2017 campaign. J. Vet. Diagn. Investig. 2020, 32, 589–593. [Google Scholar] [CrossRef]
  34. Cirone, F.; Maggiolino, A.; Cirilli, M.; Sposato, A.; De Palo, P.; Ciappetta, G.; Pratelli, A. Small ruminant lentiviruses in goats in southern Italy: Serological evidence, risk factors and implementation of control programs. Vet. Microbiol. 2019, 228, 143–146. [Google Scholar] [CrossRef] [PubMed]
  35. Arcangeli, C.; Lucarelli, D.; Torricelli, M.; Sebastiani, C.; Ciullo, M.; Pellegrini, C.; Felici, A.; Costarelli, S.; Giammarioli, M.; Feliziani, F.; et al. First Survey of SNPs in TMEM154, TLR9, MYD88 and CCR5 Genes in Sheep Reared in Italy and Their Association with Resistance to SRLVs Infection. Viruses 2021, 13, 1290. [Google Scholar] [CrossRef] [PubMed]
  36. Hall, T.A. BioEdit: A user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT. Nucleic Acids Symp. Ser. 1999, 41, 95–98. [Google Scholar]
  37. Kumar, S.; Stecher, G.; Li, M.; Knyaz, C.; Tamura, K. MEGA X: Molecular Evolutionary Genetics Analysis across Computing Platforms. Mol. Biol. Evol. 2018, 35, 1547–1549. [Google Scholar] [CrossRef]
  38. Nei, M.; Kumar, S. Molecular Evolution and Phylogenetics; Oxford University Press: New York, NY, USA, 2000. [Google Scholar]
  39. 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]
  40. Ronquist, F.; Teslenko, M.; van der Mark, P.; Ayres, D.L.; Darling, A.; Höhna, S.; Larget, B.; Liu, L.; Suchard, M.A.; Huelsenbeck, J.P. MrBayes 3.2: Efficient Bayesian Phylogenetic Inference and Model Choice Across a Large Model Space. Syst. Biol. 2012, 61, 539–542. [Google Scholar] [CrossRef] [Green Version]
  41. Huson, D.H.; Bryant, D. Application of Phylogenetic Networks in Evolutionary Studies. Mol. Biol. Evol. 2006, 23, 254–267. [Google Scholar] [CrossRef]
  42. Crooks, G.E.; Hon, G.; Chandonia, J.-M.; Brenner, S.E. WebLogo: A Sequence Logo Generator: Figure 1. Genome Res. 2004, 14, 1188–1190. [Google Scholar] [CrossRef] [Green Version]
  43. Schneider, T.D.; Stephens, R.M. Sequence logos: A new way to display consensus sequences. Nucl. Acids Res. 1990, 18, 6097–6100. [Google Scholar] [CrossRef] [PubMed]
  44. Kryazhimskiy, S.; Plotkin, J.B. The Population Genetics of dN/dS. PLoS Genet. 2008, 4, e1000304. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Korber, B. HIV Signature and Sequence Variation Analysis. In Computational Analysis of HIV Molecular Sequences; Korber, B., Rodrigo, A.G., Learn, G.H., Eds.; Kluwer Academic Publishers: Dordrecht, The Netherlands, 2000; Chapter 4; pp. 55–72. [Google Scholar]
  46. Echeverría, I.; De Miguel, R.; De Pablo-Maiso, L.; Glaria, I.; Benito, A.A.; De Blas, I.; De Andrés, D.; Luján, L.; Reina, R. Multi-Platform Detection of Small Ruminant Lentivirus Antibodies and Provirus as Biomarkers of Production Losses. Front. Vet. Sci. 2020, 7, 182. [Google Scholar] [CrossRef]
  47. Kampen, A.H.; Åkerstedt, J.; Klevar, S. The Surveillance Programme for Small Ruminant Lentivirus Infections in Sheep and Goats in Norway 2020; Norwegian Veterinary Institute: Oslo, Norway, 2021; ISSN 1890-3290. [Google Scholar]
  48. De Martin, E.; Golomingi, A.; Zahno, M.; Cachim, J.; Di Labio, E.; Perler, L.; Abril, C.; Zanoni, R.; Bertoni, G. Diagnostic response to a cross-border challenge for the Swiss caprine arthritis encephalitis virus eradication program. Schweiz. Arch. Tierheilkd. 2019, 16, 93–104. [Google Scholar] [CrossRef] [PubMed]
  49. Adjadj, N.R.; Vicca, J.; Michiels, R.; De Regge, N. (Non-)Sense of Milk Testing in Small Ruminant Lentivirus Control Programs in Goats. Comparative Analysis of Antibody Detection and Molecular Diagnosis in Blood and Milk. Viruses 2019, 12, 3. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  50. 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] [Green Version]
  51. Tavella, A.; Capello, K.; Bertoni, G.; Bettini, A. Risk Factors Associated with the Alpine Multispecies Farming System in the Eradication of CAEV in South Tyrol, Italy. Viruses 2021, 13, 1959. [Google Scholar] [CrossRef]
  52. Cvijović, I.; Good, B.H.; Desai, M.M. The Effect of Strong Purifying Selection on Genetic Diversity. Genetics 2018, 209, 1235–1278. [Google Scholar] [CrossRef] [Green Version]
  53. 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]
  54. Tanaka, M.; Robinson, B.A.; Chutiraka, K.; Geary, C.D.; Reed, J.C.; Lingappa, J.R. Mutations of Conserved Residues in the Major Homology Region Arrest Assembling HIV-1 Gag as a Membrane-Targeted Intermediate Containing Genomic RNA and Cellular Proteins. J. Virol. 2016, 90, 1944–1963. [Google Scholar] [CrossRef] [Green Version]
  55. Chu, H.-H.; Chang, Y.-F.; Wang, C.-T. Mutations in the α-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] [Green Version]
  56. Minardi da Cruz, J.; Singh, D.; 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] [Green Version]
  57. Barros, S.C.; Ramos, F.; Duarte, M.; Fagulha, T.; Cruz, B.; Fevereiro, M. Genomic Characterization of a Slow/Low Maedi Visna Virus. Virus Genes 2004, 29, 199–210. [Google Scholar] [CrossRef] [PubMed]
  58. Staskus, K.A.; Retzel, E.F.; Lewis, E.D.; Silsby, J.L.; Cyr, S.S.; Rank, J.M.; Wietgrefe, S.W.; Haase, A.T.; Cook, R.; Fast, D.; et al. Isolation of replication-competent molecular clones of visna virus. Virology 1991, 181, 228–240. [Google Scholar] [CrossRef]
  59. Braun, M.J.; Clements, J.E.; Gonda, M.A. The visna virus genome: Evidence for a hypervariable site in the env gene and sequence homology among lentivirus envelope proteins. J. Virol. 1987, 61, 4046–4054. [Google Scholar] [CrossRef] [Green Version]
  60. Grego, E.; Bertolotti, L.; Carrozza, M.; Profiti, M.; Mazzei, M.; Tolari, F.; Rosati, S. Genetic and antigenic characterization of the matrix protein of two genetically distinct ovine lentiviruses. Veter. Microbiol. 2005, 106, 179–185. [Google Scholar] [CrossRef]
  61. Herrmann-Hoesing, L.M.; Broughton-Neiswanger, L.E.; Gouine, K.C.; White, S.N.; Mousel, M.R.; Lewis, G.S.; Marshall, K.L.; Knowles, D.P. Evaluation of a Caprine Arthritis-Encephalitis Virus/Maedi-Visna Virus Indirect Enzyme-Linked Immunosorbent Assay in the Serological Diagnosis of Ovine Progressive Pneumonia Virus in U.S. Sheep. Clin. Vaccine Immunol. 2010, 17, 307–310. [Google Scholar] [CrossRef] [Green Version]
  62. Huang, J.; Sun, Y.; Liu, Y.; Xiao, H.; Zhuang, S. Development of a loop-mediated isothermal amplification method for rapid detection of caprine arthritis-encephalitis virus proviral DNA. 2012, 157, 1463–1469. Arch. Virol. 2012, 157, 1463–1469. [Google Scholar] [CrossRef]
  63. Saltarelli, M.; Querat, G.; Konings, D.A.; Vigne, R.; Clements, J.E. Nucleotide sequence and transcriptional analysis of molecular clones of CAEV which generate infectious virus. Virology 1990, 179, 347–364. [Google Scholar] [CrossRef]
  64. 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]
  65. 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]
  66. Reina, R.; Grego, E.; Bertolotti, L.; De Meneghi, D.; Rosati, S. Genome Analysis of Small-Ruminant Lentivirus Genotype E: A Caprine Lentivirus with Natural Deletions of the dUTPase Subunit, vpr -Like Accessory Gene, and 70-Base-Pair Repeat of the U3 Region. J. Virol. 2009, 83, 1152–1155. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  67. Glaria, I.; Reina, R.; Ramírez, H.; de Andrés, X.; Crespo, H.; Jáuregui, P.; Salazar, E.; Luján, L.; Pérez, 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]
Figure 1. Maximum likelihood phylogenetic tree based on the alignment of 642 nt from gag-pol region of 287 sequences: 133 analysed in this study (labeled by a different color for each subtype) and 154 reference strains available in GenBank (labeled by black color). Bar: number of substitutions per site. Correspondence between sample names and accession numbers are reported in Table S1.
Figure 1. Maximum likelihood phylogenetic tree based on the alignment of 642 nt from gag-pol region of 287 sequences: 133 analysed in this study (labeled by a different color for each subtype) and 154 reference strains available in GenBank (labeled by black color). Bar: number of substitutions per site. Correspondence between sample names and accession numbers are reported in Table S1.
Viruses 14 00686 g001
Figure 2. Amino acid sequence multiple alignment of SRLVs deduced from the gag-pol fragment. Each subtype has been aligned with the respective reference sequence. Immunodominant epitopes 2 and 3, major homology region (MHR) and variable region 2 (VR-2) are reported. Dots represent the same amino acid residue. Correspondence between sample names and accession numbers are reported in Table S1. (a) Subtypes A9, A11, A24, B1, B2; (b) Subtype B3.
Figure 2. Amino acid sequence multiple alignment of SRLVs deduced from the gag-pol fragment. Each subtype has been aligned with the respective reference sequence. Immunodominant epitopes 2 and 3, major homology region (MHR) and variable region 2 (VR-2) are reported. Dots represent the same amino acid residue. Correspondence between sample names and accession numbers are reported in Table S1. (a) Subtypes A9, A11, A24, B1, B2; (b) Subtype B3.
Viruses 14 00686 g002aViruses 14 00686 g002b
Figure 3. Graphical representation of the relative amino acid frequency of the partial SRLV gag-pol protein obtained by WebLogo 3 software. The height of the letter corresponding to each amino acid indicates its relative frequency at that specific position. Different colors indicate the physiochemical characteristics of the amino acid (black: non-polar, green: polar, red: aromatic, blue: positively charged, purple: negatively charged).
Figure 3. Graphical representation of the relative amino acid frequency of the partial SRLV gag-pol protein obtained by WebLogo 3 software. The height of the letter corresponding to each amino acid indicates its relative frequency at that specific position. Different colors indicate the physiochemical characteristics of the amino acid (black: non-polar, green: polar, red: aromatic, blue: positively charged, purple: negatively charged).
Viruses 14 00686 g003
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

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. https://doi.org/10.3390/v14040686

AMA Style

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(4):686. https://doi.org/10.3390/v14040686

Chicago/Turabian Style

Arcangeli, Chiara, Martina Torricelli, Carla Sebastiani, Daniele Lucarelli, Marcella Ciullo, Fabrizio Passamonti, Monica Giammarioli, and Massimo Biagetti. 2022. "Genetic Characterization of Small Ruminant Lentiviruses (SRLVs) Circulating in Naturally Infected Sheep in Central Italy" Viruses 14, no. 4: 686. https://doi.org/10.3390/v14040686

APA Style

Arcangeli, C., Torricelli, M., Sebastiani, C., Lucarelli, D., Ciullo, M., Passamonti, F., Giammarioli, M., & Biagetti, M. (2022). Genetic Characterization of Small Ruminant Lentiviruses (SRLVs) Circulating in Naturally Infected Sheep in Central Italy. Viruses, 14(4), 686. https://doi.org/10.3390/v14040686

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