Genome assembly, resequencing and genome-wide association analyses provide novel insights into the origin, evolution and flower colour variations of flowering cherry

New insights into genome assembly at the chromosome-level of Prunus tomentosa in evolution and cold tolerance

This study assembled a high-quality chromosome-level genome of Prunus tomentosa, offering a vital resource for elucidating its genetic architecture, evolutionary relationships, and facilitating genome-assisted breeding efforts. Multi-omics integration revealed PtIMP3 and PtMIOX1L as key factors in cold tolerance of P. tomentosa. PtIMP3 drives the conversion of glucose-6-phosphate to myo-inositol, while PtMIOX1L catalyzes myo-inositol to d-glucuronic acid. Specifically, the high expression abundance of PtIMP3 and low expression abundance of PtMIOX1L resulted in high endogenous inositol levels in P. tomentosa. The application of myo-inositol enhanced the cold tolerance of cherry rootstocks by modulating reactive oxygen species concentrations and maintaining a stable relative water content. This finding supports the superior performance of P. tomentosa in adapting to extreme low-temperatures environmental conditions. These insights advance strategies for improving cold tolerance in horticultural crops, bridging fundamental research with practical applications in developing climate-resilient crops.

Novel insight of the SVP gene involved in pedicel length based on genomics analysis in cherry

KEY MESSAGE

PcSVP was identified based on Prunus conradinae genome and was further overexpressed in A. thaliana to comfirm it was a key factor in flower development, causing the pedicels elongation. Prunus conradinae is an endemic plant resource in China with high ornamental and economic values. To generate useful genomic resources for expanding insights into the evolutionary history of this important plant, the chromosome-level genome and organelle genomes of P. conradinae are de novo assembled and functionally annotated. The chromosome-level haploid genome of autotetraploid P. conradinae was assembled with 262.79 Mb with 27,802 protein-coding genes annotated. The complete chloroplast and mitochondrial genome of P. conradinae are found to be 157,715 bp and 434,334 bp, respectively. According to evolutionary analysis, P. conradinae was closely related to P. serrulata and P. yedoensis, and they diverged from their common ancestor approximately 6.0 million years ago. There were 108 gene families that significantly expanded during P. conradinae evolution and 56 shared positively selected genes. Selective sweep analysis based on the whole-genome resequencing of wild cherries from Fujian and Zhejiang indicated that genes involved in flower development and stress responses were potentially under selection. Pedicel length varied greatly among Prunus species and was a significant identifying characteristic. Ectopic overexpression of PcSVP in Arabidopsis thaliana suggested that it was a key factor in flower development, causing the sepals curling and pedicels elongation. These findings will contribute to the discovery of key functional genes involved in the agronomic or biological traits of P. conradinae, as well as the future development, utilisation and germplasm conservation of wild cherries.

The telomere-to-telomere genome of flowering cherry (Prunus campanulata) reveals genomic evolution of the subgenus Cerasus

BACKGROUND

Prunus campanulata, a species of ornamental cherry, holds significant genetic and horticultural value. Despite the availability of various cherry genomes, a fully resolved telomere-to-telomere (T2T) assembly for this species has been lacking. Recent advancements in long-read sequencing technologies have made it possible to generate gap-free genome assemblies, providing comprehensive insights into genomic structures that were previously inaccessible.

FINDINGS

We present the first T2T genome assembly for P. campanulata "Lianmeiren" (v2.0), achieved through the integration of PacBio HiFi, ultra-long Oxford Nanopore Technologies, Illumina, and Hi-C sequencing. The assembly resulted in a highly contiguous genome with a total size of 266.23 Mb and a contig N50 of 31.6 Mb. The genome exhibits remarkable completeness (98.9% BUSCO) and high accuracy (quality value of 48.75). Additionally, 13 telomeres and putative centromere regions were successfully identified across the 8 pseudochromosomes. Comparative analysis with the previous v1.0 assembly revealed 336,943 single nucleotide polymorphisms, 107,521 indels, and 1,413 structural variations, along with the annotation of 1,402 new genes.

