Whole-Genome association analysis of susceptibility to paratuberculosis in Holstein cattle

Genetic parameters and trends for Johne's disease in US 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. The data were extracted from a JD data repository, maintained at the Council on Dairy Cattle Breeding, and initially supplied by 2 dairy record 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 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 EBVs 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. The correlations of bulls' PTA with economically important traits such as milk yield, milk protein, milk fat, somatic cell score, and mastitis were low, indicating a nonoverlapping 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.

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

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

Early-stage findings in an experimental calf model infected with Argentinean isolates of Mycobacterium avium subsp. paratuberculosis

Mycobacterium avium subsp. paratuberculosis (MAP) is an important pathogen that causes granulomatous enteritis known as Johne's disease or paratuberculosis (PTB). In this study an experimental model of calves infected with Argentinean isolates of MAP for 180 days was used to provide more data of the early PTB stages. Calves were challenged by oral route with MAP strain _IS900-RFLPA (MA; n = 3), MAP strain _IS900-RFLPC (MC; n = 2) or mock infected (MI; n = 2), and response to infection was evaluated through peripheral cytokine expression, MAP tissue distribution and histopathological early-stage findings. Specific and varied levels of IFN-γ were only detected at 80 days post-infection in infected calves. These data indicate that specific IFN-γ is not a useful indicator for early detection of MAP infection in our calf model. At 110 days post-infection, TNF-α expression was higher than IL-10 in 4 of the 5 infected animals and a significant decrease of TNF-α expression was detected in infected vs. non-infected calves. All calves challenged were identified as infected by mesenteric lymph node tissue culture and real time IS900 PCR. In addition, for lymph nodes samples, the agreement between these techniques was almost perfect (κ = 0.86). Colonization of tissues and levels of tissue infection varied between individuals. Evidence of early MAP dissemination to extraintestinal tissues such as the liver was detected by culture in one animal (MAP strain _IS900-RFLPA). In both groups microgranulomatous lesions were observed predominantly in the lymph nodes, with giant cells present only in the MA group. In summary, the findings described herein may indicate that local MAP strains induced specific immune responses with particularities that could suggest differences in their biological behavior. Further studies should be carried out in order to obtain an in-depth understanding of the influence of MAP strains in host-pathogen interactions and the outcome of disease.

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

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

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

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

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

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

Johne's Disease in Dairy Cattle: An Immunogenetic Perspective

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

Identification of loci associated with susceptibility to Mycobacterium avium subsp. paratuberculosis infection in Holstein cattle using combinations of diagnostic tests and imputed whole-genome sequence data

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

The Paratuberculosis Paradigm Examined: A Review of Host Genetic Resistance and Innate Immune Fitness in Mycobacterium avium subsp. Paratuberculosis Infection

Paratuberculosis, or Johne's Disease (JD) is a debilitating chronic enteritis mainly affecting ruminants caused by Mycobacterium avium subsp. paratuberculosis (MAP). This organism causes worldwide economic losses to the livestock industry, and is of public health importance due to the potential zoonotic risk between MAP and Crohn's disease (CD) in humans. Without economical treatments, or a vaccine capable of preventing infection without causing cross-reactions with bovine tuberculosis, test-and-cull methods for disease control are imperative. Unfortunately, difficulties in diagnostics and long subclinical stage hinder adequate control and is further complicated by variation in MAP exposure outcome. Interestingly, the majority of infections result in asymptomatic presentation and never progress to clinical disease. One contributing factor is host genetics, where polymorphisms in innate immune genes have been found to influence resistance and susceptibility to disease. Candidate genes identified across studies overlap with those found in CD and tuberculosis including; Solute carrier family 11 member 1 gene (SLC11A1), Nucleotide-binding-oligomerization domain containing gene 2 (NOD2), Major histocompatibility complex type II (MHC-II), and Toll-like receptor (TLR) genes. This review will highlight evidence supporting the vital role of these genes in MAP infection outcome, associated challenges, and implications for the future of JD research.

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

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

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

Background

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

Results

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

Conclusion

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

