Multiple quantitative trait locus analysis of bovine chromosome 6 in the Israeli Holstein population by a daughter design

Nine Israeli Holstein sire families with 2978 daughters were analyzed for quantitative trait loci effects on chromosome 6 for five milk production traits by a daughter design. All animals were genotyped for 2 markers. The three families with significant effects were genotyped for up to 10 additional markers spanning positions 0-122 cM of BTA6. Two sires were segregating for a locus affecting protein and fat percentage near position 55 cM with an estimated substitution effect of 0.18% protein, which is equivalent to one phenotypic standard deviation. This locus was localized to a confidence interval of 4 cM. One of these sires was also heterozygous for a locus affecting milk, fat, and protein production near the centromere. The hypothesis of two segregating loci was verified by multiple regression analysis. A third sire was heterozygous for a locus affecting milk and protein percentage near the telomeric end of the chromosome. Possible candidates for the major quantitative gene near position 55 cM were determined by comparative mapping. IBSP and SSP1 were used as anchors for the orthologous region on human chromosome 4. Twelve genes were detected within a 2-Mbp sequence. None of these genes have been previously associated with lactogenesis.     

Structural variation around prolactin gene linked to quantitative traits in an elite Holstein sire family

Summary Digestion of genomic DNA with the restriction endonuclease Avail disclosed a probable insertion deletion of approximately 200 base pairs (bp) near the prolactin gene. Two alleles were apparent as three distinct hybridization patterns. These alleles were statistically associated with quantitative trait loci among sons of one elite Holstein sire family. The favorable genotype was correlated with the presence of a 1.15-kb hybridization band inherited from the sire when genomic DNA was probed with a full-length cDNA for prolactin. Pedigree estimates of genetic merit among genotypes were similar, differing by only 19.3 kg for milk in ancestor merit. Comparisons of genetic estimates for quantitative yield traits in offspring of this heterozygous sire showed significant (P<0.05) differences between homozygous genotypes for predicted difference milk (PDM), predicted difference dollars (PD$), cheese yield dollars, and protein dollars. The estimated differences between homozygous genotypes for USDA Transmitting Abilities of PDM, PD$, Cheese Yield $ and Protein $ were 282.93 kg, $74.35, $48.58 and $53.67, respectively. However, the estimated breeding values from progeny ranged over 900 kg in transmitting ability for milk. Frequency of the favorable marker allele was estimated to be 0.231 in the elite cow population used as dams of sons. These results demonstrate the potential of molecular biological techniques to discriminate between individuals within a family and to predict breeding values for selection schemes.     

Linkage between amylase-1 locus and a major gene for milk fat content in cattle

Summary
Linkage between the amylase-1 ( Am-1 ) locus and a quantitative trait locus influencing fat content in milk was studied in offspring from heterozygous sires of the Swedish Red and White dairy breed. The effect on bull breeding values for fat content was estimated as interactions between sire and paternal Am-1 allele using a model eliminating the direct effects of sire and Am-1 allele. There were strong indications of linkage, confirming results of previous studies. The interaction was caused by strong associations in 7 out of 14 sire families. A test for within-family variance heterogeneity performed on the whole population of breeding bulls also supported the presence of a major gene for fat content in milk. The results indicate that there is genetic linkage between the Am-1 locus and a locus with large effect on milk fat content.     

Prediction of informativeness for microsatellite markers among progeny of sires used for detection of economic trait loci in dairy cattle

Individual loci affecting economic traits can be located using genetic linkage. Application of either daughter or granddaughter designs requires determination of allele origin in the progeny. If only the sires and their progeny are genotyped, the paternal allele origin of progeny having the same genotype as the sire cannot be determined. The expected frequency of informative sons can be predicted for each sire and genetic marker from the allele frequencies in the population. The accuracy of a predictor of the frequency of informative progeny was tested on 103 grandsire x microsatellite combinations. Number of sons per grandsire varied from 24 to 129. Allele frequencies in the population were estimated by genotyping seven sires. The regression of the frequency of informative sons on the predicted frequency was 1.04 with a zero intercept model. Thus, considering the large number of genetic markers available for analysis, predicted informative frequency is a useful criterion for selection of genetic markers.     

