Additive Genetic Correlations Research Articles

Structural Equation Models to Infer Relationships between Energy-related Blood Metabolites and Milk Daily Energy Output in Holstein Cows.

During lactation, high-yielding cows experience metabolic disturbances due to milk production. Metabolic monitoring offers valuable insights into how cows manage these challenges throughout the lactation period, making it a topic of considerable interest to breeders. In this study, we used Bayesian networks to uncover potential dependencies among various energy-related blood metabolites, i.e. glucose, urea, Beta-hydroxybutyrate (BHB), non-esterified fatty acids (NEFA), cholesterol (CHOL), and daily milk energy output (dMEO) in 1,254 Holstein cows. The inferred causal structure was then incorporated into structural equation models (SEM) to estimate heritabilities and additive genetic correlations among these phenotypes using both pedigree and genotypes from a 100k chip. Dependencies among traits were determined using the Hill Climbing algorithm, implemented with the posterior distribution of the residuals obtained from standard multiple trait model (MTM). These identified relationships were then used to construct the SEM, considering both direct and indirect relationships. The relevant dependencies and path coefficients obtained, expressed in units of measurement variation of 1σ, were as follows: dMEO → CHOL (0.181), dMEO → BHB (-0.149), dMEO → urea (0.038), glucose → BHB (-0.55), glucose → urea (-0.194), CHOL → urea (0.175), BHB → urea (-0.049), and NEFA → urea (-0.097). Heritabilities for traits of concern obtained with SEM ranged from 0.09 to 0.2. Genetic correlations with a minimum 95 % probability (P) of the posterior mean being > 0 for positive means or < 0 for negative means include those between dMEO and glucose (-0.583, P = 100), dMEO and BHB (0.349, P = 99), glucose and CHOL (0.325, P = 100), glucose and NEFA (-0.388, P = 100), and NEFA and BHB (0.759, P = 100). The results of this analysis revealed the existence of recursive relationships among the energy-related blood metabolites and dMEO. Understanding these connections is paramount for establishing effective genetic selection strategies, enhancing production and animal welfare.

The Genetic Architecture of Recombination Rates is Polygenic and Differs Between the Sexes in Wild House Sparrows (Passer domesticus).

Meiotic recombination through chromosomal crossing-over is a fundamental feature of sex and an important driver of genomic diversity. It ensures proper disjunction, allows increased selection responses, and prevents mutation accumulation; however, it is also mutagenic and can break up favorable haplotypes. This cost-benefit dynamic is likely to vary depending on mechanistic and evolutionary contexts, and indeed, recombination rates show huge variation in nature. Identifying the genetic architecture of this variation is key to understanding its causes and consequences. Here, we investigate individual recombination rate variation in wild house sparrows (Passer domesticus). We integrate genomic and pedigree data to identify autosomal crossover counts (ACCs) and intrachromosomal allelic shuffling (r¯intra) in 13,056 gametes transmitted from 2,653 individuals to their offspring. Females had 1.37 times higher ACC, and 1.55 times higher r¯intra than males. ACC and r¯intra were heritable in females and males (ACC h2 = 0.23 and 0.11; r¯intra h2 = 0.12 and 0.14), but cross-sex additive genetic correlations were low (rA = 0.29 and 0.32 for ACC and r¯intra). Conditional bivariate analyses showed that all measures remained heritable after accounting for genetic values in the opposite sex, indicating that sex-specific ACC and r¯intra can evolve somewhat independently. Genome-wide models showed that ACC and r¯intra are polygenic and driven by many small-effect loci, many of which are likely to act in trans as global recombination modifiers. Our findings show that recombination rates of females and males can have different evolutionary potential in wild birds, providing a compelling mechanism for the evolution of sexual dimorphism in recombination.

Parent of origin genetic effects on milk production traits in a population of Iranian Holstein cows.