CONCLUSIONS

This T2T genome assembly of P. campanulata "Lianmeiren" provides a critical reference for understanding the genetic architecture of the species. It enhances our ability to study structural variations, gene function, and evolutionary biology within the Prunus genus.

Functional Analysis of the PoSERK-Interacting Protein PorbcL in the Embryogenic Callus Formation of Tree Peony (Paeonia ostii T. Hong et J. X. Zhang)

SERK is a marker gene for early somatic embryogenesis. We screened and functionally verified a SERK-interacting protein to gain insights into tree-peony somatic embryogenesis. Using PoSERK as bait, we identified PorbcL (i.e., the large subunit of Rubisco) as a SERK-interacting protein from a yeast two-hybrid (Y2H) library of cDNA from developing tree-peony somatic embryos. The interaction between PorbcL and PoSERK was verified by Y2H and bimolecular fluorescence complementation analyses. PorbcL encodes a 586-amino-acid acidic non-secreted hydrophobic non-transmembrane protein that is mainly localized in the chloroplast and plasma membrane. PorbcL was highly expressed in tree-peony roots and flowers and was up-regulated during zygotic embryo development. PorbcL overexpression caused the up-regulation of PoSERK (encoding somatic embryogenesis receptor-like kinase), PoAGL15 (encoding agamous-like 15), and PoGPT1 (encoding glucose-6-phosphate translocator), while it caused the down-regulation of PoLEC1 (encoding leafy cotyledon 1) in tree-peony callus. PorbcL overexpression led to increased indole-3-acetic acid (IAA) content but decreasing contents of abscisic acid (ABA) and 6-benzyladenosine (BAPR). The changes in gene expression, high IAA levels, and increased ratio of IAA to ABA, BAPR, 1-Aminocyclopropanecarboxylic acid (ACC), 5-Deoxystrigol (5DS), and brassinolide (BL) promoted embryogenesis. These results provide a foundation for establishing a tree-peony embryogenic callus system.

Haplotype-resolved genome assembly for tetraploid Chinese cherry (Prunus pseudocerasus) offers insights into fruit firmness

Chinese cherry (Prunus pseudocerasus) holds considerable importance as one of the primary stone fruit crops in China. However, artificially improving its traits and genetic analysis are challenging due to lack of high-quality genomic resources, which mainly result from difficulties associated with resolving its tetraploid and highly heterozygous genome. Herein, we assembled a chromosome-level, haplotype-resolved genome of the cultivar 'Zhuji Duanbing', comprising 993.69 Mb assembled into 32 pseudochromosomes using PacBio HiFi, Oxford Nanopore, and Hi-C. Intra-haplotype comparative analyses revealed extensive intra-genomic sequence and expression consistency. Phylogenetic and comparative genomic analyses demonstrated that P. pseudocerasus was a stable autotetraploid species, closely related to wild P. pusilliflora, with the two diverging ~18.34 million years ago. Similar to other Prunus species, P. pseudocerasus underwent a common whole-genome duplication event that occurred ~139.96 million years ago. Because of its low fruit firmness, P. pseudocerasus is unsuitable for long-distance transportation, thereby restricting its rapid development throughout China. At the ripe fruit stage, P. pseudocerasus cv. 'Zhuji Duanbing' was significantly less firm than P. avium cv. 'Heizhenzhu'. The difference in firmness is attributed to the degree of alteration in pectin, cellulose, and hemicellulose contents. In addition, comparative transcriptomic analyses identified GalAK-like and Stv1, two genes involved in pectin biosynthesis, which potentially caused the difference in firmness between 'Zhuji Duanbing' and 'Heizhenzhu'. Transient transformations of PpsGalAK-like and PpsStv1 increase protopectin content and thereby enhance fruit firmness. Our study lays a solid foundation for functional genomic studies and the enhancement of important horticultural traits in Chinese cherries.

Chromosome-level genome assembly provides insights into the genetic diversity, evolution, and flower development of Prunus conradinae

