Draft genome sequence of wild Prunus yedoensis reveals massive inter-specific hybridization between sympatric flowering cherries

Identifying potential anthocyanin biosynthesis regulator in Chinese cherry by comprehensive genome-wide characterization of the R2R3-MYB transcription factor gene family

BACKGROUND

Chinese cherry [Cerasus pseudocerasus (Lindl.) G.Don] (syn. Prunus pseudocerasus Lindl.) is an economically important fruiting cherry species with a diverse range of attractive colors, spanning from the lightest yellow to the darkest black purple. However, the MYB transcription factors involved in anthocyanin biosynthesis underlying fruit color variation in Chinese cherry remain unknown.

RESULTS

In this study, we characterized the R2R3-MYB gene family of Chinese cherry by genome-wide identification and compared it with those of 10 Rosaceae relatives and Arabidopsis thaliana. A total of 1490 R2R3-MYBs were classified into 43 subfamilies, which included 29 subfamilies containing both Rosaceae MYBs and AtMYBs. One subfamily (S45) contained only Rosaceae MYBs, while three subfamilies (S12, S75, and S77) contained only AtMYBs. The variation in gene numbers within identical subfamilies among different species and the absence of certain subfamilies in some species indicated the species-specific expansion within MYB gene family in Chinese cherry and its relatives. Segmental and tandem duplication events primarily contributed to the expansion of Chinese cherry R2R3-CpMYBs. The duplicated gene pairs underwent purifying selection during evolution after duplication events. Phylogenetic relationships and transcript profiling revealed that CpMYB10 and CpMYB4 are involved in the regulation of anthocyanin biosynthesis in Chinese cherry fruits. Expression patterns, transient overexpression and VIGS results confirmed that CpMYB10 promotes anthocyanin accumulation in the fruit skin, while CpMYB4 acts as a repressor, inhibiting anthocyanin biosynthesis of Chinese cherry.

CONCLUSIONS

This study provides a comprehensive and systematic analysis of R2R3-MYB gene family in Chinese cherry and Rosaceae relatives, and identifies two regulators, CpMYB10 and CpMYB4, involved in anthocyanin biosynthesis in Chinese cherry. These results help to develop and utilize the potential functions of anthocyanins in Chinese cherry.

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.

The pan-plastome of Prunus mume: insights into Prunus diversity, phylogeny, and domestication history

Backgrounds

Prunus mume in the Rosaceae and commonly referred to as mei or Chinese plum is widely used as a traditional ornamental flowering plant and fruit tree in China. Although some population and genetic analyses have been conducted for this species, no extensive comparisons of genetic variation from plastomes have yet been investigated.

Methods

We de novo assembled a total of 322 complete P. mume plastomes in this study and did a series of comparative analyses to better resolve pan-plastomic patterns of P. mume. To determine the phylogeny and domestication history of this species, we reconstructed the phylogenetic tree of Prunus genus, and resolved the population structure of P. mume. We also examined the nucleotide variation of P. mume to find potential DNA barcodes.

Results

The assembled plastomes exhibited a typical quadripartite structure and ranged from 157,871 bp to 158,213 bp in total size with a GC content ranging from 36.73 to 36.75%. A total of 112 unique genes were identified. Single nucleotide variants (SNVs) were the most common variants found among the plastomes, followed by nucleotide insertions/deletions (InDels), and block substitutions with the intergenic spacer (IGS) regions containing the greatest number of variants. From the pan-plastome data six well-supported genetic clusters were resolved using multiple different population structure analyses. The different cultivars were unevenly distributed among multiple clades. We also reconstructed a phylogeny for multiple species of Prunus to better understand genus level diversity and history from which a complex introgressive relationship between mei and other apricots/plums was resolved.

Conclusion

This study constructed the pan-plastome of P. mume, which indicated the domestication of P. mume involved multiple genetic origins and possible matrilineal introgression from other species. The phylogenetic analysis in Prunus and the population structure of P. mume provide an important maternal history for Prunus and the groundwork for future studies on intergenomic sequence transfers, cytonuclear incompatibility, and conservation genetics.

A phased genome of the highly heterozygous 'Texas' almond uncovers patterns of allele-specific expression linked to heterozygous structural variants

The vast majority of traditional almond varieties are self-incompatible, and the level of variability of the species is very high, resulting in a high-heterozygosity genome. Therefore, information on the different haplotypes is particularly relevant to understand the genetic basis of trait variability in this species. However, although reference genomes for several almond varieties exist, none of them is phased and has genome information at the haplotype level. Here, we present a phased assembly of genome of the almond cv. Texas. This new assembly has 13% more assembled sequence than the previous version of the Texas genome and has an increased contiguity, in particular in repetitive regions such as the centromeres. Our analysis shows that the 'Texas' genome has a high degree of heterozygosity, both at SNPs, short indels, and structural variants level. Many of the SVs are the result of heterozygous transposable element insertions, and in many cases, they also contain genic sequences. In addition to the direct consequences of this genic variability on the presence/absence of genes, our results show that variants located close to genes are often associated with allele-specific gene expression, which highlights the importance of heterozygous SVs in almond.

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.

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.

Designing of future ornamental crops: a biotechnological driven perspective

With a basis in human appreciation of beauty and aesthetic values, the new era of ornamental crops is based on implementing innovative technologies and transforming symbols into tangible assets. Recent advances in plant biotechnology have attracted considerable scientific and industrial interest, particularly in terms of modifying desired plant traits and developing future ornamental crops. By utilizing omics approaches, genomic data, genetic engineering, and gene editing tools, scientists have successively explored the underlying molecular mechanism and potential gene(s) behind trait regulation such as floral induction, plant architecture, stress resistance, plasticity, adaptation, and phytoremediation in ornamental crop species. These signs of progress lay a theoretical and practical foundation for designing and enhancing the efficiency of ornamental plants for a wide range of applications. In this review, we briefly summarized the existing literature and advances in biotechnological approaches for the improvement of vital traits in ornamental plants. The future ornamental plants, such as light-emitting plants, biotic/abiotic stress detectors, and pollution abatement, and the introduction of new ornamental varieties via domestication of wild species are also discussed.

Selection signatures of Qinchuan cattle based on whole-genome sequences

