GSEA-SNP identifies genes associated with Johne's disease in cattle

Summary-data based Mendelian randomization identifies gene expression regulatory polymorphisms associated with bovine paratuberculosis by modulation of the nuclear factor Kappa β (NF-κß)-mediated inflammatory response

Genome-wide association studies (GWAS) have identified host genetic variants associated with paratuberculosis (PTB) susceptibility. Most of the GWAS-identified SNPs are in non-coding regions. Connecting these non-coding variants and downstream affected genes is a challenge and, up to date, only a few functional mutations or expression quantitative loci (cis-eQTLs) associated with PTB susceptibility have been identified. In the current study, the associations between imputed whole-genome sequence genotypes and whole RNA-Sequencing data from peripheral blood (PB) and ileocecal valve (ICV) samples of Spanish Holstein cows (N = 16) were analyzed with TensorQTL. This approach allowed the identification of 88 and 37 cis-eQTLs regulating the expression levels of 90 and 37 genes in PB and ICV samples, respectively (False discorey rate, FDR ≤ 0.05). Next, we applied summary-based data Mendelian randomization (SMR) to integrate the cis-eQTL dataset with GWAS data obtained from a cohort of 813 culled cattle that were classified according to the presence or absence of PTB-associated histopathological lesions in gut tissues. After multiple testing corrections (FDR ≤ 0.05), we identified two novel cis-eQTLs affecting the expression of the early growth response factor 4 (EGR4) and the bovine neuroblastoma breakpoint family member 6-like protein isoform 2 (MGC134040) that showed pleiotropic associations with the presence of multifocal and diffuse lesions in gut tissues; P = 0.002 and P = 0.017, respectively. While EGR4 acts as a brake on T-cell proliferation and cytokine production through interaction with the nuclear factor Kappa β (NF-κß), MGC134040 is a target gene of NF-κß. Our findings provide a better understanding of the genetic factors influencing PTB outcomes, confirm that the multifocal lesions are localized/confined lesions that have different underlying host genetics than the diffuse lesions, and highlight regulatory SNPs and regulated-gene targets to design future functional studies.

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.

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.

Selection signatures of Fuzhong Buffalo based on whole-genome sequences

BACKGROUND

Fuzhong buffalo, a native breed of Guangxi Zhuang Autonomous Region, is traditionally used as a draft animal to provide farm power in the rice cultivation. In addition, the Fuzhong buffalo also prepared for the bullfighting festival organized by the locals. The detection of the selective signatures in its genome can help in elucidating the selection mechanisms in its stamina and muscle development of a draft animal.

RESULTS

In this study, we analyzed 27 whole genomes of buffalo (including 15 Fuzhong buffalo genomes and 12 published buffalo genomes from Upper Yangtze region). The ZHp, ZFst, π-Ratio, and XP-EHH statistics were used to identify the candidate signatures of positive selection in Fuzhong buffalo. Our results detected a set of candidate genes involving in the pathways and GO terms associated with the response to exercise (e.g., ALDOA, STAT3, AKT2, EIF4E2, CACNA2D2, TCF4, CDH2), immunity (e.g., PTPN22, NKX2-3, PIK3R1, ITK, TMEM173), nervous system (e.g., PTPN21, ROBO1, HOMER1, MAGI2, SLC1A3, NRG3, SNAP47, CTNNA2, ADGRL3). In addition, we also identified several genes related to production and growth traits (e.g., PHLPP1, PRKN, MACF1, UCN3, RALGAPA1, PHKB, PKD1L). Our results depicted several pathways, GO terms, and candidate genes to be associated with response to exercise, immunity, nervous system, and growth traits.

CONCLUSIONS

The selective sweep analysis of the Fuzhong buffalo demonstrated positive selection pressure on potential target genes involved in behavior, immunity, and growth traits, etc. Our findings provided a valuable resource for future research on buffalo breeding and an insight into the mechanisms of artificial selection.

