Whole-Genome Analysis of Diversity and SNP-Major Gene Association in Peach Germplasm

Genome-wide presence/absence variation discovery and its application in Peach (Prunus persica)

Gene presence/absence variation (PAV) is an important contributor to the studies of genetic diversity, gene identification, and molecular marker development in plants. In the present study, 100 peach (Prunus persica) accessions were used for genome resequencing to identify PAVs. Alignmentwith a reference genome yielded a total of 2.52 Mb non-reference sequences and 923 novel genes were identified. The dispensable PAVs were enriched in resistance, perhaps reflecting their roles in plant adaptation to various environments. Furthermore, selection sweeps associated with peach domestication and improvement were identified based on PAV data. Only 4.3% and 13.4% of domestication and improvement sweeps, respectively, were identified simultaneously using single nucleotide polymorphism (SNP) data, suggesting flexible identification between the different methods. To further verify the applicability of PAV identification, a genome-wide association study was conducted using 21 agronomic traits. Some of the identified loci were consistent with those reported in previous studies, while some were mapped for the first time; the latter included petiole length, petiole gland shape, and petiole gland number. Through tissue-specific expression analysis and gene transformation experiments, a novel gene, evm.model.Contig322_A94.1, was identified and found to be involved in chilling requirements. We speculated that this novel gene might regulate the trait by participating in the ABA signaling pathway. The PAVs identified in P. persica provide valuable resources for mapping the entire gene set and identifying optional markers for molecular selection in future studies.

Multiple-statistical genome-wide association analysis and genomic prediction of fruit aroma and agronomic traits in peaches

'Chinese Cling' is an important founder in peach breeding history due to the pleasant flavor. Genome-wide association studies (GWAS) combined with genomic selection are promising tools in fruit tree breeding, as there is a considerable time lapse between crossing and release of a cultivar. In this study, 242 peaches from Shanghai germplasm were genotyped with 145 456 single-nucleotide polymorphisms (SNPs). The six agronomic traits of fruit flesh color, fruit shape, fruit hairiness, flower type, pollen sterility, and soluble solids content, along with 14 key volatile odor compounds (VOCs), were recorded for multiple-statistical GWAS. Except the reported candidate genes, six novel genes were identified as associated with these traits. Thirty-nine significant SNPs were associated with eight VOCs. The putative candidate genes were confirmed for VOCs by RNA-seq, including three genes in the biosynthesis pathway found to be associated with linalool, soluble solids content, and cis-3-hexenyl acetate. Multiple-trait genomic prediction enhanced the predictive ability for γ-decalactone to 0.7415 compared with the single-trait model value of 0.1017. One PTS1-SSR marker was designed to predict the linalool content, and the favorable genotype 187/187 was confirmed, mainly existing in the 'Shanghai Shuimi' landrace. Overall, our findings will be helpful in determining peach accessions with the ideal phenotype and show the potential of multiple-trait genomic prediction to improve accuracy for highly correlated genetic traits. The diagnostic marker will be valuable for the breeder to bridge the gap between quantitative trait loci and marker-assisted selection for developing strong-aroma cultivars.

Absence of major epigenetic and transcriptomic changes accompanying an interspecific cross between peach and almond

Hybridization has been widely used in breeding of cultivated species showing low genetic variability, such as peach (Prunus persica). The merging of two different genomes in a hybrid often triggers a so-called "genomic shock" with changes in DNA methylation and in the induction of transposable element expression and mobilization. Here, we analysed the DNA methylation and transcription levels of transposable elements and genes in leaves of Prunus persica and Prunus dulcis and in an F1 hybrid using high-throughput sequencing technologies. Contrary to the "genomic shock" expectations, we found that the overall levels of DNA methylation in the transposable elements in the hybrid are not significantly altered compared with those of the parental genomes. We also observed that the levels of transcription of the transposable elements in the hybrid are in most cases intermediate as compared with that of the parental species and we have not detected cases of higher transcription in the hybrid. We also found that the proportion of genes whose expression is altered in the hybrid compared with the parental species is low. The expression of genes potentially involved in the regulation of the activity of the transposable elements is not altered. We can conclude that the merging of the two parental genomes in this Prunus persica x Prunus dulcis hybrid does not result in a "genomic shock" with significant changes in the DNA methylation or in the transcription. The absence of major changes may facilitate using interspecific peach x almond crosses for peach improvement.

Construction of a collection of introgression lines of "Texas" almond DNA fragments in the "Earlygold" peach genetic background

Peach [Prunus persica L. Batsch] is one of the major temperate fruit tree species, the commercial materials of which have a low level of genetic variability. Almond [P. dulcis (Mill) DA Webb], a close relative of peach cultivated for its kernels, has a much higher level of diversity. The species are inter-compatible and often produce fertile hybrids, almond being a possible source of new genes for peach that could provide biotic and abiotic stress tolerance traits. In this paper we describe the development of a collection of peach-almond introgression lines (ILs) having a single fragment of almond (cv. Texas) in the peach background (cv. Earlygold). Lines with few introgressions were selected with markers from successive generations from a "Texas" × "Earlygold" F1 hybrid, initially using a set of SSRs and later with the 18 k peach SNP chip, allowing for the final extraction of 67 lines, 39 with almond heterozygous introgressions covering 99% of the genome, and 28 with homozygous introgressions covering 83% of the genome. As a proof of concept, four major genes and four quantitative characters were examined in the selected ILs giving results generally consistent with previous information on the genetics of these characters. This collection is the first of its kind produced in a woody perennial species and promises to be a valuable tool for genetic analyses, including dissection of quantitative traits, positional cloning, epistasis and as prebreeding material to introgress almond genes of interest into the peach commercial gene pool.

Genetic Diversity and Genome-Wide Association Study of Morphological and Quality Traits in Peach Using Two Spanish Peach Germplasm Collections

Peach [Prunus persica (L.) Batsch] is one of the most important stone fruits species in world production. Spanish peach production is currently the second largest in the world and the available cultivars in Spain includes a great source of genetic diversity with variability in fruit quality traits and postharvest disorders tolerance. In order to explore the genetic diversity and single nucleotide polymorphism (SNP)-trait associations in the Spanish germplasm, the new peach 18K SNP v2 array was used to genotype 287 accessions belonging to the two National Peach Germplasm Collections placed at the Agrifood Research and Technology Centre of Aragon (CITA) and at the Experimental Station of Aula Dei (EEAD)-CSIC. The high density of the new SNP array allowed the identification of 30 groups of synonymies, which had not been identified before using low-density markers. In addition, a possible large-scale molecular event in 'Starcrest', a sport of 'Springcrest', was detected showing a possible chromosome replacement of a 13.5 Mb region. Previous suggestions about Spanish diversification regions agreed with our genetic diversity and linkage disequilibrium (LD) decay results using high-density markers. A genome-wide association study (GWAS) detected 34 significant SNP-trait association with the type of leaf glands (TLG), fruit hairiness (FH), and flesh texture (FT). The impact of the significant SNPs was studied with SnpEff. Candidate genes encode several important family proteins involved in trichome formation and powdery mildew resistance (linked to TLG in peach). The genetic distance among cultivars obtained, together with SNP-trait associations found, provide new knowledge for marker-assisted selection and crossing approaches in peach breeding programmes.

Applications of Genomic Tools in Plant Breeding: Crop Biofortification

Crop breeding has mainly been focused on increasing productivity, either directly or by decreasing the losses caused by biotic and abiotic stresses (that is, incorporating resistance to diseases and enhancing tolerance to adverse conditions, respectively). Quite the opposite, little attention has been paid to improve the nutritional value of crops. It has not been until recently that crop biofortification has become an objective within breeding programs, through either conventional methods or genetic engineering. There are many steps along this long path, from the initial evaluation of germplasm for the content of nutrients and health-promoting compounds to the development of biofortified varieties, with the available and future genomic tools assisting scientists and breeders in reaching their objectives as well as speeding up the process. This review offers a compendium of the genomic technologies used to explore and create biodiversity, to associate the traits of interest to the genome, and to transfer the genomic regions responsible for the desirable characteristics into potential new varieties. Finally, a glimpse of future perspectives and challenges in this emerging area is offered by taking the present scenario and the slow progress of the regulatory framework as the starting point.

Insertion of a mMoshan transposable element in PpLMI1, is associated with the absence or globose phenotype of extrafloral nectaries in peach [Prunus persica (L.) Batsch]

