Genetic analysis of rumination time based on an analysis of 77,697 Israeli dairy cows.

Reduction of methane emission may become necessary for sustainable milk production. Several studies indicate a relationship between rumination time and the level of methane emission. The objectives of the current study were to estimate environmental factors affecting daily rumination time in high-yielding dairy cattle, genetic parameters for rumination time across parities, and environmental and genetic correlations between rumination time and economic traits, and to predict the consequence of inclusion of this trait in the Israeli breeding index. The data included more than 30 million daily records from 77,697 Israeli Holstein cows for rumination time and milk production. A lactation measure of daily rumination time per cow was computed as the mean of the residuals from a linear model analysis with rumination time as the dependent variable. The independent variables were parity and the square root, linear, quadradic and inverse of DIM by parity. Because of the shape of the lactation curve for rumination time, separate linear model analyses were performed for records up to 40 DIM and records with >40 DIM. The phenotypic correlation between first- and second-parity lactations for rumination time was almost 0.8, and close to 0.7 for milk. The heritability of lactation rumination time was close to 0.44 for parities 1 to 3. Heritability for milk production decreased from 0.5 in first parity to 0.3 in third parity. For both traits, genetic correlations among parities were all >0.9. Thus, for routine genetic analysis of rumination time, records in the different parities can be considered the same trait. The genetic correlation between rumination time and milk on first parity was 0.25 and increased slightly with increase in parity. Genetic correlations between rumination time, based on the first 40 DIM, were economically unfavorable with retained placenta but economically favorable with metritis, ketosis, and displaced abomasum. Genetic correlations between rumination time and the 9 traits included in the Israeli breeding index (milk, fat, and protein production; SCS; female fertility; herd-life; milk production persistency; calving ease; and calf mortality) were all economically favorable, except for the correlation of 0.17 with SCS. With the current index, daily rumination time with a current mean of 536 min and SD of 90 min is expected to increase by 11 min/d after 10 yr of selection. Inclusion of this trait with a positive index weight equivalent to 10% of the index should increase rumination time by 19 min. All changes in expected gains due to inclusion of rumination time in the index were economically positive, except for fat and SCS. Inclusion of rumination time in the index should result in 1 kg less gain in fat, a miniscule gain of 0.03 for SCS; and gains of 1.5 kg protein, 0.3% female fertility, and 5 d herd-life. Even though the case for a genetic correlation between rumination time and methane emission is still weak, inclusion of this trait in the commercial index may be justified, considering that equipment is now commercially available for routine recording at reasonable cost.

Genetic and genomic analysis of age at first insemination in Israeli dairy cattle.

We performed a genetic analysis of age at first insemination, including estimation of the heritability and genetic correlations with other economic traits, and the consequences of including this trait in the Israeli selection index. The genetic factors affecting age at first insemination were determined via GWAS. Five data sets were analyzed. Data sets 1, 2, and 3 were used to compute variance components among age at first insemination, first calving age, days from first insemination to calving, and the 9 traits included in the Israel breeding index. Heritabilities for age at first insemination, calving age, and days from first insemination to calving in Israeli Holsteins as computed by REML individual animal model analyses of 273,239 Israeli Holstein cows were 0.072, 0.042, and 0.014. The estimated genetic correlation between the first 2 traits was 0.88. In addition to the fact that heritability of age at first insemination is 1.7 times the heritability for calving, the former trait has the advantage that the number of records is greater, and the records are generated earlier. Absolute values of the genetic and residual correlations between age at first insemination and the 9 traits included in the Israeli index were all less than 0.2. Data set 4 included first insemination dates of 1,181,600 calves born from 1985 through 2018. Genetic evaluations were computed by a single trait animal model. Annual phenotypic and genetic trends for age at first calving for calves born since 1985 were "positive," that is, economically negative, at 0.320 ± 0.003 and 0.169 ± 0.005 d, respectively. Applying the GCTA-GREML software, 54% of variance in the transmitting ability of 1,585 sires could be explained by considering all 40,498 markers included in the GWAS analysis. The significant markers were mainly associated with milk production genes. The SNP UA-IFASA-8854 on chromosome 11 had the lowest probability value, 1.2 × 10-24. This marker is located between the genes RETSAT and ELMOD3, both of which are overexpressed in human mammary glands. The gene RETSAT is reported to be essential for lipid accumulation and adipogenesis promotion. Gene enrichment analysis found that genes in the genomic region flanking significant markers are associated with vasopressin receptor activity, which was shown to mediate puberty in humans. If age at first insemination is included in the index with a weighting to account for 9% of the index, reductions of 2.8 and 2.6 d for age at first insemination and first calving age after 10 yr of selection are predicted, as compared with reductions of 1.4 and 1.1 d with the current index. Gains for the other index traits are only marginally affected. We suggest selection on age at first insemination as an alternative to selection for early calving.

