Transcriptome sequencing data provides a solid base to understand the phylogenetic relationships, biogeography and reticulated evolution of the genus Zamia L. (Cycadales: Zamiaceae).

BACKGROUND AND AIMS: Cycad is a key lineage to understand the early evolution of seed plants and their response to past environmental changes. However, tracing the evolutionary trajectory of cycad species is challenging when the robust relationships at inter- or infrageneric level are not well resolved. METHODS: Here, using 2,901 single-copy nuclear genes, we explored the species relationships and gene flow within the second largest genus of cycads, i.e., Zamia, based on phylotranscriptomic analyses of 90% extant Zamia species. Based on a well-resolved phylogenetic framework, we performed gene flow analyses, molecular dating, and biogeographical reconstruction to examine the spatiotemporal evolution of Zamia. We also performed ancestral state reconstruction (ASR) of a total of 62 traits of the genus to comprehensively investigate its morphological evolution. KEY RESULTS: Zamia is comprised of seven major clades corresponding to seven distinct distribution areas in the Americas, with at least three reticulation nodes revealed in this genus. Extant lineages of Zamia initially diversified around 18.4-32.6 (29.14) million years ago (MA) in the Mega-Mexico, and then expanded eastward into the Caribbean and southward into Central and South America. ASR revealed homoplasy in most of the morphological characters. CONCLUSIONS: This study revealed congruent phylogenetic relationships from comparative methods/datasets, with some conflicts being the result of incomplete lineage sorting and ancient/recent hybridization events. The strong association between the clades and the biogeographic areas suggested that ancient dispersal events shaped the modern distribution pattern, and regional climatic factors may have resulted in the following in-situ diversification. Climate cooling starting during the mid Miocene is associated with the global expansion of Zamia to the tropical South America that have dramatically driven lineage diversification in the New World flora, as well as the extinction of cycad species in the nowadays cooler regions of both hemispheres as indicated by the fossil records.

Eco-evolutionary evidence for the global diversity pattern of Cycas (Cycadaceae).

The evolution of the latitudinal diversity gradient (LDG), characterized by a peak in diversity toward the tropics, has captured significant attention in evolutionary biology and ecology. However, the inverse LDG (i-LDG) mechanism, wherein species richness increases toward the poles, remains inadequately explored. Cycads are among one of the oldest lineages of extant seed plants and have undergone extensive diversification in the tropics. Intriguingly, the extant cycad abundance exhibits an i-LDG pattern, and the underlying causes for this phenomenon remain largely elusive. Here, using 1,843 nuclear genes from a nearly complete sampling, we conducted comprehensive phylogenomic analyses to establish a robust species-level phylogeny for Cycas, the largest genus within cycads. We then reconstructed the spatial-temporal dynamics and integrated global environmental data to evaluate the roles of species ages, diversification rates, contemporary environment, and conservatism to ancestral niches in shaping the i-LDG pattern. We found Cycas experienced decreased diversification rates, coupled with the cooling temperature since its origin in the Eocene from continental Asia. Different regions have distinctively contributed to the formation of i-LDG for Cycas, with the northern hemisphere acting as evolutionary museums and the southern hemisphere serving as cradles. Moreover, water-related climate variables, specifically precipitation seasonality and potential evapotranspiration, were identified as paramount factors constraining Cycas species richness in the rainforest biome near the equator. Notably, the adherence to ancestral monsoonal climates emerges as a critical factor in sustaining the diversity pattern. This study underscores the imperative of integrating both evolutionary and ecological approaches to comprehensively unravel the mechanisms underpinning global biodiversity patterns.

Chromosome-scale genomes of commercially important mahoganies, Swietenia macrophylla and Khaya senegalensis.

Mahogany species (family Meliaceae) are highly valued for their aesthetic and durable wood. Despite their economic and ecological importance, genomic resources for mahogany species are limited, hindering genetic improvement and conservation efforts. Here we perform chromosome-scale genome assemblies of two commercially important mahogany species: Swietenia macrophylla and Khaya senegalensis. By combining 10X sequencing and Hi-C data, we assemble high-quality genomes of 274.49 Mb (S. macrophylla) and 406.50 Mb (K. senegalensis), with scaffold N50 lengths of 8.51 Mb and 7.85 Mb, respectively. A total of 99.38% and 98.05% of the assembled sequences are anchored to 28 pseudo-chromosomes in S. macrophylla and K. senegalensis, respectively. We predict 34,129 and 31,908 protein-coding genes in S. macrophylla and K. senegalensis, respectively, of which 97.44% and 98.49% are functionally annotated. The chromosome-scale genome assemblies of these mahogany species could serve as a vital genetic resource, especially in understanding the properties of non-model woody plants. These high-quality genomes could support the development of molecular markers for breeding programs, conservation efforts, and the sustainable management of these valuable forest resources.

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.

