Quantitative trait loci mapping in dairy cattle: review and meta-analysis

From an extensive review of public domain information on dairy cattle quantitative trait loci (QTL), we have prepared a draft online QTL map for dairy production traits. Most publications (45 out of 55 reviewed) reported QTL for the major milk production traits (milk, fat and protein yield, and fat and protein concentration (%)) and somatic cell score. Relatively few QTL studies have been reported for more complex traits such as mastitis, fertility and health. The collated QTL map shows some chromosomal regions with a high density of QTL, as well as a substantial number of QTL at single chromosomal locations. To extract the most information from these published records, a meta-analysis was conducted to obtain consensus on QTL location and allelic substitution effect of these QTL. This required modification and development of statistical methodologies. The meta-analysis indicated a number of consensus regions, the most striking being two distinct regions affecting milk yield on chromosome 6 at 49 cM and 87 cM explaining 4.2 and 3.6 percent of the genetic variance of milk yield, respectively. The first of these regions (near marker BM143) affects five separate milk production traits (protein yield, protein percent, fat yield, fat percent, as well as milk yield).     

Analysis of chromosomal localization of loci controlling milk production traits in cattle

Analysis of the pattern of the chromosomal localization of quantitative trait loci (QTLs) is necessary for comprehensively understanding their functions. The chromosomal localization of QTLs controlling milk production traits has been studied in cattle chromosomes. The distribution of QTLs between chromosomes has proved to be binomial. Their distribution along each chromosome was, in general, uniform, except for the QTLs controlling the somatic cell score (SCS), which tended towards telomeric location. However, there are chromosomes either enriched with or particularly poor in QTLs. The QTL distribution patters are the most similar for the milk yield (M) and milk protein yield (P) and for milk fat yield (F) and milk fat content (%F). The pattern of the SCS QTLs stands out among those of other QTLs. The distance between the QTLs of contrasting traits is the shortest for M and P QTLs, longer for M and milk protein content (%P) QTLs, and still longer for M and %F QTLs, which may be explained by QTL pleiotropy, a common phenomenon in cattle.     

A genome map of divergent artificial selection between Bos taurus dairy cattle and Bos taurus beef cattle

A number of cattle breeds have become highly specialized for milk or beef production, following strong artificial selection for these traits. In this paper, we compare allele frequencies from 9323 single nucleotide polymorphism (SNP) markers genotyped in dairy and beef cattle breeds averaged in sliding windows across the genome, with the aim of identifying divergently selected regions of the genome between the production types. The value of the method for identifying selection signatures was validated by four sources of evidence. First, differences in allele frequencies between dairy and beef cattle at individual SNPs were correlated with the effects of those SNPs on production traits. Secondly, large differences in allele frequencies generally occurred in the same location for two independent data sets (correlation 0.45) between sliding window averages. Thirdly, the largest differences in sliding window average difference in allele frequencies were found on chromosome 20 in the region of the growth hormone receptor gene, which carries a mutation known to have an effect on milk production traits in a number of dairy populations. Finally, for the chromosome tested, the location of selection signatures between dairy and beef cattle was correlated with the location of selection signatures within dairy cattle.     

An expanded comparative map of bovine chromosome 27 targeting dairy form QTL regions

At present, the density of genes on the bovine maps is extremely limited and current resolution of the human-bovine comparative map is insufficient for selection of candidate genes controlling many economic traits of interest in dairy cattle. This study describes the chromosomal mapping of 10 selected gene-associated markers to bovine linkage and radiation hybrid maps to improve the breakpoint resolution in the human-bovine comparative map near two previously identified quantitative trait loci for the linear type trait, dairy form. Two regions of conserved synteny not previously described are reported between the telomeric region of bovine chromosome 27 (BTA27) and human chromosome 3 (HSA3) p24 region and between the HSA4q34.1 region and BTA8. These data increase the number of genes positioned on the bovine gene maps, refine the human-bovine comparative map, and should improve the efficiency of candidate gene selection for the dairy form trait in cattle.     

Quantitative trait loci underlying milk production traits in sheep

