Genomic selection for fruit quality traits in apple (Malus×domestica Borkh.)

Genomic prediction of resistance to Hymenoscyphus fraxineus in common ash (Fraxinus excelsior L.) populations

The increase in introduced insect pests and pathogens due to anthropogenic environmental changes has become a major concern for tree species worldwide. Common ash (Fraxinus excelsior L.) is one of such species facing a significant threat from the invasive fungal pathogen Hymenoscyphus fraxineus. Some studies have indicated that the susceptibility of ash to the pathogen is genetically determined, providing some hope for accelerated breeding programs that are aimed at increasing the resistance of ash populations. To address this challenge, we used a genomic selection strategy to identify potential genetic markers that are associated with resistance to the pathogen causing ash dieback. Through genome-wide association studies (GWAS) of 300 common ash individuals from 30 populations across Poland (ddRAD, dataset A), we identified six significant SNP loci with a p-value ≤1 × 10-4 associated with health status. To further evaluate the effectiveness of GWAS markers in predicting health status, we considered two genomic prediction scenarios. Firstly, we conducted cross-validation on dataset A. Secondly, we trained markers on dataset A and tested them on dataset B, which involved whole-genome sequencing of 20 individuals from two populations. Genomic prediction analysis revealed that the top SNPs identified via GWAS exhibited notably higher prediction accuracies compared to randomly selected SNPs, particularly with a larger number of SNPs. Cross-validation analyses using dataset A showcased high genomic prediction accuracy, predicting tree health status with over 90% accuracy across the top SNP sets ranging from 500 to 10,000 SNPs from the GWAS datasets. However, no significant results emerged for health status when the model trained on dataset A was tested on dataset B. Our findings illuminate potential genetic markers associated with resistance to ash dieback, offering support for future breeding programs in Poland aimed at combating ash dieback and bolstering conservation efforts for this invaluable tree species.

Advancing tree genomics to future proof next generation orchard production

The challenges facing tree orchard production in the coming years will be largely driven by changes in the climate affecting the sustainability of farming practices in specific geographical regions. Identifying key traits that enable tree crops to modify their growth to varying environmental conditions and taking advantage of new crop improvement opportunities and technologies will ensure the tree crop industry remains viable and profitable into the future. In this review article we 1) outline climate and sustainability challenges relevant to horticultural tree crop industries, 2) describe key tree crop traits targeted for improvement in agroecosystem productivity and resilience to environmental change, and 3) discuss existing and emerging genomic technologies that provide opportunities for industries to future proof the next generation of orchards.

Worldwide productivity and research trend on fruit quality: a bibliometric study

Introduction

As one of the important sources of food for human beings, fruits have been extensively studied. To better guide basic and applied research, it is urgent to conduct a systematic analysis of these studies based on extensive literature collection.

Methods

Based on the Web of Science Core Collection database, this study uses R language and CiteSpace to conduct bibliometric analysis and data mining on the literatures related to fruit quality from January 2013 to June 2023.

Results

The results indicated that among various fruits, tomatoes have been most frequently studied with special interests in photosynthesis, fruit development, and molecular breeding. The research direction primarily focused on fruit resistance and storage characteristics. Among the indicators related to fruit quality, antioxidant activity has the highest co-occurrence with other indicators of fruit quality, especially with nutrients such as anthocyanins, phenolic substances, sugars, and fruit firmness.

Discussion

Currently, adaptation to stress and antioxidant activity are recognized as prominent research focal points in this field. Fruit morphology, particularly fruit size, irrigation methods, application of molecular technology, and infection prevention, represent potential areas of interests in future research on fruit quality.

Mate selection: A useful approach to maximize genetic gain and control inbreeding in genomic and conventional oil palm (Elaeis guineensis Jacq.) hybrid breeding

Genomic selection (GS) is an effective method for the genetic improvement of complex traits in plants and animals. Optimization approaches could be used in conjunction with GS to further increase its efficiency and to limit inbreeding, which can increase faster with GS. Mate selection (MS) typically uses a metaheuristic optimization algorithm, simulated annealing, to optimize the selection of individuals and their matings. However, in species with long breeding cycles, this cannot be studied empirically. Here, we investigated this aspect with forward genetic simulations on a high-performance computing cluster and massively parallel computing, considering the oil palm hybrid breeding example. We compared MS and simple methods of inbreeding management (limitation of the number of individuals selected per family, prohibition of self-fertilization and combination of these two methods), in terms of parental inbreeding and genetic progress over four generations of genomic selection and phenotypic selection. The results showed that, compared to the conventional method without optimization, MS could lead to significant decreases in inbreeding and increases in annual genetic progress, with the magnitude of the effect depending on MS parameters and breeding scenarios. The optimal solution retained by MS differed by five breeding characteristics from the conventional solution: selected individuals covering a broader range of genetic values, fewer individuals selected per full-sib family, decreased percentage of selfings, selfings preferentially made on the best individuals and unbalanced number of crosses among selected individuals, with the better an individual, the higher the number of times he is mated. Stronger slowing-down in inbreeding could be achieved with other methods but they were associated with a decreased genetic progress. We recommend that breeders use MS, with preliminary analyses to identify the proper parameters to reach the goals of the breeding program in terms of inbreeding and genetic gain.

Marker-Assisted Selection in Breeding for Fruit Trait Improvement: A Review

Breeding fruit species is time-consuming and expensive. With few exceptions, trees are likely the worst species to work with in terms of genetics and breeding. Most are characterized by large trees, long juvenile periods, and intensive agricultural practice, and environmental variability plays an important role in the heritability evaluations of every single important trait. Although vegetative propagation allows for the production of a significant number of clonal replicates for the evaluation of environmental effects and genotype × environment interactions, the spaces required for plant cultivation and the intensity of work necessary for phenotypic surveys slow down the work of researchers. Fruit breeders are very often interested in fruit traits: size, weight, sugar and acid content, ripening time, fruit storability, and post-harvest practices, among other traits relevant to each individual species. The translation of trait loci and whole-genome sequences into diagnostic genetic markers that are effective and affordable for use by breeders, who must choose genetically superior parents and subsequently choose genetically superior individuals among their progeny, is one of the most difficult tasks still facing tree fruit geneticists. The availability of updated sequencing techniques and powerful software tools offered the opportunity to mine tens of fruit genomes to find out sequence variants potentially useful as molecular markers. This review is devoted to analysing what has been the role of molecular markers in assisting breeders in selection processes, with an emphasis on the fruit traits of the most important fruit crops for which examples of trustworthy molecular markers have been developed, such as the MDo.chr9.4 marker for red skin colour in apples, the CCD4-based marker CPRFC1, and LG3_13.146 marker for flesh colour in peaches, papayas, and cherries, respectively.

Deducing genotypes for loci of interest from SNP array data via haplotype sharing, demonstrated for apple and cherry

Breeders, collection curators, and other germplasm users require genetic information, both genome-wide and locus-specific, to effectively manage their genetically diverse plant material. SNP arrays have become the preferred platform to provide genome-wide genetic profiles for elite germplasm and could also provide locus-specific genotypic information. However, genotypic information for loci of interest such as those within PCR-based DNA fingerprinting panels and trait-predictive DNA tests is not readily extracted from SNP array data, thus creating a disconnect between historic and new data sets. This study aimed to establish a method for deducing genotypes at loci of interest from their associated SNP haplotypes, demonstrated for two fruit crops and three locus types: quantitative trait loci Ma and Ma3 for acidity in apple, apple fingerprinting microsatellite marker GD12, and Mendelian trait locus Rf for sweet cherry fruit color. Using phased data from an apple 8K SNP array and sweet cherry 6K SNP array, unique haplotypes spanning each target locus were associated with alleles of important breeding parents. These haplotypes were compared via identity-by-descent (IBD) or identity-by-state (IBS) to haplotypes present in germplasm important to U.S. apple and cherry breeding programs to deduce target locus alleles in this germplasm. While IBD segments were confidently tracked through pedigrees, confidence in allele identity among IBS segments used a shared length threshold. At least one allele per locus was deduced for 64-93% of the 181 individuals. Successful validation compared deduced Rf and GD12 genotypes with reported and newly obtained genotypes. Our approach can efficiently merge and expand genotypic data sets, deducing missing data and identifying errors, and is appropriate for any crop with SNP array data and historic genotypic data sets, especially where linkage disequilibrium is high. Locus-specific genotypic information extracted from genome-wide SNP data is expected to enhance confidence in management of genetic resources.

Fruit growth and development in apple: a molecular, genomics and epigenetics perspective

