The Tetracentron genome provides insight into the early evolution of eudicots and the formation of vessel elements

Chromosome-level genome assembly of the threatened resource plant Cinnamomum chago

Cinnamomum chago is a tree species endemic to Yunnan province, China, with potential economic value, phylogenetic importance, and conservation priority. We assembled the genome of C. chago using multiple sequencing technologies, resulting in a high-quality, chromosomal-level genome with annotation information. The assembled genome size is approximately 1.06 Gb, with a contig N50 length of 92.10 Mb. About 99.92% of the assembled sequences could be anchored to 12 pseudo-chromosomes, with only one gap, and 63.73% of the assembled genome consists of repeat sequences. In total, 30,497 genes were recognized according to annotation, including 28,681 protein-coding genes. This high-quality chromosome-level assembly and annotation of C. chago will assist us in the conservation and utilization of this valuable resource, while also providing crucial data for studying the evolutionary relationships within the Cinnamomum genus, offering opportunities for further research and exploration of its diverse applications.

Genomic data and ecological niche modeling reveal an unusually slow rate of molecular evolution in the Cretaceous Eupteleaceae

Living fossils are evidence of long-term sustained ecological success. However, whether living fossils have little molecular changes remains poorly known, particularly in plants. Here, we have introduced a novel method that integrates phylogenomic, comparative genomic, and ecological niche modeling analyses to investigate the rate of molecular evolution of Eupteleaceae, a Cretaceous relict angiosperm family endemic to East Asia. We assembled a high-quality chromosome-level nuclear genome, and the chloroplast and mitochondrial genomes of a member of Eupteleaceae (Euptelea pleiosperma). Our results show that Eupteleaceae is most basal in Ranunculales, the earliest-diverging order in eudicots, and shares an ancient whole-genome duplication event with the other Ranunculales. We document that Eupteleaceae has the slowest rate of molecular changes in the observed angiosperms. The unusually low rate of molecular evolution of Eupteleaceae across all three independent inherited genomes and genes within each of the three genomes is in association with its conserved genome architecture, ancestral woody habit, and conserved niche requirements. Our findings reveal the evolution and adaptation of living fossil plants through large-scale environmental change and also provide new insights into early eudicot diversification.

Cryptic divergences and repeated hybridizations within the endangered "living fossil" dove tree (Davidia involucrata) revealed by whole genome resequencing

The identification and understanding of cryptic intraspecific evolutionary units (lineages) are crucial for planning effective conservation strategies aimed at preserving genetic diversity in endangered species. However, the factors driving the evolution and maintenance of these intraspecific lineages in most endangered species remain poorly understood. In this study, we conducted resequencing of 77 individuals from 22 natural populations of Davidia involucrata, a "living fossil" dove tree endemic to central and southwest China. Our analysis revealed the presence of three distinct local lineages within this endangered species, which emerged approximately 3.09 and 0.32 million years ago. These divergence events align well with the geographic and climatic oscillations that occurred across the distributional range. Additionally, we observed frequent hybridization events between the three lineages, resulting in the formation of hybrid populations in their adjacent as well as disjunct regions. These hybridizations likely arose from climate-driven population expansion and/or long-distance gene flow. Furthermore, we identified numerous environment-correlated gene variants across the total and many other genes that exhibited signals of positive evolution during the maintenance of two major local lineages. Our findings shed light on the highly dynamic evolution underlying the remarkably similar phenotype of this endangered species. Importantly, these results not only provide guidance for the development of conservation plans but also enhance our understanding of evolutionary past for this and other endangered species with similar histories.

Wood of trees: Cellular structure, molecular formation, and genetic engineering

Wood is an invaluable asset to human society due to its renewable nature, making it suitable for both sustainable energy production and material manufacturing. Additionally, wood derived from forest trees plays a crucial role in sequestering a significant portion of the carbon dioxide fixed during photosynthesis by terrestrial plants. Nevertheless, with the expansion of the global population and ongoing industrialization, forest coverage has been substantially decreased, resulting in significant challenges for wood production and supply. Wood production practices have changed away from natural forests toward plantation forests. Thus, understanding the underlying genetic mechanisms of wood formation is the foundation for developing high-quality, fast-growing plantation trees. Breeding ideal forest trees for wood production using genetic technologies has attracted the interest of many. Tremendous studies have been carried out in recent years on the molecular, genetic, and cell-biological mechanisms of wood formation, and considerable progress and findings have been achieved. These studies and findings indicate enormous possibilities and prospects for tree improvement. This review will outline and assess the cellular and molecular mechanisms of wood formation, as well as studies on genetically improving forest trees, and address future development prospects.

