Molecular breeding of flower load related traits in dioecious autotetraploid Actinidia arguta.

Flowering plants exhibit a wide range of sexual reproduction systems, with the majority being hermaphroditic. However, some plants, such as Actinidia arguta (kiwiberry), have evolved into dioecious species with distinct female and male vines. In this study, we investigated the flower load and growth habits of female kiwiberry genotypes to identify the genetic basis of high yield with low maintenance requirements. Owing to the different selection approaches between female and male genotypes, we further extended our study to male kiwiberry genotypes. By combining both investigations, we present a novel breeding tool for dioecious crops. A population of A. arguta seedlings was phenotyped for flower load traits, in particular the proportion of non-floral shoots, proportion of floral shoots, and average number of flowers per floral shoot. Quantitative trait locus (QTL) mapping was used to analyse the genetic basis of these traits. We identified putative QTLs on chromosome 3 associated with flower-load traits. A pleiotropic effect of the male-specific region of the Y chromosome (MSY) on chromosome 3 affecting flower load-related traits between female and male vines was observed in an A. arguta breeding population. Furthermore, we utilized Genomic Best Linear Unbiased Prediction (GBLUP) to predict breeding values for the quantitative traits by leveraging genomic data. This approach allowed us to identify and select superior genotypes. Our findings contribute to the understanding of flowering and fruiting dynamics in Actinidia species, providing insights for kiwiberry breeding programs aiming to improve yield through the utilization of genomic methods and trait mapping. Supplementary Information: The online version contains supplementary material available at 10.1007/s11032-024-01476-7.

AGHmatrix: genetic relationship matrices in R.

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

A multi-omics framework reveals strawberry flavor genes and their regulatory elements.

Flavor is essential to consumer preference of foods and is an increasing focus of plant breeding programs. In fruit crops, identifying genes underlying volatile organic compounds has great promise to accelerate flavor improvement, but polyploidy and heterozygosity in many species have slowed progress. Here we use octoploid cultivated strawberry to demonstrate how genomic heterozygosity, transcriptomic intricacy and fruit metabolomic diversity can be treated as strengths and leveraged to uncover fruit flavor genes and their regulatory elements. Multi-omics datasets were generated including an expression quantitative trait loci map with 196 diverse breeding lines, haplotype-phased genomes of a highly-flavored breeding selection, a genome-wide structural variant map using five haplotypes, and volatile genome-wide association study (GWAS) with > 300 individuals. Overlaying regulatory elements, structural variants and GWAS-linked allele-specific expression of numerous genes to variation in volatile compounds important to flavor. In one example, the functional role of anthranilate synthase alpha subunit 1 in methyl anthranilate biosynthesis was supported via fruit transient gene expression assays. These results demonstrate a framework for flavor gene discovery in fruit crops and a pathway to molecular breeding of cultivars with complex and desirable flavor.

Terpene volatiles mediates the chemical basis of blueberry aroma and consumer acceptability.

Flavor is among the most important traits valued by consumers of fresh fruits. Human perception of flavor occurs primarily through two main sensory inputs, taste and aroma. Through retronasal olfaction, volatile organic compounds (VOCs) emitted by the fruit are able to produce the sensation of aroma which when combined with gustatory inputs from the tongue together underly our perception of the thousands of flavors we experience throughout our lives. In blueberry, breeders have observed that some genotypes possess berries with unique 'floral' and 'sweet' flavor and aroma notes. The potential impact these characteristics might have on consumer acceptability is largely unknown and represents an opportunity to better understand how aroma attributes affect the perception of blueberry flavor. In this study, we dissected the main components of blueberry aroma and associated it with consumer predilections by pairing metabolomics with sensory analysis. Our contribution in this study is four-fold: (i) first, we differentiated genotypes with floral and sweet aroma notes and confirmed that such characteristics are preferred by consumers; (ii) at the chemical level, we showed that a group of eight terpenoid volatiles (p-cymene, myrtenal, linalool, L-carvenol, geranyl acetone, geranyl acetate, D-limonene and ß-myrcene) constitute the primary metabolic group associated with these aroma sensations; (iii) we demonstrated that aromatic genotypes can be classified using metabolomics; and finally, (iv) we combined pedigree and metabolomic information and showed the importance of metabolomic data for flavor-assisted selection. Our findings open new avenues to explore the phenomenon of flavor in blueberries and also allow us to present an emerging view about flavor and provide a detailed blueprint of how this targeted trait could be addressed in fruit and vegetable breeding.

