Bayesian inference of the inbreeding load variance for fertility traits in Brown Swiss cattle.

Our study investigated the inbreeding load for fertility traits in the Italian Brown Swiss dairy cattle breed. Fertility traits included continuous traits (i.e., interval from calving to first service, days open, and calving interval) and categorical traits (i.e., calving rate at first insemination and nonreturn date at d 56). We included only records of the first 3 parities of cows that calved between 2010 and 2018. We traced up the pedigree of the cows with records as far as possible, ending up with a total of 73,246 animals. The final data set consisted of 59,864 records from 34,921 cows. We analyzed all models using a Bayesian approach that included a covariate with total inbreeding in addition to systematic, permanent environment, additive genetic, and inbreeding load effects. We then evaluated the trends in heritabilities and ratios of the inbreeding load using a continuum of partial inbreeding coefficients from 0.001 to 0.100 as reference. Posterior estimates of heritabilities tended to decrease across the continuum, whereas ratios of the inbreeding load tended to increase, more noticeably in categorical traits (calving rate at first insemination and nonreturn date at d 56). From the results obtained, we confirmed the presence of heterogeneity in inbreeding depression. We then predicted the inbreeding load effects, which had a low reliability of prediction, explained by having only 513 ancestors generating inbreeding. However, reliability of prediction was high enough for some of the individuals, obtaining a favorable prediction of inbreeding load for a relevant percentage, which improved the phenotypic performance of their inbred descendants. These results make it feasible to implement breeding and management strategies that select ancestors with a favorable inbreeding load prediction. In addition, it opens the possibility to define a global index for the expected consequences of the inbreeding generated by each individual.

Evaluation of bull fertility in Italian Brown Swiss dairy cattle using cow field data.

Dairy bull fertility is traditionally evaluated using semen production and quality traits; however, these attributes explain only part of the differences observed in fertility among bulls. Alternatively, bull fertility can be directly evaluated using cow field data. The main objective of this study was to investigate bull fertility in the Italian Brown Swiss dairy cattle population using confirmed pregnancy records. The data set included a total of 397,926 breeding records from 1,228 bulls and 129,858 lactating cows between first and fifth lactation from 2000 to 2019. We first evaluated cow pregnancy success, including factors related to the bull under evaluation, such as bull age, bull inbreeding, and AI organization, and factors associated with the cow that receives the dose of semen, including herd-year-season, cow age, parity, and milk yield. We then estimated sire conception rate using only factors related to the bull. Model predictive ability was evaluated using 10-fold cross-validation with 10 replicates. Interestingly, our analyses revealed that there is a substantial variation in conception rate among Brown Swiss bulls, with more than 20% conception rate difference between high-fertility and low-fertility bulls. We also showed that the prediction of bull fertility is feasible as our cross-validation analyses achieved predictive correlations equal to 0.30 for sire conception rate. Improving reproduction performance is one of the major challenges of the dairy industry worldwide, and for this, it is essential to have accurate predictions of service sire fertility. This study represents the foundation for the development of novel tools that will allow dairy producers, breeders, and artificial insemination companies to make enhanced management and selection decisions on Brown Swiss male fertility.

The genetic heritage of Alpine local cattle breeds using genomic SNP data.

BACKGROUND: Assessment of genetic diversity and population structure provides important control metrics to avoid genetic erosion, inbreeding depression and crossbreeding between exotic and locally-adapted cattle breeds since these events can have deleterious consequences and eventually lead to extinction. Historically, the Alpine Arc represents an important pocket of cattle biodiversity with a large number of autochthonous breeds that provide a fundamental source of income for the entire regional economy. By using genotype data from medium-density single nucleotide polymorphism (SNP) arrays, we performed a genome-wide comparative study of 23 cattle populations from the Alpine Arc and three cosmopolitan breeds. RESULTS: After filtering, we obtained a final genotyping dataset consisting of 30,176 SNPs for 711 individuals. The local breeds showed high or intermediate values of genetic diversity compared to the highly selected cosmopolitan breeds. Patterns of genetic differentiation, multidimensional scaling, admixture analysis and the constructed phylogenetic tree showed convergence, which indicates the presence of gene flow among the breeds according to both geographic origin and historical background. Among the most differentiated breeds, we identified the modern Brown cattle. In spite of admixture events, several local breeds have preserved distinctive characteristics, which is probably due to differences in genetic origin and geographic location. CONCLUSIONS: This study represents one of the most comprehensive genome-wide analysis of the Alpine cattle breeds to date. Using such a large dataset that includes the majority of the local breeds found in this region, allowed us to expand knowledge on the evaluation and status of Alpine cattle biodiversity. Our results indicate that although many of the analyzed local breeds are listed as endangered, they still harbor a large amount of genetic diversity, even when compared to some cosmopolitan breeds. This finding, together with the reconstruction of the phylogeny and the relationships between these Alpine Arc cattle breeds, provide crucial insights not only into the improvement of genetic stocks but also into the implementation of future conservation strategies.

