MicroRNAs modulate immunological and inflammatory responses in Holstein cattle naturally infected with Mycobacterium avium subsp. paratuberculosis.

MicroRNAs (miRNAs) regulate the post-transcriptional expression of genes by binding to their target mRNAs. In this study, whole miRNA sequencing was used to compare the expression of miRNAs in ileocecal valve (ICV) and peripheral blood (PB) samples of cows with focal or diffuse paratuberculosis (PTB)-associated lesions in gut tissues versus (vs) control cows without lesions. Among the eight miRNAs differentially expressed in the PB samples from cows with diffuse lesions vs controls, three (miR-19a, miR-144, miR32) were also down-regulated in cows with diffuse vs focal lesions. In the ICV samples, we identified a total of 4, 5, and 18 miRNAs differentially expressed in cows with focal lesions vs controls, diffuse lesions vs controls, and diffuse vs focal lesions, respectively. The differential expression of five microRNAs (miR-19a, miR-144, miR-2425-3p, miR-139, miR-101) was confirmed by RT-qPCR. Next, mRNA target prediction was performed for each differentially expressed miRNA. A functional analysis using the predicted gene targets revealed a significant enrichment of the RNA polymerase and MAPK signaling pathways in the comparison of cows with focal vs no lesions and with diffuse vs focal lesions, respectively. The identified miRNAs could be used for the development of novel diagnostic and therapeutical tools for PTB control.

Proteomic Serum Profiling of Holstein Friesian Cows with Different Pathological Forms of Bovine Paratuberculosis Reveals Changes in the Acute-Phase Response and Lipid Metabolism.

The lack of sensitive diagnostic methods to detect Mycobacterium avium subsp. paratuberculosis (Map) subclinical infections has hindered the control of paratuberculosis (PTB). The serum proteomic profiles of naturally infected cows presenting focal and diffuse pathological forms of PTB and negative controls (n = 4 per group) were analyzed using TMT-6plex quantitative proteomics. Focal and diffuse are the most frequent pathological forms in subclinical and clinical stages of PTB, respectively. One (focal versus (vs.) control), eight (diffuse vs. control), and four (focal vs. diffuse) differentially abundant (DA) proteins (q-value < 0.05) were identified. Ingenuity pathway analysis of the DA proteins revealed changes in the acute-phase response and lipid metabolism. Six candidate biomarkers were selected for further validation by specific ELISA using serum from animals with focal, multifocal, and diffuse PTB-associated lesions (n = 108) and controls (n = 56). Overall, the trends of the serum expression levels of the selected proteins were consistent with the proteomic results. Alpha-1-acid glycoprotein (ORM1)-based ELISA, insulin-like growth factor-binding protein 2 (IGFBP2)-based ELISA, and the anti-Map ELISA had the best diagnostic performance for detection of animals with focal, multifocal, and diffuse lesions, respectively. Our findings identify potential biomarkers that improve diagnostic sensitivity of PTB and help to elucidate the mechanisms involved in PTB pathogenesis.

Summary-data based Mendelian randomization identifies gene expression regulatory polymorphisms associated with bovine paratuberculosis by modulation of the nuclear factor Kappa ß (NF-κß)-mediated inflammatory response.

