Chloranthus genome provides insights into the early diversification of angiosperms

Evolution of Whirly1 in the angiosperms: sequence, splicing, and expression in a clade of early transitional mycoheterotrophic orchids

The plastid-targeted transcription factor Whirly1 (WHY1) has been implicated in chloroplast biogenesis, plastid genome stability, and fungal defense response, which together represent characteristics of interest for the study of autotrophic losses across the angiosperms. While gene loss in the plastid and nuclear genomes has been well studied in mycoheterotrophic plants, the evolution of the molecular mechanisms impacting genome stability is completely unknown. Here, we characterize the evolution of WHY1 in four early transitional mycoheterotrophic orchid species in the genus Corallorhiza by synthesizing the results of phylogenetic, transcriptomic, and comparative genomic analyses with WHY1 genomic sequences sampled from 21 orders of angiosperms. We found an increased number of non-canonical WHY1 isoforms assembled from all but the greenest Corallorhiza species, including intron retention in some isoforms. Within Corallorhiza, phylotranscriptomic analyses revealed the presence of tissue-specific differential expression of WHY1 in only the most photosynthetically capable species and a coincident increase in the number of non-canonical WHY1 isoforms assembled from fully mycoheterotrophic species. Gene- and codon-level tests of WHY1 selective regimes did not infer significant signal of either relaxed selection or episodic diversifying selection in Corallorhiza but did so for relaxed selection in the late-stage full mycoheterotrophic orchids Epipogium aphyllum and Gastrodia elata. Additionally, nucleotide substitutions that most likely impact the function of WHY1, such as nonsense mutations, were only observed in late-stage mycoheterotrophs. We propose that our findings suggest that splicing and expression changes may precede the selective shifts we inferred for late-stage mycoheterotrophic species, which therefore does not support a primary role for WHY1 in the transition to mycoheterotrophy in the Orchidaceae. Taken together, this study provides the most comprehensive view of WHY1 evolution across the angiosperms to date.

Evolutionary radiation strategy revealed in the Scarabaeidae with evidence of continuous spatiotemporal morphology and phylogenesis

Evolutionary biology faces the important challenge of determining how to interpret the relationship between selection pressures and evolutionary radiation. The lack of morphological evidence on cross-species research adds to difficulty of this challenge. We proposed a new paradigm for evaluating the evolution of branches through changes in characters on continuous spatiotemporal scales, for better interpreting the impact of biotic/abiotic drivers on the evolutionary radiation. It reveals a causal link between morphological changes and selective pressures: consistent deformation signals for all tested characters on timeline, which provided strong support for the evolutionary hypothesis of relationship between scarabs and biotic/abiotic drivers; the evolutionary strategies under niche differentiation, which were manifested in the responsiveness degree of functional morphological characters with different selection pressure. This morphological information-driven integrative approach sheds light on the mechanism of macroevolution under different selection pressures and is applicable to more biodiversity research.

De novo assembling a high-quality genome sequence of Amur grape (Vitis amurensis Rupr.) gives insight into Vitis divergence and sex determination

To date, there has been no high-quality sequence for genomes of the East Asian grape species, hindering biological and breeding efforts to improve grape cultivars. This study presents ~522 Mb of the Vitis amurensis (Va) genome sequence containing 27 635 coding genes. Phylogenetic analysis indicated that Vitis riparia (Vr) may have first split from the other two species, Va and Vitis vinifera (Vv). Divergent numbers of duplicated genes reserved among grapes suggests that the core eudicot-common hexaploidy (ECH) and the subsequent genome instability still play a non-negligible role in species divergence and biological innovation. Prominent accumulation of sequence variants might have improved cold resistance in Va, resulting in a more robust network of regulatory cold resistance genes, explaining why it is extremely cold-tolerant compared with Vv and Vr. In contrast, Va has preserved many fewer nucleotide binding site (NBS) disease resistance genes than the other grapes. Notably, multi-omics analysis identified one trans-cinnamate 4-monooxygenase gene positively correlated to the resveratrol accumulated during Va berry development. A selective sweep analysis revealed a hypothetical Va sex-determination region (SDR). Besides, a PPR-containing protein-coding gene in the hypothetical SDR may be related to sex determination in Va. The content and arrangement order of genes in the putative SDR of female Va were similar to those of female Vv. However, the putative SDR of female Va has lost one flavin-containing monooxygenase (FMO) gene and contains one extra protein-coding gene uncharacterized so far. These findings will improve the understanding of Vitis biology and contribute to the improvement of grape breeding.

The chromosome-level genomes of the herbal magnoliids Warburgia ugandensis and Saururus chinensis

Warburgia ugandensis and Saururus chinensis are two of the most important medicinal plants in magnoliids and are widely utilized in traditional Kenya and Chinese medicine, respectively. The absence of higher-quality reference genomes has hindered research on the medicinal compound biosynthesis mechanisms of these plants. We report the chromosome-level genome assemblies of W. ugandensis and S. chinensis, and generated 1.13 Gb and 0.53 Gb genomes from 74 and 27 scaffolds, respectively, using BGI-DIPSEQ, Nanopore, and Hi-C sequencing. The scaffold N50 lengths were 82.97 Mb and 48.53 Mb, and the assemblies were anchored to 14 and 11 chromosomes of W. ugandensis and S. chinensis, respectively. In total, 24,739 and 20,561 genes were annotated, and 98.5% and 98% of the BUSCO genes were fully represented, respectively. The chromosome-level genomes of W. ugandensis and S. chinensis will be valuable resources for understanding the genetics of these medicinal plants, studying the evolution of magnoliids and angiosperms and conserving plant genetic resources.

A high-quality haplotype genome of Michelia alba DC reveals differences in methylation patterns and flower characteristics

