Pathways and Genes Associated with Immune Dysfunction in Sheep Paratuberculosis

Retrospective Single Nucleotide Polymorphism Analysis of Host Resistance and Susceptibility to Ovine Johne's Disease Using Restored FFPE DNA

Johne's disease (JD), also known as paratuberculosis, is a chronic, untreatable gastroenteritis of ruminants caused by Mycobacterium avium subsp. paratuberculosis (MAP) infection. Evidence for host genetic resistance to disease progression exists, although it is limited due to the extended incubation period (years) and diagnostic challenges. To overcome this, previously restored formalin-fixed paraffin embedded tissue (FFPE) DNA from archived FFPE tissue cassettes was utilized for a novel retrospective case-control genome-wide association study (GWAS) on ovine JD. Samples from known MAP-infected flocks with ante- and postmortem diagnostic data were used. Cases (N = 9) had evidence of tissue infection, compared to controls (N = 25) without evidence of tissue infection despite positive antemortem diagnostics. A genome-wide efficient mixed model analysis (GEMMA) to conduct a GWAS using restored FFPE DNA SNP results from the Illumina Ovine SNP50 Bead Chip, identified 10 SNPs reaching genome-wide significance of p < 1 × 10-6 on chromosomes 1, 3, 4, 24, and 26. Pathway analysis using PANTHER and the Kyoto Encyclopedia of Genes and Genomes (KEGG) was completed on 45 genes found within 1 Mb of significant SNPs. Our work provides a framework for the novel use of archived FFPE tissues for animal genetic studies in complex diseases and further evidence for a genetic association in JD.

Cathepsin S Is More Abundant in Serum of Mycobacterium avium subsp. paratuberculosis-Infected Dairy Cows

Mycobacterium avium subsp. paratuberculosis (MAP) is the causative agent of bovine paratuberculosis, a chronic granulomatous enteritis leading to economic losses and posing a risk to human health due to its zoonotic potential. The pathogen cannot reliably be detected by standard methods, and immunological procedures during the infection are not well understood. Therefore, the aim of our study was to explore host-pathogen interactions in MAP-infected dairy cows and to improve diagnostic tests. Serum proteomics analysis using quantitative label-free LC-MS/MS revealed 60 differentially abundant proteins in MAP-infected dairy cows compared to healthy controls from the same infected herd and 90 differentially abundant proteins in comparison to another control group from an uninfected herd. Pathway enrichment analysis provided new insights into the immune response to MAP and susceptibility to the infection. Furthermore, we found a higher abundance of Cathepsin S (CTSS) in the serum of MAP-infected dairy cows, which is involved in multiple enriched pathways associated with the immune system. Confirmed with Western blotting, we identified CTSS as a potential biomarker for bovine paratuberculosis. This study enabled a better understanding of procedures in the host-pathogen response to MAP and improved detection of paratuberculosis-diseased cattle.

Advances in understanding the genetic architecture of antibody response to paratuberculosis in sheep by heritability estimate and LDLA mapping analyses and investigation of candidate regions using sequence-based data

BACKGROUND

Paratuberculosis is a contagious and incurable disease that is caused by Mycobacterium avium subsp. paratuberculosis (MAP) with significant negative effects on animal welfare and farm profitability. Based on a large naturally infected flock over 12 years, we analyzed repeated enzyme-linked immunosorbent assay tests (ELISA), OvineSNP50 BeadChip genotypes and whole-genome sequences imputed from 56 influential animals. The main goals were to estimate the genetic parameters of proxy traits for resistance to MAP, identify genomic regions associated with the host's immune response against MAP and search for candidate genes and causative mutations through association and functional annotation analyses of polymorphisms identified by sequencing.

RESULTS

Two variables were derived from ELISA tests. The first, a binary variable, assessed the infection status of each animal over the entire productive life, while the second considered the level of antibody recorded over time. Very similar results were obtained for both variables. Heritability estimates of about 0.20 were found and a significant region capturing 18% and 13% of the genetic variance was detected on ovine chromosome 20 by linkage disequilibrium and linkage analysis on OvineSNP50 positions. Functional annotation and association analyses on the imputed sequence polymorphisms that were identified in this region were carried out. No significant variants showed a functional effect on the genes that mapped to this region, most of which belong to the major histocompatibility complex class II (MHC II). However, the conditional analysis led to the identification of two significant polymorphisms that can explain the genetic variance associated with the investigated genomic region.

CONCLUSIONS

