Evolutionary history and pan-genome dynamics of strawberry (Fragaria spp.)

Uncovering fruit flavor and genetic diversity across diploid wild Fragaria species via comparative metabolomics profiling

Wild Fragaria resources exhibit extensive genetic diversity and desirable edible traits, such as high soluble solid content and flavor compounds. However, specific metabolites in different wild strawberry fruits remain unknown. In this study, we characterized 1008 metabolites covering 11 subclasses among 13 wild diploid resources representing eight species, including F. vesca, F. nilgerrensis, F. viridis, F. nubicola, F. pentaphylla, F. mandschurica, F. chinensis, and F. emeiensis. Fifteen potential metabolite biomarkers were identified to distinguish fruit flavors among the 13 diploid wild Fragaria accessions. A total of nine distinct modules were employed to explore key metabolites related to fruit quality through weighted gene co-expression module analysis, with significant enrichment in amino acid biosynthesis pathway. Notably, the identified significantly different key metabolites highlighted the close association of amino acids, sugars, and anthocyanins with flavor formation. These findings offer valuable resources for improving fruit quality through metabolome-assisted breeding.

Graph-Based Pangenome of Actinidia chinensis Reveals Structural Variations Mediating Fruit Degreening

Fruit ripening is associated with the degreening process (loss of chlorophyll) that occurs in most fruit species. Kiwifruit is one of the special species whose fruits may maintain green flesh by accumulating a large amount of chlorophyll even after ripening. However, little is known about the genetic variations related to the fruit degreening process. Here, a graph-based kiwifruit pangenome by analyzing 14 chromosome-scale haplotype-resolved genome assemblies from seven representative cultivars or lines in Actinidia chinensis is built. A total of 49,770 non-redundant gene families are identified, with core genes constituting 46.6%, and dispensable genes constituting 53.4%. A total of 84,591 non-redundant structural variations (SVs) are identified. The pangenome graph integrating both reference genome sequences and variant information facilitates the identification of SVs related to fruit color. The SV in the promoter of the AcBCM gene determines its high expression in the late developmental stage of fruits, which causes chlorophyll accumulation in the green-flesh fruits by post-translationally regulating AcSGR2, a key enzyme of chlorophyll catabolism. Taken together, a high-quality pangenome is constructed, unraveled numerous genetic variations, and identified a novel SV mediating fruit coloration and fruit quality, providing valuable information for further investigating genome evolution and domestication, QTL genes function, and genomics-assisted breeding.

Chromosome-level genome assembly of a cliff plant Taihangia rupestris var. ciliata provides insights into its adaptation and demographic history

BACKGROUND

Cliffs are recognized as one of the most challenging environments for plants, characterized by harsh conditions such as drought, infertile soil, and steep terrain. However, they surprisingly host ancient and diverse plant communities and play a crucial role in protecting biodiversity. The Taihang Mountains, which act as a natural boundary in eastern China, support a rich variety of plant species, including many unique to cliff habitats. However, it is little known how cliff plants adapt to harsh habitats and the demographic history in this region.

RESULTS

To better understand the demographic history and adaptation of cliff plants in this area, we analyzed the chromosome-level genome of a representative cliff plant, T. rupestris var. ciliata, which has a genome size of 769.5 Mb, with a scaffold N50 of 104.92 Mb. The rapid expansion of transposable elements may have contributed to the increasing genome and its ability to adapt to unique and challenging cliff habitats. Comparative analysis of the genome evolution between Taihangia and non-cliff plants in Rosaceae revealed a significant expansion of gene families associated with oxidative phosphorylation, which is likely a response to the abiotic stresses faced by cliff plants. This expansion may explain the long-term adaptation of Taihangia to harsh cliff environments. The effective population size of the two varieties has continuously decreased due to climatic fluctuations during the Quaternary period. Furthermore, significant differences in gene expression between the two varieties may explain the varied leaf phenotypes and adaptations to harsh conditions in different natural distributions.

CONCLUSION

Our study highlights the extraordinary adaptation of T. rupestris var. ciliata, shedding light on the evolution of cliff plants worldwide.

Genomic analyses provide insights into sex differentiation of tetraploid strawberry (Fragaria moupinensis)

The strawberry genus, Fragaria, exhibits a wide range of sexual systems and natural ploidy variation. Nearly, all polyploid strawberry species exhibit separate sexes (dioecy). Research has identified the sex-determining sequences as roughly conserved but with repeatedly changed genomic locations across octoploid strawberries. However, it remains unclear whether tetraploid wild strawberries evolved dioecy independently or shared a common origin with octoploid strawberries. In this study, we investigated the sex determinants of F. moupinensis, a dioecious plant with heterogametic females (ZW). Utilizing a combination of haplotype-resolved genome sequencing of the female F. moupinensis, k-mer-based and coverage-based genome-wide association studies (GWAS), and transcriptomic analysis, we discovered a non-recombining, approximately 33.6 kb W-specific region on chromosome 2a. Within this region, only one candidate sex-determining gene (FmoAFT) was identified. Furthermore, an extensive resequencing of the entire Fragaria genus indicated that the W-specific region displays conservative female specificity across all tetraploid species. This observation suggests that dioecy evolved independently in tetraploid and octoploid strawberries. Moreover, employing virus-induced gene silencing (VIGS), we knocked down the expression of the FmoAFT homologue transcript in cultivated strawberries, revealing its potential role in promoting female functions during early carpel development. We also applied DNA affinity purification sequencing (DAP-seq) and yeast one-hybrid assays to identify potential direct targets of FmoAFT. These insights shed new light on the genetic basis and evolutionary history of sex determination in strawberries, thereby facilitating the formulation of strategies to manipulate sex determination in breeding programs.

Novel Allelic Gene Variations in CmCLAVATA3 (CmCLV3) Were Identified in a Genetic Population of Melon (Cucumis melo L.)