Broad phenotypic impact of the effects of transgenerational heat stress in dairy cattle: a study of four consecutive generations.

BACKGROUND: Global warming has increased the frequency of heat stress in livestock. Although heat stress directly leads to negative effects on production and reproduction traits in dairy cattle, the transgenerational transition of these changes is poorly understood. We hypothesized that heat stress in pregnant cows might induce epigenetic modifications in the developing embryo germ cells, which, in turn, might lead to phenotypic effects in the offspring. Here, we examined whether transgenerational effects of heat stress contribute to the phenotypic expression of economic traits in Israel dairy cattle. Since heat stress in Israel occurs specifically between June and October, first we examined the association of the month of birth of F1 cows (pregnancy of the F0 dam) with the performance of the F2 and F3 female offspring. Then, we calculated an annual heat stress index and examined the association of the heat stress index during the pregnancy of the F0 dam with the performance of her F2 and F3 offspring. Finally, we examined intergenerational interactions of heat stress by comparing the performance of F3 cows according to the pregnancy seasons of the F0 and F1 animals. RESULTS: We found a significant association of the month of birth, season of pregnancy, and heat stress index of F0 females, with the performance of their F2 and F3 progenies, which suggests a true transgenerational effect. The most significant transgenerational effects were on fat yield and concentration, dystocia, still-birth, and maturation. CONCLUSIONS: These findings suggest that heat stress during pregnancy affects the performance of offspring, regardless of life circumstances in at least the last three generations. Therefore, heat stress can reduce selection efficiency in breeding programs and may have economic significance in livestock.

Genetic and environmental analysis of female calf survival in the Israel Holstein cattle population.

The objectives were to estimate the effects of various environmental factors on female calf survival of Israeli Holsteins, to estimate the economic value of calf survival under Israeli conditions, to estimate the genetic and environmental variance components for calf and cow survival using the individual animal model, to perform GWAS analyses of survival to first calving and herd life after first calving, to estimate the genetic and environmental trends for calf survival since 1985, to estimate genetic correlations of calf survival with the traits included in the current Israeli breeding index, and to estimate the consequences of inclusion of calf survival in the national selection index. Mean calf survival rate of Israeli Holsteins from 2001 through 2008 was 0.85, and the mean economic value of survival to first calving was $526. Birth month, gestation length, dystocia, and twin birth significantly affected calf survival rate. Dystocia and twin birth each reduced survival rate by 0.034. Survival rate was highest for calves born in October and lowest for calves born in February. The difference between these months was 3.4%. Maximum survival was at a gestation length of 276 d, the mean gestation length for this population. Survival rate was reduced to 0.76 for calves born after a gestation length of 260 d. The individual animal model was applied for all the genetic analyses. Heritability for calf survival to first calving, as estimated by REML, was 0.009, whereas heritability of herd life from first calving was 0.15. The complete data set for genetic analysis of survival to first calving included 1,235,815 calves born between 1985 and 2017. Annual genetic and phenotypic trends for calf survival were 0.019 and 0.015%, respectively. Correlations of transmitting abilities of 226 sires born since 2010 for calf survival with the traits included in the Israeli breeding index were significant only for the maternal effects of dystocia and stillbirth. The GWAS analysis was based on the transmitting abilities of 1,493 bulls with genotypes and reliabilities >0.5 for calf survival and cow herd life. There were 7 single nucleotide polymorphisms with coefficients of determination >0.03 for calf survival and 12 single nucleotide polymorphisms with coefficients of determination >0.05 for cow survival. There was no overlap between the genome-wide significant markers for the GWAS analyses of calf survival and cow herd life. This corresponds to the conclusion from the REML results and the low correlations between the sire evaluations that the genetic control of the 2 traits are not similar. Inclusion of calf survival in the Israeli breeding would result in a 0.5% increase in calf survival over 10 yr but reduce progress for the other traits by 8%.

