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

Lining the small intestine with mycobacteriophages protects from Mycobacterium avium subsp. paratuberculosis and eliminates fecal shedding

Johne's disease (JD), a chronic, infectious enteritis caused by Mycobacterium avium subsp. paratuberculosis (MAP), affects wild and domestic ruminants. There is no cure or effective prevention, and current vaccines have substantial limitations, leaving this disease widespread in all substantial dairy industries causing economic, and animal welfare implications. Mycobacteriophages (MPs) have been gaining interest in recent years and are proposed as a promising solution to curtailing MAP infection. Using a well-validated infection model, we have demonstrated the preventative potential of MPs to protect dairy calves against MAP infection. Calves were supplemented daily with a phage cocktail from birth till weaning at 2 m of age and inoculated with MAP at 2 wk of age. Infection status was measured for 4.5 mo through blood, fecal, and postmortem tissue samples. Our findings highlight the remarkable efficacy of orally administered MPs. Notably, fecal shedding of MAP was entirely eliminated within 10 wk, in contrast to the infected control group where shedding continued for the entirety of the trial period. Postmortem tissue culture analysis further supported the effectiveness of MPs, with only 1 out of 6 animals in the phage-treated group testing positive for MAP colonized tissues compared to 6 out of 6 animals in the infected control group. Additionally, plaque assay results demonstrated the ability of phages to persist within the intestinal tract. Collectively, these results underscore the potential of orally administered MP cocktails as a highly effective intervention strategy to combat JD in dairy calves and by extension in the dairy industry.

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.

The role of interleukin-10 receptor alpha (IL10Rα) in Mycobacterium avium subsp. paratuberculosis infection of a mammary epithelial cell line

BACKGROUND

Johne's disease is a chronic wasting disease caused by the bacterium Mycobacterium avium subspecies paratuberculosis (MAP). Johne's disease is highly contagious and MAP infection in dairy cattle can eventually lead to death. With no available treatment for Johne's disease, genetic selection and improvements in management practices could help reduce its prevalence. In a previous study, the gene coding interleukin-10 receptor subunit alpha (IL10Rα) was associated with Johne's disease in dairy cattle. Our objective was to determine how IL10Rα affects the pathogenesis of MAP by examining the effect of a live MAP challenge on a mammary epithelial cell line (MAC-T) that had IL10Rα knocked out using CRISPR/cas9. The wild type and the IL10Rα knockout MAC-T cell lines were exposed to live MAP bacteria for 72 h. Thereafter, mRNA was extracted from infected and uninfected cells. Differentially expressed genes were compared between the wild type and the IL10Rα knockout cell lines. Gene ontology was performed based on the differentially expressed genes to determine which biological pathways were involved.

RESULTS

Immune system processes pathways were targeted to determine the effect of IL10Rα on the response to MAP infection. There was a difference in immune response between the wild type and IL10Rα knockout MAC-T cell lines, and less difference in immune response between infected and not infected IL10Rα knockout MAC-T cells, indicating IL10Rα plays an important role in the progression of MAP infection. Additionally, these comparisons allowed us to identify other genes involved in inflammation-mediated chemokine and cytokine signalling, interleukin signalling and toll-like receptor pathways.

CONCLUSIONS

Identifying differentially expressed genes in wild type and ILR10α knockout MAC-T cells infected with live MAP bacteria provided further evidence that IL10Rα contributes to mounting an immune response to MAP infection and allowed us to identify additional potential candidate genes involved in this process. We found there was a complex immune response during MAP infection that is controlled by many genes.

No Evidence of Neutrophil Response Modulation in Goats after Immunization against Paratuberculosis with a Heat-Inactivated Vaccine

Neutrophils are believed to play a role in the initial stages of paratuberculosis, and it has recently been demonstrated that vaccination can modulate their function via priming or through epigenetic and metabolic reprogramming (training). Modulation of the neutrophil response against Mycobacterium avium subspecies paratuberculosis (Map) through vaccination has been demonstrated in a rabbit model but not in ruminants. Therefore, in the present work, the effect of vaccination on the response of caprine neutrophils against Map was studied. Neutrophils were isolated from non-vaccinated (n = 7) and Gudair®-vaccinated goat kids (n = 7), before vaccination and 30 days post-vaccination. Then, several neutrophil functions were quantified ex vivo: cell-free and anchored neutrophil extracellular trap (NET) release, phagocytosis, and the differential expression of several cytokines and TLR2. The induction of cell-free NETosis and TLR2 expression by Map is reported for the first time. However, vaccination showed no significant effect on any of the functions studied. This suggests that the protection conferred by Gudair® vaccination is based on mechanisms that are independent of the neutrophil function modulation. Further research into the impact of alternative vaccination strategies or the paratuberculosis infection stage on ruminant neutrophil function could provide valuable insights into its role in paratuberculosis.

Screening of mRNA markers in early bovine tuberculosis blood samples

Bovine tuberculosis (bTB) is a chronic zoonotic disease caused by Mycobacterium bovis. A large number of cattle are infected with bTB every year, resulting in huge economic losses. How to control bTB is an important issue in the current global livestock economy. In this study, the original transcriptome sequences related to this study were obtained from the dataset GSE192537 by searching the Gene Expression Omnibus (GEO) database. Our differential gene analysis showed that there were obvious biological activities related to immune activation and immune regulation in the early stage of bTB. Immune-related biological processes were more active in the early stage of bTB than in the late. There were obvious immune activation and immune cell recruitment in the early stage of bTB. Regulations in immune receptors are associated with pathophysiological processes of the early stage of bTB. A gene module consisting of 236 genes significantly related to the early stage of bTB was obtained by weighted gene co-expression network analysis, and 18 hub genes were further identified as potential biomarkers or therapeutic targets. Finally, by random forest algorithm and logistic regression modeling, FCRL1 was identified as a representative mRNA marker in early bTB blood. FCRL1 has the potential to be a diagnostic biomarker in early bTB.

