Quantitative trait locus mapping in dairy cattle by means of selective milk DNA pooling using dinucleotide microsatellite markers: analysis of milk protein percentage

Application of Microsatellites to Trace the Dairy Products Back to the Farm of Origin

The increasing number of food frauds, mainly targeting high quality products, is a rising concern among producers and authorities appointed to food controls. Therefore, the development or implementation of methods to reveal frauds is desired. The genetic traceability of traditional or high-quality dairy products (i.e., products of protected designation of origin, PDO) represents a challenging issue due to the technical problems that arise. The aim of the study was to set up a genetic tool for the origin traceability of dairy products. We investigated the use of Short Tandem Repeats (STRs) to assign milk and cheese to the corresponding producer. Two farms were included in the study, and the blood of the cows, bulk milk, and derived cheese were sampled monthly for one year. Twenty STRs were selected and Polymerase Chain Reactions for each locus were carried out. The results showed that bulk milk and derived cheese express an STR profile composed of a subset of STRs of the lactating animals. A bioinformatics tool was used for the exclusion analysis. The study allowed the identification of a panel of 20 markers useful for the traceability of milk and cheeses, and its effectiveness in the traceability of dairy products obtained from small producers was demonstrated.

Deciphering the Genetic Diversity of Landraces With High-Throughput SNP Genotyping of DNA Bulks: Methodology and Application to the Maize 50k Array

Genebanks harbor original landraces carrying many original favorable alleles for mitigating biotic and abiotic stresses. Their genetic diversity remains, however, poorly characterized due to their large within genetic diversity. We developed a high-throughput, cheap and labor saving DNA bulk approach based on single-nucleotide polymorphism (SNP) Illumina Infinium HD array to genotype landraces. Samples were gathered for each landrace by mixing equal weights from young leaves, from which DNA was extracted. We then estimated allelic frequencies in each DNA bulk based on fluorescent intensity ratio (FIR) between two alleles at each SNP using a two step-approach. We first tested either whether the DNA bulk was monomorphic or polymorphic according to the two FIR distributions of individuals homozygous for allele A or B, respectively. If the DNA bulk was polymorphic, we estimated its allelic frequency by using a predictive equation calibrated on FIR from DNA bulks with known allelic frequencies. Our approach: (i) gives accurate allelic frequency estimations that are highly reproducible across laboratories, (ii) protects against false detection of allele fixation within landraces. We estimated allelic frequencies of 23,412 SNPs in 156 landraces representing American and European maize diversity. Modified Roger's genetic Distance between 156 landraces estimated from 23,412 SNPs and 17 simple sequence repeats using the same DNA bulks were highly correlated, suggesting that the ascertainment bias is low. Our approach is affordable, easy to implement and does not require specific bioinformatics support and laboratory equipment, and therefore should be highly relevant for large-scale characterization of genebanks for a wide range of species.

Systematic Profiling of Short Tandem Repeats in the Cattle Genome

Short tandem repeats (STRs), or microsatellites, are genetic variants with repetitive 2–6 base pair motifs in many mammalian genomes. Using high-throughput sequencing and experimental validations, we systematically profiled STRs in five Holsteins. We identified a total of 60,106 microsatellites and generated the first high-resolution STR map, representing a substantial pool of polymorphism in dairy cattle. We observed significant STRs overlap with functional genes and quantitative trait loci (QTL). We performed evolutionary and population genetic analyses using over 20,000 common dinucleotide STRs. Besides corroborating the well-established positive correlation between allele size and variance in allele size, these analyses also identified dozens of outlier STRs based on two anomalous relationships that counter expected characteristics of neutral evolution. And one STR locus overlaps with a significant region of a summary statistic designed to detect STR-related selection. Additionally, our results showed that only 57.1% of STRs located within SNP-based linkage disequilibrium (LD) blocks whereas the other 42.9% were out of blocks. Therefore, a substantial number of STRs are not tagged by SNPs in the cattle genome, likely due to STR's distinct mutation mechanism and elevated polymorphism. This study provides the foundation for future STR-based studies of cattle genome evolution and selection.

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.

Stories and Challenges of Genome Wide Association Studies in Livestock - A Review

