Computing procedures for genetic evaluation including phenotypic, full pedigree, and genomic information

ExAutoGP: Enhancing Genomic Prediction Stability and Interpretability with Automated Machine Learning and SHAP

Machine learning has attracted much attention in the field of genomic prediction due to its powerful predictive capabilities, yet the lack of an explanatory nature in modeling decisions remains a major challenge. In this study, we propose a novel machine learning method, ExAutoGP, which aims to improve the accuracy of genomic prediction and enhance the transparency of the model by combining automated machine learning (AutoML) with SHapley Additive exPlanations (SHAP). To evaluate ExAutoGP's effectiveness, we designed a comparative experiment consisting of a simulated dataset and two real animal datasets. For each dataset, we applied ExAutoGP and five baseline models-Genomic Best Linear Unbiased Prediction (GBLUP), BayesB, Support Vector Regression (SVR), Kernel Ridge Regression (KRR), and Random Forest (RF). All models were trained and evaluated using five repeated five-fold cross-validation, and their performance was assessed based on both predictive accuracy and computational efficiency. The results show that ExAutoGP exhibits robust and excellent prediction performance on all datasets. In addition, the SHAP method not only effectively reveals the decision-making process of ExAutoGP and enhances its interpretability, but also identifies genetic markers closely related to the traits. This study demonstrates the strong potential of AutoML in genomic prediction, while the introduction of SHAP provides actionable biological insights. The synergy of high prediction accuracy and interpretability offers new perspectives for optimizing genomic selection strategies in livestock and poultry breeding.

Deep learning and genomic best linear unbiased prediction integration: An approach to identify potential nonlinear genetic relationships between traits

Genomic prediction (GP) aims to predict the breeding values of multiple complex traits, usually assumed to be multivariate normally distributed by the largely used statistical methods, thus imposing linear genetic relationships between traits. Although these methods are valuable for GP they do not account for potential nonlinear genetic relationships between traits in scenarios. For individual traits, this oversight may minimally affect prediction accuracy, but it can limit genetic progress when selection involves multiple traits. Deep learning (DL) offers a promising alternative for capturing nonlinear genetic relationships due to its ability to identify complex patterns without prior assumptions about the data structure. We proposed a novel hybrid DLGBLUP model which uses the output of the traditional GBLUP, and enhances its predicted genetic values (PGV) by accounting for nonlinear genetic relationships between traits using DL. We simulated data with linear and nonlinear genetic relationships between traits in order to verify whether DLGBLUP was able to identify nonlinearity when present and avoid inducing it when absent. We found that DLGBLUP consistently provided more accurate PGV for traits simulated with strong nonlinear genetic relationships, accurately identifying these relationships. Over 7 generations of selection, a greater genetic progress was achieved with PGV that accounted for nonlinear relationships (DLGBLUP), compared with GBLUP. When applied to a real dataset from the French Holstein dairy cattle population, DLGBLUP detected nonlinear genetic relationships between pairs of traits, such as conception rate and protein content, and somatic cell count and fat yield, although, no significant increase in prediction accuracy was observed. The integration of DL into GP enabled the modeling of nonlinear genetic relationships between traits, a possibility not previously discussed, given the linear nature of GBLUP. The detection of nonlinear genetic relationships between traits in the French Holstein population when using DLGBLUP indicates the presence of such relationships in real breeding data, suggesting that it may be relevant to further explore nonlinear relationships. This possibility of nonlinear genetic relationships between traits offers a different perspective into multitrait evaluations, with potential to further improve selection strategies in commercial livestock breeding programs. This is particularly relevant when integrating new traits into multitrait evaluations or incorporating new subpopulations, which may introduce different forms of nonlinearity. Finally, it is shown that DL can be used as a complement to the statistical methods deployed in routine genetic evaluations, rather than as an alternative, by enhancing their performance.

Genome-wide association studies of dairy cattle resistance to digital dermatitis recorded at four distinct lactation stages

Digital dermatitis (DD) is an endemic infectious hoof disease causing lameness in dairy cattle. The aim of the present study was to investigate the genetic profile of DD development using phenotypic and genotypic data on 2192 Holstein cows. The feet of each cow were clinically examined four times: pre-calving, shortly after calving, near peak of milk production, and in late lactation. Presence or absence of disease and proportion of healthy feet per cow constituted two DD phenotypes of study. For each phenotype and timepoint of clinical examination, we conducted single-step genome-wide association analyses to identify individual markers and genomic regions linked to DD. We focused on the ten 1-Mb windows that explained the largest proportion of the total genetic variance as well as windows that enclosed significant markers. Functional enrichment analysis was also applied to determine functional candidate genes for DD. Significant (P < 0.05) genomic heritability estimates were derived ranging from 0.21 to 0.25. Results revealed two markers on chromosomes 7 and 15 that were related to both disease phenotypes. Furthermore, we identified three genomic windows on chromosome 14 and one window on chromosome 7 each explaining more than 1% of the trait additive genetic variance. Functional enrichment analysis revealed multiple promising candidate genes implicated in hoof health, wound healing, and inflammatory skin diseases. Collectively, our results provide novel insights into the biological mechanism of host resistance to DD development in dairy cattle and support genomic selection towards improving foot health.

Genomic selection strategies to overcome genotype by environment interactions in biosecurity-based aquaculture breeding programs

BACKGROUND

Family-based selective breeding programs typically employ both between-family and within-family selection in aquaculture. However, these programs may exhibit a reduced genetic gain in the presence of a genotype by environment interactions (G × E) when employing biosecurity-based breeding schemes (BS), compared to non-biosecurity-based breeding schemes (NBS). Fortunately, genomic selection shows promise in improving genetic gain by taking within-family variance into account. Stochastic simulation was employed to evaluate genetic gain and G × E trends in BS for improving the body weight of L. vannamei, considering selective genotyping strategies for test group (TG) at a commercial farm environment (CE), the number individuals of the selection group (SG) genotyped at nucleus breeding center (NE), and varying levels of G × E.

RESULTS

The loss of genetic gain in BS ranged from 9.4 to 38.9% in pedigree-based selection and was more pronounced when G × E was stronger, as quantified by a lower genetic correlation for body weight between NE and CE. Genomic selection, particularly with selective genotyping of TG individuals with extreme performance, effectively offset the loss of genetic gain. With a genetic correlation of 0.8, genotyping 20 SG individuals in each candidate family achieved 93.2% of the genetic gain observed for NBS. However, when the genetic correlation fell below 0.5, the number of genotyped SG individuals per family had to be increased to 50 or more. Genetic gain improved by on average 9.4% when the number of genotyped SG individuals rose from 20 to 50, but the increase in genetic gain averaged only 2.4% when expanding from 50 to 80 individuals genotyped. In addition, the genetic correlation decreased by on average 0.13 over 30 generations of selection when performing BS and the genetic correlation fluctuated across generations.

CONCLUSIONS

Genomic selection can effectively compensate for the loss of genetic gain in BS due to G × E. However, the number of genotyped SG individuals and the level of G × E significantly affected the extra genetic gain from genomic selection. A family-based BS selective breeding program should monitor the level of G × E and genotyping 50 SG individuals per candidate family to minimize the loss of genetic gain due to G × E, unless the level of G × E is confirmed to be low.

Weighted single-step genome-wide association study identified genomic regions and candidate genes for growth and reproductive traits in Wenchang chicken

Wenchang chicken is a native meat-type chicken breed in Hainan Province, China, known for its tender meat, which is in high demand both domestically and internationally. Improving important economic traits in Wenchang chicken is of significant importance to the poultry market. In this study, we genotyped 3737 individuals from two generations of Hainan Wenchang chicken using the Jingxin No. 1 55K SNP chip, and subsequently imputed the chip data. After quality control, 41,086 high-quality SNPs were obtained. We performed weighted single-step genome-wide association study (wssGWAS) on six traits using a single-trait model, and searched for candidate genes within QTL regions explaining more than 1 % of the genetic variance. Based on the contribution to genetic variation from QTL regions, we identified 4, 8, 5, 6, 9, and 8 candidate regions for EP, 35w-EW, BW, FT, SL, and AFE, respectively. The top three significant windows cumulatively explained 5.03 %, 5.18 %, 19.05 %, 5.87 %, 6.15 %, and 3.86 % of the genetic variation for these traits, respectively. Within these regions, we identified 14 genes, including LDB2 and GIT1, which had previously been reported to be associated with ovarian development, body weight, and fat differentiation. The biological processes and pathways involving these genes include cartilage development, lipid molecule synthesis and metabolism, angiogenesis, and cell proliferation and migration. LDB2 was proposed as a strong candidate gene due to its significant function and high consistency with other studies regarding various carcass traits in chickens. Our findings provide valuable information for broiler breeding based on molecular selection.

Invited review: Phenotyping strategies and genetic background of dairy cattle behavior in intensive production systems-From trait definition to genomic selection

Understanding and assessing dairy cattle behavior is critical for developing sustainable breeding programs and management practices. The behavior of individual animals can provide valuable information on their health and welfare status, improve reproductive management, and predict efficiency traits such as feed efficiency and milking efficiency. Routine genetic evaluations of animal behavior traits can contribute to optimizing breeding and management strategies for dairy cattle but require the identification of traits that capture the most important biological processes involved in behavioral responses. These traits should be heritable, repeatable, and measured in noninvasive and cost-effective ways in many individuals from the breeding populations or related reference populations. Although behavior traits are heritable in dairy cattle populations, they are highly polygenic, with no known major genes influencing their phenotypic expression. Genetically selecting dairy cattle based on their behavior can be advantageous because of their relationship with other key traits such as animal health, welfare, and productive efficiency, as well as animal and handler safety. Trait definition and longitudinal data collection are still key challenges for breeding for behavioral responses in dairy cattle. However, the more recent developments and adoption of precision technologies in dairy farms provide avenues for more objective phenotyping and genetic selection of behavior traits. Furthermore, there is still a need to standardize phenotyping protocols for existing traits and develop guidelines for recording novel behavioral traits and integrating multiple data sources. This review gives an overview of the most common indicators of dairy cattle behavior, summarizes the main methods used for analyzing animal behavior in commercial settings, describes the genetic and genomic background of previously defined behavioral traits, and discusses strategies for breeding and improving behavior traits coupled with future opportunities for genetic selection for improved behavioral responses.