Subgenome phasing for complex allopolyploidy: case-based benchmarking and recommendations

Accurate subgenome phasing is crucial for understanding the origin, evolution and adaptive potential of polyploid genomes. SubPhaser and WGDI software are two common methodologies for subgenome phasing in allopolyploids, particularly in scenarios lacking known diploid progenitors. Triggered by a recent debate over the subgenomic origins of the cultivated octoploid strawberry, we examined four well-documented complex allopolyploidy cases as benchmarks, to evaluate and compare the accuracy of the two software. Our analysis demonstrates that the subgenomic structure phased by both software is in line with prior research, effectively tracing complex allopolyploid evolutionary trajectories despite the limitations of each software. Furthermore, using these validated methodologies, we revisited the controversial issue regarding the progenitors of the octoploid strawberry. The results of both methodologies reaffirm Fragaria vesca and Fragaria iinumae as progenitors of the octoploid strawberry. Finally, we propose recommendations for enhancing the accuracy of subgenome phasing in future studies, recognizing the potential of integrated tools for advanced complex allopolyploidy research and offering a new roadmap for robust subgenome-based phylogenetic analysis.

Telomere-to-telomere assembly of cassava genome reveals the evolution of cassava and divergence of allelic expression

Cassava is a crucial crop that makes a significant contribution to ensuring human food security. However, high-quality telomere-to-telomere cassava genomes have not been available up to now, which has restricted the progress of haploid molecular breeding for cassava. In this study, we constructed two nearly complete haploid resolved genomes and an integrated, telomere-to-telomere gap-free reference genome of an excellent cassava variety, 'Xinxuan 048', thereby providing a new high-quality genomic resource. Furthermore, the evolutionary history of several species within the Euphorbiaceae family was revealed. Through comparative analysis of haploid genomes, it was found that two haploid genomes had extensive differences in linear structure, transcriptome features, and epigenetic characteristics. Genes located within the highly divergent regions and differentially expressed alleles are enriched in the functions of auxin response and the starch synthesis pathway. The high heterozygosity of cassava 'Xinxuan 048' leads to rapid trait segregation in the first selfed generation. This study provides a theoretical basis and genomic resource for molecular breeding of cassava haploids.

Chromosome-level genome assembly and population genetic analysis of a near-threatened rosewood species (Dalbergia cultrata Pierre Graham ex Benth) provide insights into its evolutionary and cold stress responses

Dalbergia cultrata Pierre Graham ex Benth (D. cultrata) is a precious rosewood tree species that grows in the tropical and subtropical regions of Asia. In this study, we used PacBio long-reading sequencing technology and Hi-C assistance to sequence and assemble the reference genome of D. cultrata. We generated 171.47 Gb PacBio long reads and 72.43 Gb Hi-C data and yielded an assembly of 10 pseudochromosomes with a total size of 690.99 Mb and Scaffold N50 of 65.76 Mb. The analysis of specific genes revealed that the triterpenoids represented by lupeol may play an important role in D. cultrata's potential medicinal value. Using the new reference genome, we analyzed the resequencing of 19 Dalbergia accessions and found that D. cultrata and D. cochinchinensis have the latest genetic relationship. Transcriptome sequencing of D. cultrata leaves grown under cold stress revealed that MYB transcription factor and E3 ubiquitin ligase may be playing an important role in the cold response of D. cultrata. Genome resources and identified genetic variation, especially those genes related to the biosynthesis of phytochemicals and cold stress response, will be helpful for the introduction, domestication, utilization, and further breeding of Dalbergia species.

Chromosome-scale genomes of commercial timber trees (Ochroma pyramidale, Mesua ferrea, and Tectona grandis)