Phylotranscriptomics Shed Light on Intrageneric Relationships and Historical Biogeography of Ceratozamia (Cycadales).

Ceratozamia Brongn. is one of the species-rich genera of Cycadales comprising 38 species that are mainly distributed in Mexico, with a few species reported from neighboring regions. Phylogenetic relationships within the genus need detailed investigation based on extensive datasets and reliable systematic approaches. Therefore, we used 30 of the known 38 species to reconstruct the phylogeny based on transcriptome data of 3954 single-copy nuclear genes (SCGs) via coalescent and concatenated approaches and three comparative datasets (nt/nt12/aa). Based on all these methods, Ceratozamia is divided into six phylogenetic subclades within three major clades. There were a few discrepancies regarding phylogenetic position of some species within these subclades. Using these phylogenetic trees, biogeographic history and morphological diversity of the genus are explored. Ceratozamia originated from ancestors in southern Mexico since the mid-Miocene. There is a distinct distribution pattern of species through the Trans-Mexican Volcanic Belt (TMVB), that act as a barrier for the species dispersal at TMVB and its southern and northern part. Limited dispersal events occurred during the late Miocene, and maximum diversification happened during the Pliocene epoch. Our study provides a new insight into phylogenetic relationships, the origin and dispersal routes, and morphological diversity of the genus Ceratozamia. We also explain how past climatic changes affected the diversification of this Mesoamerica-native genus.

Identification of potential light deficiency response regulators in endangered species Magnolia sinostellata.

Magnolia sinostellata is one of the endangered species in China and largely suffers light deficiency stress in the understory of forest. However, the weak light response molecular mechanism remains unclear. More importantly, hub genes in the molecular network have not been pinpointed. To explore potential regulators in the mechanism, weighted gene co-expression network analysis (WGCNA) was performed to analysis the trancriptome data of M. sinostellata leaves subjected to weak light with different time points. Gene co-expression analysis illustrated that module 1, 2 and 3 were closely associated with light deficiency treatment, which. Gene ontology and KEGG analyses showed that genes in module 1 mainly participated in amino and nucleotide metabolism, module 2 mostly involved in carbon fixation and module 3 mostly regulated photosynthesis related pathways, among which 6, 7 and 8 hub genes were identified, respectively. Hub genes isoform_107196 in module 1 and isoform_55976 in module 2 were unique to M. sinostellata. This study found that light deficiency inhibited photosynthesis and stress tolerance, while improved carbon metabolism and flowering related pathways in M. sinostellata, which can impact its accumulation reserves of growth and reproduction in the next season. In addition, key shade response regulators identified in this study have laid a firm foundation for further investigation of shade response molecular mechanism and protection of other shade sensitive plants.

Phylotranscriptomics reveal the spatio-temporal distribution and morphological evolution of Macrozamia, an Australian endemic genus of Cycadales.

BACKGROUND AND AIMS: Cycads are regarded as an ancient lineage of living seed plants, and hold important clues to understand the early evolutionary trends of seed plants. The molecular phylogeny and spatio-temporal diversification of one of the species-rich genera of cycads, Macrozamia, have not been well reconstructed. METHODS: We analysed a transcriptome dataset of 4740 single-copy nuclear genes (SCGs) of 39 Macrozamia species and two outgroup taxa. Based on concatenated (maximum parsimony, maximum likelihood) and multispecies coalescent analyses, we first establish a well-resolved phylogenetic tree of Macrozamia. To identify cyto-nuclear incongruence, the plastid protein coding genes (PCGs) from transcriptome data are extracted using the software HybPiper. Furthermore, we explore the biogeographical history of the genus and shed light on the pattern of floristic exchange between three distinct areas of Australia. Six key diagnostic characters are traced on the phylogenetic framework using two comparative methods, and infra-generic classification is investigated. KEY RESULTS: The tree topologies of concatenated and multi-species coalescent analyses of SCGs are mostly congruent with a few conflicting nodes, while those from plastid PCGs show poorly supported relationships. The genus contains three major clades that correspond to their distinct distributional areas in Australia. The crown group of Macrozamia is estimated to around 11.80 Ma, with a major expansion in the last 5-6 Myr. Six morphological characters show homoplasy, and the traditional phenetic sectional division of the genus is inconsistent with this current phylogeny. CONCLUSIONS: This first detailed phylogenetic investigation of Macrozamia demonstrates promising prospects of SCGs in resolving phylogenetic relationships within cycads. Our study suggests that Macrozamia, once widely distributed in Australia, underwent major extinctions because of fluctuating climatic conditions such as cooling and mesic biome disappearance in the past. The current close placement of morphologically distinct species in the phylogenetic tree may be related to neotenic events that occurred in the genus.

