Next Article in Journal
Special Issue “Gamma Delta T Cells in Immune Response against Viruses”
Next Article in Special Issue
Expression of the Heterotrimeric GP2/GP3/GP4 Spike of an Arterivirus in Mammalian Cells
Previous Article in Journal
Efficacy of Liming Forest Soil in the Context of African Swine Fever Virus
Previous Article in Special Issue
Kinetics of the Equid Herpesvirus 2 and 5 Infections among Mares and Foals from Three Polish National Studs
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Viruses
Volume 14
Issue 4
10.3390/v14040735
viruses-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

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

Molecular Characterization of Small Ruminant Lentiviruses Isolated from Polish Goats with Arthritis

by
Monika Olech
1,2,*,
Anna Kycko
3 and
Jacek Kuźmak
2
1
Department of Swine Diseases, National Veterinary Research Institute, 24-100 Puławy, Poland
2
Department of Biochemistry, National Veterinary Research Institute, 24-100 Puławy, Poland
3
Department of Pathology, National Veterinary Research Institute, 24-100 Puławy, Poland
*
Author to whom correspondence should be addressed.
Viruses 2022, 14(4), 735; https://doi.org/10.3390/v14040735
Submission received: 25 February 2022 / Revised: 29 March 2022 / Accepted: 30 March 2022 / Published: 31 March 2022
(This article belongs to the Special Issue State-of-the-Art Animal Virus Research in Poland)
Browse Figures
Review Reports Versions Notes

Abstract

:
Previous studies revealed that the small ruminant lentivirus (SRLV) population in Poland is highly heterogeneous. All SRLVs detected from Polish sheep and goats so far have belonged to subtypes B1, B2, A1, A5, A12, A13, A16, A17, A18, A23 and A24. However, all characterized strains originated from asymptomatic animals. This is the first study that characterizes the molecular properties of SRLVs isolated from different organs of six arthritic goats. Segments from three genomic regions (gag, LTR and env) were analyzed. In addition, we quantified the SRLV proviral load in the blood and different organs and examined its association with different degrees of histopathological lesions. All sequences obtained from the goats involved in this study were homogeneous, showing an average degree of variability of 4.8%, 3.7% and 8.8% for gag, LTR and env, respectively. Phylogenetic analysis revealed that the sequences from the analyzed goats were clustered within SRLVs group A and formed a new subtype within this group, tentatively named A27. The histopathological examination of the lung, mammary gland, synovial membranes of joints and brain of the analyzed goats revealed evidence of inflammatory processes associated with SRLV infection, which was confirmed by positive immunohistochemistry assays. No significant correlation was observed between histological features and alterations in the sequences from different tissues. No tissue-specific signature pattern was identified. It was shown that animals with a higher proviral load showed more lesion severity in various SRLV-affected tissues, indicating a positive association between these two parameters. Our results also revealed differences in the SRLV load between animals even though the sequences derived from all of the goats were closely related, suggesting that host factors may restrict and control viral replication. This study provides new information about SRLV variants isolated from arthritic goats; however, more studies, including the isolation and characterization of biological properties of these viruses, should be performed to evaluate their pathogenic potential.

1. Introduction

Maedi visna virus (MVV) and Caprine arthritis encephalitis virus (CAEV) are two single-stranded RNA viruses belonging to the Lentivirus genus of the Retroviridae family that cause a lifelong persistent infection in sheep and goats. Initially, MVV and CAEV were considered pathogens specific to sheep and goats, respectively. Recently, due to their ability to cross the species barrier and their genomic and structural similarities, these two viruses are considered a single group called small ruminant lentiviruses (SRLVs). Transmission occurs mainly via the respiratory route and through the ingestion of infected colostrum/milk [1]. The main target cells of the SRLVs are monocytes/macrophages and dendritic cells, and viral replication is restricted until the maturation of monocytes to macrophages [2,3]. SRLVs cause progressive inflammatory syndromes in different host organs, including the lungs, udders, joints and central nervous system. On infection, animals may remain asymptomatic carriers for life, and only a small proportion of infected animals develop clinical signs such as interstitial pneumonia, indurative mastitis (“hard udder”), arthritis, dyspnea and, more rarely, encephalitis, ataxia or paralysis. In goats, arthritis, indurative mastitis and, more rarely, encephalomyelitis in young kids occur. Occasionally, the infected ones may develop subclinical interstitial pneumonia and progressive dyspnea [1]. Both asymptomatic and symptomatic animals can transmit the virus. Disease progression is usually slow, and the clinical course and outcome of infection are associated with the genetic background of the animals and the virulence of the infecting strain [4,5]. Histopathologically, the different SRLV disease forms involve lesions with infiltration of mononuclear cells (macrophages and lymphocytes), resulting in inflammation, fibrosis and thickening of the tissues. The severity of the histological lesions is shown to be more intense in tissues with higher viral and proviral loads [6,7]; however, the mechanism that regulates viral replication and leads to increased viral load and more severe disease outcomes is still unclear.
SRLVs are genetically heterogeneous, which may account, in part, for the different patterns of clinical diseases. SRLVs are currently divided into five phylogenetic groups, A–E, and further divided into subtypes [8,9]. The mechanism contributing to the distinct tissue specificity of SRLVs remains poorly understood. However, two areas that may contribute to tropism include long terminal repeats (LTRs), particularly the U3 promoter region, and the C-terminus of the Env protein.
The SRLV genome, which is integrated into host cells in the form of a provirus, contains three structural genes (gag, pol and env) and three regulatory genes (vpr-like, rev and vif) flanked by non-coding long terminal repeat (LTR) regions. The gag gene encodes the capsid protein, which is very conserved and contains linear epitopes, and for this reason, it is used for phylogeny and serological diagnosis. The LTR is divided into three regions, U3, R, and U5, and the U3 region contains a promoter/enhancer sequence and transcription factor binding sites, such as activator proteins (AP)-1, AP-4, AML (vis), γ-interferon activated site (GAS) and TNF-activated site (TAS) [10,11]. Several mutations, deletions or duplications in the LTR of SRLVs were shown to be associated with virulence, cellular tropism and pathogenesis [12,13,14,15,16]. Alternatively, envelope glycoproteins may also play a role in cellular host range, infectivity and disease progression. It was shown that hypervariable regions of the env gene involved in compartmentalization are implicated in tropism. In the SRLV genome, five variable (V1–V5) regions were previously identified [17]. The presence of different SRLV sequences in the V4 region in different organs suggests that SRLVs can undergo variations in this region during early infection, giving rise to different viral subpopulations. A study on the compartmentalization of SRLVs revealed the presence of different sequences in the blood, lungs, mammary glands, central nervous system and colostrum [18,19,20].
Previous studies revealed that the population of SRLVs in Poland is highly heterogeneous. The SRLVs detected so far from polish sheep and goats belonged to subtypes B1, B2, A1, A5, A12, A13, A16, A17, A18, A23 and A24 [8]. However, all characterized strains originated from asymptomatic animals. In this study, we describe SRLV sequences isolated from different organs of arthritic goats for the first time. We determined partial gag, env and LTR sequences, investigated their phylogenetic relationships with sequences of other known SRLV strains, quantified the SRLVs’ proviral load in the blood and different organs and examined its association with different degrees of histopathological lesions.

2. Materials and Methods

2.1. Animals, Serology and Tissue Sample Collection

Six naturally SRLV-infected adult goats originating from northern Poland were used for sample collection. The animals belonged to a flock with a high seroprevalence of SRLV infection as were assessed through the use of commercial test ID Screen MVV/CAEV Indirect Screening (IDvet, Grabels, France). The goats presented as thin bodied (enlargement of the carpal joint, leg weakness, ataxia) and were humanely euthanized at the owner’s request due to deterioration in their quality of life. Prior to euthanasia, peripheral blood mononuclear cells (PBMCs) were isolated from EDTA-anticoagulated blood on a Ficoll-Hypaque gradient. Blood collection was approved (no. 37/2016) by the Local Ethical Committee on Animal Testing at the University of Life Sciences in Lublin (Poland). A complete set of tissues, including joints, lungs, mammary glands, brain and choroid plexus, were collected post-mortem in a thermos flask (+4 °C) in individual containers filled with Dulbecco’s phosphate-buffered saline (DPBS) for viral isolation. Additionally, samples from the brain, including choroid plexus, lungs and mammary glands, were collected for nucleic acid extraction and histopathological analysis. In addition, synovial fluids and synovial membranes obtained from joints were used for virus isolation and molecular characterization, respectively.

2.2. Histopathology and Immunohistochemistry