Carpel number (CN) is an important trait affecting the fruit size and shape of melon, which plays a crucial role in determining the overall appearance and market value. A unique non-synonymous single nucleotide polymorphism (SNP) in CmCLAVATA3 (CmCLV3) is responsible for the variation of CN in C. melo ssp. agrestis (hereafter agrestis), but it has been unclear in C. melo ssp. melo (hereafter melo). In this study, one major locus controlling the polymorphism of 5-CN (multi-CN) and 3-CN (normal-CN) in melo was identified using bulked segregant analysis (BSA-seq). This locus was then fine-mapped to an interval of 1.8 Mb on chromosome 12 using a segregating population containing 1451 progeny. CmCLV3 is still present in the candidate region. A new allele of CmCLV3, which contains five other nucleotide polymorphisms, including a non-synonymous SNP in coding sequence (CDS), except the SNP reported in agrestis, was identified in melo. A cis-trans test confirmed that the candidate gene, CmCLV3, contributes to the variation of CNs in melo. The qRT-PCR results indicate that there is no significant difference in the expression level of CmCLV3 in the apical stem between the multi-CN plants and the normal-CN plants. Overall, this study provides a genetic resource for melon fruit development research and molecular breeding. Additionally, it suggests that melo has undergone similar genetic selection but evolved into an independent allele.

Genomic signatures of strawberry domestication and diversification

Cultivated strawberry (Fragaria × ananassa) has a brief history of less than 300 yr, beginning with the hybridization of octoploids Fragaria chiloensis and Fragaria virginiana. Here we explored the genomic signatures of early domestication and subsequent diversification for different climates using whole-genome sequences of 289 wild, heirloom, and modern varieties from two major breeding programs in the United States. Four nonadmixed wild octoploid populations were identified, with recurrent introgression among the sympatric populations. The proportion of F. virginiana ancestry increased by 20% in modern varieties over initial hybrids, and the proportion of F. chiloensis subsp. pacifica rose from 0% to 3.4%. Effective population size rapidly declined during early breeding. Meanwhile, divergent selection for distinct environments reshaped wild allelic origins in 21 out of 28 chromosomes. Overlapping divergent selective sweeps in natural and domesticated populations revealed 16 convergent genomic signatures that may be important for climatic adaptation. Despite 20 breeding cycles since initial hybridization, more than half of loci underlying yield and fruit size are still not under artificial selection. These insights add clarity to the domestication and breeding history of what is now the most widely cultivated fruit in the world.

Advances in genomics and genome editing for improving strawberry (Fragaria ×ananassa)

The cultivated strawberry, Fragaria ×ananassa, is a recently domesticated fruit species of economic interest worldwide. As such, there is significant interest in continuous varietal improvement. Genomics-assisted improvement, including the use of DNA markers and genomic selection have facilitated significant improvements of numerous key traits during strawberry breeding. CRISPR/Cas-mediated genome editing allows targeted mutations and precision nucleotide substitutions in the target genome, revolutionizing functional genomics and crop improvement. Genome editing is beginning to gain traction in the more challenging polyploid crops, including allo-octoploid strawberry. The release of high-quality reference genomes and comprehensive subgenome-specific genotyping and gene expression profiling data in octoploid strawberry will lead to a surge in trait discovery and modification by using CRISPR/Cas. Genome editing has already been successfully applied for modification of several strawberry genes, including anthocyanin content, fruit firmness and tolerance to post-harvest disease. However, reports on many other important breeding characteristics associated with fruit quality and production are still lacking, indicating a need for streamlined genome editing approaches and tools in Fragaria ×ananassa. In this review, we present an overview of the latest advancements in knowledge and breeding efforts involving CRISPR/Cas genome editing for the enhancement of strawberry varieties. Furthermore, we explore potential applications of this technology for improving other Rosaceous plant species.

Nuclear phylogenomics of angiosperms and insights into their relationships and evolution

Angiosperms (flowering plants) are by far the most diverse land plant group with over 300,000 species. The sudden appearance of diverse angiosperms in the fossil record was referred to by Darwin as the "abominable mystery," hence contributing to the heightened interest in angiosperm evolution. Angiosperms display wide ranges of morphological, physiological, and ecological characters, some of which have probably influenced their species richness. The evolutionary analyses of these characteristics help to address questions of angiosperm diversification and require well resolved phylogeny. Following the great successes of phylogenetic analyses using plastid sequences, dozens to thousands of nuclear genes from next-generation sequencing have been used in angiosperm phylogenomic analyses, providing well resolved phylogenies and new insights into the evolution of angiosperms. In this review we focus on recent nuclear phylogenomic analyses of large angiosperm clades, orders, families, and subdivisions of some families and provide a summarized Nuclear Phylogenetic Tree of Angiosperm Families. The newly established nuclear phylogenetic relationships are highlighted and compared with previous phylogenetic results. The sequenced genomes of Amborella, Nymphaea, Chloranthus, Ceratophyllum, and species of monocots, Magnoliids, and basal eudicots, have facilitated the phylogenomics of relationships among five major angiosperms clades. All but one of the 64 angiosperm orders were included in nuclear phylogenomics with well resolved relationships except the placements of several orders. Most families have been included with robust and highly supported placements, especially for relationships within several large and important orders and families. Additionally, we examine the divergence time estimation and biogeographic analyses of angiosperm on the basis of the nuclear phylogenomic frameworks and discuss the differences compared with previous analyses. Furthermore, we discuss the implications of nuclear phylogenomic analyses on ancestral reconstruction of morphological, physiological, and ecological characters of angiosperm groups, limitations of current nuclear phylogenomic studies, and the taxa that require future attention.

Maternal dominance contributes to subgenome differentiation in allopolyploid fishes