Haplotype reconstruction in connected tetraploid F1 populations.

In diploid species, many multiparental populations have been developed to increase genetic diversity and quantitative trait loci (QTL) mapping resolution. In these populations, haplotype reconstruction has been used as a standard practice to increase the power of QTL detection in comparison with the marker-based association analysis. However, such software tools for polyploid species are few and limited to a single biparental F1 population. In this study, a statistical framework for haplotype reconstruction has been developed and implemented in the software PolyOrigin for connected tetraploid F1 populations with shared parents, regardless of the number of parents or mating design. Given a genetic or physical map of markers, PolyOrigin first phases parental genotypes, then refines the input marker map, and finally reconstructs offspring haplotypes. PolyOrigin can utilize single nucleotide polymorphism (SNP) data coming from arrays or from sequence-based genotyping; in the latter case, bi-allelic read counts can be used (and are preferred) as input data to minimize the influence of genotype calling errors at low depth. With extensive simulation we show that PolyOrigin is robust to the errors in the input genotypic data and marker map. It works well for various population designs with ≥30 offspring per parent and for sequences with read depth as low as 10x. PolyOrigin was further evaluated using an autotetraploid potato dataset with a 3 × 3 half-diallel mating design. In conclusion, PolyOrigin opens up exciting new possibilities for haplotype analysis in tetraploid breeding populations.

QTL mapping in outbred tetraploid (and diploid) diallel populations.

Over the last decade, multiparental populations have become a mainstay of genetics research in diploid species. Our goal was to extend this paradigm to autotetraploids by developing software for quantitative trait locus (QTL) mapping in connected F1 populations derived from a set of shared parents. For QTL discovery, phenotypes are regressed on the dosage of parental haplotypes to estimate additive effects. Statistical properties of the model were explored by simulating half-diallel diploid and tetraploid populations with different population sizes and numbers of parents. Across scenarios, the number of progeny per parental haplotype (pph) largely determined the statistical power for QTL detection and accuracy of the estimated haplotype effects. Multiallelic QTL with heritability 0.2 were detected with 90% probability at 25 pph and genome-wide significance level 0.05, and the additive haplotype effects were estimated with over 90% accuracy. Following QTL discovery, the software enables a comparison of models with multiple QTL and nonadditive effects. To illustrate, we analyzed potato tuber shape in a half-diallel population with three tetraploid parents. A well-known QTL on chromosome 10 was detected, for which the inclusion of digenic dominance lowered the Deviance Information Criterion (DIC) by 17 points compared to the additive model. The final model also contained a minor QTL on chromosome 1, but higher-order dominance and epistatic effects were excluded based on the DIC. In terms of practical impacts, the software is already being used to select offspring based on the effect and dosage of particular haplotypes in breeding programs.

Genomic Selection in an Outcrossing Autotetraploid Fruit Crop: Lessons From Blueberry Breeding.

Blueberry (Vaccinium corymbosum and hybrids) is a specialty crop with expanding production and consumption worldwide. The blueberry breeding program at the University of Florida (UF) has greatly contributed to expanding production areas by developing low-chilling cultivars better adapted to subtropical and Mediterranean climates of the globe. The breeding program has historically focused on recurrent phenotypic selection. As an autopolyploid, outcrossing, perennial, long juvenile phase crop, blueberry breeding cycles are costly and time consuming, which results in low genetic gains per unit of time. Motivated by applying molecular markers for a more accurate selection in the early stages of breeding, we performed pioneering genomic selection studies and optimization for its implementation in the blueberry breeding program. We have also addressed some complexities of sequence-based genotyping and model parametrization for an autopolyploid crop, providing empirical contributions that can be extended to other polyploid species. We herein revisited some of our previous genomic selection studies and showed for the first time its application in an independent validation set. In this paper, our contribution is three-fold: (i) summarize previous results on the relevance of model parametrizations, such as diploid or polyploid methods, and inclusion of dominance effects; (ii) assess the importance of sequence depth of coverage and genotype dosage calling steps; (iii) demonstrate the real impact of genomic selection on leveraging breeding decisions by using an independent validation set. Altogether, we propose a strategy for using genomic selection in blueberry, with the potential to be applied to other polyploid species of a similar background.