Most commercial peach [Prunus persica (L.) Batsch] cultivars have leaves with extrafloral nectaries (EFNs). Breeders have selected this character over time, as they observed that the eglandular phenotype resulted in high susceptibility to peach powdery mildew, a major disease of peach trees. EFNs are controlled by a Mendelian locus (E), mapped on chromosome 7. However, the genetic factor underlying E was unknown. In order to address this point, we developed a mapping population of 833 individuals derived from the selfing of "Malo Konare", a Bulgarian peach cultivar, heterozygous for the trait. This progeny was used to investigate the E-locus region, along with additional resources including peach genomic resequencing data, and 271 individuals from various origins used for validation. High-resolution mapping delimited a 40.6 kbp interval including the E-locus and four genes. Moreover, three double-recombinants allowed identifying Prupe.7G121100, a LMI1-like homeodomain leucine zipper (HD-Zip) transcription factor, as a likely candidate for the trait. By comparing peach genomic resequencing data from individuals with contrasted phenotypes, a MITE-like transposable element of the hAT superfamily (mMoshan) was identified in the third exon of Prupe.7G121100. It was associated with the absence or globose phenotype of EFNs. The insertion of the transposon was positively correlated with enhanced expression of Prupe.7G121100. Furthermore, a PCR marker designed from the sequence-variants, allowed to properly assign the phenotypes of all the individuals studied. These findings provide valuable information on the genetic control of a trait poorly known so far although selected for a long time in peach.

Genetic Diversity and Population Structure Analysis of the USDA Olive Germplasm Using Genotyping-By-Sequencing (GBS)

Olives are one of the most important fruit and woody oil trees cultivated in many parts of the world. Olive oil is a critical component of the Mediterranean diet due to its importance in heart health. Olives are believed to have been brought to the United States from the Mediterranean countries in the 18th century. Despite the increase in demand and production areas, only a few selected olive varieties are grown in most traditional or new growing regions in the US. By understanding the genetic background, new sources of genetic diversity can be incorporated into the olive breeding programs to develop regionally adapted varieties for the US market. This study aimed to explore the genetic diversity and population structure of 90 olive accessions from the USDA repository along with six popular varieties using genotyping-by-sequencing (GBS)-generated SNP markers. After quality filtering, 54,075 SNP markers were retained for the genetic diversity analysis. The average gene diversity (GD) and polymorphic information content (PIC) values of the SNPs were 0.244 and 0.206, respectively, indicating a moderate genetic diversity for the US olive germplasm evaluated in this study. The structure analysis showed that the USDA collection was distributed across seven subpopulations; 63% of the accessions were grouped into an identifiable subpopulation. The phylogenetic and principal coordinate analysis (PCoA) showed that the subpopulations did not align with the geographical origins or climatic zones. An analysis of the molecular variance revealed that the major genetic variation sources were within populations. These findings provide critical information for future olive breeding programs to select genetically distant parents and facilitate future gene identification using genome-wide association studies (GWAS) or a marker-assisted selection (MAS) to develop varieties suited to production in the US.

Genetic and phenotypic analyses reveal major quantitative loci associated to fruit size and shape traits in a non-flat peach collection (P. persica L. Batsch)

Fruit size and shape are critical agronomical and pomological attributes and prime targets in peach breeding programs. Apart from the flat peach type, a Mendelian trait well-characterized at the genetic level, ample diversity of fruit size and shapes is present across peach germplasms. Nevertheless, knowledge of the underlying genomic loci remains limited. In this work, fruit size and shape were assessed in a collection of non-flat peach accessions and selections, under controlled fruit load conditions. The architecture of these traits was then dissected by combining association and linkage mapping, revealing a major locus on the proximal end of chromosome 6 (qSHL/Fs6.1) explaining a large proportion of phenotypic variability for longitudinal shape and also affecting fruit size. A second major locus for fruit longitudinal shape (qSHL5.1), probably also affecting fruit size, was found co-localizing at locus G, suggesting pleiotropic effects of peach/nectarine traits. An additional QTL for fruit longitudinal shape (qSHL6.2) was identified in the distal end of chromosome 6 in a cross with an ornamental double-flower peach and co-localized with the Di2 locus, controlling flower morphology. Besides assisting breeding activities, knowledge of loci controlling fruit size and shape paves the way for more in-depth studies aimed at the identification of underlying genetic variant(s).

Chromosome-level genome assemblies of five Prunus species and genome-wide association studies for key agronomic traits in peach

Prunus species include many important perennial fruit crops, such as peach, plum, apricot, and related wild species. Here, we report de novo genome assemblies for five species, including the cultivated species peach (Prunus persica), plum (Prunus salicina), and apricot (Prunus armeniaca), and the wild peach species Tibetan peach (Prunus mira) and Chinese wild peach (Prunus davidiana). The genomes ranged from 240 to 276 Mb in size, with contig N50 values of 2.27-8.30 Mb and 25,333-27,826 protein-coding gene models. As the phylogenetic tree shows, plum diverged from its common ancestor with peach, wild peach species, and apricot ~7 million years ago (MYA). We analyzed whole-genome resequencing data of 417 peach accessions, called 3,749,618 high-quality SNPs, 577,154 small indels, 31,800 deletions, duplications, and inversions, and 32,338 insertions, and performed a structural variant-based genome-wide association study (GWAS) of key agricultural traits. From our GWAS data, we identified a locus associated with a fruit shape corresponding to the OVATE transcription factor, where a large inversion event correlates with higher OVATE expression in flat-shaped accessions. Furthermore, a GWAS revealed a NAC transcription factor associated with fruit developmental timing that is linked to a tandem repeat variant and elevated NAC expression in early-ripening accessions. We also identified a locus encoding microRNA172d, where insertion of a transposable element into its promoter was found in double-flower accessions. Thus, our efforts have suggested roles for OVATE, a NAC transcription factor, and microRNA172d in fruit shape, fruit development period, and floral morphology, respectively, that can be connected to traits in other crops, thereby demonstrating the importance of parallel evolution in the diversification of several commercially important domesticated species. In general, these genomic resources will facilitate functional genomics, evolutionary research, and agronomic improvement of these five and other Prunus species. We believe that structural variant-based GWASs can also be used in other plants, animal species, and humans and be combined with deep sequencing GWASs to precisely identify candidate genes and genetic architecture components.

New high-quality peach (Prunus persica L. Batsch) genome assembly to analyze the molecular evolutionary mechanism of volatile compounds in peach fruits

Peach (Prunus persica L. Batsch) is an economically important fruit crop worldwide. Although a high-quality peach genome has previously been published, Sanger sequencing was used for its assembly, which generated short contigs. Here, we report a chromosome-level genome assembly and sequence analysis of Chinese Cling, an important founder cultivar for peach breeding programs worldwide. The assembled genome contained 247.33 Mb with a contig N50 of 4.13 Mb and a scaffold N50 of 29.68 Mb, representing 99.8% of the estimated genome. Comparisons between this genome and the recently published one (Lovell peach) uncovered 685 407 single nucleotide polymorphisms, 162 655 insertions and deletions, and 16 248 structural variants. Gene family analysis highlighted the contraction of the gene families involved in flavone, flavonol, flavonoid, and monoterpenoid biosynthesis. Subsequently, the volatile compounds of 256 peach varieties were quantitated in mature fruits in 2015 and 2016 to perform a genome-wide association analysis. A comparison with the identified domestication genomic regions allowed us to identify 25 quantitative trait loci, associated with seven volatile compounds, in the domestication region, which is consistent with the differences in volatile compounds between wild and cultivated peaches. Finally, a gene encoding terpene synthase, located within a previously reported quantitative trait loci region, was identified to be associated with linalool synthesis. Such findings highlight the importance of this new assembly for the analysis of evolutionary mechanisms and gene identification in peach species. Furthermore, this high-quality peach genome provides valuable information for future fruit improvement.

Population-scale peach genome analyses unravel selection patterns and biochemical basis underlying fruit flavor

A narrow genetic basis in modern cultivars and strong linkage disequilibrium in peach (Prunus persica) has restricted resolution power for association studies in this model fruit species, thereby limiting our understanding of economically important quality traits including fruit flavor. Here, we present a high-quality genome assembly for a Chinese landrace, Longhua Shui Mi (LHSM), a representative of the Chinese Cling peaches that have been central in global peach genetic improvement. We also map the resequencing data for 564 peach accessions to this LHSM assembly at an average depth of 26.34× per accession. Population genomic analyses reveal a fascinating history of convergent selection for sweetness yet divergent selection for acidity in eastern vs. western modern cultivars. Molecular-genetics and biochemical analyses establish that PpALMT1 (aluminum-activated malate transporter 1) contributes to their difference of malate content and that increases fructose content accounts for the increased sweetness of modern peach fruits, as regulated by PpERDL16 (early response to dehydration 6-like 16). Our study illustrates the strong utility of the genomics resources for both basic and applied efforts to understand and exploit the genetic basis of fruit quality in peach.

