Genomic inbreeding coefficients using imputation genotypes: assessing the effect of ancestral genotyping in Holstein-Friesian dairy cows.

The objective of this study was to assess the effect of using or not the genotypes of the parents of a cow for imputing single nucleotide polymorphisms (SNP), on the estimation of genomic inbreeding coefficients of cows. Imputation (i.e., genotyped plus imputed) genotypes from 68,127 Italian Holstein dairy cows registered in the Italian National Association of Holstein, Brown and Jersey Breeders (ANAFIBJ) were analyzed. Cows were genotyped with the HD Illumina Infinium BovineHD BeadChip and GeneSeek Genomic Profiler HD-150K, and the MD GeneSeek Genomic Profiler 3, GeneSeek Genomic Profiler 4, GeneSeek MD and the Labogena MD. To assess differences among estimators genomic inbreeding coefficients were estimated with 4 PLINK v1.9 estimators (F, Fhat1, 2, 3), 2 genomic relationship matrix (grm) based estimators (Fgrm and Fgrm2; with the latter including also pedigree information) and one estimator of runs of homozygosity (ROH; FROH). Assuming that the correct genomic inbreeding coefficients should be those estimated from genotyped SNP, a comparison of the genomic inbreeding coefficients estimated either with the genotyped SNP or the SNP after imputation was made. Information on the presence or absence of genotypic information from sire, dam and maternal grandsire during the imputation was investigated. Genomic inbreeding coefficients estimated with genotyped SNP or SNP after imputation were consistent for F, Fhat3, Fgrm2 and FROH, when at least one of the parents was genotyped. Biased (mainly higher) genomic inbreeding coefficients of imputation SNP were observed in cows that were genotyped with MD SNP panels whose SNP were poorly represented in the selected imputation SNP data set and also did not have their parents genotyped compared with what expected based on actual genotype data. For cows genotyped with MD the estimators Fhat1, Fhat2 and Fgrm provided higher genomic inbreeding coefficients of imputation SNP even with both parents and the maternal grandsire genotyped. Overall, FROH was the most robust estimator, followed by F and Fhat3. Our findings suggest that SNP selection, parental genotyping and estimator should be considered for designing imputation strategies in dairy cattle for estimating genomic inbreeding with imputation SNP. For computing genomic inbreeding coefficients, it is recommendable to have at least one parent genotyped and use an ROH based estimator.

Fecal microbiota characterization of an Italian local horse breed.

The Bardigiano horse is a traditional native Italian breed with a rich history and peculiar characteristics. Local breeds are proven to have unique genetic traits developed over generations to adapt to defined geographical regions and/or conditions. The specific microbial communities that coexist within these animals are unraveled by studying their microbiota, which permits a further step in the characterization of local heritage. This work aimed to characterize Bardigiano horse fecal microbiota composition. The data obtained were then compared with published data of a mix of athlete breeds to evaluate potential differences among local and specialized breeds. The study involved 11 Bardigiano mares between 3 and 4 years of age, from which stool was sampled for the study. Samples were processed for 16S rRNA sequencing. Data obtained were analyzed and plotted using R, RStudio, and FastTree software. The samples analyzed were similar to what literature has reported on horses of other breeds and attitudes at higher taxonomic levels (from phylum to genera). While at lower taxonomic levels, the difference was more marked highlighting specific families found in the Bardigiano breed only. Weight, province of origin, and breeding sites significantly affected microbiota composition (p-value ≤0.02, p-value ≤0.04, and p-value ≤0.05, respectively). The comparison with athlete breed showed a significant difference confirming that animal and environmental factors are crucial in determining fecal microbiota composition (p-value <0.001). Understanding the microbiota composition in local breeds like the Bardigiano horse is crucial for preserving biodiversity, managing animal health, and promoting sustainable farming practices.

Microsatellite imputation using SNP data for parentage verification in four Italian sheep breeds.

