Additive Genetic Variance Research Articles

Metabolic syndrome traits exhibit genotype-by-environment interaction in relation to socioeconomic status in the Mexican American family heart study.

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 metabolic disease. We analyzed data from the Mexican American Family Studies (MAFS) to evaluate the hypothesis that genotype-by-environment interaction (GxE) is an essential determinant of variation in risk factors for metabolic syndrome (MS). Methods: We employed a maximum likelihood estimation of the decomposition of variance components to detect GxE interaction. After excluding individuals with diabetes and individuals on medication for diabetes, hypertension, or dyslipidemia, we analyzed 12 MS risk factors: fasting glucose (FG), fasting insulin (FI), 2-h glucose (2G), 2-h insulin (2I), body mass index (BMI), waist circumference (WC), leptin (LP), high-density lipoprotein-cholesterol (HDL-C), triglycerides (TG), total serum cholesterol (TSC), systolic blood pressure (SBP), and diastolic blood pressure (DBP). Our SES variable used a combined score of Duncan's socioeconomic index and education years. Heterogeneity in the additive genetic variance across the SES continuum and a departure from unity in the genetic correlation coefficient were taken as evidence of GxE interaction. Hypothesis tests were conducted using standard likelihood ratio tests. Results: We found evidence of GxE for fasting glucose, 2-h glucose, 2-h insulin, BMI, and triglycerides. The genetic effects underlying the insulin/glucose metabolism component of MS are upregulated at the lower end of the SES spectrum. We also determined that the household variance for systolic blood pressure decreased with increasing SES. Conclusion: These results show a significant change in the GxE interaction underlying the major components of MS in response to changes in socioeconomic status. Further mRNA sequencing studies will identify genes and canonical gene pathways to support our molecular-level hypotheses.

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

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

Re-sequencing of the casein genes in Swedish Red cattle giving milk with diverse protein profiles and extreme rennet coagulation properties

Impaired rennet coagulation properties in milk could lead to prolonged processing times and production losses. Heritability for milk coagulation has previously been estimated to be 0.28 to 0.45, indicating that genetic selection can be used to manipulate this trait. The CN proteins are expressed by the genes CSN1S1, CSN2, CSN1S2, and CSN3 and are located on bovine chromosome 6. To better understand the effect of genetic variation in the CN genes on milk coagulation, blood and milk samples from 30 Swedish Red Dairy Cattle (RDC) with divergent coagulation properties were investigated. DNA from the 30 cows was sequenced for the CN genes to determine the theoretical AA sequence and to look for genetic variation in the untranslated regions. The aim is to confirm the protein genetic variants previously reported, while searching for additional genetic variation in the CN genes of 30 RDC. We observed genetic variation in 116 SNPs in the known CN genes where 10% of the SNPs are exon variants and the remaining 90% are intron variants. A total of 2.5% of the SNPs are found in the 5′- or 3′-untranslated region (UTR) regions of the exons; 2% are synonymous variants and 6% are missense variants that concurred with the known protein variants for CSN1S1, CSN2, and CSN3. Furthermore, 6% of the SNPs are splice polypyrimidine tract intron variants. The 2 genetic variants in the 5′- and 3′-UTR in CSN1S1 and CSN3 are found with protein variants CSN1S1C and CSN3B. Because both UTR variants are associated with gain and loss of micro RNA and transcription factors, this could explain differences in expression of the genetic protein variants. Preliminary chi-squared analysis and comparison with previous GWAS studies showed potential connections between the identified SNPs and coagulation properties of milk. By advancing the knowledge of the connection between the DNA sequence and the functional properties of the CN proteins, we hope to learn more about the cheese coagulation properties of milk from RDC.

Predicting bovine daily milk yield by leveraging genomic breeding values