Genetic and genomic analysis of long insemination interval in Israeli dairy cattle as an indicator of early abortions.

One of the causes of observed low fertility is embryo loss after fertilization. Previous findings suggested that more than half of fertilizations result in embryo loss before pregnancy is detected. We proposed reinsemination between 49 and 100 d after the first insemination as an indicator trait for early abortion (EA) in dairy cattle based on the mean estrus interval of 21 d. This trait was compared with conception rate from first insemination and conception status, computed as the inverse of the number of inseminations to conception. Animal model variance components were estimated by REML, including parents and grandparents of cows with records. First-parity heritability for first insemination conception rate was 3%. In the multitrait analysis of parities 1 to 3 for putative EA, heritabilities ranged from 8.9% for first parity to 10.4% for second parity. All genetic correlations were >0.9, whereas all environmental correlations were <0.12. The variance component for the service sire effect for putative EA rate was less than half the variance component for conception rate. Thus, genetic control of the 2 traits is clearly different, and analysis of EA rate by a single-trait animal model is justified. Genetic evaluation for putative EA was computed using this model, including all first- through third-parity cows with freshening dates from January 1, 1985, through December 31, 2016, that either became pregnant on first insemination or were reinseminated between 49 and 100 d after the first insemination. All known parents and grandparents of cows with records were included in the analysis. The regression of the breeding value for non-abortion rate on the cows' birth year was 0.083%/yr. The genetic correlation between first-parity EA and conception status was 0.995. The genetic correlations between first-parity EA and milk, fat, and protein production were all negative, whereas the genetic correlation between EA and herd life was 0.33. Inclusion of putative EA in the selection index instead of conception status resulted in 10 to 20% greater genetic gain for both fertility traits. In a genome-wide association study based on 1,200 dairy bulls with reliabilities >50% for abortion rate genotyped for 41,000 markers, 6 markers were found with nominal probabilities of <10-12 to reject the null hypothesis of no effect on EA rate. The markers with the lowest probabilities for EA rate were also included among the markers with the lowest probabilities for female fertility, but not vice versa. The marker explaining the most variance for abortion rate is located within the ABCA9 gene, which is found within an ATP-binding cassette (ABC) genes cluster. The ABC family is the major class of primary active transporters in the placenta.

Bos taurus-indicus hybridization correlates with intralocus sexual-conflict effects of PRDM9 on male and female fertility in Holstein cattle.