Improvement of milk production traits in dairy sheep is required to increase the competitiveness of the industry and to maintain the production of high quality cheese in regions of Mediterranean countries with less favourable conditions. Additional improvement over classical selection could be reached if genes with significant effects on the relevant traits were specifically targeted by selection. However, so far, few studies have been undertaken to detect quantitative trait loci (QTL) in dairy sheep. In this study, we present a complete genome scan performed in a commercial population of Spanish Churra sheep to identify chromosomal regions associated with phenotypic variation observed in milk production traits. Eleven half-sib families, including a total of 1213 ewes, were analysed following a daughter design. Genome-wise multi-marker regression analysis revealed a genome-wise significant QTL for milk protein percentage on chromosome 3. Eight other regions, localized on chromosomes 1, 2, 20, 23 and 25, showed suggestive significant linkage associations with some of the analysed traits. To our knowledge, this study represents the first complete genome scan for milk production traits reported in dairy sheep. The experiment described here shows that analysis of commercial dairy sheep populations has the potential to increase our understanding of the genetic determinants of complex production-related traits.     

The origin of selection signatures on bovine chromosome 6

The extent and pattern of linkage disequilibrium (LD) between closely spaced markers contain information about population history, including past population size and selection history. Selection signatures can be identified by comparing the LD surrounding a putative selected allele at a locus to the putative non-selected allele. In livestock populations, locations of selection signatures identified in this way should be correlated with QTL affecting production traits, as the populations have been under strong artificial selection for these traits. We used a dense SNP map of bovine chromosome 6 to characterize the pattern of LD on this chromosome in Norwegian Red cattle, a breed which has been strongly selected for milk production. The pattern of LD was generally consistent with strong selection in regions containing QTL affecting milk production traits, including a strong selection signature in a region containing a mutation known to affect milk production. The results demonstrate that in livestock populations, the origin of selection signatures will often be QTL for livestock production traits, and illustrate the value of selection signatures in uncovering new mutations with potential effects on quantitative traits.     

Genome wide association studies for milk production traits in Chinese Holstein population

Genome-wide association studies (GWAS) based on high throughput SNP genotyping technologies open a broad avenue for exploring genes associated with milk production traits in dairy cattle. Motivated by pinpointing novel quantitative trait nucleotide (QTN) across Bos Taurus genome, the present study is to perform GWAS to identify genes affecting milk production traits using current state-of-the-art SNP genotyping technology, i.e., the Illumina BovineSNP50 BeadChip. In the analyses, the five most commonly evaluated milk production traits are involved, including milk yield (MY), milk fat yield (FY), milk protein yield (PY), milk fat percentage (FP) and milk protein percentage (PP). Estimated breeding values (EBVs) of 2,093 daughters from 14 paternal half-sib families are considered as phenotypes within the framework of a daughter design. Association tests between each trait and the 54K SNPs are achieved via two different analysis approaches, a paternal transmission disequilibrium test (TDT)-based approach (L1-TDT) and a mixed model based regression analysis (MMRA). In total, 105 SNPs were detected to be significantly associated genome-wise with one or multiple milk production traits. Of the 105 SNPs, 38 were commonly detected by both methods, while four and 63 were solely detected by L1-TDT and MMRA, respectively. The majority (86 out of 105) of the significant SNPs is located within the reported QTL regions and some are within or close to the reported candidate genes. In particular, two SNPs, ARS-BFGL-NGS-4939 and BFGL-NGS-118998, are located close to the DGAT1 gene (160bp apart) and within the GHR gene, respectively. Our findings herein not only provide confirmatory evidences for previously findings, but also explore a suite of novel SNPs associated with milk production traits, and thus form a solid basis for eventually unraveling the causal mutations for milk production traits in dairy cattle.     

Effects of marked chromosome sections on milk performance in cattle

Summary The investigations of paternal half sibs start with the assumption that the transfer of an allele from a father to an offspring also indicates the inheritance of a distinct section of two homologous chromosomes from the father concerned. With the help of 20 gene systems, the transfers of single chromosome sections were marked and tested with regard to influences on milk performance traits. 1,457 German Friesian cattle, registered as daughters of three sires, were used. Some of the chromosome sections showed significant effects on the traits considered. Since especially those chromosomes which bear genes for milk proteins were involved, it was assumed that groups of linked loci influence the genetic variance of milk production. Possibilities for applying the results to the practical breeding situation and their significance are discussed.     

Genome-wide identification of runs of homozygosity islands and associated genes in local dairy cattle breeds