High-Resolution Linkage Map and QTL Analyses of Fruit Firmness in Autotetraploid Blueberry.

Blueberry (Vaccinium corymbosum and hybrids) is an autotetraploid crop whose commercial relevance has been growing steadily during the last 20 years. However, the ever-increasing cost of labor for hand-picking blueberry is one main constraint in competitive marketing of the fruit. Machine harvestability is, therefore, a key trait for the blueberry industry. Understanding the genetic architecture of traits related to machine harvestability through Quantitative Trait Loci (QTL) mapping is the first step toward implementation of molecular breeding for faster genetic gains. Despite recent advances in software development for autotetraploid genetic mapping, a high-resolution map is still not available for blueberry. In this study, we crafted a map for autotetraploid low-chill highbush blueberry containing 11,292 SNP markers and a total size of 1,953.97 cM (average density of 5.78 markers/cM). This map was subsequently used to perform QTL analyses in 2-year field trials for a trait crucial to machine harvesting: fruit firmness. Preliminary insights were also sought for single evaluations of firmness retention after cold storage, and fruit detachment force traits. Significant QTL peaks were identified for all the traits and overlapping QTL intervals were detected for firmness across the years. We found low-to-moderate QTL effects explaining the phenotypic variance, which suggest a quantitative nature of these traits. The QTL intervals were further speculated for putative gene repertoire. Altogether, our findings provide the basis for future fine-mapping and molecular breeding efforts for machine harvesting in blueberry.

Estimation of Molecular Pairwise Relatedness in Autopolyploid Crops.

A suitable pairwise relatedness estimation is key to genetic studies. Several methods are proposed to compute relatedness in autopolyploids based on molecular data. However, unlike diploids, autopolyploids still need further studies considering scenarios with many linked molecular markers with known dosage. In this study, we provide guidelines for plant geneticists and breeders to access trustworthy pairwise relatedness estimates. To this end, we simulated populations considering different ploidy levels, meiotic pairings patterns, number of loci and alleles, and inbreeding levels. Analysis were performed to access the accuracy of distinct methods and to demonstrate the usefulness of molecular marker in practical situations. Overall, our results suggest that at least 100 effective biallelic molecular markers are required to have good pairwise relatedness estimation if methods based on correlation is used. For this number of loci, current methods based on multiallelic markers show lower performance than biallelic ones. To estimate relatedness in cases of inbreeding or close relationships (as parent-offspring, full-sibs, or half-sibs) is more challenging. Methods to estimate pairwise relatedness based on molecular markers, for different ploidy levels or pedigrees were implemented in the AGHmatrix R package.

Signatures of selection for resistance to Haemonchus contortus in sheep and goats.

BACKGROUND: Gastrointestinal nematode infection (GNI) is the most important disease affecting the small ruminant industry in U.S. The environmental conditions in the southern United States are ideal for the survival of the most pathogenic gastrointestinal nematode, Haemonchus contortus. Host genetic variation for resistance to H. contortus allows selective breeding for increased resistance of animals. This selection process increases the prevalence of particular alleles in sheep and goats and creates unique genetic patterns in the genome of these species. The aim of this study was to identify loci with divergent allelic frequencies in a candidate gene panel of 100 genes using two different approaches (frequentist and Bayesian) to estimate Fst outliers in three different breeds of sheep and goats exposed to H. contortus. RESULTS: Our results for sheep populations showed SNPs under selection in C3AR1, CSF3, SOCS2, NOS2, STAT5B, TGFB2 and IL2RA genes using frequentist and Bayesian approaches. For goats, SNPs in CD1D, ITGA9, IL12A, IL13RA1, CD86 and TGFB2 genes were under selection. Common signatures of selection in both species were observed in NOS2, TGFB2 and TLR4 genes. Directional selection was present in all SNPs evaluated in the present study. CONCLUSIONS: A total of 13 SNPs within 7 genes of our candidate gene panel related to H. contortus exposure were identified under selection in sheep populations. For goats, 11 SNPs within 7 genes were identified under selection. Results from this study support the hypothesis that resistance to H. contortus is likely to be controlled by many loci. Shared signatures of selection related to mechanisms of immune protection against H. contortus infection in sheep and goats could be useful targets in breeding programs aimed to produce resistant animals with low FEC.