Peach allergen Pru p 1 content is generally low in fruit but with large variation in different varieties

Background

Pru p 1 is a major allergen in peach and nectarine, and the different content in varieties may affect the degree of allergic reactions. This study aimed to quantify Pru p 1 levels in representative peach varieties and select hypoallergenic Pru p 1 varieties.

Methods

To obtain monoclonal and polyclonal antibodies, mice and rabbits, respectively, were immunized with recombinant Pru p 1.01 and Pru p 1.02. The Pru p 1 levels in fruits from 83 representative peach varieties was quantified by sandwich enzyme-linked immunosorbent assay (sELISA). nPru p 1 was obtained through specific monoclonal antibody affinity purification and confirmed by Western blot and mass spectrometry. The variable Pru p 1 content of selected varieties was evaluated by Western blot and the expression level of encoding Pru p 1 genes by quantitative polymerase chain reaction.

Results

A sELISA method with monoclonal and polyclonal antibodies was built for quantifying Pru p 1 levels in peach. Pru p 1 was mainly concentrated in the peel (0.20-73.44 μg/g, fresh weight), being very low in the pulp (0.05-9.62 μg/g) and not detected in wild peach. For the 78 peach and nectarine varieties, Pru p 1 content varied widely from 0.12 to 6.45 μg/g in whole fruit. We verified that natural Pru p 1 is composed of 1.01 and 1.02 isoallergens, and the Pru p 1 expression level and Pru p 1 band intensity in the immunoblots were in agreement with protein quantity determined by ELISA for some tested varieties. In some cases, the reduced levels of Pru p 1 did not coincide with low Pru p 3 in the same variety in whole fruit, while some ancient wild peach and nectarines contained low levels of both allergens, and late-ripening yellow flesh varieties were usually highly allergenic.

Conclusion

Pru p 1 content is generally low in peach compared to Pru p 3. Several hypoallergenic Pru p 1 and Pru p 3 varieties, "Zi Xue Tao," "Wu Yue Xian," and "May Fire," were identified, which could be useful in trials for peach allergy patients.

The Di2/pet Variant in the PETALOSA Gene Underlies a Major Heat Requirement-Related QTL for Blooming Date in Peach [Prunus persica (L.) Batsch]

Environmental adaptation of deciduous fruit trees largely depends on their ability to synchronize growth and development with seasonal climate change. Winter dormancy of flower buds is a key process to prevent frost damage and ensure reproductive success. Temperature is a crucial environmental stimulus largely influencing the timing of flowering, only occurring after fulfillment of certain temperature requirements. Nevertheless, genetic variation affecting chilling or heat-dependent dormancy release still remains largely unknown. In this study, a major QTL able to delay blooming date in peach by increasing heat requirement was finely mapped in three segregating progenies, revealing a strict association with a genetic variant (petDEL) in a PETALOSA gene, previously shown to also affect flower morphology. Analysis of segregating genome-edited tobacco plants provided further evidence of the potential ability of PET variations to delay flowering time. Potential applications of the petDEL variant for improving phenological traits in peach are discussed.

Fine-tuning the performance of ddRAD-seq in the peach genome

The advance of Next Generation Sequencing (NGS) technologies allows high-throughput genotyping at a reasonable cost, although, in the case of peach, this technology has been scarcely developed. To date, only a standard Genotyping by Sequencing approach (GBS), based on a single restriction with ApeKI to reduce genome complexity, has been applied in peach. In this work, we assessed the performance of the double-digest RADseq approach (ddRADseq), by testing 6 double restrictions with the restriction profile generated with ApeKI. The enzyme pair PstI/MboI retained the highest number of loci in concordance with the in silico analysis. Under this condition, the analysis of a diverse germplasm collection (191 peach genotypes) yielded 200,759,000 paired-end (2 × 250 bp) reads that allowed the identification of 113,411 SNP, 13,661 InDel and 2133 SSR. We take advantage of a wide sample set to describe technical scope of the platform. The novel platform presented here represents a useful tool for genomic-based breeding for peach.

Genome-Wide Association Study of Brown Rot (Monilinia spp.) Tolerance in Peach

Brown rot, caused by Monilinia spp., is one of the most important diseases on stone fruit worldwide. Severe yield loss can be caused by pre- and post-harvest fruit decay. Although some degree of tolerance has been reported in peach and almond, the genetic resistance in peach cultivars is still lacking. To date, only few genomic regions associated with brown rot response in fruit skin and flesh have been detected in peach. Previous studies suggested brown rot tolerance in peach being a polygenic quantitative trait. More information is needed to uncover the genetics behind brown rot tolerance in peach. To identify the genomic regions in peach associated with this trait, 26 cultivars and progeny from 9 crosses with 'Bolinha' sources of tolerance, were phenotyped across two seasons (2015 and 2016) for brown rot disease severity index in wounded and non-wounded fruits and genotyped using a newly developed 9+9K peach SNP array. Genome wide association study using single- and multi-locus methods by GAPIT version 3, mrMLM 4.0, GAPIT and G Model, revealed 14 reliable SNPs significantly associated with brown rot infection responses in peach skin (10) and flesh (4) across whole genome except for chromosome 3. Candidate gene analysis within the haplotype regions of the detected markers identified 25 predicted genes associated with pathogen infection response/resistance. Results presented here facilitate further understanding of genetics behind brown rot tolerance in peach and provide an important foundation for DNA-assisted breeding.

Multi-Locus Genome-Wide Association Studies Reveal Fruit Quality Hotspots in Peach Genome

Peach is one of the most important fruit crops in the world, with the global annual production about 24.6 million tons. The United States is the fourth-largest producer after China, Spain, and Italy. Peach consumption has decreased over the last decade, most likely due to inconsistent quality of the fruit on the market. Thus, marker-assisted selection for fruit quality traits is highly desired in fresh market peach breeding programs and one of the major goals of the RosBREED project. The ability to use DNA information to select for desirable traits would enable peach breeders to efficiently plan crosses and select seedlings with desired quality traits early in the selection process before fruiting. Therefore, we assembled a multi-locus genome wide association study (GWAS) of 620 individuals from three public fresh market peach breeding programs (Arkansas, Texas, and South Carolina). The material was genotyped using 9K SNP array and the traits were phenotyped for three phenological (bloom date, ripening date, and days after bloom) and 11 fruit quality-related traits (blush, fruit diameter, fruit weight, adherence, fruit firmness, redness around pit, fruit texture, pit weight, soluble solid concentration, titratable acidity, and pH) over three seasons (2010, 2011, and 2012). Multi-locus association analyses, carried out using mrMLM 4.0 and FarmCPU R packages, revealed a total of 967 and 180 quantitative trait nucleotides (QTNs), respectively. Among the 88 consistently reliable QTNs detected using multiple multi-locus GWAS methods and/or at least two seasons, 44 were detected for the first time. Fruit quality hotspots were identified on chromosomes 1, 3, 4, 5, 6, and 8. Out of 566 candidate genes detected in the genomic regions harboring the QTN clusters, 435 were functionally annotated. Gene enrichment analyses revealed 68 different gene ontology (GO) terms associated with fruit quality traits. Data reported here advance our understanding of genetic mechanisms underlying important fruit quality traits and further support the development of DNA tools for breeding.

Recent Large-Scale Genotyping and Phenotyping of Plant Genetic Resources of Vegetatively Propagated Crops

Several recent national and international projects have focused on large-scale genotyping of plant genetic resources in vegetatively propagated crops like fruit and berries, potatoes and woody ornamentals. The primary goal is usually to identify true-to-type plant material, detect possible synonyms, and investigate genetic diversity and relatedness among accessions. A secondary goal may be to create sustainable databases that can be utilized in research and breeding for several years ahead. Commonly applied DNA markers (like microsatellite DNA and SNPs) and next-generation sequencing each have their pros and cons for these purposes. Methods for large-scale phenotyping have lagged behind, which is unfortunate since many commercially important traits (yield, growth habit, storability, and disease resistance) are difficult to score. Nevertheless, the analysis of gene action and development of robust DNA markers depends on environmentally controlled screening of very large sets of plant material. Although more time-consuming, co-operative projects with broad-scale data collection are likely to produce more reliable results. In this review, we will describe some of the approaches taken in genotyping and/or phenotyping projects concerning a wide variety of vegetatively propagated crops.

A candidate PpRPH gene of the D locus controlling fruit acidity in peach