New Insights in Cheese Yield Capacity of the Milk of Italian Brown and Italian Friesian Cattle in the Production of High-Moisture Mozzarella.

The aim of the present study is to investigate the effect of κ-casein B content in milk on the yield of high-moisture mozzarella cheese. The study was carried out by monitoring the production of eight mozzarella cheese batches at four cheese making factories. At each factory, two cheese making trials were performed in parallel: one using bulk milk from Italian Brown cattle and the other using bulk milk from Italian Friesian cattle. The average κ-casein B content was 0.04 g per 100 g in the Italian Friesian cows' milk, whereas it was four time higher in the Italian Brown cows' milk, reaching values of 0.16 g per 100 g. Both the κ-casein content and κ-casein B to casein ratio were positively correlated with actual cheese yield. Both parameters showed correlation coefficient values over 0.9, higher than for any other protein fraction. The influence of the level of κ-casein on the increase of the yield is probably due to smaller and more homogeneous micelles, with more efficient rennet coagulation. Consequently, milk with higher κ-casein B content produces a more elastic curd that withstands better the technological treatments and limits losses during curd mincing and stretching. In conclusion, the Italian Brown cows' milk used, characterized by higher κ-casein content than the Italian Friesian's one, allowed a yield increase of about 2.65%, which is a very relevant result for both farms and cheese making factories.

A practical approach to detect ancestral haplotypes in livestock populations.

BACKGROUND: The effects of different evolutionary forces are expected to lead to the conservation, over many generations, of particular genomic regions (haplotypes) due to the development of linkage disequilibrium (LD). The detection and identification of early (ancestral) haplotypes can be used to clarify the evolutionary dynamics of different populations as well as identify selection signatures and genomic regions of interest to be used both in conservation and breeding programs. The aims of this study were to develop a simple procedure to identify ancestral haplotypes segregating across several generations both within and between populations with genetic links based on whole-genome scanning. This procedure was tested with simulated and then applied to real data from different genotyped populations of Spanish, Fleckvieh, Simmental and Brown-Swiss cattle. RESULTS: The identification of ancestral haplotypes has shown coincident patterns of selection across different breeds, allowing the detection of common regions of interest on different bovine chromosomes and mirroring the evolutionary dynamics of the studied populations. These regions, mainly located on chromosomes BTA5, BTA6, BTA7 and BTA21 are related with certain animal traits such as coat colour and milk protein and fat content. CONCLUSION: In agreement with previous studies, the detection of ancestral haplotypes provides useful information for the development and comparison of breeding and conservation programs both through the identification of selection signatures and other regions of interest, and as indicator of the general genetic status of the populations.

Predicting haplotype carriers from SNP genotypes in Bos taurus through linear discriminant analysis.

BACKGROUND: SNP (single nucleotide polymorphisms) genotype data are increasingly available in cattle populations and, among other things, can be used to predict carriers of specific haplotypes. It is therefore convenient to have a practical statistical method for the accurate classification of individuals into carriers and non-carriers. In this paper, we present a procedure combining variable selection (i.e. the selection of predictive SNPs) and linear discriminant analysis for the identification of carriers of a haplotype on BTA19 (Bos taurus autosome 19) known to be associated with reduced cow fertility. A population of 3645 Brown Swiss cows and bulls genotyped with the 54K SNP-chip was available for the analysis. RESULTS: The overall error rate for the prediction of haplotype carriers was on average very low (∼≤1%). The error rate was found to depend on the number of SNPs in the model and their density around the region of the haplotype on BTA19. The minimum set of SNPs to still achieve accurate predictions was 5, with a total test error rate of 1.59. CONCLUSIONS: The paper describes a procedure to accurately identify haplotype carriers from SNP genotypes in cattle populations. Very few misclassifications were observed, which indicates that this is a very reliable approach for potential applications in cattle breeding.

