Identification of loci associated with tolerance to Johne's disease in Holstein cattle

Genome-wide association study identifies genomic regions associated with key reproductive traits in Korean Hanwoo cows

BACKGROUND

Conducting genome-wide association studies (GWAS) for reproductive traits in Hanwoo cattle, including age at first calving (AFC), calving interval (CI), gestation length (GL), and number of artificial inseminations per conception (NAIPC), is of paramount significance. These analyses provided a thorough exploration of the genetic basis of these traits, facilitating the identification of key markers for targeted trait improvement. Breeders can optimize their selection strategies, leading to more efficient and sustainable breeding programs, by incorporating genetic insights. This impact extends beyond individual traits and contributes to the overall productivity and profitability of the Hanwoo beef cattle industry. Ultimately, GWAS is essential in ensuring the long-term genetic resilience and adaptability of Hanwoo cattle populations. The primary goal of this study was to identify significant single nucleotide polymorphisms (SNPs) or quantitative trait loci (QTLs) associated with the studied reproductive traits and subsequently map the underlying genes that hold promise for trait improvement.

RESULTS

A genome-wide association study of reproductive traits identified 68 significant single nucleotide polymorphisms (SNPs) distributed across 29 Bos taurus autosomes (BTA). Among them, BTA14 exhibited the highest number of identified SNPs (25), whereas BTA6, BTA7, BTA8, BTA10, BTA13, BTA17, and BTA20 exhibited 8, 5, 5, 3, 8, 2, and 12 significant SNPs, respectively. Annotation of candidate genes within a 500 kb region surrounding the significant SNPs led to the identification of ten candidate genes relevant to age at first calving. These genes were: FANCG, UNC13B, TESK1, TLN1, and CREB3 on BTA8; FAM110B, UBXN2B, SDCBP, and TOX on BTA14; and MAP3K1 on BTA20. Additionally, APBA3, TCF12, and ZFR2, located on BTA7 and BTA10, were associated with the calving interval; PAX1, SGCD, and HAND1, located on BTA7 and BTA13, were linked to gestation length; and RBM47, UBE2K, and GPX8, located on BTA6 and BTA20, were linked to the number of artificial inseminations per conception in Hanwoo cows.

CONCLUSIONS

The findings of this study enhance our knowledge of the genetic factors that influence reproductive traits in Hanwoo cattle populations and provide a foundation for future breeding strategies focused on improving desirable traits in beef cattle. This research offers new evidence and insights into the genetic variants and genome regions associated with reproductive traits and contributes valuable information to guide future efforts in cattle breeding.

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.

Prevalence of paratuberculosis in cattle based on gross and microscopic lesions in Ethiopia

BACKGROUND

Paratuberculosis, caused by Mycobacterium avium subsp. paratuberculosis (MAP), is a chronic progressive granulomatous enteritis mainly affecting domestic and wild ruminants worldwide. Although paratuberculosis could be prevail in Ethiopia, there is a scarcity of epidemiological data on paratuberculosis in the country. Thus, this study was conducted to estimate the prevalence of paratuberculosis based on gross and microscopic lesions in cattle slaughtered at ELFORA Abattoir, central Ethiopia. Small intestines and associated lymph nodes of 400 apparently healthy cattle which were slaughtered at ELFORA export abattoir were examined for gross and microscopic lesions of paratuberculosis. The microscopic lesions were classified into four grades (I-IV) based on the type and number of cells infiltrated into the lesion. The prevalence of paratuberculosis was estimated on the basis of gross as well as microscopic lesion of paratuberculosis.

RESULTS

The prevalence of paratuberculosis was 11.25% (95% Confidence interval, CI = 0.083-0.148) on the basis of gross lesion. However, relatively lower prevalence (2.0%, 95% CI = 0.01, 0.039) was recorded based on microscopic lesion. The gross lesions were characterized by intestinal thickening, mucosal corrugations and enlargement of associated mesenteric lymph nodes. On the other hand, the microscopic lesions were characterized by granuloma of different grades ranging from grade I to grade III lesions.

CONCLUSIONS

The present study indicated the occurrence of paratuberculosis in cattle of Ethiopia based on the detection of gross and microscopic lesions consistent with the lesion of paratuberculosis. The result of this study could be used as baseline information for future studies on the epidemiology and economic significance of paratuberculosis.

