The advent of genomic selection paved the way for an unprecedented acceleration in genetic progress. The increased ability to select superior individuals has been coupled with a drastic reduction in ...the generation interval for most dairy populations, representing both an opportunity and a challenge. Homozygosity is now rapidly accumulating in dairy populations. Currently, inbreeding depression is managed mostly by culling at the farm level and by controlling the overall accumulation of homozygosity at the population level. A better understanding of how homozygosity and recessive load are related will guarantee continued genetic improvement while curtailing the accumulation of harmful recessives and maintaining enough genetic variability to ensure the possibility of selection in the face of changing environmental conditions. In this review, we present a snapshot of the current dairy selection structure as it relates to response to selection and accumulation of homozygosity, briefly outline the main approaches currently used to manage inbreeding and overall variability, and present some approaches that can be used in the short term to control accumulation of harmful recessives while maintaining sustained selection pressure.
•Natural and artificial selection apply selective pressure to specific genomic regions.•Heterozygous-Rich Regions help to understand the relationship between genetic diversity and fitness.•More than ...ten thousand cattle belonging to four breeds were included in the study.•Heterozygosity patterns of Local and cosmopolitan cattle breeds were compared.•Genes linked to fertility, survival, immune response and growth were identified.
The analysis of livestock heterozygosity is less common compared to the study of homozygous patterns. Heterozygous-Rich Regions (HRRs) may harbor significant loci for functional traits such as immune response, survival rate, and fertility. For this reason, this study was conducted to investigate and characterize the heterozygosity patterns of four beef cattle breeds, which included two cosmopolitan breeds (Limousine and Charolaise) and two local breeds (Sarda and Sardo Bruna). Our analysis identified regions with a high degree of heterozygosity using a consecutive runs approach, the Tajima D test, nucleotide diversity estimation, and Hardy Weinberg equilibrium test. These regions exhibited recurrent heterozygosity peaks and were consistently found on specific chromosomes across all breeds, specifically autosomes 15, 16, 20, and 23. The cosmopolitan and Sardo Bruna breeds also displayed peaks on autosomes 2 and 21, respectively. Thirty-five top runs shared by more than 25% of the populations were identified. These genomic fragments encompassed 18 genes, two of which are directly linked to male fertility, while four are associated with lactation. Two other genes play roles in survival and immune response. Our study also detected a region related to growth and carcass traits in Limousine breed. Our analysis of heterozygosity-rich regions revealed particular segments of the cattle genome linked to various functional traits. It appears that balancing selection is occurring in specific regions within the four examined breeds, and unexpectedly, they are common across cosmopolitan and local breeds. The genes identified hold potential for applications in breeding programs and conservation studies to investigate the phenotypes associated with these heterozygous genotypes. In addition, Tajima D test, Nucleotide diversity, and Hardy Weinberg equilibrium test confirmed the presence of heterozygous fragments found with Heterozygous-Rich Regions analysis.
Emphasizing increased profit through increased dairy cow production has revealed a negative relationship of production with fitness and health traits. Decreased cow health can affect herd ...profitability through increased rates of involuntary culling and decreased or lost milk sales. The development of genomic selection methodologies, with accompanying substantial gains in reliability for low-heritability traits, may dramatically improve the feasibility of genetic improvement of dairy cow health. Producer-recorded health information may provide a wealth of information for improvement of dairy cow health, thus improving profitability. The principal objective of this study was to use health data collected from on-farm computer systems in the United States to estimate variance components and heritability for health traits commonly experienced by dairy cows. A single-step analysis was conducted to estimate genomic variance components and heritabilities for health events, including cystic ovaries, displaced abomasum, ketosis, lameness, mastitis, metritis, and retained placenta. A blended H matrix was constructed for a threshold model with fixed effects of parity and year-season and random effects of herd-year and sire. The single-step genomic analysis produced heritability estimates that ranged from 0.02 (standard deviation = 0.005) for lameness to 0.36 (standard deviation = 0.08) for retained placenta. Significant genetic correlations were found between lameness and cystic ovaries, displaced abomasum and ketosis, displaced abomasum and metritis, and retained placenta and metritis. Sire reliabilities increased, on average, approximately 30% with the incorporation of genomic data. From the results of these analyses, it was concluded that genetic selection for health traits using producer-recorded data are feasible in the United States, and that the inclusion of genomic data substantially improves reliabilities for these traits.
The aim of this study was to assess the usefulness of measures derived from milk leukocyte differential (MLD) in practices that improve fresh cow mastitis monitoring and decrease mastitis incidence. ...Quarter milk samples were collected from Holstein and Jersey cows on d 4 and 11 postcalving. Samples were analyzed using MLD, whereby cell counts and quarter infection diagnosis were obtained. Measures derived from MLD included cell scores (total leukocyte, neutrophil, macrophage, and lymphocyte scores), cell proportions (neutrophil, macrophage, and lymphocyte percentages), cell thresholds (total leukocyte, neutrophil, macrophage, and lymphocyte thresholds), and MLD diagnosis at different threshold settings (A, B, and C). Microbiological culturing of milk samples was used to determine infection status to compare the MLD diagnosis and serve as an indicator of infection. Measures derived from the microbiological analysis included occurrence of major pathogens, minor pathogens, and infection. Data analysis was based on a linear mixed model, which was used on all measures for the estimation of the fixed effects of breed, lactation number, day of sample collection, time of sampling, and quarter location, and the random effects of animal and week of sampling. All the fixed effects studied were significant for one or more of the analyzed measures. The results of this study showed that MLD-derived measures justify further study on their use for management practices for mastitis screening and prevention in early lactation.
Genotype by environment interaction (G × E) in dairy cattle productive traits has been shown to exist, but current genetic evaluation methods do not take this component into account. As several ...environmental descriptors (e.g., climate, farming system) are known to vary within the United States, not accounting for the G × E could lead to reranking of bulls and loss in genetic gain. Using test-day records on milk yield, somatic cell score, fat, and protein percentage from all over the United States, we computed within herd-year-season daughter yield deviations for 1,087 Holstein bulls and regressed them on genetic and environmental information to estimate variance components and to assess prediction accuracy. Genomic information was obtained from a 50k SNP marker panel. Environmental effect inputs included herd (160 levels), geographical region (7 levels), geographical location (2 variables), climate information (7 variables), and management conditions of the herds (16 total variables divided in 4 subgroups). For each set of environmental descriptors, environmental, genomic, and G × E components were sequentially fitted. Variance components estimates confirmed the presence of G × E on milk yield, with its effect being larger than main genetic effect and the environmental effect for some models. Conversely, G × E was moderate for somatic cell score and small for milk composition. Genotype by environment interaction, when included, partially eroded the genomic effect (as compared with the models where G × E was not included), suggesting that the genomic variance could at least in part be attributed to G × E not appropriately accounted for. Model predictive ability was assessed using 3 cross-validation schemes (new bulls, incomplete progeny test, and new environmental conditions), and performance was compared with a reference model including only the main genomic effect. In each scenario, at least 1 of the models including G × E was able to perform better than the reference model, although it was not possible to find the overall best-performing model that included the same set of environmental descriptors. In general, the methodology used is promising in accounting for G × E in genomic predictions, but challenges exist in identifying a unique set of covariates capable of describing the entire variety of environments.
The efficiency of producing salable products in the pork industry is largely determined by costs associated with feed and by the amount and quality of lean meat produced. The objectives of this paper ...were 1) to explore heritability and genetic correlations for growth, feed efficiency, and real-time ultrasound traits using both pedigree and marker information and 2) to assess accuracy of genomic prediction for those traits using Bayes A prediction models in a Duroc terminal sire population. Body weight at birth (BW at birth) and weaning (BW at weaning) and real-time ultrasound traits, including back fat thickness (BF), muscle depth (MD), and intramuscular fat content (IMF), were collected on the basis of farm protocol. Individual feed intake and serial BW records of 1,563 boars obtained from feed intake recording equipment (FIRE; Osborne Industries Inc., Osborne, KS) were edited to obtain growth, feed intake, and feed efficiency traits, including ADG, ADFI, feed conversion ratio (FCR), and residual feed intake (RFI). Correspondingly, 1,047 boars were genotyped using the Illumina PorcineSNP60 BeadChip. The remaining 516 boars, as an independent sample, were genotyped with a low-density GGP-Porcine BeadChip and imputed to 60K. Magnitudes of heritability from pedigree analysis were moderate for growth, feed intake, and ultrasound traits (ranging from 0.44 ± 0.11 for ADG to 0.58 ± 0.09 for BF); heritability estimates were 0.32 ± 0.09 for FCR but only 0.10 ± 0.05 for RFI. Comparatively, heritability estimates using marker information by Bayes A models were about half of those from pedigree analysis, suggesting "missing heritability." Moderate positive genetic correlations between growth and feed intake (0.32 ± 0.05) and back fat (0.22 ± 0.04), as well as negative genetic correlations between growth and feed efficiency traits (-0.21 ± 0.08, -0.05 ± 0.07), indicate selection solely on growth traits may lead to an undesirable increase in feed intake, back fat, and reduced feed efficiency. Genetic correlations among growth, feed intake, and FCR assessed by a multiple-trait Bayes A model resulted in increased genetic correlation between ADG and ADFI, a negative correlation between ADFI and FCR, and a positive correlation between ADG and FCR. Accuracies of genomic prediction for the traits investigated, ranging from 9.4% for RFI to 36.5% for BF, were reported that might provide new insight into pig breeding and future selection programs using genomic information.
Efficient use of feed resources has become a clear challenge for the U.S. pork industry as feed costs continue to be the largest variable expense. The availability of the Illumina Porcine60K BeadChip ...has greatly facilitated whole-genome association studies to identify chromosomal regions harboring genes influencing those traits. The current study aimed at identifying genomic regions associated with variation in feed efficiency and several production traits in a Duroc terminal sire population, including ADFI, ADG, feed conversion ratio, residual feed intake (RFI), real-time ultrasound back fat thickness (BF), ultrasound muscle depth, intramuscular fat content (IMF), birth weight (BW at birth), and weaning weight (BW at weaning). Single trait association analyses were performed using Bayes B models with 35,140 SNP on 18 autosomes after quality control. Significance of nonoverlapping 1-Mb length windows (n = 2,380) were tested across 3 QTL inference methods: posterior distribution of windows variances from Monte Carlo Markov Chain, naive Bayes factor, and nonparametric bootstrapping. Genes within the informative QTL regions for the traits were annotated. A region ranging from166 to 140 Mb (4-Mb length) on SSC 1, approximately 8 Mb upstream of the MC4R gene, was significantly associated with ADFI, ADG, and BF, where SOCS6 and DOK6 are proposed as the most likely candidate genes. Another region affecting BW at weaning was identified on SSC 4 (84-85 Mb), harboring genes previously found to influence both human and cattle height: PLAG1, CHCHD7, RDHE2 (or SDR16C5), MOS, RPS20, LYN, and PENK. No QTL were identified for RFI, IMF, and BW at birth. In conclusion, we have identified several genomic regions associated with traits affecting nutrient utilization that could be considered for future genomic prediction to improve feed utilization.
Utilization of feed in livestock species consists of a wide range of biological processes, and therefore, its efficiency can be expressed in various ways, including direct measurement, such as daily ...feed intake, as well as indicator measures, such as feeding behavior. Measuring feed efficiency is important to the swine industry, and its accuracy can be enhanced by using automated feeding systems, which record feed intake and associated feeding behavior of individual animals. Each automated feeder space is often shared among several pigs and therefore raises concerns about social interactions among pen mates with regard to feeding behavior. The study herein used a data set of 14,901 Duroc boars with individual records on feed intake, feeding behavior, and other off-test traits. These traits were modeled with and without the random spatial effect of Pen_Room, a concatenation of room and pen, or random social interaction among pen mates. The nonheritable spatial effect of common Pen-Room was observed for traits directly measuring feed intake and accounted for up to 13% of the total phenotypic variance in the average daily feeding rate. The social interaction effect explained larger proportions of phenotypic variation in all the traits studied, with the highest being 59% for ADFI in the group of feeding behaviors, 73% for residual feed intake (RFI; RFI4 and RFI6) in the feed efficiency traits, and 69% for intramuscular fat percentage in the off-test traits. After accounting for the social interaction effect, residual BW gain and RFI and BW gain (RIG) were found to have the heritability of 0.38 and 0.18, respectively, and had strong genetic correlations with growth and off-test traits. Feeding behavior traits were found to be moderately heritable, ranging from 0.14 (ADFI) to 0.52 (average daily occupation time), and some of them were strongly correlated with feed efficiency measures; for example, there was a genetic correlation of 0.88 between ADFI and RFI6. Our work suggested that accounting for the social common pen effect was important for estimating genetic parameters of traits recorded by the automated feeding system. Residual BW gain and RIG appeared to be two robust measures of feed efficiency. Feeding behavior measures are worth further investigation as indicators of feed efficiency.
The objectives of this study were to estimate genetic parameters for fertility of Brown Swiss cattle, considering reproductive measures in different parities as different traits, and to estimate ...relationships between production traits of first lactation and fertility of heifers and first-parity and second-parity cows. Reproductive indicators were interval from parturition to first service, interval from first service to conception, interval from parturition to conception, number of inseminations to conception, conception rate at first service, and nonreturn rate at 56 d after first service. Production traits were peak milk yield (pMY), lactation milk yield, and lactation length (LL). Data included 37,546 records on heifers, and 24,098 and 15,653 records on first- and second-parity cows, respectively. Cows were reared in 2,035 herds, calved from 1999 to 2007, and were progeny of 527 AI bulls. Gibbs sampling was implemented to obtain (co)variance components using both univariate and bivariate threshold and censored linear sire models. Estimates of heritability for reproductive traits in heifers (0.016 to 0.026) were lower than those in first-parity (0.017 to 0.142) and second-parity (0.026 to 0.115) cows. Genetic correlations for fertility in first- and second-parity cows were very high (>0.920), whereas those between heifers and lactating cows were moderate (0.348 to 0.709). The latter result indicates that fertility in heifers is a different trait than fertility in lactating cows, and hence it cannot be used as robust indicator of cow fertility. Heifer fertility was not related to production traits in first lactation (genetic correlations between −0.215 and 0.251). Peak milk yield exerted a moderate and unfavorable effect on the interval from parturition to first service (genetic correlations of 0.414 and 0.353 after first and second calving, respectively), and a low and unfavorable effect on other fertility traits (genetic correlations between −0.281 and 0.295). Infertility after first calving caused a strong elongation of the lactation, and LL was negatively correlated with fertility of cows after second calving, so that LL can itself be regarded as a measure of fertility. Lactation milk yield depends on both pMY and LL, and, as such, is a cause and consequence of (in)fertility.
Holstein (HH), Jersey (JJ), and crosses of these breeds were mated to HH or JJ bulls to form purebreds, reciprocal crosses, backcrosses, and other crosses in a rotational mating system. The herd was ...located at the Center for Environmental Farming Systems in Goldsboro, North Carolina. Data for calf birth weight (CBW), calving ease (0 for unassisted, n=1,135, and 1 for assisted, n=96), and neonatal calf mortality (0 for alive, n=1,150, and 1 for abortions recorded after mid-gestation, stillborn, and dead within 48h, n=81) of calves (n=1,231) were recorded over 9 calving seasons from 2003 through 2011. Gestation length (GL) was calculated as the number of days from last insemination to calving. Linear mixed models for CBW and GL included fixed effects of sex, parity (first vs. later parities), twin status, and 6 genetic groups: HH, JJ, reciprocal F1 crosses (HJ, JH), crosses >50% Holsteins (HX) and crosses >50% Jerseys (JX), where sire breed is listed first. The CBW model also included GL as a covariate. Logistic regression for calving ease and neonatal calf mortality included fixed effects of sex, parity, and genetic group. Genetic groups were replaced by linear regression using percentage of HH genes as coefficients on the above models and included as covariates to determine various genetic effects. Year and dam were included as random effects in all models. Female calves (27.57±0.54kg), twins (26.39±1.0kg), and calves born to first-parity cows (27.67±0.56kg) had lower CBW than respective male calves (29.53±0.53kg), single births (30.71±0.19kg), or calves born to multiparous cows (29.43±0.52kg). Differences in genetic groups were observed for CBW and GL. Increased HH percentage in the calf increased CBW (+9.3±0.57kg for HH vs. JJ calves), and increased HH percentage in the dams increased CBW (+1.71±0.53kg for calves from HH dams vs. JJ dams); JH calves weighed 1.33kg more than reciprocal HJ calves. Shorter GL was observed for twin births (272.6±1.1 d), female calves (273.9±0.6 d), and for first-parity dams (273.8±0.6 d). Direct genetic effects of HH alleles shortened GL (−3.5±0.7 d), whereas maternal HH alleles increased GL (2.7±0.6 d). Female calves had lower odds ratio (0.32, confidence interval=0.10–0.99) for neonatal calf mortality in second and later parities than did male calves. Maternal heterosis in crossbred primiparous dams was associated with reduced calf mortality.