Live attenuated Mycobacterium bovis strains combined with the encapsulated H65 antigen as a vaccine strategy against bovine tuberculosis in a mouse model

Mycobacterium bovis is an etiological agent of bovine tuberculosis (bTB) that also infects other mammals, including humans. The lack of an effective vaccine for the control of bTB highlights the need for developing new vaccines. In this study, we developed and evaluated an M. bovis strain deleted in the virulence genes phoP, esxA and esxB as a vaccine candidate against bTB in BALBc mice. The evaluated strains were the new live vaccine and BCG, alone or in combination with ncH65vD. The immunogen ncH65vD is a fusion protein H65, encapsulated together with vitamin D3, within the oily body of a nanocapsule composed of an antigen-loading polymeric shell. All vaccines conferred protection against the M. bovis challenge. However, no significant differences were detected among the vaccinated groups regarding bacterial loads in lungs and spleen. Mice vaccinated with the mutant strain plus ncH65vD showed negative Ziehl Neelsen staining of mycobacteria in their lungs, which suggests better control of bacteria replication according to this protection parameter. Consistently, this vaccination scheme showed the highest proportion of CD4 + T cells expressing the protection markers PD-1 and CXCR3 among the vaccinated groups. Correlation studies showed that PD-1 and CXCR3 expression levels in lung-resident CD4 T cells negatively correlated with the number of colony forming units of M. bovis in the lungs of mice. Therefore, the results suggest a link between the presence of PD-1 + and CXCR3 + cells at the site of the immune response against mycobacteria and the level of mycobacterial loads.

Alternative splicing of pre-mRNA modulates the immune response in Holstein cattle naturally infected with Mycobacterium avium subsp. paratuberculosis

Little is known about the role of alternative splicing (AS) in regulating gene expression in Mycobacteria-infected individuals in distinct stages of infection. Pre-mRNA AS consists of the removal of introns and the assembly of exons contained in eukaryotic genes. AS events can influence transcript stability or structure with important physiological consequences. Using RNA-Seq data from peripheral blood (PB) and ileocecal valve (ICV) samples collected from Holstein cattle with focal and diffuse paratuberculosis (PTB)-associated histopathological lesions in gut tissues and without lesions (controls), we detected differential AS profiles between the infected and control groups. Four of the identified AS events were experimentally validated by reverse transcription-digital droplet PCR (RT-ddPCR). AS events in several genes correlated with changes in gene expression. In the ICV of animals with diffuse lesions, for instance, alternatively spliced genes correlated with changes in the expression of genes involved in endocytosis, antigen processing and presentation, complement activation, and several inflammatory and autoimmune diseases in humans. Taken together, our results identified common mechanisms of AS involvement in the pathogenesis of PTB and human diseases and shed light on novel diagnostic and therapeutic interventions to control these diseases.

Heritability estimates and genome-wide association study of methane emission traits in Nellore cattle

The objectives of the present study were to estimate the heritability for daily methane emission (CH4) and residual daily methane emission (CH4res) in Nellore cattle, as well as to perform genome-wide association studies (GWAS) to identify genomic regions and candidate genes influencing the genetic variation of CH4 and CH4res. Methane emission phenotypes of 743 Nellore animals belonging to 3 breeding programs were evaluated. CH4 was measured using the sulfur hexafluoride (SF6) tracer technique (which involves an SF6 permeation tube introduced into the rumen, and an appropriate apparatus on each animal), and CH4res was obtained as the difference between observed CH4 and CH4 adjusted for dry matter intake. A total of 6,252 genotyped individuals were used for genomic analyses. Data were analyzed with a univariate animal model by the single-step GBLUP method using the average information restricted maximum likelihood (AIREML) algorithm. The effects of single nucleotide polymorphisms (SNPs) were obtained using a single-step GWAS approach. Candidate genes were identified based on genomic windows associated with quantitative trait loci (QTLs) related to the 2 traits. Annotation of QTLs and identification of candidate genes were based on the initial and final coordinates of each genomic window considering the bovine genome ARS-UCD1.2 assembly. Heritability estimates were of moderate to high magnitude, being 0.42 ± 0.09 for CH4 and 0.21 ± 0.09 for CH4res, indicating that these traits will respond rapidly to genetic selection. GWAS revealed 11 and 15 SNPs that were significantly associated (P < 10-6) with genetic variation of CH4 and CH4res, respectively. QTLs associated with feed efficiency, residual feed intake, body weight, and height overlapped with significant markers for the traits evaluated. Ten candidate genes were present in the regions of significant SNPs; 3 were associated with CH4 and 7 with CH4res. The identified genes are related to different functions such as modulation of the rumen microbiota, fatty acid production, and lipid metabolism. CH4 and CH4res presented sufficient genetic variation and may respond rapidly to selection. Therefore, these traits can be included in animal breeding programs aimed at reducing enteric methane emissions across generations.

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.

Modeling Paratuberculosis in Laboratory Animals, Cells, or Tissues: A Focus on Their Applications for Pathogenesis, Diagnosis, Vaccines, and Therapy Studies

Paratuberculosis is a chronic granulomatous enteritis caused by Mycobacterium avium subsp. Paratuberculosis that affects a wide variety of domestic and wild animals. It is considered one of the diseases with the highest economic impact on the ruminant industry. Despite many efforts and intensive research, paratuberculosis control still remains controversial, and the existing diagnostic and immunoprophylactic tools have great limitations. Thus, models play a crucial role in understanding the pathogenesis of infection and disease, and in testing novel vaccine candidates. Ruminant animal models can be restricted by several reasons, related to space requirements, the cost of the animals, and the maintenance of the facilities. Therefore, we review the potential and limitations of the different experimental approaches currently used in paratuberculosis research, focusing on laboratory animals and cell-based models. The aim of this review is to offer a vision of the models that have been used, and what has been achieved or discovered with each one, so that the reader can choose the best model to answer their scientific questions and prove their hypotheses. Also, we bring forward new approaches that we consider worth exploring in the near future.