Genome-Wide Association Studies for Body Conformation Traits in Korean Holstein Population

The objective of this study was to identify quantitative trait loci (QTL) and nearby candidate genes that influence body conformation traits. Phenotypic data for 24 body conformation traits were collected from a population of 2329 Korean Holstein cattle, and all animals were genotyped using the 50 K Illumina bovine SNP chip. A total of 24 genome-wide significant SNPs associated with 24 body conformation traits were identified by genome-wide association analysis. The selection of the most promising candidate genes was based on gene ontology (GO) terms and the previously identified functions that influence various body conformation traits as determined in our study. These genes include KCNA1, RYBP, PTH1R, TMIE, and GNAI3 for body traits; ANGPT1 for rump traits; MALRD1, INHBA, and HOXA13 for feet and leg traits; and CDK1, RHOBTB1, and SLC17A1 for udder traits, respectively. These findings contribute to our understanding of the genetic basis of body conformation traits in this population and pave the way for future breeding strategies aimed at enhancing desirable traits in dairy cattle.

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.

New insights into the genetic resistance to paratuberculosis in Holstein cattle via single-step genomic evaluation

Background

Bovine paratuberculosis, or Johne’s disease (JD), is a contagious and incurable disease caused by Mycobacterium avium subsp. paratuberculosis (MAP). It has adverse effects on animal welfare and is very difficult to control, leading to serious economic consequences. An important line of defense to this disease is host genetic resistance to MAP, which, when it will be more fully understood, could be improved through selective breeding. Using a large dataset of Holstein cows (161,253 animals including 56,766 cows with ELISA serological phenotypes and 12,431 animals with genotypes), we applied a single-step single nucleotide polymorphism (SNP) best linear unbiased prediction approach to investigate the genetic determinism underlying resistance to this disease (heritability estimate and identification of relevant genomic regions) and estimated genetic trends, reliability, and relative risk factors associated with genomic predictions.

Results

Resistance to JD was moderately heritable (0.14) and 16 genomic regions were detected that accounted for at least 0.05% of the breeding values variance (GV) in resistance to JD, and were located on chromosomes 1, 3, 5, 6, 7, 19, 20, 21, 23, 25, and 27, with the highest percentage of variance explained by regions on chromosomes 23 (0.36% GV), 5 (0.22% GV), 1 (0.14% GV), and 3 (0.13% GV). When estimated for the whole chromosomes, the autosomes with the largest overall contributions were chromosomes 3 (5.3% GV), 10 (4.8%), 23 (4.7%), 1 (3.6%), 7 (3.4%), 5 (2.9%), 12 (2.5%), 11 (2.2%), and 13 (2%). We estimated a slightly favorable genetic trend in resistance to JD over the last two decades, which can be explained by a low positive genetic correlation between resistance to JD and total merit index (+ 0.06). Finally, in a validation population of 907 cows, relatively reliable genomic predictions (reliability = 0.55) were obtained, which allowed the identification of cows at high risk of infection.

Conclusions

This study provides new insights into the genetic determinism of resistance to JD and shows that this trait can be predicted from SNP genotypes. It has led to the implementation of a single-step genomic evaluation that should rapidly become an effective tool for controlling paratuberculosis on French Holstein farms.

Can Infectious Disease Control Be Achieved without Antibiotics by Exploiting Mechanisms of Disease Tolerance?

Antimicrobial use in animal agriculture may be contributing to the emerging public health crisis of antimicrobial resistance. The sustained prevalence of infectious diseases driving antimicrobial use industry-wide suggests that traditional methods of bolstering disease resistance are, for some diseases, ineffective. A paradigm shift in our approach to infectious disease control is needed to reduce antimicrobial use and sustain animal and human health and the global economy. Targeting the defensive mechanisms that promote the health of an infected host without impacting pathogen fitness, termed “disease tolerance,” is a novel disease control approach ripe for discovery. This article presents examples of disease tolerance dictating clinical outcomes for several infectious diseases in humans, reveals evidence suggesting a similarly critical role of disease tolerance in the progression of infectious diseases plaguing animal agriculture, and thus substantiates the assertion that exploiting disease tolerance mechanisms can positively impact animal and human health.

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.

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.