Fruit growth and development are physiological processes controlled by several internal and external factors. This complex regulatory mechanism comprises a series of events occurring in a chronological order over a growing season. Understanding the underlying mechanism of fruit development events, however, requires consideration of the events occurring prior to fruit development such as flowering, pollination, fertilization, and fruit set. Such events are interrelated and occur in a sequential order. Recent advances in high-throughput sequencing technology in conjunction with improved statistical and computational methods have empowered science to identify some of the major molecular components and mechanisms involved in the regulation of fruit growth and have supplied encouraging successes in associating genotypic differentiation with phenotypic observations. As a result, multiple approaches have been developed to dissect such complex regulatory machinery and understand the genetic basis controlling these processes. These methods include transcriptomic analysis, quantitative trait loci (QTLs) mapping, whole-genome approach, and epigenetics analyses. This review offers a comprehensive overview of the molecular, genomic and epigenetics perspective of apple fruit growth and development that defines the final fruit size and provides a detailed analysis of the mechanisms by which fruit growth and development are controlled. Though the main emphasis of this article is on the molecular, genomic and epigenetics aspects of fruit growth and development, we will also deliver a brief overview on events occurring prior to fruit growth.

A single QTL harboring multiple genetic variations leads to complicated phenotypic segregation in apple flesh firmness and crispness

Within a QTL, the genetic recombination and interactions among five and two functional variations at MdbHLH25 and MdWDR5A caused much complicated phenotype segregation in apple FFR and FCR. The storability of climacteric fruit like apple is a quantitative trait. We previously identified 62 quantitative trait loci (QTLs) associating flesh firmness retainability (FFR) and flesh crispness retainability (FCR), but only a few functional genetic variations were identified and validated. The genetic variation network controlling fruit storability is far to be understood and diagnostic markers are needed for molecular breeding. We previously identified overlapped QTLs F16.1/H16.2 for FFR and FCR using an F1 population derived from 'Zisai Pearl' × 'Red Fuji'. In this study, five and two single-nucleotide polymorphisms (SNPs) were identified on the candidate genes MdbHLH25 and MdWDR5A within the QTL region. The SNP1 A allele at MdbHLH25 promoter reduced the expression and SNP2 T allele and/or SNP4/5 GT alleles at the exons attenuated the function of MdbHLH25 by downregulating the expression of the target genes MdACS1, which in turn led to a reduction in ethylene production and maintenance of higher flesh crispness. The SNPs did not alter the protein-protein interaction between MdbHLH25 and MdWDR5A. The joint effect of SNP genotype combinations by the SNPs on MdbHLH25 (SNP1, SNP2, and SNP4) and MdWDR5A (SNPi and SNPii) led to a much broad spectrum of phenotypic segregation in FFR and FCR. Together, the dissection of these genetic variations contributes to understanding the complicated effects of a QTL and provides good potential for marker development in molecular breeding.

The influence of genetic structure on phenotypic diversity in the Australian mango (Mangifera indica) gene pool

Genomic selection is a promising breeding technique for tree crops to accelerate the development of new cultivars. However, factors such as genetic structure can create spurious associations between genotype and phenotype due to the shared history between populations with different trait values. Genetic structure can therefore reduce the accuracy of the genotype to phenotype map, a fundamental requirement of genomic selection models. Here, we employed 272 single nucleotide polymorphisms from 208 Mangifera indica accessions to explore whether the genetic structure of the Australian mango gene pool explained variation in trunk circumference, fruit blush colour and intensity. Multiple population genetic analyses indicate the presence of four genetic clusters and show that the most genetically differentiated cluster contains accessions imported from Southeast Asia (mainly those from Thailand). We find that genetic structure was strongly associated with three traits: trunk circumference, fruit blush colour and intensity in M. indica. This suggests that the history of these accessions could drive spurious associations between loci and key mango phenotypes in the Australian mango gene pool. Incorporating such genetic structure in associations between genotype and phenotype can improve the accuracy of genomic selection, which can assist the future development of new cultivars.

Breeding and genetics of disease resistance in temperate fruit trees: challenges and new opportunities

Climate change, large monocultures of disease-susceptible cultivars, overuse of pesticides, and the emergence of new pathogens or pathogenic strains causing economic losses are all major threats to our environment, health, food, and nutritional supply. Temperate tree fruit crops belonging to the Rosaceae family are the most economically important and widely grown fruit crops. These long-lived crops are under attack from many different pathogens, incurring major economic losses. Multiple chemical sprays to control various diseases annually is a common practice, resulting in significant input costs, as well as environmental and health concerns. Breeding for disease resistance has been undertaken primarily in pome fruit crops (apples and pears) for a few fungal and bacterial diseases, and to a lesser extent in some stone fruit crops. These breeding efforts have taken multiple decades due to the biological constraints and complex genetics of these tree fruit crops. Over the past couple of decades, major advances have been made in genetic and physical mapping, genomics, biotechnology, genome sequencing, and phenomics, along with accumulation of large germplasm collections in repositories. These valuable resources offer opportunities to make significant advances in greatly reducing the time needed to either develop new cultivars or modify existing economic cultivars for enhanced resistance to multiple diseases. This review will cover current knowledge, challenges, and opportunities in breeding for disease resistance in temperate tree fruit crops.

Multi-environment genomic prediction for soluble solids content in peach (Prunus persica)

Genotype-by-environment interaction (G × E) is a common phenomenon influencing genetic improvement in plants, and a good understanding of this phenomenon is important for breeding and cultivar deployment strategies. However, there is little information on G × E in horticultural tree crops, mostly due to evaluation costs, leading to a focus on the development and deployment of locally adapted germplasm. Using sweetness (measured as soluble solids content, SSC) in peach/nectarine assessed at four trials from three US peach-breeding programs as a case study, we evaluated the hypotheses that (i) complex data from multiple breeding programs can be connected using GBLUP models to improve the knowledge of G × E for breeding and deployment and (ii) accounting for a known large-effect quantitative trait locus (QTL) improves the prediction accuracy. Following a structured strategy using univariate and multivariate models containing additive and dominance genomic effects on SSC, a model that included a previously detected QTL and background genomic effects was a significantly better fit than a genome-wide model with completely anonymous markers. Estimates of an individual's narrow-sense and broad-sense heritability for SSC were high (0.57-0.73 and 0.66-0.80, respectively), with 19-32% of total genomic variance explained by the QTL. Genome-wide dominance effects and QTL effects were stable across environments. Significant G × E was detected for background genome effects, mostly due to the low correlation of these effects across seasons within a particular trial. The expected prediction accuracy, estimated from the linear model, was higher than the realised prediction accuracy estimated by cross-validation, suggesting that these two parameters measure different qualities of the prediction models. While prediction accuracy was improved in some cases by combining data across trials, particularly when phenotypic data for untested individuals were available from other trials, this improvement was not consistent. This study confirms that complex data can be combined into a single analysis using GBLUP methods to improve understanding of G × E and also incorporate known QTL effects. In addition, the study generated baseline information to account for population structure in genomic prediction models in horticultural crop improvement.

GWAS provides new insights into the genetic mechanisms of phytochemicals production and red skin colour in apple

Understanding the genetic architecture of apple phytochemicals, and their interplay with conventional selection traits, is critical for the development of new apple cultivars with enhanced health benefits. Apple accessions (n = 344) used for this genome-wide association study (GWAS) represented the wide diversity of metabolic profiles in the domesticated and wild Malus genepools. Fruit samples were phenotyped for 34 metabolites, including a stable vitamin C glycoside "ascorbic acid 2-β-glucoside" (AA-2βG), and the accessions were genotyped using the Apple 20 K SNP Array. Several fruit quality traits, including red skin over-colour (OCOL), were also assessed. Wild Malus accessions showed at least 2-fold higher average content of several metabolites (e.g. ascorbic acid, chlorogenic acid, phloridzin, and trilobatin) than Malus domestica accessions. Several new genomic regions and potential candidate genes underpinning the genetic diversity of apple phytochemicals were identified. The percentage of phenotypic variance explained by the best SNP ranged between 3% and 21% for the different metabolites. Novel association signals for OCOL in the syntenic regions on chromosomes 13 and 16 suggested that whole genome duplication has played a role in the evolution of apple red skin colour. Genetic correlations between phytochemicals and sensory traits were moderate. This study will assist in the selection of Malus accessions with specific phytochemical profiles to establish innovative genomics-based breeding strategies for the development of apple cultivars with enhanced nutritional value.

Construction of Genetic Linkage Map and Mapping QTL Specific to Leaf Anthocyanin Colouration in Mapping Population 'Allahabad Safeda' × 'Purple Guava (Local)' of Guava (Psidium guajava L.)