Tissue samples were promptly fixed in 10% buffered formalin, then routinely processed, embedded in paraffin and cut on a microtome for histopathology and immunohistochemistry (IHC) examination. For histopathology, 4 μm thick slides were stained with hematoxylin and eosin (HE). The intensity of the lesions in the target organs was estimated independently by two pathologists. The severity of arthritis/synovitis among the tissues from the examined goats was scored according to the grading system proposed by Cheevers et al. [21]. The severity of the inflammatory infiltration in the examined sections of the lungs and mammary glands was scored as follows: 0—No changes, 1—Mild diffuse inflammatory infiltration, 2—Moderate diffuse inflammatory infiltration, 3—Diffuse inflammatory infiltration with the formation of lymphoid follicles. The severity of lesions observed in the brain was scored as follows: 0—No changes, 1—Mild lymphoplasmacytic infiltration within choroid plexus, 2—Moderate lymphoplasmacytic infiltration within choroid plexus and focal perivascular cuffings, mild gliosis.
Select tissues with lesions histologically consistent with SRLV-associated inflammation were further evaluated by an immunohistochemistry assay to detect the SRLV p28 Gag antigen. IHC was performed using the DAB/HRP Envision system (DAKO, Glostrup, Denmark). Briefly, 4 μm thick paraffin-embedded tissue sections were cut and mounted on positively charged slides. After overnight air-drying, the slides were deparaffinized, rehydrated using xylene and graded alcohols and then treated with a 3% hydrogen peroxide solution for 10 min to block the activity of endogenous peroxidases. Antigen retrieval was carried out by treatment in a pressure cooker (110 °C) with citrate buffer pH 6 for 20 min. The slides were then incubated with the CAEV10A1 primary antibody (VMRD) diluted 1:100 for 60 min at room temperature (RT), followed by incubation with a peroxidase-conjugated polymer as a secondary antibody (Dako REAL EnVision Detection System, K5007, DAKO, Glostrup, Denmark) for 30 min. After color development with diaminobenzidine solution, sections were counterstained with Mayer’s hematoxylin, dehydrated and mounted. Sections incubated with PBS instead of the primary antibody were used to confirm the specificity of the staining. The tissues were analyzed under a light microscope (Axiolab, Zeiss, Germany) for the presence of positive labeling.

2.3. DNA Extraction, Amplification and Sequencing

Genomic DNA was extracted from tissue samples using a NucleoSpin Blood Quick Pure Kit (Macherey-Nagel GmbH & Co KG, Dueren, Germany), according to the manufacturer’s instructions. The quality and quantity of the DNA were checked in a Nanophotometer (Implen GmbH, Munich, Germany). The CA (625 bp) fragment of the gag gene fragment, the V4V5 fragment of the env gene and the U3-R fragment of the LTR region were amplified by nested PCRs, as previously described [8]. A water template negative control was run parallel with each PCR reaction set. PCR products were purified using NucleoSpin Gel and PCR Clean-up (Marcherey-Nagel, GmbH 7 Co., Hamburg, Germany) and directly sequenced on a 3730 xl DNA Analyzer (Applied Biosystems, Foster City, CA, USA) using a BigDye Terminator v3.1 Cycle Sequencing kit. The obtained SRLV sequences were edited and aligned against known reference sequences using the MUSCLE multiple-alignment software provided as part of Geneious Pro 5.3 (Biomatters Ltd., Auckland, New Zealand). All novel sequences reported in this study were submitted to the Gen-Bank database under accession numbers: OM517045-OM517079 for the gag sequences, OM517080-OM517094 for the env sequences and OM517095-OM517127 for the LTR. The best-fitting nucleotide substitution models were estimated, and the GTR, TN93 and K2 with the gamma distribution rates were used to infer a phylogenetic tree using the maximum likelihood (ML) method for the env, gag and LTR sequences, respectively. The reliability of the phylogenetic relationships was evaluated by nonparametric bootstrap analysis with 1000 iterations. Alignment, model testing, and tree building were performed using the MEGA 6 application [22]. Nucleotide and amino acid sequence percent identity (percentage of bases/residues which are identical) was estimated using the Geneious software, while pairwise genetic distances were calculated with the MEGA 6 software.

2.4. Proviral Load Quantification

Proviral DNA was quantified by real-time PCR using Rotor-Gene Q Series ver. 2.0.3 (Qiagen, Hilden, Germany) with primers and a probe specifically designed based on the gag sequences obtained for the SRLV subtype described in this study. The sequence of the forward and reverse primers and probe were CAPBWF (5′GAAGCAGAAAGGTGGGTAAGA 3′), CAPBWR (5′ TTGCGATGCCTGTTGATTTG3′) and CAPBWP (5′ 6-FAM-ATCCACTGTGAGGACATTTGGCCC -BHQ-1 3′), respectively. The 50 µL reaction mixture contained 25 μL of 2× QuantiTect Multiplex NoROX PCR buffer (Qiagen, Hilden, Germany), 400 nM of each primer, 200 nM of the TaqMan probe and 500 ng of extracted genomic DNA. The cycling protocol involved an initial incubation for polymerase activation at 95 °C for 15 min, followed by 45 cycles of denaturation at 94 °C for 60 s and annealing/extension at 60 °C for 60 s. A non-template control (DEPC H2O) was included in each run. Standard curves were generated using a 10-fold serial dilution from 108 to 10 copies of the reference plasmid encompassing the target gag region. All samples were tested in duplicate, and the results were expressed as a mean copy number of provirus per 500 ng of genomic DNA of each goat.

3. Results

3.1. Histopathology and Immunohistochemistry (IHC)

All of the goats presented macroscopic lesions compatible with chronic proliferative arthritis that consisted of enlarged carpal joints. Histopathological changes observed in the joints involved synoviocyte proliferation and synovial membrane thickening with the formation of papillary Villar projections, accompanied by lymphoplasmacytic infiltrations of the subsynovial connective tissue and perivascular spaces (Figure 1E) as well as diffuse fibrosis. Necrosis with focal calcification was occasionally observed within the fibrous capsule. The severity of the lesions observed in the joints is summarized in Table 1.
In the lungs, mild changes characteristic of interstitial pneumonia, including interstitial lymphocytic infiltrations and thickening of the alveolar septa, were observed (Figure 1C). In the mammary glands, around the acini and the lactiferous ducts, diffuse mononuclear inflammatory infiltrates of variable intensity were present, accompanied by multifocal inspissation and the occasional mineralization of the duct contents (Figure 1G). Finally, lesions in the brain were scant and seen only in three out of six analyzed goats. They were characterized by mild lymphoplasmacytic infiltration of the choroid plexus and, rarely, focal perivascular cuffings and gliosis in the brainstem (Figure 1A,B). The scored severity of the inflammatory infiltration in the examined sections of the lungs, mammary glands and brains of the goats is shown in Table 1.
Immunohistochemical examination showed positive staining for SRLV p28 in the articular tissues in several inflammatory cells, in the lungs in several bronchiolar epithelial cells and single mononuclear cells, as well as in single macrophage-like cells in the mammary glands (Figure 1). IHC performed on the brain tissues was negative. There was no positive immunolabeling detected in any of the examined brain sections.

3.2. Phylogenetic Analysis

Total DNA isolated from the peripheral blood mononuclear cell (PBMC), mammary gland, lung, synovial membrane and brain of all goats were used as templates for the PCR amplification of the gag gene encoding the capsid protein. The band corresponding to the expected 625 bp PCR product was detected in all tissue samples of the animals except the lung-derived samples originating from goat 5. The obtained sequences were aligned with other published sequences representing SRLV genotypes described to date and subjected to phylogenetic analysis. The maximum-likelihood (ML) tree revealed that all analyzed sequences were closely related to each other (mean genetic distance of 4.8%) and that the sequences derived from different tissues of the same animal clustered together. The mean nucleotide distance of sequences derived from different tissues of each goat did not exceed 1.3% (ranging from 0.2% to 1.3%). Interestingly, all of the analyzed sequences formed a new cluster within group A, which could be tentatively named A27. The affiliation of the new cluster, A27, was supported with high bootstrap values ≥85 (Figure 2). The mean genetic distances between sequences forming the new subtype and sequences representing other subtypes representative of genotype A varied from 10.0% to 21.7% (Table 2). Polish SRLV A27 nucleotide sequences showed the highest nucleotide sequence identity, with strains representative of subtype A13. To evaluate the robustness of our analysis, we also performed a phylogenetic analysis using the neighbor-joining method (Supplementary Materials Figure S1), which resulted in the same classification, supporting the existence of the new subtype.
In order to confirm the genetic assessment of the sequences obtained from the analyzed goats, we amplified and sequenced a genomic fragment encompassing a highly variable region of the envelope gene (env) encoding for the carboxy-terminus portion its surface subunit (SU). Unfortunately, amplification of this fragment succeeded only for the lung-derived samples of goat 4 and 6, brain-derived samples of goat 4 and 5, choroid plexus-derived samples of goat 5, mammary gland-derived samples of goat 4, 5 and 6, synovial membrane-derived samples of goat 2, 4, 5, 6 and 7 and PBMC-derived samples of goat 2 and 5. The phylogenetic tree (Figure 3) confirmed that the sequences originating from the goats analyzed in this study belonged to the new subtype A27. The mean nucleotide distance between sequences belonging to the A27 and other subtypes within group A ranged from 23.7% to 31.8%. The mean genetic distance of the env sequences originating from all goats was 8.8% and ranged from 4.7% to 12.6%. Within each animal, the mean nucleotide distance of sequences from different tissues did not exceed 7.5% (ranging from 0.8% to 7.5%).