Recent advances on phylogenomics of gymnosperms and a new classification.

Living gymnosperms comprise four major groups: cycads, Ginkgo, conifers, and gnetophytes. Relationships among/within these lineages have not been fully resolved. Next generation sequencing has made available a large number of sequences, including both plastomes and single-copy nuclear genes, for reconstruction of solid phylogenetic trees. Recent advances in gymnosperm phylogenomic studies have updated our knowledge of gymnosperm systematics. Here, we review major advances of gymnosperm phylogeny over the past 10 years and propose an updated classification of extant gymnosperms. This new classification includes three classes (Cycadopsida, Ginkgoopsida, and Pinopsida), five subclasses (Cycadidae, Ginkgoidae, Cupressidae, Pinidae, and Gnetidae), eight orders (Cycadales, Ginkgoales, Araucariales, Cupressales, Pinales, Ephedrales, Gnetales, and Welwitschiales), 13 families, and 86 genera. We also described six new tribes including Acmopyleae Y. Yang, Austrocedreae Y. Yang, Chamaecyparideae Y. Yang, Microcachrydeae Y. Yang, Papuacedreae Y. Yang, and Prumnopityeae Y. Yang, and made 27 new combinations in the genus Sabina.

Major episodes of horizontal gene transfer drove the evolution of land plants.

How horizontal gene transfer (HGT) has contributed to the evolution of animals and plants remains a major puzzle. Despite recent progress, defining the overall scale and pattern of HGT events in land plants has been largely elusive. In this study, we performed systematic analyses for acquired genes in different plant groups and throughout land plant evolution. We found that relatively recent HGT events occurred in charophytes and all major land plant groups, but their frequency declined rapidly in seed plants. Two major episodes of HGT events occurred in land plant evolution, corresponding to the early evolution of streptophytes and the origin of land plants, respectively. Importantly, a vast majority of the genes acquired in the two episodes have been retained in descendant groups, affecting numerous activities and processes of land plants. We analyzed some of the acquired genes involved in stress responses, ion and metabolite transport, growth and development, and specialized metabolism, and further assessed the cumulative effects of HGT in land plants.

Genomes shed light on the evolution of Begonia, a mega-diverse genus.

Clarifying the evolutionary processes underlying species diversification and adaptation is a key focus of evolutionary biology. Begonia (Begoniaceae) is one of the most species-rich angiosperm genera with c. 2000 species, most of which are shade-adapted. Here, we present chromosome-scale genome assemblies for four species of Begonia (B. loranthoides, B. masoniana, B. darthvaderiana and B. peltatifolia), and whole genome shotgun data for an additional 74 Begonia representatives to investigate lineage evolution and shade adaptation of the genus. The four genome assemblies range in size from 331.75 Mb (B. peltatifolia) to 799.83 Mb (B. masoniana), and harbor 22 059-23 444 protein-coding genes. Synteny analysis revealed a lineage-specific whole-genome duplication (WGD) that occurred just before the diversification of Begonia. Functional enrichment of gene families retained after WGD highlights the significance of modified carbohydrate metabolism and photosynthesis possibly linked to shade adaptation in the genus, which is further supported by expansions of gene families involved in light perception and harvesting. Phylogenomic reconstructions and genomics studies indicate that genomic introgression has also played a role in the evolution of Begonia. Overall, this study provides valuable genomic resources for Begonia and suggests potential drivers underlying the diversity and adaptive evolution of this mega-diverse clade.

Chloranthus genome provides insights into the early diversification of angiosperms.

