Genetic loci involved in antibody response to Mycobacterium avium ssp. paratuberculosis in cattle

Genetic Parameters and Genetic Trends for Johne's Disease in U.S. Holsteins: An Updated Study

Johne's disease (JD) is an infectious enteric disease in ruminants, causing substantial economic loss annually worldwide. This work aimed to estimate JD's genetic parameters and the phenotypic and genetic trends by incorporating recent data. It also explores the feasibility of a national genetic evaluation for JD susceptibility in Holstein cattle in the United States (US). The data were extracted from a Johne's disease data repository, maintained at the Council on Dairy Cattle Breeding (CDCB), and initially supplied by 2 dairy records processing centers. The data comprised 365,980 Holstein cows from 1,048 herds participating in a voluntary control program for JD. Two protocol kits, Idexx Paratuberculosis Screening Ab Test (IDX) and Parachek® 2 (PCK), were used to analyze milk samples with the Enzyme-Linked Immunosorbent Assay (ELISA) technique. Test results from the first 5 parities were considered. An animal was considered infected if it had at least one positive outcome. The overall average of JD incidence was 4.72% in these US Holstein cattle. Genotypes of 78,964 SNP markers were used for 25,000 animals randomly selected from the phenotyped population. Variance components and genetic parameters were estimated based on 3 models, namely, a pedigree-only threshold model (THR), a single-step threshold model (ssTHR), and a single-step linear model (ssLR). The posterior heritability estimates of JD susceptibility were low to moderate: 0.11 to 0.16 based on the 2 threshold models and 0.05 to 0.09 based on the linear model. The average reliability of estimated breeding values of JD susceptibility using single-step analysis for animals with or without phenotypes varied from 0.18 (THR) to 0.22 (ssLR) for IDX and from 0.14 (THR) to 0.18 (ssTHR and ssLR) for PCK. Despite no prior direct genetic selection against JD, the estimated genetic trends of JD susceptibility were negative and highly significant (P-value <0.01). The correlations of bulls' predicted transmitting abilities with economically important traits such as milk yield, milk protein, milk fat, somatic cell score, and mastitis were low, indicating a non-overlapping genetic selection process with traits in current genetic evaluations. Our results suggest the feasibility of reducing the JD incidence rate by incorporating it into the national genetic evaluation programs.

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

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

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.

Analysis of runs of homozygosity of cattle living in different climate zones

Aim of this study was to analyze the distribution and characteristics of runs of homozygosity in Bos taurus taurus and Bos taurus indicus breeds, as well as their crosses, farmed all around the world. With this aim in view, we used single-nucleotide polymorphisms (SNP) genotypes for 3,263 cattle belonging to 204 different breeds. After quality control, 23,311 SNPs were retained for the analysis. Animals were divided into seven different groups: 1) continental taurus, 2) temperate taurus, 3) temperate indicus, 4) temperate composite, 5) tropical taurus, 6) tropical indicus, and 7) tropical composite. The climatic zones were created according to the latitude of the breeds' country of origin: i) continental, latitude ≥ 45°; ii) temperate, 45°< Latitude >23.26°; iii) tropics, latitude ≤ 23.26°. Runs of homozygosity were computed as 15 SNPs spanning in at least 2 Mb; number of ROH per animal (nROH), average ROH length (meanMb), and ROH-based inbreeding coefficients (FROH) were also computed. Temperate indicus showed the largest nROH, whereas Temperate taurus the lowest value. Moreover, the largest meanMb was observed for Temperate taurus, whereas the lowest value for Tropics indicus. Temperate indicus breeds showed the largest FROH values. Genes mapped in the identified ROH were reported to be associated with the environmental adaptation, disease resistance, coat color determinism, and production traits. Results of the present study confirmed that runs of homozygosity could be used to identify genomic signatures due to both artificial and natural selection.

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.

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.

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.

Test positivity for Maedi-Visna virus and Mycobacterium avium ssp. paratuberculosis in Sarda ewes: Effects on milk composition and coagulation traits and heritability estimates for susceptibility

Maedi-Visna virus (MVV) and Mycobacterium avium ssp. paratuberculosis (MAP) are two pathogens that cause chronic, production-limiting diseases in dairy sheep. Although they are present worldwide, there are no detailed reports on their actual effects on milk traits in the literature. This study was designed to investigate the effects of test positivity to MVV and MAP on ovine milk yield, composition and coagulation properties, and curd-firming over time (CF_t) variables in clinically healthy animals at the field level. The additive genetic variation and heritabilities of MVV and MAP positivity were also estimated. Milk samples were collected from 1,079 Sarda sheep kept on 23 farms, and pedigree information was obtained from the flock book. Milk yield was also recorded on the sampling date. Positivity for MVV and MAP was determined from milk samples using indirect ELISA test kits. Milk composition traits were measured by spectroscopy, milk coagulation properties were measured with a Formagraph (Foss Italia, Padua, Italy), and CF_t traits were calculated using the data from the Formagraph diagram. The effects of MVV and MAP positivity on milk traits were determined through a set of mixed linear models, which took into account various sources of variation, such as days in milk, parity, and flock effects, and included the effects (positive or negative) of the 2 pathogens. A Bayesian threshold sire model with sire relationship was used to estimate genetic variation and heritability. The overall animal prevalence of MVV-positive ewes was 43.6%; on only 1 farm of the 23 tested were all sampled ewes negative. An overall animal prevalence of 10.6% was recorded for MAP, with 4 farms at 0%. Positivity for MVV significantly affected the logarithmic score of the bacterial count, curd firmness after 30 min and 45 min, and the curd-firming instant rate constant. We found significant effects of MAP infection on milk composition, pH, and rennet coagulation time. The mean of the posterior distributions of heritability estimates on the liability scale was 0.15 for MAP and 0.07 for MVV. Our results demonstrate that only a few traits are negatively affected by MVV and MAP positivity, and that there is exploitable genetic variation in MVV and MAP susceptibility in dairy sheep.

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.