Lateral-flow assays for bovine paratuberculosis diagnosis

Mycobacterium avium subsp. paratuberculosis (MAP) causes bovine paratuberculosis (PTB). PTB is responsible for significant economic losses in dairy herds around the word. PTB control programs that rely on testing and culling of test-positive cows have been developed. Current diagnostics, such as ELISA for detecting MAP antibodies in serum samples and PCR detecting MAP DNA in feces, have inadequate sensitivity for detecting subclinical animals. Innovative "omics" technologies such as next-generation sequencing (NGS) technology-based RNA-sequencing (RNA-Seq), proteomics and metabolomics can be used to find host biomarkers. The discovered biomarkers (RNA, microRNAs, proteins, metabolites) can then be used to develop new and more sensitive approaches for PTB diagnosis. Traditional approaches for measuring host antibodies and biomarkers, such as ELISAs, northern blotting, quantitative reverse-transcriptase polymerase chain reaction (RT-qPCR), cDNA microarrays, and mass spectrometry are time-consuming, expensive, and sometimes exhibit poor sensitivity. With the rapid development of nanotechnology, low-cost monitoring devices for measuring antibodies against MAP proteins in point-of-care (POC) settings have been developed. Lateral flow assays (LFAs), in particular, are thought to be appropriate for the on-site detection of antibodies to MAP antigens and/or host biomarkers. This review aims to summarize LFAs that have recently been developed to accurately detect antibodies against MAP antigens, as well as the benefits that host biomarkers linked with MAP infection give to PTB diagnosis. The identification of these novel biomarkers could be the basis for the development of new LFAs. The dairy industry and producers are likely to benefit from reliable and rapid technologies capable of detecting MAP infection in situ to establish a quick and sensitive PTB diagnosis.

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.

Dysbiosis of intestinal homeostasis contribute to Whitmania pigra edema disease

Whitmania pigra is widely used in traditional Chinese medicine. However, W. pigra is being threatened by an edema disease with unknown causes (WPE). In this study, a comprehensive exploration of virome, microbiome, and metabolome aberrations in the intestine of W. pigra was performed to address the aetiology of WPE. Virome analysis indicated that eukaryotic viruses did not contribute to WPE, whereas an expansion of Caudovirales was observed in WPE. Compared to the control, the microbial richness and diversity in diseased W. pigra decreased remarkably. Nine genera, including Aeromonas, Anaerotruncus, Vibrio, Proteocatella, Acinetobacter, and Brachyspira were overrepresented in WPE, whereas eleven genera, including Bifidobacterium, Phascolarctobacterium, Lactobacillus, Bacillus and AF12, were enriched in healthy individuals. Furthermore, certain metabolites, especially amino acids, short-chain fatty acids, and bile acids, were found to be linked to intestinal microbiota alterations in WPE. An integration of the microbiome and metabolome in WPE found that dysbiosis of the gut microbiota or metabolites caused WPE. Notably, W. pigra accepted intestinal microbiota transplantation from WPE donors developed WPE clinical signs eventually, and the dysbiotic intestinal microbiota can be recharacterized in this recipient W. pigra. Strikingly, pathological features of metanephridium and uraemic toxin enrichment in the gut indicated a putative interconnection between the gut and metanephridium in WPE, which represents the prototype of the gut-kidney axis in mammals. These finding exemplify the conservation of "microecological Koch's postulates" from annelids to insects and other vertebrates, which provides a direction of prevention and treatment for WPE and opens a new insight into the pathogenesis of aquatic animal diseases from an ecological perspective.

Paneth cells in farm animals: current status and future direction

A healthy intestine plays an important role in the growth and development of farm animals. In small intestine, Paneth cells are well known for their regulation of intestinal microbiota and intestinal stem cells (ISCs). Although there has been a lot of studies and reviews on human and murine Paneth cells under intestinal homeostasis or disorders, little is known about Paneth cells in farm animals. Most farm animals possess Paneth cells in their small intestine, as identified by various staining methods, and Paneth cells of various livestock species exhibit noticeable differences in cell shape, granule number, and intestinal distribution. Paneth cells in farm animals and their antimicrobial peptides (AMPs) are susceptible to multiple factors such as dietary nutrients and intestinal infection. Thus, the comprehensive understanding of Paneth cells in different livestock species will contribute to the improvement of intestinal health. This review first summarizes the current status of Paneth cells in pig, cattle, sheep, horse, chicken and rabbit, and points out future directions for the investigation of Paneth cells in the reviewed animals.

Evaluation of the innate immune response of caprine neutrophils against Mycobacterium avium subspecies paratuberculosis in vitro

