Lifetime milk production of Holstein cattle in the humid tropics compared to Holstein-Gyr and Holstein-Brahman crosses.

Crossbred cattle are commonly used for milk production in the tropics, combining the potential benefits of pure breeds with the heterosis effects of the offspring. However, no comprehensive assessment of lifetime productivity for crossbred versus purebred cattle in low-altitude tropical environments has been carried out. The present study compares the lifetime productivity of purebred Holstein (HO, n = 17,269), Gyr (GY4, n = 435), and Brahman (BR4, n = 622) with crossbreds Gyr × Holstein (GY × HO, n = 5521) and Brahman×Holstein (BR × HO, n = 5429) cows from dairy farms located in low and medium altitude tropical regions in Costa Rica. The production traits of interest were age at first calving (AFC), days open (DO), milk production per lactation (TMP), lactation length (LLEN), age at culling (ACUL), and number of lactations (NLAC). Estimates of heterosis were also calculated. The AFC for GY × HO crosses (33-34 months) was not significantly different (p > .05) from HO (33.8 months). For BR × HO crosses, a significant (p < .05) decrease in AFC (BR3HO1 35.6 months, BR2HO2 34.5 months, and BR1H03 33.3 months) was observed as the fraction of HO breed increased. Estimates of heterosis for AFC were favourable for both crosses, of a magnitude close to 3%. The DO for F1 crosses (GY2HO2 94 days; BR2HO2 96 days) was significantly (p < .05) lower than HO (123 days). Estimates of heterosis for DO were also favourable and above 15% for both crosses. The TMP and LLEN were higher for HO (TMP = 5003 kg; LLEN = 324 days) compared with GY × HO (TMP = 4428 to 4773 kg; LLEN = 298 to 312 days) and BR × HO (TMP = 3950 to 4761 kg; LLEN = 273 to 313 days) crosses. Heterosis for TMP was favourable but low for both crosses, with a magnitude below 3.0%. The NLAC for HO (4.6 lactations) was significantly (p < .05) lower than F1 (GY2HO2, 5.8 lactations; BR2HO2, 5.4 lactations). Heterosis for NLAC was above 6.0% for both crosses. Overall, estimates of lifetime income over feed costs per cow on average were USD 2637 (30.3%) and USD 734 (8.4%) higher in F1 GY × HO and BR × HO, respectively, compared to HO. In conclusion, crossbred animals, specifically those with Gyr and Brahman genetics, extend the productive lifespan, increasing economic returns.

Genomic regions, candidate genes, and pleiotropic variants associated with physiological and anatomical indicators of heat stress response in lactating sows.

BACKGROUND: Heat stress (HS) poses significant threats to the sustainability of livestock production. Genetically improving heat tolerance could enhance animal welfare and minimize production losses during HS events. Measuring phenotypic indicators of HS response and understanding their genetic background are crucial steps to optimize breeding schemes for improved climatic resilience. The identification of genomic regions and candidate genes influencing the traits of interest, including variants with pleiotropic effects, enables the refinement of genotyping panels used to perform genomic prediction of breeding values and contributes to unraveling the biological mechanisms influencing heat stress response. Therefore, the main objectives of this study were to identify genomic regions, candidate genes, and potential pleiotropic variants significantly associated with indicators of HS response in lactating sows using imputed whole-genome sequence (WGS) data. Phenotypic records for 18 traits and genomic information from 1,645 lactating sows were available for the study. The genotypes from the PorcineSNP50K panel containing 50,703 single nucleotide polymorphisms (SNPs) were imputed to WGS and after quality control, 1,622 animals and 7,065,922 SNPs were included in the analyses. RESULTS: A total of 1,388 unique SNPs located on sixteen chromosomes were found to be associated with 11 traits. Twenty gene ontology terms and 11 biological pathways were shown to be associated with variability in ear skin temperature, shoulder skin temperature, rump skin temperature, tail skin temperature, respiration rate, panting score, vaginal temperature automatically measured every 10 min, vaginal temperature measured at 0800 h, hair density score, body condition score, and ear area. Seven, five, six, two, seven, 15, and 14 genes with potential pleiotropic effects were identified for indicators of skin temperature, vaginal temperature, animal temperature, respiration rate, thermoregulatory traits, anatomical traits, and all traits, respectively. CONCLUSIONS: Physiological and anatomical indicators of HS response in lactating sows are heritable but highly polygenic. The candidate genes found are associated with important gene ontology terms and biological pathways related to heat shock protein activities, immune response, and cellular oxidative stress. Many of the candidate genes with pleiotropic effects are involved in catalytic activities to reduce cell damage from oxidative stress and cellular mechanisms related to immune response.