Undoubtedly livestock is one of the major contributors to the economy of any country. The economic value of livestock includes meat, dairy products, fiber, fertilizer etc. Understanding and identifying the associations of quantitative trait loci (QTL) with the economically important traits is believed to substantially benefit the livestock industry. The past two decades have seen a flurry of interest in mapping the QTL associated with traits of economic importance on the genome. With the availability of single nucleotide polymorphism chip of various densities it is possible to identify regions, QTL and genes on the genome that explain the association and its effect on the phenotype under consideration. Remarkable advancement has been seen in genome wide association studies (GWAS) since its inception till the present day. In this review we describe the progress and challenges of GWAS in various livestock species.

Targeted Recombinant Progeny: a design for ultra-high resolution mapping of Quantitative Trait Loci in crosses between inbred or pure lines

BACKGROUND

High-resolution mapping of the loci (QTN) responsible for genetic variation in quantitative traits is essential for positional cloning of candidate genes, and for effective marker assisted selection. The confidence interval (QTL) flanking the point estimate of QTN-location is proportional to the number of individuals in the mapping population carrying chromosomes recombinant in the given interval. Consequently, many designs for high resolution QTN mapping are based on increasing the proportion of recombinants in the mapping population. The "Targeted Recombinant Progeny" (TRP) design is a new design for high resolution mapping of a target QTN in crosses between pure, or inbred lines. It is a three-generation procedure generating a large number of recombinant individuals within a QTL previously shown to contain a QTN. This is achieved by having individuals that carry chromosomes recombinant across the target QTL interval as parents of a large mapping population; most of whom will therefore carry recombinant chromosomes targeted to the given QTL. The TRP design is particularly useful for high resolution mapping of QTN that differentiate inbred or pure lines, and hence are not amenable to high resolution mapping by genome-wide association tests.

RESULTS

In the absence of residual polygenic variation, population sizes required for achieving given mapping resolution by the TRP-F2 design relative to a standard F2 design ranged from 0.289 for a QTN with standardized allele substitution effect = 0.2, mapped to an initial QTL of 0.2 Morgan to 0.041 for equivalent QTN mapped to an initial QTL of 0.02 M. In the presence of residual polygenic variation, the relative effectiveness of the TRP design ranges from 1.068 to 0.151 for the same initial QTL intervals and QTN effect. Thus even in the presence of polygenic variation, the TRP can still provide major savings. Simulation showed that mapping by TRP should be based on 30-50 markers spanning the initial interval; and on at least 50 or more G2 families representing this number of recombination points,.

CONCLUSIONS

The TRP design can be an effective procedure for achieving high and ultra-high mapping resolution of a target QTN previously mapped to a known confidence interval (QTL).

If a Bull Were a Cow, How Much Milk Would He Give?

I became enamored of genetics at an early age. The desire to participate in the rebuilding of my people in our ancient homeland led to dairy science at Rutgers University, to Animal Breeding Plans by J.L. Lush, and to the realization that I could combine genetics and dairy science in animal breeding. It is to my mother-in-law that I owe the felicitous phrasing of the titular scientific question that has occupied my professional life: If a bull were a cow, how much milk would he give? Following my PhD (in 1956), I joined the Volcani Institute in Israel and, later (in 1972), the Applied Genetics group at the Hebrew University. The Applied Genetics group had an active marker lab, and this and a paper by Spickett & Thoday led me to explore genetic markers for quantitative trait loci mapping and marker-assisted selection. A chance encounter with Jacques Beckmann in 1980 opened my eyes to the potential of DNA-level markers for these purposes, and the rest followed.

Signatures of selection identify loci associated with milk yield in sheep

Background

Identification of genomic regions that have been targets of selection for phenotypic traits is one of the most challenging areas of research in animal genetics, particularly in livestock where few annotated genes are available. In this study a genome-wide scan using the Illumina SNP50K Beadchip was performed in the attempt to identify genomic regions associated with milk productivity in sheep. The ovine genomic regions encoding putative candidate genes were compared with the corresponding areas in Bos taurus, as the taurine genome is better annotated.

Results

A total of 100 dairy sheep were genotyped on the Illumina OvineSNP50K Beadchip. The Fisher’s exact test of significance of differences of allele frequency between each pair of the two tails of the distribution of top/worse milk yielders was performed for each marker. The genomic regions where highly divergent milk yielders showed different allele frequencies at consecutive markers was extracted from the OAR v3.1 Ovine (Texel) Genome Assembly, and was compared to the corresponding areas in Bos taurus, allowing the detection of two genes, the Palmdelphin and the Ring finger protein 145. These genes encoded non-synonymous mutations correlated with the marker alleles.

Conclusion