Neutrophils constitute an essential component of the innate immune response, readily killing most bacteria through phagocytosis, degranulation, and the release of neutrophil extracellular traps (NETs) among other mechanisms. These cells play an unclear role in mycobacterial infections such as Mycobacterium avium subspecies paratuberculosis (Map), the etiological agent of paratuberculosis, and its response is particularly understudied in ruminants. Herein, a wide set of techniques were adapted, or newly developed, to study the in vitro response of caprine neutrophils after Map infection. Immunofluorescence was used to demonstrate, simultaneously, chemotaxis, phagocytosis, degranulation, and NETs. The quantification of neutrophil phagocytic activity against Map at a 1:10 multiplicity of infection (MOI), through flow cytometry, showed values that varied from 4.54 to 5.63% of phagocyting neutrophils. By immunofluorescence, a 73.3 ± 14.5% of the fields showed NETs, and the mean release of DNA, attributable to NETosis, calculated through a fluorometric method, was 16.2 ± 3.5%. In addition, the RNA expression of TGF-β, TNF and IL-1β cytokines, measured through reverse transcription qPCR, was significantly higher in the two latter. Overall, neutrophil response was proportional to the number of bacteria. This work confirms that the simultaneous study of several neutrophil mechanisms, and the combination of different methodologies, are essential to reach a comprehensive understanding of neutrophil response against pathogens, demonstrates that, in vitro, caprine neutrophils display a strong innate response against Map, using their entire repertoire of effector functions, and sets the basis for further in vitro and in vivo studies on the role of neutrophils in paratuberculosis.

Assessment of tuberculosis biomarkers in paratuberculosis-infected cattle

Introduction

Mycobacterium bovis and Mycobacterium avium subsp. paratuberculosis, respectively the causative agents of bovine tuberculosis (bTB) and bovine paratuberculosis (PTB), share a high number of antigenic proteins. This characteristics makes the differential diagnosis of the diseases difficult. The interferon gamma (IFN-γ), C-X-C motif chemokine ligand 10 (CXCL10), matrix metallopeptidase 9 (MMP9), interleukin 22 (IL-22) and thrombospondin 1 (THBS1) bovine genes have already been shown to be accurate transcriptional biomarkers of bTB. In order to improve the diagnosis of bTB and PTB, in the present study we evaluated the risk of false positivity of these bTB biomarkers in cattle with PTB.

Material and Methods

The transcription of these genes was studied in 13 PTB-infected cattle, using Mycobacterium avium subsp. paratuberculosis (MAP)-stimulated peripheral blood mononuclear cells (PBMC).

Results

Overall, the levels of IFN-γ, CXCL10, MMP9 and IL-22 transcripts in MAP-stimulated PBMC failed to differentiate animals with PTB from healthy animals. However, as bTB-afflicted cattle do, the MAP-infected group also displayed a lower level of THBS1 transcription than the non-infected animals.

Conclusion

The results of this study add new specificity attributes to the levels of transcription of IFN-γ, CXCL10, MMP9 and IL-22 as biomarkers for bTB.

Control of Mycobacterium avium subsp. paratuberculosis load within infected bovine monocyte-derived macrophages is associated with host genetics

The genetic loci influencing individual resistance to Mycobacterium avium subsp. paratuberculosis (MAP) infection are still largely unknown. In the current study, we searched for genetic loci associated with resistance to MAP infection by evaluating the performance of monocyte-derived macrophages (MDMs) isolated from the peripheral blood of 75 healthy Holsteins cows and infected ex vivo with MAP. Bacterial load (log colony-forming units, log CFUs) within MDMs was quantified at 2 h and 7 days p. i. using a BACTEC MGIT 960 instrument. In addition, the expression levels of some genes with important roles in the innate immune response including epiregulin (EREG), complement component C3 (C3), galectin-9 (Gal9), and nitric oxide (NO^-) were measured in the supernatant of the infected cells. DNA from peripheral blood samples of the animals included in the study was isolated and genotyped with the EuroG MD bead Chip (44,779 single nucleotide-polymorphisms, SNPs). Linear mixed models were used to calculate the heritability (h² ) estimates for each indicator of MDM performance, MAP load within MDMs and EREG, C3, Gal9, and NO^-expression. After performing a genome-wide association study, the only phenotypes that showed SNPs with a significant association were the bacterial load within MDMs at 2 h (h² = 0. 87) and 7 days (h² = 0.83) p.i. A total of 6 SNPs, 5 candidate genes, and one microRNA on the Bos taurus chromosomes BTA2, BTA17, BTA18, and BTA21 were associated with MAP load at 2 h p.i. Overlap was seen in two SNPs associated with the log CFUs at 2 h and 7 d p.i. The identified SNPs had negative regression coefficients, and were, therefore, associated with a low bacterial load within MDMs. Some of the identified SNPs were located within QTLs previously associated with longevity, reproductive, and udder health traits. Some of the identified candidate genes; Oxysterol Binding Protein Like 6, Cysteine and Serine Rich Nuclear Protein 3, and the Coiled-Coil Domain Containing 92 regulate cellular cholesterol trafficking and efflux, apoptosis, and interferon production, respectively. Taken together, our results define a heritable and distinct immunogenetic profile in MAP-infected macrophages designed to limit bacterial load early after infection.

The Upregulation of Cathepsin G Is Associated with Resistance to Bovine Paratuberculosis

An in silico genomic-transcriptomic combined approach allowed the identification of a polymorphism (cis-eQTL-rs41976219) in the Bos taurus genome associated with the CTSG mRNA expression in bovine blood samples, which suggests that individual genetic variation might modulate the CTSG transcriptional response. In the current study, a sandwich ELISA is used to measure the CTSG protein levels in supernatants of monocyte-derived macrophages (MDMs) isolated from cows with the AA (n = 5) and AC (n = 11) genotypes for the rs41976219 and infected ex vivo with MAP. Cows with the AC genotype have significantly higher CTSG protein levels (1.85 ng/mL) in the supernatants of enriched CD14+-MDMs after 2 h of infection when compared with infected CD14+-MDMs from cows with the AA genotype (1.68 ng/mL). Statistically significant differences in the intracellular MAP load at 7 d p.i. are observed between animals with the AA (2.16 log CFUs) and AC (1.44 log CFUs) genotypes. Finally, the association between the rs41976219 allelic variants and resistance to PTB is tested in a larger cattle population (n = 943) classified according to the presence (n = 442) or absence (n = 501) of PTB-associated lesions. The presence of the two minor alleles in the rs41976219 (CC) is more frequent among healthy cows than in cows with PTB-associated lesions in gut tissues (2.2% vs. 1.4%, OR = 0.61). In agreement with this, the CTSG levels in plasma samples of cows without lesions in gut tissues and with the CC (n = 8) genotype are significantly higher than in the plasmas of cows with the AA + AC (n = 36) genotypes.

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.