Variation of milk components in the Italian Brown cattle.

The aim of this study was to evaluate the variations of protein, casein, saturated (SFA), unsaturated (UFA), monounsaturated (MUFA), polyunsaturated (PUFA) fatty acids contents and cheese yield in the milk of two groups of Italian Brown cows conventionally reared in indoor period of housing or consuming pasture during the summer months in 2008 and 2013. Milk components were obtained from samples collected during the national routine (conventionally reared) and 'extraordinary' (pasture period) milk recording scheme in herds located near Sondrio (Lombardia, Italy). Milk samples were processed with the MilkoScanTM FT6000 for the identification of milk casein, SFA, UFA, MUFA and PUFA composition. The groups were analysed separately per year and the environmental factors affecting milk protein, casein, and fatty acids contents (pasture/indoor, parity, data of sampling, days in milk, days from collection to analysis) were included in the MIXED procedure of SAS 9.3. A total of 778 milk samples were available, including 234 records from indoor and 544 observations from pasture feeding. Pasture intake affected the content of casein (%) and the proportion of fat in milk (g/100 g), enhancing milk casein levels (from 2.90 to 3) and reducing the concentration of milk SFA in milk from grazing cows (from 2.29 to 1.92). Additionally, the cheese yield was calculated as 'kg of cheese per 100 kg of milk' and resulted to be 10.4 and 12 in 2008 from milk of cows reared indoor and with pasture based diet, respectively. The dairy industry should take advantage of the milk production during grazing periods from which high quality products may be obtained.

Identification of runs of homozygosity associated with male fertility in Italian Brown Swiss cattle.

Intensive selection for improved productivity has been accompanied by an increase in inbreeding rates and a reduction in genetic diversity. The increase in inbreeding tends to impact performance, especially fitness-related traits such as male fertility. Inbreeding can be monitored using runs of homozygosity (ROH), defined as contiguous lengths of homozygous genotypes observed in an individual's chromosome. The goal of this study was to evaluate the presence of ROH in Italian Brown Swiss cattle and assess its association with bull fertility. First, we evaluated the association between ROH and male fertility using 1,102 Italian Brown Swiss bulls with sire conception rate records and 572 K SNPs spanning the entire genome. Second, we split the entire population into 100 high-fertility and 100 low-fertility bulls to investigate the potential enrichment of ROH segments in the low-fertility group. Finally, we mapped the significant ROH regions to the bovine genome to identify candidate genes associated with sperm biology and male fertility. Notably, there was a negative association between bull fertility and the amount of homozygosity. Four different ROH regions located in chromosomes 6, 10, 11, and 24 were significantly overrepresented in low-fertility bulls (Fisher's exact test, p-value <0.01). Remarkably, these four genomic regions harbor many genes such as WDR19, RPL9, LIAS, UBE2K, DPF3, 5S-rRNA, 7SK, U6, and WDR7 that are related to sperm biology and male fertility. Overall, our findings suggest that inbreeding and increased homozygosity have a negative impact on male fertility in Italian Brown Swiss cattle. The quantification of ROH can contribute to minimizing the inbreeding rate and avoid its negative effect on fitness-related traits, such as male fertility.

Estimates of marker effects for measures of milk flow in the Italian brown Swiss dairy cattle population.

BACKGROUND: Milkability is a complex trait that is characterized by milk flow traits including average milk flow rate, maximum milk flow rate and total milking time. Milkability has long been recognized as an economically important trait that can be improved through selection. By improving milkability, management costs of milking decrease through reduced labor and improved efficiency of the automatic milking system, which has been identified as an important factor affecting net profit. The objective of this study was to identify markers associated with electronically measured milk flow traits, in the Italian Brown Swiss population that could potentially improve selection based on genomic predictions. RESULTS: Sires (n = 1351) of cows with milk flow information were genotyped for 33,074 single nucleotide polymorphism (SNP) markers distributed across 29 Bos taurus autosomes (BTA). Among the six milk flow traits collected, ascending time, time of plateau, descending time, total milking time, maximum milk flow and average milk flow, there were 6,929 (time of plateau) to 14,585 (maximum milk flow) significant SNP markers identified for each trait across all BTA. Unique regions were found for each of the 6 traits providing evidence that each individual milk flow trait offers distinct genetic information about milk flow. This study was also successful in identifying functional processes and genes associated with SNPs that influences milk flow. CONCLUSIONS: In addition to verifying the presence of previously identified milking speed quantitative trait loci (QTL) within the Italian Brown Swiss population, this study revealed a number of genomic regions associated with milk flow traits that have never been reported as milking speed QTL. While several of these regions were not associated with a known gene or QTL, a number of regions were associated with QTL that have been formerly reported as regions associated with somatic cell count, somatic cell score and udder morphometrics. This provides further evidence of the complexity of milk flow traits and the underlying relationship it has with other economically important traits for dairy cattle. Improved understanding of the overall milking pattern will aid in identification of cows with lower management costs and improved udder health.