In the present investigation, F1 hybrids were developed in guava (Psidium guajava L.) by crossing high leaf-anthocyanin reflective-index (ARI1) content cultivars purple guava (local) 'PG', 'CISH G-1' and low leaf-ARI1 content cultivar Seedless 'SL' with Allahabad Safeda 'AS'. On the basis of phenotypic observations, high ARI1 content was observed in the cross 'AS' × 'PG' (0.214). Further, an SSR-markers-based genetic linkage map was developed from a mapping population of 238 F1 individuals derived from cross 'AS' × 'PG'. The linkage map comprised 11 linkage groups (LGs), spanning 1601.7 cM with an average marker interval distance of 29.61 cM between adjacent markers. Five anthocyanin-content related gene-specific markers from apple were tested for parental polymorphism in the genotypes 'AS' and 'PG'. Subsequently, a marker, viz., 'MdMYB10F1', revealed a strong association with leaf anthocyanin content in the guava mapping population. QTL (qARI-6-1) on LG6 explains much of the variation (PVE = 11.51% with LOD = 4.67) in levels of leaf anthocyanin colouration. This is the first report of amplification/utilization of apple anthocyanin-related genes in guava. The genotypic data generated from the genetic map can be further exploited in future for the enrichment of linkage maps and for identification of complex quantitative trait loci (QTLs) governing economically important fruit quality traits in guava.

Research Progress on Genetic Basis of Fruit Quality Traits in Apple (Malus × domestica)

Identifying the genetic variation characteristics of phenotypic traits is important for fruit tree breeding. During the long-term evolution of fruit trees, gene recombination and natural mutation have resulted in a high degree of heterozygosity. Apple (Malus × domestica Borkh.) shows strong ecological adaptability and is widely cultivated, and is among the most economically important fruit crops worldwide. However, the high level of heterozygosity and large genome of apple, in combination with its perennial life history and long juvenile phase, complicate investigation of the genetic basis of fruit quality traits. With continuing augmentation in the apple genomic resources available, in recent years important progress has been achieved in research on the genetic variation of fruit quality traits. This review focuses on summarizing recent genetic studies on apple fruit quality traits, including appearance, flavor, nutritional, ripening, and storage qualities. In addition, we discuss the mapping of quantitative trait loci, screening of molecular markers, and mining of major genes associated with fruit quality traits. The overall aim of this review is to provide valuable insights into the mechanisms of genetic variation and molecular breeding of important fruit quality traits in apple.

Genetic architecture and genomic predictive ability of apple quantitative traits across environments

Implementation of genomic tools is desirable to increase the efficiency of apple breeding. Recently, the multi-environment apple reference population (apple REFPOP) proved useful for rediscovering loci, estimating genomic predictive ability, and studying genotype by environment interactions (G × E). So far, only two phenological traits were investigated using the apple REFPOP, although the population may be valuable when dissecting genetic architecture and reporting predictive abilities for additional key traits in apple breeding. Here we show contrasting genetic architecture and genomic predictive abilities for 30 quantitative traits across up to six European locations using the apple REFPOP. A total of 59 stable and 277 location-specific associations were found using GWAS, 69.2% of which are novel when compared with 41 reviewed publications. Average genomic predictive abilities of 0.18-0.88 were estimated using main-effect univariate, main-effect multivariate, multi-environment univariate, and multi-environment multivariate models. The G × E accounted for up to 24% of the phenotypic variability. This most comprehensive genomic study in apple in terms of trait-environment combinations provided knowledge of trait biology and prediction models that can be readily applied for marker-assisted or genomic selection, thus facilitating increased breeding efficiency.

Genomic Selection: A Tool for Accelerating the Efficiency of Molecular Breeding for Development of Climate-Resilient Crops

Since the inception of the theory and conceptual framework of genomic selection (GS), extensive research has been done on evaluating its efficiency for utilization in crop improvement. Though, the marker-assisted selection has proven its potential for improvement of qualitative traits controlled by one to few genes with large effects. Its role in improving quantitative traits controlled by several genes with small effects is limited. In this regard, GS that utilizes genomic-estimated breeding values of individuals obtained from genome-wide markers to choose candidates for the next breeding cycle is a powerful approach to improve quantitative traits. In the last two decades, GS has been widely adopted in animal breeding programs globally because of its potential to improve selection accuracy, minimize phenotyping, reduce cycle time, and increase genetic gains. In addition, given the promising initial evaluation outcomes of GS for the improvement of yield, biotic and abiotic stress tolerance, and quality in cereal crops like wheat, maize, and rice, prospects of integrating it in breeding crops are also being explored. Improved statistical models that leverage the genomic information to increase the prediction accuracies are critical for the effectiveness of GS-enabled breeding programs. Study on genetic architecture under drought and heat stress helps in developing production markers that can significantly accelerate the development of stress-resilient crop varieties through GS. This review focuses on the transition from traditional selection methods to GS, underlying statistical methods and tools used for this purpose, current status of GS studies in crop plants, and perspectives for its successful implementation in the development of climate-resilient crops.

Identifying and testing marker-trait associations for growth and phenology in three pine species: Implications for genomic prediction

In tree species, genomic prediction offers the potential to forecast mature trait values in early growth stages, if robust marker-trait associations can be identified. Here we apply a novel multispecies approach using genotypes from a new genotyping array, based on 20,795 single nucleotide polymorphisms (SNPs) from three closely related pine species (Pinus sylvestris, Pinus uncinata and Pinus mugo), to test for associations with growth and phenology data from a common garden study. Predictive models constructed using significantly associated SNPs were then tested and applied to an independent multisite field trial of P. sylvestris and the capability to predict trait values was evaluated. One hundred and eighteen SNPs showed significant associations with the traits in the pine species. Common SNPs (MAF > 0.05) associated with bud set were only found in genes putatively involved in growth and development, whereas those associated with growth and budburst were also located in genes putatively involved in response to environment and, to a lesser extent, reproduction. At one of the two independent sites, the model we developed produced highly significant correlations between predicted values and observed height data (YA, height 2020: r = 0.376, p < 0.001). Predicted values estimated with our budburst model were weakly but positively correlated with duration of budburst at one of the sites (GS, 2015: r = 0.204, p = 0.034; 2018: r = 0.205, p = 0.034-0.037) and negatively associated with budburst timing at the other (YA: r = -0.202, p = 0.046). Genomic prediction resulted in the selection of sets of trees whose mean height was taller than the average for each site. Our results provide tentative support for the capability of prediction models to forecast trait values in trees, while highlighting the need for caution in applying them to trees grown in different environments.

Perspectives and recent progress of genome-wide association studies (GWAS) in fruits

BACKGROUND

Earlier next-generation sequencing technologies are being vastly used to explore, administer, and investigate the gene space with accurate profiling of nucleotide variations in the germplasm. OVERVIEW AND PROGRESS: Recently, novel advancements in high-throughput sequencing technologies allow a genotyping-by-sequencing approach that has opened up new horizons for extensive genotyping exploiting single-nucleotide-polymorphisms (SNPs). This method acts as a bridge to support and minimize a genotype to phenotype gap allowing genetic selection at the genome-wide level, named genomic selection that could facilitate the selection of traits also in the pomology sector. In addition to this, genome-wide genotyping is a prerequisite for genome-wide association studies that have been used successfully to discover the genes, which control polygenic traits including the genetic loci, associated with the trait of interest in fruit crops.

AIMS AND PROSPECTS

This review article emphasizes the role of genome-wide approaches to unlock and explore the genetic potential along with the detection of SNPs affecting the phenotype of fruit crops and highlights the prospects of genome-wide association studies in fruits.

Inheritance of Yield Components and Morphological Traits in Avocado cv. Hass From "Criollo" "Elite Trees" via Half-Sib Seedling Rootstocks

Grafting induces precocity and maintains clonal integrity in fruit tree crops. However, the complex rootstock × scion interaction often precludes understanding how the tree phenotype is shaped, limiting the potential to select optimum rootstocks. Therefore, it is necessary to assess (1) how seedling progenies inherit trait variation from elite 'plus trees', and (2) whether such family superiority may be transferred after grafting to the clonal scion. To bridge this gap, we quantified additive genetic parameters (i.e., narrow sense heritability-h ², and genetic-estimated breeding values-GEBVs) across landraces, "criollo", "plus trees" of the super-food fruit tree crop avocado (Persea americana Mill.), and their open-pollinated (OP) half-sib seedling families. Specifically, we used a genomic best linear unbiased prediction (G-BLUP) model to merge phenotypic characterization of 17 morpho-agronomic traits with genetic screening of 13 highly polymorphic SSR markers in a diverse panel of 104 avocado "criollo" "plus trees." Estimated additive genetic parameters were validated at a 5-year-old common garden trial (i.e., provenance test), in which 22 OP half-sib seedlings from 82 elite "plus trees" served as rootstocks for the cv. Hass clone. Heritability (h ²) scores in the "criollo" "plus trees" ranged from 0.28 to 0.51. The highest h ² values were observed for ribbed petiole and adaxial veins with 0.47 (CI 95%0.2-0.8) and 0.51 (CI 0.2-0.8), respectively. The h ² scores for the agronomic traits ranged from 0.34 (CI 0.2-0.6) to 0.39 (CI 0.2-0.6) for seed weight, fruit weight, and total volume, respectively. When inspecting yield variation across 5-year-old grafted avocado cv. Hass trees with elite OP half-sib seedling rootstocks, the traits total number of fruits and fruits' weight, respectively, exhibited h ² scores of 0.36 (± 0.23) and 0.11 (± 0.09). Our results indicate that elite "criollo" "plus trees" may serve as promissory donors of seedling rootstocks for avocado cv. Hass orchards due to the inheritance of their outstanding trait values. This reinforces the feasibility to leverage natural variation from "plus trees" via OP half-sib seedling rootstock families. By jointly estimating half-sib family effects and rootstock-mediated heritability, this study promises boosting seedling rootstock breeding programs, while better discerning the consequences of grafting in fruit tree crops.