3.3. Comparative Analysis of Immunodominant Regions

Nucleotide sequences of the gag and env genes were translated into amino acid (aa) sequences and aligned with the aa sequences of reference strains: Cork and K1514. The pairwise percent identity of the gag amino acid sequences of the SRLVs obtained from the tested goats was high and ranged from 91.7% to 100%. The pairwise percent identity of sequences originating from the same animal ranged from 97.4% to 100%. Furthermore, the analyzed sequences shared 84.2–86.8% and 93.1–96.1% amino acid sequence identity with the Cork and K1514 strains, respectively. Analysis revealed that all CA sequences originating from all tested goats had the asparagine-valine (NV) motif, which is specific to MVV-like (group A) sequences (Figure 4). Epitope 2 was fully conserved in all sequences originating from different tissues of all goats. Sequences of epitope 3 were also quite conserved. Only glutamic acid (E) replaced lysine (K) in sequences derived from the choroid plexus and brain of goat 6, and arginine (R) replaced lysine (K) in all sequences derived from goat 4. More alterations were found in the major homology region (MHR), which is highly conserved in all retroviruses. No tissue-specific signature pattern was identified in the gag sequences. Most of the substitutions were specific for each animal. In particular, all sequences from goat 3 had the substitution of lysine (K) instead of arginine (R) at position 137 compared to the sequences originating from the other goats, aspartic acid (D) instead of glutamic acid (E) and histidine (H) instead of glutamine (Q) in positions 7 and 98, respectively. All sequences originating from different tissues of goat 2 had the unique substitutions S41N, K129R and K139R, while all sequences from goat 4 had the unique substitutions K41R, E46D, A101S, N161T, A177S and S178T. The sequences from goat 7 had the unique substitutions M118T, N141S, V188T and T191N, while the sequences from goats 5 and 6 had the unique substitution I148V. Furthermore, all sequences from goats 2 and 3 had the specific substitution isoleucine (I) instead of threonine (T) at position 10 compared to the sequences from the other goats.
Next, the deduced amino acid sequences of the immunodominant regions of the surface glycoproteins (SU) were aligned with the corresponding reference sequences of parental strains Cork and K1514 (Figure 5). The env aa sequences obtained from the goats were more heterogeneous, showing 77.6–98.3% similarity to each other and 65.8–71.0% and 77.0–79.0% to the Cork and K1514 strains, respectively. Most of the variation was observed in the variable region (V4), especially in the highly variable region (HV2) (Figure 5). Sequences of epitope SU5 were more conserved among the sequences derived from the analyzed goats and differed prominently from the SU5 sequences of the Cork and K1514 strains and other strains representing known subtypes (Supplementary Materials Figure S2). Only the sequence originating from the choroid plexus of goat 5 differed from the sequences originating from the other goats analyzed in this study. All sequences had a perfectly conserved motif located at the N-terminus part (VRAYTYGV) and more the variable motif at the C-part.

3.4. Analysis of LTR Sequences

LTR sequences from the lungs, choroid plexus, mammary glands, synovial membrane and PBMC of all six goats and from the brain of goats 4, 5 and 7 were successfully amplified and were aligned with the corresponding sequences of SRLV strains with defined virulence, including sequences of the virulent B1 strain Cork [24], A2/A3 strain 697 isolated from sheep with neurological symptoms [25], neurovirulent A1 strain Kv1772 [13], B2 strain 496 isolated from arthritic sheep [26] as well as to the attenuated A4 strain 6221 [27].
The LTR sequences originating from all goats were closely related to each other. The mean genetic distance of the LTR sequences was 3.9%. Only the sequence derived from the choroid plexus of goat 4 was quite divergent from the rest of the sequences. The genetic distance of this sequence, compared to the other sequence from different tissues of all analyzed goats, ranged from 14.5% to 16.8%. The mean nucleotide distance of the sequences from different tissues of each goat did not exceed 4.3%. Only in the case of goat 4 the mean nucleotide distance of sequences from different tissues was 6.5% and ranged from 0.5% to 16%. Phylogenetic analysis revealed that sequences derived from different tissues of the same animal generally clustered together (Figure 6). Only goat 4 had distinct sequences isolated from the choroid plexus and synovial membrane; goat 5 had a distinct sequence isolated from the lungs, while sequences isolated from the PBMC and mammary gland of goat 6 were quite different from the sequence isolated from the choroid plexus, synovial membrane and lungs.
Multiple conserved motifs were identified within the LTR sequences, including AP-1 sites, AML, AP-4, TATA box, poly A, CAAAT, GAS and TAS. Sequences corresponding to the TATA box, AP-4, GAS and polyadenylation signal (poly A) were quite conserved among all sequences (Figure 7). Almost all sequences originating from the analyzed goats had a unique T to A substitution in the fifth position of the TATA box. Five different sites of AP-1 (TGACAC, G/TAGTCA, TCATGTA, TTAG/AGTCA, ATA/GAT/CTA/GT/G) were observed in sequences derived from all goats analyzed in this study. Furthermore, all tissue sequences from all goats presented two highly conserved AML(vis) sites. Some point mutations were observed in the TAS region; however, the sequences originating from different tissues of the same animal were very homologous. Two deletions in the U3 region of the A5 strains were noted; one (10 nt) located near the 5′ end of the U3 region, and the second (21 nt) located in the central region of the U3. The presence of these deletions was also shown in attenuated strain 6221 as well as in virulent strain 697. The sequence derived from the choroid plexus of goat 4 had a 12-nt deletion in the R region, while all sequences from goat 3 had a 9-nt deletion in this track. Similar deletions were present in virulent strain 697 as well as attenuated strain 6221. The number of repeats of the CAAAT sequence was different depending on the strains. Strain Kv1772 had three copies of the CAAAT sequence, while Cork, strain 496 (B2) and 497 (A2/A3) had only one copy of the CAAAT, and strain 6221, representing subtype A4, had no copy. The number of CAAAT sequences derived from different tissues of the analyzed goats varied from 1 to 3. Only the sequences derived from the mammary gland and brain of goat 5 had three copies of CAAAT.

3.5. Provirus Detection and Quantification

The proviral load in the blood and different organs of infected animals was quantified by real-time PCR using primers and probes specifically designed based on the gag gene sequences of SRLVs described in this study. The number of provirus copies was determined by comparing the Ct values obtained by amplifying gag fragment in a particular sample to the standard curve generated after amplification of scale-down dilutions of plasmid DNA containing the targeted region. A new standard curve was generated for every experiment, showing a mean efficiency value close to 100% and a mean R-value of R = 0.998. The provirus copy numbers differed between the animals and the tissue analyzed (Figure 8). The highest proviral load was detected in the synovial membrane of goat 7 (4618 copies per 500 ng of genomic DNA). In goats 2, 3, 4 and 7, the highest proviral load was found in the synovial membrane cells. The copy number detected in these samples ranged from 47 to 3466. In goats 6 and 5, the highest proviral load was found in the lung and choroid plexus, respectively. Brain samples were positive only in four out of six goats and showed very low copy numbers (>4 copies per 500 ng of genomic DNA). The proviral load (>8 copies per 500 ng of genomic DNA) detected in the choroid plexus also was very low in 4 out of 6 analyzed goats, while in goats 5 and 7, the number of SRLV provirus copies detected in this tissue was 2233 and 65, respectively. In the mammary gland, PBMC and lung samples, the SRLV copy number varied from 44 to 351, from 9 to 1158 and from 9 to 820, respectively.

4. Discussion

