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

Early Growth Response Factor 4 (EGR4) Expression in Gut Tissues and Regional Lymph Nodes of Cattle with Different Types of Paratuberculosis-Associated Lesions: Potential Role of EGR4 in Resilience to Paratuberculosis

Summary-data-based Mendelian randomization (SMR) analysis identified a novel cis-expression quantitative loci (cis-eQTL) associated with the upregulation of the expression of the early growth response factor 4 (EGR4) gene in animals with paratuberculosis (PTB)-associated multifocal lesions, which has been suggested to be modulating the NF-kβ-induced proinflammatory immune response to Mycobacterium avium subsp. paratuberculosis (Map) infection. To confirm these findings and to study the role of EGR4 expression in PTB resilience, the number of EGR4-expressing cells were analysed in paraffin-fixed gut tissues and regional lymph nodes of naturally Map-infected Holstein Friesian cows with focal, multifocal (subclinical and clinical), and diffuse lesions (intermediate and multibacillary), and in controls without lesions by quantitative anti-EGR4 immunohistochemistry. Subclinical animals with multifocal lesions showed a significantly higher number of EGR4-positive cells and were sacrificed at a significantly older average age than the remaining groups (p < 0.001 in all cases). We hypothesize that EGR4 could be mitigating the negative impact of Map infection on host clinical status through its involvement in three molecular mechanisms that promote resilience: (i) limiting NF-kβ-mediated proinflammatory responses, (ii) controlling tissue damage, acting as a brake on T-cell proliferation and cytokine production, and (iii) favouring tissue repair through interaction with epidermal growth factor receptor (EGFR).

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.

Genome-Wide Association Study Reveals Quantitative Trait Loci and Candidate Genes Associated with High Interferon-gamma Production in Holstein Cattle Naturally Infected with Mycobacterium Bovis

Mycobacterium bovis (Mb) is the causative agent of bovine tuberculosis (bTb). Genetic selection aiming to identify less susceptible animals has been proposed as a complementary measure in ongoing programs toward controlling Mb infection. However, individual animal phenotypes for bTb based on interferon-gamma (IFNɣ) and its use in bovine selective breeding programs have not been explored. In the current study, IFNɣ production was measured using a specific IFNɣ ELISA kit in bovine purified protein derivative (bPPD)-stimulated blood samples collected from Holstein cattle. DNA isolated from the peripheral blood samples collected from the animals included in the study was genotyped with the EuroG Medium Density bead Chip, and the genotypes were imputed to whole-genome sequences. A genome-wide association analysis (GWAS) revealed that the IFNɣ in response to bPPD was associated with a specific genetic profile (heritability = 0.23) and allowed the identification of 163 SNPs, 72 quantitative trait loci (QTLs), 197 candidate genes, and 8 microRNAs (miRNAs) associated with this phenotype. No negative correlations between this phenotype and other phenotypes and traits included in the Spanish breeding program were observed. Taken together, our results define a heritable and distinct immunogenetic profile associated with strong production of IFNɣ in response to Mb.

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.

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.

Association between High Interferon-Gamma Production in Avian Tuberculin-Stimulated Blood from Mycobacterium avium subsp. paratuberculosis-Infected Cattle and Candidate Genes Implicated in Necroptosis

The mechanisms underlying host resistance to Mycobacterium avium subsp. paratuberculosis (MAP) infection are largely unknown. In the current study, we hypothesize that cows with an ability to produce higher levels of interferon-gamma (IFNɣ) might control MAP infection more successfully. To test this hypothesis, IFNɣ production was measured using a specific IFNɣ ELISA kit in avian purified protein derivative (aPPD)-stimulated blood samples collected from 152 Holstein cattle. DNA isolated from peripheral blood samples of the animals included in the study was genotyped with the EuroG Medium-Density Bead Chip, and the genotypes were imputed to whole-genome sequencing. A genome-wide association analysis (GWAS) revealed that high levels of IFNɣ in response to the aPPD were associated with a specific genetic profile (heritability = 0.64) and allowed the identification of 71 SNPs, 40 quantitative trait loci (QTL), and 104 candidate genes. A functional analysis using the 104 candidate genes revealed a significant enrichment of genes involved in the innate immune response and, more specifically, in necroptosis. Taken together, our results define a heritable and distinct immunogenetic profile associated with the production of high IFNɣ levels and with the capacity of the host to lyse MAP-infected macrophages by necroptosis.