The aim was to estimate the relative contribution of imprinting effects from both paternal and maternal sides to phenotypic variation in milk production traits including 305 days milk yield (MY), average daily milk production (ADM), fat percentage (F%), protein percentage (P%), 305 days fat yield (FY), 305 days protein yield (PY), ratio of fat percentage to protein percentage (F:P) and somatic cell score (SCS) in Iranian Holstein cows. To do this, each trait was analysed with a series of four animal models, which were identical for fixed and additive genetic effects but differed for combinations of paternal and maternal imprinting effects. The log-likelihood ratio test (LRT) and Akaike's information criteria (AIC) were used to select the best model for each trait. Correlations between traits due to additive and imprinting effects were estimated by bivariate analyses. For all traits studied, fitting the imprinting effect led to a better data fit. Also, it resulted in a noticeable decrease in additive genetic variance from 8% (SCS) to 28% (F:P). A significant maternal imprinting effect was detected on all traits studied. Estimates of maternal imprinting heritability ( ) were 0.07 ± 0.02, 0.04 ± 0.01, 0.06 ± 0.01, 0.05 ± 0.01, 0.5 ± 0.01, 0.09 ± 0.02, 0.07 ± 0.02 and 0.06 ± 0.01 for MY, ADM, F%, P%, FY, PY, F:P and SCS, respectively. For F:P, in addition to the maternal imprinting effect, a significant paternal imprinting component was also detected with a 7% contribution to phenotypic variance of F:P. Estimates of direct heritability ( ) were 0.29 ± 0.02, 0.17 ± 0.01, 0.22 ± 0.02, 0.11 ± 0.01, 0.18 ± 0.02, 0.22 ± 0.02, 0.15 ± 0.04 and 0.06 ± 0.01 for MY, ADM, F%, P%, FY, PY, F:P and SCS, respectively. Maternal imprinting correlations (rmi) were in a wide range between -0.75 ± 0.15 (P%-SCS) and 0.95 ± 0.11 (MY-ADM). Additive genetic correlations (ra) ranged between -0.54 ± 0.05 (MY-P%) and 0.99 ± 0.01 (MY-ADM) and phenotypic correlations (rp) ranged from -0.30 ± 0.01 (MY-F%) to 0.93 ± 0.01 (MY-ADM). The Spearman's correlation between additive breeding values including and excluding imprinting effects deviated from unity especially for top-ranked animals implying re-ranking of top animals following the inclusion of imprinting effects in the model. Since including imprinting effects in the model resulted in better data fit and re-ranking of top animals, including these effects in the genetic evaluation models for milk production traits was recommended.

Genetic parameter estimates for growth traits in Texel sheep of Brazil

The present study aimed to estimate the genetic parameters for direct and maternal influences on birth weight (BW), weaning weight (WW) and postweaning weight (PW) of Texel sheep, raised in an extensive system from Brazil. Data file used in this study consisted of 834 body weight records of 834 Texel lambs (376 males and 458 females) born from 2020 to 2021, offspring of 475 ewes and 45 rams. Multivariate animal model was conducted to estimate the variance components, heritability and genetic correlations using the restricted maximum likelihood method in ASReml-R software, while the effects of fixed factors were analyzed by the general linear model procedure of statistical analysis system. Five different animal models, with or without maternal genetic effect, maternal permanent environment effect and covariance between additive genetic and maternal genetic effects, were fitted on the growth traits data of the Texel sheep population. AIC and BIC test results were used to evaluate the best fit models for growth traits. The effect of lamb's sex was significant on BW, WW and PW while birth type and year of lambing effects were on all traits (P < 0.05). Model 4, incorporating direct additive genetic, maternal additive genetic and maternal permanent environmental effects showed the minimum AIC estimate for all growth traits. The direct additive genetic heritability for BW, WW and PW were 0.11 ± 0.03, 0.37 ± 0.05 and 0.31 ± 0.08, respectively and maternal heritability for these traits were 0.00 ± 0.00, 0.04 ± 0.02 and 0.00 ± 0.00, respectively. The total heritability estimates for BW, WW and PW were 0.11 ± 0.03, 0.39 ± 0.05 and 0.31 ± 0.08, respectively. The estimates of the direct additive genetic correlation between BW-WW, BW-PW and WW-PW were 0.68 ± 0.12, 0.66 ± 0.12 and 0.78 ± 0.07, respectively. The maternal permanent environmental correlation between BW-WW, BW-PW and WW-PW were 0.38 ± 0.27, 0.29 ± 0.20 and 0.74 ± 0.16, respectively. Phenotypic correlation ranged from 0.43 ± 0.03 to 0.75 ± 0.02. It is concluded that due to the existence of genetic variation for BB, WW and PW and strong genetic correlation among them, genetic progress for these traits is possible by selection in Texel sheep from Brazil.

Plateau-linear reaction norm model analysis of number born alive in purebred Landrace pigs using meteorological data in Japan.