Genome-Wide Association Study in Mexican Holstein Cattle Reveals Novel Quantitative Trait Loci Regions and Confirms Mapped Loci for Resistance to Bovine Tuberculosis

Simple Summary

Bovine tuberculosis is an infectious disease of cattle caused by Mycobacterium bovis characterized by the formation of tubercles in any organ or tissue. Bovine tuberculosis represents a significant veterinary and public health problem in many parts of the world. It is zoonotic, transmitted to humans through consumption of infected milk and other cattle products. Although many factors influence infection and progression of the disease, there must be an important host genetic component that explains why some animals get sick and others remain healty. We present evidence of genetic variants associated with resistance to tuberculosis in Mexican Holstein dairy cattle using a case-control approach with a selective DNA pooling. Here, we identified novel quantitative trait loci regions harboring genes involved in Mycobacterium spp. immune response. This is a first screening about resistance to tuberculosis infection on Mexican Holstein cattle based on a dense single nucleotide polymorphism chip. The identified genes belong to both, the already known, and the undisclosed quantitative trait loci regions.

Abstract

Bovine tuberculosis (bTB) is a disease of cattle that represents a risk to public health and causes severe economic losses to the livestock industry. Recently, genetic studies, like genome-wide association studies (GWAS) have greatly improved the investigation of complex diseases identifying thousands of disease-associated genomic variants. Here, we present evidence of genetic variants associated with resistance to TB in Mexican dairy cattle using a case-control approach with a selective DNA pooling experimental design. A total of 154 QTLRs (quantitative trait loci regions) at 10% PFP (proportion of false positives), 42 at 5% PFP and 5 at 1% PFP have been identified, which harbored 172 annotated genes. On BTA13, five new QTLRs were identified in the MACROD2 and KIF16B genes, supporting their involvement in resistance to bTB. Six QTLRs harbor seven annotated genes that have been previously reported as involved in immune response against Mycobacterium spp: BTA (Bos taurus autosome) 1 (CD80), BTA3 (CTSS), BTA 3 (FCGR1A), BTA 23 (HFE), BTA 25 (IL21R), and BTA 29 (ANO9 and SIGIRR). We identified novel QTLRs harboring genes involved in Mycobacterium spp. immune response. This is a first screening for resistance to TB infection on Mexican dairy cattle based on a dense SNP (Single Nucleotide Polymorphism) chip.

Genome-wide association study of Mycobacterium avium subspecies Paratuberculosis infection in Chinese Holstein

Background

Paratuberculosis is a contagious, chronic and enteric disease in ruminants, which is caused by Mycobacterium avium subspecies paratuberculosis (MAP) infection, resulting in enormous economic losses worldwide. There is currently no effective cure for MAP infection or a vaccine, it is thus important to explore the genetic variants that contribute to host susceptibility to infection by MAP, which may provide a better understanding of the mechanisms of paratuberculosis and benefit animal genetic improvement. Herein we performed a genome-wide association study (GWAS) to identify genomic regions and candidate genes associated with susceptibility to MAP infection in dairy cattle.

Results

Using Illumina Bovine 50 K (54,609 SNPs) and GeneSeek HD (138,893 SNPs) chips, two analytical approaches were performed, GRAMMAR-GC and ROADTRIPS in 937 Chinese Holstein cows, among which individuals genotyped by the 50 K chip were imputed to HD SNPs with Beagle software. Consequently, 15 and 11 significant SNPs (P < 5 × 10^− 5) were identified with GRAMMAR-GC and ROADTDRIPS, respectively. A total of 10 functional genes were in proximity to (i.e., within 1 Mb) these SNPs, including IL4, IL5, IL13, IRF1, MyD88, PACSIN1, DEF6, TDP2, ZAP70 and CSF2. Functional enrichment analysis showed that these genes were involved in immune related pathways, such as interleukin, T cell receptor signaling pathways and inflammatory bowel disease (IBD), implying their potential associations with susceptibility to MAP infection. In addition, by examining the publicly available cattle QTLdb, a previous QTL for MAP was found to be overlapped with one of regions detected currently at 32.5 Mb on BTA23, where the TDP2 gene was anchored.

Conclusions

In conclusion, we identified 26 SNPs located on 15 chromosomes in the Chinese Holstein population using two GWAS strategies with high density SNPs. Integrated analysis of GWAS, biological functions and the reported QTL information helps to detect positional candidate genes and the identification of regions associated with susceptibility to MAP traits in dairy cattle.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5385-3) contains supplementary material, which is available to authorized users.

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.

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