Identification of Alternative Splicing Events Associated with Paratuberculosis in Dairy Cattle Using Multi-Tissue RNA Sequencing Data

Paratuberculosis is a major endemic disease caused by Mycobacterium avium subspecies paratuberculosis (MAP) infection and leads to huge economic loss in the dairy sector worldwide. Alternative splicing (AS) events, playing indispensable regulatory roles in many protein functions and biological pathways, are shown to be associated with complex traits and diseases. In this study, by integrating the RNA sequencing (RNA-seq) data of 24 samples from three tissues (peripheral blood, jejunum and salivary gland) of Holstein cows, we obtained 2,706,541,696 uniquely mapped reads in total that represented 12,870 expressed genes, and detected 4285 differentially expressed genes (DEGs) between MAP-infected and healthy cows (p < 0.05). Of them, 92 differentially expressed splicing factors (DESFs) were included. Further, 119, 150 and 68 differential alternative splicing (DAS) events between MAP-infected and healthy cows were identified in peripheral blood, jejunum and salivary glands, respectively. Of note, six DAS events were highly and significantly correlated with the DESFs (R2 > 0.9; p < 0.01), and their corresponding genes (COPI coat complex subunit gamma 2gene (COPG2), kinesin family member 2C gene (KIF2C), exocyst complex component 7 (EXOC7), Rab9 effector protein with kelch motifs gene (RABEPK), deoxyribonuclease 1 gene (DNASE1) and early endosome antigen 1gene (EEA1)) were significantly enriched in immune response such as vesicle-mediated transport, regulation of acute inflammatory response and tuberculosis through gene ontology (GO) and KEGG analysis. KS test showed that the DAS events in the EXOC7 and KIF2C genes indeed displayed differences between MAP-infected cows and healthy cows. The DAS in EXOC7 might produce a new protein sequence with lack of 23 amino acids, and the DAS in KIF2C induced a stop codon of premature occurrence and resulted in a lack of functional domain. In summary, this study identified the DAS events and corresponding genes related to MAP-infection base on the RNA-seq data from multiple tissues of Holstein cows, providing novel insights into the regulatory mechanisms underpinning paratuberculosis in dairy cattle.

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.

Regionally Distinct Immune and Metabolic Transcriptional Responses in the Bovine Small Intestine and Draining Lymph Nodes During a Subclinical Mycobacterium avium subsp. paratuberculosis Infection

Mycobacterium avium subsp. paratuberculosis (MAP) is the causative infectious agent of Johne’s disease (JD), an incurable granulomatous enteritis affecting domestic livestock and other ruminants around the world. Chronic MAP infections usually begin in calves with MAP uptake by Peyer’s patches (PP) located in the jejunum (JE) and ileum (IL). Determining host responses at these intestinal sites can provide a more complete understanding of how MAP manipulates the local microenvironment to support its long-term survival. We selected naturally infected (MAPinf, n=4) and naive (MAPneg, n=3) cows and transcriptionally profiled the JE and IL regions of the small intestine and draining mesenteric lymph nodes (LN). Differentially expressed (DE) genes associated with MAP infection were identified in the IL (585), JE (218), jejunum lymph node (JELN) (205), and ileum lymph node (ILLN) (117). Three DE genes (CD14, LOC616364 and ENSBTAG00000027033) were common to all MAPinf versus MAPneg tissues. Functional enrichment analysis revealed immune/disease related biological processes gene ontology (GO) terms and pathways predominated in IL tissue, indicative of an activated immune response state. Enriched GO terms and pathways in JE revealed a distinct set of host responses from those detected in IL. Regional differences were also identified between the mesenteric LNs draining each intestinal site. More down-regulated genes (52%) and fewer immune/disease pathways (n=5) were found in the ILLN compared to a higher number of up-regulated DE genes (56%) and enriched immune/disease pathways (n=13) in the JELN. Immunohistochemical staining validated myeloid cell transcriptional changes with increased CD172-positive myeloid cells in IL and JE tissues and draining LNs of MAPinf versus MAPneg cows. Several genes, GO terms, and pathways related to metabolism were significantly DE in IL and JE, but to a lesser extent (comparatively fewer enriched metabolic GO terms and pathways) in JELN suggesting distinct regional metabolic changes in IL compared to JE and JELN in response to MAP infection. These unique tissue- and regional-specific differences provides novel insight into the dichotomy in host responses to MAP infection that occur throughout the small intestine and mesenteric LN of chronically MAP infected cows.

Whole Genome Methylation Analysis Reveals Role of DNA Methylation in Cow's Ileal and Ileal Lymph Node Responses to Mycobacterium avium subsp. paratuberculosis Infection