Michelia alba DC is a highly valuable ornamental plant of the Magnoliaceae family. This evergreen tropical tree commonly grows in Southeast Asia and is adored for its delightful fragrance. Our study assembled the M. alba haplotype genome MC and MM by utilizing Nanopore ultralong reads, Pacbio Hifi long reads and parental second-generation data. Moreover, the first methylation map of Magnoliaceae was constructed based on the methylation site data obtained using Nanopore data. Metabolomic datasets were generated from the flowers of three different species to assess variations in pigment and volatile compound accumulation. Finally, transcriptome data were generated to link genomic, methylation, and morphological patterns to reveal the reasons underlying the differences between M. alba and its parental lines in petal color, flower shape, and fragrance. We found that the AP1 and AP2 genes are crucial in M. alba petal formation, while the 4CL, PAL, and C4H genes control petal color. The data generated in this study serve as a foundation for future physiological and biochemical research on M. alba, facilitate the targeted improvement of M. alba varieties, and offer a theoretical basis for molecular research on Michelia L.

Chromosome-level Alstonia scholaris genome unveils evolutionary insights into biosynthesis of monoterpenoid indole alkaloids

Alstonia scholaris of the Apocynaceae family is a medicinal plant with a rich source of bioactive monoterpenoid indole alkaloids (MIAs), which possess anti-cancer activity like vinca alkaloids. To gain genomic insights into MIA biosynthesis, we assembled a high-quality chromosome-level genome for A. scholaris using nanopore and Hi-C data. The 444.95 Mb genome contained 35,488 protein-coding genes. A total of 20 chromosomes were assembled with a scaffold N50 of 21.75 Mb. The genome contained a cluster of strictosidine synthases and tryptophan decarboxylases with synteny to other species and a saccharide-terpene cluster involved in the monoterpenoid biosynthesis pathway of the MIA upstream pathway. The multi-omics data of A. scholaris provide a valuable resource for understanding the evolutionary origins of MIAs and for discovering biosynthetic pathways and synthetic biology efforts for producing pharmaceutically useful alkaloids.

wgd v2: a suite of tools to uncover and date ancient polyploidy and whole-genome duplication

MOTIVATION

Major improvements in sequencing technologies and genome sequence assembly have led to a huge increase in the number of available genome sequences. In turn, these genome sequences form an invaluable source for evolutionary, ecological, and comparative studies. One kind of analysis that has become routine is the search for traces of ancient polyploidy, particularly for plant genomes, where whole-genome duplication (WGD) is rampant.

RESULTS

Here, we present a major update of a previously developed tool wgd, namely wgd v2, to look for remnants of ancient polyploidy, or WGD. We implemented novel and improved previously developed tools to (a) construct KS age distributions for the whole-paranome (collection of all duplicated genes in a genome), (b) unravel intragenomic and intergenomic collinearity resulting from WGDs, (c) fit mixture models to age distributions of gene duplicates, (d) correct substitution rate variation for phylogenetic placement of WGDs, and (e) date ancient WGDs via phylogenetic dating of WGD-retained gene duplicates. The applicability and feasibility of wgd v2 for the identification and the relative and absolute dating of ancient WGDs is demonstrated using different plant genomes.

AVAILABILITY AND IMPLEMENTATION

wgd v2 is open source and available at https://github.com/heche-psb/wgd.

Phylogenomics and the rise of the angiosperms

Angiosperms are the cornerstone of most terrestrial ecosystems and human livelihoods1,2. A robust understanding of angiosperm evolution is required to explain their rise to ecological dominance. So far, the angiosperm tree of life has been determined primarily by means of analyses of the plastid genome3,4. Many studies have drawn on this foundational work, such as classification and first insights into angiosperm diversification since their Mesozoic origins5-7. However, the limited and biased sampling of both taxa and genomes undermines confidence in the tree and its implications. Here, we build the tree of life for almost 8,000 (about 60%) angiosperm genera using a standardized set of 353 nuclear genes8. This 15-fold increase in genus-level sampling relative to comparable nuclear studies9 provides a critical test of earlier results and brings notable change to key groups, especially in rosids, while substantiating many previously predicted relationships. Scaling this tree to time using 200 fossils, we discovered that early angiosperm evolution was characterized by high gene tree conflict and explosive diversification, giving rise to more than 80% of extant angiosperm orders. Steady diversification ensued through the remaining Mesozoic Era until rates resurged in the Cenozoic Era, concurrent with decreasing global temperatures and tightly linked with gene tree conflict. Taken together, our extensive sampling combined with advanced phylogenomic methods shows the deep history and full complexity in the evolution of a megadiverse clade.

Plastid phylogenomics and fossil evidence provide new insights into the evolutionary complexity of the 'woody clade' in Saxifragales

BACKGROUND

The "woody clade" in Saxifragales (WCS), encompassing four woody families (Altingiaceae, Cercidiphyllaceae, Daphniphyllaceae, and Hamamelidaceae), is a phylogenetically recalcitrant node in the angiosperm tree of life, as the interfamilial relationships of the WCS remain contentious. Based on a comprehensive sampling of WCS genera, this study aims to recover a robust maternal backbone phylogeny of the WCS by analyzing plastid genome (plastome) sequence data using Bayesian inference (BI), maximum likelihood (ML), and maximum parsimony (MP) methods, and to explore the possible causes of the phylogenetic recalcitrance with respect to deep relationships within the WCS, in combination with molecular and fossil evidence.

RESULTS

Although the four WCS families were identically resolved as monophyletic, the MP analysis recovered different tree topologies for the relationships among Altingiaceae, Cercidiphyllaceae, and Daphniphyllaceae from the ML and BI phylogenies. The fossil-calibrated plastome phylogeny showed that the WCS underwent a rapid divergence of crown groups in the early Cretaceous (between 104.79 and 100.23 Ma), leading to the origin of the stem lineage ancestors of Altingiaceae, Cercidiphyllaceae, Daphniphyllaceae, and Hamamelidaceae within a very short time span (∼4.56 Ma). Compared with the tree topology recovered in a previous study based on nuclear genome data, cytonuclear discordance regarding the interfamilial relationships of the WCS was detected.

CONCLUSIONS

Molecular and fossil evidence imply that the early divergence of the WCS might have experienced radiative diversification of crown groups, extensive extinctions at the genus and species levels around the Cretaceous/Paleocene boundary, and ancient hybridization. Such evolutionarily complex events may introduce biases in topological estimations within the WCS due to incomplete lineage sorting, cytonuclear discordance, and long-branch attraction, potentially impacting the accurate reconstruction of deep relationships.