Combining genetic resources and elite material populations to improve the accuracy of genomic prediction in apple

Genomic selection is an attractive strategy for apple breeding that could reduce the length of breeding cycles. A possible limitation to the practical implementation of this approach lies in the creation of a training set large and diverse enough to ensure accurate predictions. In this study, we investigated the potential of combining two available populations, i.e., genetic resources and elite material, in order to obtain a large training set with a high genetic diversity. We compared the predictive ability of genomic predictions within-population, across-population or when combining both populations, and tested a model accounting for population-specific marker effects in this last case. The obtained predictive abilities were moderate to high according to the studied trait and small increases in predictive ability could be obtained for some traits when the two populations were combined into a unique training set. We also investigated the potential of such a training set to predict hybrids resulting from crosses between the two populations, with a focus on the method to design the training set and the best proportion of each population to optimize predictions. The measured predictive abilities were very similar for all the proportions, except for the extreme cases where only one of the two populations was used in the training set, in which case predictive abilities could be lower than when using both populations. Using an optimization algorithm to choose the genotypes in the training set also led to higher predictive abilities than when the genotypes were chosen at random. Our results provide guidelines to initiate breeding programs that use genomic selection when the implementation of the training set is a limitation.

Harnessing multivariate, penalized regression methods for genomic prediction and QTL detection of drought-related traits in grapevine

Viticulture has to cope with climate change and to decrease pesticide inputs, while maintaining yield and wine quality. Breeding is a key lever to meet this challenge, and genomic prediction a promising tool to accelerate breeding programs. Multivariate methods are potentially more accurate than univariate ones. Moreover, some prediction methods also provide marker selection, thus allowing quantitative trait loci (QTLs) detection and the identification of positional candidate genes. To study both genomic prediction and QTL detection for drought-related traits in grapevine, we applied several methods, interval mapping (IM) as well as univariate and multivariate penalized regression, in a bi-parental progeny. With a dense genetic map, we simulated two traits under four QTL configurations. The penalized regression method Elastic Net (EN) for genomic prediction, and controlling the marginal False Discovery Rate on EN selected markers to prioritize the QTLs. Indeed, penalized methods were more powerful than IM for QTL detection across various genetic architectures. Multivariate prediction did not perform better than its univariate counterpart, despite strong genetic correlation between traits. Using 14 traits measured in semi-controlled conditions under different watering conditions, penalized regression methods proved very efficient for intra-population prediction whatever the genetic architecture of the trait, with predictive abilities reaching 0.68. Compared to a previous study on the same traits, these methods applied on a denser map found new QTLs controlling traits linked to drought tolerance and provided relevant candidate genes. Overall, these findings provide a strong evidence base for implementing genomic prediction in grapevine breeding.

Genomic prediction in family bulks using different traits and cross-validations in pine

Genomic prediction integrates statistical, genomic, and computational tools to improve the estimation of breeding values and increase genetic gain. Due to the broad diversity in mating systems, breeding schemes, propagation methods, and unit of selection, no universal genomic prediction approach can be applied in all crops. In a genome-wide family prediction (GWFP) approach, the family is the basic unit of selection. We tested GWFP in two loblolly pine (Pinus taeda L.) datasets: a breeding population composed of 63 full-sib families (5–20 individuals per family), and a simulated population with the same pedigree structure. In both populations, phenotypic and genomic data was pooled at the family level in silico. Marker effects were estimated to compute genomic estimated breeding values (GEBV) at the individual and family (GWFP) levels. Less than six individuals per family produced inaccurate estimates of family phenotypic performance and allele frequency. Tested across different scenarios, GWFP predictive ability was higher than those for GEBV in both populations. Validation sets composed of families with similar phenotypic mean and variance as the training population yielded predictions consistently higher and more accurate than other validation sets. Results revealed potential for applying GWFP in breeding programs whose selection unit are family, and for systems where family can serve as training sets. The GWFP approach is well suited for crops that are routinely genotyped and phenotyped at the plot-level, but it can be extended to other breeding programs. Higher predictive ability obtained with GWFP would motivate the application of genomic prediction in these situations.

Why and How to Switch to Genomic Selection: Lessons From Plant and Animal Breeding Experience

The present study is a transversal analysis of the interest in genomic selection for plant and animal species. It focuses on the arguments that may convince breeders to switch to genomic selection. The arguments are classified into three different “bricks.” The first brick considers the addition of genotyping to improve the accuracy of the prediction of breeding values. The second consists of saving costs and/or shortening the breeding cycle by replacing all or a portion of the phenotyping effort with genotyping. The third concerns population management to improve the choice of parents to either optimize crossbreeding or maintain genetic diversity. We analyse the relevance of these different bricks for a wide range of animal and plant species and sought to explain the differences between species according to their biological specificities and the organization of breeding programs.

Genetic variation in the promoter of an R2R3-MYB transcription factor determines fruit malate content in apple (Malus domestica Borkh.)

Deciphering the mechanism of malate accumulation in apple (Malus domestica Borkh.) fruits can help to improve their flavor quality and enhance their benefits for human health. Here, we analyzed malate content as a quantitative trait that is determined mainly by genetic effects. In a previous study, we identified an R2R3-MYB transcription factor named MdMYB44 that was a candidate gene in qtl08.1 (quantitative trait locus mapped to chromosome 8) of fruit malate content. In the present study, we established that MdMYB44 negatively regulates fruit malate accumulation by repressing the promoter activity of the malate-associated genes Ma1 (Al-Activated Malate Transporter 9), Ma10 (P-type ATPase 10), MdVHA-A3 (V-type ATPase A3), and MdVHA-D2 (V-type ATPase D2). Two single-nucleotide polymorphisms (SNPs) in the MdMYB44 promoter, SNP A/G and SNP T/-, were experimentally shown to associate with fruit malate content. The TATA-box in the MdMYB44 promoter in the presence of SNP A enhances the basal activity of the MdMYB44 promoter. The binding of a basic-helix-loop-helix transcription factor MdbHLH49 to the MdMYB44 promoter was enhanced by the presence of SNP T, leading to increased MdMYB44 transcript levels and reduced malate accumulation. Furthermore, MdbHLH49 interacts with MdMYB44 and enhances MdMYB44 activity. The two SNPs could be used in combination to select for sour or non-sour apples, providing a valuable tool for the selection of fruit acidity by the apple breeding industry. This research is important for understanding the complex molecular mechanisms of fruit malate accumulation and accelerating the development of germplasm innovation in apple species and cultivars.

Genomic selection and genetic gain for nut yield in an Australian macadamia breeding population

Background

Improving yield prediction and selection efficiency is critical for tree breeding. This is vital for macadamia trees with the time from crossing to production of new cultivars being almost a quarter of a century. Genomic selection (GS) is a useful tool in plant breeding, particularly with perennial trees, contributing to an increased rate of genetic gain and reducing the length of the breeding cycle. We investigated the potential of using GS methods to increase genetic gain and accelerate selection efficiency in the Australian macadamia breeding program with comparison to traditional breeding methods. This study evaluated the prediction accuracy of GS in a macadamia breeding population of 295 full-sib progeny from 32 families (29 parents, reciprocals combined), along with a subset of parents. Historical yield data for tree ages 5 to 8years were used in the study, along with a set of 4113 SNP markers. The traits of focus were average nut yield from tree ages 5 to 8years and yield stability, measured as the standard deviation of yield over these 4 years. GBLUP GS models were used to obtain genomic estimated breeding values for each genotype, with a five-fold cross-validation method and two techniques: prediction across related populations and prediction across unrelated populations.

Results

Narrow-sense heritability of yield and yield stability was low (h²=0.30 and 0.04, respectively). Prediction accuracy for yield was 0.57 for predictions across related populations and 0.14 when predicted across unrelated populations. Accuracy of prediction of yield stability was high (r=0.79) for predictions across related populations. Predicted genetic gain of yield using GS in related populations was 474g/year, more than double that of traditional breeding methods (226g/year), due to the halving of generation length from 8 to 4years.

Conclusions

The results of this study indicate that the incorporation of GS for yield into the Australian macadamia breeding program may accelerate genetic gain due to reduction in generation length, though the cost of genotyping appears to be a constraint at present.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07694-z.

Role of MdERF3 and MdERF118 natural variations in apple flesh firmness/crispness retainability and development of QTL-based genomics-assisted prediction