The innovation of this study was to show that the DNA genotyping with the Illumina SNP50K Beadchip allowed to detect genes, and mutations in the genes, which have not yet been annotated in the livestock under investigation.

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.

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.

Mapping QTL affecting resistance to Marek's disease in an F6 advanced intercross population of commercial layer chickens

BACKGROUND

Marek's disease (MD) is a T-cell lymphoma of chickens caused by the Marek's disease virus (MDV), an oncogenic avian herpesvirus. MD is a major cause of economic loss to the poultry industry and the most serious and persistent infectious disease concern. A full-sib intercross population, consisting of five independent families was generated by crossing and repeated intercrossing of two partially inbred commercial White Leghorn layer lines known to differ in genetic resistance to MD. At the F6 generation, a total of 1615 chicks were produced (98 to 248 per family) and phenotyped for MD resistance measured as survival time in days after challenge with a very virulent plus (vv+) strain of MDV.

RESULTS

QTL affecting MD resistance were identified by selective DNA pooling using a panel of 15 SNPs and 217 microsatellite markers. Since MHC blood type (BT) is known to affect MD resistance, a total of 18 independent pool pairs were constructed according to family x BT combination, with some combinations represented twice for technical reasons. Twenty-one QTL regions (QTLR) affecting post-challenge survival time were identified, distributed among 11 chromosomes (GGA1, 2, 3, 4, 5, 8, 9, 15, 18, 26 and Z), with about two-thirds of the MD resistance alleles derived from the more MD resistant parental line. Eight of the QTLR associated with MD resistance, were previously identified in a backcross (BC) mapping study with the same parental lines. Of these, 7 originated from the more resistant line, and one from the less resistant line.

CONCLUSION

There was considerable evidence suggesting that MD resistance alleles tend to be recessive. The width of the QTLR for these QTL appeared to be reduced about two-fold in the F6 as compared to that found in the previous BC study. These results provide a firm basis for high-resolution linkage disequilibrium mapping and positional cloning of the resistance genes.

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.

Chromosomal regions underlying noncoagulation of milk in Finnish Ayrshire cows

About 10% of Finnish Ayrshire cows produce noncoagulating milk, i.e., milk that does not form a curd in a standard 30-min testing time and is thus a poor raw material for cheese dairies. This phenomenon is associated with peak and midlactation, but some cows produce noncoagulating milk persistently. A genomewide scan under a selective DNA pooling method was carried out to locate genomic regions associated with the noncoagulation of milk. On the basis of the hypothesis of the same historical mutation, we pooled the data across sires. Before testing pools for homogeneity, allele intensities were corrected for PCR artifacts, i.e., shadow bands and differential amplification. Results indicating association were verified using daughter design and selective genotyping within families. Data consisted of 18 sire families with 477 genotyped daughters in total, i.e., 12% of each tail of the milk coagulation ability. Data were analyzed using interval mapping under maximum-likelihood and nonparametric methods. BMS1126 on chromosome 2 and BMS1355 on chromosome 18 were associated with noncoagulation of milk across families on an experimentwise 0.1% significance level. By scanning gene databases, we found two potential candidate genes: LOC538897, a nonspecific serine/threonine kinase on chromosome 2, and SIAT4B, a sialyltransferase catalyzing the last step of glycosylation of kappa-casein on chromosome 18. Further studies to determine the role of the candidates in the noncoagulation of milk are clearly needed.

Detection of QTL for growth rate in the blacklip abalone (Haliotis rubra Leach) using selective DNA pooling

The objective of this study was to identify QTL for growth rate in the blacklip abalone Haliotis rubra using selective DNA pooling. Three full-sibling families of H. rubra derived from crosses of wild broodstock were used. DNA was extracted from the largest and smallest 10% of progeny and combined into two pools for each phenotypic tail. The DNA pools were typed with 139 microsatellites, and markers showing significant differences between the peak height ratios of alleles inherited from the parents were individually genotyped and analysed by interval mapping. A strong correlation (r = 0.94, P < 0.001) was found between the t-values from the analysis of pools and the t-values from the analysis of individual genotypes. Based on the interval mapping analysis, QTL were detected on nine linkage groups at a chromosome-wide P < 0.01 and one linkage group at a chromosome-wide P < 0.05. The study demonstrated that selective DNA pooling is efficient and effective as a first-pass screen for the discovery of QTL in an aquaculture species.

Mapping quantitative trait loci affecting susceptibility to Marek's disease virus in a backcross population of layer chickens