Microsatellite markers (MS) have been widely used for parentage verification in most of the livestock species over the past decades mainly due to their high polymorphic information content. In the genomic era, the spread of genotype information as single-nucleotide polymorphism (SNP) has raised the question to effectively use SNPs also for parentage testing. Despite the clear advantages of SNP panels in terms of cost, accuracy, and automation, the transition from MS to SNP markers for parentage verification is still very slow and, so far, only routinely applied in cattle. A major difficulty during this transition period is the need of SNP data for parents and offspring, which in most cases is not yet feasible due to the genotyping cost. To overcome the unavailability of same genotyping platform during the transition period, in this study we aimed to assess the feasibility of a MS imputation pipeline from SNPs in four native sheep dairy breeds: Comisana (N = 331), Massese (N = 210), Delle Langhe (N = 59) and Sarda (N = 1003). Those sheep were genotyped for 11 MS and with the Ovine SNP50 Bead Chip. Prior to imputation, a quality control (QC) was performed, and SNPs located within a window of 2 Mb from each MS were selected. The core of the developed pipeline was made up of three steps: (a) storing of both MS and SNP data in a Variant Call Format file, (b) masking MS information in a random sample of individuals (10%), (c) imputing masked MS based on non-missing individuals (90%) using an imputation program. The feasability of the proposed methodology was assessed also among different training - testing split ratio, population size, number of flanking SNPs as well as within and among breeds. The accuracy of the MS imputation was assessed based on the genotype concordance as well as at parentage verification level in a subset of animals in which assigned parents' MS were available. A total of 8 MS passed the QC, and 505 SNPs were located within the ±2 Mb window from each MS, with an average of 63 SNPs per MS. The results were encouraging since when excluding the worst imputed MS (OARAE129), and regardless on the analyses performed (within and across breeds) for all breeds, we achieved an overall concordance rate over 94%. In addition, on average, the imputed offspring MS resulted in equivalent parentage outcome in 94% of the cases when compared to verification using original MS, highlighting both the feasibility and the eventual practical advantage of using this imputation pipeline.

The use of milk Fourier-transform infrared spectra for predicting cheesemaking traits in Grana Padano Protected Designation of Origin cheese.

The prediction of the cheese yield (%CY) traits for curd, solids, and retained water and the amount of fat, protein, solids, and energy recovered from the milk into the curd (%REC) by Bayesian models, using Fourier-transform infrared spectroscopy (FTIR), can be of significant economic interest to the dairy industry and can contribute to the improvement of the cheese process efficiency. The yields give a quantitative measure of the ratio between weights of the input and output of the process, whereas the nutrient recovery allows to assess the quantitative transfer of a component from milk to cheese (expressed in % of the initial weight). The aims of this study were: (1) to investigate the feasibility of using bulk milk spectra to predict %CY and %REC traits, and (2) to quantify the effect of the dairy industry and the contribution of single-spectrum wavelengths on the prediction accuracy of these traits using vat milk samples destined to the production of Grana Padano Protected Designation of Origin cheese. Information from 72 cheesemaking days (in total, 216 vats) from 3 dairy industries were collected. For each vat, the milk was weighed and analyzed for composition (total solids [TS], lactose, protein, and fat). After 48 h from cheesemaking, each cheese was weighed, and the resulting whey was sampled for composition as well (TS, lactose, protein, and fat). Two spectra from each milk sample were collected in the range between 5,011 and 925 cm-1 and averaged before the data analysis. The calibration models were developed via a Bayesian approach by using the BGLR (Bayesian Generalized Linear Regression) package of R software. The performance of the models was assessed by the coefficient of determination (R2VAL) and the root mean squared error (RMSEVAL) of validation. Random cross-validation (CVL) was applied [80% calibration and 20% validation set] with 10 replicates. Then, a stratified cross-validation (SCV) was performed to assess the effect of the dairy industry on prediction accuracy. The study was repeated using a selection of informative wavelengths to assess the necessity of using whole spectra to optimize prediction accuracy. Results showed the feasibility of using FTIR spectra and Bayesian models to predict cheesemaking traits. The R2VAL values obtained with the CVL procedure were promising in particular for the %CY and %REC for protein, ranging from 0.44 to 0.66 with very low RMSEVAL (from 0.16 to 0.53). Prediction accuracy obtained with the SCV was strongly influenced by the dairy factory industry. The general low values gained with the SCV do not permit a practical application of this approach, but they highlight the importance of building calibration models with a dataset covering the largest possible sample variability. This study also demonstrated that the use of the full FTIR spectra may be redundant for the prediction of the cheesemaking traits and that a specific selection of the most informative wavelengths led to improved prediction accuracy. This could lead to the development of dedicated spectrometers using selected wavelengths with built-in calibrations for the online prediction of these innovative traits.

The role of inbreeding depression on productive performance in the Italian Holstein breed.