Association between the severity of histopathological lesions and Mycobacterium avium subspecies paratuberculosis (MAP) molecular diversity in cattle in southern Chile

The objective was to evaluate the association between the severity of histopathological lesions caused by Mycobacterium avium subspecies paratuberculosis (MAP) infection and the molecular diversity of this pathogen. Blood, ileum, and mesenteric lymph node samples were collected at slaughter, from 1,352 adult cattle [source population 1 (SP1)]. In addition, 42 dairy herds (n = 4,963 cows) were followed for 2 years, and samples from compatible paratuberculosis clinical cases [source population 2 (SP2)] were collected. MAP infection was confirmed using an ELISA test, liquid media culture, and PCR. Isolates were genotyped using five MIRU-VNTR markers. Tissues from confirmed samples were subjected to a histopathological examination. A histopathological severity score (HSS) system was developed and used to grade (0 to 5) the magnitude of lesions caused by MAP. In general, the HSS system assesses the number of foci and degree of macrophage infiltration, together with the presence of multinucleated giant cells (MGCs) and acid-fast bacilli (AFB), in addition to the fusion of the intestinal villi and hyperplasia of the crypts. Despite the large sampling effort, only 79 MAP isolates were successfully genotyped, where 19 different haplotypes were described. A mixed-effect Poisson regression model was used to assess the relationship between haplotypes and HSS values. The model was controlled by animal age, and the farm was used as a random effect. Haplotypes were grouped based on their relative frequency: the most frequent haplotype (group i, 49.4%), the second most frequent haplotype (group ii, 12.7%), and all other haplotypes (group iii, 37.9%). Model outputs indicated that group i had significantly higher HSS values than group iii. In addition, group i was also associated with higher optical density (OD) values of the ELISA test. These results support the existence of differences in pathogenicity between MAP haplotypes. However, results were based on a relatively small sample size; thus, these should be taken with caution. Despite this, study findings suggest that haplotypes would be associated with differences in disease progression, where the dominant haplotype tends to generate more severe lesions, which could be linked to a greater shed of MAP cells than non-dominant haplotypes, increasing their chances of transmission.

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.

Vitamin D3 alters macrophage phenotype and endosomal trafficking markers in dairy cattle naturally infected with Mycobacterium avium subsp. paratuberculosis

Macrophages are important host defense cells in ruminant paratuberculosis (Johne’s Disease; JD), a chronic enteritis caused by Mycobacterium avium subsp. paratuberculosis (MAP). Classical macrophage functions of pathogen trafficking, degradation, and antigen presentation are interrupted in mycobacterial infection. Immunologic stimulation by 25-hydroxyvitamin D₃ (25(OH)D₃) and 1,25-dihydroxyvitamin D₃ (1,25(OH)₂D₃) enhances bovine macrophage function. The present study aimed to investigate the role of vitamin D₃ on macrophage phenotype and endosomal trafficking of MAP in monocyte-derived macrophages (MDMs) cultured from JD-, JD+ subclinical, and JD+ clinically infected cattle. MDMs were pre-treated 100 ng/ml 25(OH)D₃ or 4 ng/ml 1,25(OH)₂D₃ and incubated 24 hrs with MAP at 10:1 multiplicity of infection (MOI). In vitro MAP infection upregulated pro-inflammatory (M1) CD80 and downregulated resolution/repair (M2) CD163. Vitamin D₃ generally decreased CD80 and increased CD163 expression. Furthermore, early endosomal marker Rab5 was upregulated 140× across all stages of paratuberculosis infection following in vitro MAP infection; however, Rab5 was reduced in MAP-activated MDMs from JD+ subclinical and JD+ clinical cows compared to healthy controls. Rab7 expression decreased in control and clinical cows following MDM infection with MAP. Both forms of vitamin D₃ reduced Rab5 expression in infected MDMs from JD- control cows, while 1,25(OH)₂D₃ decreased Rab7 expression in JD- and JD+ subclinical animals regardless of MAP infection in vitro. Vitamin D₃ promoted phagocytosis in MDMs from JD- and JD+ clinical cows treated with either vitamin D₃ analog. Results from this study show exogenous vitamin D₃ influences macrophage M1/M2 polarization and Rab GTPase expression within MDM culture.

