Estimated Breeding Values Research Articles

Genetic parameters for reproductive traits in the Italian Mediterranean buffalo using milk yield as a correlated trait

Until now, the genetic evaluation of the Italian Mediterranean Buffalo has been mainly focused on production traits. However, female fertility affects the efficiency of the dairy industry as it is essential to maintain the profitability of dairy farms. Indeed, the estimation of its genetic component is crucial for its improvement. In this study, 3 measures of buffalo's fertility were analyzed: the age at first calving (AFC), the interval between first and second calving (CI1), and the interval between second and successive calvings (CI2_12). Milk yield at 270 d (MY270) was used as a correlated trait. First, genetic parameters were estimated using 7,915 buffalo cows with first calving from 1991 to 2018, then breeding values were calculated from 236,087 buffalo cows. Genetic parameters were estimated by Bayesian inference fitting a multiple-trait animal model using the GIBBS1F90 program while BLUPF90 was used for estimation of breeding value. The heritability and repeatability estimates of fertility traits were low. The genetic correlations among fertility traits ranged from 0.10 (AFC-CI1) to 0.92 (CI1-CI2_12). Genetic correlation between MY270 and fertility traits was unfavorable, ranging from 0.23 to 0.48. The results from this study can be used as a basis for the future genetic improvement of fertility traits in the Italian Mediterranean Buffaloes.

Genome-wide association study identifies functional genomic variants associated with young stock survival in Nordic Red Dairy Cattle

Identifying quantitative trait loci (QTL) associated with calf survival is essential for both reducing economic loss in cattle industry and understanding the genetic basis of the trait. To identify mutations and genes underlying young stock survival (YSS), we performed GWAS using de-regressed estimated breeding values of a YSS index and its component traits defined by sex and age in 3,077 Nordic Red Dairy Cattle (RDC) bulls and 2 stillbirth traits (first lactation and later lactations) in 5,141 RDC bulls. Two associated QTL regions on Bos taurus autosome (BTA) 4 and 6 were identified for the YSS index. The results of 4 YSS component traits indicate that same QTL regions were associated with bull and heifer calf mortality, but the effects were different over the growing period and suggested an additional QTL on BTA23. The GWAS on stillbirth identified 3 additional QTL regions on BTA5, 14, and 24 compared with YSS and its component traits. The conditional test of BTA6 showed at least 2 closely located QTL segregating for YSS component traits and stillbirth. We found 2 independent QTL for stillbirth on BTA23. The post-GWAS revealed LCORL, PPM1K, SSP1, MED28, and LAP3 are putative causal genes on BTA6, and a frame shift variant within LCORL, BTA6:37401770 (rs384548488) could be the putative causal variant. On BTA4, the GRB10 gene is the putative causal gene and BTA4:5296018 is the putative causal variant. In addition, NDUFA9 and FGF23 on BTA5, LYN on BTA14, and KCNK5 on BTA23 are putative causal genes for QTL for stillbirth. The gene analysis also proposed several candidate genes. Our findings shed new light on the candidate genes affecting calf survival, and the knowledge could be utilized to reduce calf mortality and thereby enhance welfare of dairy cattle.

Accurate haplotype construction and detection of selection signatures enabled by high quality pig genome sequences

High-quality whole-genome resequencing in large-scale pig populations with pedigree structure and multiple breeds would enable accurate construction of haplotype and robust selection-signature detection. Here, we sequence 740 pigs, combine with 149 of our previously published resequencing data, retrieve 207 resequencing datasets, and form a panel of worldwide distributed wild boars, aboriginal and highly selected pigs with pedigree structures, amounting to 1096 genomes from 43 breeds. Combining with their haplotype-informative reads and pedigree structure, we accurately construct a panel of 1874 haploid genomes with 41,964,356 genetic variants. We further demonstrate its valuable applications in GWAS by identifying five novel loci for intramuscular fat content, and in genomic selection by increasing the accuracy of estimated breeding value by 36.7%. In evolutionary selection, we detect MUC13 gene under a long-term balancing selection, as well as NPR3 gene under positive selection for pig stature. Our study provides abundant genomic variations for robust selection-signature detection and accurate haplotypes for deciphering complex traits in pigs.

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

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

Genome-wide association study as a powerful tool for dissecting competitive traits in legumes.

