Multi-Platform Detection of Small Ruminant Lentivirus Antibodies and Provirus as Biomarkers of Production Losses

Small ruminant lentiviruses (SRLVs) are endemic in most areas of Europe, causing a chronic infection and a multisystemic disease affecting the udder, carpal joints, lungs, and central nervous system. Due to the lack of treatments and protective vaccination strategies, infection control is focused on the identification of infected animals through serological or molecular techniques. However, antigenic and genetic heterogeneity of SRLVs represent a clear drawback for diagnosis. Infected animals may present lower animal production parameters such as birth weight or milk production and quality, depending on productive systems considered and, likely, to the diagnostic method applied. In this study, four sheep flocks dedicated to dairy or meat production were evaluated using three different ELISA and two PCR strategies to classify animal population according to SRLV infection status. Productive parameters were recorded along one whole lactation or reproductive period and compared between positive and negative animals. SRLV was present in 19% of the total population, being unequally distributed in the different flocks. Less than half of the infected animals were detected by a single diagnostic method, highlighting the importance of combining different diagnostic techniques. Statistical analysis employing animal classification using all the diagnostic methods associated lambing size, lamb weight at birth, and daily weight gain with SRLV infection status in meat flocks. Milk production, somatic cell count, fat, and protein content in the milk were associated with SRLV infection in dairy flocks, to a greater extent in the flock showing higher seroprevalence. A multi-platform SRLV diagnostic strategy was useful for ensuring correct animal classification, thus validating downstream studies investigating production traits.

Small ruminant lentiviruses (SRLVs) are endemic in most areas of Europe, causing a chronic infection and a multisystemic disease affecting the udder, carpal joints, lungs, and central nervous system. Due to the lack of treatments and protective vaccination strategies, infection control is focused on the identification of infected animals through serological or molecular techniques. However, antigenic and genetic heterogeneity of SRLVs represent a clear drawback for diagnosis. Infected animals may present lower animal production parameters such as birth weight or milk production and quality, depending on productive systems considered and, likely, to the diagnostic method applied. In this study, four sheep flocks dedicated to dairy or meat production were evaluated using three different ELISA and two PCR strategies to classify animal population according to SRLV infection status. Productive parameters were recorded along one whole lactation or reproductive period and compared between positive and negative animals. SRLV was present in 19% of the total population, being unequally distributed in the different flocks. Less than half of the infected animals were detected by a single diagnostic method, highlighting the importance of combining different diagnostic techniques. Statistical analysis employing animal classification using all the diagnostic methods associated lambing size, lamb weight at birth, and daily weight gain with SRLV infection status in meat flocks. Milk production, somatic cell count, fat, and protein content in the milk were associated with SRLV infection in dairy flocks, to a greater extent in the flock showing higher seroprevalence. A multi-platform SRLV diagnostic strategy was useful for ensuring correct animal classification, thus validating downstream studies investigating production traits.

INTRODUCTION
Small ruminant lentiviruses (SRLVs) cause chronic infection in sheep and goats that results in the development of a multisystemic disease that may affect animal production depending on a myriad of factors including breed susceptibility (1)(2)(3), virulence of circulating strains (4,5), or production systems (6). Antibody production in response to infection can be detected by serological tests, while integrated provirus in circulating monocytes can be targeted by primers in PCR strategies. Despite the important contribution of ELISA approaches in control programs established so far (7), a description of novel genotypes that enlarge antigenic heterogeneity within SRLV jeopardizes ELISA performance, leading to diagnostic failures (8)(9)(10)(11).
Quantification of productive losses due to SRLV infection remains controversial, while some studies claim for a role of SRLV infection in decreasing quantity and quality of animal productions in both dairy and meat farms (12)(13)(14), others have revealed no differences between seronegative and seropositive animals (15)(16)(17)(18)(19). Slow disease development is a key feature of lentiviral infections and is the main cause of the underestimated losses in terms of animal production. However, differences in the production system, breed resistance, flock management, and parameters evaluated (20) may also explain the different associations between SRLV infection and productive traits.
Furthermore, diagnostics sensitivity may have significantly influenced studies aiming at evaluating production losses derived from SRLV infection. After initial assessment of antibody production for the detection of infected animals by agar gel immunodiffusion (AGID), later studies have demonstrated that ELISA methods are more sensitive (21). Beyond this, PCR and histopathology strategies can further improve the detection of virus in seronegative animals, indicating a benefit for combining molecular, and serological strategies (22). Antigenic and genetic heterogeneity of SRLVs are on the basis of the serologic and PCR test drawbacks, respectively. Indeed, antibodies against circulating strains are better detected using homologous antigens (8,9,23), and primer design is critical when developing sensitive and specific PCRs (24)(25)(26).
In this study, four sheep flocks belonging to two different production systems (dairy and meat) and breeds (Raza Navarra and Latxa Navarra) were classified as SRLV infected or uninfected using three different commercial ELISA methods, a home-made PCR kit, and a commercial PCR kit. The different commercial ELISA methods globally detected a similar number of infected animals; however, the combination of the three methods identified a significantly greater infected population in all the flocks analyzed. Furthermore, commercial PCR was more sensitive than ELISA in some cases and clearly added value to SRLV diagnosis and animal classification. Different production parameters in meat and dairy flocks were recorded during one lactation or reproductive period and, after final classification according to the different tests, were negatively (birth weight and weaning weight) or positively [somatic cell count (SCC)] associated with SRLV infection.
These results highlight the importance of using a multiplatform strategy to detect the humoral response to, as well as the presence of, different circulating strains in order to unequivocally identify infected and uninfected individuals, thereby influencing downstream studies such as production losses estimation, accreditation schemes, or control programs.