Johne's Disease (JD), caused by Mycobacterium avium subsp paratuberculosis (MAP), is an incurable disease of ruminants and other animal species and is characterized by an imbalance of gut immunity. The role of MAP infection on the epigenetic modeling of gut immunity during the progression of JD is still unknown. This study investigated the DNA methylation patterns in ileal (IL) and ileal lymph node (ILLN) tissues from cows diagnosed with persistent subclinical MAP infection over a one to 4 years period. DNA samples from IL and ILLN tissues from cows negative (MAPneg) (n = 3) or positive for MAP infection (MAPinf) (n = 4) were subjected to whole genome bisulfite sequencing. A total of 11,263 and 62,459 differentially methylated cytosines (DMCs), and 1259 and 8086 differentially methylated regions (DMRs) (FDR<0.1) were found between MAPinf and MAPneg IL and ILLN tissues, respectively. The DMRs were found on 394 genes (denoted DMR genes) in the IL and on 1305 genes in the ILLN. DMR genes with hypermethylated promoters/5'UTR [3 (IL) and 88 (ILLN)] or hypomethylated promoters/5'UTR [10 (IL) and 25 (ILLN)] and having multiple functions including response to stimulus/immune response (BLK, BTC, CCL21, AVPR1A, CHRNG, GABRA4, TDGF1), cellular processes (H2AC20, TEX101, GLA, NCKAP5L, RBM27, SLC18A1, H2AC20BARHL2, NLGN3, SUV39H1, GABRA4, PPA1, UBE2D2) and metabolic processes (GSTO2, H2AC20, SUV39H1, PPA1, UBE2D2) are potential DNA methylation candidate genes of MAP infection. The ILLN DMR genes were enriched for more biological process (BP) gene ontology (GO) terms (n = 374), most of which were related to cellular processes (27.6%), biological regulation (16.6%), metabolic processes (15.4%) and response to stimulus/immune response (8.2%) compared to 75 BP GO terms (related to cellular processes, metabolic processes and transport, and system development) enriched for IL DMR genes. ILLN DMR genes were enriched for more pathways (n = 47) including 13 disease pathways compared with 36 enriched pathways, including 7 disease/immune pathways for IL DMR genes. In conclusion, the results show tissue specific responses to MAP infection with more epigenetic changes (DMCs and DMRs) in the ILLN than in the IL tissue, suggesting that the ILLN and immune processes were more responsive to regulation by methylation of DNA relative to IL tissue. Our data is the first to demonstrate a potential role for DNA methylation in the pathogenesis of MAP infection in dairy cattle.

Multipopulational transcriptome analysis of post-weaned beef cattle at arrival further validates candidate biomarkers for predicting clinical bovine respiratory disease

Bovine respiratory disease (BRD) remains the leading infectious disease in post-weaned beef cattle. The objective of this investigation was to contrast the at-arrival blood transcriptomes from cattle derived from two distinct populations that developed BRD in the 28 days following arrival versus cattle that did not. Forty-eight blood samples from two populations were selected for mRNA sequencing based on even distribution of development (n = 24) or lack of (n = 24) clinical BRD within 28 days following arrival; cattle which developed BRD were further stratified into BRD severity cohorts based on frequency of antimicrobial treatment: treated once (treated_1) or treated twice or more and/or died (treated_2+). Sequenced reads (~ 50 M/sample, 150 bp paired-end) were aligned to the ARS-UCD1.2 bovine genome assembly. One hundred and thirty-two unique differentially expressed genes (DEGs) were identified between groups stratified by disease severity (healthy, n = 24; treated_1, n = 13; treated_2+, n = 11) with edgeR (FDR ≤ 0.05). Differentially expressed genes in treated_1 relative to both healthy and treated_2+ were predicted to increase neutrophil activation, cellular cornification/keratinization, and antimicrobial peptide production. Differentially expressed genes in treated_2+ relative to both healthy and treated_1 were predicted to increase alternative complement activation, decrease leukocyte activity, and increase nitric oxide production. Receiver operating characteristic (ROC) curves generated from expression data for six DEGs identified in our current and previous studies (MARCO, CFB, MCF2L, ALOX15, LOC100335828 (aka CD200R1), and SLC18A2) demonstrated good-to-excellent (AUC: 0.800–0.899; ≥ 0.900) predictability for classifying disease occurrence and severity. This investigation identifies candidate biomarkers and functional mechanisms in at arrival blood that predicted development and severity of BRD.

Transcriptional Profiling of Early and Late Phases of Bovine Tuberculosis

Bovine tuberculosis, caused by Mycobacterium tuberculosis var. bovis (M. bovis), is an important enzootic disease affecting mainly cattle, worldwide. Despite the implementation of national campaigns to eliminate the disease, bovine tuberculosis remains recalcitrant to eradication in several countries. Characterizing the host response to M. bovis infection is crucial for understanding the immunopathogenesis of the disease and for developing better control strategies. To profile the host responses to M. bovis infection, we analyzed the transcriptome of whole blood cells collected from experimentally infected calves with a virulent strain of M. bovis using RNA transcriptome sequencing (RNAseq). Comparative analysis of calf transcriptomes at early (8 weeks) vs. late (20 weeks) aerosol infection with M. bovis revealed divergent and unique profile for each stage of infection. Notably, at the early time point, transcriptional upregulation was observed among several of the top-ranking canonical pathways involved in T-cell chemotaxis. At the late time point, enrichment in the cell mediated cytotoxicity (e.g. Granzyme B) was the predominant host response. These results showed significant change in bovine transcriptional profiles and identified networks of chemokine receptors and monocyte chemoattractant protein (CCL) co-regulated genes that underline the host-mycobacterial interactions during progression of bovine tuberculosis in cattle. Further analysis of the transcriptomic profiles identified potential biomarker targets for early and late phases of tuberculosis in cattle. Overall, the identified profiles better characterized identified novel immunomodulatory mechanisms and provided a list of targets for further development of potential diagnostics for tuberculosis in cattle.

Identification of loci associated with pathological outcomes in Holstein cattle infected with Mycobacterium avium subsp. paratuberculosis using whole-genome sequence data

