A relationship matrix including full pedigree and genomic information

Genome-wide association study for stayability at different calvings in Nellore beef cattle

BACKGROUNDING

Stayability, which may be defined as the probability of a cow remaining in the herd until a reference age or at a specific number of calvings, is usually measured late in the animal's life. Thus, if used as selection criteria, it will increase the generation interval and consequently might decrease the annual genetic gain. Measuring stayability at an earlier age could be a reasonable strategy to avoid this problem. In this sense, a better understanding of the genetic architecture of this trait at different ages and/or at different calvings is important. This study was conducted to identify possible regions with major effects on stayability measured considering different numbers of calvings in Nellore cattle as well as pathways that can be involved in its expression throughout the female's productive life.

RESULTS

The top 10 most important SNP windows explained, on average, 17.60% of the genetic additive variance for stayability, varying between 13.70% (at the eighth calving) and 21% (at the fifth calving). These SNP windows were located on 17 chromosomes (1, 2, 4, 6, 7, 8, 9, 10, 11, 12, 13, 14, 18, 19, 20, 27, and 28), and they harbored a total of 176 annotated genes. The functional analyses of these genes, in general, indicate that the expression of stayability from the second to the sixth calving is mainly affected by genetic factors related to reproductive performance, and nervous and immune systems. At the seventh and eighth calvings, genes and pathways related to animal health, such as density bone and cancer, might be more relevant.

CONCLUSION

Our results indicate that part of the target genomic regions in selecting for stayability at earlier ages (from the 2th to the 6th calving) would be different than selecting for this trait at later ages (7th and 8th calvings). While the expression of stayability at earlier ages appeared to be more influenced by genetic factors linked to reproductive performance together with an overall health/immunity, at later ages genetic factors related to an overall animal health gain relevance. These results support that selecting for stayability at earlier ages (perhaps at the second calving) could be applied, having practical implications in breeding programs since it could drastically reduce the generation interval, accelerating the genetic progress.

Validation of reaction norm breeding values for robustness in Australian sheep

BACKGROUND

There can be variation between animals in how stable their genetic merit is across different environments due to genotype-by-environment (G×E) interactions. This variation could be used in breeding programs to select robust genotypes that combine high overall performance with stable genetic ranking across environments. There have been few attempts to validate breeding values for robustness in livestock, although this is a necessary step towards their implementation in selection decisions. The objective of this study was to validate breeding values for the robustness of body weight across different growth environments that were estimated using reaction norm models in sheep data.

RESULTS

Using threefold cross-validation for the progeny of 337 sires, the average correlation between single-step breeding values for the reaction norm slope and the realised robustness of progeny across different growth environments was 0.21. The correlation between breeding values for the reaction slope estimated independently in two different datasets linked by common sires was close to the expected correlation based on theory.

CONCLUSIONS

Slope estimated breeding values (EBV) obtained using reaction norm models were predictive of the phenotypic robustness of progeny across different environments and were consistent for sires with progeny in two different datasets. Selection based on reaction norm EBV could be used to increase the robustness of a population to environmental variation.

Genetic evaluation of productive longevity in a multibreed beef cattle population

Genetic selection for traits that have a direct impact on profitability, such as productive longevity (PL), which blends cow longevity with regular reproductive performance, is fundamental for the economic success of beef cow-calf operations. The purpose of this study was to develop a data screening strategy and a statistical model to predict genetic merit for PL in a multibreed beef cattle population. Pedigree (n = 1,352,765) and phenotype (n = 978,382) information were provided by Leachman Cattle of Colorado, and genotypes (n = 26,342) were provided by the Zoetis commercial genotyping laboratory. A repeatability model (REP) including the systematic effects of age at first calving, year-season of progeny birth, pedigree-based retained heterosis, and parity number, as well as the random effects of the additive genetic, permanent environment, contemporary group, and residual were fitted to adjust PL. In addition, a random regression model (RRM) was fitted to investigate PL considering the same effects, with the difference that random effects were regressed on parity. Estimated breeding values (EBV) were obtained by single-step GBLUP (ssGBLUP) and transformed to predict differences in the number of calves through linear regression. Predictive performance was assessed in a group of 7,268 cows born in 2010. Heritability estimates for PL were relatively low, with values of 0.109 for REP and a decreasing trend for RRM with values ranging from 0.16 to 0.04. Repeatability for PL was of moderate magnitude, with values of 0.415 for REP and from 0.29 to 0.57 for RRM. Heritability estimates suggest that most of the phenotypic variation was accounted for by environmental factors, but long-term genetic selection could still be effective. REP was more efficient than RRM, showing the lower number of iterations and time to reach convergence with comparable solutions to RRM. Validation results showed that correlations between EBV and phenotypes (observed/precorrected) increased over the years ranging from 0.04 to 0.92. Repeatability values and the validation approach suggested that using a cow's first record (second parity success or failure) is a reasonably good indicator of posterior performance for PL. Therefore, the inclusion of PL in a multibreed genetic evaluation program and incorporation into selection indexes with existing economic traits can enable more profitable selection and breeding decisions in beef cattle herds.

Quantifying phenotypic and genetic variation for cow fertility phenotypes in American Simmental using total herd reporting data