Cretaceous chloranthoids: early prominence, extinct diversity and missing links

BACKGROUND

The Chloranthaceae comprise four extant genera (Hedyosmum, Ascarina, Chloranthus and Sarcandra), all with simple flowers. Molecular phylogenetics indicates that the Chloranthaceae diverged very early in angiosperm evolution, although how they are related to eudicots, magnoliids, monocots and Ceratophyllum is uncertain. Fossil pollen similar to that of Ascarina and Hedyosmum has long been recognized in the Early Cretaceous, but over the last four decades evidence of extinct Chloranthaceae based on other types of fossils has expanded dramatically and contributes significantly to understanding the evolution of the family.

SCOPE

Studies of fossils from the Cretaceous, especially mesofossils of Early Cretaceous age from Portugal and eastern North America, recognized diverse flowers, fruits, seeds, staminate inflorescences and stamens of extinct chloranthoids. These early chloranthoids include forms related to extant Hedyosmum and also to the Ascarina, Chloranthus and Sarcandra clade. In the Late Cretaceous there are several occurrences of distinctive fossil androecia related to extant Chloranthus. The rich and still expanding Cretaceous record of Chloranthaceae contrasts with a very sparse Cenozoic record, emphasizing that the four extant genera are likely to be relictual, although speciation within the genera might have occurred in relatively recent times. In this study, we describe three new genera of Early Cretaceous chloranthoids and summarize current knowledge on the extinct diversity of the group.

CONCLUSIONS

The evolutionary lineage that includes extant Chloranthaceae is diverse and abundantly represented in Early Cretaceous mesofossil floras that provide some of the earliest evidence of angiosperm reproductive structures. Extinct chloranthoids, some of which are clearly in the Chloranthaceae crown group, fill some of the morphological gaps that currently separate the extant genera, help to illuminate how some of the unusual features of extant Chloranthaceae evolved and suggest that Chloranthaceae are of disproportionate importance for a more refined understanding of ecology and phylogeny of early angiosperm diversification.

Chromosome-level genome assembly provides insights into the genome evolution and functional importance of the phenylpropanoid-flavonoid pathway in Thymus mongolicus

BACKGROUND

Thymus mongolicus (family Lamiaceae) is a Thyme subshrub with strong aroma and remarkable environmental adaptability. Limited genomic information limits the use of this plant.

RESULTS

Chromosome-level 605.2 Mb genome of T. mongolicus was generated, with 96.28% anchored to 12 pseudochromosomes. The repetitive sequences were dominant, accounting for 70.98%, and 32,593 protein-coding genes were predicted. Synteny analysis revealed that Lamiaceae species generally underwent two rounds of whole genome duplication; moreover, species-specific genome duplication was identified. A recent LTR retrotransposon burst and tandem duplication might play important roles in the formation of the Thymus genome. Using comparative genomic analysis, phylogenetic tree of seven Lamiaceae species was constructed, which revealed that Thyme plants evolved recently in the family. Under the phylogenetic framework, we performed functional enrichment analysis of the genes on nodes that contained the most gene duplication events (> 50% support) and of relevant significant expanded gene families. These genes were highly associated with environmental adaptation and biosynthesis of secondary metabolites. Combined transcriptome and metabolome analyses revealed that Peroxidases, Hydroxycinnamoyl-CoA shikimate/quinate hydroxycinnamoyl transferases, and 4-coumarate-CoA ligases genes were the essential regulators of the phenylpropanoid-flavonoid pathway. Their catalytic products (e.g., apigenin, naringenin chalcone, and several apigenin-related compounds) might be responsible for the environmental tolerance and aromatic properties of T. mongolicus.

CONCLUSION

This study enhanced the understanding of the genomic evolution of T. mongolicus, enabling further exploration of its unique traits and applications, and contributed to the understanding of Lamiaceae genomics and evolutionary biology.

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.

Widespread incomplete lineage sorting and introgression shaped adaptive radiation in the Gossypium genus

Cotton (Gossypium) stands as a crucial economic crop, serving as the primary source of natural fiber for the textile sector. However, the evolutionary mechanisms driving speciation within the Gossypium genus remain unresolved. In this investigation, we leveraged 25 Gossypium genomes and introduced four novel assemblies-G. harknessii, G. gossypioides, G. trilobum, and G. klotzschianum (Gklo)-to delve into the speciation history of this genus. Notably, we encountered intricate phylogenies potentially stemming from introgression. These complexities are further compounded by incomplete lineage sorting (ILS), a factor likely to have been instrumental in shaping the swift diversification of cotton. Our focus subsequently shifted to the rapid radiation episode during a concise period in Gossypium evolution. For a recently diverged lineage comprising G. davidsonii, Gklo, and G. raimondii, we constructed a finely detailed ILS map. Intriguingly, this analysis revealed the non-random distribution of ILS regions across the reference Gklo genome. Moreover, we identified signs of robust natural selection influencing specific ILS regions. Noteworthy variations pertaining to speciation emerged between the closely related sister species Gklo and G. davidsonii. Approximately 15.74% of speciation structural variation genes and 12.04% of speciation-associated genes were estimated to intersect with ILS signatures. These findings enrich our understanding of the role of ILS in adaptive radiation, shedding fresh light on the intricate speciation history of the Gossypium genus.

Dysfunction of duplicated pair rice histone acetyltransferases causes segregation distortion and an interspecific reproductive barrier

Postzygotic reproductive isolation, which results in the irreversible divergence of species, is commonly accompanied by hybrid sterility, necrosis/weakness, or lethality in the F1 or other offspring generations. Here we show that the loss of function of HWS1 and HWS2, a couple of duplicated paralogs, together confer complete interspecific incompatibility between Asian and African rice. Both of these non-Mendelian determinants encode the putative Esa1-associated factor 6 (EAF6) protein, which functions as a characteristic subunit of the histone H4 acetyltransferase complex regulating transcriptional activation via genome-wide histone modification. The proliferating tapetum and inappropriate polar nuclei arrangement cause defective pollen and seeds in F2 hybrid offspring due to the recombinant HWS1/2-mediated misregulation of vitamin (biotin and thiamine) metabolism and lipid synthesis. Evolutionary analysis of HWS1/2 suggests that this gene pair has undergone incomplete lineage sorting (ILS) and multiple gene duplication events during speciation. Our findings have not only uncovered a pair of speciation genes that control hybrid breakdown but also illustrate a passive mechanism that could be scaled up and used in the guidance and optimization of hybrid breeding applications for distant hybridization.