Inbreeding depression has become an urgent issue in cosmopolitan breeds where the massive genetic progress achieved in the latest generations is counterbalanced by a dramatic loss of genetic diversity causing increased health issues. Thus, the aim of this study was to estimate inbreeding depression on productive traits in Holstein dairy cattle. More precisely, we aimed to i) determine the level of inbreeding in 27,735 Italian Holstein dairy cows using pedigree and genotype data, ii) quantify the effect of inbreeding on 305-d in milk yield (MY; kg), fat yield (FY; kg), and protein yield (PY; kg) based on different statistical approaches, iii) determine if recent inbreeding has a more harmful impact than ancestral ones, and iv) quantify chromosomal homozygosity effect on productive traits. Quality control was performed on the autosomal chromosomes resulting in a final dataset of 84,443 single nucleotide polymorphisms. Four statistical models were used to evaluate the presence of inbreeding depression, which included linear regression analysis and division of FPED and FROH into percentile classes. Moreover, FROH was partitioned into i) length classes to assess the role of recent and ancestral inbreeding and ii) chromosome-specific contributions (FROH-CHR). Results evidenced that inbreeding negatively impacted the productive performance of Italian Holstein Friesian cows. However, differences between the estimated FPED and FROH coefficients resulted in different estimates of inbreeding depression. For instance, a 1% increase in FPED and FROH was associated with a decrease in MY of about 44 and 61 kg (P < 0.01). Further, when considering the extreme inbreeding percentile classes moving from the 5th lowest to the 95th highest, there was a reduction of -263 kg and -561 kg per lactation for FPED and FROH. Increased inbreeding, estimated by FPED and FROH, had also a negative effect on PY and FY, either fit as a regressor or percentile classes. When evaluating the impact of inbreeding based on runs of homozygosity (ROH) length classes, longer ROH (over 8 Mb) had a negative effect in all traits, indicating that recent inbreeding might be more harmful than the ancestral one. Finally, results within chromosome homozygosity highlighted specific chromosomes with a more deleterious effect on productive traits.

Inbreeding depression is a reduction in performance or health due to the mating of closely related individuals. The overall aim of this study was to investigate the level of inbreeding in the Italian Holstein dairy cow breed and quantify its negative effect on productive performances. The level of inbreeding was estimated by pedigree (FPED) and genomic data by looking at stretches of homozygosity (FROH). Both methods revealed a reduction in milk yield, fat yield, and protein yield when inbreeding increased. Moreover, the study demonstrated that FROH was able to capture more inbreeding depression compared to FPED. In addition, the more recent inbreeding had a stronger negative impact on productive performances compared to ancestral ones. Then, since the amount of runs of homozygosity can vary across the chromosomes of an individual, the effect of each chromosomal homozygosity region on productive traits was also evaluated. The chromosome-level results might be included in breeding programs to limit the accumulation of homozygosity in particular regions that appear to have a more detrimental effect on productive traits. Overall, this study highlights the importance of avoiding inbreeding in animal breeding programs to keep productive animals in the long term.

Heritability and genetic correlations of total and differential somatic cell count with milk yield and composition traits in Italian Simmental cows.

Costs of production have deeply increased each year in the last decades, breeders are continuously looking for more cost effective and more efficient ways to produce milk. Despite the major signs of progress in productivity, it is fundamental to optimize rather than maximize the performances of the dairy cows. Mastitis is still a highly prevalent disease in the dairy sector which causes several economic losses and environmental effect. Its accurate and early diagnosis is crucial to improve profitability of dairy cows and contribute to a more sustainable dairy industry. Among mastitis reduction strategies, there is the urgent need to implement breeding objectives to select cows displaying mastitis resistance by investigating the genetic mechanisms at the base of the inflammatory response. Therefore, in this study we aimed to further understand the genetic background of the differential somatic cell count (DSCC), which provides thorough insights on the actual inflammatory status of the mammary glands. The objectives of this study were to estimate on a cohort of 20,215 Italian Simmental cows over a 3-yr period: (1) the heritability and repeatability values of somatic cell score (SCS) and DSCC, (2) the genetic and phenotypic correlations between these 2 traits and milk production and milk composition traits, (3) the heritability and repeatability values of SCS and DSCC within class of udder health status. Heritability was low both for SCS (0.06) and DSCC (0.08), whereas the repeatability values for these traits were 0.43 and 0.36, suggesting that the magnitude of cow permanent environmental effect for these traits is remarkable. The genetic and phenotypic correlation of SCS with DSCC was 0.612 and 0.605, respectively. Because both significantly differed from the unit, we must consider those traits as different ones. This latter aspect corroborates the need to consider the DSCC as a further indicator of inflammatory status which might be implemented in the Simmental breed genetic evaluation. It is worthy to mention that heritability estimates for SCS and DSCC were the highest in healthy cows compared with the other udder health classes. This implies that when the udder health status changes, it is most likely due to environmental factors rather than aspects related to the animal's genetics. In contrast, the highest additive genetic variance and heritability found for SCS and DSCC in the healthy group might reveal the potential to further implement breeding strategies to select for healthier animals.

