Fine mapping and identification of a candidate gene for a major locus controlling maturity date in peach

Tandem transcription factors PpNAC1 and PpNAC5 synergistically activate the transcription of the PpPGF to regulate peach softening during fruit ripening

Peach fruit rapidly soften after harvest, a significant challenge for producers and marketers as it results in rotting fruit and significantly reduces shelf life. In this study, we identified two tandem genes, PpNAC1 and PpNAC5, within the sr (slow ripening) locus. Phylogenetic analysis showed that NAC1 and NAC5 are highly conserved in dicots and that PpNAC1 is the orthologous gene of Non-ripening (NOR) in tomato. PpNAC1 and PpNAC5 were highly expressed in peach fruit, with their transcript levels up-regulated at the onset of ripening. Yeast two-hybrid and bimolecular fluorescence complementation assays showed PpNAC1 interacting with PpNAC5 and this interaction occurs with the tomato and apple orthologues. Transient gene silencing experiments showed that PpNAC1 and PpNAC5 positively regulate peach fruit softening. Yeast one-hybrid and dual luciferase assays and LUC bioluminescence imaging proved that PpNAC1 and PpNAC5 directly bind to the PpPGF promoter and activate its transcription. Co-expression of PpNAC1 and PpNAC5 showed higher levels of PpPGF activation than expression of PpNAC1 or PpNAC5 alone. In summary, our findings demonstrate that the tandem transcription factors PpNAC1 and PpNAC5 synergistically activate the transcription of PpPGF to regulate fruit softening during peach fruit ripening.

Ppe.RPT/SSC-1: from QTL mapping to a predictive KASP test for ripening time and soluble solids concentration in peach

Genomic regions associated with ripening time (RPT) and soluble solids concentration (SSC) were mapped using a pedigreed population including multiple F1 and F2 families from the Clemson University peach breeding program (CUPBP). RPT and SSC QTLs were consistently identified in two seasons (2011 and 2012) and the average datasets (average of two seasons). A target region spanning 10,981,971-11,298,736 bp on chromosome 4 of peach reference genome used for haplotype analysis revealed four haplotypes with significant differences in trait values among different diplotype combinations. Favorable alleles at the target region for both RPT and SSC were determined and a DNA test for predicting RPT and SSC was developed. Two Kompetitive Allele Specific PCR (KASP) assays were validated on 84 peach cultivars and 163 seedlings from the CUPBP, with only one assay (Ppe.RPT/SSC-1) needed to predict between early and late-season ripening cultivars and low and high SSC. These results advance our understanding of the genetic basis of RPT and SSC and facilitate selection of new peach cultivars with the desired RPT and SSC.

Genome variation and LTR-RT analyses of an ancient peach landrace reveal mechanism of blood-flesh fruit color formation and fruit maturity date advancement

Peach (Prunus persica) landrace has typical regional characteristics, strong environmental adaptability, and contains many valuable genes that provide the foundation for breeding excellent varieties. Therefore, it is necessary to assemble the genomes of specific landraces to facilitate the localization and utilization of these genes. Here, we de novo assembled a high-quality genome from an ancient blood-fleshed Chinese landrace Tianjin ShuiMi (TJSM) that originated from the China North Plain. The assembled genome size was 243.5 Mb with a contig N50 of 23.7 Mb and a scaffold N50 of 28.6 Mb. Compared with the reported peach genomes, our assembled TJSM genome had the largest number of specific structural variants (SVs) and long terminal repeat-retrotransposons (LTR-RTs). Among the LTR-RTs with the potential to regulate their host genes, we identified a 6688 bp LTR-RT (named it blood TE) in the promoter of NAC transcription factor-encoding PpBL, a gene regulating peach blood-flesh formation. The blood TE was not only co-separated with the blood-flesh phenotype but also associated with fruit maturity date advancement and different intensities of blood-flesh color formation. Our findings provide new insights into the mechanism underlying the development of the blood-flesh color and determination of fruit maturity date and highlight the potential of the TJSM genome to mine more variations related to agronomic traits in peach fruit.

Two adjacent NAC transcription factors regulate fruit maturity date and flavor in peach

Although maturity date (MD) is an essential factor affecting fresh fruit marketing and has a pleiotropic effect on fruit taste qualities, the underlying mechanisms remain largely unclear. In this study, we functionally characterized two adjacent NAM-ATAF1/2-CUC2 (NAC) transcription factors (TFs), PpNAC1 and PpNAC5, both of which were associated with fruit MD in peach. PpNAC1 and PpNAC5 were found capable of activating transcription of genes associated with cell elongation, cell wall degradation and ethylene biosynthesis, suggesting their regulatory roles in fruit enlargement and ripening. Furthermore, PpNAC1 and PpNAC5 had pleiotropic effects on fruit taste due to their ability to activate transcription of genes for sugar accumulation and organic acid degradation. Interestingly, both PpNAC1 and PpNAC5 orthologues were found in fruit-producing angiosperms and adjacently arranged in all 91 tested dicots but absent in fruitless gymnosperms, suggesting their important roles in fruit development. Our results provide insight into the regulatory roles of NAC TFs in MD and fruit taste.

Discovery of a key gene associated with fruit maturity date and analysis of its regulatory pathway in peach

Fruit maturity is an important agronomic trait of fruit crops. Although in previous studies, several molecular markers are developed for the trait, the knowledge about its candidate genes is particularly limited. In this study, a total of 357 peach accessions were re-sequenced to obtain 949,638 SNPs. Combing with 3-year fruit maturity dates, a genome-wide association analysis was performed, and 5, 8, and 9 association loci were identified. To screen the candidate genes for those year-stable loci on chromosomes 4 and 5, two maturity date mutants were used for transcriptome sequencing. Gene expression analysis indicated that Prupe.4G186800 and Prupe.4G187100 on chromosome 4 were essential to fruit ripening in peaches. However, the expression analysis of different tissues showed that the first gene has no tissue-specific character, but transgenic studies showed that the latter is more likely to be a key candidate gene than the first for the maturity date in peach. The yeast two-hybrid assay showed that the proteins encoded by the two genes interacted and then regulated fruit ripening. Moreover, the previously identified 9bp insertion in Prupe.4G186800 may affect their interaction ability. This research is of great significance for understanding the molecular mechanism of peach fruit ripening and developing practical molecular markers in a breeding program.

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.

Inheritance of Fruit Red-Flesh Patterns in Peach

Fruit color is an important trait in peach from the point of view of consumer preference, nutritional content, and diversification of fruit typologies. Several genes and phenotypes have been described for peach flesh and skin color, and although peach color knowledge has increased in the last few years, some fruit color patterns observed in peach breeding programs have not been carefully described. In this work, we first describe some peach mesocarp color patterns that have not yet been described in a collection of commercial peach cultivars, and we also study the genetic inheritance of the red dots present in the flesh (RDF) and red color around the stone (CAS) in several intra- and interspecific segregating populations for both traits. For RDF, we identified a QTL at the beginning of G5 in two intraspecific populations, and for CAS we identified a major QTL in G4 in both an intraspecific and an interspecific population between almond and peach. Finally, we discuss the interaction between these QTLs and some other genes previously identified in peach, such as dominant blood flesh (DBF), color around the stone (Cs), subacid (D) and the maturity date (MD), and the implications for peach breeding. The results obtained here will help peach germplasm curators and breeders to better characterize their plant materials and to develop an integrated system of molecular markers to select these traits.