Runs of homozygosity (ROH) are widely used as predictors of whole-genome inbreeding levels in cattle. They identify regions that have an unfavorable effect on a phenotype when homozygous, but also identify the genes associated with traits of economic interest present in these regions. Here, the distribution of ROH islands and enriched genes within these regions in four dairy cattle breeds were investigated. Cinisara (71), Modicana (72), Reggiana (168) and Italian Holstein (96) individuals were genotyped using the 50K v2 Illumina BeadChip. The genomic regions most commonly associated with ROHs were identified by selecting the top 1% of the single nucleotide polymorphisms (SNPs) most commonly observed in the ROH of each breed. In total, 11 genomic regions were identified in Cinisara and Italian Holstein, and eight in Modicana and Reggiana, indicating an increased ROH frequency level. Generally, ROH islands differed between breeds. The most homozygous region (>45% of individuals with ROH) was found in Modicana on chromosome 6 within a quantitative trail locus affecting milk fat and protein concentrations. We identified between 126 and 347 genes within ROH islands, which are involved in multiple signaling and signal transduction pathways in a wide variety of biological processes. The gene ontology enrichment provided information on possible molecular functions, biological processes and cellular components under selection related to milk production, reproduction, immune response and resistance/susceptibility to infection and diseases. Thus, scanning the genome for ROH could be an alternative strategy to detect genomic regions and genes related to important economic traits.     

Application of a canonical transformation to detection of quantitative trait loci with the aid of genetic markers in a multi-trait experiment

Effects of individual quantitative trait loci (QTLs) can be isolated with the aid of linked genetic markers. Most studies have analyzed each marker or pair of linked markers separately for each trait included in the analysis. Thus, the number of contrasts tested can be quite large. The experimentwise type-I error can be readily derived from the nominal type-I error if all contrasts are statistically independent, but different traits are generally correlated. A new set of uncorrelated traits can be derived by application of a canonical transformation. The total number of effective traits will generally be less than the original set. An example is presented for DNA microsatellite D21S4, which is used as a marker for milk production traits of Israeli dairy cattle. This locus had significant effects on milk and protein production but not on fat. It had a significant effect on only one of the canonical variables that was highly correlated with both milk and protein, and this variable explained 82% of the total variance. Thus, it can be concluded that a single QTL is affecting both traits. The effects on the original traits could be derived by a reverse transformation of the effects on the canonical variable.     

Estimation of (co)variances for genomic regions of flexible sizes: application to complex infectious udder diseases in dairy cattle

Background

Multi-trait genomic models in a Bayesian context can be used to estimate genomic (co)variances, either for a complete genome or for genomic regions (e.g. per chromosome) for the purpose of multi-trait genomic selection or to gain further insight into the genomic architecture of related traits such as mammary disease traits in dairy cattle.

Methods

Data on progeny means of six traits related to mastitis resistance in dairy cattle (general mastitis resistance and five pathogen-specific mastitis resistance traits) were analyzed using a bivariate Bayesian SNP-based genomic model with a common prior distribution for the marker allele substitution effects and estimation of the hyperparameters in this prior distribution from the progeny means data. From the Markov chain Monte Carlo samples of the allele substitution effects, genomic (co)variances were calculated on a whole-genome level, per chromosome, and in regions of 100 SNP on a chromosome.

Results

Genomic proportions of the total variance differed between traits. Genomic correlations were lower than pedigree-based genetic correlations and they were highest between general mastitis and pathogen-specific traits because of the part-whole relationship between these traits. The chromosome-wise genomic proportions of the total variance differed between traits, with some chromosomes explaining higher or lower values than expected in relation to chromosome size. Few chromosomes showed pleiotropic effects and only chromosome 19 had a clear effect on all traits, indicating the presence of QTL with a general effect on mastitis resistance. The region-wise patterns of genomic variances differed between traits. Peaks indicating QTL were identified but were not very distinctive because a common prior for the marker effects was used. There was a clear difference in the region-wise patterns of genomic correlation among combinations of traits, with distinctive peaks indicating the presence of pleiotropic QTL.

Conclusions

The results show that it is possible to estimate, genome-wide and region-wise genomic (co)variances of mastitis resistance traits in dairy cattle using multivariate genomic models.     

High-density mapping and comparative analysis of agronomically important traits on wheat chromosome 3A