Marek's disease (MD), caused by the oncogenic MD avian herpes virus (MDV), is a major source of economic losses to the poultry industry. A reciprocal backcross (BC) population (total 2052 individuals) was generated by crossing two partially inbred commercial Leghorn layer lines known to differ in MDV resistance, measured as survival time after challenge with a (vv+) MDV. QTL affecting resistance were identified by selective DNA pooling using a panel of 198 microsatellite markers covering two-thirds of the chicken genome. Data for each BC were analyzed separately, and as a combined data set. Markers showing significant association with resistance generally appeared in blocks of two or three, separated by blocks of nonsignificant markers. Defined this way, 15 chromosomal regions (QTLR) affecting MDV resistance, distributed among 10 chromosomes (GGA 1, 2, 3, 4, 5, 7, 8, 9, 15, and Z), were identified. The identified QTLR include one gene and three QTL associated with resistance in previous studies of other lines, and three additional QTL associated with resistance in previous studies of the present lines. These QTL could be used in marker-assisted selection (MAS) programs for MDV resistance and as a platform for high-resolution mapping and positional cloning of the resistance genes.

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.

Mapping, fine mapping, and molecular dissection of quantitative trait Loci in domestic animals

Artificial selection has created myriad breeds of domestic animals, each characterized by unique phenotypes pertaining to behavior, morphology, physiology, and disease. Most domestic animal populations share features with isolated founder populations, making them well suited for positional cloning. Genome sequences are now available for most domestic species, and with them a panoply of tools including high-density single-nucleotide polymorphism panels. As a result, domestic animal populations are becoming invaluable resources for studying the molecular architecture of complex traits and of adaptation. Here we review recent progress and issues in the positional identification of genes underlying complex traits in domestic animals. As many phenotypes studied in animals are quantitative, we focus on mapping, fine mapping, and cloning of quantitative trait loci.

Identification of a QTL on BTA20 affecting susceptibility to Mycobacterium avium ssp. paratuberculosis infection in US Holsteins

The objective of this study was to identify QTL affecting susceptibility to Mycobacterium paratuberculosis infection in US Holsteins. Twelve paternal half-sib families were selected for the study based on large numbers of daughters in production and limited relationships among sires. Serum and faecal samples from 4350 daughters of these 12 sires were obtained for disease testing. Case definition for an infected cow was an ELISA sample-to-positive ratio >/=0.25, a positive faecal culture or both. Three families were selected for genotyping based on a high apparent prevalence (6.8-10.4% infected cows), high faecal culture prevalence (46.2-52.9% positive faecal cultures) and large numbers of daughters tested for disease (264-585). DNA pooling was used to genotype cows, with an average of 159 microsatellites within each sire family. Infected cows (the positive pool) were matched with two of their non-infected herdmates in the same lactation (the negative pool) to control for herd and age effects. Eight chromosomal regions putatively linked with susceptibility to M. paratuberculosis infection were identified using a Z-test (P < 0.01). Significant results were more rigorously tested by individually genotyping cows with three to five informative microsatellites within 15 cM of the significant markers identified with the DNA pools. Probability of infection based on both diagnostic tests was estimated for each individual and used as the dependent variable for interval mapping. Based on this analysis, evidence for the presence of a QTL segregating within families on BTA20 was found (chromosome-wide P-value = 0.0319).

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.

The slick hair coat locus maps to chromosome 20 in Senepol-derived cattle

The ability to maintain normal temperatures during heat stress is an important attribute for cattle in the subtropics and tropics. Previous studies have shown that Senepol cattle and their crosses with Holstein, Charolais and Angus animals are as heat tolerant as Brahman cattle. This has been attributed to the slick hair coat of Senepol cattle, which is thought to be controlled by a single dominant gene. In this study, a genome scan using a DNA-pooling strategy indicated that the slick locus is most likely on bovine chromosome 20 (BTA20). Interval mapping confirmed the BTA20 assignment and refined the location of the locus. In total, 14 microsatellite markers were individually genotyped in two pedigrees consisting of slick and normal-haired cattle (n = 36), representing both dairy and beef breeds. The maximum LOD score was 9.4 for a 4.4-cM support interval between markers DIK2416 and BM4107. By using additional microsatellite markers in this region, and genotyping in six more pedigrees (n = 86), the slick locus was further localized to the DIK4835 - DIK2930 interval.

A simplified QTL mapping approach for screening and mapping of novel AFLP markers associated with beef marbling