Animals and Samples
A total of 1,497 animals belonging to four different flocks dedicated to dairy or meat production were used ( Table 1). Flocks A and B (meat flocks, focused on lamb production) and Flocks C and D (dairy flocks) were located in the north of Navarra (humid climate) and bred in a semi-extensive system with periods of housing, especially during lactation. In addition, Flock D was bred in an ecological production scheme. None of the studied animals presented clinical signs of SRLV disease.
Whole blood was obtained by jugular puncture and placed in EDTA-K 3+ tubes. After centrifugation, plasma samples were stored at −20 • C until use in ELISA. Buffy coats were washed, erythrocytes lysed, resuspended in PBS, and stored at −20 • C until DNA extraction.

ELISA Tests
Plasma samples were tested for the presence of SRLV antibodies with three commercial ELISA kits: ELISA#1 uses an Eradikit TM SRLV screening test (In3Diagnostic, Torino, Italia) that includes capsid recombinant proteins from different genotypes; ELISA#2 uses ELITEST TM MVV/CAEV (Hyphen Biomed, Neuville-sur-Oise, France) that uses recombinant a capsid protein and synthetic peptide of the TM region; and ELISA#3 uses INgezim Maedi screening TM (Ingenasa, Eurofins Technologies, Spain) that includes synthetic peptides from the Env protein. All tests were performed following manufacturers' instructions.
Data were analyzed by considering each ELISA individually and the combinations of diagnostic kits. "Total ELISA" results were built by summarizing all positive samples to at least one of the ELISAs tested. Diagnostic efficiency was determined for each of the ELISAs in samples analyzed by the three methods, in comparison with the Total ELISA result.

DNA Extraction and PCR
Genomic DNA was extracted from buffy coat samples with and E.Z.N.A. tissue/blood kit (OMEGA, Bio-Tek, Norcross, GA, USA). DNA was quantified (NanoDrop OneC, Thermo Scientific R , Waltham, MA, USA), and real-time PCRs were performed using 250 ng of DNA (AriaMx Real Time PCR System) using the commercial kit EXOone Maedi Visna CAEV oneMix (Exopol, Spain) following manufacturer's instructions.
A total of 341 samples randomly distributed across all the flocks were analyzed by Gag PCR using primers described elsewhere (27)(28)(29)(30). Amplicons from 28 positive reactions were sequenced for molecular characterization purposes.
An animal was finally considered as infected when at least one commercial ELISA test or one PCR method revealed a positive result ("Total Infected").

Productive Parameters
In meat production farms, the following factors were studied during one production period: lambing size, weaning date, birth and weaning weight, and daily weight gain (DWG).
In dairy production farms, total volume of produced milk, fat, and protein content, and SCC were parameters quantified in each monthly quality control during a whole lactation.

Statistical Analysis
A generalized linear model (GLM) was used to assess the effect of each diagnostic strategy on each productive parameter recorded in meat flocks, taking into account the effect of the flock and the age of the mother included as covariables. Statistical results were shown as mean and standard error (SE). Production parameters were analyzed in dairy flocks by Mann Whitney's test due to significant heterogeneity in terms of SRLV prevalence, animal management, and production values between the flocks. Medians and interquartile range (IR) are shown. The kappa coefficient was calculated to assess the agreement between tests using the methods of Cohen-Fleiss (weighted) and Pearson-Fleiss.
Statistical analyses were carried out with IBM SPSS Statistics 19.0 for Windows, and alpha error was set at 0.050.