Legumes are extremely valuable because of their high protein content and several other nutritional components. The major challenge lies in maintaining the quantity and quality of protein and other nutritional compounds in view of climate change conditions. The global need for plant-based proteins has increased the demand for seeds with a high protein content that includes essential amino acids. Genome-wide association studies (GWAS) have evolved as a standard approach in agricultural genetics for examining such intricate characters. Recent development in machine learning methods shows promising applications for dimensionality reduction, which is a major challenge in GWAS. With the advancement in biotechnology, sequencing, and bioinformatics tools, estimation of linkage disequilibrium (LD) based associations between a genome-wide collection of single-nucleotide polymorphisms (SNPs) and desired phenotypic traits has become accessible. The markers from GWAS could be utilized for genomic selection (GS) to predict superior lines by calculating genomic estimated breeding values (GEBVs). For prediction accuracy, an assortment of statistical models could be utilized, such as ridge regression best linear unbiased prediction (rrBLUP), genomic best linear unbiased predictor (gBLUP), Bayesian, and random forest (RF). Both naturally diverse germplasm panels and family-based breeding populations can be used for association mapping based on the nature of the breeding system (inbred or outbred) in the plant species. MAGIC, MCILs, RIAILs, NAM, and ROAM are being used for association mapping in several crops. Several modifications of NAM, such as doubled haploid NAM (DH-NAM), backcross NAM (BC-NAM), and advanced backcross NAM (AB-NAM), have also been used in crops like rice, wheat, maize, barley mustard, etc. for reliable marker-trait associations (MTAs), phenotyping accuracy is equally important as genotyping. Highthroughput genotyping, phenomics, and computational techniques have advanced during the past few years, making it possible to explore such enormous datasets. Each population has unique virtues and flaws at the genomics and phenomics levels, which will be covered in more detail in this review study. The current investigation includes utilizing elite breeding lines as association mapping population, optimizing the choice of GWAS selection, population size, and hurdles in phenotyping, and statistical methods which will analyze competitive traits in legume breeding.

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.

Insights from the early generations of the Swedish rainbow trout (Oncorhynchus mykiss) breeding program

ABSTRACT This study aimed to gain insights about the status of the Swedish breeding program through studying key phenotypic traits. In total, 133 and 73 full-sib families were formed during the reproductive seasons of 2016 and 2019, respectively. Growth-related recordings were available from two occasions: ∼ 9 and 24–27 months post-hatch. A preliminary descriptive analysis of the fecundity and early embryo survival identified substantial differences in favor of G0 (year class 2016) which was partly explained by the fact that first-time 3-year-old spawners were used in G1 (year class 2019). Moderate to high heritability values (0.23–0.49) were obtained for body length and weight across both time points and generations. Furthermore, the prediction accuracy of the estimated breeding values was ∼0.65, while the expected genetic gain was 36 g per year. Overall, our results demonstrate the positive performance of the Swedish rainbow trout breeding program.

Marker-assisted and Genomic Selection of Beef Cattle

This article provides an overview of modern genetic technologies for improving production traits and predicting breeding value in beef cattle. In particular, in marker-assisted selection the most promising is the selectionby desirable genotypes in the genes of myostatin (MSTN), calpain (CAPN), calpastatin (CAST), growth hormone (GH), leptin (LEP), thyroglobulin (TG), fatty acid binding protein (FABP), retinoic acid C-receptor (RORC), diacyl-glycerol acyltransferase (DGATI), sterol-Co desaturase (SCD). A modern and much more advanced approach is the Single Step Genomic Best Linear Unbiased Predictions (ssGBLUP) method, which calculates a Genomic Estimated Breeding Value (GEBV) using DNA chip genotyping, phenotype and pedigree data. Genome-wide association studies (GWAS), based on the use of genetic markers distributed throughout the genome and in non-equilibrium linkage with at least one of the quantitative traits, are currently recognised as more informative for finding new genes for beef cattle productivity. New genes associated with live weight at different stages of ontogenesis, average daily live weight gain, residual feed intake, carcass weight and flesh content have been identified. Most of the identified genes control cell division, lipid and carbohydrate metabolism. The accumulated data on full-genome association studies and exome sequencing led to new methods of genetic analysis – gene ontology and gene networks. The use of gene networks provided the first detailed understanding of the genetic basis for the formation of complex phenotypic traits based on the complex interaction of regulatory networks of «major» and «peripheral» genes controlling the development of a particular trait.

