Identification of loci associated with susceptibility to subspecies () tissue infection in cattle

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.

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.

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.

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

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

Johne's Disease in Dairy Cattle: An Immunogenetic Perspective

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

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

Bovine paratuberculosis (PTB) is a chronic inflammatory disease caused by Mycobacterium avium susbp. paratuberculosis (MAP). Genome-wide association studies (GWAS) have identified single nucleotide polymorphisms (SNPs) significantly associated with susceptibility to bovine PTB. The main objective of this study was to identify quantitative trait loci (QTLs) associated with MAP infection in Spanish Holstein cows (N = 983) using combinations of diagnostic tests and imputed whole-genome sequence (WGS) data. The infection status of these animals was defined by three diagnostic methods including ELISA for MAP-antibodies detection, and tissue culture and PCR for MAP detection. The 983 cows included in this study were genotyped with the Bovine MD SNP50 Bead Chip, and the corresponding genotypes were imputed to WGS using the 1,000 Bull genomes reference population. In total, 33.77 million SNP variants per animal were identified across the genome. Linear mixed models were used to calculate the heritability (h²) estimates for each diagnostic test and test combinations. Next, we performed a case-control GWAS using the imputed WGS datasets and the phenotypes and combinations of phenotypes with h² estimates > 0.080. After performing the GWAS, the test combinations that showed SNPs with a significant association (P_FDR ≤ 0.05), were the ELISA-tissue PCR-tissue culture, ELISA-tissue culture, and ELISA-tissue PCR. A total of twelve quantitative trait loci (QTLs) highly associated with MAP infection status were identified on the Bos taurus autosomes (BTA) 4, BTA5, BTA11, BTA12, BTA14, BTA23, BTA24, and BTA28, and some of these QTLs were linked to immune-modulating genes. The identified QTLs on BTA23 spanning from 18.81 to 22.95 Mb of the Bos taurus genome overlapped with several QTLs previously found to be associated with PTB susceptibility, bovine tuberculosis susceptibility, and clinical mastitis. The results from this study provide more clues regarding the molecular mechanisms underlying susceptibility to PTB infection in cattle and might be used to develop national genetic evaluations for PTB in Spain.

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.

MAP, Johne's disease and the microbiome; current knowledge and future considerations

Mycobacterium avium subsp. paratuberculosis is the causative agent of Johne's disease in ruminants. As an infectious disease that causes reduced milk yields, effects fertility and, eventually, the loss of the animal, it is a huge financial burden for associated industries. Efforts to control MAP infection and Johne's disease are complicated due to difficulties of diagnosis in the early stages of infection and challenges relating to the specificity and sensitivity of current testing methods. The methods that are available contribute to widely used test and cull strategies, vaccination programmes also in place in some countries. Next generation sequencing technologies have opened up new avenues for the discovery of novel biomarkers for disease prediction within MAP genomes and within ruminant microbiomes. Controlling Johne's disease in herds can lead to improved animal health and welfare, in turn leading to increased productivity. With current climate change bills, such as the European Green Deal, targeting livestock production systems for more sustainable practices, managing animal health is now more important than ever before. This review provides an overview of the current knowledge on genomics and detection of MAP as it pertains to Johne's disease.

Functional Variants Surrounding Endothelin 2 Are Associated With Mycobacterium avium Subspecies paratuberculosis Infection