Johne's disease is a contagious bacterial infection of cattle caused by ssp. (). A previous genome-wide association analysis (GWAA) in Holstein cattle identified QTL on BTA3 and BTA9 that were highly associated (P < 5 × 10) and on BTA1, BTA16, and BTA21 that were moderately associated (P < 5 × 10) with Map tissue infection. The objectives of this study were to validate previous GWAA results in Jersey cattle ( = 57), Holstein cattle from the Pacific Northwest (PNW, = 205) and a combined Holstein population from the PNW and the Northeast (PNW + NE, = 423), and also identify new loci associated with tissue infection. DNA was genotyped using the Illumina BovineSNP50 BeadChip, and the PNW + NE data was also imputed to whole genome sequence level using Run4 of the 1000 Bull Genomes project with Beagle v 4.1 and FImpute. Cases were ileocecal node positive and controls were negative for by quantitative PCR (qPCR). Individuals were removed for SNP call rate < 90%, and SNP were removed for genotype call rate < 90% or minor allele frequency < 1%. For the Jersey, PNW, and PNW + NE, GWAA were conducted using an allelic dosage model. For the PNW and the PNW + NE, an additional efficient mixed-model association eXpedited (EMMAX) analysis was performed using additive, dominance and recessive models. Seven QTL on BTA22 were identified in the Jersey population with the most significant ( = 4.45 × 10) located at 21.7 megabases (Mb). Six QTL were associated in the PNW and the PNW + NE analyses, including a QTL previously identified on BTA16 in the NE population. The most significant locus for the PNW was located on BTA21 at 61 Mb ( = 8.61 × 10) while the most significant locus for the PNW + NE was on BTA12 at 90 Mb ( = 2.33 × 10). No additional QTL were identified with the imputed GWAA. Putative positional candidate genes were identified within 50 kb 5' and 3' of each QTL. Two positional candidate genes were identified in Jersey cattle, 1 identified in the PNW and 8 in the PNW + NE populations. Many identified positional candidate genes are involved in signal transduction, have immunological functions, or have putative functional relevance in entry into host cells. This study supported 2 previously identified SNP within a QTL on BTA16 and identified 16 new QTL, including 2 found in the PNW and the PNW+NE, associated with tissue infection.

Responses of Bovine Innate Immunity to Mycobacterium avium subsp. paratuberculosis Infection Revealed by Changes in Gene Expression and Levels of MicroRNA

Paratuberculosis in cattle is a chronic granulomatous gastroenteritis caused by Mycobacterium avium subsp. paratubercolosis (MAP) which is endemic worldwide. In dairy herds, it is responsible for huge economic losses. However, current diagnostic methods do not detect subclinical infection making control of the disease difficult. The identification of MAP infected animals during the sub-clinical phase of infection would play a key role in preventing the dissemination of the pathogen and in reducing transmission. Gene expression and circulating microRNA (miRNA) signatures have been proposed as biomarkers of disease both in the human and veterinary medicine. In this paper, gene expression and related miRNA levels were investigated in cows positive for MAP, by ELISA and culture, in order to identify potential biomarkers to improve diagnosis of MAP infection. Three groups, each of 5 animals, were used to compare the results of gene expression from positive, exposed and negative cows. Overall 258 differentially expressed genes were identified between unexposed, exposed, but ELISA negative and positive groups which were involved in biological functions related to inflammatory response, lipid metabolism and small molecule biochemistry. Differentially expressed miRNA was also found among the three groups: 7 miRNAs were at a lower level and 2 at a higher level in positive animals vs unexposed animals, while 5 and 3 miRNAs were respectively reduced and increased in the exposed group compared to the unexposed group. Among the differentially expressed miRNAs 6 have been previously described as immune-response related and two were novel miRNAs. Analysis of the miRNA levels showed correlation with expression of their target genes, known to be involved in the immune process. This study suggests that miRNA expression is affected by MAP infection and play a key role in tuning the host response to infection. The miRNA and gene expression profiles may be biomarkers of infection and potential diagnostic of MAP infection earlier than the current ELISA based diagnostic tests.

An alternative experimental case-control design for genetic association studies on bovine mastitis

The possibility of using genetic control strategies to increase disease resistance to infectious diseases relies on the identification of markers to include in the breeding plans. Possible incomplete exposure of mastitis-free (control) animals, however, is a major issue to find relevant markers in genetic association studies for infectious diseases. Usually, designs based on elite dairy sires are used in association studies, but an epidemiological case-control strategy, based on cows repeatedly field-tested could be an alternative for disease traits. To test this hypothesis, genetic association results obtained in the present work from a cohort of Italian Holstein cows tested for mastitis over time were compared with those from a previous genome-wide scan on Italian Holstein sires genotyped with 50k single nucleotide polymorphisms for de-regressed estimated breeding values for somatic cell counts (SCCs) on Bos taurus autosome (BTA6) and BTA14. A total of 1121 cows were selected for the case-control approach (cases=550, controls=571), on a combination of herd level of SCC incidence and of within herd individual level of SCC. The association study was conducted on nine previously identified markers, six on BTA6 and four on BTA14, using the R statistical environment with the 'qtscore' function of the GenABEL package, on high/low adjusted linear score as a binomial trait. The results obtained in the cow cohort selected on epidemiological information were in agreement with those obtained from the previous sire genome-wide association study (GWAS). Six out of the nine markers showed significant association, four on BTA14 (rs109146371, rs109234250, rs109421300, rs109162116) and two on BTA6 (rs110527224 and rs42766480). Most importantly, using mastitis as a case study, the current work further validated the alternative use of historical field disease data in case-control designs for genetic analysis of infectious diseases in livestock.

Signatures of selection identify loci associated with fat tail in sheep

Identification of genomic regions that have been targets of selection for phenotypic traits is one of the most challenging applications of dense marker panels in animal genetics. In this study, a genomewide scan using approximately 50,000 SNP was performed in an attempt to identify genomic regions associated with fat deposition in sheep, the importance of this not only being limited to livestock facing future climate changes but also for elucidating the physiology of lipid metabolism. The genotyping results obtained with the Ovine SNP50K BeadChip in 2 fat tail breeds were compared with those obtained in 13 thin tail breeds. Direct sequencing of the annotated genes located in proximity to the markers with opposite allele frequency in thin tail vs. fat tail sheep gave additional SNP of interest. To further confirm the results of the genomewide scan, we genotyped the SNP within these genes in the 2 groups of sheep. A missense mutation in the gene, with different allele frequency in the 2 groups, was detected. The results indicated and as the most probable genes involved in the fat tail phenotype.

A genome-wide association study for genetic susceptibility to Mycobacterium bovis infection in dairy cattle identifies a susceptibility QTL on chromosome 23

Background

Bovine tuberculosis (bTB) infection in cattle is a significant economic concern in many countries, with annual costs to the UK and Irish governments of approximately €190 million and €63 million, respectively, for bTB control. The existence of host additive and non-additive genetic components to bTB susceptibility has been established.