Chloranthales remain the last major mesangiosperm lineage without a nuclear genome assembly. We therefore assemble a high-quality chromosome-level genome of Chloranthus spicatus to resolve enigmatic evolutionary relationships, as well as explore patterns of genome evolution among the major lineages of mesangiosperms (eudicots, monocots, magnoliids, Chloranthales, and Ceratophyllales). We find that synteny is highly conserved between genomic regions of Amborella, Vitis, and Chloranthus. We identify an ancient single whole-genome duplication (WGD) (κ) prior to the divergence of extant Chloranthales. Phylogenetic inference shows Chloranthales as sister to magnoliids. Furthermore, our analyses indicate that ancient hybridization may account for the incongruent phylogenetic placement of Chloranthales + magnoliids relative to monocots and eudicots in nuclear and chloroplast trees. Long genes and long introns are found to be prevalent in both Chloranthales and magnoliids compared to other angiosperms. Overall, our findings provide an improved context for understanding mesangiosperm relationships and evolution and contribute a valuable genomic resource for future investigations.

The complete mitochondrial genome of Cycas debaoensis revealed unexpected static evolution in gymnosperm species.

Mitochondrial genomes of vascular plants are well known for their liability in architecture evolution. However, the evolutionary features of mitogenomes at intra-generic level are seldom studied in vascular plants, especially among gymnosperms. Here we present the complete mitogenome of Cycas debaoensis, an endemic cycad species to the Guangxi region in southern China. In addition to assemblage of draft mitochondrial genome, we test the conservation of gene content and mitogenomic stability by comparing it to the previously published mitogenome of Cycas taitungensis. Furthermore, we explored the factors such as structural rearrangements and nuclear surveillance of double-strand break repair (DSBR) proteins in Cycas in comparison to other vascular plant groups. The C. debaoensis mitogenome is 413,715 bp in size and encodes 69 unique genes, including 40 protein coding genes, 26 tRNAs, and 3 rRNA genes, similar to that of C. taitungensis. Cycas mitogenomes maintained the ancestral intron content of seed plants (26 introns), which is reduced in other lineages of gymnosperms, such as Ginkgo biloba, Taxus cuspidata and Welwitschia mirabilis due to selective pressure or retroprocessing events. C. debaoensis mitogenome holds 1,569 repeated sequences (> 50 bp), which partially account for fairly large intron size (1200 bp in average) of Cycas mitogenome. The comparison of RNA-editing sites revealed 267 shared non-silent editing site among predicted vs. empirically observed editing events. Another 33 silent editing sites from empirical data increase the total number of editing sites in Cycas debaoensis mitochondrial protein coding genes to 300. Our study revealed unexpected conserved evolution between the two Cycas species. Furthermore, we found strict collinearity of the gene order along with the identical set of genomic content in Cycas mt genomes. The stability of Cycas mt genomes is surprising despite the existence of large number of repeats. This structural stability may be related to the relative expansion of three DSBR protein families (i.e., RecA, OSB, and RecG) in Cycas nuclear genome, which inhibit the homologous recombinations, by monitoring the accuracy of mitochondrial chromosome repair.

Plastid genomes and phylogenomics of liverworts (Marchantiophyta): Conserved genome structure but highest relative plastid substitution rate in land plants.

With some 7300 species of small nonvascular spore-producing plants, liverworts represent one of the major lineages of land plants. Although multi-locus molecular phylogenetic studies have elucidated relationships of liverworts at different taxonomic categories, the backbone phylogeny of liverworts is still to be fully resolved, especially for the placement of Ptilidiales and the relationships within Jungermanniales and Marchantiales. Here, we provided phylogenomic inferences of liverworts based on 42 newly sequenced and 24 published liverwort plastid genomes representing all but two orders of liverworts, and characterized the evolution of the plastome in liverworts. The structure of the plastid genome is overall conserved across the phylogeny of liverworts, with only two structural variants detected from simple thalloids, besides 18 out of 43 liverwort genera showing intron variations in their plastomes. Complex thalloid liverworts maintain the most plastid genes, and seem to undergo fewer gene deletions and pseudogenization events than other liverworts. Plastid phylogenetic inferences yielded mostly robustly supported relationships, and consistently resolved Ptilidiales as the sister to Porellales. The relative ratio of silent substitutions across the three genetic compartments (i.e., 1:15:10, for mitochondrial:plastid:nuclear) suggests that liverwort plastid genes have the potential to evolve faster than their nuclear counterparts, unlike in any other major land plant lineages where the mutation rate of nuclear genes overwhelm those of their plastid and mitochondrial counterparts.