Characterization of the Complete Mitochondrial Genome of Wintersweet (Chimonanthus praecox) and Comparative Analysis within Magnoliids

Mitochondrial genome sequencing is a valuable tool for investigating mitogenome evolution, species phylogeny, and population genetics. Chimonanthus praecox (L.) Link, also known as "La Mei" in Chinese, is a famous ornamental and medical shrub belonging to the order Laurales of the Calycanthaceae family. Although the nuclear genomes and chloroplast genomes of certain Laurales representatives, such as Lindera glauca, Laurus nobilis, and Piper nigrum, have been sequenced, the mitochondrial genome of Laurales members remains unknown. Here, we reported the first complete mitogenome of C. praecox. The mitogenome was 972,347 bp in length and comprised 60 unique coding genes, including 40 protein-coding genes (PCGs), 17 tRNA genes, and three rRNA genes. The skewness of the PCGs showed that the AT skew (-0.0096233) was negative, while the GC skew (0.031656) was positive, indicating higher contents of T's and G's in the mitochondrial genome of C. praecox. The Ka/Ks ratio analysis showed that the Ka/Ks values of most genes were less than one, suggesting that these genes were under purifying selection. Furthermore, there is a substantial abundance of dispersed repeats in C. praecox, constituting 16.98% of the total mitochondrial genome. A total of 731 SSR repeats were identified in the mitogenome, the highest number among the eleven available magnoliids mitogenomes. The mitochondrial phylogenetic analysis based on 29 conserved PCGs placed the C. praecox in Lauraceae, and supported the sister relationship of Laurales with Magnoliales, which was congruent with the nuclear genome evidence. The present study enriches the mitogenome data of C. praecox and promotes further studies on phylogeny and plastid evolution.

Chemistry and bioactivity of lindenane sesquiterpenoids and their oligomers

Covering: 1925 to July 2023Among the sesquiterpenoids with rich structural diversity and potential bioactivities, lindenane sesquiterpenoids (LSs) possess a characteristic cis, trans-3,5,6-carbocyclic skeleton and mainly exist as monomers and diverse oligomers in plants from the Lindera genus and Chloranthaceae family. Since the first identification of lindeneol from Lindera strychnifolia in 1925, 354 natural LSs and their oligomers with anti-inflammatory, antitumor, and anti-infective activities have been discovered. Structurally, two-thirds of LSs exist as oligomers with interesting skeletons through diverse polymeric patterns, especially Diels-Alder [4 + 2] cycloaddition. Fascinated by their diverse bioactivities and intriguing polycyclic architectures, synthetic chemists have engaged in the total synthesis of natural LSs in recent decades. In this review, the research achievements related to LSs from 1925 to July of 2023 are systematically and comprehensively summarized, focusing on the classification of their structures, chemical synthesis, and bioactivities, which will be helpful for further research on LSs and their oligomers.

The chromosome-scale genome of Magnolia sinica (Magnoliaceae) provides insights into the conservation of plant species with extremely small populations (PSESP)

Magnolia sinica (Magnoliaceae) is a highly threatened tree endemic to southeast Yunnan, China. In this study, we generated for the first time a high-quality chromosome-scale genome sequence from M. sinica, by combining Illumina and ONT data with Hi-C scaffolding methods. The final assembled genome size of M. sinica was 1.84 Gb, with a contig N50 of ca. 45 Mb and scaffold N50 of 92 Mb. Identified repeats constituted approximately 57% of the genome, and 43,473 protein-coding genes were predicted. Phylogenetic analysis shows that the magnolias form a sister clade with the eudicots and the order Ceratophyllales, while the monocots are sister to the other core angiosperms. In our study, a total of 21 individuals from the 5 remnant populations of M. sinica, as well as 22 specimens belonging to 8 related Magnoliaceae species, were resequenced. The results showed that M. sinica had higher genetic diversity (θw = 0.01126 and θπ = 0.01158) than other related species in the Magnoliaceae. However, population structure analysis suggested that the genetic differentiation among the 5 M. sinica populations was very low. Analyses of the demographic history of the species using different models consistently revealed that 2 bottleneck events occurred. The contemporary effective population size of M. sinica was estimated to be 10.9. The different patterns of genetic loads (inbreeding and numbers of deleterious mutations) suggested constructive strategies for the conservation of these 5 different populations of M. sinica. Overall, this high-quality genome will be a valuable genomic resource for conservation of M. sinica.

Chromosome-scale genomes of five Hongmu species in Leguminosae

The Legume family (Leguminosae or Fabaceae), is one of the largest and economically important flowering plants. Heartwood, the core of a tree trunk or branch, is a valuable and renewable resource employed for centuries in constructing sturdy and sustainable structures. Hongmu refers to a category of precious timber trees in China, encompassing 29 woody species, primarily from the legume genus. Due to the lack of genome data, detailed studies on their economic and ecological importance are limited. Therefore, this study generates chromosome-scale assemblies of five Hongmu species in Leguminosae: Pterocarpus santalinus, Pterocarpus macrocarpus, Dalbergia cochinchinensis, Dalbergia cultrata, and Senna siamea, using a combination of short-reads, long-read nanopore, and Hi-C data. We obtained 623.86 Mb, 634.58 Mb, 700.60 Mb, 645.98 Mb, and 437.29 Mb of pseudochromosome level assemblies with the scaffold N50 lengths of 63.1 Mb, 63.7 Mb, 70.4 Mb, 61.1 Mb and 32.2 Mb for P. santalinus, P. macrocarpus, D. cochinchinensis, D. cultrata and S. siamea, respectively. These genome data will serve as a valuable resource for studying crucial traits, like wood quality, disease resistance, and environmental adaptation in Hongmu.

Scaphopoda is the sister taxon to Bivalvia: Evidence of ancient incomplete lineage sorting