Genetic evaluations of dairy goats with few pedigree data: different approaches to use molecular information.

One of the limitations of implementing animal breeding programs in small-scale or extensive production systems is the lack of production records and genealogical records. In this context, molecular markers could help to gain information for the breeding program. This study addresses the inclusion of molecular data into traditional genetic evaluation models as a random effect by molecular pedigree reconstruction and as a fixed effect by Bayesian clustering. The methods were tested for lactation curve traits in 14 dairy goat herds with incomplete phenotypic data and pedigree information. The results showed an increment of 37.3% of the relationships regarding the originals with MOLCOAN and clustering into five genetic groups. Data leads to estimating additive variance, error variance, and heritability with four different models, including pedigree and molecular information. Deviance Information Criterion (DIC) values demonstrate a greater fitting of the models that include molecular information either as fixed (genetic clusters) or as random (molecular matrix) effects. The molecular information of simple markers can complement genetic improvement strategies in populations with little information.

Genetic association between fat-to-protein ratio and traits of economic interest in early lactation Holstein cows in Brazil.

The aims of this study were to estimate the genetic parameters for fat-to-protein ratio (F:P) within the first 90 days of lactation and to examine their genetic associations with daily milk yield (MY), somatic cell score (SCS), and calving interval between the first and second calving (IFSC) and between the second and third calving (ISTC) during the first three lactations of Holstein cows. We utilized 200,626 production-related data officially recorded from 77,436 cows milked two or three times a day from 2012 to 2022, sourced from the Holstein Cattle Breeders Association of Paraná State, Brazil. The (co)variance components were estimated using animal models, adopting the restricted maximum likelihood (REML) method with single-trait analysis (for heritability and repeatability) and two-trait analysis (for genetic and phenotypic correlations), per lactation. Regardless of lactation number, heritability estimates were relatively low, ranging from 0.08 ± 0.005 to 0.10 ± 0.003 for F:P; 0.08 ± 0.01 to 0.18 ± 0.005 for MY; 0.04 ± 0.01 to 0.07 ± 0.004 for SCS; and 0.03 ± 0.01 for both IFSC and ISTC. Repeatability estimates within the same lactation were low for F:P (ranging from 0.17 ± 0.002 to 0.19 ± 0.03), high for MY (between 0.50 ± 0.003 and 0.53 ± 0.002), and moderate to high for SCS (between 0.39 ± 0.003 and 0.44 ± 0.004). Genetic correlations between F:P and MY ranged from -0.26 ± 0.03 to -0.15 ± 0.02; F:P and SCS, from -0.06 ± 0.03 to -0.03 ± 0.08; F:P and IFSC, 0.31 ± 0.01; F:P and ISTC, 0.20 ± 0.01; MY and IFSC, 0.24 ± 0.05; and MY and ISTC, 0.13 ± 0.08. The fat-to-protein ratio during early lactation showed low genetic variability, regardless of lactation number. Furthermore, it was genetically correlated with MY, IFSC, and ISTC, although there is an antagonistic and unfavorable correlation between traits that can limit genetic progress.

Genetic parameters and genome-wide association studies for mozzarella and milk production traits, lactation length, and lactation persistency in Murrah buffaloes.