The main goal of this work, conducted on a herd of 502 Holstein cows situated in Italy, is to propose a machine learning-based approach to forecast the individual bovine daily milk production by explicitly leveraging genotypic information. As part of our study, we also evaluated the importance in the prediction of genotypic and phenotypic variables usually available within herd. The methodology we propose is based on two consecutive models: a genomic prediction one to calculate the animal’s genomic breeding value from marker data, followed by a feed forward neural network combining such additive genetic effect and the environmental features (parity, days in milk, age at calving in months, month of calving) for milk yield forecasting. In particular, we both assess the inclusion of genomic breeding values calculated within herd or provided by the breeders’ association, discovering that the latter ones allow for better final predictions. The results of our model outperform the ones by a linear mixed model with the same inputs on average, day by day and at the individual level. Moreover, we propose a problem formulation that also leverages additional factors partially controllable by breeders: in this case, features such as the number of milkings and the concentrate consumption inside the automatic milking system prove to highly impact on the final prediction, and hence on milk production. To the best of our knowledge, the proposed problem formulation based on genomic breeding values is a novelty in the individual bovine milk yield machine learning forecasting literature. Given the low genotyping costs and the availability of a larger number of environmental features in farms equipped with a wide range of sensors, as automatic milking systems, our solution can support breeders’ herd management and animal monitoring, thanks to the possibility to forecast the full lactation curve in advance even for primiparous bovines and newborn calves. With this work, we successfully achieve our objectives of including genomic information in bovine milk yield machine learning-based forecasting, thus improving the performance on this task, and of evaluating the impact on prediction of common genotypic and phenotypic information available to breeders.

Adult body mass is heritable, positively genetically correlated and under selection of differing shapes between the sexes in a bird with little apparent sexual dimorphism.

In most animals, body mass varies with ecological conditions and is expected to reflect how much energy can be allocated to reproduction and survival. Because the sexes often differ in their resource acquisition and allocation strategies, variations in adult body mass and their consequences on fitness can differ between the sexes. Assessing the relative contributions of environmental and genetic effects (i.e. heritability)-and whether these effects and their fitness consequences are sex-specific-is essential to gain insights into the evolution of sexual dimorphism and sexual conflicts. We used 20+ years of data to study the sources of variation in adult body mass and associated fitness consequences in a bird with biparental care, the Alpine swift (Tachymarptis melba). Swifts appear monomorphic to human observers, though subtle dimorphisms are present. We first investigated the effects of weather conditions on adult body mass using a sliding window analysis approach. We report a positive effect of temperature and a negative effect of rainfall on adult body mass, as expected for an aerial insectivorous bird. We then quantified the additive genetic variance and heritability of body mass in both sexes and assessed the importance of genetic constraints on mass evolution by estimating the cross-sex genetic correlation. Heritability was different from zero in both sexes at ~0.30. The positive cross-sex genetic correlation and comparable additive genetic variance between the sexes suggest the possibility for evolutionary constraints when it comes to body mass. Finally, we assessed the sex-specific selection on adjusted body mass using multiple fitness components. We report directional positive and negative selection trending towards stabilizing and diversifying selection on females and males respectively in relation to the weighted proportion of surviving fledglings. Overall, these results suggest that while body mass may be able to respond to environmental conditions and evolve, genetic constraints would result in similar changes in both sexes or an overall absence of response to selection. It remains unclear whether the weak (1%) dimorphism in Alpine swift body mass we report is simply a result of the similar fitness peaks between the sexes or of genetic constraints.

Directional epistasis is common in morphological divergence.

Epistasis is often portrayed as unimportant in evolution. While random patterns of epistasis may have limited effects on the response to selection, systematic directional epistasis can have substantial effects on evolutionary dynamics. Directional epistasis occurs when allele substitutions that change a trait also modify the effects of allele substitution at other loci in a systematic direction. In this case, trait evolution may induce correlated changes in allelic effects and additive genetic variance (evolvability) that modify further evolution. Although theory thus suggests a potentially important role for directional epistasis in evolution, we still lack empirical evidence about its prevalence and magnitude. Using a new framework to estimate systematic patterns of epistasis from line-crosses experiments, we quantify its effects on 197 size-related traits from diverging natural populations in 24 animal and 17 plant species. We show that directional epistasis is common and tends to become stronger with increasing morphological divergence. In animals, most traits displayed negative directionality toward larger size, suggesting that epistatic constraints reducing evolvability toward larger size. Dominance was also common but did not systematically alter the effects of epistasis.