We performed a plateau-linear reaction norm model (RNM) analysis of number born alive (NBA) in purebred Landrace pigs, where breeding value changes according to maximum temperature at mating day, using public meteorological observation data in Japan. We analysed 52,668 NBA records obtained from 10,320 Landrace sows. Pedigree data contained 99,201 animals. Off-farm daily temperature data at the nearest weather station from each of the farms were downloaded from the Japan Meteorological Agency website. A plateau-linear RNM analysis based on daily maximum temperature on mating day (threshold temperature of 16.6°C) was performed. The percentage of the records with daily maximum temperatures at mating days of ≤16.6, ≥25.0, ≥30.0 and ≥35.0°C were 34.3%, 33.6%, 14.0% and 0.8%, respectively. The value of Akaike's information criterion for the plateau-linear RNM was lower than that for a simple repeatability model (RM). With the plateau-linear RNM, estimated value of heritability ranged from 0.14 to 0.15, while that from the RM analysis was 0.15. Additive genetic correlation between intercept and slope terms was estimated to be -0.52 from the plateau-linear RNM analysis. Estimated additive genetic correlations were >0.9 between NBA at different temperatures ranging from 16.6 to 37.6°C. For the 10,320 sows, average values of prediction reliability of the intercept and slope terms for breeding values in the plateau-linear RNM were 0.47 and 0.16, respectively. Increasing weight for slope term in linear selection index could bring positive genetic gain in the slope part, but prediction accuracy would decrease. Our results imply that genetically improving heat tolerance in sows reared in Japan focusing on NBA using RNM is possible, while RNM is more complex to implement and interpret. Therefore, further study should be encouraged to make genetic improvement for heat tolerance in sows more efficient.

Estimation of Additive and Dominance Effects in an Acacia crassicarpa Multi-Environment Progeny Trial Using Genomic Pedigree Reconstruction

Abstract Acacia crassicarpa is an important tree species in Southeast Asia, where hundreds of thousands of hectares of planted forests are supported by advancements in silviculture and genetic improvement. Although possible, controlled pollination is impractical for advancing breeding populations, requiring an unreasonable effort to produce more than a few crosses per year. For this reason, breeding populations often are bred by open pollination. This study used large-scale pedigree reconstruction in multi-environment trials to assess full-sib families to model the genetics of the quantitative traits survival, straightness, height, diameter at breast height, tree volume, mean annual increment (MAI), and basic density. The traits were predominantly controlled by additive effects, with heritabilities between 0.09 for survival and 0.45 for basic density. The genetic correlation across sites was high for all traits, showing the low impact of genotype-by-environment interaction. The trait-trait correlation showed that straightness was independent of any other traits, survival was only correlated with MAI, and growth traits were highly correlated among themselves. Basic density was positively correlated with growth traits and MAI. Study Implications: Parentage analysis using an informative single nucleotide polymorphism panel was used to reconstruct pedigree and allow a full-sib family model to estimate additive and dominance effects and genetic correlations across sites and among important traits in an open-pollinated population. The genetic control of all traits assessed in this study was mainly additive. In this scenario, the recommended breeding strategy is forward selection of outstanding progeny for advanced generation breeding and backward selection of outstanding parents to produce seed for deployment via family forestry. Full-sib families can be identified by pedigree reconstruction at a seedling stage, followed by tissue culture multiplication, rooted cutting propagation, and plantation establishment.

Delineating maternal influence in regulation of variance in major economic traits of White Leghorns: Bayesian insights.

Proper variance partitioning and estimation of genetic parameters at appropriate time interval is crucial for understanding the dynamics of trait variance and genetic correlations and for deciding the future breeding strategy of the population. This study was conducted on the same premise to estimate genetic parameters of major economic traits in a White Leghorn strain IWH using Bayesian approach and to identify the role of maternal effects in the regulation of trait variance. Three different models incorporating the direct additive effect (Model 1), direct additive and maternal genetic effect (Model 2) and direct additive, maternal genetic and maternal permanent environmental effects (Model 3) were tried to estimate the genetic parameters for body weight traits (birth weight, body weight at 16, 20, 40 and 52 weeks), Age at sexual maturity (ASM), egg production traits (egg production up to 24, 28, 40, 52, 64 and 72 weeks) and egg weight traits (egg weight at 28, 40 and 52 weeks). Model 2 and Model 3 with maternal effects were found to be the best having the highest accuracy for almost all the traits. The direct additive genetic heritability was moderate for ASM, moderate to high for body weight traits and egg weight traits and low to moderate for egg production traits. Though the maternal heritability (h2mat) and permanent environmental effect (c2mpe) was low (<0.1) for most of the traits, they formed an important component of trait variance. Traits like egg weight at 28 weeks (0.14±0.06) and egg production at 72 weeks (0.13±0.07) reported comparatively higher values for c2mpe and h2mat respectively. Additive genetic correlation was high and positive between body weight traits, between egg weight traits, between consecutive egg production traits and between body weight and egg weight traits. However, a negative genetic correlation existed between egg production and egg weight traits, egg production and body weight traits, ASM and early egg production traits. Overall, a moderate positive genetic correlation was estimated between ASM and body weight traits and ASM and egg weight traits. Based on our findings, we can deduce that maternal effects constitute an important source of variation for all the major economic traits in White Leghorn and should be necessarily considered in genetic evaluation programs.

