Search

Twinning is a complex trait with negative impacts on health and reproduction, which cause economic loss in dairy production. Several twinning rate quantitative trait loci (QTL) have been detected in previous studies, but confidence intervals for QTL location are broad and many QTL are unreplicated. To identify genomic regions or genes associated with twinning rate, QTL analysis based on linkage combined with linkage disequilibrium (LLD) and individual marker associations was conducted across the genome using high-throughput single nucleotide polymorphism (SNP) genotypes. A total of 9919 SNP markers were genotyped with 200 sires and sons in 19 half-sib North American Holstein dairy cattle families. After SNPs were genotyped, informative markers were selected for genome-wide association tests and QTL searches. Evidence for twinning rate QTL was found throughout the genome. Thirteen markers significantly associated with twinning rate were detected on chromosomes 2, 5 and 14 (P < 2.3 x 10(-5)). Twenty-six regions on fourteen chromosomes were identified by LLD analysis at P < 0.0007. Seven previously reported ovulation or twinning rate QTL were supported by results of single marker association or LLD analyses. Single marker association analysis and LLD mapping were complementary tools for the identification of putative QTL in this genome scan.

Linkage disequilibrium was estimated using 7119 single nucleotide polymorphism markers across the genome and 200 animals from the North American Holstein cattle population. The analysis of maternally inherited haplotypes revealed strong linkage disequilibrium (r(2) > 0.8) in genomic regions of approximately 50 kb or less. While linkage disequilibrium decays as a function of genomic distance, genomic regions within genes showed greater linkage disequilibrium and greater variation in linkage disequilibrium compared with intergenic regions. Identification of haplotype blocks could characterize the most common haplotypes. Although maximum haplotype block size was over 1 Mb, mean block size was 26-113 kb by various definitions, which was larger than that observed in humans ( approximately 10 kb). Effective population size of the dairy cattle population was estimated from linkage disequilibrium between single nucleotide polymorphism marker pairs in various haplotype ranges. Rapid reduction of effective population size of dairy cattle was inferred from linkage disequilibrium in recent generations. This result implies a loss of genetic diversity because of the high rate of inbreeding and high selection intensity in dairy cattle. The pattern observed in this study indicated linkage disequilibrium in the current dairy cattle population could be exploited to refine mapping resolution. Changes in effective population size during past generations imply a necessity of plans to maintain polymorphism in the Holstein population.

Twinning in cattle is a complex trait that is associated with economic loss and health issues such as abortion, dystocia, and reduced calf survival. Twinning-rate QTL have been detected previously on BTA5 in the North American Holstein and Norwegian dairy cattle populations and in a USDA herd selected for high twinning rate. In previous work with the North American Holstein population, the strongest evidence for a QTL was obtained from analysis of an extended, multiple-generation family. Using additional animals, an increased density of SNP marker association tests, and a combined linkage and linkage disequilibrium mapping method, we refined the position of this QTL in the North American Holstein population. Two sets of twinning-rate predicted transmitting abilities estimated during 2 different time periods in the North American dairy cattle population were used to provide validation of results. A total of 106 SNP and 3 microsatellites were used to scan the genomic region between 5 and 80 Mb on BTA5. Combined linkage-linkage disequilibrium analysis identified significant evidence for QTL within the 25- to 35-Mb and 64- to 70-Mb regions of BTA5. The IGF-1 gene (IGF1) was examined as a positional candidate gene and an SNP in intron 2 of IGF1 was significantly associated with twinning rate by using both data sets (P = 0.003 and P = 1.05 x 10(-6)). Replication of this association in other cattle populations will be required to examine the extent of linkage disequilibrium with the underlying quantitative trait nucleotide across breeds.