Genetic and genomic analyses of longitudinal traits related to milk production efficiency are paramount for optimizing water buffaloes breeding schemes. Therefore, this study aimed to (1) compare single-trait random regression models under a single-step genomic BLUP setting based on alternative covariance functions (i.e., Wood, Wilmink, and Ali and Schaeffer) to describe milk (MY), fat (FY), protein (PY), and mozzarella (MZY) yields, fat-to-protein ratio (FPR), somatic cell score (SCS), lactation length (LL), and lactation persistency (LP) in Murrah dairy buffaloes (Bubalus bubalis); (2) combine the best functions for each trait under a multiple-trait framework; (3) estimate time-dependent SNP effects for all the studied longitudinal traits; and (4) identify the most likely candidate genes associated with the traits. A total of 323,140 test-day records from the first lactation of 4,588 Murrah buffaloes were made available for the study. The model included the average curve of the population nested within herd-year-season of calving, systematic effects of number of milkings per day, and age at first calving as linear and quadratic covariates, and additive genetic, permanent environment, and residual as random effects. The Wood model had the best goodness of fit based on the deviance information criterion and posterior model probabilities for all traits. Moderate heritabilities were estimated over time for most traits (0.30 ± 0.02 for MY; 0.26 ± 0.03 for FY; 0.45 ± 0.04 for PY; 0.28 ± 0.05 for MZY; 0.13 ± 0.02 for FPR; and 0.15 ± 0.03 for SCS). The heritability estimates for LP ranged from 0.38 ± 0.02 to 0.65 ± 0.03 depending on the trait definition used. Similarly, heritabilities estimated for LL ranged from 0.10 ± 0.01 to 0.14 ± 0.03. The genetic correlation estimates across days in milk (DIM) for all traits ranged from -0.06 (186-215 DIM for MY-SCS) to 0.78 (66-95 DIM for PY-MZY). The SNP effects calculated for the random regression model coefficients were used to estimate the SNP effects throughout the lactation curve (from 5 to 305 d). Numerous relevant genomic regions and candidate genes were identified for all traits, confirming their polygenic nature. The candidate genes identified contribute to a better understanding of the genetic background of milk-related traits in Murrah buffaloes and reinforce the value of incorporating genomic information in their breeding programs.

Genomic-based genetic parameters for milkability traits derived from automatic milking systems in North American Holstein cattle.

The number of dairy farms adopting automatic milking systems (AMS) has considerably increased around the world aiming to reduce labor costs, improve cow welfare, increase overall performance, and generate a large amount of daily data, including production, behavior, health, and milk quality records. In this context, this study aimed to (1) estimate genomic-based variance components for milkability traits derived from AMS in North American Holstein cattle based on random regression models; and (2) derive and estimate genetic parameters for novel behavioral indicators based on AMS-derived data. A total of 1,752,713 daily records collected using 36 milking robot stations and 70,958 test-day records from 4,118 genotyped Holstein cows were used in this study. A total of 57,600 SNP remained after quality control. The daily-measured traits evaluated were milk yield (MY, kg), somatic cell score (SCS, score unit), milk electrical conductivity (EC, mS), milking efficiency (ME, kg/min), average milk flow rate (FR, kg/min), maximum milk flow rate (FRM, kg/min), milking time (MT, min), milking failures (MFAIL), and milking refusals (MREF). Variance components and genetic parameters for MY, SCS, ME, FR, FRM, MT, and EC were estimated using the AIREMLF90 software under a random regression model fitting a third-order Legendre orthogonal polynomial. A threshold Bayesian model using the THRGIBBS1F90 software was used for genetically evaluating MFAIL and MREF. The daily heritability estimates across days in milk (DIM) ranged from 0.07 to 0.28 for MY, 0.02 to 0.08 for SCS, 0.38 to 0.49 for EC, 0.45 to 0.56 for ME, 0.43 to 0.52 for FR, 0.47 to 0.58 for FRM, and 0.22 to 0.28 for MT. The estimates of heritability (± SD) for MFAIL and MREF were 0.02 ± 0.01 and 0.09 ± 0.01, respectively. Slight differences in the genetic correlations were observed across DIM for each trait. Strong and positive genetic correlations were observed among ME, FR, and FRM, with estimates ranging from 0.94 to 0.99. Also, moderate to high and negative genetic correlations (ranging from -0.48 to -0.86) were observed between MT and other traits such as SCS, ME, FR, and FRM. The genetic correlation (± SD) between MFAIL and MREF was 0.25 ± 0.02, indicating that both traits are influenced by different sets of genes. High and negative genetic correlations were observed between MFAIL and FR (-0.58 ± 0.02) and MFAIL and FRM (-0.56 ± 0.02), indicating that cows with more MFAIL are those with lower FR. The use of random regression models is a useful alternative for genetically evaluating AMS-derived traits measured throughout the lactation. All the milkability traits evaluated in this study are heritable and have demonstrated selective potential, suggesting that their use in dairy cattle breeding programs can improve dairy production efficiency in AMS.

An insight into population structure and genetic progress of Argentinean Holstein cattle.

