Single-Step Genomic Evaluations from Theory to Practice: Using SNP Chips and Sequence Data in BLUPF90

Estimation of genetic parameters and genome-wide association study for enzootic bovine leukosis resistance in Canadian Holstein cattle

Bovine leukemia virus (BLV) is the causative agent of enzootic bovine leukosis (leukosis) frequently observed in North American dairy herds. Infection with BLV can lead to persistent lymphocytosis and tumors, and is associated with decreased production, immunity and fertility. With no available treatment or vaccine, reducing the prevalence of leukosis through management and culling has not yet been successful. Genetic selection could contribute to permanent improvement in dairy cattle resistance to leukosis. This study aimed to examine the prevalence and impact of leukosis in Canada, and to assess the potential for including leukosis resistance in Canadian national genetic evaluations by characterizing the genetic architecture of leukosis resistance using pedigree and genomic information. A total of 117,349 milk enzyme-linked immunosorbent assay test records on 96,779 Holstein cows from 950 Canadian herds taken between 2007 and 2021 were provided by Lactanet Canada (Guelph, ON, Canada). Each cow was classified as test-positive for leukosis or test-negative for leukosis. Leukosis was present in approximately 77% of herds tested; within those herds, an average of 39% of cows tested were test positive for leukosis. Heritabilities of 0.10 (SE = 0.001) and 0.07 (SE <0.001) were estimated for leukosis resistance using a linear animal model and BLUP or single-step GBLUP methodology, respectively. Breeding value correlations were estimated between leukosis resistance and economically important and phenotypically relevant traits. Most correlations between leukosis resistance and traits already included in Canadian genetic evaluations were favorable, with the exception of somatic cell score. The candidate genes for leukosis resistance identified using a genome-wide association study, were on chromosome 23, with some being part of the major histocompatibility complex. This study showed that genetic evaluation for leukosis resistance is possible, and could be considered for inclusion in Canadian national selection indices.

Impact of truncating historical data on prediction ability of dairy sheep selection candidates

Along the last decades, the genetic evaluation methodology has evolved, improving breeding value estimates. Many breeding programmes have historical phenotypic records and large number of generations, but to make use of them could result in more inconveniences than benefits. In this study, the prediction ability of genotyped young animals was assessed by simultaneously evaluating the removal of historical data, two pedigree deepness and two methodologies (traditional BLUP and single-step genomic BLUP or ssGBLUP), using milk yield records of 40 years of three Latxa dairy sheep populations. The linear regression method was used to compare predictions of young rams before and after progeny testing, with six cut-off points, by intervals of 4 years (from 1992 to 2012), and statistics of ratio of accuracies, bias, and dispersion were calculated. The prediction accuracy of selection candidates, when genomic information was included, was the highest in all Latxa populations (between 0.54 and 0.69 with full data set). Nevertheless, the deletion of historical phenotypic data resulted on moderate accuracy gain in the bigger data size populations (mean gain 2.5%), and the smaller population took advantage of a moderate data deletion (2.7% gain by removing data until 2004), reducing accuracy when more records were removed. The bias of validation individuals was lower when the breeding value was predicted based on genomic information (between 2.1 and 13.9), being lower when the biggest amount of data was deleted in the bigger data size populations (5.2% reduction), and the smaller population was benefited from data deletion between 1996 and 2008 (3.8% bias reduction). Meanwhile, the slope of estimated genetic trend was lower when less data were included, and an overestimation of the unknown parent group estimates was observed. The results indicated that ssGBLUP evaluations were outstanding, compared with traditional BLUP evaluations, while the depth of pedigree had a very small influence, and deletion of historical phenotypic data was beneficial. Thus, Latxa routine genetic evaluations would benefit from truncating phenotypic records between 2000 and 2004, the use of two pedigree generations and the implementation of ssGBLUP methodology.

Breed of origin analysis in genome-wide association studies: enhancing SNP-based insights into production traits in a commercial Brangus population

BACKGROUND

Carcass weight (HCW) and marbling (MARB) are critical for meat quality and market value in beef cattle. In composite breeds like Brangus, which meld the genetics of Angus and Brahman, SNP-based analyses have illuminated some genetic influences on these traits, but they fall short in fully capturing the nuanced effects of breed of origin alleles (BOA) on these traits. Focus on the impacts of BOA on phenotypic features within Brangus populations can result in a more profound understanding of the specific influences of Angus and Brahman genetics. Moreover, the consideration of BOA becomes particularly significant when evaluating dominance effects contributing to heterosis in crossbred populations. BOA provides a more comprehensive measure of heterosis due to its ability to differentiate the distinct genetic contributions originating from each parent breed. This detailed understanding of genetic effects is essential for making informed breeding decisions to optimize the benefits of heterosis in composite breeds like Brangus.

OBJECTIVE

This study aims to identify quantitative trait loci (QTL) influencing HCW and MARB by utilizing SNP and BOA information, incorporating additive, dominance, and overdominance effects within a multi-generational Brangus commercial herd.

METHODS

We analyzed phenotypic data from 1,066 genotyped Brangus steers. BOA inference was performed using LAMP-LD software using Angus and Brahman reference sets. SNP-based and BOA-based GWAS were then conducted considering additive, dominance, and overdominance models.

RESULTS

The study identified numerous QTLs for HCW and MARB. A notable QTL for HCW was associated to the SGCB gene, pivotal for muscle growth, and was identified solely in the BOA GWAS. Several BOA GWAS QTLs exhibited a dominance effect underscoring their importance in estimating heterosis.

CONCLUSIONS

Our findings demonstrate that SNP-based methods may not detect all genetic variation affecting economically important traits in composite breeds. BOA inclusion in genomic evaluations is crucial for identifying genetic regions contributing to trait variation and for understanding the dominance value underpinning heterosis. By considering BOA, we gain a deeper understanding of genetic interactions and heterosis, which is integral to advancing breeding programs. The incorporation of BOA is recommended for comprehensive genomic evaluations to optimize trait improvements in crossbred cattle populations.

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

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

Machine learning methods for genomic prediction of cow behavioral traits measured by automatic milking systems in North American Holstein cattle

Identifying genome-enabled methods that provide more accurate genomic prediction is crucial when evaluating complex traits such as dairy cow behavior. In this study, we aimed to compare the predictive performance of traditional genomic prediction methods and deep learning algorithms for genomic prediction of milking refusals (MREF) and milking failures (MFAIL) in North American Holstein cows measured by automatic milking systems (milking robots). A total of 1,993,509 daily records from 4,511 genotyped Holstein cows were collected by 36 milking robot stations. After quality control, 57,600 SNPs were available for the analyses. Four genomic prediction methods were considered: Bayesian least absolute shrinkage and selection operator (LASSO), multiple layer perceptron (MLP), convolutional neural network (CNN), and GBLUP. We implemented the first 3 methods using the Keras and TensorFlow libraries in Python (v.3.9) but the GBLUP method was implemented using the BLUPF90+ family programs. The accuracy of genomic prediction (mean square error) for MREF and MFAIL was 0.34 (0.08) and 0.27 (0.08) based on LASSO, 0.36 (0.09) and 0.32 (0.09) for MLP, 0.37 (0.08) and 0.30 (0.09) for CNN, and 0.35 (0.09) and 0.31(0.09) based on GBLUP, respectively. Additionally, we observed a lower reranking of top selected individuals based on the MLP versus CNN methods compared with the other approaches for both MREF and MFAIL. Although the deep learning methods showed slightly higher accuracies than GBLUP, the results may not be sufficient to justify their use over traditional methods due to their higher computational demand and the difficulty of performing genomic prediction for nongenotyped individuals using deep learning procedures. Overall, this study provides insights into the potential feasibility of using deep learning methods to enhance genomic prediction accuracy for behavioral traits in livestock. Further research is needed to determine their practical applicability to large dairy cattle breeding programs.

Modeling missing pedigree with metafounders and validating single-step genomic predictions in a small dairy cattle population with a great influence of foreign genetics