Genome-wide association studies (GWAS) have identified host genetic variants associated with paratuberculosis (PTB) susceptibility. Most of the GWAS-identified SNPs are in non-coding regions. Connecting these non-coding variants and downstream affected genes is a challenge and, up to date, only a few functional mutations or expression quantitative loci (cis-eQTLs) associated with PTB susceptibility have been identified. In the current study, the associations between imputed whole-genome sequence genotypes and whole RNA-Sequencing data from peripheral blood (PB) and ileocecal valve (ICV) samples of Spanish Holstein cows (N = 16) were analyzed with TensorQTL. This approach allowed the identification of 88 and 37 cis-eQTLs regulating the expression levels of 90 and 37 genes in PB and ICV samples, respectively (False discorey rate, FDR ≤ 0.05). Next, we applied summary-based data Mendelian randomization (SMR) to integrate the cis-eQTL dataset with GWAS data obtained from a cohort of 813 culled cattle that were classified according to the presence or absence of PTB-associated histopathological lesions in gut tissues. After multiple testing corrections (FDR ≤ 0.05), we identified two novel cis-eQTLs affecting the expression of the early growth response factor 4 (EGR4) and the bovine neuroblastoma breakpoint family member 6-like protein isoform 2 (MGC134040) that showed pleiotropic associations with the presence of multifocal and diffuse lesions in gut tissues; P = 0.002 and P = 0.017, respectively. While EGR4 acts as a brake on T-cell proliferation and cytokine production through interaction with the nuclear factor Kappa ß (NF-κß), MGC134040 is a target gene of NF-κß. Our findings provide a better understanding of the genetic factors influencing PTB outcomes, confirm that the multifocal lesions are localized/confined lesions that have different underlying host genetics than the diffuse lesions, and highlight regulatory SNPs and regulated-gene targets to design future functional studies.

Efficacy of heat-inactivated Mycobacterium bovis vaccine delivered to European badgers (Meles meles) through edible bait.

Badgers (Meles meles) are a major tuberculosis (TB) reservoir in Europe, with the potential to transmit infection to cattle. Here we assessed whether a recently described oral tuberculosis vaccine based on heat-inactivated Mycobacterium bovis (HIMB), delivered as edible baits, can protect badgers from infection. Eight badgers were given individually five baits, each one consisting of a ball of peanut butter, natural peanut and oat flakes including a dose of the vaccine containing 5 × 107 colony-forming units. In parallel, a control group of seven badgers did not receive the vaccine. One month and a half later a second dose of the vaccine was offered to the vaccinated group. Ninety-four days after the second dose, all badgers were challenged with M. bovis (103 colony-forming units per animal) delivered endobronchially to the right middle lung lobe. Clinical, immunological, pathological and bacteriological variables were measured throughout the whole study to assess the efficacy of the vaccine. Two vaccinated animals showed high bacterial load of M. bovis and worsening of pathological lesions of TB. Conversely, the other six vaccinated animals showed slight improvement in bacterial load and pathology with respect to the control group. These results suggest that delivering the TB vaccine via food bait can partially protect wild badger populations, although vaccination can lead to either protection or tolerization, likely depending on the animal's immune status and general condition at the time of vaccination. Further optimization of the vaccination trial/strategy is needed to reduce the rate of tolerization, such as altering vaccine dose, number of doses, type of bait, use of adjuvants or route of administration.

Evaluation of a droplet digital PCR assay for quantification of Mycobacterium avium subsp. paratuberculosis DNA in whole-blood and fecal samples from MAP-infected Holstein cattle.

Bovine paratuberculosis (PTB) is an infectious disease that affects ruminants worldwide and is a burden on the dairy industry. PTB control measures include culling of Mycobacterium avium subsp. paratuberculosis (MAP)-infected animals from the herd and the enhancement of farm-biosecurity measures. Diagnostics tools for the direct detection of MAP are fecal real-time qPCR and bacteriological culture, the last one being considered the gold standard. However, both show limitations for detecting subclinical MAP-infected cattle with low bacterial load in feces and gut tissues. Droplet digital polymerase chain reaction (ddPCR) is a third-generation PCR method that shows high reproducibility for the quantification of low DNA copy numbers. The objective of this study was to design a ddPCR assay to detect and quantify a fragment of the F57 MAP-specific sequence in samples of naturally MAP-infected Holstein cattle. DNA was isolated from whole-blood and fecal samples from control cows with a negative ELISA and qPCR result (N = 75) and from cows with PTB-associated focal (N = 32), multifocal (N = 21), and diffuse lesions (N = 17) in gut tissues. After ddPCR, the DNA extracted from fecal samples of cows with diffuse lesions showed higher mean copies per microliter (13,791.2 copies/µl) than samples from cows with multifocal lesions (78.8 copies/µl), focal lesions (177.1 copies/µl) or control cows (4.8 copies/µl) (P ≤ 0.05). Significant differences in mean DNA copies/µl were also observed in the blood samples from cows with focal lesions (47.7 copies/µl) when compared with cows with multifocal and diffuse lesions; 18.1 and 12.4 copies/µl, respectively. Using a principal component analysis, the results of the fecal ddPCR clustered together with the results of a commercial ELISA for the specific detection of MAP antibodies, fecal and tissue qPCR, and bacteriological culture results. In contrast, blood ddPCR results clustered together with the results of an ELISA for the detection of a biomarker of subclinical PTB, the ABCA13 transporter. Blood ddPCR was the most sensitive tool (sensitivity 71%, specificity 100%) of all the quantitative methods used in the study for the detection of subclinical cows with focal lesions.