The genome of Magnolia biondii Pamp. provides insights into the evolution of Magnoliales and biosynthesis of terpenoids.

Magnolia biondii Pamp. (Magnoliaceae, magnoliids) is a phylogenetically, economically, and medicinally important ornamental tree species widely grown and cultivated in the north-temperate regions of China. Determining the genome sequence of M. biondii would help resolve the phylogenetic uncertainty of magnoliids and improve the understanding of individual trait evolution within the Magnolia genus. We assembled a chromosome-level reference genome of M. biondii using ~67, ~175, and ~154 Gb of raw DNA sequences generated via Pacific Biosciences single-molecule real-time sequencing, 10X Genomics Chromium, and Hi-C scaffolding strategies, respectively. The final genome assembly was ~2.22 Gb, with a contig N50 value of 269.11 kb and a BUSCO complete gene percentage of 91.90%. Approximately 89.17% of the genome was organized into 19 chromosomes, resulting in a scaffold N50 of 92.86 Mb. The genome contained 47,547 protein-coding genes, accounting for 23.47% of the genome length, whereas 66.48% of the genome length consisted of repetitive elements. We confirmed a WGD event that occurred very close to the time of the split between the Magnoliales and Laurales. Functional enrichment of the Magnolia-specific and expanded gene families highlighted genes involved in the biosynthesis of secondary metabolites, plant-pathogen interactions, and responses to stimuli, which may improve the ecological fitness and biological adaptability of the lineage. Phylogenomic analyses revealed a sister relationship of magnoliids and Chloranthaceae, which are sister to a clade comprising monocots and eudicots. The genome sequence of M. biondii could lead to trait improvement, germplasm conservation, and evolutionary studies on the rapid radiation of early angiosperms.

Dipterisshenzhenensis, a new endangered species of Dipteridaceae from Shenzhen, southern China.

Dipterisshenzhenensis, a new species of ferns from Shenzhen, Guangdong, southern China, is identified and described. It closely resembles D.chinensis but possesses several unique traits, such as long rhizome scales, castaneous stipe, and abaxially pale fronds with two fan-shaped fronds connected by a broad wing. Molecular evidence showed that D.shenzhenensis is allied to D.conjugata, whereas it has morphologically significant differences (P < 0.05) on the basis of quantitative trait statistical analysis. Overall, the morphological evidence, taken together with the result of cpDNA indicated that D.shenzhenensis is a distinct species.

Gnetum chinense, a new species of Gnetaceae from southwestern China.

Gnetum chinense sp. nov., a new lianoid species of Gnetaceae, is described from southwestern China. The new species is morphologically similar to G. montanum Markgr. in its oblong elliptic leaves and the ovoid to ellipsoid chlamydosperm, but differs from the latter by its shorter male spikes having fewer involucral collars (7-10 vs. 13-18 in G. montanum). We also did a new molecular analysis using one nuclear marker (i.e. nrITS) and four chloroplast markers (i.e. matK gene, rpoC1 intron, psbB-rps12 IGS, and trnF-trnV IGS). The result suggests that this specific clade is sister to a large clade consisting of all other known Chinese lianoid species of Gnetum except G. parvifolium (Warb.) W.C. Cheng.

Severe Plastid Genome Size Reduction in a Mycoheterotrophic Orchid, Danxiaorchis singchiana, Reveals Heavy Gene Loss and Gene Relocations.

Danxiaorchis singchiana (Orchidaceae) is a leafless mycoheterotrophic orchid in the subfamily Epidendroideae. We sequenced the complete plastome of D. singchiana. The plastome has a reduced size of 87,931 bp, which includes a pair of inverted repeat (IR) regions of 13,762 bp each that are separated by a large single copy (LSC) region of 42,575 bp and a small single copy (SSC) region of 17,831 bp. When compared to its sister taxa, Cremastra appendiculata and Corallorhiza striata var. involuta, D. singchiana showed an inverted gene block in the LSC and SSC regions. A total of 61 genes were predicted, including 21 tRNA, 4 rRNA, and 36 protein-coding genes. While most of the housekeeping genes were still intact and seem to be protein-coding, only four photosynthesis-related genes appeared presumably intact. The majority of the presumably intact protein-coding genes seem to have undergone purifying selection (dN/dS < 1), and only the psaC gene was positively selected (dN/dS > 1) when compared to that in Cr. appendiculata. Phylogenetic analysis of 26 complete plastome sequences from 24 species of the tribe Epidendreae had revealed that D. singchiana diverged after Cr. appendiculata and is sister to the genus Corallorhiza with strong bootstrap support (100%).