Reproduction plays a major role in the production efficiency of livestock species. However, cow-centric reproductive traits tend to be lowly heritable and are not expressed until later in an animal's lifetime, making phenotypic selection alone inefficient at generating genetic gain. Genetic progress can be accelerated by focusing selection on the predicted genetic component of reproductive traits using Expected Progeny Differences. We used the American Simmental Association's performance and Total Herd Enrollment data, made up of 533,155 calving records from 303,158 females (132,403 cows and 170,755 heifers), 33,732 of which are genotyped, to explore three continuous and two discrete phenotypes focused on quantifying early and sustained fertility in beef cows. We analyzed calving date (CD) (cow's CD relative to the start of the calving season), calving interval (CI) (days between calves), first calving interval (FCI) (CI observation between the first and second calving record for a female), heifer pregnancy (HP) (did the animal calve as 2-yr-old), and discrete early calving (DEC) (did animal calve in the first 30 d of the calving season) as distinct, but correlated measures of fertility. This dataset provides insight into population-wide trends related to cow attrition, calving season lengths, and phenotypic variation in fertility. We used pedigree and genomic REML to estimate these six phenotypes' genetic, permanent environment, and temporary environmental variance components. Pedigree-estimated heritabilities were 0.06 (± 0.000011) for CD, (0.04 ± 0.000005) for CI, 0.07 (± 0.000016) for DEC, 0.05 ± 0.000041 for FCI, and 0.23 (± 0.000099) for HP, consistent with other fertility traits across beef and dairy cattle. The incorporation of genomics increased the heritability estimate for HP (0.24 ± 0.000098) and decreased the estimate for FCI (0.04 ± 0.000029). Positive phenotypic and genetic correlations were found among these phenotypes (rG = 0.01 to 0.96). These results call for further work in optimizing genetic predictions and exploration of the genetic architecture through genome-wide studies. Whole herd reporting date frameworks represent opportunities for measuring new reproductive phenotypes, but their utility in genetic evaluations will rely on novel trait initiatives and consistent recording that captures more detailed data.

Combining large broiler populations into a single genomic evaluation: dealing with genetic divergence

Combining breeding populations that have diverged at some point is a conventional practice, particularly in the poultry industry, where generation intervals are short and genetic evaluations should be frequently available. This study aimed to assess the feasibility of combining large, distantly genetically connected broiler populations into a single genomic evaluation within the single-step GBLUP framework. The pedigree data for broiler lines 1 and 2 consisted of 428,790 and 477,488 animals, being 156,088 and 186,387 genotyped, respectively. Phenotypic data for body weight (kg), carcass yield (%), mortality (1 to 2), and feet health (1 to 7) were collected for 397,974 animals in line 1 and 458,881 in line 2. A 4-trait model was employed for the analyses, and genetic differences between the populations were addressed through different approaches: introducing an additional fixed effect accounting for the line of origin (M2) or making each fixed effect origin-specific (M3). Those models were compared against a conventional model (M1) that did not account for animal origin in the evaluation. Unknown parent groups (UPG) and Metafounders (MF) were fit to account for the genetic differences in M1, M2, and M3; they were set based on the animal's line of origin and sex. Accuracy, bias, and dispersion were used to assess the performances of the models using the Linear Regression method. Validations were performed separately within individual lines and collectively after combining the 2 lines to better assess the advantages of combining the 2 populations. Overall, the accuracy increased when the 2 populations were combined compared to the accuracies obtained from evaluating each line individually. Notably, there were no apparent differences among the models regarding accuracy and dispersion. Regarding bias, using models M2 or M3 with UPG yields the least biased estimates in the combined evaluation. Thus, when combining different populations into a single genomic evaluation, accounting for the genetic and non-genetic differences among the lines ensures accurate and less biased predictions.

Genomic evaluation for breeding and genetic management in Cordia africana, a multipurpose tropical tree species

BACKGROUND

Planting tested forest reproductive material is crucial to ensure the increased resilience of intensively managed productive stands for timber and wood product markets under climate change scenarios. Single-step Genomic Best Linear Unbiased Prediction (ssGBLUP) analysis is a cost-effective option for using genomic tools to enhance the accuracy of predicted breeding values and genetic parameter estimation in forest tree species. Here, we tested the efficiency of ssGBLUP in a tropical multipurpose tree species, Cordia africana, by partial population genotyping. A total of 8070 trees from three breeding seedling orchards (BSOs) were phenotyped for height. We genotyped 6.1% of the phenotyped individuals with 4373 single nucleotide polymorphisms. The results of ssGBLUP were compared with pedigree-based best linear unbiased prediction (ABLUP) and genomic best linear unbiased prediction (GBLUP), based on genetic parameters, theoretical accuracy of breeding values, selection candidate ranking, genetic gain, and predictive accuracy and prediction bias.

RESULTS

Genotyping a subset of the study population provided insights into the level of relatedness in BSOs, allowing better genetic management. Due to the inbreeding detected within the genotyped provenances, we estimated genetic parameters both with and without accounting for inbreeding. The ssGBLUP model showed improved performance in terms of additive genetic variance and theoretical breeding value accuracy. Similarly, ssGBLUP showed improved predictive accuracy and lower bias than the pedigree-based relationship matrix (ABLUP).

CONCLUSIONS

This study of C. africana, a species in decline due to deforestation and selective logging, revealed inbreeding depression. The provenance exhibiting the highest level of inbreeding had the poorest overall performance. The use of different relationship matrices and accounting for inbreeding did not substantially affect the ranking of candidate individuals. This is the first study of this approach in a tropical multipurpose tree species, and the analysed BSOs represent the primary effort to breed C. africana.

Construction of a strawberry breeding core collection to capture and exploit genetic variation

BACKGROUND

Genetic diversity is crucial for the success of plant breeding programs and core collections are important resources to capture this diversity. Many core collections have already been constructed by gene banks, whose main goal is to obtain a panel of a limited number of genotypes to simplify management practices and to improve shareability while retaining as much diversity as possible. However, as gene banks have a different composition and goal than plant breeding programs, constructing a core collection for a plant breeding program should consider different aspects.

RESULTS

In this study, we present a novel approach for constructing a core collection by integrating both genomic and pedigree information to maximize the representation of the breeding germplasm in a minimum subset of genotypes while accounting for future genetic variation within a strawberry breeding program. Our stepwise approach starts with selecting the most important crossing parents of advanced selections and genotypes included for specific traits, to represent also future genetic variation. We then use pedigree-genomic-based relationship coefficients combined with the 'accession to nearest entry' criterion to complement the core collection and maximize its representativeness of the current breeding program. Combined pedigree-genomic-based relationship coefficients allow for accurate relationship estimation without the need to genotype every individual in the breeding program.

CONCLUSIONS

This stepwise construction of a core collection in a strawberry breeding program can be applied in other plant breeding programs to construct core collections for various purposes.

Genetic analysis of prolificacy, milk production and composition traits in Egyptian Zaraibi goats

This study was conducted to genetically and environmentally characterize prolificacy (litter size and weight at birth; LSB and LWB and litter size and weight at weaning; LSW and LWW, respectively), milk yield at the 7th (MY7), 15th (MY15), 30th (MY30), 60th (MY60), 90th (MY90) day of lactation and monthly milk yield (MMY) and milk composition traits in Egyptian Zaraibi goats. A total of 443 and 421 records produced by 121 Zaraibi lactating goats were used to assess prolificacy and milk production traits, respectively. The milk composition traits were measured using 371 milk samples obtained at random from 53 goats. The fourth parity showed the highest values for LWB, LWW, and MMY (3.62, 18.15, and 28.99 kg, respectively). Milk composition traits revealed an inverse tendency, decreasing until the second month and then increasing until the seventh month. The heritability estimates ranged from 0.07 to 0.13, from 0.04 to 0.39, and from 0.07 to 0.33 for prolificacy, milk yield, and milk composition traits, respectively. Negatively high genetic correlations between MMY and all milk composition traits were found. MMY had the highest estimate of heritability (0.39 ± 0.07), this means that the genetic improvement of this trait could be achieved through direct selection.

Additive and non-additive genetic variance in juvenile Sitka spruce (Picea sitchensis Bong. Carr)

Many quantitative genetic models assume that all genetic variation is additive because of a lack of data with sufficient structure and quality to determine the relative contribution of additive and non-additive variation. Here the fractions of additive (fa) and non-additive (fd) genetic variation were estimated in Sitka spruce for height, bud burst and pilodyn penetration depth. Approximately 1500 offspring were produced in each of three sib families and clonally replicated across three geographically diverse sites. Genotypes from 1525 offspring from all three families were obtained by RADseq, followed by imputation using 1630 loci segregating in all families and mapped using the newly developed linkage map of Sitka spruce. The analyses employed a new approach for estimating fa and fd, which combined all available genotypic and phenotypic data with spatial modelling for each trait and site. The consensus estimate for fa increased with age for height from 0.58 at 2 years to 0.75 at 11 years, with only small overlap in 95% support intervals (I95). The estimated fa for bud burst was 0.83 (I95=[0.78, 0.90]) and 0.84 (I95=[0.77, 0.92]) for pilodyn depth. Overall, there was no evidence of family heterogeneity for height or bud burst, or site heterogeneity for pilodyn depth, and no evidence of inbreeding depression associated with genomic homozygosity, expected if dominance variance was the major component of non-additive variance. The results offer no support for the development of sublines for crossing within the species. The models give new opportunities to assess more accurately the scale of non-additive variation.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11295-023-01627-5.

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.

Investigation of potential genetic factors for growth traits in yellow-feather broilers using weighted single-step genome-wide association study

Yellow-feather broilers take a large portion of poultry industry in China due to its meat characteristics. Improving the growth traits of yellow-feathered broilers will have great significance for the Chinese poultry market. The current study was designed to investigate the potential genetic factors using the weighted single-step genome-wide association study (wssGWAS) method, which takes consideration of more factors including pedigree, sex, environment and has more accuracy than traditional GWAS. The yellow-feather dwarf chickens from Wens Nanfang Poultry Breeding Co. Ltd. were revolved to recode 9 growth traits: Average daily gain (ADG), body weight (BW) at 45 d, 49 d, 56 d, 63 d, 70 d, 77 d, 84 d, 91 d for analysis. For the results, the region 4.63 to 5.03 Mb on chromosome 15, which was the QTL overlapped in BW45, BW49, BW56, BW63, BW84, might be the crucial genetic region for growth traits. Seven GO terms and 3 KEGG pathways, GO:0005200, GO:0005882, GO:0045111, GO:0099513, GO:0099081, GO:0099512, GO:0099080, KEGG:gga04020, KEGG:gga04540, KEGG:gga04210, were detected to relevant with growth traits. The genes enriched in these biological processes (NRAS, TUBB1, ADORA2B, NTRK3, NGF, TNNC2, F-KER, LOC429492, LOC431325, LOC431324, LOC396480) might have the function in growth of yellow-feather broilers.

Effects of prenatal heat stress on birth weight and birth weight genetic parameters in German Holstein calves

The aim of this study was to infer the effects of heat stress (HS) during late gestation of dams on phenotypes and on direct and maternal genetic parameters for birth weight (BiW). We considered 171,221 Holstein calves kept in 56 large-scale co-operator herds. For a clear separation of maternal effects, only calves from dams with at least 3 offspring were included in the analyses. The genotype data set comprised 41,143 SNPs from 1,883 Holstein bulls. Temperature-humidity indices (THI) during the last 8 wk of gestation were calculated in each herd to reflect prenatal HS. A further prenatal HS descriptor was the first principal component (PC1) from principal component analysis considering the daily THI during the last 56 d of gestation. Regression coefficients of BiW on prenatal THI during the last 12 wk of gestation and PC1 were estimated in 13 consecutive phenotypic analyses. The strongest BiW decline was -0.63 kg per standardized THI, identified during 50 to 56 d before birth. A reaction norm model with weekly prenatal THI or PC1 nested within maternal genetic and maternal permanent environmental effects was defined to infer maternal sensitivity in response to prenatal THI alterations. Direct BiW heritabilities were close to 0.33 in the course of prenatal THI. Maternal BiW heritabilities marginally increased from 0.07 to 0.08 with increasing THI. Genetic correlations between maternal genetic effects at maximum HS levels and remaining THI were larger than 0.95, indicating the absence of genotype by time-lagged HS interactions. In contrast, maternal permanent environmental correlations between BiW at prenatal THI indicating HS with BiW at remaining THI substantially declined with increasing THI distances. Hence, from a herd management perspective, avoiding HS during the dry period of the dams will contribute to a slight increase in fetus growth.

Animal board invited review: Practical applications of genomic information in livestock

Access to high-dimensional genomic information in many livestock species is accelerating. This has been greatly aided not only by continual reductions in genotyping costs but also an expansion in the services available that leverage genomic information to create a greater return-on-investment. Genomic information on individual animals has many uses including (1) parentage verification and discovery, (2) traceability, (3) karyotyping, (4) sex determination, (5) reporting and monitoring of mutations conferring major effects or congenital defects, (6) better estimating inbreeding of individuals and coancestry among individuals, (7) mating advice, (8) determining breed composition, (9) enabling precision management, and (10) genomic evaluations; genomic evaluations exploit genome-wide genotype information to improve the accuracy of predicting an animal's (and by extension its progeny's) genetic merit. Genomic data also provide a huge resource for research, albeit the outcome from this research, if successful, should eventually be realised through one of the ten applications already mentioned. The process for generating a genotype all the way from sample procurement to identifying erroneous genotypes is described, as are the steps that should be considered when developing a bespoke genotyping panel for practical application.