Genomic inbreeding coefficients using imputed genotypes: assessing differences among SNP panels in Holstein-Friesian dairy cows.

The objective of this study was to evaluate the effect of imputation of single nucleotide polymorphisms (SNP) on the estimation of genomic inbreeding coefficients. Imputed genotypes of 68,127 Italian Holstein dairy cows were analyzed. Cows were initially genotyped with two high density (HD) SNP panels, namely the Illumina Infinium BovineHD BeadChip (678 cows; 777,962 SNP) and the Genomic Profiler HD-150K (641 cows; 139,914 SNP), and four medium density (MD): GeneSeek Genomic Profiler 3 (10,679 cows; 26,151 SNP), GeneSeek Genomic Profiler 4 (33,394 cows; 30,113 SNP), GeneSeek MD (12,030 cows; 47,850 SNP) and the Labogena MD (10,705 cows; 41,911 SNP). After imputation, all cows had genomic information on 84,445 SNP. Seven genomic inbreeding estimators were tested: (i) four PLINK v1.9 estimators (F, Fhat1,2,3), (ii) two genomic relationship matrix (grm) estimators [VanRaden's 1st method, but with observed allele frequencies (Fgrm) and VanRaden's 3rd method that is allelic free and pedigree dependent (Fgrm2)], and (iii) a runs of homozygosity (roh) - based estimator (Froh). Genomic inbreeding coefficients of each SNP panel were compared with genomic inbreeding coefficients derived from the 84,445 imputation SNP. Coefficients of the HD SNP panels were consistent between genotyped-imputed SNP (Pearson correlations ~99%), while variability across SNP panels and estimators was observed in the MD SNP panels, with Labogena MD providing, on average, more consistent estimates. The robustness of Labogena MD, can be partly explained by the fact that 97.85% of the SNP of this panel is included in the 84,445 SNP selected by ANAFIBJ for routine genomic imputations, while this percentage for the other MD SNP panels varied between 55 and 60%. Runs of homozygosity was the most robust estimator. Genomic inbreeding estimates using imputation SNP are influenced by the SNP number of the SNP panel that are included in the imputed SNP, and performance of genomic inbreeding estimators depends on the imputation.

Distribution of Calcium, Phosphorus and Magnesium in Yak (Bos grunniens) Milk from the Qinghai Plateau in China.

This research was aimed to assess the distribution of calcium, phosphorus and magnesium within the casein micelles of yak milk. To this aim, nine bulk yak milk samples (Y-milk), collected in three yak farms located in the Chinese province of Qinghai, were compared to nine bulk cow milk samples used as a reference. A quite similar content of colloidal calcium (0.80 vs. 0.77 mmol/g of casein; p > 0.05), a higher content of magnesium (0.05 vs. 0.04 mmol/g of casein; p ≤ 0.01) and a lower content of colloidal phosphorus (0.48 vs. 0.56 mmol/g of casein; p ≤ 0.01) between yak and cow casein micelles were found. Moreover, the yak casein micelles showed a lower value of prosthetic phosphorus (0.20 vs. 0.26 mmol/g of casein; p ≤ 0.05) compared to the cow micelles. The lower values of colloidal and prosthetic phosphorus in yak casein micelles suggest that the yak casein is less phosphorylated than the cow one.

Fatty acids composition and lipolysis of Parmigiano Reggiano PDO cheese: effect of the milk cooling temperature at the farm.

OBJECTIVE: The aim was to study the influence of cooling milk at 9°C at the farm versus keeping it at 20°C on Parmigiano Reggiano cheese lipolysis. METHODS: A total of six cheesemaking trials (3 in winter and 3 in summer) were performed. In each trial, milk was divided continuously into two identical aliquots, one of which was kept at 9°C (MC9) and the other at 20°C (MC20). For each trial and milk temperature, vat milk (V-milk) and the resulting 21 month ripened cheese were analysed. RESULTS: Fat and dry matter and fat/casein ratio were lower in MC9 V-milk (p≤0.05) than in MC20. Total bacteria, mesophilic lactic acid and psychrotrophic and lipolytic bacteria showed significant differences (p≤0.05) between the two V-milks. Regarding cheese, fat content resulted lower and crude protein higher (p≤0.05) both in outer (OZ) and in inner zone (IZ) of the MC9 cheese wheels. Concerning total fatty acids, the MC9 OZ had a lower concentration of butyric, capric (p≤0.05) and medium chain fatty acids (p≤0.05), while the MC9 IZ had lower content of butyric (p≤0.05), caproic (p≤0.01) and short chain fatty acids (p≤0.05). The levels of short chain and medium chain free fatty acids (p≤0.05) were lower and that of long chain fatty acids (p≤0.05) was higher in MC9 OZ cheese. The principal component analysis of total and free fatty acids resulted in a clear separation among samples by seasons, whereas slight differences were observed between the two different milk temperatures. CONCLUSION: Storing milk at 9°C at the herd affects the chemical composition of Parmigiano Reggiano, with repercussion on lipolysis. However, the changes are not very relevant, and since the cheese can present a high variability among the different cheese factories, such changes should be considered within the "normal variations" of Parmigiano Reggiano chemical characteristics.

