Quantitative trait loci for clinical mastitis on chromosomes 2, 6, 14 and 20 in Norwegian Red cattle

Meta-analysis of six dairy cattle breeds reveals biologically relevant candidate genes for mastitis resistance

BACKGROUND

Mastitis is a disease that incurs significant costs in the dairy industry. A promising approach to mitigate its negative effects is to genetically improve the resistance of dairy cattle to mastitis. A meta-analysis of genome-wide association studies (GWAS) across multiple breeds for clinical mastitis (CM) and its indicator trait, somatic cell score (SCS), is a powerful method to identify functional genetic variants that impact mastitis resistance.

RESULTS

We conducted meta-analyses of eight and fourteen GWAS on CM and SCS, respectively, using 30,689 and 119,438 animals from six dairy cattle breeds. Methods for the meta-analyses were selected to properly account for the multi-breed structure of the GWAS data. Our study revealed 58 lead markers that were associated with mastitis incidence, including 16 loci that did not overlap with previously identified quantitative trait loci (QTL), as curated at the Animal QTLdb. Post-GWAS analysis techniques such as gene-based analysis and genomic feature enrichment analysis enabled prioritization of 31 candidate genes and 14 credible candidate causal variants that affect mastitis.

CONCLUSIONS

Our list of candidate genes can help to elucidate the genetic architecture underlying mastitis resistance and provide better tools for the prevention or treatment of mastitis, ultimately contributing to more sustainable animal production.

Genome-wide association and functional genomic analyses for various hoof health traits in North American Holstein cattle

Hoof diseases are a major welfare and economic issue in the global dairy cattle production industry, which can be minimized through improved management and breeding practices. Optimal genetic improvement of hoof health could benefit from a deep understanding of the genetic background and biological underpinning of indicators of hoof health. Therefore, the primary objectives of this study were to perform genome-wide association studies, using imputed high-density genetic markers data from North American Holstein cattle, for 8 hoof-related traits: digital dermatitis, sole ulcer, sole hemorrhage, white line lesion, heel horn erosion, interdigital dermatitis, interdigital hyperplasia, and toe ulcer, and a hoof health index. De-regressed estimated breeding values from 25,580 Holstein animals were used as pseudo-phenotypes for the association analyses. The genomic quality control, genotype phasing, and genotype imputation were performed using the PLINK (version 1.9), Eagle (version 2.4.1), and Minimac4 software, respectively. The functional genomic analyses were performed using the GALLO R package and the DAVID platform. We identified 22, 34, 14, 22, 28, 33, 24, 43, and 15 significant markers for digital dermatitis, heel horn erosion, interdigital dermatitis, interdigital hyperplasia, sole hemorrhage, sole ulcer, toe ulcer, white line lesion disease, and the hoof health index, respectively. The significant markers were located across all autosomes, except BTA10, BTA12, BTA20, BTA26, BTA27, and BTA28. Moreover, the genomic regions identified overlap with various previously reported quantitative trait loci for exterior, health, meat and carcass, milk, production, and reproduction traits. The enrichment analyses identified 44 significant gene ontology terms. These enriched genomic regions harbor various candidate genes previously associated with bone development, metabolism, and infectious and immunological diseases. These findings indicate that hoof health traits are highly polygenic and influenced by a wide range of biological processes.

Impact of temporal variations and risk factors associated with udder inflammation in Hardhenu cattle (Bos taurus × Bos indicus)

This study aimed to explore the occurrence and risk factors associated with clinical mastitis within the Hardhenu cattle herd over a span of 14 years (2008-2021). A comprehensive analysis of 1515 lactation records was conducted to ascertain the incidence of clinical mastitis. The investigation determined an overall incidence rate of 26.80% in the studied population. A significant relationship between the year and clinical mastitis incidence was established through Chi-square analysis (p < .05). Temporal variations in clinical mastitis odds were apparent, with the highest odds (ranging from 0.91 to 1.00) observed during the initial years of 2008-2009 and 2009-2010. Logistic regression revealed that odds values for clinical mastitis incidence were highest in 2008-2009 (1.00), succeeded by 2009-2010 (0.91), 2012-2013 (0.88), 2018-2019 (0.67) and reaching the lowest in 2021-2022 (0.35). Subsequent rankings included 2010-2011 (0.39), 2014-2015 (0.43) and 2019-2020 (0.45). Parity was found to be significantly associated with clinical mastitis occurrence. When compared to Parity 3, both Parity 1 (odds ratio: 1.516, 95% confidence interval: 0.881-2.612) and Parity 2 (odds ratio: 2.626, 95% confidence interval: 1.568-4.398) exhibited higher odds values for clinical mastitis incidence. While the period of calving did not exert a significant influence on clinical mastitis incidence, a heightened occurrence was observed during the rainy season within the targeted population. These findings offer valuable insights into the patterns of incidence, temporal fluctuations, and non-genetic determinants impacting clinical mastitis within the Hardhenu cattle. The implications of this study can facilitate the development of targeted strategies and management protocols aimed at enhancing udder health and overall productivity in dairy cattle.

Candidate genes for mastitis resistance in dairy cattle: a data integration approach

BACKGROUND

Inflammation of the mammary tissue (mastitis) is one of the most detrimental health conditions in dairy ruminants and is considered the most economically important infectious disease of the dairy sector. Improving mastitis resistance is becoming an important goal in dairy ruminant breeding programmes. However, mastitis resistance is a complex trait and identification of mastitis-associated alleles in livestock is difficult. Currently, the only applicable approach to identify candidate loci for complex traits in large farm animals is to combine different information that supports the functionality of the identified genomic regions with respect to a complex trait.

METHODS

To identify the most promising candidate loci for mastitis resistance we integrated heterogeneous data from multiple sources and compiled the information into a comprehensive database of mastitis-associated candidate loci. Mastitis-associated candidate genes reported in association, expression, and mouse model studies were collected by searching the relevant literature and databases. The collected data were integrated into a single database, screened for overlaps, and used for gene set enrichment analysis.

RESULTS

The database contains candidate genes from association and expression studies and relevant transgenic mouse models. The 2448 collected candidate loci are evenly distributed across bovine chromosomes. Data integration and analysis revealed overlaps between different studies and/or with mastitis-associated QTL, revealing promising candidate genes for mastitis resistance.

CONCLUSION

Mastitis resistance is a complex trait influenced by numerous alleles. Based on the number of independent studies, we were able to prioritise candidate genes and propose a list of the 22 most promising. To our knowledge this is the most comprehensive database of mastitis associated candidate genes and could be helpful in selecting genes for functional validation studies.

New Insights on Nucleotide Sequence Variants and mRNA Levels of Candidate Genes Assessing Resistance/Susceptibility to Mastitis in Holstein and Montbéliarde Dairy Cows

A major factor in the propagation of an infectious disease is host genetics. In this study, 180 dairy cows (90 of each breed: Holstein and Montbéliarde) were used. Each breed's tested dairy cows were divided into two groups of comparable size (45 cows each), mastitis-free and mastitis-affected groups. Each cow's jugular vein was punctured to obtain blood samples for DNA and RNA extraction. In the examined Holstein and Montbéliarde dairy cows, single nucleotide polymorphisms (SNPs) related with mastitis resistance/susceptibility were found in the RASGRP1, NFkB, CHL1, MARCH3, PDGFD, MAST3, EPS15L1, C1QTNF3, CD46, COX18, NEURL1, PPIE, and PTX3 genes. Chi-square analysis of identified SNPs revealed a significant difference in gene frequency between mastitic and healthy cows. Except for CHL1, mastitic dairy cows of two breeds had considerably higher mRNA levels of the examined genes than did healthy ones. Marker-assisted selection and monitoring of dairy cows' susceptibility to mastitis may be accomplished through the use of discovered SNPs and changes in the gene expression profile of the studied genes. These findings also point to a possible method for reducing mastitis in dairy cows through selective breeding of animals using genetic markers linked to an animal's ability to resist infection.

Multi-Omics Integration and Network Analysis Reveal Potential Hub Genes and Genetic Mechanisms Regulating Bovine Mastitis