Realized genomic selection across generations in a reciprocal recurrent selection breeding program of Eucalyptus hybrids

Introduction

Genomic selection (GS) experiments in forest trees have largely reported estimates of predictive abilities from cross-validation among individuals in the same breeding generation. In such conditions, no effects of recombination, selection, drift, and environmental changes are accounted for. Here, we assessed the effectively realized predictive ability (RPA) for volume growth at harvest age by GS across generations in an operational reciprocal recurrent selection (RRS) program of hybrid Eucalyptus.

Methods

Genomic best linear unbiased prediction with additive (GBLUP_G), additive plus dominance (GBLUP_G+D), and additive single-step (HBLUP) models were trained with different combinations of growth data of hybrids and pure species individuals (N = 17,462) of the G1 generation, 1,944 of which were genotyped with ~16,000 SNPs from SNP arrays. The hybrid G2 progeny trial (HPT267) was the GS target, with 1,400 selection candidates, 197 of which were genotyped still at the seedling stage, and genomically predicted for their breeding and genotypic values at the operational harvest age (6 years). Seedlings were then grown to harvest and measured, and their pedigree-based breeding and genotypic values were compared to their originally predicted genomic counterparts.

Results

Genomic RPAs ≥0.80 were obtained as the genetic relatedness between G1 and G2 increased, especially when the direct parents of selection candidates were used in training. GBLUP_G+D reached RPAs ≥0.70 only when hybrid or pure species data of G1 were included in training. HBLUP was only marginally better than GBLUP. Correlations ≥0.80 were obtained between pedigree and genomic individual ranks. Rank coincidence of the top 2.5% selections was the highest for GBLUP_G (45% to 60%) compared to GBLUP_G+D. To advance the pure species RRS populations, GS models were best when trained on pure species than hybrid data, and HBLUP yielded ~20% higher predictive abilities than GBLUP, but was not better than ABLUP for ungenotyped trees.

Discussion

We demonstrate that genomic data effectively enable accurate ranking of eucalypt hybrid seedlings for their yet-to-be observed volume growth at harvest age. Our results support a two-stage GS approach involving family selection by average genomic breeding value, followed by within-top-families individual GS, significantly increasing selection intensity, optimizing genotyping costs, and accelerating RRS breeding.

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.

Predictor bias in genomic and phenomic selection

KEY MESSAGE

NIRS of wheat grains as phenomic predictors for grain yield show inflated prediction ability and are biased toward grain protein content. Estimating the breeding value of individuals using genome-wide marker data (genomic prediction) is currently one of the most important drivers of breeding progress in major crops. Recently, phenomic technologies, including remote sensing and aerial hyperspectral imaging of plant canopies, have made it feasible to predict the breeding value of individuals in the absence of genetic marker data. This is commonly referred to as phenomic prediction. Hyperspectral measurements in the form of near-infrared spectroscopy have been used since the 1980 s to predict compositional parameters of harvest products. Moreover, in recent studies NIRS from grains was used to predict grain yield. The same studies showed that phenomic prediction can outperform genomic prediction for grain yield. The genome is static and not environment dependent, thereby limiting genomic prediction ability. Gene expression is tissue specific and differs under environmental influences, leading to a tissue- and environment-specific phenome, potentially explaining the higher predictive ability of phenomic prediction. Here, we compare genomic prediction and phenomic prediction from hyperspectral measurements of wheat grains for the prediction of a variety of traits including grain yield. We show that phenomic predictions outperform genomic prediction for some traits. However, phenomic predictions are biased toward the information present in the predictor. Future studies on this topic should investigate whether population parameters are retained in phenomic prediction as they are in genomic prediction. Furthermore, we find that unbiased phenomic prediction abilities are considerably lower than previously reported and recommend a method to circumvent this issue.

Genome-wide association study for growth traits in Blanco Orejinegro and Romosinuano cattle

Growth traits are economically important characteristics for the genetic improvement of local cattle breeds. Genome-wide association studies (GWAS) provide valuable information to enhance the understanding on the genetics of complex traits. The aim of this study was to perform a GWAS to identify genomic regions and genes associated to birth weight, weaning weight adjusted for 240 days, 16 months, and 24 months weight in Romosinuano (ROMO) and Blanco Orejinegro (BON) cattle. A single-step genomic-BLUP was implemented using 596 BON and 569 ROMO individuals that were genotyped with an Illumina BovineSNP50 BeadChip. There were 25 regions of interest identified on different chromosomes, with few of them simultaneously associated with two or more growth traits and some were common to both breeds. The gene mapping allowed to find 173 annotations on these regions, from which 49 represent potential candidate genes with known growth-related functions in cattle and other species. Among the regions that were associated with several growth traits, that at 24 - 27 MB of BTA14, has important candidate genes such as LYPLA1, XKR4, TMEM68 and PLAG1. Another region of interest at 0.40-0.77 Mb of BTA23 was identified in both breeds, containing KHDRBS2 as a potential candidate gene influencing body weight. Future studies targeting these regions could provide more knowledge to uncover the genetic architecture underlying growth traits in BON and ROMO cattle. The genomic regions and genes identified in this study could be used to improve the prediction of genetic merit for growth traits in these creole cattle breeds.

Genomic Scanning of Inbreeding Depression for Litter Size in Two Varieties of Iberian Pigs

Inbreeding depression is expected to be more pronounced in fitness-related traits, such as pig litter size. Recent studies have suggested that the genetic determinism of inbreeding depression may be heterogeneous across the genome. Therefore, the objective of this study was to conduct a genomic scan of the whole pig autosomal genome to detect the genomic regions that control inbreeding depression for litter size in two varieties of Iberian pigs (Entrepelado and Retinto). The datasets consisted of 2069 (338 sows) and 2028 (327 sows) records of litter size (Total Number Born and Number Born Alive) for the Entrepelado and Retinto varieties. All sows were genotyped using the Geneseek GGP PorcineHD 70 K chip. We employed the Unfavorable Haplotype Finder software to extract runs of homozygosity (ROHs) and conducted a mixed-model analysis to identify highly significant differences between homozygous and heterozygous sows for each specific ROH. A total of eight genomic regions located on SSC2, SSC5, SSC7, SSC8, and SSC13 were significantly associated with inbreeding depression, housing some relevant genes such as FSHR, LHCGR, CORIN, AQP6, and CEP120.

Interpreting single-step genomic evaluation as a neural network of three layers: pedigree, genotypes, and phenotypes