The almost simultaneous emergence of major animal phyla during the early Cambrian shaped modern animal biodiversity. Reconstructing evolutionary relationships among such closely spaced branches in the animal tree of life has proven to be a major challenge, hindering understanding of early animal evolution and the fossil record. This is particularly true in the species-rich and highly varied Mollusca where dramatic inconsistency among paleontological, morphological, and molecular evidence has led to a long-standing debate about the group's phylogeny and the nature of dozens of enigmatic fossil taxa. A critical step needed to overcome this issue is to supplement available genomic data, which is plentiful for well-studied lineages, with genomes from rare but key lineages, such as Scaphopoda. Here, by presenting chromosome-level genomes from both extant scaphopod orders and leveraging complete genomes spanning Mollusca, we provide strong support for Scaphopoda as the sister taxon of Bivalvia, revitalizing the morphology-based Diasoma hypothesis originally proposed 50 years ago. Our molecular clock analysis confidently dates the split between Bivalvia and Scaphopoda at ~520 Ma, prompting a reinterpretation of controversial laterally compressed Early Cambrian fossils, including Anabarella, Watsonella, and Mellopegma, as stem diasomes. Moreover, we show that incongruence in the phylogenetic placement of Scaphopoda in previous phylogenomic studies was due to ancient incomplete lineage sorting (ILS) that occurred during the rapid radiation of Conchifera. Our findings highlight the need to consider ILS as a potential source of error in deep phylogeny reconstruction, especially in the context of the unique nature of the Cambrian Explosion.

Is there silicon in flowers and what does it tell us?

The emergence of flowers marked an important development in plant evolution. Flowers in many species evolved to attract animal pollinators to increase fertilisation chances. In leaves, silicon (Si) discourages herbivores, for example by wearing down mouthparts. Flowers are essentially modified leaves and hence may also have the capacity to accumulate Si. If Si in flowers discourages animal visitors as it does in leaves, Si accumulation may be disadvantageous for pollination. Whether flowers accumulate Si, and what the implications may be, was not known for many species. We analysed leaves and flowers of different taxa, separated into their different anatomical parts. Flowers mostly have low Si concentrations in all parts (mean ± SE of BSi in mg g-1 was 0.22 ± 0.04 in petals, 0.59 ± 0.24 in sepals, 0.14 ± 0.03 in stamens, 0.15 ± 0.04 in styles and stigmas and 0.37 ± 0.19 in ovaries for a subset of 56 species). In most cases, less Si was accumulated in flowers than in leaves (mean ± SE of BSi in mg g-1 was 1.51 ± 0.55 in whole flowers vs. 2.97 ± 0.57 in leaves in 104 species) though intriguing exceptions are found, with some species accumulating more Si in flowers than leaves. The large variation in concentration among flowers across the taxa examined, with a particularly high concentration in grass inflorescences, tantalisingly suggests differences in the use of Si for flowers across plant groups. We conclude that the study of the functions of Si for flowers warrants more attention, with pollination strategy a potential contributing factor.

Comparing genomes of Fructus Amomi-producing species reveals genetic basis of volatile terpenoid divergence

Wurfbainia longiligularis and Wurfbainia villosa are both rich in volatile terpenoids and are 2 primary plant sources of Fructus Amomi used for curing gastrointestinal diseases. Metabolomic profiling has demonstrated that bornyl diphosphate (BPP)-related terpenoids are more abundant in the W. villosa seeds and have a wider tissue distribution in W. longiligularis. To explore the genetic mechanisms underlying the volatile terpenoid divergence, a high-quality chromosome-level genome of W. longiligularis (2.29 Gb, contig N50 of 80.39 Mb) was assembled. Functional characterization of 17 terpene synthases (WlTPSs) revealed that WlBPPS, along with WlTPS 24/26/28 with bornyl diphosphate synthase (BPPS) activity, contributes to the wider tissue distribution of BPP-related terpenoids in W. longiligularis compared to W. villosa. Furthermore, transgenic Nicotiana tabacum showed that the GCN4-motif element positively regulates seed expression of WvBPPS and thus promotes the enrichment of BPP-related terpenoids in W. villosa seeds. Systematic identification and analysis of candidate TPS in 29 monocot plants from 16 families indicated that substantial expansion of TPS-a and TPS-b subfamily genes in Zingiberaceae may have driven increased diversity and production of volatile terpenoids. Evolutionary analysis and functional identification of BPPS genes showed that BPP-related terpenoids may be distributed only in the Zingiberaceae of monocot plants. This research provides valuable genomic resources for breeding and improving Fructus Amomi with medicinal and edible value and sheds light on the evolution of terpenoid biosynthesis in Zingiberaceae.

Bioinformatics Analysis of MSH1 Genes of Green Plants: Multiple Parallel Length Expansions, Intron Gains and Losses, Partial Gene Duplications, and Alternative Splicing

MutS homolog 1 (MSH1) is involved in the recombining and repairing of organelle genomes and is essential for maintaining their stability. Previous studies indicated that the length of the gene varied greatly among species and detected species-specific partial gene duplications in Physcomitrella patens. However, there are critical gaps in the understanding of the gene size expansion, and the extent of the partial gene duplication of MSH1 remains unclear. Here, we screened MSH1 genes in 85 selected species with genome sequences representing the main clades of green plants (Viridiplantae). We identified the MSH1 gene in all lineages of green plants, except for nine incomplete species, for bioinformatics analysis. The gene is a singleton gene in most of the selected species with conserved amino acids and protein domains. Gene length varies greatly among the species, ranging from 3234 bp in Ostreococcus tauri to 805,861 bp in Cycas panzhihuaensis. The expansion of MSH1 repeatedly occurred in multiple clades, especially in Gymnosperms, Orchidaceae, and Chloranthus spicatus. MSH1 has exceptionally long introns in certain species due to the gene length expansion, and the longest intron even reaches 101,025 bp. And the gene length is positively correlated with the proportion of the transposable elements (TEs) in the introns. In addition, gene structure analysis indicated that the MSH1 of green plants had undergone parallel intron gains and losses in all major lineages. However, the intron number of seed plants (gymnosperm and angiosperm) is relatively stable. All the selected gymnosperms contain 22 introns except for Gnetum montanum and Welwitschia mirabilis, while all the selected angiosperm species preserve 21 introns except for the ANA grade. Notably, the coding region of MSH1 in algae presents an exceptionally high GC content (47.7% to 75.5%). Moreover, over one-third of the selected species contain species-specific partial gene duplications of MSH1, except for the conserved mosses-specific partial gene duplication. Additionally, we found conserved alternatively spliced MSH1 transcripts in five species. The study of MSH1 sheds light on the evolution of the long genes of green plants.