Single-step genomic BLUP (ssGBLUP) effectively models small cattle populations: lessons from the Israeli-Holstein Herdbook

BACKGROUND

Routine genomic-estimated breeding values (gEBVs) are computed for the Israeli dairy cattle population by a two-step methodology in combination with the much larger Dutch population. Only sire genotypes are included. This work evaluated the contribution of cow genotypes obtained from the Israeli Holstein population to enhance gEBVs predictions via single-step genomic best-linear unbiased prediction (ssGBLUP). The gEBV values of 141 bulls with daughter information and high reliabilities for 305-day lactation yield of milk, fat, and protein were compared with the bulls' predicted ssGBLUP-gEBVs using a truncated dataset omitting production data of the last five years. We investigated how these sire gEBVs were affected by varying polygenic weights in the genomic relationship matrices and by deleting old phenotypic or genotypic records.

RESULTS

The correlations of the predicted gEBVs for milk, fat and protein computed from the truncated data with the current gEBVs based also on daughter records of the last five years were 0.64, 0.57, and 0.56, respectively, for a polygenic weight of 0.5, similar to the values achieved by the current two-step methodology. The regressions of the current gEBVs on the predicted values were 0.9 for milk and 0.7 for fat and protein. Genotyping of 1.8-5 cows had the approximate statistical power of one additional bull depending on the trait. Omitting phenotype records earlier than 2000 resulted in similar gEBV values. Omitting genotypes before 1995 improved the regression coefficients. For all experiments, varying the polygenic weights over the range of 0.1 to 0.9 resulted in a trade-off between correlations and overestimation of gEBVs for young bulls.

CONCLUSIONS

The model suffers from overestimation of the predicted values for young bulls. The time interval used for inclusion of genotypic and phenotypic records and adjustment of the polygenic weight can improve gEBV predictions and should be tuned to fit the tested population. For relatively small populations, genotyping of cows can significantly increase the reliability of gEBVs computed by single-step methodology. By extrapolation of our results, records of ~ 13,000 genotyped cows should provide a sufficiently large training population to obtain reliable estimates of gEBVs using ssGBLUP.

The Effect of Preselection on the Level of Bias and Accuracy in a Broiler Breeder Population, a Simulation Study

Many breeding programmes have to perform preselection, as genotyping and phenotyping all potential breeder candidates is often not a feasible option. There is need to understand how preselection affects the quality of the genomic estimated breeding values (EBVs) at final selection and thereby can affect genetic progress. This simulation study evaluated nine different preselection strategies in a broiler breeder programme and their effect on the quality of the (genomic) EBVs and genetic progress for three different traits: body weight (Body Weight), residual feed intake (RFI) and body weight gain (Gain). All birds have Body Weight recorded at preselection, but only the preselected birds were phenotyped for RFI and Gain and genotyped. The following criteria and intensities were studied: preselection based on phenotypic Body Weight (P), on a BLUP index (B) or on an ssGBLUP Index (G). Additionally, all criteria were studied with three different selection intensities, 10% of the males and 30% of the females (P10, B10, G10), 15% of the males and 45% of the females (P15, B15, G15) and 20% of the males and 60% of the females (P20, B20, G20). The accuracy at preselection with G10 was more accurate than B10 for both RFI and Gain (0.71 vs. 0.58 and 0.65 vs. 0.55 respectively), and also G15 was more accurate than B15 for both RFI and Gain (0.72 vs. 0.63 and 0.67 vs. 0.64 respectively); thus, the difference in accuracy reduces with an increasing number of birds being preselected. Differences in accuracy at final selection were mostly notable in the RFI trait between P10, B10 and G10, where G10 showed the highest accuracy (0.82 vs. 0.84 vs. 0.86 respectively). This difference in accuracy for RFI disappeared when more animals were preselected. For Body Weight and Gain, the BLUP preselection resulted in the highest accuracy at final selection when selection intensity decreased. The dispersion bias of EBVs at final selection was most pronounced in the P10 and P15 for Body Weight (0.81 and 0.92) but disappeared at P20 (0.97). The dispersion bias for all other criteria and traits was relatively small. Genetic progress was mostly affected when P10 or P15 was used at preselection, where the progress in Body Weight was noticeably higher, but prominently lower for RFI and Gain. The BLUP and ssGBLUP preselection had very similar genetic progress across traits and showed comparable improvements in the selection index. In conclusion, with high preselection intensity, the use of ssGBLUP at preselection might be favoured as there is an improvement in genetic progress across traits in all scenarios, which is due to the increased preselection accuracy. When preselection intensity decreases, the benefit of using ssGBLUP over BLUP at preselection disappears.

Weighted single step GWAS reveals genomic regions associated with economic traits in Murrah buffaloes

The objective of the present study was to identify genomic regions influencing economic traits in Murrah buffaloes using weighted single step Genome Wide Association Analysis (WssGWAS). Data on 2000 animals, out of which 120 were genotyped using a double digest Restriction site Associated DNA (ddRAD) sequencing approach. The phenotypic data were collected from NDRI, India, on growth traits, viz., body weight at 6M (month), 12M, 18M and 24M, production traits like 305D (day) milk yield, lactation length (LL) and dry period (DP) and reproduction traits like age at first calving (AFC), calving interval (CI) and first service period (FSP). The biallelic genotypic data consisted of 49353 markers post-quality check. The heritability estimates were moderate to high, low to moderate, low for growth, production, reproduction traits, respectively. Important genomic regions explaining more than 0.5% of the total additive genetic variance explained by 30 adjacent SNPs were selected for further analysis of candidate genes. In this study, 105 genomic regions were associated with growth, 35 genomic regions with production and 42 window regions with reproduction traits. Different candidate genes were identified in these genomic regions, of which important are OSBPL8, NAP1L1 for growth, CNTNAP2 for production and ILDR2, TADA1 and POGK for reproduction traits.

Enhancing animal breeding through quality control in genomic data - a review

High-throughput genotyping and sequencing has revolutionized animal breeding by providing access to vast amounts of genomic data to facilitate precise selection for desirable traits. This shift from traditional methods to genomic selection provides dense marker information for predicting genetic variants. However, the success of genomic selection heavily depends on the accuracy and quality of the genomic data. Inaccurate or low-quality data can lead to flawed predictions, compromising breeding programs and reducing genetic gains. Therefore, stringent quality control (QC) measures are essential at every stage of data processing. QC in genomic data involves managing single nucleotide polymorphism (SNP) quality, assessing call rates, and filtering based on minor allele frequency (MAF) and Hardy-Weinberg equilibrium (HWE). High-quality SNP data is crucial because genotyping errors can bias the estimates of breeding values. Cost-effective low-density genotyping platforms often require imputation to deduce missing genotypes. QC is vital for genomic selection, genome-wide association studies (GWAS), and population genetics analyses because it ensures data accuracy and reliability. This paper reviews QC strategies for genomic data and emphasizes their applications in animal breeding programs. By examining various QC tools and methods, this review highlights the importance of data integrity in achieving successful outcomes in genomic selection, GWAS, and population analyses. Furthermore, this review covers the critical role of robust QC measures in enhancing the reliability of genomic predictions and advancing animal breeding practices.

Genomic prediction of male fertility in Brown Swiss cattle

Bull fertility has been recognized as an important factor affecting dairy herd fertility. The objective of this study was to assess the feasibility of predicting male fertility in Brown Swiss cattle using genomic data. The dataset consisted of 1,102 Italian Brown Swiss bulls with sire conception rate (SCR) records and genotype data for roughly 480k SNP. The analyses included the use of linear kernel-based regression models fitting all SNPs or incorporating markers with large effect. Predictive performance was evaluated in 5-fold cross-validation using the correlation between observed and predicted SCR values and mean squared error of prediction. The entire SNP set exhibited predictive correlations around 0.19. Interestingly, the inclusion of 2 markers with large effect yielded predictive correlations around 0.32. Overall, using linear kernel-based models fitting markers with large effect is a promising approach. Our findings could help Brown Swiss breeders make enhanced genome-guided management and selection decisions on male fertility.

Genomic selection in forest trees comes to life: unraveling its potential in an advanced four-generation Eucalyptus grandis population

Genomic Selection (GS) in tree breeding optimizes genetic gains by leveraging genomic data to enable early selection of seedlings without phenotypic data reducing breeding cycle and increasing selection intensity. Traditional assessments of the potential of GS in forest trees have typically focused on model performance using cross-validation within the same generation but evaluating effectively realized predictive ability (RPA) across generations is crucial. This study estimated RPAs for volume growth (VOL), wood density (WD), and pulp yield (PY) across four generations breeding of Eucalyptus grandis. The training set spanned three generations, including 34,461 trees with three-year growth data, 6,014 trees with wood quality trait data, and 1,918 trees with 12,695 SNPs (single nucleotide polymorphisms) data. Employing single-step genomic BLUP, we compared the genomic predictions of breeding values (GEBVs) for 1,153 fourth-generation full-sib seedlings in the greenhouse with their later-collected phenotypic estimated breeding values (EBVs) at age three years. RPAs were estimated using three GS targets (individual trees, trees within families, and families), two selection criteria (single- and multiple-trait), and training populations of either all 1,918 genotyped trees or the 67 direct ancestors of the selection candidates. RPAs were higher for wood quality traits (0.33 to 0.59) compared to VOL (0.14 to 0.19) and improved for wood traits (0.42 to 0.75) but not for VOL when trained only with direct ancestors, highlighting the challenges in accurately predicting growth traits. GS was more effective at excluding bottom-ranked candidates than selecting top-ranked ones. The between-family GS approach outperformed individual-tree selection for VOL (0.11 to 0.16) and PY (0.72 to 0.75), but not for WD (0.43 vs. 0.42). Furthermore, higher levels of relatedness and lower genotype by environment (G × E) interaction between training and testing populations enhanced RPAs for VOL (0.39). In summary, despite limited effectiveness in ranking top VOL individuals, GS effectively identified low-performing individuals and families. These multi-generational findings underscore GS's potential in tree breeding, stressing the importance of considering relatedness and G × E interaction for optimal performance.

Comparison of genomic prediction accuracy using different models for egg production traits in Taiwan country chicken