Bovine paratuberculosis, caused by Mycobacterium avium subspecies paratuberculosis (MAP), continues to impact the dairy industry through increased morbidity, mortality, and lost production. Although genome-wide association analyses (GWAAs) have identified loci associated with susceptibility to MAP, limited progress has been made in identifying mutations that cause disease susceptibility. A 235-kb region on Bos taurus chromosome 3 (BTA3), containing a 70-kb haplotype block surrounding endothelin 2 (EDN2), has previously been associated with the risk of MAP infection. EDN2 is highly expressed in the gut and is involved in intracellular calcium signaling and a wide array of biological processes. The objective of this study was to identify putative causal mutations for disease susceptibility in the region surrounding EDN2 in Holstein and Jersey cattle. Using sequence data from 10 Holstein and 10 Jersey cattle, common variants within the 70-kb region containing EDN2 were identified. A custom SNP genotyping array fine-mapped the region using 221 Holstein and 51 Jersey cattle and identified 17 putative causal variants (P < 0.01) located in the 5′ region of EDN2 and a SNP in the 3′ UTR (P = 0.00009) associated with MAP infection. MicroRNA interference assays, mRNA stability assays, and electrophoretic mobility shift assays were performed to determine if allelic changes at each SNP resulted in differences in EDN2 stability or expression. Two SNPs [rs109651404 (G/A) and rs110287192 (G/T)] located within the promoter region of EDN2 displayed differential binding affinity for transcription factors in binding sequences harboring the alternate SNP alleles. The luciferase reporter assay revealed that the transcriptional activity of the EDN2 promoter was increased (P < 0.05) with the A allele for rs109651404 and the G allele for rs110287192. These results suggest that the variants rs109651404 and rs110287192 are mutations that alter transcription and thus may alter susceptibility to MAP infection in Holstein and Jersey cattle.

Genome-wide association analysis identified both RNA-seq and DNA variants associated to paratuberculosis in Canadian Holstein cattle 'in vitro' experimentally infected macrophages

Background

Mycobacterium avium ssp. paratuberculosis (MAP) is the causative agent of paratuberculosis, or Johne’s disease (JD), an incurable bovine disease. The evidence for susceptibility to MAP disease points to multiple interacting factors, including the genetic predisposition to a dysregulation of the immune system. The endemic situation in cattle populations can be in part explained by a genetic susceptibility to MAP infection. In order to identify the best genetic improvement strategy that will lead to a significant reduction of JD in the population, we need to understand the link between genetic variability and the biological systems that MAP targets in its assault to dominate macrophages. MAP survives in macrophages where it disseminates. We used next-generation RNA (RNA-Seq) sequencing to study of the transcriptome in response to MAP infection of the macrophages from cows that have been naturally infected and identified as positive for JD (JD (+); n = 22) or negative for JD (healthy/resistant, JD (−); n = 28). In addition to identifying genetic variants from RNA-seq data, SNP variants were also identified using the Bovine SNP50 DNA chip.

Results

The complementary strategy allowed the identification of 1,356,248 genetic variants, including 814,168 RNA-seq and 591,220 DNA chip variants. Annotation using SnpEff predicted that the 2435 RNA-seq genetic variants would produce high functional effect on known genes in comparison to the 33 DNA chip variants. Significant variants from JD(+/−) macrophages were identified by genome-wide association study and revealed two quantitative traits loci: BTA4 and 11 at (P < 5 × 10^− 7). Using BovineMine, gene expression levels together with significant genomic variants revealed pathways that potentially influence JD susceptibility, notably the energy-dependent regulation of mTOR by LKB1-AMPK and the metabolism of lipids.

Conclusion

In the present study, we succeeded in identifying genetic variants in regulatory pathways of the macrophages that may affect the susceptibility of cows that are healthy/resistant to MAP infection. RNA-seq provides an unprecedented opportunity to investigate gene expression and to link the genetic variations to biological pathways that MAP normally manipulate during the process of killing macrophages. A strategy incorporating functional markers into genetic selection may have a considerable impact in improving resistance to an incurable disease. Integrating the findings of this research into the conventional genetic selection program may allow faster and more lasting improvement in resistance to bovine paratuberculosis in dairy cattle.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07487-4.

Loci Associated With Antibody Response in Feral Swine (Sus scrofa) Infected With Brucella suis