Genomic Selection with Allele Dosage in Panicum maximum Jacq.

Genomic selection is an efficient approach to get shorter breeding cycles in recurrent selection programs and greater genetic gains with selection of superior individuals. Despite advances in genotyping techniques, genetic studies for polyploid species have been limited to a rough approximation of studies in diploid species. The major challenge is to distinguish the different types of heterozygotes present in polyploid populations. In this work, we evaluated different genomic prediction models applied to a recurrent selection population of 530 genotypes of Panicum maximum, an autotetraploid forage grass. We also investigated the effect of the allele dosage in the prediction, i.e., considering tetraploid (GS-TD) or diploid (GS-DD) allele dosage. A longitudinal linear mixed model was fitted for each one of the six phenotypic traits, considering different covariance matrices for genetic and residual effects. A total of 41,424 genotyping-by-sequencing markers were obtained using 96-plex and Pst1 restriction enzyme, and quantitative genotype calling was performed. Six predictive models were generalized to tetraploid species and predictive ability was estimated by a replicated fivefold cross-validation process. GS-TD and GS-DD models were performed considering 1,223 informative markers. Overall, GS-TD data yielded higher predictive abilities than with GS-DD data. However, different predictive models had similar predictive ability performance. In this work, we provide bioinformatic and modeling guidelines to consider tetraploid dosage and observed that genomic selection may lead to additional gains in recurrent selection program of P. maximum.

How can a high-quality genome assembly help plant breeders?

The decreasing costs of next-generation sequencing and the improvements in de novo sequence assemblers have made it possible to obtain reference genomes for most eukaryotes, including minor crops such as the blueberry (Vaccinium corymbosum). Nevertheless, these genomes are at various levels of completeness and few have been anchored to chromosome scale and/or are haplotype-phased. We highlight the impact of a high-quality genome assembly for plant breeding and genetic research by showing how it affects our understanding of the genetic architecture of important traits and aids marker selection and candidate gene detection. We compared the results of genome-wide association studies and genomic selection that were already published using a blueberry draft genome as reference with the results using the recent released chromosome-scale and haplotype-phased blueberry genome. We believe that the benefits shown herein reinforce the importance of genome assembly projects for other non-model species.

Genomic Prediction of Autotetraploids; Influence of Relationship Matrices, Allele Dosage, and Continuous Genotyping Calls in Phenotype Prediction.

Estimation of allele dosage, using genomic data, in autopolyploids is challenging and current methods often result in the misclassification of genotypes. Some progress has been made when using SNP arrays, but the major challenge is when using next generation sequencing data. Here we compare the use of read depth as continuous parameterization with ploidy parameterizations in the context of genomic selection (GS). Additionally, different sources of information to build relationship matrices were compared. A real breeding population of the autotetraploid species blueberry (Vaccinium corybosum), composed of 1,847 individuals was phenotyped for eight yield and fruit quality traits over two years. Continuous genotypic based models performed as well as the best models. This approach also reduces the computational time and avoids problems associated with misclassification of genotypic classes when assigning dosage in polyploid species. This approach could be very valuable for species with higher ploidy levels or for emerging crops where ploidy is not well understood. To our knowledge, this work constitutes the first study of genomic selection in blueberry. Accuracies are encouraging for application of GS for blueberry breeding. GS could reduce the time for cultivar release by three years, increasing the genetic gain per cycle by 86% on average when compared to phenotypic selection, and 32% when compared with pedigree-based selection. Finally, the genotypic and phenotypic data used in this study are made available for comparative analysis of dosage calling and genomic selection prediction models in the context of autopolyploids.

High-Resolution Genetic Map and QTL Analysis of Growth-Related Traits of Hevea brasiliensis Cultivated Under Suboptimal Temperature and Humidity Conditions.