The single-step approach has become the most widely-used methodology for genomic evaluations when only a subset of phenotyped individuals in the pedigree are genotyped, where the genotypes for non-genotyped individuals are imputed based on gene contents (i.e., genotypes) of genotyped individuals through their pedigree relationships. We proposed a new method named single-step neural network with mixed models (NNMM) to represent single-step genomic evaluations as a neural network of three sequential layers: pedigree, genotypes, and phenotypes. These three sequential layers of information create a unified network instead of two separate steps, allowing the unobserved gene contents of non-genotyped individuals to be sampled based on pedigree, observed genotypes of genotyped individuals, and phenotypes. In addition to imputation of genotypes using all three sources of information, including phenotypes, genotypes, and pedigree, single-step NNMM provides a more flexible framework to allow nonlinear relationships between genotypes and phenotypes, and for individuals to be genotyped with different single-nucleotide polymorphism (SNP) panels. The single-step NNMM has been implemented in the software package "JWAS'.

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.

Factors affecting the accuracy of genomic prediction in joint pig populations

Genomic prediction (GP) has greatly advanced animal and plant breeding over the past two decades. GP in joint populations is a feasible method to improve the accuracy of genomic estimated breeding values in small populations. However, there is still a need to understand the factors that influence GP in joint populations. This study used simulated data and real data from Duroc pig populations to examine the impact of linkage disequilibrium (LD), causal variants effect sizes (CVESs), and minor allele frequencies (MAF) of SNPs on the accuracy of genomic prediction in joint populations. Three prediction methods were used: genomic best linear unbiased prediction (GBLUP), single-step GBLUP and multi-trait GBLUP. Results from the simulated datasets showed that the accuracies of GP in joint populations were always higher than those in a single population when only LD inconsistencies existed. However, single-step GBLUP accuracy in joint populations decreased as the correlation of MAF between populations decreased, while the accuracy of GBLUP is consistently higher in joint populations than in a single population. As the correlation of CVES between populations decreased, the accuracy of both GBLUP and single-step GBLUP in joint populations declined. Analysis of real Duroc populations showed low genetic correlation, similar to the simulated relationship between the most distant populations. In most cases in Duroc populations, GP have higher accuracies in joint populations than in individual population. In conclusion, the consistency of CVES plays a more important role in multi-population GP. The genetic relatedness of the Duroc populations is so weak that the prediction accuracy of GP in joint populations is reduced in some traits. Multi-trait GBLUP is a competitive method for the joint breeding evaluation.

Single-step SNPBLUP evaluation in six German beef cattle breeds

The implementation of genomic selection for six German beef cattle populations was evaluated. Although the multiple-step implementation of genomic selection is the status quo in most national dairy cattle evaluations, the breeding structure of German beef cattle, coupled with the shortcoming and complexity of the multiple-step method, makes single step a more attractive option to implement genomic selection in German beef cattle populations. Our objective was to develop a national beef cattle single-step genomic evaluation in five economically important traits in six German beef cattle populations and investigate its impact on the accuracy and bias of genomic evaluations relative to the current pedigree-based evaluation. Across the six breeds in our study, 461,929 phenotyped and 14,321 genotyped animals were evaluated with a multi-trait single-step model. To validate the single-step model, phenotype data in the last 2 years were removed in a forward validation study. For the conventional and single-step approaches, the genomic estimated breeding values of validation animals and other animals were compared between the truncated and the full evaluations. The correlation of the GEBVs between the full and truncated evaluations in the validation animals was slightly higher in the single-step evaluation. The regression of the full GEBVs on truncated GEBVs was close to the optimal value of 1 for both the pedigree-based and the single-step evaluations. The SNP effect estimates from the truncated evaluation were highly correlated with those from the full evaluation, with values ranging from 0.79 to 0.94. The correlation of the SNP effect was influenced by the number of genotyped animals shared between the full and truncated evaluations. The regression coefficients of the SNP effect of the full evaluation on the truncated evaluation were all close to the expected value of 1, indicating unbiased estimates of the SNP markers for the production traits. The Manhattan plot of the SNP effect estimates identified chromosomal regions harbouring major genes for muscling and body weight in breeds of French origin. Based on the regression intercept and slope of the GEBVs of validation animals, the single-step evaluation was neither inflated nor deflated across the six breeds. Overall, the single-step model resulted in a more accurate and stable evaluation. However, due to the small number of genotyped individuals, the single-step method only provided slightly better results when compared to the pedigree-based method.

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.

Feet and legs malformation in Nellore cattle: genetic analysis and prioritization of GWAS results

Beef cattle affected by feet and legs malformations (FLM) cannot perform their productive and reproductive functions satisfactorily, resulting in significant economic losses. Accelerated weight gain in young animals due to increased fat deposition can lead to ligaments, tendon and joint strain and promote gene expression patterns that lead to changes in the normal architecture of the feet and legs. The possible correlated response in the FLM due to yearling weight (YW) selection suggest that this second trait could be used as an indirect selection criterion. Therefore, FLM breeding values and the genetic correlation between FLM and yearling weight (YW) were estimated for 295,031 Nellore animals by fitting a linear-threshold model in a Bayesian approach. A genome-wide association study was performed to identify genomic windows and positional candidate genes associated with FLM. The effects of single nucleotide polymorphisms (SNPs) on FLM phenotypes (affected or unaffected) were estimated using the weighted single-step genomic BLUP method, based on genotypes of 12,537 animals for 461,057 SNPs. Twelve non-overlapping windows of 20 adjacent SNPs explaining more than 1% of the additive genetic variance were selected for candidate gene annotation. Functional and gene prioritization analysis of candidate genes identified six genes (ATG7, EXT1, ITGA1, PPARD, SCUBE3, and SHOX) that may play a role in FLM expression due to their known role in skeletal muscle development, aberrant bone growth, lipid metabolism, intramuscular fat deposition and skeletogenesis. Identifying genes linked to foot and leg malformations enables selective breeding for healthier herds by reducing the occurrence of these conditions. Genetic markers can be used to develop tests that identify carriers of these mutations, assisting breeders in making informed breeding decisions to minimize the incidence of malformations in future generations, resulting in greater productivity and animal welfare.

Using residual regressions to quantify and map signal leakage in genomic prediction

BACKGROUND

Most genomic prediction applications in animal breeding use genotypes with tens of thousands of single nucleotide polymorphisms (SNPs). However, modern sequencing technologies and imputation algorithms can generate ultra-high-density genotypes (including millions of SNPs) at an affordable cost. Empirical studies have not produced clear evidence that using ultra-high-density genotypes can significantly improve prediction accuracy. However, (whole-genome) prediction accuracy is not very informative about the ability of a model to capture the genetic signals from specific genomic regions. To address this problem, we propose a simple methodology that detects chromosome regions for which a specific model (e.g., single-step genomic best linear unbiased prediction (ssGBLUP)) may fail to fully capture the genetic signal present in such segments-a phenomenon that we refer to as signal leakage. We propose to detect regions with evidence of signal leakage by testing the association of residuals from a pedigree or a genomic model with SNP genotypes. We discuss how this approach can be used to map regions with signals that are poorly captured by a model and to identify strategies to fix those problems (e.g., using a different prior or increasing marker density). Finally, we explored the proposed approach to scan for signal leakage of different models (pedigree-based, ssGBLUP, and various Bayesian models) applied to growth-related phenotypes (average daily gain and backfat thickness) in pigs.

RESULTS

We report widespread evidence of signal leakage for pedigree-based models. Including a percentage of animals with SNP data in ssGBLUP reduced the extent of signal leakage. However, local peaks of missed signals remained in some regions, even when all animals were genotyped. Using variable selection priors solves leakage points that are caused by excessive shrinkage of marker effects. Nevertheless, these models still miss signals in some regions due to low linkage disequilibrium between the SNPs on the array used and causal variants. Thus, we discuss how such problems could be addressed by adding sequence SNPs from those regions to the prediction model.

CONCLUSIONS

Residual single-marker regression analysis is a simple approach that can be used to detect regional genomic signals that are poorly captured by a model and to indicate ways to fix such problems.

A promising resilience parameter for breeding: the use of weight and feed trajectories in growing pigs

BACKGROUND

Increasing resilience is a priority in modern pig breeding. Recent research shows that general resilience can be quantified via variability in longitudinal data. The collection of such longitudinal data on weight, feed intake and feeding behaviour in pigs has been facilitated by the development of technologies such as automated feeding stations. The goal of this study was to investigate resilience traits, which were estimated as deviations from longitudinal weight, feed intake and feeding behaviour data during the finishing phase. A dataset with 324,207 records between the age of 95 and 155 days on 5,939 Piétrain pigs with known pedigree and genomic information was used. We provided guidelines for a rigid quality control of longitudinal body weight data, as we found that outliers can significantly affect results. Gompertz growth curve analysis, linear modelling and trajectory analyses were used for quantifying resilience traits.

RESULTS