BACKGROUND: Crossover localization during meiotic recombination is mediated by the fast-evolving zinc-finger (ZnF) domain of gene PRDM9. To study its impact on dairy cattle performance, we compared its genetic variation between the relatively small Israeli (IL) Holsteins and the North American (US) Holsteins that count millions. RESULTS: Initially, we analyzed the major BTA1 haplotypes present in IL Holsteins based on the 10 most telomeric SNPs of the BovineSNP50 BeadChip. Sequencing of representative haplotype carriers indicated that for all frequent haplotypes (> 6%), the variable PRDM9 ZnF array consisted of seven tandem ZnF repeats. Two rare haplotypes (frequency < 4%) carried an indicine PRDM9, whereas all others were variants of the taurine type. These two haplotypes included the minor SNP allele, which was perfectly linked with a previously described PRDM9 allele known to induce unique localization of recombination hotspots. One of them had a significant (p = 0.03) negative effect on IL sire fertility. This haplotype combined the rare minor alleles of the only SNPs with significant (p < 0.05) negative substitution effects on US sire fertility (SCR). Analysis of telomeric SNPs indicated general agreement of allele frequencies (R = 0.95) and of the substitution effects on sire fertility (SCR, R = 0.6) between the US and IL samples. Surprisingly, the alleles that had a negative impact on male fertility had the most positive substitution effects on female fertility traits (DPR, CCR and HCR). CONCLUSIONS: A negative genetic correlation between male and female fertility is encoded within the BTA1 telomere. Cloning the taurine PRDM9 gene, which is the common form carried by Holsteins, we encountered the infiltration of an indicine PRDM9 variant into this population. During meiosis, in heterozygous males, the indicine PRDM9 variant may induce incompatibility of recombination hotspots and male infertility. However, this variant is associated with favorable female fertility, which would explain its survival and the general negative correlation (R = - 0.3) observed between male and female fertility in US Holsteins. Further research is needed to explain the mechanism underlying this positive effect and to devise a methodology to unlink it from the negative effect on male fertility during breeding.

Genetic and environmental analysis of diseases with major economic impact in Israeli Holsteins.

Incidences of ketosis, metritis, mastitis, and retained placenta were studied in Israeli Holstein cows calving between 2008 and 2017. These diseases were selected based on their economic impact. Ketosis, metritis, and retained placenta were scored dichotomously. Mastitis was scored as absent, a single occurrence during the lactation, or more than 1 occurrence. Ketosis and metritis were recorded during the first 21 d after calving, retained placenta during the first 5 d after calving, and mastitis up to 305 d in milk. The effects of herd-year-season, calving age, month of calving, gestation length, and occurrence of dystocia were included in the first-parity analysis models. All effects were significant for metritis and retained placenta. For ketosis, all effects were significant, except for gestation length. For mastitis, only the effects of herd-year-season and calving age were significant. Variance components were computed by the multitrait animal model. The 4 diseases were analyzed jointly based on first-parity records, and each disease was analyzed separately for parities 1 to 3 with the different parities considered separate traits. The 4 disease traits in first parity were also analyzed jointly with the 6 major traits included in the Israeli breeding index: milk, fat, and protein production; somatic cell score; female fertility; and longevity. Heritability was highest for metritis and lowest for mastitis, but all heritabilities were <0.07, similar to previous studies. For all 4 diseases, genetic correlations among the first 3 parities were >0.65, and all residual correlations were <0.07. Selection of herd-years assumed to have more accurate recording of mastitis did not result in higher heritability estimates. Genetic correlations between the disease traits and milk, fat, and protein production were economically unfavorable, while correlations between the disease traits and somatic cell score, female fertility, and longevity were economically favorable. Expected genetic changes in the disease traits after 10 yr of selection with the current Israeli breeding index were all <1%, except for ketosis, which was predicted to increase by 1.5%. Inclusion of these traits in a proposed index with the disease traits constituting 7% of the index would result in only marginal improvements for the disease traits and adversely affect genetic gain for fat and protein production. Thus, inclusion of these traits in the breeding index cannot be justified economically.

Short communication: Selection for persistency by changing the shape of lactation curves with geometric morphometrics.