Comparative network analysis reveals the dynamics of organic acid diversity during fruit ripening in peach (Prunus persica L. Batsch)

BACKGROUND

Organic acids are important components that determine the fruit flavor of peach (Prunus persica L. Batsch). However, the dynamics of organic acid diversity during fruit ripening and the key genes that modulate the organic acids metabolism remain largely unknown in this kind of fruit tree which yield ranks sixth in the world.

RESULTS

In this study, we used 3D transcriptome data containing three dimensions of information, namely time, phenotype and gene expression, from 5 different varieties of peach to construct gene co-expression networks throughout fruit ripening of peach. With the network inferred, the time-ordered network comparative analysis was performed to select high-acid specific gene co-expression network and then clarify the regulatory factors controlling organic acid accumulation. As a result, network modules related to organic acid synthesis and metabolism under high-acid and low-acid comparison conditions were identified for our following research. In addition, we obtained 20 candidate genes as regulatory factors related to organic acid metabolism in peach.

CONCLUSIONS

The study provides new insights into the dynamics of organic acid accumulation during fruit ripening, complements the results of classical co-expression network analysis and establishes a foundation for key genes discovery from time-series multiple species transcriptome data.

Transcriptome and Gene Regulatory Network Analyses Reveal New Transcription Factors in Mature Fruit Associated with Harvest Date in Prunus persica

Harvest date is a critical parameter for producers and consumers regarding agro-industrial performance. It involves a pleiotropic effect controlling the development of other fruit quality traits through finely controlling regulatory mechanisms. Fruit ripening is a process in which various signals and biological events co-occur and are regulated by hormone signaling that produces the accumulation/degradation of multiple compounds. However, the regulatory mechanisms that control the hormone signaling involved in fruit development and ripening are still unclear. To investigate the issue, we used individuals with early, middle and late harvest dates from a peach segregating population to identify regulatory candidate genes controlling fruit quality traits at the harvest stage and validate them in contrasting peach varieties for this trait. We identified 467 and 654 differentially expressed genes for early and late harvest through a transcriptomic approach. In addition, using the Arabidopsis DAP-seq database and network analysis, six transcription factors were selected. Our results suggest significant hormonal balance and cell wall composition/structure differences between early and late harvest samples. Thus, we propose that higher expression levels of the transcription factors HB7, ERF017 and WRKY70 in early harvest individuals would induce the expression of genes associated with the jasmonic acid pathway, photosynthesis and gibberellins inhibition. While on the other hand, the high expression levels of LHY, CDF3 and NAC083 in late harvest individuals would promote the induction of genes associated with abscisic acid biosynthesis, auxins and cell wall remodeling.

Genetic dissection of fruit maturity date in apricot (P. armeniaca L.) through a Single Primer Enrichment Technology (SPET) approach

Background

Single primer enrichment technology (SPET) is an emerging and increasingly popular solution for high-throughput targeted genotyping in plants. Although SPET requires a priori identification of polymorphisms for probe design, this technology has potentially higher reproducibility and transferability compared to other reduced representation sequencing (RRS) approaches, also enabling the discovery of closely linked polymorphisms surrounding the target one.

Results

The potential for SPET application in fruit trees was evaluated by developing a 25K target SNPs assay to genotype a panel of apricot accessions and progenies. A total of 32,492 polymorphic sites were genotyped in 128 accessions (including 8,188 accessory non-target SNPs) with extremely low levels of missing data and a significant correlation of allelic frequencies compared to whole-genome sequencing data used for array design. Assay performance was further validated by estimating genotyping errors in two biparental progenies, resulting in an overall 1.8% rate. SPET genotyping data were used to infer population structure and to dissect the architecture of fruit maturity date (MD), a quantitative reproductive phenological trait of great agronomical interest in apricot species. Depending on the year, GWAS revealed loci associated to MD on several chromosomes. The QTLs on chromosomes 1 and 4 (the latter explaining most of the phenotypic variability in the panel) were the most consistent over years and were further confirmed by linkage mapping in two segregating progenies.

Conclusions

Besides the utility for marker assisted selection and for paving the way to in-depth studies to clarify the molecular bases of MD trait variation in apricot, the results provide an overview of the performance and reliability of SPET for fruit tree genetics.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08901-1.

An apple somatic mutation of delayed fruit maturation date is primarily caused by a retrotransposon insertion-associated large deletion

Somatic mutations may alter important traits in tree fruits, such as fruit color, size and maturation date. Autumn Gala (AGala), a somatic mutation from apple cultivar Gala, matures 4 weeks later than Gala. To understand the mechanisms underlying the delayed maturation, RNA-seq analyses were conducted with fruit sampled at 13 (Gala) and 16 (AGala) time-points during their growth and development. Weighted gene co-expression network analysis (WGCNA) of 23 372 differentially expressed genes resulted in 25 WGCNA modules. Of these, modules 1 (r = -0.98, P = 2E-21) and 2 (r = -0.52, P = 0.004), which were suppressed in AGala, were correlated with fruit maturation date. Surprisingly, 77 of the 152 member genes in module 1 were harbored in a 2.8-Mb genomic region on chromosome 6 that was deleted and replaced by a 10.7-kb gypsy-like retrotransposon (Gy-36) from chromosome 7 in AGala. Among the 77 member genes, MdACT7 was the most suppressed (by 10.5-fold) in AGala due to a disruptive 2.5-kb insertion in coding sequence. Moreover, MdACT7 is the exclusive apple counterpart of Arabidopsis ACT7 known of essential roles in plant development, and the functional allele MdACT7, which was lost to the deletion in AGala, was associated with early fruit maturation in 268 apple accessions. Overexpressing alleles MdACT7 and Mdact7 in an Arabidopsis act7 line showed that MdACT7 largely rescued its stunted growth and delayed initial flowering while Mdact7 did not. Therefore, the 2.8-Mb hemizygous deletion is largely genetically causal for fruit maturation delay in AGala, and the total loss of MdACT7 might have contributed to the phenotype.

QTLs Identification for Iron Chlorosis in a Segregating Peach-Almond Progeny Through Double-Digest Sequence-Based Genotyping (SBG)