A candidate gene, designate PpRPH, in the D locus was identified to control fruit acidity in peach. Fruit acidity has a strong impact on organoleptic quality of fruit. Peach fruit acidity is controlled by a large-effect D locus on chromosome 5. In this study, the D locus was mapped to a 509-kb interval, with two markers, 5dC720 and 5C1019, co-segregating with the non-acid/acid trait of peach fruit. Within this interval, a candidate gene encoding a putative small protein, designated PpRPH, showed a consistency between gene expression and fruit acidity, with up- and down-regulation in non-acidic and acidic fruits, respectively. Transient ectopic expression of PpRPH in tobacco leaves caused an increase of pH by approximately 40% compared to the control transformed with empty vector. Whereas, the concentrations of citrate and malate decreased significantly by 22% and 37%, respectively, with respect to the empty vector control. All these results suggest that PpRPH is a strong candidate gene of the D locus. These findings contribute to our overall understanding of the complex mechanism underlying fruit acidity in peach as well as that in other fruit crops.

Preliminary Identification of Key Genes Controlling Peach Pollen Fertility Using Genome-Wide Association Study

Previous genetic mapping helped detect a ~7.52 Mb putative genomic region for the pollen fertility trait on peach Chromosome 06 (Chr.06), which was too long for candidate gene characterization. In this study, using the whole-genome re-sequencing data of 201 peach accessions, we performed a genome-wide association study to identify key genes related to peach pollen fertility trait. The significant association peak was detected at Chr.06: 2,116,368 bp, which was in accordance with the previous genetic mapping results, but displayed largely improved precision, allowing for the identification of nine candidate genes. Among these candidates, gene PpABCG26, encoding an ATP-binding cassette G (ABCG) transporter and harboring the most significantly associated SNP (Single Nucleotide Polymorphism) marker in its coding region, was hypothesized to control peach pollen fertility/sterility based on the results of gene function comparison, gene relative expression, and nucleotide sequence analysis. The obtained results will help us to understand the genetic basis of peach pollen fertility trait, and to discover applicable markers for pre-selection in peach.

Genome structure variation analyses of peach reveal population dynamics and a 1.67 Mb causal inversion for fruit shape

BACKGROUND

Structural variations (SVs), a major resource of genomic variation, can have profound consequences on phenotypic variation, yet the impacts of SVs remain largely unexplored in crops.

RESULTS

Here, we generate a high-quality de novo genome assembly for a flat-fruit peach cultivar and produce a comprehensive SV map for peach, as a high proportion of genomic sequence is occupied by heterozygous SVs in the peach genome. We conduct population-level analyses that indicate SVs have undergone strong purifying selection during peach domestication, and find evidence of positive selection, with a significant preference for upstream and intronic regions during later peach improvement. We perform a SV-based GWAS that identifies a large 1.67-Mb heterozygous inversion that segregates perfectly with flat-fruit shape. Mechanistically, this derived allele alters the expression of the PpOFP2 gene positioned near the proximal breakpoint of the inversion, and we confirm in transgenic tomatoes that PpOFP2 is causal for flat-fruit shape.

CONCLUSIONS

Thus, beyond introducing new genomics resources for peach research, our study illustrates how focusing on SV data can drive basic functional discoveries in plant science.

Single Nucleotide Polymorphism Detection for Peach Gummosis Disease Resistance by Genome-Wide Association Study

Peach gummosis is one of the most widespread and destructive diseases. It causes growth stunting, yield loss, branch, trunk, and tree death, and is becoming a restrictive factor in healthy and sustainable development of peach production. Although a locus has been identified based on bi-parental quantitative trait locus (QTL) mapping, selection of gummosis-resistant cultivars remains challenging due to the lack of resistant parents and of the complexity of an inducing factor. In this study, an integrated approach of genome-wide association study (GWAS) and comparative transcriptome was used to elucidate the genetic architecture associated with the disease using 195 accessions and 145,456 genome-wide single nucleotide polymorphisms (SNPs). The broad-sense and narrow-sense heritabilities were estimated using 2-year phenotypic data and genotypic data, which gave high values of 70 and 73%, respectively. Evaluation of population structure by neighbor-joining and principal components analysis (PCA) clustered all accessions into three major groups and six subgroups, mainly according to fruit shape, hairy vs. glabrous fruit skin, pedigree, geographic origin, and domestication history. Five SNPs were found to be significantly associated with gummosis disease resistance, of which SNPrs285957, located on chromosome6 across 28 Mb, was detected by both the BLINK and the FarmCPU model. Six candidate genes flanked by or harboring the significant SNPs, previously implicated in biotic stress tolerance, were significantly associated with this resistance. Two highly resistant accessions were identified with low disease severity, which could be potential sources of resistance genes for breeding. Our results provide a fresh insight into the genetic control of peach gummosis disease.

The Multisite PeachRefPop Collection: A True Cultural Heritage and International Scientific Tool for Fruit Trees

Plants have evolved a range of adaptive mechanisms that adjust their development and physiology to variable external conditions, particularly in perennial species subjected to long-term interplay with the environment. Exploiting the allelic diversity within available germplasm and leveraging the knowledge of the mechanisms regulating genotype interaction with the environment are crucial to address climatic challenges and assist the breeding of novel cultivars with improved resilience. The development of multisite collections is of utmost importance for the conservation and utilization of genetic materials and will greatly facilitate the dissection of genotype-by-environment interaction. Such resources are still lacking for perennial trees, especially with the intrinsic difficulties of successful propagation, material exchange, and living collection maintenance. This work describes the concept, design, and realization of the first multisite peach (Prunus persica) reference collection (PeachRefPop) located across different European countries and sharing the same experimental design. Other than an invaluable tool for scientific studies in perennial species, PeachRefPop provides a milestone in an international collaborative project for the conservation and exploitation of European peach germplasm resources and, ultimately, as a true heritage for future generations.

Construction of a high-density SNP-based genetic map and identification of fruit-related QTLs and candidate genes in peach [Prunus persica (L.) Batsch]

BACKGROUND

High-density genetic mapping is a valuable tool for mapping loci that control specific traits for perennial fruit trees. Peach is an economically important fruit tree and a model Rosaceae species for genomic and genetic research. In peach, even though many molecular markers, genetic maps and QTL mappings have been reported, further research on the improvement of marker numbers, map densities, QTL accuracy and candidate gene identification is still warranted.

RESULTS

A high-density single nucleotide polymorphism (SNP)-based peach linkage map was constructed using specific locus amplified fragment sequencing (SLAF-seq). This genetic map consisted of 7998 SLAF markers, spanning 1098.79 cM with an average distance of 0.17 cM between adjacent markers. A total of 40 QTLs and 885 annotated candidate genes were detected for 10 fruit-related traits, including fruit weight (FW), fruit diameter (FD), percentage of red skin colour (PSC), eating quality (EQ), fruit flavour (FV), red in flesh (RF), red around pit (RP), adherence to pit (AP), fruit development period (FDP) and fruit fibre content (FFC). Eighteen QTLs for soluble solid content (SSC) were identified along LGs 1, 4, 5, and 6 in 2015 and 2016, and 540 genes were annotated in QTL intervals. Thirty-two QTLs for fruit acidity content (FA) were detected on LG1, and 2, 4, 5, 6, and 1232 candidate genes were identified. The expression profiles of 2 candidate genes for SSC and 4 for FA were analysed in parents and their offspring.

CONCLUSIONS

We constructed a high-density genetic map in peach based on SLAF-seq, which may contribute to the identification of important agronomic trait loci. Ninety QTLs for 12 fruit-related traits were identified, most of which overlapped with previous reports, and some new QTLs were obtained. A large number of candidate genes for fruit-related traits were screened and identified. These results may improve our understanding of the genetic control of fruit quality traits and provide useful information in marker-assisted selection for fruit quality in peach breeding programmes.

Resynthesis: Marker-Based Partial Reconstruction of Elite Genotypes in Clonally-Reproducing Plant Species

We propose a method for marker-based selection of cultivars of clonally-reproducing plant species which keeps the basic genetic architecture of a top-performing cultivar (usually a partly heterozygous genotype), with the addition of some agronomically relevant differences (such as production time, product appearance or quality), providing added value to the product or cultivation process. The method is based on selecting a) two complementary nearly-inbred lines from successive selfing generations (ideally only F2 and F3) of large size, that may generate individuals with most of their genome identical to the original cultivar but being homozygous for either of the two component haplotypes in the rest, and b) individuals with such characteristics already occurring in the F2. Option a) allows for introgressing genes from other individuals in one or both of these nearly-inbred lines. Peach, a woody-perennial, clonally-reproduced species, was chosen as a model for a proof of concept of the Resynthesis process due to its biological characteristics: self-compatibility, compact and genetically well-known genome, low recombination rates and relatively short intergeneration time (3-4 years). From 416 F2 seedlings from cultivar Sweet Dream (SD), we obtained seven individuals with 76-94% identity with SD, and selected five pairs of complementary lines with average homozygosity of the two parents ≥0.70 such that crossing would produce some individuals highly similar to SD. The application of this scheme to other species with more complex genomes or biological features, including its generalization to F1 hybrids, is discussed.