Supplementary Information

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

Survey of Candidate Single-Nucleotide Polymorphisms in SLC11A1, TLR4, NOD2, PGLYRP1, and IFNγ in Ankole Longhorn Cattle in Central Region of Uganda to Determine Their Role in Mycobacterium avium Subspecies paratuberculosis Infection Outcome

Mycobacterium avium ssp. paratuberculosis (MAP) is the cause of Johne's disease (JD) in a wide range of domestic and wild ruminants. Single-nucleotide polymorphisms (SNPs) in several genes including solute-like carrier 11A1 (SLC11A1), interferon gamma (IFNγ), Toll-like receptor 4 (TLR4), nucleotide-binding oligomerization domain 2 gene (NOD2), and bovine peptidoglycan recognition protein 1 (PGLYRP1) have been implicated in influencing the infection outcome of MAP in cattle. We have carried out a survey in a population of Ankole cattle from three districts in the central region of Uganda including Isingiro, Lyantonde, and Rakai to determine the role played by several SNPs on the above genes in the infection outcome of local cattle in Uganda. Nine hundred fifty-five heads of cattle obtained from 93 herds were tested using ELISA. Thirty-five ELISA-positive cattle and 35 negative herd mates from a total of 955 cattle tested for MAP were genotyped using iPLEX MassARRAY genotyping systems to detect the presence of a total of 13 SNPS in five different genes (SLC11A1, IFNγ, TLR4, NOD2, and PGLYRP1). The cow-level prevalence of MAP infection in Ankole Longhorn cattle in the three districts was 3.98% (35/955), while the herd-level prevalence was 27.9% and within-herd prevalence was 12 ± 1.5% (95% CI = 9.1-14.8%). The genotypes and allele frequencies of the MAP-positive cattle were compared with those of their ELISA-negative herd mates to determine the significance of the polymorphisms. The results showed that SNPs rs109915208, rs110514940, and rs110905610 on SLC11A1, c.480G>A and c.625C>A on PGLYRP1, and c.2021C>T on TLR4 were monomorphic in both seropositive and seronegative cattle and therefore had no influence on the infection outcome. The remaining SNPs studied in the five genes [SLC11A1: rs109614179; TLR4: rs29017188 (c.226G>C), c.2021C>T; NOD2: rs110536091, rs111009394; PGLYRP1: c.102G>C, c.480G>A, c.625C>A; IFNγ: rs110853455] were polymorphic, but their allele and genotype frequencies did not show any significant difference between the seropositive and seronegative cattle. No significant difference was observed for any haplotype at the gene level.

Functional analysis of bovine interleukin-10 receptor alpha in response to Mycobacterium avium subsp. paratuberculosis lysate using CRISPR/Cas9

Background

The interleukin-10 receptor alpha (IL10RA) gene codes for the alpha chain of the IL-10 receptor which binds the cytokine IL-10. IL-10 is an anti-inflammatory cytokine with immunoregulatory function during the pathogenesis of many inflammatory disorders in livestock, including Johne’s disease (JD). JD is a chronic enteritis in cattle caused by Mycobacterium avium subsp. paratuberculosis (MAP) and is responsible for significant economic losses to the dairy industry. Several candidate genes including IL10RA have been found to be associated with JD. The aim of this study was to better understand the functional significance of IL10RA in the context of immune stimulation with MAP cell wall lysate.

Results

An IL10RA knock out (KO) bovine mammary epithelial cell (MAC-T) line was generated using the CRISPR/cas9 (Clustered Regularly Interspaced Short Palindromic Repeats/CRISPR-associated protein 9) gene editing system. These IL10RA KO cells were stimulated with the immune stimulant MAP lysate +/− IL-10, or with LPS as a positive control. In comparison to unedited cells, relative quantification of immune-related genes after stimulation revealed that knocking out IL10RA resulted in upregulation of pro-inflammatory cytokine gene expression (TNFA, IL1A, IL1B and IL6) and downregulation of suppressor of cytokine signaling 3 (SOCS3), a negative regulator of pro-inflammatory cytokine signaling. At the protein level knocking out IL10RA also resulted in upregulation of inflammatory cytokines - TNF-α and IL-6 and chemokines - IL-8, CCL2 and CCL4, relative to unedited cells.