To our knowledge, this is the first study comparing resilience traits from longitudinal body weight, feed intake and feeding behaviour data in pigs. We demonstrated that the resilience traits are lowly to moderately heritable for deviations in body weight (h² = 2.9%-20.2%), in feed intake (9.4%-23.3%) and in feeding behaviour (16.2%-28.3%). Additionally, these traits have good predictive abilities in cross-validation analyses. Deviations in individual body weight and feed intake trajectories are highly correlated (r_g = 0.78) with low to moderate favourable genetic correlations with feed conversion ratio (r_g = 0.39-0.49). Lastly, we showed that some resilience traits, such as the natural logarithm of variances of observed versus predicted body weights (lnvar_weight), are more robust to lower observation frequencies and are repeatable over three different time periods of the finishing phase.

CONCLUSIONS

Our results will help future studies investigating resilience traits and resilience-related traits. Moreover, our study provides first results on standardization of quality control and efficient data sampling from automated feeding station data. Our findings will be valuable for breeding organizations as they offer evidence that pigs' general resilience can be selected on with good accuracy. Moreover, this methodology might be extended to other species to quantify resilience based on longitudinal data.

First Genomic Prediction of Single-Step Models in Large Yellow Croaker

Genome selection is mainly used in disease-resistant traits of aquatic species; however, its implementation is hindered by a high cost of genotype and phenotype data collection. Single-step genomic best linear unbiased prediction (SSGBLUP) can integrate phenotypes, genetic markers, and pedigree records into simultaneous prediction without increasing genotyping costs. The objective of this study is to investigate the performance of SSGBLUP in large yellow croaker and to evaluate the effects of the number of phenotypic records and genotyping per family on the predictive ability of SSGBLUP. A large yellow croaker population consists of 6898 individuals from 14 families with survival time resistant against Cryptocaryon irritans (C. irritans), body weight (BW), and body length (BL) traits were collected, of which 669 individuals were genotyped. Results showed that the mean predictive ability of all traits in the individuals randomly sampling for SSGBLUP, GBLUP, and BLUP was 0.738, 0.738, and 0.736, respectively. Moreover, the predictive ability of SSGBLUP and BLUP models did not increase with the extra phenotypic records per family, in which the predictive ability of SSGBLUP and BLUP in survival time was 0.853 and 0.851 for only genotyped data (N = 0) used, and 0.852, 0.845 for all phenotypic records (N = 600) used, respectively. However, with the increase in the genotype number of training set, the prediction ability of SSGBLUP and GBLUP model was increased and the highest predictive ability was gained when the genotype number per family was 40 or 45. In addition, the prediction ability of SSGBLUP model was higher than that of GBLUP. Our study showed that the SSGBLUP model still has great potential and advantages in genomic breeding of large yellow croakers. It is recommended that each family provide 100 phenotypic individuals, of which 40 individuals with genotyping data for SSGBLUP model prediction and family resistance evaluation.

The trend of breeding value research in animal science: bibliometric analysis

This study aims to identify trends and hot topics in breeding value to support researchers in finding new directions for future research in that area. The data of this study consist of 7072 academic studies on breeding value in the Web of Science database. Network visualizations and in-depth bibliometric analysis were performed on cited references, authors, countries, institutions, journals, and keywords through CiteSpace. VanRaden (2008) is the most cited work and has an essential place in the field. The most prolific writer is Ignacy Misztal. While the most productive country in breeding value studies is the United States, the People's Republic of China is an influential country that has experienced a strong citation burst in the last 3 years. The National Institute for Agricultural Research and Wageningen University are important institutions that play a critical role in connecting other institutions. Also, these two institutions have the highest centrality values. "Genomic prediction" is the outstanding sub-study field in the active clusters appearing in the analysis results. We have summarized the literature on breeding value, including publication information, country, institution, author, and journal. We can say that hot topics today are "genome-wide association", "feed efficiency", and "genomic prediction". While the studies conducted in the past years have focused on economic value and accuracy, the studies conducted in recent years have started to be studies that consider technological developments and changing world conditions such as global warming and carbon emission.

Invited review: Good practices in genome-wide association studies to identify candidate sequence variants in dairy cattle

Genotype data from dairy cattle selection programs have greatly facilitated GWAS to identify variants related to economic traits. Results can enhance the accuracy of genomic prediction, analyze more complex models that go beyond additive effects, elucidate the genetic architecture of a trait, and finally, decipher the underlying biology of traits. The entire process, comprising data generation, quality control, statistical analyses, interpretation of association results, and linking results to biology should be designed and executed to minimize the generation of false-positive and false-negative associations and misleading links to biological processes. This review aims to provide general guidelines for data analysis that address data quality control, association tests, adjustment for population stratification, and significance evaluation to improve the reliability of conclusions. We also provide guidance on post-GWAS strategy and the interpretation of results. These guidelines are tailored to dairy cattle, which are characterized by long-range linkage disequilibrium, large half-sib families, and routinely collected phenotypes, requiring different approaches than those applied in human GWAS. We discuss common limitations and challenges that have been overlooked in the analysis and interpretation of GWAS to identify candidate sequence variants in dairy cattle.

Integrated analysis of genome-wide association studies and 3D epigenomic characteristics reveal the BMP2 gene regulating loin muscle depth in Yorkshire pigs

BACKGROUND

The lack of integrated analysis of genome-wide association studies (GWAS) and 3D epigenomics restricts a deep understanding of the genetic mechanisms of meat-related traits. With the application of techniques as ChIP-seq and Hi-C, the annotations of cis-regulatory elements in the pig genome have been established, which offers a new opportunity to elucidate the genetic mechanisms and identify major genetic variants and candidate genes that are significantly associated with important economic traits. Among these traits, loin muscle depth (LMD) is an important one as it impacts the lean meat content. In this study, we integrated cis-regulatory elements and genome-wide association studies (GWAS) to identify candidate genes and genetic variants regulating LMD.

RESULTS

Five single nucleotide polymorphisms (SNPs) located on porcine chromosome 17 were significantly associated with LMD in Yorkshire pigs. A 10 kb quantitative trait locus (QTL) was identified as a candidate functional genomic region through the integration of linkage disequilibrium and linkage analysis (LDLA) and high-throughput chromosome conformation capture (Hi-C) analysis. The BMP2 gene was identified as a candidate gene for LMD based on the integrated results of GWAS, Hi-C meta-analysis, and cis-regulatory element data. The identified QTL region was further verified through target region sequencing. Furthermore, through using dual-luciferase assays and electrophoretic mobility shift assays (EMSA), two SNPs, including SNP rs321846600, located in the enhancer region, and SNP rs1111440035, located in the promoter region, were identified as candidate SNPs that may be functionally related to the LMD.

CONCLUSIONS

Based on the results of GWAS, Hi-C, and cis-regulatory elements, the BMP2 gene was identified as an important candidate gene regulating variation in LMD. The SNPs rs321846600 and rs1111440035 were identified as candidate SNPs that are functionally related to the LMD of Yorkshire pigs. Our results shed light on the advantages of integrating GWAS with 3D epigenomics in identifying candidate genes for quantitative traits. This study is a pioneering work for the identification of candidate genes and related genetic variants regulating one key production trait (LMD) in pigs by integrating genome-wide association studies and 3D epigenomics.

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.

Efficient large-scale single-step evaluations and indirect genomic prediction of genotyped selection candidates

BACKGROUND

Single-step genomic best linear unbiased prediction (ssGBLUP) models allow the combination of genomic, pedigree, and phenotypic data into a single model, which is computationally challenging for large genotyped populations. In practice, genotypes of animals without their own phenotype and progeny, so-called genotyped selection candidates, can become available after genomic breeding values have been estimated by ssGBLUP. In some breeding programmes, genomic estimated breeding values (GEBV) for these animals should be known shortly after obtaining genotype information but recomputing GEBV using the full ssGBLUP takes too much time. In this study, first we compare two equivalent formulations of ssGBLUP models, i.e. one that is based on the Woodbury matrix identity applied to the inverse of the genomic relationship matrix, and one that is based on marker equations. Second, we present computationally-fast approaches to indirectly compute GEBV for genotyped selection candidates, without the need to do the full ssGBLUP evaluation.

RESULTS