In local chickens targeted for niche markets, genotyping costs are relatively high due to the small population size and diverse breeding goals. The single-step genomic best linear unbiased prediction (ssGBLUP) model, which combines pedigree and genomic information, has been introduced to increase the accuracy of genomic estimated breeding value (GEBV). Therefore, this model may be more beneficial than the genomic BLUP (GBLUP) model for genomic selection in local chickens. Additionally, the single-step genome-wide association study (ssGWAS) can be used to extend the ssGBLUP model results to animals with available phenotypic information but without genotypic data. In this study, we compared the accuracy of (G)EBVs using the pedigree-based BLUP (PBLUP), GBLUP, and ssGBLUP models. Moreover, we conducted single-SNP GWAS (SNP-GWAS), GBLUP-GWAS, and ssGWAS methods to identify genes associated with egg production traits in the NCHU-G101 chicken to understand the feasibility of using genomic selection in a small population. The average prediction accuracy of (G)EBV for egg production traits using the PBLUP, GBLUP, and ssGBLUP models is 0.536, 0.531, and 0.555, respectively. In total, 22 suggestive- and 5% Bonferroni genome-wide significant-level SNPs for total egg number (EN), average laying rate (LR), average clutch length, and total clutch number are detected using 3 GWAS methods. These SNPs are mapped onto Gallus gallus chromosomes (GGA) 4, 6, 10, 18, and 25 in NCHU-G101 chicken. Furthermore, through SNP-GWAS and ssGWAS methods, we identify 2 genes on GGA4 associated with EN and LR: ENSGALG00000023172 and PPARGC1A. In conclusion, the ssGBLUP model demonstrates superior prediction accuracy, performing on average 3.41% than the PBLUP model. The implications of our gene results may guide future selection strategies for Taiwan Country chickens. Our results highlight the applicability of the ssGBLUP model for egg production traits selection in a small population, specifically NCHU-G101 chicken in Taiwan.

Ovarian function and endocrine phenotypes of lactating dairy cows during the estrous cycle are associated with genomic-enhanced predictions of fertility potential

The objectives of this prospective cohort study were to characterize associations among genomic merit for fertility with ovarian and endocrine function and the estrous behavior of dairy cows during an entire nonhormonally manipulated estrous cycle. Lactating Holstein cows entering their first (n = 82) or second (n = 37) lactation had ear-notch tissue samples collected for genotyping using a commercial genomic test. Based on genomic predicted transmitting ability values for daughter pregnancy rate (gDPR), cows were classified into high (Hi-Fert; gDPR > 0.6, n = 36), medium (Med-Fert; gDPR -1.3 to 0.6, n = 45), and low fertility (Lo-Fert; gDPR < -1.3, n = 38) groups. At 33 to 39 DIM, cohorts of cows were enrolled in the Presynch-Ovsynch protocol for synchronization of ovulation and initiation of a new estrous cycle. Thereafter, the ovarian function and endocrine dynamics were monitored daily until the next ovulation by transrectal ultrasonography and concentrations of progesterone (P4), estradiol, and FSH. Estrous behavior was monitored with an ear-attached automated estrus detection system that recorded physical activity and rumination time. Overall, we observed an association between fertility group and the ovarian and hormonal phenotype of dairy cows during the estrous cycle. Cows in the Hi-Fert group had greater circulating concentrations of P4 than cows in the Lo-Fert group from d 4 to 13 after induction of ovulation and from day -3 to -1 before the onset of luteolysis. The frequency of atypical estrous cycles was 3-fold greater for cows in the Lo-Fert than the Hi-Fert group. We also observed other modest associations between genomic merit for fertility with the follicular dynamics and estrous behavior. We found several associations between milk yield and parity with ovarian, endocrine, and estrous behavior phenotypes as cows with greater milk yield and in the second lactation were more likely to have unfavorable phenotypes. These results demonstrate that differences in reproductive performance between cows of different genomic merit for fertility classified based on gDPR may be partially associated with circulating concentrations of P4, the incidence of atypical phenotypes during the estrous cycles, and, to a lesser extent, the follicular wave dynamics. The observed physiological and endocrine phenotypes might help explain part of the differences in reproductive performance between cows of superior and inferior genomic merit for fertility.

An overview of recent technological developments in bovine genomics

Cattle are regarded as highly valuable animals because of their milk, beef, dung, fur, and ability to draft. The scientific community has tried a number of strategies to improve the genetic makeup of bovine germplasm. To ensure higher returns for the dairy and beef industries, researchers face their greatest challenge in improving commercially important traits. One of the biggest developments in the last few decades in the creation of instruments for cattle genetic improvement is the discovery of the genome. Breeding livestock is being revolutionized by genomic selection made possible by the availability of medium- and high-density single nucleotide polymorphism (SNP) arrays coupled with sophisticated statistical techniques. It is becoming easier to access high-dimensional genomic data in cattle. Continuously declining genotyping costs and an increase in services that use genomic data to increase return on investment have both made a significant contribution to this. The field of genomics has come a long way thanks to groundbreaking discoveries such as radiation-hybrid mapping, in situ hybridization, synteny analysis, somatic cell genetics, cytogenetic maps, molecular markers, association studies for quantitative trait loci, high-throughput SNP genotyping, whole-genome shotgun sequencing to whole-genome mapping, and genome editing. These advancements have had a significant positive impact on the field of cattle genomics. This manuscript aimed to review recent advances in genomic technologies for cattle breeding and future prospects in this field.

Investigating the genomic background of calving-related traits in Canadian Jersey cattle

Traits related to calving have a significant impact on animal welfare and farm profitability in dairy production systems. Identifying genomic regions associated with calving traits could contribute to refining dairy cattle breeding programs and management practices in the dairy industry. Therefore, the primary objectives of this study were to estimate genetic parameters and perform genome-wide association studies (GWAS) and functional enrichment analyses for stillbirth, gestation length, calf size, and calving ease traits in North American Jersey cattle. A total of 40,503 animals with phenotypic records and 5,398 animals genotyped for 45,101 single nucleotide polymorphisms (SNPs) were included in the analyses. Genetic parameters were estimated based on animal models and Bayesian methods. The effects of SNPs were estimated using the Single-step Genomic Best Linear Unbiased Prediction (ssGBLUP) method. The heritability (standard error) estimates ranged from 0.01 (0.01) for stillbirths (SB) in heifers to 0.11 (0.01) for gestation length (GL) in cows. The genetic correlations ranged from -0.58 (0.11) between calving ease (CE) and SB in heifers to 0.44 (0.14) between calving ease and calf size (CZ) in cows. CE showed the highest genetic correlation between heifers and cows, 0.8 (0.22) respectively. The candidate genes identified, including MTHFR, SERPINA5, IGFBP3, and ZRANB1, are involved in key biological processes and metabolic pathways related to the studied traits. Reducing environmental variation and identifying novel indicators of reproduction traits in the Jersey breed are needed given the low heritability estimates for most traits evaluated in this study. In conclusion, this study provides a characterization of the genetic background of calving-related traits in Jersey cattle. The estimates obtained can be used to improve or build selection indexes in Jersey cattle breeding programs in North America.

Residual network improves the prediction accuracy of genomic selection

Genetic improvement of complex traits in animal and plant breeding depends on the efficient and accurate estimation of breeding values. Deep learning methods have been shown to be not superior over traditional genomic selection (GS) methods, partially due to the degradation problem (i.e. with the increase of the model depth, the performance of the deeper model deteriorates). Since the deep learning method residual network (ResNet) is designed to solve gradient degradation, we examined its performance and factors related to its prediction accuracy in GS. Here we compared the prediction accuracy of conventional genomic best linear unbiased prediction, Bayesian methods (BayesA, BayesB, BayesC, and Bayesian Lasso), and two deep learning methods, convolutional neural network and ResNet, on three datasets (wheat, simulated and real pig data). ResNet outperformed other methods in both Pearson's correlation coefficient (PCC) and mean squared error (MSE) on the wheat and simulated data. For the pig backfat depth trait, ResNet still had the lowest MSE, whereas Bayesian Lasso had the highest PCC. We further clustered the pig data into four groups and, on one separated group, ResNet had the highest prediction accuracy (both PCC and MSE). Transfer learning was adopted and capable of enhancing the performance of both convolutional neural network and ResNet. Taken together, our findings indicate that ResNet could improve GS prediction accuracy, affected potentially by factors such as the genetic architecture of complex traits, data volume, and heterogeneity.

Genetic background of walking ability and its relationship with leg defects, mortality, and performance traits in turkeys (Meleagris gallopavo)

This study aimed to explore the genetic basis of walking ability and potentially related performance traits in turkey purebred populations. Phenotypic, pedigree, and genomic datasets from 2 turkey lines hatched between 2010 and 2023 were included in the study. Walking ability data, defined based on a scoring system ranging from 1 (worst) to 6 (best), were collected on 192,019 animals of a female line and 235,461 animals of a male line. Genomic information was obtained for 46,427 turkeys (22,302 from a female line and 24,125 from a male line) using a 65K single nucleotide polymorphism (SNP) panel. Variance components and heritability for walking ability were estimated. Furthermore, genetic and phenotypic correlations among walking ability, mortality and disorders, and performance traits were calculated. A genome-wide association study (GWAS) was also conducted to identify SNPs associated with walking ability. Walking ability is moderately heritable (0.23 ± 0.01) in both turkey lines. The genetic correlations between walking ability and the other evaluated traits ranged from -0.02 to -0.78, with leg defects exhibiting the strongest negative correlation with walking ability. In the female line, 31 SNPs were associated with walking ability and overlapped with 116 genes. These positional genes are linked to 6 gene ontology (GO) terms. Notably, genes such as CSRP2, DDX1, RHBDL1, SEZ6L, and CTSK are involved in growth, development, locomotion, and bone disorders. GO terms, including fibronectin binding (GO:0001968), peptide cross-linking (GO:0018149), and catabolic process (GO:0009057), are directly linked with mobility. In the male line, 66 markers associated with walking ability were identified and overlapped with 281 genes. These genes are linked to 12 GO terms. Genes such as RB1CC1, TNNI1, MSTN, FN1, SIK3, PADI2, ERBB4, B3GNT2, and BACE1 are associated with cell growth, myostatin development, and disorders. GO terms in the male line are predominantly related to lipid metabolism. In conclusion, walking ability is moderately heritable in both populations. Furthermore, walking ability can be enhanced through targeted genetic selection, emphasizing its relevance to both animal welfare and productivity.