Linkage maps are highly appreciated tools for cultivar and rootstock breeding programs because they are suitable for genetic and genomic studies. In this study, we report on using sequence-based genotyping (SBG) approach to simultaneously discover and genotype SNPs from two peach-based rootstocks ("Adafuel" and "Flordaguard") and their progeny (n = 118): from a initial mapping population composed of 131 seedlings. The plant material was developed at the EEAD-CSIC Prunus rootstocks breeding program, aiming to obtain a segregating progeny for a range of characters of agronomical interest to rootstock breeding (iron-chlorosis and root-asphyxia tolerance, nematode resistance, vigor traits, and other effects on scion cultivars). Sequence reads obtained from double-digest SBG were aligned to the P. persica reference genome (Peach v2.0). While eight linkage groups were constructed for "Adafuel," only four linkage groups were constructed for "Flordaguard," given the low heterozygosity of this last genotype. High synteny and co-linearity were observed between obtained maps and Peach v2.0. On the other hand, this work aimed to elucidate the genetic basis of leaf chlorosis tolerance using the phenotypic segregation of the progeny to iron-chlorosis tolerance, along with the QTLs responsible for leaf chlorosis. The F1 mapping population, composed initially of 131 seedlings, was growing in four field trials established on calcareous soils at the experimental field of the EEAD-CSIC in Zaragoza, Spain. From the initial mapping population, 131 individuals were selected for their phenotypical characterization with SPAD measurements of plants grown in the field, exhibiting a great variability. Significant QTLs associated with tolerance to iron chlorosis were found in LG1, LG5, LG7, and LG8. The significant QTLs detected in LG5 and LG7 have not been associated with this abiotic stress before in Prunus. Several candidate genes such as Prupe.1G541100, predicted as glutamyl-tRNA reductase 1, Prupe.1G468200, encoding a 2-oxoglutarate (2OG), and Fe(II)-dependent oxygenase superfamily protein or Prupe.1G577000 (ppa011050.m), a NIFU-like protein 2 (NIFU2) were detected. The exact biological function of some of these genes should be verified for the future development of marker-assisted selection for peach iron chlorosis tolerance.

De novo chromosome-level genome of a semi-dwarf cultivar of Prunus persica identifies the aquaporin PpTIP2 as responsible for temperature-sensitive semi-dwarf trait and PpB3-1 for flower type and size

Peach (Prunus persica) is one of the most important fruit crops globally, but its cultivation can be hindered by large tree size. 'Zhongyoutao 14' (CN14) is a temperature-sensitive semi-dwarf (TSSD) cultivar which might be useful as breeding stock. The genome of CN14 was sequenced and assembled de novo using single-molecule real-time sequencing and chromosome conformation capture assembly. A high-quality genome was assembled and annotated, with 228.82 Mb mapped to eight chromosomes. Eighty-six re-sequenced F₁ individuals and 334 previously re-sequenced accessions were used to identify candidate genes controlling TSSD and flower type and size. An aquaporin tonoplast intrinsic protein (PpTIP2) was a strong candidate gene for control of TSSD. Sequence variations in the upstream regulatory region of PpTIP2 correlated with different transcriptional activity at different temperatures. PpB3-1, a candidate gene for flower type (SH) and flower size, contributed to petal development and promoted petal enlargement. The locus of another 12 agronomic traits was identified through genome-wide association study. Most of these loci exhibited consistent and precise association signals, except for flesh texture and flesh adhesion. A 6015-bp insertion in exon 3 and a 26-bp insertion upstream of PpMYB25 were associated with fruit hairless. Along with a 70.5-Kb gap at the F-M locus in CN14, another two new alleles were identified in peach accessions. Our findings will not only promote genomic research and agronomic breeding in peach but also provide a foundation for the peach pan-genome.

De novo chromosome-level genome of a semi-dwarf cultivar of Prunus persica identifies the aquaporin PpTIP2 as responsible for temperature-sensitive semi-dwarf trait and PpB3-1 for flower type and size

Peach (Prunus persica) is one of the most important fruit crops globally, but its cultivation can be hindered by large tree size. 'Zhongyoutao 14' (CN14) is a temperature-sensitive semi-dwarf (TSSD) cultivar which might be useful as breeding stock. The genome of CN14 was sequenced and assembled de novo using single-molecule real-time sequencing and chromosome conformation capture assembly. A high-quality genome was assembled and annotated, with 228.82 Mb mapped to eight chromosomes. Eighty-six re-sequenced F₁ individuals and 334 previously re-sequenced accessions were used to identify candidate genes controlling TSSD and flower type and size. An aquaporin tonoplast intrinsic protein (PpTIP2) was a strong candidate gene for control of TSSD. Sequence variations in the upstream regulatory region of PpTIP2 correlated with different transcriptional activity at different temperatures. PpB3-1, a candidate gene for flower type (SH) and flower size, contributed to petal development and promoted petal enlargement. The locus of another 12 agronomic traits was identified through genome-wide association study. Most of these loci exhibited consistent and precise association signals, except for flesh texture and flesh adhesion. A 6015-bp insertion in exon 3 and a 26-bp insertion upstream of PpMYB25 were associated with fruit hairless. Along with a 70.5-Kb gap at the F-M locus in CN14, another two new alleles were identified in peach accessions. Our findings will not only promote genomic research and agronomic breeding in peach but also provide a foundation for the peach pan-genome.

Monitoring Apricot (Prunus armeniaca L.) Ripening Progression through Candidate Gene Expression Analysis

This study aimed at the monitoring of the apricot (Prunus armeniaca L.) ripening progression through the expression analysis of 25 genes related to fruit quality traits in nine cultivars with great differences in fruit color and ripening date. The level of pigment compounds, such as anthocyanins and carotenoids, is a key factor in food taste, and is responsible for the reddish blush color or orange skin and flesh color in apricot fruit, which are desirable quality traits in apricot breeding programs. The construction of multiple linear regression models to predict anthocyanins and carotenoids content from gene expression allows us to evaluate which genes have the strongest influence over fruit color, as these candidate genes are key during biosynthetic pathways or gene expression regulation, and are responsible for the final fruit phenotype. We propose the gene CHS as the main predictor for anthocyanins content, CCD4 and ZDS for carotenoids content, and LOX2 and MADS-box for the beginning and end of the ripening process in apricot fruit. All these genes could be applied as RNA markers to monitoring the ripening stage and estimate the anthocyanins and carotenoids content in apricot fruit during the ripening process.

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.

An Upgraded, Highly Saturated Linkage Map of Japanese Plum (Prunus salicina Lindl.), and Identification of a New Major Locus Controlling the Flavan-3-ol Composition in Fruits

Japanese plum fruits are rich in phenolic compounds, such as anthocyanins and flavan-3-ols, whose contents vary significantly among cultivars. Catechin (C) and epicatechin (EC) are flavan-3-ol monomers described in the fruits of this species and are associated with bitterness, astringency, antioxidant capacity, and susceptibility to enzymatic mesocarp browning. In this study, we aimed to identify quantitative trait loci (QTL) associated with the content of flavan-3-ol in Japanese plum fruits. We evaluated the content of C and EC in the mesocarp and exocarp of samples from 79 and 64 seedlings of an F1 progeny (<'98-99' × 'Angeleno'>) in the first and second seasons, respectively. We also constructed improved versions of linkage maps from '98-99' and 'Angeleno,' presently called single-nucleotide polymorphisms (SNPs) after mapping the already available GBS reads to Prunus salicina Lindl. cv. 'Sanyueli' v2.0 reference genome. These data allowed for describing a cluster of QTLs in the cultivar, 'Angeleno,' associated with the flavan-3-ol composition of mesocarp and exocarp, which explain up to 100% of the C/EC ratio. Additionally, we developed a C/EC metabolic marker, which was mapped between the markers with the highest log of odds (LOD) scores detected by the QTL analysis. The C/EC locus was located in the LG1, at an interval spanning 0.70 cM at 108.30-108.90 cM. Our results suggest the presence of a novel major gene controlling the preferential synthesis of C or EC in the Japanese plum fruits. This study is a significant advance in understanding the regulation of synthesizing compounds associated with fruit quality, postharvest, and human health promotion.