The indirect approaches use information from the latest ssGBLUP evaluation and rely on the decomposition of GEBV into its components. The two equivalent ssGBLUP models and indirect approaches were tested on a six-trait calving difficulty model using Irish dairy and beef cattle data that include 2.6 million genotyped animals of which about 500,000 were considered as genotyped selection candidates. When using the same computational approaches, the solving phase of the two equivalent ssGBLUP models showed similar requirements for memory and time per iteration. The computational differences between them were due to the preprocessing phase of the genomic information. Regarding the indirect approaches, compared to GEBV obtained from single-step evaluations including all genotypes, indirect GEBV had correlations higher than 0.99 for all traits while showing little dispersion and level bias.

CONCLUSIONS

In conclusion, ssGBLUP predictions for the genotyped selection candidates were accurately approximated using the presented indirect approaches, which are more memory efficient and computationally fast, compared to solving a full ssGBLUP evaluation. Thus, indirect approaches can be used even on a weekly basis to estimate GEBV for newly genotyped animals, while the full single-step evaluation is done only a few times within a year.

Restriction site-associated DNA sequencing technologies as an alternative to low-density SNP chips for genomic selection: a simulation study in layer chickens

BACKGROUND

To reduce the cost of genomic selection, a low-density (LD) single nucleotide polymorphism (SNP) chip can be used in combination with imputation for genotyping selection candidates instead of using a high-density (HD) SNP chip. Next-generation sequencing (NGS) techniques have been increasingly used in livestock species but remain expensive for routine use for genomic selection. An alternative and cost-efficient solution is to use restriction site-associated DNA sequencing (RADseq) techniques to sequence only a fraction of the genome using restriction enzymes. From this perspective, use of RADseq techniques followed by an imputation step on HD chip as alternatives to LD chips for genomic selection was studied in a pure layer line.

RESULTS

Genome reduction and sequencing fragments were identified on reference genome using four restriction enzymes (EcoRI, TaqI, AvaII and PstI) and a double-digest RADseq (ddRADseq) method (TaqI-PstI). The SNPs contained in these fragments were detected from the 20X sequence data of the individuals in our population. Imputation accuracy on HD chip with these genotypes was assessed as the mean correlation between true and imputed genotypes. Several production traits were evaluated using single-step GBLUP methodology. The impact of imputation errors on the ranking of the selection candidates was assessed by comparing a genomic evaluation based on ancestry using true HD or imputed HD genotyping. The relative accuracy of genomic estimated breeding values (GEBVs) was investigated by considering the GEBVs estimated on offspring as a reference. With AvaII or PstI and ddRADseq with TaqI and PstI, more than 10 K SNPs were detected in common with the HD SNP chip, resulting in an imputation accuracy greater than 0.97. The impact of imputation errors on genomic evaluation of the breeders was reduced, with a Spearman correlation greater than 0.99. Finally, the relative accuracy of GEBVs was equivalent.

CONCLUSIONS

RADseq approaches can be interesting alternatives to low-density SNP chips for genomic selection. With more than 10 K SNPs in common with the SNPs of the HD SNP chip, good imputation and genomic evaluation results can be obtained. However, with real data, heterogeneity between individuals with missing data must be considered.

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.

Multi-line ssGBLUP evaluation using preselected markers from whole-genome sequence data in pigs

Genomic evaluations in pigs could benefit from using multi-line data along with whole-genome sequencing (WGS) if the data are large enough to represent the variability across populations. The objective of this study was to investigate strategies to combine large-scale data from different terminal pig lines in a multi-line genomic evaluation (MLE) through single-step GBLUP (ssGBLUP) models while including variants preselected from whole-genome sequence (WGS) data. We investigated single-line and multi-line evaluations for five traits recorded in three terminal lines. The number of sequenced animals in each line ranged from 731 to 1,865, with 60k to 104k imputed to WGS. Unknown parent groups (UPG) and metafounders (MF) were explored to account for genetic differences among the lines and improve the compatibility between pedigree and genomic relationships in the MLE. Sequence variants were preselected based on multi-line genome-wide association studies (GWAS) or linkage disequilibrium (LD) pruning. These preselected variant sets were used for ssGBLUP predictions without and with weights from BayesR, and the performances were compared to that of a commercial porcine single-nucleotide polymorphisms (SNP) chip. Using UPG and MF in MLE showed small to no gain in prediction accuracy (up to 0.02), depending on the lines and traits, compared to the single-line genomic evaluation (SLE). Likewise, adding selected variants from the GWAS to the commercial SNP chip resulted in a maximum increase of 0.02 in the prediction accuracy, only for average daily feed intake in the most numerous lines. In addition, no benefits were observed when using preselected sequence variants in multi-line genomic predictions. Weights from BayesR did not help improve the performance of ssGBLUP. This study revealed limited benefits of using preselected whole-genome sequence variants for multi-line genomic predictions, even when tens of thousands of animals had imputed sequence data. Correctly accounting for line differences with UPG or MF in MLE is essential to obtain predictions similar to SLE; however, the only observed benefit of an MLE is to have comparable predictions across lines. Further investigation into the amount of data and novel methods to preselect whole-genome causative variants in combined populations would be of significant interest.

Assessment of long-term trends in genetic mean and variance after the introduction of genomic selection in layers: a simulation study

Nucleus-based breeding programs are characterized by intense selection that results in high genetic gain, which inevitably means reduction of genetic variation in the breeding population. Therefore, genetic variation in such breeding systems is typically managed systematically, for example, by avoiding mating the closest relatives to limit progeny inbreeding. However, intense selection requires maximum effort to make such breeding programs sustainable in the long-term. The objective of this study was to use simulation to evaluate the long-term impact of genomic selection on genetic mean and variance in an intense layer chicken breeding program. We developed a large-scale stochastic simulation of an intense layer chicken breeding program to compare conventional truncation selection to genomic truncation selection optimized with either minimization of progeny inbreeding or full-scale optimal contribution selection. We compared the programs in terms of genetic mean, genic variance, conversion efficiency, rate of inbreeding, effective population size, and accuracy of selection. Our results confirmed that genomic truncation selection has immediate benefits compared to conventional truncation selection in all specified metrics. A simple minimization of progeny inbreeding after genomic truncation selection did not provide any significant improvements. Optimal contribution selection was successful in having better conversion efficiency and effective population size compared to genomic truncation selection, but it must be fine-tuned for balance between loss of genetic variance and genetic gain. In our simulation, we measured this balance using trigonometric penalty degrees between truncation selection and a balanced solution and concluded that the best results were between 45° and 65°. This balance is specific to the breeding program and depends on how much immediate genetic gain a breeding program may risk vs. save for the future. Furthermore, our results show that the persistence of accuracy is better with optimal contribution selection compared to truncation selection. In general, our results show that optimal contribution selection can ensure long-term success in intensive breeding programs using genomic selection.

FSF-GA: A Feature Selection Framework for Phenotype Prediction Using Genetic Algorithms

(1) Background: Phenotype prediction is a pivotal task in genetics in order to identify how genetic factors contribute to phenotypic differences. This field has seen extensive research, with numerous methods proposed for predicting phenotypes. Nevertheless, the intricate relationship between genotypes and complex phenotypes, including common diseases, has resulted in an ongoing challenge to accurately decipher the genetic contribution. (2) Results: In this study, we propose a novel feature selection framework for phenotype prediction utilizing a genetic algorithm (FSF-GA) that effectively reduces the feature space to identify genotypes contributing to phenotype prediction. We provide a comprehensive vignette of our method and conduct extensive experiments using a widely used yeast dataset. (3) Conclusions: Our experimental results show that our proposed FSF-GA method delivers comparable phenotype prediction performance as compared to baseline methods, while providing features selected for predicting phenotypes. These selected feature sets can be used to interpret the underlying genetic architecture that contributes to phenotypic variation.

Validation of single-step genomic BLUP random regression test-day models and SNP effects analysis on milk yield in French Saanen goats

The shape of the lactation curve is linked to an animal's health, feed requirements, and milk production throughout the year. Random regression models (RRM) are widely used for genetic evaluation of total milk production throughout the lactation and for milk yield persistency. Genomic information used with the single-step genomic BLUP method (ssGBLUP) substantially improves the accuracy of genomic prediction of breeding values in the main dairy cattle breeds. The aim of this study was to implement an RRM using ssGBLUP for milk yield in Saanen dairy goats in France. The data set consisted of 7,904,246 test-day records from 1,308,307 lactations of Saanen goats collected in France between 2000 and 2017. The performance of this type of evaluation was assessed by applying a validation step with data targeting candidate bucks. The model was compared with a nongenomic evaluation and a traditional evaluation that use cumulated performance throughout the lactation model (LM). The incorporation of genomic information increased correlations between daughter yield deviations (DYD) and estimated breeding values (EBV) obtained with a partial data set for candidate bucks. The LM and the RRM had similar correlation between DYD and EBV. However, the RRM reduced overestimation of EBV and improved the slope of the regression of DYD on EBV obtained at birth. This study shows that a genomic evaluation from a ssGBLUP RRM is possible in dairy goats in France and that RRM performance is comparable to a LM but with the additional benefit of a genomic evaluation of persistency. Variance of adjacent SNPs was studied with LM and RRM following the ssGBLUP. Both approaches converged on approximately the same regions explaining more than 1% of total variance. Regions associated with persistency were also found.