Teleost fishes, which are the largest and most diverse group of living vertebrates, have a rich history of ancient and recent polyploidy. Previous studies of allotetraploid common carp and goldfish (cyprinids) reported a dominant subgenome, which is more expressed and exhibits biased gene retention. However, the underlying mechanisms contributing to observed 'subgenome dominance' remains poorly understood. Here we report high-quality genomes of twenty-one cyprinids to investigate the origin and subsequent subgenome evolution patterns following three independent allopolyploidy events. We identify the closest extant relatives of the diploid progenitor species, investigate genetic and epigenetic differences among subgenomes, and conclude that observed subgenome dominance patterns are likely due to a combination of maternal dominance and transposable element densities in each polyploid. These findings provide an important foundation to understanding subgenome dominance patterns observed in teleost fishes, and ultimately the role of polyploidy in contributing to evolutionary innovations.

Comparative Analysis of Transposable Elements in Strawberry Genomes of Different Ploidy Levels

Transposable elements (TEs) make up a large portion of plant genomes and play a vital role in genome structure, function, and evolution. Cultivated strawberry (Fragaria x ananassa) is one of the most important fruit crops, and its octoploid genome was formed through several rounds of genome duplications from diploid ancestors. Here, we built a pan-genome TE library for the Fragaria genus using ten published strawberry genomes at different ploidy levels, including seven diploids, one tetraploid, and two octoploids, and performed comparative analysis of TE content in these genomes. The TEs comprise 51.83% (F. viridis) to 60.07% (F. nilgerrensis) of the genomes. Long terminal repeat retrotransposons (LTR-RTs) are the predominant TE type in the Fragaria genomes (20.16% to 34.94%), particularly in F. iinumae (34.94%). Estimating TE content and LTR-RT insertion times revealed that species-specific TEs have shaped each strawberry genome. Additionally, the copy number of different LTR-RT families inserted in the last one million years reflects the genetic distance between Fragaria species. Comparing cultivated strawberry subgenomes to extant diploid ancestors showed that F. vesca and F. iinumae are likely the diploid ancestors of the cultivated strawberry, but not F. viridis. These findings provide new insights into the TE variations in the strawberry genomes and their roles in strawberry genome evolution.

Blueberry and cranberry pangenomes as a resource for future genetic studies and breeding efforts

Domestication of cranberry and blueberry began in the United States in the early 1800s and 1900s, respectively, and in part owing to their flavors and health-promoting benefits are now cultivated and consumed worldwide. The industry continues to face a wide variety of production challenges (e.g. disease pressures), as well as a demand for higher-yielding cultivars with improved fruit quality characteristics. Unfortunately, molecular tools to help guide breeding efforts for these species have been relatively limited compared with those for other high-value crops. Here, we describe the construction and analysis of the first pangenome for both blueberry and cranberry. Our analysis of these pangenomes revealed both crops exhibit great genetic diversity, including the presence-absence variation of 48.4% genes in highbush blueberry and 47.0% genes in cranberry. Auxiliary genes, those not shared by all cultivars, are significantly enriched with molecular functions associated with disease resistance and the biosynthesis of specialized metabolites, including compounds previously associated with improving fruit quality traits. The discovery of thousands of genes, not present in the previous reference genomes for blueberry and cranberry, will serve as the basis of future research and as potential targets for future breeding efforts. The pangenome, as a multiple-sequence alignment, as well as individual annotated genomes, are publicly available for analysis on the Genome Database for Vaccinium-a curated and integrated web-based relational database. Lastly, the core-gene predictions from the pangenomes will serve useful to develop a community genotyping platform to guide future molecular breeding efforts across the family.

Deciphering the genetic architecture of fruit color in strawberry

Fruits of Fragaria species usually have an appealing bright red color due to the accumulation of anthocyanins, water-soluble flavonoid pigments. Octoploid cultivated strawberry (Fragaria × ananassa) is a major horticultural crop for which fruit color and associated nutritional value are main breeding targets. Great diversity in fruit color intensity and pattern is observed not only in cultivated strawberry but also in wild relatives such as its octoploid progenitor F. chiloensis or the diploid woodland strawberry F. vesca, a model for fruit species in the Rosaceae. This review examines our understanding of fruit color formation in strawberry and how ongoing developments will advance it. Natural variations of fruit color as well as color changes during fruit development or in response to several cues have been used to explore the anthocyanin biosynthetic pathway and its regulation. So far, the successful identification of causal genetic variants has been largely driven by the availability of high-throughput genotyping tools and high-quality reference genomes of F. vesca and F. × ananassa. The current completion of haplotype-resolved genomes of F. × ananassa combined with QTL mapping will accelerate the exploitation of the untapped genetic diversity of fruit color and help translate the findings into strawberry improvement.

An arginine-to-histidine mutation in flavanone-3-hydroxylase results in pink strawberry fruits

Fruit color is a very important external commodity factor for consumers. Compared to the most typical red octoploid strawberry (Fragaria × ananassa), the pink strawberry often sells for a more expensive price and has a higher economic benefit due to its outstanding color. However, few studies have examined the molecular basis of pink-colored strawberry fruit. Through an EMS mutagenesis of woodland strawberry (Fragaria vesca), we identified a mutant with pink fruits and green petioles. Bulked-segregant analysis sequencing analysis and gene function verification confirmed that the responsible mutation resides in a gene encoding flavanone-3-hydroxylase (F3H) in the anthocyanin synthesis pathway. This nonsynonymous mutation results in an arginine-to-histidine change at position 130 of F3H. Molecular docking experiments showed that the arginine-to-histidine mutation results in a reduction of intermolecular force-hydrogen bonding between the F3H protein and its substrates. Enzymatic experiments showed a greatly reduced ability of the mutated F3H protein to catalyze the conversion of the substrates and hence a blockage of the anthocyanin synthesis pathway. The discovery of a key residue in the F3H gene controlling anthocyanin synthesis provides a clear target of modification for the molecular breeding of strawberry varieties with pink-colored fruits, which may be of great commercial value.