New insights into flowering date in Prunus: fine mapping of a major QTL in sweet cherry

Flowering date is an important trait in Prunus fruit species, especially for their adaptation in a global warming context. Numerous quantitative trait loci (QTLs) have been identified and a major one was previously located on LG4. The objectives of this study were to fine-map this QTL in sweet cherry, to identify robust candidate genes by using the new sweet cherry genome sequence of the cultivar 'Regina' and to define markers usable in marker-assisted selection (MAS). We performed QTL analyses on two populations derived from crosses using cultivars 'Regina' and 'Garnet' as parents. The first one (n = 117) was phenotyped over ten years, while the second one (n = 1386) was evaluated during three years. Kompetitive allele specific PCR (KASP) markers located within the QTL region on LG4 were developed and mapped within this region, consisting in the first fine mapping in sweet cherry. The QTL interval was narrowed from 380 kb to 68 kb and candidate genes were identified by using the genome sequence of 'Regina'. Their expression was analyzed from bud dormancy period to flowering in cultivars 'Regina' and 'Garnet'. Several genes, such as PavBOI-E3, PavSR45a and PavSAUR71, were differentially expressed in these two cultivars and could be then considered as promising candidate genes. Two KASP markers were validated using a population derived from a cross between cultivars 'Regina' and 'Lapins' and two collections, including landraces and modern cultivars. Thanks to the high synteny within the Prunus genus, these results give new insights into the control of flowering date in Prunus species and pave the way for the development of molecular breeding strategies.

Synergistic Effects of Plant Hormones on Spontaneous Late-ripening Mutant of ‘Jinghong’ Peach Detected by Transcriptome Analysis

Objectives

Peach (Prunus persica L.) is an ancient fruit tree that originated from China. It is the climacteric fruit which belongs to genus Prunus in family Rosaceae. Ethylene, which is produced during ripening, accelerates fruit softening, and therefore peaches cannot be stored for a long time.

Materials and Methods

To study the mechanism of fruit late ripening, transcriptome analysis of the fruit of a late-ripening mutant of ‘Jinghong’ peach was performed to identify genes and pathways involved in fruit late ripening.

Results

A total of 1,805, 1,511, and 2,309 genes were found to be differentially expressed in W2_vs_M1, W3_vs_M2, and W3_vs_M3, respectively. Functional enrichment analysis of the differentially expressed genes showed they were related to carotenoid biosynthesis, starch and sucrose metabolism plant hormone signal transduction, flavonoid biosynthesis, and photosynthesis. The expression trends of ripening-related genes that encode transcription factors and plant hormone signal transduction-related genes that encode enzymes were similar.

Conclusions

It will help to elucidate the transcriptional regulatory network of fruit development in the spontaneous late-ripening mutant of ‘Jinghong’ peach and provide a theoretical basis for understanding the molecular regulatory mechanism of fruit ripening.

NAC Transcription Factor Family Regulation of Fruit Ripening and Quality: A Review

The NAC transcription factor (TF) family is one of the largest plant-specific TF families and its members are involved in the regulation of many vital biological processes during plant growth and development. Recent studies have found that NAC TFs play important roles during the ripening of fleshy fruits and the development of quality attributes. This review focuses on the advances in our understanding of the function of NAC TFs in different fruits and their involvement in the biosynthesis and signal transduction of plant hormones, fruit textural changes, color transformation, accumulation of flavor compounds, seed development and fruit senescence. We discuss the theoretical basis and potential regulatory models for NAC TFs action and provide a comprehensive view of their multiple roles in modulating different aspects of fruit ripening and quality.

Knock-Out of CmNAC-NOR Affects Melon Climacteric Fruit Ripening

Fruit ripening is an important process that affects fruit quality. A QTL in melon, ETHQV6.3, involved in climacteric ripening regulation, has been found to be encoded by CmNAC-NOR, a homologue of the tomato NOR gene. To further investigate CmNAC-NOR function, we obtained two CRISPR/Cas9-mediated mutants (nor-3 and nor-1) in the climacteric Védrantais background. nor-3, containing a 3-bp deletion altering the NAC domain A, resulted in ~8 days delay in ripening without affecting fruit quality. In contrast, the 1-bp deletion in nor-1 resulted in a fully disrupted NAC domain, which completely blocked climacteric ripening. The nor-1 fruits did not produce ethylene, no abscission layer was formed and there was no external color change. Additionally, volatile components were dramatically altered, seeds were not well developed and flesh firmness was also altered. There was a delay in fruit ripening with the nor-1 allele in heterozygosis of ~20 days. Our results provide new information regarding the function of CmNAC-NOR in melon fruit ripening, suggesting that it is a potential target for modulating shelf life in commercial climacteric melon varieties.

Unravelling the Molecular Regulation Mechanisms of Slow Ripening Trait in Prunus persica

Fruit development is a complex process that involves the interplay of cell division, expansion, and differentiation. As a model to study fruit development, nectarines incapable of ripening were described as slow ripening. Slow ripening fruits remained firm and exhibited no rise in CO₂ or ethylene production rates for one month or more at 20 °C. Different studies suggest that this trait is controlled by a single gene (NAC072). Transcriptome analysis between normal and slow ripening fruits showed a total of 157, 269, 976, and 5.224 differentially expressed genes in each fruit developmental stage analyzed (T1, T2, T3, and T7, respectively), and no expression of NAC072 was found in the slow ripening individuals. Using this transcriptomic information, we identified a correlation of NAC072 with auxin-related genes and two genes associated with terpene biosynthesis. On the other hand, significant differences were observed in hormonal biosynthetic pathways during fruit development between the normal and slow ripening individuals (gibberellin, ethylene, jasmonic acid and abscisic acid). These results suggest that the absence of NAC072 by the direct or indirect expression control of auxins or terpene-related genes prevents normal peach fruit development.

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).

Molecular Insights of Fruit Quality Traits in Peaches, Prunus persica

Fleshy fruits are the most demanded fruits because of their organoleptic qualities and nutritional values. The genus Prunus is a rich source of diversified stone/drupe fruits such as almonds, apricots, plums, sweet cherries, peaches, and nectarines. The fruit-ripening process in Prunus involves coordinated biochemical and physiological changes resulting in changes in fruit texture, aroma gain, color change in the pericarp, sugar/organic acid balance, fruit growth, and weight gain. There are different varieties of peaches with unique palatable qualities and gaining knowledge in the genetics behind these quality traits helps in seedling selection for breeding programs. In addition, peaches have shorter post-harvest life due to excessive softening, resulting in fruit quality reduction and market loss. Many studies have been executed to understand the softening process at the molecular level to find the genetic basis. To summarize, this review focused on the molecular aspects of peach fruit quality attributes and their related genetics to understand the underlying mechanisms.

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.

Genome-wide association study of individual sugar content in fruit of Japanese pear (Pyrus spp.)

BACKGROUND