The number of loci that affect milk production traits in dairy cattle

We have used the results of an experiment mapping quantitative trait loci (QTL) affecting milk yield and composition to estimate the total number of QTL affecting these traits. We did this by estimating the number of segregating QTL within a half-sib daughter design using logic similar to that used to estimate the "false discovery rate" (FDR). In a half-sib daughter design with six sire families we estimate that the average sire was heterozygous for approximately 5 QTL per trait. Also, in most cases only one sire was heterozygous for any one QTL; therefore at least 30 QTL were likely to be segregating for these milk production traits in this Holstein population.     

A QTL affecting milk yield and composition maps to bovine chromosome 20: a confirmation

As part of a whole genome scan undertaken to detect quantitative trait loci (QTL) affecting milk yield and composition, we have genotyped a granddaughter design comprising 1152 sons for six microsatellite markers spanning bovine chromosome 20. An analysis performed across families provided strong evidence (experiment-wise P -values < 0·01) for the presence of a QTL affecting primarily protein percentage towards the telomeric end of the chromosome. A founder sire, shown in a previous study to segregate for a similar QTL in the corresponding chromosome region, was characterized by 29 and 57 sons and maternal grandsons, respectively, in the present design. Sorting corresponding sons and grandsons by paternal or grandpaternal allele provided significant evidence for the segregation of a QTL on chromosome 20. Altogether these results confirm the location of a QTL affecting milk production on bovine chromosome 20.     

Fine mapping quantitative trait loci affecting milk production traits on bovine chromosome 6 in a Chinese Holstein population

To fine map the previously detected quantitative trait loci （QTLs） affecting milk production traits on bovine chromosome 6 （BTA6）, 15 microsatellite markers situated within an interval of 14.3 cM spanning from BMS690 to BM4528 were selected and 918 daughters of 8 sires were genotyped. Two mapping approaches, haplotype sharing based LD mapping and single marker regression mapping, were used to analyze the data. Both approaches revealed a quantitative trait locus （QTL） with significant effects on milk yield, fat yield and protein yield located in the segment flanked by markers BMS483 and MNB209, which spans a genetic distance of 0.6 cM and a physical distance of 1.5 Mb. In addition, the single marker regression mapping also revealed a QTL affecting fat percentage and protein percentage at marker DIK2291. Our fine mapping work will facilitate the cloning of candidate genes underlying the QTLs for milk production traits.     

Associations between blood groups, blood protein polymorphisms and breeding values for production traits in Swedish Red and White Dairy bulls

The relationships of nine blood group systems and two blood protein polymorphisms with breeding values for several production traits were examined in dairy cattle of the Swedish Red and White (SRB) breed. The material consisted of 2212 bulls; the bulls were performance tested for growth rate and their breeding values for milk yield, fat and protein content in milk were estimated from progeny tests. The direct effect of marker alleles or marker phenotypes was analysed in a multiple regression model. Several significant associations were found; many supported earlier findings. However, the contribution of the markers to the total variation of the breeding values was very small. Linkage between marker loci and production loci was studied in offspring from heterozygous sires by estimating the interactions between sire and marker alle using a model eliminating the direct effects of sire and marker alle. There were strong indications of linkage between some marker loci (e.g. B, J and Am-1 loci) and loci with large effects on production traits.     

Estimation of quantitative trait locus allele frequency via a modified granddaughter design

A method is described on the basis of a modification of the granddaughter design to obtain estimates of quantitative trait loci (QTL) allele frequencies in dairy cattle populations and to determine QTL genotypes for both homozygous and heterozygous grandsires. The method is based on determining the QTL allele passed from grandsires to their maternal granddaughters using haplotypes consisting of several closely linked genetic markers. This method was applied to simulated data of 10 grandsire families, each with 500 granddaughters, and a QTL with a substitution effect of 0.4 phenotypic standard deviations and to actual data for a previously analyzed QTL in the center of chromosome 6, with substitution effect of 1 phenotypic standard deviation on protein percentage. In the simulated data the standard error for the estimated QTL substitution effect with four closely linked multiallelic markers was only 7% greater than the expected standard error with completely correct identification of QTL allele origin. The method estimated the population QTL allelic frequency as 0.64 +/- 0.07, compared to the simulated value of 0.7. In the actual data, the frequency of the allele that increases protein percentage was estimated as 0.63 +/- 0.06. In both data sets the hypothesis of equal allelic frequencies was rejected at P < 0.05.     

