A primary assembly of a bovine haplotype block map based on a 15,036-single-nucleotide polymorphism panel genotyped in holstein-friesian cattle

High density linkage disequilibrium maps of chromosome 14 in Holstein and Angus cattle

BACKGROUND

Linkage disequilibrium (LD) maps can provide a wealth of information on specific marker-phenotype relationships, especially in areas of the genome where positional candidate genes with similar functions are located. A recently published high resolution radiation hybrid map of bovine chromosome 14 (BTA14) together with the bovine physical map have enabled the creation of more accurate LD maps for BTA14 in both dairy and beef cattle.

RESULTS

Over 500 Single Nucleotide Polymorphism (SNP) markers from both Angus and Holstein animals had their phased haplotypes estimated using GENOPROB and their pairwise r2 values compared. For both breeds, results showed that average LD extends at moderate levels up to 100 kilo base pairs (kbp) and falls to background levels after 500 kbp. Haplotype block structure analysis using HAPLOVIEW under the four gamete rule identified 122 haplotype blocks for both Angus and Holstein. In addition, SNP tagging analysis identified 410 SNPs and 420 SNPs in Holstein and Angus, respectively, for future whole genome association studies on BTA14. Correlation analysis for marker pairs common to these two breeds confirmed that there are no substantial correlations between r-values at distances over 10 kbp. Comparison of extended haplotype homozygosity (EHH), which calculates the LD decay away from a core haplotype, shows that in Holstein there is long range LD decay away from the DGAT1 region consistent with the selection for milk fat % in this population. Comparison of EHH values for Angus in the same region shows very little long range LD.

CONCLUSION

Overall, the results presented here can be applied in future single or haplotype association analysis for both populations, aiding in confirming or excluding potential polymorphisms as causative mutations, especially around Quantitative Trait Loci regions. In addition, knowledge of specific LD information among markers will aid the research community in selecting appropriate markers for whole genome association studies.

Extent of linkage disequilibrium in Holstein cattle in North America

The success of fine-scale mapping and genomic selection depends mainly on the strength of linkage disequilibrium (LD) between markers and causal mutations. With Lewontin's measure of LD (known as D'), high levels of LD that extend over several million base pairs have been reported in livestock. However, this measure of LD can be strongly biased upward by small samples and by low allele frequencies. The aim of this study was to characterize the level and extent of LD in Holstein cattle in North America (Canada and the United States for purposes of this study) by using the squared correlation of the alleles at 2 loci (r(2)). The Affymetrix MegAllele GeneChip Bovine Mapping 10K single nucleotide polymorphism (SNP) array was used to genotype 821 bulls, from which 497 were used in the analysis of the extent of LD. A total of 5,564 SNP were used after filtering out SNP with more than 5% of Mendelian inconsistencies, with more than 20% missing genotypes, or with a minor allele frequency of less than 10%. Analysis of syntenic pairs revealed that useful LD (measured as r(2) > 0.3) occurred at distances shorter than 100 kb. Linkage disequilibrium decayed very rapidly, within a few hundred kilobase pairs. In addition, no substantial LD between unlinked loci was found. Using a sliding window analysis, we observed an irregular pattern of LD across the genome. These findings suggest that to capture useful LD, which is required for whole-genome fine mapping and genomic selection, a denser SNP map would be needed.

Fine mapping of quantitative trait loci affecting female fertility in dairy cattle on BTA03 using a dense single-nucleotide polymorphism map

Fertility quantitative trait loci (QTL) are of high interest in dairy cattle since insemination failure has dramatically increased in some breeds such as Holstein. High-throughput SNP analysis and SNP microarrays give the opportunity to genotype many animals for hundreds SNPs per chromosome. In this study, due to these techniques a dense SNP marker map was used to fine map a QTL underlying nonreturn rate measured 90 days after artificial insemination previously detected with a low-density microsatellite marker map. A granddaughter design with 17 Holstein half-sib families (926 offspring) was genotyped for a set of 437 SNPs mapping to BTA3. Linkage analysis was performed by both regression and variance components analysis. An additional analysis combining both linkage analysis and linkage-disequilibrium information was applied. This method first estimated identity-by-descent probabilities among base haplotypes. These probabilities were then used to group the base haplotypes in different clusters. A QTL explaining 14% of the genetic variance was found with high significance (P < 0.001) at position 19 cM with the linkage analysis and four sires were estimated to be heterozygous (P < 0.05). Addition of linkage-disequilibrium information refined the QTL position to a set of narrow peaks. The use of the haplotypes of heterozygous sires offered the possibility to give confidence in some peaks while others could be discarded. Two peaks with high likelihood-ratio test values in the region of which heterozygous sires shared a common haplotype appeared particularly interesting. Despite the fact that the analysis did not fine map the QTL in a unique narrow region, the method proved to be able to handle efficiently and automatically a large amount of information and to refine the QTL position to a small set of narrow intervals. In addition, the QTL identified was confirmed to have a large effect (explaining 13.8% of the genetic variance) on dairy cow fertility as estimated by nonreturn rate at 90 days.