Retention of flesh texture attributes during cold storage is critical for the long-term maintenance of fruit quality. The genetic variations determining flesh firmness and crispness retainability are not well understood. The objectives of this study are to identify gene markers based on quantitative trait loci (QTLs) and to develop genomics-assisted prediction (GAP) models for apple flesh firmness and crispness retainability. Phenotype data of 2664 hybrids derived from three Malus domestica cultivars and a M. asiatica cultivar were collected in 2016 and 2017. The phenotype segregated considerably with high broad-sense heritability of 83.85% and 83.64% for flesh firmness and crispness retainability, respectively. Fifty-six candidate genes were predicted from the 62 QTLs identified using bulked segregant analysis and RNA-seq. The genotype effects of the markers designed on each candidate gene were estimated. The genomics-predicted values were obtained using pyramiding marker genotype effects and overall mean phenotype values. Fivefold cross-validation revealed that the prediction accuracy was 0.5541 and 0.6018 for retainability of flesh firmness and crispness, respectively. An 8-bp deletion in the MdERF3 promoter disrupted MdDOF5.3 binding, reduced MdERF3 expression, relieved the inhibition on MdPGLR3, MdPME2, and MdACO4 expression, and ultimately decreased flesh firmness and crispness retainability. A 3-bp deletion in the MdERF118 promoter decreased its expression by disrupting the binding of MdRAVL1, which increased MdPGLR3 and MdACO4 expression and reduced flesh firmness and crispness retainability. These results provide insights regarding the genetic variation network regulating flesh firmness and crispness retainability, and the GAP models can assist in apple breeding.

Lipid nutritional quality of marine and freshwater bivalves and their aquaculture potential

Omega-3 Long-chain polyunsaturated fatty acids (n-3 LC-PUFA) are beneficial to human health. Since the industrial revolution, with the tremendous increase of human population, the supply of natural n-3 LC-PUFA is far lower than the nutritional need of n-3 LC-PUFA. Therefore, a new alternative source of natural n-3 LC-PUFA is urgently needed to reduce the supply and demand gap of n-3 LC-PUFA. Mollusks, mainly bivalves, are rich in n-3 LC-PUFA, but the information of bivalves' lipid profile is not well organized. Therefore, this study aims to analyze the published fatty acid profiles of bivalves and reveal the potential of bivalve aquaculture in meeting the nutritional needs of human for n-3 LC-PUFA. There are growing evidence show that the nutritional quality of bivalve lipid is not only species-specific, but also geographical specific. To date, bivalve aquaculture has not been evenly practiced across the globe. It can be seen that aquaculture is predominant in Asia, especially China. Unlike fish aquaculture, bivalve aquaculture does not rely on fishmeal and fish oil inputs, so it has better room for expansion. In order to unleash the full potential of bivalve aquaculture, there are some challenges need to be addressed, including recurrent mass mortalities of farmed bivalves, food safety and food security issues. The information of this article is very useful to provide an overview of lipid nutritional quality of bivalves, and reveal the potential of bivalve aquaculture in meeting the growing demand of human for n-3 LC-PUFA.

Integration of Infinium and Axiom SNP array data in the outcrossing species Malus × domestica and causes for seemingly incompatible calls

Background

Single nucleotide polymorphism (SNP) array technology has been increasingly used to generate large quantities of SNP data for use in genetic studies. As new arrays are developed to take advantage of new technology and of improved probe design using new genome sequence and panel data, a need to integrate data from different arrays and array platforms has arisen. This study was undertaken in view of our need for an integrated high-quality dataset of Illumina Infinium® 20 K and Affymetrix Axiom® 480 K SNP array data in apple (Malus × domestica). In this study, we qualify and quantify the compatibility of SNP calling, defined as SNP calls that are both accurate and concordant, across both arrays by two approaches. First, the concordance of SNP calls was evaluated using a set of 417 duplicate individuals genotyped on both arrays starting from a set of 10,295 robust SNPs on the Infinium array. Next, the accuracy of the SNP calls was evaluated on additional germplasm (n = 3141) from both arrays using Mendelian inconsistent and consistent errors across thousands of pedigree links. While performing this work, we took the opportunity to evaluate reasons for probe failure and observed discordant SNP calls.

Results

Concordance among the duplicate individuals was on average of 97.1% across 10,295 SNPs. Of these SNPs, 35% had discordant call(s) that were further curated, leading to a final set of 8412 (81.7%) SNPs that were deemed compatible. Compatibility was highly influenced by the presence of alternate probe binding locations and secondary polymorphisms. The impact of the latter was highly influenced by their number and proximity to the 3′ end of the probe.

Conclusions

The Infinium and Axiom SNP array data were mostly compatible. However, data integration required intense data filtering and curation. This work resulted in a workflow and information that may be of use in other data integration efforts. Such an in-depth analysis of array concordance and accuracy as ours has not been previously described in the literature and will be useful in future work on SNP array data integration and interpretation, and in probe/platform development.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07565-7.

Tracing founder haplotypes of Japanese apple varieties: application in genomic prediction and genome-wide association study

Haplotypes provide useful information for genomics-based approaches, genomic prediction, and genome-wide association study. As a small number of superior founders have contributed largely to the breeding history of fruit trees, the information of founder haplotypes may be relevant for performing the genomics-based approaches in these plants. In this study, we proposed a method to estimate 14 haplotypes from 7 founders and automatically trace the haplotypes forward to apple parental (185 varieties) and breeding (659 F1 individuals from 16 full-sib families) populations based on 11,786 single-nucleotide polymorphisms, by combining multiple algorithms. Overall, 92% of the single-nucleotide polymorphisms information in the parental and breeding populations was characterized by the 14 founder haplotypes. The use of founder haplotype information improved the accuracy of genomic prediction in 7 traits and the resolution of genome-wide association study in 13 out of 27 fruit quality traits analyzed in this study. We also visualized the significant propagation of the founder haplotype with the largest genetic effect in genome-wide association study over the pedigree tree of the parental population. These results suggest that the information of founder haplotypes can be useful for not only genetic improvement of fruit quality traits in apples but also for understanding the selection history of founder haplotypes in the breeding program of Japanese apple varieties.

Integrated genomic selection for rapid improvement of crops

An increase in the rate of crop improvement is essential for achieving sustained food production and other needs of ever-increasing population. Genomic selection (GS) is a potential breeding tool that has been successfully employed in animal breeding and is being incorporated into plant breeding. GS promises accelerated breeding cycles through a rapid selection of superior genotypes. Numerous empirical and simulation studies on GS and realized impacts on improvement in the crop yields are recently being reported. For a holistic understanding of the technology, we briefly discuss the concept of genetic gain, GS methodology, its current status, advantages of GS over other breeding methods, prediction models, and the factors controlling prediction accuracy in GS. Also, integration of speed breeding and other novel technologies viz. high throughput genotyping and phenotyping technologies for enhancing the efficiency and pace of GS, followed by its prospective applications in varietal development programs is reviewed.

Inheritance of Rootstock Effects in Avocado (Persea americana Mill.) cv. Hass

Grafting is typically utilized to merge adapted seedling rootstocks with highly productive clonal scions. This process implies the interaction of multiple genomes to produce a unique tree phenotype. However, the interconnection of both genotypes obscures individual contributions to phenotypic variation (rootstock-mediated heritability), hampering tree breeding. Therefore, our goal was to quantify the inheritance of seedling rootstock effects on scion traits using avocado (Persea americana Mill.) cv. Hass as a model fruit tree. We characterized 240 diverse rootstocks from 8 avocado cv. Hass orchards with similar management in three regions of the province of Antioquia, northwest Andes of Colombia, using 13 microsatellite markers simple sequence repeats (SSRs). Parallel to this, we recorded 20 phenotypic traits (including morphological, biomass/reproductive, and fruit yield and quality traits) in the scions for 3 years (2015-2017). Relatedness among rootstocks was inferred through the genetic markers and inputted in a "genetic prediction" model to calculate narrow-sense heritabilities (h 2) on scion traits. We used three different randomization tests to highlight traits with consistently significant heritability estimates. This strategy allowed us to capture five traits with significant heritability values that ranged from 0.33 to 0.45 and model fits (r) that oscillated between 0.58 and 0.73 across orchards. The results showed significance in the rootstock effects for four complex harvest and quality traits (i.e., total number of fruits, number of fruits with exportation quality, and number of fruits discarded because of low weight or thrips damage), whereas the only morphological trait that had a significant heritability value was overall trunk height (an emergent property of the rootstock-scion interaction). These findings suggest the inheritance of rootstock effects, beyond root phenotype, on a surprisingly wide spectrum of scion traits in "Hass" avocado. They also reinforce the utility of polymorphic SSRs for relatedness reconstruction and genetic prediction of complex traits. This research is, up to date, the most cohesive evidence of narrow-sense inheritance of rootstock effects in a tropical fruit tree crop. Ultimately, our work highlights the importance of considering the rootstock-scion interaction to broaden the genetic basis of fruit tree breeding programs while enhancing our understanding of the consequences of grafting.