Qinchuan cattle has gradually improved in body shape and growth rate in the long-term breeding process from the draft cattle to beef cattle. As the head of the five local yellow cattle in China, the Qinchuan cattle has been designated as a specialized beef cattle breed. We investigated the selection signatures using whole genome sequencing data in Qinchuan cattle. Based on Fst, we detected hundreds of candidate genes under selection across Qinchuan, Red Angus, and Japanese Black cattle. Through protein-protein interaction analysis and functional annotation of candidate genes, the results revealed that KMT2E, LTBP1 and NIPBL were related to brain size, body characteristics, and limb development, respectively, suggesting that these potential genes may affect the growth and development traits in Qinchuan cattle. ARIH2, DACT1 and DNM2, et al. are related to meat quality. Meanwhile, TBXA2R can be used as a gene associated with reproductive function, and USH2A affect coat color. This provided a glimpse into the formation of breeds and molecular genetic breeding. Our findings will promote genome-assisted breeding to improve animal production and health.

Genomic insights into selection for heterozygous alleles and woody traits in Populus tomentosa

Heterozygous alleles are widespread in outcrossing and clonally propagated woody plants. The variation in heterozygosity that underlies population adaptive evolution and phenotypic variation, however, remains largely unknown. Here, we describe a de novo chromosome-level genome assembly of Populus tomentosa, an economic and ecologically important native tree in northern China. By resequencing 302 natural accessions, we determined that the South subpopulation (Pop_S) encompasses the ancestral strains of P. tomentosa, while the Northwest subpopulation (Pop_NW) and Northeast subpopulation (Pop_NE) experienced different selection pressures during population evolution, resulting in significant population differentiation and a decrease in the extent of heterozygosity. Analysis of heterozygous selective sweep regions (HSSR) suggested that selection for lower heterozygosity contributed to the local adaptation of P. tomentosa by dwindling gene expression and genetic load in the Pop_NW and Pop_NE subpopulations. Genome-wide association studies (GWAS) revealed that 88 single nucleotide polymorphisms (SNPs) within 63 genes are associated with nine wood composition traits. Among them, the selection for the homozygous AA allele in PtoARF8 is associated with reductions in cellulose and hemicellulose contents by attenuating PtoARF8 expression, and the increase in lignin content is attributable to the selection for decreases in exon heterozygosity in PtoLOX3 during adaptive evolution of natural populations. This study provides novel insights into allelic variations in heterozygosity associated with adaptive evolution of P. tomentosa in response to the local environment and identifies a series of key genes for wood component traits, thereby facilitating genomic-based breeding of important traits in perennial woody plants.

Chromosome-level reference genome assembly provides insights into the evolution of Pennisetum alopecuroides

Pennisetum alopecuroides is an important forage grass resource, which plays a vital role in ecological environment improvement. Therefore, the acquisition of P. alopecuroides genome resources is conducive to the study of the adaptability of Pennisetum species in ecological remediation and forage breeding development. Here we assembled a P. alopecuroides cv. 'Liqiu' genome at the chromosome level with a size of approximately 845.71 Mb, contig N50 of 84.83Mb, and genome integrity of 99.13% as assessed by CEGMA. A total of 833.41-Mb sequences were mounted on nine chromosomes by Hi-C technology. In total, 60.66% of the repetitive sequences and 34,312 genes were predicted. The genomic evolution analysis showed that P. alopecuroides cv. 'Liqiu' was isolated from Setaria 7.53-13.80 million years ago and from Cenchrus 5.33-8.99 million years ago, respectively. The whole-genome event analysis showed that P. alopecuroides cv. 'Liqiu' underwent two whole-genome duplication (WGD) events in the evolution process, and the duplication events occurred at a similar time to that of Oryza sativa and Setaria viridis. The completion of the genome sequencing of P. alopecuroides cv. 'Liqiu' provides data support for mining high-quality genetic resources of P. alopecuroides and provides a theoretical basis for the origin and evolutionary characteristics of Pennisetum.

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.

Genome of tetraploid sour cherry (Prunus cerasus L.) 'Montmorency' identifies three distinct ancestral Prunus genomes