Molecular bases of strawberry fruit quality traits: Advances, challenges, and opportunities

The strawberry is one of the world's most popular fruits, providing humans with vitamins, fibers, and antioxidants. Cultivated strawberry (Fragaria × ananassa) is an allo-octoploid and highly heterozygous, making it a challenge for breeding, quantitative trait locus (QTL) mapping, and gene discovery. Some wild strawberry relatives, such as Fragaria vesca, have diploid genomes and are becoming laboratory models for the cultivated strawberry. Recent advances in genome sequencing and CRISPR-mediated genome editing have greatly improved the understanding of various aspects of strawberry growth and development in both cultivated and wild strawberries. This review focuses on fruit quality traits that are most relevant to the consumers, including fruit aroma, sweetness, color, firmness, and shape. Recently available phased-haplotype genomes, single nucleotide polymorphism (SNP) arrays, extensive fruit transcriptomes, and other big data have made it possible to locate key genomic regions or pinpoint specific genes that underlie volatile synthesis, anthocyanin accumulation for fruit color, and sweetness intensity or perception. These new advances will greatly facilitate marker-assisted breeding, the introgression of missing genes into modern varieties, and precise genome editing of selected genes and pathways. Strawberries are poised to benefit from these recent advances, providing consumers with fruit that is tastier, longer-lasting, healthier, and more beautiful.

Exploring Pan-Genomes: An Overview of Resources and Tools for Unraveling Structure, Function, and Evolution of Crop Genes and Genomes

The availability of multiple sequenced genomes from a single species made it possible to explore intra- and inter-specific genomic comparisons at higher resolution and build clade-specific pan-genomes of several crops. The pan-genomes of crops constructed from various cultivars, accessions, landraces, and wild ancestral species represent a compendium of genes and structural variations and allow researchers to search for the novel genes and alleles that were inadvertently lost in domesticated crops during the historical process of crop domestication or in the process of extensive plant breeding. Fortunately, many valuable genes and alleles associated with desirable traits like disease resistance, abiotic stress tolerance, plant architecture, and nutrition qualities exist in landraces, ancestral species, and crop wild relatives. The novel genes from the wild ancestors and landraces can be introduced back to high-yielding varieties of modern crops by implementing classical plant breeding, genomic selection, and transgenic/gene editing approaches. Thus, pan-genomic represents a great leap in plant research and offers new avenues for targeted breeding to mitigate the impact of global climate change. Here, we summarize the tools used for pan-genome assembly and annotations, web-portals hosting plant pan-genomes, etc. Furthermore, we highlight a few discoveries made in crops using the pan-genomic approach and future potential of this emerging field of study.

Application of Multi-Perspectives in Tea Breeding and the Main Directions

Tea plants are an economically important crop and conducting research on tea breeding contributes to enhancing the yield and quality of tea leaves as well as breeding traits that satisfy the requirements of the public. This study reviews the current status of tea plants germplasm resources and their utilization, which has provided genetic material for the application of multi-omics, including genomics and transcriptomics in breeding. Various molecular markers for breeding were designed based on multi-omics, and available approaches in the direction of high yield, quality and resistance in tea plants breeding are proposed. Additionally, future breeding of tea plants based on single-cellomics, pangenomics, plant-microbe interactions and epigenetics are proposed and provided as references. This study aims to provide inspiration and guidance for advancing the development of genetic breeding in tea plants, as well as providing implications for breeding research in other crops.

Blueberry and cranberry pangenomes as a resource for future genetic studies and breeding efforts

Domestication of cranberry and blueberry began in the United States in the early 1800s and 1900s, respectively, and in part owing to their flavors and health-promoting benefits are now cultivated and consumed worldwide. The industry continues to face a wide variety of production challenges (e.g. disease pressures) as well as a demand for higher-yielding cultivars with improved fruit quality characteristics. Unfortunately, molecular tools to help guide breeding efforts for these species have been relatively limited compared with those for other high-value crops. Here, we describe the construction and analysis of the first pangenome for both blueberry and cranberry. Our analysis of these pangenomes revealed both crops exhibit great genetic diversity, including the presence-absence variation of 48.4% genes in highbush blueberry and 47.0% genes in cranberry. Auxiliary genes, those not shared by all cultivars, are significantly enriched with molecular functions associated with disease resistance and the biosynthesis of specialized metabolites, including compounds previously associated with improving fruit quality traits. The discovery of thousands of genes, not present in the previous reference genomes for blueberry and cranberry, will serve as the basis of future research and as potential targets for future breeding efforts. The pangenome, as a multiple-sequence alignment, as well as individual annotated genomes, are publicly available for analysis on the Genome Database for Vaccinium - a curated and integrated web-based relational database. Lastly, the core-gene predictions from the pangenomes will serve useful to develop a community genotyping platform to guide future molecular breeding efforts across the family.

Haplotype-resolved genomes of wild octoploid progenitors illuminate genomic diversifications from wild relatives to cultivated strawberry

Strawberry is an emerging model for studying polyploid genome evolution and rapid domestication of fruit crops. Here we report haplotype-resolved genomes of two wild octoploids (Fragaria chiloensis and Fragaria virginiana), the progenitor species of cultivated strawberry. Substantial variation is identified between species and between haplotypes. We redefine the four subgenomes and track the genetic contributions of diploid species by additional sequencing of the diploid F. nipponica genome. We provide multiple lines of evidence that F. vesca and F. iinumae, rather than other described extant species, are the closest living relatives of these wild and cultivated octoploids. In response to coexistence with quadruplicate gene copies, the octoploid strawberries have experienced subgenome dominance, homoeologous exchanges and coordinated expression of homoeologous genes. However, some homoeologues have substantially altered expression bias after speciation and during domestication. These findings enhance our understanding of the origin, genome evolution and domestication of strawberries.

