Genetic correlations between type and test-day milk yield in small dual-purpose cattle populations: The Aosta Red Pied breed as a case study

Uncovering genetic parameters and environmental influences on fertility, milk production, and quality in autochthonous Reggiana cattle

Reggiana is a local cattle breed from northern Italy known for its rusticity and profitability, due to the production of branded Parmigiano Reggiano cheese. To ensure the persistence of such profitability in the long term, an adequate breeding program is required. To this aim, in the present study we estimate the genetic parameters of the main productive and reproductive traits, and we evaluate the effect of genotype by environment interaction (GxE) on these traits using 2 environmental covariates: (1) productivity and (2) temperature-humidity index (THI). Milk, fat, protein, and casein yield were considered as daily production traits, whereas protein, fat, casein percentage, casein index, and somatic cell score were considered as milk quality traits. Finally, reproductive traits such as the number of inseminations, days open, calving interval, and calving-to-first-insemination interval were evaluated. Reggiana cattle produce an average of 19 kg of milk per day with 3.7% fat and 3.4% protein content and have excellent fertility parameters. Compared with other breeds, they have slightly lower heritability for production and quality for production traits (e.g., 0.12 [0.09; 0.15] for milk yield), but similar heritability for fertility traits. Milk, protein, and fat daily yields are highly correlated but negatively correlated with the percentage of protein, fat, and casein, whereas fertility traits have an unfavorable genetic correlation with daily production traits. When considering productivity, a consistent amount of variability due to GxE was observed for all daily production traits, somatic cell count, and casein index. A modest amount of GxE was observed for fertility parameters, while the percentage of solid content showed almost no GxE effect. A similar situation occurred when considering the THI, but no GxE interaction was observed for reproduction traits. In conclusion, this study provides useful information for the implementation of accurate selection plans in this local breed, accounting for environmental plasticity measured through the consistent GxE interaction observed.

Genetic correlations of direct and indirect genetic components of social dominance with fitness and morphology traits in cattle

BACKGROUND

Within the same species, individuals show marked variation in their social dominance. Studies on a handful of populations have indicated heritable genetic variation for this trait, which is determined by both the genetic background of the individual (direct genetic effect) and of its opponent (indirect genetic effect). However, the evolutionary consequences of selection for this trait are largely speculative, as it is not a usual target of selection in livestock populations. Moreover, studying social dominance presents the challenge of working with a phenotype with a mean value that cannot change in the population, as for every winner of an agonistic interaction there will necessarily be a loser. Thus, to investigate what could be the evolutionary response to selection for social dominance, it is necessary to focus on traits that might be correlated with it. This study investigated the genetic correlations of social dominance, both direct and indirect, with several morphology and fitness traits. We used a dataset of agonistic contests involving cattle (Bos taurus): during these contests, pairs of cows compete in ritualized interactions to assess social dominance. The outcomes of 37,996 dominance interactions performed by 8789 cows over 20 years were combined with individual data for fertility, mammary health, milk yield and morphology and analysed using bivariate animal models including indirect genetic effects.

RESULTS

We found that winning agonistic interactions has a positive genetic correlation with more developed frontal muscle mass, lower fertility, and poorer udder health. We also discovered that the trends of changes in the estimated breeding values of social dominance, udder health and more developed muscle mass were consistent with selection for social dominance in the population.

CONCLUSIONS

We present evidence that social dominance is genetically correlated with fitness traits, as well as empirical evidence of the possible evolutionary trade-offs between these traits. We show that it is feasible to estimate genetic correlations involving dyadic social traits.

Genotype by Environment Interaction and Selection Response for Milk Yield Traits and Conformation in a Local Cattle Breed Using a Reaction Norm Approach

Local breeds are often reared in various environmental conditions (EC), suggesting that genotype by environment interaction (GxE) could influence genetic progress. This study aimed at investigating GxE and response to selection (R) in Rendena cattle under diverse EC. Traits included milk, fat, and protein yields, fat and protein percentage, and somatic cell score, three-factor scores and 24 linear type traits. The traits belonged to 11,085 cows (615 sires). Variance components were estimated in a two-step reaction norm model (RNM). A single trait animal model was run to obtain the solutions of herd-EC effect, then included in a random regression sire model. A multivariate response to selection (R) in different EC was computed for traits under selection including beef traits from a performance test. GxE accounted on average for 10% of phenotypic variance, and an average rank correlation of over 0.97 was found between bull estimated breeding values (EBVs) by either including or not including GxE, with changing top ranks. For various traits, significantly greater genetic components and R were observed in plain farms, loose housing rearing system, feeding total mixed ration, and without summer pasture. Conversely, for beef traits, a greater R was found for mountain farms, loose housing, hay-based feeding and summer pasture.