Herd and animal factors affect the variability of total and differential somatic cell count in bovine milk.

The aim of this study was to quantify some environmental (individual herds, herd productivity, milking system, and season) and animal factors [individual animals, breed, days in milk (DIM) and parity] on the variability of the log-10 transformation of somatic cell count (LSCC) and differential somatic cell count (DSCC) on individual bovine milk. A total of 159,360 test-day records related to milk production and composition were extracted from 12,849 Holstein-Friesian and 9,275 Simmental cows distributed across 223 herds. Herds were classified into high and low productivity, defined according to the average daily milk net energy output (DMEO) yielded by the cows. Data included daily milk yield (DYM; kg/d), milk fat, protein, lactose, SCC, and DSCC, and information on herds (i.e., productivity, milking system). The daily production of total and differential somatic cells in milk was calculated and then log-10 transformed, obtaining DLSCC and DLDSCC, respectively. Data were analyzed using a mixed model including the effects of individual herd, animal, repeated measurements intra animal as random, and herd productivity, milking system, season, breed, DIM, parity, DIM × parity, breed × season, DIM × milking system and parity × milking system as fixed factors. Herds with a high DMEO were characterized by a lower content of LSCC and DSCC, and higher DLSCC and DLDSCC, compared to the low DMEO herds. The association between milking system and somatic cell traits suggested that the use of the automatic milking systems would not allow for a rapid intervention on the cow, as evidenced by the higher content of all somatic cell traits compared to the other milking systems. Season was an important source of variation, as evidenced by high LSCC and DSCC content in milk during summer. Breed of cow had a large influence, with Holstein-Friesian having greater LSCC, DSCC, DLSCC, and DLDSCC compared to Simmental. With regard to DIM, the variability of LSCC was mostly related to that of DSCC, showing an increase from calving to the end of lactation, and suggesting the higher occurrence of chronic mastitis in cows toward the end of lactation. All the somatic cell traits increased across number of parities, possibly because older cows may have increased susceptibility to intramammary infections.

This study investigated factors affecting the variability of somatic cell traits in bovine milk. Animal had greater influence on somatic cell score (SCS) and differential somatic cell count (DSCC) compared to herd factors. Herds producing high average of daily milk energy were characterized by lower SCS and DSCC compared to the low average daily milk energy herds. The SCS and DSCC were higher in Holstein-Friesian than in Simmental, and during summer with respect to the other seasons. Older cows at the end of lactation showed the highest content of somatic cell traits. These results are helpful for the management of somatic cell traits at herd and animal levels.

Prediction of fresh and ripened cheese yield using detailed milk composition and udder health indicators from individual Brown Swiss cows.

The composition of raw milk is of major importance for dairy products, especially fat, protein, and casein (CN) contents, which are used worldwide in breeding programs for dairy species because of their role in human nutrition and in determining cheese yield (%CY). The aim of the study was to develop formulas based on detailed milk composition to disentangle the role of each milk component on %CY traits. To this end, 1,271 individual milk samples (1.5 L/cow) from Brown Swiss cows were processed according to a laboratory model cheese-making procedure. Fresh %CY (%CYCURD), total solids and water retained in the fresh cheese (%CYSOLIDS and %CYWATER), and 60-days ripened cheese (%CYRIPENED) were the reference traits and were used as response variables. Training-testing linear regression modeling was performed: 80% of observations were randomly assigned to the training set, 20% to the validation set, and the procedure was repeated 10 times. Four groups of predictive equations were identified, in which different combinations of predictors were tested separately to predict %CY traits: (i) basic composition, i.e., fat, protein, and CN, tested individually and in combination; (ii) udder health indicators (UHI), i.e., fat + protein or CN + lactose and/or somatic cell score (SCS); (iii) detailed protein profile, i.e., fat + protein fractions [CN fractions, whey proteins, and nonprotein nitrogen (NPN) compounds]; (iv) detailed protein profile + UHI, i.e., fat + protein fractions + NPN compounds and/or UHI. Aside from the positive effect of fat, protein, and total casein on %CY, our results allowed us to disentangle the role of each casein fraction and whey protein, confirming the central role of ß-CN and κ-CN, but also showing α-lactalbumin (α-LA) to have a favorable effect, and ß-lactoglobulin (ß-LG) a negative effect. Replacing protein or casein with individual milk protein and NPN fractions in the statistical models appreciably increased the validation accuracy of the equations. The cheese industry would benefit from an improvement, through genetic selection, of traits related to cheese yield and this study offers new insights into the quantification of the influence of milk components in composite selection indices with the aim of directly enhancing cheese production.