Health assessment of Holstein calves born after in vitro fertilization, biopsy-based genotyping at the blastocyst stage and subsequent embryo transfer

Establishing methods for evaluating genomic estimated breeding values of bovine embryos can potentially increase the efficiency of breeding programs by transferring only embryos with a high genomic estimated breeding value. This may be achieved by analyzing DNA from trophectoderm biopsies. However, manipulation of bovine embryos is associated with a risk of impaired conceptus health. More knowledge on the health implications of embryonic handling procedures is required. In this study, we followed pregnancies after transfer of in vitro-produced (IVP) embryos and assessed the health of the offspring during the first 2 weeks of life. Three groups of calves were studied: i) freshly transferred non-biopsied embryos (39 transfers, 17 calves; Group B-/C-); ii) biopsied and freshly transferred IVP embryos (42 transfers, 21 calves; Group B+/C-); iii) biopsied and cryopreserved IVP embryos (17 transfers, 6 calves; Group B+/C+). Blood biochemical and hematologic values were compared between groups and to a control group of 13 calves produced by conventional artificial insemination. The pregnancy rate on day 50 and the calving rate did not differ among the groups, but the average gestation length of the B+/C+ group was significantly shorter and with wider variation than the two other groups. There was a tendency toward a higher average body weight at birth in group B+/C+ (45.1 kg) and the standard deviation in body weight was larger (11.7 kg) compared to the B-/C- (39.5 kg; 3.2 kg) and B+/C- (41.8 kg; 6 kg) groups. Body weight on day 14 was higher in the B+/C+ calves compared to the other groups. There was no difference in the biochemical and hematological values at birth between the groups and these were within the normal range. However, when compared to a group of calves produced by standard artificial insemination, significantly higher concentrations were found for the hepatic-related enzymes ALAT, ASAT, ALP, and GGT in group B-/C-and B+/C-, while only higher ALP concentrations were found in B+/C+ calves. The biochemical findings indicate higher heterogeneity in IVP calves compared to calves produced by artificial insemination. The more manipulated IVP embryos also showed increased heterogeneity in body weight at birth, with a shift toward heavier calves, which calls for closer attendance at parturition to handle dystocia in a timely manner and minimize fetal losses.

A statistical package for evaluation of hybrid performance in plant breeding via genomic selection

Hybrid breeding employs heterosis, which could potentially improve the yield and quality of a crop. Genomic selection (GS) is a promising approach for the selection of quantitative traits in plant breeding. The main objectives of this study are to (i) propose a GS-based approach to identify potential parental lines and superior hybrid combinations from a breeding population, which is composed of hybrids produced by a half diallel mating design; (ii) develop a software package for users to carry out the proposed approach. An R package, designated EHPGS, was generated to facilitate the employment of the genomic best linear unbiased model considering additive plus dominance marker effects for the hybrid performance evaluation. The R package contains a Bayesian statistical algorithm for calculating genomic estimated breeding value (GEBVs), GEBV-based specific combining ability, general combining ability, mid-parent heterosis, and better-parent heterosis. Three datasets that have been published in literature, including pumpkin (Cucurbita maxima), maize (Zea mays), and wheat (Triticum aestivum L.), were reanalyzed to illustrate the use of EHPGS.

Using social interaction models for genetic analysis of skin damage in gilts.

Skin damage is a trait of economic and welfare importance that results from social interactions between animals. These interactions may produce wound signs on the gilt's skin as a result of damage behavior (i.e., fighting), biting syndromes (i.e., tail, vulva, or ear biting), and swine inflammation and necrosis syndrome. Although current selection for traits that are affected by social interactions primarily focuses on improving direct genetic effects, combined selection on direct and social genetic effects could increase genetic gain and avoid a negative response to selection in cases of competitive behavior. The objectives of this study were to (1) estimate variance components for combined skin damage (CSD), with or without accounting for social genetic effects, (2) investigate the impact of including genomic information on the prediction accuracy, bias, and dispersion of CSD estimated breeding values, and (3) perform a single-step genome-wide association study (ssGWAS) of CSD under a classical and a social interaction model. Our results show that CSD is heritable and affected by social genetic effects. Modeling CSD with social interaction models increased the total heritable variance relative to the phenotypic variance by three-fold compared to the classical model. Including genomic information increased the prediction accuracy of direct, social, and total estimated breeding values for purebred sires by at least 21.2%. Bias and dispersion of estimated breeding values were reduced by including genomic information in classical and social interaction models but remained present. The ssGWAS did not identify any single nucleotide polymorphism that was significantly associated with social or direct genetic effects for CSD. Combined skin damage is heritable, and genetic selection against this trait will increase the welfare of animals in the long term. Combined skin damage is affected by social genetic effects, and modeling this trait with a social interaction model increases the potential for genetic improvement. Including genomic information increases the prediction accuracy of estimated breeding values and reduces their bias and dispersion, although some biases persist. The results of the genome-wide association study indicate that CSD has a polygenic architecture and no major quantitative trait locus was detected.