The objective was to simulate the response to selection by modifying the shape of the lactation curve to improve persistency in Israeli first-parity Holstein cattle using geometric morphometrics and multivariate quantitative genetics. The 2-dimensional Cartesian coordinates representing the shape of the lactation curve were characterized using geometric morphometrics. Genetic parameters of daily milk yield were estimated for the individual animal model using REML. Results of selection were simulated using the multivariate breeders' equation. We defined 6 selection regimens: (1) decrease in the third month in milk (MIM), (2) increase in the ninth MIM, (3) decrease in the third MIM and increase in the ninth MIM, (4) decrease in the first MIM, (5) decrease in the second MIM, and (6) decrease in the fourth MIM. The first 3 selection criteria targeted changes in the MIM of previous measurements of persistency in the Israeli Holstein population, and the last 3 targeted the months with the highest variability. All selection regimens attempted to modify the shape of the lactation curve while keeping total milk yield constant. The first 3 selection regimens showed greater total responses, with the highest response for decrease in the third MIM and increase in the ninth MIM. Magnitudes of correlated responses consistently surpassed direct responses in all selection criteria, suggesting that selection on specific parts of the lactation curve would produce correlated responses in other parts. Predicted responses to selections were depicted as shape changes from the mean shape in the data set. Most predicted responses demonstrated a fairly similar shape compared with the lactation curve. The predicted responses of lactation shapes of the last 3 selection regimens were overall lower in the first 3 MIM and decreased less in the last MIM, demonstrating a much flatter shape and thus improved persistency.

Using geometric morphometrics for the genetics analysis of shape and size of lactation curves in Israeli first-parity Holstein cattle.

Our objective was to combine the methods of geometric morphometrics and multivariate quantitative genetics to genetic evaluation of the size and shape of lactation curves of milk of 3,492 Israeli first-parity Holstein cattle. Lactation records were treated as morphological data, for which 2 different lactation shape functions were evaluated, one depicted by a line graph and the other by an orbital graph. The lactation curves from both shape functions were represented by 2-dimensional Cartesian landmark coordinates. The 2 sets of landmarks were then analyzed individually for each shape function with geometric morphometrics to separate variation into components of size and shape. The analysis yielded 2 size measures and 2 sets of shape variables, and they were the inputs to estimate variance components using the MTC REML individual animal model program. Variance components were also estimated for the 305-d lactation production as a reference. Shape variables showed negligible correlation with 305-d production, providing evidence of size and shape of lactation curve as separate characters. The size measure derived from the orbital-depicted lactation curve had equal heritability (0.39 ± 0.01; ± standard error) and complete genetic and environmental correlations with 305-d production, whereas the size measure derived from the line-depicted lactation curve showed low heritability (0.09 ± 0.01) and environmental correlation (0.02 ± 0.004) and relative high genetic correlation with 305-d production (0.48 ± 0.04). This may validate both the orbital graph to depict lactation records and the use of geometric morphometrics to split variation of lactation curve into size and shape components. The maximal heritability for shape of lactation curve was 0.55 for orbital- and 0.56 for line-depicted lactation curves. The respective patterns of variations were visualized as shape changes from the mean shape in the data set. Geometric morphometrics are well grounded within the theory of shape analysis and can be paired with conventional methods in the field to characterize the patterns of phenotypic and genetic variation of shape and size of lactation curve in dairy cattle.

Genetic and Genomic Analysis of Cow Mortality in the Israeli Holstein Population.

"Livability" was defined as the inverse of the probability of death. The objectives of this study were to estimate the heritability, genetic and phenotypic trends for the livability of Israeli Holstein cows; estimate the genetic and environmental correlations between livability and the nine traits included in the Israeli breeding index; estimate the effect of the inclusion of livability in the Israeli breeding index on expected genetic gains; and compute a genome-wide association study (GWAS) for livability. Seven data sets were analyzed. All data were derived from the database of the Israeli dairy cattle herd-book. The mean livability for the complete data set of 523,954 cows born from 2000 through 2016 was 89.6%. Pregnancy reduced livability by 15%. Livability generally increased with parity and days in milk within parity. Heritability of livability was 0.0082. Phenotypic and genetic trends over the 14-year period from 2000 through 2013 were -0.42% and -0.22% per year. If livability is included in the Israeli breeding index, accounting for 9% of the index, livability would increase by 1.3% and protein production would decrease by 11 kg over the next decade, as compared to the current index. A marker in proximity to the oxytocin-vasopressin locus had the greatest effect in the GWAS. Oxytocin activity in cattle affects calving-associated pathologies and maternal death. Inclusion of livability in the Israeli breeding index is not recommended.