Genomic inbreeding coefficients using imputed genotypes: Assessing different estimators in Holstein-Friesian dairy cows.

The objective of this study was to estimate inbreeding coefficients in Holstein dairy cattle using imputed SNPs data. A data set of 95,540 Italian Holstein dairy cows from the routine genomic evaluations of the Italian National Association of Holstein, Brown, and Jersey Breeders were analyzed, with 84,445 imputed SNP. Ten widely used genomic inbreeding estimators were tested, including 4 PLINK v1.9 estimators (F, FHAT1, FHAT2, FHAT3), 3 genomic relationship matrix (GRM)-based methods [VanRaden's first method with observed allele frequencies (FGRM) or with fixed frequencies at 0.5 (FGRM05), VanRaden's third method, allelic frequency free and pedigree regressed (FGRM2)], runs of homozygosity (ROH)-based estimators in a complete (FROH) and simplified version (FROH2), and proportion of homozygous SNP (FPH). Pairwise comparisons among them were made, including the comparison with traditional pedigree-based inbreeding coefficients (FPED). Our results showed variability among the genomic inbreeding estimators. Coefficients of FGRM and FHAT3 were >1, meaning that more variability has been lost than the variability that existed in the base population. Regarding the remaining ones, FGRM05, FROH, FROH2, and FPH provided coefficients within the [0,1] space and are considered comparable to FPED. Not comparable to FPED, yet with an interpretable value, can be considered the coefficients of F, FHAT2, and FGRM2. Estimators based on ROH had the highest correlation with pedigree-based coefficients (0.59-0.66), among all estimators tested. In this study, Spearman correlations were shown to possibly provide a clearer estimation of the strength of the relationship between estimators. We hypothesize that imputation might cause extreme genomic inbreeding values that deserves further investigation.

Predictive formulas for different measures of cheese yield using milk composition from individual goat samples.

The objective of this study was to develop formulas based on milk composition of individual goat samples for predicting cheese yield (%CY) traits (fresh curd, milk solids, and water retained in the curd). The specific aims were to assess and quantify (1) the contribution of major milk components (fat, protein, and casein) and udder health indicators (lactose, somatic cell count, pH, and bacterial count) on %CY traits (fresh curd, milk solids, and water retained in the curd); (2) the cheese-making method; and (3) goat breed effects on prediction accuracy of the %CY formulas. The %CY traits were analyzed in duplicate from 600 goats, using an individual laboratory cheese-making procedure (9-MilCA method; 9 mL of milk per observation) for a total of 1,200 observations. Goats were reared in 36 herds and belonged to 6 breeds (Saanen, Murciano-Granadina, Camosciata delle Alpi, Maltese, Sarda, and Sarda Primitiva). Fresh %CY (%CYCURD), total solids (%CYSOLIDS), and water retained (%CYWATER) in the curd were used as response variables. Single and multiple linear regression models were tested via different combinations of standard milk components (fat, protein, casein) and indirect udder health indicators (UHI; lactose, somatic cell count, pH, and bacterial count). The 2 %CY observations within animal were averaged, and a cross-validation (CrV) scheme was adopted, in which 80% of observations were randomly assigned to the calibration (CAL) set and 20% to the validation (VAL) set. The procedure was repeated 10 times to account for sampling variability. Further, the model presenting the best prediction accuracy in CrV (i.e., comprehensive formula) was used in a secondary analysis to assess the accuracy of the %CY predictive formulas as part of the laboratory cheese-making procedure (within-animal validation, WAV), in which the first %CY observation within animal was assigned to CAL, and the second to the VAL set. Finally, a stratified CrV (SCrV) was adopted to assess the %CY traits prediction accuracy across goat breeds, again using the best model, in which 5 breeds were included in CAL and the remaining one in the VAL set. Fitting statistics of the formulas were assessed by coefficient of determination of validation (R2VAL) and the root mean square error of validation (RMSEVAL). In CrV, the formula with the best prediction accuracy for all %CY traits included fat, casein, and UHI (R2VAL = 0.65, 0.96, and 0.23 for %CYCURD, %CYSOLIDS, and %CYWATER, respectively). The WAV procedure showed R2VAL higher than those obtained in CrV, evidencing a low effect of the 9-MilCA method and, indirectly, its high repeatability. In the SCrV, large differences for %CYCURD and %CYWATER among breeds evidenced that the breed is a fundamental factor to consider in %CY predictive formulas. These results may be useful to monitor milk composition and quantify the influence of milk traits in the composite selection indices of specific breeds, and for the direct genetic improvement of cheese production.