Roles and evolution of four LEAFY homologs in floral patterning and leaf development in woodland strawberry

The plant-specific transcription factor LEAFY (LFY), generally maintained as a single-copy gene in most angiosperm species, plays critical roles in flower development. The woodland strawberry (Fragaria vesca) possesses four LFY homologs in the genome; however, their respective functions and evolution remain unknown. Here, we identified and validated that mutations in one of the four LFY homologs, FveLFYa, cause homeotic conversion of floral organs and reiterative outgrowth of ectopic flowers. In contrast to FveLFYa, FveLFYb/c/d appear dispensable under normal growth conditions, as fvelfyc mutants are indistinguishable from wild type and FveLFYb and FveLFYd are barely expressed. Transgenic analysis and yeast one-hybrid assay showed that FveLFYa and FveLFYb, but not FveLFYc and FveLFYd, are functionally conserved with AtLFY in Arabidopsis (Arabidopsis thaliana). Unexpectedly, LFY-binding site prediction and yeast one-hybrid assay revealed that the transcriptional links between LFY and the APETALA1 (AP1) promoter/the large AGAMOUS (AG) intron are missing in F. vesca, which is due to the loss of LFY-binding sites. The data indicate that mutations in cis-regulatory elements could contribute to LFY evolution. Moreover, we showed that FveLFYa is involved in leaf development, as approximately 30% of mature leaves have smaller or fewer leaflets in fvelfya. Phylogenetic analysis indicated that LFY homologs in Fragaria species may arise from recent duplication events in their common ancestor and are undergoing convergent gene loss. Together, these results provide insight into the role of LFY in flower and leaf development in strawberry and have important implications for the evolution of LFY.

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.

Rapid evolution of T2/S-RNase genes in Fragaria linked to multiple transitions from self-incompatibility to self-compatibility

The T2/RNase gene family is widespread in eukaryotes, and particular members of this family play critical roles in the gametophytic self-incompatibility (GSI) system in plants. Wild diploid strawberry (Fragaria) species have diversified their sexual systems via self-incompatible and self-compatible traits, yet how these traits evolved in Fragaria remains elusive. By integrating the published and de novo assembled genomes and the newly generated RNA-seq data, members of the RNase T2 gene family were systematically identified in six Fragaria species, including three self-incompatible species (Fragaria nipponica, Fragaria nubicola, and Fragaria viridis) and three self-compatible species (Fragaria nilgerrensis, Fragaria vesca, and Fragaria iinumae). In total, 115 RNase T2 genes were identified in the six Fragaria genomes and can be classified into three classes (I-III) according to phylogenetic analysis. The identified RNase T2 genes could be divided into 22 homologous gene sets according to amino acid sequence similarity and phylogenetic and syntenic relationships. We found that extensive gene loss and pseudogenization coupled with small-scale duplications mainly accounted for variations in the RNase T2 gene numbers in Fragaria. Multiple copies of homologous genes were mainly generated from tandem and segmental duplication events. Furthermore, we newly identified five S-RNase genes in three self-incompatible Fragaria genomes, including two in F. nipponica, two in F. viridis, and one in F. nubicola, which fit for typical features of a pistil determinant, including highly pistil-specific expression, highly polymorphic proteins and alkaline isoelectric point (pI), while no S-RNase genes were found in all three self-compatible Fragaria species. Surprisingly, these T2/S-RNase genes contain at least one large intron (>10 kb). This study revealed that the rapid evolution of T2/S-RNase genes within the Fragaria genus could be associated with its sexual mode, and repeated evolution of the self-compatible traits in Fragaria was convergent via losses of S-RNase.

The population genomic analyses of chloroplast genomes shed new insights on the complicated ploidy and evolutionary history in Fragaria

The genus Fragaria consists of a rich diversity of ploidy levels with diploid (2x), tetraploid (4x), pentaploid (5x), hexaploidy (6x), octoploid (8x) and decaploid (10x) species. Only a few studies have explored the origin of diploid and octoploid strawberry, and little is known about the roles of tetraploidy and hexaploidy during the evolution of octoploid strawberry. The chloroplast genome is usually a stable circular genome and is widely used in investigating the evolution and matrilineal identification. Here, we assembled the chloroplast genomes of F. x ananassa cv. 'Benihoppe' (8x) using Illumina and HiFi data seperately. The genome alignment results showed that more InDels were located in the chloroplast genomes based on the PacBio HiFi data than Illumina data. We obtain highly accurate chloroplast genomes assembled through GetOrganelle using Illumina reads. We assembled 200 chloroplast genomes including 198 Fragaria (21 species) and 2 Potentilla samples. Analyses of sequence variation, phylogenetic and PCA analyses showed that Fragaria was divided into five groups. F. iinumae, F. nilgerrensis and all octoploid accessions formed Group A, C and E separately. Species native to western China were clustered into Group B. Group D consisted of F. virdis, F. orientalis, F. moschata, and F. vesca. STRUCTURE and haplotype network confirmed that the diploid F. vesca subsp. bracteata was the last maternal donator of octoploid strawberry. The dN/dS ratio estimated for the protein-coding genes revealed that genes involved in ATP synthase and photosystem function were under positive selection. These findings demonstrate the phylogeny of totally 21 Fragaria species and the origin of octoploid species. F. vesca was the last female donator of octoploid, which confirms the hypothesis that the hexaploid species F. moschata may be an evolutionary intermediate between the diploids and wild octoploid species.

Going broad and deep: sequencing-driven insights into plant physiology, evolution, and crop domestication