Genetic study of quantitative traits supports the use of Guzerá as dual-purpose cattle

Objective

The aim of this study was to estimate genetic parameters for 305-day cumulative milk yield and components, growth, and reproductive traits in Guzerá cattle.

Methods

The evaluated traits were 305-day first-lactation cumulative yields (kg) of milk (MY305), fat (FY305), protein (PY305), lactose (LY305), and total solids (SY305); age at first calving (AFC) in days; adjusted scrotal perimeter (cm) at the ages of 365 (SP365) and 450 (SP450) days; and adjusted body weight (kg) at the ages of 210 (W210), 365 (W365), and 450 (W450) days. The (co)variance components were estimated using the restricted maximum likelihood method for single-trait, bi-trait and tri-trait analyses. Contemporary groups and additive genetic effects were included in the general mixed model. Maternal genetic and permanent environmental effects were also included for W210.

Results

The direct heritability estimates ranged from 0.16 (W210) to 0.32 (MY305). The maternal heritability estimate for W210 was 0.03. Genetic correlation estimates among milk production traits and growth traits ranged from 0.92 to 0.99 and from 0.92 to 0.99, respectively. For milk production and growth traits, the genetic correlations ranged from 0.33 to 0.56. The genetic correlations among AFC and all other traits were negative (−0.43 to −0.27). Scrotal perimeter traits and body weights showed genetic correlations ranging from 0.41 to 0.46, and scrotal perimeter and milk production traits showed genetic correlations ranging from 0.11 to 0.30. The phenotypic correlations were similar in direction (same sign) and lower than the corresponding genetic correlations.

Conclusion

These results suggest the viability and potential of joint selection for dairy and beef traits in Guzerá cattle, taking into account reproductive traits.

Genome Wide Association Study of Beef Traits in Local Alpine Breed Reveals the Diversity of the Pathways Involved and the Role of Time Stratification

Knowledge of the genetic architecture of key growth and beef traits in livestock species has greatly improved worldwide thanks to genome-wide association studies (GWAS), which allow to link target phenotypes to Single Nucleotide Polymorphisms (SNPs) across the genome. Local dual-purpose breeds have rarely been the focus of such studies; recently, however, their value as a possible alternative to intensively farmed breeds has become clear, especially for their greater adaptability to environmental change and potential for survival in less productive areas. We performed single-step GWAS and post-GWAS analysis for body weight (BW), average daily gain (ADG), carcass fleshiness (CF) and dressing percentage (DP) in 1,690 individuals of local alpine cattle breed, Rendena. This breed is typical of alpine pastures, with a marked dual-purpose attitude and good genetic diversity. Moreover, we considered two of the target phenotypes (BW and ADG) at different times in the individuals' life, a potentially important aspect in the study of the traits' genetic architecture. We identified 8 significant and 47 suggestively associated SNPs, located in 14 autosomal chromosomes (BTA). Among the strongest signals, 3 significant and 16 suggestive SNPs were associated with ADG and were located on BTA10 (50-60 Mb), while the hotspot associated with CF and DP was on BTA18 (55-62 MB). Among the significant SNPs some were mapped within genes, such as SLC12A1, CGNL1, PRTG (ADG), LOC513941 (CF), NLRP2 (CF and DP), CDC155 (DP). Pathway analysis showed great diversity in the biological pathways linked to the different traits; several were associated with neurogenesis and synaptic transmission, but actin-related and transmembrane transport pathways were also represented. Time-stratification highlighted how the genetic architectures of the same traits were markedly different between different ages. The results from our GWAS of beef traits in Rendena led to the detection of a variety of genes both well-known and novel. We argue that our results show that expanding genomic research to local breeds can reveal hitherto undetected genetic architectures in livestock worldwide. This could greatly help efforts to map genomic complexity of the traits of interest and to make appropriate breeding decisions.