Use of ATP-Binding Cassette Subfamily A Member 13 (ABCA13) for Sensitive Detection of Focal Pathological Forms of Subclinical Bovine Paratuberculosis.

Bovine paratuberculosis (PTB) is a chronic enteritis caused by Mycobacterium avium subspecies paratuberculosis (Map) that causes a heavy economic impact worldwide. Map infected animals can remain asymptomatic for years while transmitting the mycobacteria to other members of the herd. Therefore, accurate detection of subclinically infected animals is crucial for disease control. In a previous RNA-Seq study, we identified several mRNAs that were overexpressed in whole blood of cows with different PTB-associated histological lesions compared with control animals without detected lesions. The proteins encoded by two of these mRNAs, ATP binding cassette subfamily A member 13 (ABCA13) and Matrix Metallopeptidase 8 (MMP8) were significantly overexpressed in whole blood of animals with focal histological lesions, the most frequent pathological form in the subclinical stages of the disease. In the current study, the potential of sensitive early diagnostic tools of commercial ELISAs, based on the detection of these two biomarkers, was evaluated in serum samples of 704 Holstein Friesian cows (566 infected animals and 138 control animals from PTB-free farms). For this evaluation, infected animals were classified into three groups, according to the type of histological lesions present in their gut tissues: focal (n = 447), multifocal (n = 59), and diffuse (n = 60). The ELISA based on the detection of ABCA13 was successfully validated showing good discriminatory power between animals with focal lesions and control animals (sensitivity 82.99% and specificity 80.43%). Conversely, the MMP8-based ELISA showed a poor discriminatory power between the different histological groups and non-infected controls. The ABCA13-based ELISA showed a higher diagnostic value (0.822) than the IDEXX ELISA (0.517), the fecal bacterial isolation (0.523) and the real-time PCR (0.531) for the detection of animals with focal lesions. Overall, our results indicate that this ABCA13 ELISA greatly improves the identification of subclinically infected animals with focal lesions that are undetectable using current diagnostic methods.

Bovine Intelectin 2 Expression as a Biomarker of Paratuberculosis Disease Progression.

Paratuberculosis (PTB), a chronic granulomatous enteritis caused by Mycobacterium avium subsp. paratuberculosis (MAP), is responsible for important economic losses in the dairy industry. Our previous RNA-sequencing (RNA-Seq) analysis showed that bovine intelectin 2 (ITLN2) precursor gene was overexpressed in ileocecal valve (ICV) samples of animals with focal (log2 fold-change = 10.6) and diffuse (log2 fold-change = 6.8) PTB-associated lesions compared to animals without lesions. This study analyzes the potential use of ITLN2, a protein that has been described as fundamental in the innate immune response to infections, as a biomarker of MAP infection. The presence of ITLN2 was investigated by quantitative immunohistochemical analysis of ICV samples of 20 Holstein Friesian cows showing focal (n = 5), multifocal (n = 5), diffuse (n = 5) and no histological lesions (n = 5). Significant differences were observed in the mean number of ITLN2 immunostained goblet and Paneth cells between the three histopathological types and the control. The number of immunolabelled cells was higher in the focal histopathological type (116.9 ± 113.9) followed by the multifocal (108.7 ± 140.5), diffuse (76.5 ± 97.8) and control types (41.0 ± 81.3). These results validate ITLN2 as a post-mortem biomarker of disease progression.