The objective of the present study was to investigate the most critical issues associated with the limited genetic progress evidenced in the Argentinean Holstein ("Holando Argentino") breed in the last 20 years (only 26% of the phenotypic trend in milk yield was due to genetics). The study comprised the analysis of population structure, realized genetic selection differentials, genetic progress and partition of genetic trends by sex and country of origin from 1936 to 2019 (1,045,582 records; 24,680 sires and 619,322 dams in the pedigree). Average inbreeding steadily increased in the last 15 generations (ΔF = 0.6%, which translates to Ne = 75). Partition of genetic trends revealed that local genetics made a negligible contribution to genetic progress, which for most traits was highly dependent on imported genetics (>80%). Mean generation intervals were fairly constant until 2009 (8-9 years for males and 5-6 years for females, respectively) and then decreased, especially in the paths of sires of bulls and dams of bulls (to 5 and 4 years, respectively) mostly due to the influence of imported sires. The reduction in generation intervals was counterbalanced by a marked deterioration of realized selection differentials, particularly in the path of sires of bulls that nevertheless made the largest contribution to genetic progress. In the last 20 years, realized selection differentials in this path went from 533.6 to 170.8 kg for milk yield and from 16.7 to 13.3 kg for protein yield (1.7-0.5 and 1.6-1.3 standard deviation units, respectively). Among all considered traits (milk yield, fat yield, protein yield, stature, final score and daughter pregnancy rate) in the analysed period, annual genetic gain was negative for milk yield, fairly constant for composition and conformation traits, and positive only in the case of daughter pregnancy rate. Considered together, these results suggest that limited genetic progress is due to the absence of a sound breeding programme that includes genomic selection and a carefully defined selection objective, together with the absence of stronger regulations in germplasm importation; however, other factors such as potential genetics by environment interactions cannot be ruled out.

Genotype-environment interaction for milk production of Gyr cattle in Brazil and Colombia.

BACKGROUND: Genotype by environment interactions (G × E) can play an important role in cattle populations and should be included in breeding programs in order to select the best animals for different environments. OBJECTIVE: The aim of this study was to investigate the G × E for milk production of Gyr cattle in Brazil and Colombia by applying a reaction norm model used genomics information, and to identify genomic regions associated with milk production in the two countries. METHODS: The Brazilian and Colombian database included 464 animals (273 cows and 33 sires from Brazil and 158 cows from Colombia) and 27,505 SNPs. A two-trait animal model was used for milk yield adjusted to 305 days in Brazil and Colombia as a function of country of origin, which included genomic information obtained with a single-step genomic reaction norm model. The GIBBS3F90 and POSTGSf90 programs were used. RESULTS: The results obtained indicate G × E based on the reranking of bulls between Brazil and Colombia, demonstrating environmental differences between the two countries. The findings highlight the importance of considering the environment when choosing breeding animals in order to ensure the adequate performance of their progeny. Within this context, the reranking of bulls and the different SNPs associated with milk production in the two countries suggest that G × E is an important effect that should be included in the genetic evaluation of Dairy Gyr cattle in Brazil and Colombia. CONCLUSION: The Gyr breeding program can be optimized by choosing a selection environment that will allow maximum genetic progress in milk production in different environments within and between countries.

Heat stress on breeding value prediction for milk yield and composition of a Brazilian Holstein cattle population.

Due to the high milk production of Holstein cows, many countries have chosen to import semen to improve local dairy herds. This strategy would be more effective if this semen was used in the same environment conditions in which the bulls were selected. If the effect of genotype by environment (G × E) interaction is not considered, the estimated breeding values (EBVs) may vary, potentially reducing the selection response. We evaluate the impact of heat stress on selection for milk yield and composition of Holstein cows using random regression models. To verify the interference of heat stress in milk yield (MY) and composition traits (fat, protein, total saturated, and total unsaturated fatty acids content in milk), temperature-humidity index (THI) on test-day milk records was used. The threshold value to divide the environments using test-day information from Brazilian Holstein cows was 72 units of THI, i.e., < 72 represented no heat stress and > 72 represented heat stress. Legendre polynomials of second-order (Leg 2) model and two lactation points (33 and 122 DIM) were used to estimate heritabilities and EBVs for five important dairy traits. The heritabilities of milk components and fatty acids were low (0.09-0.29), regardless of lactation period and degree of heat stress, with the exception of protein content (0.30-0.35). Fat content was the only milk component that was reduced according to the degree of heat stress and lactation period. The EBVs tended to decrease in heat stress conditions, thus animals with high genetic potential demonstrated evidence of G × E interaction. However, acclimatization of dairy cows to heat stress in the farm production systems may have been responsible for the low differences among genetic parameters and EBVs with and without heat stress found in this study.