Selection for Favorable Health Traits: A Potential Approach to Cope with Diseases in Farm Animals

Disease is a global problem for animal farming industries causing tremendous economic losses (>USD 220 billion over the last decade) and serious animal welfare issues. The limitations and deficiencies of current non-selection disease control methods (e.g., vaccination, treatment, eradication strategy, genome editing, and probiotics) make it difficult to effectively, economically, and permanently eliminate the adverse influences of disease in the farm animals. These limitations and deficiencies drive animal breeders to be more concerned and committed to dealing with health problems in farm animals by selecting animals with favorable health traits. Both genetic selection and genomic selection contribute to improving the health of farm animals by selecting certain health traits (e.g., disease tolerance, disease resistance, and immune response), although both of them face some challenges. The objective of this review was to comprehensively review the potential of selecting health traits in coping with issues caused by diseases in farm animals. Within this review, we highlighted that selecting health traits can be applied as a method of disease control to help animal agriculture industries to cope with the adverse influences caused by diseases in farm animals. Certainly, the genetic/genomic selection solution cannot solve all the disease problems in farm animals. Therefore, management, vaccination, culling, medical treatment, and other measures must accompany selection solution to reduce the adverse impact of farm animal diseases on profitability and animal welfare.

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.

Integrative Analyses of Long Non-coding RNA and mRNA Involved in Piglet Ileum Immune Response to Clostridium perfringens Type C Infection

Long non-coding RNAs (lncRNAs) have been shown to play important roles in regulating host immune and inflammatory responses to bacterial infection. Infection with Clostridium perfringens (C. perfringens), a food-borne zoonotic pathogen, can lead to a series of inflammatory diseases in human and piglet, greatly challenging the healthy development of global pig industry. However, the roles of lncRNAs involved in piglet immune response against C. perfringens type C infection remain unknown. In this study, the regulatory functions of ileum lncRNAs and mRNAs were investigated in piglet immune response to C. perfringens type C infection among resistance (IR), susceptibility (IS) and sham-inoculation (control, IC) groups. A total of 480 lncRNAs and 3,669 mRNAs were significantly differentially expressed, the differentially expressed lncRNAs and mRNAs in the IR and IS groups were enriched in various pathways of ABC transporters, olfactory transduction, PPAR signaling pathway, chemokine signaling pathway and Toll-like receptor signaling pathway, involving in regulating piglet immune responses and resistance during infection. There were 212 lncRNAs and 505 target mRNAs found to have important association with C. perfringens infectious diseases, furthermore, 25 dysregulated lncRNAs corresponding to 13 immune-related target mRNAs were identified to play potential roles in defense against bacterial infection. In conclusion, the results improve our understanding on the characteristics of lncRNAs and mRNAs on regulating host immune response against C. perfringens type C infection, which will provide a reference for future research into exploring C. perfringens-related diseases in human.

Novel Insights reveal Anti-microbial Gene Regulation of Piglet Intestine Immune in response to Clostridium perfringens Infection

LncRNA play important roles in regulation of host immune and inflammation responses in defending bacterial infection. Clostridium perfringens (C. perfringens) type C is one of primary bacteria leading to piglet diarrhea and other intestinal inflammatory diseases. For the differences of host immune capacity, individuals usually show resistance and susceptibility to bacterial infection. However, whether and how lncRNAs involved in modulating host immune resistance have not been reported. We have investigated the expression patterns of ileum lncRNAs of 7-day-old piglets infected by C. perfringens type C through RNA sequencing. A total of 16 lncRNAs and 126 mRNAs were significantly differentially expressed in resistance (IR) and susceptibility (IS) groups. Many lncRNAs and mRNAs were identified to regulate resistance and susceptibility of piglets through immune related pathways. Five lncRNAs may have potential function on regulating the expressions of cytokines, these lncRNAs and cytokines work together to co-regulated piglet immune response to C. perfringens, affecting host resistance and susceptibility. These results provide valuable information for understanding the functions of lncRNA and mRNA in affecting piglet diarrhea resistance of defensing to C. perfringens type C, these lncRNAs and mRNAs may be used as the important biomarkers for decreasing C. perfringens spread and diseases in human and piglets.

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.

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

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

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.