Dissecting the genetic architecture of quantitative traits using genome-wide identity-by-descent sharing.

Additive and dominance genetic variances underlying the expression of quantitative traits are important quantities for predicting short-term responses to selection, but they are notoriously challenging to estimate in most non-model wild populations. Specifically, large-sized or panmictic populations may be characterized by low variance in genetic relatedness among individuals which, in turn, can prevent accurate estimation of quantitative genetic parameters. We used estimates of genome-wide identity-by-descent (IBD) sharing from autosomal SNP loci to estimate quantitative genetic parameters for ecologically important traits in nine-spined sticklebacks (Pungitius pungitius) from a large, outbred population. Using empirical and simulated datasets, with varying sample sizes and pedigree complexity, we assessed the performance of different crossing schemes in estimating additive genetic variance and heritability for all traits. We found that low variance in relatedness characteristic of wild outbred populations with high migration rate can impair the estimation of quantitative genetic parameters and bias heritability estimates downwards. On the other hand, the use of a half-sib/full-sib design allowed precise estimation of genetic variance components and revealed significant additive variance and heritability for all measured traits, with negligible dominance contributions. Genome-partitioning and QTL mapping analyses revealed that most traits had a polygenic basis and were controlled by genes at multiple chromosomes. Furthermore, different QTL contributed to variation in the same traits in different populations suggesting heterogeneous underpinnings of parallel evolution at the phenotypic level. Our results provide important guidelines for future studies aimed at estimating adaptive potential in the wild, particularly for those conducted in outbred large-sized populations.

Genetic and environmental influence on alcohol intent and alcohol sips among U.S. children-Effects across sex, race, and ethnicity.

Alcohol intent (the susceptibility to initiating alcohol use) and alcohol sips (the initiation of alcohol) in youth are a multifactorial puzzle with many components. This research aims to examine the connection between genetic and environmental factors across sex, race and ethnicity. Data was obtained from the twin hub of the Adolescent Brain Cognitive Development (ABCD) study at baseline (2016-2018). Variance component models were conducted to dissect the additive genetic (A), common (C) and unique environmental (E) effects on alcohol traits. The proportion of the total alcohol phenotypic variation attributable to additive genetic factors is reported as heritability (h2). The sample (n = 1,772) included an approximately equal male-female distribution. The 886 same-sex twin pairs were 60.4% dizygotic (DZ), 39.6% monozygotic (MZ), 65.4% non-Hispanic Whites, 13.9% non-Hispanic Blacks, 10.8% of Hispanics with a mean age of 121.2 months. Overall, genetic predisposition was moderate for alcohol intent (h2 = 28%, p = .006) and low for alcohol initiation (h2 = 4%, p = 0.83). Hispanics (h2 = 53%, p < .0001) and Blacks (h2 = 48%, p < .0001) demonstrated higher alcohol intent due to additive genetic factors than Whites (h2 = 34%, p < .0001). Common environmental factors explained more variation in alcohol sips in females (c2 = 63%, p = .001) than in males (c2 = 55%, p = .003). Unique environmental factors largely attributed to alcohol intent, while common environmental factors explained the substantial variation in alcohol initiation. Sex and racial/ethnic disparities in genetic and environmental risk factors for susceptibility to alcohol initiation can lead to significant health disparities. Certain populations may be at greater risk for alcohol use due to their genetic and ecological factors at an early age.

Investigating pedigree- and SNP-associated components of heritability in a wild population of Soay sheep