Genetic parameters and trends for Johne's disease in US Holsteins: An updated study

Johne's disease (JD) is an infectious enteric disease in ruminants, causing substantial economic loss annually worldwide. This work aimed to estimate JD's genetic parameters and the phenotypic and genetic trends by incorporating recent data. It also explores the feasibility of a national genetic evaluation for JD susceptibility in Holstein cattle in the United States. The data were extracted from a JD data repository, maintained at the Council on Dairy Cattle Breeding, and initially supplied by 2 dairy record processing centers. The data comprised 365,980 Holstein cows from 1,048 herds participating in a voluntary control program for JD. Two protocol kits, IDEXX Paratuberculosis Screening Ab Test (IDX) and Parachek 2 (PCK), were used to analyze milk samples with the ELISA technique. Test results from the first 5 parities were considered. An animal was considered infected if it had at least one positive outcome. The overall average of JD incidence was 4.72% in these US Holstein cattle. Genotypes of 78,964 SNP markers were used for 25,000 animals randomly selected from the phenotyped population. Variance components and genetic parameters were estimated based on 3 models, namely, a pedigree-only threshold model (THR), a single-step threshold model (ssTHR), and a single-step linear model (ssLR). The posterior heritability estimates of JD susceptibility were low to moderate: 0.11 to 0.16 based on the 2 threshold models and 0.05 to 0.09 based on the linear model. The average reliability of EBVs of JD susceptibility using single-step analysis for animals with or without phenotypes varied from 0.18 (THR) to 0.22 (ssLR) for IDX and from 0.14 (THR) to 0.18 (ssTHR and ssLR) for PCK. Despite no prior direct genetic selection against JD, the estimated genetic trends of JD susceptibility were negative and highly significant. The correlations of bulls' PTA with economically important traits such as milk yield, milk protein, milk fat, somatic cell score, and mastitis were low, indicating a nonoverlapping genetic selection process with traits in current genetic evaluations. Our results suggest the feasibility of reducing the JD incidence rate by incorporating it into the national genetic evaluation programs.

Reproductive physiological outcomes of dairy cows with different genomic merit for fertility: biomarkers, uterine health, endocrine status, estrus features, and response to ovarian synchronization

Our overarching objective was to characterize associations between genomic merit for fertility and the reproductive function of lactating dairy cows in a prospective cohort study. In this manuscript, we present results of the association between genomic merit for fertility and indicators of metabolic status and inflammation, uterine health, endocrine status, response to synchronization, and estrous behavior in dairy cows. Lactating Holstein cows entering their first (n = 82) or second (n = 37) lactation were enrolled at parturition and fitted with an ear-attached sensor for automated detection of estrus. Ear-notch tissue samples were collected from all cows and submitted for genotyping using a commercial genomic test. Based on genomic predicted transmitting ability values for daughter pregnancy rate (gDPR) cows were classified into a high (Hi-Fert; gDPR > 0.6; n = 36), medium (Med-Fert; gDPR -1.3 to 0.6; n = 45), and low (Lo-Fert; gDPR < -1.3; n = 38) group. At 33 to 39 d in milk (DIM), cohorts of cows were enrolled in the Presynch-Ovsynch protocol for synchronization of estrus and ovulation. Body weights, body condition scores (BCS), and uterine health measurements (i.e., vaginal discharge, uterine cytology) were collected from parturition to 60 DIM and milk yield was collected through 90 DIM. Blood samples were collected weekly through 3 wk of lactation for analysis of β-hydroxybutyrate, nonesterified fatty acids, and haptoglobin plasma concentrations. Body weight, BCS, NEFA, BHB, and Haptoglobin were not associated with fertility groups from 1 to 9 wk after parturition. The proportion of cows classified as having endometritis at 33 to 36 DIM tended to be greater for the Lo-Fert than the Hi-Fert group. The proportion of cows that resumed cyclicity did not differ at any time point evaluated and there were no significant associations between probability or duration and intensity of estrus with fertility group. Cows of superior genetic merit for fertility were more likely to ovulate, have a functional CL, have greater circulating P4, and have larger ovulatory size than cows of inferior fertility potential at key time points during synchronization of estrus and ovulation. Despite observing numerical differences with potential performance consequences for the proportion of cows that responded to synchronization of ovulation and were both cyclic and responded to the Ovsynch portion of the synchronization protocol, we did not observe significant differences between fertility groups. Although not consistent and modest in magnitude, the collective physiological and endocrine differences observed suggested that cows of superior genetic fertility potential might have improved reproductive performance, at least in part, because of modestly improved endocrine status, uterine health, and ability to ovulate.

Genetic factors underlying host resistance to Rhipicephalus microplus tick infestation in Braford cattle: a systems biology perspective

Approximately 80% of the world's cattle are raised in regions with a high risk of tick-borne diseases, resulting in significant economic losses due to parasitism by Rhipicephalus (Boophilus) microplus. However, the lack of a systemic biology approach hampers a comprehensive understanding of tick-host interactions that mediate tick resistance phenotypes. Here, we conducted a genome-wide association study (GWAS) of 2933 Braford cattle and found 340 single-nucleotide polymorphisms (SNPs) associated with tick counts. Gene expression analyses were performed on skin samples obtained from previously tick-exposed heifers with extremely high or low estimated breeding values for R. microplus counts. Evaluations were performed both before and after artificial infestation with ticks. Differentially expressed genes were found within 1-Mb windows centered at significant SNPs from GWAS. A total of 330 genes were related to the breakdown of homeostasis that was induced by larval attachment to bovine skin. Enrichment analysis pointed to a key role of proteolysis and signal transduction via JAK/STAT, NFKB and WNT/beta catenin signaling pathways. Integrative analysis on matrixEQTL revealed two cis-eQTLs and four significant SNPs in the genes peptidyl arginine deiminase type IV (PADI4) and LOC11449251. The integration of genomic data from QTL maps and transcriptome analyses has identified a set of twelve key genes that show significant associations with tick loads. These genes could be key candidates to improve the accuracy of genomic predictions for tick resistance in Braford cattle.

Accuracy of Genomic prediction for fleece traits in Inner Mongolia Cashmere goats

The fleece traits are important economic traits of goats. With the reduction of sequencing and genotyping cost and the improvement of related technologies, genomic selection for goats has become possible. The research collect pedigree, phenotype and genotype information of 2299 Inner Mongolia Cashmere goats (IMCGs) individuals. We estimate fixed effects, and compare the estimates of variance components, heritability and genomic predictive ability of fleece traits in IMCGs when using the pedigree based Best Linear Unbiased Prediction (ABLUP), Genomic BLUP (GBLUP) or single-step GBLUP (ssGBLUP). The fleece traits considered are cashmere production (CP), cashmere diameter (CD), cashmere length (CL) and fiber length (FL). It was found that year of production, sex, herd and individual ages had highly significant effects on the four fleece traits (P < 0.01). All of these factors should be considered when the genetic parameters of fleece traits in IMCGs are evaluated. The heritabilities of FL, CL, CP and CD with ABLUP, GBLUP and ssGBLUP methods were 0.26 ~ 0.31, 0.05 ~ 0.08, 0.15 ~ 0.20 and 0.22 ~ 0.28, respectively. Therefore, it can be inferred that the genetic progress of CL is relatively slow. The predictive ability of fleece traits in IMCGs with GBLUP (56.18% to 69.06%) and ssGBLUP methods (66.82% to 73.70%) was significantly higher than that of ABLUP (36.73% to 41.25%). For the ssGBLUP method is significantly (29% ~ 33%) higher than that with ABLUP, and which is slightly (4% ~ 14%) higher than that of GBLUP. The ssGBLUP will be as an superiors method for using genomic selection of fleece traits in Inner Mongolia Cashmere goats.

Genetic parameters for calf feeding traits derived from automated milk feeding machines and number of bovine respiratory disease treatments in North American Holstein calves

Precision livestock farming technologies, such as automatic milk feeding machines, have increased the availability of on-farm data collected from dairy operations. We analyzed feeding records from automatic milk feeding machines to evaluate the genetic background of milk feeding traits and bovine respiratory disease (BRD) in North American Holstein calves. Data from 10,076 preweaning female Holstein calves were collected daily over a period of 6 yr (3 yr included per-visit data), and daily milk consumption (DMC), per-visit milk consumption (PVMC), daily sum of drinking duration (DSDD), drinking duration per-visit, daily number of rewarded visits (DNRV), and total number of visits per day were recorded over a 60-d preweaning period. Additional traits were derived from these variables, including total consumption and duration variance (TCV and TDV), feeding interval, drinking speed (DS), and preweaning stayability. A single BRD-related trait was evaluated, which was the number of times a calf was treated for BRD (NTT). The NTT was determined by counting the number of BRD incidences before 60 d of age. All traits were analyzed using single-step genomic BLUP mixed-model equations and fitting either repeatability or random regression models in the BLUPF90+ suite of programs. A total of 10,076 calves with phenotypic records and genotypic information for 57,019 SNP after the quality control were included in the analyses. Feeding traits had low heritability estimates based on repeatability models (0.006 ± 0.0009 to 0.08 ± 0.004). However, total variance traits using an animal model had greater heritabilities of 0.21 ± 0.023 and 0.23 ± 0.024, for TCV and TDV, respectively. The heritability estimates increased with the repeatability model when using only the first 32 d preweaning (e.g., PVMC = 0.040 ± 0.003, DMC = 0.090 ± 0.009, DSDD = 0.100 ± 0.005, DS = 0.150 ± 0.007, DNRV = 0.020 ± 0.002). When fitting random regression models (RRM) using the full dataset (60-d period), greater heritability estimates were obtained (e.g., PVMC = 0.070 [range: 0.020, 0.110], DMC = 0.460 [range: 0.050, 0.680], DSDD = 0.180 [range: 0.010, 0.340], DS = 0.19 [range: 0.070, 0.430], DNRV = 0.120 [range: 0.030, 0.450]) for the majority of the traits, suggesting that RRM capture more genetic variability than the repeatability model with better fit being found for RRM. Moderate negative genetic correlations of -0.59 between DMC and NTT were observed, suggesting that automatic milk feeding machines records have the potential to be used for genetically improving disease resilience in Holstein calves. The results from this study provide key insights of the genetic background of early in-life traits in dairy cattle, which can be used for selecting animals with improved health outcomes and performance.