Selection Response Due to Different Combination of Antagonistic Milk, Beef, and Morphological Traits in the Alpine Grey Cattle Breed

Selection in local dual-purpose breeds requires great carefulness because of the need to preserve peculiar traits and also guarantee the positive genetic progress for milk and beef production to maintain economic competitiveness. A specific breeding plan accounting for milk, beef, and functional traits is required by breeders of the Alpine Grey cattle (AG), a local dual-purpose breed of the Italian Alps. Hereditability and genetic correlations among all traits have been analyzed for this purpose. After that, different selection indexes were proposed to identify the most suitable for this breed. Firstly, a genetic parameters analysis was carried out with different datasets. The milk dataset contained 406,918 test day records of milk, protein, and fat yields and somatic cells (expressed as SCS). The beef dataset included performance test data conducted on 749 young bulls. Average daily gain, in vivo estimated carcass yields, and carcass conformation (SEUROP) were the phenotypes obtained from the performance tests. The morphological dataset included 21 linear type evaluations of 11,320 first party cows. Linear type traits were aggregated through factor analysis and three factors were retained, while head typicality (HT) and rear muscularity (RM) were analyzed as single traits. Heritability estimates (h2) for milk traits ranged from 0.125 to 0.219. Analysis of beef traits showed h2 greater than milk traits, ranging from 0.282 to 0.501. Type traits showed a medium value of h2 ranging from 0.238 to 0.374. Regarding genetic correlation, SCS and milk traits were strongly positively correlated. Milk traits had a negative genetic correlation with the factor accounting for udder conformations (-0.40) and with all performance test traits and RM. These latter traits showed also a negative genetic correlation with udder volume (-0.28). The HT and the factor accounting for rear legs traits were not correlated with milk traits, but negatively correlated with beef traits (-0.32 with RM). We argue that the consequence of these results is that the use of the current selection index, which is mainly focused on milk attitude, will lead to a deterioration of all other traits. In this study, we propose more appropriate selection indexes that account for genetic relationships among traits, including functional traits.

Genetic relationship between litter size traits at birth and body measurement and production traits in purebred Duroc pigs

Heritabilities of litter size traits at birth (total number born (TNB), number born alive (NBA), and number still born (NSB)) and their genetic correlations with body measurement (body height, body length, front width (FW), chest width (CW), hind width, chest depth, chest girth, front cannon circumference, and rear cannon circumference) and production traits (ages at the start and end of performance testing (D30 and D105), average daily gain (ADG), backfat thickness, and loin muscle area) in purebred Duroc pigs were estimated. Records of performance testing for 2,835 animals and farrowing records of 1,168 litters from 437 dams were used. Genetic parameters were estimated using single-trait and two-trait animal models. Permanent environment effect was considered for litter size traits and common litter environmental effect was considered for body measurement and production traits. The estimated heritability was 0.10 ± 0.06 for TNB, 0.16 ± 0.06 for NBA, and 0.08 ± 0.05 for NSB. Positive genetic correlation of NBA was estimated with D30, D105, and ADG (0.51, 0.11, and 0.39). The estimated genetic correlation of NBA was 0.47 ± 0.17 with FW and 0.55 ± 0.18 with CW, implying that FW and CW could be promising indicator traits for efficiently improving NBA.

Multivariate analysis of body morphometric traits in conjunction with performance of reproduction and milk traits in crossbred progeny of Murrah × Jafarabadi buffalo (Bubalus bubalis) in North-Eastern Brazil