Random regression models using B-splines functions provide more accurate genomic breeding values for milk yield and lactation persistence in Murrah buffaloes.

There is a great worldwide demand for cheese made with buffalo milk, due to its flavour and nutritional properties. In this context, there is a need for increasing the efficiency of buffalo milk production (including lactation persistence), which can be achieved through genomic selection. The most used methods for the genetic evaluation of longitudinal data, such as milk-related traits, are based on random regression models (RRM). The choice of the best covariance functions and polynomial order for modelling the random effects is an important step to properly fit RRM. To our best knowledge, there are no studies evaluating the impact of the order and covariance function (Legendre polynomials-LEG and B-splines-BSP) used to fit RRM for genomic prediction of breeding values in dairy buffaloes. Therefore, the main objectives of this study were to estimate variance components and evaluate the performance of LEG and BSP functions of different orders on the predictive ability of genomic breeding values for the first three lactations of milk yield (MY1, MY2, and MY3) and lactation persistence (LP1, LP2, and LP3) of Brazilian Murrah. Twenty-two models for each lactation were contrasted based on goodness of fit, genetic parameter estimates, and predictive ability. Overall, the models of higher orders of LEG or BSP had a better performance based on the deviance information criterion (DIC). The daily heritability estimates ranged from 0.01 to 0.30 for MY1, 0.08 to 0.42 for MY2, and from 0.05 to 0.47 for MY3. For lactation persistence (LP), the heritability estimates ranged from 0.09 to 0.32 for LP1, from 0.15 to 0.33 for LP2, and from 0.06 to 0.32 for LP3. In general, the curves plotted for variance components and heritability estimates based on BSP models presented lower oscillation along the lactation trajectory. Similar predictive ability was observed among the models. Considering a balance between the complexity of the model, goodness of fit, and credibility of the results, RRM using quadratic B-splines functions based on four or five segments to model the systematic, additive genetic, and permanent environment curves provide better fit with no significant differences between genetic variances estimates, heritabilities, and predictive ability for the genomic evaluation of dairy buffaloes.

Impact of different numbers of milk test-day records during lactation on the reliability of estimated breeding values.

The objective of this study was to evaluate the impact of the reduction in the number of test-day records per lactation on genetic parameters of test-day milk yield (TDMY) and the reliability of estimated breeding values for 305-day milk yield in Holstein cattle. Estimates of genetic parameters and breeding values were performed using the animal model of random regression and adjustment of the Legendre polynomial (fourth order). When comparing sires with the same number of daughters, greater reliability was found in the subpopulations with the highest number of milk test-day records per lactation per cow. It was also found that the elimination of at least one test-day record affected the reliability of estimated breeding value for 305-day milk yield in the sires, regardless of the class of number of daughters per sire. When selecting the 5% best sires and 20% best cows, the lowest order correlations were observed between the population with 10 test-day records per lactation (complete lactation) and the other subpopulations with incomplete lactations (4, 5, 6, 7, 8, 9 test-day records per lactation). The reduction in the number of test-day milk records per lactation interferes in the reliability of the estimated breeding value for sires and cows, negatively impacting the precision of selecting genetically superior animals.

Random regression test-day models to describe milk production and fatty acid traits in first lactation Walloon Holstein cows.

We investigated the use of different Legendre polynomial orders to estimate genetic parameters for milk production and fatty acid (FA) traits in the first lactation Walloon Holstein cows. The data set comprised 302,684 test-day records of milk yield, fat and protein contents, and FAs generated by mid-infrared (MIR) spectroscopy, C16:0 (palmitic acid), C18:1 cis-9 (oleic acid), LCFAs (long-chain FAs), SFAs (saturated FAs) and UFAs (unsaturated FAs) were studied. The models included random regression coefficients for herd-year of calving (h), additive genetic (a) and permanent environment (p) effects. The selection of the best random regression model (RRM) was based on the deviance information criterion (DIC), and genetic parameters were estimated via a Bayesian approach. For all analysed random effects, DIC values decreased as the order of the Legendre polynomials increased. Best-fit models had fifth-order (degree 4) for the p effect and ranged from second- to fifth-order (degree 1-4) for the a and h effects (LEGhap: LEG555 for milk yield and protein content; LEG335 for fat content and SFA; LEG545 for C16:0 and UFA; and LEG535 for C18:1 cis-9 and LCFA). Based on the best-fit models, an effect of overcorrection was observed in early lactation (5-35 days in milk [DIM]). On the contrary, third-order (LEG333; degree 2) models showed flat residual trajectories throughout lactation. In general, the estimates of genetic variance tended to increase over DIM, for all traits. Heritabilities for milk production traits ranged from 0.11 to 0.58. Milk FA heritabilities ranged from low-to-high magnitude (0.03-0.56). High Spearman correlations (>0.90 for all bulls and >0.97 for top 100) were found among breeding values for 155 and 305 DIM between the best RRM and LEG333 model. Therefore, third-order Legendre polynomials seem to be most parsimonious and sufficient to describe milk production and FA traits in Walloon Holstein cows.