Genetic improvement in small countries rely heavily on foreign genetics. In an importing country such as Uruguay, consideration of unknown parent groups (UPG) for foreign sires is essential. However, the use of UPG in genomic model evaluations may lead to bias in genomic estimated breeding values (GEBV). The objective of this study was to study different models including UPG or metafounders (MF) in the Uruguayan Holstein evaluation and to analyze bias, dispersion, and accuracy of GEBV predictions in BLUP and single-step genomic BLUP (ssGBLUP). A gamma matrix (Γ) was estimated either by using base allele population frequencies obtained by bounded linear regression (MFbounded), or by using 2 values to design Γ (i.e., a single value for the diagonal and a different value for the off-diagonal [MFrobust]). Both Γ estimators performed well in terms of GEBV predictions, but MFbounded was the best option. There is, however, some bias whose origin was not completely understood. UPG or MF seem to model correctly genetic progress for unknown parents except for the very first groups (earlier time period). As for validation bulls, bias was observed across all models, whereas for validation cows it was only observed with UPG in BLUP. Overdispersion was found in all models, but it was mostly detected in validation bulls. Ratio of accuracies indicated that ssGBLUP gave better predictions than BLUP.

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

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

Study on Changes in Gut Microbiota and Microbiability in Rabbits at Different Developmental Stages

This study used feces from 0-day-old (36 rabbits), 10-day-old (119 rabbits), and 60-day-old (119 rabbits) offspring rabbits and their corresponding female rabbits (36 rabbits) as experimental materials. Using 16s rRNA sequencing, the study analyzed the types and changes of gut microbiota in rabbits at different growth and development stages, as well as the correlation between gut microbiota composition and the weight of 60-day-old rabbits. All experimental rabbits were placed in the same rabbit shed. Juvenile rabbits were fed solid feed at 18 days of age and weaned at 35 days of age. In addition to identifying the dominant bacterial phyla of gut microbiota in rabbits at different age stages, it was found that the abundance of Clostridium tertium and Clostridium paraputrificum in all suckling rabbits (10-day-old) was significantly higher than that in rabbits fed with whole feed (60-day-old) (p < 0.05), while the abundance of Gram-negative bacterium cTPY13 was significantly lower (p < 0.05). In addition, Fast Expected Maximum Microbial Source Tracing (FEAST) analysis showed that the contribution of female rabbits' gut microbiota to the colonization of offspring rabbits' gut microbiota was significantly higher than that of unrelated rabbits' gut microbiota (p < 0.05). The contribution of female rabbits' gut microbiota to the colonization of gut microbiota in 0-day-old rabbits was significantly higher than that to the colonization of gut microbiota in the 10- and 60-day-old rabbits (p < 0.05). Finally, the correlation between gut microbiota composition and body weight of 60-day-old rabbits was analyzed based on a mixed linear model, and six ASVs significantly affecting body weight were screened. The above results provide important theoretical and practical guidance for maintaining gut health, improving growth and development performance, and feeding formulation in rabbits.

Genetic evaluations and genome-wide association studies for specific digital dermatitis diagnoses in dairy cows considering genotype × housing system interactions

The present study aimed to use detailed phenotyping for the claw disorder digital dermatitis (DD) considering specific DD stages in 2 housing systems (conventional cubicle barns [CON] and compost-bedded pack barns [CBPB]) to infer possible genotype × housing system interactions. The DD stages included 2,980 observations for the 3 traits DD-sick, DD-acute, and DD-chronic from 1,311 Holstein-Friesian and 399 Fleckvieh-Simmental cows. Selection of the 5 CBPB and 5 CON herds was based on a specific protocol to achieve a high level of herd similarity with regard to climate, feeding, milking system, and location, but with pronounced housing-system differences. Five other farms had a "mixed system" with 2 subherds, one representing CBPB and the other one CON. The CBPB system was represented by 899 cows (1,530 observations), and 811 cows (1,450 observations) represented the CON system. The average disease prevalence was 20.47% for DD-sick, 13.88% for DD-acute, and 5.34% for DD-chronic, with a higher prevalence in CON than in CBPB. After quality control of 50K genotypes, 38,495 SNPs from 926 cows remained for the ongoing genomic analyses. Genetic parameters for DD-sick, DD-acute, and DD-chronic were estimated by applying single-step approaches for single-trait repeatability animal models considering the whole dataset, and separately for the CON and CBPB subsets. Genetic correlations between same DD traits from different housing systems, and between DD-sick, DD-chronic, and DD-acute, were estimated via bivariate animal models. Heritabilities based on the whole dataset were 0.16 for DD-sick, 0.14 for DD-acute, and 0.11 for DD-chronic. A slight increase of heritabilities and genetic variances was observed in CON compared with the "well-being" CBPB system, indicating a stronger genetic differentiation of diseases in a more challenging environment. Genetic correlations between same DD traits recorded in CON or CBPB were close to 0.80, disproving obvious genotype × housing system interactions. Genetic correlations among DD-sick, DD-acute and DD-chronic ranged from 0.58 to 0.81. SNP main effects and SNP × housing system interaction effects were estimated simultaneously via GWAS, considering only the phenotypes from genotyped cows. Ongoing annotations of potential candidate genes focused on chromosomal segments 100 kb upstream and downstream from the significantly associated candidate SNP. GWAS for main effects indicated heterogeneous Manhattan plots especially for DD-acute and DD-chronic, indicating particularities in disease pathogenesis. Nevertheless, a few shared annotated potential candidate genes, that is, METTL25, AFF3, PRKG1, and TENM4 for DD-sick and DD-acute, were identified. These genes have direct or indirect effects on disease resistance or immunology. For the SNP × housing system interaction, the annotated genes ASXL1 and NOL4L on BTA 13 were relevant for DD-sick and DD-acute. Overall, the very similar genetic parameters for the same traits in different environments and negligible genotype × housing system interactions indicate only minor effects on genetic evaluations for DD due to housing-system particularities.

Unknown parent groups and truncated pedigree in single-step genomic evaluations of Murrah buffaloes

Missing pedigree may produce bias in genomic evaluations. Thus, strategies to deal with this problem have been proposed as using unknown parent groups (UPG) or truncated pedigrees. The aim of this study was to investigate the impact of modeling missing pedigree under ssGBLUP evaluations for productive and reproductive traits in dairy buffalos using different approaches: 1) traditional BLUP without UPG (BLUP), 2) traditional BLUP including UPG (BLUP/UPG), 3) ssGBLUP without UPG (ssGBLUP), 4) ssGBLUP including UPG in the A and A22 matrices (ssGBLUP/A_UPG), 5) ssGBLUP including UPG in all elements of the H matrix (ssGBLUP/H_UPG), 6) BLUP with pedigree truncation for the last 3 generations (BLUP/truncated), and 7) ssGBLUP with pedigree truncation for the last 3 generations (ssGBLUP/ truncated). UPGs were not used in the scenarios with truncated pedigree. A total of 3,717, 4,126 and 3,823 records of the first lactation for accumulated 305 d milk yield (MY), age at first calving (AFC) and lactation length (LL), respectively were used. Accuracies ranged from 0.27 for LL (BLUP) to 0.46 for MY (BLUP), bias ranged from -0.62 for MY (ssGBLUP) to 0.0002 for AFC (BLUP/truncated), and dispersion ranged from 0.88 for MY (BLUP/ A_UPG) to 1.13 for LL (BLUP). Genetic trend showed genetic gains for all traits across 20 years of selection and the impact of including either genomic information, UPG or pedigree truncation under GEBV accuracies ranged among the evaluated traits. Overall, methods using UPGs, truncation pedigree and genomic information exhibited potential to improve GEBV accuracies, bias and dispersion for all traits compared with other methods. Truncated scenarios promoted high genetic gains. In small populations with few genotyped animals, combining truncated pedigree or UPG with genomic information is a feasible approach to deal with missing pedigrees.

Genome-wide association and functional annotation analysis for the calving interval in Nellore cattle

