Genetic dissection of milk yield traits and mastitis resistance quantitative trait loci on chromosome 20 in dairy cattle.

Intense selection to increase milk yield has had negative consequences for mastitis incidence in dairy cattle. Due to low heritability of mastitis resistance and an unfavorable genetic correlation with milk yield, a reduction in mastitis through traditional breeding has been difficult to achieve. Here, we examined quantitative trait loci (QTL) that segregate for clinical mastitis and milk yield on Bos taurus autosome 20 (BTA20) to determine whether both traits are affected by a single polymorphism (pleiotropy) or by multiple closely linked polymorphisms. In the latter but not the former situation, undesirable genetic correlation could potentially be broken by selecting animals that have favorable variants for both traits. First, we performed a within-breed association study using a haplotype-based method in Danish Holstein cattle (HOL). Next, we analyzed Nordic Red dairy cattle (RDC) and Danish Jersey cattle (JER) with the goal of determining whether these QTL identified in Holsteins were segregating across breeds. Genotypes for 12,566 animals (5,966 HOL, 5,458 RDC, and 1,142 JER) were determined by using the Illumina Bovine SNP50 BeadChip (50K; Illumina, San Diego, CA), which identifies 1,568 single nucleotide polymorphisms on BTA20. Data were combined, phased, and clustered into haplotype states, followed by within- and across-breed haplotype-based association analyses using a linear mixed model. Association signals for both clinical mastitis and milk yield peaked in the 26- to 40-Mb region on BTA20 in HOL. Single-variant association analyses were carried out in the QTL region using whole sequence level variants imputed from references of 2,036 HD genotypes (BovineHD BeadChip; Illumina) and 242 whole-genome sequences. The milk QTL were also segregating in RDC and JER on the BTA20-targeted region; however, an indication of differences in the causal factor(s) was observed across breeds. A previously reported F279Y mutation (rs385640152) within the growth hormone receptor gene showed strong association with milk, fat, and protein yields. In HOL, the highest peaks for milk yield and susceptibility to mastitis were separated by over 3.5 Mb (3.8 Mb by haplotype analysis, 3.6 Mb by single nucleotide polymorphism analysis), suggesting separate genetic variants for the traits. Further analysis yielded 2 candidate mutations for the mastitis QTL, at 33,642,072 bp (rs378947583) in an intronic region of the caspase recruitment domain protein 6 gene and 35,969,994 bp (rs133596506) in an intronic region of the leukemia-inhibitory factor receptor gene. These findings suggest that it may be possible to separate these beneficial and detrimental genetic factors through targeted selective breeding.

Kadri NK ，Guldbrandtsen B ，Lund MS ，Sahana G ... -  《-》

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

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

Cai Z ，Dusza M ，Guldbrandtsen B ，Lund MS ，Sahana G ... -  《-》

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.

Sahana G ，Guldbrandtsen B ，Thomsen B ，Holm LE ，Panitz F ，Brøndum RF ，Bendixen C ，Lund MS ... -  《-》

Molecular characterization of a long range haplotype affecting protein yield and mastitis susceptibility in Norwegian Red cattle.

Previous fine mapping studies in Norwegian Red cattle (NRC) in the region 86-90.4 Mb on Bos taurus chromosome 6 (BTA6) has revealed a quantitative trait locus (QTL) for protein yield (PY) around 88 Mb and a QTL for clinical mastitis (CM) around 90 Mb. The close proximity of these QTLs may partly explain the unfavorable genetic correlation between these two traits in NRC. A long range haplotype covering this region was introduced into the NRC population through the importation of a Holstein-Friesian bull (1606 Frasse) from Sweden in the 1970s. It has been suggested that this haplotype has a favorable effect on milk protein content but an unfavorable effect on mastitis susceptibility. Selective breeding for milk production traits is likely to have increased the frequency of this haplotype in the NRC population. Association mapping for PY and CM in NRC was performed using genotypes from 556 SNPs throughout the region 86-97 Mb on BTA6 and daughter-yield-deviations (DYDs) from 2601 bulls made available from the Norwegian dairy herd recording system. Highest test scores for PY were found for single-nucleotide polymorphisms (SNPs) within and surrounding the genes CSN2 and CSN1S2, coding for the β-casein and α(S2)-casein proteins. High coverage re-sequencing by high throughput sequencing technology enabled molecular characterization of a long range haplotype from 1606 Frasse encompassing these two genes. Haplotype analysis of a large number of descendants from this bull indicated that the haplotype was not markedly disrupted by recombination in this region. The haplotype was associated with both increased milk protein content and increased susceptibility to mastitis, which might explain parts of the observed genetic correlation between PY and CM in NRC. Plausible causal polymorphisms affecting PY were detected in the promoter region and in the 5'-flanking UTR of CSN1S2. These polymorphisms could affect transcription or translation of CSN1S2 and thereby affect the amount of α(S2)-casein in milk. Highest test scores for CM were found in the region 89-91 Mb on BTA6, very close to a cluster of genes coding for CXC chemokines. Expression levels of some of these CXC chemokines have previously been shown to increase in bovine mammary gland cell lines after exposure to bacterial cell wall components. Molecular characterization of the long range haplotype from the Holstein-Friesian bull 1606 Frasse, imported into NRC in the 1970s, revealed polymorphisms that could affect transcription or translation of the casein gene CSN1S2. Sires with this haplotype had daughters with significantly elevated milk protein content and selection for milk production traits is likely to have increased the frequency of this haplotype in the NRC population. The haplotype was also associated with increased mastitis susceptibility, which might explain parts of the genetic correlation between PY and CM in NRC.

Sodeland M ，Grove H ，Kent M ，Taylor S ，Svendsen M ，Hayes BJ ，Lien S ... -  《BMC GENETICS》

Association analysis for udder index and milking speed with imputed whole-genome sequence variants in Nordic Holstein cattle.

Genome-wide association testing facilitates the identification of genetic variants associated with complex traits. Mapping genes that promote genetic resistance to mastitis could reduce the cost of antibiotic use and enhance animal welfare and milk production by improving outcomes of breeding for udder health. Using imputed whole-genome sequence variants, we carried out association studies for 2 traits related to udder health, udder index, and milking speed in Nordic Holstein cattle. A total of 4,921 bulls genotyped with the BovineSNP50 BeadChip array were imputed to high-density genotypes (Illumina BovineHD BeadChip, Illumina, San Diego, CA) and, subsequently, to whole-genome sequence variants. An association analysis was carried out using a linear mixed model. Phenotypes used in the association analyses were deregressed breeding values. Multitrait meta-analysis was carried out for these 2 traits. We identified 10 and 8 chromosomes harboring markers that were significantly associated with udder index and milking speed, respectively. Strongest association signals were observed on chromosome 20 for udder index and chromosome 19 for milking speed. Multitrait meta-analysis identified 13 chromosomes harboring associated markers for the combination of udder index and milking speed. The associated region on chromosome 20 overlapped with earlier reported quantitative trait loci for similar traits in other cattle populations. Moreover, this region was located close to the FYB gene, which is involved in platelet activation and controls IL-2 expression; FYB is a strong candidate gene for udder health and worthy of further investigation.

Jardim JG ，Guldbrandtsen B ，Lund MS ，Sahana G ... -  《-》