Our results confirm the involvement of the host's genetics in susceptibility to MAP in sheep and suggest that selective breeding may be an option to limit the infection. The estimated heritability is moderate with a relevant portion being due to a highly significant region on ovine chromosome 20. The results of the combined use of sequence-based data and functional analyses suggest several genes belonging to the MHC II as the most likely candidates, although no mutations in their coding regions showed a significant association. Nevertheless, information from genotypes of two highly significant polymorphisms in the region can enhance the efficiency of selective breeding programs.

Landscape genomics reveals adaptive divergence of indigenous sheep in different ecological environments of Xinjiang, China

Sheep are important livestock animals that have evolved under various ecological pressures. Xinjiang is a region with diverse and harsh environments that have shaped many local sheep breeds with unique characteristics and environmental adaptability. However, these breeds are losing ecological flexibility due to the promotion of intensive farming practices. Here we sequenced 14 local sheep breeds from Xinjiang and analyzed their genetic structure and gene flow with other sheep breeds from neighboring regions. The Tibetan Plateau was the geographic origin of Xinjiang native sheep evolution. We performed genome-environment association analysis and identified Bio9: Mean Temperature of Driest Quarter and Bio15: Precipitation Seasonality as the key environmental factors affecting Xinjiang local sheep and the key genes involved in their survival and adaptation. We classified Xinjiang native sheep breeds into six groups based on their differential genes by pairwise selective sweep analysis and Community Network Analysis. We analyzed transcriptome expression data of 832 sheep tissues and detected tissue-specific enrichment of six group-specific genes in different biological systems. Our results revealed the genetic basis of year-round estrus, drought tolerance, hypoxia resistance, and cold tolerance traits of Xinjiang sheep breeds. Moreover, we proposed conservation strategies for Xinjiang local sheep breeds and provided theoretical guidance for breeding new sheep breeds under global extreme environments.

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.

Transcriptomic meta-analysis reveals unannotated long non-coding RNAs related to the immune response in sheep

Long non-coding RNAs (lncRNAs) are involved in several biological processes, including the immune system response to pathogens and vaccines. The annotation and functional characterization of lncRNAs is more advanced in humans than in livestock species. Here, we take advantage of the increasing number of high-throughput functional experiments deposited in public databases in order to uniformly analyse, profile unannotated lncRNAs and integrate 422 ovine RNA-seq samples from the ovine immune system. We identified 12302 unannotated lncRNA genes with support from independent CAGE-seq and histone modification ChIP-seq assays. Unannotated lncRNAs showed low expression levels and sequence conservation across other mammal species. There were differences in expression levels depending on the genomic location-based lncRNA classification. Differential expression analyses between unstimulated and samples stimulated with pathogen infection or vaccination resulted in hundreds of lncRNAs with changed expression. Gene co-expression analyses revealed immune gene-enriched clusters associated with immune system activation and related to interferon signalling, antiviral response or endoplasmic reticulum stress. Besides, differential co-expression networks were constructed in order to find condition-specific relationships between coding genes and lncRNAs. Overall, using a diverse set of immune system samples and bioinformatic approaches we identify several ovine lncRNAs associated with the response to an external stimulus. These findings help in the improvement of the ovine lncRNA catalogue and provide sheep-specific evidence for the implication in the general immune response for several lncRNAs.

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.

The Paratuberculosis Paradigm Examined: A Review of Host Genetic Resistance and Innate Immune Fitness in Mycobacterium avium subsp. Paratuberculosis Infection

Paratuberculosis, or Johne's Disease (JD) is a debilitating chronic enteritis mainly affecting ruminants caused by Mycobacterium avium subsp. paratuberculosis (MAP). This organism causes worldwide economic losses to the livestock industry, and is of public health importance due to the potential zoonotic risk between MAP and Crohn's disease (CD) in humans. Without economical treatments, or a vaccine capable of preventing infection without causing cross-reactions with bovine tuberculosis, test-and-cull methods for disease control are imperative. Unfortunately, difficulties in diagnostics and long subclinical stage hinder adequate control and is further complicated by variation in MAP exposure outcome. Interestingly, the majority of infections result in asymptomatic presentation and never progress to clinical disease. One contributing factor is host genetics, where polymorphisms in innate immune genes have been found to influence resistance and susceptibility to disease. Candidate genes identified across studies overlap with those found in CD and tuberculosis including; Solute carrier family 11 member 1 gene (SLC11A1), Nucleotide-binding-oligomerization domain containing gene 2 (NOD2), Major histocompatibility complex type II (MHC-II), and Toll-like receptor (TLR) genes. This review will highlight evidence supporting the vital role of these genes in MAP infection outcome, associated challenges, and implications for the future of JD research.

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.