Linkage disequilibrium structures in cattle and their application to breed identification testing

We examined the extent of linkage disequilibrium (LD) block lengths in four breed populations: Japanese Black, Angus, Hereford and Holstein. Three chromosomal regions in which QTL were previously mapped in Japanese Black populations were scanned with 84 microsatellite markers. The estimated LD lengths in these four purebred populations varied from 535 to 683 kb, which is much shorter than the values reported previously. Our findings suggest that QTL can be mapped in sub-centimorgan regions in these populations using an LD-mapping method. We also developed breed identification methods to distinguish Japanese Black from Angus, Hereford, Holstein and F(1) animals (Japanese Black x Holstein) respectively using the haplotypic frequencies of a pair of markers in the breed populations. After assessing the distributions of posterior probabilities to be Japanese Black, we obtained several pairs of markers that completely distinguished Japanese Black from the other breeds. We also obtained several combinations of six markers that completely distinguished Japanese Black animals from F(1) animals.

Extent of genome-wide linkage disequilibrium in Australian Holstein-Friesian cattle based on a high-density SNP panel

BACKGROUND

The extent of linkage disequilibrium (LD) within a population determines the number of markers that will be required for successful association mapping and marker-assisted selection. Most studies on LD in cattle reported to date are based on microsatellite markers or small numbers of single nucleotide polymorphisms (SNPs) covering one or only a few chromosomes. This is the first comprehensive study on the extent of LD in cattle by analyzing data on 1,546 Holstein-Friesian bulls genotyped for 15,036 SNP markers covering all regions of all autosomes. Furthermore, most studies in cattle have used relatively small sample sizes and, consequently, may have had biased estimates of measures commonly used to describe LD. We examine minimum sample sizes required to estimate LD without bias and loss in accuracy. Finally, relatively little information is available on comparative LD structures including other mammalian species such as human and mouse, and we compare LD structure in cattle with public-domain data from both human and mouse.

RESULTS

We computed three LD estimates, D', Dvol and r2, for 1,566,890 syntenic SNP pairs and a sample of 365,400 non-syntenic pairs. Mean D' is 0.189 among syntenic SNPs, and 0.105 among non-syntenic SNPs; mean r2 is 0.024 among syntenic SNPs and 0.0032 among non-syntenic SNPs. All three measures of LD for syntenic pairs decline with distance; the decline is much steeper for r2 than for D' and Dvol. The value of D' and Dvol are quite similar. Significant LD in cattle extends to 40 kb (when estimated as r2) and 8.2 Mb (when estimated as D'). The mean values for LD at large physical distances are close to those for non-syntenic SNPs. Minor allelic frequency threshold affects the distribution and extent of LD. For unbiased and accurate estimates of LD across marker intervals spanning < 1 kb to > 50 Mb, minimum sample sizes of 400 (for D') and 75 (for r2) are required. The bias due to small samples sizes increases with inter-marker interval. LD in cattle is much less extensive than in a mouse population created from crossing inbred lines, and more extensive than in humans.

CONCLUSION

For association mapping in Holstein-Friesian cattle, for a given design, at least one SNP is required for each 40 kb, giving a total requirement of at least 75,000 SNPs for a low power whole-genome scan (median r2 > 0.19) and up to 300,000 markers at 10 kb intervals for a high power genome scan (median r2 > 0.62). For estimation of LD by D' and Dvol with sufficient precision, a sample size of at least 400 is required, whereas for r2 a minimum sample of 75 is adequate.

SNiPORK - a microarray of SNPs in candidate genes potentially associated with pork yield and quality - development and validation in commercial breeds

SNiPORK is an oligonucleotide microarray based on the arrayed primer extension (APEX) technique, allowing genotyping of single nucleotide polymorphisms (SNPs) in genes of interest for pork yield and quality traits. APEX consists of a sequencing reaction primed by an oligonucleotide anchored with its 5' end to a glass slide and terminating one nucleotide before the polymorphic site. Extension with one fluorescently labeled dideoxynucleotide complementary to the template reveals the polymorphism. Ninety SNPs were selected from those associated directly or potentially with pork traits. Of the 90 SNPs, 5 did not produce a positive signal. For 85 SNPs, 100% repeatiblity was proved by double genotyping of 13 randomly chosen boars. In addition, the accuracy of genotyping was verified in 2 sib-families by a Mendelian inheritance of 49-50 homozygous genotypes from sire to sons. Three genotype discrepancies were found (97% accuracy rate). All inaccurities were confirmed by an alternative method (sequencing and PCR-RFLP assays). Moreover, the exclusion power of the chip was evalueted by an SNP inheritance analysis of unrelated boars within each sib-family. In the validation step, 88 boars (13 Pietrain, 31 Landrace, 16 Large White, 8 Duroc, 7 Hampshire x Pietrain crosses, and 13 other hybrid lines) were screened to validate SNPs. Among the 85 selected SNPs, 12 were found to be monoallelic, the rest showing at least two genotypes for the entire population under study. The primary application of the SNiPORK chip is the simultaneous genotyping of dozens of SNPs to study gene interaction and consequently better understand the genetic background of pork yield and quality. The chip may prospectively be used for evolutionary studies, evaluation of genetic distances between wild and domestic pig breeds, traceability tests, as well as the starting point for developing a platform for identification and paternity analysis.