PyAGH: a python package to fast construct kinship matrices based on different levels of omic data

BACKGROUND

Construction of kinship matrices among individuals is an important step for both association studies and prediction studies based on different levels of omic data. Methods for constructing kinship matrices are becoming diverse and different methods have their specific appropriate scenes. However, software that can comprehensively calculate kinship matrices for a variety of scenarios is still in an urgent demand.

RESULTS

In this study, we developed an efficient and user-friendly python module, PyAGH, that can accomplish (1) conventional additive kinship matrces construction based on pedigree, genotypes, abundance data from transcriptome or microbiome; (2) genomic kinship matrices construction in combined population; (3) dominant and epistatic effects kinship matrices construction; (4) pedigree selection, tracing, detection and visualization; (5) visualization of cluster, heatmap and PCA analysis based on kinship matrices. The output from PyAGH can be easily integrated in other mainstream software based on users' purposes. Compared with other softwares, PyAGH integrates multiple methods for calculating the kinship matrix and has advantages in terms of speed and data size compared to other software. PyAGH is developed in python and C +  + and can be easily installed by pip tool. Installation instructions and a manual document can be freely available from https://github.com/zhaow-01/PyAGH .

CONCLUSION

PyAGH is a fast and user-friendly Python package for calculating kinship matrices using pedigree, genotype, microbiome and transcriptome data as well as processing, analyzing and visualizing data and results. This package makes it easier to perform predictions and association studies processes based on different levels of omic data.

A potent allele marker related to low bull conception rate in Japanese Black bulls

Over the years, there has been considerable variation in the bull conception rate (BCR) of Japanese Black cattle; moreover, several Japanese Black bulls with a low BCR of ≤10% have been identified. However, the alleles responsible for the low BCR are not determined yet. Therefore, in this study, we aimed to identify single-nucleotide polymorphisms (SNPs) for predicting low BCR. To this end, the genome of Japanese Black bulls was comprehensively examined by a genome-wide association study with whole-exome sequencing (WES), and the effect of the identified marker regions on BCR was determined. The WES analysis of six sub-fertile bulls with a BCR of ≤10% and 73 normal bulls with a BCR of ≥40% identified a homozygous genotype for low BCR in Bos taurus autosome 5 in the region between 116.2 and 117.9 Mb. The g.116408653G > A SNP in this region had the most significant effect on the BCR (P-value = 1.0 × 10^-23), and the GG (55.4 ± 11.2%) and AG (54.4 ± 9.4%) genotypes in the SNP had a higher phenotype than the AA (9.5 ± 6.1%) genotype for the BCR. The mixed model analysis revealed that g.116408653G > A was related to approximately 43% of the total genetic variance. In conclusion, the AA genotype of g.116408653G > A is a useful index for identifying sub-fertile Japanese Black bulls. Some positive and negative effects of SNP on the BCR were presumed to identify the causative mutations, which can help evaluate bull fertility.

Fully efficient, two-stage analysis of multi-environment trials with directional dominance and multi-trait genomic selection

R/StageWise enables fully efficient, two-stage analysis of multi-environment, multi-trait datasets for genomic selection, including support for dominance heterosis and polyploidy. Plant breeders interested in genomic selection often face challenges to fully utilizing multi-trait, multi-environment datasets. R package StageWise was developed to go beyond the capabilities of most specialized software for genomic prediction, without requiring the programming skills needed for more general-purpose software for mixed models. As the name suggests, one of the core features is a fully efficient, two-stage analysis for multiple environments, in which the full variance-covariance matrix of the Stage 1 genotype means is used in Stage 2. Another feature is directional dominance, including for polyploids, to account for inbreeding depression in outbred crops. StageWise enables selection with multi-trait indices, including restricted indices with one or more traits constrained to have zero response. For a potato dataset with 943 genotypes evaluated over 6 years, including the Stage 1 errors in Stage 2 reduced the Akaike Information Criterion (AIC) by 29, 67, and 104 for maturity, yield, and fry color, respectively. The proportion of variation explained by heterosis was largest for yield but still only 0.03, likely because of limited variation for the genomic inbreeding coefficient. Due to the large additive genetic correlation (0.57) between yield and maturity, naïve selection on an index combining yield and fry color led to an undesirable response for later maturity. The restricted index coefficients to maximize genetic merit without delaying maturity were identified. The software and three vignettes are available at https://github.com/jendelman/StageWise .

Increasing genomic prediction accuracy for unphenotyped full-sib families by modeling additive and dominance effects with large datasets in white spruce

Introduction

Genomic selection is becoming a standard technique in plant breeding and is now being introduced into forest tree breeding. Despite promising results to predict the genetic merit of superior material based on their additive breeding values, many studies and operational programs still neglect non-additive effects and their potential for enhancing genetic gains.

Methods

Using two large comprehensive datasets totaling 4,066 trees from 146 full-sib families of white spruce (Picea glauca (Moench) Voss), we evaluated the effect of the inclusion of dominance on the precision of genetic parameter estimates and on the accuracy of conventional pedigree-based (ABLUP-AD) and genomic-based (GBLUP-AD) models.

Results

While wood quality traits were mostly additively inherited, considerable non-additive effects and lower heritabilities were detected for growth traits. For growth, GBLUP-AD better partitioned the additive and dominance effects into roughly equal variances, while ABLUP-AD strongly overestimated dominance. The predictive abilities of breeding and total genetic value estimates were similar between ABLUP-AD and GBLUP-AD when predicting individuals from the same families as those included in the training dataset. However, GBLUP-AD outperformed ABLUP-AD when predicting for new unphenotyped families that were not represented in the training dataset, with, on average, 22% and 53% higher predictive ability of breeding and genetic values, respectively. Resampling simulations showed that GBLUP-AD required smaller sample sizes than ABLUP-AD to produce precise estimates of genetic variances and accurate predictions of genetic values. Still, regardless of the method used, large training datasets were needed to estimate additive and non-additive genetic variances precisely.

Discussion

This study highlights the different quantitative genetic architectures between growth and wood traits. Furthermore, the usefulness of genomic additive-dominance models for predicting new families should allow practicing mating allocation to maximize the total genetic values for the propagation of elite material.

Genomic evaluation for two-way crossbred performance in cattle

BACKGROUND

Dairy cattle production systems are mostly based on purebreds, but recently the use of crossbreeding has received increased interest. For genetic evaluations including crossbreds, several methods based on single-step genomic best linear unbiased prediction (ssGBLUP) have been proposed, including metafounder ssGBLUP (MF-ssGBLUP) and breed-specific ssGBLUP (BS-ssGBLUP). Ideally, models that account for breed effects should perform better than simple models, but knowledge on the performance of these methods is lacking for two-way crossbred cattle. In addition, the differences in the estimates of genetic parameters (such as the genetic variance component and heritability) between these methods have rarely been investigated. Therefore, the aims of this study were to (1) compare the estimates of genetic parameters for average daily gain (ADG) and feed conversion ratio (FCR) between these methods; and (2) evaluate the impact of these methods on the predictive ability for crossbred performance.

METHODS

Bivariate models using standard ssGBLUP, MF-ssGBLUP and BS-ssGBLUP for the genetic evaluation of ADG and FCR were investigated. To measure the predictive ability of these three methods, we estimated four estimators, bias, dispersion, population accuracy and ratio of population accuracies, using the linear regression (LR) method.

RESULTS

The results show that, for both ADG and FCR, the heritabilities were low with the three methods. For FCR, the differences in the estimated genetic parameters were small between the three methods, while for ADG, those estimated with BS-ssGBLUP deviated largely from those estimated with the other two methods. Bias and dispersion were similar across the three methods. Population accuracies for both ADG and FCR were always higher with MF-ssGBLUP than with ssGBLUP, while with BS-ssGBLUP the population accuracy was highest for FCR and lowest for ADG.

CONCLUSIONS

Our results indicate that in the genetic evaluation for crossbred performance in a two-way crossbred cattle production system, the predictive ability of MF-ssGBLUP and BS-ssGBLUP is greater than that of ssGBLUP, when the estimated variance components are consistent across the three methods. Compared with BS-ssGBLUP, MF-ssGBLUP is more robust in its superiority over ssGBLUP.