Rubber tree (Hevea brasiliensis) cultivation is the main source of natural rubber worldwide and has been extended to areas with suboptimal climates and lengthy drought periods; this transition affects growth and latex production. High-density genetic maps with reliable markers support precise mapping of quantitative trait loci (QTL), which can help reveal the complex genome of the species, provide tools to enhance molecular breeding, and shorten the breeding cycle. In this study, QTL mapping of the stem diameter, tree height, and number of whorls was performed for a full-sibling population derived from a GT1 and RRIM701 cross. A total of 225 simple sequence repeats (SSRs) and 186 single-nucleotide polymorphism (SNP) markers were used to construct a base map with 18 linkage groups and to anchor 671 SNPs from genotyping by sequencing (GBS) to produce a very dense linkage map with small intervals between loci. The final map was composed of 1,079 markers, spanned 3,779.7 cM with an average marker density of 3.5 cM, and showed collinearity between markers from previous studies. Significant variation in phenotypic characteristics was found over a 59-month evaluation period with a total of 38 QTLs being identified through a composite interval mapping method. Linkage group 4 showed the greatest number of QTLs (7), with phenotypic explained values varying from 7.67 to 14.07%. Additionally, we estimated segregation patterns, dominance, and additive effects for each QTL. A total of 53 significant effects for stem diameter were observed, and these effects were mostly related to additivity in the GT1 clone. Associating accurate genome assemblies and genetic maps represents a promising strategy for identifying the genetic basis of phenotypic traits in rubber trees. Then, further research can benefit from the QTLs identified herein, providing a better understanding of the key determinant genes associated with growth of Hevea brasiliensis under limiting water conditions.

Metabolite Profiles of Sugarcane Culm Reveal the Relationship Among Metabolism and Axillary Bud Outgrowth in Genetically Related Sugarcane Commercial Cultivars.

Metabolic composition is known to exert influence on several important agronomic traits, and metabolomics, which represents the chemical composition in a cell, has long been recognized as a powerful tool for bridging phenotype-genotype interactions. In this work, sixteen truly representative sugarcane Brazilian varieties were selected to explore the metabolic networks in buds and culms, the tissues involved in the vegetative propagation of this species. Due to the fact that bud sprouting is a key trait determining crop establishment in the field, the sprouting potential among the genotypes was evaluated. The use of partial least square discriminant analysis indicated only mild differences on bud outgrowth potential under controlled environmental conditions. However, primary metabolite profiling provided information on the variability of metabolic features even under a narrow genetic background, typical for modern sugarcane cultivars. Metabolite-metabolite correlations within and between tissues revealed more complex patterns for culms in relation to buds, and enabled the recognition of key metabolites (e.g., sucrose, putrescine, glutamate, serine, and myo-inositol) affecting sprouting ability. Finally, those results were associated with the genetic background of each cultivar, showing that metabolites can be potentially used as indicators for the genetic background.

AGHmatrix: R Package to Construct Relationship Matrices for Autotetraploid and Diploid Species: A Blueberry Example.

Progress in the rate of improvement in autopolyploid species has been limited compared with diploids, mainly because software and methods to apply advanced prediction and selection methodologies in autopolyploids are lacking. The objectives of this research were to (i) develop an R package for autopolyploids to construct the relationship matrix derived from pedigree information that accounts for autopolyploidy and double reduction and (ii) use the package to estimate the level and effect of double reduction in an autotetraploid blueberry breeding population with extensive pedigree information. The package is unique, as it can create A-matrices for different levels of ploidy and double reduction, which can then be used by breeders to fit mixed models in the context of predicting breeding values (BVs). Using the data from this blueberry population, we found for all the traits that tetrasomic inheritance creates a better fit than disomic inheritance. In one of the five traits studied, the level of double reduction was different from zero, decreasing the estimated heritability, but it did not affect the prediction of BVs. We also discovered that different depths of pedigree would have significant implications on the estimation of double reduction using this approach. This freely available R package is available for autopolyploid breeders to estimate the level of double reduction present in their populations and the impact in the estimation of genetic parameters as well as to use advanced methods of prediction and selection.