Methods

Two approaches i.e. single-SNP (single nucleotide polymorphism) regression and a Bayesian method were applied to genome-wide association studies (GWAS) using high-density SNP genotypes (n = 597,144 SNPs) from 841 dairy artificial insemination (AI) sires. Deregressed estimated breeding values for bTB susceptibility were used as the quantitative dependent variable. Network analysis was performed using the quantitative trait loci (QTL) that were identified as significant in the single-SNP regression and Bayesian analyses separately. In addition, an identity-by-descent analysis was performed on a subset of the most prolific sires in the dataset that showed contrasting prevalences of bTB infection in daughters.

Results

A significant QTL region was identified on BTA23 (P value >1 × 10⁻⁵, Bayes factor >10) across all analyses. Sires with the minor allele (minor allele frequency = 0.136) for this QTL on BTA23 had estimated breeding values that conferred a greater susceptibility to bTB infection than those that were homozygous for the major allele. Imputation of the regions that flank this QTL on BTA23 to full sequence indicated that the most significant associations were located within introns of the FKBP5 gene.

Conclusions

A genomic region on BTA23 that is strongly associated with host susceptibility to bTB infection was identified. This region contained FKBP5, a gene involved in the TNFα/NFκ-B signalling pathway, which is a major biological pathway associated with immune response. Although there is no study that validates this region in the literature, our approach represents one of the most powerful studies for the analysis of bTB susceptibility to date.

Electronic supplementary material

The online version of this article (doi:10.1186/s12711-016-0197-x) contains supplementary material, which is available to authorized users.

Genomic scan for identifying candidate genes for paratuberculosis resistance in sheep

Breeding objectives relating to health, functional traits and welfare need to receive priority in the research programs and selection schemes, but very few reports are available on natural resistant genotypes in livestock, where some important diseases cause severe economic losses and pose serious zoonotic threats. In this study, diagnosis of paratuberculosis was performed on 759 adult sheep, from a single flock, with the serum antibody enzyme-linked immunosorbent assay; 100 sheep were selected among the extreme divergent animals for the S/P ratio obtained from the test, and were genotyped on the Illumina Ovine SNP50K BeadChip. A genome-wide scan was then performed on the individual marker genotypes, in the attempt to identify genomic regions associated with disease resistance in sheep. For each marker, the allelic substitution effect was calculated by regressing the S/P value on the number of copies of the reference allele. The position on the OARv3.1 Genome Assembly was searched for 32 markers, which showed a statistically significant allelic substitution effect (Raw P < 0.0006 and FDR P < 0.09). All markers were located within, or close to, annotated genes. Five of these genes, SEMA3, CD109, PCP4, PRDM2 and ITFG2 are referred in literature to play a role in either disease resistance or cell-mediated immune response.

CARD15 Gene Polymorphisms Are Associated with Tuberculosis Susceptibility in Chinese Holstein Cows

Bovine tuberculosis (BTB) is a significant veterinary and financial problem in many parts of the world. Associations between specific host genes and susceptibility to mycobacterial infections, such as tuberculosis, have been reported in several species. The objective of this study was to identify and evaluate the relationship of single-nucleotide polymorphisms (SNPs) in the CARD15 gene with susceptibility to BTB in Chinese Holstein cows. DNA samples from 201 Chinese Holstein cows (103 cases and 98 controls) were collected from Kunming City, Yuxi City, and Dali City in China. SNPs in the CARD15 gene were assessed using polymerase chain reaction (PCR) and restriction fragment length polymorphism-polymerase chain reaction (RFLP-PCR). Case-control association testing and statistical analysis identified six SNPs associated with susceptibility to BTB in Chinese Holstein cows. The frequency of genotypes C/T, A/G, A/G, A/G, C/T, and A/G in E4 (-37), 208, 1644, 1648, 1799, and E10 (+107), respectively, was significantly higher in cases than in controls, and also the alleles C, A, A, G, T, and A, respectively, were associated with a greater relative risk in cases than in controls. The distribution of two haplotypes, TGGACA and CAGACA, was significantly different between cases and controls. Overall, this case-control study suggested that E4 (-37)(C/T), 208(A/G), 1644(A/G), 1648(A/G), 1799(C/T), and E10 (+107)(A/G) in the CARD15 gene were significantly associated with susceptibility to BTB in Chinese Holstein cows and that haplotypes TGGACA and CAGACA could be used as genetic markers in marker-assisted breeding programs for breeding cows with high resistance to BTB.

Phenotype definition is a main point in genome-wide association studies for bovine Mycobacterium avium ssp. paratuberculosis infection status

Paratuberculosis caused by Mycobacterium avium ssp. paratuberculosis (MAP) causes economic losses and is present in dairy herds worldwide. Different studies used different diagnostic tests to detect infection status and are the basis of genome-wide association (GWA) studies with inconsistent results. Therefore, the aim of this study was to identify and compare genomic regions associated with MAP susceptibility in the same cohort of cattle using different diagnostic tests. The GWA study was performed in German Holsteins within a case-control assay using 305 cows tested for MAP by fecal culture and additional with four different commercial ELISA-tests. Genotyping was performed with the Illumina Bovine SNP50 BeadChip. The results using fecal culture or ELISA test led to the identification of different genetic loci. Two single-nucleotide polymorphisms showed significant association with the ELISA-status. However, no significant association for MAP infection could be confirmed. Our results show that the definition of the MAP-phenotype has an important impact on the outcome of GWA studies for paratuberculosis.

Wingless-type MMTV integration site family member 2 (WNT2) gene is associated with resistance to MAP in faecal culture and antibody response in Holstein cattle

