- 1Department of Animal Health, NEIKER-Basque Institute for Agricultural Research and Development, Basque Research and Technology Alliance (BRTA), Derio, Spain
- 2Departamento de Mejora Genética Animal, Instituto Nacional de Investigación y Tecnología Agraria y Alimentaria, CSIC, Madrid, Spain
- 3Departamento de Producción Agraria, Escuela Técnica Superior de Ingeniería Agronómica, Alimentaria y de Biosistemas, Universidad Politécnica de Madrid, Ciudad Universitaria, Madrid, Spain
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.
Introduction
Paratuberculosis (PTB) or Johne´s disease is a granulomatous enteritis of domestic and wild ruminants caused by Mycobacterium avium subsp. paratuberculosis (MAP). PTB is a major problem for animal health and must be reported to the World Organization for Animal Health. In Europe and North America, PTB is endemic in dairy cattle, with herd prevalence estimates higher than 50% (1). PTB causes substantial economic losses to the dairy industry mainly due to decreased milk production, weight loss, and premature culling or death of the infected animals (2). Animals are infected early in life through the fecal-oral route but clinical onset appears when animals are 18 months or older. Current diagnostic tests, ELISA and fecal PCR, often fail to detect MAP-infected cattle early in the course of infection (3, 4). According to their extension, cellular infiltrate, and amount of MAP, PTB-associated histological lesions have been classified into the following categories: focal, multifocal, and diffuse (diffuse paucibacillary, diffuse intermediate, and diffuse multibacillary) (5). High bacterial burden, clinical signs, and gross lesions are associated with the presence of diffuse lesions. MAP has been associated with Crohn´s disease (CD) in humans and has been detected in samples of patients with CD, ulcerative colitis, and idiopathic inflammatory bowel disease (IBD)-associated colorectal cancer (6–8). MAP is also a possible trigger factor in some human autoimmune diseases such as rheumatoid arthritis, multiple sclerosis, and type I diabetes (9–12). Currently, no efficacious vaccine exists that confers protection against MAP. The current control measures include the removal of MAP-infected animals from the herd and the enhancement of farm-biosecurity measures. Genetic selection to enhance the resistance of dairy cattle to PTB and other diseases is extensively being explored (13, 14).
In previous studies, we identified a total of 380 single nucleotide polymorphisms (SNPs) associated (FDR ≤ 0.05, P < 5 × 10-7) with PTB susceptibility using serum ELISA for the detection of humoral responses against MAP and post-mortem diagnostic definitions (tissue PCR, tissue bacteriological culture, and histopathological analysis) in a common set of Spanish Holstein cattle (N= 983) (15, 16). These studies showed that PTB susceptibility is under complex genetic control, with heritability (h2) estimates ranging from 0.139 for the combination of ELISA-PCR-culture positive results to 0.189 for the risk of developing diffuse lesions. Using these phenotypes, we were, however, unable to identify SNPs associated with PTB resistance. While test positive cows are very likely infected, animals can test negative due to the lack of sensitivity of the diagnostic method or lack of exposure to MAP. Host defense strategies against infectious diseases are comprised of both host resistance and disease tolerance. Resistance is defined as “the ability of the host to prevent invasion (i.e., absence of a target receptor) or to clear the pathogen at the early stage of infection by mounting a protective immune response” (17). Disease tolerance is defined as “the mechanisms that decrease host susceptibility to tissue damage, or other fitness costs caused by pathogens or by the immune response” (18). Unlike resistance, disease tolerance does not necessarily imply direct effects on the pathogen load. Initially, it was observed that plants with a distinct genetic profile could tolerate an infection without affecting the pathogen load, which was termed “disease tolerance” (19, 20). More recently, the crucial role of disease tolerance in invertebrates and vertebrates against infectious diseases has been recognized (21, 22). The term disease tolerance is not to be confused with immunological tolerance, which is defined as unresponsiveness to self-antigens (18). Evidence of tolerance can potentially be observed in the host genome, with genetic variants that confer disease tolerance displaying positive or advantageous selection (23). However, the genetic loci influencing tolerance induction and the primary molecular and cellular mechanisms underlying disease tolerance in asymptomatic individuals are still largely unknown.
Mycobacterium tuberculosis and MAP have coevolved with humans and animals infrequently compromising host survival. Although this has been attributed to host resistance, our understanding of the mechanisms associated with disease tolerance in the infected individuals who become asymptomatic and disease-free is very limited. Since disease tolerance does not affect the pathogen burden as resistance is expected to do, it does not provoke the pathogen to develop new strategies to evade the host immune response. In addition, pathogen prevalence induces a continuous selection in favor of beneficial genetic variants for disease tolerance. While improving host resistance could lead to disease eradication, this is unlikely if hosts are tolerant as they can harbor the pathogen without showing clinical symptoms and lesions (24). For endemic diseases that result in a high number of asymptomatic animals, it might be more beneficial to improve tolerance than resistance. Furthermore, whereas resistance mechanisms are often pathogen-specific (e.g. mobilization of specific immune cells), tolerance mechanisms that prevent or repair damage may offer protection for a wide range of pathogens. Ultimately, tolerant mechanisms determine the overall health of the host and its longevity. Mycobacteria have developed several strategies to evade the resistance mechanisms of the host and persist for long periods. Clinical signs associated with Mycobacterial infections emerge as the host homeostasis becomes compromised due to tissue damage. In this context, the induction and modulation of tolerance pathways are critical to control inflammation and minimize pathology.
The response to MAP infection is complex and heritable leading to differences between individuals which can be categorized as (i) susceptible host, individuals who progress to severe forms of the disease, (ii) resistant host, individuals able to eliminate the bacteria by inducing protective innate immune responses at the early stage of infection, (iii) tolerant host, if innate immunity is unable to eliminate MAP, the host initiates tolerance mechanisms to contain MAP infection and prevent tissue damage. In the current study, we searched for genetic loci associated with tolerance to PTB by using WGS data from infected Spanish Holstein cows with MAP detected by PCR and bacteriological culture but without lesions in gut tissues and regional lymph nodes (N=24). A PTB tolerant animal was PCR and culture-positive (infected) but displayed no lesions in gut tissues (no disease). Using this case population, three case-control approaches were designed including as control population: 1) cows with negative PCR and culture results and with PTB-associated focal lesions (N= 253), 2) cows with PTB-associated focal lesions irrespective of their PCR and culture result (N= 288), and 3) cows with negative ELISA-PCR and culture results and without lesions in gut tissues (N= 373). Variance components, heritability (h2) estimates, SNPs, quantitative trait loci (QTL), and candidate genes were calculated for the three case-control approaches. Subsequently, functional pathways were identified for the case-control study with the highest h2 estimate, the case-control 1 (h2 = 0.55). Estimated breeding values (EBVs) for this population were calculated and cross-validated. Since there is evidence that some allelic variants may contribute to tolerance to multiple diseases, the identified QTLs and candidate genes were compared with QTLs and candidate genes for bovine diseases, and with human candidate genes for CD, IBD, colorectal cancer, multiple sclerosis, type I diabetes mellitus, and rheumatoid arthritis.
Materials and Methods
Ethics Statement
Animals used in this study were not submitted to any in vivo experimentation and, therefore, no specific ethics committee authorization was required. The cows were slaughtered in the Bilbao and Donostia municipal slaughterhouses (Basque Country, Spain) from March 2007 to May 2010 under the pertinent Basque (Basque Government Decree 454/1994), Spanish (Spanish Government Law 32/2007 and Royal decree 731/2007), and European (Council Regulation No 1099/2009) legislation on animal welfare.
Animals
The study population consisted of 685 culled Spanish Holstein cattle from several herds located on seven regions: Basque Country (N= 286, 41.75%), Catalonia (N= 157, 22.92%), Navarre (N= 147, 21.46%), Cantabria (N= 49, 7.15%), Aragon (N= 22, 3.21%), Castile and Leon (N= 16, 2.34%), La Rioja (N= 5, 0.73%) and Asturias (N=3, 0.44%). All cows were 18 months years or older (5.6 years mean age).
Histopathological Examination, Tissue PCR and Bacteriological Culture, and ELISA
Post-mortem tissue sampling, PCR, bacteriological culture, and ELISA were previously performed (25). For the histopathological analysis of gut tissues standardized protocols were used (5). Briefly, samples from the ileocecal lymph node, jejunal lymph node, ileocecal valve (ICV), jejunum, and terminal ileum were collected from each animal and fixed in 10% neutral buffered formalin for 72 h, dehydrated, embedded in paraffin, and cut into 4 μm sections. Several sections were mounted on microscope slides, stained with hematoxylin and eosin (HE) and with Ziehl-Neelsen (ZN), and examined. PTB-associated lesions were classified in focal, multifocal, and diffuse lesions as previously described (5). More specifically, PTB-associated focal lesions were defined according to their morphology and localization. They appear mainly in the lymph nodes (ileal and jejunal), characterize by the presence of small granulomas formed by 5 to 30 macrophages with undetected or low MAP load. Since similar focal lesions have been reported in experimental MAP infections and in the early stages of natural MAP infections in cattle (26–28), they were considered to represent initial forms of PTB. A pool (2 gr) of ICV, distal ileum and the jejunal caudal lymph nodes was decontaminated with 38 ml of hexadecyl pyridinium chloride at a final concentration of 0.75% (Sigma, St. Louis, MO) and homogenized in a stomacher blender. After 30 min of incubation at room temperature, 15 ml of the suspension was transferred to a new tube and incubated overnight for decontamination and sedimentation. Approximately, 200 µl of the suspension was taken and inoculated in duplicate Herrold´s egg yolk (Becton Dickinson, Franklin Lakes, NJ, USA) and Lowenstein-Jensen media (Difco, Detroit, MI, USA), both supplemented with 2 mg/L of mycobactin J (Allied Monitor, Fayette, MO, USA). If one or more MAP colonies were observed in any of the four medium slants the tissue culture was considered positive. Bacterial load in tissues was classified as low (<10 CFU/g), medium (between 10 to 50 CFU/g), or heavy (>50 CFU/gr). A second aliquot from the same tissue homogenates was used for DNA isolation and MAP IS900 amplification (Adiagene, Saint Brieuc, France). PCR amplifications were performed on an ABI Prism 7000 Sequence Detection System (Applied Biosystems, Foster City, CA, US). Samples with a threshold cycle (Ct) below 40.0 were considered positive. Antibody production against MAP in serum samples was tested using a two-step ELISA PTB antibody screening and verification kit (IDEXX Laboratories, Inc., Westbrook, ME, USA) according to the manufacturer´s instructions. The results were categorized as positive according to the sample-to-positive control ratio defined by the manufacturer.
Genotyping and Imputation
DNA was extracted from peripheral blood samples collected from the study population at the time of slaughter and genotyped with the EuroG MD Bead Chip at the molecular genetic laboratory service of the Spanish Federation of Holstein Cattle (CONAFE). The InfiniumTM iScan software (Illumina, San Diego, CA) was used for allele assignation. The individual genotypes were previously phased and imputed to WGS level as previously described (15). Briefly, individual genotypes were phased using Eagle 2.4 (29) and imputed with minimac4 (30) to the Bovine HD Bead Chip using a reference panel of 1,278 Bos taurus from Run7.0 of the 1,000 Bull Genomes project and 581,712 SNPs (ASR-UCD1.2). Imputation to WGS level was then undertaken using a reference population of 2,333 Bos taurus from Run6.0 of the 1,000 Bull Genomes project (31). All the SNPs passed a call rate > 0.80. After filtering out SNPs with minimum allele frequency (MAF) < 0.01, we obtained 13,881,067 SNPs.
Case-Control Studies, Variance Components, Standard Errors (SE), and h2 Estimates
A PTB tolerant animal was defined as an animal naturally infected with MAP that could harbor the pathogen without showing lesions in gut tissues that could be observed after the histopathological analysis of gut tissues. More specifically, animals with positive PCR and culture results but without PTB-associated lesions in gut tissues (N= 24). Using this case population, three case-control approaches were designed including as control population: 1) cows with negative PCR and culture results and with PTB-associated focal lesions (N= 253), 2) cows with PTB-associated focal lesions irrespective of their PCR and culture result (N= 288), and 3) cows with negative ELISA-PCR and culture results and without lesions in gut tissues (N= 373). The variance components, standard errors (SE), and h2 estimates explained by all the SNPs of each case-control population were calculated using the genome-wide complex trait analysis (GCTA) software 1.93.2, according to the following formula:
where is the variance explained by all the SNPs and is the residual variance (32).
Genome-Wide Association Study (GWAS)
In the current study, imputed genotypes (15) and the results of the diagnostic tests (25) were analyzed using the mixed linear model association analysis of the GCTA 1.93.2 software (32). Briefly, the model is y = a + bx + g + e, where y is the phenotype, a is the mean term, b is the additive effect (fixed effect) of the candidate SNP to be tested for association, x is the SNP genotype indicator variable coded as 0, 1 or 2, g is the polygenic effect (random effect) assumed to be distributed as N~(0, ), and e is the residual assumed to be distributed as N~(0, ). Age was included as a covariate in the analysis. After the GWAS, the SNPs with R2 values higher than 70% were retained. To account for multiple testing, a 5% genome-wide false discovery rate (FDR) was used. As strong evidence for association, a level of significance (P < 5 × 10-7) appropriate for genome-wide analysis was established according to the Welcome Trust Case Control Consortium (33). The inflation factor (λ) and quantile-quantile plots were calculated to compare the observed distributions of –log (P-values) to the expected distribution under the no association model. The regression coefficients (b) calculated using the GCTA 1.93.2 software represents “the estimated increase in the log odds of the outcome per unit increase in the value of the exposure. In other words, the exponential function of b is the odds ratio (OR) associated with a one-unit increase in exposure” (34). In addition, the OR and their 95% confidence intervals (CI) for the SNPs associated with PTB tolerance (P < 5 × 10-7) were calculated using logistic regression analysis with the WGassociation function of SNPassoc 1.9.2 (35).
SNPs, QTLs, and Candidate Genes Identification
Localization of SNPs and QTLs and identification of candidate genes was performed using the ARS-UCD1.2 reference genome as previously described (15, 16). The genomic localization of the identified SNPs was determined using the Ensembl Variant Effect Predictor (VEP). QTLs were defined based on SNPs on linkage disequilibrium (LD) patterns with SNPs that surpassed the threshold (P < 5 × 10-7) in a given chromosome. Pairwise LD (r2 values) among all SNPs within 500,000 base pairs upstream and downstream of the suggestive SNP were calculated using PLINK 1.9. The beginning and end of each QTL were defined by the SNPs that were furthest upstream and downstream of the suggestive SNP and that had an r2 ≥ 0.9 (36, 37). Overlapping QTL regions were merged and considered as a single QTL. The defined QTL regions we further investigated for the presence of candidate genes within 50,000 base pairs to each side of the defined QTL using Ensembl (https://www.ensembl.org). The identified QTLs and candidate genes were compared with reported QTLs for bovine diseases (http://www.animalgenome.org) and with human candidate genes previously identified for CD, IBD, colorectal cancer, multiple sclerosis, type I diabetes mellitus, and rheumatoid arthritis (https://www.ebi.ac.uk/gwas). To further investigate the function of the identified candidate genes, we uploaded the 98 genes into the Innate DB database for innate immunity genes (https://innatedb.com).
Gene Ontology (GO) and Pathway Enrichment Analysis
Candidate genes were investigated for significant enrichment of GO categories and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways using the Cluster Profiler Bioconductor package 3.10.1 (38). To account for multiple testing, P-adjust values were estimated.
Estimated Breeding Values (EBVs) and Cross-Validation of the Genomic Predictions
The GBLUP (genomic best linear unbiased prediction) model (39) of GCTA 1. 93. 2. was used to estimate the breeding value of each individual of the study population attributed to the aggregative effect of the SNPs with evidence of association with the phenotype; P < 5 × 10-5 and P < 5 × 10-7, respectively. The accuracy of the genomic predictions was cross-validated in a validation population containing animals with a tissue culture positive result and without lesions in gut tissues (cases, N= 10) and animals with PTB-associated lesions in gut tissues (controls, N= 190). Cross-validation was performed using a 5-fold strategy. Briefly, phenotypes were masked (i.e., set to missing) for all the animals in a random 20% of the validation population and both phenotypes and genotypes for all other animals (i.e. training data set) were used to predict EBVs of individuals from the omitted group (i.e. testing dataset). This cross-validation process was iterated five times, such that each animal´s phenotype was masked at least once to generate its EBV. The ability of the group to determine whether or not each masked animal has a PTB tolerant phenotype was measured by the area under the receiver operating curve (AUC) for each model in each simulation using the ROCR package in R (40). The ROC curve is a graphical plot of the sensitivity, or true vs. false positive rate (1-specificity). The AUC would be 0.5 if its discrimination changes based on random chance, whereas if the model is perfectly able to determine the phenotype of each animal, the area under the curve would be 1.0 (41).
Results
Data, Variance Components, and h2 Estimates
The Spanish Holstein cows included in our study were previously tested by tissue PCR and bacteriological culture, ELISA, and histopathological analysis of gut tissues (25). Genotypes from the study population were previously obtained and imputed to WGS (15). In the current study, the associations between genotypes and PTB tolerance were investigated by using the genotypes from infected cows with a positive PCR and culture result and without lesions in gut tissues and regional lymph nodes (N= 24). All these animals had a low bacterial load in gut tissues (< 10 CFU/g) and a negative ELISA result. Using this case population, three case-control approaches were designed using as control population 1) cows with negative PCR and culture results and with PTB-associated focal lesions (N= 253), 2) cows with PTB-associated focal lesions irrespective of their PCR and culture result (N= 288), and 3) cows with negative ELISA-PCR and culture results and without lesions in gut tissues (N= 373). Table 1 summarizes the number of cases and controls, variance components, and h2 estimates, for each of the three case-control approaches. The h2 estimates were 0.55, 0.44 and 0.18 for case-control studies 1, 2, and 3, respectively (Table 1). The small differences in h2 estimates (0.55 and 0.44) between the case-control studies 1 and 2 suggest that there would be little bias due to the use of either sample as the control group. The highest h2 estimate (h2 = 0.55) was obtained when cows with negative PCR and culture results and with PTB-associated focal lesions (N= 253) were the control population, highlighting the importance of the selection of control individuals from the lower extremity of the relevant trait distribution to increase sensitivity. The b-values ranged from 0.57 to 0.94 for the first case-control study and the range was 0.27-0.59 for the second case-control study providing more evidence that the most powerful approach was the first case-control study.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/820965/fimmu-13-820965-HTML/image_m/fimmu-13-820965-t001.jpg)
Table 1 Number (N) of cases and controls, variance components, standard errors (SE), and h2 estimates.
GWAS
GWAS using the imputed WGS datasets (15) and the phenotypes defined in the three case-control approaches were performed. The number of SNPs that surpassed the P < 5 x 10-7 threshold is presented in Table 1. Three SNPs on BTA25 (11.89-12.93 Mb) overlapped with the three case-control approaches suggesting an important link of this BTA25 region with the PTB tolerant animals. Twenty-eight SNPs were associated with the case-control studies 1 and 2, and four with the case-control studies 1 and 3.
A total of 40 (P < 5 x 10-7) SNPs were associated with PTB tolerance defined using the case-control approach with the highest hereditability (h2 = 0.55). In this approach, animals were considered cases when they were PCR and culture positives but they did not show lesions in gut tissues (N= 24). Controls were PCR and culture negatives with focal PTB-associated lesions in gut tissues (N= 253). A Manhattan plot showing –log10 (P-values) of association between every single SNP and PTB tolerance is presented in Figure 1A. The 40 SNPs associated with tolerance (FDR ≤ 0.05, P < 5 × 10-7) were located on the Bos taurus chromosomes BTA4, BTA9, BTA16, BTA25, and BTA26. Most of these SNPs (59%) were located in intronic regions (Figure 1B). A quantile-quantile plot was generated (Figure 1C). The plot showed a distribution close to the expected line (λmedian= 1.03), indicating that significant values were not overestimated due to population stratification or cryptic relatedness. The regression coefficients (b-values) for the 40 SNPs associated with PTB tolerance (P < 5 × 10-7) were all positive (data not shown). In addition, the OR and 95% confidence intervals (CI) calculated for the 40 SNPs were > 1 indicating that the animals with minor alleles had a higher probability of being tolerant (Table 2). The two SNPs with the highest OR (25.89 and 20.71) were located on BTA4 and BTA9. Interestingly, three of the 12 SNPs with OR= 15.86 and located on BTA25 (11.89-12.93 Mb) were common to the three case-control approaches.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/820965/fimmu-13-820965-HTML/image_m/fimmu-13-820965-g001.jpg)
Figure 1 GWAS for PTB tolerance. (A) Manhattan plots showing –log10 (P-values) of association between every single SNP and phenotype. The associations between host genetic and tolerance were investigated using the genotypes from 277 Spanish Holstein cows with two distinct phenotypes: cases) positive PCR and culture result without lesions (N = 24), and controls) negative PCR and culture result with PTB-associated lesions (N = 253). Each dot represents the result from the test association for a single SNP. The horizontal red line shows –log10 (5 × 10-7) and the blue line –log10 (5 × 10-5). (B) Genomic distribution of the SNPs associated (FDR≤ 0,05, P < 5 × 10-7) with tolerance according to the Ensembl Variant Effect Predictor (VEP). (C) Quantile-quantile plot compared the observed distribution of –log (P-values) to the expected values under the null hypothesis.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/820965/fimmu-13-820965-HTML/image_m/fimmu-13-820965-t002.jpg)
Table 2 Frequency of cases and controls and odds ratio (OR) under the additive genetic model for the 15-top SNPs associated with tolerance (FDR ≤ 0.05, P < 5 × 10-7).
SNP Variants, QTLs, and Candidate Genes Identification
The SNPs surpassing the significance threshold (P < 5 × 10-7), P-values, and candidate genes located within the defined QTLs are presented in Table 3 for the case-control 1, and in Supplementary Tables 1, 2 for the case-control studies 2 and 3. Five of the nine QTLs identified in the case-control study 1 and located on BTA4, BTA9, BTA25 (8.90-9.99 Mb, 11.89-12.93 Mb), and BTA26 were associated with the case-control study 2 as well. Case-control studies 1 and 3 shared two QTLs on BTA16 (71.51-72.51 Mb) and BTA25 (11.89-12.93 Mb). According to the Welcome Trust Consortium “a study design using hypercontrols (that is, selection of control individuals from the lower extremity of the relevant trait distribution) would generally be the most powerful approach in a study focusing on one disease” (33). Therefore, we focused on the results of the case-control study 1 with the highest h2 estimate (0.55) and b-values (0.57-0.90). No differences in age between the cases (74.73 ± 38.13 months) and controls (68.29 ± 25.79 months) (P= 0.2659) included in the case-control study 1 were observed.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/820965/fimmu-13-820965-HTML/image_m/fimmu-13-820965-t003.jpg)
Table 3 QTLs surpassing the significance threshold (P < 5 × 10-7) for evidence of an association with tolerance.
The 40 SNPs associated with PTB tolerance resided within 9 QTLs on 5 chromosomes including BTA4, BTA9, BTA16, BTA25, and BTA26 (Table 3). Twenty-five SNPs (62.5%) were located on BTA25 (8.90-12.93 Mb) highlighting the important association of this region with the phenotype. The BTA25 and BTA26 harbored 4 of the 9 identified QTLs. The QTL that harbored the most genome-wide significantly-associated SNP (P= 1.58E-07) was located on BTA16 (21.30-22.30 Mb). By examining the available cattle QTL database (Table 4), we observed that this QTL together with a QTL located in BTA26 (33.90-34.90 Mb) overlap with two QTLs previously reported in the literature that were associated with PTB susceptibility; QTL:14860 and QTL:14880, respectively (42). The largest QTL was located on BTA25 (8.90-9.99 Mb) and contained 13 SNPs and 16 candidate genes. This QTL overlap with QTLs associated with clinical mastitis (QTL:1751) (43) and somatic cell score (QTL:4993, QTL:52494) (44, 45). A second QTL identified on BTA25 (11.89-12.93 Mb) overlap with QTLs associated with bovine tuberculosis susceptibility (QTL:96564) (46), clinical mastitis (QTL:1751) (43), and somatic cell score (QTL:177968, QTL:32484, QTL:32485, QTL4993, QTL:52521) (45, 47, 48). The two QTLs identified in the current study on BTA16 and one of the QTLs identified on BTA26 (20.05-21.05 Mb) overlap with QTLs associated with length of productive life highlighting the association of this trait with PTB and health (45). Interestingly, a QTL identified in our study on BTA26 (33.90-34.90 Mb) overlaps with a QTL associated with bovine diarrhea virus persistent infection (QTL:66096) (49).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/820965/fimmu-13-820965-HTML/image_m/fimmu-13-820965-t004.jpg)
Table 4 Overlapping of the identified QTLs with regions associated with other cattle diseases and relevant traits.
The candidate genes associated with MAP tolerance that were identified in the current study had not been associated with bovine PTB, tuberculosis, or mastitis before. Table 3 presents a list of the 98 candidate genes. We identified candidate genes in all the QTLs. Only one candidate gene, the Excision Repair 4 (ERCC4), was identified on one of the QTLs located on BTA25 (11.89-12.93 Mb). ERCC4 is a DNA repair endonuclease (XPF). Together with ERCC1, XPF forms the ERCC1-XPF enzyme complex that participates in the nucleotide excision DNA repair (NER) pathway, a highly conserved mechanism that recognizes and removes helical distortions through the genome.
Functional GO and Metabolic Pathways
We identified 7 GOs and 2 metabolic pathways significantly enriched in the PTB tolerant animals. As seen in Table 5, two biological processes (BPs) related to DNA packaging and DNA conformation change were significantly enriched (P ≤ 0.05). We observed that five cellular compartments (CC) associated with DNA packaging such as nucleosome, chromatin, chromosome, DNA packaging complex, and protein-DNA complex were also enriched. Gene members of the Transition protein 2 (TNP2) and Protamine (PRM) family such as PRM1, PRM2, and PRM3 were common to these seven GOs. TNP2 plays a key role in the replacement of histones to PRMs during spermatogenesis and PRMs compact DNA into a highly condensed, stable, and inactive complex. Pathway analysis on the 98 candidate genes revealed a significant enrichment of the toxoplasmosis (bta05145; TGFβ2/CHUK/CIITA/SOCS1) and TNF-signaling (bta04668; TRAF5/CREB5/CASP7/CHUK) pathways with four proteins matching each route. Interestingly, the nuclear Factor NF-κβ Inhibitor Kinase Alpha (CHUK or IkBKA), a key molecule in the NF- κβ signaling pathway, was enriched in both pathways (P ≤ 0.01).
Estimated Breeding Values (EBVs) and Cross-Validation of the Genomic Predictions
The GBLUP model (39) was used to estimate breeding values of each individual in the study population attributed to the aggregative effect of the 142 and 40 SNPs with evidence of association with the phenotype; P < 5 × 10-5 and P < 5 × 10-7, respectively. The genomic predictions were cross-validated using a 5-fold strategy as described in Materials and Methods. The efficacy of the models was examined by ROC analysis, which provides a measure of the probability of the model´s successful classification of samples (Figure 2). The AUC (± SD) of the models generated using the 142 (AUC= 0.657 ± 0.017) and 40 (AUC= 0.676 ± 0.014) SNPs were very similar which suggested that the model developed with the 40 most significant SNPs is a good classifier. The cut-off point for the EBVs in the study population was obtained with the Statistical Analysis in Epidemiology (Epi) R package (cut-off= 0.0356).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/820965/fimmu-13-820965-HTML/image_m/fimmu-13-820965-g002.jpg)
Figure 2 Receiver Operating Characteristics (ROC) curve for application of genomic predictions from the study population to the validation population. GBLUP models were developed using the 142 (A) and 40 (B) SNPs with evidence of association with the case-control population 1; P < 5 × 10-5 and P < 5 × 10-7, respectively. Genetic predictions were cross-validated in the validation population using a five-fold strategy as described in Materials and Methods. Each curve represents one of the five generated models and models might occasionally overlap. The average AUC and standard deviation (SD) across the five cross-validations are included in each figure.
Discussion
Understanding tolerance pathways is a significant challenge due to the difficulty in determining the tolerance phenotype from field data. The identification of genetic variants underlying variation in tolerance to pathogens requires both great genome coverage and measuring disease outcomes that sometimes are not readily measured. In this study, we aim to elucidate cellular and molecular mechanisms responsible for disease tolerance using cattle naturally infected with MAP that can harbor the pathogen without showing lesions in gut tissues. Although a previous GWAS identified loci associated with tolerance to Johne´s disease (50), our GWAS is the first to complete an analysis of the genetic markers associated with tolerance to PTB using WGS data and epidemiological data based on three diagnostic tests. Results from the current study do not corroborate those of Zanella et al. (50) who identified a tolerance-associated region of 6.5 Kb on BTA15 which we did not identify here. This might be due to differences in the phenotypes used. In the study by Zanella et al, tolerance measured the relationship between infection intensity (level of MAP in the tissues) and fitness (level of MAP fecal shedding), but to define tolerance is critical to measure disease outcomes, not just infection burden. In our study, PTB tolerant cows are infected animals with positive PCR and bacteriological culture results but without lesions in gut tissues. Histopathological analysis usually detects more infected animals than bacteriological culture and PCR (5). However, in the study by Gonzalez et al, 3/116 (2.5%) of the animals were PCR and culture positives but without lesions in gut tissues. Similarly, in our study population (N=983) (15), 24 animals without lesions in gut tissues were PCR and culture positives (2.4%). The highest h2 estimate (h2 = 0.55) was obtained when cows with negative PCR and culture results and with focal lesions (N= 253) were the control population, highlighting the importance of the selection of control individuals from the lower extremity of the relevant trait distribution to increase sensitivity. In a large number of PTB-associated focal lesions, MAP could not be detected due to the low numbers of bacteria that cannot be easily detected by bacteriological culture or PCR (25–28, 51, 52). In humans, inflammatory bowel disease is solely defined by pathological findings since no causal pathogen has been unanimously recognized. Despite the small sample size of cases (N=24), the Manhattan plot showed loci with associations with the underlying SNPs. A total of 40 novel SNPs (FDR ≤ 0.05 and P < 5 × 10-7) defining 9 QTLs and 98 candidate genes were associated with the tolerant cattle. All these SNPs had positive b-values and OR > 1 and were, therefore, associated with PTB tolerance. EBVs were estimated for the study population using the SNPs with a significant association with the phenotype (FDR ≤ 0.05) and the genetic predictions were subsequently cross-validated in a different population.
The 98 candidate genes identified in our study were compared with candidate genes previously reported in CD, IBD, colorectal cancer, multiple sclerosis, type I diabetes mellitus, and rheumatoid arthritis. Seven genes including Glutamic-Oxaloacetic Transaminase 1 (GOT1), NK2 Homeobox 3 transcription factor (NKX2-3), Solute Carrier Family 25 Member 8 (SLC25A28), Fibroblast Growth Factor Receptor 1 Oncogene Partner (CPE43), Chemokine Receptor-Like 3 (CCR6), C-Type Lectin Domain Containing 16A (CLEC16A), and RecQ Mediated Genome Instability 2 (RMI2), were previously associated to rheumatoid arthritis. The last two genes, CLEC16A and RMI2, were also associated with CD, IBD, multiple sclerosis, and type I diabetes mellitus. CLEC16A is a regulator of mitophagy/autophagy and mitochondrial health and RMI2 plays a critical role in homologous recombination-dependent DNA repair (53). NKX2-3 and SLC25A28 were associated with CD, IBD, colorectal cancer, type I diabetes mellitus, and rheumatoid arthritis. Five of the 98 candidate genes have relevant innate immunity functions including CCR6, the nuclear Factor NF-κβ Inhibitor Kinase Alpha (CHUK or IkBKA), STAT-Induced STAT Inhibitor 1 or suppressor of cytokine Signaling 1 (SOCS1), TNF Receptor Associated Factor 5 (TRAF5), and caspase 7 (CASP7), a caspase involved in the activation cascade of caspases responsible for apoptosis.
Although poorly understood, the mechanisms underlying disease tolerance appear to involve conserved stress and damage responses (54). When these responses fail, the default response becomes programmed cell death of infected cells which is coupled to their immediate phagocytosis by macrophages and subsequent polarization. This macrophage polarization is also associated with the production of cytokines such as transforming growth factor β1 (TGFβ1) which can promote tissue regeneration. Consistent with these results, our pathways analysis with the 98 identified candidate genes revealed a significant enrichment of the DNA packaging (TNP2/PRM1/PRM2/PRM3), TNF-signaling pathway (TRAF5/CREB5/CASP7/CHUK), and toxoplasmosis (TGFβ2/CHUK/CIITA/SOCS1) in the tolerant animals (Figure 3). Upstream of the TNP2, PRM1, PRM2, and PRM3 genes located on BTA25, we identified the RMI2 gene which is critical for DNA damage repair and genome stability (53). Moreover, we identified 12 genetic variants on BTA25 (11.89-12.93 Mb) upstream of the Excision Repair 4 (ERCC4) gene. ERCC4 forms the ERCC1-XPF enzyme complex that participates in the nucleotide excision DNA repair (NER) pathway. Altogether, these findings provide evidence for a link between PTB tolerance and genetic variants in candidate genes involved in DNA packaging, DNA damage repair response, innate immune responses, and pathogen persistence.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/820965/fimmu-13-820965-HTML/image_m/fimmu-13-820965-g003.jpg)
Figure 3 Schematic diagram illustrating the identified PTB-tolerance mechanisms. Candidate genes involved in DNA packaging, TNF signaling pathway, and toxoplasmosis were associated with the PTB tolerant individuals. These mechanisms play conserved and important roles in disease tolerance.
In most eukaryotic cells, nucleosomes are the basic structural units of chromatin formed by a histone octamer core. A unique chromatin remodeling process occurs extensively during spermatogenesis, which leads to more than 90% of core histones being replaced by transition proteins (TNP), and then PRMs (55, 56). The number of histones that are replaced by PRMs and the proportion of PRM1/PRM2 are essential to keep chromatin structure. So if this proportion changes, the DNA is poorly packed and susceptible to damage (57). The specialized structure of the sperm chromatin-PRM has a dual function; first, to protect the DNA from damage during transport to the oocyte, and then to enable a rapid and complete unpacking of the undamaged paternal genome in the ooplasm (58). In the ooplasm, maternal factors rapidly access the paternal chromatin and begin the replacement of PRM with maternal histones (59). Any damage to the sperm DNA during the transition from the testicle to the oocyte cannot be repaired until the DNA is accessible for DNA repair systems in the ooplasm (58). On the other hand, oocytes are especially sensitive to DNA damage when they resume meiosis during the menstrual cycle. These DNA lesions, if not repaired correctly, can endanger the fitness and viability of the cell. Since DNA damage is one of the most important factors that can compromise the viability of the cell and/or organism and contribute to diseases, it is well recognized that individuals with enhanced DNA packaging and DNA damage responses have improved homeostasis and prolonged healthy lifespan (60).
Schneider et al., reported that mutants in the fly TNF–related molecule “Eiger” when challenged with a bacterial species, lived longer with similar pathogen titers (61). Similarly, we found a significant enrichment of candidate genes involved in the TNF-signaling pathway (TRAF5/CREB5/CASP7/CHUK) in the PTB tolerant animals. This pathway can trigger programmed cell death while also activating pro-survival responses via activation of NF-κβ. TRAF5, CASP7, and CHUK or IkBKA, have well-recognized roles in the innate immune responses. CHUK is a serine kinase that plays an essential role in the canonical NF-κβ signaling pathway. CHUK acts as part of the canonical IKK complex phosphorylating inhibitors of the NF-κβ on serine residues. In turn, free NF-κβ is translocated into the nucleus and activates the transcription of genes involved in immune response, growth control, or protection against apoptosis. CASP7 is responsible for apoptosis, which is regarded as an important mechanism through which the host can restrict pathogen replication in infected cells and maximize disease tolerance. TRAF5 is a member of the TNF receptor-associated factor (TRAF) family that mediates TNF-induced activation. TNF-signaling activates NF-κβ and JNK pathways.
Although the TNF-signaling pathway constitutes a critical host defense against Mycobacteria, its excess is also implicated in immunopathogenesis and activation of programmed necrosis pathway via induction of mitochondrial reactive oxygen species (ROS) that kills both bacteria and infected macrophages (62). Therefore, the innate immune response, while important to clear the pathogen, can also cause tissue damage. In addition to the TNF pathway, we detected enrichment of candidate genes involved in the toxoplasmosis (TGFβ2/CHUK/CIITA/SOCS1) pathway in the tolerant animals. Interestingly, CHUK was one of the candidate genes in the toxoplasmosis pathway which highlights the role of CHUK as an important modulator of the inflammatory response. As part of the non-canonical pathway of NF-κβ activation, CHUK participates also in the negative feedback of NF-κβ signaling pathway preventing TNF-mediated RIPK1-dependent cell death. While the canonical NF-κβ activation is associated with rapid and acute production of diverse pro-inflammatory mediators, the non-canonical NF-κβ signaling pathway is tightly regulated and associated with less acute, chronic inflammation, and with maintaining homeostasis. Consequently, the role of CHUK in controlling inflammation by activating or limiting the activation of NF-κβ may be timing and/or environmental dependent. Overall, the activation of the toxoplasmosis pathway in the tolerant animals might result in a lifelong and persistent MAP infection as seen in Toxoplama gondii infected hosts. Toxoplasma gondii and other intracellular parasites which typically establish themselves as long-term residents in mammalian hosts are often associated with elevated host TGFβ expression (63). Direct activation of the signal transducer and activator of transcription 3 (STAT3) by the parasite enhances the anti-inflammatory function of TGFβ causing lifelong infection by establishing an anti-apoptotic environment. In addition, the role of TGFβ2 in tissue healing and repair is well established promoting disease tolerance at later time points (64). We also found two novel candidate genes associated with PTB tolerance, the Class II Major Histocompatibility Complex Transactivator (CIITA) and the suppressor of cytokine signaling (SOCS1). Both genes are located on BTA25 upstream of the RMI2/TNP2/PRIM1/PRM2/PRM3 region and adjacent to the CLEC16A gene, a regulator of phagosome-lysosome fusion during late mitophagy. CIITA is a positive regulator of class II major histocompatibility complex gene transcription and is the “master control factor” for the expression of these genes and T-cell activation. SOCS1 is a member of the STAT-induced STAT inhibitor (SSI) family, with important roles in immune cell homeostasis and regulation of inflammation which is crucial to avoid an excessive inflammatory response and to minimize tissue damage.
Conclusions
In this study, we explored the concept of disease tolerance in the context of a chronic intestinal infection. We were able to show that there is genetic variation associated with PTB tolerance. We revealed the genetic basis for PTB tolerance and the molecular pathways involved. Pathways analysis using the identified candidate genes revealed mechanisms underlying DNA packaging, innate immunity, and establishment of persistent infection that support tolerance to PTB. Strategies designed to improve diseases tolerance are expected to improve diseases control and to reduce the dependence on antimicrobials to kill pathogens, which in turn would reduce antimicrobial resistance. In livestock, the identification of genetic variants that regulate disease tolerance could help to breed for more tolerant cattle to a variety of diseases and likely to other environmental stressors as well.
Data Availability Statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.
Ethics Statement
Animals used in this study were not submitted to any in vivo experimentation and, therefore, no specific ethics committee authorization was required.
Author Contributions
MC performed the GWAS and wrote the manuscript. PV, JG, and RJ participated in samples collection and phenotypic characterization of the animals. GB-B performed the functional analysis, estimation of breeding values and cross-validation, and helped with the figures. OG-R and AF supervised the imputation to WGS, the GWAS, and the cross-validation of the genomic predictions. MA-H conceived and coordinated the project, supervised the data analysis, and wrote the manuscript. All authors revised and approved the final manuscript.
Funding
Financial support for this study was provided by a grant (RTI2018-094192-R-C21) funded by MCIN/AEI/10.13039/501100011033 and by FEDER, “Una manera de hacer Europa”. MC and GB-B have been awarded fellowships from INIA and MCIN/AEI/10.13039/501100011033 and “FSE Invierte en tu futuro”; grants FPI2016-00041 and PRE2019-090562, respectively.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Acknowledgments
This work has been possible thanks to the support of the computing infrastructure of the i2BASQUE Research and Academic Network. We gratefully acknowledge the Bull Genomes Consortium for providing accessibility to the 1,000 Bull genomes database.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fimmu.2022.820965/full#supplementary-material
References
1. Nielsen SS, Toft N. A Review of Prevalences of Paratuberculosis in Farmed Animals in Europe. Prev Vet Med (2009) 88:1–14. doi: 10.1016/j.prevetmed.2008.07.003
2. Hasonova L, Pavlik I. Economic Impact of Paratuberculosis in Dairy Cattle Herds: A Review. Vet Med (2006) 51:193–211. doi: 10.17221/5539-VETMED
3. Blanco Vázquez C, Alonso-Hearn M, Juste RA, Canive M, Iglesias T, Iglesias N, et al. Detection of Latent Forms of Mycobacterium Avium Subsp. Paratuberculosis Infection Using Host Biomarker-Based ELISAs Greatly Improves Paratuberculosis Diagnostic Sensitivity. PloS One (2020) 15:1–21. doi: 10.1371/journal.pone.0236336
4. Blanco Vázquez C, Balseiro A, Alonso-Hearn M, Juste RA, Iglesias N, Canive M, et al. Bovine Intelectin 2 Expression as a Biomarker of Paratuberculosis Disease Progression. Animals (2021) 11:1–13. doi: 10.3390/ani11051370
5. González J, Geijo MV, García-Pariente C, Verna A, Corpa JM, Reyes LE, et al. Histopathological Classification of Lesions Associated With Natural Paratuberculosis Infection in Cattle. J Comp Pathol (2005) 133:184–96. doi: 10.1016/j.jcpa.2005.04.007
6. Feller M, Huwiler K, Stephan R, Altpeter E, Shang A, Furrer H, et al. Mycobacterium Avium Subspecies Paratuberculosis and Crohn’s Disease: A Systematic Review and Meta-Analysis. Lancet Infect Dis (2007) 7:607–13. doi: 10.1016/S1473-3099(07)70211-6
7. Juste RA, Elguezabal N, Garrido JM, Pavon A, Geijo MV, Sevilla I, et al. On the Prevalence of M. Avium Subspecies Paratuberculosis DNA in the Blood of Healthy Individuals and Patients With Inflammatory Bowel Disease. PloS One (2008) 3:3–8. doi: 10.1371/journal.pone.0002537
8. Jeyanathan M, Boutros-Tadros O, Radhi J, Semret M, Bitton A, Behr MA. Visualization of Mycobacterium Avium in Crohn’s Tissue by Oil-Immersion Microscopy. Microbes Infect (2007) 9:1567–73. doi: 10.1016/j.micinf.2007.09.001
9. Bo M, Erre GL, Niegowska M, Piras M, Taras L, Longu MG, et al. Interferon Regulatory Factor 5 Is a Potential Target of Autoimmune Response Triggered by Epstein-Barr Virus and Mycobacterium Avium Subsp. Paratuberculosis in Rheumatoid Arthritis: Investigating a Mechanism of Molecular Mimicry. Clin Exp Rheumatol (2018) 36:376–81.
10. Yokoyama K, Cossu D, Hoshino Y, Tomizawa Y, Momotani E, Hattori N. Anti-Mycobacterial Antibodies in Paired Cerebrospinal Fluid and Serum Samples From Japanese Patients With Multiple Sclerosis or Neuromyelitis Optica Spectrum Disorder. J Clin Med (2018) 7:522. doi: 10.3390/jcm7120522
11. Mameli G, Cocco E, Frau J, Marrosu MG, Sechi LA. Epstein Barr Virus and Mycobacterium Avium Subsp. Paratuberculosis Peptides Are Recognized in Sera and Cerebrospinal Fluid of MS Patients. Sci Rep (2016) 6:1–10. doi: 10.1038/srep22401
12. Niegowska M, Rapini N, Piccinini S, Mameli G, Caggiu E, Manca Bitti ML, et al. Type 1 Diabetes at-Risk Children Highly Recognize Mycobacterium Avium Subspecies Paratuberculosis Epitopes Homologous to Human Znt8 and Proinsulin. Sci Rep (2016) 6:2–11. doi: 10.1038/srep22266
13. Mallikarjunappa S, Brito LF, Pant SD, Schenkel FS, Meade KG, Karrow NA. Johne’s Disease in Dairy Cattle: An Immunogenetic Perspective. Front Vet Sci (2021) 8:718987. doi: 10.3389/fvets.2021.718987
14. Kravitz A, Pelzer K, Nammalwar S. The Paratuberculosis Paradigm Examined: A Review of Host Genetic Resistance and Innate Immune Fitness in Mycobacterium Avium Subsp. Paratuberculosis Infection. Front Vet Sci (2021) 8:721706. doi: 10.3389/fvets.2021.721706
15. Canive M, González-Recio O, Fernández A, Vázquez P, Badia-Bringué G, Lavín JL, et al. (a) Identification of Loci Associated With Susceptibility to Mycobacterium Avium Subsp. Paratuberculosis Infection in Holstein Cattle Using Combinations of Diagnostic Tests and Imputed Whole-Genome Sequence Data. PloS One (2021) 16:1–17. doi: 10.1371/journal.pone.0256091
16. Canive M, Badia-Bringué G, Vázquez P, González-Recio O, Fernández A, Garrido JM, et al. (B) Identification of Loci Associated With Pathological Outcomes in Holstein Cattle Infected With Mycobacterium Avium Subsp. Paratuberculosis Using Whole-Genome Sequence Data. Sci Rep (2021) 11:1–13. doi: 10.1038/s41598-021-99672-4
17. Divangahi M, Aaby P, Khader SA, Barreiro LB, Bekkering S, Chavakis T, et al. Trained Immunity, Tolerance, Priming and Differentiation: Distinct Immunological Processes. Nat Immunol (2021) 22:2–6. doi: 10.1038/s41590-021-00960-y
18. Medzhitov R, Schneider DS, Soares MP. Disease Tolerance as a Defense Strategy. Science (2012) 335(6071):936–41. doi: 10.1126/science.1214935
19. Caldwell RM, Schafer JF, Compton LE, Patterson FL. Tolerance to Cereal Leaf Rusts. Science (1958) 128(3326):714–5. doi: 10.1126/science.128.3326.714
20. Schafer JF. Tolerance to Plant Disease. Annu Rev Phytopathol (1971) 9:235–52. doi: 10.1146/annurev.py.09.090171.001315
21. Raberg L, Sim D, Read AF. Disentangling Genetic Variation for Resistance and Tolerance to Infectious Diseases in Animals. Science (2007) 318(5851):812–4. doi: 10.1126/science.1148526
22. Read AF, Graham AL, Råberg L. Animal Defenses Against Infectious Agents: Is Damage Control More Important Than Pathogen Control? PloS Biol (2008) 6:2638–41. doi: 10.1371/journal.pbio.1000004
23. Barreiro LB, Quintana-Murci L. From Evolutionary Genetics to Human Immunology: How Selection Shapes Host Defence Genes. Nat Rev Genet (2010) 11:17–30. doi: 10.1038/nrg2698
24. Doeschl-Wilson AB, Kyriazakis I. Should We Aim for Genetic Improvement in Host Resistance or Tolerance to Infectious Pathogens? Front Genet (2012) 14:272. doi: 10.3389/fgene.2012.00272
25. Vázquez P, Ruiz-Larrañaga O, Garrido JM, Iriondo M, Manzano C, Agirre M, et al. Genetic Association Analysis of Paratuberculosis Forms in Holstein-Friesian Cattle. Vet Med Int (2014) 2014:321327. doi: 10.1155/2014/321327
26. Larsen AB, Merkal RS, Cutlip RC. Age of Cattle as Related to Resistance to Infection With Mycobacterium Paratuberculosis. Am J Vet Res (1975) 36:255–7.
27. Sigurethardóttir OG, Valheim M, Press CM. Establishment of Mycobacterium Avium Subsp. Paratuberculosis Infection in the Intestine of Ruminants. Adv Drug Deliv Rev (2004) 56:819–34. doi: 10.1016/j.addr.2003.10.032
28. Kurade NP, Tripathi BN, Rajukumar K, Parihar NS. Sequential Development of Histologic Lesions and Their Relationship With Bacterial Isolation, Fecal Shedding, and Immune Responses During Progressive Stages of Experimental Infection of Lambs With Mycobacterium Avium Subsp. Paratuberculosis. Vet Pathol (2004) 41:378–87. doi: 10.1354/vp.41-4-378
29. Loh PR, Palamara PF, Price AL. Fast and Accurate Long-Range Phasing in a UK Biobank Cohort. Nat Genet (2016) 48:811–6. doi: 10.1038/ng.3571
30. Das S, Forer L, Schönherr S, Sidore C, Locke AE, Kwong A, et al. Next-Generation Genotype Imputation Service and Methods. Nat Genet (2016) 48:1284–7. doi: 10.1038/ng.3656
31. Hayes BJ, Daetwyler HD. 1000 Bull Genome Project to Map Simple and Complex Genetic Traits in Cattle: Applications and Outcomes. Annu Rev Anim. Biosci (2019) 7:89–102. doi: 10.1146/annurev-animal-020518-115024
32. Yang J, Lee SH, Goddard ME, Visscher PM. GCTA: A Tool for Genome-Wide Complex Trait Analysis. Am J Hum Genet (2011) 88:76–82. doi: 10.1016/j.ajhg.2010.11.011
33. The Welcome Trust Case Control Consortium. Genome-Wide Association Study of 14,000 Cases of Seven Common Disease and 3,000 Shared Controls. Nature (2007) 447:661–77. doi: 10.1038/nature05911
35. González JR, Armengol L, Solé X, Guinó E, Mercader JP, Estivill X, et al. SNPassoc: An R Package to Perform Whole Genome Association Studies. Bioinformatics (2007) 23:644–5. doi: 10.1093/bioinformatics/btm025
36. Sargolzaei M, Schenkel FS, Jansen GB, Schaeffer LR. Extent of Linkage Disequilibrium in Holstein Cattle in North America. J Dairy Sci (2008) 91:2106–17. doi: 10.3168/jds.2007-0553
37. Kim E-S, Kirkpatrick BW. Linkage Disequilibrium in the North American Holstein Population. Anim Genet (2009) 40:279–88. doi: 10.1111/j.1365-2052.2008.01831.x
38. Yu G, Wang LG, Han Y, He QY. ClusterProfiler: An R Package for Comparing Biological Themes Among Gene Clusters. OMICS (2012) 16:284–77. doi: 10.1089/omi.2011.0118
39. VanRaden PM, Sullivan PG. International Genomic Evaluation Methods for Dairy Cattle. Genet Sel. Evol (2010) 42:1–9. doi: 10.1186/1297-9686-42-7
40. Sing T, Sander O, Beerenwinkel N, Lengauer T. ROCR: Visualizing Classifier Performance in R. Bioinformatics (2005) 21:3940–1. doi: 10.1093/bioinformatics/bti623
41. Kirkwood BR, Sterne JAC. Essential Medical Statistics. 2nd. Oxford, UK: Blackwell Publishing (2010).
42. Kirkpatrick BW, Shi X, Shook GE, Collins MT. Whole-Genome Association Analysis of Susceptibility to Paratuberculosis in Holstein Cattle. Anim Genet (2011) 42:149–60. doi: 10.1111/j.1365-2052.2010.02097.x
43. Holmberg M, Andersson-Eklund L. Quantitative Trait Loci Affecting Health Traits in Swedish Dairy Cattle. J Dairy Sci (2004) 87:2653–9. doi: 10.3168/jds.S0022-0302(04)73391-3
44. Lund MS, Guldbrandtsen B, Buitenhuis AJ, Thomsen B, Bendixen C. Detection of Quantitative Trait Loci in Danish Holstein Cattle Affecting Clinical Mastitis, Somatic Cell Score, Udder Conformation Traits, and Assessment of Associated Effects on Milk Yield. J Dairy Sci (2008) 91:4028–36. doi: 10.3168/jds.2007-0290
45. Cole JB, Wiggans GR, Ma L, Sonstegard TS, Lawlor TJ Jr, Crooker BA, et al. Genome-Wide Association Analysis of Thirty One Production, Health, Reproduction and Body Conformation Traits in Contemporary U.S. Holstein Cows. BMC Genomics (2011) 12:1–17. doi: 10.1186/1471-2164-12-408
46. Richardson IW, Berry DP, Wiencko HL, Higgins IM, More SJ, McClure J, et al. A Genome-Wide Association Study for Genetic Susceptibility to Mycobacterium Bovis Infection in Dairy Cattle Identifies a Susceptibility QTL on Chromosome 23. Genet Sel Evol (2016) 48:1–13. doi: 10.1186/s12711-016-0197-x
47. Strillacci MG, Frigo E, Schiavini F, Samoré AB, Canavesi F, Vevey M, et al. Genome-Wide Association Study for Somatic Cell Score in Valdostana Red Pied Cattle Breed Using Pooled DNA. BMC Genet (2014) 15:106. doi: 10.1186/s12863-014-0106-7
48. Oliveira HR, Cant JP, Brito LF, Feitosa FLB, Chud TCS, Fonseca PAS, et al. Genome-Wide Association for Milk Production Traits and Somatic Cell Score in Different Lactation Stages of Ayrshire, Holstein, and Jersey Dairy Cattle. J Dairy Sci (2019) 102:8159–74. doi: 10.3168/jds.2019-16451
49. Neibergs H, Zanella R, Casas E, Snowder GD, Wenz J, Neibergs JS, et al. Loci on Bos Taurus Chromosome 2 and Bos Taurus Chromosome 26 Are Linked With Bovine Respiratory Disease and Associated With Persistent Infection of Bovine Viral Diarrhea Virus. J Anim Sci (2011) 89:907–15. doi: 10.2527/jas.2010-3330
50. Zanella R, Whitlock RH, Neibergs HL. Short Communication: Refinement of Genetic Regions Associated With Mycobacterium Avium Subspecies Paratuberculosis Tissue Infection and Tolerance to Johne’s Disease. J Dairy Sci (2011) 94:4230–6. doi: 10.3168/jds.2011-4192
51. Buergelt CD, Hall C, McEntee K, Duncan JR. Pathological Evaluation of Paratuberculosis in Naturally Infected Cattle. Vet Pathol (1978) 15:196–207. doi: 10.1177/030098587801500206
52. Huda A, Jensen HE. Comparison of Histopathology, Cultivation of Tissues and Rectal Contents, and Interferon-Gamma and Serum Antibody Responses for the Diagnosis of Bovine Paratuberculosis. J Comp Pathol (2003) 129:259–67. doi: 10.1016/S0021-9975(03)00042-2
53. Xu D, Guo R, Sobeck A, Bachrati CZ, Yang J, Enomoto T, et al. RMI, a New OB-Fold Complex Essential for Bloom Syndrome Protein to Maintain Genome Stability. Genes Dev (2008) 22:2843–55. doi: 10.1101/gad.1708608
54. Soares MP, Teixeira L, Moita LF. Disease Tolerance and Immunity in Host Protection Against Infection. Nat Rev Immunol (2017) 7(2):83–96. doi: 10.1038/nri.2016.136
55. Balhorn R. The Protamine Family of Sperm Nuclear Proteins. Genome Biol (2007) 8:227. doi: 10.1186/gb-2007-8-9-227
56. Rousseaux S, Reynoird N, Escoffier E, Thevenon J, Caron C, Khochbin S. Epigenetic Reprogramming of the Male Genome During Gametogenesis and in the Zygote. Reprod BioMed (2008) 16:492–503. doi: 10.1016/S1472-6483(10)60456-7
57. Simon L, Castillo J, Oliva R, Lewis SEM. Relationships Between Human Sperm Protamines, DNA Damage and Assisted Reproduction Outcomes. Reprod BioMed (2011) 23:724–34. doi: 10.1016/j.rbmo.2011.08.010
58. Björndahl L, Kvist U. Structure of Chromatin in Spermatozoa. Adv Exp Med Biol (2014) 791:1–11. doi: 10.1007/978-1-4614-7783-9_1
59. Miller D, Brinkworth M, Iles D. Paternal DNA Packaging in Spermatozoa: More Than the Sum of Its Parts? DNA, Histones, Protamines and Epigenetics. Reproduction (2010) 139:287–301. doi: 10.1530/REP-09-0281
60. Sheldon IM, Cronin JG, Bromfield JJ. Tolerance and Innate Immunity Shape the Development of Postpartum Uterine Disease and the Impact of Endometritis in Dairy Cattle. Annu Rev Anim Biosci (2019) 7:361–84. doi: 10.1146/annurev-animal-020518-115227
61. Schneider DS, Ayres JS, Brandt SM, Costa A, Dionne MS, Gordon MD, et al. Drosophila Eiger Mutants Are Sensitive to Extracellular Pathogens. PlosPathogens (2007) 3:e41. doi: 10.1371/journal.ppat.0030041
62. Roca FJ, Ramakrishnan L. TNF Dually Mediates Resistance and Susceptibility to Mycobacteria via Mitochondrial Reactive Oxygen Species. Cell (2013) 153:521–34. doi: 10.1016/j.cell.2013.03.022
63. Johnson PTJ, Wood CL, Joseph MB, Preston DL, Haas SE, Springer YP, et al. Habitat Heterogeneity Drives the Host-Diversity-Begets-Parasite-Diversity Relationship: Evidence From Experimental and Field Studies. Ecol Lett (2016) 19:752–61. doi: 10.1111/ele.12609
Keywords: disease tolerance, paratuberculosis, DNA packaging, DNA damage repair, innate immunity, tissue damage control, pathogen persistence
Citation: Canive M, Badia-Bringué G, Vázquez P, Garrido JM, Juste RA, Fernandez A, González-Recio O and Alonso-Hearn M (2022) A Genome-Wide Association Study for Tolerance to Paratuberculosis Identifies Candidate Genes Involved in DNA Packaging, DNA Damage Repair, Innate Immunity, and Pathogen Persistence. Front. Immunol. 13:820965. doi: 10.3389/fimmu.2022.820965
Received: 15 December 2021; Accepted: 14 March 2022;
Published: 06 April 2022.
Edited by:
Philip John Griebel, University of Saskatchewan, CanadaReviewed by:
Holly Louise Neibergs, Washington State University, United StatesPetr Horin, University of Veterinary and Pharmaceutical Sciences Brno, Czechia
Copyright © 2022 Canive, Badia-Bringué, Vázquez, Garrido, Juste, Fernandez, González-Recio and Alonso-Hearn. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Marta Alonso-Hearn, malonso@neiker.eus