Comparison of genomic and traditional BLUP-estimated breeding value accuracy and selection response under alternative trait and genomic parameters

Accuracy of prediction of estimated breeding values based on genome-wide markers (GEBV) and selection based on GEBV as compared with traditional Best Linear Unbiased Prediction (BLUP) was examined for a number of alternatives, including low heritability, number of generations of training, marker density, initial distributions, and effective population size (Ne). Results show that the more the generations of data in which both genotypes and phenotypes were collected, termed training generations (TG), the better the accuracy and persistency of accuracy based on GEBV. GEBV excelled for traits of low heritability regardless of initial equilibrium conditions, as opposed to traditional marker-assisted selection, which is not useful for traits of low heritability. Effective population size is critical for populations starting in Hardy-Weinberg equilibrium but not for populations started from mutation-drift equilibrium. In comparison with traditional BLUP, GEBV can exceed the accuracy of BLUP provided enough TG are included. Unfortunately selection rapidly reduces the accuracy of GEBV. In all cases examined, classic BLUP selection exceeds what was possible for GEBV selection. Even still, GEBV could have an advantage over traditional BLUP in cases such as sex-limited traits, traits that are expensive to measure, or can only be measured on relatives. A combined approach, utilizing a mixed model with a second random effect to account for quantitative trait loci in linkage equilibrium (the polygenic effect) was suggested as a way to capitalize on both methodologies.

Genomewide analysis of secretory activation in mouse models

Mouse models have been widely used to elucidate the biology of mammary gland development and secretory activation. Recent advances in the availability of genomic resources for mice will generate a renewed effort to define the genetic basis of lactation phenotypes and help identify candidate gene pathways. Specific aspects of these advances are relevant to the dairy industry and may provide a rationale for improving milk production in the dairy cow. Differences are evident in mammary gland morphology and various characteristics of milk production of inbred mouse strains, but few studies have undertaken any systematic phenotypic analysis of the different inbred strains of mice for lactation performance. Whole genome association analysis using recent strain-specific genotype data and detailed phenotype measurements from available inbred strains, along with transcript profiling of divergent inbred strains for lactation performance, provides a valuable approach to identify putative candidate genes and associated pathways underlying dairy QTL intervals. Here we discuss the utility of integrating mouse phenomic and genomic resources for understanding secretory activation in the mammary gland.

Whole genome linkage disequilibrium maps in cattle

Background

Bovine whole genome linkage disequilibrium maps were constructed for eight breeds of cattle. These data provide fundamental information concerning bovine genome organization which will allow the design of studies to associate genetic variation with economically important traits and also provides background information concerning the extent of long range linkage disequilibrium in cattle.

Results

Linkage disequilibrium was assessed using r² among all pairs of syntenic markers within eight breeds of cattle from the Bos taurus and Bos indicus subspecies. Bos taurus breeds included Angus, Charolais, Dutch Black and White Dairy, Holstein, Japanese Black and Limousin while Bos indicus breeds included Brahman and Nelore. Approximately 2670 markers spanning the entire bovine autosomal genome were used to estimate pairwise r² values. We found that the extent of linkage disequilibrium is no more than 0.5 Mb in these eight breeds of cattle.

Conclusion

Linkage disequilibrium in cattle has previously been reported to extend several tens of centimorgans. Our results, based on a much larger sample of marker loci and across eight breeds of cattle indicate that in cattle linkage disequilibrium persists over much more limited distances. Our findings suggest that 30,000–50,000 loci will be needed to conduct whole genome association studies in cattle.

Genetic and haplotypic structure in 14 European and African cattle breeds

To evaluate and compare the extent of LD in cattle, 1536 SNPs, mostly localized on BTA03, were detected in silico from available sequence data using two different methods and genotyped on samples from 14 distinct breeds originating from Europe and Africa. Only 696 SNPs could be validated, confirming the importance of trace-quality information for the in silico detection. Most of the validated SNPs were informative in several breeds and were used for a detailed description of their genetic structure and relationships. Results obtained were in agreement with previous studies performed on microsatellite markers and using larger samples. In addition, the majority of the validated SNPs could be mapped precisely, reaching an average density of one marker every 311 kb. This allowed us to analyze the extent of LD in the different breeds. Decrease of LD with physical distance across breeds revealed footprints of ancestral LD at short distances (<10 kb). As suggested by the haplotype block structure, these ancestral blocks are organized, within a breed, into larger blocks of a few hundred kilobases. In practice, such a structure similar to that already reported in dogs makes it possible to develop a chip of <300,000 SNPs, which should be efficient for mapping purposes in most cattle breeds.