PSXII-5 Genetic Parameters for Dry Matter Intake in Beef Cattle as a Function of Water Restriction: Insights into Environmental Sensitivity and Genetic-By-Environment Interactions

Abstract Genetic selection for animals more tolerant or less sensitive to inadequate water availability may become critical to beef production, particularly as water resources become scarce or represent a meaningful cost to producers. The objectives of this study were 1) to explore the genetic component and environmental sensitivity (ES) of dry matter intake (DMI) to water restriction (WR), and 2) to evaluate genetic-by-environment interactions (G×E) with respect to WR using a reaction norm. A group of 788 crossbred steers arranged across seven cohorts were provided feed and water whereby daily water intake (WI) and DMI values were recorded automatically on an individual basis. After an initial acclimation period, data were collected for 70 days in which ad libitum water was provided and baseline WI was determined on each steer. Next, steers were restricted by approximately 10% of their baseline WI every seven days until 50% WR was achieved and subsequently maintained for 42 days. The reaction norm was modeled with a linear random regression allowing for additive genetic animal-specific intercepts and slope parameters as a function of the WR environmental covariate. The relationship covariance matrix for the additive genetic animal parameters was estimated using genotypes. The model was fit using the THR3GIBBSf90 program in the BLUPf90 software. The additive genetic variances for animal-specific intercepts and WR-driven slopes, as well as their covariance, were estimated. The estimated (posterior median and 95% highest posterior density) ES to WR, represented by the regression coefficient on the WR covariate, was -0.410 (-0.422, -0.399) kg daily DMI per 10% WR, indicating evidence for population-level DMI sensitivity to WR. The estimated additive genetic correlation between animal-specific intercept (i.e., ad libitum WI) and animal-specific slope (i.e., sensitivity to WR) was -0.75 (-0.88, -0.59), supporting an association between selection for increased DMI at ad libitum WI and ES. The additive genetic correlation between DMI at 0% and 25% WR was estimated at 0.96 (0.94, 0.99) and at 0% and 50% WR was 0.78 (0.65, 0.91), thereby indicating G×E had limited influence on the overall EBV ranking of animals as WR increased. The results of this study do not indicate major G×E concerns in DMI with respect to WR in the overall population; however, ES to WR may be important to producers as selection on DMI at ad libitum WI could increase ES.

Genotype-by-socioeconomic status interaction influences heart disease risk scores and carotid artery thickness in Mexican Americans: the predominant role of education in comparison to household income and socioeconomic index.

Background: Socioeconomic status (SES) is a potent environmental determinant of health. To our knowledge, no assessment of genotype-environment interaction has been conducted to consider the joint effects of socioeconomic status and genetics on risk for cardiovascular disease (CVD). We analyzed Mexican American Family Studies (MAFS) data to evaluate the hypothesis that genotype-by-environment interaction (GxE) is an important determinant of variation in CVD risk factors. Methods: We employed a linear mixed model to investigate GxE in Mexican American extended families. We studied two proxies for CVD [Pooled Cohort Equation Risk Scores/Framingham Risk Scores (FRS/PCRS) and carotid artery intima-media thickness (CA-IMT)] in relation to socioeconomic status as determined by Duncan's Socioeconomic Index (SEI), years of education, and household income. Results: We calculated heritability for FRS/PCRS and carotid artery intima-media thickness. There was evidence of GxE due to additive genetic variance heterogeneity and genetic correlation for FRS, PCRS, and CA-IMT measures for education (environment) but not for household income or SEI. Conclusion: The genetic effects underlying CVD are dynamically modulated at the lower end of the SES spectrum. There is a significant change in the genetic architecture underlying the major components of CVD in response to changes in education.

Estimation of genetic parameters for feed efficiency traits using random regression models in dairy cattle