The draft mitochondrial genome of Magnolia biondii and mitochondrial phylogenomics of angiosperms.

The mitochondrial genomes of flowering plants are well known for their large size, variable coding-gene set and fluid genome structure. The available mitochondrial genomes of the early angiosperms show extreme genetic diversity in genome size, structure, and sequences, such as rampant HGTs in Amborella mt genome, numerous repeated sequences in Nymphaea mt genome, and conserved gene evolution in Liriodendron mt genome. However, currently available early angiosperm mt genomes are still limited, hampering us from obtaining an overall picture of the mitogenomic evolution in angiosperms. Here we sequenced and assembled the draft mitochondrial genome of Magnolia biondii Pamp. from Magnoliaceae (magnoliids) using Oxford Nanopore sequencing technology. We recovered a single linear mitochondrial contig of 967,100 bp with an average read coverage of 122 × and a GC content of 46.6%. This draft mitochondrial genome contains a rich 64-gene set, similar to those of Liriodendron and Nymphaea, including 41 protein-coding genes, 20 tRNAs, and 3 rRNAs. Twenty cis-spliced and five trans-spliced introns break ten protein-coding genes in the Magnolia mt genome. Repeated sequences account for 27% of the draft genome, with 17 out of the 1,145 repeats showing recombination evidence. Although partially assembled, the approximately 1-Mb mt genome of Magnolia is still among the largest in angiosperms, which is possibly due to the expansion of repeated sequences, retention of ancestral mtDNAs, and the incorporation of nuclear genome sequences. Mitochondrial phylogenomic analysis of the concatenated datasets of 38 conserved protein-coding genes from 91 representatives of angiosperm species supports the sister relationship of magnoliids with monocots and eudicots, which is congruent with plastid evidence.

The hornwort genome and early land plant evolution.

Hornworts, liverworts and mosses are three early diverging clades of land plants, and together comprise the bryophytes. Here, we report the draft genome sequence of the hornwort Anthoceros angustus. Phylogenomic inferences confirm the monophyly of bryophytes, with hornworts sister to liverworts and mosses. The simple morphology of hornworts correlates with low genetic redundancy in plant body plan, while the basic transcriptional regulation toolkit for plant development has already been established in this early land plant lineage. Although the Anthoceros genome is small and characterized by minimal redundancy, expansions are observed in gene families related to RNA editing, UV protection and desiccation tolerance. The genome of A. angustus bears the signatures of horizontally transferred genes from bacteria and fungi, in particular of genes operating in stress-response and metabolic pathways. Our study provides insight into the unique features of hornworts and their molecular adaptations to live on land.

The water lily genome and the early evolution of flowering plants.

Water lilies belong to the angiosperm order Nymphaeales. Amborellales, Nymphaeales and Austrobaileyales together form the so-called ANA-grade of angiosperms, which are extant representatives of lineages that diverged the earliest from the lineage leading to the extant mesangiosperms1-3. Here we report the 409-megabase genome sequence of the blue-petal water lily (Nymphaea colorata). Our phylogenomic analyses support Amborellales and Nymphaeales as successive sister lineages to all other extant angiosperms. The N. colorata genome and 19 other water lily transcriptomes reveal a Nymphaealean whole-genome duplication event, which is shared by Nymphaeaceae and possibly Cabombaceae. Among the genes retained from this whole-genome duplication are homologues of genes that regulate flowering transition and flower development. The broad expression of homologues of floral ABCE genes in N. colorata might support a similarly broadly active ancestral ABCE model of floral organ determination in early angiosperms. Water lilies have evolved attractive floral scents and colours, which are features shared with mesangiosperms, and we identified their putative biosynthetic genes in N. colorata. The chemical compounds and biosynthetic genes behind floral scents suggest that they have evolved in parallel to those in mesangiosperms. Because of its unique phylogenetic position, the N. colorata genome sheds light on the early evolution of angiosperms.