Wood is the most important natural and endlessly renewable source of energy. Despite the ecological and economic importance of wood, many aspects of its formation have not yet been investigated. We performed chromosome-scale genome assemblies of three timber trees (Ochroma pyramidale, Mesua ferrea, and Tectona grandis) which exhibit different wood properties such as wood density, hardness, growth rate, and fiber cell wall thickness. The combination of 10X, stLFR, Hi-Fi sequencing and HiC data led us to assemble high-quality genomes evident by scaffold N50 length of 55.97 Mb (O. pyramidale), 22.37 Mb (M. ferrea) and 14.55 Mb (T. grandis) with >97% BUSCO completeness of the assemblies. A total of 35774, 24027, and 44813 protein-coding genes were identified in M. ferrea, T. grandis and O. pyramidale, respectively. The data generated in this study is anticipated to serve as a valuable genetic resource and will promote comparative genomic analyses, and it is of practical importance in gaining a further understanding of the wood properties in non-model woody species.

Manipulating microRNA miR408 enhances both biomass yield and saccharification efficiency in poplar

The conversion of lignocellulosic feedstocks to fermentable sugar for biofuel production is inefficient, and most strategies to enhance efficiency directly target lignin biosynthesis, with associated negative growth impacts. Here we demonstrate, for both laboratory- and field-grown plants, that expression of Pag-miR408 in poplar (Populus alba × P. glandulosa) significantly enhances saccharification, with no requirement for acid-pretreatment, while promoting plant growth. The overexpression plants show increased accessibility of cell walls to cellulase and scaffoldin cellulose-binding modules. Conversely, Pag-miR408 loss-of-function poplar shows decreased cell wall accessibility. Overexpression of Pag-miR408 targets three Pag-LACCASES, delays lignification, and modestly reduces lignin content, S/G ratio and degree of lignin polymerization. Meanwhile, the LACCASE loss of function mutants exhibit significantly increased growth and cell wall accessibility in xylem. Our study shows how Pag-miR408 regulates lignification and secondary growth, and suggest an effective approach towards enhancing biomass yield and saccharification efficiency in a major bioenergy crop.

A method for analyzing programmed cell death in xylem development by flow cytometry

Programmed cell death (PCD) is a genetically regulated developmental process leading to the death of specific types of plant cells, which plays important roles in plant development and growth such as wood formation. However, an efficient method needs to be established to study PCD in woody plants. Flow cytometry is widely utilized to evaluate apoptosis in mammalian cells, while it is rarely used to detect PCD in plants, especially in woody plants. Here, we reported that the xylem cell protoplasts from poplar stem were stained with a combination of fluorescein annexin V-FITC and propidium iodide (PI) and then sorted by flow cytometry. As expected, living cells (annexin V-FITC negative/PI negative), early PCD cells (annexin V-FITC positive/PI negative), and late PCD cells (annexin V-FITC positive/PI positive) could be finely distinguished through this method and then subjected for quantitative analysis. The expression of cell-type- and developmental stages-specific marker genes was consistent with the cell morphological observation. Therefore, the newly developed fluorescence-activated cell sorting (FACS) method can be used to study PCD in woody plants, which will be beneficial for studying the molecular mechanisms of wood formation.

A high-quality Buxus austro-yunnanensis (Buxales) genome provides new insights into karyotype evolution in early eudicots

Background

Eudicots are the most diverse group of flowering plants that compromise five well-defined lineages: core eudicots, Ranunculales, Proteales, Trochodendrales, and Buxales. However, the phylogenetic relationships between these five lineages and their chromosomal evolutions remain unclear, and a lack of high-quality genome analyses for Buxales has hindered many efforts to address this knowledge gap.

Results