Feed efficiency has become an increasingly important research topic in recent years. As feed costs rise and the environmental impacts of agriculture become more apparent, improving the efficiency with which dairy cows convert feed to milk is increasingly important. However, feed intake is expensive to measure accurately on large populations, making the inclusion of this trait in breeding programs difficult. Understanding how the genetic parameters of feed efficiency and traits related to feed efficiency vary throughout the lactation period is valuable to gain understanding into the genetic nature of feed efficiency. This study used 121,226 dry matter intake (DMI) records, 120,500 energy-corrected milk (ECM) records, and 98,975 metabolic body weight (MBW) records, collected on 7,440 first-lactation Holstein cows from 6 countries (Canada, Denmark, Germany, Spain, Switzerland, and the United States), from January 2003 to February 2022. Genetic parameters were estimated using a multiple-trait random regression model with a fourth-order Legendre polynomial for all traits. Weekly phenotypes for DMI were re-parameterized using linear regressions of DMI on ECM and MBW, creating a measure of feed efficiency that was genetically corrected for ECM and MBW, referred to as genomic residual feed intake (gRFI). Heritability (SE) estimates varied from 0.15 (0.03) to 0.29 (0.02) for DMI, 0.24 (0.01) to 0.29 (0.03) for ECM, 0.55 (0.03) to 0.83 (0.05) for MBW, and 0.12 (0.03) to 0.22 (0.06) for gRFI. In general, heritability estimates were lower in the first stage of lactation compared with the later stages of lactation. Additive genetic correlations between weeks of lactation varied, with stronger correlations between weeks of lactation that were close together. The results of this study contribute to a better understanding of the change in genetic parameters across the first lactation, providing insight into potential selection strategies to include feed efficiency in breeding programs.

Estimation of heritabilities and additive genetic correlations for reproduction traits in swine: insights for tropical commercial production systems using multiple trait animal models.

This study was to estimate heritabilities, additive genetic correlations, and phenotypic correlations between number of piglets born alive (NBA), litter birth weight (LTBW), number of piglets weaned (NPW) and litter weaning weight (LTWW) in different parities of Landrace (L), Yorkshire (Y), Landrace×Yorkshire (LY), and Yorkshire×Landrace (YL) sows in a commercial swine operation in Northern Thailand. Two models were utilized, a single trait repeatability model (RM) and a multiple trait animal model (MTM). The RM assumed reproductive records from different parities to be repeated values of the same trait, whereas the MTM assumed these records to be different traits. The two models accounted for the fixed effects of farrowing year-season, genetic group of the sow, heterosis, and age at first farrowing, and the random effects of sow, boar, and residual. Heritability estimates from RM were 0.02±0.01 for NBA, 0.10±0.01 for LTBW, 0.04±0.01 for NPW, and 0.11±0.01 for LTWW. Heritability estimates from MTM fluctuated across parities, ranging from 0.04±0.01 in parity 2 to 0.09±0.02 in parity 4 for NBA, 0.07±0.02 in parity 2 to 0.16±0.02 in parity 3 for LTBW, 0.04±0.02 in parity 4 to 0.08±0.01 in parity 1 for NPW, and 0.16±0.02 in parity 1 to 0.20±0.02 in parity 2 for LTWW. Additive genetic correlation estimates from MTM were also variable, ranging from 0.29±0.24 between NBA in parity 1 and NBA in parity 2 to 0.99±0.05 between LTWW in parity 3 and LTWW in parity 4. The findings of this study highlight the advantage of using MTM for the genetic improvement of reproductive traits in swine and contribute to the development of sustainable swine breeding programs in Thailand.

An assessment of the interaction between sucrose content and seed quality traits in soybeans

AbstractSoybean (Glycine max [L.] Merr.) cultivars with increased protein and sucrose content are desirable to enhance soybean meal quality; however, the complex interactions between soybean traits, mainly trade‐offs, make the plant breeders’ job challenging because traits cannot simultaneously be enhanced. Information about the interactions of seed carbohydrate composition with other seed traits remains inconclusive and contradictory; thus, this study explored the interactions among soluble carbohydrate content with other seed traits, flowering time, and maturity in a diverse panel of 1096 soybean accessions from the USDA Soybean Germplasm Collection. Trait interactions were explored through phenotypic, additive genetic, and residual correlations. We did not identify traits that can be useful in the indirect selection of sucrose. Additive genetic correlations of sucrose with flowering and maturity time suggest that obtaining cultivars with high sucrose content might be easier in earlier flowering and later maturing cultivars. Sucrose correlations with other seed traits were primarily driven by their additive genetic correlations, which make them more stable, whereas correlations of raffinose with oil and maturity were population and environment specific. The two main obstacles to enhancing soymeal quality are the trade‐offs of protein with sucrose and oil.

Genetic and Environmental Influences on Blood Pressure and Serum Lipids Across Age-Groups.