Mycobacterium avium subspecies paratuberculosis (MAP) is a pathogenic bacterium responsible for the lethal Johne's disease in cattle. So far, several genome-wide association studies (GWAS) have been carried out to identify chromosomal regions highly associated with Johne's disease. The aim of this study was to investigate the genetic variability within a pool of seven genes (LAMB1, DLD, WNT2, PRDM1, SOCS5, PTGER4 and IL10) indicated by former GWAS/RNA-Seq studies as putatively associated with MAP infections and to achieve a confirmation study of association with paratuberculosis susceptibility in a population of 324 German Holstein cattle (162 cases MAP positive and 162 controls MAP negative) using ELISA and fecal cultural tests. SNP validation and genotyping information are provided, quick methods for allelic discrimination were set up and transcription factor binding analyses were performed. The rs43390642:G>TSNP in the WNT2 promoter region is associated with paratuberculosis susceptibility (P = 0.013), suggesting a protective role of the T allele (P = 0.043; odds ratio 0.50 [0.25-0.97]). The linkage disequilibrium with the DLD rs134692583:A>T might suggest a combined mechanism of action of these neighboring genes in resistance to MAP infection, which is also supported by a significant effect shown by the haplotype DLD(T) /WNT2(T) (P = 0.047). In silico analysis predicted rs43390642:G>T and rs134692583:A>T as essential parts of binding sites for the transcription factors GR, C/EBPβ and GATA-1, hence suggesting a potential influence on WNT2 and DLD gene expression. This study confirmed the region on BTA 4 (UMD 3.1: 50639460-51397892) as involved in tolerance/resistance to Johne's disease. In addition, this study clarifies the involvement of the investigated genes in MAP infection and contributes to the understanding of genetic variability involved in Johne's disease susceptibility.

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.

A genome-wide association study of immune response traits in Canadian Holstein cattle

Background

Breeding for enhanced immune response (IR) has been suggested as a tool to improve inherent animal health. Dairy cows with superior antibody-mediated (AMIR) and cell-mediated immune responses (CMIR) have been demonstrated to have a lower occurrence of many diseases including mastitis. Adaptive immune response traits are heritable, and it is, therefore, possible to breed for improved IR, decreasing the occurrence of disease. The objective of this study was to perform genome-wide association studies to determine differences in genetic profiles among Holstein cows classified as High or Low for AMIR and CMIR. From a total of 680 cows with immune response phenotypes, 163 cows for AMIR (81 High and 82 Low) and 140 for CMIR (75 High and 65 Low) were selectively genotyped using the Illumina Bovine SNP50 BeadChip. Results were validated using an unrelated population of 164 Holstein bulls IR phenotyped for AMIR and 146 for CMIR.

Results

A generalized quasi likelihood score method was used to determine single nucleotide polymorphisms (SNP) and chromosomal regions associated with immune response. After applying a 5% chromosomal false discovery rate, 186 SNPs were significantly associated with AMIR. The majority (93%) of significant markers were on chromosome 23, with a similar peak found in the bull population. For CMIR, 21 SNP markers remained significant. Candidate genes within 250,000 base pairs of significant SNPs were identified to determine biological pathways associated with AMIR and CMIR. Various pathways were identified, including the antigen processing and presentation pathway, important in host defense. Candidate genes included those within the bovine Major Histocompatability Complex such as BoLA-DQ, BoLA-DR and the non-classical BoLA-NC1 for AMIR and BoLA-DQ for CMIR, the complement system including C2 and C4 for AMIR and C1q for CMIR, and cytokines including IL-17A, IL17F for AMIR and IL-17RA for CMIR and tumor necrosis factor for both AMIR and CMIR. Additional genes associated with CMIR included galectins 1, 2 and 3, BCL2 and β-defensin.

Conclusions

The significant genetic variation associated with AMIR and CMIR in this study may imply feasibility to include immune response in genomic breeding indices as an approach to improve inherent animal health.

Genome-wide association analysis and genomic prediction of Mycobacterium avium subspecies paratuberculosis infection in US Jersey cattle

Paratuberculosis (Johne’s disease), an enteric disorder in ruminants caused by Mycobacterium avium subspecies paratuberculosis (MAP), causes economic losses in excess of $200 million annually to the US dairy industry. To identify genomic regions underlying susceptibility to MAP infection in Jersey cattle, a case-control genome-wide association study (GWAS) was performed. Blood and fecal samples were collected from ∼5,000 mature cows in 30 commercial Jersey herds from across the US. Discovery data consisted of 450 cases and 439 controls genotyped with the Illumina BovineSNP50 BeadChip. Cases were animals with positive ELISA and fecal culture (FC) results. Controls were animals negative to both ELISA and FC tests that matched cases on birth date and herd. Validation data consisted of 180 animals including 90 cases (positive to FC) and 90 controls (negative to ELISA and FC), selected from discovery herds and genotyped by Illumina BovineLD BeadChip (∼7K SNPs). Two analytical approaches were used: single-marker GWAS using the GRAMMAR-GC method and Bayesian variable selection (Bayes C) using GenSel software. GRAMMAR-GC identified one SNP on BTA7 at 68 megabases (Mb) surpassing a significance threshold of 5×10⁻⁵. ARS-BFGL-NGS-11887 on BTA23 (27.7 Mb) accounted for the highest percentage of genetic variance (3.3%) in the Bayes C analysis. SNPs identified in common by GRAMMAR-GC and Bayes C in both discovery and combined data were mapped to BTA23 (27, 29 and 44 Mb), 3 (100, 101, 106 and 107 Mb) and 17 (57 Mb). Correspondence between results of GRAMMAR-GC and Bayes C was high (70–80% of most significant SNPs in common). These SNPs could potentially be associated with causal variants underlying susceptibility to MAP infection in Jersey cattle. Predictive performance of the model developed by Bayes C for prediction of infection status of animals in validation set was low (55% probability of correct ranking of paired case and control samples).