In this study, we revealed, for the first time, the gag, LTR and env sequences of SRLVs detected in Polish goats with clinical signs of arthritis and determined their phylogenetic relationship in the context of known SRLV sequences. In addition, we examined the potential association between the degree of histopathological changes and provirus concentration in different organs of the infected goats.
All sequences obtained from the goats involved in this study were homogeneous, showing an average degree of variability of 4.8%, 3.7% and 8.8% for the gag, LTR and env sequences, respectively. Phylogenetic analysis on the basis of the gag and env sequences revealed that the sequences from the analyzed goats clustered within SRLV group A and formed a new subtype within this group, tentatively named A27. The SRLVs isolated until now from asymptomatic sheep and goats in Poland belonged to the well-known subtypes B1, B2, A1, A5 and A16, as well as subtypes A12, A13, A17, A18, A23 and A24, detected only in Poland [8,28,29,30,31]. The sequences obtained in this study were grouped in a cluster with a significant divergence from the other strains isolated from goats and sheep from Poland. The sequences were closely related to sequences belonging to subtype A13, which was previously detected only in sheep from Poland, suggesting their common origin. It was not known if the examined goats had contact with sheep, but infection with A subtypes in goats was reported [9,27], also in Poland, where the circulation of subtypes A1, A5, A12, A16, A17 and A23 [8,28,29,31] was detected in goats.
Indirect serological assays usually use the capsid protein as the antigen, so analysis of these sequences is useful for detecting antigenic variability, which may affect the sensitivity of serological tests. The primers used in this study allowed the sequencing of two major immunodominant epitopes identified by Rosati et al. [32] and the major homology region (MHR). Immunodominant epitope 2 was fully conserved in all sequences originating from different tissues of all goats, and only minor variations were found in sequences of epitope 3. This high conservation is important for the maintenance of cross-reactivity in serological tests based on Gag-derived antigen [17,33]. More alterations were observed in the major homology region (MHR), suggesting that this region is quite variable. Therefore, our results are in line with other studies demonstrating the variability of the MHR region [8,34]. No tissue-specific signature pattern was identified in the gag sequences. Most of the substitutions were specific to each animal, suggesting that these changes could have arisen as a result of long-term host–virus adaptation and evolution.
The env aa sequences obtained from the analyzed goats were more heterogeneous. Most amino acid substitution was observed in variable region V4, especially in highly variable region HV2, previously proposed to be part of a variable, conformational neutralization epitope [35,36]. It is proposed that the V4 region forms a highly constrained and surface-exposed domain and that a cysteine loop may have an analogous function to the V3 domain of HIV-1. Recently, a ‘signature pattern’ related to the SRLV genotype and different clinical statuses in sheep and goats was found in the V4 region. Possible signature patterns were detected at positions 54, 78, 79 and 80 in the sequence loop of the V4 region. However, the most important statistical significance was found at position 54, where residue N occurred only in sequences derived from arthritic animals infected with genotype B, while residue T and G occurred only in sequences derived from asymptomatic animals infected with genotype A and B, respectively [37]. Our results showed the presence of residue T at position 54 in all sequences derived from the different tissues of the analyzed goats. This may suggest that this position is not associated with clinical signs but may be associated with SRLV genotypes as all sequences with residue T at position 54 are derived only from animals infected with genotype A. We did not find a tissue-specific signature pattern in the env sequences. However, we were unable to obtain all sequences derived from the different tissues of the analyzed goats. PCR amplification detected env proviral DNA mostly in synovial membrane tissues, whereas samples from other tissues yielded limited positive results. This indicates a high variability of env sequences in different tissues, which may arise from compartmentalization and rapid sequence evolution during SRLV infection [18,19]. A study on SRLV compartmentalization revealed the presence of different sequences in the blood, lungs, mammary glands, central nervous system and colostrum [18,19,20]. Our study revealed limited genetic variability within the host. Only the sequence originating from the choroid plexus of goat 5 differed from the sequences originating from the other tissues of this goat as well as sequences originating from other goats. Moreover, the fact that the highest SRLV proviral load and moderate lesions were only observed in the choroid plexus of this goat may suggest that env sequences may be involved in tissue tropism, viral compartmentalization and pathology.
The U3 region of lentiviruses contains promoter sequences and transcription factor binding sites important for virus transcription and replication. It was demonstrated that variability in the LTR sequence is associated with tissue tropism and disease outcomes in a number of retroviruses [13,38,39]. However, in this study, we did not identify any tissue-specific motifs. Sequences isolated from the same tissue in different goats did not aggregate into a phylogenetic cluster. Our results revealed a rather coherent grouping of sequences obtained from different anatomical compartments of the same animals, indicating that the LTR sequences are not good predictors of viral tissue tropism. Sequences corresponding to the TATA box, AP-4, AML (vis) and polyadenylation signal (poly A) were conserved among all sequences derived from different tissues of the analyzed goats. Furthermore, our results revealed that almost all sequences derived from different tissues of the analyzed goats had a unique T to A substitution in the fifth position of the TATA box. The same mutation was also noted for Polish strains belonging to subtype A17 [8], but the meaning of this mutation is unknown. More alterations were observed in sequences of AP-1 sites. Additionally, GAS and TAS response elements for some cytokines were identified in the U3 region of all sequences. It was demonstrated that the pro-inflammatory cytokines, IFN-γ and tumor necrosis factor (TNF)-α, which are increased in goats with arthritis, activate viral transcription through the GAS and TAS sites [40]. In our study, GAS was identical in all sequences derived from the goats. These findings are compatible with previous studies showing absolute conservation of the GAS region in 40 of 41 viral promoter sequences, suggesting that this motif is important for viral pathogenesis [16]. The TAS region was less conserved; however, sequences originating from different tissues of the same animal were very homologous. Generally, our results are consistent with previous studies demonstrating only minimal variation in U3 sequences derived from different anatomical sites of a single host and that the TAS and AP-1 sites are the least conserved transcriptional features of the U3 region [10,15,16,41].
Previous studies also showed that duplication in the LTR U3 region of Icelandic neurological strains is a determinant of cell tropism. It was demonstrated that deletion in the duplicated sequence in this region caused deficient virus replication in choroid plexus cells, indicating that this duplication may be associated with neurovirulence [13]. Our results revealed the absence of 10-nt and 21-nt insertions in the U3 region of all sequences originating from the goats analyzed in this study, insertions that were present only in the Kv1772 strain. The lack of a 21-nt repeat resulted in a reduction of potential transcription factor binding sites (AP-1 and AML(vis)) and the CAAAT sequence located in this repeat. No putative transcription factor binding domains were identified within the 10-nt deletion. However, the presence of these two deletions was also shown in attenuated strain 6221 as well as in neuropathogenic strain 697, suggesting that these insertions are not related to a particular host tropism. Furthermore, our results revealed that the sequence derived from the choroid plexus of goat 4 had a 12-nt deletion in the R region, while all sequences from different tissues of goat 3 had a 9-nt deletion in this track. This deletion was previously described by Angelopoulos et al. [41], who suggested that the presence of this specific deletion may be associated with lower pathogenicity of SRLV strains. However, no substantial difference was found in the LTR region of small ruminants with and without clinical signs of disease. Our results revealed that a similar deletion was present in sequences of attenuated strain 6221 and strain 697 isolated from sheep from a neurological outbreak. Moreover, this deletion was found previously in the R region of the Polish strains isolated from clinically healthy sheep and goats [8] and clinically affected sheep from Spain [26], suggesting that factors besides this deletion are involved in avoiding the appearance of clinical symptoms. Generally, all of these observations are in line with the observations of Blatti-Cardinaux et al., who indicated that the role of LTR in tissue tropism is still unclear, suggesting that other viral sequences might be involved in tissue tropism and pathology [42].
SRLV proviral DNA was detected in all examined goats confirming the persistent nature of SRLV infection. In four out of six goats, the highest proviral load was found in the synovial membrane cells, confirming the tropism of isolated viruses to the joints. The relatively high proviral load was also observed in the mammary glands of all goats and the lungs of three goats. An increased proviral load in these affected organs was expected since they are the primary target organs of SRLVs and are directly involved in the efficient transmission of SRLVs through the colostrum, milk and respiratory exudates [43,44]. Several types of cells in the mammary glands, including macrophages and epithelial cells, harbor the virus playing an important role in its transmission to offspring. Additionally, SRLVs belonging to genotypes A4 and E1, which express low pathogenicity for goats, were abundantly presented in the mammary glands, confirming that lactogenic transmission is the main route of infection of SRLVs in this animal species [43,44]. It was also shown that SRLV subtype A is more efficiently transmitted to offspring, suggesting that A group-derived subtypes are particularly efficient in lactogenic transmission [45]. The lowest proviral load was detected in DNA extracted from the choroid plexus and brain cells, suggesting poor neurotropism of these viruses. The highest proviral load was only found in the choroid plexus cells of goat 5; however, no neurological symptoms were observed. High proviral loads in the absence of overt clinical signs in the infected animals are similar to the situation observed in simian immunodeficiency virus (SIV)-infected natural hosts and are probably related to the long coevolution between this lentivirus and their respective hosts [43,46]. Studies on neurovirulent variants of SIV showed that macrophage tropism is not sufficient for neurological disease [46].
Many studies describe a phylogenetic analysis of SRLVs, but only a few of them describe the histopathological lesions observed in the target organs of the infected goats/sheep [6,34,41,43,47]. The pathogenesis of the SRLV strains found in arthritic goats in Poland was not characterized in any case. The histopathological examination of the lungs, mammary glands, synovial membranes of joints and brain of the analyzed goats revealed evidence of inflammatory processes associated with SRLV infection, which was confirmed by positive immunohistochemistry assays, amplification and sequencing. It was shown that animals with a higher proviral load showed more lesion severity in various SRLV-affected tissues, indicating that the proviral load is positively correlated with the presence and severity of clinical disease symptoms [6,7]. Therefore, the goal of this work was to examine SRLVs’ target organs for their proviral load and the potential presence of histopathological lesions. Comparative analysis of the distribution of lesions according to the lesion degree and proviral load performed in this study indicated a positive association between these two parameters. This was most clearly seen in the case of the synovial membrane and brain tissues. The more severe lesions were observed in the synovial membrane of three goats, which indeed showed very high proviral loads. Furthermore, moderate lymphoplasmacytic infiltration within the choroid plexus and focal perivascular cuffings with mild gliosis was only observed in two goats, which showed higher proviral loads. Our results also revealed differences in the SRLV loads between animals even though the sequences derived from all of the goats were closely related. Goats 3, 4 and 6 had relatively lower proviral loads in all tissue samples, while goats 2, 5 and 7 had higher overall proviral loads. This different reactivity to SRLV infection suggests that host factors may restrict and control viral replication. We may speculate that goats with a lower proviral load mounted an immune response that was more efficient, which resulted in a reduced viral load. Several works identified immune response loci that could influence resistance/susceptibility to SRLV infection and disease, providing evidence that genetic factors might modulate the outcome [48,49,50,51,52,53,54,55,56]. It was demonstrated that host control of infection with SRLVs, including the provirus level, may have a genetic basis. Molecular studies in humans infected with HIV pointed to the role of SNPs in the Toll-like receptor (TLR) gene, especially in TLR7, and their association with viral load. In infection with SRLVs, it was also demonstrated that some polymorphisms were significantly associated with the SRLV provirus concentration. The obtained results suggest that SNPs of TLR7/8 genes may induce differential innate immune responses toward SRLVs, affecting the proviral concentration and, thereby, disease progression [57].
In conclusion, this is the first study investigating the association between the SRLV sequences analysis and histopathological lesions observed in different tissues of clinically affected goats from Poland. No significant correlation was observed between histological features and alterations in the sequences from different tissues. No tissue-specific signature pattern was identified. It was shown that the provirus load is positively correlated with lesion severity in various SRLV-affected tissues. The results of the phylogenetic analysis revealed the existence of a putative new subtype which should lead to the consideration of an update of the current SRLV classification. The A27 subtype in goats induces histopathological changes that are typical of SRLVs in general and contribute toward ongoing interspecies transmission and viral persistence. Our results indicated that the relationship between small ruminants and SRLVs is complex and that the severity of infection varies depending on the viral strain and the genetic background of the host. More extensive studies, including the isolation and characterization of the biological properties of these viruses, should be performed to evaluate their pathogenic potential.
The findings of this study have to be seen in light of some limitations. The number of tested animals was small (n = 6); all goats originated from the same flock and all animals were infected with the same SRLV subtype (A27). Furthermore, a limitation of this study is undoubtedly a lack of all env gene sequences.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/v14040735/s1, Figure S1: Neighbor-joining phylogenetic tree based on the alignment of CA fragment of gag gene; Figure S2: Alignment of deduced amino acid sequences of the immunodominant epitope of ENV protein and variable region V4 of the SRLVs obtained in this study and reference strains representing different subtypes.

