Equivalence of variance components between standard and recursive genetic models using LDL' transformations.

BACKGROUND: Recursive models are a category of structural equation models that propose a causal relationship between traits. These models are more parameterized than multiple trait models, and they require imposing restrictions on the parameter space to ensure statistical identification. Nevertheless, in certain situations, the likelihood of recursive models and multiple trait models are equivalent. Consequently, the estimates of variance components derived from the multiple trait mixed model can be converted into estimates under several recursive models through LDL' or block-LDL' transformations. RESULTS: The procedure was employed on a dataset comprising five traits (birth weight-BW, weight at 90 days-W90, weight at 210 days-W210, cold carcass weight-CCW and conformation-CON) from the Pirenaica beef cattle breed. These phenotypic records were unequally distributed among 149,029 individuals and had a high percentage of missing data. The pedigree used consisted of 343,753 individuals. A Bayesian approach involving a multiple-trait mixed model was applied using a Gibbs sampler. The variance components obtained at each iteration of the Gibbs sampler were subsequently used to estimate the variance components within three distinct recursive models. CONCLUSIONS: The LDL' or block-LDL' transformations applied to the variance component estimates achieved from a multiple trait mixed model enabled inference across multiple sets of recursive models, with the sole prerequisite of being likelihood equivalent. Furthermore, the aforementioned transformations simplify the handling of missing data when conducting inference within the realm of recursive models.

Integrating perspectives in developing mycorrhizal trait databases: a call for inclusive and collaborative continental efforts.

Global assessments of mycorrhizal symbiosis present large sampling gaps in rich biodiversity regions. Filling these gaps is necessary to build large-scale, unbiased mycorrhizal databases to obtain reliable analyses and prevent misleading generalizations. Underrepresented regions in mycorrhizal research are mainly in Africa, Asia, and South America. Despite the high biodiversity and endemism in these regions, many groups of organisms remain understudied, especially mycorrhizal fungi. In this Viewpoint, we emphasize the importance of inclusive and collaborative continental efforts in integrating perspectives for comprehensive trait database development and propose a conceptual framework that can help build large mycorrhizal databases in underrepresented regions. Based on the four Vs of big data (volume, variety, veracity, and velocity), we identify the main challenges of constructing a large mycorrhizal dataset and propose solutions for each challenge. We share our collaborative methodology, which involves employing open calls and working groups to engage all mycorrhizal researchers in the region to build a South American Mycorrhizal Database. By fostering interdisciplinary collaborations and embracing a continental-scale approach, we can create robust mycorrhizal trait databases that provide valuable insights into the evolution, ecology, and functioning of mycorrhizal associations, reducing the geographical biases that are so common in large-scale ecological studies.

Mining genomic regions associated with agronomic and biochemical traits in quinoa through GWAS.

Quinoa (Chenopodium quinoa Willd.), an Andean crop, is a facultative halophyte food crop recognized globally for its high nutritional value and plasticity to adapt to harsh conditions. We conducted a genome-wide association study on a diverse set of quinoa germplasm accessions. These accessions were evaluated for the following agronomic and biochemical traits: days to 50% flowering (DTF), plant height (PH), panicle length (PL), stem diameter (SD), seed yield (SY), grain diameter (GD), and thousand-grain weight (TGW). These accessions underwent genotyping-by-sequencing using the DNBSeq-G400R platform. Among all evaluated traits, TGW represented maximum broad-sense heritability. Our study revealed average SNP density of ≈ 3.11 SNPs/10 kb for the whole genome, with the lowest and highest on chromosomes Cq1B and Cq9A, respectively. Principal component analysis clustered the quinoa population in three main clusters, one clearly representing lowland Chilean accessions, whereas the other two groups corresponded to germplasm from the highlands of Peru and Bolivia. In our germplasm set, we estimated linkage disequilibrium decay to be ≈ 118.5 kb. Marker-trait analyses revealed major and consistent effect associations for DTF on chromosomes 3A, 4B, 5B, 6A, 7A, 7B and 8B, with phenotypic variance explained (PVE) as high as 19.15%. Nine associations across eight chromosomes were also found for saponin content with 20% PVE by qSPN5A.1. More QTLs were identified for PL and TGW on multiple chromosomal locations. We identified putative candidate genes in the genomic regions associated with DTF and saponin content. The consistent and major-effect genomic associations can be used in fast-tracking quinoa breeding for wider adaptation across marginal environments.