Here, we present a high-quality chromosome-level genome of Buxus austro-yunnanensis (Buxales). Our phylogenomic analyses revealed that Buxales and Trochodendrales are genetically similar and classified as sisters. Additionally, both are sisters to the core eudicots, while Ranunculales was found to be the first lineage to diverge from these groups. Incomplete lineage sorting and hybridization were identified as the main contributors to phylogenetic discordance (34.33%) between the lineages. In fact, B. austro-yunnanensis underwent only one whole-genome duplication event, and collinear gene phylogeny analyses suggested that separate independent polyploidizations occurred in the five eudicot lineages. Using representative genomes from these five lineages, we reconstructed the ancestral eudicot karyotype (AEK) and generated a nearly gapless karyotype projection for each eudicot species. Within core eudicots, we recovered one common chromosome fusion event in asterids and malvids, respectively. Further, we also found that the previously reported fused AEKs in Aquilegia (Ranunculales) and Vitis (core eudicots) have different fusion positions, which indicates that these two species have different karyotype evolution histories.

Conclusions

Based on our phylogenomic and karyotype evolution analyses, we revealed the likely relationships and evolutionary histories of early eudicots. Ultimately, our study expands genomic resources for early-diverging eudicots.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-022-01420-1.

Pervasive genome duplications across the plant tree of life and their links to major evolutionary innovations and transitions

Whole-genome duplication (WGD) has occurred repeatedly during plant evolution and diversification, providing genetic layers for evolving new functions and phenotypes. Advances in long-read sequencing technologies have enabled sequencing and assembly of over 1000 plant genomes spanning nearly 800 species, in which a large set of ancient WGDs has been uncovered. Here, we review the recently reported WGDs that occurred in major plant lineages and key evolutionary positions, and highlight their contributions to morphological innovation and adaptive evolution. Current gaps and challenges in integrating enormous volumes of sequenced plant genomes, accurately inferring WGDs, and developing web-based analysis tools are emphasized. Looking to the future, ambitious genome sequencing projects and global efforts may substantially recapitulate the plant tree of life based on broader sampling of phylogenetic diversity, reveal much of the timetable of ancient WGDs, and address the biological significance of WGDs in plant adaptation and radiation.

Paleopolyploidies and Genomic Fractionation in Major Eudicot Clades

Eudicots account for ~75% of living angiosperms, containing important food and energy crops. Recently, high-quality genome sequences of several eudicots including Aquilegia coerulea and Nelumbo nucifera have become available, providing an opportunity to investigate the early evolutionary characteristics of eudicots. We performed genomic hierarchical and event-related alignments to infer homology within and between representative species of eudicots. The results provide strong evidence for multiple independent polyploidization events during the early diversification of eudicots, three of which are likely to be allopolyploids: The core eudicot-common hexaploidy (ECH), Nelumbo-specific tetraploidy (NST), and Ranunculales-common tetraploidy (RCT). Using different genomes as references, we constructed genomic alignment to list the orthologous and paralogous genes produced by polyploidization and speciation. This could provide a fundamental framework for studying other eudicot genomes and gene(s) evolution. Further, we revealed significantly divergent evolutionary rates among these species. By performing evolutionary rate correction, we dated RCT to be ~118-134 million years ago (Mya), after Ranunculales diverged with core eudicots at ~123-139 Mya. Moreover, we characterized genomic fractionation resulting from gene loss and retention after polyploidizations. Notably, we revealed a high degree of divergence between subgenomes. In particular, synonymous nucleotide substitutions at synonymous sites (Ks) and phylogenomic analyses implied that A. coerulea might provide the subgenome(s) for the gamma-hexaploid hybridization.

Buxus and Tetracentron genomes help resolve eudicot genome history

Ancient whole-genome duplications (WGDs) characterize many large angiosperm lineages, including angiosperms themselves. Prominently, the core eudicot lineage accommodates 70% of all angiosperms and shares ancestral hexaploidy, termed gamma. Gamma arose via two WGDs that occurred early in eudicot history; however, the relative timing of these is unclear, largely due to the lack of high-quality genomes among early-diverging eudicots. Here, we provide complete genomes for Buxus sinica (Buxales) and Tetracentron sinense (Trochodendrales), representing the lineages most closely related to core eudicots. We show that Buxus and Tetracentron are both characterized by independent WGDs, resolve relationships among early-diverging eudicots and their respective genomes, and use the RACCROCHE pipeline to reconstruct ancestral genome structure at three key phylogenetic nodes of eudicot diversification. Our reconstructions indicate genome structure remained relatively stable during early eudicot diversification, and reject hypotheses of gamma arising via inter-lineage hybridization between ancestral eudicot lineages, involving, instead, only stem lineage core eudicot ancestors.