Understanding mechanisms of sugar accumulation and composition is essential to determining fruit quality and maintaining a desirable balance of sugars in plant storage organs. The major sugars in mature Rosaceae fruits are sucrose, fructose, glucose, and sorbitol. Among these, sucrose and fructose have high sweetness, whereas glucose and sorbitol have low sweetness. Japanese pear has extensive variation in individual sugar contents in mature fruit. Increasing total sugar content and that of individual high-sweetness sugars is a major target of breeding programs. The objective of this study was to identify quantitative trait loci (QTLs) associated with fruit traits including individual sugar accumulation, to infer the candidate genes underlying the QTLs, and to assess the potential of genomic selection for breeding pear fruit traits.

RESULTS

We evaluated 10 fruit traits and conducted genome-wide association studies (GWAS) for 106 cultivars and 17 breeding populations (1112 F1 individuals) using 3484 tag single-nucleotide polymorphisms (SNPs). By implementing a mixed linear model and a Bayesian multiple-QTL model in GWAS, 56 SNPs associated with fruit traits were identified. In particular, a SNP located close to acid invertase gene PPAIV3 on chromosome 7 and a newly identified SNP on chromosome 11 had quite large effects on accumulation of sucrose and glucose, respectively. We used 'Golden Delicious' doubled haploid 13 (GDDH13), an apple reference genome, to infer the candidate genes for the identified SNPs. In the region flanking the SNP on chromosome 11, there is a tandem repeat of early responsive to dehydration (ERD6)-like sugar transporter genes that might play a role in the phenotypes observed.

CONCLUSIONS

SNPs associated with individual sugar accumulation were newly identified at several loci, and candidate genes underlying QTLs were inferred using advanced apple genome information. The candidate genes for the QTLs are conserved across Pyrinae genomes, which will be useful for further fruit quality studies in Rosaceae. The accuracies of genomic selection for sucrose, fructose, and glucose with genomic best linear unbiased prediction (GBLUP) were relatively high (0.67-0.75), suggesting that it would be possible to select individuals having high-sweetness fruit with high sucrose and fructose contents and low glucose content.

Apple Ripening Is Controlled by a NAC Transcription Factor

Softening is a hallmark of ripening in fleshy fruits, and has both desirable and undesirable implications for texture and postharvest stability. Accordingly, the timing and extent of pre-harvest ripening and associated textural changes following harvest are key targets for improving fruit quality through breeding. Previously, we identified a large effect locus associated with harvest date and firmness in apple (Malus domestica) using genome-wide association studies (GWAS). Here, we present additional evidence that polymorphisms in or around a transcription factor gene, NAC18.1, may cause variation in these traits. First, we confirmed our previous findings with new phenotype and genotype data from ∼800 apple accessions. In this population, we compared a genetic marker within NAC18.1 to markers targeting three other firmness-related genes currently used by breeders (ACS1, ACO1, and PG1), and found that the NAC18.1 marker was the strongest predictor of both firmness at harvest and firmness after 3 months of cold storage. By sequencing NAC18.1 across 18 accessions, we revealed two predominant haplotypes containing the single nucleotide polymorphism (SNP) previously identified using GWAS, as well as dozens of additional SNPs and indels in both the coding and promoter sequences. NAC18.1 encodes a protein that is orthogolous to the NON-RIPENING (NOR) transcription factor, a regulator of ripening in tomato (Solanum lycopersicum). We introduced both NAC18.1 transgene haplotypes into the tomato nor mutant and showed that both haplotypes complement the nor ripening deficiency. Taken together, these results indicate that polymorphisms in NAC18.1 may underlie substantial variation in apple firmness through modulation of a conserved ripening program.

Mapping and characterization QTLs for phenological traits in seven pedigree-connected peach families

Background

Environmental adaptation and expanding harvest seasons are primary goals of most peach [Prunus persica (L.) Batsch] breeding programs. Breeding perennial crops is a challenging task due to their long breeding cycles and large tree size. Pedigree-based analysis using pedigreed families followed by haplotype construction creates a platform for QTL and marker identification, validation, and the use of marker-assisted selection in breeding programs.

Results

Phenotypic data of seven F₁ low to medium chill full-sib families were collected over 2 years at two locations and genotyped using the 9 K SNP Illumina array. Three QTLs were discovered for bloom date (BD) and mapped on linkage group 1 (LG1) (172–182 cM), LG4 (48–54 cM), and LG7 (62–70 cM), explaining 17–54%, 11–55%, and 11–18% of the phenotypic variance, respectively. The QTL for ripening date (RD) and fruit development period (FDP) on LG4 was co-localized at the central part of LG4 (40–46 cM) and explained between 40 and 75% of the phenotypic variance. Haplotype analyses revealed SNP haplotypes and predictive SNP marker(s) associated with desired QTL alleles and the presence of multiple functional alleles with different effects for a single locus for RD and FDP.

Conclusions

A multiple pedigree-linked families approach validated major QTLs for the three key phenological traits which were reported in previous studies across diverse materials, geographical distributions, and QTL mapping methods. Haplotype characterization of these genomic regions differentiates this study from the previous QTL studies. Our results will provide the peach breeder with the haplotypes for three BD QTLs and one RD/FDP QTL to create predictive DNA-based molecular marker tests to select parents and/or seedlings that have desired QTL alleles and cull unwanted genotypes in early seedling stages.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07483-8.

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.

Molecular Bases of Fruit Quality in Prunus Species: An Integrated Genomic, Transcriptomic, and Metabolic Review with a Breeding Perspective

In plants, fruit ripening is a coordinated developmental process that requires the change in expression of hundreds to thousands of genes to modify many biochemical and physiological signal cascades such as carbohydrate and organic acid metabolism, cell wall restructuring, ethylene production, stress response, and organoleptic compound formation. In Prunus species (including peaches, apricots, plums, and cherries), fruit ripening leads to the breakdown of complex carbohydrates into sugars, fruit firmness reductions (softening by cell wall degradation and cuticle properties alteration), color changes (loss of green color by chlorophylls degradation and increase in non-photosynthetic pigments like anthocyanins and carotenoids), acidity decreases, and aroma increases (the production and release of organic volatile compounds). Actually, the level of information of molecular events at the transcriptional, biochemical, hormonal, and metabolite levels underlying ripening in Prunus fruits has increased considerably. However, we still poorly understand the molecular switch that occurs during the transition from unripe to ripe fruits. The objective of this review was to analyze of the molecular bases of fruit quality in Prunus species through an integrated metabolic, genomic, transcriptomic, and epigenetic approach to better understand the molecular switch involved in the ripening process with important consequences from a breeding point of view.

An integrated peach genome structural variation map uncovers genes associated with fruit traits

BACKGROUND

Genome structural variations (SVs) have been associated with key traits in a wide range of agronomically important species; however, SV profiles of peach and their functional impacts remain largely unexplored.

RESULTS

Here, we present an integrated map of 202,273 SVs from 336 peach genomes. A substantial number of SVs have been selected during peach domestication and improvement, which together affect 2268 genes. Genome-wide association studies of 26 agronomic traits using these SVs identify a number of candidate causal variants. A 9-bp insertion in Prupe.4G186800, which encodes a NAC transcription factor, is shown to be associated with early fruit maturity, and a 487-bp deletion in the promoter of PpMYB10.1 is associated with flesh color around the stone. In addition, a 1.67 Mb inversion is highly associated with fruit shape, and a gene adjacent to the inversion breakpoint, PpOFP1, regulates flat shape formation.