Calving interval (CI) measures the number of days between two consecutive calves of the same cow, and previous studies based on phenotype and pedigree data reported low heritability for this trait. However, the genetic architecture of CI in the Nellore breed was not evaluated based on genomic data. Thus, this study aimed to estimate the heritability based on genomic data and carry out a genome-wide association study (GWAS) for CI in the Nellore breed, using 12,599 pedigree records, 5078 CI records, and 3818 animals genotyped with 50k SNPchip panel. Both quality control and GWAS were performed in BLUPF90 family packages, which use the single-step genomic best linear unbiased predictor (ssGBLUP) method. The average CI was 427.6 days, with a standard deviation of 106.9 and a total range of 270-730 days. The heritability estimate was 0.04 ± 0.04. The p-values of GWAS analysis resulted in a genomic inflation factor (lambda) of 1.08. The only significant SNP (rs136725686) at the genome-wide level (p-value = 1.53E-06) was located on BTA13. Other 19 SNPs were significant at the chromosome-wide level, distributed on BTA1, 2, 3, 6, 10, 13, 14, 17, 18, 22, and 26. Functional annotation analysis found thirty-six protein-coding genes, including genes related to cell cycle (RAD21, BCAR3), oocyte function (LHX8, CLPX, UTP23), immune system (TXK, TEC, NFATC2), endocrine function (LRRFIP2, GPR158), estrous cycle (SLC38A7), and female fertility (CCK, LYZL4, TRAK1, FOXP1, STAC). Therefore, CI is a complex trait with small heritability in Nellore cattle, and various biological processes may be involved with the genetic architecture of CI in Nellore cattle.

Genetic gains underpinning a little-known strawberry Green Revolution

The annual production of strawberry has increased by one million tonnes in the US and 8.4 million tonnes worldwide since 1960. Here we show that the US expansion was driven by genetic gains from Green Revolution breeding and production advances that increased yields by 2,755%. Using a California population with a century-long breeding history and phenotypes of hybrids observed in coastal California environments, we estimate that breeding has increased fruit yields by 2,974-6,636%, counts by 1,454-3,940%, weights by 228-504%, and firmness by 239-769%. Using genomic prediction approaches, we pinpoint the origin of the Green Revolution to the early 1950s and uncover significant increases in additive genetic variation caused by transgressive segregation and phenotypic diversification. Lastly, we show that the most consequential Green Revolution breeding breakthrough was the introduction of photoperiod-insensitive, PERPETUAL FLOWERING hybrids in the 1970s that doubled yields and drove the dramatic expansion of strawberry production in California.

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

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

Improving Breeding Value Reliability with Genomic Data in Breeding Groups of Charolais

The aim of this study was to assess the impact of incorporating genomic data using the single-step genomic best linear unbiased prediction (ssGBLUP) method compared to the best linear unbiased prediction (BLUP) method on the reliability of breeding values for age at first calving, calving interval, and productive longevity at 78 months in Charolais cattle. The study included 48,590 purebred Charolais individuals classified into four subgroups based on genotyping and performance records. The results showed that considering genotypes significantly improved genomic estimated breeding values (GEBV) reliability across all categories except nongenotyped individuals. For young genotyped individuals, the increase in reliability was up to 27% for both sexes. The highest average reliability was achieved for genotyped proven bulls and cows with performance records, and the inclusion of genomic data further improved the reliability by up to 22% and 21% for cows and bulls, respectively. The gain in reliability was observed mainly during the first three calvings, and then the differences decreased. The imported individuals showed lower estimated breeding values (EBV) and GEBV reliabilities than the domestic population, probably due to the weak genetic connection with the domestic population. However, when the progeny of imported heifers were sired by domestic bulls, the reliability increased by up to 24%. For nongenotyped individuals, only a slight increase in reliability was observed; however, the number of genotyped individuals in the population was still relatively small.

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.

Validating genomic prediction for nitrogen efficiency index and its composition traits of Holstein cows in early lactation

Nitrogen (N) use efficiency (NUE) is an economically important trait for dairy cows. Recently, we proposed a new N efficiency index (NEI), that simultaneously considers both NUE and N pollution. This study aimed to validate the genomic prediction for NEI and its composition traits and investigate the relationship between SNP effects estimated directly from NEI and indirectly from its composition traits. The NEI composition included genomic estimated breeding value of N intake (NINT), milk true protein N (MTPN) and milk urea N yield. The edited data were 132,899 records on 52,064 cows distributed in 773 herds. The pedigree contained 122,368 animals. Genotypic data of 566,294 SNP was available for 4514 individuals. A total of 4413 cows (including 181 genotyped) and 56 bulls (including 32 genotyped) were selected as the validation populations. The linear regression method was used to validate the genomic prediction of NEI and its composition traits using best linear unbiased prediction (BLUP) and single-step genomic BLUP (ssGBLUP). The mean theoretical accuracies of validation populations obtained from ssGBLUP were higher than those obtained from BLUP for both NEI and its composition traits, ranging from 0.57 (MTPN) to 0.72 (NINT). The highest mean prediction accuracies for NEI and its composition traits were observed for the genotyped cows estimated under ssGBLUP, ranging from 0.48 (MTPN) to 0.66 (NINT). Furthermore, the SNP effects estimated from NEI composition traits, multiplied by the relative weight were the same as those estimated directly from NEI. This study preliminary showed that genomic prediction can be used for NEI, however, we acknowledge the need for further validation of this result in a larger dataset. Moreover, the SNP effects of NEI can be indirectly calculated using the SNP effects estimated from its composition traits. This study provided a basis for adding genomic information to establish NEI as part of future routine genomic evaluation programs.

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 pre-selected variants from large-scale whole-genome sequence data for single-step genomic predictions in pigs

BACKGROUND

Whole-genome sequence (WGS) data harbor causative variants that may not be present in standard single nucleotide polymorphism (SNP) chip data. The objective of this study was to investigate the impact of using preselected variants from WGS for single-step genomic predictions in maternal and terminal pig lines with up to 1.8k sequenced and 104k sequence imputed animals per line.

METHODS

Two maternal and four terminal lines were investigated for eight and seven traits, respectively. The number of sequenced animals ranged from 1365 to 1491 for the maternal lines and 381 to 1865 for the terminal lines. Imputation to sequence occurred within each line for 66k to 76k animals for the maternal lines and 29k to 104k animals for the terminal lines. Two preselected SNP sets were generated based on a genome-wide association study (GWAS). Top40k included the SNPs with the lowest p-value in each of the 40k genomic windows, and ChipPlusSign included significant variants integrated into the porcine SNP chip used for routine genotyping. We compared the performance of single-step genomic predictions between using preselected SNP sets assuming equal or different variances and the standard porcine SNP chip.

RESULTS

In the maternal lines, ChipPlusSign and Top40k showed an average increase in accuracy of 0.6 and 4.9%, respectively, compared to the regular porcine SNP chip. The greatest increase was obtained with Top40k, particularly for fertility traits, for which the initial accuracy based on the standard SNP chip was low. However, in the terminal lines, Top40k resulted in an average loss of accuracy of 1%. ChipPlusSign provided a positive, although small, gain in accuracy (0.9%). Assigning different variances for the SNPs slightly improved accuracies when using variances obtained from BayesR. However, increases were inconsistent across the lines and traits.

CONCLUSIONS

The benefit of using sequence data depends on the line, the size of the genotyped population, and how the WGS variants are preselected. When WGS data are available on hundreds of thousands of animals, using sequence data presents an advantage but this remains limited in pigs.

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.

Genome-wide association study and functional annotation analyses for nitrogen efficiency index and its composition traits in dairy cattle

