Detection of quantitative trait loci affecting twinning rate in Israeli Holsteins by the daughter design

Twinning rate was analyzed in the Israeli Holstein dairy cattle population by the multiple-trait animal model, and a daughter design genome scan for quantitative trait loci was performed. Each parity was considered a separate trait. Heritabilities of twinning rate were very low, but increased by parity from 0.01 in first parity to 0.03 in fifth parity. All genetic correlations among parities were >0.77, but all environmental correlations were <0.07. Genetic correlations between twinning rate and female fertility (measured as the inverse of the number of inseminations to conception) in the first 3 parities were negative for all 9 parity-by-trait combinations. All environmental correlations were very small, but generally negative. The overall genetic trend since 1985 was positive at 0.02% twinning/yr, whereas the phenotypic trends were positive for parities 3 and 4 and negative for the other parities, but all trends were quite small. A total of 5,221 cows, daughters of 11 sires, were genotyped for 73 markers spanning all 29 autosomes. There were 9 markers with significant effects on twinning rate at P < 0.05, for a false discovery rate of 0.4; thus, about 5 of these probably represent true effects. Significant effects were found on chromosomes 1, 6, 7, 8, 14, 15, and 23. Of these, 3 effects were significant at P < 0.01, for a false discovery rate of 0.24. All 11 families were analyzed by interval mapping of chromosome 7. Only 2 families showed nominally significant effects, but chromosome-wise significance at P < 0.05 was not obtained for either family. Suggestive evidence of quantitative trait loci near the beginning of the chromosome and near position 50 cM were found in both families. Sire 3070 also had a significant effect for female fertility near the beginning of the chromosome. There was also evidence for a third quantitative trait loci at the end of the chromosome for sire 2357.