Genome-wide association studies for the identification of cattle susceptible and resilient to paratuberculosis

Mycobacterium avium subsp. paratuberculosis (MAP) causes Johne's disease or paratuberculosis (PTB), with important animal health and economic implications. There are no therapeutic strategies to control this disease, and vaccination with inactivated vaccines is limited in many countries because it can interfere with the intradermal test used for bovine tuberculosis detection. Thus, infected animals either get culled after a positive ELISA or fecal PCR result or die due to clinical disease. In this study, we review recent studies aimed to discover genetic markers which could help to identify and select cattle less susceptible and more resilient to PTB. In recent years, the genotyping and subsequent imputation to whole-genome sequence (WGS) has allowed the identification of single-nucleotide polymorphisms (SNPs), quantitative trait loci (QTL), and candidate genes in the Bos taurus genome associated with susceptibility to MAP infection. In most of these genome-wide association studies (GWAS), phenotypes were based on ante-mortem test results including serum ELISA, milk ELISA, and detection of MAP by fecal PCR and bacteriological culture. Cattle infected with MAP display lesions with distinct severity but the associations between host genetics and PTB-associated pathology had not been explored until very recently. On the contrary, the understanding of the mechanisms and genetic loci influencing pathogen resistance, and disease tolerance in asymptomatic individuals is currently very limited. The identification of long-time asymptomatic cattle that is able to resist the infection and/or tolerate the disease without having their health and milk production compromised is important for disease control and breeding purposes.

A Genome-Wide Association Study for Tolerance to Paratuberculosis Identifies Candidate Genes Involved in DNA Packaging, DNA Damage Repair, Innate Immunity, and Pathogen Persistence

Although the genetic susceptibility to diseases has been extensively studied, the genetic loci and the primary molecular and cellular mechanisms that control disease tolerance are still largely unknown. Bovine paratuberculosis (PTB) is an enteritis caused by Mycobacterium avium subsp. paratuberculosis (MAP). PTB affects cattle worldwide and represents a major issue on animal health. In this study, the associations between host genetic and PTB tolerance were investigated using the genotypes from 277 Spanish Holstein cows with two distinct phenotypes: cases) infected animals with positive PCR and bacteriological culture results but without lesions in gut tissues (N= 24), and controls) animals with negative PCR and culture results but with PTB-associated lesions (N= 253). DNA from peripheral blood of the study population was genotyped with the Bovine EuroG MD Bead Chip, and the corresponding genotypes were imputed to whole-genome sequencing (WGS) data. A genome-wide association study was performed using the WGS data and the defined phenotypes in a case-control approach. A total of 142 single nucleotide polymorphisms (SNPs) were associated (false discovery rate ≤ 0.05, P values between 1.5 × 10^-7 and 5.7 × 10^-7) with tolerance (heritability= 0.55). The 40 SNPs with P-values < 5 × 10^-7 defined 9 QTLs and 98 candidate genes located on BTA4, BTA9, BTA16, BTA25, and BTA26. Some of the QTLs identified in this study overlap with QTLs previously associated with PTB, bovine tuberculosis, mastitis, somatic cell score, bovine diarrhea virus persistent infection, tick resistance, and length of productive life. Two candidate genes with important roles in DNA damage response (ERCC4 and RMI2) were identified on BTA25. Functional analysis using the 98 candidate genes revealed a significant enrichment of the DNA packaging process (TNP2/PRMI1/PRM2/PRM3). In addition, the TNF-signaling (bta04668; TRAF5/CREB5/CASP7/CHUK) and the toxoplasmosis (bta05145; TGFβ2/CHUK/CIITA/SOCS1) pathways were significantly enriched. Interestingly, the nuclear Factor NF-κβ Inhibitor Kinase Alpha (CHUK), a key molecule in the regulation of the NF-κB pathway, was enriched in both pathways. Taken together, our results define a distinct immunogenetic profile in the PTB-tolerant animals designed to control bacterial growth, modulate inflammation, limit tissue damage and increase repair, thus reducing the severity of the disease.