Economic Aspects of Bovine Ephemeral Fever (BEF) Outbreaks in Dairy Cattle Herds.

Bovine ephemeral fever virus (BEFV) is an arthropod-borne virus (arbovirus) transmitted by blood-feeding insects (mosquitoes and Culicoides biting midges). While the dispersal of arboviral diseases such as bovine ephemeral fever (BEF) into naive areas is often the result of globalization and animal movement, the endemization and local outbreaks of these diseases are mainly influenced by environmental changes. Climate change affects the activity, distribution, dynamics, and life cycles of these vectors (arthropods), the replication of viruses within their vectors, and weakens animal's immune systems. Although BEF does not currently occur in the Americas and Europe (other than in the western regions of Turkey), the risk of BEFV emergence, spread, and endemization in Europe is real. Over the past two decades, arboviruses such as the bluetongue virus (BTV) and Schmallenberg virus (SBV) have emerged in Europe without warning and caused significant losses to the dairy and meat industries. Since the European cattle population has never been exposed to BEFV, the economic losses to dairy and beef production in this continent due to the reduction in milk production, loss of valuable cows, and abortion, should BEF emerge, would probably be considerable. Moreover, arboviruses can also cause substantial financial damage due to restrictions on animal trade and transportation, like the current EHDV-8 outbreak in the Mediterranean basin. In this study, we used national data stored in the Israeli herd book to examine the economic aspects of BEF outbreaks in affected dairy cattle farms countrywide. Our results demonstrate that BEF outbreaks can have immediate and delayed effects, causing severe economic losses due to culling (loss of valuable cows) and a reduction in milk production that affects dairy farm income for months after clinical diagnosis. To our knowledge, this is the first extensive study on the impact of a BEF outbreak at a population level, enabling to conduct accurate risk assessments in future cases of BEFV emergence and re-emergence.

Estimating the Effect of the Kappa Casein Genotype on Milk Coagulation Properties in Israeli Holstein Cows.

In Israel, about 26% of produced milk is used to produce hard cheeses and 29% for soft cheeses. Milk with preferred coagulation properties requires a shorter coagulation time and yields a higher curd firmness than milk with inferior coagulation properties. Studies have shown that milk from cows with the B allele of kappa casein (κ-CN) produces more cheese than milk from those with A and E alleles. There is evidence that milk from AE or EE genotype cows is unsuitable for cheese production. In the early 1990s, the proportion of the B allele in Israeli Holstein cattle was about 17%, similar to its prevalence in the Holstein population worldwide. In recent years, however, its proportion has increased to about 40%. We analyzed milk coagulation properties as a function of the cow's κ-CN genotype, including time in minutes until the beginning of coagulation and curd firmness after 60 min-measured in volts via an optigraph device and scored on a scale of 0-4 by a laboratory technician. Cow selection was based on their sire's genotype, so that there would be sufficient genotypes that include the rare E allele. A total of 359 cows were sampled from 15 farms: 64 with genotype AA, 142 with AB, 41 with AE, 65 with BB, and 47 with BE. Data were analyzed via the general linear model procedure of SAS. We found the following: (a) There were significant differences between genotypes for optigraph-measured curd firmness. In a multi-comparison test, the BB genotype gave the highest curd firmness, and AB and BE showed a significant advantage compared to AA and AE (9.4, 8.6, 8.4, 6.9, 6.8 V, respectively). Assuming a frequency of about 55% for the A allele, about 30% of the milk delivered to dairy plants comes from AA cows. (b) There was a significant difference between the genotypes in technician-observed curd firmness, with BB scoring significantly higher than AA and AE. (c) The optigraph-measured curd firmness was significantly higher for milk from primiparous cows as compared to milk from second, third, or fourth lactation cows (8.9, 7.8, 7.9, 7.7 V, respectively). The technician-observed curd firmness was significantly higher for primiparous vs. multiparous cows. There was a clear advantage in curd firmness for genotypes that included the B allele compared to those with AA and AE genotypes. We can increase the proportion of the B allele in the population by insemination of cows using bulls with the genotypes AB and BB. This factor should therefore be included in the selection index.