Candidate genes for longitudinal traits under sequential sampling in beef cattle

Both the measurement age of a longitudinal trait and the common pre-sampling procedures used in beef cattle herds may affect the identification of a functional candidate gene (FCG) that is potentially associated with a trait. To identify the FCG that takes part in the genetic control of body weight at five different ages in a beef cattle population with and without sequential sampling, the animals were weighed at different measurement events, around 330, 385, 440, 495 and 550 days old. Genetic parameters were estimated for body weight at each age using a single trait (STM) and a random regression model (RRM). In addition, two different databases were used to estimate the genetic parameters: the first (DB100) was formed by all animals that were weighed in the five measurement events, and the second (DB70) has records of the same population, considering that 70% of the heaviest animals were selected after each measurement event. For DB100, genome-wide association studies (GWAS) were performed with 21,667 SNP markers to identify genomic windows that explained at least 1% of the genetic variance. Additionally, prioritization analyses were performed and FCGs were selected. We associated seven different FCGs with body weight at different ages. Among them, the gene DUSP10 was suggested as FCG in all five ages evaluated. Genetic parameters estimated for body weight using DB100 were similar when STM and RRM were applied. However, when DB70 was used as phenotypic data, there were differences between the two models. When the STM was applied, there were differences between the genetic parameters estimated for body weight when DB100 or DB70 were used as sources of phenotypes, but not for the estimates obtained with RRM. The importance of each gene for animal growth can change at different ages, and different genes may be more relevant to body weight at each different growth stage for beef cattle. Besides, sequential sampling can affect the GWAS results of a longitudinal trait. The age of the animal when a longitudinal trait is measured and pre-sampling can also contribute to inconsistencies in GWAS results for body weight in beef cattle, depending on the time when that data were collected, and consequently on the identification of FCG between studies, even when models that consider a covariance structure are used.

Expected values for the accuracy of predicted breeding values accounting for genetic differences between reference and target populations

BACKGROUND

Genetic merit, or breeding values as referred to in livestock and crop breeding programs, is one of the keys to the successful selection of animals in commercial farming systems. The developments in statistical methods during the twentieth century and single nucleotide polymorphism (SNP) chip technologies in the twenty-first century have revolutionized agricultural production, by allowing highly accurate predictions of breeding values for selection candidates at a very early age. Nonetheless, for many breeding populations, realized accuracies of predicted breeding values (PBV) remain below the theoretical maximum, even when the reference population is sufficiently large, and SNPs included in the model are in sufficient linkage disequilibrium (LD) with the quantitative trait locus (QTL). This is particularly noticeable over generations, as we observe the so-called erosion of the effects of SNPs due to recombinations, accompanied by the erosion of the accuracy of prediction. While accurately quantifying the erosion at the individual SNP level is a difficult and unresolved task, quantifying the erosion of the accuracy of prediction is a more tractable problem. In this paper, we describe a method that uses the relationship between reference and target populations to calculate expected values for the accuracies of predicted breeding values for non-phenotyped individuals accounting for erosion. The accuracy of the expected values was evaluated through simulations, and a further evaluation was performed on real data.

RESULTS

Using simulations, we empirically confirmed that our expected values for the accuracy of PBV accounting for erosion were able to correctly determine the prediction accuracy of breeding values for non-phenotyped individuals. When comparing the expected to the realized accuracies of PBV with real data, only one out of the four traits evaluated presented accuracies that were significantly higher than the expected, approachingh 2 .

CONCLUSIONS

We defined an index of genetic correlation between reference and target populations, which summarizes the expected overall erosion due to differences in allele frequencies and LD patterns between populations. We used this correlation along with a trait's heritability to derive expected values for the accuracy ( R ) of PBV accounting for the erosion, and demonstrated that our derived E R | erosion is a reliable metric.

Review of sheep breeding and genetic research in Türkiye

The livestock industry in Türkiye is vital to the country's agricultural sector and economy. In particular, sheep products are an important source of income and livelihood for many Turkish smallholder farmers in semi-arid and highland areas. Türkiye is one of the largest sheep producers in the world and its sheep production system is heavily dependent on indigenous breeds. Given the importance of the sheep industry in Türkiye, a systematic literature review on sheep breeding and genetic improvement in the country is needed for the development and optimization of sheep breeding programs using modern approaches, such as genomic selection. Therefore, we conducted a comprehensive literature review on the current characteristics of sheep populations and farms based on the most up-to-date census data and breeding and genetic studies obtained from scientific articles. The number of sheep has increased in recent years, mainly due to the state's policy of supporting livestock farming and the increase in consumer demand for sheep dairy products with high nutritional and health benefits. Most of the genetic studies on indigenous Turkish sheep have been limited to specific traits and breeds. The use of genomics was found to be incipient, with genomic analysis applied to only two major breeds for heritability or genome-wide association studies. The scope of heritability and genome-wide association studies should be expanded to include traits and breeds that have received little or no attention. It is also worth revisiting genetic diversity studies using genome-wide single nucleotide polymorphism markers. Although there was no report of genomic selection in Turkish sheep to date, genomics could contribute to overcoming the difficulties of implementing traditional pedigree-based breeding programs that require accurate pedigree recording. As indigenous sheep breeds are better adapted to the local environmental conditions, the proper use of breeding strategies will contribute to increased income, food security, and reduced environmental footprint in a sustainable manner.

Genetic background of juniper (Juniperus spp.) consumption predicted by fecal near-infrared spectroscopy in divergently selected goats raised in harsh rangeland environments

BACKGROUND

Junipers (Juniperus spp.) are woody native, invasive plants that have caused encroachment problems in the U.S. western rangelands, decreasing forage productivity and biodiversity. A potential solution to this issue is using goats in targeted grazing programs. However, junipers, which grow in dry and harsh environmental conditions, use chemical defense mechanisms to deter herbivores. Therefore, genetically selecting goats for increased juniper consumption is of great interest for regenerative rangeland management. In this context, the primary objectives of this study were to: 1) estimate variance components and genetic parameters for predicted juniper consumption in divergently selected Angora (ANG) and composite Boer x Spanish (BS) goat populations grazing on Western U.S. rangelands; and 2) to identify genomic regions, candidate genes, and biological pathways associated with juniper consumption in these goat populations.

RESULTS

The average juniper consumption was 22.4% (± 18.7%) and 7.01% (± 12.1%) in the BS and ANG populations, respectively. The heritability estimates (realized heritability within parenthesis) for juniper consumption were 0.43 ± 0.02 (0.34 ± 0.06) and 0.19 ± 0.03 (0.13 ± 0.03) in BS and ANG, respectively, indicating that juniper consumption can be increased through genetic selection. The repeatability values of predicted juniper consumption were 0.45 for BS and 0.28 for ANG. A total of 571 significant SNP located within or close to 231 genes in BS, and 116 SNP related to 183 genes in ANG were identified based on the genome-wide association analyses. These genes are primarily associated with biological pathways and gene ontology terms related to olfactory receptors, intestinal absorption, and immunity response.

CONCLUSIONS

These findings suggest that juniper consumption is a heritable trait of polygenic inheritance influenced by multiple genes of small effects. The genetic parameters calculated indicate that juniper consumption can be genetically improved in both goat populations.

Developing flexible models for genetic evaluations in smallholder crossbred dairy farms

The productivity of smallholder dairy farms is very low in developing countries. Important genetic gains could be realized using genomic selection, but genetic evaluations need to be tailored for lack of pedigree information and very small farm sizes. To accommodate this situation, we propose a flexible Bayesian model for the genetic evaluation of milk yield, which allows us to simultaneously account for nongenetic random effects for farms and varying SNP variance (BayesR model). First, we used simulations based on real genotype data from Indian crossbred dairy cattle to demonstrate that the proposed model can separate the true genetic and nongenetic parameters even for small farm sizes (2 cows on average) although with high standard errors in scenarios with low heritability. The accuracy of genomic genetic evaluation increased until farm size was approximately 5. We then applied the model to real data from 4,655 crossbred cows with 106,109 monthly test day milk records and 689,750 autosomal SNPs. We estimated a heritability of 0.16 (0.04) for milk yield and using cross-validation, a genomic estimated breeding value (GEBV) accuracy of 0.45 and bias (regression of phenotype on GEBV) of 1.04 (0.26). Estimated genetic parameters were very similar using BayesR, BayesC, and genomic BLUP approaches. Candidate genes near the top variants, IMMP2L and ARHGEF2, have been previously associated with milk protein composition, mastitis resistance, and milk cholesterol content. The estimated heritability and GEBV accuracy for milk yield are much lower than those from intensive or pasture-based systems in many countries. Further increases in the number of phenotyped and genotyped animals in farms with at least 2 cows (preferably 3-5, to allow for dropout of cows) are needed to improve the estimation of genetic effects in these smallholder dairy farms.

Relationship between age, scrotal circumference, postweaning weight and semen quality in Nellore and Caracu bulls: a cross sectional study

The aim of this study was to investigate the relationship between age, scrotal circumference, postweaning weight and semen quality in Nellore and Caracu bulls selected for postweaning weight. Data from the andrological evaluation of 836 bulls born between 2000 and 2019, including 583 Nellore animals (Bos indicus) and 253 Caracu animals (Bos taurus), were used. The bulls were divided into categories of age at the time of assessment: category 1 consisted of animals aged 20 to 23 months (22 ± 0.76 months, 518 ± 94.17 kg), category 2 consisted of animals aged 24 to 35 months (30 ± 4.42 months, 679 ± 137.19 kg), and category 3 consisted of animals ≥ 36 months (60 ± 14.12 months, 907 ± 161.73 kg). The statistical model included the effects of breed, age category, date of semen collection, and breed x age category interaction. Heritability estimates for scrotal circumference at 13 months of age (SC1year) and semen quality traits were obtained for the sample of Nellore animals. Most semen quality traits improved with increasing age in both Nellore and Caracu animals. High heritability was observed for SC1year (0.45), while sperm motility, vigor, turbulence, and major, minor and total sperm defects exhibited low heritability (0.11, 0.019, 0.047, 0.017, 0.017 and 0.019, respectively). Spearman correlations of breeding values for postweaning weight (W378) and SC1year with the semen quality traits were low. Nellore and Caracu bulls have similar semen quality that improves with increasing age. In the Nellore breed, the heritability of SC is high, while semen quality traits exhibit low heritability. Selection for higher postweaning weight does not phenotypically affect the semen quality of bulls at breeding age.