Serodiagnosis
Serological diagnosis strongly depended on the flock and the test used. While Flocks A, B, and C presented low seroprevalence rates of up to 10%, Flock D showed a moderate rate of seropositive animals of around 30%.
Seroprevalence values in meat flocks ranged from 2.5 to 7.5% depending on the ELISA used, showing an average value of 5% with similar ELISA efficiencies ( Table 2).
However, seroprevalence in dairy flocks differed not only depending on the ELISA used but also depending on the flock.
While Flock C showed low seroprevalence values, up to 10% (range 1.9-10%), seropositive animals in Flock D ranged from 24 to 34.1% depending on the test used. Average values considering both flocks ranged from 14.9 to 22%. ELISA efficiencies differed between flocks and ELISA tests, ranging from 10 to 84% ( Table 2).
Overall, ELISA#1 showed better results in terms of reactivity, taking into account the whole population, followed by ELISA#3 and ELISA#2. However, when considering data within flocks, ELISA#1 detected a higher number of animals in Flocks A and C, whereas ELISA#2 and ELISA#3 were more accurate in detecting infected animals in Flocks D and B, respectively. Considering samples positive to any of the kits tested (Total ELISA), a single ELISA only detected from 19 to 57% of the seropositive animals in Flock A and from 22 to 54% in Flock B; and the range in Flock C was from 10 to 53%. Single ELISA reactivity in Flock D ranged from 61 to 84.1% in comparison with Total ELISA results ( Table 2).
Indeed, agreement between serological tests calculated in samples analyzed with the three methods was very good between ELISA#1 and ELISA#2 only in Flock D. Values in Flocks A and C were not calculated since no samples were found to be simultaneously positive to more than one ELISA. In addition, association was minimal or low in Flock B (data not shown).

PCR
PCR analysis included a commercial kit that employs real-time PCR to detect and quantify different genotypes of SRLV and a home-made PCR based on previous publications. Commercial PCR detected more positive animals than did the classical homemade procedure, and both detected a total percentage of positive animals of 0, 36.8, 3.3, and 38.5% in Flocks A, B, C, and D, respectively. These values were higher than those obtained by serological methods in all flocks, except for Flock A in which the number of samples analyzed by PCR was significantly lower. In addition, an average of 10% of the seronegative animals rendered PCR-positive reactions. Consistently, PCR and ELISA agreement was minimal ( Table 3).
Sequence characterization of the Gag region in Flocks B, C, and D showed a mixed virus population between genotypes A and B2, the latter being prevalent, with high similarity in respect to previous deposited sequences. Gag immunodominant region alignment revealed no significant changes compared to sequences described (data not shown).

Evaluation of Production in Meat Farms
Lambing size, birth body weight, weaning date, and weight gain per day were affected depending on animal classification using single or combined ELISA results as well as PCRs. Lambing size was not affected when animal classification was conducted with single ELISAs; however, a combination of the three ELISAs and Total ELISA data revealed a significant positive association. A combination of PCRs also associated SRLV infection status with higher lambing size. However, a combination of serological and molecular tests associated SRLV infection with lower lambing size ( Table 4).
Birth weight was also associated when animals were classified according to ELISA combinations but not when considering the Total Infected result. SRLV infection influenced weaning date when considering a combination of the serological tests or the Total Infected data. Finally, gain per day was not associated when animals were classified according to single considered methods; however, the Total Infected result showed significant association with lower DWG values ( Table 4).
Total ELISA-positive animals were further classified into PCR negative or positive, and production losses were evaluated ( Table 5). Lower birth weight, weaning date, weight, and DWG values were observed in double-positive animals.

Evaluation of Production in Dairy Farms
Classification of animals according to serological tests, individually, or collectively considered, revealed differences in SCCs between seropositive and seronegative animals ( Table 6). Animal classification with both PCRs also identified elevated SCC in infected animals. Considering Total Infected animals (ELISA and PCR), SCCs were elevated in positive animals to an extent of 62% compared to uninfected individuals (p < 0.01). Milk yield was reduced by 6% in infected animals classified by single or combined ELISA results, as well as by commercial PCR (p < 0.05). Fat and protein contents were also related to SRLV serodiagnosis, while fat percentage was higher in milk from infected animals, protein content was reduced ( Table 6).    Total ELISA-positive animals were again further classified into PCR-negative or PCR-positive animals ( Table 7). While milk yield and protein content were reduced in ELISA-and PCRpositive animals, milk fat was inversely affected. A tendency to higher SCC was also observed in double-positive animals.