Spatial and Temporal Distribution of Mycobacterium tuberculosis Complex Infection in Eurasian Badger (Meles meles) and Cattle in Asturias, Spain.

The present work investigated the prevalence, spatial distribution, and temporal distribution of tuberculosis (TB) in free-ranging Eurasian badgers (Meles meles) and cattle in Asturias (Atlantic Spain) during a 13-year follow-up. The study objective was to assess the role of badgers as a TB reservoir for cattle and other sympatric wild species in the region. Between 2008 and 2020, 673 badgers (98 trapped and 575 killed in road traffic accidents) in Asturias were necropsied, and their tissue samples were cultured for the Mycobacterium tuberculosis complex (MTC) isolation. Serum samples were tested in an in-house indirect P22 ELISA to detect antibodies against the MTC. In parallel, data on MTC isolation and single intradermal tuberculin test results were extracted for cattle that were tested and culled as part of the Spanish National Program for the Eradication of Bovine TB. A total of 27/639 badgers (4.23%) were positive for MTC based on bacterial isolation, while 160/673 badgers (23.77%) were found to be positive with the P22 ELISA. The rate of seropositivity was higher among adult badgers than subadults. Badger TB status was spatially and temporally associated with cattle TB status. Our results cannot determine the direction of possible interspecies transmission, but they are consistent with the idea that the two hosts may exert infection pressure on each other. This study highlights the importance of the wildlife monitoring of infection and disease during epidemiological interventions in order to optimize outcomes.

Identification of loci associated with susceptibility to bovine paratuberculosis and with the dysregulation of the MECOM, eEF1A2, and U1 spliceosomal RNA expression.

Although genome-wide association studies have identified single nucleotide polymorphisms (SNPs) associated with the susceptibility to Mycobacterium avium subsp. paratuberculosis (MAP) infection, only a few functional mutations for bovine paratuberculosis (PTB) have been characterized. Expression quantitative trait loci (eQTLs) are genetic variants typically located in gene regulatory regions that alter gene expression in an allele-specific manner. eQTLs can be considered as functional links between genomic variants, gene expression, and ultimately phenotype. In the current study, peripheral blood (PB) and ileocecal valve (ICV) gene expression was quantified by RNA-Seq from fourteen Holstein cattle with no lesions and with PTB-associated histopathological lesions in gut tissues. Genotypes were generated from the Illumina LD EuroG10K BeadChip. The associations between gene expression levels (normalized read counts) and genetic variants were analyzed by a linear regression analysis using R Matrix eQTL 2.2. This approach allowed the identification of 192 and 48 cis-eQTLs associated with the expression of 145 and 43 genes in the PB and ICV samples, respectively. To investigate potential relationships between these cis-eQTLs and MAP infection, a case-control study was performed using the genotypes for all the identified cis-eQTLs and phenotypical data (histopathology, ELISA for MAP-antibodies detection, tissue PCR, and bacteriological culture) of 986 culled cows. Our results suggested that the heterozygous genotype in the cis-eQTL-rs43744169 (T/C) was associated with the up-regulation of the MDS1 and EVI1 complex (MECOM) expression, with positive ELISA, PCR, and bacteriological culture results, and with increased risk of progression to clinical PTB. As supporting evidence, the presence of the minor allele was associated with higher MECOM levels in plasma samples from infected cows and with increased MAP survival in an ex-vivo macrophage killing assay. Moreover, the presence of the two minor alleles in the cis-eQTL-rs110345285 (C/C) was associated with the dysregulation of the eukaryotic elongation factor 1-α2 (eEF1A2) expression and with increased ELISA (OD) values. Finally, the presence of the minor allele in the cis-eQTL rs109859270 (C/T) was associated with the up-regulation of the U1 spliceosomal RNA expression and with an increased risk of progression to clinical PTB. The introduction of these novel functional variants into marker-assisted breeding programs is expected to have a relevant effect on PTB control.

Detection of latent forms of Mycobacterium avium subsp. paratuberculosis infection using host biomarker-based ELISAs greatly improves paratuberculosis diagnostic sensitivity.