Genome screening of quantitative trait loci (QTL) for a complex trait is usually costly and highly laborious, as it requires a large number of markers spanning the whole genome. Here we present a simplified approach for screening and mapping of QTL-linked markers for beef marbling using a WagyuxLimousin F(2) reference population. This simplified approach involves integration of the amplified fragment length polymorphism (AFLP) with DNA pooling and selective genotyping and comparative bioinformatics tools. AFLP analysis on two high and two low marbling DNA pools yielded ten visually different markers. Among them, four were confirmed based on individual AFLP validation. Sequencing and in silico characterization assigned two of these AFLP markers to bovine chromosomes 1 (BTA1) and 13 (BTA13), which are orthologous to human chromosomes HSA21q22.2 and HSA10p11.23 with both regions harboring QTL for obesity-related phenotypes. Both AFLP markers showed significantly large additive genetic effects (0.28+/-0.11 on BTA1 and 0.54+/-0.21 on BTA13) on beef-marbling score (BMS) (P<0.05). Overall, this approach is less time consuming, inexpensive and in particular, suitable for screening and mapping QTL-linked markers when targeting one or a few complex traits.

Rapid identification of single nucleotide polymorphisms and estimation of allele frequencies using sequence traces from DNA pools

Millions of single nucleotide polymorphisms (SNP) are now known in the chicken, but the ones that segregate may differ among populations. The current study characterized a simple method to simultaneously identify SNP in populations and to estimate allele frequencies by using relative peak heights of nucleotides in sequencing traces from DNA pools. In the first experiment, known concentration gradients of alleles were produced by measured admixture of DNA from highly inbred individual chickens. In the second experiment, pools of DNA from 100 individuals of each of 3 populations were produced, and the estimated frequencies were compared with those from genotyping of individuals sampled from the same populations. Linear relationships (slope = 0.93 to 1.02) were demonstrated between relative peak heights of nucleotide sequence traces and allelic frequency from both types of experimental samples. The correlation of frequency estimates from pools and individual genotyping was 0.917 to 0.997. Therefore, relative peak heights of nucleotides in sequence traces from DNA isolated from blood pools provide accurate estimates of SNP frequencies of the population from which the pool individuals were sampled and can, therefore, increase efficiency of selecting SNP to genotype in specific populations.

Microsatellite markers associated with resistance to Marek's disease in commercial layer chickens

The objective of the current study was to identify QTL conferring resistance to Marek's disease (MD) in commercial layer chickens. To generate the resource population, 2 partially inbred lines that differed in MD-caused mortality were intermated to produce 5 backcross families. Vaccinated chicks were challenged with very virulent plus (vv+) MD virus strain 648A at 6 d and monitored for MD symptoms. A recent field isolate of the MD virus was used because the lines were resistant to commonly used older laboratory strains. Selective genotyping was employed using 81 microsatellites selected based on prior results with selective DNA pooling. Linear regression and Cox proportional hazard models were used to detect associations between marker genotypes and survival. Significance thresholds were validated by simulation. Seven and 6 markers were significant based on proportion of false positive and false discovery rate thresholds less than 0.2, respectively. Seventeen markers were associated with MD survival considering a comparison-wise error rate of 0.10, which is about twice the number expected by chance, indicating that at least some of the associations represent true effects. Thus, the present study shows that loci affecting MD resistance can be mapped in commercial layer lines. More comprehensive studies are under way to confirm and extend these results.

Family-based association tests for different family structures using pooled DNA

DNA pooling is a cost-effective strategy for genomewise association studies to identify disease genes. In the context of family-based association studies, Risch & Teng (1998) mainly considered families of identical structures to detect associations between genetic markers and disease, and suggested possible approaches to incorporating different family types without a thorough study of their properties. However, families collected in real genetic studies often have different structures and, more importantly, the informativeness of each family structure depends on the disease model which is generally unknown. So there is a need to develop and investigate statistical methods to combine information from diverse family types. In this article, we propose a general strategy to incorporate different family types by assigning each family an "optimal" weight in association tests. In addition, we consider measurement errors in our analysis. When we evaluate our approach under different disease models and measurement errors, we find that our weighting scheme may lead to a substantial reduction in sample size required over the approach suggested by Risch & Teng (1998), and measurement errors may have significant impact on the required sample size when the error rates are not negligible.

Evaluation of DNA pooling for the estimation of microsatellite allele frequencies: a case study using striped bass (Morone saxatilis)