Estimates of narrow sense heritability derived from genomic data that contain related individuals may be biased due to the within-family effects such as dominance, epistasis and common environmental factors. However, for many wild populations, removal of related individuals from the data would result in small sample sizes. In 2013, Zaitlen et al. proposed a method to estimate heritability in populations that include close relatives by simultaneously fitting an identity-by-state (IBS) genomic relatedness matrix (GRM) and an identity-by-descent (IBD) GRM. The IBD GRM is identical to the IBS GRM, except relatedness estimates below a specified threshold are set to 0. We applied this method to a sample of 8557 wild Soay sheep from St. Kilda, with genotypic information for 419,281 single nucleotide polymorphisms. We aimed to see how this method would partition heritability into population-level (IBS) and family-associated (IBD) variance for a range of genetic architectures, and so we focused on a mixture of polygenic and monogenic traits. We also implemented a variant of the model in which the IBD GRM was replaced by a GRM constructed from SNPs with low minor allele frequency to examine whether any additive genetic variance is captured by rare alleles. Whilst the inclusion of the IBD GRM did not significantly improve the fit of the model for the monogenic traits, it improved the fit for some of the polygenic traits, suggesting that dominance, epistasis and/or common environment not already captured by the non-genetic random effects fitted in our models may influence these traits.

Sex-specific genetic parameter estimates of body weight in Mazandaran indigenous chickens.

Body weight is an economically important trait in poultry that shows sexual dimorphism (SD). In the present study, variation in SD in Mazandaran native chickens was investigated in terms of the (Co) variance components and genetic parameters of body weight between males and females. Studied traits were body weights at hatch (BW1), 8 weeks (BW8) and 12 weeks of age (BW12). Also, for weight at sexual maturity (WSM) covariance components were only estimated in females. Cross-sex direct and maternal correlations were also estimated for studied traits except for WSM. For this purpose, a deep 21-generation pedigree and body weight data (57,576 BW1, 72,925 BW8, 62,727 BW12 and, 42,496 WSM) were used. Evaluation of SD of body weight was performed using six bivariate animal models with and without considering the genetic and permanent maternal environmental effects under the restricted maximum likelihood method in WOMBAT software. Model with direct additive genetic effects and maternal genetic effects without covariance between them was identified as the best model for BW1 and BW8. The Model including direct additive genetic effects and permanent maternal environmental effects was the best model for BW12 and WSM. Direct heritability (h2 ) estimates for BW1, BW8 and, BW12 were, respectively, 0.05 ± 0.013, 0.17 ± 0.02 and, 0.25 ± 0.03 in males and, 0.05 ± 0.012, 0.15 ± 0.01 and 0.21 ± 0.01 in females. Also, the direct heritability of WSM based on univariate analysis in females was estimated to be 0.40 ± 0.01. Maternal heritability ( ) varied from 0.39 ± 0.01 (BW1) to 0.04 ± 0.009 (BW8) in males, and 0.36 ± 0.10 (BW1) to 0.04 ± 0.006 (BW8) in females. The correlation between direct genetic effects between males and females for BW1 was not significantly different from one. The direct genetic correlation between the two sexes for BW8 and BW12 was significantly different from 1 concluding that these traits are dimorphic in terms of direct genetic effects and therefore independent selection in both sexes is possible.

The estimation of additive genetic variance of body size in a wild passerine is sensitive to the method used to estimate relatedness among the individuals.

Assessing additive genetic variance is a crucial step in predicting the evolutionary response of a target trait. However, the estimated genetic variance may be sensitive to the methodology used, e.g., the way relatedness is assessed among the individuals, especially in wild populations where social pedigrees can be inaccurate. To investigate this possibility, we investigated the additive genetic variance in tarsus length, a major proxy of skeletal body size in birds. The model species was the collared flycatcher (Ficedula albicollis), a socially monogamous but genetically polygamous migratory passerine. We used two relatedness matrices to estimate the genetic variance: (1) based solely on social links and (2) a genetic similarity matrix based on a large array of single-nucleotide polymorphisms (SNPs). Depending on the relatedness matrix considered, we found moderate to high additive genetic variance and heritability estimates for tarsus length. In particular, the heritability estimates were higher when obtained with the genetic similarity matrix instead of the social pedigree. Our results confirm the potential for this crucial trait to respond to selection and highlight methodological concerns when calculating additive genetic variance and heritability in phenotypic traits. We conclude that using a social pedigree instead of a genetic similarity matrix to estimate relatedness among individuals in a genetically polygamous wild population may significantly deflate the estimates of additive genetic variation.