Deep sequencing is a term that has become embedded in the plant genomic literature in recent years and with good reason. A torrent of (largely) high-quality genomic and transcriptomic data has been collected and most of this has been publicly released. Indeed, almost 1000 plant genomes have been reported (www.plabipd.de) and the 2000 Plant Transcriptomes Project has long been completed. The EarthBioGenome project will dwarf even these milestones. That said, massive progress in understanding plant physiology, evolution, and crop domestication has been made by sequencing broadly (across a species) as well as deeply (within a single individual). We will outline the current state of the art in genome and transcriptome sequencing before we briefly review the most visible of these broad approaches, namely genome-wide association and transcriptome-wide association studies, as well as the compilation of pangenomes. This will include both (i) the most commonly used methods reliant on single nucleotide polymorphisms and short InDels and (ii) more recent examples which consider structural variants. We will subsequently present case studies exemplifying how their application has brought insight into either plant physiology or evolution and crop domestication. Finally, we will provide conclusions and an outlook as to the perspective for the extension of such approaches to different species, tissues, and biological processes.

Hexokinase1: A glucose sensor involved in drought stress response and sugar metabolism depending on its kinase activity in strawberry

Hexokinase1 (HXK1) is a bifunctional enzyme that plays indispensable roles in plant growth, nitrogen utilization, and stress resistance. However, information on the HXK family members of strawberries and their functions in glucose sensing and metabolic regulation is scarce. In the present study, four HXKs were firstly identified in the genome of Fragaria vesca and F. pentaphylla. The conserved domains of the HXK1s were confirmed, and a site-directed mutation (S177A) was introduced into the FpHXK1. FpHXK1, which shares the highest identity with the AtHXK1 was able to restore the glucose sensitivity and developmental defects of the Arabidopsis gin2-1 mutant, but not its kinase-activity-impaired mutant (FpHXK1^S177A ). The transcription of FpHXK1 was dramatically up-regulated under PEG-simulated drought stress conditions. The inhibition of the HXK kinase activity delayed the strawberry plant's responses to drought stress. Transient overexpression of the FpHXK1 and its kinase-impaired mutant differentially affected the level of glucose, sucrose, anthocyanins, and total phenols in strawberry fruits. All these results indicated that the FpHXK1, acting as a glucose sensor, was involved in drought stress response and sugar metabolism depending on its kinase activity.

Graph-based pan-genomes: increased opportunities in plant genomics

Due to the development of sequencing technology and the great reduction in sequencing costs, an increasing number of plant genomes have been assembled, and numerous genomes have revealed large amounts of variations. However, a single reference genome does not allow the exploration of species diversity, and therefore the concept of pan-genome was developed. A pan-genome is a collection of all sequences available for a species, including a large number of consensus sequences, large structural variations, and small variations including single nucleotide polymorphisms and insertions/deletions. A simple linear pan-genome does not allow these structural variations to be intuitively characterized, so graph-based pan-genomes have been developed. These pan-genomes store sequence and structural variation information in the form of nodes and paths to store and display species variation information in a more intuitive manner. The key role of graph-based pan-genomes is to expand the coordinate system of the linear reference genome to accommodate more regions of genetic diversity. Here, we review the origin and development of graph-based pan-genomes, explore their application in plant research, and further highlight the application of graph-based pan-genomes for future plant breeding.

Integrated transcriptome and methylome analyses reveal the molecular regulation of drought stress in wild strawberry (Fragaria nilgerrensis)

BACKGROUND

Fragaria nilgerrensis, which is a diploid wild strawberry with excellent drought-resistance, would provide useful candidate genes for improving drought resistance of cultivated strawberry. So far, its molecular regulatory networks involved in drought stress are unclear. We therefore investigated the drought response regulatory networks of F. nilgerrensis based on the integrated analysis of DNA methylation, transcriptome and physiological traits during four time points under drought stress.  RESULTS: The most differentially expressed genes and the physiological changes were found at 8 days (T8) compared with 0 day (T0, control). Methylome analysis revealed slight dynamic changes in genome-wide mC levels under drought conditions, while the most hypomethylated and hypermethylated regions were identified at T4 and T8. Association analysis of the methylome and transcriptome revealed that unexpressed genes exhibited expected hypermethylation levels in mCHG and mCHH contexts, and highly expressed genes exhibited corresponding hypomethylation levels in the gene body, but mCG contexts showed the opposite trend. Then, 835 differentially methylated and expressed genes were identified and grouped into four clustering patterns to characterize their functions. The genes with either negative or positive correlation between methylation and gene expression were mainly associated with kinases, Reactive Oxygen Species (ROS) synthesis, scavenging, and the abscisic acid (ABA) signal pathway. Consistently, weighted gene co-expression network analysis (WGCNA) revealed Hub genes including NCED, CYP707A2, PP2Cs and others that play important roles in the ABA signaling pathway.

CONCLUSION

F. nilgerrensis drought is dominated by ABA-dependent pathways, possibly accompanied by ABA-independent crosstalk. DNA methylation may affect gene expression, but their correlation was more subtle and multiple types of association exist. Maintaining the balance between ROS regeneration and scavenging is an important factor in drought resistance in F. nilgerrensis. These results deepen our understanding of drought resistance and its application in breeding in strawberry plants.

Genome and pan-genome assembly of asparagus bean (Vigna unguiculata ssp. sesquipedialis) reveal the genetic basis of cold adaptation