Prunus conradinae, a valuable flowering cherry belonging to the Rosaceae family subgenus Cerasus and endemic to China, has high economic and ornamental value. However, a high-quality P. conradinae genome is unavailable, which hinders our understanding of its genetic relationships and phylogenesis, and ultimately, the possibility of mining of key genes for important traits. Herein, we have successfully assembled a chromosome-scale P. conradinae genome, identifying 31,134 protein-coding genes, with 98.22% of them functionally annotated. Furthermore, we determined that repetitive sequences constitute 46.23% of the genome. Structural variation detection revealed some syntenic regions, inversions, translocations, and duplications, highlighting the genetic diversity and complexity of Cerasus. Phylogenetic analysis demonstrated that P. conradinae is most closely related to P. campanulata, from which it diverged ~ 19.1 million years ago (Mya). P. avium diverged earlier than P. cerasus and P. conradinae. Similar to the other Prunus species, P. conradinae underwent a common whole-genome duplication event at ~ 138.60 Mya. Furthermore, 79 MADS-box members were identified in P. conradinae, accompanied by the expansion of the SHORT VEGETATIVE PHASE subfamily. Our findings shed light on the complex genetic relationships, and genome evolution of P. conradinae and will facilitate research on the molecular breeding and functions of key genes related to important horticultural and economic characteristics of subgenus Cerasus.

Haplotype-resolved genome assembly of poplar line NL895 provides a valuable tree genomic resource

Poplar line NL895 can potentially become a model plant for poplar study as it is a widely cultivated elite line. However, the lack of genome resources hindered the use of NL895 as the major plant material in poplar. In this study, we provided a high-quality genome assembly for poplar line NL895 with PacBio single molecule real-time (SMRT) sequencing and High-throughput chromosome conformation capture (Hi-C) technology. The raw assembly of NL895 for the diploid genome included 606 contigs with a total size of ~815 Mb, and the monoploid genome included 246 contigs with a total size of ~412 Mb. The haplotype-resolved chromosomes in the diploid genomes were also generated. All the monoploid, diploid, and haplotype-resolved genomes showed more than 97% completeness and they can largely improve the mapping efficiency in RNA-Seq analysis. By comprehensively comparing the two haplotype genomes we found the heterozygosity of NL895 is much higher than other poplar lines. We also found that NL895 harbors more genomic variants and more gene diversity. The haplotype-specific genes showed higher variable gene expression patterns. These characters would be attributed to the high heterosis of poplar line NL895. The allele-specific expression (ASE) was also investigated and lots of alleles showed biased expressions in different tissues or environmental conditions. Taken together, the genome sequence for NL895 is a valuable tree genomic resource and it would greatly facilitate studies in poplar.

Characterization of Key Odorants in Lushan Yunwu Tea in Response to Intercropping with Flowering Cherry

Lushan Yunwu tea (LSYWT) is a famous green tea in China. However, the effects of intercropping tea with flowering cherry on the overall aroma of tea have not been well understood. In this study, headspace solid-phase microextraction (HS-SPME) coupled with gas chromatography-mass spectrometry (GC-MS) was used for analysis. A total of 54 volatile compounds from eight chemical classes were identified in tea samples from both the intercropping and pure-tea-plantation groups. Principal component analysis (PCA), orthogonal partial least-squares discriminant analysis (OPLS-DA), and odor activity value (OAV) methods combined with sensory evaluation identified cis-jasmone, nonanal, and linalool as the key aroma compounds in the intercropping group. Benzaldehyde, α-farnesene, and methyl benzene were identified as the main volatile compounds in the flowering cherry using headspace solid-phase microextraction/gas chromatography-mass spectrometry (HS-SPME/GC-MS). These findings will enrich the research on tea aroma chemistry and offer new insights into the product development and quality improvement of LSYWT.

In vitro regeneration, photomorphogenesis and light signaling gene expression in Hydrangea quercifolia cv. 'Harmony' under different LED environments

MAIN CONCLUSION