Deciphering the genetic basis of male fertility in Italian Brown Swiss dairy cattle.

Improving reproductive performance remains a major goal in dairy cattle worldwide. Service sire has been recognized as an important factor affecting herd fertility. The main objective of this study was to reveal the genetic basis of male fertility in Italian Brown Swiss dairy cattle. Dataset included 1102 Italian Brown Swiss bulls with sire conception rate records genotyped with 454k single nucleotide polymorphisms. The analysis included whole-genome scans and gene-set analyses to identify genomic regions, individual genes and genetic mechanisms affecting Brown Swiss bull fertility. One genomic region on BTA1 showed significant additive effects. This region harbors gene RABL3 which is implicated cell proliferation and motility. Two genomic regions, located on BTA6 and BTA26, showed marked non-additive effects. These regions harbor genes, such as WDR19 and ADGRA1, that are directly involved in male fertility, including sperm motility, acrosome reaction, and embryonic development. The gene-set analysis revealed functional terms related to cell adhesion, cellular signaling, cellular transport, immune system, and embryonic development. Remarkably, a gene-set analysis also including Holstein and Jersey data, revealed significant processes that are common to the three dairy breeds, including cell migration, cell-cell interaction, GTPase activity, and the immune function. Overall, this comprehensive study contributes to a better understanding of the genetic basis of male fertility in cattle. In addition, our findings may guide the development of novel genomic strategies for improving service sire fertility in Brown Swiss cattle.

Integration of Wet-Lab Measures, Milk Infrared Spectra, and Genomics to Improve Difficult-to-Measure Traits in Dairy Cattle Populations.

The objective of this study was to evaluate the contribution of Fourier-transformed infrared spectroscopy (FTIR) data for dairy cattle breeding through two different approaches: (i) estimating the genetic parameters for 30 measured milk traits and their FTIR predictions and investigating the additive genetic correlation between them and (ii) evaluating the effectiveness of FTIR-derived phenotyping to replicate a candidate bull's progeny testing or breeding value prediction at birth. Records were available from 1,123 cows phenotyped using gold standard laboratory methodologies (LAB data). This included phenotypes related to fine milk composition and milk technological characteristics, milk acidity, and milk protein fractions. The dataset used to generate FTIR predictions comprised 729,202 test-day records from 51,059 Brown Swiss cows (FIELD data). A first approach consisted of estimating genetic parameters for phenotypes available from LAB and FIELD datasets. To do so, a set of bivariate animal models were run, and genetic correlations between LAB and FIELD phenotypes were estimated using FIELD information obtained at the population level. Heritability estimates were generally higher for FIELD predictions than for the corresponding LAB measures. The additive genetic correlations (r a ) between LAB and FIELD phenotypes had different magnitudes across traits but were generally strong. Overall, these results demonstrated the potential of using FIELD information as indicator traits for the indirect genetic improvement of LAB measures. In the second approach, we included genotype information for 1,011 cows from the LAB dataset, 1,493 cows from the FIELD dataset, 181 sires with daughters in both LAB and FIELD datasets, and 540 sires with daughters in the FIELD dataset only. Predictions were obtained using the single-step GBLUP method. A four fold cross-validation was used to assess the predictive ability of the different models, assessed as the ability to predict masked LAB records from daughters of progeny testing bulls. The correlation between observed and predicted LAB measures in validation was averaged over the four training-validation sets. Different sets of phenotypic information were used sequentially in cross-validation schemes: (i) LAB cows from the training set; (ii) FIELD cows from the training set; and (iii) FIELD cows from the validation set. Models that included FIELD records showed an improvement for the majority of traits. This study suggests that breeding programs for difficult-to-measure traits could be implemented using FTIR information. While these programs should use progeny testing, acceptable values of accuracy can be achieved also for bulls without phenotyped progeny. Robust calibration equations are, deemed as essential.