Genetic and phenotypic time trends of litter size, piglet mortality, and birth weight in pigs

IntroductionLitter size in pigs has increased steadily since 1990. Because of unfavorable genetic correlations with piglet mortality, breeding goals should include survival traits next to litter size. Unbalanced breeding programs that neglect this requirement have produced increased mortality levels, attracting negative public attention. Balanced breeding does not have this disadvantage, but the general public is largely unaware of this.MethodsWe present long-term time trends as realized in commercial breeding. The data includes (i) phenotypes of litter size, piglet birth weight, and piglet mortality, as used in routine breeding value estimation; and (ii) the genomic Best Linear Unbiased Prediction (gBLUP) estimated breeding values thus obtained. Piglet mortality (2001–2022) and birth weight (2009–2022) phenotypes were related to litter size by recording year. Estimated breeding values (EBVs) for the mortality traits were regressed on those for litter size by birth year (2012–2022).ResultsAverage litter size is very weakly correlated to the mortality (R2 ≤ 0.06) and birth weight (0.07 ≤ R2 ≤ 0.26) traits, and those correlations are unfavorable (antagonistic) within each year. However, all traits analyzed here show favorable simultaneous phenotypic and genetic trends over time: the antagonisms are neutralized by balanced breeding. Above the annual mean litter size level, farrowing and lactation mortality rates increased with increasing litter size in every year (unfavorable), but the annual intercepts and the slopes decreased from 2001 to 2022 (favorable). Average litter birth weight decreased with litter size in every year (unfavorable), but the annual intercepts increased and the slopes decreased from 2009 to 2022 (favorable). The within-litter birth weight variation coefficient increased with litter size in every year (unfavorable), but the annual intercepts decreased from 2009 to 2022 (favorable). The proportion of low birth weights (i.e.,< 0.9 kg) for a given litter size is decreasing over time, and the critical birth weight level (below which lactation mortality increases strongly) is clearly population dependent and changes over time too.DiscussionThe increases in litter size and piglet survival rates due to balanced breeding policies lead to reduced total numbers of dead piglets (i.e., per country, per year) coinciding with a certain pig production volume (i.e., with a certain total number of weaned piglets).

Genotype by Prenatal Environment Interaction for Postnatal Growth of Nelore Beef Cattle Raised under Tropical Grazing Conditions.

The prenatal environment is recognized as crucial for the postnatal performance in cattle. In tropical regions, pregnant beef cows commonly experience nutritional restriction during the second half of the gestation period. Thus, the present study was designed to analyze the genotype by prenatal environment interaction (G × Epn) and to identify genomic regions associated with the level and response in growth and reproduction-related traits of beef cattle to changes in the prenatal environment. A reaction norm model was applied to data from two Nelore herds using the solutions of contemporary groups for birth weight as a descriptor variable of the gestational environment quality. A better gestational environment favored weights until weaning, scrotal circumference at yearling, and days to first calving of the offspring. The G × Epn was strong enough to result in heterogeneity of variance components and genetic parameters in addition to reranking of estimated breeding values and SNPs effects. Several genomic regions associated with the level of performance and specific responses of the animals to variations in the gestational environment were revealed, which harbor QTLs and can be exploited for selection purposes. Therefore, genetic evaluation models considering G × Epn and special management and nutrition care for pregnant cows are recommended.

Heterogeneity of variance and genetic parameters for milk production in cattle, using Bayesian inference.