Effects of the Cooling Temperature at the Farm on Milk Maturation and Cheesemaking Process in the Manufacture of Parmigiano Reggiano PDO Cheese.

Parmigiano Reggiano is a hard PDO cheese made from bovine raw milk, whose microbiological characteristics have important repercussions on cheese quality. According to the EU official production protocol, milk temperature at the farm must not drop below 18 °C. The present research aimed to study the effect of cooling milk at the farm at 9 °C on the characteristics of milk and on the cheesemaking process and losses during manufacture. Six cheesemaking trials were performed in two different dairies. In each of them, two cheesemakings were made in parallel: one with milk kept at 9 °C (TM9) and the other with milk kept at 20 °C (TM20). TM9 milk, in comparison with TM20, showed after the creaming process a significant reduction not only of total bacterial count but also of psychrotrophic and lipolytic bacteria. At the same time, TM9 milk showed a higher creaming capacity and, consequently, a lower fat content than TM20. TM9 vat milk had worst coagulation properties than TM20, which caused slightly higher loss of fat and curd fines into the whey. Nevertheless, these changes were too small to influence the efficiency of the cheesemaking process; conversely, maintaining milk at the farm at 9 °C led to a reduction of the number of spoilage bacteria.

Genome-wide scan reveals genetic divergence in Italian Holstein cows bred within PDO cheese production chains.

Dairy cattle breeds have been exposed to intense artificial selection for milk production traits over the last fifty years. In Italy, where over 80% of milk is processed into cheese, selection has also focused on cheese-making traits. Due to a deep-rooted tradition in cheese-making, currently fifty Italian cheeses are marked with the Protected Designation of Origin (PDO) label as they proved traditional land of origin and procedures for milk transformation. This study aimed to explore from a genetic point of view if the presence of such diverse productive contexts in Italy have shaped in a different manner the genome of animals originally belonging to a same breed. We analyzed high density genotype data from 1000 Italian Holstein cows born between 2014 and 2018. Those animals were either farmed in one of four Italian PDO consortia or used for drinkable milk production only. Runs of Homozygosity, Bayesian Information Criterion and Discriminant Analysis of Principal Components were used to evaluate potential signs of genetic divergence within the breed. We showed that the analyzed Italian Holstein cows have genomic inbreeding level above 5% in all subgroups, reflecting the presence of ongoing artificial selection in the breed. Our study provided a comprehensive representation of the genetic structure of the Italian Holstein breed, highlighting the presence of potential genetic subgroups due to divergent dairy farming systems. This study can be used to further investigate genetic variants underlying adaptation traits in these subgroups, which in turn might be used to design more specialized breeding programs.

The mineral profile affects the coagulation pattern and cheese-making efficiency of bovine milk.

Natural variations in milk minerals, their relationships, and their associations with the coagulation process and cheese-making traits present an opportunity for the differentiation of milk destined for high-quality natural products, such as traditional specialties or Protected Designation of Origin (PDO) cheeses. The aim of this study was to quantify the effects of the native contents of Ca, P, Na, K, and Mg on 18 traits describing traditional milk coagulation properties (MCP), curd firming over time (CFt) equation parameters, cheese yield (CY) measures, and nutrient recoveries in the curd (REC) using models that either included or omitted the simultaneous effects of milk fat and casein contents. The results showed that, by including milk fat and casein and the minerals in the statistical model, we were able to determine the specific effects of each mineral on coagulation and cheese-making efficiency. In general, about two-thirds of the apparent effects of the minerals on MCP and the CFt equation parameters are actually mediated by their association with milk composition, especially casein content, whereas only one-third of the effects are direct and independent of milk composition. In the case of cheese-making traits, the effects of the minerals were mediated only negligibly by their association with milk composition. High Ca content had a positive effect on the coagulation pattern and cheese-making traits, favoring water retention in the curd in particular. Phosphorus positively affected the cheese-making traits in that it was associated with an increase in CY in terms of curd solids, and in all the nutrient recovery traits. However, a very high P content in milk was associated with lower fat recovery in the curd. The variation in the Na content in milk only mildly affected coagulation, whereas with regard to cheese-making, protein recovery was negatively associated with high concentrations of this mineral. Potassium seemed not to be actively involved in coagulation and the cheese-making process. Magnesium content tended to slow coagulation and reduce CY measures. Further studies on the relationships of minerals with casein and protein fractions could deepen our knowledge of the role of all minerals in coagulation and the cheese-making process.