Gamma triplication arises via two whole-genome duplications early in eudicot history, but the relative timing of these is unclear. Here, the authors report the genomes of Buxales and Trochodendrales and reject the hypothesis of gamma arising via inter-lineage hybridization between ancestral eudicot lineages.

The genome of the endangered Macadamia jansenii displays little diversity but represents an important genetic resource for plant breeding

Macadamia, a recently domesticated expanding nut crop in the tropical and subtropical regions of the world, is one of the most economically important genera in the diverse and widely adapted Proteaceae family. All four species of Macadamia are rare in the wild with the most recently discovered, M. jansenii, being endangered. The M. jansenii genome has been used as a model for testing sequencing methods using a wide range of long read sequencing techniques. Here, we report a chromosome level genome assembly, generated using a combination of Pacific Biosciences sequencing and Hi-C, comprising 14 pseudo-molecules, with a N50 of 52 Mb and a total genome assembly size of 758 Mb of which 56% is repetitive. Completeness assessment revealed that the assembly covered -97.1% of the conserved single copy genes. Annotation predicted 31,591 protein coding genes and allowed the characterization of genes encoding biosynthesis of cyanogenic glycosides, fatty acid metabolism, and anti-microbial proteins. Re-sequencing of seven other genotypes confirmed low diversity and low heterozygosity within this endangered species. Important morphological characteristics of this species such as small tree size and high kernel recovery suggest that M. jansenii is an important source of these commercial traits for breeding. As a member of a small group of families that are sister to the core eudicots, this high-quality genome also provides a key resource for evolutionary and comparative genomics studies.

The genome of the endangered Macadamia jansenii displays little diversity but represents an important genetic resource for plant breeding

Macadamia, a recently domesticated expanding nut crop in the tropical and subtropical regions of the world, is one of the most economically important genera in the diverse and widely adapted Proteaceae family. All four species of Macadamia are rare in the wild with the most recently discovered, M. jansenii, being endangered. The M. jansenii genome has been used as a model for testing sequencing methods using a wide range of long read sequencing techniques. Here, we report a chromosome level genome assembly, generated using a combination of Pacific Biosciences sequencing and Hi-C, comprising 14 pseudo-molecules, with a N50 of 52 Mb and a total genome assembly size of 758 Mb of which 56% is repetitive. Completeness assessment revealed that the assembly covered -97.1% of the conserved single copy genes. Annotation predicted 31,591 protein coding genes and allowed the characterization of genes encoding biosynthesis of cyanogenic glycosides, fatty acid metabolism, and anti-microbial proteins. Re-sequencing of seven other genotypes confirmed low diversity and low heterozygosity within this endangered species. Important morphological characteristics of this species such as small tree size and high kernel recovery suggest that M. jansenii is an important source of these commercial traits for breeding. As a member of a small group of families that are sister to the core eudicots, this high-quality genome also provides a key resource for evolutionary and comparative genomics studies.

The Chloranthus sessilifolius genome provides insight into early diversification of angiosperms

Most extant angiosperms belong to Mesangiospermae, which comprises eudicots, monocots, magnoliids, Chloranthales and Ceratophyllales. However, phylogenetic relationships between these five lineages remain unclear. Here, we report the high-quality genome of a member of the Chloranthales lineage (Chloranthus sessilifolius). We detect only one whole genome duplication within this species and find that polyploidization events in different Mesangiospermae lineage are mutually independent. We also find that the members of all floral development-related gene lineages are present in C. sessilifolius despite its extremely simplified flower. The AP1 and PI genes, however, show a weak floral tissue-specialized expression. Our phylogenomic analyses suggest that Chloranthales and magnoliids are sister groups, and both are together sister to the clade comprising Ceratophyllales and eudicots, while the monocot lineage is sister to all other Mesangiospermae. Our findings suggest that in addition to hybridization, incomplete lineage sorting may largely account for phylogenetic inconsistencies between the observed gene trees.