Feral swine (Sus scrofa) are a destructive invasive species widespread throughout the United States that disrupt ecosystems, damage crops, and carry pathogens of concern for the health of domestic stock and humans including Brucella suis-the causative organism for swine brucellosis. In domestic swine, brucellosis results in reproductive failure due to abortions and infertility. Contact with infected feral swine poses spillover risks to domestic pigs as well as humans, companion animals, wildlife, and other livestock. Genetic factors influence the outcome of infectious diseases; therefore, genome wide association studies (GWAS) of differential immune responses among feral swine can provide an understanding of disease dynamics and inform management to prevent the spillover of brucellosis from feral swine to domestic pigs. We sought to identify loci associated with differential antibody responses among feral swine naturally infected with B. suis using a case-control GWAS. Tissue, serum, and genotype data (68,516 bi-allelic single nucleotide polymorphisms) collected from 47 feral swine were analyzed in this study. The 47 feral swine were culture positive for Brucella spp. Of these 47, 16 were antibody positive (cases) whereas 31 were antibody negative (controls). Single-locus GWAS were performed using efficient mixed-model association eXpedited (EMMAX) methodology with three genetic models: additive, dominant, and recessive. Eight loci associated with seroconversion were identified on chromosome 4, 8, 9, 10, 12, and 18. Subsequent bioinformatic analyses revealed nine putative candidate genes related to immune function, most notably phagocytosis and induction of an inflammatory response. Identified loci and putative candidate genes may play an important role in host immune responses to B. suis infection, characterized by a detectable bacterial presence yet a differential antibody response. Given that antibody tests are used to evaluate brucellosis infection in domestic pigs and for disease surveillance in invasive feral swine, additional studies are needed to fully understand the genetic component of the response to B. suis infection and to more effectively translate estimates of Brucella spp. antibody prevalence among feral swine to disease control management action.

Association of genetic polymorphisms related to Johne's disease with estimated breeding values of Holstein sires for milk ELISA test scores

BACKGROUND

Johne's disease (JD) is a chronic intestinal inflammatory disease caused by Mycobacterium avium subsp. paratuberculosis (MAP) infection in ruminants. Since there are currently no effective vaccine or treatment options available to control JD, genetic selection may be an alternative strategy to enhance JD resistance. Numerous Single Nucleotide Polymorphisms (SNPs) have been reported to be associated with MAP infection status based on published genome-wide association and candidate gene studies. The main objective of this study was to validate these SNPs that were previously identified to be associated with JD by testing their effect on Holstein bulls' estimated breeding values (EBVs) for milk ELISA test scores, an indirect indicator of MAP infection status in cattle.

RESULTS

Three SNPs, rs41810662, rs41617133 and rs110225854, located on Bos taurus autosomes (BTA) 16, 23 and 26, respectively, were confirmed as significantly associated with Holstein bulls' EBVs for milk ELISA test score (FDR < 0.01) based on General Quasi Likelihood Scoring analysis (GQLS) analysis. Single-SNP regression analysis identified four SNPs that were associated with sire EBVs (FDR < 0.05). This includes two SNPs that were common with GQLS (rs41810662 and rs41617133), with the other two SNPs being rs110494981 and rs136182707, located on BTA9 and BTA16, respectively.

CONCLUSIONS

The findings of this study validate the association of SNPs with JD MAP infection status and highlight the need to further investigate the genomic regions harboring these SNPs.

Identification of the ABCC4, IER3, and CBFA2T2 candidate genes for resistance to paratuberculosis from sequence-based GWAS in Holstein and Normande dairy cattle

Background

Bovine paratuberculosis is a contagious disease, caused by Mycobacterium avium subsp. paratuberculosis (MAP), with adverse effects on animal welfare and serious economic consequences. Published results on host genetic resistance to MAP are inconsistent, mainly because of difficulties in characterizing the infection status of cows. The objectives of this study were to identify quantitative trait loci (QTL) for resistance to MAP in Holstein and Normande cows with an accurately defined status for MAP.

Results

