The genome of Chinese flowering cherry (Cerasus serrulata) provides new insights into Cerasus species

Haplotype-resolved genome of Prunus zhengheensis provides insight into its evolution and low temperature adaptation in apricot

Prunus zhengheensis, an extremely rare population of apricots, originated in warm South-East China and is an excellent material for genetic breeding. However, most apricots and two related species (P. sibirica, P. mandshurica) are found in the cold northern regions in China and the mechanism of their distribution is still unclear. In addition, the classification status of P. zhengheensis is controversial. Thus, we generated a high-quality haplotype-resolved genome for P. zhengheensis, exploring key genetic variations in its adaptation and the causes of phylogenetic incongruence. We found extensive phylogenetic discordances between the nuclear and organelle phylogenies of P. zhengheensis, which could be explained by incomplete lineage sorting. A 242.22-Mb pan-genome of the Armeniaca section was developed with 13 chromosomal genomes. Importantly, we identified a 566-bp insertion in the promoter of the HSFA1d gene in apricot and showed that the activity of the HSFA1d promoter increased under low temperatures. In addition, HSFA1d overexpression in Arabidopsis thaliana indicated that HSFA1d positively regulated plant growth under chilling. Therefore, we hypothesized that the insertion in the promoter of HSFA1d in apricot improved its low-temperature adaptation, allowing it to thrive in relatively cold locations. The findings help explain the weather adaptability of Armeniaca plants.

Molecular characterization and evolutionary relationships of DOFs in four cherry species and functional analysis in sweet cherry

The DOF (DNA binding with one finger) has multiple functions in plants. However, it has received little attention in the research field of cherries. In this study, the evolutionary relationship and molecular characterization of DOF in four cherry species were analyzed, revealing its expression pattern in sweet cherry. There are 23 members in Prunus avium cv. 'Tieton', 88 in Prunus cerasus, 53 in Cerasus × yedoensis, and 27 in Cerasus serrulata. Most of these genes are intron-less or non-intron, with a conserved C2-C2 domain. Due to heterozygosity and chromosomal ploidy, whole-genome duplication (WGD) events occur to varying degrees, and DOF genes are contracted during evolution. Furthermore, these genes are affected by purifying selection pressure. Under low-temperature treatment, the expression of PavDOF2 and PavDOF18 were significantly up-regulated, while that of PavDOF16 is significantly down-regulated. The expression of PavDOF9, PavDOF12, PavDOF14, PavDOF16, PavDOF17, PavDOF18, and PavDOF19 exhibits an increasing trend during flower development and varies during sweet cherry fruit development. PavDOF1, PavDOF8, PavDOF9, and PavDOF15 are localized in the nucleus but is not transcriptionally active. The findings systemically demonstrate the molecular characteristics of DOF in different cherry varieties, providing a basis for further research on the functions of these genes.

Population Variation and Phylogeography of Cherry Blossom (Prunus conradinae) in China