Sour cherry (Prunus cerasus L.) is a valuable fruit crop in the Rosaceae family and a hybrid between progenitors closely related to extant Prunus fruticosa (ground cherry) and Prunus avium (sweet cherry). Here we report a chromosome-scale genome assembly for sour cherry cultivar Montmorency, the predominant cultivar grown in the USA. We also generated a draft assembly of P. fruticosa to use alongside a published P. avium sequence for syntelog-based subgenome assignments for 'Montmorency' and provide compelling evidence P. fruticosa is also an allotetraploid. Using hierarchal k-mer clustering and phylogenomics, we show 'Montmorency' is trigenomic, containing two distinct subgenomes inherited from a P. fruticosa-like ancestor (A and A') and two copies of the same subgenome inherited from a P. avium-like ancestor (BB). The genome composition of 'Montmorency' is AA'BB and little-to-no recombination has occurred between progenitor subgenomes (A/A' and B). In Prunus, two known classes of genes are important to breeding strategies: the self-incompatibility loci (S-alleles), which determine compatible crosses, successful fertilization, and fruit set, and the Dormancy Associated MADS-box genes (DAMs), which strongly affect dormancy transitions and flowering time. The S-alleles and DAMs in 'Montmorency' and P. fruticosa were manually annotated and support subgenome assignments. Lastly, the hybridization event 'Montmorency' is descended from was estimated to have occurred less than 1.61 million years ago, making sour cherry a relatively recent allotetraploid. The 'Montmorency' genome highlights the evolutionary complexity of the genus Prunus and will inform future breeding strategies for sour cherry, comparative genomics in the Rosaceae, and questions regarding neopolyploidy.

Transcriptome Analysis of the Late-Acting Self-Incompatibility Associated with RNase T2 Family in Camellia oleifera

The Camellia oil tree (Camellia oleifera Abel.) is an important nonwood forest species in China, and the majority of its cultivars are late-acting self-incompatibility (LSI) types. Although several studies have examined the mechanism of LSI, the process is quite complicated and unclear. In this study, pollen tube growth and fruit setting of two Camellia oil tree cultivars Huashuo (HS) and Huajin (HJ) were investigated after non and self-pollination, and transcriptomic analysis of the ovaries was performed 48 h after self-pollination to identify the potential genes implicated in the LSI of Camellia oil trees. The results showed that the fruit set of HS was significantly higher than that of HJ after self-pollination. Transcriptomic analysis revealed that plant hormone signal transduction, the phosphatidylinositol signaling system, ATP-binding cassette (ABC) transporters, reactive oxygen species (ROS) metabolism, and Ca2+ signaling were mainly contributed in the LSI of reaction of Camellia oil tree. Moreover, nine RNase T2 genes were identified from the transcriptome analysis, which also showed that CoRNase7 participated in the self-incompatibility reaction in HS. Based on phylogenetic analysis, CoRNase6 was closely related to S-RNase from coffee, and CoRNase7 and CoRNase8 were closely related to S-RNase from Camellia sinensis. The 9 RNase T2 genes successfully produced proteins in prokaryotes. Subcellular localization indicated that CoRNase1 and CoRNase5 were cytoplasmic proteins, while CoRNase7 was a plasma membrane protein. These results screened the main metabolic pathways closely related to LSI in Camellia oil tree, and SI signal transduction might be regulated by a large molecular regulatory network. The discovery of T2 RNases provided evidence that Camellia oil tree might be under RNase-based gametophytic self-incompatibility.

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.

Deep reticulation: the long legacy of hybridization in vascular plant evolution

Hybridization has long been recognized as a fundamental evolutionary process in plants but, until recently, our understanding of its phylogenetic distribution and biological significance across deep evolutionary scales has been largely obscure. Over the past decade, genomic and phylogenomic datasets have revealed, perhaps not surprisingly, that hybridization, often associated with polyploidy, has been common throughout the evolutionary history of plants, particularly in various lineages of flowering plants. However, phylogenomic studies have also highlighted the challenges of disentangling signals of ancient hybridization from other sources of genomic conflict (in particular, incomplete lineage sorting). Here, we provide a critical review of ancient hybridization in vascular plants, outlining well-documented cases of ancient hybridization across plant phylogeny, as well as the challenges unique to documenting ancient versus recent hybridization. We provide a definition for ancient hybridization, which, to our knowledge, has not been explicitly attempted before. Further documenting the extent of deep reticulation in plants should remain an important research focus, especially because published examples likely represent the tip of the iceberg in terms of the total extent of ancient hybridization. However, future research should increasingly explore the macroevolutionary significance of this process, in terms of its impact on evolutionary trajectories (e.g. how does hybridization influence trait evolution or the generation of biodiversity over long time scales?), as well as how life history and ecological factors shape, or have shaped, the frequency of hybridization across geologic time and plant phylogeny. Finally, we consider the implications of ubiquitous ancient hybridization for how we conceptualize, analyze, and classify plant phylogeny. Networks, as opposed to bifurcating trees, represent more accurate representations of evolutionary history in many cases, although our ability to infer, visualize, and use networks for comparative analyses is highly limited. Developing improved methods for the generation, visualization, and use of networks represents a critical future direction for plant biology. Current classification systems also do not generally allow for the recognition of reticulate lineages, and our classifications themselves are largely based on evidence from the chloroplast genome. Updating plant classification to better reflect nuclear phylogenies, as well as considering whether and how to recognize hybridization in classification systems, will represent an important challenge for the plant systematics community.

Genome-Wide Identification and Gene Expression Analysis of Sweet Cherry Aquaporins (Prunus avium L.) under Abiotic Stresses

Aquaporins (AQPs) are integral transmembrane proteins well known as channels involved in the mobilization of water, small uncharged molecules and gases. In this work, the main objective was to carry out a comprehensive study of AQP encoding genes in Prunus avium (cv. Mazzard F12/1) on a genome-wide scale and describe their transcriptional behaviors in organs and in response to different abiotic stresses. A total of 28 non-redundant AQP genes were identified in Prunus spp. Genomes, which were phylogenetically grouped into five subfamilies (seven PIPs, eight NIPs, eight TIPs, three SIPs and two XIPs). Bioinformatic analyses revealed a high synteny and remarkable conservation of structural features among orthologs of different Prunus genomes. Several cis-acting regulatory elements (CREs) related to stress regulation were detected (ARE, WRE3, WUN, STRE, LTR, MBS, DRE, AT-rich and TC-rich). The above could be accounting for the expression variations associated with plant organs and, especially, each abiotic stress analyzed. Gene expressions of different PruavAQPs were shown to be preferentially associated with different stresses. PruavXIP2;1 and PruavXIP1;1 were up-regulated in roots at 6 h and 72 h of hypoxia, and in PruavXIP2;1 a slight induction of expression was also detected in leaves. Drought treatment strongly down-regulated PruavTIP4;1 but only in roots. Salt stress exhibited little or no variation in roots, except for PruavNIP4;1 and PruavNIP7;1, which showed remarkable gene repression and induction, respectively. Interestingly, PruavNIP4;1, the AQP most expressed in cherry roots subjected to cold temperatures, also showed this pattern in roots under high salinity. Similarly, PruavNIP4;2 consistently was up-regulated at 72 h of heat and drought treatments. From our evidence is possible to propose candidate genes for the development of molecular markers for selection processes in breeding programs for rootstocks and/or varieties of 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.

Spontaneous, Artificial, and Genome Editing-Mediated Mutations in Prunus

Mutation is a source of genetic diversity widely used in breeding programs for the acquisition of agronomically interesting characters in commercial varieties of the Prunus species, as well as in the rest of crop species. Mutation can occur in nature at a very low frequency or can be induced artificially. Spontaneous or bud sport mutations in somatic cells can be vegetatively propagated to get an individual with the mutant phenotype. Unlike animals, plants have unlimited growth and totipotent cells that let somatic mutations to be transmitted to the progeny. On the other hand, in vitro tissue culture makes it possible to induce mutation in plant material and perform large screenings for mutant’s selection and cleaning of chimeras. Finally, targeted mutagenesis has been boosted by the application of CRISPR (Clustered Regularly Interspaced Short Palindromic Repeats)/Cas9 and Transcription activator-like effector nuclease (TALEN) editing technologies. Over the last few decades, environmental stressors such as global warming have been threatening the supply of global demand for food based on population growth in the near future. For this purpose, the release of new varieties adapted to such changes is a requisite, and selected or generated Prunus mutants by properly regulated mechanisms could be helpful to this task. In this work, we reviewed the most relevant mutations for breeding traits in Prunus species such as flowering time, self-compatibility, fruit quality, and disease tolerance, including new molecular perspectives in the present postgenomic era including CRISPR/Cas9 and TALEN editing technologies.

Botany and national identities: The Tokyo Cherry

When Japan faced the world after the collapse of its feudal system, it had to invent its own modern identity in which the Tokyo Cherry became the National Flower. Despite being a garden plant, it received a Latin scientific species name as if it was an endemic species. After Japan's colonial conquest of Korea, exploring the flora of the peninsula became part of imperial knowledge practices of Japan. In the wild, a different cherry was discovered in Korea that was proposed as the endemic parent of the Tokyo Cherry, supporting imperialist policies. Following Japan's defeat after the Pacific War, South Korea in turn entered its search for cultural identity. The supposed parent of the Tokyo Cherry was now successfully acclaimed as the parent species of the colonial oppressor's Tokyo Cherry and named the King Cherry. Such scientific practice into cherries smoothly intertwined with nationalism and its legacy continues to interfere with research today.

The chromosome-level genome provides insight into the molecular mechanism underlying the tortuous-branch phenotype of Prunus mume

Plant with naturally twisted branches is referred to as a tortuous-branch plant, which have extremely high ornamental value due to their zigzag shape and the natural twisting of their branches. Prunus mume is an important woody ornamental plant. However, the molecular mechanism underlying this unique trait in Prunus genus is unknown. Here, we present a chromosome-level genome assembly of the cultivated P. mume var. tortuosa created using Oxford Nanopore combined with Hi-C scaffolding, which resulted in a 237.8 Mb genome assembly being anchored onto eight pseudochromosomes. Molecular dating indicated that P. mume is the most recently differentiated species in Prunus. Genes associated with cell division, development and plant hormones play essential roles in the formation of tortuous branch trait. A putative regulatory pathway for the tortuous branch trait was constructed based on gene expression levels. Furthermore, after transferring candidate PmCYCD genes into Arabidopsis thaliana, we found that seedlings overexpressing these genes exhibited curled rosette leaves. Our results provide insights into the evolutionary history of recently differentiated species in Prunus genus, the molecular basis of stem morphology, and the molecular mechanism underlying the tortuous branch trait and highlight the utility of multi-omics in deciphering the properties of P. mume plant architecture.

Genome Assembly and Analysis of the Flavonoid and Phenylpropanoid Biosynthetic Pathways in Fingerroot Ginger (Boesenbergia rotunda)

Boesenbergia rotunda (Zingiberaceae), is a high-value culinary and ethno-medicinal plant of Southeast Asia. The rhizomes of this herb have a high flavanone and chalcone content. Here we report the genome analysis of B. rotunda together with a complete genome sequence as a hybrid assembly. B. rotunda has an estimated genome size of 2.4 Gb which is assembled as 27,491 contigs with an N50 size of 12.386 Mb. The highly heterozygous genome encodes 71,072 protein-coding genes and has a 72% repeat content, with class I TEs occupying ~67% of the assembled genome. Fluorescence in situ hybridization of the 18 chromosome pairs at the metaphase showed six sites of 45S rDNA and two sites of 5S rDNA. An SSR analysis identified 238,441 gSSRs and 4604 EST-SSRs with 49 SSR markers common among related species. Genome-wide methylation percentages ranged from 73% CpG, 36% CHG and 34% CHH in the leaf to 53% CpG, 18% CHG and 25% CHH in the embryogenic callus. Panduratin A biosynthetic unigenes were most highly expressed in the watery callus. B rotunda has a relatively large genome with a high heterozygosity and TE content. This assembly and data (PRJNA71294) comprise a source for further research on the functional genomics of B. rotunda, the evolution of the ginger plant family and the potential genetic selection or improvement of gingers.

Chromosome-level genome assemblies of four wild peach species provide insights into genome evolution and genetic basis of stress resistance

Background

Peach (Prunus persica) is an economically important stone fruit crop in Rosaceae and widely cultivated in temperate and subtropical regions, emerging as an excellent material to study the interaction between plant and environment. During its genus, there are four wild species of peach, all living in harsh environments. For example, one of the wild species, P. mira, originates from the Qinghai-Tibet Plateau (QTP) and exhibits strong cold/ultraviolet ray environmental adaptations. Although remarkable progresses in the gene discovery of fruit quality-related traits in peach using previous assembled genome were obtained, genomic basis of the response of these wild species to different geographical environments remains unclear.

Results

To uncover key genes regulating adaptability in different species and analyze the role of genetic variations in resistance formation, we performed de novo genome assembling of four wild relatives of peach (P. persica), P. mira, P. davidiana, P. kansuensis, and P. ferganensis and resequenced 175 peach varieties. The phylogenetic tree showed that the divergence time of P. mira and other wild relatives of peach was 11.5 million years ago, which was consistent with the drastic crustal movement of QTP. Abundant genetic variations were identified in four wild species when compared to P. persica, and the results showed that plant-pathogen interaction pathways were enriched in genes containing small insertions and deletions and copy number variations in all four wild relatives of peach. Then, the data were used to identify new genes and variations regulating resistance. For example, presence/absence variations which result from a hybridization event that occurred between P. mira and P. dulcis enhanced the resistance of their putative hybrid, P. davidiana. Using bulked segregant analysis, we located the nematode resistance locus of P. kansuensis in chromosome 2. Within the mapping region, a deletion in the promoter of one NBS-LRR gene was found to involve the resistance by regulating gene expression. Furthermore, combined with RNA-seq and selective sweeps analysis, we proposed that a deletion in the promoter of one CBF gene was essential for high-altitude adaptation of P. mira through increasing its resistance to low temperature.

Conclusions

In general, the reference genomes assembled in the study facilitate our understanding of resistance mechanism of perennial fruit crops, and provide valuable resources for future breeding and improvement.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-022-01342-y.

Genome-wide identification of the SWEET gene family mediating the cold stress response in Prunus mume

The Sugars Will Eventually be Exported Transporter (SWEET) gene family encodes a family of sugar transporters that play essential roles in plant growth, reproduction, and biotic and abiotic stresses. Prunus mume is a considerable ornamental wood plant with high edible and medicinal values; however, its lack of tolerance to low temperature has severely limited its geographical distribution. To investigate whether this gene family mediates the response of P. mume to cold stress, we identified that the P. mume gene family consists of 17 members and divided the family members into four groups. Sixteen of these genes were anchored on six chromosomes, and one gene was anchored on the scaffold with four pairs of segmental gene duplications and two pairs of tandem gene duplications. Cis-acting regulatory element analysis indicated that the PmSWEET genes are potentially involved in P. mume development, including potentially regulating roles in procedure, such as circadian control, abscisic acid-response and light-response, and responses to numerous stresses, such as low-temperature and drought. We performed low-temperature treatment in the cold-tolerant cultivar 'Songchun' and cold-sensitive cultivar 'Zaolve' and found that the expression of four of 17 PmSWEETs was either upregulated or downregulated with prolonged treatment times. This finding indicates that these family members may potentially play a role in cold stress responses in P. mume. Our study provides a basis for further investigation of the role of SWEET proteins in the development of P. mume and its responses to cold stress.

An improved Raphanus sativus cv. WK10039 genome localizes centromeres, uncovers variation of DNA methylation and resolves arrangement of the ancestral Brassica genome blocks in radish chromosomes

This study presents an improved genome of Raphanus sativus cv. WK10039 uncovering centromeres and differentially methylated regions of radish chromosomes. Comprehensive genome comparison of radish and diploid Brassica species of U's triangle reveals that R. sativus arose from the Brassica B genome lineage and is a sibling species of B. nigra. Radish (Raphanus sativus L.) is a key root vegetable crop closely related to the Brassica crop species of the family Brassicaceae. We reported a draft genome of R. sativus cv. WK10039 (Rs1.0), which had 54.6 Mb gaps. To study the radish genome and explore previously unknown regions, we generated an improved genome assembly (Rs2.0) by long-read sequencing and high-resolution genome-wide mapping of chromatin interactions. Rs2.0 was 434.9 Mb in size with 0.27 Mb gaps, and the N50 scaffold length was 37.3 Mb (40-fold larger assembly compared to Rs1.0). Approximately 38% of Rs2.0 was comprised of repetitive sequences, and 52,768 protein-coding genes and 4845 non-protein-coding genes were predicted and annotated. The improved contiguity and coverage of Rs2.0, along with the detection of highly methylated regions, enabled localization of centromeres where R. sativus-specific centromere-associated repeats, full-length OTA and CRM LTR-Gypsy retrotransposons, hAT-Ac, CMC-EnSpm and Helitron DNA transposons, and sequences highly homologous to B. nigra centromere-specific CENH3-associated CL sequences were enriched. Whole-genome bisulfite sequencing combined with mRNA sequencing identified differential epigenetic marks in the radish genome related to tissue development. Synteny comparison and genomic distance analysis of radish and three diploid Brassica species of U's triangle suggested that the radish genome arose from the Brassica B genome lineage through unique rearrangement of the triplicated ancestral Brassica genome after splitting of the Brassica A/C and B genomes.

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.

The chromosome-scale genomes of Dipterocarpus turbinatus and Hopea hainanensis (Dipterocarpaceae) provide insights into fragrant oleoresin biosynthesis and hardwood formation

Dipterocarpaceae are typical tropical plants (dipterocarp forests) that are famous for their high economic value because of their production of fragrant oleoresins, top-quality timber and usage in traditional Chinese medicine. Currently, the lack of Dipterocarpaceae genomes has been a limiting factor to decipher the fragrant oleoresin biosynthesis and gain evolutionary insights into high-quality wood formation in Dipterocarpaceae. We generated chromosome-level genome assemblies for two representative Dipterocarpaceae species viz. Dipterocarpus turbinatus Gaertn. f. and Hopea hainanensis Merr. et Chun. Our whole-genome duplication (WGD) analysis revealed that Dipterocarpaceae underwent a shared WGD event, which showed significant impacts on increased copy numbers of genes related to the biosynthesis of terpene, BAHD acyltransferases, fatty acid and benzenoid/phenylpropanoid, which probably confer to the formation of their characteristic fragrant oleoresin. Additionally, compared with common soft wood plants, the expansion of gene families was also found to be associated with wood formation, such as in CESA (cellulose synthase), CSLE (cellulose synthase-like protein E), laccase and peroxidase in Dipterocarpaceae genomes, which might also contribute to the formation of harder, stronger and high-density timbers. Finally, an integrative analysis on a combination of genomic, transcriptomic and metabolic data from different tissues provided further insights into the molecular basis of fragrant oleoresins biosynthesis and high-quality wood formation of Dipterocarpaceae. Our study contributes the first two representative genomes for Dipterocarpaceae, which are valuable genetic resources for further researches on the fragrant oleoresins and superior-quality timber, genome-assisted breeding and improvement, and conservation biology of this family.

Cherry Blossom Forecast Based on Transcriptome of Floral Organs Approaching Blooming in the Flowering Cherry (Cerasus × yedoensis) Cultivar 'Somei-Yoshino'

To gain insights into the genetic mechanisms underlying blooming and petal movement in flowering cherry (Cerasus × yedoensis), we performed time-course RNA-seq analysis of the floral buds and open-flowers of the most popular flowering cherry cultivar, 'Somei-Yoshino.' Independent biological duplicate samples of floral buds and open-flowers were collected from 'Somei-Yoshino' trees grown at three different locations in Japan. RNA-seq reads obtained from floral bud and open-flower samples collected in the current study (in 2019) and in a previous study (in 2017) were aligned against the genome sequence of 'Somei-Yoshino' to quantify gene transcript levels. Clustering analysis of RNA-seq reads revealed dynamic changes in the transcriptome, with genes in seven modules predominantly expressed at specific time points, ranging from 5 weeks before flowering to 2 weeks after flowering. Based on the identified gene modules and Gene Ontology (GO) terms enriched at different floral stages, we speculate that the genetic mechanisms underlying petal movement and flower opening in cherry involve the processes of development, cell wall organization, reproduction, and metabolism, which are executed by genes encoding transcription factors, phytohormones, transporters, and polysaccharide metabolic enzymes. Furthermore, we established a statistical model for cherry bloom forecasting, based on gene expression levels as RNA markers at different time points before flowering.

Different scales of gene duplications occurring at different times have jointly shaped the NBS-LRR genes in Prunus species

In this study, genome-wide identification, phylogenetic relationships, duplication time and selective pressure of the NBS-LRR genes, an important group of plant disease-resistance genes (R genes), were performed to uncover their genetic evolutionary patterns in the six Prunus species. A total of 1946 NBS-LRR genes were identified; specifically, 589, 361, 284, 281, 318, and 113 were identified in Prunus yedoensis, P. domestica, P. avium, P. dulcis, P. persica and P. yedoensis var. nudiflora, respectively. Two NBS-LRR gene subclasses, TIR-NBS-LRR (TNL) and non-TIR-NBS-LRR (non-TNL), were also discovered. In total, 435 TNL and 1511 non-TNL genes were identified and could be classified into 30/55/75 and 103/158/191 multi-gene families, respectively, according to three different criteria. Higher Ks and Ka/Ks values were detected in TNL gene families than in non-TNL gene families. These results indicated that the TNL genes had more members involved in relatively ancient duplications and were affected by stronger selection pressure than the non-TNL genes. In general, the NBS-LRR genes were shaped by species-specific duplications, and lineage-specific duplications occurred at recent and relatively ancient periods among the six Prunus species. Therefore, different duplicated copies of NBS-LRRs can resist specific pathogens and will provide an R-gene library for resistance breeding in Prunus species.

Recent Duplications Dominate VQ and WRKY Gene Expansions in Six Prunus Species

Genes encoding VQ motif-containing (VQ) transcriptional regulators and WRKY transcription factors can participate separately or jointly in plant growth, development, and abiotic and biotic stress responses. In this study, 222 VQ and 645 WRKY genes were identified in six Prunus species. Based on phylogenetic tree topologies, the VQ and WRKY genes were classified into 13 and 32 clades, respectively. Therefore, at least 13 VQ gene copies and 32 WRKY gene copies were present in the genome of the common ancestor of the six Prunus species. Similar small Ks value peaks for the VQ and WRKY genes suggest that the two gene families underwent recent duplications in the six studied species. The majority of the Ka/Ks ratios were less than 1, implying that most of the VQ and WRKY genes had undergone purifying selection. Pi values were significantly higher in the VQ genes than in the WRKY genes, and the VQ genes therefore exhibited greater nucleotide diversity in the six species. Forty-one of the Prunus VQ genes were predicted to interact with 44 of the WRKY genes, and the expression levels of some predicted VQ-WRKY interacting pairs were significantly correlated. Differential expression patterns of the VQ and WRKY genes suggested that some might be involved in regulating aphid resistance in P. persica and fruit development in P. avium.

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.

Resequencing of 388 cassava accessions identifies valuable loci and selection for variation in heterozygosity

Background

Heterozygous genomes are widespread in outcrossing and clonally propagated crops. However, the variation in heterozygosity underlying key agronomic traits and crop domestication remains largely unknown. Cassava is a staple crop in Africa and other tropical regions and has a highly heterozygous genome.

Results

We describe a genomic variation map from 388 resequenced genomes of cassava cultivars and wild accessions. We identify 52 loci for 23 agronomic traits through a genome-wide association study. Eighteen allelic variations in heterozygosity for nine candidate genes are significantly associated with seven key agronomic traits. We detect 81 selective sweeps with decreasing heterozygosity and nucleotide diversity, harboring 548 genes, which are enriched in multiple biological processes including growth, development, hormone metabolisms and responses, and immune-related processes. Artificial selection for decreased heterozygosity has contributed to the domestication of the large starchy storage root of cassava. Selection for homozygous GG allele in MeTIR1 during domestication contributes to increased starch content. Selection of homozygous AA allele in MeAHL17 is associated with increased storage root weight and cassava bacterial blight (CBB) susceptibility. We have verified the positive roles of MeTIR1 in increasing starch content and MeAHL17 in resistance to CBB by transient overexpression and silencing analysis. The allelic combinations in MeTIR1 and MeAHL17 may result in high starch content and resistance to CBB.

Conclusions

This study provides insights into allelic variation in heterozygosity associated with key agronomic traits and cassava domestication. It also offers valuable resources for the improvement of cassava and other highly heterozygous crops.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13059-021-02524-7.

The draft chromosome-level genome assembly of tetraploid ground cherry (Prunus fruticosa Pall.) from long reads

Cherries are stone fruits and belong to the economically important plant family of Rosaceae with worldwide cultivation of different species. The ground cherry, Prunus fruticosa Pall., is an ancestor of cultivated sour cherry, an important tetraploid cherry species. Here, we present a long read chromosome-level draft genome assembly and related plastid sequences using the Oxford Nanopore Technology PromethION platform and R10.3 pore type. We generated a final consensus genome sequence of 366 Mb comprising eight chromosomes. The N50 scaffold was ~44 Mb with the longest chromosome being 66.5 Mb. The chloroplast and mitochondrial genomes were 158,217 bp and 383,281 bp long, which is in accordance with previously published plastid sequences. This is the first report of the genome of ground cherry (P. fruticosa) sequenced by long read technology only. The datasets obtained from this study provide a foundation for future breeding, molecular and evolutionary analysis in Prunus studies.

A draft genome, resequencing, and metabolomes reveal the genetic background and molecular basis of the nutritional and medicinal properties of loquat (Eriobotrya japonica (Thunb.) Lindl)

Loquat (Eriobotrya japonica) is a popular fruit and medicinal plant. Here, a high-quality draft genome of the E. japonica 'Big Five-pointed Star' cultivar that covers ~98% (733.32 Mb) of the estimated genome size (749.25 Mb) and contains a total of 45,492 protein-coding genes is reported. Comparative genomic analysis suggests that the loquat genome has evolved a unique genetic mechanism of chromosome repair. Resequencing data from 52 loquat cultivars, including 16 white-fleshed and 36 yellow-fleshed variants, were analyzed, and the flower, leaf, and root metabolomes of 'Big Five-pointed Star' were determined using a UPLC-ESI-MS/M system. A genome-wide association study identified several candidate genes associated with flesh color in E. japonica, linking these phenotypes to sugar metabolism. A total of 577 metabolites, including 98 phenolic acids, 95 flavonoids, and 28 terpenoids, were found, and 191 metabolites, including 46 phenolic acids, 33 flavonoids, and 7 terpenoids, showed no differences in concentration among the leaves, roots, and flowers. Candidate genes related to the biosynthesis of various medicinal ingredients, such as phenolics, flavonoids, terpenoids, and polysaccharides, were identified. Some of these genes were confirmed to be members of expanding gene families, suggesting that the high concentrations of beneficial metabolites in loquat may be associated with the number of biosynthetic genes in this plant. In summary, this study provides fundamental molecular insights into the nutritional and medical properties of E. japonica.

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.

Frequent germplasm exchanges drive the high genetic diversity of Chinese-cultivated common apricot germplasm

The genetic diversity of germplasm is critical for exploring genetic and phenotypic resources and has important implications for crop-breeding sustainability and improvement. However, little is known about the factors that shape and maintain genetic diversity. Here, we assembled a high-quality chromosome-level reference of the Chinese common apricot 'Yinxiangbai', and we resequenced 180 apricot accessions that cover four major ecogeographical groups in China and other accessions from occidental countries. We concluded that Chinese-cultivated common apricot germplasms possessed much higher genetic diversity than those cultivated in Western countries. We also detected seven migration events among different apricot groups, where 27% of the genome was identified as being introgressed. Remarkably, we demonstrated that these introgressed regions drove the current high level of germplasm diversity in Chinese-cultivated common apricots by introducing different genes related to distinct phenotypes from different cultivated groups. Our results highlight the consideration that introgressed regions may provide an important reservoir of genetic resources that can be used to sustain modern breeding programs.

The chromosome-scale reference genome of Rubus chingii Hu provides insight into the biosynthetic pathway of hydrolyzable tannins

Rubus chingii Hu (Fu-Pen-Zi), a perennial woody plant in the Rosaceae family, is a characteristic traditional Chinese medicinal plant because of its unique pharmacological effects. There are abundant hydrolyzable tannin (HT) components in R. chingii that provide health benefits. Here, an R. chingii chromosome-scale genome and related functional analysis provide insights into the biosynthetic pathway of HTs. In total, sequence data of 231.21 Mb (155 scaffolds with an N50 of 8.2 Mb) were assembled into seven chromosomes with an average length of 31.4 Mb, and 33 130 protein-coding genes were predicted, 89.28% of which were functionally annotated. Evolutionary analysis showed that R. chingii was most closely related to Rubus occidentalis, from which it was predicted to have diverged 22.46 million years ago (Table S8). Comparative genomic analysis showed that there was a tandem gene cluster of UGT, carboxylesterase (CXE) and SCPL genes on chromosome 02 of R. chingii, including 11 CXE, eight UGT, and six SCPL genes, which may be critical for the synthesis of HTs. In vitro enzyme assays indicated that the proteins encoded by the CXE (LG02.4273) and UGT (LG02.4102) genes have tannin hydrolase and gallic acid glycosyltransferase functions, respectively. The genomic sequence of R. chingii will be a valuable resource for comparative genomic analysis within the Rosaceae family and will be useful for understanding the biosynthesis of HTs.

The first draft genome of Picrorhiza kurrooa, an endangered medicinal herb from Himalayas

Picrorhiza kurrooa is an endangered medicinal herb which is distributed across the Himalayan region at an altitude between 3000–5000 m above mean sea level. The medicinal properties of P. kurrooa are attributed to monoterpenoid picrosides present in leaf, rhizome and root of the plant. However, no genomic information is currently available for P. kurrooa, which limits our understanding about its molecular systems and associated responses. The present study brings the first assembled draft genome of P. kurrooa by using 227 Gb of raw data generated by Illumina and PacBio RS II sequencing platforms. The assembled genome has a size of n = ~ 1.7 Gb with 12,924 scaffolds. Four pronged assembly quality validations studies, including experimentally reported ESTs mapping and directed sequencing of the assembled contigs, confirmed high reliability of the assembly. About 76% of the genome is covered by complex repeats alone. Annotation revealed 24,798 protein coding and 9789 non-coding genes. Using the assembled genome, a total of 710 miRNAs were discovered, many of which were found responsible for molecular response against temperature changes. The miRNAs and targets were validated experimentally. The availability of draft genome sequence will aid in genetic improvement and conservation of P. kurrooa. Also, this study provided an efficient approach for assembling complex genomes while dealing with repeats when regular assemblers failed to progress due to repeats.

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.

An Atlas of Genomic Resources for Studying Rosaceae Fruits and Ornamentals

Rosaceae, a large plant family of more than 3,000 species, consists of many economically important fruit and ornamental crops, including peach, apple, strawberry, raspberry, cherry, and rose. These horticultural crops are not only important economic drivers in many regions of the world, but also major sources of human nutrition. Additionally, due to the diversity of fruit types in Rosaceae, this plant family offers excellent opportunities for investigations into fleshy fruit diversity, evolution, and development. With the development of high-throughput sequencing technologies and computational tools, an increasing number of high-quality genomes and transcriptomes of Rosaceae species have become available and will greatly facilitate Rosaceae research and breeding. This review summarizes major genomic resources and genome research progress in Rosaceae, highlights important databases, and suggests areas for further improvement. The availability of these big data resources will greatly accelerate research progress and enhance the agricultural productivity of Rosaceae.

The Draft Genome of a Flat Peach (Prunus persica L. cv. '124 Pan') Provides Insights into Its Good Fruit Flavor Traits

The flat peach has become more and more popular worldwide for its fruit quality with relatively low acidity, high sugar content and rich flavor. However, the draft genome assembly of flat peach is still unavailable and the genetic basis for its fruit flavor remains unclear. In this study, the draft genome of a flat peach cultivar '124 Pan' was assembled by using a hybrid assembly algorithm. The final assembly resulted in a total size of 206 Mb with a N50 of 26.3 Mb containing eight chromosomes and seven scaffolds. Genome annotation revealed that a total of 25,233 protein-coding genes were predicted with comparable gene abundance among the sequenced peach species. The phylogenetic tree and divergence times inferred from 572 single copy genes of 13 plant species confirmed that Prunus ferganensis was the ancestor of the domesticated peach. By comparing with the genomes of Prunus persica (Lovell) and Prunus ferganensis, the expansion of genes encoding enzymes involved in terpene biosynthesis was found, which might contribute to the good fruit flavor traits of '124 Pan'. The flat peach draft genome assembly obtained in this study will provide a valuable genomic resource for peach improvement and molecular breeding.

Genome-wide characterization of PEBP family genes in nine Rosaceae tree species and their expression analysis in P. mume

BACKGROUND

Phosphatidylethanolamine-binding proteins (PEBPs) constitute a common gene family found among animals, plants and microbes. Plant PEBP proteins play an important role in regulating flowering time, plant architecture as well as seed dormancy. Though PEBP family genes have been well studied in Arabidopsis and other model species, less is known about these genes in perennial trees.

RESULTS

To understand the evolution of PEBP genes and their functional roles in flowering control, we identified 56 PEBP members belonging to three gene clades (MFT-like, FT-like, and TFL1-like) and five lineages (FT, BFT, CEN, TFL1, and MFT) across nine Rosaceae perennial species. Structural analysis revealed highly conserved gene structure and protein motifs among Rosaceae PEBP proteins. Codon usage analysis showed slightly biased codon usage across five gene lineages. With selection pressure analysis, we detected strong purifying selection constraining divergence within most lineages, while positive selection driving the divergence of FT-like and TFL1-like genes from the MFT-like gene clade. Spatial and temporal expression analyses revealed the essential role of FT in regulating floral bud breaking and blooming in P. mume. By employing a weighted gene co-expression network approach, we inferred a putative FT regulatory module required for dormancy release and blooming in P. mume.

CONCLUSIONS

We have characterized the PEBP family genes in nine Rosaceae species and examined their phylogeny, genomic syntenic relationship, duplication pattern, and expression profiles during flowering process. These results revealed the evolutionary history of PEBP genes and their functions in regulating floral bud development and blooming among Rosaceae tree species.

Pilot RNA-seq data from 24 species of vascular plants at Harvard Forest

PREMISE

Large-scale projects such as the National Ecological Observatory Network (NEON) collect ecological data on entire biomes to track climate change. NEON provides an opportunity to launch community transcriptomic projects that ask integrative questions in ecology and evolution. We conducted a pilot study to investigate the challenges of collecting RNA-seq data from diverse plant communities.

METHODS

We generated >650 Gbp of RNA-seq for 24 vascular plant species representing 12 genera and nine families at the Harvard Forest NEON site. Each species was sampled twice in 2016 (July and August). We assessed transcriptome quality and content with TransRate, BUSCO, and Gene Ontology annotations.

RESULTS

Only modest differences in assembly quality were observed across multiple k-mers. On average, transcriptomes contained hits to >70% of loci in the BUSCO database. We found no significant difference in the number of assembled and annotated transcripts between diploid and polyploid transcriptomes.

DISCUSSION

We provide new RNA-seq data sets for 24 species of vascular plants in Harvard Forest. Challenges associated with this type of study included recovery of high-quality RNA from diverse species and access to NEON sites for genomic sampling. Overcoming these challenges offers opportunities for large-scale studies at the intersection of ecology and genomics.

Characteristics and Expression Analyses of Trehalose-6-Phosphate Synthase Family in Prunus mume Reveal Genes Involved in Trehalose Biosynthesis and Drought Response

Trehalose and its key synthase (trehalose-6-phosphate synthase, TPS) can improve the drought tolerance of plants. However, little is known about the roles of trehalose and the TPS family in Prunus mume response to drought. In our study, we discovered that the trehalose content in leaf, root, and stem tissues significantly increased in P. mume in response to drought. Therefore, the characteristics and functions of the TPS family are worth investigating in P. mume. We identified nine TPS family members in P. mume, which were divided into two sub-families and characterized by gene structure, promoter elements, protein conserved domains, and protein motifs. We found that the Hydrolase_3 domain and several motifs were highly conserved in Group II instead of Group I. The distinctions between the two groups may result from selective constraints, which we estimated by the d_N/d_S (ω) ratio. The ω values of all the PmTPS family gene pairs were evaluated as less than 1, indicating that purity selection facilitated their divergence. A phylogenetic tree was constructed using 92 TPSs from 10 Rosaceae species, which were further divided into five clusters. Based on evolutionary analyses, the five clusters of TPS family proteins mainly underwent varied purity selection. The expression patterns of PmTPSs under drought suggested that the TPS family played an important role in the drought tolerance of P. mume. Combining the expression patterns of PmTPSs and the trehalose content changes in leaf, stem, and root tissues under normal conditions and drought stress, we found that the PmTPS2 and PmTPS6 mainly function in the trehalose biosynthesis in P. mume. Our findings not only provide valuable information about the functions of trehalose and TPSs in the drought response of P. mume, but they also contribute to the future drought breeding of P. mume.

Chromosome-scale genome assembly of sweet cherry (Prunus avium L.) cv. Tieton obtained using long-read and Hi-C sequencing

Sweet cherry (Prunus avium) is an economically significant fruit species in the genus Prunus. However, in contrast to other important fruit trees in this genus, only one draft genome assembly is available for sweet cherry, which was assembled using only Illumina short-read sequences. The incompleteness and low quality of the current sweet cherry draft genome limit its use in genetic and genomic studies. A high-quality chromosome-scale sweet cherry reference genome assembly is therefore needed. A total of 65.05 Gb of Oxford Nanopore long reads and 46.24 Gb of Illumina short reads were generated, representing ~190x and 136x coverage, respectively, of the sweet cherry genome. The final de novo assembly resulted in a phased haplotype assembly of 344.29 Mb with a contig N50 of 3.25 Mb. Hi-C scaffolding of the genome resulted in eight pseudochromosomes containing 99.59% of the bases in the assembled genome. Genome annotation revealed that more than half of the genome (59.40%) was composed of repetitive sequences, and 40,338 protein-coding genes were predicted, 75.40% of which were functionally annotated. With the chromosome-scale assembly, we revealed that gene duplication events contributed to the expansion of gene families for salicylic acid/jasmonic acid carboxyl methyltransferase and ankyrin repeat-containing proteins in the genome of sweet cherry. Four auxin-responsive genes (two GH3s and two SAURs) were induced in the late stage of fruit development, indicating that auxin is crucial for the sweet cherry ripening process. In addition, 772 resistance genes were identified and functionally predicted in the sweet cherry genome. The high-quality genome assembly of sweet cherry obtained in this study will provide valuable genomic resources for sweet cherry improvement and molecular breeding.

A de novo assembly of the sweet cherry (Prunus avium cv. Tieton) genome using linked-read sequencing technology

The sweet cherry (Prunus avium) is one of the most economically important fruit species in the world. However, there is a limited amount of genetic information available for this species, which hinders breeding efforts at a molecular level. We were able to describe a high-quality reference genome assembly and annotation of the diploid sweet cherry (2n = 2x = 16) cv. Tieton using linked-read sequencing technology. We generated over 750 million clean reads, representing 112.63 GB of raw sequencing data. The Supernova assembler produced a more highly-ordered and continuous genome sequence than the current P. avium draft genome, with a contig N50 of 63.65 KB and a scaffold N50 of 2.48 MB. The final scaffold assembly was 280.33 MB in length, representing 82.12% of the estimated Tieton genome. Eight chromosome-scale pseudomolecules were constructed, completing a 214 MB sequence of the final scaffold assembly. De novo, homology-based, and RNA-seq methods were used together to predict 30,975 protein-coding loci. 98.39% of core eukaryotic genes and 97.43% of single copy orthologues were identified in the embryo plant, indicating the completeness of the assembly. Linked-read sequencing technology was effective in constructing a high-quality reference genome of the sweet cherry, which will benefit the molecular breeding and cultivar identification in this species.

Developmental Mechanisms of Fleshy Fruit Diversity in Rosaceae

Rosaceae (the rose family) is an economically important family that includes species prized for high-value fruits and ornamentals. The family also exhibits diverse fruit types, including drupe (peach), pome (apple), drupetum (raspberry), and achenetum (strawberry). Phylogenetic analysis and ancestral fruit-type reconstruction suggest independent evolutionary paths of multiple fleshy fruit types from dry fruits. A recent whole genome duplication in the Maleae/Pyreae tribe (with apple, pear, hawthorn, and close relatives; referred to as Maleae here) may have contributed to the evolution of pome fruit. MADS-box genes, known to regulate floral organ identity, are emerging as important regulators of fruit development. The differential competence of floral organs to respond to fertilization signals may explain the different abilities of floral organs to form fleshy fruit. Future comparative genomics and functional studies in closely related Rosaceae species with distinct fruit types will test hypotheses and provide insights into mechanisms of fleshy fruit diversity. These efforts will be facilitated by the wealth of genome data and resources in Rosaceae.