Mastitis, inflammation of the mammary gland, is the most prevalent disease in dairy cattle that has a potential impact on profitability and animal welfare. Specifically designed multi-omics studies can be used to prioritize candidate genes and identify biomarkers and the molecular mechanisms underlying mastitis in dairy cattle. Hence, the present study aimed to explore the genetic basis of bovine mastitis by integrating microarray and RNA-Seq data containing healthy and mastitic samples in comparative transcriptome analysis with the results of published genome-wide association studies (GWAS) using a literature mining approach. The integration of different information sources resulted in the identification of 33 common and relevant genes associated with bovine mastitis. Among these, seven genes—CXCR1, HCK, IL1RN, MMP9, S100A9, GRO1, and SOCS3—were identified as the hub genes (highly connected genes) for mastitis susceptibility and resistance, and were subjected to protein-protein interaction (PPI) network and gene regulatory network construction. Gene ontology annotation and enrichment analysis revealed 23, 7, and 4 GO terms related to mastitis in the biological process, molecular function, and cellular component categories, respectively. Moreover, the main metabolic-signalling pathways responsible for the regulation of immune or inflammatory responses were significantly enriched in cytokine–cytokine-receptor interaction, the IL-17 signaling pathway, viral protein interaction with cytokines and cytokine receptors, and the chemokine signaling pathway. Consequently, the identification of these genes, pathways, and their respective functions could contribute to a better understanding of the genetics and mechanisms regulating mastitis and can be considered a starting point for future studies on bovine mastitis.

Suitability of GWAS as a Tool to Discover SNPs Associated with Tick Resistance in Cattle: A Review

Understanding the biological mechanisms underlying tick resistance in cattle holds the potential to facilitate genetic improvement through selective breeding. Genome wide association studies (GWAS) are popular in research on unraveling genetic determinants underlying complex traits such as tick resistance. To date, various studies have been published on single nucleotide polymorphisms (SNPs) associated with tick resistance in cattle. The discovery of SNPs related to tick resistance has led to the mapping of associated candidate genes. Despite the success of these studies, information on genetic determinants associated with tick resistance in cattle is still limited. This warrants the need for more studies to be conducted. In Africa, the cost of genotyping is still relatively expensive; thus, conducting GWAS is a challenge, as the minimum number of animals recommended cannot be genotyped. These population size and genotype cost challenges may be overcome through the establishment of collaborations. Thus, the current review discusses GWAS as a tool to uncover SNPs associated with tick resistance, by focusing on the study design, association analysis, factors influencing the success of GWAS, and the progress on cattle tick resistance studies.

A 12 kb multi-allelic copy number variation encompassing a GC gene enhancer is associated with mastitis resistance in dairy cattle

Clinical mastitis (CM) is an inflammatory disease occurring in the mammary glands of lactating cows. CM is under genetic control, and a prominent CM resistance QTL located on chromosome 6 was reported in various dairy cattle breeds. Nevertheless, the biological mechanism underpinning this QTL has been lacking. Herein, we mapped, fine-mapped, and discovered the putative causal variant underlying this CM resistance QTL in the Dutch dairy cattle population. We identified a ~12 kb multi-allelic copy number variant (CNV), that is in perfect linkage disequilibrium with a lead SNP, as a promising candidate variant. By implementing a fine-mapping and through expression QTL mapping, we showed that the group-specific component gene (GC), a gene encoding a vitamin D binding protein, is an excellent candidate causal gene for the QTL. The multiplicated alleles are associated with increased GC expression and low CM resistance. Ample evidence from functional genomics data supports the presence of an enhancer within this CNV, which would exert cis-regulatory effect on GC. We observed that strong positive selection swept the region near the CNV, and haplotypes associated with the multiplicated allele were strongly selected for. Moreover, the multiplicated allele showed pleiotropic effects for increased milk yield and reduced fertility, hinting that a shared underlying biology for these effects may revolve around the vitamin D pathway. These findings together suggest a putative causal variant of a CM resistance QTL, where a cis-regulatory element located within a CNV can alter gene expression and affect multiple economically important traits.

Improvement of Disease Resistance in Livestock: Application of Immunogenomics and CRISPR/Cas9 Technology

Disease occurrence adversely affects livestock production and animal welfare, and have an impact on both human health and public perception of food-animals production. Combined efforts from farmers, animal scientists, and veterinarians have been continuing to explore the effective disease control approaches for the production of safe animal-originated food. Implementing the immunogenomics, along with genome editing technology, has been considering as the key approach for safe food-animal production through the improvement of the host genetic resistance. Next-generation sequencing, as a cutting-edge technique, enables the production of high throughput transcriptomic and genomic profiles resulted from host-pathogen interactions. Immunogenomics combine the transcriptomic and genomic data that links to host resistance to disease, and predict the potential candidate genes and their genomic locations. Genome editing, which involves insertion, deletion, or modification of one or more genes in the DNA sequence, is advancing rapidly and may be poised to become a commercial reality faster than it has thought. The clustered regulatory interspaced short palindromic repeats (CRISPR)/CRISPR-associated protein 9 (Cas9) [CRISPR/Cas9] system has recently emerged as a powerful tool for genome editing in agricultural food production including livestock disease management. CRISPR/Cas9 mediated insertion of NRAMP1 gene for producing tuberculosis resistant cattle, and deletion of CD163 gene for producing porcine reproductive and respiratory syndrome (PRRS) resistant pigs are two groundbreaking applications of genome editing in livestock. In this review, we have highlighted the technological advances of livestock immunogenomics and the principles and scopes of application of CRISPR/Cas9-mediated targeted genome editing in animal breeding for disease resistance.

Genome-Wide Association Studies for the Concentration of Albumin in Colostrum and Serum in Chinese Holstein

Albumin can be of particular benefit in fighting infections for newborn calves due to its anti-inflammatory and anti-oxidative stress properties. To identify the candidate genes related to the concentration of albumin in colostrum and serum, we collected the colostrum and blood samples from 572 Chinese Holstein cows within 24 h after calving and measured the concentration of albumin in the colostrum and serum using the ELISA methods. The cows were genotyped with GeneSeek 150 K chips (containing 140,668 single nucleotide polymorphisms; SNPs). After quality control, we performed GWASs via GCTA software with 91,620 SNPs and 563 cows. Consequently, 9 and 7 genome-wide significant SNPs (false discovery rate (FDR) at 1%) were identified. Correspondingly, 42 and 206 functional genes that contained or were approximate to (±1 Mbp) the significant SNPs were acquired. Integrating the biological process of these genes and the reported QTLs for immune and inflammation traits in cattle, 3 and 12 genes were identified as candidates for the concentration of colostrum and serum albumin, respectively; these are RUNX1, CBR1, OTULIN,CDK6, SHARPIN, CYC1, EXOSC4, PARP10, NRBP2, GFUS, PYCR3, EEF1D, GSDMD, PYCR2 and CXCL12. Our findings provide important information for revealing the genetic mechanism behind albumin concentration and for molecular breeding of disease-resistance traits in dairy cattle.

Time- and population-dependent genetic patterns underlie bovine milk somatic cell count