Author Contributions

Conceptualization, M.O. and J.K.; methodology, M.O. and A.K.; formal analysis, M.O.; investigation, M.O. and A.K.; resources, M.O. and J.K.; writing—original draft preparation, M.O.; writing—review and editing, M.O., A.K. and J.K.; visualization, M.O. and A.K.; supervision, M.O.; funding acquisition, M.O. All authors have read and agreed to the published version of the manuscript.

Funding

The study was funded by the National Veterinary Research Institute (NVRI), Pulawy, Poland, within the “NVRI Research Fund” awarded to project no. F/104.

Institutional Review Board Statement

The study was conducted in accordance with the relevant guidelines and regulations. Specifically, blood collection was approved (no. 37/2016) by the Local Ethical Committee on Animal Testing at the University of Life Sciences in Lublin (Poland).

Informed Consent Statement

Not applicable.

Data Availability Statement

All data generated and analyzed in this study are included in this article.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. 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] [PubMed]
  2. Narayan, O.; Wolinsky, J.S.; Clements, J.E.; Strandberg, J.D.; Grin, D.E.; Cork, L.C. Slow virus replication: The role of macrophages in the persistence and expression of visna viruses of sheep and goats. J. Gen. Virol. 1982, 59, 345–356. [Google Scholar] [CrossRef] [PubMed]
  3. Ryan, S.; Tiley, L.; McConnell, I.; Blacklaws, B. Infection of dendritic cells by the Maedi-Visna lentivirus. J. Virol. 2000, 74, 10096–10103. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Heaton, M.P.; Clawson, M.L.; Chitko-Mckown, C.G.; Leymaster, K.A.; Smith, T.P.; Harhay, G.P.; White, S.N.; Herrmann-Hoesing, L.M.; Mousel, M.R.; Lewis, G.S.; et al. Reduced lentivirus susceptibility in sheep with TMEM154 mutations. PLoS Genet. 2012, 8, e1002467. [Google Scholar] [CrossRef] [Green Version]
  5. Sider, L.H.; Heaton, M.P.; Chitko-McKown, C.G.; Harhay, G.P.; Smith, T.P.; Leymaster, K.A.; Laegreid, W.W.; Clawson, M.L. Small ruminant lentivirus genetic subgroups associate with sheep TMEM154 genotypes. Vet. Res. 2013, 44, 64. [Google Scholar] [CrossRef] [Green Version]
  6. Ravazzolo, A.P.; Nenci, C.; Vogt, H.R.; Waldvogel, A.; Obexer-Ruff, G.; Peterhans, E.; Bertoni, G. Viral load, organ distribution, histopathological lesions, and cytokine mRNA expression in goats infected with a molecular clone of the caprine arthritis encephalitis virus. Virology 2006, 350, 116–127. [Google Scholar] [CrossRef]
  7. Herrmann-Hoesing, L.M.; Noh, S.M.; White, S.N.; Snekvik, K.R.; Truscott, T.; Knowles, D.P. Peripheral ovine progressive pneumonia provirus levels correlate with and predict histological tissue lesion severity in naturally infected sheep. Clin. Vaccine Immunol. 2009, 16, 551–557. [Google Scholar] [CrossRef] [Green Version]
  8. 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]
  9. 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]
  10. Gomez-Lucia, E.; Rowe, J.; Collar, C.; Murphy, B. Diversity of caprine arthritis-encephalitis virus promoters isolated from goat milk and passaged in vitro. Vet. J. 2013, 196, 431–438. [Google Scholar] [CrossRef]
  11. Gomez-Lucia, E.; Barquero, N.; Domenech, A. Maedi-Visna virus: Current perspectives. Vet. Med. 2018, 9, 11–21. [Google Scholar] [CrossRef] [Green Version]
  12. 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]
  13. Oskarsson, T.; Hreggvidsdottir, H.S.; Agnarsdottir, G.; Matthiasdottir, S.; Ogmundsdottir, M.H.; Jónsson, S.R.; Georgsson, G.; Ingvarsson, S.; Andrésson, O.S.; Andrésdóttir, V. Duplicated sequence motif in the long terminal repeat of maedi-visna virus extends cell tropism and is associated with neurovirulence. J. Virol. 2007, 81, 4052–4057. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Agnarsdóttir, G.; Thorsteinsdóttir, H.; Óskarsson, T.; Matthíasdóttir, S.; St Haflidadóttir, B.; Andrésson, Ó.S.; Andrésdóttir, V. The long terminal repeat is a determinant of cell tropism of maedi-visna virus. J. Gen. Virol. 2000, 81, 1901–1905. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Adedeji, A.O.; Barr, B.; Gomez-Lucia, E.; Murphy, B. A polytropic caprine arthritis encephalitis virus promoter isolated from multiple tissues from a sheep with multisystemic lentivirus-associated inflammatory disease. Viruses 2013, 5, 2005–2018. [Google Scholar] [CrossRef] [Green Version]
  16. Murphy, B.; McElliott, V.; Vapniarsky, N.; Oliver, A.; Rowe, J. Tissue tropism and promoter sequence variation in caprine arthritis encephalitis virus infected goats. Virus Res. 2010, 151, 177–184. [Google Scholar] [CrossRef]
  17. Valas, S.; Benoit, C.; Baudry, C.; Perrin, G.; Mamoun, R.Z. Variability and immunogenicity of caprine arthritis encephalitis virus surface glycoprotein. J. Virol. 2000, 74, 6178–6185. [Google Scholar] [CrossRef] [Green Version]
  18. Ramírez, H.; Reina, R.; Bertolotti, L.; Cenoz, A.; Hernández, M.M.; San Román, B.; Glaria, I.; de Andrés, X.; Crespo, H.; Jáuregui, P.; et al. Study of compartmentalization in the visna clinical form of small ruminant lentivirus infection in sheep. BMC Vet. Res. 2012, 8, 8. [Google Scholar] [CrossRef] [Green Version]
  19. Pisoni, G.; Moroni, P.; Turin, L.; Bertoni, G. Compartmentalization of small ruminant lentivirus between blood and colostrum in infected goats. Virology 2007, 369, 119–130. [Google Scholar] [CrossRef] [Green Version]
  20. Olech, M.; Kuźmak, J. Compartmentalization of Subtype A17 of Small Ruminant Lentiviruses between Blood and Colostrum in Infected Goats Is Not Exclusively Associated to the env Gene. Viruses 2019, 11, 270. [Google Scholar] [CrossRef] [Green Version]
  21. Cheevers, W.P.; Snekvik, K.R.; Trujillo, J.D.; Kumpula-McWhirter, N.M.; Top, K.P.O.; Knowles, D.P. Prime-boost vaccination with plasmid DNA encoding caprine-arthritis encephalitis lentivirus env and viral SU suppresses challenge virus and development of arthritis. Virology 2003, 306, 116–125. [Google Scholar] [CrossRef]
  22. Tamura, K.; Stecher, G.; Peterson, D.; Filipski, A.; Kumar, S. MEGA6: Molecular Evolutionary Genetics Analysis version 6.0. Mol. Biol. Evol. 2013, 30, 2725–2729. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Colitti, B.; Coradduzza, E.; Puggioni, G.; Capucchio, M.T.; Reina, R.; Bertolotti, L.; Rosati, S. A new approach for Small Ruminant Lentivirus full genome characterization revealed the circulation of divergent strains. PLoS ONE 2019, 14, e0212585. [Google Scholar] [CrossRef] [PubMed]
  24. Hess, J.L.; Pyper, J.M.; Clements, J.E. Nucleotide sequence and transcriptional activity of the caprine arthritis-encephalitis virus long terminal repeat. J. Virol. 1986, 60, 385–393. [Google Scholar] [CrossRef] [Green Version]
  25. Glaria, I.; Reina, R.; Ramírez, H.; de Andrés, X.; Crespo, H.; Jauregui, P.; Salazar, E.; Luján, L.; Pérez, M.M.; Benavides, J.; et al. Visna/Maedi virus genetic characterization and serological diagnosis of infection in sheep from a neurological outbreak. Vet. Microbiol. 2012, 155, 137–146. [Google Scholar] [CrossRef]
  26. 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] [Green Version]
  27. Cardinaux, L.; Zahno, M.L.; Deubelbeiss, M.; Zanoni, R.; Vogt, H.R.; Bertoni, G. Virological and phylogenetic characterization of attenuated small ruminant lentivirus isolates eluding eluding efficient serological detection. Vet. Microbiol. 2013, 162, 572–581. [Google Scholar] [CrossRef]
  28. 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]
  29. 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, e193892. [Google Scholar] [CrossRef] [Green Version]
  30. Olech, M.; Murawski, M.; Kuźmak, J. Molecular analysis of small-ruminant lentiviruses in Polish flocks reveals the existence of a novel subtype in sheep. Arch. Virol. 2019, 164, 1193–1198. [Google Scholar] [CrossRef] [Green Version]
  31. 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] [PubMed]
  32. 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]
  33. 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] [PubMed] [Green Version]
  34. Gayo, E.; Cuteri, V.; Polledo, L.; Rossi, G.; García Marín, J.F.; Preziuso, S. Genetic Characterization and Phylogenetic Analysis of Small Ruminant Lentiviruses Detected in Spanish Assaf Sheep with Different Mammary Lesions. Viruses 2018, 10, 315. [Google Scholar] [CrossRef] [Green Version]
  35. Hötzel, I.; Kumpula-McWhirter, N.; Cheevers, W.P. Rapid evolution of two discrete regions of the caprine arthritis-encephalitis virus envelope surface glycoprotein during persistent infection. Virus Res. 2002, 84, 17–25. [Google Scholar] [CrossRef]
  36. Skraban, R.; Matthíasdóttir, S.; Torsteinsdóttir, S.; Agnarsdóttir, G.; Gudmundsson, B.; Georgsson, G.; Meloen, R.H.; Andrésson, Ó.S.; Staskus, K.A.; Thormar, H.; et al. Naturally occurring mutations within 39 amino acids in the envelope glycoprotein of maedi-visna virus alter the neutralization phenotype. J. Virol. 1999, 73, 8064–8072. [Google Scholar] [CrossRef] [Green Version]
  37. González Méndez, A.S.; Cerón Téllez, F.; Tórtora Pérez, J.L.; Martínez Rodríguez, H.A.; García Flores, M.M.; Ramírez Álvarez, H. Signature patterns in region V4 of small ruminant lentivirus surface protein in sheep and goats. Virus Res. 2020, 280, 197900. [Google Scholar] [CrossRef]
  38. Ramírez, H.; Reina, R.; Amorena, B.; de Andrés, D.; Martínez, H.A. Small ruminant lentiviruses: Genetic variability, tropism and diagnosis. Viruses 2013, 5, 1175–1207. [Google Scholar] [CrossRef] [Green Version]
  39. Maury, W.; Thompson, R.J.; Jones, Q.; Bradley, S.; Denke, T.; Baccam, P.; Smazik, M.; Oaks, J.L. Evolution of the equine infectious anemia virus long terminal repeat during the alteration of cell tropism. J. Virol. 2005, 79, 5653–5664. [Google Scholar] [CrossRef] [Green Version]
  40. Murphy, B.; Hillman, C.; Castillo, D.; Vapniarsky, N.; Rowe, J. The presence or absence of the gamma-activated site determines IFN gamma-mediated transcriptional activation in CAEV promoters cloned from the mammary gland and joint synovium of a single CAEV-infected goat. Virus Res. 2012, 163, 537–545. [Google Scholar] [CrossRef]
  41. Angelopoulou, K.; Poutahidis, T.; Brellou, G.D.; Greenland, T.; Vlemmas, I. A deletion in the R region of long terminal repeats in small ruminant lentiviruses is associated with decreased pathology in the lung. Vet. J. 2008, 175, 346–355. [Google Scholar] [CrossRef] [PubMed]
  42. Blatti-Cardinaux, L.; Sanjosé, L.; Zahno, M.L.; Zanoni, R.; Reina, R.; Bertoni, G. Detailed analysis of the promoter activity of an attenuated lentivirus. J. Gen. Virol. 2016, 97, 1699–1708. [Google Scholar] [CrossRef] [PubMed]
  43. Deubelbeiss, M.; Blatti-Cardinaux, L.; Zahno, M.L.; Zanoni, R.; Vogt, H.R.; Posthaus, H.; Bertoni, G. Characterization of small ruminant lentivirus A4 subtype isolates and assessment of their pathogenic potential in naturally infected goats. Virol. J. 2014, 11, 65. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Grego, E.; Reina, R.; Lanfredini, S.; Tursi, M.; Favole, A.; Profiti, M.; Lungu, M.M.; Perona, G.; Gay, L.; Stella, M.C.; et al. Viral load, tissue distribution and histopathological lesions in goats naturally and experimentally infected with the Small Ruminant Lentivirus Genotype E (subtype E1 Roccaverano strain). Res. Vet. Sci. 2018, 118, 107–114. [Google Scholar] [CrossRef]
  45. Pisoni, G.; Bertoni, G.; Manarolla, G.; Vogt, H.R.; Scaccabarozzi, L.; Locatelli, C.; Moroni, P. Genetic analysis of small ruminant lentiviruses following lactogenic transmission. Virology 2010, 407, 91–99. [Google Scholar] [CrossRef] [Green Version]
  46. Mankowski, J.L.; Flaherty, M.T.; Spelman, J.P.; Hauer, D.A.; Didier, P.J.; Amedee, A.M.; Murphey-Corb, M.; Kirstein, L.M.; Muñoz, A.; Clements, J.E.; et al. Pathogenesis of simian immunodeficiency virus encephalitis: Viral determinants of neurovirulence. J. Virol. 1997, 71, 6055–6060. [Google Scholar] [CrossRef] [Green Version]
  47. Pérez, M.; Biescas, E.; Reina, R.; Glaria, I.; Marín, B.; Marquina, A.; Salazar, E.; Álvarez, N.; de Andrés, D.; Fantova, E.; et al. Small ruminant lentivirus-induced arthritis: Clinicopathologic findings in sheep infected by a highly replicative SRLV B2 genotype. Vet. Pathol. 2015, 52, 132–139. [Google Scholar] [CrossRef] [Green Version]
  48. Larruskain, A.; Jugo, B.M. Retroviral infections in sheep and goats: Small ruminant lentiviruses and host interaction. Viruses 2013, 5, 2043–2061. [Google Scholar] [CrossRef] [Green Version]
  49. Herrmann-Hoesing, L.M.; White, S.N.; Mousel, M.R.; Lewis, G.S.; Knowles, D.P. Ovine progressive pneumonia provirus levels associate with breed and Ovar-DRB1. Immunogenetics 2008, 60, 749–758. [Google Scholar] [CrossRef]
  50. Larruskain, A.; Bernales, I.; Luján, L.; de Andrés, D.; Amorena, B.; Jugo, B.M. Expression analysis of 13 ovine immune response candidate genes in Visna/Maedi disease progression. Comp. Immunol. Microbiol. Infect. Dis. 2013, 36, 405–413. [Google Scholar] [CrossRef] [Green Version]
  51. White, S.N.; Knowles, D.P. Expanding possibilities for intervention against small ruminant lentiviruses through genetic markerassisted selective breeding. Viruses 2013, 5, 1466–1499. [Google Scholar] [CrossRef] [PubMed]
  52. Cecchi, F.; Dadousis, C.; Bozzi, R.; Fratini, F.; Russo, C.; Bandecchi, P.; Cantile, C.; Mazzei, M. Genome scan for the possibility of identifying candidate resistance genes for goat lentiviral infections in the Italian Garfagnina goat breed. Trop. Anim. Health Prod. 2019, 51, 729–733. [Google Scholar] [CrossRef] [PubMed]
  53. White, S.N.; Mousel, M.R.; Herrmann-Hoesing, L.M.; Reynolds, J.O.; Leymaster, K.A.; Neibergs, H.L.; Lewis, G.S.; Knowles, D.P. Genome-wide association identifies multiple genomic regions associated with susceptibility to and control of ovine lentivirus. PLoS ONE 2012, 7, e47829. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  54. Colussi, S.; Desiato, R.; Beltramo, C.; Peletto, S.; Modesto, P.; Maniaci, M.G.; Campia, V.; Quasso, A.; Rosati, S.; Bertolotti, L.; et al. A single nucleotide variant in the promoter region of the CCR5 gene increases susceptibility to arthritis encephalitis virus in goats. BMC Vet. Res. 2019, 15, 230. [Google Scholar] [CrossRef] [PubMed]
  55. Molaee, V.; Eltanany, M.; Lühken, G. First survey on association of TMEM154 and CCR5 variants with serological maedi-visna status of sheep in German flocks. Vet. Res. 2018, 49, 36. [Google Scholar] [CrossRef] [Green Version]
  56. Molaee, V.; Otarod, V.; Abdollahi, D.; Lühken, G. Lentivirus Susceptibility in Iranian and German Sheep Assessed by Determination of TMEM154 E35K. Animals 2019, 9, 685. [Google Scholar] [CrossRef] [Green Version]
  57. Olech, M.; Ropka-Molik, K.; Szmatoła, T.; Piórkowska, K.; Kuźmak, J. Single Nucleotide Polymorphisms in Genes Encoding Toll-Like Receptors 7 and 8 and Their Association with Proviral Load of SRLVs in Goats of Polish Carpathian Breed. Animals 2021, 11, 1908. [Google Scholar] [CrossRef]