Bovine paratuberculosis (PTB), caused by Mycobacterium avium subsp. paratuberculosis (MAP), is a chronic granulomatous enteritis that affects cattle worldwide. According to their severity and extension, PTB-associated histological lesions have been classified into the following groups; focal, multifocal, and diffuse. It is unknown whether these lesions represent sequential stages or divergent outcomes. In the current study, the associations between host genetic and pathology were explored by genotyping 813 Spanish Holstein cows with no visible lesions (N = 373) and with focal (N = 371), multifocal (N = 33), and diffuse (N = 33) lesions in gut tissues and regional lymph nodes. DNA from peripheral blood samples of these animals was genotyped with the bovine EuroG MD Bead Chip, and the corresponding genotypes were imputed to whole-genome sequencing (WGS) data using the 1000 Bull genomes reference population. A genome-wide association study (GWAS) was performed using the WGS data and the presence or absence of each type of histological lesion in a case–control approach. A total of 192 and 92 single nucleotide polymorphisms (SNPs) defining 13 and 9 distinct quantitative trait loci (QTLs) were highly-associated (P ≤ 5 × 10⁻⁷) with the multifocal (heritability = 0.075) and the diffuse (heritability = 0.189) lesions, respectively. No overlap was seen in the SNPs controlling these distinct pathological outcomes. The identified QTLs overlapped with some QTLs previously associated with PTB susceptibility, bovine tuberculosis susceptibility, clinical mastitis, somatic cell score, bovine respiratory disease susceptibility, tick resistance, IgG level, and length of productive life. Pathway analysis with candidate genes overlapping the identified QTLs revealed a significant enrichment of the keratinization pathway and cholesterol metabolism in the animals with multifocal and diffuse lesions, respectively. To test whether the enrichment of SNP variants in candidate genes involved in the cholesterol metabolism was associated with the diffuse lesions; the levels of total cholesterol were measured in plasma samples of cattle with focal, multifocal, or diffuse lesions or with no visible lesions. Our results showed reduced levels of plasma cholesterol in cattle with diffuse lesions. Taken together, our findings suggested that the variation in MAP-associated pathological outcomes might be, in part, genetically determined and indicative of distinct host responses.

Johne's Disease in Dairy Cattle: An Immunogenetic Perspective

Johne's disease (JD), also known as paratuberculosis, is a severe production-limiting disease with significant economic and welfare implications for the global cattle industry. Caused by infection with Mycobacterium avium subspecies paratuberculosis (MAP), JD manifests as chronic enteritis in infected cattle. In addition to the economic losses and animal welfare issues associated with JD, MAP has attracted public health concerns with potential association with Crohn's disease, a human inflammatory bowel disease. The lack of effective treatment options, such as a vaccine, has hampered JD control resulting in its increasing global prevalence. The disease was first reported in 1895, but in recognition of its growing economic impact, extensive recent research facilitated by a revolution in technological approaches has led to significantly enhanced understanding of the immunological, genetic, and pathogen factors influencing disease pathogenesis. This knowledge has been derived from a variety of diverse models to elucidate host-pathogen interactions including in vivo and in vitro experimental infection models, studies measuring immune parameters in naturally-infected animals, and by studies conducted at the population level to enable the estimation of genetic parameters, and the identification of genetic markers and quantitative trait loci (QTL) putatively associated with susceptibility or resistance to JD. The main objectives of this review are to summarize these recent developments from an immunogenetics perspective and attempt to extract the principal and common findings emerging from this wealth of recent information. Based on these analyses, and in light of emerging technologies such as gene-editing, we conclude by discussing potential future avenues for effectively mitigating JD in cattle.

Oral vaccination stimulates neutrophil functionality and exerts protection in a Mycobacterium avium subsp. paratuberculosis infection model

Mycobacterium avium subsp. paratuberculosis (Map) causes paratuberculosis (PTB), a granulomatous enteritis in ruminants that exerts high economic impact on the dairy industry worldwide. Current vaccines have shown to be cost-effective against Map and in some cases confer beneficial non-specific effects against other pathogens suggesting the existence of trained immunity. Although Map infection is mainly transmitted by the fecal-oral route, oral vaccination has not been deeply studied. Therefore, the aim of this study was to compare the oral route with a set of mycobacterial and non-mycobacterial vaccines with a subcutaneously administered commercially available vaccine. Training effects on polymorphonuclear neutrophils (PMNs) and homologous and heterologous in vivo protection against Map were investigated in the rabbit infection model. Oral vaccination with inactivated or live vaccines was able to activate mucosal immunity as seen by elevation of serum IgA and the expression of IL4 in peripheral blood mononuclear cells (PBMCs). In addition, peripheral PMN phagocytosis against Map was enhanced by vaccination and extracellular trap release against Map and non-related pathogens was modified by both, vaccination and Map-challenge, indicating trained immunity. Finally, PBMCs from vaccinated animals stimulated in vitro with Map antigens showed a rapid innate activation cytokine profile. In conclusion, our data show that oral vaccination against PTB can stimulate neutrophil activity and both innate and adaptive immune responses that correlate with protection.

RNA-Seq Transcriptome Analysis of Peripheral Blood From Cattle Infected With Mycobacterium bovis Across an Experimental Time Course

Bovine tuberculosis, caused by infection with members of the Mycobacterium tuberculosis complex, particularly Mycobacterium bovis, is a major endemic disease affecting cattle populations worldwide, despite the implementation of stringent surveillance and control programs in many countries. The development of high-throughput functional genomics technologies, including RNA sequencing, has enabled detailed analysis of the host transcriptome to M. bovis infection, particularly at the macrophage and peripheral blood level. In the present study, we have analysed the transcriptome of bovine whole peripheral blood samples collected at −1 week pre-infection and +1, +2, +6, +10, and +12 weeks post-infection time points. Differentially expressed genes were catalogued and evaluated at each post-infection time point relative to the −1 week pre-infection time point and used for the identification of putative candidate host transcriptional biomarkers for M. bovis infection. Differentially expressed gene sets were also used for examination of cellular pathways associated with the host response to M. bovis infection, construction of de novo gene interaction networks enriched for host differentially expressed genes, and time-series analyses to identify functionally important groups of genes displaying similar patterns of expression across the infection time course. A notable outcome of these analyses was identification of a 19-gene transcriptional biosignature of infection consisting of genes increased in expression across the time course from +1 week to +12 weeks post-infection.