The EDN2 rs110287192 gene polymorphism is associated with paratuberculosis susceptibility in multibreed cattle population

Paratuberculosis (pTB), also known as Johne's disease (JD), is a contagious, chronic, and granulomatous inflammatory disease of the intestines of ruminants which is caused by Mycobacterium avium subsp. paratuberculosis (MAP) infection, resulting in billions of dollars in economic losses worldwide. Since, currently, no effective cure is available for MAP infection, it is important to explore the genetic variants that affect the host MAP susceptibility. The aim of this study was to analyze a potential association between EDN2 synonymous gene mutations (rs110287192, rs109651404 and rs136707411), that modifies susceptibility to pTB. EDN2 rs110287192, rs109651404 and rs136707411 mutations were genotyped in 68 infected and 753 healthy animals from East Anatolian Red crossbred, Anatolian Black crossbred and Holstein breed cattle by using Custom TaqMan SNP Genotyping Assays. For pTB status, serum antibody levels S/P ≥ 1.0 were assessed in carriers of the different EDN2 genotypes. EDN2 rs110287192 mutation showed a significant association with bovine pTB (adj. p < 0.05). For rs110287192 locus, the odd ratios for GG and TG genotypes versus TT genotypes were 1.73; (95% CI = 0.34–8.59) and 0.53 (95% CI = 0.12–2.37) respectively, which indicated that proportion of TG heterozygotes were significantly higher in control animals as compared to pTB animals. On the other hand, while rs136707411 mutation showed a suggestive association with pTB status in the examined cattle population (nominal p < 0.05); no association was detected between rs109651404 genotypes and pTB status. Selecting animals against rs110287192-GG genotype may decrease the risk of pTB in cattle of the Bos taurus taurus subspecies.

A missense mutation in the coding region of the toll-like receptor 4 gene affects milk traits in Barki sheep

OBJECTIVE

Milk production is one of the most desirable traits in livestock. Recently, the toll-like receptor (TLR) has been identified as a candidate gene for milk traits in cows. So far, there is no information concerning the contribution of this gene in milk traits in sheep. This study was designed to investigate the TLR 4 gene polymorphisms in Barki ewes in Egypt and then correlate that with milk traits in order to identify potential single nucleotide polymorphisms (SNPs) for these traits in sheep.

METHODS

A part of the ovine TLR 4 gene was amplified in Barki ewes, to identify the SNPs. Consequently; Barki ewes were genotyped using polymerase chain reaction-single strand conformation polymorphism protocol. These genotypes were correlated with milk traits, which were the daily milk yield (DMY), protein percentage (PP), fat percentage (FP), lactose percentage, and total solid percentage (TSP).

RESULTS

Age and parity of the ewe had a significant effect (p<0.05 or p<0.01) on DMY, FP, and TSP. The direct sequencing identified a missense mutation located in the coding sequence of the gene (rs592076818; c.1710C>A) and was predicted to change the amino acid sequence of the resulted protein (p.Asn570Lys). The association analyses suggested a significant effect (p<0.05) of the TLR genotype on the FP and PP, while the DMY tended to be influenced as well (p = 0.07). Interestingly, the presence of the G allele tended to increase the DMY (+40.5 g/d) and significantly (p<0.05 or p<0.01) decreased the FP (-1.11%), PP (-1.21%), and TSP (-7.98%).

CONCLUSION

The results of this study suggested the toll-like receptor 4 (TLR4) as a candidate gene to improve milk traits in sheep worldwide, which will enhance the ability to understand the genetic architecture of genes underlying SNPs that affect such traits.

Candidate genes associated with the heritable humoral response to Mycobacterium avium ssp. paratuberculosis in dairy cows have factors in common with gastrointestinal diseases in humans