CONCLUSIONS

The integrated peach SV map and the identified candidate genes and variants represent valuable resources for future genomic research and breeding in peach.

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.

A draft genome of sweet cherry (Prunus avium L.) reveals genome-wide and local effects of domestication

Sweet cherry (Prunus avium L.) trees are both economically important fruit crops but also important components of natural forest ecosystems in Europe, Asia and Africa. Wild and domesticated trees currently coexist in the same geographic areas with important questions arising on their historical relationships. Little is known about the effects of the domestication process on the evolution of the sweet cherry genome. We assembled and annotated the genome of the cultivated variety "Big Star*" and assessed the genetic diversity among 97 sweet cherry accessions representing three different stages in the domestication and breeding process (wild trees, landraces and modern varieties). The genetic diversity analysis revealed significant genome-wide losses of variation among the three stages and supports a clear distinction between wild and domesticated trees, with only limited gene flow being detected between wild trees and domesticated landraces. We identified 11 domestication sweeps and five breeding sweeps covering, respectively, 11.0 and 2.4 Mb of the P. avium genome. A considerable fraction of the domestication sweeps overlaps with those detected in the related species, Prunus persica (peach), indicating that artificial selection during domestication may have acted independently on the same regions and genes in the two species. We detected 104 candidate genes in sweep regions involved in different processes, such as the determination of fruit texture, the regulation of flowering and fruit ripening and the resistance to pathogens. The signatures of selection identified will enable future evolutionary studies and provide a valuable resource for genetic improvement and conservation programs in sweet cherry.

Multiple-population QTL mapping of maturity and fruit-quality traits reveals LG4 region as a breeding target in sweet cherry (Prunus avium L.)

Sweet cherry maturity date and fruit quality are relevant traits for its marketability, transport, and consumer acceptance. In this work, sweet cherry fruit development time, maturity date, and commercial fruit-quality traits (size, weight, firmness, soluble solid content, and titratable acidity) were investigated to improve the knowledge of their genetic control, and to identify alleles of breeding interest. Six sweet cherry populations segregating for these traits were used for QTL analyses. These populations descend from cross- and self-pollinations of local Spanish sweet cherries 'Ambrunés' and 'Cristobalina', and breed cultivars ('Brooks', 'Lambert', or 'Vic'). The six populations (n = 411), previously genotyped with RosBREED Cherry 6 K SNP array, were phenotyped for 2 years. QTL analyses were conducted using a multifamily approach implemented by FlexQTL™. Fruit development time, soluble solid content, and titratable acidity QTLs are first reported in sweet cherry in this work. Significant QTLs were detected for all the traits. Eighteen were more stable as they were detected for 2 years. Of these, nine are first reported in this work. The major QTLs for fruit development time, maturity date, firmness, and soluble solid content were identified on the same narrow region of linkage group 4. These traits also showed significant positive correlation (long fruit development time associated with late maturity, high firmness, and high SSC). NAC transcription factor genes identified on this LG4 region may be candidate genes for the regulation of these traits in sweet cherry, as previously described in syntenic regions of other Rosaceae species. Haplotypes of breeding interest on this LG4 genomic region were identified and will be useful for sweet cherry breeding from this and related plant material.

Identification and characterization of QTLs for fruit quality traits in peach through a multi-family approach

BACKGROUND

Fruit quality traits have a significant effect on consumer acceptance and subsequently on peach (Prunus persica (L.) Batsch) consumption. Determining the genetic bases of key fruit quality traits is essential for the industry to improve fruit quality and increase consumption. Pedigree-based analysis across multiple peach pedigrees can identify the genomic basis of complex traits for direct implementation in marker-assisted selection. This strategy provides breeders with better-informed decisions and improves selection efficiency and, subsequently, saves resources and time.

RESULTS

Phenotypic data of seven F1 low to medium chill full-sib families were collected over 2 years at two locations and genotyped using the 9 K SNP Illumina array. One major QTL for fruit blush was found on linkage group 4 (LG4) at 40-46 cM that explained from 20 to 32% of the total phenotypic variance and showed three QTL alleles of different effects. For soluble solids concentration (SSC), one QTL was mapped on LG5 at 60-72 cM and explained from 17 to 39% of the phenotypic variance. A major QTL for titratable acidity (TA) co-localized with the major locus for low-acid fruit (D-locus). It was mapped at the proximal end of LG5 and explained 35 to 80% of the phenotypic variance. The new QTL for TA on the distal end of LG5 explained 14 to 22% of the phenotypic variance. This QTL co-localized with the QTL for SSC and affected TA only when the first QTL is homozygous for high acidity (epistasis). Haplotype analyses revealed SNP haplotypes and predictive SNP marker(s) associated with desired QTL alleles.

CONCLUSIONS

A multi-family-based QTL discovery approach enhanced the ability to discover a new TA QTL at the distal end of LG5 and validated other QTLs which were reported in previous studies. Haplotype characterization of the mapped QTLs distinguishes this work from the previous QTL studies. Identified predictive SNPs and their original sources will facilitate the selection of parents and/or seedlings that have desired QTL alleles. Our findings will help peach breeders develop new predictive, DNA-based molecular marker tests for routine use in marker-assisted breeding.

RNA-sequencing analysis revealed genes associated drought stress responses of different durations in hexaploid sweet potato

Purple-fleshed sweet potato (PFSP) is an important food crop, as it is a rich source of nutrients and anthocyanin pigments. Drought has become a major threat to sustainable sweetpotato production, resulting in huge yield losses. Therefore, the present study was conducted to identify drought stress-responsive genes using next-generation (NGS) and third-generation sequencing (TGS) techniques. Five cDNA libraries were constructed from seedling leaf segments treated with a 30% solution of polyethylene glycol (PEG-6000) for 0, 1, 6, 12, and 48 h for second-generation sequencing. Leaf samples taken from upper third of sweet potato seedlings after 1, 6, 12, and 48 h of drought stress were used for the construction of cDNA libraries for third-generation sequencing; however, leaf samples from untreated plants were collected as controls. A total of 184,259,679 clean reads were obtained using second and third-generation sequencing and then assembled into 17,508 unigenes with an average length of 1,783 base pairs. Out of 17,508 unigenes, 642 (3.6%) unigenes failed to hit any homologs in any databases, which might be considered novel genes. A total of 2, 920, 1578, and 2,418 up-regulated unigenes and 3,834, 2,131, and 3,337 down-regulated unigenes from 1 h, 6 h, 12 h, and 48 h library were identified, respectively in drought stress versus control. In addition, after 6, 12, and 48 h of drought stress, 540 up-regulated unigenes, 486 down-regulated unigenes and 414 significantly differentially expressed unigenes were detected. It was found that several gene families including Basic Helix-loop-helix (bHLH), basic leucine zipper (bZIP), Cystein2/Histidine2 (C2H2), C3H, Ethylene-responsive transcription factor (ERF), Homo domain-leucine zipper (HD-ZIP), MYB, NAC (NAM, ATAF1/2, and CUC2), Thiol specific antioxidant and WRKY showed responses to drought stress. In total, 17,472 simple sequence repeats and 510,617 single nucleotide polymorphisms were identified based on transcriptome sequencing of the PFSP. About 96.55% of the obtained sequences are not available online in sweet potato genomics resources. Therefore, it will enrich annotated sweet potato gene sequences and enhance understanding of the mechanisms of drought tolerance through genetic manipulation. Moreover, it represents a sequence resource for genetic and genomic studies of sweet potato.