Bovine Intelectin 2 Expression as a Biomarker of Paratuberculosis Disease Progression

Simple Summary

The potential of the bovine intelectin 2 as a biomarker of Mycobacterium avium subsp. paratuberculosis infection was investigated using quantitative immunohistochemical analysis of ileocecal valve samples of animals with increasing degrees of lesion severity (focal, multifocal and diffuse histological lesions) and control animals without detected lesions. Significant differences were observed in the mean number of intelectin 2 immunolabelled cells between the three histopathological types and the control. Specifically, the mean number of intelectin 2 labelled cells was indicative of disease progression as the focal group had the highest number of intelectin 2 secreting cells followed by the multifocal, diffuse and control groups indicating that intelectin 2 is a good biomarker for the different stages of Mycobacterium avium subsp. paratuberculosis infection. Quantification of bovine intelectin 2 secreting cells could constitute a good post-mortem tool, complementary to histopathology, to improve detection of Mycobacterium avium subsp. Paratuberculosis infections, especially latent forms of infection.

Abstract

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.

Bovine Neutrophils Release Extracellular Traps and Cooperate With Macrophages in Mycobacterium avium subsp. paratuberculosis clearance In Vitro

Mycobacterium avium subsp. paratuberculosis (Map) is the underlying pathogen causing bovine paratuberculosis (PTB), an enteric granulomatous disease that mainly affects ruminants and for which an effective treatment is needed. Macrophages are the primary target cells for Map, which survives and replicates intracellularly by inhibiting phagosome maturation. Neutrophils are present at disease sites during the early stages of the infection, but seem to be absent in the late stage, in contrast to healthy tissue. Although neutrophil activity has been reported to be impaired following Map infection, their role in PTB pathogenesis has not been fully defined. Neutrophils are capable of releasing extracellular traps consisting of extruded DNA and proteins that immobilize and kill microorganisms, but this mechanism has not been evaluated against Map. Our main objective was to study the interaction of neutrophils with macrophages during an in vitro mycobacterial infection. For this purpose, neutrophils and macrophages from the same animal were cultured alone or together in the presence of Map or Mycobacterium bovis Bacillus-Calmette-Guérin (BCG). Extracellular trap release, mycobacteria killing as well as IL-1β and IL-8 release were assessed. Neutrophils released extracellular traps against mycobacteria when cultured alone and in the presence of macrophages without direct cell contact, but resulted inhibited in direct contact. Macrophages were extremely efficient at killing BCG, but ineffective at killing Map. In contrast, neutrophils showed similar killing rates for both mycobacteria. Co-cultures infected with Map showed the expected killing effect of combining both cell types, whereas co-cultures infected with BCG showed a potentiated killing effect beyond the expected one, indicating a potential synergistic cooperation. In both cases, IL-1β and IL-8 levels were lower in co-cultures, suggestive of a reduced inflammatory reaction. These data indicate that cooperation of both cell types can be beneficial in terms of decreasing the inflammatory reaction while the effective elimination of Map can be compromised. These results suggest that neutrophils are effective at Map killing and can exert protective mechanisms against Map that seem to fail during PTB disease after the arrival of macrophages at the infection site.

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.

Biomarkers as diagnostic tools for mycobacterial infections in cattle

Mycobacterial infections are widely distributed in animals and cause considerable economic losses, especially in livestock animals. Bovine paratuberculosis and bovine tuberculosis, which are representative mycobacterial infections in cattle, are difficult to diagnose using current-generation diagnostics due to their relatively long incubation periods. Thus, alternative diagnostic tools are needed for the detection of mycobacterial infections in cattle. A biomarker is an indicator present in biological fluids that reflects the biological state of an individual during the progression of a specific disease. Therefore, biomarkers are considered a potential diagnostic tool for various diseases. Recently, the number of studies investigating biomarkers as tools for diagnosing mycobacterial infections has increased. In human medicine, many diagnostic biomarkers have been developed and applied in clinical practice. In veterinary medicine, however, many such developments are still in the early stages. In this review, we summarize the current progress in biomarker research related to the development of diagnostic biomarkers for mycobacterial infections in cattle.

Epithelial processed Mycobacterium avium subsp. paratuberculosis induced prolonged Th17 response and suppression of phagocytic maturation in bovine peripheral blood mononuclear cells

Johne’s disease (JD) caused by Mycobacterium avium subsp. paratuberculosis (MAP) is a chronic, wasting infectious disease in ruminants that causes enormous economic losses to the dairy and beef cattle industries. Understanding the mechanism of persistency of MAP is key to produce novel ideas for the development of new diagnostic methods or prevention techniques. We sought interactions between the host and MAP using epithelial passage model, which mimic initial stage of infection. From the transcriptomic analysis of bovine immune cells (PBMCs), it was suggested that infection through the epithelial cells elicited prolonged Th17-derived immune response, as indicated by upregulation of IL-17A, IL-17F and RORC until 120 h p.i., compared to directly infected PBMCs. Global downregulation of gene expression was observed after 72 h p.i., especially for genes encoding cell surface receptors of phagocytic cells, such as Toll-like receptors and MHC class II molecules. In addition, the cholesterol efflux transporters ABCA1, ABCG1, and APOE, which are regulated by the LXR/RXR pathway, were downregulated. In summary, it would be suggested that the host initiate immune response to activate Th17-derived cytokines, and MAP survives persistently by altering the host adaptive immune response by suppressing surface receptors and manipulating lipid metabolism in phagocytic cells.

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.