Breeding Dairy Cattle for Female Fertility and Production in the Age of Genomics.

Phenotypic and genetic changes for female fertility and production traits in the Israeli Holstein population over the last three decades were studied in order to determine if long term selection has resulted in reduced heritability and negative genetic correlations. Annual means for conception status, defined as the inverse of the number of inseminations to conception in percent, decreased from 55.6 for cows born in 1983 to 46.5 for cows born in 2018. Mean estimated breeding values increased by 1.8% for cow born in 1981 to cows born in 2018. Phenotypic records from 1988 to 2016 for the nine Israeli breeding index traits were divided into three time periods for multi-trait REML analysis by the individual animal model. For all traits, heritabilities increased or stayed the same for the later time periods. Heritability for conception status was 0.05. The first parity genetic correlation between conception status and protein yield was -0.38. Heritabilities decreased with the increase in parity for protein but remained the same for conception status. Realized genetic trends were greater than expected for cows born from 2008 through 2016 for somatic cell score, conception status and herd-life, and lower than expected for the production traits.

Genetic and Genome-Wide Association Analysis of Yearling Weight Gain in Israel Holstein Dairy Calves.

Yearling weight gain in male and female Israeli Holstein calves, defined as 365 × ((weight - 35)/age at weight) + 35, was analyzed from 814,729 records on 368,255 animals from 740 herds recorded between 1994 and 2021. The variance components were calculated based on valid records from 2008 through 2017 for each sex separately and both sexes jointly by a single-trait individual animal model analysis, which accounted for repeat records on animals. The analysis model also included the square root, linear, and quadratic effects of age at weight. Heritability and repeatability were 0.35 and 0.71 in the analysis of both sexes and similar in the single sex analyses. The regression of yearling weight gain on birth date in the complete data set was -0.96 kg/year. The complete data set was also analyzed by the same model as the variance component analysis, including both sexes and accounting for differing variance components for each sex. The genetic trend for yearling weight gain, including both sexes, was 1.02 kg/year. Genetic evaluations for yearling weight gain was positively correlated with genetic evaluations for milk, fat, protein production, and cow survival but negatively correlated with female fertility. Yearling weight gain was also correlated with the direct effect on dystocia, and increased yearling weight gain resulted in greater frequency of dystocia. Of the 1749 Israeli Holstein bulls genotyped with reliabilities >50%, 1445 had genetic evaluations. As genotyping of these bulls was performed using several single nucleotide polymorhphism (SNP) chip platforms, we included only those markers that were genotyped in >90% of the tested cohort. A total of 40,498 SNPs were retained. More than 400 markers had significant effects after permutation and correction for multiple testing (pnominal < 1 × 10-8). Considering all SNPs simultaneously, 0.69 of variance among the sires' transmitting ability was explained. There were 24 markers with coefficients of determination for yearling weight gain >0.04. One marker, BTA-75458-no-rs on chromosome 5, explained ≈6% of the variance among the estimated breeding values for yearling weight gain. ARS-BFGL-NGS-39379 had the fifth largest coefficient of determination in the current study and was also found to have a significant effect on weight at an age of 13-14 months in a previous study on Holsteins. Significant genomic effects on yearling weight gain were mainly associated with milk production quantitative trait loci, specifically with kappa casein metabolism.

Analysis of copy loss and gain variations in Holstein cattle autosomes using BeadChip SNPs.