The aims of this study were (1) to identify genomic regions associated with a N efficiency index (NEI) and its composition traits and (2) to analyze the functional annotation of identified genomic regions. The NEI included N intake (NINT1), milk true protein N (MTPN1), milk urea N yield (MUNY1) in primiparous cattle, and N intake (NINT2+), milk true protein N (MTPN2+), and milk urea N yield (MUNY2+) in multiparous cattle (2 to 5 parities). The edited data included 1,043,171 records on 342,847 cows distributed in 1,931 herds. The pedigree consisted of 505,125 animals (17,797 males). Data of 565,049 SNPs were available for 6,998 animals included in the pedigree (5,251 females and 1,747 males). The SNP effects were estimated using a single-step genomic BLUP approach. The proportion of the total additive genetic variance explained by windows of 50 consecutive SNPs (with an average size of about 240 kb) was calculated. The top 3 genomic regions explaining the largest rate of the total additive genetic variance of the NEI and its composition traits were selected for candidate gene identification and quantitative trait loci (QTL) annotation. The selected genomic regions explained from 0.17% (MTPN2+) to 0.58% (NEI) of the total additive genetic variance. The largest explanatory genomic regions of NEI, NINT1, NINT2+, MTPN1, MTPN2+, MUNY1, and MUNY2+ were Bos taurus autosome 14 (1.52-2.09 Mb), 26 (9.24-9.66 Mb), 16 (75.41-75.51 Mb), 6 (8.73-88.92 Mb), 6 (8.73-88.92 Mb), 11 (103.26-103.41 Mb), 11 (103.26-103.41 Mb). Based on the literature, gene ontology, Kyoto Encyclopedia of Genes and Genomes, and protein-protein interaction, 16 key candidate genes were identified for NEI and its composition traits, which are mainly expressed in the milk cell, mammary, and liver tissues. The number of enriched QTL related to NEI, NINT1, NINT2+, MTPN1, and MTPN2+ were 41, 6, 4, 11, 36, 32, and 32, respectively, and most of them were related to the milk, health, and production classes. In conclusion, this study identified genomic regions associated with NEI and its composition traits, and identified key candidate genes describing the genetic mechanisms of N use efficiency-related traits. Furthermore, the NEI reflects not only its composition traits but also the interactions among them.

Testicular hypoplasia in Nellore Cattle: Genetic analysis and functional analysis of genome-wide association study results

Characterized by the incomplete development of the germinal epithelium of the seminiferous tubules, Testicular hypoplasia (TH) leads to decreased sperm concentration, increased morphological changes in sperm and azoospermia. Economic losses resulting from the disposal of affected bulls reduce the efficiency of meat production systems. A genome-wide association study and functional analysis were performed to identify genomic windows and the underlying positional candidate genes associated with TH in Nellore cattle. Phenotypic and pedigree data from 207,195 animals and genotypes (461,057 single nucleotide polymorphism, SNP) from 17,326 sires were used in this study. TH was evaluated as a binary trait measured at 18 months of age. A possible correlated response on TH resulting from the selection for scrotal circumference was evaluated by using a two-trait analysis. Thus, estimated breeding values were calculated by fitting a linear-threshold animal model in a Bayesian approach. The SNP effects were estimated using the weighted single-step genomic BLUP method. Twelve non-overlapping windows of 20 adjacent SNP that explained more than 1% of the additive genetic variance were selected for candidate gene annotation. Functional and gene prioritization analysis of the candidate genes identified genes (KHDRBS3, GPX5, STAR, ERLIN2), which might play an important role in the expression of TH due to their known roles in the spermatogenesis process, synthesis of steroids and lipid metabolism.

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

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

Mathematical Modeling and Software Tools for Breeding Value Estimation Based on Phenotypic, Pedigree and Genomic Information of Holstein Friesian Cattle in Serbia

In this paper, six univariate and two multivariate best linear unbiased prediction (BLUP) models were tested for the estimation of breeding values (BV) in Holstein Friesian cattle in Serbia. Two univariate models were formed using the numerator relationship matrix (NRM), four using the genomic relationship matrix (GRM). Multivariate models contained only an NRM. Two cases were studied, the first when only first lactations were observed, and the second when all lactations were observed using a repeatability model. A total of 6041 animals were included, and of them, 2565 had data on milk yield (MY), milk fat yield (FY), milk fat content (FC), milk protein yield (PY) and milk protein content (PC). Finally, out of those 2565 cows, 1491 were genotyped. A higher accuracy of BV was obtained when using a combination of NRM and GRM compared to NRM alone in univariate analysis, while multivariate analysis with repeated measures gave the highest accuracy with all 6041 animals. When only genotyped animals were observed, the highest accuracy of the estimated BV was calculated by the ssGBLUPp model, and the lowest by the univariate BLUP model. In conclusion, the current breeding programs in Serbia should be changed to use multivariate analysis with repeated measurements until the optimal size of the reference population, which must include genotyping data on both bulls and cows, is reached.

Genome-wide association study for milk production and conformation traits in Canadian Alpine and Saanen dairy goats

Increasing the productivity of Canadian dairy goats is critical to the competitiveness of the sector; however, little is known about the underlying genetic architecture of economically important traits in these populations. Consequently, the objectives of this study were as follows: (1) to perform a single-step GWAS for milk production traits (milk, protein, and fat yields, and protein and fat percentages in first and later lactations) and conformation traits (body capacity, dairy character, feet and legs, fore udder, general appearance, rear udder, suspensory ligament, and teats) in the Canadian Alpine and Saanen breeds; and (2) to identify positional and functional candidate genes related to these traits. The data available for analysis included 305-d milk production records for 6,409 Alpine and 3,434 Saanen does in first lactation and 5,827 Alpine and 2,632 Saanen does in later lactations; as well as linear type conformation records for 5,158 Alpine and 2,342 Saanen does. Genotypes were available for 833 Alpine and 874 Saanen animals. Both single-breed and multiple-breed GWAS were performed using single-trait animal models. Positional and functional candidate genes were then identified in downstream analyses. The GWAS identified 189 unique SNP that were significant at the chromosomal level, corresponding to 271 unique positional candidate genes within 50 kb up- and downstream, across breeds and traits. This study provides evidence for the economic importance of several candidate genes (e.g., CSN1S1, CSN2, CSN1S2, CSN3, DGAT1, and ZNF16) in the Canadian Alpine and Saanen populations that have been previously reported in other dairy goat populations. Moreover, several novel positional and functional candidate genes (e.g., RPL8, DCK, and MOB1B) were also identified. Overall, the results of this study have provided greater insight into the genetic architecture of milk production and conformation traits in the Canadian Alpine and Saanen populations. Greater understanding of these traits will help to improve dairy goat breeding programs.

Implication of the order of blending and tuning when computing the genomic relationship matrix in single-step GBLUP

Single-step genomic BLUP (ssGBLUP) relies on the combination of the genomic ( G $$ \mathbf{G} $$ ) and pedigree relationship matrices for all ( A $$ \mathbf{A} $$ ) and genotyped ( A 22 $$ {\mathbf{A}}_{22} $$ ) animals. The procedure ensures G $$ \mathbf{G} $$ and A 22 $$ {\mathbf{A}}_{22} $$ are compatible so that both matrices refer to the same genetic base ('tuning'). Then G $$ \mathbf{G} $$ is combined with a proportion of A 22 $$ {\mathbf{A}}_{22} $$ ('blending') to avoid singularity problems and to account for the polygenic component not accounted for by markers. This computational procedure has been implemented in the reverse order (blending before tuning) following the sequential research developments. However, blending before tuning may result in less optimal tuning because the blended matrix already contains a proportion of A 22 $$ {\mathbf{A}}_{22} $$ . In this study, the impact of 'tuning before blending' was compared with 'blending before tuning' on genomic estimated breeding values (GEBV), single nucleotide polymorphism (SNP) effects and indirect predictions (IP) from ssGBLUP using American Angus Association and Holstein Association USA, Inc. data. Two slightly different tuning methods were used; one that adjusts the mean diagonals and off-diagonals of G $$ \mathbf{G} $$ to be similar to those in A 22 $$ {\mathbf{A}}_{22} $$ and another one that adjusts based on the average difference between all elements of G $$ \mathbf{G} $$ and A 22 $$ {\mathbf{A}}_{22} $$ . Over 6 million Angus growth records and 5.9 million Holstein udder depth records were available. Genomic information was available on 51,478 Angus and 105,116 Holstein animals. Average realized relationship estimates among groups of animals were similar across scenarios. Scatterplots show that GEBV, SNP effects and IP did not noticeably change for all animals in the evaluation regardless of the order of computations and when using blending parameter of 0.05. Formulas were derived to determine the blending parameter that maximizes changes in the genomic relationship matrix and GEBV when changing the order of blending and tuning. Algebraically, the change is maximized when the blending parameter is equal to 0.5. Overall, tuning G $$ \mathbf{G} $$ before blending, regardless of blending parameter used, had a negligible impact on genomic predictions and SNP effects in this study.