Bread wheat chromosome 3A has been shown to contain genes/QTLs controlling grain yield and other agronomic traits. The objectives of this study were to generate high-density physical and genetic-linkage maps of wheat homoeologous group 3 chromosomes and reveal the physical locations of genes/QTLs controlling yield and its component traits, as well as agronomic traits, to obtain a precise estimate of recombination for the corresponding regions and to enrich the QTL-containing regions with markers. Physical mapping was accomplished by 179 DNA markers mostly representing expressed genes using 41 single-break deletion lines. Polymorphism survey of cultivars Cheyenne (CNN) and Wichita (WI), and a substitution line of CNN carrying chromosome 3A from WI [CNN(WI3A)], with 142 RFLP probes and 55 SSR markers revealed that the extent of polymorphism is different among various group 3 chromosomal regions as well as among the homoeologs. A genetic-linkage map for chromosome 3A was developed by mapping 17 QTLs for seven agronomic traits relative to 26 RFLP and 15 SSR chromosome 3A-specific markers on 95 single-chromosome recombinant inbred lines. Comparison of the physical maps with the 3A genetic-linkage map localized the QTLs to gene-containing regions and accounted for only about 36% of the chromosome. Two chromosomal regions containing 9 of the 17 QTLs encompassed less than 10% of chromosome 3A but accounted for almost all of the arm recombination. To identify rice chromosomal regions corresponding to the particular QTL-containing wheat regions, 650 physically mapped wheat group 3 sequences were compared with rice genomic sequences. At an E value of E < or = 10(-5), 82% of the wheat group 3 sequences identified rice homologs, of which 54% were on rice chromosome 1. The rice chromosome 1 region collinear with the two wheat regions that contained 9 QTLs was about 6.5 Mb.     

Mapping of fertility traits in Finnish Ayrshire by genome-wide association analysis

A whole-genome scan using single marker association was used to detect chromosome regions associated with seven female fertility traits in Finnish Ayrshire dairy cattle. The phenotypic data consisted of de-regressed estimated breeding values for 340 bulls which were estimated using a single trait model. Genotypes were obtained with the Illumina BovineSNP50 panel and a total of 35 630 informative, high-quality single nucleotide polymorphism (SNP) markers were used. The association analysis was performed using a mixed-model approach which fitted a fixed effect for each SNP and a random polygenic effect. We detected eleven genome-wide significant associations on eight different chromosomes. With at least chromosome-wise significance after Bonferroni correction, sixteen SNPs on nine chromosomes showed significant associations with one or more fertility traits. The results confirmed quantitative trait loci on three chromosomes (1, 2 and 20) for fertility traits previously reported for the same breed and one on chromosome four previously detected in Holstein cattle.     

The role of the bovine growth hormone receptor and prolactin receptor genes in milk, fat and protein production in Finnish Ayrshire dairy cattle

We herein report new evidence that the QTL effect on chromosome 20 in Finnish Ayrshire can be explained by variation in two distinct genes, growth hormone receptor (GHR) and prolactin receptor (PRLR). In a previous study in Holstein-Friesian dairy cattle an F279Y polymorphism in the transmembrane domain of GHR was found to be associated with an effect on milk yield and composition. The result of our multimarker regression analysis suggests that in Finnish Ayrshire two QTL segregate on the chromosomal region including GHR and PRLR. By sequencing the coding sequences of GHR and PRLR and the sequence of three GHR promoters from the pooled samples of individuals of known QTL genotype, we identified two substitutions that were associated with milk production traits: the previously reported F-to-Y substitution in the transmembrane domain of GHR and an S-to-N substitution in the signal peptide of PRLR. The results provide strong evidence that the effect of PRLR S18N polymorphism is distinct from the GHR F279Y effect. In particular, the GHR F279Y has the highest influence on protein percentage and fat percentage while PRLR S18N markedly influences protein and fat yield. Furthermore, an interaction between the two loci is suggested.     

Genetic mapping of loci responsible for milk quality parameters in dairy cattle

The review presents a definition of loci controlling quantitative traits (quantitative trait loci, QTLs) and localization of all currently known QTLs responsible for milk production traits in dairy cattle. The QTL number and chromosome localization are verified, with special reference to chromosomes 1, 3, 6, 14, 20, and 23. In a number of cases, close location of QTLs for mastitis and for milk production traits was found. Some aspects of QTL pleiotropy and epistasis are discussed and mapping methods of major QTLs are listed.     

Polymorphism within the 3' flanking region of the bovine growth hormone receptor gene

Growth hormone receptor (GHR) has a major role in the regulation of growth hormone action, and thus, is an obvious candidate gene associated with milk production traits in mammals. The present authors have sequenced 273 bp of the 3' flanking region of the bovine GHR , and found three length variants and one base substitution polymorphism in this region. Allele frequencies of the length variants differ between Finnish native and commercial dairy cattle breeds. The chromosomal localization of GHR was confirmed to bovine chromosome 20 by synteny mapping and linkage analysis.     