Conclusions

The findings of this study illustrate the broad and significant effects of knocking out the IL10RA gene in enhancing pro-inflammatory cytokine expression and further support the immunoregulatory role of IL10RA in eliciting an anti-inflammatory response as well as its potential functional involvement during the immune response associated with JD.

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

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

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

Background

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

Results

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

Conclusions

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

Divergent Antigen-Specific Cellular Immune Responses during Asymptomatic Subclinical and Clinical States of Disease in Cows Naturally Infected with Mycobacterium avium subsp. paratuberculosis

Infection of the host with Mycobacterium avium subsp. paratuberculosis results in chronic and progressive enteritis that traverses both subclinical and clinical stages. The mechanism(s) for the shift from an asymptomatic subclinical disease state to advanced clinical disease is not fully understood. In the present study, naturally infected dairy cattle were divided into subclinical and clinical infection groups, along with noninfected control cows of similar parity, to study host immune responses in different stages of infection. Both infection groups had higher levels of secretion of gamma interferon (IFN-γ), tumor necrosis factor alpha (TNF-α), and interleukin-2 (IL-2) than control cows, whereas only clinical cows had increased secretion of IL-10, IL-12, and IL-18 upon stimulation of peripheral blood mononuclear cells (PBMCs) with antigen. Conversely, secretion of IL-17Α was decreased for clinical cows compared to subclinical and control cows. Proinflammatory cytokine genes were upregulated only for subclinical cows, whereas increased IL-10 and IL-17 gene expression levels were observed for both infection groups. Increased CD4⁺, CD8⁺, and γδ T cell receptor-positive (TCR⁺) T cells were observed for subclinical cows compared to clinical cows. Although clinical cows expressed antigen-specific immune responses, the profile for subclinical cows was one of a dominant proinflammatory response to infection. We reason that a complex coordination of immune responses occurs during M. avium subsp. paratuberculosis infection, with these responses shifting as the host transitions through the different stages of infection and disease (subclinical to clinical). A further understanding of the series of events characterized by Th1/Th2/Th17 responses will provide mechanisms for disease progression and may direct insightful intervention strategies.

Variance components for bovine tuberculosis infection and multi-breed genome-wide association analysis using imputed whole genome sequence data

Bovine tuberculosis (bTB) is an infectious disease of cattle generally caused by Mycobacterium bovis, a bacterium that can elicit disease humans. Since the 1950s, the objective of the national bTB eradication program in Republic of Ireland was the biological extinction of bTB; that purpose has yet to be achieved. Objectives of the present study were to develop the statistical methodology and variance components to undertake routine genetic evaluations for resistance to bTB; also of interest was the detection of regions of the bovine genome putatively associated with bTB infection in dairy and beef breeds. The novelty of the present study, in terms of research on bTB infection, was the use of beef breeds in the genome-wide association and the utilization of imputed whole genome sequence data. Phenotypic bTB data on 781,270 animals together with imputed whole genome sequence data on 7,346 of these animals' sires were available. Linear mixed models were used to quantify variance components for bTB and EBVs were validated. Within-breed and multi-breed genome-wide associations were undertaken using a single-SNP regression approach. The estimated genetic standard deviation (0.09), heritability (0.12), and repeatability (0.30) substantiate that genetic selection help to eradicate bTB. The multi-breed genome-wide association analysis identified 38 SNPs and 64 QTL regions associated with bTB infection; two QTL regions (both on BTA23) identified in the multi-breed analysis overlapped with the within-breed analyses of Charolais, Limousin, and Holstein-Friesian. Results from the association analysis, coupled with previous studies, suggest bTB is controlled by an infinitely large number of loci, each having a small effect. The methodology and results from the present study will be used to develop national genetic evaluations for bTB in the Republic of Ireland. In addition, results can also be used to help uncover the biological architecture underlying resistance to bTB infection in cattle.

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.