Prediction of Johne's disease state based on quantification of T cell markers and their interaction with macrophages in the bovine intestine

Cell-mediated immune responses to Mycobacterium avium subsp. paratuberculosis (MAP) are regulated by various types of T lymphocytes. The aim of this study was to quantitate T cell subsets in the mid-ileum of cows naturally infected with MAP to identify differences during different stages of infection, and to determine whether these subsets could be used as predictors of disease state. Immunofluorescent labeling of T cell subsets and macrophages was performed on frozen mid-ileal tissue sections archived from naturally infected dairy cows in either subclinical or clinical disease status, and noninfected control cows. Comprehensive IF staining for CD4, CD8α, TcR1-N24 (gamma delta), FoxP3, CXCR3 and CCR9 served to define T cell subsets and was correlated with macrophages present. Clinically affected cows demonstrated significantly higher numbers of CXCR3⁺ (Th1-type) and CCR9⁺ (total small intestinal lymphocytes) cells at the site of infection compared to the subclinical cows and noninfected controls. Further, predictive modeling indicated a significant interaction between CXCR3⁺ and AM3K⁺ (macrophages) cells, suggesting that progression to clinical disease state aligns with increased numbers of these cell types at the site of infection. The ability to predict disease state with this model was improved from previous modeling using immunofluorescent macrophage data. Predictive modelling indicated an interaction between CXCR3⁺ and AM3K⁺ cells, which could more sensitively detect subclinical cows compared to clinical cows. It may be possible to use this knowledge to improve and develop an assay to detect subclinically infected animals with more confidence during the early stages of the disease.

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.

Paratuberculosis in Captive Scimitar-Horned Oryxes (Oryx dammah)

Paratuberculosis, a chronic disease caused by Mycobacterium avium subsp. paratuberculosis (MAP), in ten scimitar-horned oryxes (SHOs) hosted in an Italian zoological park and originating from a Slovakian flock, was documented by pathology, molecular, cultural, and serological testing. The infection origin in this threatened species was also investigated by genomic analyses. Following the death of six of the 10 SHOs, serial investigations of dead and alive animals were performed. Necropsy, carried out on five out of six animals, identified intestinal thickening and mesenteric lymphadenomegaly in one of the animals. Histopathology (5/6) revealed lepromatous (2/5) and tuberculoid (2/5) intestinal forms or lack of lesions (1/5). Ziehl-Neelsen and immunohistochemistry stains identified two multibacillary, two paucibacillary forms, and one negative case. MAP was identified by quantitative PCR (qPCR) in tissue samples in five out of five SHOs and was microbiologically isolated from two of the three animals whose fresh tissue samples were available. Fecal samples were collected in four of the six dead animals: all four resulted positive to qPCR and in MAP was isolated in three. ELISA identified MAP-specific antibodies in three of the five dead animals whose serum was available. qPCR identified MAP in the freshly deposited feces of two out of the four alive animals. From the feces of these two animals, MAP was microbiologically isolated in one case. All isolates were classified as MAP type C and profiled as INMV2 and MVS27 by molecular analysis. Genomic analysis of a field isolate revealed clusterization with a European clade but was more similar to Italian than East European isolates. Our findings underline that paratuberculosis should always be considered in zoological parks in which endangered species are hosted. Infection can be subclinical, and multiple combined testing techniques may be necessary.

Biomarkers for Detecting Resilience against Mycobacterial Disease in Animals

Paratuberculosis and bovine tuberculosis are two mycobacterial diseases of ruminants which have a considerable impact on livestock health, welfare, and production. These are chronic "iceberg" diseases which take years to manifest and in which many subclinical cases remain undetected. Suggested biomarkers to detect infected or diseased animals are numerous and include cytokines, peptides, and expression of specific genes; however, these do not provide a strong correlation to disease. Despite these advances, disease detection still relies heavily on dated methods such as detection of pathogen shedding, skin tests, or serology. Here we review the evidence for suitable biomarkers and their mechanisms of action, with a focus on identifying animals that are resilient to disease. A better understanding of these factors will help establish new strategies to control the spread of these diseases.

A Novel Approach to Deliver a Mycobacterium avium subsp. paratuberculosis Antigen in Eukaryotic Cells