Mitigating type 1 error inflation and power loss in GxE PRS: Genotype-environment interaction in polygenic risk score models.

The use of polygenic risk score (PRS) models has transformed the field of genetics by enabling the prediction of complex traits and diseases based on an individual's genetic profile. However, the impact of genotype-environment interaction (GxE) on the performance and applicability of PRS models remains a crucial aspect to be explored. Currently, existing genotype-environment interaction polygenic risk score (GxE PRS) models are often inappropriately used, which can result in inflated type 1 error rates and compromised results. In this study, we propose novel GxE PRS models that jointly incorporate additive and interaction genetic effects although also including an additional quadratic term for nongenetic covariates, enhancing their robustness against model misspecification. Through extensive simulations, we demonstrate that our proposed models outperform existing models in terms of controlling type 1 error rates and enhancing statistical power. Furthermore, we apply the proposed models to real data, and report significant GxE effects. Specifically, we highlight the impact of our models on both quantitative and binary traits. For quantitative traits, we uncover the GxE modulation of genetic effects on body mass index by alcohol intake frequency. In the case of binary traits, we identify the GxE modulation of genetic effects on hypertension by waist-to-hip ratio. These findings underscore the importance of employing a robust model that effectively controls type 1 error rates, thus preventing the occurrence of spurious GxE signals. To facilitate the implementation of our approach, we have developed an innovative R software package called GxEprs, specifically designed to detect and estimate GxE effects. Overall, our study highlights the importance of accurate GxE modeling and its implications for genetic risk prediction, although providing a practical tool to support further research in this area.

The effect of repeated measurements and within-individual variance on the estimation of heritability: a simulation study

The estimation of heritability is a common practice in the field of ecology and evolution. Heritability of the traits is often estimated using one single measurement per individual, although many traits (especially behavioural and physiological traits) are characterized by large within-individual variance, and ideally a large number of within individual measurements can be obtained. Importantly, the effect of the within-individual variance and the rate at which this variance is sampled on the estimation of heritability has not been thoroughly tested. We fill this gap of knowledge with a simulation study, and assess the effect of within- and between-individual sample size, and the true value of the variance components on the estimation of heritability. In line with previous studies we found that the accuracy and precision of heritability estimation increased with sample size and accuracy with higher values of additive genetic variance. When the sample size was above 500 accuracy and power of heritability estimates increased in the models including repeated measurements, especially when within-individual variance was high. We thus suggest to use a sample of more than 100 individuals and to include more than two repeated measurements per individual in the models to improve estimation when investigating heritability of labile traits.Significance statementHeritability reflects the part of the trait’s phenotypic variation underlined by genetic variation. Despite the difficulties of heritability calculation (high number of individuals is needed with known relatedness), it is a widely used measure in evolutionary studies. However, not every factor potentially affecting the quality of heritability estimation is well understood. We thus investigated with a comprehensive simulation study how the number of repeated measurements per individuals and the amount of within-individual variation influence the goodness of heritability estimation. We found that although the previously described effect of the number of studied individuals was the most important, including repeated measurements also improved the reliability of the heritability estimates, especially when within-individual variation was high. Our results thus highlight the importance of including repeated measurements when investigating the heritability of highly plastic traits, such as behavioural or physiological traits.

Field Insights into the Genetic Symphony: Decoding Gene Dynamics for Enhanced Sorghum Yield (Sorghum bicolor [L.] Moench)