From MAP-infected herds, cows without clinical signs of disease were subjected to at least four repeated serum ELISA and fecal PCR tests over time to determine both infected and non-infected statuses. Clinical cases were confirmed using PCR. Only cows that had concordant results for all tests were included in further analyses. Positive and control cows were matched within herd according to their birth date to ensure a same level of exposure to MAP. Cows with accurate phenotypes, i.e. unaffected (control) or affected (clinical or non-clinical cases), were genotyped with the Illumina BovineSNP50 BeadChip. Genotypes were imputed to whole-genome sequences using the 1000 Bull Genomes reference population (run6). A genome-wide association study (GWAS) of MAP status of 1644 Holstein and 649 Normande cows, using either two (controls versus cases) or three classes of phenotype (controls, non-clinical and clinical cases), revealed three regions, on Bos taurus (BTA) chromosomes 12, 13, and 23, presenting significant effects in Holstein cows, while only one of those was identified in Normande cows (BTA23). The most significant effect was found on BTA13, in a short 8.5-kb region. Conditional analyses revealed that only one causal variant may be responsible for the effects observed on each chromosome with the ABCC4 (BTA12), CBFA2T2 (BTA13), and IER3 (BTA23) genes as good functional candidates.

Conclusions

A sequence-based GWAS on cows for which resistance to MAP was accurately defined, was able to identify candidate variants located in genes that were functionally related to resistance to MAP; these explained up to 28% of the genetic variance of the trait. These results are very encouraging for efforts towards implementation of a breeding strategy aimed at improving resistance to paratuberculosis in Holstein cows.

Genome-wide pleiotropy and shared biological pathways for resistance to bovine pathogens

Host genetic architecture is a major factor in resistance to pathogens and parasites. The collection and analysis of sufficient data on both disease resistance and host genetics has, however, been a major obstacle to dissection the genetics of resistance to single or multiple pathogens. A severe challenge in the estimation of heritabilities and genetic correlations from pedigree-based studies has been the confounding effects of the common environment shared among relatives which are difficult to model in pedigree analyses, especially for health traits with low incidence rates. To circumvent this problem we used genome-wide single-nucleotide polymorphism data and implemented the Genomic-Restricted Maximum Likelihood (G-REML) method to estimate the heritabilities and genetic correlations for resistance to 23 different infectious pathogens in calves and cows in populations undergoing natural pathogen challenge. Furthermore, we conducted gene-based analysis and generalized gene-set analysis to understand the biological background of resistance to infectious diseases. The results showed relatively higher heritabilities of resistance in calves than in cows and significant pleiotropy (both positive and negative) among some calf and cow resistance traits. We also found significant pleiotropy between resistance and performance in both calves and cows. Finally, we confirmed the role of the B-lymphocyte pathway as one of the most important biological pathways associated with resistance to all pathogens. These results both illustrate the potential power of these approaches to illuminate the genetics of pathogen resistance in cattle and provide foundational information for future genomic selection aimed at improving the overall production fitness of cattle.

Identification of genes associated with susceptibility to Mycobacterium avium ssp. paratuberculosis (Map) tissue infection in Holstein cattle using gene set enrichment analysis-SNP

Multiple genome-wide association analyses have investigated susceptibility to bovine paratuberculosis, but few loci have been identified across independent cattle populations. A SNP-based gene set enrichment analysis (GSEA-SNP) allows expanded identification of genes with moderate effects on a trait through the enrichment of gene sets instead of identifying only few loci with large effects. Therefore, the objective of this study was to identify genes that were moderately associated with Mycobacterium avium ssp. paratuberculosis (Map) tissue infection using GSEA-SNP in Holstein cattle from the Pacific Northwest (PNW; n = 205) and from the PNW and Northeast (PNW+NE; n = 245) which were previously genotyped with the Illumina BovineSNP50 BeadChip. The GSEA-SNP utilized 4389 gene sets from five databases. For each annotated gene in the UMD3.1 assembly (n = 19,723), the most significant SNP within each gene and its surrounding region (10 kb up- and downstream) was selected as a proxy for that gene. Any gene set with a normalized enrichment score > 2.5 was considered enriched. Thirteen gene sets (8 PNW GSEA-SNP; 5 PNW+NE) were enriched in these analyses and all have functions that relate to nuclear factor kappa beta. Nuclear factor kappa beta is critical to gut immune responses, implicated in host immune responses to other mycobacterial diseases, and has established roles in inflammation as well as cancer. Gene sets and genes moderately associated with Map infection could be used in genomic selection to allow producers to select for less susceptible cattle, lower the prevalence of the disease, and reduce economic losses.