Single-step genomic predictions for heat tolerance of production yields in US Holsteins and Jerseys

The physiological stress caused by excessive heat affects dairy cattle health and production. This study sought to investigate the effect of heat stress on test-day yields in US Holstein and Jersey cows and develop single-step genomic predictions to identify heat tolerant animals. Data included 12.8 million and 2.1 million test-day records, respectively, for 923,026 Holstein and 153,710 Jersey cows in 27 US states. From 2015 through 2021, test-day records from the first 5 lactations included milk, fat, and protein yields (kg). Cow records were included if they had at least 5 test-day records per lactation. Heat stress was quantified by analyzing the effect of a 5-d hourly average temperature-humidity index (THI5d¯) on observed test-day yields. Using a multiple trait repeatability model, a heat threshold (THI threshold) was determined fowr each breed based on the point that the average adjusted yields started to decrease, which was 69 for Holsteins and 72 for Jerseys. An additive genetic component of general production and heat tolerance production were estimated using a multiple trait reaction norm model and single-step genomic BLUP methodology. Random effects were regressed on a function of 5-d hourly average (THI5d¯) and THI threshold. The proportion of test-day records that occurred on or above the respective heat thresholds was 15% for Holstein and 10% for Jersey. Heritability of milk, fat, and protein yields under heat stress for Holsteins increased, with a small standard error, indicating that the additive genetic component for heat tolerance of these traits was observed. This was not as evident in Jersey traits. For Jersey, the permanent environment explained the same or more of the variation in fat and protein yield under heat stress indicating that nongenetic factors may determine heat tolerance for these Jersey traits. Correlations between the general genetic merit of production (in the absence of heat stress) and heat tolerance genetic merit of production traits were moderate in strength and negative. This indicated that selecting for general genetic merit without consideration of heat tolerance genetic merit of production may result in less favorable performance in hot and humid climates. A general genomic estimated breeding value for genetic merit and a heat tolerance genomic estimated breeding value were calculated for each animal. This study contributes to the investigation of the impact of heat stress on US dairy cattle production yields and offers a basis for the implementation of genomic selection. The results indicate that genomic selection for heat tolerance of production yields is possible for US Holsteins and Jerseys, but a study to validate the genomic predictions should be explored.

Combined single-step evaluation of functional longevity of dairy cows including correlated traits

BACKGROUND

For years, multiple trait genetic evaluations have been used to increase the accuracy of estimated breeding values (EBV) using information from correlated traits. In France, accurate approximations of multiple trait evaluations were implemented for traits that are described by different models by combining the results of univariate best linear unbiased prediction (BLUP) evaluations. Functional longevity (FL) is the trait that has most benefited from this approach. Currently, with many single-step (SS) evaluations, only univariate FL evaluations can be run. The aim of this study was to implement a "combined" SS (CSS) evaluation that extends the "combined" BLUP evaluation to obtain more accurate genomic (G) EBV for FL when information from five correlated traits (somatic cell score, clinical mastitis, conception rate for heifers and cows, and udder depth) is added.

RESULTS

GEBV obtained from univariate SS (USS) evaluations and from a CSS evaluation were compared. The correlations between these GEBV showed the benefits of including information from correlated traits. Indeed, a CSS evaluation run without any performances on FL showed that the indirect information from correlated traits to evaluate FL was substantial. USS and CSS evaluations that mimic SS evaluations with data available in 2016 were compared. For each evaluation separately, the GEBV were sorted and then split into 10 consecutive groups (deciles). Survival curves were calculated for each group, based on the observed productive life of these cows as known in 2021. Regardless of their genotyping status, the worst group of heifers based on their GEBV in 2016 was well identified in the CSS evaluation and they had a substantially shorter herd life, while those in the best heifer group had a longer herd life. The gaps between groups were more important for the genotyped than the ungenotyped heifers, which indicates better prediction of future survival.

CONCLUSIONS

A CSS evaluation is an efficient tool to improve FL. It allows a proper combination of information on functional traits that influence culling. In contrast, because of the strong selection intensity on young bulls for functional traits, the benefit of such a "combined" evaluation of functional traits is more modest for these males.

Effects of a mastitis J5 bacterin vaccination on the productive performance of dairy cows: An observational study using propensity score matching techniques

Inferring causal effects between variables when utilizing observational data is challenging due to confounding factors not controlled through a randomized experiment. Propensity score matching can decrease confounding in observational studies and offers insights about potential causal effects of prophylactic management interventions such as vaccinations. The objective of this study was to determine potential causality and impact of vaccination with an Escherichia coli J5 bacterin on the productive performance of dairy cows applying propensity score matching techniques to farm-recorded (e.g., observational) data. Traits of interest included 305-d milk yield (MY305), 305-d fat yield (FY305), 305-d protein yield (PY305), and somatic cell score (SCS). Records from 6,418 lactations generated by 5,121 animals were available for the analysis. Vaccination status of each animal was obtained from producer-recorded information. Confounding variables considered were herd-year-season groups (56 levels), parity (5 levels: 1, 2, 3, 4, and ≥5), and genetic quartile groups (4 levels: top 25% through bottom 25%) derived from genetic predictions for MY305, FY305, PY305, and SCS, as well as for the genetic susceptibility to mastitis. A logistic regression model was applied to estimate the propensity score (PS) for each cow. Subsequently, PS values were used to form pairs of animals (1 vaccinated with 1 unvaccinated control), depending on their PS similarities (difference in PS values of cows within a match required to be <20% of 1 standard deviation of the logit of PS). After the matching process, 2,091 pairs of animals (4,182 records) remained available to infer the causal effects of vaccinating dairy cows with the E. coli J5 bacterin. Causal effects estimation was performed using 2 approaches: simple matching and a bias-corrected matching. According to the PS methodology, causal effects of vaccinating dairy cows with a J5 bacterin on their productive performance were identified for MY305. The simple matched estimator suggested that vaccinated cows produced 163.89 kg more milk over an entire lactation when compared with nonvaccinated counterparts, whereas the bias-corrected estimator suggested that such increment in milk production was of 150.48 kg. Conversely, no causal effects of immunizing dairy cows with a J5 bacterin were identified for FY305, PY305, or SCS. In conclusion, the utilization of PS matching techniques applied to farm-recorded data was feasible and allowed us to identify that vaccination with an E. coli J5 bacterin relates to an overall milk production increment without compromising milk quality.

Appropriate Genetic Approaches for Heat Tolerance and Maintaining Good Productivity in Tropical Poultry Production: A Review

Heat stress is a major environmental threat to poultry production systems, especially in tropical areas. The effects of heat stress have been discovered in several areas, including reduced growth rate, reduced egg production, low feed efficiency, impaired immunological responses, changes in intestinal microflora, metabolic changes, and deterioration of meat quality. Although several methods have been used to address the heat stress problem, it persists. The answer to this problem can be remedied sustainably if genetic improvement approaches are available. Therefore, the purpose of this review article was to present the application of different approaches to genetic improvement in poultry in the hope that users will find suitable solutions for their poultry population and be able to plan future poultry breeding programs.

Effect of reproductive management programs that prioritized artificial insemination at detected estrus or timed artificial insemination on the reproductive performance of primiparous Holstein cows of different genetic merit for fertility

Our objective was to compare reproductive outcomes of primiparous lactating Holstein cows of different genetic merit for fertility submitted for insemination with management programs that prioritized artificial insemination (AI) at detected estrus (AIE) or timed AI (TAI). Moreover, we aimed to determine whether subgroups of cows with different fertility potential would present a distinct response to the reproductive management strategies compared. Lactating primiparous Holstein cows (n = 6 commercial farms) were stratified into high (Hi-Fert), medium (Med-Fert), and low (Lo-Fert) genetic fertility groups (FG) based on a Reproduction Index value calculated from multiple genomic-enhanced predicted transmitting abilities. Within herd and FG, cows were randomly assigned either to a program that prioritized TAI and had an extended voluntary waiting period (P-TAI; n = 1,338) or another that prioritized AIE (P-AIE; n = 1,416) and used TAI for cows, not AIE. Cows in P-TAI received first service by TAI at 84 ± 3 d in milk (DIM) after a Double-Ovsynch protocol, were AIE if detected in estrus after a previous AI, and received TAI after an Ovsynch-56 protocol at 35 ± 3 d after a previous AI if a corpus luteum (CL) was visualized at nonpregnancy diagnosis (NPD) 32 ± 3 d after AI. Cows with no CL visualized at NPD received TAI at 42 ± 3 d after AI after an Ovsynch-56 protocol with progesterone supplementation (P4-Ovsynch). Cows in P-AIE were eligible for AIE after a PGF2α treatment at 53 ± 3 DIM and after a previous AI. Cows not AIE by 74 ± 3 DIM or by NPD 32 ± 3 d after AI received P4-Ovsynch for TAI at 74 ± 3 DIM or 42 ± 3 d after AI. Binary data were analyzed with logistic regression, count data with Poisson regression, continuous data by ANOVA, and time to event data by Cox's proportional hazard regression. Pregnancy per AI (P/AI) to first service was greater for cows in the Hi-Fert (59.8%) than the Med-Fert (53.6%) and Lo-Fert (47.7%) groups, and for the P-TAI (58.7%) than the P-AIE (48.7%) treatment. Overall, P/AI for all second and subsequent AI combined did not differ by treatment (P-TAI = 45.2%; P-AIE = 44.5%) or FG (Hi-Fert = 46.1%; Med-Fert = 46.0%; Lo-Fert = 42.4%). The hazard of pregnancy after calving was greater for the P-AIE than the P-TAI treatment [hazard ratio (HR) = 1.27, 95% CI: 1.17 to 1.37)], and for the Hi-Fert than the Med-Fert (HR = 1.16, 95% CI: 1.05 to 1.28) and Lo-Fert (HR = 1.34, 95% CI: 1.20 to 1.49) groups. More cows in the Hi-Fert (91.2%) than the Med-Fert (88.4%) and Lo-Fert (85.8%) groups were pregnant at 200 DIM. Within FG, the hazard of pregnancy was greater for the P-AIE than the P-TAI treatment for the Hi-Fert (HR = 1.41, 95% CI: 1.22 to 1.64) and Med-Fert (HR = 1.28, 95% CI: 1.12 to 1.46) groups but not for the Lo-Fert group (HR = 1.13, 95% CI: 0.98 to 1.31). We conclude that primiparous Holstein cows of superior genetic merit for fertility had better reproductive performance than cows of inferior genetic merit for fertility, regardless of the type of reproductive management used. In addition, the effect of programs that prioritized AIE or TAI on reproductive performance for cows of superior or inferior genetic merit for fertility depended on the outcomes evaluated. Thus, programs that prioritize AIE or TAI could be used to affect certain outcomes of reproductive performance or management.