Genetic effects of heat stress on milk fatty acids in Brazilian Holstein cattle.

The present study aimed to estimate covariance components of milk fatty acids (FA) and to compare the genomic estimated breeding values under general and heat-stress effects. Data consisted of 38,762 test-day records from 6,344 Holstein cows obtained from May 2012 through January 2018 on 4 dairy herds from Brazil. Single-trait repeatability test-day models with random regressions as a function of temperature-humidity index values were used for genetic analyses. The models included contemporary groups, parity order (1-6), and days in milk classes as fixed effects, and general and thermotolerance additive genetic and permanent environmental as random effects. Notably, differences in heritability estimates between environments (general and heat stress) increased (0.03 to 0.06) for unsaturated FA traits, such as unsaturated, monounsaturated, and polyunsaturated, at higher heat-stress levels. In contrast, heritability estimated between environments for saturated FA traits, including saturated FA, palmitic acid (C16:0), and stearic acid (C18:0) did not observe significant differences between environments. In addition, our study revealed negative genetic correlations between general and heat-stress additive genetic effects (antagonistic effect) for the saturated FA, C16:0, C18:0, and C18:1, which ranged from -0.007 to -0.32. Spearman's ranking correlation between genomic estimated breeding values ranged from -0.27 to 0.99. Results indicated a moderate to strong interaction of genotype by the environment for most FA traits comparing a heat-stress environment with thermoneutral conditions. Our findings point out novel opportunities to explore the use of FA milk profile and heat-stress models.

Genomic study of the resilience of buffalo cows to a negative energy balance.

The research article was carried out with the objective of studying the genetic variation on the resilience of buffaloes to negative energy balance-NEB (measured by changes in body weight in early lactation)-as well as investigating genomic regions of interest for this trait. A model of reaction norms was used, considering milk production as the trait to be analyzed and solutions of the contemporary groups to weight changes as environmental gradient. In this methodology, the genetic value of the slope represents the measure of resilience of the animals. After the estimation step, a genome-wide association analysis was performed for the slope of the reaction norms model, to obtain a list of windows and associated genes. The heritability estimates for milk production over the resilience gradient ranged from 0.13 to 0.28, with lower values in the intermediate environmental groups. Regarding the productive resilience of dairy buffalo cows to NEB, the genomic windows with the highest contribution to the genetic variance were detected on chromosomes BBU 1, 2, 3, 4, 9, 12, 19, and 21. A functional analysis of the genes described in the selected windows indicated association with metabolic routes related to growth and immunity of the animals, with an emphasis on the STAT6 gene. The results presented indicate that there is for this trait genetic variation to be used as selection criteria, in addition to genomic regions that can increase the precision of the selection.

Rumination time, activity index, and productive performance of Holstein and crossbred Holstein × jersey cows exposed to different temperature-humidity indexes.

The aim of this study was to evaluate the rumination time, activity index, milk yield (MY), physicochemical milk characteristics, and physiological indicators of pure Holstein (H) and first (F1- ½ Holstein × Jersey) and second-generation (R1- ¾ Holstein × » Jersey) of crossbred cows in a pasture-based system under the influence of different temperature-humidity indexes (THI). Twenty-two multiparous cows (H = 7, F1 = 5, and R1 = 10) were evaluated throughout 1-year period. Daily information on rumination time and activity index was obtained using the software HealthyCow24® and the daily THI was determined from data logger information. Weekly, MY, physicochemical milk characteristics, body weight, body condition score, and physiological indicators were evaluated. THI were grouped into six classes from safe (< 68) to emergency (≥ 84). Variance analyses were performed. There was no interaction between the genetic group and THI classes for any trait. The rumination time was longer for crossbred R1 and F1 cows than for Holstein cows. The activity index increased, and rumination time decreased as the THI classes increased. There was no difference in MY between the genetic groups, but crossbred cows showed higher milk fat and protein content. Our study indicated that the high THI increases the activity index and decreases the rumination time, MY, and the physicochemical quality of milk, BW, BCS, RF, and RT of the three genetic groups (H, F1, and R1).