Asparagus bean (Vigna unguiculata ssp. sesquipedialis) is an important cowpea subspecies. We assembled the genomes of Ningjiang 3 (NJ, 550.31 Mb) and Dubai bean (DB, 564.12 Mb) for comparative genomics analysis. The whole-genome duplication events of DB and NJ occurred at 64.55 and 64.81 Mya, respectively, while the divergence between soybean and Vigna occurred in the Paleogene period. NJ genes underwent positive selection and amplification in response to temperature and abiotic stress. In species-specific gene families, NJ is mainly enriched in response to abiotic stress, while DB is primarily enriched in respiration and photosynthesis. We established the pan-genomes of four accessions (NJ, DB, IT97K-499-35 and Xiabao II) and identified 20,336 (70.5%) core genes present in all the accessions, 6,507 (55.56%) variable genes in two individuals, and 2,004 (6.95%) unique genes. The final pan genome is 616.35 Mb, and the core genome is 399.78 Mb. The variable genes are manifested mainly in stress response functions, ABC transporters, seed storage, and dormancy control. In the pan-genome sequence variation analysis, genes affected by presence/absence variants were enriched in biological processes associated with defense responses, immune system processes, signal transduction, and agronomic traits. The results of the present study provide genetic data that could facilitate efficient asparagus bean genetic improvement, especially in producing cold-adapted asparagus bean.

Fragaria mitogenomes evolve rapidly in structure but slowly in sequence and incur frequent multinucleotide mutations mediated by microinversions

Plant mitochondrial DNA has been described as evolving rapidly in structure but slowly in sequence. However, many of the noncoding portions of plant mitogenomes are not homologous among species, raising questions about the rate and spectrum of mutations in noncoding regions. Recent studies have suggested that the lack of homology in noncoding regions could be due to increased sequence divergence. We compared 30 kb of coding and 200 kb of noncoding DNA from 13 sequenced Fragaria mitogenomes, followed by analysis of the rate of sequence divergence, microinversion events and structural variations. Substitution rates in synonymous sites and nongenic sites are nearly identical, suggesting that the genome-wide point mutation rate is generally consistent. A surprisingly high number of large multinucleotide substitutions were detected in Fragaria mitogenomes, which may have resulted from microinversion events and could affect phylogenetic signal and local rate estimates. Fragaria mitogenomes preferentially accumulate deletions relative to insertions and substantial genomic arrangements, whereas mutation rates could positively associate with these sequence and structural changes among species. Together, these observations suggest that plant mitogenomes exhibit low point mutations genome-wide but exceptionally high structural variations, and our results favour a gain-and-loss model for the rapid loss of homology among plant mitogenomes.

Chromosome-level genome assembly of Fragaria pentaphylla using PacBio and Hi-C technologies

Fragaria pentaphylla, a wild diploid quinquefoliolate species of Fragaria, is native to Southwest China. It has two morphs of red and white fruit color in nature and has characteristics of unique fragrance and resistance, which made it not only a valuable breeding material but also a potential model plant for molecular function researches. Here, we generate a high-quality chromosome-level genome assembly of a F. pentaphylla accession, BAAFS-FP039 employing a combination of PacBio Long-Read Sequencing, Illumina Short-Read Sequencing, and Hi-C Sequencing. The assembled genome contained 256.74 Mb and a contig N50 length of 32.38 Mb, accounting for 99.9% of the estimated genome (256.77 Mb). Based on Hi-C data, seven pseudo-chromosomes of F. pentaphylla-FP039 genome were assembled, covering 99.39% of the genome assembly. The genome was composed of 44.61% repetitive sequences and 29,623 protein-coding genes, 97.62% of protein-coding genes could be functionally annotated. Phylogenetic and chromosome syntenic analysis revealed that F. pentaphylla-FP039 was closely related to F. nubicola. This high-quality genome could provides fundamental molecular resources for evolutionary studies, breeding efforts, and exploring the unique biological characteristics of F. pentaphylla.

Towards smart and sustainable development of modern berry cultivars in Europe

Fresh berries are a popular and important component of the human diet. The demand for high-quality berries and sustainable production methods is increasing globally, challenging breeders to develop modern berry cultivars that fulfill all desired characteristics. Since 1994, research projects have characterized genetic resources, developed modern tools for high-throughput screening, and published data in publicly available repositories. However, the key findings of different disciplines are rarely linked together, and only a limited range of traits and genotypes has been investigated. The Horizon2020 project BreedingValue will address these challenges by studying a broader panel of strawberry, raspberry and blueberry genotypes in detail, in order to recover the lost genetic diversity that has limited the aroma and flavor intensity of recent cultivars. We will combine metabolic analysis with sensory panel tests and surveys to identify the key components of taste, flavor and aroma in berries across Europe, leading to a high-resolution map of quality requirements for future berry cultivars. Traits linked to berry yields and the effect of environmental stress will be investigated using modern image analysis methods and modeling. We will also use genetic analysis to determine the genetic basis of complex traits for the development and optimization of modern breeding technologies, such as molecular marker arrays, genomic selection and genome-wide association studies. Finally, the results, raw data and metadata will be made publicly available on the open platform Germinate in order to meet FAIR data principles and provide the basis for sustainable research in the future.

Rearrangement and domestication as drivers of Rosaceae mitogenome plasticity

Background

The mitochondrion is an important cellular component in plants and that functions in producing vital energy for the cell. However, the evolution and structure of mitochondrial genomes (mitogenomes) remain unclear in the Rosaceae family. In this study, we assembled 34 Rosaceae mitogenomes and characterized genome variation, rearrangement rate, and selection signal variation within these mitogenomes.

Results

Comparative analysis of six genera from the Amygdaloideae and five from the Rosoideae subfamilies of Rosaceae revealed that three protein-coding genes were absent from the mitogenomes of five Rosoideae genera. Positive correlations between genome size and repeat content were identified in 38 Rosaceae mitogenomes. Twenty repeats with high recombination frequency (> 50%) provided evidence for predominant substoichiometric conformation of the mitogenomes. Variations in rearrangement rates were identified between eleven genera, and within the Pyrus, Malus, Prunus, and Fragaria genera. Based on population data, phylogenetic inferences from Pyrus mitogenomes supported two distinct maternal lineages of Asian cultivated pears. A Pyrus-specific deletion (DEL-D) in selective sweeps was identified based on the assembled genomes and population data. After the DEL-D sequence fragments originally arose, they may have experienced a subsequent doubling event via homologous recombination and sequence transfer in the Amygdaloideae; afterwards, this variant sequence may have significantly expanded to cultivated groups, thereby improving adaptation during the domestication process.