Accelerated matrix-vector multiplications for matrices involving genotype covariates with applications in genomic prediction

In the last decade, a number of methods have been suggested to deal with large amounts of genetic data in genomic predictions. Yet, steadily growing population sizes and the suboptimal use of computational resources are pushing the practical application of these approaches to their limits. As an extension to the C/CUDA library miraculix, we have developed tailored solutions for the computation of genotype matrix multiplications which is a critical bottleneck in the empirical evaluation of many statistical models. We demonstrate the benefits of our solutions at the example of single-step models which make repeated use of this kind of multiplication. Targeting modern Nvidia® GPUs as well as a broad range of CPU architectures, our implementation significantly reduces the time required for the estimation of breeding values in large population sizes. miraculix is released under the Apache 2.0 license and is freely available at https://github.com/alexfreudenberg/miraculix.

Comparison of accuracy of breeding value for cow from three methods in Hanwoo (Korean cattle) population

In Korea, Korea Proven Bulls (KPN) program has been well-developed. Breeding and evaluation of cows are also an essential factor to increase earnings and genetic gain. This study aimed to evaluate the accuracy of cow breeding value by using three methods (pedigree index [PI], pedigree-based best linear unbiased prediction [PBLUP], and genomic-BLUP [GBLUP]). The reference population (n = 16,971) was used to estimate breeding values for 481 females as a test population. The accuracy of GBLUP was 0.63, 0.66, 0.62 and 0.63 for carcass weight (CWT), eye muscle area (EMA), back-fat thickness (BFT), and marbling score (MS), respectively. As for the PBLUP method, accuracy of prediction was 0.43 for CWT, 0.45 for EMA, 0.43 for MS, and 0.44 for BFT. Accuracy of PI method was the lowest (0.28 to 0.29 for carcass traits). The increase by approximate 20% in accuracy of GBLUP method than other methods could be because genomic information may explain Mendelian sampling error that pedigree information cannot detect. Bias can cause reducing accuracy of estimated breeding value (EBV) for selected animals. Regression coefficient between true breeding value (TBV) and GBLUP EBV, PBLUP EBV, and PI EBV were 0.78, 0.625, and 0.35, respectively for CWT. This showed that genomic EBV (GEBV) is less biased than PBLUP and PI EBV in this study. In addition, number of effective chromosome segments (Me) statistic that indicates the independent loci is one of the important factors affecting the accuracy of BLUP. The correlation between Me and the accuracy of GBLUP is related to the genetic relationship between reference and test population. The correlations between Me and accuracy were -0.74 in CWT, -0.75 in EMA, -0.73 in MS, and -0.75 in BF, which were strongly negative. These results proved that the estimation of genetic ability using genomic data is the most effective, and the smaller the Me, the higher the accuracy of EBV.

An effective hyper-parameter can increase the prediction accuracy in a single-step genetic evaluation

The H-matrix best linear unbiased prediction (HBLUP) method has been widely used in livestock breeding programs. It can integrate all information, including pedigree, genotypes, and phenotypes on both genotyped and non-genotyped individuals into one single evaluation that can provide reliable predictions of breeding values. The existing HBLUP method requires hyper-parameters that should be adequately optimised as otherwise the genomic prediction accuracy may decrease. In this study, we assess the performance of HBLUP using various hyper-parameters such as blending, tuning, and scale factor in simulated and real data on Hanwoo cattle. In both simulated and cattle data, we show that blending is not necessary, indicating that the prediction accuracy decreases when using a blending hyper-parameter <1. The tuning process (adjusting genomic relationships accounting for base allele frequencies) improves prediction accuracy in the simulated data, confirming previous studies, although the improvement is not statistically significant in the Hanwoo cattle data. We also demonstrate that a scale factor, α, which determines the relationship between allele frequency and per-allele effect size, can improve the HBLUP accuracy in both simulated and real data. Our findings suggest that an optimal scale factor should be considered to increase prediction accuracy, in addition to blending and tuning processes, when using HBLUP.

Single-Step Genomic Prediction of Superovulatory Response Traits in Japanese Black Donor Cows

We assessed the performance of single-step genomic prediction of breeding values for superovulatory response traits in Japanese Black donor cows. A total of 25,332 records of the total number of embryos and oocytes (TNE) and the number of good embryos (NGE) per flush for 1874 Japanese Black donor cows were collected during 2008 and 2022. Genotype information on 36,426 autosomal single-nucleotide polymorphisms (SNPs) for 575 out of the 1,874 cows was used. Breeding values were predicted exploiting a two-trait repeatability animal model. Two genetic relationship matrices were used, one based on pedigree information (A matrix) and the other considering both pedigree and SNP marker genotype information (H matrix). Estimated heritabilities of TNE and NGE were 0.18 and 0.11, respectively, when using the H matrix, which were both slightly lower than when using the A matrix (0.26 for TNE and 0.16 for NGE). Estimated genetic correlations between the traits were 0.61 and 0.66 when using H and A matrices, respectively. When the variance components were the same in breeding value prediction, the mean reliability was greater when using the H matrix than when using the A matrix. This advantage seems more prominent for cows with low reliability when using the A matrix. The results imply that introducing single-step genomic prediction could boost the rate of genetic improvement of superovulatory response traits, but efforts should be made to maintain genetic diversity when performing selection.

Selective genotyping to implement genomic selection in beef cattle breeding

Genomic selection (GS) plays an essential role in livestock genetic improvement programs. In dairy cattle, the method is already a recognized tool to estimate the breeding values of young animals and reduce generation intervals. Due to the different breeding structures of beef cattle, the implementation of GS is still a challenge and has been adopted to a much lesser extent than dairy cattle. This study aimed to evaluate genotyping strategies in terms of prediction accuracy as the first step in the implementation of GS in beef while some restrictions were assumed for the availability of phenotypic and genomic information. For this purpose, a multi-breed population of beef cattle was simulated by imitating the practical system of beef cattle genetic evaluation. Four genotyping scenarios were compared to traditional pedigree-based evaluation. Results showed an improvement in prediction accuracy, albeit a limited number of animals being genotyped (i.e., 3% of total animals in genetic evaluation). The comparison of genotyping scenarios revealed that selective genotyping should be on animals from both ancestral and younger generations. In addition, as genetic evaluation in practice covers traits that are expressed in either sex, it is recommended that genotyping covers animals from both sexes.

Comparing Methods to Select Candidates for Re-Genotyping to Impute Higher-Density Genotype Data in a Japanese Black Cattle Population: A Case Study

As optimization methods to identify the best animals for dense genotyping to construct a reference population for genotype imputation, the MCA and MCG methods, which use the pedigree-based additive genetic relationship matrix (A matrix) and the genomic relationship matrix (G matrix), respectively, have been proposed. We assessed the performance of MCA and MCG methods using 575 Japanese Black cows. Pedigree data were provided to trace back up to five generations to construct the A matrix with changing the pedigree depth from 1 to 5 (five MCA methods). Genotype information on 36,426 single-nucleotide polymorphisms was used to calculate the G matrix based on VanRaden's methods 1 and 2 (two MCG methods). The MCG always selected one cow per iteration, while MCA sometimes selected multiple cows. The number of commonly selected cows between the MCA and MCG methods was generally lower than that between different MCA methods or between different MCG methods. For the studied population, MCG appeared to be more reasonable than MCA in selecting cows as a reference population for higher-density genotype imputation to perform genomic prediction and a genome-wide association study.

Computer vision system using deep learning to predict rib and loin yield in the fish Colossoma macropomum

Computer vision system (CVSs) are effective tools that enable large-scale phenotyping with a low-cost and non-invasive method, which avoids animal stress. Economically important traits, such as rib and loin yield, are difficult to measure; therefore, the use of CVS is crucial to accurately predict several measures to allow their inclusion in breeding goals by indirect predictors. Therefore, this study aimed (1) to validate CVS by a deep learning approach and to automatically predict morphometric measurements in tambaqui and (2) to estimate genetic parameters for growth traits and body yield. Data from 365 individuals belonging to 11 full-sib families were evaluated. Seven growth traits were measured. After biometrics, each fish was processed in the following body regions: head, rib, loin, R + L (rib + loin). For deep learning image segmentation, we adopted a method based on the instance segmentation of the Mask R-CNN (Region-based Convolutional Neural Networks) model. Pearson's correlation values between measurements predicted manually and automatically by the CVS were high and positive. Regarding the classification performance, visible differences were detected in only about 3% of the images. Heritability estimates for growth and body yield traits ranged from low to high. The genetic correlations between the percentage of body parts and morphometric characteristics were favorable and highly correlated, except for percentage head, whose correlations were unfavorable. In conclusion, the CVS validated in this image dataset proved to be resilient and can be used for large-scale phenotyping in tambaqui. The weight of the rib and loin are traits under moderate genetic control and should respond to selection. In addition, standard length and pelvis length can be used as an efficient and indirect selection criterion for body yield in this tambaqui population.

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.

Genomic prediction in Brazilian Guzerá cattle: application of a single-step approach to productive and reproductive traits