Detection and application of genome-wide variations in peach for association and genetic relationship analysis

BACKGROUND

Peach (Prunus persica L.) is a diploid species and model plant of the Rosaceae family. In the past decade, significant progress has been made in peach genetic research via DNA markers, but the number of these markers remains limited.

RESULTS

In this study, we performed a genome-wide DNA markers detection based on sequencing data of six distantly related peach accessions. A total of 650,693~1,053,547 single nucleotide polymorphisms (SNPs), 114,227~178,968 small insertion/deletions (InDels), 8386~12,298 structure variants (SVs), 2111~2581 copy number variants (CNVs) and 229,357~346,940 simple sequence repeats (SSRs) were detected and annotated. To demonstrate the application of DNA markers, 944 SNPs were filtered for association study of fruit ripening time and 15 highly polymorphic SSRs were selected to analyze the genetic relationship among 221 accessions.

CONCLUSIONS

The results showed that the use of high-throughput sequencing to develop DNA markers is fast and effective. Comprehensive identification of DNA markers, including SVs and SSRs, would be of benefit to genetic diversity evaluation, genetic mapping, and molecular breeding of peach.

Refining the Genomic Region Containing a Major Locus Controlling Fruit Maturity in Peach

Maturity date (MD), defined as the duration between the first calendar day of the year and maturity, and fruit development period (FDP), defined as the duration between full bloom and maturity, are highly variable in peach [Prunus persica (L.) Batsch]. There is a need to discover molecular markers associated with these traits in order to enhance the efficiency and reliability of breeding for extending the harvest season in peach. An association mapping population consisting of 132 peach accessions was phenotypically evaluated for MD and FDP, and genotypically characterized using the genotyping-by-sequencing (GBS) approach. The phenotypic and genotypic data collected were used to conduct a genome-wide association study (GWAS). The GWAS identified three SNPs on chromosome 4 that are significantly associated with both FDP and MD. These three SNPs covered a region of 43,067 bp; we referred to this region as the MD/FDP locus. Seven genes were identified in the MD/FDP locus. One or more of these genes is believed to regulate some aspect of maturity in peach. The data reported here is expected to aid in marker-assisted seedling selection (MASS) targeted towards widening peach germplasm for maturity, particularly early maturity.

A fruit firmness QTL identified on linkage group 4 in sweet cherry (Prunus avium L.) is associated with domesticated and bred germplasm

Fruit firmness is an important market driven trait in sweet cherry (Prunus avium L.) where the desirable increase in fruit firmness is associated with landrace and bred cultivars. The aim of this work was to investigate the genetic basis of fruit firmness using plant materials that include wild cherry (syn. mazzard), landrace and bred sweet cherry germplasm. A major QTL for fruit firmness, named qP-FF4.1, that had not previously been reported, was identified in three sweet cherry populations. Thirteen haplotypes (alleles) associated with either soft or firm fruit were identified for qP-FF4.1 in the sweet cherry germplasm, and the “soft” alleles were dominant over the “firm” alleles. The finding that sweet cherry individuals that are homozygous for the “soft” alleles for qP-FF4.1 are exclusively mazzards and that the vast majority of the bred cultivars are homozygous for “firm” alleles suggests that this locus is a signature of selection. Candidate genes related to plant cell wall modification and various plant hormone signaling pathways were identified, with an expansin gene being the most promising candidate. These results advance our understanding of the genetic basis of fruit firmness and will help to enable the use of DNA informed breeding for this trait in sweet cherry breeding programs.

Expression QTL (eQTLs) Analyses Reveal Candidate Genes Associated With Fruit Flesh Softening Rate in Peach [Prunus persica (L.) Batsch]

Significant differences in softening rate have been reported between melting flesh in peach and nectarine varieties. This trait seems to be controlled by several genes. We aimed to identify candidate genes involved in fruit softening rate by integrating quantitative trait loci (QTL) and expression QTL (eQTL) analyses, comparing siblings with contrasting softening rates. We used a segregating population derived from nectarine cv. 'Venus' selfing, which was phenotyped for softening rate during three seasons. Six siblings with high (HSR) and six with low softening rate (LSR) were sequenced using RNA-Seq. A group of 5,041 differentially expressed genes was identified. Also, we found a QTL with a LOD (logarithm of odds) score of 9.7 on LG4 in all analyzed seasons. Furthermore, we detected 1,062 eQTLs, of which 133 were found co-localizing with the identified QTL. Gene Ontology (GO) analysis showed 'Response to auxin' as one the main over-represented categories. Our findings suggest over-expression of auxin biosynthetic related genes in the HSR group, which implies a higher expression and/or accumulation of auxin, thereby triggering fast softening. Conversely, the LSR phenotype might be explained by an altered auxin-homeostasis associated with low auxin levels. This work will contribute to unraveling the genetic mechanisms responsible for the softening rate in peaches and nectarines and lead to the development of molecular markers.

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.

High-density multi-population consensus genetic linkage map for peach

Highly saturated genetic linkage maps are extremely helpful to breeders and are an essential prerequisite for many biological applications such as the identification of marker-trait associations, mapping quantitative trait loci (QTL), candidate gene identification, development of molecular markers for marker-assisted selection (MAS) and comparative genetic studies. Several high-density genetic maps, constructed using the 9K SNP peach array, are available for peach. However, each of these maps is based on a single mapping population and has limited use for QTL discovery and comparative studies. A consensus genetic linkage map developed from multiple populations provides not only a higher marker density and a greater genome coverage when compared to the individual maps, but also serves as a valuable tool for estimating genetic positions of unmapped markers. In this study, a previously developed linkage map from the cross between two peach cultivars 'Zin Dai' and 'Crimson Lady' (ZC2) was improved by genotyping additional progenies. In addition, a peach consensus map was developed based on the combination of the improved ZC2 genetic linkage map with three existing high-density genetic maps of peach and a reference map of Prunus. A total of 1,476 SNPs representing 351 unique marker positions were mapped across eight linkage groups on the ZC2 genetic map. The ZC2 linkage map spans 483.3 cM with an average distance between markers of 1.38 cM/marker. The MergeMap and LPmerge tools were used for the construction of a consensus map based on markers shared across five genetic linkage maps. The consensus linkage map contains a total of 3,092 molecular markers, consisting of 2,975 SNPs, 116 SSRs and 1 morphological marker associated with slow ripening in peach (SR). The consensus map provides valuable information on marker order and genetic position for QTL identification in peach and other genetic studies within Prunus and Rosaceae.

SlAREB1 transcriptional activation of NOR is involved in abscisic acid-modulated ethylene biosynthesis during tomato fruit ripening