Aging plays a crucial role in the mechanisms of the impacts of genetic and environmental factors on blood pressure and serum lipids. However, to our knowledge, how the influence of genetic and environmental factors on the correlation between blood pressure and serum lipids changes with age remains to be determined. In this study, data from the Chinese National Twin Registry (CNTR) were used. Resting blood pressure, including systolic and diastolic blood pressure (SBP and DBP), and fasting serum lipids, including total cholesterol (TC), high-density lipoprotein cholesterol (HDL-C), low-density lipoprotein cholesterol (LDL-C) and triglycerides (TGs) were measured in 2378 participants (1189 twin pairs). Univariate and bivariate structural equation models examined the genetic and environmental influences on blood pressure and serum lipids among three age groups. All phenotypes showed moderate to high heritability (0.37-0.59) and moderate unique environmental variance (0.30-0.44). The heritability of all phenotypes showed a decreasing trend with age. Among all phenotypes, SBP and DBP showed a significant monotonic decreasing trend. For phenotype-phenotype pairs, the phenotypic correlation (Rph) of each pair ranged from -0.04 to 0.23, and the additive genetic correlation (Ra) ranged from 0.00 to 0.36. For TC&SBP, TC&DBP, TG&SBP and TGs&DBP, both the Rph and Ra declined with age, and the Ra difference between the young group and the older adult group is statistically significant (p < .05). The unique environmental correlation (Re) of each pair did not follow any pattern with age and remained relatively stable with age. In summary, we observed that the heritability of blood pressure was affected by age. Moreover, blood pressure and serum lipids shared common genetic backgrounds, and age had an impact on the phenotypic correlation and genetic correlations.

Tests for associations between sexual dimorphism and patterns of quantitative genetic variation in the water strider, Aquarius remigis

The evolution of sexual dimorphisms requires divergence between sexes in the evolutionary trajectories of the traits involved. Discerning how genetic architecture could facilitate such divergence has proven challenging because of the difficulty in estimating non-additive and sex-linked genetic variances using traditional quantitative genetic designs. Here we use a three-generation, double-first-cousin pedigree design to estimate additive, sex-linked and dominance (co)variances for 12 traits in the water strider, Aquarius remigis. Comparisons among these traits, which have size ratios ranging from 1 to 5 (larger/smaller), allow us to ask if sexual dimorphisms are associated with characteristic patterns of quantitative genetic variation. We frame our analysis around three main questions, derived from existing theory and empirical evidence: Are sexual dimorphisms associated with (1) lower additive inter-sex genetic correlations, (2) higher proportions of sex-linked variance, or (3) differences between sexes in autosomal additive and dominance genetic variances? For questions (1) and (2), we find weak and non-significant trends in the expected directions, which preclude definitive conclusions. However, in answer to question (3), we find strong evidence for a positive relationship between sexual dimorphism and differences between sexes in proportions of autosomal dominance variance. We also find strong interactions among the three genetic components indicating that their relative influence differs among traits and between sexes. These results highlight the need to include all three components of genetic (co)variance in both theoretical evolutionary models and empirical estimations of the genetic architecture of dimorphic traits.

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

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

Genetic Parameter Estimates from a Polymix Breeding Population of Pinus taeda L.

Abstract Estimation of genetic parameters from progeny testing is essential for many important decisions in forest tree breeding. In this study, we estimated heritabilities, trait-trait genetic correlations, and genotype by environment (GxE) interactions using a large multi-environmental data set of Pinus taeda L. in the southern United States. In the study, 284 parents were pollinated with a pollen mix and were field tested in four different test series. A total of twenty tests (two to seven tests per series) were established using randomized complete block design with single tree plots. Half-sib family-mean heritability estimates within individual test series varied considerably for all the traits. Pooled estimates of half-sib family-mean heritability estimates across series were 0.81 for fusiform rust incidence; they were relatively lower for tree height (0.68), diameter at breast height (DBH; 0.66), and stem straightness (0.58). Genetic correlations between pairs of traits were low except for height and diameter, which had a high genetic correlation (0.79). Additive genetic correlations between pairs of sites within test series (as a measure of GxE) varied. The pooled estimates ranged between 0.65 (DBH) to 0.80 (fusiform rust). The results suggest that polymix mating is efficient for parental selection in Pinus taeda L. and for deployment in the seed orchards. Study Implications: This study focused on genetic variation in the Piedmont breeding population of Pinus taeda L. in the United States. The Piedmont population was developed for more inland and colder regions of the species. Fusiform rust disease incidence of pine species of southeastern United States is a major threat to pine plantations. Understanding the level of genetic control on disease outcome is vital to develop sound breeding strategies and recommend deployment practices to mitigate the economic loss of landowners. The results provide valuable information to breeders to select disease-resistant genotypes to breed and put into seed orchards to mass-produce stock for planting.

Nonparallel genome changes within subpopulations over time contributed to genetic diversity within the US Holstein population