Whole genome re-sequencing of sweet cherry (Prunus avium L.) yields insights into genomic diversity of a fruit species

Sweet cherries, Prunus avium L. (Rosaceae), are gaining importance due to their perenniallity and nutritional attributes beneficial for human health. Interestingly, sweet cherry cultivars exhibit a wide range of phenotypic diversity in important agronomic traits, such as flowering time and defense reactions against pathogens. In this study, whole-genome resequencing (WGRS) was employed to characterize genetic variation, population structure and allelic variants in a panel of 20 sweet cherry and one wild cherry genotypes, embodying the majority of cultivated Greek germplasm and a representative of a local wild cherry elite phenotype. The 21 genotypes were sequenced in an average depth of coverage of 33.91×. and effective mapping depth, to the genomic reference sequence of 'Satonishiki' cultivar, between 22.21× to 36.62×. Discriminant analysis of principal components (DAPC) with SNPs revealed two clusters of genotypes. There was a rapid linkage disequilibrium decay, as the majority of SNP pairs with r2 in near complete disequilibrium (>0.8) were found at physical distances less than 10 kb. Functional analysis of the variants showed that the genomic ratio of non-synonymous/synonymous (dN/dS) changes was 1.78. The higher dN frequency in the Greek cohort of sweet cherry could be the result of artificial selection pressure imposed by breeding, in combination with the vegetative propagation of domesticated cultivars through grafting. The majority of SNPs with high impact (e.g., stop codon gaining, frameshift), were identified in genes involved in flowering time, dormancy and defense reactions against pathogens, providing promising resources for future breeding programs. Our study has established the foundation for further large scale characterization of sweet cherry germplasm, enabling breeders to incorporate diverse germplasm and allelic variants to fine tune flowering and maturity time and disease resistance in sweet cherry cultivars.

Maximising recombination across macadamia populations to generate linkage maps for genome anchoring

The Proteaceae genus Macadamia has a recent history of domestication as a commercial nut crop. We aimed to establish the first sequence-based haploid-correlated reference genetic linkage maps for this primarily outcrossing perennial tree crop, with marker density suitable for genome anchoring. Four first generation populations were used to maximise the segregation patterns available within full-sib, biparental and self-pollinated progeny. This allowed us to combine segregation data from overlapping subsets of >4,000 informative sequence-tagged markers to increase the effective coverage of the karyotype represented by the recombinant crossover events detected. All maps had 14 linkage groups, corresponding to the Macadamia haploid chromosome number, and enabled the anchoring and orientation of sequence scaffolds to construct a pseudo-chromosomal genome assembly for macadamia. Comparison of individual maps indicated a high level of congruence, with minor discrepancies satisfactorily resolved within the integrated maps. The combined set of maps significantly improved marker density and the proportion (70%) of the genome sequence assembly anchored. Overall, increasing our understanding of the genetic landscape and genome for this nut crop represents a substantial advance in macadamia genetics and genomics. The set of maps, large number of sequence-based markers and the reconstructed genome provide a toolkit to underpin future breeding that should help to extend the macadamia industry as well as provide resources for the long term conservation of natural populations in eastern Australia of this unique genus.

A High Resolution Melting Analysis-Based Genotyping Toolkit for the Peach (Prunus persica) Chilling Requirement

The chilling requirement (CR) is the main factor controlling the peach floral bud break and subsequent reproductive growth. To date, several peach CR quantitative trait loci (QTLs) have been identified. To improve the accessibility and convenience of this genetic information for peach breeders, the aim of this study was to establish an easy-to-use genotype screening system using peach CR molecular markers as a toolkit for marker-assisted selection. Here, we integrated 22 CR-associated markers from three published QTLs and positioned them on the Prunus persica physical map. Then, we built a PCR-based genotyping platform by using high-resolution melting (HRM) analysis with specific primers and trained this platform with 27 peach cultivars. Due to ambiguous variant calls from a commercial HRM software, we developed an R-based pipeline using principal component analysis (PCA) to accurately differentiate genotypes. Based on the PCA results, this toolkit was able to determine the genotypes at the CR-related single nucleotide polymorphisms (SNPs) in all tested peach cultivars. In this study, we showed that this HRM-PCA pipeline served as a low-cost, high-throughput, and non-gel genotyping solution. This system has great potential to accelerate CR-focused peach breeding.

Genome-Wide SNP discovery and genomic characterization in avocado (Persea americana Mill.)

Modern crop breeding is based on the use of genetically and phenotypically diverse plant material and, consequently, a proper understanding of population structure and genetic diversity is essential for the effective development of breeding programs. An example is avocado, a woody perennial fruit crop native to Mesoamerica with an increasing popularity worldwide. Despite its commercial success, there are important gaps in the molecular tools available to support on-going avocado breeding programs. In order to fill this gap, in this study, an avocado 'Hass' draft assembly was developed and used as reference to study 71 avocado accessions which represent the three traditionally recognized avocado horticultural races or subspecies (Mexican, Guatemalan and West Indian). An average of 5.72 M reads per individual and a total of 7,108 single nucleotide polymorphism (SNP) markers were produced for the 71 accessions analyzed. These molecular markers were used in a study of genetic diversity and population structure. The results broadly separate the accessions studied according to their botanical race in four main groups: Mexican, Guatemalan, West Indian and an additional group of Guatemalan × Mexican hybrids. The high number of SNP markers developed in this study will be a useful genomic resource for the avocado community.

Comparative population genomics identified genomic regions and candidate genes associated with fruit domestication traits in peach

Crop evolution is a long-term process involving selection by natural evolutionary forces and anthropogenic influences; however, the genetic mechanisms underlying the domestication and improvement of fruit crops have not been well studied to date. Here, we performed a population structure analysis in peach (Prunus persica) based on the genome-wide resequencing of 418 accessions and confirmed the presence of an obvious domestication event during evolution. We identified 132 and 106 selective sweeps associated with domestication and improvement, respectively. Analysis of their tissue-specific expression patterns indicated that the up-regulation of selection genes during domestication occurred mostly in fruit and seeds as opposed to other organs. However, during the improvement stage, more up-regulated selection genes were identified in leaves and seeds than in the other organs. Genome-wide association studies (GWAS) using 4.24 million single nucleotide polymorphisms (SNPs) revealed 171 loci associated with 26 fruit domestication traits. Among these loci, three candidate genes were highly associated with fruit weight and the sorbitol and catechin content in fruit. We demonstrated that as the allele frequency of the SNPs associated with high polyphenol composition decreased during peach evolution, alleles associated with high sugar content increased significantly. This indicates that there is genetic potential for the breeding of more nutritious fruit with enhanced bioactive polyphenols without disturbing a harmonious sugar and acid balance by crossing with wild species. This study also describes the development of the genomic resources necessary for evolutionary research in peach and provides the large-scale characterization of key agronomic traits in this crop species.

A new genomic tool for walnut (Juglans regia L.): development and validation of the high-density Axiom™ J. regia 700K SNP genotyping array

Over the last 20 years, global production of Persian walnut (Juglans regia L.) has grown enormously, likely reflecting increased consumption due to its numerous benefits to human health. However, advances in genome-wide association (GWA) studies and genomic selection (GS) for agronomically important traits in walnut remain limited due to the lack of powerful genomic tools. Here, we present the development and validation of a high-density 700K single nucleotide polymorphism (SNP) array in Persian walnut. Over 609K high-quality SNPs have been thoroughly selected from a set of 9.6 m genome-wide variants, previously identified from the high-depth re-sequencing of 27 founders of the Walnut Improvement Program (WIP) of University of California, Davis. To validate the effectiveness of the array, we genotyped a collection of 1284 walnut trees, including 1167 progeny of 48 WIP families and 26 walnut cultivars. More than half of the SNPs (55.7%) fell in the highest quality class of 'Poly High Resolution' (PHR) polymorphisms, which were used to assess the WIP pedigree integrity. We identified 151 new parent-offspring relationships, all confirmed with the Mendelian inheritance test. In addition, we explored the genetic variability among cultivars of different origin, revealing how the varieties from Europe and California were differentiated from Asian accessions. Both the reconstruction of the WIP pedigree and population structure analysis confirmed the effectiveness of the Applied Biosystems™ Axiom™ J. regia 700K SNP array, which initiates a novel genomic and advanced phase in walnut genetics and breeding.