Genome-Wide Association Study for Body Conformation Traits and Fitness in Czech Holsteins

The aim of this study was a genome-wide association study (GWAS) on conformation traits using 25,486 genotyped Czech Holsteins, with 35,227 common SNPs for each genotype. Linear trait records were collected between 1995 and 2020. The Interbull information from Multiple Across Country Evaluation (MACE) was included for bulls that mostly had daughter records in a foreign country. When using the Bonferroni correction, the number of SNPs that were either significant or approached the significance threshold was low-dairy capacity composite on BTA4, feet and legs composite BTA21, total score BTA10, stature BTA24, body depth BTA6, angularity BTA20, fore udder attachment BTA10. Without the Bonferroni correction, the total number of significant or near of significance SNPs was 32. The SNPs were localized on BTA1,2,4,5,6,7,8,18,22,25,26,28 for dairy capacity composite, BTA15,21 for feet and legs composite, BTA10 for total score, BTA24 stature, BTA6,23 body depth, BTA20 angularity, BTA2 rump angle, BTA9,10 rear legs rear view, BTA2,19 rear legs side view, BTA10 fore udder attachment, BTA2 udder depth, BTA10 rear udder height, BTA12 central alignment, BTA24 rear teat placement, BTA8,29 rear udder width. The results provide biological information for the improvement of body conformation and fitness in the Holstein population.

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

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

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

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

Genomic selection in tropical perennial crops and plantation trees: a review

To overcome the multiple challenges currently faced by agriculture, such as climate change and soil deterioration, more efficient plant breeding strategies are required. Genomic selection (GS) is crucial for the genetic improvement of quantitative traits, as it can increase selection intensity, shorten the generation interval, and improve selection accuracy for traits that are difficult to phenotype. Tropical perennial crops and plantation trees are of major economic importance and have consequently been the subject of many GS articles. In this review, we discuss the factors that affect GS accuracy (statistical models, linkage disequilibrium, information concerning markers, relatedness between training and target populations, the size of the training population, and trait heritability) and the genetic gain expected in these species. The impact of GS will be particularly strong in tropical perennial crops and plantation trees as they have long breeding cycles and constrained selection intensity. Future GS prospects are also discussed. High-throughput phenotyping will allow constructing of large training populations and implementing of phenomic selection. Optimized modeling is needed for longitudinal traits and multi-environment trials. The use of multi-omics, haploblocks, and structural variants will enable going beyond single-locus genotype data. Innovative statistical approaches, like artificial neural networks, are expected to efficiently handle the increasing amounts of heterogeneous multi-scale data. Targeted recombinations on sites identified from profiles of marker effects have the potential to further increase genetic gain. GS can also aid re-domestication and introgression breeding. Finally, GS consortia will play an important role in making the best of these opportunities.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-022-01326-4.

Factor Analysis of Genetic Parameters for Body Conformation Traits in Dual-Purpose Simmental Cattle

Simple Summary

Body conformation traits are closely related to economically important characteristics and should be considered in cattle breeding programs. A variety of body conformation traits recorded by classifiers can complicate the analysis process. Factor analysis can reduce the number of variables by combining two or more variables into a single factor, which has biological significance. The results of this study could be used by breeders to define conformation indexes and implement genetic assessments for conformation traits in dual-purpose breeds.

Abstract

In this study, we estimated the genetic parameters for 6 composite traits and 27 body conformation traits of 1016 dual-purpose Simmental cattle reared in northwestern China from 2010 to 2019 using a linear animal mixed model. To integrate these traits, a variety of methods were used as follows: (1) genetic parameters estimates for composite and individual body conformation traits based on the pedigree relationship matrix (A) and combined genomic-pedigree relationship matrix (H); (2) factor analysis to explore the relationships among body conformation traits; and (3) genetic parameters of factor scores estimated using A and H, and the correlations of EBVs of the factor scores and EBVs of the composite traits. Heritability estimates of the composite traits using A and H were low to medium (0.07–0.47). The 24 common latent factors explained 96.13% of the total variance. Among factors with eigenvalues ≥ 1, F1 was mainly related to body frame, muscularity, and rump; F2 was related to feet and legs; F3, F4, F5, and F6 were related to teat placement, teat size, udder size, and udder conformation; and F7 was related to body frame. Single-trait analysis of factor scores yielded heritability estimates that were low to moderate (0.008–0.43 based on A and 0.04–0.43 based on H). Spearman and Pearson correlations, derived from the best linear unbiased prediction analysis of composite traits and factor scores, showed a similar pattern. Thus, incorporating factor analysis into the morphological evaluation to simplify the assessment of body conformation traits may improve the genetics of dual-purpose Simmental cattle.

Genomic prediction and genetic correlations estimated for milk production and fatty acid traits in Walloon Holstein cattle using random regression models

The aims of this study were to: (1) estimate genetic correlation for milk production traits (milk, fat and protein yields and fat and protein contents) and fatty acids (FA: C16:0, C18:1 cis-9, LCFA, SFA, and UFA) over days in milk, (2) investigate the performance of genomic predictions using single-step GBLUP (ssGBLUP) based on random regression models (RRM), and (3) identify the optimal scaling and weighting factors to be used in the construction of the H matrix. A total of 302 684 test-day records of 63.875 first lactation Walloon Holstein cows were used. Positive genetic correlations were found between milk yield and fat and protein yield (rg from 0.46 to 0.85) and between fat yield and milk FA (rg from 0.17 to 0.47). On the other hand, negative correlations were estimated between fat and protein contents (rg from -0.22 to -0.59), between milk yield and milk FA (rg from -0.22 to -0.62), and between protein yield and milk FA (rg from -0.11 to -0.19). The selection for high fat content increases milk FA throughout lactation (rg from 0.61 to 0.98). The test-day ssGBLUP approach showed considerably higher prediction reliability than the parent average for all milk production and FA traits, even when no scaling and weighting factors were used in the H matrix. The highest validation reliabilities (r2 from 0.09 to 0.38) and less biased predictions (b1 from 0.76 to 0.92) were obtained using the optimal parameters (i.e., ω = 0.7 and α = 0.6) for the genomic evaluation of milk production traits. For milk FA, the optimal parameters were ω = 0.6 and α = 0.6. However, biased predictions were still observed (b1 from 0.32 to 0.81). The findings suggest that using ssGBLUP based on RRM is feasible for the genomic prediction of daily milk production and FA traits in Walloon Holstein dairy cattle.

Impact of blending the genomic relationship matrix with different levels of pedigree relationships or the identity matrix on genetic evaluations

Evaluations using single-step genomic BLUP require blending the genomic relationship matrix (G) with a positive definite matrix to ensure nonsingularity for solving the mixed model equations. Many organizations blend G with a proportion of the numerator relationship matrix for genotyped animals (A ₂₂) to improve stability and possibly add a residual polygenic effect. However, when nearly all the polygenic variance is explained by G, blending with A ₂₂ may cause inflation and add excess computing time; thus, blending with an identity matrix (I) multiplied by a small value may be a better solution. The objective of this study was to evaluate changes in reliability and inflation of genomic estimated breeding values, convergence rate, elapsed wall-clock time for blending G with different levels of A ₂₂ or I, and develop a more time-efficient blending method. A US Holstein cattle data set was used with 9.7 million animals in the pedigree, 569,404 animals with genotypes, and 10.1 million stature phenotypes. Blending G by adding a small value to the diagonal elements had comparable performance to A ₂₂ with fewer rounds to convergence required to solve the system of equations. Reliability and inflation of genomic estimated breeding values ranged from 0.63 to 0.68 and 0.86 to 0.89 for all blending scenarios tested. The current blending default in the BLUPF90 software is to replace G with (1 - β)G + βA ₂₂, where β equals 0.05. In this study, β values of 0.30, 0.20, 0.05, 0.01, 0.005, and 0.001 were evaluated with A ₂₂ and I. Negligible differences in elapsed computing time between the blending types and levels were observed. Subsequently, the current blending algorithm used in the BLUPF90 family of programs was optimized, reducing the blending time from approximately 2 h to 5 min for A ₂₂ and less than 1 s for I. The new time difference between blending with A ₂₂ or I is negligible and not computationally critical. The results indicate that blending G with A ₂₂ does not have clear advantages over blending with a small proportion of I.

