1. Introduction
Johne’s disease (JD) is a global chronic gastroenteric condition that affects ruminants, causing the dairy industry serious economic concern [
1,
2]. JD is triggered by
Mycobacterium avium subsp. paratuberculosis (MAP), an acid-fast bacterium that belongs to the Mycobacteriaceae family [
2]. Many studies have emphasized the association of MAP with other gastrointestinal conditions such as Crohn’s disease (CD) in humans due to the symptomatic similarities that CD has in common with JD [
3,
4,
5,
6]. Animals that are infected with this bacterium are predisposed to other infections and conditions such as lameness, mastitis, and pneumonia due to JD causing an incurable immunodeficiency [
2]. Johne’s disease has three main stages—silent, subclinical, and clinical. The symptom manifestation usually begins after a long latency period (around 2 to 10 years). Therefore, the infected animals remain as reservoirs for the rest of their lives or transmit the disease to healthy individuals, either by direct contact and shedding of the MAP into the environment (horizontal transmission) or by reproducing newborn lambs (vertical transmission) [
2,
7].
Different MAP strains may infect different ruminants [
1]. Sheep are susceptible to both S and C MAP types, whereas cattle are resistant to the S type [
8]. These types differ not only by genotypic characteristics, but also by phenotypic divergence in pigmentation, iron metabolism, cytokine induction, and transcriptional patterns, as shown in experiments performed on a human monocyte cell line model (THP-1) [
1]. MAP strains belonging to the S type are notorious for their fastidious and slow-growing profiles, making the detection of this bacterium by some routine laboratory methods more laborious [
1,
9].
In Sardinia, unpasteurized sheep’s milk is usually used for cheesemaking purposes, which raises concerns over the consumption of these products and increases the possibility of some MAP-relevant diseases occurring in people in this region. This underlines the importance of using an efficient and confirmatory monitoring assay that can recognize the infected animals prior to using their milk in the cheese production cycle.
MAP has been reported at the herd level by testing various sample types, such as tissue, feces, and milk [
10,
11]. As milk can feasibly be manipulated, while blood and fecal samples cannot, milk is considered a potent specimen for diagnosing diseases relevant to ruminants [
12,
13].
Interestingly, the bulk tank milk (BTM) test is one the most economical and efficient herd-level analyses for screening of JD, which provides a general overview of the health status of the animals in each herd, especially as the presence of asymptomatic animals can lead to underestimation of the sensitivity and specificity of the diagnostic assays [
14,
15]. Evidence of MAP was discovered in milk and feces cultured from both symptomatic and asymptomatic carriers [
16,
17].
Culturing is one of the most reliable and specific techniques commonly used as a gold-standard reference in MAP diagnostic studies of milk and fecal samples. However, it lacks sensitivity (30–50%), partly due to the prolonged incubation time of between 7 and 16 weeks on solid media [
9,
10,
18,
19]. In addition, the concentrations of MAP (especially viable ones) in milk and colostrum are naturally lower than in other samples, such as feces [
9]. Although pretreatment of milk samples could augment the load of the bacterium in such samples, some preparation steps, such as chemical decontamination and incubation with antimicrobial brews, sacrifice the sensitivity of culture-based assessments for more specificity [
9]. This highlights the importance of applying alternative gold-standard approaches that increase the sensitivity and specificity of the diagnostic assay simultaneously.
The advent of the quantitative polymerase chain reaction (qPCR) technique was a great revolution in the area of molecular diagnosis that led to remarkable improvements in the detection and quantification of MAP DNA in milk and other clinical samples. Among the advantages of the qPCR technique are its availability, sensitivity, reproducibility, rapidity, and cost-effectiveness. Until now, some MAP sequences such as F57, IS-MAP02, and hspX have been targeted for qPCR analysis; however, the insertion sequence IS900, due to its high sensitivity (12–18 copies in the total MAP genome) and specificity, has been adopted in MAP screening programs [
9,
20,
21,
22].
On the other hand, the enzyme-linked immunosorbent assay (ELISA) is commonly used to assess the presence of sera antibodies directed against MAP in infected animal serum [
10]. Although some studies have demonstrated that MAP-infected animals develop titers of antibodies against MAP in milk as well as sera [
23,
24], some studies believe that milk and serum ELISA act imprecisely before the clinical stage of JD [
11,
14].
However, the sensitivity and specificity of the diagnostic assay (i.e., qPCR, ELISA) could be affected by the type of gold standard assigned for the detection. In addition, the efficiency of a reference test is also influenced by some factors, such as age of the animals, stage of the ailment, and source of the specimen [
25]. In contrast to some worldwide studies that used only a single test, such as fecal culture, as the gold standard [
25], the following study aimed to compare the efficiency of various gold-standard reference models of serum ELISA (SELISA) and fecal PCR (FPCR) at the individual level, and milk qPCR (MqPCR) and milk ELISA (MELISA) at the herd level (in isolation and in combination), in order to evaluate the sensitivity and specificity of each milk assay (MqPCR and MELISA) in predicting real MAP-infected cases and prove whether or not milk is a top-of-the-line sample over other specimens to detect MAP in sheep flocks. Moreover, in this study, a high-throughput IS900-based qPCR assay was developed using a modified DNA extraction protocol from low-quantity milk samples (5 mL), and its efficiency was compared with milk ELISA in a longitudinal and low-scale seasonal study as well.
2. Materials and Methods
2.1. Bacterial Strains
Mycobacterium avium subsp. paratuberculosis strain 1515 (ATCC 43015) and Mycobacterium smegmatis strain MC2155 (ATCC 700084) were used as positive and negative controls through the study (RIVM, Bilthoven, The Netherlands). They were grown in Middlebrooke 7H9 broth (Sigma-Aldrich, Milan, Italy) supplemented with 10% Oleic Albumin Dextrose Catalase (OADC; Sigma-Aldrich, Milan, Italy) and 2 mg mycobactin J (Allied Monitor, Fayette, MO, USA) and incubated at 37 °C for 3 days to 4 weeks (depending on the mycobacterial species).
2.2. Sample Collection
Raw sheep milk is commonly used in the cheese production process in Sardinia. A total of 201 unpasteurized sheep milk samples, including 73 bulk tank milk (BTM) samples from 59 dairy herds and 128 individual milk samples, were used in this study. The 128 individual samples were certified by the Veterinary Department of the University of Sassari and collected from a MAP-infected flock (MIF) of Northern Sardinia (Italy) (
Table 1). BTM sampling started from July 2018 and continued in different intervals until July 2019. Of the 59 herds, 14 herds were selected for a longitudinal study and monitored twice by BTM analyses (qPCR and ELISA) on different dates. The samples were taken in a sterile condition, kept at 4 °C during transportation, divided into aliquots immediately after transferring to the diagnostic laboratory, and stored at −80 °C without using any preservatives.
2.3. DNA Extraction from BTMs and MIF Milk Samples (MIFMs)
A modified DNA extraction protocol was established based on minimizing the quantity of milk samples (5 mL) required for DNA extraction and adding some milk preparation steps prior to DNA extraction as follows: dilution with phosphate-buffered saline (PBS; pH = 7.4; up to 30 mL), treatment with 0.75% (
w/v) hexadecylpyridinium chloride (HPC; Sigma-Aldrich, Milano, Italy), and homogenization by ribolyzer using glass beads (3 mm). Briefly, milk samples were allowed to stand at room temperature (RT) for one hour. Five milliliters of each milk sample was diluted in 25 mL 1× PBS (pH = 7.4) and centrifuged at 5000×
g rpm and 4 °C for 30 min. Then, the whey phase was decanted, the cream and pellet layers were collected, and DNA was extracted from both layers. At the next step, the cream and pellet were resuspended in 6 mL of hexadecylpyridinium chloride (HPC: Sigma-Aldrich, Milan, Italy; 0.75%
w/v), vigorously agitated by vortex, and incubated at RT for 1 h. The aim of this step was to release any MAP cells that might be surrounded by a lipid or cream layer in the milk samples. Then, samples were centrifuged at 5000×
g rpm at 15 °C for 20 min and supernatant including cream residuals and HPC was discarded. Subsequently, pellets were resuspended in 1 mL 1× PBS (pH = 7.4), transferred into specific microtubes, and homogenized by ribolyzer using glass beads (Sigma-Aldrich, Milan, Italy; diameter of 3 mm) for 4 cycles of 45 s at 4 m/s. This step aimed to homogenize the samples and release any MAP cells that might be trapped in somatic cells [
26,
27]. Then, this suspension was centrifuged at 10,000×
g rpm and 4 °C for 15 min, and DNA was extracted from the pellet using an RTP Mycobacteria kit (Stratec kit, Stratec Molecular GmbH, Berlin, Germany). The RTP Mycobacteria protocol was modified as follows: NAC buffer (NAC buffer comes with the other RTP-kit reagents and it probably contains N-acetylcysteine that is used for reducing the sputum viscosity) was replaced with 1× PBS (pH = 7.4); incubation times at 95 °C (15 min) and 65 °C (10 min) were increased to 30 and 20 min, respectively; and the column centrifugation steps were extended by 2 min based on the density of samples.
2.4. Milk qPCR (MqPCR)
qPCR has been performed to quantify the presence of IS900 in milk samples. Two mycobacterial species of Mycobacterium avium subsp. paratuberculosis strain 1515 and Mycobacterium smegmatis strain MC2155 were used as the positive and negative controls, respectively, through the qPCR analysis. Autoclaved Milli-Q water and SYBR Select Master Mix (Thermofisher Scientific, Applied Biosystems, Milan, Italy) were also included as the second and third negative controls, respectively. The MAP-specific primers AV1 (5′-ATGTGGTTGCTGTGTTGGATGG-3′) and AV2 (5′-CCGCCGCAATCAACTCCAG-3′) (Sigma-Aldrich, Milan, Italy) were used for amplification of IS900. Two to five microliters of extracted DNA (20–100 ng) were added into each reaction consisting of 10 μL SYBR Select Master Mix, 0.4 μL (10 μM) of each of the primers AV1 and AV2, and water added up to the final volume of 20 μL. The amplification conditions were as follows: initial denaturation at 95 °C for 3 min, followed by 50 cycles of denaturation at 95 °C for 40 s, annealing at 68 °C for 40 s, and extension at 72 °C for 40 s, with a final extension at 72 °C for 5 min.
2.5. Milk ELISA (MELISA)
In order to assess the presence of antibodies against MAP in milk samples, an indirect commercial milk ELISA kit (IDEXX Laboratories, Westbrook, ME, USA) was purchased.
As a sample preparation step prior to ELISA, inhibitors such as lipids (cream) and somatic cells (pellet) were removed from milk samples. To do this, samples were centrifuged at 10,000× g and 4 °C for 2 min, and then the whey phase, the liquid between the cream and pellet, was aspirated into a new Eppendorf tube (1.5 mL) and stored at −28 °C.
MELISA was performed according to the manufacturer’s (IDEXX) instructions. Briefly, samples (BTMs-MIFMs) were gently inverted few times and diluted (1:2) in dilution buffer N.12,
Mycobacterium phlei extract (using an uncoated 96-well microplate), in order to reduce cross-reactivity between nonspecific mycobacterial antibodies (abs; other than MAP-abs) and MAP epitopes (coating peptide). The positive and negative controls were also included in the assessment and diluted (1:20) in the same buffer (N.12) as milk samples. Then, they were homogenized using an orbital shaker (100 rpm/5 min) and incubated at RT for 2 h. Afterwards, samples (100 μL) were transferred into the wells of epitope-coated plates and incubated at RT for around 1 h. This step was followed by a washing step (300 μL × 5 times), adding conjugate (1:100, diluted in dilution buffer N.1), incubation at RT for 30 min, a further washing step (300 μL × 3 times), and adding 100 μL substrate 3,3′,5,5′-Tetramethylbenzidine (TMB) to each well. Finally, the optical densities (ODs) were read at the wavelength of 450 nm, converted to sample-to-positive ratios (S/P%) and the results were interpreted accordingly. Samples with S/P% of below 20, between 20–30, and 30 and above were classified as negative, suspect, and positive, respectively.
3. Statistical Analysis
Statistical analysis was carried out and then authenticated using R and Graphpad Prism 8, respectively. A number of samples tested by MELISA or SELISA produced suspect results. This transformed the dataset features from binary to ternary and might affect the ease of analysis. In order to overcome this data discrepancy, the tests related to these samples were repeated. MELISA-suspect BTMs and MIFMs were subjected to another trial. Regarding the SELISA-suspect MIFMs, samples that had at least two negative test results among MqPCR, MELISA, and FPCR were considered negative.
In order to perform a receiver operating characteristic (ROC) curve test, three different binary-based reference models (0,1) were suggested for both the MAP-infected flock (MIF) and BTM groups. In the MIF group, a binary dataset (0,1) was created based on the results of both SELISA and FPCR (individuals with at least one positive test were considered positive) and adjusted as the first gold-standard model reflecting the positive/negative status of each animal [
12,
28]. SELISA and FPCR datasets were separately considered as the second and third reference models to present the MAP status of each animal. In the BTM group, since there was no information on the clinical status of each herd animal, three different reference models were suggested based on the results of MqPCR and MELISA. In the first approach, the results of MqPCR and MELISA were interpreted as a binary-format dataset (0,1) and adjusted as a gold standard demonstrating the clinical status of each herd (herds with only one milk assay (MqPCR or MELISA) considered positive) [
12]. In the second and third approaches, each of the MqPCR and MELISA datasets were separately converted into binary numbers (0,1) and used as a gold standard to test the efficiency of the other milk assay (MELISA or MqPCR). Then, ROC curve and area under the curve (AUC) analyses were carried out as described before, and the sensitivity and specificity of the reference model was assessed with different thresholds (cutoffs were defined according to the results of MqPCR or MELISA that were continuous variables). In another analysis, the level of dependency between different test results was measured by Chi-square analysis using Graphpad Prism 8. For this purpose in the MIF group, continuous variables (MqPCR-MELISA) were converted to categorical variables (binary: (0,1)) and together with SELISA and FPCR datasets (that were originally categorical variables) used for Chi-square analysis. Regarding BTMs, a Chi-square test was executed on the categorical format of the MqPCR and MELISA datasets.
5. Discussion
Our findings suggest that milk is a potent sample for screening paratuberculosis in sheep flocks. However, the sensitivity and specificity of milk tests, i.e., MqPCR and MELISA, might be affected by some factors, such as the type of milk samples (BTM or individual samples), the disease status of animals participating in the survey, and the selected gold standard for statistical analysis [
28,
29]. Our BTM analyses (MqPCR and MELISA) demonstrated that MAP positivity follows a seasonal pattern, which was proven through a longitudinal study on 14 sheep herds.
Our MIF-relevant studies showed that MELISA and SELISA had the highest levels of concordance among the tests (89.8%). Such agreements were previously recorded by two other studies about the evaluation of the kinetics of antibodies directed against MAP during the lactation period (R
2 = 0.5358) [
30] and the assessment of the efficiency of a multiplex bead-based immunoassay in the detection of MAP-immunogenic antigens in animals with JD (R
2 = 0.572 to 0.756) [
31].
We compared the efficiency of different gold-standard models to reach the highest sensitivity and specificity in detection of the true incidence rate of MAP in milk samples. We diagnosed that SELISA together with FPCR is one of the best practical reference models in individual-level milk assessments. However, each of them as an independent gold standard favors a specific milk assay (MqPCR or MELISA). The results suggest remarkable sensitivity and specificity of MELISA when the binary reference models were based on SELISA in isolation (SP: 0.94, SN: 0.75;
p < 0.0001) and in combination with FPCR (SP: 0.95, SN: 0.61;
p < 0.0001). In comparison, MqPCR had its highest specificity (SP: 0.87; SN: 0.4;
p = 0.047) and sensitivity (SP: 0.84, SN: 0.47;
p = 0.044) when the gold standards were based on SELISA + FPCR and FPCR, respectively. However, when FPCR was the standard model, the specificity of MELISA dropped slightly (SP: 0.72, SN: 0.67;
p = 0.129). Our result is in accordance with another study that was conducted on BTMs from 21 dairy sheep flocks, evaluating the SP and SN of a modified milk ELISA in the detection of antibodies directed against MAP. A recent study that applied each of fecal PCR and serum ELISA as an independent reference test depicted that MELISA introduced a notable SP and SN when the reference test was based on SELISA (SP: 100%, SN: 72.7%) compared to when FPCR was the gold standard (SP: 46.7%, SN: 83.3%) [
29]. Interestingly, in another comparative work on the evaluation of the efficiency of a high-yield fecal qPCR (YHDEqPCR) assay in the detection of MAP, the poorest level of agreement was seen between YHDEqPCR and milk ELISA (1%;
p = 0.739), although the correlation between YHDEqPCR and milk qPCR was also insignificant (36%;
p = 0.591) [
14].
At the herd level (BTM), we developed a binary gold standard from both MqPCR and MELISA results that remarkably enhanced the sensitivity and specificity of each milk assay, i.e., MqPCR or MELISA. In order to assess whether or not this gold standard favors a specific milk assay, each of MqPCR and MELISA was used as an independent gold standard to evaluate the efficiency of the other diagnostic assay (MELISA or MqPCR). MqPCR offered a significantly greater sensitivity and specificity (SP: 1, SN: 0.95;
p < 0.0001) than MELISA (SP: 0.89, SN: 0.59;
p = 0.00035) when both MqPCR and MELISA were the reference models. We reached a similar SP but higher SN in MqPCR analysis than another BTM-MAP-detection study conducted on 21 sheep flocks under the condition that any of serum ELISA, milk ELISA, fecal PCR, and fecal culture was assigned as the gold standard (SP: 100%; SN: 25%) [
29].
We found that MqPCR as a reference test induced a remarkable specificity (SP: 0.95, SN: 0.34; p = 0.037) to MELISA compared to when MELISA was the reference test and MqPCR was the diagnostic assay (SP: 0.79, SN: 0.5; p = 0.098).
Our results revealed that the incidence of MAP and antibodies against it were lower in the MIF level than the herd groups (BTM) [
29,
32,
33,
34]. MqPCR and MELISA results were concordant with each other by 75% and 50% in MIF and BTM levels, respectively. Accordingly, the positivity rates of MqPCR and MELISA were estimated to be 19.53–21.09% in the MIF group and 72.6–30.14% in the BTM group. Thus, 47% of BTM and 11.7% of MIFM cases were detected as being positive by MqPCR without showing a sufficient antibody titer for positivity by MELISA. The significant disagreement ratio (47%) between MqPCR and MELISA in the BTM group was probably due to the high proportion of intermittent MAP shedders in each herd, lack of environmental hygiene practices that inhibit the dissemination of MAP through the environment, the inequity of the animals that participated in the survey in terms of immune status, stages of disease, age, and genetic predisposition [
16,
28,
35,
36]. Furthermore, MqPCR IS900 is one of the most sensitive MAP-detection assays and can detect the lowest concentration of MAP in the MAP-shedder animals in the initial stages of JD. As an American survey on the assessment of the sensitivity of various MAP-detection approaches (fecal–milk culture and qPCR) demonstrated that qPCR IS900 can distinguish MAP-shedder cows via milk, colostrum, and feces more sensitively than the culture-based assays [
9].
The results of BTM analyses (MqPCR and MELISA) on 73 samples from 59 herds showed an association between MAP-positivity rhythm and seasonal changes. We found that MqPCR and MELISA positivity have different seasonal patterns, in which the dominant number of positive MqPCR cases (45.21%) occurred in the spring. However, the highest number of MELISA-positive subjects (17.81%) was seen in the summer. This suggests that the number of MAP shedders increases in the spring due to seasonal breeding, possible sexual transmission of MAP, and the presence of animals shedding MAP via feces, which all corresponded to higher percentages of MqPCR-positive cases in the spring. However, humoral immunity due to MAP positivity was notable in the summer (MELISA). As a study on the detection of MAP in reproductive tissue showed that horizontal sex transmission increases the risk of MAP dissemination in infected rams, even though the lesions related to Johne’s disease are not developed in the reproductive tissues, MAP can be disseminated to the reproductive tissues [
37,
38,
39,
40,
41,
42]. Also, an Italian study on MAP seroprevalence in dairy sheep flocks evidenced that the rate of MAP seropositivity was enhanced in the spring and autumn (
p < 0.071) [
43].
However, our longitudinal study on 14 herds (BTMs) showed that the intensity of MAP-positivity, corresponding to lower TC values in MqPCR (ranged from 15 to 34 cycles), was more noticeable in the summer. A similar Canadian study on detecting MAP in milk and colostrum detected a higher proportion of MAP shedders in the summer by milk and colostrum qPCR than fecal and milk cultures [
9]. Our longitudinal study also determined an increase in the number of positive MELISA cases in the summer. Even though some cases were MELISA-negative both in the spring and summer, the antibody titers against MAP had an upward trend in the summer and changes were remarkable in some cases (S/P% increased by between 1.3% and 60.4%). This finding is in accordance with the results of another study about the association of season of sampling with change of MAP antibody titers in BTM. This study also confirmed that the MAP antibody titer follows a seasonal pattern in milk samples: peaking in the summer and drastically dropping in the winter [
23].
Our findings are limited by some factors, such as sample size, lack of communication with the herd’s owners regarding the health status of animals that participated in this survey and hygiene practices, and limitations in following up all herds in each season.
Our experiments on the efficiency of different reference standard models suggest that the type of gold-standard test remarkably affects the sensitivity and specificity of both MqPCR and MELISA tests. In MIF study, SELISA in isolation and together with FPCR introduced the best sensitivity and specificity to MqPCR and MELISA tests. Furthermore, MELISA and SELISA results were concordant in 75% of cases (p < 0.0001), and this proved that MELISA could be used as a predicting test in the MAP screening programs at the MIF level. Of course, this conclusion directly depends on some factors, such as the age, immune status, stage of disease and so forth of the animals that participated in the survey.
Our herd-level studies demonstrated that the results of MqPCR and MELISA together could form a practical binary reference model (gold standard) under the condition that there is no information on the health status of the herd’s animals. The results of both milk assays (MqPCR and MELISA) could represent how much potential each herd has to disseminate MAP through the environment and transmit the disease to other herds. The BTM-level studies depicted that MqPCR is more sensitive and specific (p < 0.0001) than MELISA (p = 0.00035). However, the agreement between the results of both milk assays was statistically significant when the p-value was adjusted to 0.1 (p = 0.083); this was equal to 50% concordance. This suggests that a high proportion of animals in each herd were MAP shedders.
Our studies on 73 sheep herds indicated that MAP positivity follows a seasonal pattern. The number of positive MqPCR cases (45.21%) peaked in the spring, while the trend of MELISA positivity was notable (17.81%) in the summer. However, the result of the longitudinal study on 14 herds showed that the lower TC values, corresponding to stronger positivity, are most obvious in the summer. This may suggest that a large number of MqPCR-positive cases in the spring could be MAP shedders.
The results obtained in this study highlight the risk of the transmission of this pathogenic mycobacterium to the community. A higher exposition to MAP, which causes a persistent infection in its host [
44], may lead to trigger not only inflammatory bowel diseases including Crohn’s disease [
44,
45], but also different autoimmune diseases, as it has been previously associated with in several studies [
46,
47,
48,
49,
50,
51,
52,
53,
54,
55]. This should motivate reflection by legislators and encourage the adoption of the best policies to lower the leakage of this bacterium into communities.