In this investigation using line × tester analysis, fifteen parents (3 Cytoplasmic Male Sterility females and 12 males) and 36 hybrids, including the standard check CSH-25, underwent evaluation in a randomized block design with three replications. Notably, significant general combining ability variances for 1000-seed weight and highly significant specific combining ability variances for all traits were observed. The research highlights the prevalence of additive genetic variance in 1000-seed weight and the combined influence of additive and nonadditive gene effects in panicle length ratio. Recognized good general combiners, such as SR 2991 , SR 2950 , SR 2872, SR 2984 and 691A, provide valuable insights for breeding strategies. Subsequent evaluations of promising hybrids, particularly 691 A × SR 2991, 691 A × SR 2872 and 691 A × SR 2950 hold potential for optimizing grain yield, panicle length, and dwarfness. These findings provide nuanced insights for future Sorghum breeding programs, enhancing understanding of gene dynamics for improved productivity and emphasizing the pivotal role of hybrid sorghum in contemporary agricultural contexts.

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

BackgroundingStayability, which may be defined as the probability of a cow remaining in the herd until a reference age or at a specific number of calvings, is usually measured late in the animal’s life. Thus, if used as selection criteria, it will increase the generation interval and consequently might decrease the annual genetic gain. Measuring stayability at an earlier age could be a reasonable strategy to avoid this problem. In this sense, a better understanding of the genetic architecture of this trait at different ages and/or at different calvings is important. This study was conducted to identify possible regions with major effects on stayability measured considering different numbers of calvings in Nellore cattle as well as pathways that can be involved in its expression throughout the female’s productive life.ResultsThe top 10 most important SNP windows explained, on average, 17.60% of the genetic additive variance for stayability, varying between 13.70% (at the eighth calving) and 21% (at the fifth calving). These SNP windows were located on 17 chromosomes (1, 2, 4, 6, 7, 8, 9, 10, 11, 12, 13, 14, 18, 19, 20, 27, and 28), and they harbored a total of 176 annotated genes. The functional analyses of these genes, in general, indicate that the expression of stayability from the second to the sixth calving is mainly affected by genetic factors related to reproductive performance, and nervous and immune systems. At the seventh and eighth calvings, genes and pathways related to animal health, such as density bone and cancer, might be more relevant.ConclusionOur results indicate that part of the target genomic regions in selecting for stayability at earlier ages (from the 2th to the 6th calving) would be different than selecting for this trait at later ages (7th and 8th calvings). While the expression of stayability at earlier ages appeared to be more influenced by genetic factors linked to reproductive performance together with an overall health/immunity, at later ages genetic factors related to an overall animal health gain relevance. These results support that selecting for stayability at earlier ages (perhaps at the second calving) could be applied, having practical implications in breeding programs since it could drastically reduce the generation interval, accelerating the genetic progress.

Genetic parameters for novel mastitis traits defined by combining test-day somatic cell score and differential somatic cell count in the first lactation of Japanese Holsteins