Maintaining genetic variation in a population is important for long-term genetic gain. The existence of subpopulations within a breed helps maintain genetic variation and diversity. The 20,990 genotyped animals, representing the breeding animals in the year 2014, were identified as the sires of animals born after 2010 with at least 25 progenies, and females measured for type traits within the last 2 yr of data. K-means clustering with 5 clusters (C1, C2, C3, C4, and C5) was applied to the genomic relationship matrix based on 58,990 SNP markers to stratify the selected candidates into subpopulations. The general higher inbreeding resulting from within-cluster mating than across-cluster mating suggests the successful stratification into genetically different groups. The largest cluster (C4) contained animals that were less related to each animal within and across clusters. The average fixation index was 0.03, indicating that the populations were differentiated, and allele differences across the subpopulations were not due to drift alone. Starting with the selected candidates within each cluster, a family unit was identified by tracing back through the pedigree, identifying the genotyped ancestors, and assigning them to a pseudogeneration. Each of the 5 families (F1, F2, F3, F4, and F5) was traced back for 10 generations, allowing for changes in frequency of individual SNPs over time to be observed, which we call allele frequencies change. Alternative procedures were used to identify SNPs changing in a parallel or nonparallel way across families. For example, markers that have changed the most in the whole population, markers that have changed differently across families, and genes previously identified as those that have changed in allele frequency. The genomic trajectory taken by each family involves selective sweeps, polygenic changes, hitchhiking, and epistasis. The replicate frequency spectrum was used to measure the similarity of change across families and showed that populations have changed differently. The proportion of markers that reversed direction in allele frequency change varied from 0.00 to 0.02 if the rate of change was greater than 0.02 per generation, or from 0.14 to 0.24 if the rate of change was greater than 0.005 per generation within each family. Cluster-specific SNP effects for stature were estimated using only females and applied to obtain indirect genomic predictions for males. Reranking occurs depending on SNP effects used. Additive genetic correlations between clusters show possible differences in populations. Further research is required to determine how this knowledge can be applied to maintain diversity and optimize selection decisions in the future.

Predicting milk protein fractions using infrared spectroscopy and a gradient boosting machine for breeding purposes in Holstein cattle

In recent years, increasing attention has been focused on the genetic evaluation of protein fractions in cow milk with the aim of improving milk quality and technological characteristics. In this context, advances in high-throughput phenotyping by Fourier transform infrared (FTIR) spectroscopy offer the opportunity for large-scale, efficient measurement of novel traits that can be exploited in breeding programs as indicator traits. We took milk samples from 2,558 Holstein cows belonging to 38 herds in northern Italy, operating under different production systems. Fourier transform infrared spectra were collected on the same day as milk sampling and stored for subsequent analysis. Two sets of data (i.e., phenotypes and FTIR spectra) collected in 2 different years (2013 and 2019-2020) were compiled. The following traits were assessed using HPLC: true protein, major casein fractions [αS1-casein (CN), αS2-CN, β-CN, κ-CN, and glycosylated-κ-CN], and major whey proteins (β-lactoglobulin and α-lactalbumin), all of which were measured both in grams per liter (g/L) and proportion of total nitrogen (% N). The FTIR predictions were calculated using the gradient boosting machine technique and tested by 3 different cross-validation (CRV) methods. We used the following CRV scenarios: (1) random 10-fold, which randomly split the whole into 10-folds of equal size (9-folds for training and 1-fold for validation); (2) herd/date-out CRV, which assigned 80% of herd/date as the training set with independence of 20% of herd/date assigned as the validation set; (3) forward/backward CRV, which split the data set in training and validation set according with the year of milk sampling (FTIR and gold standard data assessed in 2013 or 2019-2020) using the "old" and "new" databases for training and validation, and vice-versa with independence among them; (4) the CRV for genetic parameters (CRV-gen), where animals without pedigree as assigned as a fixed training population and animals with pedigree information was split in 5-folds, in which 1-fold was assigned to the fixed training population, and 4-folds were assigned to the validation set (independent from the training set). The results (i.e., measures and predictions) of CRV-gen were used to infer the genetic parameters for gold standard laboratory measurements (i.e., proteins assessed with HPLC) and FTIR-based predictions considering the CRV-gen scenario from a bi-trait animal model using single-step genomic BLUP. We found that the prediction accuracies of the gradient boosting machine equations differed according to the way in which the proteins were expressed, achieving higher accuracy when expressed in g/L than when expressed as % N in all CRV scenarios. Concerning the reproducibility of the equations over the different years, the results showed no relevant differences in predictive ability between using "old" data as the training set and "new" data as the validation set and vice-versa. Comparing the additive genetic variance estimates for milk protein fractions between the FTIR predicted and HPLC measures, we found reductions of -19.7% for milk protein fractions expressed in g/L, and -21.19% expressed as % N. Although we found reductions in the heritability estimates, they were small, with values ranging from -1.9 to -7.25% for g/L, and -1.6 to -7.9% for % N. The posterior distributions of the additive genetic correlations (ra) between the FTIR predictions and the laboratory measurements were generally high (>0.8), even when the milk protein fractions were expressed as % N. Our results show the potential of using FTIR predictions in breeding programs as indicator traits for the selection of animals to enhance milk protein fraction contents. We expect acceptable responses to selection due to the high genetic correlations between HPLC measurements and FTIR predictions.