Bovine CLEC7A genetic variants and their association with seropositivity in Johne's disease ELISA

Mycobacterium avium ssp. paratuberculosis (MAP) infection in cattle causes significant economic losses to the dairy and beef industries resulting from reduced productivity, premature culling and mortality. Bovine Dectin-1, an important pattern recognition molecule that is able to generate a proinflammatory response by acting alongside Toll like receptor (TLR) signaling, is known to co-operate with TLR2 to specifically activate a macrophage proinflammatory response against mycobacterial infections. Therefore, the goal of this study was to identify single nucleotide polymorphisms (SNPs) in the gene encoding bovine Dectin-1 (CLEC7A) and to assess their association with susceptibility to MAP infection in dairy cattle. Blood and milk samples, collected from commercial dairy operations, were tested for MAP infection using blood and milk ELISAs and a resource population consisting of 197 infected and 242 healthy cattle was constructed. Pooled DNA was used for sequencing and eight single nucleotide polymorphisms (SNPs) were identified. Identified SNPs were genotyped on the resource population using the iPLEX MassARRAY system and statistical analysis was performed using logistic regression fitting the additive and dominance effects of each SNP in the model. Out of a total of eight identified SNPs, five were successfully genotyped, and three out of these five SNPs were found to be in complete linkage. Statistical analysis revealed a strong association between a non-synonymous SNP c.589A>G (p=0.008), and MAP infection status of the resource population inferred by seropositivity in MAP antibody specific ELISAs. This SNP c.589A>G was located in the geneic region that encodes the carbohydrate recognition domain of bovine Dectin-1. Therefore, further investigation of its functional relevance is warranted.

Genome-wide association study to identify chromosomal regions associated with antibody response to Mycobacterium avium subspecies paratuberculosis in milk of Dutch Holstein-Friesians

Heritability of susceptibility to Johne's disease in cattle has been shown to vary from 0.041 to 0.159. Although the presence of genetic variation involved in susceptibility to Johne's disease has been demonstrated, the understanding of genes contributing to the genetic variance is far from complete. The objective of this study was to contribute to further understanding of genetic variation involved in susceptibility to Johne's disease by identifying associated chromosomal regions using a genome-wide association approach. Log-transformed ELISA test results of 265,290 individual Holstein-Friesian cows from 3,927 herds from the Netherlands were analyzed to obtain sire estimated breeding values for Mycobacterium avium subspecies paratuberculosis (MAP)-specific antibody response in milk using a sire-maternal grandsire model with fixed effects for parity, year of birth, lactation stage, and herd; a covariate for milk yield on test day; and random effects for sire, maternal grandsire, and error. For 192 sires with estimated breeding values with a minimum reliability of 70%, single nucleotide polymorphism (SNP) typing was conducted by a multiple SNP analysis with a random polygenic effect fitting 37,869 SNP simultaneously. Five SNP associated with MAP-specific antibody response in milk were identified distributed over 4 chromosomal regions (chromosome 4, 15, 18, and 28). Thirteen putative SNP associated with MAP-specific antibody response in milk were identified distributed over 10 chromosomes (chromosome 4, 14, 16, 18, 19, 20, 21, 26, 27, and 29). This knowledge contributes to the current understanding of genetic variation involved in Johne's disease susceptibility and facilitates control of Johne's disease and improvement of health status by breeding.

New phenotypes for new breeding goals in dairy cattle

Cattle production faces new challenges regarding sustainability with its three pillars - economic, societal and environmental. The following three main factors will drive dairy cattle selection in the future: (1) During a long period, intensive selection for enhanced productivity has deteriorated most functional traits, some reaching a critical point and needing to be restored. This is especially the case for the Holstein breed and for female fertility, mastitis resistance, longevity and metabolic diseases. (2) Genomic selection offers two new opportunities: as the potential genetic gain can be almost doubled, more traits can be efficiently selected; phenotype recording can be decoupled from selection and limited to several thousand animals. (3) Additional information from other traits can be used, either from existing traditional recording systems at the farm level or from the recent and rapid development of new technologies and precision farming. Milk composition (i.e. mainly fatty acids) should be adapted to better meet human nutritional requirements. Fatty acids can be measured through a new interpretation of the usual medium infrared spectra. Milk composition can also provide additional information about reproduction and health. Modern milk recorders also provide new information, that is, on milking speed or on the shape of milking curves. Electronic devices measuring physiological or activity parameters can predict physiological status like estrus or diseases, and can record behavioral traits. Slaughterhouse data may permit effective selection on carcass traits. Efficient observatories should be set up for early detection of new emerging genetic defects. In the near future, social acceptance of cattle production could depend on its capacity to decrease its ecological footprint. The first solution consists in increasing survival and longevity to reduce replacement needs and the number of nonproductive animals. At the individual level, selection on rumen activity may lead to decreased methane production and concomitantly to improved feed efficiency. A major effort should be dedicated to this new field of research and particularly to rumen flora metagenomics. Low input in cattle production is very important and tomorrow's cow will need to adapt to a less intensive production environment, particularly lower feed quality and limited care. Finally, global climate change will increase pathogen pressure, thus more accurate predictors for disease resistance will be required.

Genes controlling vaccine responses and disease resistance to respiratory viral pathogens in cattle