BACKGROUND: Copy number variation (CNV) has been recently identified in human and other mammalian genomes, and there is a growing awareness of CNV's potential as a major source for heritable variation in complex traits. Genomic selection is a newly developed tool based on the estimation of breeding values for quantitative traits through the use of genome-wide genotyping of SNPs. Over 30,000 Holstein bulls have been genotyped with the Illumina BovineSNP50 BeadChip, which includes 54,001 SNPs (~SNP/50,000 bp), some of which fall within CNV regions. RESULTS: We used the BeadChip data obtained for 912 Israeli bulls to investigate the effects of CNV on SNP calls. For each of the SNPs, we estimated the frequencies of occurrence of loss of heterozygosity (LOH) and of gain, based either on deviation from the expected Hardy-Weinberg equilibrium (HWE) or on signal intensity (SI) using the PennCNV "detect" option. Correlations between LOH/CNV frequencies predicted by the two methods were low (up to r = 0.08). Nevertheless, 418 locations displayed significantly high frequencies by both methods. Efficiency of designating large genomic clusters of olfactory receptors as CNVs was 29%. Frequency values for copy loss were distinguishable in non-autosomal regions, indicating misplacement of a region in the current BTA7 map. Analysis of BTA18 placed major quantitative trait loci affecting net merit in the US Holstein population in regions rich in segmental duplications and CNVs. Enrichment of transporters in CNV loci suggested their potential effect on milk-production traits. CONCLUSIONS: Expansion of HWE and PennCNV analyses allowed estimating LOH/CNV frequencies, and combining the two methods yielded more sensitive detection of inherited CNVs and better estimation of their possible effects on cattle genetics. Although this approach was more effective than methodologies previously applied in cattle, it has severe limitations. Thus the number of CNVs reported here for the Holstein breed may represent as little as one-tenth of inherited common structural variation.

Factors affecting fertility of dairy cows in Israel.

UNLABELLED: The objectives of this review are to describe the reproductive parameters monitored in Israeli dairy herds and to evaluate their changes in recent years. Eighty percent of the cows and 70% of the farms use the Israel Cattle Breeders' Association Herdbook and about 50% of them use pedometry systems. Intensive herd medicine is practiced in 80% of the herds by Hachaklait Veterinary Services Ltd. Herd-health reports monitor calving, production and reproduction. Causal analysis explains the effects and interactions of various risk factors involved. The average of 305 days of milk production per cow increased between 2004 and 2008 from 11,200 to 11,903 kg. At the same time the first A.I. conception rate (C.R) dropped from 43.0 to 40.7% and from 35.6 to 30.5% in primiparous cows (PC) and multiparous cows (MC), respectively. The waiting period (WP) was shortened from 106.2 to 93.4 days in PC and from 99.9 to 87.3 days in MC. The undetected heat rate per herd increased from 30.3 to 38.9% and from 33.9 to 43.9% in PC and MC, respectively. The average of days open per herd dropped from 127 to 118.4 and from 127.5 to 120.5 in PC and MC, respectively. The rate of cows open by 150 days in lactation dropped from 42% (+/- 10.2) to 34.2% (+/- 8.1) and 47.1% (+/- 8.8) to 39.5% (+/- 7.1) in PC and MC, respectively. The ratio between summer inseminations and winter inseminations increased from 0.81 to 1.04 from 2000 to 2008. The calving interval (CI) average fluctuated around 424.5 (+/- 2.0) days and 417.5 (+/- 1.7) days in PC and MC, respectively. The average duration of the dry period in 2008 was 60.7 (+/- 4.7, 47-72) days. From 2004 to 2008, the average herd rate of endometritis increased from 38.1 to 46.0% and from 25.5 to 30.1% in PC and MC, respectively. The milk fat to protein ratio in the first test day of lactation has remained steady during the past 5 years. Genetic trends in the breeding values of fertility and milk showed consistent improvement from 2000 to 2006. CONCLUSIONS: In recent years there has been a small decline in some reproductive parameters, while at the same time others have remained unchanged. The farmer's economical viewpoint and management practices have contributed to the changes.