Effects of Milk Storage Temperature at the Farm on the Characteristics of Parmigiano Reggiano Cheese: Chemical Composition and Proteolysis.

Parmigiano Reggiano is a Protected Designation of Origin (PDO) cheese whose official production protocol provides that milk cannot be stored at less than 18 °C at the farm. The possibility of refrigerating milk at the farm is highly debated, since it should allow for the limiting of bacterial growth, thus improving the quality of the cheese. The present research aimed to study the influence of storing the milk at 9 °C on the chemical composition and proteolysis during the ripening of Parmigiano Reggiano cheese. The experimentation considered six cheese-making trials, in which both evening and morning milks were subdivided into two parts that were maintained at 9 and 20 °C. After Parmigiano Reggiano cheese-making, one of the twin wheels obtained was analyzed after 21 months of ripening. From each cheese, two different samples were taken, one from the inner zone, and the other from the outer zone. The results of the chemical analyses evidenced that milk storage at 9 °C significantly (p ≤ 0.05) influenced fat, crude protein, soluble nitrogen and peptone nitrogen contents. Nevertheless, the differences observed with respect to the cheese obtained with milk stored under standard condition were very small and should be considered within the "normal variations" of Parmigiano Reggiano chemical characteristics.

Genetic Diversity in the Italian Holstein Dairy Cattle Based on Pedigree and SNP Data Prior and After Genomic Selection.

Genetic diversity has become an urgent matter not only in small local breeds but also in more specialized ones. While the use of genomic data in livestock breeding programs increased genetic gain, there is increasing evidence that this benefit may be counterbalanced by the potential loss of genetic variability. Thus, in this study, we aimed to investigate the genetic diversity in the Italian Holstein dairy cattle using pedigree and genomic data from cows born between 2002 and 2020. We estimated variation in inbreeding, effective population size, and generation interval and compared those aspects prior to and after the introduction of genomic selection in the breed. The dataset contained 84,443 single-nucleotide polymorphisms (SNPs), and 74,485 cows were analyzed. Pedigree depth based on complete generation equivalent was equal to 10.67. A run of homozygosity (ROH) analysis was adopted to estimate SNP-based inbreeding (FROH). The average pedigree inbreeding was 0.07, while the average FROH was more than double, being equal to 0.17. The pattern of the effective population size based on pedigree and SNP data was similar although different in scale, with a constant decrease within the last five generations. The overall inbreeding rate (ΔF) per year was equal to +0.27% and +0.44% for Fped and FROH throughout the studied period, which corresponded to about +1.35% and +2.2% per generation, respectively. A significant increase in the ΔF was found since the introduction of genomic selection in the breed. This study in the Italian Holstein dairy cattle showed the importance of controlling the loss of genetic diversity to ensure the long-term sustainability of this breed, as well as to guarantee future market demands.

Quantification of the Effect of the Cattle Breed on Milk Cheese Yield: Comparison between Italian Brown Swiss and Italian Friesian.

Milk from different cattle breeds can present different casein and fat contents, which are reflected in different cheese yields (CY). However, CY is also related to some breed-related molecular characteristics. The aim of the present work was to quantify the effect of these characteristics by comparing a series of Parmigiano Reggiano (PR) cheese-making trials made with milks from Italian Brown (IB) and Italian Friesian (IF) cattle herds. Twelve trials were carried out in a cheese factory in one year (one trial per month), each one consisting of four vats processed in parallel: three vats contained milk from three different IF cattle herds (IF1, IF2 and IF3) and one contained milk from a single IB cattle herd. A 24-h CY prediction formula was developed with data from IF1, IF2 and IF3 trials (calibration) and successively validated by applying it to 12 PR trials made with IF milk in six different cheese factories (external validation). The predicted values of 24-h CY were no different to the actual ones in both calibration and external validation. Finally, the formula was tested on trials made with IB milk. In this case, the predicted values were lower than the actual ones. The quantity of IF milk casein necessary to give the same CY of IB milk was 0.20 g/100 g.

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.