Farm animals remain at risk of endemic, exotic and newly emerging viruses. Vaccination is often promoted as the best possible solution, and yet for many pathogens, either there are no appropriate vaccines or those that are available are far from ideal. A complementary approach to disease control may be to identify genes and chromosomal regions that underlie genetic variation in disease resistance and response to vaccination. However, identification of the causal polymorphisms is not straightforward as it generally requires large numbers of animals with linked phenotypes and genotypes. Investigation of genes underlying complex traits such as resistance or response to viral pathogens requires several genetic approaches including candidate genes deduced from knowledge about the cellular pathways leading to protection or pathology, or unbiased whole genome scans using markers spread across the genome. Evidence for host genetic variation exists for a number of viral diseases in cattle including bovine respiratory disease and anecdotally, foot and mouth disease virus (FMDV). We immunised and vaccinated a cattle cross herd with a 40-mer peptide derived from FMDV and a vaccine against bovine respiratory syncytial virus (BRSV). Genetic variation has been quantified. A candidate gene approach has grouped high and low antibody and T cell responders by common motifs in the peptide binding pockets of the bovine major histocompatibility complex (BoLA) DRB3 gene. This suggests that vaccines with a minimal number of epitopes that are recognised by most cattle could be designed. Whole genome scans using microsatellite and single nucleotide polymorphism (SNP) markers has revealed many novel quantitative trait loci (QTL) and SNP markers controlling both humoral and cell-mediated immunity, some of which are in genes of known immunological relevance including the toll-like receptors (TLRs). The sequencing, assembly and annotation of livestock genomes and is continuing apace. In addition, provision of high-density SNP chips should make it possible to link phenotypes with genotypes in field populations without the need for structured populations or pedigree information. This will hopefully enable fine mapping of QTL and ultimate identification of the causal gene(s). The research could lead to selection of animals that are more resistant to disease and new ways to improve vaccine efficacy.

Genome-wide association study of porcine hematological parameters in a Large White × Minzhu F2 resource population

Hematological traits, which are important indicators of immune function in animals, have been commonly examined as biomarkers of disease and disease severity in humans and animals. Genome-wide significant quantitative trait loci (QTLs) provide important information for use in breeding programs of animals such as pigs. QTLs for hematological parameters (hematological traits) have been detected in pig chromosomes, although these are often mapped by linkage analysis to large intervals making identification of the underlying mutation problematic. Single nucleotide polymorphisms (SNPs) are the common form of genetic variation among individuals and are thought to account for the majority of inherited traits. In this study, a genome-wide association study (GWAS) was performed to detect regions of association with hematological traits in a three-generation resource population produced by intercrossing Large White boars and Minzhu sows during the period from 2007 to 2011. Illumina PorcineSNP60 BeadChip technology was used to genotype each animal and seven hematological parameters were measured (hematocrit (HCT), hemoglobin (HGB), mean corpuscular hemoglobin (MCH), mean corpuscular hemoglobin concentration (MCHC), mean corpuscular volume (MCV), red blood cell count (RBC) and red blood cell volume distribution width (RDW)). Data were analyzed in a three step Genome-wide Rapid Association using the Mixed Model and Regression-Genomic Control (GRAMMAR-GC) method. A total of 62 genome-wide significant and three chromosome-wide significant SNPs associated with hematological parameters were detected in this GWAS. Seven and five SNPs were associated with HCT and HGB, respectively. These SNPs were all located within the region of 34.6-36.5 Mb on SSC7. Four SNPs within the region of 43.7-47.0 Mb and fifty-five SNPs within the region of 42.2-73.8 Mb on SSC8 showed significant association with MCH and MCV, respectively. At chromosome-wide significant level, one SNP at 29.2 Mb on SSC1 and two SNPs within the region of 26.0-26.2 Mb were found to be significantly associated with RBC and RDW, respectively. Many of the SNPs were located within previously reported QTL regions and appeared to narrow down the regions compared with previously described QTL intervals. In current research, a total of seven significant SNPs were found within six candidate genes SCUBE3, KDR, TDO, IGFBP7, ADAMTS3 and AFP. In addition, the KIT gene, which has been previously reported to relate to hematological parameters, was located within the region significantly associated with MCH and MCV and could be a candidate gene. These results of this study may lead to a better understanding of the molecular mechanisms of hematological parameters in pigs.

Crohn's disease and the mycobacterioses: a quarter century later. Causation or simple association?

It has been more than 25 years since Mycobacterium paratuberculosis was first proposed as an etiologic agent in Crohn's disease based on the isolation of this organism from several patients. Since that time, a great deal of information has been accumulated that clearly establishes an association between M. paratuberculosis and Crohn's disease. However, data are conflicting and difficult to interpret and the field has become divided into committed advocates and confirmed skeptics. This review is an attempt to provide a thorough and objective summary of current knowledge from both basic and clinical research from the views and interpretations of both the antagonists and proponents. The reader is left to draw his or her own conclusions related to the validity of the issues and claims made by the opposing views and data interpretations. Whether M. paratuberculosis is a causative agent in some cases or simply represents an incidental association remains a controversial topic, but current evidence suggests that the notion should not be so readily dismissed. Remaining questions that need to be addressed in defining the role of M. paratuberculosis in Crohn's disease and future implications are discussed.

Genome-wide association analysis of meat quality traits in a porcine Large White × Minzhu intercross population

Pork quality is an economically important trait and one of the main selection criteria for breeding in the swine industry. In this genome-wide association study (GWAS), 455 pigs from a porcine Large White × Minzhu intercross population were genotyped using the Illumina PorcineSNP60K Beadchip, and phenotyped for intramuscular fat content (IMF), marbling, moisture, color L*, color a*, color b* and color score in the longissimus muscle (LM). Association tests between each trait and the SNPs were performed via the Genome Wide Rapid Association using the Mixed Model and Regression-Genomic Control (GRAMMAR-GC) approach. From the Ensembl porcine database, SNP annotation was implemented using Sus scrofa Build 9. A total of 45 SNPs showed significant association with one or multiple meat quality traits. Of the 45 SNPs, 36 were located on SSC12. These significantly associated SNPs aligned to or were in close approximation to previously reported quantitative trait loci (QTL) and some were located within introns of previously reported candidate genes. Two haplotype blocks ASGA0100525-ASGA0055225-ALGA0067099-MARC0004712-DIAS0000861, and ASGA0085522-H3GA0056170 were detected in the significant region. The first block contained the genes MYH1, MYH2 and MYH4. A SNP (ASGA0094812) within an intron of the USP43 gene was significantly associated with five meat quality traits. The present results effectively narrowed down the associated regions compared to previous QTL studies and revealed haplotypes and candidate genes on SSC12 for meat quality traits in pigs.