Marker Density and Models to Improve the Accuracy of Genomic Selection for Growth and Slaughter Traits in Meat Rabbits.

The selection and breeding of good meat rabbit breeds are fundamental to their industrial development, and genomic selection (GS) can employ genomic information to make up for the shortcomings of traditional phenotype-based breeding methods. For the practical implementation of GS in meat rabbit breeding, it is necessary to assess different marker densities and GS models. Here, we obtained low-coverage whole-genome sequencing (lcWGS) data from 1515 meat rabbits (including parent herd and half-sibling offspring). The specific objectives were (1) to derive a baseline for heritability estimates and genomic predictions based on randomly selected marker densities and (2) to assess the accuracy of genomic predictions for single- and multiple-trait linear mixed models. We found that a marker density of 50 K can be used as a baseline for heritability estimation and genomic prediction. For GS, the multi-trait genomic best linear unbiased prediction (GBLUP) model results in more accurate predictions for virtually all traits compared to the single-trait model, with improvements greater than 15% for all of them, which may be attributed to the use of information on genetically related traits. In addition, we discovered a positive correlation between the performance of the multi-trait GBLUP and the genetic correlation between the traits. We anticipate that this approach will provide solutions for GS, as well as optimize breeding programs, in meat rabbits.

Country-wide, multi-location trials of Green Super Rice lines for yield performance and stability analysis using genetic and stability parameters.

Rice (Oryza sativa L.) is an important member of the family Poaceae and more than half of world population depend for their dietary nutrition on rice. Rice cultivars with higher yield, resilience to stress and wider adaptability are essential to ensure production stability and food security. The fundamental objective of this study was to identify higher-yielding rice genotypes with stable performance and wider adaptability in a rice growing areas of Pakistan. A triplicate RCBD design experiment with 20 Green Super Rice (GSR) advanced lines was conducted at 12 rice growing ecologies in four Provinces of Pakistan. Grain yield stability performance was assessed by using different univariate and multivariate statistics. Analysis of variance revealed significant differences among genotypes, locations, and G x E interaction for mean squares (p < 0.05) of major yield contributing traits. All the studied traits except for number of tillers per plant revealed higher genotypic variance than environmental variance. Broad sense heritability was estimated in the range of 44.36% to 98.60%. Based on ASV, ASI, bi, Wi2, σ2i and WAAS statistics, the genotypes G1, G4, G5, G8, G11 and G12 revealed lowest values for parametric statistics and considered more stable genotypes based on  paddy yield. The additive main effects and multiplicative interaction (AMMI) model revealed significant variation (p < 0.05) for genotypes, non-signification for environment and highly significant for G × E interaction. The variation proportion of PC1 and PC2 from interaction revealed 67.2% variability for paddy yield. Based on 'mean verses stability analysis of GGE biplot', 'Which-won-where' GGE Biplot, 'discriminativeness vs. representativeness' pattern of stability, 'IPCA and WAASB/GY' ratio-based stability Heat-map, and ranking of genotypes, the genotypes G1, G2, G3, G5, G8, G10, G11 and G13 were observed ideal genotypes with yield potential more than 8 tons ha-1. Discriminativeness vs. representativeness' pattern of stability identifies two environments, E5 (D.I Khan, KPK) and E6 (Usta Muhammad, Baluchistan) were best suited for evaluating genotypic yield performance. Based on these findings we have concluded that the genotypes G1, G2, G3, G5, G8, G10, G11 and G13 could be included in the commercial varietal development process and future breeding program.

Better together against genetic heterogeneity: A sex-combined joint main and interaction analysis of 290 quantitative traits in the UK Biobank.

Genetic effects can be sex-specific, particularly for traits such as testosterone, a sex hormone. While sex-stratified analysis provides easily interpretable sex-specific effect size estimates, the presence of sex-differences in SNP effect implies a SNP×sex interaction. This suggests the usage of the often overlooked joint test, testing for an SNP's main and SNP×sex interaction effects simultaneously. Notably, even without individual-level data, the joint test statistic can be derived from sex-stratified summary statistics through an omnibus meta-analysis. Utilizing the available sex-stratified summary statistics of the UK Biobank, we performed such omnibus meta-analyses for 290 quantitative traits. Results revealed that this approach is robust to genetic effect heterogeneity and can outperform the traditional sex-stratified or sex-combined main effect-only tests. Therefore, we advocate using the omnibus meta-analysis that captures both the main and interaction effects. Subsequent sex-stratified analysis should be conducted for sex-specific effect size estimation and interpretation.