This study aimed to investigate the feasibility of genomic prediction for productive and reproductive traits in Guzerá cattle using single-step genomic best linear unbiased prediction (ssGBLUP). Evaluations included the 305-day cumulative yields (first lactation, in kg) of milk, lactose, protein, fat, and total solids; adjusted body weight (kg) at the ages of 450, 365, and 210 days; and age at first calving (in days), from a database containing 197,283 measurements from Guzerá males and females born between 1954 and 2018. The pedigree included 433,823 animals spanning up to 14 overlapping generations. A total of 1618 animals were genotyped. The analyses were performed using ssGBLUP and traditional BLUP methods. Predictive ability and bias were accessed using cross-validation: predictive ability was similar between the methods and ranged from 0.27 to 0.47 for the genomic-based model and from 0.30 to 0.45 for the pedigree-based model; the bias was also similar between the methods, ranging from 0.88 to 1.35 in the genomic-based model and from 0.96 to 1.41 in the pedigree-based model. The individual accuracies of breeding values were evidently increased in the genomic evaluation, with values ranging from 0.41 to 0.56 in the genomic-based model and from 0.26 to 0.54 in the pedigree-based model. Even based on a small number of genotyped animals and a small database for some traits, the results suggest that ssGBLUP is feasible and may be applied to national genetic evaluation of the breed to increase the accuracy of breeding values without greatly impacting predictive ability and bias.

Unraveling the phenotypic and genomic background of behavioral plasticity and temperament in North American Angus cattle

BACKGROUND

Longitudinal records of temperament can be used for assessing behavioral plasticity, such as aptness to learn, memorize, or change behavioral responses based on affective state. In this study, we evaluated the phenotypic and genomic background of North American Angus cow temperament measured throughout their lifetime around the weaning season, including the development of a new indicator trait termed docility-based learning and behavioral plasticity. The analyses included 273,695 and 153,898 records for yearling (YT) and cow at weaning (CT) temperament, respectively, 723,248 animals in the pedigree, and 8784 genotyped animals. Both YT and CT were measured when the animal was loading into/exiting the chute. Moreover, CT was measured around the time in which the cow was separated from her calf. A random regression model fitting a first-order Legendre orthogonal polynomial was used to model the covariance structure of temperament and to assess the learning and behavioral plasticity (i.e., slope of the regression) of individual cows. This study provides, for the first time, a longitudinal perspective of the genetic and genomic mechanisms underlying temperament, learning, and behavioral plasticity in beef cattle.

RESULTS

CT measured across years is heritable (0.38-0.53). Positive and strong genetic correlations (0.91-1.00) were observed among all CT age-group pairs and between CT and YT (0.84). Over 90% of the candidate genes identified overlapped among CT age-groups and the estimated effect of genomic markers located within important candidate genes changed over time. A small but significant genetic component was observed for learning and behavioral plasticity (heritability = 0.02 ± 0.002). Various candidate genes were identified, revealing the polygenic nature of the traits evaluated. The pathways and candidate genes identified are associated with steroid and glucocorticoid hormones, development delay, cognitive development, and behavioral changes in cattle and other species.

CONCLUSIONS

Cow temperament is highly heritable and repeatable. The changes in temperament can be genetically improved by selecting animals with favorable learning and behavioral plasticity (i.e., habituation). Furthermore, the environment explains a large part of the variation in learning and behavioral plasticity, leading to opportunities to also improve the overall temperament by refining management practices. Moreover, behavioral plasticity offers opportunities to improve the long-term animal and handler welfare through habituation.

Derivation of indirect predictions using genomic recursions across generations in a broiler population

Genomic estimated breeding values (GEBV) of animals without phenotypes can be indirectly predicted using recursions on GEBV of a subset. To maximize predictive ability of indirect predictions (IP), the subset must represent the independent chromosome segments segregating in the population. We aimed to 1) determine the number of animals needed in recursions to maximize predictive ability, 2) evaluate equivalency IP-GEBV, and 3) investigate trends in predictive ability of IP derived from recent vs. distant generations or accumulating phenotypes from recent to past generations. Data comprised pedigree of 825K birds hatched over 12 overlapping generations, phenotypes for body weight (BW; 820K), residual feed intake (RF; 200K) and weight gain during a trial period (WG; 200K), and breast meat percent (BP; 43K). A total of 154K birds (last six generations) had genotypes. The number of animals that maximize predictive ability was assessed based on the number of largest eigenvalues explaining 99% of variation in the genomic relationship matrix (1Me = 7,131), twice (2Me), or a fraction of this number (i.e., 0.75, 0.50, or 0.25Me). Equivalency between IP and GEBV was measured by correlating these two sets of predictions. GEBV were obtained as if generation 12 (validation animals) was part of the evaluation. IP were derived from GEBV of animals from generations 8 to 11 or generations 11, 10, 9, or 8. IP predictive ability was defined as the correlation between IP and adjusted phenotypes. The IP predictive ability increased from 0.25Me to 1Me (11%, on average); the change from 1Me to 2Me was negligible (0.6%). The correlation IP-GEBV was the same when IP were derived from a subset of 1Me animals chosen randomly across generations (8 to 11) or from generation 11 (0.98 for BW, 0.99 for RF, WG, and BP). A marginal decline in the correlation was observed when IP were based on GEBV of animals from generation 8 (0.95 for BW, 0.98 for RF, WG, and BP). Predictive ability had a similar trend; from generation 11 to 8, it changed from 0.32 to 0.31 for BW, from 0.39 to 0.38 for BP, and was constant at 0.33(0.22) for RF(WG). Predictive ability had a slight to moderate increase accumulating up to four generations of phenotypes. 1Me animals provide accurate IP, equivalent to GEBV. A minimum decay in predictive ability is observed when IP are derived from GEBV of animals from four generations back, possibly because of strong selection or the model not being completely additive.

An efficient genomic prediction method without the direct inverse of the genomic relationship matrix

GBLUP, the most widely used genomic prediction (GP) method, consumes large and increasing amounts of computational resources as the training population size increases due to the inverse of the genomic relationship matrix (GRM). Therefore, in this study, we developed a new genomic prediction method (RHEPCG) that avoids the direct inverse of the GRM by combining randomized Haseman-Elston (HE) regression (RHE-reg) and a preconditioned conjugate gradient (PCG). The simulation results demonstrate that RHEPCG, in most cases, not only achieves similar predictive accuracy with GBLUP but also significantly reduces computational time. As for the real data, RHEPCG shows similar or better predictive accuracy for seven traits of the Arabidopsis thaliana F2 population and four traits of the Sorghum bicolor RIL population compared with GBLUP. This indicates that RHEPCG is a practical alternative to GBLUP and has better computational efficiency.

Candidate genes for marker-assisted selection for growth, carcass and meat quality traits in rabbits

Growth and meat production are the most relevant targets for animal breeders, there are strong relationships between animal growth regulation, body composition and meat quality. Therefore, it is essential to identify the genetic factors that are controlling growth, carcass, and meat quality traits and to explore the correlations between identified genes of those traits. Identification of candidate genes may shift rabbit breeding from classical to modern approaches, which offer great potential to accelerate genetic improvement plans, especially in developing countries. The current work reviews several genes and mutations affecting growth, carcass and meat quality traits. These candidate genes and mutations can be incorporated into MAS programs to improve rabbit breeds especially local breeds, provided that a reasonable proportion of trait additive genetic variance is explained by the significant marker. Furthermore, we highlighted the indispensable need for more researches investigating candidate genes for different traits.

The prediction accuracies of linear-type traits in Czech Holstein cattle when using ssGBLUP or wssGBLUP

The aim of this study was to assess the contribution of the weighted single-step genomic best linear unbiased prediction (wssGBLUP) method compared to the single-step genomic best linear unbiased prediction (ssGBLUP) method for genomic evaluation of 25 linear-type traits in the Czech Holstein cattle population. The nationwide database of linear-type traits with 6,99,681 records combined with deregressed proofs from Interbull (MACE method) was used as the input data. Genomic breeding values (GEBVs) were predicted based on these phenotypes using ssGBLUP and wssGBLUP methods using the BLUPF90 software. The bull validation test was employed which was based on comparing GEBVs of young bulls (N = 334) with no progeny in 2016. A minimum of 50 daughters with their own performance in 2020 was chosen to verify the contribution to the GEBV prediction, GEBV reliability, validation reliabilities (R2), and regression coefficients (b1). The results showed that the differences between the two methods were negligible. The low benefit of wssGBLUP may be due to the inclusion of a small number of SNPs; therefore, most predictions rely on polygenic relationships between animals. Nevertheless, the benefits of wssGBLUP analysis should be assessed with respect to specific population structures and given traits.

A Note on the Conditioning of the -1 Matrix Used in Single-Step GBLUP

The single-step genomic BLUP (ssGBLUP) is used worldwide for the simultaneous genetic evaluation of genotyped and non-genotyped animals. It is easily extendible to all BLUP models by replacing the pedigree-based additive genetic relationship matrix (A) with an augmented pedigree-genomic relationship matrix (H). Theoretically, H does not introduce any artificially inflated variance. However, inflated genetic variances have been observed due to the incomparability between the genomic relationship matrix (G) and A used in H. Usually, G is blended and tuned with A22 (the block of A for genotyped animals) to improve its numerical condition and compatibility. If deflation/inflation is still needed, a common approach is weighting G-1-A22-1 in the form of τG-1-ωA22-1, added to A-1 to form H-1. In some situations, this can violate the conditional properties upon which H is built. Different ways of weighting the H-1 components (A-1, G-1, A22-1, and H-1 itself) were studied to avoid/minimise the violations of the conditional properties of H. Data were simulated on ten populations and twenty generations. Responses to weighting different components of H-1 were measured in terms of the regression of phenotypes on the estimated breeding values (the lower the slope, the higher the inflation) and the correlation between phenotypes and the estimated breeding values (predictive ability). Increasing the weight on H-1 increased the inflation. The responses to weighting G-1 were similar to those for H-1. Increasing the weight on A-1 (together with A22-1) was not influential and slightly increased the inflation. Predictive ability is a direct function of the slope of the regression line and followed similar trends. Responses to weighting G-1-A22-1 depend on the inflation/deflation of evaluations from A-1 to H-1 and the compatibility of the two matrices with the heritability used in the model. One possibility is a combination of weighting G-1-A22-1 and weighting H-1. Given recent advances in ssGBLUP, conditioning H-1 might become an interim solution from the past and then not be needed in the future.