The objective of this study was to determine whether genetic regulation of bovine milk somatic cell count (SCC) varied throughout the course of an individual lactation and to identify quantitative trait loci (QTL) that may differentiate populations of chronically mastitic and robustly healthy cows. Milk SCC has long been a proxy for clinical mastitis diagnosis in management and genetic improvement strategies to control the disease. Cows (n = 471) were genotyped on the Illumina BovineHD 777K BeadChip (Illumina Inc., San Diego, CA), and composite milk samples were collected for SCC at 0-1 d in milk (DIM), 3-5 DIM, 10-14 DIM, 90-110 DIM, and 210-230 DIM, with each time span representing key physiological transitions for the cow. Median lactation somatic cell score (SCS) and area under the SCS curve were calculated from farm test data. A total of 8 genome-wide associations were performed and 167 SNP spanning the genome were significantly associated (false discovery rate <0.05). Of these associated regions, 27 of 48 associated QTL were novel for clinical mastitis or SCC. The linkage disequilibrium block surrounding the associated QTL or a 1-Mb window in the absence of linkage disequilibrium was interrogated for candidate genes, and many of those identified were related to multiple arms of the immune system, including toll-like receptor signaling, macrophage activation, B-cell maturation, T-cell recruitment, and the complement pathway. These genes included EXOC4, BAMBI, ITSN2, IL34, FCN3, CD8A, and CD8B. In addition, we identified populations of robustly healthy (SCS ≤4 from 10-14 DIM until study end), chronically mastitic (SCS >4 from 10-14 DIM until study end), and average cows with fluctuating SCS, and calculated fixation indices to identify regions of the genome differentiating these 3 populations. A total of 12 SNP were identified that showed moderate allelic differentiation (Wright's F statistic, F_ST ≥ 0.4) between the "chronic," "healthy," and "average" populations of cows. Candidate genes in the region surrounding differentiated QTL were related to cell signaling and immune response, such as JAKMIP1 and MADCAM1. The wide range of significantly associated QTL spanning the genome and the diversity of gene functions reinforces that mastitis is a complex trait and suggests that selection based on lactation stage-specific SCS rather than a generalized score may lead to greater success in breeding mastitis-resistant cows.

Ancestral Haplotype Mapping for GWAS and Detection of Signatures of Selection in Admixed Dairy Cattle of Kenya

Understanding the genetic structure of adaptation and productivity in challenging environments is necessary for designing breeding programs that suit such conditions. Crossbred dairy cattle in East Africa resulting from over 60 years of crossing exotic dairy breeds with indigenous cattle plus inter se matings form a highly variable admixed population. This population has been subject to natural selection in response to environmental stresses, such as harsh climate, low-quality feeds, poor management, and strong disease challenge. Here, we combine two complementary sets of analyses, genome-wide association (GWA) and signatures of selection (SoS), to identify genomic regions that contribute to variation in milk yield and/or contribute to adaptation in admixed dairy cattle of Kenya. Our GWA separates SNP effects due to ancestral origin of alleles from effects due to within-population linkage disequilibrium. The results indicate that many genomic regions contributed to the high milk production potential of modern dairy breeds with no region having an exceptional effect. For SoS, we used two haplotype-based tests to compare haplotype length variation within admixed and between admixed and East African Shorthorn Zebu cattle populations. The integrated haplotype score (iHS) analysis identified 16 candidate regions for positive selection in the admixed cattle while the between population Rsb test detected 24 divergently selected regions in the admixed cattle compared to East African Shorthorn Zebu. We compare the results from GWA and SoS in an attempt to validate the most significant SoS results. Only four candidate regions for SoS intersect with GWA regions using a low stringency test. The identified SoS candidate regions harbored genes in several enriched annotation clusters and overlapped with previously found QTLs and associations for different traits in cattle. If validated, the GWA and SoS results indicate potential for SNP-based genomic selection for genetic improvement of smallholder crossbred cattle.

Distinguishing pleiotropy from linked QTL between milk production traits and mastitis resistance in Nordic Holstein cattle

BACKGROUND

Production and health traits are central in cattle breeding. Advances in next-generation sequencing technologies and genotype imputation have increased the resolution of gene mapping based on genome-wide association studies (GWAS). Thus, numerous candidate genes that affect milk yield, milk composition, and mastitis resistance in dairy cattle are reported in the literature. Effect-bearing variants often affect multiple traits. Because the detection of overlapping quantitative trait loci (QTL) regions from single-trait GWAS is too inaccurate and subjective, multi-trait analysis is a better approach to detect pleiotropic effects of variants in candidate genes. However, large sample sizes are required to achieve sufficient power. Multi-trait meta-analysis is one approach to deal with this problem. Thus, we performed two multi-trait meta-analyses, one for three milk production traits (milk yield, protein yield and fat yield), and one for milk yield and mastitis resistance.

RESULTS

For highly correlated traits, the power to detect pleiotropy was increased by multi-trait meta-analysis compared with the subjective assessment of overlapping of single-trait QTL confidence intervals. Pleiotropic effects of lead single nucleotide polymorphisms (SNPs) that were detected from the multi-trait meta-analysis were confirmed by bivariate association analysis. The previously reported pleiotropic effects of variants within the DGAT1 and MGST1 genes on three milk production traits, and pleiotropic effects of variants in GHR on milk yield and fat yield were confirmed. Furthermore, our results suggested that variants in KCTD16, KCNK18 and ENSBTAG00000023629 had pleiotropic effects on milk production traits. For milk yield and mastitis resistance, we identified possible pleiotropic effects of variants in two genes, GC and DGAT1.

CONCLUSIONS

Multi-trait meta-analysis improves our ability to detect pleiotropic interactions between milk production traits and identifies variants with pleiotropic effects on milk production traits and mastitis resistance. In particular, this should contribute to better understand the biological mechanisms that underlie the unfavorable genetic correlation between milk yield and mastitis.

Identification of candidate genes affecting chronic subclinical mastitis in Norwegian Red cattle: combining genome-wide association study, topologically associated domains and pathway enrichment analysis

The aim of this study was to identify genes associated with chronic subclinical mastitis (SCM) in Norwegian Red (NR) cattle. Twelve SCM traits defined based on fixed threshold for test-day somatic cell count (SCC) were, together with lactation-average somatic cell score (LSCS) used for association and pathway enrichment analyses. A GWAS was performed on 3795 genotyped NR bulls with 777K SNP data and phenotypic information from 7 300 847 test-day SCC observations from 3 543 764 cows. At 5% chromosome-wide significance level 36 unique SNP were detected to be associated with one or more of the traits. These SNPs were analysed for linked genes using genomic positions of topologically associated domains (TAD). For the SCM traits with SCC >50 000 and >100 000 cells/ml on two test-days in a row and LSCS, the same top significant genes were identified - checkpoint clamp loader component (RAD17) and cyclin B1 (CCNB1). The SCM traits with SCC >250 000, 300 000, 350 000 or 400 000 cells/ml on two test-days in a row and D400 (number of days before the first case with SCC >400 000 cells/ml) displayed similar top significant genes: acyl-CoA thioesterase 2 and 4 (ACOT2; ACOT4). For the traits SCM200_3 (SCC >200 000 cells/ml on three test-days in a row) and SCM150, SCM200 (SCC >150 000; 200 000 cells/ml on two test-days in a row) a group of chemokine (C-X-C motif) ligand genes and the Fos proto-oncogene, AP-1 transcription factor subunit (FOS) gene, were identified. Further functional studies of these identified candidate genes are necessary to clarify their actual role in development of chronic SCM in NR cattle.

Genome-wide association study of milk and reproductive traits in dual-purpose Xinjiang Brown cattle

Background

Dual-purpose cattle are more adaptive to environmental challenges than single-purpose dairy or beef cattle. Balance among milk, reproductive, and mastitis resistance traits in breeding programs is therefore more critical for dual-purpose cattle to increase net income and maintain well-being. With dual-purpose Xinjiang Brown cattle adapted to the Xinjiang Region in northwestern China, we conducted genome-wide association studies (GWAS) to dissect the genetic architecture related to milk, reproductive, and mastitis resistance traits. Phenotypic data were collected for 2410 individuals measured during 1995–2017. By adding another 445 ancestors, a total of 2855 related individuals were used to derive estimated breeding values for all individuals, including the 2410 individuals with phenotypes. Among phenotyped individuals, we genotyped 403 cows with the Illumina 150 K Bovine BeadChip.

Results

GWAS were conducted with the FarmCPU (Fixed and random model circulating probability unification) method. We identified 12 markers significantly associated with six of the 10 traits under the threshold of 5% after a Bonferroni multiple test correction. Seven of these SNPs were in QTL regions previously identified to be associated with related traits. One identified SNP, BovineHD1600006691, was significantly associated with both age at first service and age at first calving. This SNP directly overlapped a QTL previously reported to be associated with calving ease. Within 160 Kb upstream and downstream of each significant SNP identified, we speculated candidate genes based on functionality. Four of the SNPs were located within four candidate genes, including CDH2, which is linked to milk fat percentage, and GABRG2, which is associated with milk protein yield.

Conclusions

These findings are beneficial not only for breeding through marker-assisted selection, but also for genome editing underlying the related traits to enhance the overall performance of dual-purpose cattle.

Three Novel Players: PTK2B, SYK, and TNFRSF21 Were Identified to Be Involved in the Regulation of Bovine Mastitis Susceptibility via GWAS and Post-transcriptional Analysis

Bovine mastitis is a common inflammatory disease caused by multiple factors in early lactation or dry period. Genome wide association studies (GWAS) can provide a convenient and effective strategy for understanding the biological basis of mastitis and better prevention. 2b-RADseq is a high-throughput sequencing technique that offers a powerful method for genome-wide genetic marker development and genotyping. In this study, single nucleotide polymorphisms (SNPs) of the immune-regulated gene correlative with mastitis were screened and identified by two stage association analysis via GWAS-2b-RADseq in Chinese Holstein cows. We have screened 10,058 high quality SNPs from 7,957,920 tags and calculated their allele frequencies. Twenty-seven significant SNPs were co-labeled in two GWAS analysis models [Bayesian (P < 0.001) and Logistic regression (P < 0.01)], and only three SNPs (rs75762330, C > T, PIC = 0.2999; rs88640083, A > G, PIC = 0.1676; rs20438858, G > A, PIC = 0.3366) were annotated to immune-regulated genes (PTK2B, SYK, and TNFRSF21). Identified three SNPs are located in non-coding regions with low or moderate genetic polymorphisms. However, independent sample population validation (Case-control study) data showed that three important SNPs (rs75762330, P < 0.025, OR > 1; rs88640083, P < 0.005, OR > 1; rs20438858, P < 0.001, OR < 1) were significantly associated with clinical mastitis trait. Importantly, PTK2B and SYK expression was down-regulated in both peripheral blood leukocytes (PBLs) of clinical mastitis cows and in vitro LPS (E. coli)-stimulated bovine mammary epithelial cells, while TNFRSF21 was up-regulated. Under the same conditions, expression of Toll-like receptor 4 (TLR4), AKT1, and pro-inflammatory factors (IL-1β and IL-8) were also up-regulated. Interestingly, network analysis indicated that PTK2B and SYK are co-expressed in innate immune signaling pathway of Chinese Holstein. Taken together, these results provided strong evidence for the study of SNPs in bovine mastitis, and revealed the role of SYK, PTK2B, and TNFRSF21 in bovine mastitis susceptibility/tolerance.

The genetic background of clinical mastitis in Holstein-Friesian cattle

Mastitis is an inflammatory disease of the mammary gland, which has a significant economic impact and is an animal welfare concern. This work examined the association between single nucleotide polymorphisms (SNPs) and copy number variations (CNVs) with the incidence of clinical mastitis (CM). Using information from 16 half-sib pairs of Holstein-Friesian cows (32 animals in total) we searched for genomic regions that differed between a healthy (no incidence of CM) and a mastitis-prone (multiple incidences of CM) half-sib. Three cows with average sequence depth of coverage below 10 were excluded, which left 13 half-sib pairs available for comparisons. In total, 191 CNV regions were identified, which were deleted in a mastitis-prone cow, but present in its healthy half-sib and overlapped in at least nine half-sib pairs. These regions overlapped with exons of 46 genes, among which APP (BTA1), FOXL2 (BTA1), SSFA2 (BTA2), OTUD3 (BTA2), ADORA2A (BTA17), TXNRD2 (BTA17) and NDUFS6 (BTA20) have been reported to influence CM. Moreover, two duplicated CNV regions present in nine healthy individuals and absent in their mastitis-affected half-sibs overlapped with exons of a cholinergic receptor nicotinic α 10 subunit on BTA15 and a novel gene (ENSBTAG00000008519) on BTA27. One CNV region deleted in nine mastitis-affected sibs overlapped with two neighbouring long non-coding RNA sequences located on BTA12. Single nucleotide polymorphisms with differential genotypes between a healthy and a mastitis-affected sib included 17 polymorphisms with alternate alleles in eight affected and healthy half-sib families. Three of these SNPs were located introns of genes: MET (BTA04), RNF122 (BTA27) and WRN (BTA27). In summary, structural polymorphisms in form of CNVs, putatively play a role in susceptibility to CM. Specifically, sequence deletions have a greater effect on reducing resistance against mastitis, than sequence duplications have on increasing resistance against the disease.

Prioritizing candidate genes post-GWAS using multiple sources of data for mastitis resistance in dairy cattle

Background

Improving resistance to mastitis, one of the costliest diseases in dairy production, has become an important objective in dairy cattle breeding. However, mastitis resistance is influenced by many genes involved in multiple processes, including the response to infection, inflammation, and post-infection healing. Low genetic heritability, environmental variations, and farm management differences further complicate the identification of links between genetic variants and mastitis resistance. Consequently, studies of the genetics of variation in mastitis resistance in dairy cattle lack agreement about the responsible genes.

Results

We associated 15,552,968 imputed whole-genome sequencing markers for 5147 Nordic Holstein cattle with mastitis resistance in a genome-wide association study (GWAS). Next, we augmented P-values for markers in genes in the associated regions using Gene Ontology terms, Kyoto Encyclopedia of Genes and Genomes pathway analysis, and mammalian phenotype database. To confirm results of gene-based analyses, we used gene expression data from E. coli-challenged cow udders. We identified 22 independent quantitative trait loci (QTL) that collectively explained 14% of the variance in breeding values for resistance to clinical mastitis (CM). Using association test statistics with multiple pieces of independent information on gene function and differential expression during bacterial infection, we suggested putative causal genes with biological relevance for 12 QTL affecting resistance to CM in dairy cattle.

Conclusion

Combining information on the nearest positional genes, gene-based analyses, and differential gene expression data from RNA-seq, we identified putative causal genes (candidate genes with biological evidence) in QTL for mastitis resistance in Nordic Holstein cattle. The same strategy can be applied for other traits.

Electronic supplementary material

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

Genome-Wide Association Study for Susceptibility to and Recoverability From Mastitis in Danish Holstein Cows

Because mastitis is very frequent and unavoidable, adding recovery information into the analysis for genetic evaluation of mastitis is of great interest from economical and animal welfare point of view. Here we have performed genome-wide association studies (GWAS) to identify associated single nucleotide polymorphisms (SNPs) and investigate the genetic background not only for susceptibility to – but also for recoverability from mastitis. Somatic cell count records from 993 Danish Holstein cows genotyped for a total of 39378 autosomal SNP markers were used for the association analysis. Single SNP regression analysis was performed using the statistical software package DMU. Substitution effect of each SNP was tested with a t-test and a genome-wide significance level of P-value < 10^-4 was used to declare significant SNP-trait association. A number of significant SNP variants were identified for both traits. Many of the SNP variants associated either with susceptibility to – or recoverability from mastitis were located in or very near to genes that have been reported for their role in the immune system. Genes involved in lymphocyte developments (e.g., MAST3 and STAB2) and genes involved in macrophage recruitment and regulation of inflammations (PDGFD and PTX3) were suggested as possible causal genes for susceptibility to – and recoverability from mastitis, respectively. However, this is the first GWAS study for recoverability from mastitis and our results need to be validated. The findings in the current study are, therefore, a starting point for further investigations in identifying causal genetic variants or chromosomal regions for both susceptibility to – and recoverability from mastitis.

Genomic breeding values, SNP effects and gene identification for disease traits in cow training sets

Holstein Friesian cow training sets were created according to disease incidences. The different datasets were used to investigate the impact of random forest (RF) and genomic BLUP (GBLUP) methodology on genomic prediction accuracies. In addition, for further verifications of some specific scenarios, single-step genomic BLUP was applied. Disease traits included the overall trait categories of (i) claw disorders, (ii) clinical mastitis and (iii) infertility from 80 741 first lactation Holstein cows kept in 58 large-scale herds. A subset of 6744 cows was genotyped (50K SNP panel). Response variables for all scenarios were de-regressed proofs (DRPs) and pre-corrected phenotypes (PCPs). Initially, all sick cows were allocated to the testing set, and healthy cows represented the training set. For the ongoing cow allocation schemes, the number of sick cows in the training set increased stepwise by moving 10% of the sick cows from the testing to the training set in each step. The size of training and testing sets was kept constant by replacing the same number of cows in the testing set with (randomly selected) healthy cows from the training set. For both the RF and GBLUP methods, prediction accuracies were larger for DRPs compared to PCPs. For PCPs as a response variable, the largest prediction accuracies were observed when the disease incidences in training sets reflected the disease incidence in the whole population. A further increase in prediction accuracies for some selected cow allocation schemes (i.e. larger prediction accuracies compared to corresponding scenarios with RF or GBLUB) was achieved via single-step GBLUP applications. Correlations between genome-wide association study SNP effects and RF importance criteria for single SNPs were in a moderate range, from 0.42 to 0.57, when considering SNPs from all chromosomes or from specific chromosome segments. RF identified significant SNPs close to potential positional candidate genes: GAS1, GPAT3 and CYP2R1 for clinical mastitis; SPINK5 and SLC26A2 for laminitis; and FGF12 for endometritis.

Genome-wide association study of conformation and milk yield in mixed-breed dairy goats

Identification of genetic markers that affect economically important traits is of high value from a biological point of view, enabling the targeting of candidate genes and providing practical benefits for the industry such as wide-scale genomic selection. This study is one of the first to investigate the genetic background of economically important traits in dairy goats using the caprine 50K single nucleotide polymorphism (SNP) chip. The aim of the project was to perform a genome-wide association study for milk yield and conformation of udder, teat, and feet and legs. A total of 137,235 milk yield records on 4,563 goats each scored for 10 conformation traits were available. Out of these, 2,381 goats were genotyped with the Illumina Caprine 50K BeadChip (Illumina Inc., San Diego, CA). A range of pseudo-phenotypes were used including deregressed breeding values and pseudo-estimated breeding values. Genome-wide association studies were performed using the multi-locus mixed model (MLMM) algorithm implemented in SNP & Variation Suite v7.7.8 (Golden Helix Inc., Bozeman, MT). A genome-wise significant [-log₁₀(P-value) > 5.95] SNP for milk yield was identified on chromosome 19, with additional chromosome-wise significant (-log₁₀(P-value) > 4.46] SNP on chromosomes 4, 8, 14, and 29. Three genome-wise significant SNP for conformation of udder attachment, udder depth, and front legs were identified on chromosome 19, and chromosome-wise SNP were found on chromosomes 4, 5, 6, 10, 11, 12, 13, 14, 15, 16, 17, 18, 21, 23, and 27. The proportion of variance explained by the significant SNP was between 0.4 and 7.0% for milk yield and between 0.1 and 13.8% for conformation traits. This study is the first attempt to identify SNP associated with milk yield and conformation in dairy goats. Two genome-wise significant SNP for milk yield and 3 SNP for conformation of udder attachment, udder depth, and front legs were found. Our results suggest that conformation traits have a polygenic background because, for most of them, we did not identify any quantitative trait loci with major effect.

Population structure and genomic inbreeding in nine Swiss dairy cattle populations

BACKGROUND

Domestication, breed formation and intensive selection have resulted in divergent cattle breeds that likely exhibit their own genomic signatures. In this study, we used genotypes from 27,612 autosomal single nucleotide polymorphisms to characterize population structure based on 9214 sires representing nine Swiss dairy cattle populations: Brown Swiss (BS), Braunvieh (BV), Original Braunvieh (OB), Holstein (HO), Red Holstein (RH), Swiss Fleckvieh (SF), Simmental (SI), Eringer (ER) and Evolèner (EV). Genomic inbreeding (F ROH) and signatures of selection were determined by calculating runs of homozygosity (ROH). The results build the basis for a better understanding of the genetic development of Swiss dairy cattle populations and highlight differences between the original populations (i.e. OB, SI, ER and EV) and those that have become more popular in Switzerland as currently reflected by their larger populations (i.e. BS, BV, HO, RH and SF).

RESULTS

The levels of genetic diversity were highest and lowest in the SF and BS breeds, respectively. Based on F ST values, we conclude that, among all pairwise comparisons, BS and HO (0.156) differ more than the other pairs of populations. The original Swiss cattle populations OB, SI, ER, and EV are clearly genetically separated from the Swiss cattle populations that are now more common and represented by larger numbers of cows. Mean levels of F ROH ranged from 0.027 (ER) to 0.091 (BS). Three of the original Swiss cattle populations, ER (F ROH: 0.027), OB (F ROH: 0.029), and SI (F ROH: 0.039), showed low levels of genomic inbreeding, whereas it was much higher in EV (F ROH: 0.074). Private signatures of selection for the original Swiss cattle populations are reported for BTA4, 5, 11 and 26.

CONCLUSIONS

The low levels of genomic inbreeding observed in the original Swiss cattle populations ER, OB and SI compared to the other breeds are explained by a lesser use of artificial insemination and greater use of natural service. Natural service results in more sires having progeny at each generation and thus this breeding practice is likely the major reason for the remarkable levels of genetic diversity retained within these populations. The fact that the EV population is regionally restricted and its small census size of herd-book cows explain its high level of genomic inbreeding.

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.

Fine mapping of a QTL on bovine chromosome 6 using imputed full sequence data suggests a key role for the group-specific component (GC) gene in clinical mastitis and milk production

BACKGROUND

Clinical mastitis is an inflammation of the mammary gland and causes significant costs to dairy production. It is unfavourably genetically correlated to milk production, and, thus, knowledge of the mechanisms that underlie these traits would be valuable to improve both of them simultaneously through breeding. A quantitative trait locus (QTL) that affects both clinical mastitis and milk production has recently been fine-mapped to around 89 Mb on bovine chromosome 6 (BTA6), but identification of the gene that underlies this QTL was not possible due to the strong linkage disequilibrium between single nucleotide polymorphisms (SNPs) within this region. Our aim was to identify the gene and, if possible, the causal polymorphism(s) responsible for this QTL through association analysis of high-density SNPs and imputed full sequence data in combination with analyses of transcript and protein levels of the identified candidate gene.

RESULTS

Associations between SNPs and the studied traits were strongest for SNPs that were located within and immediately upstream of the group-specific component (GC) gene. This gene encodes the vitamin D-binding protein (DBP) and has multiple roles in immune defense and milk production. A 12-kb duplication that was identified downstream of this gene covered its last exon and segregated with the QTL allele that is associated with increased mastitis susceptibility and milk production. However, analyses of GC mRNA levels on the available samples revealed no differences in expression between animals having or lacking this duplication. Moreover, we detected no differences in the concentrations of DBP and its ligand vitamin D between the animals with different GC genotypes that were available for this study.

CONCLUSIONS

Our results suggest GC as the gene that underlies the QTL for clinical mastitis and milk production. However, since only healthy animals were sampled for transcription and expression analyses, we could not draw any final conclusion on the absence of quantitative differences between animals with different genotypes. Future studies should investigate GC RNA expression and protein levels in cows with different genotypes during an infection.

Bovine P-selectin mediates leukocyte adhesion and is highly polymorphic in dairy breeds

Bovine P-selectin (SELP) mediates leukocyte rolling and primes leukocyte adhesion to endothelium, both essential for leukocyte recruitment to an infection site. We investigated SELP-mediated adhesion between bovine peripheral blood leukocytes (PBLs) and cultured bovine aortic endothelial cells pre-activated with lipopolysaccharide (LPS). We examined gene polymorphism for bovine selectins SELP, l-selectin (SELL) and E-selectin (SELE) and compared their SNP frequency between five dairy breeds (Holstein, Friesian, Jersey, Ayrshire and Brown Swiss). LPS treatment caused a rapid (10min) and slower (4h) enhancement of PBL adhesion (P<0.01). Antibody blocking of SELP inhibited LPS induced cell adhesion. SELP was highly polymorphic, with 9 of the 13 SNPs in its exons, whereas only three synonymous SNPs in SELL and one in SELE. The resulting amino acid changes for the three missense SELP SNP were located in the lectin domain and in two consensus repeat (CR) regions, CR2 and CR5. The Val475Met variant locus in the CR4 and CR5 linking region was very close to a predicted N-acetyl-d-glucosamine glycosylation site, which is likely to influence SELP function. The AA genotype was under-represented, only being found in 1% of 373 heifers genotyped from the 5 breeds (P=0.056), suggesting that AA homozygous animals carrying the Val475Met substitution for SELP may have compromised development. Our study thus confirmed that SELP mediates the attachment of PBL to endothelium and provides novel evidence that its high polymorphism is likely to affect biological function. This may potentially influence leukocyte migration and fertility, both key to successful performance in dairy cows.

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

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

Deciphering Transcriptome and Complex Alternative Splicing Transcripts in Mammary Gland Tissues from Cows Naturally Infected with Staphylococcus aureus Mastitis

Alternative splicing (AS) contributes to the complexity of the mammalian proteome and plays an important role in diseases, including infectious diseases. The differential AS patterns of these transcript sequences between the healthy (HS3A) and mastitic (HS8A) cows naturally infected by Staphylococcus aureus were compared to understand the molecular mechanisms underlying mastitis resistance and susceptibility. In this study, using the Illumina paired-end RNA sequencing method, 1352 differentially expressed genes (DEGs) with higher than twofold changes were found in the HS3A and HS8A mammary gland tissues. Gene ontology and KEGG pathway analyses revealed that the cytokine–cytokine receptor interaction pathway is the most significantly enriched pathway. Approximately 16k annotated unigenes were respectively identified in two libraries, based on the bovine Bos taurus UMD3.1 sequence assembly and search. A total of 52.62% and 51.24% annotated unigenes were alternatively spliced in term of exon skipping, intron retention, alternative 5′ splicing and alternative 3ʹ splicing. Additionally, 1,317 AS unigenes were HS3A-specific, whereas 1,093 AS unigenes were HS8A-specific. Some immune-related genes, such as ITGB6, MYD88, ADA, ACKR1, and TNFRSF1B, and their potential relationships with mastitis were highlighted. From Chromosome 2, 4, 6, 7, 10, 13, 14, 17, and 20, 3.66% (HS3A) and 5.4% (HS8A) novel transcripts, which harbor known quantitative trait locus associated with clinical mastitis, were identified. Many DEGs in the healthy and mastitic mammary glands are involved in immune, defense, and inflammation responses. These DEGs, which exhibit diverse and specific splicing patterns and events, can endow dairy cattle with the potential complex genetic resistance against mastitis.

Detection of quantitative trait loci for maternal traits using high-density genotypes of Blonde d'Aquitaine beef cattle

Background

The genetic determinism of the calving and suckling performance of beef cows is little known whereas these maternal traits are of major economic importance in beef cattle production systems. This paper aims to identify QTL regions and candidate genes that affect maternal performance traits in the Blonde d’Aquitaine breed. Three calving performance traits were studied: the maternal effect on calving score from field data, the calving score and pelvic opening recorded in station for primiparous cows. Three other traits related to suckling performance were also analysed: the maternal effect on weaning weight from field data, milk yield and the udder swelling score recorded in station for primiparous cows. A total of 2,505 animals were genotyped from various chip densities and imputed in high density chips for 706,791 SNP. The number of genotyped animals with phenotypes ranged from 1,151 to 2,284, depending on the trait considered.

Results

QTL detections were performed using a Bayes C approach. Evidence for a QTL was based on Bayes Factor values. Putative candidate genes were proposed for the QTL with major evidence for one of the six traits and for the QTL shared by at least two of the three traits underlying either calving or suckling performance. Nine candidate genes were proposed for calving performance among the nine highlighted QTL regions. The neuroregulin gene on chromosome 27 was notably identified as a very likely candidate gene for maternal calving performance. As for suckling abilities, seven candidate genes were identified among the 15 highlighted QTL. In particular, the Group-Specific Component gene on chromosome 6, which encodes vitamin D binding protein, is likely to have a major effect on maternal weaning weight in the Blonde d’Aquitaine breed. This gene had already been linked to milk production and clinical mastitis in dairy cattle.

Conclusion

In the near future, these QTL findings and the preliminary proposals of candidate genes which act on the maternal performance of beef cows should help to identify putative causal mutations based on sequence data from different cattle breeds.

Electronic supplementary material

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

Genome-wide association for milk production and female fertility traits in Canadian dairy Holstein cattle

BACKGROUND

Genome-wide association studies (GWAS) are a powerful tool for detecting genomic regions explaining variation in phenotype. The objectives of the present study were to identify or refine the positions of genomic regions affecting milk production, milk components and fertility traits in Canadian Holstein cattle, and to use these positions to identify genes and pathways that may influence these traits.

RESULT

Several QTL regions were detected for milk production (MILK), fat production (FAT), protein production (PROT) and fat and protein deviation (FATD, PROTD respectively). The identified QTL regions for production traits (including milk production) support previous findings and some overlap with genes with known relevant biological functions identified in earlier studies such as DGAT1 and CPSF1. A significant region on chromosome 21 overlapping with the gene FAM181A and not previous linked to fertility in dairy cattle was identified for the calving to first service interval and days open. A functional enrichment analysis of the GWAS results yielded GO terms consistent with the specific phenotypes tested, for example GO terms GO:0007595 (lactation) and GO:0043627 (response to estrogen) for milk production (MILK), GO:0051057 (positive regulation of small GTPase mediated signal transduction) for fat production (FAT), GO:0040019 (positive regulation of embryonic development) for first service to calving interval (CTFS) and GO:0043268 (positive regulation of potassium ion transport) for days open (DO). In other cases the connection between the enriched GO terms and the traits were less clear, for example GO:0003279 (cardiac septum development) for FAT and GO:0030903 (notochord development) for DO trait.

CONCLUSION

The chromosomal regions and enriched pathways identified in this study confirm several previous findings and highlight new regions and pathways that may contribute to variation in production or fertility traits in dairy cattle.

Genome-wide association study in Chinese Holstein cows reveal two candidate genes for somatic cell score as an indicator for mastitis susceptibility

Backgrounds

Bovine mastitis is a typical inflammatory disease causing seriously economic loss. Genome-wide association study (GWAS) can be a powerful method to promote marker assistant selection of this kind of complex disease. The present study aimed to analyze and identify single nucleotide polymorphisms (SNPs) and candidate genes that associated with mastitis susceptibility traits in Chinese Holstein.

Results

Forty eight SNPs were identified significantly associated with mastitis resistance traits in Chinese Holstein cows, which are mainly located on the BTA 14. A total of 41 significant SNPs were linked to 31 annotated bovine genes. Gene Ontology and pathway enrichment revealed 5 genes involved in 32 pathways, in which, TRAPPC9 and ARHGAP39 genes participate cell differentiation and developmental pathway together. The six common genome-wide significant SNPs are found located within TRAPPC9 and flanking ARHGAP39 genes.

Conclusions

Our data identified the six SNPs significantly associated with SCS EBVs, which suggest that their linked two genes (TRAPPC9 and ARHGAP39) are novel candidate genes of mastitis susceptibility in Holsteins.

Electronic supplementary material

The online version of this article (doi:10.1186/s12863-015-0263-3) contains supplementary material, which is available to authorized users.

Association analysis for udder health based on SNP-panel and sequence data in Danish Holsteins

BACKGROUND

The sensitivity of genome-wide association studies for the detection of quantitative trait loci (QTL) depends on the density of markers examined and the statistical models used. This study compares the performance of three marker densities to refine six previously detected QTL regions for mastitis traits: 54 k markers of a medium-density SNP (single nucleotide polymorphism) chip (MD), imputed 777 k markers of a high-density SNP chip (HD), and imputed whole-genome sequencing data (SEQ). Each dataset contained data for 4496 Danish Holstein cattle. Comparisons were performed using a linear mixed model (LM) and a Bayesian variable selection model (BVS).

RESULTS

After quality control, 587, 7825, and 78 856 SNPs in the six targeted regions remained for MD, HD, and SEQ data, respectively. In general, the association patterns between SNPs and traits were similar for the three marker densities when tested using the same statistical model. With the LM model, 120 (MD), 967 (HD), and 7209 (SEQ) SNPs were significantly associated with mastitis, whereas with the BVS model, 43 (MD), 131 (HD), and 1052 (SEQ) significant SNPs (Bayes factor > 3.2) were observed. A total of 26 (MD), 75 (HD), and 465 (SEQ) significant SNPs were identified by both models. In addition, one, 16, and 33 QTL peaks for MD, HD, and SEQ data were detected according to the QTL intensity profile of SNP bins by post-analysis of the BVS model.

CONCLUSIONS

The power to detect significant associations increased with increasing marker density. The BVS model resulted in clearer boundaries between linked QTL than the LM model. Using SEQ data, the six targeted regions were refined to 33 candidate QTL regions for udder health. The comparison between these candidate QTL regions and known genes suggested that NPFFR2, SLC4A4, DCK, LIFR, and EDN3 may be considered as candidate genes for mastitis susceptibility.

Reactive oxygen species generation by bovine blood neutrophils with different CXCR1 (IL8RA) genotype following Interleukin-8 incubation

Background

Associations between polymorphisms in the bovine CXCR1 gene, encoding the chemokine (C-X-C motif) receptor 1 (IL8RA), and neutrophil traits and mastitis have been described. In the present study, blood neutrophils were isolated from 20 early lactating heifers with different CXCR1 genotype at position 735 or 980. The cells were incubated with different concentrations of recombinant bovine IL-8 (rbIL-8) for 2 or 6 h and stimulated with phorbol 12-myristate 13-acetate (PMA) or opsonized zymosan particles (OZP). Potential association between CXCR1 genotype and production of reactive oxygen species (ROS) was studied.

Results

Although on single nucleotide polymorphisms (SNPs) may potentially affect CXCR1 function, SNPs c.735C > G and c.980A > G showed no association with ROS production with or without incubation of rbIL-8. Neutrophils incubated with rbIL-8 for 2 or 6 h showed higher PMA- and lower OZP-induced ROS production compared to control without rbIL-8.

Conclusions

In the present study no association could be detected between superoxide production by isolated bovine neutrophils during early lactation and CXCR1 gene polymorphism. IL-8 showed to possess inhibitory effects on ROS generation in bovine neutrophils.

A genome-wide association study for clinical mastitis in first parity US Holstein cows using single-step approach and genomic matrix re-weighting procedure

Clinical mastitis (CM) is one of the health disorders with large impacts on dairy farming profitability and animal welfare. The objective of this study was to perform a genome-wide association study (GWAS) for CM in first-lactation Holstein. Producer-recorded mastitis event information for 103,585 first-lactation cows were used, together with genotype information on 1,361 bulls from the Illumina BovineSNP50 BeadChip. Single-step genomic-BLUP methodology was used to incorporate genomic data into a threshold-liability model. Association analysis confirmed that CM follows a highly polygenic mode of inheritance. However, 10-adjacent-SNP windows showed that regions on chromosomes 2, 14 and 20 have impacts on genetic variation for CM. Some of the genes located on chromosome 14 (LY6K, LY6D, LYNX1, LYPD2, SLURP1, PSCA) are part of the lymphocyte-antigen-6 complex (LY6) known for its neutrophil regulation function linked to the major histocompatibility complex. Other genes on chromosome 2 were also involved in regulating immune response (IFIH1, LY75, and DPP4), or are themselves regulated in the presence of specific pathogens (ITGB6, NR4A2). Other genes annotated on chromosome 20 are involved in mammary gland metabolism (GHR, OXCT1), antibody production and phagocytosis of bacterial cells (C6, C7, C9, C1QTNF3), tumor suppression (DAB2), involution of mammary epithelium (OSMR) and cytokine regulation (PRLR). DAVID enrichment analysis revealed 5 KEGG pathways. The JAK-STAT signaling pathway (cell proliferation and apoptosis) and the 'Cytokine-cytokine receptor interaction' (cytokine and interleukines response to infectious agents) are co-regulated and linked to the 'ABC transporters' pathway also found here. Gene network analysis performed using GeneMania revealed a co-expression network where 665 interactions existed among 145 of the genes reported above. Clinical mastitis is a complex trait and the different genes regulating immune response are known to be pathogen-specific. Despite the lack of information in this study, candidate QTL for CM were identified in the US Holstein population.

Eight SNVs in NF-κB pathway genes and their different performances between subclinical mastitis and mixed Chinese Holstein cows

The nuclear factor-kappa B (NF-κB) pathway proteins are key players in controlling both innate and adaptive immunity. However, the information on NF-κB pathway genes is very limited in mastitis resistance and milk production of Chinese Holstein cows. In this study, we examine the association of the NF-κB pathway gene variants with milk quality traits and somatic cell score (SCS) in Chinese Holstein cows. Eight single nucleotide variants (SNVs) were identified within the bovine NF-κB pathway genes, using DNA pooled sequencing, PCR-RFLP, and forced PCR-RFLP methods. These SNVs include SNV1: g. 536 C>T (exon 10 of Rel), SNV2: g. 94 G>A (exon 20 of p100), SNV3: g. 43 T>C (intron 6 of p105), SNV4: g. 2397 T>G (intron 9 of p105), SNV5: g. 382 G>C (intron 1 of IκBδ), SNV6: g. 21 C>T (exon 5 of IκBζ), SNV7: g. 272 G>A (intron 6 of IκBζ), and SNV8: g. 18 C>T (intron 10 of IκBζ). The association analysis in mixed Chinese Holstein population showed that SNV1 was significantly or highly significantly associated (P<0.01 and P<0.05) with fat rate, protein rate and SCS. Furthermore, the SNV1-CC (wild genotype) determined serine showed the significantly lower SCS and higher milk production traits compared to TT and TC. SNV2 was significantly associated (P<0.05) with SCS; SNV3 was significantly associated (P<0.05) with fat rate; and SNV4 was significantly associated (P<0.05) with fat rate and SCS. In 199 subclinical mastitis Chinese Holstein cows, the statistical results absolutely differed from the mixed Chinese Holstein individuals. Splice-site prediction by SplicePort showed that single nucleotide difference at eight SNVs results in the acceptor score and donor score changing obviously that may lead to alternative splicing. In brief, SNV1, SNV2, SNV3 and SNV4 could be useful genetic markers for mastitis resistance selection and breeding in Chinese Holstein cows. Furthermore, whether these SNVs lead to alternative splicing need further research.

Associations between CXCR1 polymorphisms and pathogen-specific incidence rate of clinical mastitis, test-day somatic cell count, and test-day milk yield

The CXCR1 gene plays an important role in the innate immunity of the bovine mammary gland. Associations between single nucleotide polymorphisms (SNP) CXCR1c.735C>G and c.980A>G and udder health have been identified before in small populations. A fluorescent multiprobe PCR assay was designed specifically and validated to genotype both SNP simultaneously in a reliable and cost-effective manner. In total, 3,106 cows from 50 commercial Flemish dairy herds were genotyped using this assay. Associations between genotype and detailed phenotypic data, including pathogen-specific incidence rate of clinical mastitis (IRCM), test-day somatic cell count, and test-day milk yield (MY) were analyzed. Staphylococcus aureus IRCM tended to associate with SNP c.735C>G. Cows with genotype c.735GG had lower Staph. aureus IRCM compared with cows with genotype c.735CC (rate ratio = 0.35, 95% confidence interval = 0.14–0.90). Additionally, a parity-specific association between Staph. aureus IRCM and SNP c.980A>G was detected. Heifers with genotype c.980GG had a lower Staph. aureus IRCM compared with heifers with genotype c.980AG (rate ratio = 0.15, 95% confidence interval = 0.04–0.56). Differences were less pronounced in multiparous cows. Associations between CXCR1 genotype and somatic cell count were not detected. However, MY was associated with SNP c.735C>G. Cows with genotype c.735GG out-produced cows with genotype c.735CC by 0.8 kg of milk/d. Results provide a basis for further research on the relation between CXCR1 polymorphism and pathogen-specific mastitis resistance and MY.

Bovine mastitis: frontiers in immunogenetics

Mastitis is one of the most prevalent and costly diseases in the dairy industry with losses attributable to reduced milk production, discarded milk, early culling, veterinary services, and labor costs. Typically, mastitis is an inflammation of the mammary gland most often, but not limited to, bacterial infection, and is characterized by the movement of leukocytes and serum proteins from the blood to the site of infection. It contributes to compromised milk quality and the potential spread of antimicrobial resistance if antibiotic treatment is not astutely applied. Despite the implementation of management practises and genetic selection approaches, bovine mastitis control continues to be inadequate. However, some novel genetic strategies have recently been demonstrated to reduce mastitis incidence by taking advantage of a cow's natural ability to make appropriate immune responses against invading pathogens. Specifically, dairy cattle with enhanced and balanced immune responses have a lower occurrence of disease, including mastitis, and they can be identified and selected for using the high immune response (HIR) technology. Enhanced immune responsiveness is also associated with improved response to vaccination, increased milk, and colostrum quality. Since immunity is an important fitness trait, beneficial associations with longevity and reproduction are also often noted. This review highlights the genetic regulation of the bovine immune system and its vital contributions to disease resistance. Genetic selection approaches currently used in the dairy industry to reduce the incidence of disease are reviewed, including the HIR technology, genomics to improve disease resistance or immune response, as well as the Immunity(+)™ sire line. Improving the overall immune responsiveness of cattle is expected to provide superior disease resistance, increasing animal welfare and food quality while maintaining favorable production levels to feed a growing population.

Genome-wide association study using high-density single nucleotide polymorphism arrays and whole-genome sequences for clinical mastitis traits in dairy cattle

Mastitis is a mammary disease that frequently affects dairy cattle. Despite considerable research on the development of effective prevention and treatment strategies, mastitis continues to be a significant issue in bovine veterinary medicine. To identify major genes that affect mastitis in dairy cattle, 6 chromosomal regions on Bos taurus autosome (BTA) 6, 13, 16, 19, and 20 were selected from a genome scan for 9 mastitis phenotypes using imputed high-density single nucleotide polymorphism arrays. Association analyses using sequence-level variants for the 6 targeted regions were carried out to map causal variants using whole-genome sequence data from 3 breeds. The quantitative trait loci (QTL) discovery population comprised 4,992 progeny-tested Holstein bulls, and QTL were confirmed in 4,442 Nordic Red and 1,126 Jersey cattle. The targeted regions were imputed to the sequence level. The highest association signal for clinical mastitis was observed on BTA 6 at 88.97 Mb in Holstein cattle and was confirmed in Nordic Red cattle. The peak association region on BTA 6 contained 2 genes: vitamin D-binding protein precursor (GC) and neuropeptide FF receptor 2 (NPFFR2), which, based on known biological functions, are good candidates for affecting mastitis. However, strong linkage disequilibrium in this region prevented conclusive determination of the causal gene. A different QTL on BTA 6 located at 88.32 Mb in Holstein cattle affected mastitis. In addition, QTL on BTA 13 and 19 were confirmed to segregate in Nordic Red cattle and QTL on BTA 16 and 20 were confirmed in Jersey cattle. Although several candidate genes were identified in these targeted regions, it was not possible to identify a gene or polymorphism as the causal factor for any of these regions.

The early phase transcriptome of bovine monocyte-derived macrophages infected with Staphylococcus aureus in vitro

Background

In the mammary gland, local recruitment and action of macrophages is a key immunological defence mechanism against infection. Macrophages are members of the innate immune system, serve as the first line of the defence against invading pathogens and are critical effectors and regulators of inflammation. We have examined the early phase response of bovine macrophages to infection with live Staphylococcus aureus. Genome-wide transcript profiling of blood monocyte-derived macrophages from six Norwegian Red heifers infected with live S. aureus for 2 and 6 hours in vitro was performed.

Results

About 420 of the 17 000 genes on the ARK-Genomics bovine cDNA array were differentially regulated at 6 hours post infection. Approximately 70% of the responding genes had a known identity (Entrez Gene ID) and were used in the identification of overrepresented pathways and biological functions in the dataset.

Analysis of a subset of differentially regulated genes (List eQG) obtained by comparison with data from genome-wide association mapping in Norwegian Red cattle identified anti-inflammatory cytokines interleukin 4 and interleukin 13 as putative expression quantitative trait loci, suggesting that S. aureus infection triggers alternative activation of macrophages. Moreover, several classical activation pathways were found, mainly cellular immune response and cytokine signaling pathways, i.e. triggering receptor expressed on myeloid cells 1 (TREM1) and nucleotide-binding and oligomerization domain-like receptor (NLR) pathways. Tumor necrosis factor receptor superfamily member 5 (CD40 ligand) was identified as an upstream regulator which points toward CD40 likely acting as a co-stimulatory receptor during Toll-like receptor 2(TLR2)-mediated inflammatory response of bovine macrophages to S. aureus infection. Furthermore, peptidoglycan was identified as an upstream regulator in the List eQG, which indicates that this bacterial cell-wall component might be pivotal in macrophage intracellular bacterial recognition during early inflammation.

Conclusions

Here we have shown that in vitro infection of bovine macrophages with live S. aureus induced both alternative and classical activation pathways. Alternative activation of macrophages may be a mechanism contributing to intracellular persistence of S. aureus in the course of inflammation such as during mastitis in dairy cattle.

Effects of Bos taurus autosome 9-located quantitative trait loci haplotypes on the disease phenotypes of dairy cows with experimentally induced Escherichia coli mastitis

Several quantitative trait loci (QTL) affecting mastitis incidence and mastitis-related traits such as somatic cell score exist in dairy cows. Previously, QTL haplotypes associated with susceptibility to Escherichia coli mastitis in Nordic Holstein-Friesian (HF) cows were identified on Bos taurus autosome 9. In the present study, we induced experimental E. coli mastitis in Danish HF cows to investigate the effect of 2 E. coli mastitis-associated QTL haplotypes on the cows' disease phenotypes and recovery in early lactation. Thirty-two cows were divided in 2 groups bearing haplotypes with either low (HL) or high (HH) susceptibility to E. coli. In addition, biopsies (liver and udder) were collected from half of the cows (n=16), resulting in a 2 × 2 factorial design, with haplotype being one factor (HL vs. HH) and biopsy being the other factor (biopsies vs. no biopsies). Each cow was inoculated with a low E. coli dose (20 to 40 cfu) in one front quarter at time 0 h. Liver biopsies were collected at -144, 12, 24, and 192 h; udder biopsies were collected at 24h and 192 h post-E. coli inoculation. The clinical parameters: feed intake, milk yield, body temperature, heart rate, respiration rate, rumen motility; and the paraclinical parameters: bacterial counts, somatic cell count (SCC), and milk amyloid A levels in milk; and white blood cell count, polymorphonuclear neutrophilic leukocyte (PMNL) count, and serum amyloid A levels in blood were recorded at different time points post-E. coli inoculation. Escherichia coli inoculation changed the clinical and paraclinical parameters in all cows except one that was not infected. Clinically, the HH group tended to have higher body temperature and heart rate than the HL group did. Paraclinically, the HL group had faster PMNL recruitment and SCC recovery than the HH group did. However, we also found interactions between the effects of haplotype and biopsy for body temperature, heart rate, and PMNL. In conclusion, when challenged with E. coli mastitis, HF cows with the specific Bos taurus autosome 9-located QTL haplotypes were associated with differences in leukocyte kinetics, with low-susceptibility cows having faster blood PMNL recruitment and SCC recovery and a tendency for a milder clinical response than the high-susceptibility cows did.

Pathogen-group specific association between CXCR1 polymorphisms and subclinical mastitis in dairy heifers

The chemokine (C-X-C motif) receptor 1 (CXCR1) gene encodes the homonymous receptor for interleukin 8 (IL8) on polymorphonuclear neutrophilic leucocytes (PMNL). Binding causes migration from blood to milk, activation and prolonged survival of PMNL, a crucial process in the innate immune defence of the bovine mammary gland against invading mastitis-causing pathogens. The main objective of this study was to screen the entire coding region of the CXCR1 gene for polymorphisms and to analyse their association with udder health of dairy heifers. One-hundred-and-forty Belgian Holstein heifers originating from 20 commercial dairy farms were genotyped by DNA sequencing. Detailed phenotypic data on udder health was available including quarter bacteriological culture results and somatic cell count (SCC) in early lactation and composite milk SCC during first lactation. In total, 16 polymorphisms (including 8 missense mutations) were detected. Polymorphism c.980A>G was associated with pathogen-group specific IMI: heifers with genotype AG were less likely to have an IMI due to major mastitis pathogens compared with heifers with genotype GG but did not have less IMI by coagulase-negative staphylococci, so-called minor pathogens. CXCR1 genotype was neither associated with quarter SCC in early lactation nor with composite SCC during lactation. Although mastitis susceptibility is influenced by many factors, some genetic polymorphisms potentially have major effects on udder health of heifers, as was shown here. These results trigger us to further study the relationship between CXCR1 polymorphisms and mastitis susceptibility in both observational and experimental trials.

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

Genomic associations with somatic cell score in first-lactation Holstein cows

This genome-wide association study aimed to identify loci associated with lactation-average somatic cell score (LASCS) and the standard deviation of test-day somatic cell score (SCS-SD). It is one of the first studies to combine detailed phenotypic and genotypic cow data from research dairy herds located in different countries. The combined data set contained up to 52 individual test-days per lactation and thereby aimed to capture temporary increases in somatic cell score associated with infection. Phenotypic data for analysis consisted of 46,882 test-day records on 1,484 cows, and genotypic data consisted of 37,590 single nucleotide polymorphisms (SNP). Using an animal model, the associations between each individual SNP and the phenotypic data were estimated. To account for the risk of false positives, a false discovery rate threshold of 0.20 was set. The analyses showed that LASCS was significantly associated with a SNP on Bos taurus autosome (BTA) 4 and a SNP on BTA18. Likewise, SCS-SD was associated with this SNP on BTA18. In addition, SCS-SD significantly associated with a SNP on BTA6. Relatively few associations were found, suggesting that LASCS and SCS-SD are controlled by multiple loci distributed across the genome, each with a relatively small effect. Increased knowledge on genetic regulation of LASCS and SCS-SD may aid in identification of genes that play a role in mastitis resistance. Such knowledge helps us understand the genetic mechanisms leading to mastitis and in discovery of targets for mastitis therapeutics.

Combining quantitative trait loci and heterogeneous microarray data analyses reveals putative candidate pathways affecting mastitis in cattle

Mastitis is a frequent disease and considerable problem for the global dairy industry. Identification of solutions leading to the development of new control strategies is therefore of high importance. In this study, we have integrated genomic data from genome-wide association mapping in cattle with transcriptomic data from microarray studies of several mastitis pathogens and host species in vitro and in vivo. To identify significant candidate pathways directly and indirectly involved in the immune response to mastitis, ingenuity pathway analysis (ipa) and database for annotation, visualization and integrated discovery bioinformatic (david) were applied. Several candidate pathways were found. Of great interest are IL-17 and IL-8 signalling pathways, responsible for the recruitment and migration of inflammatory cells into tissue during inflammation and infection. These results may emphasize further functional studies for identification of factors contributing to resistance to mastitis pathogens in cattle.

Gene polymorphisms: the keys for marker assisted selection and unraveling core regulatory pathways for mastitis resistance

One of the most frequent mammary diseases impacting lactating animals is mastitis, an inflammation of the mammary gland most commonly caused by bacterial infection. The severity of mastitis is greatly influenced by the invading organism and the subsequent immune response which must recognize the foreign organism, recruit immune cells, eliminate the invading pathogen, and resolve the inflammatory response. The speed, strength, and duration of this response and subsequent disease susceptibility are critically tied to the genetic background of an animal. However, the genetic contribution has been difficult to identify due to the complex interactions that must occur for effective disease resistance. Recent studies have utilized polymorphisms to better define the genes and chromosomal regions that contribute to mastitis resistance. This review will examine these studies with primary emphasis in bovine systems, as the most work regarding mastitis has been conducted in this species.