Many studies have shown that abscisic acid (ABA) regulates climacteric fruits ripening by inducing ethylene production. Nevertheless, the key components involved in the crosstalk between these two phytohormones in controlling fruit ripening remain unknown. SlAREB1, a downstream transcription factor in ABA signaling pathway, has been reported to mediate ABA signaling that regulates tomato ripening through induction of ethylene biosynthetic genes. NOR, a member of NAC domain family, was proved to act upstream of ethylene and essential for ripening- and ethylene-associated genes expression. Here, we found that the expression of SlAREB1 and NOR are both ABA-inducible, and SlAREB1 transcription reaches the peak level prior to NOR during the ripening process. Yeast one-hybrid (Y1H), electrophoretic mobility shift assay (EMSA) and dual luciferase assay indicated NOR as a novel direct target of SlAREB1. Transient over-expression of SlAREB1 in tomato fruits results in elevated expression of NOR as well as a number of downstream ethylene biosynthetic genes including SlACS2, SlACS4 and SlACO1, suggesting that SlAREB1 can mediate ABA signal to activate NOR transcription and ultimately promote ethylene synthesis. Based on these data, we present a model suggesting that the SlAREB1-NOR regulation is a crucial node modulating ABA-regulated ethylene biosynthesis during tomato fruit ripening.

Genome encode analyses reveal the basis of convergent evolution of fleshy fruit ripening

Fleshy fruits using ethylene to regulate ripening have developed multiple times in the history of angiosperms, presenting a clear case of convergent evolution whose molecular basis remains largely unknown. Analysis of the fruitENCODE data consisting of 361 transcriptome, 71 accessible chromatin, 147 histone and 45 DNA methylation profiles reveals three types of transcriptional feedback circuits controlling ethylene-dependent fruit ripening. These circuits are evolved from senescence or floral organ identity pathways in the ancestral angiosperms either by neofunctionalisation or repurposing pre-existing genes. The epigenome, H3K27me3 in particular, has played a conserved role in restricting ripening genes and their orthologues in dry and ethylene-independent fleshy fruits. Our findings suggest that evolution of ripening is constrained by limited hormone molecules and genetic and epigenetic materials, and whole-genome duplications have provided opportunities for plants to successfully circumvent these limitations.

ETHQV6.3 is involved in melon climacteric fruit ripening and is encoded by a NAC domain transcription factor

Fruit ripening is divided into climacteric and non-climacteric types depending on the presence or absence of a transient rise in respiration rate and the production of autocatalytic ethylene. Melon is ideal for the study of fruit ripening, as both climacteric and non-climacteric varieties exist. Two introgressions of the non-climacteric accession PI 161375, encompassed in the QTLs ETHQB3.5 and ETHQV6.3, into the non-climacteric 'Piel de Sapo' background are able to induce climacteric ripening independently. We report that the gene underlying ETHQV6.3 is MELO3C016540 (CmNAC-NOR), encoding a NAC (NAM, ATAF1,2, CUC2) transcription factor that is closely related to the tomato NOR (non-ripening) gene. CmNAC-NOR was functionally validated through the identification of two TILLING lines carrying non-synonymous mutations in the conserved NAC domain region. In an otherwise highly climacteric genetic background, both mutations provoked a significant delay in the onset of fruit ripening and in the biosynthesis of ethylene. The PI 161375 allele of ETHQV6.3 is similar to that of climacteric lines of the cantalupensis type and, when introgressed into the non-climacteric 'Piel de Sapo', partially restores its climacteric ripening capacity. CmNAC-NOR is expressed in fruit flesh of both climacteric and non-climacteric lines, suggesting that the causal mutation may not be acting at the transcriptional level. The use of a comparative genetic approach in a species with both climacteric and non-climacteric ripening is a powerful strategy to dissect the complex mechanisms regulating the onset of fruit ripening.

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.

The Peach v2.0 release: high-resolution linkage mapping and deep resequencing improve chromosome-scale assembly and contiguity

BACKGROUND

The availability of the peach genome sequence has fostered relevant research in peach and related Prunus species enabling the identification of genes underlying important horticultural traits as well as the development of advanced tools for genetic and genomic analyses. The first release of the peach genome (Peach v1.0) represented a high-quality WGS (Whole Genome Shotgun) chromosome-scale assembly with high contiguity (contig L50 214.2 kb), large portions of mapped sequences (96%) and high base accuracy (99.96%). The aim of this work was to improve the quality of the first assembly by increasing the portion of mapped and oriented sequences, correcting misassemblies and improving the contiguity and base accuracy using high-throughput linkage mapping and deep resequencing approaches.

RESULTS

Four linkage maps with 3,576 molecular markers were used to improve the portion of mapped and oriented sequences (from 96.0% and 85.6% of Peach v1.0 to 99.2% and 98.2% of v2.0, respectively) and enabled a more detailed identification of discernible misassemblies (10.4 Mb in total). The deep resequencing approach fixed 859 homozygous SNPs (Single Nucleotide Polymorphisms) and 1347 homozygous indels. Moreover, the assembled NGS contigs enabled the closing of 212 gaps with an improvement in the contig L50 of 19.2%.

CONCLUSIONS

The improved high quality peach genome assembly (Peach v2.0) represents a valuable tool for the analysis of the genetic diversity, domestication, and as a vehicle for genetic improvement of peach and related Prunus species. Moreover, the important phylogenetic position of peach and the absence of recent whole genome duplication (WGD) events make peach a pivotal species for comparative genomics studies aiming at elucidating plant speciation and diversification processes.

The Peach v2.0 release: high-resolution linkage mapping and deep resequencing improve chromosome-scale assembly and contiguity

BACKGROUND

The availability of the peach genome sequence has fostered relevant research in peach and related Prunus species enabling the identification of genes underlying important horticultural traits as well as the development of advanced tools for genetic and genomic analyses. The first release of the peach genome (Peach v1.0) represented a high-quality WGS (Whole Genome Shotgun) chromosome-scale assembly with high contiguity (contig L50 214.2 kb), large portions of mapped sequences (96%) and high base accuracy (99.96%). The aim of this work was to improve the quality of the first assembly by increasing the portion of mapped and oriented sequences, correcting misassemblies and improving the contiguity and base accuracy using high-throughput linkage mapping and deep resequencing approaches.

RESULTS

Four linkage maps with 3,576 molecular markers were used to improve the portion of mapped and oriented sequences (from 96.0% and 85.6% of Peach v1.0 to 99.2% and 98.2% of v2.0, respectively) and enabled a more detailed identification of discernible misassemblies (10.4 Mb in total). The deep resequencing approach fixed 859 homozygous SNPs (Single Nucleotide Polymorphisms) and 1347 homozygous indels. Moreover, the assembled NGS contigs enabled the closing of 212 gaps with an improvement in the contig L50 of 19.2%.

CONCLUSIONS

The improved high quality peach genome assembly (Peach v2.0) represents a valuable tool for the analysis of the genetic diversity, domestication, and as a vehicle for genetic improvement of peach and related Prunus species. Moreover, the important phylogenetic position of peach and the absence of recent whole genome duplication (WGD) events make peach a pivotal species for comparative genomics studies aiming at elucidating plant speciation and diversification processes.