A quantitative trait locus for live weight maps to bovine Chromosome 23

A multiple-marker mapping approach was used to search for quantitative trait loci (QTLs) affecting production, health, and fertility traits in Finnish Ayrshire dairy cattle. As part of a whole-genome scan, altogether 469 bulls were genotyped for six microsatellite loci in 12 families on Chromosome (Chr) 23. Both multiple-marker interval mapping with regression and maximum-likelihood methods were applied with a granddaughter design. Eighteen traits, belonging to 11 trait groups, were included in the analysis. One QTL exceeded experiment level and one QTL genome level significance thresholds. Across-families analysis provided strong evidence (P_experiment= 0.0314) for a QTL affecting live weight. The QTL for live weight maps between markers BM1258 and BoLA DRBP1. A QTL significant at genome level (P_genome= 0.0087) was mapped for veterinary treatment, and the putative QTL probably affects susceptibility to milk fever or ketosis. In addition, three traits exceeded the chromosome 5% significance threshold: protein percentage of milk, calf mortality (sire), and milking speed. In within-family analyses, protein percentage was associated with markers in one family (LOD score = 4.5). Received: 14 December 1998 / Accepted: 28 March 1998     

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.     

The Effects of Selection on Linkage Analysis for Quantitative Traits

The effects of within-sample selection on the outcome of analyses detecting linkage between genetic markers and quantitative traits were studied. It was found that selection by truncation for the trait of interest significantly reduces the differences between marker genotype means thus reducing the power to detect linked quantitative trait loci (QTL). The size of this reduction is a function of proportion selected, the magnitude of the QTL effect, recombination rate between the marker locus and the QTL, and the allele frequency of the QTL. Proportion selected was the most influential of these factors on bias, e.g., for an allele substitution effect of one standard deviation unit, selecting the top 80%, 50% or 20% of the population required 2, 6 or 24 times the number of progeny, respectively, to offset the loss of power caused by this selection. The effect on power was approximately linear with respect to the size of gene effect, almost invariant to recombination rate, and a complex function of QTL allele frequency. It was concluded that experimental samples from animal populations which have been subjected to even minor amounts of selection will be inefficient in yielding information on linkage between markers and loci influencing the quantitative trait under selection.     

Restriction fragment length polymorphism in amplification products of the bovine butyrophilin gene: assignment of bovine butyronhilin to bovine chromosome 23

A polymorphism was identified in the bovine butyrophilin (BTN) gene by digesting poly-merase chain reaction products with the restriction enzyme Hae III. This polymorphism was segregating in a Holstein-Friesian sire selected as part of an ongoing study directed towards the identification of quantitative trait loci affecting milk composition. Screening of a half-sib family established for the heterozygous sire allowed the localization of BTN to bovine chromosome 23 (BTA23).     

Microsatellite markers associated with quantitative trait loci controlling antibody response to Escherichia coli and Salmonella enteritidis in young broilers

A unique resource population was produced to facilitate detection of microsatellite markers associated with quantitative trait loci controlling antibody (Ab) response in broiler chickens. Three F1 males were produced by mating two lines divergently selected on Ab response to Escherichia coli vaccination. Each F1 male was mated with females from four genetic backgrounds: F1, high-Ab line (HH), low-Ab line and commercial line, producing three resource families, each with four progeny types. About 1700 chicks were immunized with E. coli and Salmonella enteritidis vaccines. Selective genotyping was conducted on the individuals with highest or lowest average Ab to E. coli and S. enteritidis within each progeny type in each sire family. Twelve markers were significantly associated with Ab to E. coli and six of them were also associated with Ab to S. enteritidis, mostly exhibiting a similar low effect (approximately 0.35 phenotypic SD) in all progeny types. Four markers exhibited a highly significant and much larger effect (approximately 1.7 SD), but only in progeny of females from the HH, suggesting that a backcross to the high parental line should be preferred over the commonly used F2 population. Results from two markers suggested a quantitative trait locus on chromosome 2 around 400 cM. The marker MCW0083, significant in two sire families, is closely linked to the bone morphogenetic protein 2 (BMP2) gene, known to be associated with the control of T-cell transformation in humans.     