DISCUSSION
Serological diagnosis is currently the best choice for SRLV detection in livestock. It has been widely applied in control programs but also in downstream studies evaluating production losses (21,31) or genetic susceptibility and resistance to lentiviral infection (32). However, serological methods may fail at detecting the whole infected population due to virus antigenic diversity (10,33,34) or to delayed seroconversion (35), encouraging the update of existing serological methods to new variants and challenging the development and evaluation of molecular methods.
In this study, we analyze an ovine population of ∼1,500 individuals by ELISA methods detecting antibodies that recognize different antigen preparations. Our results demonstrate that care should be taken when ELISA tests are considered individually, since the combination of tests is able to increase the detection of seropositive animals up to 50%. In addition, we included the evaluation of a recently developed commercial qPCR that showed better detection of infected animals when compared with individual ELISA results ( Table 3). These results demonstrate that the truly infected population cannot be assessed by using a single ELISA strategy or even when applying three different commercial ELISA tests, since about 10% of infected animals remained seronegative but provirus positive and detected by PCR. SRLV seroprevalence in animals varied substantially when applying just one of the herein studied ELISA methods ranging from 2 to 34% of positive animals. A combination with PCR results, either commercial or home-made, enriched the infected population in different percentages depending on the flock studied, further modifying animal classification.
Three of the analyzed flocks, two of the Rasa Navarra breed (Flocks A and B) intended for meat production, and one milk flock of the Latxa Navarra breed (Flock C) presented low seroprevalence values. However, the remaining milk flock from the Latxa Navarra breed presented a moderate seroprevalence of around 30% excluding breed or production system influence.
Antigenic heterogeneity of SRLVs seems to be at the basis for this relatively low individual sensitivity. Consistently, the antigenic preparations included in the different ELISAs used may account for this different performance, since antigenic spectrum of the circulating strains in the studied population is also a key point to consider (36). Preliminary genetic characterization of the SRLVs in the different flocks reveals the presence of a mixed virus population including strains from different genotypes and subtypes. This may explain the better performance of the ELISA test including the highest antigenic diversity.
In spite of including antigenic preparations from different genotypes, ELISA failed at detecting a proportion of infected animals that were evidenced by molecular techniques. Infection by divergent SRLV strains in these animals is unlikely since primers used in PCR were designed based on known genotypes. Instead, low or fluctuating antibody titers may account for this discrepancy between serological and molecular techniques (36,37). Indeed, absence of serological response has been described in the periparturient period (38) as well as a result of recent infections (39). Additionally, antibody response in infected animals to viral epitopes not included in ELISAs cannot be ruled out. Therefore, the previously suggested combination between ELISA and PCR to really achieve a "gold standard" (20, 21) is reinforced from these results.
The use of more than one diagnostic technique allowed the evaluation of different animal classifications according to single ELISA, combined Total ELISA results, and the combination of these with PCRs, resulting in a Total Infected classification. Proper classification enabled determination of the real effect SRLV infection had on production traits. Of note, the use of a single ELISA may represent the detection of roughly the 61-84% of the seropositive animals as shown in Flock D of moderate seroprevalence. Furthermore, in low-seroprevalence flocks, a single ELISA may only detect 10% of the seropositive animals. ELISA performance differed among flocks, while results from Flocks A, B and C showed low seroprevalence and ELISA efficiency; in Flock D, seroprevalence and efficiency reached 32 and 84%, respectively. Breed, age, production system, nutrition, or animal management cannot be argued as important factors, since Flock C shared these features. Instead, circulation of a more prototypic SRLV in Flock D could explain this better performance.
In meat flocks, SRLV influence on animal production was clearly evidenced when serological tests were combined between them or with molecular tests. Furthermore, different interpretations could be reached, taking into account the different animal classifications. Total Infected animals, including PCR and ELISA, showed lower lambing sizes and a trend to  higher birth body weight. Despite one variable possibly being related to the other, since higher lambing size implies lower birth weights (40), previous studies relate SRLV with lower birth weights (18) or rather did not find association (16,17,19) likely due to the low epidemiologic importance of natural in utero transmission (41). In contrast, serological methods associated SRLV seropositivity with lower birth body weight and with lambing size depending on the data analysis performed. While positive animals to the three ELISAs used (ELISA#1, ELISA#2, and ELISA#3) showed higher lambing size, ELISA combination (Total ELISA) associated lower lambing size with SRLV positivity. Since Total PCR results were in accordance with higher lambing size in positive animals, inclusion of ELISA falsepositive reactions in the Total ELISA and Total Infected groups may help to explain this discrepancy. Previous studies including one of the ELISAs used, reported specificity values ranging from 98.4 to 99.8% with respect to AGID (42), further supporting this hypothesis. Despite the very low seroprevalence observed in meat flocks to single ELISAs, a moderate presence of infected animals (∼30%) was evidenced by PCR. Thus, PCR analysis has improved the results presented here due to detection of incipient infections that may mask SRLV influence. Chronic infections and especially SRLV show a long asymptomatic period in which ewe's body condition may inadvertently diminish, likely determining a reduced nutrient transfer to the fetus (43,44). Sustained immune response in these infections may also alter the metabolism to a more catabolic profile, thereby reducing disposable input for the lamb. Actually, HIV infection has severe impact on pregnancy outcomes such as low birth weight and preterm delivery (45)(46)(47).
In dairy flocks, the application of single ELISA already identified higher SCC and fat content as well as lower milk yield and protein in milk from SRLV-seropositive sheep. A combination of ELISAs and PCR further confirmed this finding. Total Infected animals showed lower milk production (up to 3%) and elevated SCCs (60% increment). Augmented SCC has been already linked to SRLV infection due to epithelial cell desquamation derived from microscopic alterations in the mammary gland (48,49) and may represent lower milk quality and, beyond, economic losses to farmers due to penalties. In the absence of clinical signs, increased SCC could be related to systemic incipient lesions that may be present in up to 20% of infected animals (50). Interestingly, recent studies show that up to 90.9% of naturally SRLV-infected animals exhibit minimal to moderate lesions in the mammary gland, this prevalence being even higher in intensive milk-producing systems (22). Increased fat content in the milk could be the simple consequence of lower production (51). Decreased protein content was found in infected animals, further pointing out SRLV influence on milk production parameters.
Among ELISA-positive animals, the PCR-negative population showed lower production losses as compared to PCR-positive animals in meat and dairy flocks. Higher viral load implies higher PCR sensitivity as well as increased disease severity (52)(53)(54). These results suggest that antibodies revealed in ELISA may play a protective role, thereby reducing clinical signs and production losses. In agreement, the presence of antibodies against SRLV in milk may reduce proviral load detection in milk cells (55).
Interestingly, lower weight at weaning presented by lambs from seropositive ewes in meat farms could be explained by the lower milk production observed in infected sheep from dairy flocks. However, milk production parameters were not evaluated in meat flocks.
Exhaustive estimation of production losses derived from infections, especially those chronic, should be evaluated after proper infection status evaluation. The multi-platform strategy applied here to classify more than 1,000 animals into SRLV infected vs. uninfected enabled the analysis of different production parameters in meat and milk-oriented semiextensive production systems. Proper diagnosis was achieved when three different ELISA methods and two different PCRs were used. Both meat and dairy flocks showed diminished production parameters in infected animals, mainly affecting birth and weaning weights as well as milk production together with an increased number of somatic cells counts. These results highlight the crucial importance of proper SRLV infection status determination in sheep production studies and help to clarify previous colliding results obtained by other authors.

DATA AVAILABILITY STATEMENT
All datasets generated for this study are included in the article/supplementary material.

ETHICS STATEMENT
Ethical review and approval was not required for the animal study because Samples were obtained for institutional control campaings and later used in this study. Written informed consent was obtained from the owners for the participation of their animals in this study.

AUTHOR CONTRIBUTIONS
IE, RR, and LD were involved in processing of the samples and DNA extraction. IE, RD, and LD performed ELISA and PCR. AB, IE, and RD contributed to real-time PCR analysis. IE, RD, and ID performed statistical analysis of data. IE, RD, and RR drafted the manuscript. ID, DD, LL, IG, and RR revised the manuscript and wrote the final version. DD, LL, and RR obtained funds and coordinated and supervised the study. All co-authors approved the final version of the manuscript.

FUNDING
This research was funded by the Spanish Ministry of Science, Innovation and Universities, grant number RTI2018-096172-B-C31, and by the Government of Navarra and European Union ERDF program, grant number PDR (210160006). LD and IE were funded by Universidad Pública de Navarra. RR was supported by the Spanish Ministry of Science and Innovation Ramón y Cajal contract. Authors acknowledge APC funding by the CSIC Open Access Publication Support Initiative through its Unit of Information Resources for Research (URICI). RD was a Ph.D. student funded by the Department of Innovation, Research and University of Aragón.

ACKNOWLEDGMENTS
Authors acknowledge INTIA of Navarra for blood samples and breeders association of the Raza Navarra breed (ARANA) and breeders association of the Latxa Navarra breed (ASLANA) for production data. Ana Mariela Dodero (INTA, Salta, Argentina) is acknowledged for technical assistance in obtaining buffy coat. Authors acknowledge the reviewers for their contributions as their comments have considerably improved the manuscript.