A high-quality Bougainvillea genome provides new insights into evolutionary history and pigment biosynthetic pathways in the Caryophyllales

Bougainvillea is a perennial ornamental shrub that is highly regarded in ornamental horticulture around the world. However, the absence of genome data limits our understanding of the pathways involved in bract coloration and breeding. Here, we report a chromosome-level assembly of the giga-genome of Bougainvillea × buttiana 'Mrs Butt', a cultivar thought to be the origin of many other Bougainvillea cultivars. The assembled genome is ~5 Gb with a scaffold N50 of 151 756 278 bp and contains 86 572 genes which have undergone recent whole-genome duplication. We confirmed that multiple rounds of whole-genome multiplication have occurred in the evolutionary history of the Caryophyllales, reconstructed the relationship in the Caryophyllales at whole genome level, and found discordance between species and gene trees as the result of complex introgression events. We investigated betalain and anthocyanin biosynthetic pathways and found instances of independent evolutionary innovations in the nine different Caryophyllales species. To explore the potential formation mechanism of diverse bract colors in Bougainvillea, we analyzed the genes involved in betalain and anthocyanin biosynthesis and found extremely low expression of ANS and DFR genes in all cultivars, which may limit anthocyanin biosynthesis. Our findings indicate that the expression pattern of the betalain biosynthetic pathway did not directly correlate with bract color, and a higher expression level in the betalain biosynthetic pathway is required for colored bracts. This improved understanding of the correlation between gene expression and bract color allows plant breeding outcomes to be predicted with greater certainty.

The genome of Acorus deciphers insights into early monocot evolution

Acorales is the sister lineage to all the other extant monocot plants. Genomic resource enhancement of this genus can help to reveal early monocot genomic architecture and evolution. Here, we assemble the genome of Acorus gramineus and reveal that it has ~45% fewer genes than the majority of monocots, although they have similar genome size. Phylogenetic analyses based on both chloroplast and nuclear genes consistently support that A. gramineus is the sister to the remaining monocots. In addition, we assemble a 2.2 Mb mitochondrial genome and observe many genes exhibit higher mutation rates than that of most angiosperms, which could be the reason leading to the controversies of nuclear genes- and mitochondrial genes-based phylogenetic trees existing in the literature. Further, Acorales did not experience tau (τ) whole-genome duplication, unlike majority of monocot clades, and no large-scale gene expansion is observed. Moreover, we identify gene contractions and expansions likely linking to plant architecture, stress resistance, light harvesting, and essential oil metabolism. These findings shed light on the evolution of early monocots and genomic footprints of wetland plant adaptations.

Diploid and tetraploid genomes of Acorus and the evolution of monocots

Monocots are a major taxon within flowering plants, have unique morphological traits, and show an extraordinary diversity in lifestyle. To improve our understanding of monocot origin and evolution, we generate chromosome-level reference genomes of the diploid Acorus gramineus and the tetraploid Ac. calamus, the only two accepted species from the family Acoraceae, which form a sister lineage to all other monocots. Comparing the genomes of Ac. gramineus and Ac. calamus, we suggest that Ac. gramineus is not a potential diploid progenitor of Ac. calamus, and Ac. calamus is an allotetraploid with two subgenomes A, and B, presenting asymmetric evolution and B subgenome dominance. Both the diploid genome of Ac. gramineus and the subgenomes A and B of Ac. calamus show clear evidence of whole-genome duplication (WGD), but Acoraceae does not seem to share an older WGD that is shared by most other monocots. We reconstruct an ancestral monocot karyotype and gene toolkit, and discuss scenarios that explain the complex history of the Acorus genome. Our analyses show that the ancestors of monocots exhibit mosaic genomic features, likely important for that appeared in early monocot evolution, providing fundamental insights into the origin, evolution, and diversification of monocots.

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.

A plastid phylogenomic framework for the palm family (Arecaceae)

BACKGROUND

Over the past decade, phylogenomics has greatly advanced our knowledge of angiosperm evolution. However, phylogenomic studies of large angiosperm families with complete species or genus-level sampling are still lacking. The palms, Arecaceae, are a large family with ca. 181 genera and 2600 species and are important components of tropical rainforests bearing great cultural and economic significance. Taxonomy and phylogeny of the family have been extensively investigated by a series of molecular phylogenetic studies in the last two decades. Nevertheless, some phylogenetic relationships within the family are not yet well-resolved, especially at the tribal and generic levels, with consequent impacts for downstream research.

RESULTS

Plastomes of 182 palm species representing 111 genera were newly sequenced. Combining these with previously published plastid DNA data, we were able to sample 98% of palm genera and conduct a plastid phylogenomic investigation of the family. Maximum likelihood analyses yielded a robustly supported phylogenetic hypothesis. Phylogenetic relationships among all five palm subfamilies and 28 tribes were well-resolved, and most inter-generic phylogenetic relationships were also resolved with strong support.

CONCLUSIONS

The inclusion of nearly complete generic-level sampling coupled with nearly complete plastid genomes strengthened our understanding of plastid-based relationships of the palms. This comprehensive plastid genome dataset complements a growing body of nuclear genomic data. Together, these datasets form a novel phylogenomic baseline for the palms and an increasingly robust framework for future comparative biological studies of this exceptionally important plant family.

The origin and early evolution of plants