Using striped bass (Morone saxatilis) and six multiplexed microsatellite markers, we evaluated procedures for estimating allele frequencies by pooling DNA from multiple individuals, a method suggested as cost-effective relative to individual genotyping. Using moment-based estimators, we estimated allele frequencies in experimental DNA pools and found that the three primary laboratory steps, DNA quantitation and pooling, PCR amplification, and electrophoresis, accounted for 23, 48, and 29%, respectively, of the technical variance of estimates in pools containing DNA from 2-24 individuals. Exact allele-frequency estimates could be made for pools of sizes 2-8, depending on the locus, by using an integer-valued estimator. Larger pools of size 12 and 24 tended to yield biased estimates; however, replicates of these estimates detected allele frequency differences among pools with different allelic compositions. We also derive an unbiased estimator of Hardy-Weinberg disequilibrium coefficients that uses multiple DNA pools and analyze the cost-efficiency of DNA pooling. DNA pooling yields the most potential cost savings when a large number of loci are employed using a large number of individuals, a situation becoming increasingly common as microsatellite loci are developed in increasing numbers of taxa.

Genome screen for twinning rate QTL in four North American Holstein families

The objective of this study was to identify twinning rate quantitative trait loci (QTL) by typing pooled samples in a preliminary screening followed by interval mapping to test QTL effects. Four elite North American Holstein half-sib sire families with high twinning rate predicted transmitting abilities (PTA) were used in this study. Chromosomes 5, 7, 19 and 23 were not genotyped as these chromosomes were scanned for QTL in these families in a previous study. DNA was extracted from phenotypically extreme sons in each sire family. Two pools were prepared from sons of sires in each phenotypic tail, two each for high and low PTA levels for twinning rates. Each pool contained DNA from 4 to 15% of all sons of the sire depending on family. A total of 268 fluorescently labelled microsatellite markers were tested for heterozygosity in sires. About 135--170 informative markers per family were genotyped using pooled DNA samples. Based on the preliminary evidence for potential twinning rate QTL from pooled typing, interval mapping was performed subsequently on 12 chromosomal regions by family combinations. Evidence of QTL for twinning rate was found in one family on BTA 21 and 29 at a chromosome-wide P<0.05 and on BTA 8, 10 and 14 with a chromosome-wide P<0.01.

Extent and consistency across generations of linkage disequilibrium in commercial layer chicken breeding populations