Low-density SNP panel for efficient imputation and genomic selection of milk production and technological traits in dairy sheep

The present study aimed to ascertain how different strategies for leveraging genomic information enhance the accuracy of estimated breeding values for milk and cheese-making traits and to evaluate the implementation of a low-density (LowD) SNP chip designed explicitly for that aim. Thus, milk samples from a total of 2,020 dairy ewes from 2 breeds (1,039 Spanish Assaf and 981 Churra) were collected and analyzed to determine 3 milk production and composition traits and 2 traits related to milk coagulation properties and cheese yield. The 2 studied populations were genotyped with a customized 50K Affymetrix SNP chip (Affymetrix Inc.) containing 55,627 SNP markers. The prediction accuracies were obtained using different multitrait methodologies, such as the BLUP model based on pedigree information, the genomic BLUP (GBLUP), and the BLUP at the SNP level (SNP-BLUP), which are based on genotypic data, and the single-step GBLUP (ssGBLUP), which combines both sources of information. All of these methods were analyzed by cross-validation, comparing predictions of the whole population with the test population sets. Additionally, we describe the design of a LowD SNP chip (3K) and its prediction accuracies through the different methods mentioned previously. Furthermore, the results obtained using the LowD SNP chip were compared with those based on the 50K SNP chip data sets. Finally, we conclude that implementing genomic selection through the ssGBLUP model in the current breeding programs would increase the accuracy of the estimated breeding values compared with the BLUP methodology in the Assaf (from 0.19 to 0.39) and Churra (from 0.27 to 0.44) dairy sheep populations. The LowD SNP chip is cost-effective and has proven to be an accurate tool for estimating genomic breeding values for milk and cheese-making traits, microsatellite imputation, and parentage verification. The results presented here suggest that the routine use of this LowD SNP chip could potentially increase the genetic gains of the breeding selection programs of the 2 Spanish dairy sheep breeds considered here.

Weighted Single-Step Genomic Best Linear Unbiased Prediction Method Application for Assessing Pigs on Meat Productivity and Reproduction Traits

Changes in the accuracy of the genomic estimates obtained by the ssGBLUP and wssGBLUP methods were evaluated using different reference groups. The weighting procedure’s reasonableness of application Pwas considered to improve the accuracy of genomic predictions for meat, fattening and reproduction traits in pigs. Six reference groups were formed to assess the genomic data quantity impact on the accuracy of predicted values (groups of genotyped animals). The datasets included 62,927 records of meat and fattening productivity (fat thickness over 6–7 ribs (BF1, mm)), muscle depth (MD, mm) and precocity up to 100 kg (age, days) and 16,070 observations of reproductive qualities (the number of all born piglets (TNB) and the number of live-born piglets (NBA), according to the results of the first farrowing). The wssGBLUP method has an advantage over ssGBLUP in terms of estimation reliability. When using a small reference group, the difference in the accuracy of ssGBLUP over BLUP AM is from −1.9 to +7.3 percent points, while for wssGBLUP, the change in accuracy varies from +18.2 to +87.3 percent points. Furthermore, the superiority of the wssGBLUP is also maintained for the largest group of genotyped animals: from +4.7 to +15.9 percent points for ssGBLUP and from +21.1 to +90.5 percent points for wssGBLUP. However, for all analyzed traits, the number of markers explaining 5% of genetic variability varied from 71 to 108, and the number of such SNPs varied depending on the size of the reference group (79–88 for BF1, 72–81 for MD, 71–108 for age). The results of the genetic variation distribution have the greatest similarity between groups of about 1000 and about 1500 individuals. Thus, the size of the reference group of more than 1000 individuals gives more stable results for the estimation based on the wssGBLUP method, while using the reference group of 500 individuals can lead to distorted results of GEBV.

Status quo of genetic improvement in local goats: a review

This review aims to summarize and synthesize the fragmented information available on the genetic improvement of local goats (criollo, indigenous, native) on the American and other continents, where populations with these goats have an important role in food security and the economy of rural communities, as well as in conservation of biodiversity and productivity improvement. Topics such as the current state of goat production globally, conservation programs, resistance to parasites and diseases, use of phenotypical characteristics and genomic information, and molecular markers for genetic improvement are addressed. The main challenges, opportunities, and limitations described in recent literature concerning local goats in the immediate future are discussed.

Using genomic prediction to detect microevolutionary change of a quantitative trait

Detecting microevolutionary responses to natural selection by observing temporal changes in individual breeding values is challenging. The collection of suitable datasets can take many years and disentangling the contributions of the environment and genetics to phenotypic change is not trivial. Furthermore, pedigree-based methods of obtaining individual breeding values have known biases. Here, we apply a genomic prediction approach to estimate breeding values of adult weight in a 35-year dataset of Soay sheep (Ovis aries). Comparisons are made with a traditional pedigree-based approach. During the study period, adult body weight decreased, but the underlying genetic component of body weight increased, at a rate that is unlikely to be attributable to genetic drift. Thus cryptic microevolution of greater adult body weight has probably occurred. Genomic and pedigree-based approaches gave largely consistent results. Thus, using genomic prediction to study microevolution in wild populations can remove the requirement for pedigree data, potentially opening up new study systems for similar research.

Genomic Prediction of Complex Traits in Perennial Plants: A Case for Forest Trees

This chapter provides an overview of the genomic selection progress in long-lived forest tree species. Factors affecting the prediction accuracy in genomic prediction are assessed with examples from empirical studies. Infrastructure and resources required for the implementation of genomic selection are evaluated. Some general guidelines are provided for the successful application of genomic selection in forest tree breeding programs.

Dissection of the practical soybean breeding pipeline by developing ZDX1, a high-throughput functional array

KEY MESSAGE

We developed the ZDX1 high-throughput functional soybean array for high accuracy evaluation and selection of both parents and progeny, which can greatly accelerate soybean breeding. Microarray technology facilitates rapid, accurate, and economical genotyping. Here, using resequencing data from 2214 representative soybean accessions, we developed the high-throughput functional array ZDX1, containing 158,959 SNPs, covering 90.92% of soybean genes and sites related to important traits. By application of the array, a total of 817 accessions were genotyped, including three subpopulations of candidate parental lines, parental lines and their progeny from practical breeding. The fixed SNPs were identified in progeny, indicating artificial selection during the breeding process. By identifying functional sites of target traits, novel soybean cyst nematode-resistant progeny and maturity-related novel sources were identified by allele combinations, demonstrating that functional sites provide an efficient method for the rapid screening of desirable traits or gene sources. Notably, we found that the breeding index (BI) was a good indicator for progeny selection. Superior progeny were derived from the combination of distantly related parents, with at least one parent having a higher BI. Furthermore, new combinations based on good performance were proposed for further breeding after excluding redundant and closely related parents. Genomic best linear unbiased prediction (GBLUP) analysis was the best analysis method and achieved the highest accuracy in predicting four traits when comparing SNPs in genic regions rather than whole genomic or intergenic SNPs. The prediction accuracy was improved by 32.1% by using progeny to expand the training population. Collectively, a versatile assay demonstrated that the functional ZDX1 array provided efficient information for the design and optimization of a breeding pipeline for accelerated soybean breeding.

Improving lodgepole pine genomic evaluation using spatial correlation structure and SNP selection with single-step GBLUP