Viruses 14 00735 g001 550
Figure 1. Example photomicrographs of tissue sections from the CAEV-infected goats. (A–C,E,G)—Sections stained with hematoxylin and eosin (HE); (D,F,H)—Immunohistochemical (IHC) labeling for CAEV p28 detection. (A) Brainstem, goat 5. Locally present mononuclear infiltrations around the blood vessel; (B) Choroid plexus, goat 5. Mononuclear infiltration within the choroid; (C) Lung, goat 7. Moderate interstitial pneumonia. The alveolar septa and peribronchial interstitium are widened by diffuse mononuclear infiltration and eosinophilic amorphic material; (D) Lung, goat 7. Positive CAEV p28 immunolabeling visible as brown particles in the cytoplasm of the bronchiolar epithelium (black arrows) and single macrophage-like cells (white arrows); (E) Carpal joint, goat 5. Hypertrophic synovial membrane-forming villous projections thickened by mononuclear infiltrates; (F) Carpal joint, goat 7. Positive CAEV p28 immunolabeling visible as brown particles in the cytoplasm of singular cells (arrows) within subsynovial mononuclear infiltrations; (G) Mammary gland, goat 5. Multifocal to diffuse mononuclear infiltrations around lactiferous ducts, blood vessels and within interstitial tissue; (H) Mammary gland, goat 5. Positive CAEV p28 immunolabeling in the cytoplasm of macrophage-like cells infiltrating the periacinar area; Scale bars = 50 µm.
Figure 1. Example photomicrographs of tissue sections from the CAEV-infected goats. (A–C,E,G)—Sections stained with hematoxylin and eosin (HE); (D,F,H)—Immunohistochemical (IHC) labeling for CAEV p28 detection. (A) Brainstem, goat 5. Locally present mononuclear infiltrations around the blood vessel; (B) Choroid plexus, goat 5. Mononuclear infiltration within the choroid; (C) Lung, goat 7. Moderate interstitial pneumonia. The alveolar septa and peribronchial interstitium are widened by diffuse mononuclear infiltration and eosinophilic amorphic material; (D) Lung, goat 7. Positive CAEV p28 immunolabeling visible as brown particles in the cytoplasm of the bronchiolar epithelium (black arrows) and single macrophage-like cells (white arrows); (E) Carpal joint, goat 5. Hypertrophic synovial membrane-forming villous projections thickened by mononuclear infiltrates; (F) Carpal joint, goat 7. Positive CAEV p28 immunolabeling visible as brown particles in the cytoplasm of singular cells (arrows) within subsynovial mononuclear infiltrations; (G) Mammary gland, goat 5. Multifocal to diffuse mononuclear infiltrations around lactiferous ducts, blood vessels and within interstitial tissue; (H) Mammary gland, goat 5. Positive CAEV p28 immunolabeling in the cytoplasm of macrophage-like cells infiltrating the periacinar area; Scale bars = 50 µm.
Viruses 14 00735 g001
Viruses 14 00735 g002 550
Figure 2. Maximum-likelihood phylogenetic tree based on the alignment of the CA fragment of the gag gene. Sequences from this study are labeled by black circles. Numbers at the branches indicate the percentage of bootstrap values obtained from 1000 replicates. Bootstrap values > 70% are shown. In the present study, the SRLV subtypes found by Colitti et al. [23] were renamed from A18 to A19 and from A19 to A20, and the subtypes found by Olech et al. [8] were renamed from A23 to A25 and from A24 to A26; goat 2—Blue, goat 3—Black, goat 4—Yellow, goat 5—Green, goat 6—Red, goat 7—Pink.
Figure 2. Maximum-likelihood phylogenetic tree based on the alignment of the CA fragment of the gag gene. Sequences from this study are labeled by black circles. Numbers at the branches indicate the percentage of bootstrap values obtained from 1000 replicates. Bootstrap values > 70% are shown. In the present study, the SRLV subtypes found by Colitti et al. [23] were renamed from A18 to A19 and from A19 to A20, and the subtypes found by Olech et al. [8] were renamed from A23 to A25 and from A24 to A26; goat 2—Blue, goat 3—Black, goat 4—Yellow, goat 5—Green, goat 6—Red, goat 7—Pink.
Viruses 14 00735 g002
Viruses 14 00735 g003 550
Figure 3. Maximum-likelihood phylogenetic tree based on the alignment of the V4V5 fragment of the env gene. Sequences from this study are labeled by black circles. Numbers at the branches indicate the percentage of bootstrap values obtained from 1000 replicates. Bootstrap values >70% are shown. In the present study, the SRLV subtypes found by Olech et al. [8] were renamed from A23 to A25 and from A24 to A26.
Figure 3. Maximum-likelihood phylogenetic tree based on the alignment of the V4V5 fragment of the env gene. Sequences from this study are labeled by black circles. Numbers at the branches indicate the percentage of bootstrap values obtained from 1000 replicates. Bootstrap values >70% are shown. In the present study, the SRLV subtypes found by Olech et al. [8] were renamed from A23 to A25 and from A24 to A26.
Viruses 14 00735 g003
Viruses 14 00735 g004 550
Figure 4. Alignment (MUSCLE) of deduced amino acid sequences of the gag-p25 capsid protein of SRLVs obtained in this study and K1514 (GenBank accession number M60609) and Cork (GenBank accession number M33677) reference strains, which are MVV (group A) and CAEV (group B) prototype strains, respectively. Identical residues are indicated by dots (.). Two immunodominant epitopes and the major homology region (MHR) are within squares. PBMC—Peripheral blood mononuclear cells, L—Lung, SM—Synovial membrane, B—Brain, CHP—Choroid plexus, MG—Mammary gland.
Figure 4. Alignment (MUSCLE) of deduced amino acid sequences of the gag-p25 capsid protein of SRLVs obtained in this study and K1514 (GenBank accession number M60609) and Cork (GenBank accession number M33677) reference strains, which are MVV (group A) and CAEV (group B) prototype strains, respectively. Identical residues are indicated by dots (.). Two immunodominant epitopes and the major homology region (MHR) are within squares. PBMC—Peripheral blood mononuclear cells, L—Lung, SM—Synovial membrane, B—Brain, CHP—Choroid plexus, MG—Mammary gland.
Viruses 14 00735 g004
Viruses 14 00735 g005 550
Figure 5. Alignment of deduced amino acid sequences of an immunodominant epitope of the ENV protein and variable region V4 of SRLVs obtained in this study and reference strains. Deletions are indicated by a dash (-), and identical residues are indicated by dots (.). PBMC—Peripheral blood mononuclear cells, L—Lung, SM—Synovial membrane, B—Brain, CHP—Choroid plexus, MG—Mammary gland.
Figure 5. Alignment of deduced amino acid sequences of an immunodominant epitope of the ENV protein and variable region V4 of SRLVs obtained in this study and reference strains. Deletions are indicated by a dash (-), and identical residues are indicated by dots (.). PBMC—Peripheral blood mononuclear cells, L—Lung, SM—Synovial membrane, B—Brain, CHP—Choroid plexus, MG—Mammary gland.
Viruses 14 00735 g005
Viruses 14 00735 g006 550
Figure 6. Maximum-likelihood phylogenetic tree based on the alignment of the LTR fragment. Sequences from this study are labeled by black circles. Numbers at the branches indicate the percentage of bootstrap values obtained from 1000 replicates. Bootstrap values >70% are shown. PBMC—Peripheral blood mononuclear cells, L—Lung, SM—Synovial membrane, B—Brain, CHP—Choroid plexus, MG—Mammary gland, goat 2—Blue, goat 3—Black, goat 4—Yellow, goat 5—Green, goat 6—Red, goat 7—Pink.
Figure 6. Maximum-likelihood phylogenetic tree based on the alignment of the LTR fragment. Sequences from this study are labeled by black circles. Numbers at the branches indicate the percentage of bootstrap values obtained from 1000 replicates. Bootstrap values >70% are shown. PBMC—Peripheral blood mononuclear cells, L—Lung, SM—Synovial membrane, B—Brain, CHP—Choroid plexus, MG—Mammary gland, goat 2—Blue, goat 3—Black, goat 4—Yellow, goat 5—Green, goat 6—Red, goat 7—Pink.
Viruses 14 00735 g006
Viruses 14 00735 g007a 550Viruses 14 00735 g007b 550
Figure 7. Alignment of U3-R sequences of LTR region from Polish SRLV strains. Sequences are aligned against prototype strains K1514 and Cork representative of SRLV groups A and B, respectively. Dots indicate identity with Cork, and dashes represent gaps. Boundaries between U3, R and U5 are indicated by straight arrows. AP-1, AP-4, AML (vis) motifs, TATA box, GAS, TAS and polyadenylation signal (poly A) are marked by boxes. Gray boxes represent CAAAT sequences. PBMC—Peripheral blood mononuclear cells, L—Lung, SM—Synovial membrane, B—Brain, CHP—Choroid plexus, MG—Mammary gland.
Figure 7. Alignment of U3-R sequences of LTR region from Polish SRLV strains. Sequences are aligned against prototype strains K1514 and Cork representative of SRLV groups A and B, respectively. Dots indicate identity with Cork, and dashes represent gaps. Boundaries between U3, R and U5 are indicated by straight arrows. AP-1, AP-4, AML (vis) motifs, TATA box, GAS, TAS and polyadenylation signal (poly A) are marked by boxes. Gray boxes represent CAAAT sequences. PBMC—Peripheral blood mononuclear cells, L—Lung, SM—Synovial membrane, B—Brain, CHP—Choroid plexus, MG—Mammary gland.
Viruses 14 00735 g007aViruses 14 00735 g007b
Viruses 14 00735 g008a 550Viruses 14 00735 g008b 550
Figure 8. Quantification of SRLV proviral load in different organs and cells of goats analyzed in this study. PBMC—Peripheral blood mononuclear cells, L—Lung, SM—Synovial membrane, B—Brain, CHP—Choroid plexus, MG—Mammary gland. Mean value of SRLV proviral copy number is shown above the bar.
Figure 8. Quantification of SRLV proviral load in different organs and cells of goats analyzed in this study. PBMC—Peripheral blood mononuclear cells, L—Lung, SM—Synovial membrane, B—Brain, CHP—Choroid plexus, MG—Mammary gland. Mean value of SRLV proviral copy number is shown above the bar.
Viruses 14 00735 g008aViruses 14 00735 g008b
Table
Table 1. Inflammation intensity scored in histopathological sections from the examined goats.
Table 1. Inflammation intensity scored in histopathological sections from the examined goats.
JointsBrainsLungMammary Glands
Goat 23 *012
Goat 32 *012
Goat 42 *013
Goat 53 *212
Goat 62 *012
Goat 73 *222
* arthritis score according to the grading system by Cheevers et al. [21].
Table
Table 2. Mean pairwise nucleotide genetic distance between subtypes of genotype A (inter-genotype) based on the CA fragment of the gag gene. In the present study, the SRLV subtypes found by Colitti et al. [23] were renamed from A18 to A19 and from A19 to A20, and the subtypes found by Olech et al. [8] were renamed from A23 to A25 and from A24 to A26.
Table 2. Mean pairwise nucleotide genetic distance between subtypes of genotype A (inter-genotype) based on the CA fragment of the gag gene. In the present study, the SRLV subtypes found by Colitti et al. [23] were renamed from A18 to A19 and from A19 to A20, and the subtypes found by Olech et al. [8] were renamed from A23 to A25 and from A24 to A26.
A1A2A3A4A5A8A9A11A12A13A16A17A18A19A20A21A22A23A24A25A26A27
A1----------------------
A216.6---------------------
A313.811.5--------------------
A414.716.115.3-------------------
A515.714.711.915.8------------------
A816.217.115.317.317.2-----------------
A915.214.512.216.514.315.1----------------
A1115.914.614.618.415.915.913.6---------------
A1216.212.710.916.613.516.715.814.5--------------
A1316.010.811.914.913.817.516.014.713.3-------------
A1617.915.917.215.716.419.620.618.016.116.3------------
A1714.211.58.815.011.015.711.515.510.013.214.7-----------
A1819.311.914.218.316.220.417.517.114.412.315.913.6----------
A1914.813.210.914.512.314.34.612.814.814.118.110.515.9---------
A2019.313.713.515.713.516.916.716.013.314.115.214.514.914.8--------
A2115.813.812.618.115.418.016.113.512.912.516.913.714.515.015.9-------
A2224.420.619.724.218.624.022.522.720.522.623.320.221.022.218.821.3------
A2319.617.216.418.817.417.015.714.916.918.821.716.820.314.515.919.421.3-----
A2415.912.612.016.014.213.211.813.914.013.318.212.815.612.516.213.121.115.0----
A2515.111.112.315.614.016.214.015.312.711.915.812.612.213.311.614.621.017.514.4---
A2617.711.012.916.613.217.215.514.413.011.116.013.010.214.811.413.118.817.612.110.3--
A2713.211.311.916.013.017.414.713.212.410.014.212.412.813.614.012.221.718.214.211.010.7-
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Olech, M.; Kycko, A.; Kuźmak, J. Molecular Characterization of Small Ruminant Lentiviruses Isolated from Polish Goats with Arthritis. Viruses 2022, 14, 735. https://doi.org/10.3390/v14040735

AMA Style

Olech M, Kycko A, Kuźmak J. Molecular Characterization of Small Ruminant Lentiviruses Isolated from Polish Goats with Arthritis. Viruses. 2022; 14(4):735. https://doi.org/10.3390/v14040735

Chicago/Turabian Style

Olech, Monika, Anna Kycko, and Jacek Kuźmak. 2022. "Molecular Characterization of Small Ruminant Lentiviruses Isolated from Polish Goats with Arthritis" Viruses 14, no. 4: 735. https://doi.org/10.3390/v14040735

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