Comparison of Marker Effects and Breeding Values at Two Levels at THI for Milk Yield and Quality Traits in Brazilian Holstein Cows.

Genomic tools can help in the selection of animals genetically resistant to heat stress, especially the genome-wide association studies (GWAS). The objective of this study was to compare the variance explained by SNPs and direct genomic breeding values (DGVs) at two levels of a temperature and humidity index (THI). Records of milk yield (MY), somatic cell score (SCS), and percentages of casein (CAS), saturated fatty acids (SFA), and unsaturated fatty acids (UFA) in milk from 1157 Holstein cows were used. Traditional breeding values (EBV) were determined in a previous study and used as pseudo-phenotypes. Two levels of THI (heat comfort zone and heat stress zone) were used as environments and were treated as "traits" in a bi-trait model. The GWAS was performed using the genomic best linear unbiased prediction (GBLUP) method. Considering the top 50 SNPs, a total of 36 SNPs were not common between environments, eight of which were located in gene regions related to the evaluated traits. Even for those SNPs that had differences in their explained variances between the two environments, the differences were very small. The animals showed virtually no rank order, with rank correlation values of 0.90, 0.88, 1.00, 0.88, and 0.97 for MY, CAS, SCS, SFA, and UFA, respectively. The small difference between the environments studied can be attributed to the small difference in the pseudo-phenotypes used between the environments, on-farm acclimation, the polygenic nature of the traits, and the THI values studied near the threshold between comfort and heat stress. It is recommended that future studies be conducted with a larger number of animals and at more extreme THI levels.

Detection of genomic regions underlying milk production traits in Valle del Belice dairy sheep using regional heritability mapping.

The aim of this study was to identify genomic regions underlying milk production traits in the Valle del Belice dairy sheep using regional heritability mapping (RHM). Repeated measurements for milk yield (MY), fat percentage and yield (F% and FY) and protein percentage and yield (P% and PY), collected over a period of 6 years (2006-2012) on 481 Valle del Belice ewes, were used for the analysis. Animals were genotyped with the Illumina 50k SNP chip. Variance components, heritabilities and repeatabilities within and across lactations were estimated, fitting parity, litter size, season of lambing and fortnights in milk, as fixed; and additive genetic, permanent environment within and across lactations, flock by test-day interaction and residual as random effects. For the RHM analysis, the model included the same fixed and random effects as before, plus an additional regional genomic additive effect (specific for the region being tested) as random. While the whole genomic additive effect was estimated using the genomic relationship matrix (GRM) constructed from all SNPs, the regional genomic additive effect was estimated from a GRM matrix constructed from the SNPs within each region. Heritability estimates ranged between 0.06 and 0.15, with repeatabilities being between 0.14 and 0.24 across lactations and between 0.23 and 0.39 within lactation for all milk production traits. A substantial effect of flock-test-day on milk production traits was also estimated. Significant genomic regions at either genome-wide (p < .05) or suggestive (i.e., one false positive per genome scan) level were identified on chromosome (OAR) 2, 3 and 20 for F% and on OAR3 for P%, with the regions on OAR3 in common between the two traits. Our results confirmed the role of LALBA and AQP genes, on OAR3, as candidate genes for milk production traits in sheep.

Genetic principal components for reproductive and productive traits in Holstein cows reared under tropical conditions.

The objective of this study was to compare the standard multi-trait model and five reduced-rank models fitted to the first principal components and genetic parameter estimates in order to determine the most appropriate method to model the covariance structure of reproductive and productive traits in Brazilian Holstein cows. Individual records of the following traits from 5217 cows were analyzed: 305-day milk yield (MY305), peak yield, milk yield per day of calving interval, days from calving to first estrus, days from calving to last service (CLS), calving interval (CI), and gestation length. Schwarz's Bayesian information criterion was used to compare the different models. The results indicated that four principal components were necessary to model the genetic (co)variance structure, reducing the number of parameters to be estimated. Analysis of genetic and phenotypic correlations showed that milk production-related traits were strongly correlated with each other (ranging from 0.74 to 0.99), while the correlation of these traits with the reproductive traits was weak (ranging from - 0.14 to 0.27). Heritability estimates for the traits ranged from 0.03 to 0.18. The reproductive traits CLS and CI and the production trait MY305 should be included as selection criteria in dairy cattle breeding programs because they are correlated with the first two principal components, retaining 91% of the genetic variability of the data.