Plant growth regulators, sucrose concentration, and light quality significantly impact in vitro regeneration of 'Harmony'. Blue light promotes photomorphogenesis by enhancing light energy utilization, adjusting transcription of light signal genes, and altering hormone levels. Hydrangea quercifolia cv. 'Harmony', celebrated for lush green foliage and clusters of white flowers, has been extensively researched for its regenerative properties. Regeneration in stem segments, leaves, and petioles is facilitated by exogenous auxin and cytokinins (CTKs), with the concentration of sucrose (SC) being a key determinant for shoot regeneration from leaves. The study also highlights the significant impact of light conditions on photomorphogenesis. With an increase in the proportion of red (R) light, there is an inhibitory effect, leading to a reduction in leaf area, a decrease in the quantum yield of PSII (ΦPSII), and an increase in non-photochemical quenching (ΦNPQ) and non-regulated energy dissipation in PSII (ΦNO). Conversely, blue (B) light enhances growth, characterized by an increase in leaf area, elevated ΦPSII, and stable ΦNPQ and ΦNO levels. Additionally, B light induces the upregulation of HqCRYs, HqHY5-like, HqXTH27-like, and HqPHYs genes, along with an increase in endogenous CTKs levels, which positively influence photomorphogenesis independent of HqHY5-like regulation. This light condition also suppresses the synthesis of endogenous gibberellins (GA) and brassinosteroids (BR), further facilitating photomorphogenesis. In essence, B light is fundamental in expediting photomorphogenesis in 'Harmony', demonstrating the vital role in plant growth and development.

A chromosome-scale assembly of the early-flowering Prunus campanulata and comparative genomics of cherries

Prunus campanulata is an important flowering cherry germplasm of high ornamental value. Given its early-flowering phenotypes, P. campanulata could be used for molecular breeding of ornamental species and fruit crops belonging to the subgenus Cerasus. Here, we report a chromosome-scale assembly of P. campanulata with a genome size of 282.6 Mb and a contig N50 length of 12.04 Mb. The genome contained 24,861 protein-coding genes, of which 24,749 genes (99.5%) were functionally annotated, and 148.20 Mb (52.4%) of the assembled sequences are repetitive sequences. A combination of genomic and population genomic analyses revealed a number of genes under positive selection or accelerated molecular evolution in P. campanulata. Our study provides a reliable genome resource, and lays a solid foundation for genetic improvement of flowering cherry germplasm.

Genome-Wide Identified MADS-Box Genes in Prunus campanulata 'Plena' and Theirs Roles in Double-Flower Development

The MADS-box gene family plays key roles in flower induction, floral initiation, and floral morphogenesis in flowering plants. To understand their functions in the double-flower formation of Prunus campanulata 'Plena' (hereafter referred to as PCP), which is an excellent flowering cherry cultivar, we performed genome-wide identification of the MADS-box gene family. In this study, 71 MADS-box genes were identified and grouped into the Mα, Mβ, Mγ and MIKC subfamilies according to their structures and phylogenetic relationships. All 71 MADS-box genes were located on eight chromosomes of PCP. Analysis of the cis-acting elements in the promoter region of MADS-box genes indicated that they were associated mainly with auxin, abscisic acid, gibberellin, MeJA (methyl jasmonate), and salicylic acid responsiveness, which may be involved in floral development and differentiation. By observing the floral organ phenotype, we found that the double-flower phenotype of PCP originated from petaloid stamens. The analysis of MIKC-type MADS-box genes in PCP vegetative and floral organs by qRT-PCR revealed six upregulated genes involved in petal development and three downregulated genes participating in stamen identity. Comparative analysis of petaloid stamens and normal stamens also indicated that the expression level of the AG gene (PcMADS40) was significantly reduced. Thus, we speculated that these upregulated and downregulated genes, especially PcMADS40, may lead to petaloid stamen formation and thus double flowers. This study lays a theoretical foundation for MADS-box gene identification and classification and studying the molecular mechanism underlying double flowers in other ornamental plants.