Bovine paratuberculosis (PTB) is a chronic granulomatous enteritis, caused by Mycobacterium avium subsp. paratuberculosis (MAP), responsible for important economic losses in the dairy industry. Current diagnostic methods have low sensitivities for detection of latent forms of MAP infection, defined by focal granulomatous lesions and scarce humoral response or MAP presence. In contrast, patent infections correspond to multifocal and diffuse types of enteritis where there is increased antibody production, and substantial mycobacterial load. Our previous RNA-Seq analysis allowed the selection of five candidate biomarkers overexpressed in peripheral blood of MAP infected Holstein cows with focal (ABCA13 and MMP8) and diffuse (FAM84A, SPARC and DES) lesions vs. control animals with no detectable PTB-associated lesions in intestine and regional lymph nodes. The aim of the current study was to assess the PTB diagnostic potential of commercial ELISAs designed for the specific detection of these biomarkers. The ability of these ELISAs to identify animals with latent and/or patent forms of MAP infection was investigated using serum from naturally infected cattle (n = 88) and non-infected control animals (n = 67). ROC analysis revealed that the ABCA13-based ELISA showed the highest diagnostic accuracy for the detection of infected animals with focal lesions (AUC 0.837, sensitivity 79.25% and specificity 88.06%) and with any type of histological lesion (AUC 0.793, sensitivity 69.41% and specificity 86.57%) improving on the diagnostic performance of the popular IDEXX ELISA and other conventional diagnostic methods. SPARC and MMP8 showed the highest diagnostic accuracy for the detection of animals with multifocal (AUC 0.852) and diffuse lesions (AUC 0.831), respectively. In conclusion, our results suggest that quantification of ABCA13, SPARC and MMP8 by ELISA has the potential for implementation as a diagnostic tool to reliably identify MAP infection, greatly improving early detection of MAP latent infections when antibody responses and fecal shedding are undetectable using conventional diagnostic methods.

Correlations between single nucleotide polymorphisms in bovine CD209, SLC11A1, SP110 and TLR2 genes and estimated breeding values for several traits in Spanish Holstein cattle.

Little is known about the correlations between the genetic susceptibility/resistance to Mycobacterium avium subsp. paratuberculosis (Map) infection and the estimated breeding values for type, production and functional traits. Previously, we identified 70 combinations of five single nucleotide polymorphisms (SNPs) in four bovine innate immune genes (SLC11A1, SP110, TLR2, CD209) that are associated with the genetic risk of paratuberculosis (PTB) or Johne's disease progression, which can be graded as low (LOWIN), latent (LATIN), or patent (PATIN) risk. Other possible combinations of these 5 SNPs were grouped in the average group (AVERIN). In the current study, differences in estimated breeding values (EBVs) for several traits were analyzed using linear models in a large cohort of Holstein cows (N = 15656) genotyped across Spain in 2016 or 2017. After the assignment of each genotyped cow to a risk group, cows within the PATIN risk group (N = 1448) had a superior combined genetic index (2797.57), type genetic index (524.62), milk yield (653.92 kg), protein yield (21.77 kg), fat yield (24.82 kg) and economic merit index (125 Euros) compared with the other three risk groups. Statistically significant differences in the longevity scores between the cows that were included in the PATIN risk group (108.85) and the LOWIN (107.82) and AVERIN (107.92) groups were also observed. The associations between the genetic risk groups and PTB diagnostic results were validated in a population of 99 cows from a Spanish farm with a high prevalence of PTB. Significant differences in ELISA readings between the PATIN (65.49 %) and the AVERIN (15.97 %), LATIN (2.11 %), and LOWIN (3.27 %) groups were observed. In addition, significant differences in Map DNA copies/gram of feces were observed between the PATIN and the other three risk groups. These results together with the substantial economic impact of PTB in dairy cattle support the selection of the animals with less susceptibility to PTB in the Spanish breeding program.

Local Lung Immune Response to Mycobacterium bovis Challenge after BCG and M. bovis Heat-Inactivated Vaccination in European Badger (Meles meles).