Plant (archaeplastid) evolution has transformed the biosphere, but we are only now beginning to learn how this took place through comparative genomics, phylogenetics, and the fossil record. This has illuminated the phylogeny of Archaeplastida, Viridiplantae, and Streptophyta, and has resolved the evolution of key characters, genes, and genomes - revealing that many key innovations evolved long before the clades with which they have been casually associated. Molecular clock analyses estimate that Streptophyta and Viridiplantae emerged in the late Mesoproterozoic to late Neoproterozoic, whereas Archaeplastida emerged in the late-mid Palaeoproterozoic. Together, these insights inform on the coevolution of plants and the Earth system that transformed ecology and global biogeochemical cycles, increased weathering, and precipitated snowball Earth events, during which they would have been key to oxygen production and net primary productivity (NPP).

The Saururus chinensis genome provides insights into the evolution of pollination strategies and herbaceousness in magnoliids

Saururus chinensis, an herbaceous magnoliid without perianth, represents a clade of early-diverging angiosperms that have gone through woodiness-herbaceousness transition and pollination obstacles: the characteristic white leaves underneath inflorescence during flowering time are considered a substitute for perianth to attract insect pollinators. Here, using the newly sequenced S. chinensis genome, we revisited the phylogenetic position of magnoliids within mesangiosperms, and recovered a sister relationship for magnoliids and Chloranthales. By considering differentially expressed genes, we identified candidate genes that are involved in the morphogenesis of the white leaves in S. chinensis. Among those genes, we verified - in a transgenic experiment with Arabidopsis - that increasing the expression of the "pseudo-etiolation in light" gene (ScPEL) can inhibit the biosynthesis of chlorophyll. ScPEL is thus likely responsible for the switches between green and white leaves, suggesting that changes in gene expression may underlie the evolution of pollination strategies. Despite being an herbaceous plant, S. chinensis still has vascular cambium and maintains the potential for secondary growth as a woody plant, because the necessary machinery, i.e., the entire gene set involved in lignin biosynthesis, is well preserved. However, similar expression levels of two key genes (CCR and CAD) between the stem and other tissues in the lignin biosynthesis pathway are possibly associated with the herbaceous nature of S. chinensis. In conclusion, the S. chinensis genome provides valuable insights into the adaptive evolution of pollination in Saururaceae and reveals a possible mechanism for the evolution of herbaceousness in magnoliids.

Phylogenetic studies of magnoliids: Advances and perspectives

Magnoliids are the largest flowering plant clades outside of the eudicots and monocots, which are distributed worldwide and have high economic, ornamental and ecological values. Eudicots, monocots and magnoliids are the three major clades of Mesangiospermae, and their phylogenetic relationship is one of the most interesting issues. In recent years, with the continuous accumulation of genomic information, the evolutionary status of magnoliids has become a hot spot in plant phylogenetic research. Although great efforts have been made to study the evolution of magnoliids using molecular data from several representative species such as nuclear genome, plastid genome, mitochondrial genome, and transcriptome, the results of current studies on the phylogenetic status of magnoliids are inconsistent. Here, we systematically describe the current understanding of the molecular research on magnoliid phylogeny and review the differences in the evolutionary state of magnoliids. Understanding the research approaches and limitations of magnoliid phylogeny can guide research strategies to further improve the study of the phylogenetic evolution of magnoliids.

Genome assembly of Musa beccarii shows extensive chromosomal rearrangements and genome expansion during evolution of Musaceae genomes

BACKGROUND

Musa beccarii (Musaceae) is a banana species native to Borneo, sometimes grown as an ornamental plant. The basic chromosome number of Musa species is x = 7, 10, or 11; however, M. beccarii has a basic chromosome number of x = 9 (2n = 2x = 18), which is the same basic chromosome number of species in the sister genera Ensete and Musella. Musa beccarii is in the section Callimusa, which is sister to the section Musa. We generated a high-quality chromosome-scale genome assembly of M. beccarii to better understand the evolution and diversity of genomes within the family Musaceae.

FINDINGS

The M. beccarii genome was assembled by long-read and Hi-C sequencing, and genes were annotated using both long Iso-seq and short RNA-seq reads. The size of M. beccarii was the largest among all known Musaceae assemblies (∼570 Mbp) due to the expansion of transposable elements and increased 45S ribosomal DNA sites. By synteny analysis, we detected extensive genome-wide chromosome fusions and fissions between M. beccarii and the other Musa and Ensete species, far beyond those expected from differences in chromosome number. Within Musaceae, M. beccarii showed a reduced number of terpenoid synthase genes, which are related to chemical defense, and enrichment in lipid metabolism genes linked to the physical defense of the cell wall. Furthermore, type III polyketide synthase was the most abundant biosynthetic gene cluster (BGC) in M. beccarii. BGCs were not conserved in Musaceae genomes.

CONCLUSIONS

The genome assembly of M. beccarii is the first chromosome-scale genome assembly in the Callimusa section in Musa, which provides an important genetic resource that aids our understanding of the evolution of Musaceae genomes and enhances our knowledge of the pangenome.

A draft genome of the medicinal plant Cremastra appendiculata (D. Don) provides insights into the colchicine biosynthetic pathway

Cremastra appendiculata (D. Don) Makino is a rare terrestrial orchid with a high market value as an ornamental and Chinese traditional medicinal herb with a wide range of pharmacological properties. The pseudobulbs of C. appendiculata are one of the primary sources of the famous traditional Chinese medicine "Shancigu", which has been clinically used for treating many diseases, especially, as the main component to treat gout. The lack of genetic research and genome data restricts the modern development and clinical use of C. appendiculata. Here, we report a 2.3 Gb chromosome-level genome of C. appendiculata. We identify a series of candidates of 35 candidate genes responsible for colchicine biosynthesis, among which O-methyltransferase (OMT) gene exhibits an important role in colchicine biosynthesis. Co-expression analysis reveal purple and green-yellow module have close relationships with pseudobulb parts and comprise most of the colchicine pathway genes. Overall, our genome data and the candidate genes reported here set the foundation to decipher the colchicine biosynthesis pathways in medicinal plants.

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.

Shoot Development in Members of an Ancient Aquatic Angiosperm Lineage, Ceratophyllaceae: A New Interpretation Facilitates Comparisons with Chloranthaceae