This study was aimed to express and deliver a Mycobacterium avium subsp. paratuberculosis antigen to macrophages using salmonella as carrier. The coding sequence of a fibronectin attachment protein which is expressed by Mycobacterium avium subsp. paratuberculosis was cloned into pcDNA3.1 (+) plasmid. The construct was introduced into the attenuated Salmonella typhimurium strain SL7207 (ΔhisG, ΔaroA) as carrier. In order to evaluate the delivery capacity of Salmonella and gene expression by antigen-presenting cells, the THP-1 derived macrophages were infected with the salmonella carrier. SDS-PAGE and western blot analysis showed the successful delivery and expression of targeted gene in THP-1 cell line. Although, in vitro stimulation of peripheral blood mononuclear cells with Salmonella containing plasmid did not trigger IFNγ production significantly. But it seems that this carrier can increase plasmid uptake and antigen expression by host intestinal antigen-presenting cells after mucosal administration. So, the construct can be used for further in vivo studies on the Salmonella carrier's efficiency in mycobacterial DNA vaccines.

Application of Transcriptomics to Enhance Early Diagnostics of Mycobacterial Infections, with an Emphasis on Mycobacterium avium ssp. paratuberculosis

Mycobacteria cause a wide variety of disease in human and animals. Species that infect ruminants include M. bovis and M. avium ssp. paratuberculosis (MAP). MAP is the causative agent of Johne’s disease in ruminants, which is a chronic granulomatous enteric infection that leads to severe economic losses worldwide. Characteristic of MAP infection is the long, latent phase in which intermittent shedding can take place, while diagnostic tests are unable to reliably detect an infection in this stage. This leads to unnoticed dissemination within herds and the presence of many undetected, silent carriers, which makes the eradication of Johne’s disease difficult. To improve the control of MAP infection, research is aimed at improving early diagnosis. Transcriptomic approaches can be applied to characterize host-pathogen interactions during infection, and to develop novel biomarkers using transcriptional profiles. Studies have focused on the identification of specific RNAs that are expressed in different infection stages, which will assist in the development and clinical implementation of early diagnostic tests.

Gene expression profiles during subclinical Mycobacterium avium subspecies paratuberculosis infection in sheep can predict disease outcome

Paratuberculosis in ruminants is caused by infection with Mycobacterium avium subspecies paratuberculosis (MAP) however exposure does not predetermine progression to clinical disease. The pathogenesis incorporates a subclinical phase during which MAP is capable of evading host immune responses through adaptation of host cellular immune mechanisms. Presented are results of transcriptomic analysis of Merino sheep experimentally exposed to MAP and repeatedly sampled over the subclinical phase, identifying genes consistently changed over time in comparison to unexposed controls and associated with different disease outcomes. MAP exposed sheep were classified as diseased 45% (n = 9) or resilient 55% (n = 11). Significant gene expression changes were identified in the white blood cells of paucibacillary (n = 116), multibacillary (n = 98) and resilient cohorts (n = 53) compared to controls. Members of several gene families were differentially regulated, including S100 calcium binding, lysozyme function, MHC class I and class II, T cell receptor and transcription factors. The microarray findings were validated by qPCR. These differentially regulated genes are presented as putative biomarkers of MAP exposure, or of the specified disease or resilience outcomes. Further, in silico functional analysis of genes suggests that experimental MAP exposure in Merino sheep results in adaptations to cellular growth, proliferation and lipid metabolism.

A Cell Proliferation and Inflammatory Signature Is Induced by Lawsonia intracellularis Infection in Swine

Lawsonia intracellularis causes porcine proliferative enteropathy. This is an enteric disease characterized by thickening of the wall of the ileum that leads to decreased growth of animals and diarrhea. In this study, we investigated the host response to L. intracellularis infection by performing transcriptomic and pathway analysis of intestinal tissue samples from groups of infected and noninfected animals at 14, 21, and 28 days postchallenge. At the peak of infection, when animals developed the most severe lesions, infected animals had higher levels of several gene transcripts involved in cellular proliferation and inflammation, including matrix metalloproteinase-7 (MMP7), transglutaminase-2 (TGM2), and oncostatin M (OSM). Histomorphology also revealed general features of intestinal inflammation. This study identified important pathways associated with the host response in developing and resolving lesions due to L. intracellularis infection.IMPORTANCE Lawsonia intracellularis is among the most important enteric pathogens of swine, and it can also infect other mammalian species. Much is still unknown regarding its pathogenesis and the host response, especially at the site of infection. In this study, we uncovered several novel genes and pathways associated with infection. Differentially expressed transcripts, in addition to histological changes in infected tissue, revealed striking similarities between L. intracellularis infection and cellular proliferation mechanisms described in some cancers and inflammatory diseases of the gastrointestinal tract. This research sheds important light into the pathogenesis of L. intracellularis and the host response associated with the lesions caused by infection.