The goal of this study was to verify the effect of heterogeneity of variance (HV) on milk production in up to 305 days of lactation (L305) of daughters of Girolando, Gir and Holstein sires, as well as in the genetic evaluation of these sires and their progenies. in Brazil. The model included contemporary groups (consisting of herd, year and calving season) as a fixed effect, cow age at calving (linear and quadratic effects) and heterozygosity (linear effect) as covariates, in addition to the random effects of direct additive genetic and environmental, permanent and residual. The first analysis consisted of the single-trait animal model, with L305 records (disregarding HV). The second considered classes of standard deviations (SD): two-trait model including low and high classes (considering HV), according to the standardized means of L305 for herd-year of calving. The low SD class was composed of herds with SD equal to or less than zero and the high class with positive SD values. Estimates of (co)variance components and breeding values were obtained separately for each scenario using Bayesian inference via Gibbs sampling. Different heritability was estimated. Higher for the high DP class in the Gir (0.20) and Holstein (0.15) breeds, not occurring the same in the Girolando breed, with a lower value among the classes for the high DP (0.10). High values of genetic correlations were also found between low and high SD classes (0.88; 0.85 and 0.79) for the Girolando, Gir and Holstein breeds, respectively. Like the order correlations (Spearman) which were also high for the three breeds analyzed (equal to or above 0.92). Thus, the presence of HV had a smaller impact for L305 and did not affect the genetic evaluation of sires.

Efficient Genomic Prediction of Yield and Dry Matter in Hybrid Potato.

There is an ongoing endeavor within the potato breeding sector to rapidly adapt potato from a clonal polyploid crop to a diploid hybrid potato crop. While hybrid breeding allows for the efficient generation and selection of parental lines, it also increases breeding program complexity and results in longer breeding cycles. Over the past two decades, genomic prediction has revolutionized hybrid crop breeding through shorter breeding cycles, lower phenotyping costs, and better population improvement, resulting in increased genetic gains for genetically complex traits. In order to accelerate the genetic gains in hybrid potato, the proper implementation of genomic prediction is a crucial milestone in the rapid improvement of this crop. The authors of this paper set out to test genomic prediction in hybrid potato using current genotyped material with two alternative models: one model that predicts the general combining ability effects (GCA) and another which predicts both the general and specific combining ability effects (GCA+SCA). Using a training set comprising 769 hybrids and 456 genotyped parental lines, we found that reasonable a prediction accuracy could be achieved for most phenotypes with both zero common parents (ρ=0.36-0.61) and one (ρ=0.50-0.68) common parent between the training and test sets. There was no benefit with the inclusion of non-additive genetic effects in the GCA+SCA model despite SCA variance contributing between 9% and 19% of the total genetic variance. Genotype-by-environment interactions, while present, did not appear to affect the prediction accuracy, though prediction errors did vary across the trial's targets. These results suggest that genomically estimated breeding values on parental lines are sufficient for hybrid yield prediction.

Association of Phenotypic Markers of Heat Tolerance with Australian Genomic Estimated Breeding Values and Dairy Cattle Selection Indices.

Dairy cattle predicted by genomic breeding values to be heat tolerant are known to have less milk production decline and lower core body temperature increases in response to elevated temperatures. In a study conducted at the University of Melbourne's Dookie Robotic Dairy Farm during summer, we identified the 20 most heat-susceptible and heat-tolerant cows in a herd of 150 Holstein Friesian lactating cows based on their phenotypic responses (changes in respiration rate, surface body temperature, panting score, and milk production). Hair samples were collected from the tip of the cows' tails following standard genotyping protocols. The results indicated variation in feed saved and HT genomic estimated breeding values (GEBVs) (p ≤ 0.05) across age, indicating a potential for their selection. As expected, the thermotolerant group had higher GEBVs for HT and feed saved but lower values for milk production. In general, younger cows had superior GEBVs for the Balanced Performance Index (BPI) and Australian Selection Index (ASI), whilst older cows were superior in fertility, feed saved (FS), and HT. This study demonstrated highly significant (p ≤ 0.001) negative correlations (-0.28 to -0.74) between HT and GEBVs for current Australian dairy cattle selection indices (BPI, ASI, HWI) and significant (p ≤ 0.05) positive correlations between HT and GEBVs for traits like FS (0.45) and fertility (0.25). Genomic selection for HT will help improve cow efficiency and sustainability of dairy production under hot summer conditions. However, a more extensive study involving more lactating cows across multiple farms is recommended to confirm the associations between the phenotypic predictors of HT and GEBVs.

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

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