In this study, we aimed to improve current udder health genetic evaluations by addressing the limitations of monthly sampled somatic cell score (SCS) for distinguishing cows with robust innate immunity from those susceptible to chronic infections. The objectives were to (1) establish novel somatic cell traits by integrating SCS and the differential somatic cell count (DSCC), which represents the combined proportion of polymorphonuclear leukocytes and lymphocytes in somatic cells and (2) estimate genetic parameters for the new traits, including their daily heritability and genetic correlations with milk production traits and SCS, using a random regression test-day model (RRTDM). We derived 3 traits, namely ML_SCS_DSCC, SCS_4_DSCC_65_binary, and ML_SCS_DSCC_binary, by using milk loss estimates at corresponding SCS and DSCC levels, thresholds established in previous studies, and a threshold established from milk loss estimates, respectively. Data consisted of test-day records collected during January 2021 through March 2022 from 265 herds in Hokkaido, Japan. From these records, we extracted records between 7 to 305 d in milk (DIM) in the first lactation to fit the RRTDM. The model included the random effect of herd-test-day, the fixed effect of year-month, fixed lactation curves nested with calving age groups, and random regressions with Legendre polynomials of order 3 for additive genetic and permanent environmental effects. The analysis was performed using Gibbs sampling with Gibbsf90+ software. The averages (ranges) of daily heritability estimates over lactation were 0.086 (0.075 to 0.095) for SCS, 0.104 (0.073 to 0.127) for ML_SCS_DSCC, 0.137 (0.014 to 0.297) for SCS_4_DSCC_65_binary, and 0.138 (0.115 to 0.185) for ML_SCS_DSCC_binary; the heritability curve for SCS_4_DSCC_65_binary was erratic. Genetic correlations within the trait decreased as the DIM interval widened, especially for those integrating DSCC, indicating that these traits should be analyzed using RRTDM rather than repeatability models. The averages (ranges) of genetic correlations with milk yield over lactation were 0.01 (-0.22 to 0.28) for SCS, -0.05 (-0.40 to 0.13) for ML_SCS_DSCC, -0.08 (-0.17 to 0.09) for SCS_4_DSCC_65_binary, and -0.08 (-0.22 to 0.27) for ML_SCS_DSCC_binary. Compared with SCS, the newly defined traits exhibited slightly stronger negative genetic correlations with milk yield. Especially in late lactation stages, the genetic correlation between ML_SCS_DSCC and milk yield was significantly below zero, with a posterior median of -0.40. Furthermore, the new traits showed positive correlations with SCS, having estimates varying from 0.68 to 0.85 for ML_SCS_DSCC, 0.14 to 0.47 for SCS_4_DSCC_65_binary, and 0.61 to 0.66 for ML_SCS_DSCC_binary, depending on DIM. Considering that ML_SCS_DSCC and ML_SCS_DSCC_binary have relatively high heritability (compared with SCS) and favorable genetic correlations with milk production traits and SCS, their incorporation into breeding programs appears promising. Nevertheless, their genetic relationships with (sub)clinical mastitis require further investigation.

The Combining Ability and Heterosis Analysis of Sweet-Waxy Corn Hybrids for Yield-Related Traits and Carotenoids.

Improving sweet-waxy corn hybrids enriched in carotenoids via a hybrid breeding approach may provide an alternative cash crop for growers and provide health benefits for consumers. This study estimates the combining ability and heterosis of sweet-waxy corn hybrids for yield-related traits and carotenoids. Eight super sweet corn and three waxy corn lines were crossed to generate 24 F1 hybrids according to the North Carolina Design II scheme, and these hybrids were evaluated across two seasons of 2021/22. The results showed that both additive and non-additive genetic effects were involved in expressing the traits, but the additive genetic effect was more predominant. Most observed traits exhibited moderate to high narrow-sense heritability. Three parental lines, namely the ILS2 and ILS7 females and the ILW1 male, showed the highest positive GCA effects on yield-related traits, making them desirable for developing high-yielding hybrids. Meanwhile, five parental lines, namely the ILS3, ILS5, and ILS7 females and the ILW1 and ILW2 males, were favorable general combiners for high carotenoids. A tested hybrid, ILS2 × ILW1, was a candidate biofortified sweet-waxy corn hybrid possessing high yields and carotenoids. Heterosis and per se performance were more positively correlated with GCAsum than SCA, indicating that GCAsum can predict heterosis for improving biofortified sweet-waxy corn hybrid enriched in carotenoids. The breeding strategies of biofortified sweet-waxy corn hybrids with high yield and carotenoid content are discussed.

Environmental effects on genetic variance are likely to constrain adaptation in novel environments.