Genomic Tools and Animal Health

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

Genome-wide association study of periweaning failure-to-thrive syndrome (PFTS) in pigs

Porcine periweaning-failure-to-thrive syndrome (PFTS) is a condition that affects newly weaned piglets. It is characterised by a progressive debilitation leading to death, in the absence of infectious, nutritional, management or environmental factors. In this study, we present the first report of PFTS in South America and the results of a genome-wide association study to identify the genetic markers associated with the appearance of this condition in a crossbred swine population. Four chromosomal regions were associated with PFTS predisposition, one located on SSCX, one on SSC8, and the two other regions on SSC14. Regions on SSC8 and SSC14 harbour important functional candidate genes involved in human depression and might have an important role in PFTS. Our findings contribute to the increasing knowledge about this syndrome, which has been investigated since 2007, and to the identification of the aetiology of this disease.

Genetic diversity analysis of two commercial breeds of pigs using genomic and pedigree data

Background

Genetic improvement in livestock populations can be achieved without significantly affecting genetic diversity if mating systems and selection decisions take genetic relationships among individuals into consideration. The objective of this study was to examine the genetic diversity of two commercial breeds of pigs. Genotypes from 1168 Landrace (LA) and 1094 Large White (LW) animals from a commercial breeding program in Brazil were obtained using the Illumina PorcineSNP60 Beadchip. Inbreeding estimates based on pedigree (F_x) and genomic information using runs of homozygosity (F_ROH) and the single nucleotide polymorphisms (SNP) by SNP inbreeding coefficient (F_SNP) were obtained. Linkage disequilibrium (LD), correlation of linkage phase (r) and effective population size (N_e) were also estimated.

Results

Estimates of inbreeding obtained with pedigree information were lower than those obtained with genomic data in both breeds. We observed that the extent of LD was slightly larger at shorter distances between SNPs in the LW population than in the LA population, which indicates that the LW population was derived from a smaller N_e. Estimates of N_e based on genomic data were equal to 53 and 40 for the current populations of LA and LW, respectively. The correlation of linkage phase between the two breeds was equal to 0.77 at distances up to 50 kb, which suggests that genome-wide association and selection should be performed within breed. Although selection intensities have been stronger in the LA breed than in the LW breed, levels of genomic and pedigree inbreeding were lower for the LA than for the LW breed.

Conclusions

The use of genomic data to evaluate population diversity in livestock animals can provide new and more precise insights about the effects of intense selection for production traits. Resulting information and knowledge can be used to effectively increase response to selection by appropriately managing the rate of inbreeding, minimizing negative effects of inbreeding depression and therefore maintaining desirable levels of genetic diversity.

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.

Unravelling the genetic components involved in the immune response of pigs vaccinated against influenza virus

A genome-wide association study for immune response to influenza vaccination in a crossbred swine population was conducted. Swine influenza is caused by influenza A virus (FLUAV) which is considered one of the most prevalent respiratory pathogens in swine worldwide. The main strategy used to control influenza in swine herds is through vaccination. However, the currently circulating FLUAV subtypes in swine are genetically and antigenically diverse and their interaction with the host genetics poses a challenge for the production of efficacious and cross-protective vaccines. In this study, 103 pigs vaccinated with an inactivated H1N1 pandemic virus were genotyped with the Illumina PorcineSNP60V2 BeadChip for the identification of genetic markers associated with immune response efficacy to influenza A virus vaccination. Immune response was measured based on the presence or absence of HA (hemagglutinin) and NP (nucleoprotein) antibodies induced by vaccination and detected in swine sera by the hemagglutination inhibition (HI) and ELISA assays, respectively. The ELISA test was also used as a measurement of antibody levels produced following the FLUAV vaccination. Associations were tested with x(2) test for a case and control data and using maximum likelihood method for the quantitative data, where a moderate association was considered if p<5×10(-5). When testing the association using the HI results, three markers with unknown location and three located on chromosomes SSCX, SSC14 and SSC18 were identified as associated with the immune response. Using the response to vaccination measured by ELISA as a qualitative and quantitative phenotype, four genomic regions were associated with immune response: one on SSC12 and three on chromosomes SSC1, SSC7, and SSC15, respectively. Those regions harbor important functional candidate genes possibly involved with the degree of immune response to vaccination. These results show an important role of host genetics in the immune response to influenza vaccination. Genetic selection for pigs with better response to FLUAV vaccination might be an alternative to reduce the impact of influenza virus infection in the swine industry. However, these results should to be validated in additional populations before its use.

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.

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.