Genomic studies of milk-related traits in water buffalo (Bubalus bubalis) based on single-step genomic best linear unbiased prediction and random regression models.

Genomic selection has been widely implemented in many livestock breeding programs, but it remains incipient in buffalo. Therefore, this study aimed to (1) estimate variance components incorporating genomic information in Murrah buffalo; (2) evaluate the performance of genomic prediction for milk-related traits using single- and multitrait random regression models (RRM) and the single-step genomic best linear unbiased prediction approach; and (3) estimate longitudinal SNP effects and candidate genes potentially associated with time-dependent variation in milk, fat, and protein yields, as well as somatic cell score (SCS) in multiple parities. The data used to estimate the genetic parameters consisted of a total of 323,140 test-day records. The average daily heritability estimates were moderate (0.35 ± 0.02 for milk yield, 0.22 ± 0.03 for fat yield, 0.42 ± 0.03 for protein yield, and 0.16 ± 0.03 for SCS). The highest heritability estimates, considering all traits studied, were observed between 20 and 280 d in milk (DIM). The genetic correlation estimates at different DIM among the evaluated traits ranged from -0.10 (156 to 185 DIM for SCS) to 0.61 (36 to 65 DIM for fat yield). In general, direct selection for any of the traits evaluated is expected to result in indirect genetic gains for milk yield, fat yield, and protein yield but also increase SCS at certain lactation stages, which is undesirable. The predicted RRM coefficients were used to derive the genomic estimated breeding values (GEBV) for each time point (from 5 to 305 DIM). In general, the tuning parameters evaluated when constructing the hybrid genomic relationship matrices had a small effect on the GEBV accuracy and a greater effect on the bias estimates. The SNP solutions were back-solved from the GEBV predicted from the Legendre random regression coefficients, which were then used to estimate the longitudinal SNP effects (from 5 to 305 DIM). The daily SNP effect for 3 different lactation stages were performed considering 3 different lactation stages for each trait and parity: from 5 to 70, from 71 to 150, and from 151 to 305 DIM. Important genomic regions related to the analyzed traits and parities that explain more than 0.50% of the total additive genetic variance were selected for further analyses of candidate genes. In general, similar potential candidate genes were found between traits, but our results suggest evidence of differential sets of candidate genes underlying the phenotypic expression of the traits across parities. These results contribute to a better understanding of the genetic architecture of milk production traits in dairy buffalo and reinforce the relevance of incorporating genomic information to genetically evaluate longitudinal traits in dairy buffalo. Furthermore, the candidate genes identified can be used as target genes in future functional genomics studies.

Micronucleated erythrocytes in peripheral blood from neonate rats fed by nursing mothers exposed to X-rays.

Most women with breast cancer can become pregnant and give birth while undergoing radiation therapy and breastfeeding is generally not contraindicated. The induction of long-lived reactive species in proteins, such as casein by X-ray radiation and DNA damage to unexposed organisms, has been shown when ingesting irradiated cheese. To determine whether exposing lactating rats to X-rays increases the number of micronucleated erythrocytes (MNEs) in peripheral blood of their unexposed or breastfeeding rat pups, 15 female Wistar rats were divided into three groups: Negative control; Experimental group exposed to X-rays, and group exposed to X-rays plus vitamin C. The mothers of groups 2 and 3 were irradiated for three consecutive days after giving birth, returning them to their respective cages each time to continue lactation. A blood sample was taken from the mothers and pups at 0, 24, and 48 hr. Blood smears were stained with acridine orange to analyze MNEs. In mother rats, the frequency of micronucleated polychromatic erythrocytes (MNPCEs) increased significantly at 24 and 48 hr in both study groups exposed to radiation. Likewise, in rat pups the MNPCE and MNE frequencies increased in both groups with radiation and radiation plus vitamin C at 24 and 48 hr, and a protection from vitamin C was observed. In conclusion, the genotoxic damage produced in rat pups that were lactated by mothers irradiated with X-rays is possibly due to the effect of long-lived reactive species that were formed in the breast milk of female Wistar rats during the irradiation process.