We studied the relationship between body morphometric traits, and their underlying association with milk production (MP), lactation length (LL), first calving interval (FCI) and subsequent calving interval (CI) of crossbred progeny of Murrah × Jafarabadi buffalo aiming to assist in selection programs. We carried out principal component analysis (PCA) of the body morphometric traits, which include breast width (BW), thigh width (THW), hip width (HW), rump width (RW), rump length (RL), body depth (BD), body length (BL), height withers (HEW), rear height (RH), shoulder width (SW), thorax width (TW), loin width (LW), distance from the head to ischium (DHI), and thoracic perimeter (TP). We determined the association of morphometric traits with milk/reproduction traits using canonical correlation analysis (CCA). The analysis revealed that the first six PCA accounted for 82.14% of the total observed variation, and the traits THW, HW, TW, LW, RW, HEW, TP, RH, and BW, accounted for almost half (48.00%) of the total variance indicating a higher contribution in body structural conformation. The first canonical function was significant (p<0.05), accounted for 72.46% of the total variance, and the canonical correlation was 0.56, indicating the dependence between both groups of traits. Higher canonical loadings were obtained for LL (0.49), FCI (0.46), BW (-0.71), BL (-0.56), DHI (-0.34), HEW, (-0.38) and TP (-0.50). These traits were most important for the derivation of canonical statistical variables, and presented a higher canonical correlations (r) between the dependent (LL/FCI) and independent (BW, BL, DHI, HEW and TP) groups. The results could suggest that the body morphometric traits THW, HW, TW, LW, RW, HEW, BD, TP, RH, and BW could play important role in body structural composition, indicating a suitable functional type, and aid designing of selection programs for buffalo breeding.

Investigation of myostatin and calpain 3 gene polymorphisms and their association with milk-production traits in Sfakia sheep

Context Genetic selection based on genetic markers for economically important traits in Sfakia sheep. Aims The aim of the present study was to investigate variation in the ovine myostatin gene (MSTN) and calpain 3 gene (CAPN3), and their association with milk-production traits. Methods Records for milk yield, milk fat content, protein content, lactose content, and non-fat solid content, pH and somatic-cell score (log), were obtained from 376 Sfakia ewes. Polymerase chain reaction–single-strand conformational polymorphism (PCR–SSCP) analyses were used to detect variation in intron 1 of MSTN and exon 10 of CAPN3. General linear models were then used to test for associations between the variation in MSTN and CAPN3, and milk-production traits. Key results The SSCP banding patterns for MSTN showed four variants (A1, A2, A3 and A4), which contained nine nucleotide sequence differences. Four SSCP banding patterns (C1, C2, C3 and C4) were observed for CAPN3 and these contained eight nucleotide-sequence differences. The MSTN variation was associated (P &lt; 0.05) with variation in milk yield and non-fat milk solid content. Variation in CAPN3 was associated with milk yield (P &lt; 0.001), fat content (P &lt; 0.05) and lactose content (P &lt; 0.05). Association analyses between the presence/absence of MSTN and CAPN3 variants and milk-production traits showed that a variant of MSTN that had previously between associated with muscle hypertrophy was associated with decreased milk yield (P &lt; 0.05) and a lower non-fat milk solid content (P &lt; 0.01). A CAPN3 variant that had previously been associated with increased sheep-carcass loin lean-meat yield was associated with a decreased milk yield (P &lt; 0.01) and a decreased milk fat content (P &lt; 0.05). Conclusions Our results have provided an insight into the effects of variation in ovine MSTN and CAPN3 on milk-production traits in sheep. Implications To preserve the dual-purpose characteristics of Sfakia sheep, breeding goals should take into account the possible antagonism between meat and milk traits.

Genetic correlations of fighting ability with somatic cells and longevity in cattle

The success in competitions may be stressful for animals and costly in terms of immune functions and longevity. Focusing on Aosta Chestnut and Aosta Black Pied cattle, selected for their fighting ability in traditional competitions, this study investigated the genetic relationships of fighting ability with udder health traits (somatic cell score and two threshold traits for somatic cells), longevity (length of productive life and number of calvings) and test-day milk, fat and protein yield. Herdbook information and phenotypic records that have been routinely collected for breeding programs in 16 years were used for the abovementioned traits. Data belonged to 9328 cows and 19 283 animals in pedigree. Single-trait animal model analyses were run using a Gibbs sampling algorithm to estimate the variance components of traits, and bivariate analyses were then performed to estimate the genetic correlations. Moderate positive genetic correlations (ra) were found for fighting ability with somatic cell score (ra=0.255), suggesting that greater fighting ability is genetically related to a detriment in udder health, in agreement with the theory. The high positive genetic correlation between fighting ability and longevity (average ra=0.669) suggests that the economic importance of fighting ability (the winning cows get an higher price at selling) had probably masked the true genetic covariances. The genetic correlation between milk yield traits and fighting ability showed large intervals, but the negative values (average ra=-0.121) agreed with previous research. This study is one of the few empirical studies on genetic correlations for the competitive success v. immune functions and longevity traits. The knowledge of the genetic correlations among productive and functional traits of interest, including fighting ability, is important in animal breeding for a sustainable genetic improvement.