Adaptive plasticity allows populations to cope with environmental variation but is expected to fail as conditions become unfamiliar. In novel conditions, populations may instead rely on rapid adaptation to increase fitness and avoid extinction. Adaptation should be fastest when both plasticity and selection occur in directions of the multivariate phenotype that contain abundant genetic variation. However, tests of this prediction from field experiments are rare. Here, we quantify how additive genetic variance in a multivariate phenotype changes across an elevational gradient, and test whether plasticity and selection align with genetic variation. We do so using two closely related, but ecologically distinct, sister species of Sicilian daisy (Senecio, Asteraceae) adapted to high and low elevations on Mt. Etna. Using a quantitative genetic breeding design, we generated and then reciprocally planted c. 19,000 seeds of both species, across an elevational gradient spanning each species' native elevation, and then quantified mortality and five leaf traits of emergent seedlings. We found that genetic variance in leaf traits changed more across elevations than between species. The high-elevation species at novel lower elevations showed changes in the distribution of genetic variance among the leaf traits, which reduced the amount of genetic variance in the directions of selection and the native phenotype. By contrast, the low-elevation species mainly showed changes in the amount of genetic variance at the novel high elevation, and genetic variance was concentrated in the direction of the native phenotype. For both species, leaf trait plasticity across elevations was in a direction of the multivariate phenotype that contained a moderate amount of genetic variance. Together, these data suggest that where plasticity is adaptive, selection on genetic variance for an initially plastic response could promote adaptation. However, large environmental effects on genetic variance are likely to reduce adaptive potential in novel environments.

Wright was right: leveraging old data and new methods to illustrate the critical role of epistasis in genetics and evolution.

Much of evolutionary theory is predicated on assumptions about the relative importance of simple additive versus complex epistatic genetic architectures. Previous work suggests traits strongly associated with fitness will lack additive genetic variation, whereas traits less strongly associated with fitness are expected to exhibit more additive genetic variation. We use a quantitative genetics method, line cross analysis, to infer genetic architectures that contribute to trait divergence. By parsing over 1,600 datasets by trait type, clade, and cross divergence, we estimated the relative importance of epistasis across the tree of life. In our comparison between life-history traits and morphological traits, we found greater epistatic contributions to life-history traits. Our comparison between plants and animals showed that animals have more epistatic contribution to trait divergence than plants. In our comparison of within-species versus between-species crosses, we found that only animals exhibit a greater epistatic contribution to trait divergence as divergence increases. While many scientists have argued that epistasis is ultimately of little importance, our results show that epistasis underlies much of trait divergence and must be accounted for in theory and practical applications like domestication, conservation breeding design, and understanding complex diseases.

Scald resistance in hybrid rye (Secale cereale): genomic prediction and GWAS.

Rye (Secale cereale L.) is an important cereal crop used for food, beverages, and feed, especially in North-Eastern Europe. While rye is generally more tolerant to biotic and abiotic stresses than other cereals, it still can be infected by several diseases, including scald caused by Rhynchosporium secalis. The aims of this study were to investigate the genetic architecture of scald resistance, to identify genetic markers associated with scald resistance, which could be used in breeding of hybrid rye and to develop a model for genomic prediction for scald resistance. Four datasets with records of scald resistance on a population of 251 hybrid winter rye lines grown in 2 years and at 3 locations were used for this study. Four genomic models were used to obtain variance components and heritabilities of scald resistance. All genomic models included additive genetic effects of the parental components of the hybrids and three of the models included additive-by-additive epistasis and/or dominance effects. All models showed moderate to high broad sense heritabilities in the range of 0.31 (SE 0.05) to 0.76 (0.02). The model without non-additive genetic effects and the model with dominance effects had moderate narrow sense heritabilities ranging from 0.24 (0.06) to 0.55 (0.08). None of the models detected significant non-additive genomic variances, likely due to a limited data size. A genome wide association study was conducted to identify markers associated with scald resistance in hybrid winter rye. In three datasets, the study identified a total of twelve markers as being significantly associated with scald resistance. Only one marker was associated with a major quantitative trait locus (QTL) influencing scald resistance. This marker explained 11-12% of the phenotypic variance in two locations. Evidence of genotype-by-environment interactions was found for scald resistance between one location and the other two locations, which suggested that scald resistance was influenced by different QTLs in different environments. Based on the results of the genomic prediction models and GWAS, scald resistance seems to be a quantitative trait controlled by many minor QTL and one major QTL, and to be influenced by genotype-by-environment interactions.