Ceratophyllum is an ancient and phylogenetically isolated angiosperm lineage. Comparisons between Ceratophyllum and other angiosperms are hampered by uncertainty in inferring organ homologies in this genus of specialized aquatics. Interpretation of shoot morphology is especially problematic in Ceratophyllum. Each node has several leaf-like appendages interpreted as verticillate leaves, modified parts of one and the same leaf or parts of two leaves under decussate phyllotaxis. Vegetative branches are axillary, but reproductive units (interpreted as flowers or inflorescences) are commonly viewed as developing from collateral accessory buds. We studied shoot development in Ceratophyllum submersum, C. tanaiticum, and C. demersum using scanning electron microscopy to clarify shoot morphology and branching patterns. Our data support the idea that the phyllotaxis is essentially decussate with appendages of stipular origin resembling leaf blades. We conclude that a leaf axil of Ceratophyllum possesses a complex of two serial buds, the lower one producing a vegetative branch and the upper one developing a reproductive unit. The reproductive unit is congenitally displaced to the subsequent node, a phenomenon known as concaulescence. Either member of the serial bud complex may be absent. There is a theory based on a synthesis of molecular and morphological data that Chloranthaceae are the closest extant relatives of Ceratophyllum. Serial buds and concaulescence are known in Hedyosmum (Chloranthaceae). Our new interpretation facilitates morphological comparisons between Hedyosmum and Ceratophyllum.

Refined Interpretation of the Pistillate Flower in Ceratophyllum Sheds Fresh Light on Gynoecium Evolution in Angiosperms

Molecular phylogenetic analyses have revealed a superclade of mesangiosperms with five extant lineages: monocots, eudicots, magnoliids, Ceratophyllum and Chloranthaceae. Both Ceratophyllum and Chloranthaceae are ancient lineages with a long fossil record; their precise placement within mesangiosperms is uncertain. Morphological studies have suggested that they form a clade together with some Cretaceous fossils, including Canrightia, Montsechia and Pseudoasterophyllites. Apart from Canrightia, members of this clade share unilocular gynoecia commonly interpreted as monomerous with ascidiate carpels. Alternatively, the gynoecium of Ceratophyllum has also been interpreted as syncarpous with a single fertile carpel (pseudomonomerous). We investigate patterns of morphological, anatomical and developmental variation in gynoecia of three Ceratophyllum species to explore the controversial interpretation of its gynoecium as either monomerous or pseudomonomerous. We use an angiosperm-wide morphological data set and contrasting tree topologies to estimate the ancestral gynoecium type in both Ceratophyllum and mesangiosperms. Gynoecia of all three Ceratophyllum species possess a small (sometimes vestigial) glandular appendage on the abaxial side and an occasionally bifurcating apex. The ovary is usually unilocular with two procambium strands, but sometimes bilocular and/or with three strands in C. demersum. None of the possible phylogenetic placements strongly suggest apocarpy in the stem lineage of Ceratophyllum. Rescoring Ceratophyllum as having two united carpels affects broader-scale reconstructions of the ancestral gynoecium in mesangiosperms. Our interpretation of the glandular appendage as a tepal or staminode homologue makes the Ceratophyllum ovary inferior, thus resembling (semi)inferior ovaries of most Chloranthaceae and potentially related fossils Canrightia and Zlatkocarpus. The entire structure of the flower of Ceratophyllum suggests strong reduction following a long and complex evolutionary history. The widely accepted notion that apocarpy is ancestral in mesangiosperms (and angiosperms) lacks robust support, regardless of which modes of carpel fusion are considered. Our study highlights the crucial importance of incorporating fossils into large-scale analyses to understand character evolution.

A High-Quality Reference Genome Sequence and Genetic Transformation System of Aralia elata

Aralia elata is a perennial woody plant of the genus Aralia in the family Araliaceae. It is rich in saponins and therefore has a wide range of pharmacological effects. Here, we report a high-quality reference genome of A. elata, with a genome size of 1.21 Gb and a contig N50 of 51.34 Mb, produced by PacBio HiFi sequencing technology. This is the first genome assembly for the genus Aralia. Through genome evolutionary analysis, we explored the phylogeny and whole genome duplication (WGD) events in the A. elata genome. The results indicated that a recent WGD event occurred in the A. elata genome. Estimation of the divergence times indicated that the WGD may be shared by Araliaceae. By analyzing the genome sequence of A. elata and combining the transcriptome data from three tissues, we discovered important genes related to triterpene saponins biosynthesis. Furthermore, based on the embryonic callus induction system of A. elata established in our laboratory, we set up the genetic transformation system of this plant. The genomic resources and genetic transformation system obtained in this study provide insights into A. elata and lays the foundation for further exploration of the A. elata regulatory mechanism.

Comparative Analyses of 3,654 Plastid Genomes Unravel Insights Into Evolutionary Dynamics and Phylogenetic Discordance of Green Plants

The plastid organelle is essential for many vital cellular processes and the growth and development of plants. The availability of a large number of complete plastid genomes could be effectively utilized to understand the evolution of the plastid genomes and phylogenetic relationships among plants. We comprehensively analyzed the plastid genomes of Viridiplantae comprising 3,654 taxa from 298 families and 111 orders and compared the genomic organizations in their plastid genomic DNA among major clades, which include gene gain/loss, gene copy number, GC content, and gene blocks. We discovered that some important genes that exhibit similar functions likely formed gene blocks, such as the psb family presumably showing co-occurrence and forming gene blocks in Viridiplantae. The inverted repeats (IRs) in plastid genomes have doubled in size across land plants, and their GC content is substantially higher than non-IR genes. By employing three different data sets [all nucleotide positions (nt123), only the first and second codon positions (nt12), and amino acids (AA)], our phylogenomic analyses revealed Chlorokybales + Mesostigmatales as the earliest-branching lineage of streptophytes. Hornworts, mosses, and liverworts forming a monophylum were identified as the sister lineage of tracheophytes. Based on nt12 and AA data sets, monocots, Chloranthales and magnoliids are successive sister lineages to the eudicots + Ceratophyllales clade. The comprehensive taxon sampling and analysis of different data sets from plastid genomes recovered well-supported relationships of green plants, thereby contributing to resolving some long-standing uncertainties in the plant phylogeny.