Evidence of passive faecal shedding of Mycobacterium avium subsp. paratuberculosis in a Limousin cattle herd

It has been suggested that passive shedding of Mycobacterium avium subsp. paratuberculosis (MAP) in faeces may occur, but reliable data are missing. Passive shedding assumes the ingestion of MAP in contaminated feed and passive passage through the gastrointestinal tract without causing infection. In this study the presence of MAP in faeces in a closed herd of Limousin cattle was monitored for 53 months using quantitative real time PCR (qPCR) and culture. The initial prevalence of MAP in the herd was determined to be 63.4% and 4.9% using qPCR and culture, respectively. After the removal of two culture- and qPCR-positive (>10(4) MAP cells/g) cows, the prevalence of MAP using qPCR decreased to 42.1% and later to 15.6% and 6.7%. The continuous removal of suspected animals from the herd during the monitoring period minimised the presence of MAP in faeces to sporadic, which may have resulted from a decrease in the environmental infectious pressure. The findings suggest that the presence of low numbers of MAP in bovine faeces may not necessarily be caused by real infection, but rather by passive passage of MAP. This phenomenon should therefore be considered when interpreting MAP qPCR data.

TLR4 single nucleotide polymorphisms (SNPs) associated with Salmonella shedding in pigs

Toll-like receptor 4 (TLR4) is a key factor in the innate immune recognition of lipopolysaccharide (LPS) from Gram-negative bacteria. Previous studies from our group identified differences in the expression profile of TLR4 and genes affected by the TLR4 signaling pathway among pigs that shed varying levels of Salmonella, a Gram-negative bacterium. Therefore, genetic variation in this gene may be involved with the host’s immune response to bacterial infections. The current study screened for single nucleotide polymorphisms (SNPs) in the TLR4 gene and tested their association with Salmonella fecal shedding. Pigs (n = 117) were intranasally challenged at 7 weeks of age with 1 × 10⁹ CFU of S. Typhimurium χ4232 and were classified as low or persistent Salmonella shedders based on the levels of Salmonella being excreted in fecal material. Salmonella fecal shedding was determined by quantitative bacteriology on days 2, 7, 14, and 20/21 post exposure, and the cumulative levels of Salmonella were calculated to identify the low (n = 20) and persistent (n = 20) Salmonella shedder pigs. From those 40 animals, the TLR4 region was sequenced, and 18 single nucleotide polymorphisms (SNPs) in TLR4 were identified. Twelve SNPs have been previously described and six are novel SNPs of which five are in the 5′ untranslated region and one is in intron 2. Single marker association test identified 13 SNPs associated with the qualitative trait of Salmonella fecal shedding, and seven of those SNPs were also associated with a quantitative measurement of fecal shedding (P < 0.05). Using a stepwise regression process, a haplotype composed of SNPs rs80787918 and rs80907449 (P ≤ 4.0 × 10⁻³) spanning a region of 4.9 Kb was identified, thereby providing additional information of the influence of those SNPs on Salmonella fecal shedding in pigs.

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.

Effects of fertility on gene expression and function of the bovine endometrium