Genetic linkage maps of rose constructed with new microsatellite markers and locating QTL controlling flowering traits

New microsatellites markers [simple sequence repeat (SSR)] have been isolated from rose and integrated into an existing amplified fragment-length polymorphism genetic map. This new map was used to identify quantitative trait locus (QTL) controlling date of flowering and number of petals. From a rose bud expressed sequence tag (EST) database of 2,556 unigenes and a rose genomic library, 44 EST-SSRs and 20 genomic-SSR markers were developed, respectively. These new rose SSRs were used to expand genetic maps of the rose interspecific F₁ progeny. In addition, SSRs from other Rosaceae genera were also tested in the mapping progeny. Genetic maps for the two parents of the progeny were constructed using pseudo-testcross mapping strategy. The maps consist of seven linkage groups of 105 markers covering 432 cM for the maternal map and 136 markers covering 438 cM for the paternal map. Homologous relationships among linkage groups between the maternal and paternal maps were established using SSR markers. Loci controlling flowering traits were localised on genetic maps as a major gene and QTL for the number of petals and a QTL for the blooming date. New SSR markers developed in this study will provide tools for the establishment of a consensus linkage map for roses that combine traits and markers in various rose genetic maps.     

A whole genome scan for quantitative trait loci affecting milk protein percentage in Israeli-Holstein cattle, by means of selective milk DNA pooling in a daughter design, using an adjusted false discovery rate criterion

Selective DNA pooling was employed in a daughter design to screen all bovine autosomes for quantitative trait loci (QTL) affecting estimated breeding value for milk protein percentage (EBVP%). Milk pools prepared from high and low daughters of each of seven sires were genotyped for 138 dinucleotide microsatellites. Shadow-corrected estimates of sire allele frequencies were compared between high and low pools. An adjusted false discovery rate (FDR) method was employed to calculate experimentwise significance levels and empirical power. Significant associations with milk protein percentage were found for 61 of the markers (adjusted FDR = 0.10; estimated power, 0.68). The significant markers appear to be linked to 19--28 QTL. Mean allele substitution effects of the putative QTL averaged 0.016 (0.009--0.028) in units of the within-sire family standard deviation of EBVP% and summed to 0.460 EBVP%. Overall QTL heterozygosity was 0.40. The identified QTL appear to account for all of the variation in EBVP% in the population. Through use of selective DNA pooling, 4400 pool data points provided the statistical power of 600,000 individual data points.     

Quantitative trait loci for fertility traits in Finnish Ayrshire cattle

A whole genome scan was carried out to detect quantitative trait loci (QTL) for fertility traits in Finnish Ayrshire cattle. The mapping population consisted of 12 bulls and 493 sons. Estimated breeding values for days open, fertility treatments, maternal calf mortality and paternal non-return rate were used as phenotypic data. In a granddaughter design, 171 markers were typed on all 29 bovine autosomes. Associations between markers and traits were analysed by multiple marker regression. Multi-trait analyses were carried out with a variance component based approach for the chromosomes and trait combinations, which were observed significant in the regression method. Twenty-two chromosome-wise significant QTL were detected. Several of the detected QTL areas were overlapping with milk production QTL previously identified in the same population. Multi-trait QTL analyses were carried out to test if these effects were due to a pleiotropic QTL affecting fertility and milk yield traits or to linked QTL causing the effects. This distinction could only be made with confidence on BTA1 where a QTL affecting milk yield is linked to a pleiotropic QTL affecting days open and fertility treatments.     

Linkage analysis of quantitative traits in an interspecific cross of tomato (lycopersicon esculentum x lycopersicon pimpinellifolium) by means of genetic markers