Study of genetic parameters for pre-weaning growth traits in inner Mongolia white Arbas cashmere goats.

Inner Mongolia Arbas white cashmere goats is a dual-purpose breed for producing cashmere and meat. In recent years, its meat has become more and more popular among consumers because of rich nutrients and delicious flavor. Therefore, it is particularly important to study the genetic and non-genetic factors affecting the early growth traits and estimate variance components of pre-weaning growth traits of Inner Mongolia Albas white cashmere goats. A total of 37487 kidding records such as birth weight (BWT), weaning weight (WWT), average daily gain from birth to weaning (ADG) and Kleiber ratio (KR) from 343 sires and 7296 dams were used in this study. The most appropriate model was chosen on the basis of likelihood ratio test by fitting six models which excluding or including maternal genetic, maternal permanent environmental effects. The parameters were estimated under the most appropriate model using AIREML method by WOMBAT software. With the best model (Model 6), heritability estimates were 0.0435, 0.0911, 0.0932 and 0.2339 for BWT, WWT, ADG and KR traits, respectively. Maternal heritability estimates were 0.0143, 0.0246, 0.0220 and 0.0186 for BWT, WWT, ADG, and KR traits respectively. The correlation between different traits was estimated with the most suitable model by using bivariate analysis method. The direct additive genetic correlation among the traits ranged from -0.026 (BWT~KR) to 0.772 (ADG~KR). The maternal permanent environment correlation is between -0.289 (BWT-KR) ~0.900 (WWT-ADG). Results indicated that maternal effects and direct-maternal genetic covariance should be considered in any program aimed at improving pre-weaning growth traits to have an accurate genetic evaluation. In addition, positive and medium to high genetic correlations generally exist among WWT, ADG and KR due to the existence of genetic variation for early growth traits. The results showed that the genetic progress of these traits could be slowly through selection except for KR.

Genetic and Environmental Factors Underlying Parallel Changes in Body Mass Index and Alcohol Consumption: A 36-Year Longitudinal Study of Adult Twins

Introduction: While the genetic and environmental underpinnings of body weight and alcohol use are fairly well-known, determinants of simultaneous changes in these traits are still poorly known. We sought to quantify the environmental and genetic components underlying parallel changes in weight and alcohol consumption and to investigate potential covariation between them. Methods: The analysis comprised 4,461 adult participants (58% women) from the Finnish Twin Cohort with four measures of alcohol consumption and body mass index (BMI) over a 36-year follow-up. Trajectories of each trait were described by growth factors, defined as intercepts (i.e., baseline) and slopes (i.e., change over follow-up), using latent growth curve modeling. Growth values were used for male (190 monozygotic pairs, 293 dizygotic pairs) and female (316 monozygotic pairs, 487 dizygotic pairs) same-sex complete twin pairs in multivariate twin modeling. The variances and covariances of growth factors were then decomposed into genetic and environmental components. Results: The baseline heritabilities were similar in men (BMI: h² = 79% [95% confidence interval: 74, 83]; alcohol consumption: h² = 49% [32, 67]) and women (h² = 77% [73, 81]; h² = 45% [29, 61]). Heritabilities of BMI change were similar in men (h² = 52% [42, 61]) and women (h² = 57% [50, 63]), but the heritability of change in alcohol consumption was significantly higher (p = 0.03) in men (h² = 45% [34, 54]) than in women (h² = 31% [22, 38]). Significant additive genetic correlations between BMI at baseline and change in alcohol consumption were observed in both men (rA = −0.17 [−0.29, −0.04]) and women (rA = −0.18 [−0.31, −0.06]). Non-shared environmental factors affecting changes in alcohol consumption and BMI were correlated in men (rE = 0.18 [0.06, 0.30]). Among women, non-shared environmental factors affecting baseline alcohol consumption and the change in BMI were inversely correlated (rE = −0.11 [−0.20, −0.01]). Conclusions: Based on genetic correlations, genetic variation underlying BMI may affect changes in alcohol consumption. Independent of genetic effects, change in BMI correlates with change in alcohol consumption in men, suggesting direct effects between them.