Adoption and Optimization of Genomic Selection To Sustain Breeding for Apricot Fruit Quality

Genomic selection (GS) is a breeding approach which exploits genome-wide information and whose unprecedented success has shaped several animal and plant breeding schemes through delivering their genetic progress. This is the first study assessing the potential of GS in apricot (Prunus armeniaca) to enhance postharvest fruit quality attributes. Genomic predictions were based on a F1 pseudo-testcross population, comprising 153 individuals with contrasting fruit quality traits. They were phenotyped for physical and biochemical fruit metrics in contrasting climatic conditions over two years. Prediction accuracy (PA) varied from 0.31 for glucose content with the Bayesian LASSO (BL) to 0.78 for ethylene production with RR-BLUP, which yielded the most accurate predictions in comparison to Bayesian models and only 10% out of 61,030 SNPs were sufficient to reach accurate predictions. Useful insights were provided on the genetic architecture of apricot fruit quality whose integration in prediction models improved their performance, notably for traits governed by major QTL. Furthermore, multivariate modeling yielded promising outcomes in terms of PA within training partitions partially phenotyped for target traits. This provides a useful framework for the implementation of indirect selection based on easy-to-measure traits. Thus, we highlighted the main levers to take into account for the implementation of GS for fruit quality in apricot, but also to improve the genetic gain in perennial species.

The apple REFPOP-a reference population for genomics-assisted breeding in apple

Breeding of apple is a long-term and costly process due to the time and space requirements for screening selection candidates. Genomics-assisted breeding utilizes genomic and phenotypic information to increase the selection efficiency in breeding programs, and measurements of phenotypes in different environments can facilitate the application of the approach under various climatic conditions. Here we present an apple reference population: the apple REFPOP, a large collection formed of 534 genotypes planted in six European countries, as a unique tool to accelerate apple breeding. The population consisted of 269 accessions and 265 progeny from 27 parental combinations, representing the diversity in cultivated apple and current European breeding material, respectively. A high-density genome-wide dataset of 303,239 SNPs was produced as a combined output of two SNP arrays of different densities using marker imputation with an imputation accuracy of 0.95. Based on the genotypic data, linkage disequilibrium was low and population structure was weak. Two well-studied phenological traits of horticultural importance were measured. We found marker-trait associations in several previously identified genomic regions and maximum predictive abilities of 0.57 and 0.75 for floral emergence and harvest date, respectively. With decreasing SNP density, the detection of significant marker-trait associations varied depending on trait architecture. Regardless of the trait, 10,000 SNPs sufficed to maximize genomic prediction ability. We confirm the suitability of the apple REFPOP design for genomics-assisted breeding, especially for breeding programs using related germplasm, and emphasize the advantages of a coordinated and multinational effort for customizing apple breeding methods in the genomics era.

Quantitative trait loci-based genomics-assisted prediction for the degree of apple fruit cover color

Apple fruit cover color is an important appearance trait determining fruit quality, high degree of fruit cover color or completely red fruit skin is also the ultimate breeding goal. MdMYB1 has repeatedly been reported as a major gene controlling apple fruit cover color. There are also multiple minor-effect genes affecting degree of fruit cover color (DFC). This study was to identify genome-wide quantitative trait loci (QTLs) and to develop genomics-assisted prediction for apple DFC. The DFC phenotype data of 9,422 hybrids from five full-sib families of Malus asiatica 'Zisai Pearl', M. domestica 'Red Fuji', 'Golden Delicious', and 'Jonathan' were collected in 2014-2017. The phenotype varied considerably among hybrids with the same MdMYB1 genotype. Ten QTLs for DFC were identified using MapQTL and bulked segregant analysis via sequencing. From these QTLs, ten candidate genes were predicted, including MdMYB1 from a year-stable QTL on chromosome 9 of 'Zisai Pearl' and 'Red Fuji'. Then, kompetitive allele-specific polymerase chain reaction (KASP) markers were designed on these candidate genes and 821 randomly selected hybrids were genotyped. The genotype effects of the markers were estimated. MdMYB1-1 (represented by marker H162) exhibited a partial dominant allelic effect on MdMYB1-2 and showed non-allelic epistasis on markers H1245 and G6. Finally, a non-additive QTL-based genomics assisted prediction model was established for DFC. The Pearson's correlation coefficient between the genomic predicted value and the observed phenotype value was 0.5690. These results can be beneficial for apple genomics-assisted breeding and may provide insights for understanding the mechanism of fruit coloration.

Genomic predictions and genome-wide association studies based on RAD-seq of quality-related metabolites for the genomics-assisted breeding of tea plants

Effectively using genomic information greatly accelerates conventional breeding and applying it to long-lived crops promotes the conversion to genomic breeding. Because tea plants are bred using conventional methods, we evaluated the potential of genomic predictions (GPs) and genome-wide association studies (GWASs) for the genetic breeding of tea quality-related metabolites using genome-wide single nucleotide polymorphisms (SNPs) detected from restriction site-associated DNA sequencing of 150 tea accessions. The present GP, based on genome-wide SNPs, and six models produced moderate prediction accuracy values (r) for the levels of most catechins, represented by ( -)-epigallocatechin gallate (r = 0.32-0.41) and caffeine (r = 0.44-0.51), but low r values for free amino acids and chlorophylls. Integrated analysis of GWAS and GP detected potential candidate genes for each metabolite using 80-160 top-ranked SNPs that resulted in the maximum cumulative prediction value. Applying GPs and GWASs to tea accession traits will contribute to genomics-assisted tea breeding.

Identification of two QTLs associated with high fruit acidity in apple using pooled genome sequencing analysis

Acidity is a critical component determining apple fruit quality. Previous studies reported two major acidity quantitative trait loci (QTLs) on linkage groups (LGs) 16 (Ma) and 8 (Ma3), respectively, and their homozygous genotypes mama and ma3ma3 usually confer low titratable acidity (TA) (<3.0 mg ml-1) to apple fruit. However, apples of genotypes Ma- (MaMa and Mama) or Ma3- (Ma3Ma3 and Ma3ma3) frequently show an acidity range spanning both regular (TA 3.0-10.0 mg ml-1) and high (TA > 10 mg ml-1) acidity levels. To date, the genetic control for high-acidity apples remains essentially unknown. In order to map QTLs associated with high acidity, two genomic DNA pools, one for high acidity and the other for regular acidity, were created in an interspecific F1 population Royal Gala (Malus domestica) × PI 613988 (M. sieversii) of 191 fruit-bearing progenies. By Illumina paired-end sequencing of the high and regular acidity pools, 1,261,640 single-nucleotide variants (SNVs) commonly present in both pools were detected. Using allele frequency directional difference and density (AFDDD) mapping approach, one region on chromosome 4 and another on chromosome 6 were identified to be putatively associated with high acidity, and were named Ma6 and Ma4, respectively. Trait association analysis of DNA markers independently developed from the Ma6 and Ma4 regions confirmed the mapping of Ma6 and Ma4. In the background of MaMa, 20.6% of acidity variation could be explained by Ma6, 28.5% by Ma4, and 50.7% by the combination of both. The effects of Ma6 and Ma4 in the background of Mama were also significant, but lower. These findings provide important genetic insight into high acidity in apple.

Genomic prediction of fruit texture and training population optimization towards the application of genomic selection in apple

Texture is a complex trait and a major component of fruit quality in apple. While the major effect of MdPG1, a gene controlling firmness, has already been exploited in elite cultivars, the genetic basis of crispness remains poorly understood. To further improve fruit texture, harnessing loci with minor effects via genomic selection is therefore necessary. In this study, we measured acoustic and mechanical features in 537 genotypes to dissect the firmness and crispness components of fruit texture. Predictions of across-year phenotypic values for these components were calculated using a model calibrated with 8,294 SNP markers. The best prediction accuracies following cross-validations within the training set of 259 genotypes were obtained for the acoustic linear distance (0.64). Predictions for biparental families using the entire training set varied from low to high accuracy, depending on the family considered. While adding siblings or half-siblings into the training set did not clearly improve predictions, we performed an optimization of the training set size and composition for each validation set. This allowed us to increase prediction accuracies by 0.17 on average, with a maximal accuracy of 0.81 when predicting firmness in the 'Gala' × 'Pink Lady' family. Our results therefore identified key genetic parameters to consider when deploying genomic selection for texture in apple. In particular, we advise to rely on a large training population, with high phenotypic variability from which a 'tailored training population' can be extracted using a priori information on genetic relatedness, in order to predict a specific target population.

Genomics-assisted prediction of salt and alkali tolerances and functional marker development in apple rootstocks

Background

Saline, alkaline, and saline-alkaline stress severely affect plant growth and development. The tolerance of plants to these stressors has long been important breeding objectives, especially for woody perennials like apple. The aims of this study were to identify quantitative trait loci (QTLs) and to develop genomics-assisted prediction models for salt, alkali, and salt-alkali tolerance in apple rootstock.