Association between combinations of genetic polymorphisms and epidemiopathogenic forms of bovine paratuberculosis

Control of major mycobacterial diseases affecting livestock is a challenging issue that requires different approaches. The use of genetic markers for improving resistance to Mycobacterium avium subsp. paratuberculosis infection in cattle has been explored as a promising population strategy We performed paratuberculosis epidemiopathogenic phenotypic and genotypic characterization involving 24 SNPs in six candidate genes (NOD2, CD209, SLC11A1, SP110, TLR2 and TLR4) on 502 slaughtered Friesian cows. In the current study, we investigate whether recently proposed paratuberculosis (PTB) epidemiopathogenic (EP) forms (apparently free-AF, latent-LAT and patent-PAT) could be associated with some combination of these 24 SNPs. Best EP form grouping was obtained using a combination of 5 SNPs in four genes (CD209: rs210748127; SLC11A1: rs110090506; SP110: rs136859213 and rs110480812; and TLR2: rs41830058). These groups were defined according to the level of infection progression risk to patent epidemiopathogenic forms and showed the following distributions: LOWIN (low) with 39 (8%) cases (94.9% AF/5.1% LAT/0% PAT); LATIN (low) with 17 (3%) cases (5.9% AF/94.1% LAT/0% PAT); AVERIN (average) with 413 (82%) cases (52.1% AF/38.5% LAT/9.4% PAT) and PATIN (patent) with 33 (7%) cases (36.4% AF/24.2% LAT/39.4% PAT). Age of slaughter was significantly higher for LATIN (88.3 months) compared to AVERIN (65.3 months; p = 0.0007) and PATIN (59.1 months; p = 0.0004), and for LOWIN (73.9 months) compared to PATIN (p = 0.0233), and nearly significant compared to AVERIN (p = 0.0572) These results suggest that some selected genetic polymorphisms have a potential use as markers of PTB EP forms and thus add a new tool for the control of this widespread infection.