Infection of cattle with bovine paratuberculosis (i.e., Johne's disease) is caused by Mycobacterium avium ssp. paratuberculosis (MAP) and results in a chronic incurable gastroenteritis. This disease, which has economic ramifications for the cattle industry, is increasing in detected prevalence globally; subclinically infected animals can silently shed the bacterium into the environment for years, exposing contemporaries and hampering disease-control programs. The objective of the present study was to first quantify the genetic parameters for humoral response to MAP in dairy cattle. This was followed by a genome-based association analysis and subsequent downstream bioinformatic analyses from imputed whole genome sequence SNP data. After edits, ELISA test records were available on 136,767 cows; analyses were also undertaken on a subset of 33,818 of these animals from herds with at least 5 MAP ELISA-positive cows, with at least 1 of those positive cows being homebred. Variance components were estimated using univariate animal and sire linear mixed models. The heritability calculated from the animal model for humoral response to MAP using alternative phenotype definitions varied from 0.02 (standard error = 0.003) to 0.05 (standard error = 0.008). The genome-based associations were undertaken within a mixed model framework using weighted deregressed estimated breeding values as a dependent variable on 1,883 phenotyped animals that were ≥87.5% Holstein-Friesian. Putative susceptibility quantitative trait loci (QTL) were identified on Bos taurus autosome 1, 3, 5, 6, 8, 9, 10, 11, 13, 14, 18, 21, 23, 25, 26, 27, and 29; mapping the most significant SNP to genes within and overlapping these QTL revealed that the most significant associations were with the 10 functional candidate genes KALRN, ZBTB20, LPP, SLA2, FI3A1, LRCH3, DNAJC6, ZDHHC14, SNX1, and HAS2. Pathway analysis failed to reveal significantly enriched biological pathways, when both bovine-specific pathway data and human ortholog data were taken into account. The existence of genetic variation for MAP susceptibility in a large data set of dairy cows signifies the potential of breeding programs for reducing MAP susceptibility. Furthermore, the identification of susceptible QTL facilitates greater biological understanding of bovine paratuberculosis and potential therapeutic targets for future investigation. The novel molecular similarities identified between bovine paratuberculosis and human inflammatory bowel disease suggest potential for human therapeutic interventions to be translated to veterinary medicine and vice versa.

Bovine Genome Database: Tools for Mining the Bos taurus Genome

The Bovine Genome Database (BGD; http://bovinegenome.org ) is a web-accessible resource that supports bovine genomics research by providing genome annotation and data mining tools. BovineMine is a tool within BGD that integrates BGD data, including the genome, genes, precomputed gene expression levels and variant consequences, with external data sources that include quantitative trait loci (QTL), orthologues, Gene Ontology, gene interactions, and pathways. BovineMine enables researchers without programming skills to create custom integrated datasets for use in downstream analyses. This chapter describes how to enhance a bovine genomics project using the Bovine Genome Database, with data mining examples demonstrating BovineMine.

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.

Which phenotypic traits of resistance should be improved in cattle to control paratuberculosis dynamics in a dairy herd: a modelling approach

Paratuberculosis is a worldwide disease causing production losses in dairy cattle herds. Variability of cattle response to exposure to Mycobacterium avium subsp. paratuberculosis (Map) has been highlighted. Such individual variability could influence Map spread at larger scale. Cattle resistance to paratuberculosis has been shown to be heritable, suggesting genetic selection could enhance disease control. Our objective was to identify which phenotypic traits characterising the individual course of infection influence Map spread in a dairy cattle herd. We used a stochastic mechanistic model. Resistance consisted in the ability to prevent infection and the ability to cope with infection. We assessed the effect of varying (alone and combined) fourteen phenotypic traits characterising the infection course. We calculated four model outputs 25 years after Map introduction in a naïve herd: cumulative incidence, infection persistence, and prevalence of infected and affected animals. A cluster analysis identified influential phenotypes of cattle resistance. An ANOVA quantified the contribution of traits to model output variance. Four phenotypic traits strongly influenced Map spread: the decay in susceptibility with age (the most effective), the quantity of Map shed in faeces by high shedders, the incubation period duration, and the required infectious dose. Interactions contributed up to 12% of output variance, highlighting the expected added-value of improving several traits simultaneously. Combinations of the four most influential traits decreased incidence to less than one newly infected animal per year in most scenarios. Future genetic selection should aim at improving simultaneously the most influential traits to reduce Map spread in cattle populations.

Genetic susceptibility to infectious diseases: Current status and future perspectives from genome-wide approaches

Genome-wide association studies (GWASs) have been widely applied to identify genetic factors that affect complex diseases or traits. Presently, the GWAS Catalog includes > 2800 human studies. Of these, only a minority have investigated the susceptibility to infectious diseases or the response to therapies for the treatment or prevention of infections. Despite their limited application in the field, GWASs have provided valuable insights by pinpointing associations to both innate and adaptive immune response loci, as well as novel unexpected risk factors for infection susceptibility. Herein, we discuss some issues and caveats of GWASs for infectious diseases, we review the most recent findings ensuing from these studies, and we provide a brief summary of selected GWASs for infections in non-human mammals. We conclude that, although the general trend in the field of complex traits is to shift from GWAS to next-generation sequencing, important knowledge on infectious disease-related traits can be still gained by GWASs, especially for those conditions that have never been investigated using this approach. We suggest that future studies will benefit from the leveraging of information from the host's and pathogen's genomes, as well as from the exploration of models that incorporate heterogeneity across populations and phenotypes. Interactions within HLA genes or among HLA variants and polymorphisms located outside the major histocompatibility complex may also play an important role in shaping the susceptibility and response to invading pathogens.

•

Relatively few GWASs for infectious diseases were performed.

•

Phenotype heterogeneity and case/control misclassification can affect GWAS power.

•

Adaptive and innate immunity loci were identified in several infectious disease GWASs.

•

Unexpected loci (e.g., lncRNAs) were also associated with infection susceptibility.

•

GWASs should integrate host and pathogen diversity and use complex association models.

Identification of additional loci associated with antibody response to Mycobacterium avium ssp. Paratuberculosis in cattle by GSEA-SNP analysis

Mycobacterium avium subsp. paratuberculosis: (MAP) causes a contagious chronic infection results in Johne's disease in a wide range of animal species, including cattle. Several genome-wide association studies (GWAS) have been carried out to identify loci putatively associated with MAP susceptibility by testing each marker separately and identifying SNPs that show a significant association with the phenotype, while SNP with modest effects are usually ignored. The objective of this study was to identify modest-effect genes associated with MAP susceptibility using a pathway-based approach. The Illumina BovineSNP50 BeadChip was used to genotype 966 Holstein cows, 483 positive and 483 negative for antibody response to MAP, data were then analyzed using novel SNP-based Gene Set Enrichment Analysis (GSEA-SNP) and validated with Adaptive Rank Truncated Product methodology. An allele-based test was carried out to estimate the statistical association for each marker with the phenotype, subsequently SNPs were mapped to the closest genes, considering for each gene the single variant with the highest value within a window of 50 kb, then pathway-statistics were tested using the GSEA-SNP method. The GO biological process "embryogenesis and morphogenesis" was most highly associated with antibody response to MAP. Within this pathway, five genes code for proteins which play a role in the immune defense relevant to response to bacterial infection. The immune response genes identified would not have been considered using a standard GWAS, thus demonstrating that the pathway approach can extend the interpretation of genome-wide association analyses and identify additional candidate genes for target traits.

Genomic Tools and Animal Health

Animals have been selected to improve their productivity in order to increase the profitability to the producer. In this scenario, not much attention was given to health traits. As a consequence of that, selection was made for animals with higher production and a shortened productive life. In addition to that, the intense production system used in livestock has forced animals to be exposed to higher pathogen loads, therefore predisposing them to infections. Infectious diseases are known to be caused by micro-organisms that are able to infect and colonize the host, affecting their physiological functions and causing problems in their production and on animal welfare. Even with the best management practices, diseases are still the most important cause of economic losses in the animal industry. In this review article we have addressed the new tools that could be used to select animals to better cope with diseases and pathogens.

Bovine Genome Database: new tools for gleaning function from the Bos taurus genome

We report an update of the Bovine Genome Database (BGD) (http://BovineGenome.org). The goal of BGD is to support bovine genomics research by providing genome annotation and data mining tools. We have developed new genome and annotation browsers using JBrowse and WebApollo for two Bos taurus genome assemblies, the reference genome assembly (UMD3.1.1) and the alternate genome assembly (Btau_4.6.1). Annotation tools have been customized to highlight priority genes for annotation, and to aid annotators in selecting gene evidence tracks from 91 tissue specific RNAseq datasets. We have also developed BovineMine, based on the InterMine data warehousing system, to integrate the bovine genome, annotation, QTL, SNP and expression data with external sources of orthology, gene ontology, gene interaction and pathway information. BovineMine provides powerful query building tools, as well as customized query templates, and allows users to analyze and download genome-wide datasets. With BovineMine, bovine researchers can use orthology to leverage the curated gene pathways of model organisms, such as human, mouse and rat. BovineMine will be especially useful for gene ontology and pathway analyses in conjunction with GWAS and QTL studies.

Genome-wide association study of susceptibility to infection by Mycobacterium avium subspecies paratuberculosis in Holstein cattle

Paratuberculosis, or Johne's disease, is a chronic, granulomatous, gastrointestinal tract disease of cattle and other ruminants caused by the bacterium Mycobacterium avium, subspecies paratuberculosis (MAP). Control of Johne's disease is based on programs of testing and culling animals positive for infection with MAP while concurrently modifying management to reduce the likelihood of infection. The current study is motivated by the hypothesis that genetic variation in host susceptibility to MAP infection can be dissected and quantifiable associations with genetic markers identified. For this purpose, a case-control, genome-wide association study was conducted using US Holstein cattle phenotyped for MAP infection using a serum ELISA and/or fecal culture test. Cases included cows positive for either serum ELISA, fecal culture or both. Controls consisted of animals negative for the serum ELISA test or both serum ELISA and fecal culture when both were available. Controls were matched by herd and proximal birth date with cases. A total of 856 cows (451 cases and 405 controls) were used in initial discovery analyses, and an additional 263 cows (159 cases and 104 controls) from the same herds were used as a validation data set. Data were analyzed in a single marker analysis controlling for relatedness of individuals (GRAMMAR-GC) and also in a Bayesian analysis in which multiple marker effects were estimated simultaneously (GenSel). For the latter, effects of non-overlapping 1 Mb marker windows across the genome were estimated. Results from the two discovery analyses were generally concordant; however, discovery results were generally not well supported in analysis of the validation data set. A combined analysis of discovery and validation data sets provided strongest support for SNPs and 1 Mb windows on chromosomes 1, 2, 6, 7, 17 and 29.

Impact of host genetics on susceptibility and resistance to Mycobacterium avium subspecies Paratuberculosis infection in domestic ruminants

Johne's disease or Paratuberculosis has emerged as major infectious disease of animals in general and domestic livestock in particular on global basis. There have been major initiatives in developed countries for the control of this incurable malady of animals and human beings alike (inflammatory bowel disease or Crohn's disease). Disease has not received similar attention due to inherent complexities of disease, diagnosis and control, in resource poor counties around the world. However, the rich genetic diverstiy of the otherwise low productive animal population offers opportunity for the control of Johne's disease and improve per animal productivity. Present review aims to gather and compile information available on genetics or resistance to Johne's disease and its future exploitation by resource poor countries rich in animal diversity. This review will also help to create awareness and share knowledge and experience on prevalence and opportunities for control of Johne's disease in the livestock population to boost per animal productivity among developing and poor countries of the world. Breeding of animals for disease resistance provides good, safe, effective and cheaper way of controlling Johne's disease in animals, with especial reference to domestic livestock of developing and poor countries. Study will help to establish better understanding of the correlation between host cell factors and resistance to MAP infection which may have ultimately help in the control of Johne's disease in future.

Genomics of complex traits

The analysis of complex genetic traits, including mapping and identification of causative genes, has long been an enigma of genetic biology, whether in the animal sciences or in medical sciences. Traits of agricultural interest and traits of medical interest are often under the influence of both environmental factors and multiple genes, each with modest contributions to the total variance in the trait. Although the number of known mutations underlying complex traits is still relatively small, advances in genomics have greatly enhanced traditional pathways to their analysis and gene mining. The candidate gene approach, linkage analysis, and association studies are all significantly more powerful with recent advances in genome mapping, sequencing, and analysis of individual variation. Avenues to gene discovery are discussed with emphasis on genome wide association studies (GWAS) and the use of single nucleotide polymorphisms (SNPs) as revealed by increasingly powerful commercially available microarrays.

Meta-analysis of two genome-wide association studies of bovine paratuberculosis

Background

Bovine paratuberculosis (ParaTB) also known as Johne's disease, is a contagious fatal disease resulting from infection by Mycobacterium avium subspecies paratuberculosis (MAP). Previous studies have identified loci associated with ParaTB using different measurements to define cases and controls. The objective of this study was to combine the data from two recent studies to identify genetic loci associated with MAP tissue infection and humoral immune response, defined by MAP ELISA-positive cattle, by comparing cases and control animals for one or both measures of infection.

Methodology/Principal Findings

The two populations used for the association analyses were a cohort of MAP tissue infected animals and control Holstein cows from the USA and the second cohort composed of ELISA-positive and ELISA-negative Holstein cows from Italy. Altogether 1190 cattle were genotyped with the Illumina BovineSNP50 BeadChip. SNP markers were removed if the minor allele frequency <0.01 or genotyping failure was >5%. Animals were removed with >5% genotyping failure. Whole genome association analyses were conducted with the GRAMMAR-CG method using two different definitions of control populations.

Conclusion/Significance

The analyses identified several loci (P<5 e-05) associated with ParaTB, defined by positive ELISA and presence of bacteria in tissue compared to ELISA and tissue negative animals, on chromosomes 1, 12 and 15 and one unassigned SNP. These results confirmed associations on chromosome 12 and the unassigned SNP with ParaTB which had been found in the Italian population alone. Furthermore, several additional genomic regions were found associated with ParaTB when ELISA and tissue positive animals were compared with tissue negative samples. These loci were on chromosomes 1, 6, 7, 13, 16, 21,23 and 25 (P<5 e-05). The results clearly indicate the importance of the phenotype definition when seeking to identify markers associated with different disease responses.

Genetic susceptibility to paratuberculosis

Multiple studies indicate that host animal genetics play a role in susceptibility to Mycobacterium avium subsp. paratuberculosis (MAP) infection. However, due to differences in methods used to define MAP-infected animals and controls and differences in methods of genetic analysis, there is as yet no clear consensus on the genes or markers to reliably define the MAP infection susceptibility of any animal species. Meta-analysis of combined studies and larger studies will help resolve the situation in the coming years.

Short communication: refinement of genetic regions associated with Mycobacterium avium subspecies paratuberculosis tissue infection and tolerance to Johne's disease

Johne's disease is a highly transmissible bacterial disease caused by Mycobacterium avium ssp. paratuberculosis (MAP). The objective of this study was to refine the locus associated with MAP tissue infection and the locus associated with tolerance to Johne's disease. Using a genome-wide association analysis, single nucleotide polymorphisms associated with MAP tissue infection and tolerance to Johne's disease on Bos taurus autosome (BTA)3 and BTA15, respectively, have previously been identified. A 235-kb region on BTA3 was evaluated with 42 single nucleotide polymorphisms, and a 193-kb region on BTA15 was evaluated with 54 single nucleotide polymorphisms in a group of 209 Holstein cows. Using a single marker association analysis and haplotype tests, we refined a region of 10.6 kb on BTA3 as being associated with MAP tissue infection and a region of 6.5 kb on BTA15 as being associated with tolerance to Johne's disease.

Evolution of the bovine TLR gene family and member associations with Mycobacterium avium subspecies paratuberculosis infection

Members of the Toll-like receptor (TLR) gene family occupy key roles in the mammalian innate immune system by functioning as sentries for the detection of invading pathogens, thereafter provoking host innate immune responses. We utilized a custom next-generation sequencing approach and allele-specific genotyping assays to detect and validate 280 biallelic variants across all 10 bovine TLR genes, including 71 nonsynonymous single nucleotide polymorphisms (SNPs) and one putative nonsense SNP. Bayesian haplotype reconstructions and median joining networks revealed haplotype sharing between Bos taurus taurus and Bos taurus indicus breeds at every locus, and specialized beef and dairy breeds could not be differentiated despite an average polymorphism density of 1 marker/158 bp. Collectively, 160 tagSNPs and two tag insertion-deletion mutations (indels) were sufficient to predict 100% of the variation at 280 variable sites for both Bos subspecies and their hybrids, whereas 118 tagSNPs and 1 tagIndel predictively captured 100% of the variation at 235 variable sites for B. t. taurus. Polyphen and SIFT analyses of amino acid (AA) replacements encoded by bovine TLR SNPs indicated that up to 32% of the AA substitutions were expected to impact protein function. Classical and newly developed tests of diversity provide strong support for balancing selection operating on TLR3 and TLR8, and purifying selection acting on TLR10. An investigation of the persistence and continuity of linkage disequilibrium (r²≥0.50) between adjacent variable sites also supported the presence of selection acting on TLR3 and TLR8. A case-control study employing validated variants from bovine TLR genes recognizing bacterial ligands revealed six SNPs potentially eliciting small effects on susceptibility to Mycobacterium avium spp paratuberculosis infection in dairy cattle. The results of this study will broadly impact domestic cattle research by providing the necessary foundation to explore several avenues of bovine translational genomics, and the potential for marker-assisted vaccination.

Candidate gene and genome-wide association studies of Mycobacterium avium subsp. paratuberculosis infection in cattle and sheep: a review

Paratuberculosis (Johne's disease), caused by Mycobacterium avium subspecies paratuberculosis, is responsible for significant economic losses in livestock industries worldwide. This organism is also of public health concern due to an unconfirmed link to Crohn's disease. Susceptibility to paratuberculosis has been suggested to have a genetic component. In livestock, a number of candidate genes have been studied, selected on their association to susceptibility in other mycobacterial diseases, their known role in disease pathogenesis or links to susceptibility of humans to Crohn's disease. These genes include solute carrier family 11 member 1 (SLC11A1, formerly NRAMP1), toll-like receptors, caspase associated recruitment domain 15 (CARD15, formerly NOD2), major histocompatibility complex (MHC) and cytokines (interleukin-10 and interferon-gamma) and their receptors. Genome wide association studies have attempted to confirm associations found and identify new genes involved in pathogenesis and susceptibility. There are a number of limitations and difficulties in these approaches, some peculiar to paratuberculosis but others generally applicable to identification of genetic associations for complex traits. The technical approaches and available information for paratuberculosis have expanded rapidly, particularly relating to sheep and cattle. Here we review the current published evidence for a genetic association with paratuberculosis susceptibility, technological advances that have progressed the field and potential avenues for future research.