Quantitative trait loci associated with the immune response to a bovine respiratory syncytial virus vaccine

Infectious disease is an important problem for animal breeders, farmers and governments worldwide. One approach to reducing disease is to breed for resistance. This linkage study used a Charolais-Holstein F2 cattle cross population (n = 501) which was genotyped for 165 microsatellite markers (covering all autosomes) to search for associations with phenotypes for Bovine Respiratory Syncytial Virus (BRSV) specific total-IgG, IgG1 and IgG2 concentrations at several time-points pre- and post-BRSV vaccination. Regions of the bovine genome which influenced the immune response induced by BRSV vaccination were identified, as well as regions associated with the clearance of maternally derived BRSV specific antibodies. Significant positive correlations were detected within traits across time, with negative correlations between the pre- and post-vaccination time points. The whole genome scan identified 27 Quantitative Trait Loci (QTL) on 13 autosomes. Many QTL were associated with the Thymus Helper 1 linked IgG2 response, especially at week 2 following vaccination. However the most significant QTL, which reached 5% genome-wide significance, was on BTA 17 for IgG1, also 2 weeks following vaccination. All animals had declining maternally derived BRSV specific antibodies prior to vaccination and the levels of BRSV specific antibody prior to vaccination were found to be under polygenic control with several QTL detected.

Heifers from the same population (n = 195) were subsequently immunised with a 40-mer Foot-and-Mouth Disease Virus peptide (FMDV) in a previous publication. Several of these QTL associated with the FMDV traits had overlapping peak positions with QTL in the current study, including the QTL on BTA23 which included the bovine Major Histocompatibility Complex (BoLA), and QTL on BTA9 and BTA24, suggesting that the genes underlying these QTL may control responses to multiple antigens. These results lay the groundwork for future investigations to identify the genes underlying the variation in clearance of maternal antibody and response to vaccination.

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.

Control of paratuberculosis in beef cattle

As with any susceptible livestock species, the key to control of paratuberculosis in beef cattle is to reduce exposure of the susceptible calves to Mycobacterium avium subsp. paratuberculosis (MAP)–contaminated feces. Because beef calves remain with mature, potentially shedding cattle until weaning, control strategies are aimed at providing an environment with the least possible fecal burden and removing MAP shedders as soon as possible. Testing and culling or separation may be more important in beef cattle than in dairy cattle. Seedstock owners have greater potential for economic loss from paratuberculosis, making control program more financially attractive to them than to commercial beef cattle producers.

Genetic variation of toll-like receptor genes and infection by Mycobacterium avium ssp. paratuberculosis in Holstein-Friesian cattle

Toll-like receptors (TLR) are membrane proteins that play a key role in innate immunity, by recognizing pathogens and subsequently activating appropriate responses. Mutations in TLR genes are associated with susceptibility to inflammatory and infectious diseases in humans. In cattle, 3 members of the TLR family, TLR1, TLR2, and TLR4, are associated with Mycobacterium avium ssp. paratuberculosis infection, although the extent of this association for the TLR1 and TLR4 receptors has not yet been determined. Moreover, the causal variant in the TLR2 gene has not yet been unequivocally established. In this study, 24 single nucleotide polymorphisms (SNP) in the bovine TLR1, TLR2, and TLR4 genes were selected from the literature, databases, and in silico searches, for a population-based genetic association study of a Spanish Holstein-Friesian sample. Whereas previous results regarding the TLR1 gene were not corroborated, a risk haplotype was detected in TLR2; however, its low frequency indicates that this detected association should be interpreted with caution. In the case of the TLR4 gene, 3 tightly linked SNP were found to be associated with susceptibility to M. avium ssp. paratuberculosis infection. Moreover, one of these SNP, the SNP c.-226G>C, which is localized in the 5'UTR region of the TLR4 gene, has been reported to be able to alter TLR4 expression, raising the possibility that this mutation may contribute to the response of the individual to infection.

Characterization and management of asymptomatic Mycobacterium infections at the Zebrafish International Resource Center

The Zebrafish International Resource Center (ZIRC) supplies wildtype, mutant, and transgenic zebrafish (Danio rerio) to the international research community. In 2005, the ZIRC halted shipment of adult Tübingen (TU) zebrafish, a popular wildtype line, after diagnosis of asymptomatic Mycobacterium chelonae infections in a high proportion of the TU stock. Mycobacterium presents a zoonotic risk to fish handlers. In addition, the presence of underlying chronic disease in a model organism is unacceptable. The TU stock was depopulated and replaced by a new import of TU with the intent of reducing disease prevalence. In the current study, we sampled the new population of TU and fish of the AB, Tupfel long-fin (TL), TAB5 and TAB14 (2 AB × TU hybrid lines), and wildtype-in-Kalkutta (WIK) lines for histologic evaluation and acid-fast staining and compared the prevalence of subclinical mycobacteriosis between these lines. Although prevalence in the new TU stock was lower than that of the original TU stock, asymptomatic infections with Mycobacterium remained high (10%) in the new TU stock held in 20-gal tanks. The prevalence was similar (10%) in the TAB5 line compared with other wildtype lines held in similar conditions. Prevalence of infections in TU can be minimized by husbandry adjustments, including tank size, population density, and cleaning method. Application of these findings has allowed us to decrease mycobacteriosis in TU zebrafish and resume shipment of TU adults to the research community.

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.