Conclusions

This study characterizes the variations in gene content, genome size, rearrangement rate, and the impact of domestication in Rosaceae mitogenomes and provides insights into their structural variation patterns and phylogenetic relationships.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-022-01383-3.

Plant pan-genomics: recent advances, new challenges, and roads ahead

Pan-genomics can encompass most of the genetic diversity of a species or population and has proved to be a powerful tool for studying genomic evolution and the origin and domestication of species, and for providing information for plant improvement. Plant genomics has greatly progressed because of improvements in sequencing technologies and the rapid reduction of sequencing costs. Nevertheless, pan-genomics still presents many challenges, including computationally intensive assembly methods, high costs with large numbers of samples, ineffective integration of big data, and difficulty in applying it to downstream multi-omics analysis and breeding research. In this review, we summarize the definition and recent achievements of plant pan-genomics, computational technologies used for pan-genome construction, and the applications of pan-genomes in plant genomics and molecular breeding. We also discuss challenges and perspectives for future pan-genomics studies and provide a detailed pipeline for sample selection, genome assembly and annotation, structural variation identification, and construction and application of graph-based pan-genomes. The aim is to provide important guidance for plant pan-genome research and a better understanding of the genetic basis of genome evolution, crop domestication, and phenotypic diversity for future studies.

Characterization of Expression and Epigenetic Features of Core Genes in Common Wheat

The availability of multiple wheat genome sequences enables us to identify core genes and characterize their genetic and epigenetic features, thereby advancing our understanding of their biological implications within individual plant species. It is, however, largely understudied in wheat. To this end, we reanalyzed genome sequences from 16 different wheat varieties and identified 62,299 core genes. We found that core and non-core genes have different roles in subgenome differentiation. Meanwhile, according to their expression profiles, these core genes can be classified into genes related to tissue development and stress responses, including 3376 genes highly expressed in both spikelets and at high temperatures. After associating with six histone marks and open chromatin, we found that these core genes can be divided into eight sub-clusters with distinct epigenomic features. Furthermore, we found that ca. 51% of the expressed transcription factors (TFs) were marked with both H3K27me3 and H3K4me3, indicative of the bivalency feature, which can be involved in tissue development through the TF-centered regulatory network. Thus, our study provides a valuable resource for the functional characterization of core genes in stress responses and tissue development in wheat.

Unraveling NPR-like Family Genes in Fragaria spp. Facilitated to Identify Putative NPR1 and NPR3/4 Orthologues Participating in Strawberry-Colletotrichum fructicola Interaction

The salicylic acid receptor NPR1 (nonexpressor of pathogenesis-related genes) and its paralogues NPR3 and NPR4 are master regulators of plant immunity. Commercial strawberry (Fragaria × ananassa) is a highly valued crop vulnerable to various pathogens. Historic confusions regarding the identity of NPR-like genes have hindered research in strawberry resistance. In this study, the comprehensive identification and phylogenic analysis unraveled this family, harboring 6, 6, 5, and 23 members in F. vesca, F. viridis, F. iinumae, and F. × ananassa, respectively. These genes were clustered into three clades, with each diploid member matching three to five homoalleles in F. × ananassa. Despite the high conservation in terms of gene structure, protein module, and functional residues/motifs/domains, substantial divergence was observed, hinting strawberry NPR proteins probably function in ways somewhat different from Arabidopsis. RT-PCR and RNAseq analysis evidenced the transcriptional responses of FveNPR1 and FxaNPR1a to Colletotrichum fructicola. Extended expression analysis for strawberry NPR-likes helped to us understand how strawberry orchestrate the NPRs-centered defense system against C. fructicola. The cThe current work supports that FveNPR1 and FxaNPR1a, as well as FveNPR31 and FxaNPR31a-c, were putative functional orthologues of AtNPR1 and AtNPR3/4, respectively. These findings set a solid basis for the molecular dissection of biological functions of strawberry NPR-like genes for improving disease resistance.

Dynamic Changes of DNA Methylation During Wild Strawberry (Fragaria nilgerrensis) Tissue Culture

Tissue culture is an important tool for asexual propagation and genetic transformation of strawberry plants. In plant tissue culture, variation of DNA methylation is a potential source of phenotypic variation in regenerated plants. However, the genome wide dynamic methylation patterns of strawberry tissue culture remain unclear. In this study, we used whole-genome bisulfite sequencing (WGBS) to study genomic DNA methylation changes of a wild strawberry Fragaria nilgerrensis at six stages: from explants of shoot tips to outplanting and acclimation. Global methylation levels showed that CG sites exhibited the highest methylation level in all stages with an average of 49.5%, followed by CHG (33.2%) and CHH (12.4%). Although CHH accounted for the lowest proportion of total cytosine methylation, it showed the most obvious methylation change and the most of these changes occurred in the transposable element regions. The overall methylation levels alternately decreased and increased during the entire tissue culture process and the distribution of DNA methylation was non-uniform among different genetic regions. Furthermore, much more differentially methylated regions (DMRs) were detected in dedifferentiation and redifferentiation stages and most of them were transposable elements, suggesting these processes involved activating or silencing of amounts of transposons. The functional enrichment of the DMR-related genes indicated that genes involved in hormone metabolic processes, plant development and the stress response changed methylation throughout the tissue culture process. Finally, the quantitative real-time PCR (qRT-PCR) was conducted to examine the association of methylation and gene expression of a set of different methylated genes. Our findings give deeper insight into the epigenetic regulation of gene expression during the plant tissue cultures process, which will be useful in the efficient control of somaclonal variations and in crop improvement.