Single-step genome-wide association analyses of claw horn lesions in Holstein cattle using linear and threshold models

BACKGROUND

Lameness in dairy cattle is primarily caused by foot lesions including the claw horn lesions (CHL) of sole haemorrhage (SH), sole ulcers (SU), and white line disease (WL). This study investigated the genetic architecture of the three CHL based on detailed animal phenotypes of CHL susceptibility and severity. Estimation of genetic parameters and breeding values, single-step genome-wide association analyses, and functional enrichment analyses were performed.

RESULTS

The studied traits were under genetic control with a low to moderate heritability. Heritability estimates of SH and SU susceptibility on the liability scale were 0.29 and 0.35, respectively. Heritability of SH and SU severity were 0.12 and 0.07, respectively. Heritability of WL was relatively lower, indicating stronger environmental influence on the presence and development of WL than the other two CHL. Genetic correlations between SH and SU were high (0.98 for lesion susceptibility and 0.59 for lesion severity), whereas genetic correlations of SH and SU with WL also tended to be positive. Candidate quantitative trait loci (QTL) were identified for all CHL, including some on Bos taurus chromosome (BTA) 3 and 18 with potential pleiotropic effects associated with multiple foot lesion traits. A genomic window of 0.65 Mb on BTA3 explained 0.41, 0.50, 0.38, and 0.49% of the genetic variance for SH susceptibility, SH severity, WL susceptibility, and WL severity, respectively. Another window on BTA18 explained 0.66, 0.41, and 0.70% of the genetic variance for SH susceptibility, SU susceptibility, and SU severity, respectively. The candidate genomic regions associated with CHL harbour annotated genes that are linked to immune system function and inflammation responses, lipid metabolism, calcium ion activities, and neuronal excitability.

CONCLUSIONS

The studied CHL are complex traits with a polygenic mode of inheritance. Most traits exhibited genetic variation suggesting that animal resistance to CHL can be improved with breeding. The CHL traits were positively correlated, which will facilitate genetic improvement for resistance to CHL as a whole. Candidate genomic regions associated with lesion susceptibility and severity of SH, SU, and WL provide insights into a global profile of the genetic background underlying CHL and inform genetic improvement programmes aiming at enhancing foot health in dairy cattle.

Selective genotyping strategies for a sib test scheme of a broiler breeder program

BACKGROUND

In broiler breeding, genotype-by-environment interaction is known to result in a genetic correlation between body weight measured in bio-secure and commercial environments that is substantially less than 1. Thus, measuring body weights on sibs of selection candidates in a commercial environment and genotyping them could increase genetic progress. Using real data, the aim of this study was to evaluate which genotyping strategy and which proportion of sibs placed in the commercial environment should be genotyped to optimize a sib-testing breeding program in broilers. Phenotypic body weight and genomic information were collected on all sibs raised in a commercial environment, which allowed to retrospectively analyze different sampling strategies and genotyping proportions.

RESULTS

Accuracies of genomic estimated breeding values (GEBV) obtained with the different genotyping strategies were assessed by computing their correlation with GEBV obtained when all sibs in the commercial environment were genotyped. Results showed that, compared to random sampling (RND), genotyping sibs with extreme phenotypes (EXT) resulted in higher GEBV accuracy across all genotyping proportions, especially for genotyping proportions of 12.5% or 25%, which resulted in correlations of 0.91 vs 0.88 for 12.5% and 0.94 vs 0.91 for 25% genotyped. Including pedigree on birds with phenotype in the commercial environment that were not genotyped increased accuracy at lower genotyping proportions, especially for the RND strategy (correlations of 0.88 vs 0.65 at 12.5% and 0.91 vs 0.80 at 25%), and a smaller but still substantial increase in accuracy for the EXT strategy (0.91 vs 0.79 for 12.5% and 0.94 vs 0.88 for 25% genotyped). Dispersion bias was virtually absent for RND if 25% or more birds were genotyped. However, GEBV were considerably inflated for EXT, especially when the proportion genotyped was low, which was further exacerbated if the pedigree of non-genotyped sibs was excluded.

CONCLUSIONS

When less than 75% of all animals placed in a commercial environment are genotyped, it is recommended to use the EXT strategy, because it yields the highest accuracy. However, caution should be taken when interpreting the resulting GEBV because they will be over-dispersed. When 75% or more of the animals are genotyped, random sampling is recommended because it yields virtually no bias of GEBV and results in similar accuracies as the EXT strategy.

Invited review: Reliability computation from the animal model era to the single-step genomic model era

The calculation of exact reliabilities involving the inversion of mixed model equations poses a heavy computational challenge when the system of equations is large. This has prompted the development of different approximation methods. We give an overview of the various methods and computational approaches in calculating reliability from the era before the animal model to the era of single-step genomic models. The different methods are discussed in terms of modeling, development, and applicability in large dairy cattle populations. The paper also describes the problems faced in reliability computation. Many details dispersed throughout the literature are presented in this paper. It is clear that a universal solution applicable to every model and input data may not be possible, but we point out several efficient and accurate algorithms developed recently for a variety of very large genomic evaluations.

Genetic parameters of sole lesion recovery in Holstein cows

Sole hemorrhage and sole ulcers, referred to as sole lesions, are important causes of lameness in dairy cattle. The objective of this study was to estimate the genetic parameters of a novel trait reflecting how well cows recovered from sole lesions and the genetic correlation of this trait with overall susceptibility to sole lesions. A cohort of Holstein dairy cows was prospectively enrolled on 4 farms and assessed at 4 timepoints: before calving, immediately after calving, in early lactation, and in late lactation. At each timepoint, sole lesions were recorded at the claw level by veterinary surgeons and used to define 2 binary traits: (1) susceptibility to sole lesions-whether animals were affected with sole lesions at least once during the study or were unaffected at every assessment, and (2) sole lesion recovery-whether sole lesions healed between early and late lactation. Animals were genotyped and pedigree details extracted from the national database. Analyses were conducted with BLUPF90 software in a single-step framework; genetic parameters were estimated from animal threshold models using Gibbs sampling. The genetic correlation between both traits was approximated as the correlation between genomic estimated breeding values, adjusting for their reliabilities. A total of 2,025 animals were used to estimate the genetic parameters of sole lesion susceptibility; 44% of animals recorded a sole lesion at least once during the study period. The heritability of sole lesion susceptibility, on the liability scale, was 0.25 (95% highest density interval = 0.16-0.34). A total of 498 animals were used to estimate the genetic parameters of sole lesion recovery; 71% of animals had recovered between the early and late lactation assessments. The heritability of sole lesion recovery, on the liability scale, was 0.27 (95% highest density interval = 0.02-0.52). The approximate genetic correlation between each trait was -0.11 (95% confidence interval = -0.20 to -0.02). Our results indicate that recovery from sole lesions is heritable. If this finding is corroborated in further studies, it may be possible to use selective breeding to reduce the frequency of chronically lame cows. As sole lesion recovery appears to be weakly genetically related to sole lesion susceptibility, successful genetic improvement of sole lesion recovery would benefit from selection on this trait directly.

Applying genetic technologies to combat infectious diseases in aquaculture

Disease and parasitism cause major welfare, environmental and economic concerns for global aquaculture. In this review, we examine the status and potential of technologies that exploit genetic variation in host resistance to tackle this problem. We argue that there is an urgent need to improve understanding of the genetic mechanisms involved, leading to the development of tools that can be applied to boost host resistance and reduce the disease burden. We draw on two pressing global disease problems as case studies-sea lice infestations in salmonids and white spot syndrome in shrimp. We review how the latest genetic technologies can be capitalised upon to determine the mechanisms underlying inter- and intra-species variation in pathogen/parasite resistance, and how the derived knowledge could be applied to boost disease resistance using selective breeding, gene editing and/or with targeted feed treatments and vaccines. Gene editing brings novel opportunities, but also implementation and dissemination challenges, and necessitates new protocols to integrate the technology into aquaculture breeding programmes. There is also an ongoing need to minimise risks of disease agents evolving to overcome genetic improvements to host resistance, and insights from epidemiological and evolutionary models of pathogen infestation in wild and cultured host populations are explored. Ethical issues around the different approaches for achieving genetic resistance are discussed. Application of genetic technologies and approaches has potential to improve fundamental knowledge of mechanisms affecting genetic resistance and provide effective pathways for implementation that could lead to more resistant aquaculture stocks, transforming global aquaculture.