Infertility and subfertility are important and pervasive reproductive problems in both domestic animals and humans. The majority of embryonic loss occurs during the first three weeks of pregnancy in cattle and women due, in part, to inadequate endometrial receptivity for support of embryo implantation. To identify heifers of contrasting fertility, serial rounds of artificial insemination (AI) were conducted in 201 synchronized crossbred beef heifers. The heifers were then fertility classified based on number of pregnancies detected on day 35 in four AI opportunities. Heifers, classified as having high fertility, subfertility or infertility, were selected for further study. The fertility-classified heifers were superovulated and flushed, and the recovered embryos were graded and then transferred to synchronized recipients. Quantity of embryos recovered per flush, embryo quality, and subsequent recipient pregnancy rates did not differ by fertility classification. Two in vivo-produced bovine embryos (stage 4 or 5, grade 1 or 2) were then transferred into each heifer on day 7 post-estrus. Pregnancy rates were greater in high fertility than lower fertility heifers when heifers were used as embryo recipients. The reproductive tracts of the classified heifers were obtained on day 14 of the estrous cycle. No obvious morphological differences in reproductive tract structures and histology of the uterus were observed in the heifers. Microarray analysis revealed differences in the endometrial transcriptome based on fertility classification. A genome-wide association study, based on SNP genotyping, detected 7 moderate associations with fertility across 6 different chromosomes. Collectively, these studies support the idea that innate differences in uterine function underlie fertility and early pregnancy loss in ruminants. Cattle with defined early pregnancy success or loss is useful to elucidate the complex biological and genetic mechanisms governing endometrial receptivity and uterine competency for pregnancy.

Genomics of complex traits

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

The molecular pathways underlying host resistance and tolerance to pathogens

Breeding livestock that are better able to withstand the onslaught of endemic- and exotic pathogens is high on the wish list of breeders and farmers world-wide. However, the defense systems in both pathogens and their hosts are complex and the degree of genetic variation in resistance and tolerance will depend on the trade-offs that they impose on host fitness as well as their life-histories. The genes and pathways underpinning resistance and tolerance traits may be distinct or intertwined as the outcome of any infection is a result of a balance between collateral damage of host tissues and control of the invading pathogen. Genes and molecular pathways associated with resistance are mainly expressed in the mucosal tract and the innate immune system and control the very early events following pathogen invasion. Resistance genes encode receptors involved in uptake of pathogens, as well as pattern recognition receptors (PRR) such as the toll-like receptor family as well as molecules involved in strong and rapid inflammatory responses which lead to rapid pathogen clearance, yet do not lead to immunopathology. In contrast tolerance genes and pathways play a role in reducing immunopathology or enhancing the host's ability to protect against pathogen associated toxins. Candidate tolerance genes may include cytosolic PRRs and unidentified sensors of pathogen growth, perturbation of host metabolism and intrinsic danger or damage associated molecules. In addition, genes controlling regulatory pathways, tissue repair and resolution are also tolerance candidates. The identities of distinct genetic loci for resistance and tolerance to infectious pathogens in livestock species remain to be determined. A better understanding of the mechanisms involved and phenotypes associated with resistance and tolerance should ultimately help to improve livestock health and welfare.

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.

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

Background

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

Methodology/Principal Findings

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

Conclusion/Significance

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

Genetic susceptibility to paratuberculosis

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

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

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

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

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

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

SNP-based gene-set enrichment analysis from single nucleotide polymorphisms, or GSEA-SNP, is a tool to identify candidate genes based on enrichment analysis of sets of genes rather than single SNP associations. The objective of this study was to identify modest-effect genes associated with Mycobacterium avium subsp. paratuberculosis (Map) tissue infection or fecal shedding using GSEA-SNP applied to KEGG pathways or Gene Ontology (GO) gene sets. The Illumina Bovine SNP50 BeadChip was used to genotype 209 Holstein cows for the GSEA-SNP analyses. For each of 13,744 annotated genes genome-wide located within 50 kb of a Bovine SNP50 SNP, the single SNP with the highest Cochran-Armitage Max statistic was used as a proxy statistic for that gene's strength of affiliation with Map. Gene-set enrichment was tested using a weighted Kolmogorov-Smirnov-like running sum statistic with data permutation to adjust for multiple testing. For tissue infection and fecal shedding, no gene sets in KEGG pathways or in GO sets for molecular function or cellular component were enriched for signal. The GO biological process gene set for positive regulation of cell motion (GO:0051272, q = 0.039, 5/11 genes contributing to the core enrichment) was enriched for Map tissue infection, while no GO biological process gene sets were enriched for fecal shedding. GSEA-SNP complements traditional SNP association approaches to identify genes of modest effects as well as genes with larger effects as demonstrated by the identification of one locus that we previously found to be associated with Map tissue infection using a SNP-by-SNP genome-wide association study.