Comparison Between Flat and Round Peaches, Genomic Evidences of Heterozygosity Events

Bud sports occur in many plant species, including fruit trees. Although they are correlated with genetic variance in somatic cells, the mechanisms responsible for bud sports are mostly unknown. In this study, a peach bud sport whose fruit shape was transformed to round from flat was identified by next generation sequencing (NGS), and we provide evidence that a long loss of heterozygosity (LOH) event may be responsible for this alteration in fruit shape. Moreover, compared to the reference genome, we identified 237,476 high quality single nucleotide polymorphisms (SNPs) in the wild-type and bud sport genomes. Using this SNP set, a long LOH event was identified at the distal end of scaffold Pp06 of the bud sport genome. Haplotypes from 155 additional peach accessions were phased, suggesting that the homozygous distal end of scaffold Pp06 of the bud sport was likely derived from only one haplotype of the wild-type flat peach. A genome-wide association study (GWAS) of 127 peach accessions was conducted to associate a SNP found at 26,924,482 bp of scaffold Pp06 to differences in fruit shape. All accessions with round-shaped fruit were found to have an A/A genotype, while those with A/T, or T/T genotypes had flat-shaped fruits. Finally, we also found that 236 peach accessions and 141 Prunus species with round-type fruit were found to have an A/A genotype at this SNP, while 22 flat peach accessions had an A/T genotype. Taken together, our results suggest that genes flanking this A/T polymorphism, and haplotyped carrying the T allele may determine flat fruit shape in this population. Furthermore, the LOH event resulting in the loss of the haplotype carrying the T allele may therefore be responsible for fruit shape alteration in wild-type flat peach.

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.

Prunus genetics and applications after de novo genome sequencing: achievements and prospects

Prior to the availability of whole-genome sequences, our understanding of the structural and functional aspects of Prunus tree genomes was limited mostly to molecular genetic mapping of important traits and development of EST resources. With public release of the peach genome and others that followed, significant advances in our knowledge of Prunus genomes and the genetic underpinnings of important traits ensued. In this review, we highlight key achievements in Prunus genetics and breeding driven by the availability of these whole-genome sequences. Within the structural and evolutionary contexts, we summarize: (1) the current status of Prunus whole-genome sequences; (2) preliminary and ongoing work on the sequence structure and diversity of the genomes; (3) the analyses of Prunus genome evolution driven by natural and man-made selection; and (4) provide insight into haploblocking genomes as a means to define genome-scale patterns of evolution that can be leveraged for trait selection in pedigree-based Prunus tree breeding programs worldwide. Functionally, we summarize recent and ongoing work that leverages whole-genome sequences to identify and characterize genes controlling 22 agronomically important Prunus traits. These include phenology, fruit quality, allergens, disease resistance, tree architecture, and self-incompatibility. Translationally, we explore the application of sequence-based marker-assisted breeding technologies and other sequence-guided biotechnological approaches for Prunus crop improvement. Finally, we present the current status of publically available Prunus genomics and genetics data housed mainly in the Genome Database for Rosaceae (GDR) and its updated functionalities for future bioinformatics-based Prunus genetics and genomics inquiry.

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.

Identification of a new allele of the Dw gene causing brachytic dwarfing in peach

Objective

Peach brachytic dwarfism determined by Dwarf gene (Dw) is an undesired trait segregating in some peach breeding programs. Recently, a single nucleotide polymorphism (SNP) mutation in the gibberellin insensitive dwarf 1 (GID1) peach gene causing brachytic dwarfism was described. In this research we wanted to validate this marker in an F₂ population of the ‘Nectavantop’ peach cultivar (Nv) to include it as a marker assisted selection tool for peach breeding programs.

Results

The observed segregation of the trait was in agreement with that of a recessive gene, the individuals homozygous for the recessive allele (dwdw) presenting the dwarf genotype. Dw was mapped to the distal part of linkage group 6 as previously described. The SNP marker based on the causal mutation previously described did not segregate in Nv F₂ population. The sequence of the GID1c gene in Nv revealed a second SNP in its coding sequence which cosegregated with the dwarf phenotype. This SNP was predicted by the SNAP² software to cause a major functional change and was validated in the dwarf peach cultivar ‘Small sunning’. These results suggest the existence of at least two independent mutations of the Dw gene causing the peach brachytic dwarf phenotype.

Electronic supplementary material

The online version of this article (10.1186/s13104-018-3490-7) contains supplementary material, which is available to authorized users.

Integrative genomics approaches validate PpYUC11-like as candidate gene for the stony hard trait in peach (P. persica L. Batsch)

BACKGROUND

Texture is one of the most important fruit quality attributes. In peach, stony hard (SH) is a recessive monogenic trait (hd/hd) that confers exceptionally prolonged firm flesh to fully ripe fruit. Previous studies have shown that the SH mutation affects the fruit ability to synthesize appropriate amounts of indol-3-acetic acid (IAA), which orchestrates the ripening processes through the activation of system 2 ethylene pathway. Allelic variation in a TC microsatellite located within the first intron of PpYUC11-like (a YUCCA-like auxin-biosynthesis gene) has been recently proposed as the causal mutation of the SH phenotype.

RESULTS

The simple genetic determinism of the SH trait has been clarified through genome-wide association and LD analyses in a diverse set of accessions, restricting the hd locus to an interval of about 1.8 Mbp in chromosome 6. The comparison of fruit transcriptome data from non-SH (melting flesh) and SH accessions provided an expression patterns overview of the annotated transcripts within the hd locus, confirming the absence of PpYUC11-like expression in SH fruits. To explore further possible associations between genomic variants at the hd locus and the SH phenotype, re-sequencing data of the SH accession 'D41-62' were compared with several SH and non-SH accessions with different genetic backgrounds. A further step of validation was provided through the evaluation of variant-trait association in two bi-parental F₂ populations issued from the SH accession 'D41-62' and a panel of advanced breeding selections, showing perfect co-segregation of the PpYUC11-like intron TC₂₀ allele and the SH phenotype.

CONCLUSIONS

In this study, we provide a multi-level validation of the genetic control of the SH trait through the integration of genome-wide association mapping, transcriptome analysis and whole-genome resequencing data for SH and non-SH accessions, and marker-trait association in a panel of advanced breeding selections and segregating progenies. Collectively, our data confirm with high confidence the role of allelic variation at PpYUC11-like locus as the genetic determinant of the SH trait, opening interesting perspectives at both biological and applied research level.

An integrated approach for increasing breeding efficiency in apple and peach in Europe

Despite the availability of whole genome sequences of apple and peach, there has been a considerable gap between genomics and breeding. To bridge the gap, the European Union funded the FruitBreedomics project (March 2011 to August 2015) involving 28 research institutes and private companies. Three complementary approaches were pursued: (i) tool and software development, (ii) deciphering genetic control of main horticultural traits taking into account allelic diversity and (iii) developing plant materials, tools and methodologies for breeders. Decisive breakthroughs were made including the making available of ready-to-go DNA diagnostic tests for Marker Assisted Breeding, development of new, dense SNP arrays in apple and peach, new phenotypic methods for some complex traits, software for gene/QTL discovery on breeding germplasm via Pedigree Based Analysis (PBA). This resulted in the discovery of highly predictive molecular markers for traits of horticultural interest via PBA and via Genome Wide Association Studies (GWAS) on several European genebank collections. FruitBreedomics also developed pre-breeding plant materials in which multiple sources of resistance were pyramided and software that can support breeders in their selection activities. Through FruitBreedomics, significant progresses were made in the field of apple and peach breeding, genetics, genomics and bioinformatics of which advantage will be made by breeders, germplasm curators and scientists. A major part of the data collected during the project has been stored in the FruitBreedomics database and has been made available to the public. This review covers the scientific discoveries made in this major endeavour, and perspective in the apple and peach breeding and genomics in Europe and beyond.

Genetic dissection of Sharka disease tolerance in peach (P. persica L. Batsch)

BACKGROUND

Plum pox virus (PPV), agent of Sharka disease, is the most important quarantine pathogen of peach (P. persica L. Batsch). Extensive evaluation of peach germplasm has highlighted the lack of resistant sources, while suggesting the presence of a quantitative disease resistance, expressed as reduction in the intensity of symptoms. Unravelling the genetic architecture of peach response to PPV infection is essential for pyramiding resistant genes and for developing more tolerant varieties. For this purpose, a genome-wide association (GWA) approach was applied in a panel of accessions phenotyped for virus susceptibility and genotyped with the IPSC peach 9 K SNP Array, and coupled with an high-coverage resequencing of the tolerant accession 'Kamarat'.