Mapping of MYF5, GlR, MYHL, TPIl, IAPP, A2MR and RNR onto sheep chromosome 3q

Five new loci, myogenic factor 5 (MYF5), complement 1 receptor (CIR), myosin-like heavy chain (MYHL), islet amyloid polypeptide (IAPP), and alpha-2-macroglobulin receptor (A2MR), were mapped onto sheep chromosome 3q by Southern hybridization to a panel of chro-mosomally characterized sheep × hamster cell hybrid lines. The location of the triose phosphate isomerase (TPI1) gene and one of the nucleolar organizer regions (RNR) on sheep 3q was confirmed by Southern analysis. This study provides further evidence for the existence of a large conserved chromosomal segment comprising much of sheep chromosome 3q, cattle chromosome 5, and human chromosome 12. The distal evolutionary breakpoint on human chromosome 12, producing the chromosomal segment U23 in cattle marked by aldehyde dehydrogenase (ALDH2), also produces a separate segment in sheep. Neither ALDH2 nor pancreatic lipase (PLA2), which is also distally located on human chromosome 12, were mapped onto sheep chromosome 3q.     

Genetics of heat tolerance for milk yield and quality in Holsteins

Tropical and sub-tropical climates are characterized by high temperature and humidity, during at least part of the year. Consequently, heat stress is common in Holstein cattle and productive and reproductive losses are frequent. Our objectives were as follows: (1) to quantify losses in production and quality of milk due to heat stress; (2) to estimate genetic correlations within and between milk yield (MY) and milk quality traits; and (3) to evaluate the trends of genetic components of tolerance to heat stress in multiple lactations of Brazilian Holstein cows. Thus, nine analyses using two-trait random regression animal models were carried out to estimate variance components and genetic parameters over temperature–humidity index (THI) values for MY and milk quality traits (three lactations: MY×fat percentage (F%), MY×protein percentage (P%) and MY×somatic cell score (SCS)) of Brazilian Holstein cattle. It was demonstrated that the effects of heat stress can be harmful for traits related to milk production and milk quality of Holstein cattle even though most herds were maintained in a modified environment, for example, with fans and sprinklers. For MY, the effect of heat stress was more detrimental in advanced lactations (−0.22 to −0.52 kg/day per increase of 1 THI unit). In general, the mean heritability estimates were higher for lower THI values and longer days in milk for all traits. In contrast, the heritability estimates for SCS increased with increasing THI values in the second and third lactation. For each trait studied, lower genetic correlations (different from unity) were observed between opposite extremes of THI (THI 47 v. THI 80) and in advanced lactations. The genetic correlations between MY and milk quality trait varied across the THI scale and lactations. The genotype×environment interaction due to heat stress was more important for MY and SCS, particularly in advanced lactations, and can affect the genetic relationship between MY and milk quality traits. Selection for higher MY, F% or P% may result in a poor response of the animals to heat stress, as a genetic antagonism was observed between the general production level and specific ability to respond to heat stress for these traits. Genetic trends confirm the adverse responses in the genetic components of heat stress over the years for milk production and quality. Consequently, the selection of Holstein cattle raised in modified environments in both tropical and sub-tropical regions should take into consideration the genetic variation in heat stress.     

The Genome Response to Artificial Selection: A Case Study in Dairy Cattle

Dairy cattle breeds have been subjected over the last fifty years to intense artificial selection towards improvement of milk production traits. In this study, we performed a whole genome scan for differentiation using 42,486 SNPs in the three major French dairy cattle breeds (Holstein, Normande and Montbéliarde) to identify the main physiological pathways and regions which were affected by this selection. After analyzing the population structure, we estimated F_ST within and across the three breeds for each SNP under a pure drift model. We further considered two different strategies to evaluate the effect of selection at the genome level. First, smoothing F_ST values over each chromosome with a local variable bandwidth kernel estimator allowed identifying 13 highly significant regions subjected to strong and/or recent positive selection. Some of them contained genes within which causal variants with strong effect on milk production traits (GHR) or coloration (MC1R) have already been reported. To go further in the interpretation of the observed signatures of selection we subsequently concentrated on the annotation of differentiated genes defined according to the F_ST value of SNPs localized close or within them. To that end we performed a comprehensive network analysis which suggested a central role of somatotropic and gonadotropic axes in the response to selection. Altogether, these observations shed light on the antagonism, at the genome level, between milk production and reproduction traits in highly producing dairy cows.