Case definition terminology for paratuberculosis (Johne's disease)

Paratuberculosis (Johne's disease) is an economically significant condition caused by Mycobacterium avium subsp. paratuberculosis. However, difficulties in diagnosis and classification of individual animals with the condition have hampered research and impeded efforts to halt its progressive spread in the global livestock industry. Descriptive terms applied to individual animals and herds such as exposed, infected, diseased, clinical, sub-clinical, infectious and resistant need to be defined so that they can be incorporated consistently into well-understood and reproducible case definitions. These allow for consistent classification of individuals in a population for the purposes of analysis based on accurate counts. The outputs might include the incidence of cases, frequency distributions of the number of cases by age class or more sophisticated analyses involving statistical comparisons of immune responses in vaccine development studies, or gene frequencies or expression data from cases and controls in genomic investigations. It is necessary to have agreed definitions in order to be able to make valid comparisons and meta-analyses of experiments conducted over time by a given researcher, in different laboratories, by different researchers, and in different countries. In this paper, terms are applied systematically in an hierarchical flow chart to enable classification of individual animals. We propose descriptive terms for different stages in the pathogenesis of paratuberculosis to enable their use in different types of studies and to enable an independent assessment of the extent to which accepted definitions for stages of disease have been applied consistently in any given study. This will assist in the general interpretation of data between studies, and will facilitate future meta-analyses.

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.

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.

An across-breed genome wide association analysis of susceptibility to paratuberculosis in dairy cattle

Paratuberculosis is a chronic disease of ruminants caused by Mycobacterium avium subspecies paratuberculosis (MAP). It occurs worldwide and causes a significant loss in the animal production industry. There is no cure for MAP infection and vaccination is problematic. Identification of genetics of susceptibility could be a useful adjunct for programs that focus on management, testing and culling of diseased animals. A case-control, genome-wide association study (GWAS) was conducted using Holstein and Jersey cattle in a combined analysis in order to identify markers and chromosomal regions associated with susceptibility to MAP infection across-breed. A mixed-model method (GRAMMAR-GC) implemented in the GenABEL R package and a Bayes C analysis implemented in GenSel software were used as alternative approaches to conduct GWAS analysis focused on single SNPs and chromosomal segments, respectively. After conducting quality control, 22 406 SNPs from 2157 individuals were available for the GRAMMAR-GC (Bayes C) analysis and 45 640 SNPs from 2199 individuals were available for the Bayes C analysis. One SNP located on BTA27 (8·6 Mb) was identified as moderately associated (P < 5 × 10-5, FDR = 0·44) in the GRAMMAR-GC analysis of the combined breed data. Nine 1 Mb windows located on BTA 2, 3 (3 windows), 6, 8, 25, 27 and 29 each explained ≥1% of the total proportion of genetic variance in the Bayes C analysis. In an analysis ignoring differences in linkage phase, two moderately significantly associated SNPs were identified; ARS-BFGL-NGS-19381 on BTA23 (32 Mb) and Hapmap40994-BTA-46361 on BTA19 (61 Mb). New common genomic regions and candidate genes have been identified from the across-breed analysis that might be involved in the immune response and susceptibility to MAP infection.

Whole genome sequencing of Guzerá cattle reveals genetic variants in candidate genes for production, disease resistance, and heat tolerance

In bovines, artificial selection has produced a large number of breeds which differ in production, environmental adaptation, and health characteristics. To investigate the genetic basis of these phenotypical differences, several bovine breeds have been sequenced. Millions of new SNVs were described at every new breed sequenced, suggesting that every breed should be sequenced. Guzerat or Guzerá is an indicine breed resistant to drought and parasites that has been the base for some important breeds such as Brahman. Here, we describe the sequence of the Guzerá genome and the in silico functional analyses of intragenic breed-specific variations. Mate-paired libraries were generated using the ABI SOLiD system. Sequences were mapped to the Bos taurus reference genome (UMD 3.1) and 87% of the reference genome was covered at a 26X. Among the variants identified, 2,676,067 SNVs and 463,158 INDELs were homozygous, not found in any database searched, and may represent true differences between Guzerá and B. taurus. Functional analyses investigated with the NGS-SNP package focused on 1069 new, non-synonymous SNVs, splice-site variants (including acceptor and donor sites, and the conserved regions at both intron borders, referred to here as splice regions) and coding INDELs (NS/SS/I). These NS/SS/I map to 935 genes belonging to cell communication, environmental adaptation, signal transduction, sensory, and immune systems pathways. These pathways have been involved in phenotypes related to health, adaptation to the environment and behavior, and particularly, disease resistance and heat tolerance. Indeed, 105 of these genes are known QTLs for milk, meat and carcass, production, reproduction, and health traits. Therefore, in addition to describing new genetic variants, our approach provided groundwork for unraveling key candidate genes and mutations.

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.

Genetic analysis of infectious diseases: estimating gene effects for susceptibility and infectivity

Background

Genetic selection of livestock against infectious diseases can complement existing interventions to control infectious diseases. Most genetic approaches that aim at reducing disease prevalence assume that individual disease status (infected/not-infected) is solely a function of its susceptibility to a particular pathogen. However, individual infectivity also affects the risk and prevalence of an infection in a population. Variation in susceptibility and infectivity between hosts affects transmission of an infection in the population, which is usually measured by the value of the basic reproduction ratio R₀. R₀ is an important epidemiological parameter that determines the risk and prevalence of infectious diseases. An individual’s breeding value for R₀ is a function of its genes that influence both susceptibility and infectivity. Thus, to estimate the effects of genes on R₀, we need to estimate the effects of genes on individual susceptibility and infectivity. To that end, we developed a generalized linear model (GLM) to estimate relative effects of genes for susceptibility and infectivity. A simulation was performed to investigate bias and precision of the estimates, the effect of R₀, the size of the effects of genes for susceptibility and infectivity, and relatedness among group mates on bias and precision. We considered two bi-allelic loci that affect, respectively, the individuals’ susceptibility only and individuals’ infectivity only.

Results

A GLM with complementary log–log link function can be used to estimate the relative effects of genes on the individual’s susceptibility and infectivity. The model was developed from an equation that describes the probability of an individual to become infected as a function of its own susceptibility genotype and infectivity genotypes of all its infected group mates. Results show that bias is smaller when R₀ ranges approximately from 1.8 to 3.1 and relatedness among group mates is higher. With larger effects, both absolute and relative standard deviations become clearly smaller, but the relative bias remains the same.

Conclusions

We developed a GLM to estimate the relative effect of genes that affect individual susceptibility and infectivity. This model can be used in genome-wide association studies that aim at identifying genes that influence the prevalence of infectious diseases.

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.

Genomics and disease resistance studies in livestock

This paper considers the application of genetic and genomic techniques to disease resistance, the interpretation of data arising from such studies and the utilisation of the research outcomes to breed animals for enhanced resistance. Resistance and tolerance are defined and contrasted, factors affecting the analysis and interpretation of field data presented, and appropriate experimental designs discussed. These general principles are then applied to two detailed case studies, infectious pancreatic necrosis in Atlantic salmon and bovine tuberculosis in dairy cattle, and the lessons learnt are considered in detail. It is concluded that the rate limiting step in disease genetic studies will generally be provision of adequate phenotypic data, and its interpretation, rather than the genomic resources. Lastly, the importance of cross-disciplinary dialogue between the animal health and animal genetics communities is stressed.

Genome wide association study identifies 20 novel promising genes associated with milk fatty acid traits in Chinese Holstein

Detecting genes associated with milk fat composition could provide valuable insights into the complex genetic networks of genes underling variation in fatty acids synthesis and point towards opportunities for changing milk fat composition via selective breeding. In this study, we conducted a genome-wide association study (GWAS) for 22 milk fatty acids in 784 Chinese Holstein cows with the PLINK software. Genotypes were obtained with the Illumina BovineSNP50 Bead chip and a total of 40,604 informative, high-quality single nucleotide polymorphisms (SNPs) were used. Totally, 83 genome-wide significant SNPs and 314 suggestive significant SNPs associated with 18 milk fatty acid traits were detected. Chromosome regions that affect milk fatty acid traits were mainly observed on BTA1, 2, 5, 6, 7, 9, 13, 14, 18, 19, 20, 21, 23, 26 and 27. Of these, 146 SNPs were associated with more than one milk fatty acid trait; most of studied fatty acid traits were significant associated with multiple SNPs, especially C18:0 (105 SNPs), C18 index (93 SNPs), and C14 index (84 SNPs); Several SNPs are close to or within the DGAT1, SCD1 and FASN genes which are well-known to affect milk composition traits of dairy cattle. Combined with the previously reported QTL regions and the biological functions of the genes, 20 novel promising candidates for C10:0, C12:0, C14:0, C14:1, C14 index, C18:0, C18:1n9c, C18 index, SFA, UFA and SFA/UFA were found, which composed of HTR1B, CPM, PRKG1, MINPP1, LIPJ, LIPK, EHHADH, MOGAT1, ECHS1, STAT1, SORBS1, NFKB2, AGPAT3, CHUK, OSBPL8, PRLR, IGF1R, ACSL3, GHR and OXCT1. Our findings provide a groundwork for unraveling the key genes and causal mutations affecting milk fatty acid traits in dairy cattle.

Genetic association analysis of paratuberculosis forms in holstein-friesian cattle

A genetic susceptibility to Mycobacterium avium subsp. paratuberculosis (MAP) infections in ruminants has been longtime suspected to exist. Recently, natural infections in cattle have been reclassified into latent and patent forms based on histopathological findings and their associations with immunological and microbiological variables. This study aims to explore whether these newly defined phenotypes are associated with twenty-four single-nucleotide polymorphisms (SNPs) in six bovine candidate genes: nucleotide-binding oligomerization domain 2 (NOD2), solute carrier family 11 member A1 (SLC11A1), nuclear body protein SP110 (SP110), toll-like receptors (TLRs) 2 and 4, and CD209 (also known as DC-SIGN, dendritic cell-specific ICAM3-grabbing nonintegrin). SNPs were genotyped for 772 Holstein-Friesian animals (52.6% apparently free; 38.1% latent; 9.3% patent) by TaqMan OpenArray technology. Genotypic-phenotypic associations were assessed by logistic regression analysis adjusted for age at slaughter, under five models (codominant, dominant, recessive, overdominant, and log-additive), and corrected for multiple testing. The rs208222804 C allele (CD209 gene) was found to be associated with latent paratuberculosis (log-additive model: P < 0.0034 after permutation procedure; OR = 0.64, 95% CI = 0.48–0.86). No significant association was detected between any SNP and the patent phenotype. Consequently, CD209 gene may play a key role in the pathogenesis of bovine paratuberculosis.

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.

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).

Genome wide association studies for body conformation traits in the Chinese Holstein cattle population

BACKGROUND

Genome-wide association study (GWAS) is a powerful tool for revealing the genetic basis of quantitative traits. However, studies using GWAS for conformation traits of cattle is comparatively less. This study aims to use GWAS to find the candidates genes for body conformation traits.

RESULTS

The Illumina BovineSNP50 BeadChip was used to identify single nucleotide polymorphisms (SNPs) that are associated with body conformation traits. A least absolute shrinkage and selection operator (LASSO) was applied to detect multiple SNPs simultaneously for 29 body conformation traits with 1,314 Chinese Holstein cattle and 52,166 SNPs. Totally, 59 genome-wide significant SNPs associated with 26 conformation traits were detected by genome-wide association analysis; five SNPs were within previously reported QTL regions (Animal Quantitative Trait Loci (QTL) database) and 11 were very close to the reported SNPs. Twenty-two SNPs were located within annotated gene regions, while the remainder were 0.6-826 kb away from known genes. Some of the genes had clear biological functions related to conformation traits. By combining information about the previously reported QTL regions and the biological functions of the genes, we identified DARC, GAS1, MTPN, HTR2A, ZNF521, PDIA6, and TMEM130 as the most promising candidate genes for capacity and body depth, chest width, foot angle, angularity, rear leg side view, teat length, and animal size traits, respectively. We also found four SNPs that affected four pairs of traits, and the genetic correlation between each pair of traits ranged from 0.35 to 0.86, suggesting that these SNPs may have a pleiotropic effect on each pair of traits.

CONCLUSIONS

A total of 59 significant SNPs associated with 26 conformation traits were identified in the Chinese Holstein population. Six promising candidate genes were suggested, and four SNPs showed genetic correlation for four pairs of traits.

Evaluation of age-dependent susceptibility in calves infected with two doses of Mycobacterium avium subspecies paratuberculosis using pathology and tissue culture

The longstanding assumption that calves of more than 6 months of age are more resistant to Mycobacterium avium subspecies paratuberculosis (MAP) infection has recently been challenged. In order to elucidate this, a challenge experiment was performed to evaluate age- and dose-dependent susceptibility to MAP infection in dairy calves. Fifty-six calves from MAP-negative dams were randomly allocated to 10 MAP challenge groups (5 animals per group) and a negative control group (6 calves). Calves were inoculated orally on 2 consecutive days at 5 ages: 2 weeks and 3, 6, 9 or 12 months. Within each age group 5 calves received either a high - or low - dose of 5 × 10⁹ CFU or 5 × 10⁷ CFU, respectively. All calves were euthanized at 17 months of age. Macroscopic and histological lesions were assessed and bacterial culture was done on numerous tissue samples. Within all 5 age groups, calves were successfully infected with either dose of MAP. Calves inoculated at < 6 months usually had more culture-positive tissue locations and higher histological lesion scores. Furthermore, those infected with a high dose had more severe scores for histologic and macroscopic lesions as well as more culture-positive tissue locations compared to calves infected with a low dose. In conclusion, calves to 1 year of age were susceptible to MAP infection and a high infection dose produced more severe lesions than a low dose.

Combining cow and bull reference populations to increase accuracy of genomic prediction and genome-wide association studies

Genomic selection holds the promise to be particularly beneficial for traits that are difficult or expensive to measure, such that access to phenotypes on large daughter groups of bulls is limited. Instead, cow reference populations can be generated, potentially supplemented with existing information from the same or (highly) correlated traits available on bull reference populations. The objective of this study, therefore, was to develop a model to perform genomic predictions and genome-wide association studies based on a combined cow and bull reference data set, with the accuracy of the phenotypes differing between the cow and bull genomic selection reference populations. The developed bivariate Bayesian stochastic search variable selection model allowed for an unbalanced design by imputing residuals in the residual updating scheme for all missing records. The performance of this model is demonstrated on a real data example, where the analyzed trait, being milk fat or protein yield, was either measured only on a cow or a bull reference population, or recorded on both. Our results were that the developed bivariate Bayesian stochastic search variable selection model was able to analyze 2 traits, even though animals had measurements on only 1 of 2 traits. The Bayesian stochastic search variable selection model yielded consistently higher accuracy for fat yield compared with a model without variable selection, both for the univariate and bivariate analyses, whereas the accuracy of both models was very similar for protein yield. The bivariate model identified several additional quantitative trait loci peaks compared with the single-trait models on either trait. In addition, the bivariate models showed a marginal increase in accuracy of genomic predictions for the cow traits (0.01-0.05), although a greater increase in accuracy is expected as the size of the bull population increases. Our results emphasize that the chosen value of priors in Bayesian genomic prediction models are especially important in small data sets.

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.

[Advances on genome-wide association study for economically important traits in dairy cattle]

Milk performance traits have been considered as the most important production traits. With the development of balance breeding sense, much attention has been paid to the functional traits such as reproduction, type, health, and longevity traits, which is put into breeding programs. Identification of major genes or genetic markers for milk production and functional traits and their applications in breeding programs of dairy cattle are expected to improve the genetic progress of these traits. Occurrence of millions of SNPs in whole genome and high throughout genotyping techniques has made genome-wide association study to be an important strategy to identify genes responsible for economic traits in domestic animals. This paper reviewed the advances on genome-wide association study for milk production traits and functional traits in dairy cattle.

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.

Predicted accuracy of and response to genomic selection for new traits in dairy cattle

Genomic selection relaxes the requirement of traditional selection tools to have phenotypic measurements on close relatives of all selection candidates. This opens up possibilities to select for traits that are difficult or expensive to measure. The objectives of this paper were to predict accuracy of and response to genomic selection for a new trait, considering that only a cow reference population of moderate size was available for the new trait, and that selection simultaneously targeted an index and this new trait. Accuracy for and response to selection were deterministically evaluated for three different breeding goals. Single trait selection for the new trait based only on a limited cow reference population of up to 10 000 cows, showed that maximum genetic responses of 0.20 and 0.28 genetic standard deviation (s.d.) per year can be achieved for traits with a heritability of 0.05 and 0.30, respectively. Adding information from the index based on a reference population of 5000 bulls, and assuming a genetic correlation of 0.5, increased genetic response for both heritability levels by up to 0.14 genetic s.d. per year. The scenario with simultaneous selection for the new trait and the index, yielded a substantially lower response for the new trait, especially when the genetic correlation with the index was negative. Despite the lower response for the index, whenever the new trait had considerable economic value, including the cow reference population considerably improved the genetic response for the new trait. For scenarios with a zero or negative genetic correlation with the index and equal economic value for the index and the new trait, a reference population of 2000 cows increased genetic response for the new trait with at least 0.10 and 0.20 genetic s.d. per year, for heritability levels of 0.05 and 0.30, respectively. We conclude that for new traits with a very small or positive genetic correlation with the index, and a high positive economic value, considerable genetic response can already be achieved based on a cow reference population with only 2000 records, even when the reliability of individual genomic breeding values is much lower than currently accepted in dairy cattle breeding programs. New traits may generally have a negative genetic correlation with the index and a small positive economic value. For such new traits, cow reference populations of at least 10 000 cows may be required to achieve acceptable levels of genetic response for the new trait and for the whole breeding goal.

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.