E. Lipkin

Quantitative trait loci affecting milk yield and protein percentage in a three-country Brown Swiss population.

Quantitative trait loci (QTL) mapping projects have been implemented mainly in the Holstein dairy cattle breed for several traits. The aim of this study is to map QTL for milk yield (MY) and milk protein percent (PP) in the Brown Swiss cattle populations of Austria, Germany, and Italy, considered in this study as a single population. A selective DNA pooling approach using milk samples was applied to map QTL in 10 paternal half-sib daughter families with offspring spanning from 1,000 to 3,600 individuals per family. Three families were sampled in Germany, 3 in Italy, 1 in Austria and 3 jointly in Austria and Italy. The pools comprised the 200 highest and 200 lowest performing daughters, ranked by dam-corrected estimated breeding value for each sire-trait combination. For each tail, 2 independent pools, each of 100 randomly chosen daughters, were constructed. Sire marker allele frequencies were obtained by densitometry and shadow correction analyses of 172 genome-wide allocated autosomal markers. Particular emphasis was placed on Bos taurus chromosomes 3, 6, 14, and 20. Marker association for MY and PP with a 10% false discovery rate resulted in nominal P-values of 0.071 and 0.073 for MY and PP, respectively. Sire marker association tested at a 20% false discovery rate (within significant markers) yielded nominal P-values of 0.031 and 0.036 for MY and PP, respectively. There were a total of 36 significant markers for MY, 33 for PP, and 24 for both traits; 75 markers were not significant for any of the traits. Of the 43 QTL regions found in the present study, 10 affected PP only, 8 affected MY only, and 25 affected MY and PP. Remarkably, all 8 QTL regions that affected only MY in the Brown Swiss, also affected MY in research reported in 3 Web-based QTL maps used for comparison with the findings of this study (http://www.vetsci.usyd.edu.au/reprogen/QTL_Map/; http://www.animalgenome.org/QTLdb/cattle.html; http://bovineqtl.tamu.edu/). Similarly, all 10 QTL regions in the Brown Swiss that affected PP only, affected only PP in the databases. Thus, many QTL appear to be common to Brown Swiss and other breeds in the databases (mainly Holstein), and an appreciable fraction of QTL appears to affect MY or PP primarily or exclusively, with little or no effect on the other trait. Although QTL information available today in the Brown Swiss population can be utilized only in a within family marker-assisted selection approach, knowledge of QTL segregating in the whole population should boost gene identification and ultimately the implementation and efficiency of an individual genomic program.

Association of Bovine Fibroblast Growth Factor 2 (FGF2) Gene with Milk Fat and Productive Life: An Example of the Ability of the Candidate Pathway Strategy to Identify Quantitative Trait Genes

Fibroblast growth factor 2 (FGF2) is expressed in the bovine mammary gland and may play a role in the development and reorganization of the mammary gland. It is also expressed by the uterine endometrium throughout the estrous cycle and early pregnancy. The FGF2 was chosen for this study because it regulates the expression of interferon-tau, a key member of the signal transduction pathway involved in milk production. In previous studies, we reported the association of several genes in this pathway with milk production and health traits in dairy cattle. The objective of this study was to examine the association of FGF2 polymorphisms with milk composition, somatic cell score, and productive life in 3 Holstein cattle populations from the United States and Israel. The pooled DNA sequencing approach was used to identify single nucleotide polymorphisms (SNP) in FGF2. Sequencing of a total of 6.4 kb including 3 exons of the gene revealed only one SNP (A/G) in intron 1 at position 11646. This SNP was investigated for association with production traits in 2,773 individuals from 3 Holstein populations: the granddaughter-design Cooperative Dairy DNA Repository and the daughter-design University of Wisconsin populations from the United States and a daughter-design population from Israel. For both the Israeli and the UW populations, FGF2 variants were associated with fat yield and percentage, somatic cell score, and productive life with significant dominance and complete dominance effects. For the Cooperative Dairy DNA Repository population, no significant associations were observed for the examined traits. Given that FGF2 was chosen for this study because of its role in the interferon-tau signal transduction pathway and was found to be associated with production traits, results suggest that the candidate pathway could be an attractive strategy to search for candidate quantitative trait genes.

A genome scan for quantitative trait loci affecting milk somatic cell score in Israeli and Italian Holstein cows by means of selective DNA pooling with single- and multiple-marker mapping.

Mastitis is an important and common dairy cattle disease affecting milk yield, quality, and consumer safety as well as cheese yields and quality. Animal welfare and residues of the antibiotics used to treat mastitis cause public concern. Considerable genetic variation may allow selection for increased resistance to mastitis. Because of high genetic correlation to milk somatic cell score (SCS), SCS can serve as a surrogate trait for mastitis resistance. The present study intended to identify quantitative trait loci (QTL) affecting SCS in Israeli and Italian Holstein dairy cattle (IsH and ItH, respectively), using selective DNA pooling with single and multiple marker mapping. Milk samples of 4,788 daughters of 6 IsH and 7 ItH sires were used to construct sire-family high- and low-tail pools, which were genotyped at 123 (IsH) and 133 (ItH) microsatellite markers. Shadow correction was used to obtain pool allele frequency estimates. Frequency difference between the tails and empirical standard error of D, SE(D), were used to obtain P-values. All markers significant by single marker mapping were also significant by multiple marker mapping, but not vice versa. Combining both populations, 22 QTL on 21 chromosomes were identified; all corresponded to previous reports in the literature. Confidence intervals were set by chi-squared drop method. Heterozygosity of QTL was estimated at 44.2%. Allele substitution effects ranged from 1,782 to 4,930 cells/mL in estimated breeding value somatic cell count units. Most (80%) of the observed variation in estimated breeding value somatic cell score could be explained by the QTL identified under the stringent criteria. The results found here can be used as a basis for further genome-wide association studies for the same trait.

Fractioned DNA Pooling: A New Cost-Effective Strategy for Fine Mapping of Quantitative Trait Loci

Selective DNA pooling (SDP) is a cost-effective means for an initial scan for linkage between marker and quantitative trait loci (QTL) in suitable populations. The method is based on scoring marker allele frequencies in DNA pools from the tails of the population trait distribution. Various analytical approaches have been proposed for QTL detection using data on multiple families with SDP analysis. This article presents a new experimental procedure, fractioned-pool design (FPD), aimed to increase the reliability of SDP mapping results, by “fractioning” the tails of the population distribution into independent subpools. FPD is a conceptual and structural modification of SDP that allows for the first time the use of permutation tests for QTL detection rather than relying on presumed asymptotic distributions of the test statistics. For situations of family and cross mapping design we propose a spectrum of new tools for QTL mapping in FPD that were previously possible only with individual genotyping. These include: joint analysis of multiple families and multiple markers across a chromosome, even when the marker loci are only partly shared among families; detection of families segregating (heterozygous) for the QTL; estimation of confidence intervals for the QTL position; and analysis of multiple-linked QTL. These new advantages are of special importance for pooling analysis with SNP chips. Combining SNP microarray analysis with DNA pooling can dramatically reduce the cost of screening large numbers of SNPs on large samples, making chip technology readily applicable for genomewide association mapping in humans and farm animals. This extension, however, will require additional, nontrivial, development of FPD analytical tools.

Variation in the ovocalyxin-32 gene in commercial egg-laying chickens and its relationship with egg production and egg quality traits.

Avian eggshell quality is an important trait for commercial egg production, as the eggshell is the primary packaging material and antimicrobial barrier for the internal food resource. Strong eggshells are essential to ensure that eggs can reach their final destination without damage. Ovocalyxin-32 (OCX32) is a matrix protein found within the outer layers of the eggshell and in the cuticle. Numerous reports in the literature have identified association between variants in the gene encoding this protein, OCX32, and various eggshell quality traits. Thus, OCX32 is a candidate gene for selection for eggshell traits in commercial poultry populations. Sequencing of exons 2-6 of the OCX32 gene in eight elite brown and white eggshell commercial egg-laying lines revealed 28 SNPs and one SNP/indel. Eighteen of these SNPs were predicted to alter the amino acid sequence of the protein. Clusters of SNPs in complete linkage disequilibrium were found in both exons 2 and 6. A total of 19 different versions or protein-sequence haplotypes of the OCX32 protein were inferred, revealing considerable variation within commercial lines. Genotypes for 13 of the SNPs were determined for 330-1819 individuals per line. Trait association studies revealed a significant effect of OCX32 on shell color in white egg lines and line-specific significant effects on albumen height, early egg weight, puncture score, and yolk weight. Three of the lines showed a significant change in OCX32 frequency over time, indicating selection pressure for certain variants of this gene during the breeding program.

Effect of Quantitative Trait Loci for Milk Protein Percentage on Milk Protein Yield and Milk Yield in Israeli Holstein Dairy Cattle

Although numerous quantitative trait loci (QTL) mapping studies involving milk protein percent (PP), milk yield (MY), and protein yield (PY) have been carried out, there has not been any systematic evaluation of the effects of individual QTL on these 3 interrelated traits. Consequently, the aim of the present study was to investigate the effects on MY and PY of QTL for PP previously mapped in various laboratories. The study, based on selective DNA pooling of milk samples, included 10 Israeli Holstein artificial insemination bulls, each the sire of 1,800 or more milk-recorded daughters. For each sire-trait combination across the 10 sires, milk samples of the highest and lowest daughters with respect to estimated breeding values for PP, PY, and MY were collected for pooling. A total of 134 dinucleotide microsatellites distributed over 25 bovine autosomes were used. An empirical standard error for marker-QTL linkage testing was calculated based on the variation among split samples within the same tail. Threshold comparison-wise error rate P-values were set to control proportion of false positives at P = 0.10 level for declaring significant effects at the marker-trait level. Estimates of the number of true null hypotheses for each trait were obtained from the histogram of marker comparison-wise error rate P-values. Based on these estimates, effective power of the experiment at the marker-trait level was estimated as 0.75, 0.41, and 0.73 for PP, PY, and MY. The proportion of heterozygosity at the QTL was estimated as 0.46, 0.39, and 0.40, respectively. After correcting for incomplete power and proportion of false positives, it was estimated that 38.7 and 37.5% of the markers affecting PP and MY, respectively, also affected PY. Of the markers affecting PY, 68.9 and 76.5%, respectively, also affected PP and MY. Apparently, none of the significant markers affected PY exclusively, and only 6.5 and 16.0%, respectively, affected PP or MY exclusively. Thus, almost all significant markers, and by inference almost all QTL, had effects on at least 2 of the 3 traits.

Extensive Long-Range and Nonsyntenic Linkage Disequilibrium in Livestock Populations: Deconstruction of a Conundrum

Great interest was aroused by reports, based on microsatellite markers, of high levels of statistically significant long-range and nonsyntenic linkage disequilibrium (LD) in livestock. Simulation studies showed that this could result from population family structure. In contrast, recent SNP-based studies of livestock populations report much lower levels of LD. In this study we show, on the basis of microsatellite data from four cattle populations, that high levels of long-range LD are indeed obtained when using the multi-allelic D′ measure of LD. Long-range and nonsyntenic LD are exceedingly low, however, when evaluated by the standardized chi-square measure of LD, which stands in relation to the predictive ability of LD. Furthermore, specially constructed study populations provided no evidence for appreciable LD resulting from family structure at the grandparent level. We propose that the high statistical significance and family structure effects observed in the earlier studies are due to the use of large sample sizes, which accord high statistical significance to even slight deviations from asymptotic expectations under the null hypothesis. Nevertheless, even after taking sample size into account, our results indicate that microsatellites testify to the presence of usable LD at considerably wider separation distances than SNPs, suggesting that use of SNP haplotypes may considerably increase the usefulness of a given fixed SNP array.

Genome-wide association study for somatic cell score in Valdostana Red Pied cattle breed using pooled DNA

BackgroundMastitis is a major disease of dairy cattle occurring in response to environmental exposure to infective agents with a great economic impact on dairy industry. Somatic cell count (SCC) and its log transformation in somatic cell score (SCS) are traits that have been used as indirect measures of resistance to mastitis for decades in selective breeding. A selective DNA pooling (SDP) approach was applied to identify Quantitative Trait Loci (QTL) for SCS in Valdostana Red Pied cattle using the Illumina Bovine HD BeadChip.ResultsA total of 171 SNPs reached the genome-wide significance for association with SCS. Fifty-two SNPs were annotated within genes, some of those involved in the immune response to mastitis. On BTAs 1, 2, 3, 4, 9, 13, 15, 17, 21 and 22 the largest number of markers in association to the trait was found. These regions identified novel genomic regions related to mastitis (1-Mb SNP windows) and confirmed those already mapped. The largest number of significant SNPs exceeding the threshold for genome-wide significant signal was found on BTA 15, located at 50.43-51.63 Mb.ConclusionsThe genomic regions identified in this study contribute to a better understanding of the genetic control of the mastitis immune response in cattle and may allow the inclusion of more detailed QTL information in selection programs.

Fine mapping and association analysis of a quantitative trait locus for milk production traits on Bos taurus autosome 4

To fine map a quantitative trait locus (QTL) affecting milk production traits previously associated with microsatellite RM188, we implemented an interval mapping analysis by using microsatellite markers in a large Israeli Holstein half-sib sire family, and linkage disequilibrium (LD) mapping in a large set of US Holstein bulls. Interval mapping located the target QTL to the near vicinity of RM188. For the LD mapping, we identified 42 single nucleotide polymorphisms (SNP) in 15 genes in a 12-Mb region on bovine chromosome 4. A total of 24 tag SNP were genotyped in 882 bulls belonging to the University of California Davis archival collection of Holstein bull DNA samples with predicted transmitted ability phenotypes. Marker-to-marker LD analysis revealed 2 LD blocks, with intrablock r(2) values of 0.10 and 0.46, respectively; outside the blocks, r(2) values ranged from 0.002 to 0.23. A standard additive/dominance model using the generalized linear model procedure of SAS and the regression module of HelixTree software were used to test marker-trait associations. Single nucleotide polymorphism 9 on ARL4A, SNP10 on XR_027435.1, SNP12 on ETV1, SNP21 on SNX13, and SNP24 were significantly associated with milk production traits. We propose the interval encompassing ARL4A and SNX13 genes as a candidate region in bovine chromosome 4 for a concordant QTL related to milk protein traits in dairy cattle. Functional studies are needed to confirm this result.

A whole genome scan for QTL affecting milk protein percentage in Italian Holstein cattle, applying selective milk DNA pooling and multiple marker mapping in a daughter design

We report on a complete genome scan for quantitative trait loci (QTL) affecting milk protein percentage (PP) in the Italian Holstein-Friesian cattle population, applying a selective DNA pooling strategy in a daughter design. Ten Holstein-Friesian sires were chosen, and for each sire, about 200 daughters, each from the high and low tails of estimated breeding value for PP, were used to construct milk DNA pools. Sires and pools were genotyped for 181 dinucleotide microsatellites covering all cattle autosomes. Sire marker allele frequencies in the pools were obtained by shadow correction of peak height in the electropherograms. After quality control, pool data from eight sires were used for all subsequent analyses. The QTL heterozygosity estimate was lower than that of similar studies in other cattle populations. Multiple marker mapping identified 19 QTL located on 14 chromosomes (BTA1, 2, 3, 4, 5, 6, 8, 9, 12, 14, 17, 20, 23 and 27). The sires were also genotyped for seven polymorphic sites in six candidate genes (ABCG2, SPP1, casein kappa, DGAT1, GHR and PRLR) located within QTL regions of BTA6, 14 and 20 found in this study. The results confirmed or excluded the involvement of some of the analysed markers as the causative polymorphic sites of the identified QTL. The QTL identified, combined with genotype data of these candidate genes, will help to identify other quantitative trait genes and clarify the complex QTL patterns observed for a few chromosomes. Overall, the results are consistent with the Italian Holstein population having been under long-term selection for high PP.