AGHmatrix: genetic relationship matrices in R.

The resemble between relatives computed from pedigree and genomic data is an important resource for geneticists and ecologists, who are interested in understanding how genes influence phenotypic variation, fitness adaptation, and population dynamics. The AGHmatrix software is an R-package focused on the construction of pedigree (A matrix) and/or molecular markers (G matrix), with the possibility of building a combined matrix of pedigree corrected by molecular markers (H matrix). Designed to estimate the relationships for any ploidy level, the software also includes auxiliary functions related to filtering molecular markers, compute linkage disequilibrium (LD) and check pedigree errors in large data sets. After computing the relationship matrices, results from the AGHmatrix can be used in different contexts, including on prediction of (genomic) estimated breeding values and genome-wide association studies (GWAS). AGHmatrix v2.1.0 is available under GPL-3 license in CRAN at https://cran.r-project.org/web/packages/AGHmatrix/index.html and also in GitHub at https://github.com/rramadeu/AGHmatrix. It has a comprehensive tutorial, and it follows with real data examples. Supplementary data are available at Bioinformatics online.

Comparative genetic analysis of frequentist and Bayesian approach for reproduction, production and life time traits showing favourable association of age at first calving in Tharparkar cattle.

The present study was aimed primarily for estimating various genetic parameters (heritability, genetic correlations) of reproduction (age at first calving [AFC], first service period [FSP]); production (first lactation milk, solid-not fat, and fat yield) and lifetime traits (lifetime milk yield, productive life [PL], herd life [HL]) in Tharparkar cattle to check the association of reproduction traits with lifetime traits through two different methods (Frequentist and Bayesian) for comparative purpose. Animal breeding data of Tharparkar cattle (n = 964) collected from Livestock farm unit of ICAR-NDRI Karnal for the period 1990 through 2019 were analyzed using a Frequentist least squares maximum likelihood method (LSML; Harvey, 1990) and a multitrait Bayesian-Gibbs sampler approach (MTGSAM) for genetic correlations estimation of all the traits. Estimated breeding values of sires was obtained by BLUP and Bayesian analysis for the production traits. Heritability estimates of most of the traits were medium to high with the LSML (0.20±0.44 to 0.49±0.71) and Bayesian approach (0.24±0.009 to 0.61±0.017), respectively. However, more reliable estimates were obtained using the Bayesian technique. A higher heritability estimate was obtained for AFC (0.61±0.017) followed by first lactation fat yield, first lactation solid-not fat yield, FSP, first lactation milk yield (FLMY), PL (0.60±0.013, 0.60±0.006, 0.57±0.024, 0.57±0.020, 0.42±0.025); while a lower estimate for HL (0.38±0.034) by MTGSAM approach. Genetic and phenotypic correlations were negative for AFC-PL, AFC-HL, FSP-PL, and FSP-HL (-0.59±0.19, -0.59±0.24, -0.38±0.101 and -0.34±0.076) by the multi-trait Bayesian analysis. Breed and traits of economic importance are important for selection decisions to ensure genetic gain in cattle breeding programs. Favourable genetic and phenotypic correlations of AFC with production and lifetime traits compared to that of FSP indicated better scope of AFC for indirect selection of life-time traits at an early age. This also indicated that the present Tharparkar cattle herd had sufficient genetic diversity through the selection of AFC for the improvement of first lactation production and lifetime traits.

Comparison of various models for the estimation of heritability and breeding values.

The exploration of additive genetic variance for the selection of animals is the central paradigm in quantitative genetics and it is important to use appropriate animal models considering important factors. This study compares various factor effects for heritability and breeding values estimations on data collected on the Corriedale. Overall, the heritability estimates were the highest for birthweight (BW). Both the additive genetic variance and the heritability values were reduced upon adding more random effects. Year-effect as random reduced the heritability except for BW. The highest average accuracy for breeding values for birth, weaning, 6, 9, and 12-month-weights for were 0.656, 0.567, 0.519, 0.486, and 0.515. These were observed for models having sex or sex and parity as fixed effects. It is concluded that fixed and random factors have a significant impact on genetic-factors and breeding-value accuracies and should be considered during selection programs and make the selection program more effective.