Recent studies report a surprisingly high degree of marker-to-marker linkage disequilibrium (LD) in ruminant livestock populations. This has important implications for QTL mapping and marker-assisted selection. This study evaluated LD between microsatellite markers in a number of breeding populations of layer chickens using the standardized chi-square (chi(2')) measure. The results show appreciable LD among markers separated by up to 5 cM, decreasing rapidly with increased separation between markers. The LD within 5 cM was strongly conserved across generations and differed among chromosomal regions. Using marker-to-marker LD as an indication for marker-QTL LD, a genome scan of markers spaced 2 cM apart at moderate power would have good chances of uncovering most QTL segregating in these populations. However, of markers showing significant trait associations, only 57% are expected to be within 5 cM of the responsible QTL, and the remainder will be up to 20 cM away. Thus, high-resolution LD mapping of QTL will require dense marker genotyping across the region of interest to allow for interval mapping of the QTL.

Accurate determination of microsatellite allele frequencies in pooled DNA samples

Pooling of DNA samples instead of individual genotyping can speed up genetic association studies. However, for microsatellite markers, the electrophoretic pattern of DNA pools can be complex, and procedures for deriving allele frequencies are often confounded by PCR-induced stutter artefacts. We have developed a mathematical procedure to remove stutter noise and accurately determine allele frequencies in pools. A stutter correction model can be reliably derived from one standard 'training set' of the same 10 individual DNA samples for each marker, which can also include heterozygous patterns with partially overlapping peaks. Compared with earlier methods, this reduces the number of genotypes needed in the training set considerably, and allows standardization of analyses for different markers. Moreover, the use of a procedure that fits all data simultaneously makes the method less sensitive to aberrant data. The model was tested with 34 markers, 18 of which were newly defined from human sequence data. Allele frequencies derived from stutter-corrected DNA pool patterns were compared with the summed individual genotyping results of all the individuals in the pools (n = 109 and n = 64). We show that the model is robust and accurately extracts allele frequencies from pooled DNA samples for 32 of the 34 microsatellite markers tested. Finally, we performed a case-control study in celiac disease and found that weakly associated disease alleles, identified by individual genotyping, were only detectable in pools after stutter correction. This efficient method for correcting stutter artefacts in microsatellite markers enables large-scale genetic association studies using DNA pools to be performed.

Novel algorithm for automated genotyping of microsatellites

Microsatellites or short tandem repeats (STRs) are abundant in the human genome with easily assayed polymorphisms, providing powerful genetic tools for mapping both Mendelian and complex traits. Microsatellite genotyping requires detection of the products of polymerase chain reaction (PCR) amplification by electrophoresis, and analysis of the peak data for discrimination of the true allele. A high-throughput genotyping approach requires computer-based automation at both the detection and analysis phases. In order to achieve this, complicated peak patterns from individual alleles must be interpreted in order to assign alleles. Previous methods consider limited types of noise peaks and cannot provide enough accuracy. By pattern recognition of various types of noise peaks, such as stutter peaks and additional peaks, we have achieved an overall average accuracy of 94% for allele calling in our actual data. Our algorithm is crucial for a high-throughput genotyping system for microsatellite markers by reducing manual editing and human errors.

Ovine alpha-amylase genes: isolation, linkage mapping and association analysis with milk traits

On the basis of comparisons between cattle and sheep genome mapping information the ovine alpha-amylase gene was examined as a possible genetic marker for milk traits in sheep. The objective of the present study was to isolate, map and determine whether this gene is a candidate gene for milk traits. DNA fragments (832 and 2360 bp) corresponding to two different AMY genes were isolated, and one SNP in intron 3 and one GTG deletion in exon 3 of the 2360 bp DNA fragment were found. The 2360 bp ovine AMY DNA fragment was located on chromosome 1 by linkage mapping using the International Mapping Flock. No association was found between estimated breeding values for milk yield, protein and fat contents and AMY genotypes in a daughter design comprising 13 Manchega families with an average of 29 daughters (12-62) per sire.

Identification of an ovulation rate QTL in cattle on BTA14 using selective DNA pooling and interval mapping

Increased twinning incidence in beef cattle has the potential to improve production efficiency. However, phenotypic selection for twinning rate is difficult because of the trait's low heritability and the long time interval necessary to collect phenotypic records. Therefore, this trait and the correlated trait of ovulation rate are ideal candidates for marker-assisted selection. The objective of this study was to complete a genome-wide search for ovulation rate quantitative trait loci (QTL) in two related sire families. The families (paternal halfsib sires 839802 and 839803) were from a population of cattle selected for ovulation rate at the USDA Meat Animal Research Center, Clay Center, Nebraska. Putative ovulation rate QTL have previously been identified in the 839802 family on chromosomes 7 and 19; however, marker coverage in the original scan was not complete. This study fills the gaps in marker coverage of the earlier study by adding approximately 60 informative microsatellites to each sire family. Each family was genotyped using selective DNA pooling. Sons and daughters were included in either the high or low pool based on their estimated breeding value deviations from the mid-parent average (EBVMD) for ovulation rate. Approximately 40% (839802) and 26% (839803) of available progeny comprised the high and low pools combined. Pooled typing revealed possible associations (nominal P < 0.05) between ovulation rate and marker genotype for 11 and 15 microsatellites in the 839802 and 839803 families, respectively. Subsequent interval mapping strengthened support for the presence of an ovulation rate QTL on BTA14 (chromosome-wise P < 0.02).

Verification of selective DNA pooling methodology through identification and estimation of the DGAT1 effect

The discovery of markers linked to genes that are responsible for traits of interest to the dairy industry might prove useful because they could aid in selection and breeding decisions. We have developed a selective DNA pooling methodology to allow us to efficiently screen the bovine genome in order to find genes responsible for production traits. Using markers on chromosome 14 as a test case, we identified a gene (DGAT1) previously known to affect three traits (fat yield, protein yield and total milk yield). Furthermore, we predicted similar effects to those previously shown for DGAT1 in a New Zealand Holstein-Friesian herd. Additionally, we showed a low error rate (1.6%) for the pooling procedure. Hence we are confident that we can apply this procedure to an entire genome scan in the search for quantitative trait loci (QTL).

High-resolution mapping of quantitative trait loci by selective recombinant genotyping

Selective recombinant genotyping (SRG) is a three-stage procedure for high-resolution mapping of a QTL that has previously been mapped to a known confidence interval (target C.I.). In stage 1, a large mapping population is accessed and phenotyped, and a proportion, P, of the high and low tails is selected. In stage 2, the selected individuals are genotyped for a pair of markers flanking the target C.I., and a group of R individuals carrying recombinant chromosomes in the target interval are identified. In stage 3, the recombinant individuals are genotyped for a set of M markers spanning the target C.I. Extensive simulations showed that: (1) Standard error of QTL location (SEQTL) decreased when QTL effect (d) or population size (N) increased, but was constant for given "power factor" (PF = d(2)N); (2) increasing the proportion selected in the tails beyond 0.25 had only a negligible effect on SEQTL; and (3) marker spacing in the target interval had a remarkably powerful effect on SEQTL, yielding a reduction of up to 10-fold in going from highest (24 cM) to lowest (0.29 cM) spacing at given population size and QTL effect. At the densest marker spacing, SEQTL of 1.0-0.06 cM were obtained at PF = 500-16,000. Two new genotyping procedures, the half-section algorithm and the golden section/half-section algorithm, allow the equivalent of complete haplotyping of the target C.I. in the recombinant individuals to be achieved with many fewer data points than would be required by complete individual genotyping.

Polymorphism of the bovine CSN1S1 promoter: linkage mapping, intragenic haplotypes, and effects on milk production traits

The bovine CSN1S1 5' flanking region (CSN1S1-5') was screened for polymorphisms in different cattle breeds. Single-strand conformation polymorphisms (SSCP) and sequence analyses revealed four alleles (1-4), two of them being new allelic forms (3 and 4). Sequences were deposited in GenBank with accession numbers AF549499-502. In alleles 1 and 4, potential transcription factor binding sites are altered by the mutations. Using SSCP analysis, all four alleles were identified in German Holsteins. Six intragenic haplo-types comprising CSN1S1-5' (alleles 1, 2, 3, 4) and exon 17 (CSN1S1*B and C) genotypes were found. Linkage mapping using half-sib families from the German QTL project positioned CSN1S1 between the markers FBN14 and CSN3, with 5.6 cM distance between CSN1S1 and CSN3. Variance analysis, using family and CSN1S1 promoter genotypes as fixed effects, of breeding values and deregressed proofs for milk production traits (milk, fat, and protein yield and also fat and protein percentage) revealed significant effects on protein percentage when all families and genotypes were considered. Contrast calculations assigned a highly significant effect to genotype 24, which was associated with highest LS-means for protein percentage breeding values. As CSN1S1 is one of the main caseins in milk, this could be an effect of mutations in regulatory elements in the promoter region. An effect on milk yield breeding values was indicated for genotype 12, but is probably caused by a linked locus.

Quantitative technologies for allele frequency estimation of SNPs in DNA pools

We have compared several genotyping methods to assess their applicability to single nucleotide polymorphism (SNP) allele frequency estimation in DNA pools. The accuracy of these methods (restriction fragment length polymorphism, real-time pyrophosphate DNA sequencing, single base extension with fluorescently labeled ddNTPs, homogeneous 5'-nuclease assay, and MALDI-TOF mass spectrometry) was tested by calculating the standard deviation among heterozygous individuals (which are natural DNA pools with 50% representation of each allele) and by estimating allele frequency in artificial pools. We show that although the methods differ in their accuracy, they can all serve for quantification of allele frequency in DNA pools with reasonable accuracy. We found that the influence of the error variance attributed to pool construction on quantification accuracy is insignificant and is SNP dependent.

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.

A genome scan for quantitative trait loci affecting milk production in Norwegian dairy cattle

An autosomal genome scan for quantitative trait loci (QTL) affecting milk production traits was carried out on the Norwegian Dairy Cattle population. Six half-sibling families with a total of 285 sons organized according to a granddaughter design were analyzed by a multiple marker regression method. Suggestive QTL for one or several of the five milk traits (milk yield, protein percentage, protein yield, fat percentage and fat yield) were detected on chromosomes 3, 5, 6, 11, 13, 18 and 20. Among these results, the findings on chromosomes 3, 6, and 20 are highly supported by literature. The most convincing result was found close to marker FBN9 on chromosome 6, where a QTL was detected with alleles that cause a marked reduction in both protein and fat percentages and an increase in milk yield. The results for fat and protein percentage were highly significant even after accounting for multiple testing across the genome. Using bootstrapping, a 95% confidence interval for the position of the QTL for the percentage traits on chromosome 6 was estimated to 16 cM.