Linkage relationships between loci affecting quantitative traits (QTL) and marker loci were examined in an interspecific cross between Lycopersicon esculentum and Lycopersicon pimpinellifolium. Parental lines differed for six morphological markers and for four electrophoretic markers. Almost 1700 F-2 plants were scored with respect to the genetic markers and also with respect to 18 quantitative traits. Major genes affecting the quantitative traits were not found, but out of 180 possible marker x trait combinations, 85 showed significant quantitative effects associated with the genetic markers. The average marker-associated main effect was on the order of 6% of the mean value of the trait. Most of the main effects were apparently due to linkage of QTL to the marker loci rather than to pleiotropy. Fourteen of the traits showed at least one highly significant effect of opposite sign to the overall difference between the parental lines, demonstrating the ability of this design to uncover cryptic genetic variation. Significant variance and skewness effects on the quantitative traits were found to be associated with the genetic markers, suggesting the possible presence of loci affecting the variance and shape of quantitative trait distribution in a population. Most marker-associated quantitative effects showed some degree of dominance, generally in the direction of the L. pimpinellifolium parent. When the significant marker-associated effects were examined in pairs, 12% showed significant interaction effects. The results of this study illustrate the potential usefulness of this type of analysis for the detailed genetic investigation of quantitative trait variation in suitably marked populations.     

Application of the a posteriori granddaughter design to the Holstein genome

An a posteriori granddaughter design was applied to estimate quantitative trait loci genotypes of sires with many sons in the US Holstein population. The results of this analysis can be used to determine concordance between specific polymorphisms and segregating quantitative trait loci. Determination of the actual polymorphisms responsible for observed genetic variation should increase the accuracy of genomic evaluations and rates of genetic gain. A total of 52 grandsire families, each with ⩾100 genotyped sons with genetic evaluations based on progeny tests, were analyzed for 33 traits (milk, fat and protein yields; fat and protein percentages; somatic cell score (SCS); productive life; daughter pregnancy rate; heifer and cow conception rates; service-sire and daughter calving ease; service-sire and daughter stillbirth rates; 18 conformation traits; and net merit). Of 617 haplotype segments spanning the entire bovine genome and each including ~5×10⁶ bp, 5 cM and 50 genes, 608 autosomal segments were analyzed. A total of 19 335 unique haplotypes were found among the 52 grandsires. There were a total of 133 chromosomal segment-by-trait combinations, for which the nominal probability of significance for the haplotype effect was <10⁻⁸, which corresponds to genome-wide significance of <10⁻⁴. The number of chromosomal regions that met this criterion by trait ranged from one for rear legs (rear view) to seven for net merit. For each of the putative quantitative trait loci, at least one grandsire family had a within-family contrast with a t-value of >3. Confidence intervals (CIs) were estimated by the nonparametric bootstrap for the largest effect for each of nine traits. The bootstrap distribution generated by 100 samples was bimodal only for net merit, which had the widest 90% CI (eight haplotype segments). This may be due to the fact that net merit is a composite trait. For all other chromosomes, the CI spanned less than a third of the chromosome. The narrowest CI (a single haplotype segment) was found for SCS. It is likely that analysis by more advanced methods could further reduce CIs at least by half. These results can be used as a first step to determine the actual polymorphisms responsible for observed quantitative variation in dairy cattle.     

DNA fingerprint bands applied to linkage analysis with quantitative trait loci in chickens

Summary
An efficient approach to detect association between quantitative traits and bands of DNA fingerprint patterns uses intra-family tail analysis, which compares fingerprints of DNA mixes from individuals at the two tails of a phenotypic distribution. In analysis of 67 paternal half-sibs of a meat-type chicken family, of 57 sire bands generated by two probes, one sire-specific band (S6–6) was associated with abdominal fat deposition. The band effect was estimated by a linear model analysis to be 0–88 standard deviations, or about 30% of the family mean. The association between band S6–6 and abdominal fat was further examined by testing progeny of paternal half-sibs of the chickens which were used in the tail analysis, establishing genetic linkage between the DNA marker and a genetic locus affecting abdominal fat deposition.