A gap‐free genome of pillar peach (<i>Prunus persica</i> L.) provides new insights into branch angle and double flower traits

Abstract

Peach (Prunus persica L.), a deciduous fruit tree in the Rosaceae family, is widely cultivated around the world. The release of the peach genome which significantly propelled the progress of gene mapping (Lian et al., 2022; Verde et al., 2013; Zhou et al., 2023). However, the existing peach genomes all contain multiple gaps (Table S1), which may lead to inaccurate gene annotation or gene mapping (Zhou et al., 2022). Here, we present ‘Zhaoshouhong’ peach (ZSH, also named ‘Terutebeni’, pillar type and double flower) gap-free genome assembly generated by combing Nanopore ultra-long and Hi-C reads. After scaffolding with the Hi-C data, the contigs were anchored to eight chromosomes (Appendix S1). Notably, Chromosome 1 (Chr1), Chr2, Chr3, Chr4, Chr6 and Chr7 were each covered by a single contig. Two gaps remaining on Chr5, and a single gap on Chr8 were successfully filled using the raw Nanopore ultra-long reads. Finally, a complete gap-free peach nuclear genome was obtained, with a total size of 239.34 Mb and contig N50 of 29.67 Mb (Figure 1a, Figure S1, Tables S2 and S3). Additionally, the mitochondrial and chloroplast genomes were assembled using the ultra-long and Illumina sequencing reads (Figures S2 and S3). A total of 112 959 366 bps was annotated as repeat elements (Table S4). The telomeric satellites were present at the 16 chromosome ends (Table S5), eight candidate centromeres were predicated in the ZSH genome (Table S6), the high Benchmarking universal single-copy orthologs (BUSCO) value (98.88%), Long terminal repeats (LTR) assembly index (31.03), integrity (99.63%) and accuracy (QV = 53.3) indicate the genome meets a high-quality level of assembly (Table S7). The gene models were adjusted manually using IGV-GSAman (v0.6.38, https://tbtools.cowtransfer.com/s/a11146181df14f, Figure S4), 24 901 protein-coding genes were obtained, and 23 253 genes were functionally annotated (Figure S5 and Table S8). Branch angle is one of the most important agronomic traits in fruit trees. To identify the major genes influencing branch angle in ZSH peach, structural variations were identified between the ZSH and the other peach genomes (standard type, Figure 1b; Figures S6, S7; Table S9), and 9100 variations involving in 3523 genes were commonly detected between the ZSH and the other peach genomes (Figure S8; Table S10). Notably, 25 genes showed high expression levels in standard types than that in pillar types (Figure S9; Table S11). Among these 25 genes, Pp02G17890.t1 (PpTAC1), a homologous gene of OsTAC1 in rice (Ku et al., 2011), was related to the tiller angle. Interestingly, compared with other peach genomes, the PpTAC1 in ZSH exhibited an 11 bp deletion in the promoter and a 4422 bp insertion in the exon (Figure 1c). Compared with standard peach, the variants were present in the coding or promoter sequence of PpTAC1 in other 10 pillar peach cultivars (Figure 1c; Table S12). These findings showed a close relationship between variation in PpTAC1 and the branch angle in peach. For single/double flower traits, a significant peak on Chr2 and a minor peak on Chr6 were found using 334 natural peach germplasms (Figure 1d; Table S13). Two new mutations (a 5033 bp insertion and a 1210 bp insertion) that disrupt the miR172d gene were identified in ZSH using comparative genomic analysis and gene cloning. The variation in miR172d was used as a molecular marker to distinguish the single- and double-flower trait (Table S14). The presence of either the 1210 bp or 5033 bp insertions detected in 27 varieties with double flower was absent in the six single-flower accessions and in other five double-flower accessions (‘No.18’, ‘HongChuizhi’, ‘1-1-4’, ‘1-2-7’ and ‘Huayulu’) (Figure S10; Table S14). To further identify other candidate genes that contribute to the single/double flower trait, a hybrid population was produced from ‘No.18’ (double-flower) and ‘Okubo’ (single flower) (Table S15). A major locus on Chr6 was identified using the two pools with single flower (pool 1) and double flower (pool 2) based on bulked-segregant analysis (BSA) (Figure 1e). A 994 bp heterozygous deletion was identified in the coding region of Pp06G22680.t1 in No. 18 peach based on combing re-sequencing data (Figure S11). Intriguingly, Pp06G22680.t1 encodes an APETALA2 transcription factor (PpAP2), which is known to play roles in flower development. PCR results demonstrated that the 994 bp deletion was present in all hybrids with double-flower but absent in hybrids with single flower (Figure S12; Table S15). This finding indicates that the 994 bp deletion within PpAP2 contributes to the double-flower trait in the ‘No. 18’ peach. Furthermore, the 994 bp deletion genotypes of PpAP2 were present in ‘Hongchuizhi’, ‘1-1-4’ and ‘1-2-7’, but it was absent in ‘Huayulu’ (Figure S13). Interestingly, a single nucleotide polymorphism (SNP) (G/T) located in the binding site of miR172d in PpAP2 was identified in ‘Huayulu’ (Figure 1f). We therefore postulated that this SNP prevents miR172d from effectively targeting and degrading PpAP2. In transiently transformed Nicotiana benthamiana leaves, the PpAP2 (G1346) group showed a dramatic decreased fluorescence signal, indicating that PpAP2 (G1346) can be targeted and degraded by miR172d (Figure 1g). Therefore, the mutation from G to T in the miR172d binding site prevents it from targeting and degrading PpAP2, resulting in the double-flower trait in ‘Huayulu’ peach (Figure 1h). Overall, miR172d and PpAP2 were identified as co-regulators of the single/double flower phenotype in peach. Other than the above, candidate loci associating with showy/nonshowy flowers, hairiness/hairless and three other important agronomic traits were pinpointed (Table S13; Figure S14). In conclusion, a complete gap-free peach genome was obtained. The gene structure was manually refined to ensure high accuracy. Using this genome, we identified the variations in PpTAC1, miR172d, and PpAP2 were associated with their corresponding traits in peach. The gap-free peach genome offers a valuable genomic resource for facilitating the genetic improvement of peach and related species. This work was supported by the National Natural Science Foundation of China (32102329), Modern Agricultural Industry Technology of Henan Province (HARS-22-09-G1) and the Special Fund for Young Talents in Henan Agricultural University (30501339). The authors declare that they have no conflicts of interest. B.T., J.F., H.Z. and X.L. designed this experiment. F.G., C.S., B.F., X.W., J.C., W.W., N.H., X.Z., X.Y. H.Z., L.Z. and J.L. conducted work. H.Z., B.T. and J.F. edited and revised the manuscript. The genome has been deposited in the Peach Genome Database (http://www.stylebio.cn/index.html) and NCBI under the accession number PRJNA1152333. Appendix S1 Methods. Figure S1 Hi-C interactions among the eight ZSH chromosomes. Figure S2 The chloroplast genome of ZSH peach was assembly and annotated. Figure S3 The mitochondrial genome of ZSH peach was assembly and annotated. Figure S4 Gene adjusted manually based on transcriptome data using GSAman software. Figure S5 Gene function was annotated using KEGG, KOBAS, GO, Pfam and Mercator databases. Figure S6 Colinearity between seven peach genomes. Figure S7 Variants number between ZSH and the other peach genomes. Figure S8 The expression patterns of the 3523 genes in two standard (HSM and ‘Okubo’) and two pillar (ZSH and SHLZ) peaches. Figure S9 Candidate genes might contribute to branch angle were identified. Figure S10 The 5033 or 1210 bp insertion were served as molecular marker to verify single and double flower trait in natural populations. Figure S11 A 994 bp heterozygous deletion in PpAP2 was identified in ‘No.18’ and ‘Okubo’ peach using re-sequencing data. Figure S12 The 994 bp deletion was served as molecular marker to verify single and double flower trait in the hybrid population from ‘No.18’ (double flower, D) and ‘Okubo’ (single flower, S). Figure S13 The 994 bp deletion was served as molecular marker to verify single and double flower trait in natural populations. Figure S14 The locus of showy/non-showy flower (a), hairiness/hairless (b), maturity data (c), red flesh (d) and yellow flesh (e) traits were identified using ZSH peach genome. Table S1 Summary of seven peach genome assembly features. Table S2 Statistics of the Prunus persica cv. ‘Zhaoshouhong’ draft genome. Table S3 The length of chromosomes of ZSH genome. Table S4 Repetitive elements in ZSH genome. Table S5 Telomere repeat at the chromosome ends in ZSH genome were identified. Table S6 The candidate centromere at the chromosomes in ZSH genome were identified. Table S7 Evaluation of ZSH peach genome assembly quality. Table S8 Gene function was annotated using five databases. Table S9 The statistics of variations between ZSH and other peach genomes. Table S10 The expression patterns of the 3523 genes in two standard (HSM and ‘Okubo’) and two pillar (ZSH and SHLZ) peaches. Table S11 Gene annotation of these genes which might involve in branch angle development. Table S12 The variation in the promoter and coding region of PpTAC1 in four standard and eleven pillar peach accessions. Table S13 Lead SNPs and candidate genes associated with six agronomic traits using GWAS based on the ZSH genome. Table S14 The phenotype of single/double flower and genotype of miR172d and PpAP2 genes. Table S15 Flower type and genotype of PpAP2 in the hybrids of ‘No.18’ and ‘Okubo’. Table S16 Primers were used in this study. Please note: The publisher is not responsible for the content or functionality of any supporting information supplied by the authors. Any queries (other than missing content) should be directed to the corresponding author for the article.

Topics & Concepts