Results

A total of 3258 hybrids derived from the apple rootstock cultivars ‘Baleng Crab’ (Malus robusta Rehd., tolerant) × ‘M9’ (M. pumila Mill., sensitive) were used to identify 17, 13, and two QTLs for injury indices of salt, alkali, and salt–alkali stress via bulked segregant analysis. The genotype effects of single nucleotide polymorphism (SNP) markers designed on candidate genes in each QTL interval were estimated. The genomic predicted value of an individual hybrid was calculated by adding the sum of all marker genotype effects to the mean phenotype value of the population. The prediction accuracy was 0.6569, 0.6695, and 0.5834 for injury indices of salt, alkali, and salt–alkali stress, respectively. SNP182G on MdRGLG3, which changes a leucine to an arginine at the vWFA-domain, conferred tolerance to salt, alkali, and salt-alkali stress. SNP761A on MdKCAB, affecting the Kv_beta domain that cooperated with the linked allelic variation SNP11, contributed to salt, alkali, and salt–alkali tolerance in apple rootstock.

Conclusions

The genomics-assisted prediction models can potentially be used in breeding saline, alkaline, and saline-alkaline tolerant apple rootstocks. The QTLs and the functional markers may provide insight for future studies into the genetic variation of plant abiotic stress tolerance.

Linkage disequilibrium vs. pedigree: Genomic selection prediction accuracy in conifer species

BACKGROUND

The presupposition of genomic selection (GS) is that predictive accuracies should be based on population-wide linkage disequilibrium (LD). However, in species with large, highly complex genomes the limitation of marker density may preclude the ability to resolve LD accurately enough for GS. Here we investigate such an effect in two conifer species with ~ 20 Gbp genomes, Douglas-fir (Pseudotsuga menziesii Mirb. (Franco)) and Interior spruce (Picea glauca (Moench) Voss x Picea engelmannii Parry ex Engelm.). Random sampling of markers was performed to obtain SNP sets with totals in the range of 200-50,000, this was replicated 10 times. Ridge Regression Best Linear Unbiased Predictor (RR-BLUP) was deployed as the GS method to test these SNP sets, and 10-fold cross-validation was performed on 1,321 Douglas-fir trees, representing 37 full-sib F1 families and on 1,126 Interior spruce trees, representing 25 open-pollinated (half-sib) families. Both trials are located on 3 sites in British Columbia, Canada.

RESULTS

As marker number increased, so did GS predictive accuracy for both conifer species. However, a plateau in the gain of accuracy became apparent around 10,000-15,000 markers for both Douglas-fir and Interior spruce. Despite random marker selection, little variation in predictive accuracy was observed across replications. On average, Douglas-fir prediction accuracies were higher than those of Interior spruce, reflecting the difference between full- and half-sib families for Douglas-fir and Interior spruce populations, respectively, as well as their respective effective population size.

CONCLUSIONS

Although possibly advantageous within an advanced breeding population, reducing marker density cannot be recommended for carrying out GS in conifers. Significant LD between markers and putative causal variants was not detected using 50,000 SNPS, and GS was enabled only through the tracking of relatedness in the populations studied. Dramatically increasing marker density would enable said markers to better track LD with causal variants in these large, genetically diverse genomes; as well as providing a model that could be used across populations, breeding programs, and traits.

Turbocharging introgression breeding of perennial fruit crops: a case study on apple

The allelic diversity of primitive germplasm of fruit crops provides a useful resource for introgressing novel genes to meet consumer preferences and environmental challenges. Pre-breeding facilitates the identification of novel genetic variation in the primitive germplasm and expedite its utilisation in cultivar breeding programmes. Several generations of pre-breeding could be required to minimise linkage drag from the donor parent and to maximise the genomic content of the recipient parent. In this study we investigated the potential of genomic selection (GS) as a tool for rapid background selection of parents for the successive generation. A diverse set of 274 accessions was genotyped using random-tag genotyping-by-sequencing, and phenotyped for eight fruit quality traits. The relationship between 'own phenotypes' of 274 accessions and their general combining ability (GCA) was also examined. Trait heritability influenced the strength of correspondence between own phenotype and the GCA. The average (across eight traits) accuracy of predicting own phenotype was 0.70, and the correlations between genomic-predicted own phenotype and GCA were similar to the observed correlations. Our results suggest that genome-assisted parental selection (GAPS) is a credible alternative to phenotypic parental selection, so could help reduce the generation interval to allow faster accumulation of favourable alleles from donor and recipient parents.

Genetic and genomic resources for Rubus breeding: a roadmap for the future

Rubus fruits are high-value crops that are sought after by consumers for their flavor, visual appeal, and health benefits. To meet this demand, production of red and black raspberries (R. idaeus L. and R. occidentalis L.), blackberries (R. subgenus Rubus), and hybrids, such as Boysenberry and marionberry, is growing worldwide. Rubus breeding programmes are continually striving to improve flavor, texture, machine harvestability, and yield, provide pest and disease resistance, improve storage and processing properties, and optimize fruits and plants for different production and harvest systems. Breeders face numerous challenges, such as polyploidy, the lack of genetic diversity in many of the elite cultivars, and until recently, the relative shortage of genetic and genomic resources available for Rubus. This review will highlight the development of continually improving genetic maps, the identification of Quantitative Trait Loci (QTL)s controlling key traits, draft genomes for red and black raspberry, and efforts to improve gene models. The development of genetic maps and markers, the molecular characterization of wild species and germplasm, and high-throughput genotyping platforms will expedite breeding of improved cultivars. Fully sequenced genomes and accurate gene models facilitate identification of genes underlying traits of interest and enable gene editing technologies such as CRISPR/Cas9.

Genome-wide association studies in apple reveal loci of large effect controlling apple polyphenols

Apples are a nutritious food source with significant amounts of polyphenols that contribute to human health and wellbeing, primarily as dietary antioxidants. Although numerous pre- and post-harvest factors can affect the composition of polyphenols in apples, genetics is presumed to play a major role because polyphenol concentration varies dramatically among apple cultivars. Here we investigated the genetic architecture of apple polyphenols by combining high performance liquid chromatography (HPLC) data with ~100,000 single nucleotide polymorphisms (SNPs) from two diverse apple populations. We found that polyphenols can vary in concentration by up to two orders of magnitude across cultivars, and that this dramatic variation was often predictable using genetic markers and frequently controlled by a small number of large effect genetic loci. Using GWAS, we identified candidate genes for the production of quercitrin, epicatechin, catechin, chlorogenic acid, 4-O-caffeoylquinic acid and procyanidins B1, B2, and C1. Our observation that a relatively simple genetic architecture underlies the dramatic variation of key polyphenols in apples suggests that breeders may be able to improve the nutritional value of apples through marker-assisted breeding or gene editing.

High-quality, genome-wide SNP genotypic data for pedigreed germplasm of the diploid outbreeding species apple, peach, and sweet cherry through a common workflow

High-quality genotypic data is a requirement for many genetic analyses. For any crop, errors in genotype calls, phasing of markers, linkage maps, pedigree records, and unnoticed variation in ploidy levels can lead to spurious marker-locus-trait associations and incorrect origin assignment of alleles to individuals. High-throughput genotyping requires automated scoring, as manual inspection of thousands of scored loci is too time-consuming. However, automated SNP scoring can result in errors that should be corrected to ensure recorded genotypic data are accurate and thereby ensure confidence in downstream genetic analyses. To enable quick identification of errors in a large genotypic data set, we have developed a comprehensive workflow. This multiple-step workflow is based on inheritance principles and on removal of markers and individuals that do not follow these principles, as demonstrated here for apple, peach, and sweet cherry. Genotypic data was obtained on pedigreed germplasm using 6-9K SNP arrays for each crop and a subset of well-performing SNPs was created using ASSIsT. Use of correct (and corrected) pedigree records readily identified violations of simple inheritance principles in the genotypic data, streamlined with FlexQTL software. Retained SNPs were grouped into haploblocks to increase the information content of single alleles and reduce computational power needed in downstream genetic analyses. Haploblock borders were defined by recombination locations detected in ancestral generations of cultivars and selections. Another round of inheritance-checking was conducted, for haploblock alleles (i.e., haplotypes). High-quality genotypic data sets were created using this workflow for pedigreed collections representing the U.S. breeding germplasm of apple, peach, and sweet cherry evaluated within the RosBREED project. These data sets contain 3855, 4005, and 1617 SNPs spread over 932, 103, and 196 haploblocks in apple, peach, and sweet cherry, respectively. The highly curated phased SNP and haplotype data sets, as well as the raw iScan data, of germplasm in the apple, peach, and sweet cherry Crop Reference Sets is available through the Genome Database for Rosaceae.

Marker-trait associations and genomic predictions of interspecific pear (Pyrus) fruit characteristics