Modeling environmental spatial heterogeneity can improve the efficiency of forest tree genomic evaluation. Furthermore, genotyping costs can be lowered by reducing the number of markers needed. We investigated the impact on variance components, breeding value accuracy, and bias of two phenotypic data adjustments (experimental design and autoregressive spatial models), and a relationship matrix calculated from a subset of markers selected for their ability to infer ancestry. Using a multiple-trait multiple-site single-step Genomic Best Linear Unbiased Prediction (ssGBLUP) approach, four scenarios (2 phenotype adjustments × 2 marker sets) were applied to diameter at breast height (DBH), height (HT), and resistance to western gall rust (WGR) in four open-pollinated progeny trials of lodgepole pine, with 1490 (out of 11,188) trees genotyped with 25,099 SNPs. As a control, we fitted the conventional ABLUP model using pedigree information. The highest heritability estimates were achieved for the ABLUP followed closely by the ssGBLUP with the full marker set and using the spatial phenotype adjustments. The highest predictive ability was obtained by using a reduced marker subset (8000 SNPs) when either the spatial (DBH: 0.429, and WGR: 0.513) or design (HT: 0.467) phenotype corrections were used. No significant difference was detected in prediction bias among the six fitted models, and all values were close to 1 (0.918-1.014). Results demonstrated that selecting informative markers, such as those capturing ancestry, can improve the predictive ability. The use of spatial correlation structure increased traits' heritability and reduced prediction bias, while increases in predictive ability were trait-dependent.

First-degree relationships and genotyping errors deciphered by a high-density SNP array in a Duroc × Iberian pig cross

BACKGROUND

Two individuals with a first-degree relationship share about 50 percent of their alleles. Parent-offspring relationships cannot be homozygous for alternative alleles (genetic exclusion).

METHODS

Applying the concept of genetic exclusion to HD arrays typed in animals for experimental purposes or genomic selection allows estimation of the rate of rejection of first-degree relationships as the rate at which two individuals typed for a large number of Single Nucleotide Polymorphisms (SNPs) do not share at least one allele. An Expectation-Maximization algorithm is applied to estimate parentage. In addition, genotyping errors are estimated in true parent-offspring relationships. Samples from nine candidate Duroc sires and 55 Iberian dams producing 214 Duroc × Iberian barrows were typed for the HD porcine Affymetrix array.

RESULTS

We were able to establish paternity and maternity of 75 and 85 piglets, respectively. Rate of rejection in true parent-offspring relationships was estimated as 0.000735. This is a lower bound of the genotyping error since rate of rejection depends on allele frequencies. After accounting for allele frequencies, our estimate of the genotyping error is 0.6%. A total of 7,744 SNPs were rejected in five or more true parent-offspring relationships facilitating identification of "problematic" SNPs with inconsistent inheritance in multiple parent-offspring relationships.

CONCLUSIONS

This study shows that animal experiments and routine genotyping in genomic selection allow to establish or to verify first-degree relationships as well as to estimate genotyping errors for each batch of animals or experiment.

Prediction ability for growth and maternal traits using SNP arrays based on different marker densities in Nellore cattle using the ssGBLUP

This study aimed to investigate the prediction ability for growth and maternal traits using different low-density customized SNP arrays selected by informativeness and distribution of markers across the genome employing single-step genomic BLUP (ssGBLUP). Phenotypic records for adjusted weight at 210 and 450 days of age were utilized. A total of 945 animals were genotyped with high-density chip, and 267 individuals born after 2008 were selected as validation population. We evaluated 11 scenarios using five customized density arrays (40 k, 20 k, 10 k, 5 k and 2 k) and the HD array was used as desirable scenario. The GEBV predictions and BIF (Beef Improvement Federation) accuracy were obtained with BLUPF90 family programs. Linear regression was used to evaluate the prediction ability, inflation, and bias of GEBV of each customized array. An overestimation of partial GEBVs in contrast with complete GEBVs and increase of BIF accuracy with the density arrays diminished were observed. For all traits, the prediction ability was higher as the array density increased and it was similar with customized arrays higher than 10 k SNPs. Level of inflation was lower as the density array increased of and was higher for MW210 effect. The bias was susceptible to overestimation of GEBVs when the density customized arrays decreased. These results revealed that the BIF accuracy is sensible to overestimation using low-density customized arrays while the prediction ability with least 10,000 informative SNPs obtained from the Illumina BovineHD BeadChip shows accurate and less biased predictions. Low-density customized arrays under ssGBLUP method could be feasible and cost-effective in genomic selection.

Determining Heat Stress Effects of Multiple Genetic Traits in Tropical Dairy Cattle Using Single-Step Genomic BLUP

Heat stress is becoming a significant problem in dairy farming, especially in tropical countries, making accurate genetic selection for heat tolerance a priority. This study investigated the effect of heat stress manifestation on genetics for milk yield, milk quality, and dairy health traits with and without genomic information using single-step genomic best linear unbiased prediction (ssGBLUP) and BLUP in Thai−Holstein crossbred cows. The dataset contained 104,150 test-day records from the first lactation of 15,380 Thai−Holstein crossbred cows. A multiple-trait random regression test-day model on a temperature−humidity index (THI) function was used to estimate the genetic parameters and genetic values. Heat stress started at a THI of 76, and the heritability estimates ranged from moderate to low. The genetic correlation between those traits and heat stress in both BLUP methods was negative. The accuracy of genomic predictions in the ssGBLUP method was higher than the BLUP method. In conclusion, heat stress negatively impacted milk production, increased the somatic cell score, and disrupted the energy balance. Therefore, in dairy cattle genetic improvement programs, heat tolerance is an important trait. The new genetic evaluation method (ssGBLUP) should replace the traditional method (BLUP) for more accurate genetic selection.

Single-step genomic evaluation of milk production traits in Canadian Alpine and Saanen dairy goats

Genomic evaluations are routine in most plant and livestock breeding programs but are used infrequently in dairy goat breeding schemes. In this context, the purpose of this study was to investigate the use of the single-step genomic BLUP method for predicting genomic breeding values for milk production traits (milk, protein, and fat yields; protein and fat percentages) in Canadian Alpine and Saanen dairy goats. There were 6,409 and 12,236 Alpine records and 3,434 and 5,008 Saanen records for each trait in first and later lactations, respectively, and a total of 1,707 genotyped animals (833 Alpine and 874 Saanen). Two validation approaches were used, forward validation (i.e., animals born after 2013 with an average estimated breeding value accuracy from the full data set ≥0.50) and forward cross-validation (i.e., subsets of all animals included in the forward validation were used in successive replications). The forward cross-validation approach resulted in similar validation accuracies (0.55 to 0.66 versus 0.54 to 0.61) and biases (-0.01 to -0.07 versus -0.03 to 0.11) to the forward validation when averaged across traits. Additionally, both single and multiple-breed analyses were compared, and similar average accuracies and biases were observed across traits. However, there was a small gain in accuracy from the use of multiple-breed models for the Saanen breed. A small gain in validation accuracy for genomically enhanced estimated breeding values (GEBV) relative to pedigree-based estimated breeding values (EBV) was observed across traits for the Alpine breed, but not for the Saanen breed, possibly due to limitations in the validation design, heritability of the traits evaluated, and size of the training populations. Trait-specific gains in theoretical accuracy of GEBV relative to EBV for the validation animals ranged from 17 to 31% in Alpine and 35 to 55% in Saanen, using the cross-validation approach. The GEBV predicted from the full data set were 12 to 16% more accurate than EBV for genotyped animals, but no gains were observed for nongenotyped animals. The largest gains were found for does without lactation records (35-41%) and bucks without daughter records (46-54%), and consequently, the implementation of genomic selection in the Canadian dairy goat population would be expected to increase selection accuracy for young breeding candidates. Overall, this study represents the first step toward implementation of genomic selection in Canadian dairy goat populations.

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

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

Additive and Dominance Genomic Analysis for Litter Size in Purebred and Crossbred Iberian Pigs