RESULTS

Genome-wide association identified three highly significant associated loci on chromosome 2 and 3, accounting for most of the reduction in PPV-M susceptibility within the analysed peach population. The exploration of associated intervals through whole-genome comparison of the tolerant accession 'Kamarat' and other susceptible accessions, including the PPV-resistant wild-related species P. davidiana, allow the identification of allelic variants in promising candidate genes, including an RTM2-like gene already characterized in A. thaliana.

CONCLUSIONS

The present study is the first effort to identify genetic factors involved in Sharka disease in peach germplasm through a GWA approach. We provide evidence of the presence of quantitative resistant loci in a collection of peach accessions, identifying major loci and highly informative SNPs that could be useful for marker assisted selection. These results could serve as reference bases for future research aimed at the comprehension of genetic mechanism regulating the complex peach-PPV interaction.

A deletion affecting an LRR-RLK gene co-segregates with the fruit flat shape trait in peach

In peach, the flat phenotype is caused by a partially dominant allele in heterozygosis (Ss), fruits from homozygous trees (SS) abort a few weeks after fruit setting. Previous research has identified a SSR marker (UDP98-412) highly associated with the trait, found suitable for marker assisted selection (MAS). Here we report a ∼10 Kb deletion affecting the gene PRUPE.6G281100, 400 Kb upstream of UDP98-412, co-segregating with the trait. This gene is a leucine-rich repeat receptor-like kinase (LRR-RLK) orthologous to the Brassinosteroid insensitive 1-associated receptor kinase 1 (BAK1) group. PCR markers suitable for MAS confirmed its strong association with the trait in a collection of 246 cultivars. They were used to evaluate the DNA from a round fruit derived from a somatic mutation of the flat variety 'UFO-4', revealing that the mutation affected the flat associated allele (S). Protein BLAST alignment identified significant hits with genes involved in different biological processes. Best protein hit occurred with AtRLP12, which may functionally complement CLAVATA2, a key regulator that controls the stem cell population size. RT-PCR analysis revealed the absence of transcription of the partially deleted allele. The data support PRUPE.6G281100 as a candidate gene for flat shape in peach.

Integrated QTL detection for key breeding traits in multiple peach progenies

Background

Peach (Prunus persica (L.) Batsch) is a major temperate fruit crop with an intense breeding activity. Breeding is facilitated by knowledge of the inheritance of the key traits that are often of a quantitative nature. QTLs have traditionally been studied using the phenotype of a single progeny (usually a full-sib progeny) and the correlation with a set of markers covering its genome. This approach has allowed the identification of various genes and QTLs but is limited by the small numbers of individuals used and by the narrow transect of the variability analyzed. In this article we propose the use of a multi-progeny mapping strategy that used pedigree information and Bayesian approaches that supports a more precise and complete survey of the available genetic variability.

Results

Seven key agronomic characters (data from 1 to 3 years) were analyzed in 18 progenies from crosses between occidental commercial genotypes and various exotic lines including accessions of other Prunus species. A total of 1467 plants from these progenies were genotyped with a 9 k SNP array. Forty-seven QTLs were identified, 22 coinciding with major genes and QTLs that have been consistently found in the same populations when studied individually and 25 were new. A substantial part of the QTLs observed (47%) would not have been detected in crosses between only commercial materials, showing the high value of exotic lines as a source of novel alleles for the commercial gene pool. Our strategy also provided estimations on the narrow sense heritability of each character, and the estimation of the QTL genotypes of each parent for the different QTLs and their breeding value.

Conclusions

The integrated strategy used provides a broader and more accurate picture of the variability available for peach breeding with the identification of many new QTLs, information on the sources of the alleles of interest and the breeding values of the potential donors of such valuable alleles. These results are first-hand information for breeders and a step forward towards the implementation of DNA-informed strategies to facilitate selection of new cultivars with improved productivity and quality.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3783-6) contains supplementary material, which is available to authorized users.

Genome-enabled predictions for fruit weight and quality from repeated records in European peach progenies

Background

Highly polygenic traits such as fruit weight, sugar content and acidity strongly influence the agroeconomic value of peach varieties. Genomic Selection (GS) can accelerate peach yield and quality gain if predictions show higher levels of accuracy compared to phenotypic selection. The available IPSC 9K SNP array V1 allows standardized and highly reliable genotyping, preparing the ground for GS in peach.

Results

A repeatability model (multiple records per individual plant) for genome-enabled predictions in eleven European peach populations is presented. The analysis included 1147 individuals derived from both commercial and non-commercial peach or peach-related accessions. Considered traits were average fruit weight (FW), sugar content (SC) and titratable acidity (TA). Plants were genotyped with the 9K IPSC array, grown in three countries (France, Italy, Spain) and phenotyped for 3–5 years. An analysis of imputation accuracy of missing genotypic data was conducted using the software Beagle, showing that two of the eleven populations were highly sensitive to increasing levels of missing data. The regression model produced, for each trait and each population, estimates of heritability (FW:0.35, SC:0.48, TA:0.53, on average) and repeatability (FW:0.56, SC:0.63, TA:0.62, on average). Predictive ability was estimated in a five-fold cross validation scheme within population as the correlation of true and predicted phenotypes. Results differed by populations and traits, but predictive abilities were in general high (FW:0.60, SC:0.72, TA:0.65, on average).

Conclusions

This study assessed the feasibility of Genomic Selection in peach for highly polygenic traits linked to yield and fruit quality. The accuracy of imputing missing genotypes was as high as 96%, and the genomic predictive ability was on average 0.65, but could be as high as 0.84 for fruit weight or 0.83 for titratable acidity. The estimated repeatability may prove very useful in the management of the typical long cycles involved in peach productions. All together, these results are very promising for the application of genomic selection to peach breeding programmes.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3781-8) contains supplementary material, which is available to authorized users.

Reduced representation genome sequencing reveals patterns of genetic diversity and selection in apple

Identifying DNA sequence variations is a fundamental step towards deciphering the genetic basis of traits of interest. Here, a total of 20 cultivated and 10 wild apples were genotyped using specific-locus amplified fragment sequencing, and 39,635 single nucleotide polymorphisms with no missing genotypes and evenly distributed along the genome were selected to investigate patterns of genome-wide genetic variations between cultivated and wild apples. Overall, wild apples displayed higher levels of genetic diversity than cultivated apples. Linkage disequilibrium (LD) decays were observed quite rapidly in cultivated and wild apples, with an r² -value below 0.2 at 440 and 280 bp, respectively. Moreover, bidirectional gene flow and different distribution patterns of LD blocks were detected between domesticated and wild apples. Most LD blocks unique to cultivated apples were located within QTL regions controlling fruit quality, thus suggesting that fruit quality had probably undergone selection during apple domestication. The genome of the earliest cultivated apple in China, Nai, was highly similar to that of Malus sieversii, and contained a small portion of genetic material from other wild apple species. This suggested that introgression could have been an important driving force during initial domestication of apple. These findings will facilitate future breeding and genetic dissection of complex traits in apple.

Deciphering the genetic control of fruit texture in apple by multiple family-based analysis and genome-wide association

Fruit texture is a complex feature composed of mechanical and acoustic properties relying on the modifications occurring in the cell wall throughout fruit development and ripening. Apple is characterized by a large variation in fruit texture behavior that directly impacts both the consumer's appreciation and post-harvest performance. To decipher the genetic control of fruit texture comprehensively, two complementing quantitative trait locus (QTL) mapping approaches were employed. The first was represented by a pedigree-based analysis (PBA) carried out on six full-sib pedigreed families, while the second was a genome-wide association study (GWAS) performed on a collection of 233 apple accessions. Both plant materials were genotyped with a 20K single nucleotide polymorphism (SNP) array and phenotyped with a sophisticated high-resolution texture analyzer. The overall QTL results indicated the fundamental role of chromosome 10 in controlling the mechanical properties, while chromosomes 2 and 14 were more associated with the acoustic response. The latter QTL, moreover, showed a consistent relationship between the QTL-estimated genotypes and the acoustic performance assessed among seedlings. The in silico annotation of these intervals revealed interesting candidate genes potentially involved in fruit texture regulation, as suggested by the gene expression profile. The joint integration of these approaches sheds light on the specific control of fruit texture, enabling important genetic information to assist in the selection of valuable fruit quality apple varieties.