Tuberculosis (TB) vaccination could be used as a key part of integrated strategies for the disease's control if an effective and safe vaccine under field conditions is obtained. Recent studies in Spain have evaluated the protective efficacy of two oral vaccines against experimental challenge with live intra-bronchial Mycobacterium bovis in captive badgers: the live-attenuated M. bovis BCG vaccine (Danish strain) and a heat-inactivated M. bovis (HIMB) vaccine. With the objective of increasing the knowledge of the cellular development progress of infection and generating further tools to discriminate between mild and severe TB lesions between and within animals, the immunopathology of tuberculous lesions was studied to characterize the local immune response (cell type profile) within lung granulomas from control (non-vaccinated), BCG vaccinated and HIMB-vaccinated experimentally infected badgers with M. bovis. Four immunohistochemical protocols, for the specific detection of macrophages, T lymphocytes, B lymphocytes and plasma cells within TB granulomas in formalin fixed sections of the right middle lung lobe (lobe targeted for the M. bovis delivery), were performed. Immunolabelled sections were scanned and five randomly selected areas were analyzed with digital image analysis software. The results were expressed as the proportion of the positively immunolabelled area within the total area of the selected site. Data was analyzed using the statistical analysis software (SAS). In the three treatment groups, macrophages were the most abundant inflammatory cells within the granulomas, followed by B lymphocytes and plasma cells. T lymphocyes were absent in those granulomas. This would suggest a predominance of a non-specific innate response mediated by phagocytic cells over an adaptative humoral immune response. The proportion of macrophages and plasma cells was higher in BCG and HIMB-vaccinated badgers, respectively, suggesting the establishment of an adaptative humoral response in HIMB-vaccinated badgers. The lower bacterial load at the lung level, as well as the volume of lesions in lungs using magnetic resonance imaging in badgers with the HIMB vaccine in relation with local immune response presented, must be highlighted, since it would be an advantage in favor of its use under field conditions in terms of reducing TB transmission and environmental contamination.

RNA-Seq analysis of ileocecal valve and peripheral blood from Holstein cattle infected with Mycobacterium avium subsp. paratuberculosis revealed dysregulation of the CXCL8/IL8 signaling pathway.

Paratuberculosis is chronic granulomatous enteritis of ruminants caused by Mycobacterium avium subsp. paratuberculosis (MAP). Whole RNA-sequencing (RNA-Seq) is a promising source of novel biomarkers for early MAP infection and disease progression in cattle. Since the blood transcriptome is widely used as a source of biomarkers, we analyzed whether it recapitulates, at least in part, the transcriptome of the ileocecal valve (ICV), the primary site of MAP colonization. Total RNA was prepared from peripheral blood (PB) and ICV samples, and RNA-Seq was used to compare gene expression between animals with focal or diffuse histopathological lesions in gut tissues versus control animals with no detectable signs of infection. Our results demonstrated both shared, and PB and ICV-specific gene expression in response to a natural MAP infection. As expected, the number of differentially expressed (DE) genes was larger in the ICV than in the PB samples. Among the DE genes in the PB and ICV samples, there were some common genes irrespective of the type of lesion including the C-X-C motif chemokine ligand 8 (CXCL8/IL8), apolipoprotein L (APOLD1), and the interferon inducible protein 27 (IFI27). The biological processes (BP) enriched in the PB gene expression profiles from the cows with diffuse lesions included the killing of cells of other organism, defense response, immune response and the regulation of neutrophil chemotaxis. Two of these BP, the defense and immune response, were also enriched in the ICV from the cows with diffuse lesions. Metabolic analysis of the DE genes revealed that the N-glycan biosynthesis, bile secretion, one-carbon pool by folate and purine metabolism were significantly enriched in the ICV from the cows with focal lesions. In the ICV from cows with diffuse lesions; the valine, leucine and isoleucine degradation route, purine metabolism, vitamin digestion and absorption and the cholesterol routes were enriched. Some of the identified DE genes, BP and metabolic pathways will be studied further to develop novel diagnostic tools, vaccines and immunotherapeutics.