Prunus conradinae (subgenus Cerasus, Rosaceae) is a significant germplasm resource of wild cherry blossom in China. To ensure the comprehensiveness of this study, we used a large sample size (12 populations comprising 244 individuals) which involved the fresh leaves of P. conradinae in Eastern, Central, and Southwestern China. We combined morphological and molecular evidence (three chloroplast DNA (cpDNA) sequences and one nuclear DNA (nr DNA) sequence) to examine the population of P. conradinae variation and differentiation. Our results revealed that Central, East, and Southwest China are important regions for the conservation of P. conradinae to ensure adequate germplasm resources in the future. We also found support for a new variant, P. conradinae var. rubrum. We observed high genetic diversity within P. conradinae (haplotype diversity [Hd] = 0.830; ribotype diversity [Rd] = 0.798), with novel genetic variation and a distinct genealogical structure among populations. There was genetic variation among populations and phylogeographic structure among populations and three geographical groups (Central, East, and Southwest China). The genetic differentiation coefficient was the lowest in the Southwest region and the gene exchange was obvious, while the differentiation was obvious in Central China. In the three geographic groups, we identified two distinct lineages: an East China lineage (Central China and East China) and a Southwest China lineage ((Central China and Southwest China) and East China). These two lineages originated approximately 4.38 million years ago (Mya) in the early Pliocene due to geographic isolation. P. conradinae expanded from Central China to East China at 3.32 Mya (95% HPD: 1.12-5.17 Mya) in the Pliocene. The population of P. conradinae spread from East China to Southwest China, and the differentiation time was 2.17 Mya (95% (HPD: 0.47-4.54 Mya), suggesting that the population of P. conradinae differentiated first in Central and East China. The population of P. conradinae experienced differentiation from Central China to Southwest China around 1.10 Mya (95% HPD: 0.11-2.85 Mya) during the early Pleistocene of the Quaternary period. The southeastern region of East China, near Mount Wuyi, likely serves as a refuge for P. conradinae. This study establishes a theoretical foundation for the classification, identification, conservation, and exploitation of germplasm resources of P. conradinae.

Genomic and transcriptomic studies on flavonoid biosynthesis in Lagerstroemia indica

BACKGROUND

Lagerstroemia indica is a widely cultivated ornamental woody shrub/tree of the family Lythraceae that is used as a traditional medicinal plant in East Asia and Egypt. However, unlike other ornamental woody plants, its genome is not well-investigated, which hindered the discovery of the key genes that regulate important traits and the synthesis of bioactive compounds.

RESULTS

In this study, the genomic sequences of L. indica were determined using several next-generation sequencing technologies. Altogether, 324.01 Mb sequences were assembled and 98.21% (318.21 Mb) of them were placed in 24 pseudo-chromosomes. The heterozygosity, repeated sequences, and GC residues occupied 1.65%, 29.17%, and 38.64% of the genome, respectively. In addition, 28,811 protein-coding gene models, 327 miRNAs, 552 tRNAs, 214 rRNAs, and 607 snRNAs were identified. The intra- and interspecies synteny and Ks analysis revealed that L. indica exhibits a hexaploidy. The co-expression profiles of the genes involved in the phenylpropanoid (PA) and flavonoid/anthocyanin (ABGs) pathways with the R2R3 MYB genes (137 members) showed that ten R2R3 MYB genes positively regulate flavonoid/anthocyanin biosynthesis. The colors of flowers with white, purple (PB), and deep purplish pink (DPB) petals were found to be determined by the levels of delphinidin-based (Dp) derivatives. However, the substrate specificities of LiDFR and LiOMT probably resulted in the different compositions of flavonoid/anthocyanin. In L. indica, two LiTTG1s (LiTTG1-1 and LiTTG1-2) were found to be the homologs of AtTTG1 (WD40). LiTTG1-1 was found to repress anthocyanin biosynthesis using the tobacco transient transfection assay.

CONCLUSIONS

This study showed that the ancestor L. indica experienced genome triplication approximately 38.5 million years ago and that LiTTG1-1 represses anthocyanin biosynthesis. Furthermore, several genes such as LiDFR, LiOMTs, and R2R3 LiMYBs are related to anthocyanin biosynthesis. Further studies are required to clarify the mechanisms and alleles responsible for flower color development.

Identification and Biological Characteristics of Mortierella alpina Associated with Chinese Flowering Cherry (Cerasus serrulata) Leaf Blight in China

The Chinese flowering cherry (Cerasus serrulata), an ornamental tree with established medicinal values, is observed to suffer from leaf blight within Xi'an's greenbelts. This disease threatens both the plant's growth and its ornamental appeal. In this study, 26 isolates were obtained from plants with typical leaf blight, and only 3 isolates (XA-10, XA-15, and XA-18) were found to be pathogenic, causing similar symptoms on the leaves of the host plant. Based on sequence alignment, the ITS and LSU sequences of the three selected isolates were consistent, respectively. Following morphological and molecular analyses, the three selected isolates were further identified as Mortierella alpina. The three selected isolates exhibited similar morphological characteristics, including wavy colonies with dense, milky-white aerial mycelia on PDA medium. Therefore, isolate XA-10 was used as a representative strain for subsequent experiments. The representative strain XA-10 was found to exhibit optimal growth at a temperature of 30 °C and a pH of 7.0. Host range infection tests further revealed that the representative strain XA-10 could also inflict comparable disease symptoms on both the leaves and fruits of three different Rosaceae species (Prunus persica, Pyrus bretschneideri, and Prunus salicina). This study reveals, for the first time, the causative agent of leaf blight disease affecting the Chinese flowering cherry. This provides a deeper understanding of the biology and etiology of M. alpina. This study lays a solid foundation for the sustainable control and management of leaf blight disease in the Chinese flowering cherry.

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.

Phylogeography of 912 Cherry Accessions Insight into Independent Origins of Fruiting Cherries and Domestication Footprints of Cultivated Chinese Cherry (Prunus pseudocerasus Lindl.)

The subgenus Cerasus (Rosaceae) contain numerous fruit trees and ornamentals with high economic values. The origin and genetic divergence among various types of fruiting cherries always remain a perplexing issue. We employed three plastom fragments and ITS sequence matrices derived from 912 cherry accessions to elucidate the phylogeographic structure and genetic relationship among fruiting cherries, as well as the origin and domestication of cultivated Chinese cherry. The integration of haplotype genealogies, Approximate Bayesian computation (ABC) approach and estimation of genetic differentiation within and between different groups and lineages has facilitated the resolution of several previously unresolved questions. Firstly, distant phylogenetic relationships between Cerasus and Microcerasus accessions, as indicated by both nuclear and chloroplast data, suggested independent origins and evolution for these two taxa. Moreover, two distinct geographic origin centers (Europe and China) have been confirmed, with significant phylogeographic signals and high genetic differentiation observed between cherries from these regions. This may be attributed to long-term geographic isolation caused by Himalaya-Hengduan Mountains. Our phylogeographic analyses and ABC analysis suggested that cherries inhabiting in China may have undergone multiple hybridization events during the glacial refugia of the eastern edge and southern Himalaya-Hengduan Mountains, followed by rapid radiation throughout their current habitats during interglacial period. The discrepancy between nuclear and chloroplast data may be attributed to hybridization events and incomplete lineage sorting. Furthermore, we speculated that the domesticated Chinese cherries were derived from wild accessions in Longmenshan Fault Zones approximately 2600 years ago. We have also traced the domestication processes and dispersal routes of cultivated Chinese cherries.

Chromosome-scale genome assembly of Prunus pusilliflora provides novel insights into genome evolution, disease resistance, and dormancy release in Cerasus L

Prunus pusilliflora is a wild cherry germplasm resource distributed mainly in Southwest China. Despite its ornamental and economic value, a high-quality assembled P. pusilliflora genome is unavailable, hindering our understanding of its genetic background, population diversity, and evolutionary processes. Here, we de novo assembled a chromosome-scale P. pusilliflora genome using Oxford Nanopore, Illumina, and chromosome conformation capture sequencing. The assembled genome size was 309.62 Mb, with 76 scaffolds anchored to eight pseudochromosomes. We predicted 33 035 protein-coding genes, functionally annotated 98.27% of them, and identified repetitive sequences covering 49.08% of the genome. We found that P. pusilliflora is closely related to Prunus serrulata and Prunus yedoensis, having diverged from them ~41.8 million years ago. A comparative genomic analysis revealed that P. pusilliflora has 643 expanded and 1128 contracted gene families. Furthermore, we found that P. pusilliflora is more resistant to Colletotrichum viniferum, Phytophthora capsici, and Pseudomonas syringae pv. tomato (Pst) DC3000 infections than cultivated Prunus avium. P. pusilliflora also has considerably more nucleotide-binding site-type resistance gene analogs than P. avium, which explains its stronger disease resistance. The cytochrome P450 and WRKY families of 263 and 61 proteins were divided into 42 and 8 subfamilies respectively in P. pusilliflora. Furthermore, 81 MADS-box genes were identified in P. pusilliflora, accompanying expansions of the SVP and AGL15 subfamilies and loss of the TM3 subfamily. Our assembly of a high-quality P. pusilliflora genome will be valuable for further research on cherries and molecular breeding.

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

Flowering cherry is a very popular species around the world. High-quality genome resources for different elite cultivars are needed, and the understanding of their origins and the regulation of key ornamental traits are limited for this tree. Here, a high-quality chromosome-scale genome of Prunus campanulata 'Plena' (PCP), which is a native and elite flowering cherry cultivar in China, was generated. The contig N50 of the genome was 18.31 Mb, and 99.98% of its contigs were anchored to eight chromosomes. Furthermore, a total of 306 accessions of flowering cherry germplasm and six lines of outgroups were collected. Resequencing of these 312 lines was performed, and 761 267 high-quality genomic variants were obtained. The origins of flowering cherry were predicted, and these 306 accessions could be classified into three clades, A, B and C. According to phylogenetic analysis, we predicted two origins of flowering cherry. Flowering cherry in clade A originated in southern China, such as in the Himalayan Mountains, while clades B and C originated in northeastern China. Finally, a genome-wide association study of flower colour was performed for all 312 accessions of flowering cherry germplasm. A total of seven quantitative trait loci (QTLs) were identified. One gene encoding glycosylate transferase was predicted as the candidate gene for one QTL. Taken together, our results provide a valuable genomic resource and novel insights into the origin, evolution and flower colour variations of flowering cherry.

Chromosome-Level Assembly of Flowering Cherry (Prunus campanulata) Provides Insight into Anthocyanin Accumulation

The flowering cherries (genus Prunus, subgenus Cerasus) are popular ornamental trees in China, Japan, Korea, and elsewhere. Prunus campanulata Maxim. is an important species of flowering cherry native to Southern China, which is also distributed in Taiwan, the Ryukyu Islands of Japan, and Vietnam. It produces bell-shaped flowers with colors ranging from bright pink to crimson during the Chinese Spring Festival from January to March each year. We selected the P. campanulata cultivar "Lianmeiren", with only 0.54% of heterozygosity, as the focus of this study, and generated a high-quality chromosome-scale genome assembly of P. campanulata by combining Pacific Biosciences (PacBio) single-molecule sequencing, 10× Genomics sequencing, and high-throughput chromosome conformation capture (Hi-C) technology. We first assembled a 300.48 Mb genome assembly with a contig N50 length of 2.02 Mb. In total, 28,319 protein-coding genes were predicted from the genome, 95.8% of which were functionally annotated. Phylogenetic analyses indicated that P. campanulata diverged from a common ancestor of cherry approximately 15.1 million years ago. Comparative genomic analyses showed that the expanded gene families were significantly involved in ribosome biogenesis, diterpenoid biosynthesis, flavonoid biosynthesis, and circadian rhythm. Furthermore, we identified 171 MYB genes from the P. campanulata genome. Based on the RNA-seq of five organs at three flowering stages, expression analyses revealed that the majority of the MYB genes exhibited tissue-specific expression patterns, and some genes were identified as being associated with anthocyanin accumulation. This reference sequence is an important resource for further studies of floral morphology and phenology, and comparative genomics of the subgenera Cerasus and Prunus.

Geographical distribution and predict potential distribution of Cerasus serrulata

Climate change is closely related to the distribution of plant resources. Cerasus serrulata is an important plant resource in China. The study on the influence of environmental factors on the distribution of suitable areas of C. serrulata is conducive to the protection and development of C. serrulata. In this paper, the distribution information of 238 Chinese Cerasus serrulata plants was processed by DIVA-GIS. The MaxEnt model was used to simulate the current and future distribution, and the ecological distribution and richness of Cerasus were analyzed. The results showed that the Cerasus serrulata was widely distributed in Hebei, Heilongjiang, Shandong, Jiangsu, Zhejiang, Anhui, Jiangxi, Hunan, and Guizhou provinces, mainly in the low and middle elevation areas of 10 to 1200 m. Based on this model, the precipitation of the warmest quarter, the precipitation of the driest month, and the mean temperature of the coldest were the most significant bioclimatic variables affecting the distribution of C. serrulata. In the future, climate change may lead to a slight increase of 2.31% in the area of suitable habitat for Cerasus serrulata, while the optimal habitat will decrease from 20.81 to 14.55%. Therefore, conservation measures should be taken to protect these precious resources.

Accurate Chromosome Identification in the Prunus Subgenus Cerasus (Prunus pseudocerasus) and its Relatives by Oligo-FISH

A precise, rapid and straightforward approach to chromosome identification is fundamental for cytogenetics studies. However, the identification of individual chromosomes was not previously possible for Chinese cherry or other Prunus species due to the small size and similar morphology of their chromosomes. To address this issue, we designed a pool of oligonucleotides distributed across specific pseudochromosome regions of Chinese cherry. This oligonucleotide pool was amplified through multiplex PCR with specific internal primers to produce probes that could recognize specific chromosomes. External primers modified with red and green fluorescence tags could produce unique signal barcoding patterns to identify each chromosome concomitantly. The same oligonucleotide pool could also discriminate all chromosomes in other Prunus species. Additionally, the 5S/45S rDNA probes and the oligo pool were applied in two sequential rounds of fluorescence in situ hybridization (FISH) localized to chromosomes and showed different distribution patterns among Prunus species. At the same time, comparative karyotype analysis revealed high conservation among P. pseudocerasus, P. avium, and P. persica. Together, these findings establish this oligonucleotide pool as the most effective tool for chromosome identification and the analysis of genome organization and evolution in the genus Prunus.

Comparative Analysis of Transposable Elements and the Identification of Candidate Centromeric Elements in the Prunus Subgenus Cerasus and Its Relatives

The subgenus Cerasus and its relatives include many crucial economic drupe fruits and ornamental plants. Repetitive elements make up a large part of complex genomes, and some of them play an important role in gene regulation that can affect phenotypic variation. However, the variation in their genomes remains poorly understood. This work conducted a comprehensive repetitive sequence identification across the draft genomes of eight taxa of the genus Prunus, including four of the Prunus subgenus Cerasus (Prunus pseudocerasus, P. avium, P. yedoensis and P. × yedoensis) as well as congeneric species (Prunus salicina, P. armeniaca, P. dulcis and P. persica). Annotation results showed high proportions of transposable elements in their genomes, ranging from 52.28% (P. armeniaca) to 61.86% (P. pseudocerasus). The most notable differences in the contents of long terminal repeat retrotransposons (LTR-RTs) and tandem repeats (TRs) were confirmed with de novo identification based on the structure of each genome, which significantly contributed to their genome size variation, especially in P. avium and P.salicina. Sequence comparisons showed many similar LTR-RTs closely related to their phylogenetic relationships, and a highly similar monomer unit of the TR sequence was conserved among species. Additionally, the predicted centromere-associated sequence was located in centromeric regions with FISH in the 12 taxa of Prunus. It presented significantly different signal intensities, even within the diverse interindividual phenotypes for Prunus tomentosa. This study provides insight into the LTR-RT and TR variation within Prunus and increases our knowledge about its role in genome evolution.

Twenty years of plant genome sequencing: achievements and challenges

Publication of the complete genome sequence of Arabidopsis thaliana, the first plant reference genome, in December 2000 heralded the beginning of the plant genome era. Over the past 20 years reference genomes have been generated for hundreds of plant species, spanning non-vascular to flowering plants. Releasing these plant genomes has dramatically advanced studies in all disciplines of plant biology. Importantly, multiple reference-level genomes have been generated for the major crops and their progenitors, enabling the creation of pan-genomes and exploration of domestication history and natural variations that can be adopted by modern crop breeding. We summarize the progress of plant genome sequencing and the challenges of sequencing more complex plant genomes and generating pan-genomes.

Evolution of Rosaceae Plastomes Highlights Unique Cerasus Diversification and Independent Origins of Fruiting Cherry

Rosaceae comprises numerous types of economically important fruits, ornamentals, and timber. The lack of plastome characteristics has blocked our understanding of the evolution of plastome and plastid genes of Rosaceae crops. Using comparative genomics and phylogenomics, we analyzed 121 Rosaceae plastomes of 54 taxa from 13 genera, predominantly including Cerasus (true cherry) and its relatives. To our knowledge, we generated the first comprehensive map of genomic variation across Rosaceae plastomes. Contraction/expansion of inverted repeat regions and sequence losses of the two single-copy regions underlie large genomic variations in size among Rosaceae plastomes. Plastid protein-coding genes were characterized with a high proportion (over 50%) of synonymous variants and insertion-deletions with multiple triplets. Five photosynthesis-related genes were specially selected in perennial woody trees. Comparative genomic analyses implied divergent evolutionary patterns between pomaceous and drupaceous trees. Across all examined plastomes, unique and divergent evolution was detected in Cerasus plastomes. Phylogenomic analyses and molecular dating highlighted the relatively distant phylogenetic relationship between Cerasus and relatives (Microcerasus, Amygdalus, Prunus, and Armeniaca), which strongly supported treating the monophyletic true cherry group as a separate genus excluding dwarf cherry. High genetic differentiation and distinct phylogenetic relationships implied independent origins and domestication between fruiting cherries, particularly between Prunus pseudocerasus (Cerasus pseudocerasus) and P. avium (C. avium). Well-resolved maternal phylogeny suggested that cultivated P. pseudocerasus originated from Longmenshan Fault zone, the eastern edge of Himalaya-Hengduan Mountains, where it was subjected to frequent genomic introgression between its presumed wild ancestors and relatives.

Genome sequencing and analysis of two early-flowering cherry (Cerasus × kanzakura) varieties, 'Kawazu-zakura' and 'Atami-zakura'

To gain genetic insights into the early-flowering phenotype of ornamental cherry, also known as sakura, we determined the genome sequences of two early-flowering cherry (Cerasus × kanzakura) varieties, ‘Kawazu-zakura’ and ‘Atami-zakura’. Because the two varieties are interspecific hybrids, likely derived from crosses between Cerasus campanulata (early-flowering species) and Cerasus speciosa, we employed the haplotype-resolved sequence assembly strategy. Genome sequence reads obtained from each variety by single-molecule real-time sequencing (SMRT) were split into two subsets, based on the genome sequence information of the two probable ancestors, and assembled to obtain haplotype-phased genome sequences. The resultant genome assembly of ‘Kawazu-zakura’ spanned 519.8 Mb with 1,544 contigs and an N50 value of 1,220.5 kb, while that of ‘Atami-zakura’ totalled 509.6 Mb with 2,180 contigs and an N50 value of 709.1 kb. A total of 72,702 and 69,528 potential protein-coding genes were predicted in the genome assemblies of ‘Kawazu-zakura’ and ‘Atami-zakura’, respectively. Gene clustering analysis identified 2,634 clusters uniquely presented in the C. campanulata haplotype sequences, which might contribute to its early-flowering phenotype. Genome sequences determined in this study provide fundamental information for elucidating the molecular and genetic mechanisms underlying the early-flowering phenotype of ornamental cherry tree varieties and their relatives.

Research advances in and prospects of ornamental plant genomics

The term 'ornamental plant' refers to all plants with ornamental value, which generally have beautiful flowers or special plant architectures. China is rich in ornamental plant resources and known as the "mother of gardens". Genomics is the science of studying genomes and is useful for carrying out research on genome evolution, genomic variations, gene regulation, and important biological mechanisms based on detailed genome sequence information. Due to the diversity of ornamental plants and high sequencing costs, the progress of genome research on ornamental plants has been slow for a long time. With the emergence of new sequencing technologies and a reduction in costs since the whole-genome sequencing of the first ornamental plant (Prunus mume) was completed in 2012, whole-genome sequencing of more than 69 ornamental plants has been completed in <10 years. In this review, whole-genome sequencing and resequencing of ornamental plants will be discussed. We provide analysis with regard to basic data from whole-genome studies of important ornamental plants, the regulation of important ornamental traits, and application prospects.