Mapping the sensory perception of apple using descriptive sensory evaluation in a genome wide association study

Breeding apples is a long-term endeavour and it is imperative that new cultivars are selected to have outstanding consumer appeal. This study has taken the approach of merging sensory science with genome wide association analyses in order to map the human perception of apple flavour and texture onto the apple genome. The goal was to identify genomic associations that could be used in breeding apples for improved fruit quality. A collection of 85 apple cultivars was examined over two years through descriptive sensory evaluation by a trained sensory panel. The trained sensory panel scored randomized sliced samples of each apple cultivar for seventeen taste, flavour and texture attributes using controlled sensory evaluation practices. In addition, the apple collection was subjected to genotyping by sequencing for marker discovery. A genome wide association analysis suggested significant genomic associations for several sensory traits including juiciness, crispness, mealiness and fresh green apple flavour. The findings include previously unreported genomic regions that could be used in apple breeding and suggest that similar sensory association mapping methods could be applied in other plants.

Comparative Transcriptome and Microscopy Analyses Provide Insights into Flat Shape Formation in Peach (Prunus persica)

Fruit shape is an important external characteristic that consumers use to select preferred fruit cultivars. In peach, the flat fruit cultivars have become more and more popular worldwide. Genetic markers closely linking to the flat fruit trait have been identified and are useful for marker-assisted breeding. However, the cellular and genetic mechanisms underpinning flat fruit formation are still poorly understood. In this study, we have revealed the differences in fruit cell number, cell size, and in gene expression pattern between the traditional round fruit and modern flat fruit cultivars. Flat peach cultivars possessed significantly lower number of cells in the vertical axis because cell division in the vertical direction stopped early in the flat fruit cultivars at 15 DAFB (day after full bloom) than in round fruit cultivars at 35 DAFB. This resulted in the reduction in vertical development in the flat fruit. Significant linear relationship was observed between fruit vertical diameter and cell number in vertical axis for the four examined peach cultivars (R² = 0.9964) at maturation stage, and was also observed between fruit vertical diameter and fruit weight (R² = 0.9605), which indicated that cell number in vertical direction contributed to the flat shape formation. Furthermore, in RNA-seq analysis, 4165 differentially expressed genes (DEGs) were detected by comparing RNA-seq data between flat and round peach cultivars at different fruit development stages. In contrast to previous studies, we discovered 28 candidate genes potentially responsible for the flat shape formation, including 19 located in the mapping site and 9 downstream genes. Our study indicates that flat and round fruit shape in peach is primarily determined by the regulation of cell production in the vertical direction during early fruit development.

Genome-Wide Association Mapping of Flowering and Ripening Periods in Apple

Deciphering the genetic control of flowering and ripening periods in apple is essential for breeding cultivars adapted to their growing environments. We implemented a large Genome-Wide Association Study (GWAS) at the European level using an association panel of 1,168 different apple genotypes distributed over six locations and phenotyped for these phenological traits. The panel was genotyped at a high-density of SNPs using the Axiom®Apple 480 K SNP array. We ran GWAS with a multi-locus mixed model (MLMM), which handles the putatively confounding effect of significant SNPs elsewhere on the genome. Genomic regions were further investigated to reveal candidate genes responsible for the phenotypic variation. At the whole population level, GWAS retained two SNPs as cofactors on chromosome 9 for flowering period, and six for ripening period (four on chromosome 3, one on chromosome 10 and one on chromosome 16) which, together accounted for 8.9 and 17.2% of the phenotypic variance, respectively. For both traits, SNPs in weak linkage disequilibrium were detected nearby, thus suggesting the existence of allelic heterogeneity. The geographic origins and relationships of apple cultivars accounted for large parts of the phenotypic variation. Variation in genotypic frequency of the SNPs associated with the two traits was connected to the geographic origin of the genotypes (grouped as North+East, West and South Europe), and indicated differential selection in different growing environments. Genes encoding transcription factors containing either NAC or MADS domains were identified as major candidates within the small confidence intervals computed for the associated genomic regions. A strong microsynteny between apple and peach was revealed in all the four confidence interval regions. This study shows how association genetics can unravel the genetic control of important horticultural traits in apple, as well as reduce the confidence intervals of the associated regions identified by linkage mapping approaches. Our findings can be used for the improvement of apple through marker-assisted breeding strategies that take advantage of the accumulating additive effects of the identified SNPs.

DNA-informed breeding of rosaceous crops: promises, progress and prospects

Crops of the Rosaceae family provide valuable contributions to rural economies and human health and enjoyment. Sustained solutions to production challenges and market demands can be met with genetically improved new cultivars. Traditional rosaceous crop breeding is expensive and time-consuming and would benefit from improvements in efficiency and accuracy. Use of DNA information is becoming conventional in rosaceous crop breeding, contributing to many decisions and operations, but only after past decades of solved challenges and generation of sufficient resources. Successes in deployment of DNA-based knowledge and tools have arisen when the 'chasm' between genomics discoveries and practical application is bridged systematically. Key steps are establishing breeder desire for use of DNA information, adapting tools to local breeding utility, identifying efficient application schemes, accessing effective services in DNA-based diagnostics and gaining experience in integrating DNA information into breeding operations and decisions. DNA-informed germplasm characterization for revealing identity and relatedness has benefitted many programs and provides a compelling entry point to reaping benefits of genomics research. DNA-informed germplasm evaluation for predicting trait performance has enabled effective reallocation of breeding resources when applied in pioneering programs. DNA-based diagnostics is now expanding from specific loci to genome-wide considerations. Realizing the full potential of this expansion will require improved accuracy of predictions, multi-trait DNA profiling capabilities, streamlined breeding information management systems, strategies that overcome plant-based features that limit breeding progress and widespread training of current and future breeding personnel and allied scientists.

Evolutionary Genomics of Peach and Almond Domestication

The domesticated almond [Prunus dulcis (L.) Batsch] and peach [P. persica (Mill.) D. A. Webb] originated on opposite sides of Asia and were independently domesticated ∼5000 yr ago. While interfertile, they possess alternate mating systems and differ in a number of morphological and physiological traits. Here, we evaluated patterns of genome-wide diversity in both almond and peach to better understand the impacts of mating system, adaptation, and domestication on the evolution of these taxa. Almond has around seven times the genetic diversity of peach, and high genome-wide FST values support their status as separate species. We estimated a divergence time of ∼8 MYA (million years ago), coinciding with an active period of uplift in the northeast Tibetan Plateau and subsequent Asian climate change. We see no evidence of a bottleneck during domestication of either species, but identify a number of regions showing signatures of selection during domestication and a significant overlap in candidate regions between peach and almond. While we expected gene expression in fruit to overlap with candidate selected regions, instead we find enrichment for loci highly differentiated between the species, consistent with recent fossil evidence suggesting fruit divergence long preceded domestication. Taken together, this study tells us how closely related tree species evolve and are domesticated, the impact of these events on their genomes, and the utility of genomic information for long-lived species. Further exploration of this data will contribute to the genetic knowledge of these species and provide information regarding targets of selection for breeding application, and further the understanding of evolution in these species.

Application of Genomic Technologies to the Breeding of Trees

The recent introduction of next generation sequencing (NGS) technologies represents a major revolution in providing new tools for identifying the genes and/or genomic intervals controlling important traits for selection in breeding programs. In perennial fruit trees with long generation times and large sizes of adult plants, the impact of these techniques is even more important. High-throughput DNA sequencing technologies have provided complete annotated sequences in many important tree species. Most of the high-throughput genotyping platforms described are being used for studies of genetic diversity and population structure. Dissection of complex traits became possible through the availability of genome sequences along with phenotypic variation data, which allow to elucidate the causative genetic differences that give rise to observed phenotypic variation. Association mapping facilitates the association between genetic markers and phenotype in unstructured and complex populations, identifying molecular markers for assisted selection and breeding. Also, genomic data provide in silico identification and characterization of genes and gene families related to important traits, enabling new tools for molecular marker assisted selection in tree breeding. Deep sequencing of transcriptomes is also a powerful tool for the analysis of precise expression levels of each gene in a sample. It consists in quantifying short cDNA reads, obtained by NGS technologies, in order to compare the entire transcriptomes between genotypes and environmental conditions. The miRNAs are non-coding short RNAs involved in the regulation of different physiological processes, which can be identified by high-throughput sequencing of RNA libraries obtained by reverse transcription of purified short RNAs, and by in silico comparison with known miRNAs from other species. All together, NGS techniques and their applications have increased the resources for plant breeding in tree species, closing the former gap of genetic tools between trees and annual species.