Interspecific pear (Pyrus spp.) hybrid populations are often used to develop novel cultivars. Pear cultivar breeding is a lengthy process because of long juvenility and the subsequent time required for reliable fruit phenotyping. Molecular techniques such as genome-wide association (GWA) and genomic selection (GS) provide an opportunity to fast-forward the development of high-value cultivars. We evaluated the genetic architecture of 10 pear fruit phenotypes (including sensory traits) and the potential of GS using genotyping-by-sequencing of 550 hybrid seedlings from nine interrelated full-sib families. Results from GWA suggested a complex polygenic nature of all 10 traits as the maximum variance explained by each marker was less than 4% of the phenotypic variance. The effect-size of SNPs for each trait suggested many genes of small effect and few of moderate effect. Some genomic regions associated with pear sensory traits were similar to those reported for apple - possibly a result of high synteny between the apple and pear genomes. The average (across nine families) GS accuracy varied from 0.32 (for crispness) to 0.62 (for sweetness), with an across-trait average of 0.42. Further efforts are needed to develop larger genotype-phenotype datasets in order to predict fruit phenotypes of untested seedlings with sufficient efficiency.

An R Package for Bayesian Analysis of Multi-environment and Multi-trait Multi-environment Data for Genome-Based Prediction

Evidence that genomic selection (GS) is a technology that is revolutionizing plant breeding continues to grow. However, it is very well documented that its success strongly depends on statistical models, which are used by GS to perform predictions of candidate genotypes that were not phenotyped. Because there is no universally better model for prediction and models for each type of response variable are needed (continuous, binary, ordinal, count, etc.), an active area of research aims to develop statistical models for the prediction of univariate and multivariate traits in GS. However, most of the models developed so far are for univariate and continuous (Gaussian) traits. Therefore, to overcome the lack of multivariate statistical models for genome-based prediction by improving the original version of the BMTME, we propose an improved Bayesian multi-trait and multi-environment (BMTME) R package for analyzing breeding data with multiple traits and multiple environments. We also introduce Bayesian multi-output regressor stacking (BMORS) functions that are considerably efficient in terms of computational resources. The package allows parameter estimation and evaluates the prediction performance of multi-trait and multi-environment data in a reliable, efficient and user-friendly way. We illustrate the use of the BMTME with real toy datasets to show all the facilities that the software offers the user. However, for large datasets, the BME() and BMTME() functions of the BMTME R package are very intense in terms of computing time; on the other hand, less intensive computing is required with BMORS functions BMORS() and BMORS_Env() that are also included in the BMTME package.

Development of a highly efficient Axiom™ 70 K SNP array for Pyrus and evaluation for high-density mapping and germplasm characterization

Background

Both a source of diversity and the development of genomic tools, such as reference genomes and molecular markers, are equally important to enable faster progress in plant breeding. Pear (Pyrus spp.) lags far behind other fruit and nut crops in terms of employment of available genetic resources for new cultivar development. To address this gap, we designed a high-density, high-efficiency and robust single nucleotide polymorphism (SNP) array for pear, with the main objectives of conducting genetic diversity and genome-wide association studies.

Results

By applying a two-step design process, which consisted of the construction of a first ‘draft’ array for the screening of a small subset of samples, we were able to identify the most robust and informative SNPs to include in the Applied Biosystems™ Axiom™ Pear 70 K Genotyping Array, currently the densest SNP array for pear. Preliminary evaluation of this 70 K array in 1416 diverse pear accessions from the USDA National Clonal Germplasm Repository (NCGR) in Corvallis, OR identified 66,616 SNPs (93% of all the tiled SNPs) as high quality and polymorphic (PolyHighResolution). We further used the Axiom Pear 70 K Genotyping Array to construct high-density linkage maps in a bi-parental population, and to make a direct comparison with available genotyping-by-sequencing (GBS) data, which suggested that the SNP array is a more robust method of screening for SNPs than restriction enzyme reduced representation sequence-based genotyping.

Conclusions

The Axiom Pear 70 K Genotyping Array, with its high efficiency in a widely diverse panel of Pyrus species and cultivars, represents a valuable resource for a multitude of molecular studies in pear. The characterization of the USDA-NCGR collection with this array will provide important information for pear geneticists and breeders, as well as for the optimization of conservation strategies for Pyrus.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5712-3) contains supplementary material, which is available to authorized users.

Apple whole genome sequences: recent advances and new prospects

In 2010, a major scientific milestone was achieved for tree fruit crops: publication of the first draft whole genome sequence (WGS) for apple (Malus domestica). This WGS, v1.0, was valuable as the initial reference for sequence information, fine mapping, gene discovery, variant discovery, and tool development. A new, high quality apple WGS, GDDH13 v1.1, was released in 2017 and now serves as the reference genome for apple. Over the past decade, these apple WGSs have had an enormous impact on our understanding of apple biological functioning, trait physiology and inheritance, leading to practical applications for improving this highly valued crop. Causal gene identities for phenotypes of fundamental and practical interest can today be discovered much more rapidly. Genome-wide polymorphisms at high genetic resolution are screened efficiently over hundreds to thousands of individuals with new insights into genetic relationships and pedigrees. High-density genetic maps are constructed efficiently and quantitative trait loci for valuable traits are readily associated with positional candidate genes and/or converted into diagnostic tests for breeders. We understand the species, geographical, and genomic origins of domesticated apple more precisely, as well as its relationship to wild relatives. The WGS has turbo-charged application of these classical research steps to crop improvement and drives innovative methods to achieve more durable, environmentally sound, productive, and consumer-desirable apple production. This review includes examples of basic and practical breakthroughs and challenges in using the apple WGSs. Recommendations for "what's next" focus on necessary upgrades to the genome sequence data pool, as well as for use of the data, to reach new frontiers in genomics-based scientific understanding of apple.

Prediction of genetic value for sweet cherry fruit maturity among environments using a 6K SNP array

The timing of fruit maturity is an important trait in sweet cherry production and breeding. Phenotypic variation for phenology of fruit maturity in sweet cherry appears to be under strong genetic control, but that control might be complicated by phenotypic instability across environments. Although such genotype-by-environment interaction (G × E) is a common phenomenon in crop plants, knowledge about it is lacking for fruit maturity timing and other sweet cherry traits. In this study, 1673 genome-wide SNP markers were used to estimate genomic relationships among 597 weakly pedigree-connected individuals evaluated over two seasons at three locations in Europe and one location in the USA, thus sampling eight 'environments'. The combined dataset enabled a single meta-analysis to investigate the environmental stability of genomic predictions. Linkage disequilibrium among marker loci declined rapidly with physical distance, and ordination of the relationship matrix suggested no strong structure among germplasm. The most parsimonious G × E model allowed heterogeneous genetic variance and pairwise covariances among environments. Narrow-sense genomic heritability was very high (0.60-0.83), as was accuracy of predicted breeding values (>0.62). Average correlation of additive effects among environments was high (0.96) and breeding values were highly correlated across locations. Results indicated that genomic models can be used in cherry to accurately predict date of fruit maturity for untested individuals in new environments. Limited G × E for this trait indicated that phenotypes of individuals will be stable across similar environments. Equivalent analyses for other sweet cherry traits, for which multiple years of data are commonly available among breeders and cultivar testers, would be informative for predicting performance of elite selections and cultivars in new environments.

Validation of SNP markers for fruit quality and disease resistance loci in apple (Malus × domestica Borkh.) using the OpenArray® platform

Genome mapping has promised much to tree fruit breeding during the last 10 years. Nevertheless, one of the greatest challenges remaining to tree fruit geneticists is the translation of trait loci and whole genome sequences into diagnostic genetic markers that are efficient and cost-effective for use by breeders, who must select genetically optimal parents and subsequently select genetically superior individuals among their progeny. To take this translational step, we designed the apple International RosBREED SNP Consortium OpenArray v1.0 (IRSCOA v1.0) assay using a set of 128 apple single nucleotide polymorphisms (SNPs) linked to fruit quality and pest and disease resistance trait loci. The Thermo Fisher Scientific OpenArray® technology enables multiplexed screening of SNP markers using a real-time PCR instrument with fluorescent probe-based Taqman® assays. We validated the apple IRSCOA v1.0 multi-trait assay by screening 240 phenotyped individuals from the Plant & Food Research apple cultivar breeding programme. This set of individuals comprised commercial and heritage cultivars, elite selections, and families segregating for traits of importance to breeders. In total, 33 SNP markers of the IRSCOA v1.0 were validated for use in marker-assisted selection (MAS) for the scab resistances Rvi2/Vh2, Rvi4/Vh4, Rvi6/Vf, fire blight resistance MR5/RLP1, powdery mildew resistance Pl2, fruit firmness, skin colour, flavour intensity, and acidity. The availability of this set of validated trait-associated SNP markers, which can be used individually on multiple genotyping platforms available to various apple breeding programmes or re-designed using the flanking sequences, represents a large translational genetics step from genomics to crop improvement of apple.