Identifying latent genetic interactions in genome-wide association studies using multiple traits.

The "missing" heritability of complex traits may be partly explained by genetic variants interacting with other genes or environments that are difficult to specify, observe, and detect. We propose a new kernel-based method called Latent Interaction Testing (LIT) to screen for genetic interactions that leverages pleiotropy from multiple related traits without requiring the interacting variable to be specified or observed. Using simulated data, we demonstrate that LIT increases power to detect latent genetic interactions compared to univariate methods. We then apply LIT to obesity-related traits in the UK Biobank and detect variants with interactive effects near known obesity-related genes (URL: https://CRAN.R-project.org/package=lit ).

Exploring the role of biotic factors in regulating the spatial variability in land surface phenology across four temperate forest sites.

Land surface phenology (LSP), the characterization of plant phenology with satellite data, is essential for understanding the effects of climate change on ecosystem functions. Considerable LSP variation is observed within local landscapes, and the role of biotic factors in regulating such variation remains underexplored. In this study, we selected four National Ecological Observatory Network terrestrial sites with minor topographic relief to investigate how biotic factors regulate intra-site LSP variability. We utilized plant functional type (PFT) maps, functional traits, and LSP data to assess the explanatory power of biotic factors for the start and end of season (SOS and EOS) variability. Our results indicate that PFTs alone explain only 0.8-23.4% of intra-site SOS and EOS variation, whereas including functional traits significantly improves explanatory power, with cross-validation correlations ranging from 0.50 to 0.85. While functional traits exhibited diverse effects on SOS and EOS across different sites, traits related to competitive ability and productivity were important for explaining both SOS and EOS variation at these sites. These findings reveal that plants exhibit diverse phenological responses to comparable environmental conditions, and functional traits significantly contribute to intra-site LSP variability, highlighting the importance of intrinsic biotic properties in regulating plant phenology.

Estimation of genetic parameters for the implementation of selective breeding in commercial insect production.

BACKGROUND: There is a burgeoning interest in using insects as a sustainable source of food and feed, particularly by capitalising on various waste materials and by-products that are typically considered of low value. Enhancing the commercial production of insects can be achieved through two main approaches: optimising environmental conditions and implementing selective breeding strategies. In order to successfully target desirable traits through selective breeding, having a thorough understanding of the genetic parameters pertaining to those traits is essential. In this study, a full-sib half-sib mating design was used to estimate variance components and heritabilities for larval size and survival at day seven of development, development time and survival from egg to adult, and to estimate correlations between these traits, within an outbred population of house flies (Musca domestica), using high-throughput phenotyping for data collection. RESULTS: The results revealed low to intermediate heritabilities and positive genetic correlations between all traits except development time and survival to day seven of development and from egg to adulthood. Surprisingly, larval size at day seven exhibited a comparatively low heritability (0.10) in contrast to development time (0.25), a trait that is believed to have a stronger association with overall fitness. A decline in family numbers resulting from low mating success and high overall mortality reduced the amount of available data which resulted in large standard errors for the estimated parameters. Environmental factors made a substantial contribution to the phenotypic variation, which was overall high for all traits. CONCLUSIONS: There is potential for genetic improvement in all studied traits and estimates of genetic correlations indicate a partly shared genetic architecture among the traits. All estimates have large standard errors. Implementing high-throughput phenotyping is imperative for the estimation of genetic parameters in fast developing insects, and facilitates age synchronisation, which is vital in a breeding population. In spite of endeavours to minimise non-genetic sources of variation, all traits demonstrated substantial influences from environmental components. This emphasises the necessity of thorough attention to the experimental design before breeding is initiated in insect populations.

Investigating genomic prediction strategies for grain carotenoid traits in a tropical/subtropical maize panel.

Vitamin A deficiency remains prevalent on a global scale, including in regions where maize constitutes a high percentage of human diets. One solution for alleviating this deficiency has been to increase grain concentrations of provitamin A carotenoids in maize (Zea mays ssp. mays L.)-an example of biofortification. The International Maize and Wheat Improvement Center (CIMMYT) developed a Carotenoid Association Mapping panel of 380 inbred lines adapted to tropical and subtropical environments that have varying grain concentrations of provitamin A and other health-beneficial carotenoids. Several major genes have been identified for these traits, 2 of which have particularly been leveraged in marker-assisted selection. This project assesses the predictive ability of several genomic prediction strategies for maize grain carotenoid traits within and between 4 environments in Mexico. Ridge Regression-Best Linear Unbiased Prediction, Elastic Net, and Reproducing Kernel Hilbert Spaces had high predictive abilities for all tested traits (ß-carotene, ß-cryptoxanthin, provitamin A, lutein, and zeaxanthin) and outperformed Least Absolute Shrinkage and Selection Operator. Furthermore, predictive abilities were higher when using genome-wide markers rather than only the markers proximal to 2 or 13 genes. These findings suggest that genomic prediction models using genome-wide markers (and assuming equal variance of marker effects) are worthwhile for these traits even though key genes have already been identified, especially if breeding for additional grain carotenoid traits alongside ß-carotene. Predictive ability was maintained for all traits except lutein in between-environment prediction. The TASSEL (Trait Analysis by aSSociation, Evolution, and Linkage) Genomic Selection plugin performed as well as other more computationally intensive methods for within-environment prediction. The findings observed herein indicate the utility of genomic prediction methods for these traits and could inform their resource-efficient implementation in biofortification breeding programs.

The genetic basis of apple shape and size unraveled by digital phenotyping.

Great diversity of shape, size, and skin color is observed among the fruits of different apple genotypes. These traits are critical for consumers and therefore interesting targets for breeding new apple varieties. However, they are difficult to phenotype and their genetic basis, especially for fruit shape and ground color, is largely unknown. We used the FruitPhenoBox to digitally phenotype 525 genotypes of the apple reference population (apple REFPOP) genotyped for 303,148 single nucleotide polymorphism (SNP) markers. From the apple images, 573 highly heritable features describing fruit shape and size as well as 17 highly heritable features for fruit skin color were extracted to explore genotype-phenotype relationships. Out of these features, seven principal components (PCs) and 16 features with the Pearson's correlation r < 0.75 (selected features) were chosen to carry out genome-wide association studies (GWAS) for fruit shape and size. Four PCs and eight selected features were used in GWAS for fruit skin color. In total, 69 SNPs scattered over all 17 apple chromosomes were significantly associated with round, conical, cylindrical, or symmetric fruit shapes and fruit size. Novel associations with major effect on round or conical fruit shapes and fruit size were identified on chromosomes 1 and 2. Additionally, 16 SNPs associated with PCs and selected features related to red overcolor as well as green and yellow ground color were found on eight chromosomes. The identified associations can be used to advance marker-assisted selection in apple fruit breeding to systematically select for desired fruit appearance.

Spatio-temporal modeling of high-throughput multispectral aerial images improves agronomic trait genomic prediction in hybrid maize.

Design randomizations and spatial corrections have increased understanding of genotypic, spatial, and residual effects in field experiments, but precisely measuring spatial heterogeneity in the field remains a challenge. To this end, our study evaluated approaches to improve spatial modeling using high-throughput phenotypes (HTP) via unoccupied aerial vehicle (UAV) imagery. The normalized difference vegetation index was measured by a multispectral MicaSense camera and processed using ImageBreed. Contrasting to baseline agronomic trait spatial correction and a baseline multitrait model, a two-stage approach was proposed. Using longitudinal normalized difference vegetation index data, plot level permanent environment effects estimated spatial patterns in the field throughout the growing season. Normalized difference vegetation index permanent environment were separated from additive genetic effects using 2D spline, separable autoregressive models, or random regression models. The Permanent environment were leveraged within agronomic trait genomic best linear unbiased prediction either modeling an empirical covariance for random effects, or by modeling fixed effects as an average of permanent environment across time or split among three growth phases. Modeling approaches were tested using simulation data and Genomes-to-Fields hybrid maize (Zea mays L.) field experiments in 2015, 2017, 2019, and 2020 for grain yield, grain moisture, and ear height. The two-stage approach improved heritability, model fit, and genotypic effect estimation compared to baseline models. Electrical conductance and elevation from a 2019 soil survey significantly improved model fit, while 2D spline permanent environment were most strongly correlated with the soil parameters. Simulation of field effects demonstrated improved specificity for random regression models. In summary, the use of longitudinal normalized difference vegetation index measurements increased experimental accuracy and understanding of field spatio-temporal heterogeneity.

Harnessing the predicted maize pan-interactome for putative gene function prediction and prioritization of candidate genes for important traits.

The recent assembly and annotation of the 26 maize nested association mapping population founder inbreds have enabled large-scale pan-genomic comparative studies. These studies have expanded our understanding of agronomically important traits by integrating pan-transcriptomic data with trait-specific gene candidates from previous association mapping results. In contrast to the availability of pan-transcriptomic data, obtaining reliable protein-protein interaction (PPI) data has remained a challenge due to its high cost and complexity. We generated predicted PPI networks for each of the 26 genomes using the established STRING database. The individual genome-interactomes were then integrated to generate core- and pan-interactomes. We deployed the PPI clustering algorithm ClusterONE to identify numerous PPI clusters that were functionally annotated using gene ontology (GO) functional enrichment, demonstrating a diverse range of enriched GO terms across different clusters. Additional cluster annotations were generated by integrating gene coexpression data and gene description annotations, providing additional useful information. We show that the functionally annotated PPI clusters establish a useful framework for protein function prediction and prioritization of candidate genes of interest. Our study not only provides a comprehensive resource of predicted PPI networks for 26 maize genomes but also offers annotated interactome clusters for predicting protein functions and prioritizing gene candidates. The source code for the Python implementation of the analysis workflow and a standalone web application for accessing the analysis results are available at https://github.com/eporetsky/PanPPI.

Variation in leaf carbon economics, energy balance, and heat tolerance traits highlights differing timescales of adaptation and acclimation.

Multivariate leaf trait correlations are hypothesized to originate from natural selection on carbon economics traits that control lifetime leaf carbon gain, and energy balance traits governing leaf temperatures, physiological rates, and heat injury. However, it is unclear whether macroevolution of leaf traits primarily reflects selection for lifetime carbon gain or energy balance, and whether photosynthetic heat tolerance is coordinated along these axes. To evaluate these hypotheses, we measured carbon economics, energy balance, and photosynthetic heat tolerance traits for 177 species (157 families) in a common garden that minimizes co-variation of taxa and climate. We observed wide variation in carbon economics, energy balance, and heat tolerance traits. Carbon economics and energy balance (but not heat tolerance) traits were phylogenetically structured, suggesting macroevolution of leaf mass per area and leaf dry matter content reflects selection on carbon gain rather than energy balance. Carbon economics and energy balance traits varied along a common axis orthogonal to heat tolerance traits. Our results highlight a fundamental mismatch in the timescales over which morphological and heat tolerance traits respond to environmental variation. Whereas carbon economics and energy balance traits are constrained by species' evolutionary histories, photosynthetic heat tolerance traits are not and can acclimate readily to leaf microclimates.

Root-associated fungi and acquisitive root traits facilitate permafrost nitrogen uptake from long-term experimentally warmed tundra.

Root-associated fungi (RAF) and root traits regulate plant acquisition of nitrogen (N), which is limiting to growth in Arctic ecosystems. With anthropogenic warming, a new N source from thawing permafrost has the potential to change vegetation composition and increase productivity, influencing climate feedbacks. Yet, the impact of warming on tundra plant root traits, RAF, and access to permafrost N is uncertain. We investigated the relationships between RAF, species-specific root traits, and uptake of N from the permafrost boundary by tundra plants experimentally warmed for nearly three decades at Toolik Lake, Alaska. Warming increased acquisitive root traits of nonmycorrhizal and mycorrhizal plants. RAF community composition of ericoid (ERM) but not ectomycorrhizal (ECM) shrubs was impacted by warming and correlated with root traits. RAF taxa in the dark septate endophyte, ERM, and ECM guilds strongly correlated with permafrost N uptake for ECM and ERM shrubs. Overall, a greater proportion of variation in permafrost N uptake was related to root traits than RAF. Our findings suggest that warming Arctic ecosystems will result in interactions between roots, RAF, and newly thawed permafrost that may strongly impact feedbacks to the climate system through mechanisms of carbon and N cycling.

Pervasive under-dominance in gene expression underlying emergent growth trajectories in Arabidopsis thaliana hybrids.

BACKGROUND: Complex traits, such as growth and fitness, are typically controlled by a very large number of variants, which can interact in both additive and non-additive fashion. In an attempt to gauge the relative importance of both types of genetic interactions, we turn to hybrids, which provide a facile means for creating many novel allele combinations. RESULTS: We focus on the interaction between alleles of the same locus, i.e., dominance, and perform a transcriptomic study involving 141 random crosses between different accessions of the plant model species Arabidopsis thaliana. Additivity is rare, consistently observed for only about 300 genes enriched for roles in stress response and cell death. Regulatory rare-allele burden affects the expression level of these genes but does not correlate with F1 rosette size. Non-additive, dominant gene expression in F1 hybrids is much more common, with the vast majority of genes (over 90%) being expressed below the parental average. Unlike in the additive genes, regulatory rare-allele burden in the dominant gene set is strongly correlated with F1 rosette size, even though it only mildly covaries with the expression level of these genes. CONCLUSIONS: Our study underscores under-dominance as the predominant gene action associated with emergence of rosette growth trajectories in the A. thaliana hybrid model. Our work lays the foundation for understanding molecular mechanisms and evolutionary forces that lead to dominance complementation of rare regulatory alleles.

Sex-specific differences in the genetic and environmental effects on cardiac phenotypic variation assessed by echocardiography.

The drivers of sexual dimorphism in heart failure phenotypes are currently poorly understood. Divergent phenotypes may result from differences in heritability and genetic versus environmental influences on the interplay of cardiac structure and function. To assess sex-specific heritability and genetic versus environmental contributions to variation and inter-relations between echocardiography traits in a large community-based cohort. We studied Framingham Heart Study participants of Offspring Cohort examination 8 (2005-2008) and Third Generation Cohort examination 1 (2002-2005). Five cardiac traits and six functional traits were measured using standardized echocardiography. Sequential Oligogenic Linkage Analysis Routines (SOLAR) software was used to perform singular and bivariate quantitative trait linkage analysis. In our study of 5674 participants (age 49 ± 15 years; 54% women), heritability for all traits was significant for both men and women. There were no significant differences in traits between men and women. Within inter-trait correlations, there were two genetic, and four environmental trait pairs with sex-based differences. Within both significant genetic trait pairs, men had a positive relation, and women had no significant relation. We observed significant sex-based differences in inter-trait genetic and environmental correlations between cardiac structure and function. These findings highlight potential pathways of sex-based divergent heart failure phenotypes.

Genomic evaluation of carcass traits of Korean beef cattle Hanwoo using a single-step marker effect model.

Hanwoo beef cattle are well known for the flavor and tenderness of their meat. Genetic improvement programs have been extremely successful over the last 40 yr. Recently, genomic selection was initiated in Hanwoo to enhance genetic progress. Routine genomic evaluation based on the single-step breeding value model was implemented in 2020 for all economically important traits. In this study, we tested a single-step marker effect model for the genomic evaluation of four carcass traits, namely, carcass weight (CW), eye muscle area, backfat thickness, and marbling score. In total, 8,023,666 animals with carcass records were jointly evaluated, including 29,965 genotyped animals. To assess the prediction stability of the single-step model, carcass data from the last 4 yr were removed in a forward validation study. The estimated genomic breeding values (GEBV) of the validation animals and other animals were compared between the truncated and full evaluations. A parallel conventional best linear unbiased prediction (BLUP) evaluation with either the full or the truncated dataset was also conducted for comparison with the single-step model. The estimates of the marker effect from the truncated evaluation were highly correlated with those from the full evaluation, ranging from 0.88 to 0.92. The regression coefficients of the estimates of the marker effect for the full and truncated evaluations were close to their expected value of 1, indicating unbiased estimates for all carcass traits. Estimates of the marker effect revealed three chromosomal regions (chromosomes 4, 6, and 14) harboring the major genes for CW in Hanwoo. For validation of cows or steers, the single-step model had a much higher R2 value for the linear regression model than the conventional BLUP model. Based on the regression intercept and slope of the validation, the single-step evaluation was neither inflated nor deflated. For genotyped animals, the estimated GEBV from the full and truncated evaluations were more correlated than the estimated breeding values from the two conventional BLUP evaluations. The single-step model provided a more accurate and stable evaluation over time.

Hanwoo beef cattle are well known for the flavor and tenderness of their meat. Genetic improvement programs have been successful over the last 40 yr. Recently, genomic selection was initiated in Hanwoo to enhance genetic progress. A routine genomic evaluation based on the single-step breeding value model was implemented in 2020 for all economically important traits. In this study, we tested a single-step marker effect model for the genomic evaluation of four carcass traits. In total, 8,023,666 cows or steers with carcass records were jointly evaluated, including 29,965 genotyped animals. To assess the prediction accuracy of the single-step model, carcass data from the last 4 yr were removed in a forward validation study. Estimated genomic breeding values (GEBV) of validation animals were compared between truncated and full evaluations. A parallel conventional best linear unbiased prediction (BLUP) evaluation with either the full or truncated dataset was conducted for comparison with the single-step model. Plots of the estimates of the marker effect showed three chromosomal regions harboring the major genes for carcass weight in Hanwoo. The single-step model yielded a more accurate and stable evaluation over time than the conventional BLUP model.

Detection of CD36 gene polymorphism associated with chicken carcass traits and skin yellowness.

Investigations into the association between chicken traits and genetic variations provide helpful breeding information to improve production performance and economic benefits in chickens. The single nucleotide polymorphism technique is an important method in agricultural molecular breeding. In this study, we detected 11 SNPs in the CD36 gene, 2 SNPs (g.-1974 A>G, g.-1888 T>C) located in the 5' flanking regions, 8 SNPs (g.23496 G>A, g.23643 C>T, g.23931 T>C, g.23937 G>A, g.31256 C>A, g.31258 C>T, g.31335 C>T, g.31534 A>C) located in the intron region, 1 SNPs (g.23743 G>T) located in the exon region and it belongs to synonymous mutation. In SNPs g.23743 G>T, the abdominal fat weight and abdominal fat weight rate of the GG genotype were lower than that of the TT genotype. In SNPs g.23931 T>C, the full-bore weight rate and half-bore weight rate of the TT genotype were higher compared with the CC genotype. And the SNPs g.-1888 T>C, g.23496 G>A, g.23643 C>T, g.31335 C>T and g.31534 A>C were significantly associated with skin yellowness traits, the cloacal skin yellowness before slaughter of the TT genotype was higher than that of the TC and CC genotype in SNPs g.-1888 T>C. Furthermore, 3 haplotypes of the above eleven SNPs were calculated and they correlated with heart weight, stomach weight, wing weight, leg skin yellowness and shin skin yellowness before slaughter. Finally, the CD36 expression profile displayed the expression pattern of CD36 mRNA variation in different tissues.

Epistatic QTLs for yield heterosis in tomato.

Controlled population development and genome-wide association studies have proven powerful in uncovering genes and alleles underlying complex traits. An underexplored dimension of such studies is the phenotypic contribution of nonadditive interactions between quantitative trait loci (QTLs). Capturing of such epistasis in a genome-wide manner requires very large populations to represent replicated combinations of loci whose interactions determine phenotypic outcomes. Here, we dissect epistasis using a densely genotyped population of 1,400 backcross inbred lines (BILs) between a modern processing tomato inbred (Solanum lycopersicum) and the Lost Accession (LA5240) of a distant, green-fruited, drought-tolerant wild species, Solanum pennellii. The homozygous BILs, each harboring an average of 11 introgressions and their hybrids with the recurrent parents, were phenotyped for tomato yield components. Population-wide mean yield of the BILs was less than 50% of that of their hybrids (BILHs). All the homozygous introgressions across the genome reduced yield relative to recurrent parent, while several QTLs of the BILHs independently improved productivity. Analysis of two QTL scans showed 61 cases of less-than-additive interactions and 19 cases of more-than-additive interactions. Strikingly, a single epistatic interaction involving S. pennellii QTLs on chromosomes 1 and 7, that independently did not affect yield, increased fruit yield by 20 to 50% in the double introgression hybrid grown in irrigated and dry fields over a period of 4 y. Our work demonstrates the power of large, interspecific controlled population development to uncover hidden QTL phenotypes and how rare epistatic interactions can improve crop productivity via heterosis.