Genetic parameter estimates for body conformation traits using composite index, principal component, and factor analysis

Information about genetic parameters is population specific and it is crucial for designing animal breeding programs and predicting response to selection. This study was carried out to estimate the genetic parameters for 23 body conformation traits of 45,517 Chinese Holstein reared in Eastern China from 1995 to 2017 with the Bayesian inference method using a linear animal mixed model. The methods to integrate these traits included (1) using the composite index from the Dairy Association of China and (2) applying principal component analysis and factor analysis to explore the relationship between the conformation traits. Estimates of heritability using the composite index were low (0.04; feet and legs) to moderate (0.23; body capacity). Strong genetic correlations were observed between the individual body conformation traits. Both principal components (1 to 7; eigenvalues ≥ 1) and latent factors (1 to 7; eigenvalues ≥ 1) explained 60.37% of total variability. Principal component 1 and factor 1 accounted for the traits that are usually associated with milk production. Moderate to low heritability were estimated through multi-trait analysis for principal components (from 0.07 to 0.21) and latent factors (from 0.07 to 0.23). Genetic correlations among the 2 multivariate techniques are typically lower compared with the one existing among the measured traits. Results from these analyses suggest the possibility of using both principal component analysis and factor analysis in morphological evaluation, simplifying the information given by the body conformation traits into new variables that could be useful for the genetic improvement of the Chinese Holstein population. This information could also be used to avoid analyzing large number of correlated traits, thereby improving precision and reducing computation burdens to analyze large and complex data.

Genetic correlations among milk yield, morphology, performance test traits and somatic cells in dual-purpose Rendena breed

Selection in native local breeds needs great carefulness due to the small population size and the risk of inbreeding. Furthermore, most breeds are dual-purpose, and milk and beef attitudes are antagonistic. For preservation purposes functional traits need to be considered. Focusing on the small local Rendena cattle, this study aimed to analyse the genetic correlations among milk, beef and udder health traits and the response to selection predicted under different scenarios. The study considered milk, fat and protein yields (MY), factor scores for udder volume (UV), conformation (UC) and muscularity obtained from type traits scored on primiparous cows, and performance test traits (PT) measured on young bulls at test station: average daily gain, in vivo SEUROP fleshiness, in vivo dressing percentage. Somatic cell score (SCS) was considered as a functional trait, with a possibility of restricting its genetic gain to zero. The study considered 281 497 MY test-day data collected on 16 974 cows, and data from linear type evaluation on 11 992 primiparous cows for factor scores. The PT data were recorded on 1428 young bulls, and SCS obtained from cell counts at milk recording. Bi-trait restricted maximum likelihood animal model analyses were performed to assess genetic parameters. Heritability varied from 0.157 (fat) to 0.442 (dressing percentage). Udder volume and MY resulted positively genetically correlated (average correlation 0.427), whereas the low-negative genetic correlation between MY and UC (-0.141) suggested a negative impact of milk gain on udder form. Beef traits of factor muscularity and PT showed medium-high favourable genetic correlations (from 0.357 to 0.984), excluding a null correlation between daily gain and muscularity. The genetic correlation MY v. muscularity was unfavourable (-0.328 on average), whereas null correlations were found in MY v. PT, apart from fat v. dressing percentage (-0.151). Somatic cell score showed low unfavourable correlations with protein (0.111) and UV (0.092), and favourable correlations with UC (-0.193). Response to selection in different scenarios indicated a good balanced gain for milk and beef when standardized economic weights of 0.66 and 0.34 are given to the two attitudes, and SCS genetic gain is restricted. Current genetic trends (MY and PT increasing, but muscularity lessening) reflect a stronger selection for milk, suggesting a slight progressive change towards a milk conformation. Aiming to preserve the dual-purpose characteristics of a breed, proper breeding policies taking into account the genetic relationships among traits and including functional traits should be applied in local dual-purpose populations.