INGA FOOD S. A., as a Spanish company that produces and commercializes fattened pigs, has produced a hybrid Iberian sow called CASTÚA by crossing the Retinto and Entrepelado varieties. The selection of the parental populations is based on selection criteria calculated from purebred information, under the assumption that the genetic correlation between purebred and crossbred performance is high; however, these correlations can be less than one because of a GxE interaction or the presence of non-additive genetic effects. This study estimated the additive and dominance variances of the purebred and crossbred populations for litter size, and calculated the additive genetic correlations between the purebred and crossbred performances. The dataset consisted of 2030 litters from the Entrepelado population, 1977 litters from the Retinto population, and 1958 litters from the crossbred population. The individuals were genotyped with a GeneSeek^® GGP Porcine70K HDchip. The model of analysis was a ‘biological’ multivariate mixed model that included additive and dominance SNP effects. The estimates of the additive genotypic variance for the total number born (TNB) were 0.248, 0.282 and 0.546 for the Entrepelado, Retinto and Crossbred populations, respectively. The estimates of the dominance genotypic variances were 0.177, 0.172 and 0.262 for the Entrepelado, Retinto and Crossbred populations. The results for the number born alive (NBA) were similar. The genetic correlations between the purebred and crossbred performance for TNB and NBA—between the brackets—were 0.663 in the Entrepelado and 0.881 in Retinto poplulations. After backsolving to obtain estimates of the SNP effects, the additive genetic variance associated with genomic regions containing 30 SNPs was estimated, and we identified four genomic regions that each explained > 2% of the additive genetic variance in chromosomes (SSC) 6, 8 and 12: one region in SSC6, two regions in SSC8, and one region in SSC12.

A Comprehensive Comparison of Haplotype-Based Single-Step Genomic Predictions in Livestock Populations With Different Genetic Diversity Levels: A Simulation Study

The level of genetic diversity in a population is inversely proportional to the linkage disequilibrium (LD) between individual single nucleotide polymorphisms (SNPs) and quantitative trait loci (QTLs), leading to lower predictive ability of genomic breeding values (GEBVs) in high genetically diverse populations. Haplotype-based predictions could outperform individual SNP predictions by better capturing the LD between SNP and QTL. Therefore, we aimed to evaluate the accuracy and bias of individual-SNP- and haplotype-based genomic predictions under the single-step-genomic best linear unbiased prediction (ssGBLUP) approach in genetically diverse populations. We simulated purebred and composite sheep populations using literature parameters for moderate and low heritability traits. The haplotypes were created based on LD thresholds of 0.1, 0.3, and 0.6. Pseudo-SNPs from unique haplotype alleles were used to create the genomic relationship matrix ( G ) in the ssGBLUP analyses. Alternative scenarios were compared in which the pseudo-SNPs were combined with non-LD clustered SNPs, only pseudo-SNPs, or haplotypes fitted in a second G (two relationship matrices). The GEBV accuracies for the moderate heritability-trait scenarios fitting individual SNPs ranged from 0.41 to 0.55 and with haplotypes from 0.17 to 0.54 in the most (Ne ≅ 450) and less (Ne < 200) genetically diverse populations, respectively, and the bias fitting individual SNPs or haplotypes ranged between -0.14 and -0.08 and from -0.62 to -0.08, respectively. For the low heritability-trait scenarios, the GEBV accuracies fitting individual SNPs ranged from 0.24 to 0.32, and for fitting haplotypes, it ranged from 0.11 to 0.32 in the more (Ne   ≅ 250) and less (Ne   ≅ 100) genetically diverse populations, respectively, and the bias ranged between -0.36 and -0.32 and from -0.78 to -0.33 fitting individual SNPs or haplotypes, respectively. The lowest accuracies and largest biases were observed fitting only pseudo-SNPs from blocks constructed with an LD threshold of 0.3 (p < 0.05), whereas the best results were obtained using only SNPs or the combination of independent SNPs and pseudo-SNPs in one or two G matrices, in both heritability levels and all populations regardless of the level of genetic diversity. In summary, haplotype-based models did not improve the performance of genomic predictions in genetically diverse populations.

Genome‑wide association study and genomic prediction for growth traits in yellow-plumage chicken using genotyping-by-sequencing

Background

Growth traits are of great importance for poultry breeding and production and have been the topic of extensive investigation, with many quantitative trait loci (QTL) detected. However, due to their complex genetic background, few causative genes have been confirmed and the underlying molecular mechanisms remain unclear, thus limiting our understanding of QTL and their potential use for the genetic improvement of poultry. Therefore, deciphering the genetic architecture is a promising avenue for optimising genomic prediction strategies and exploiting genomic information for commercial breeding. The objectives of this study were to: (1) conduct a genome-wide association study to identify key genetic factors and explore the polygenicity of chicken growth traits; (2) investigate the efficiency of genomic prediction in broilers; and (3) evaluate genomic predictions that harness genomic features.

Results

We identified five significant QTL, including one on chromosome 4 with major effects and four on chromosomes 1, 2, 17, and 27 with minor effects, accounting for 14.5 to 34.1% and 0.2 to 2.6% of the genomic additive genetic variance, respectively, and 23.3 to 46.7% and 0.6 to 4.5% of the observed predictive accuracy of breeding values, respectively. Further analysis showed that the QTL with minor effects collectively had a considerable influence, reflecting the polygenicity of the genetic background. The accuracy of genomic best linear unbiased predictions (BLUP) was improved by 22.0 to 70.3% compared to that of the conventional pedigree-based BLUP model. The genomic feature BLUP model further improved the observed prediction accuracy by 13.8 to 15.2% compared to the genomic BLUP model.

Conclusions

A major QTL and four minor QTL were identified for growth traits; the remaining variance was due to QTL effects that were too small to be detected. The genomic BLUP and genomic feature BLUP models yielded considerably higher prediction accuracy compared to the pedigree-based BLUP model. This study revealed the polygenicity of growth traits in yellow-plumage chickens and demonstrated that the predictive ability can be greatly improved by using genomic information and related features.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12711-021-00672-9.

Genome-wide association study on milk production and somatic cell score for Thai dairy cattle using weighted single-step approach with random regression test-day model

Genome-wide association studies are a powerful tool to identify genomic regions and variants associated with phenotypes. However, only limited mutual confirmation from different studies is available. The objectives of this study were to identify genomic regions as well as genes and pathways associated with the first-lactation milk, fat, protein, and total solid yields; fat, protein, and total solid percentage; and somatic cell score (SCS) in a Thai dairy cattle population. Effects of SNPs were estimated by a weighted single-step GWAS, which back-solved the genomic breeding values predicted using single-step genomic BLUP (ssGBLUP) fitting a single-trait random regression test-day model. Genomic regions that explained at least 0.5% of the total genetic variance were selected for further analyses of candidate genes. Despite the small number of genotyped animals, genomic predictions led to an improvement in the accuracy over the traditional BLUP. Genomic predictions using weighted ssGBLUP were slightly better than the ssGBLUP. The genomic regions associated with milk production traits contained 210 candidate genes on 19 chromosomes [Bos taurus autosome (BTA) 1 to 7, 9, 11 to 16, 20 to 21, 26 to 27 and 29], whereas 21 candidate genes on 3 chromosomes (BTA 11, 16, and 21) were associated with SCS. Many genomic regions explained a small fraction of the genetic variance, indicating polygenic inheritance of the studied traits. Several candidate genes coincided with previous reports for milk production traits in Holstein cattle, especially a large region of genes on BTA14. We identified 141 and 5 novel genes related to milk production and SCS, respectively. These novel genes were also found to be functionally related to heat tolerance (e.g., SLC45A2, IRAG1, and LOC101902172), longevity (e.g., SYT10 and LOC101903327), and fertility (e.g., PAG1). These findings may be attributed to indirect selection in our population. Identified biological networks including intracellular cell transportation and protein catabolism implicate milk production, whereas the immunological pathways such as lymphocyte activation are closely related to SCS. Further studies are required to validate our findings before exploiting them in genomic selection.