The water lily genome and the early evolution of flowering plants

Underwater CAM photosynthesis elucidated by Isoetes genome

To conserve water in arid environments, numerous plant lineages have independently evolved Crassulacean Acid Metabolism (CAM). Interestingly, Isoetes, an aquatic lycophyte, can also perform CAM as an adaptation to low CO2 availability underwater. However, little is known about the evolution of CAM in aquatic plants and the lack of genomic data has hindered comparison between aquatic and terrestrial CAM. Here, we investigate underwater CAM in Isoetes taiwanensis by generating a high-quality genome assembly and RNA-seq time course. Despite broad similarities between CAM in Isoetes and terrestrial angiosperms, we identify several key differences. Notably, Isoetes may have recruited the lesser-known 'bacterial-type' PEPC, along with the 'plant-type' exclusively used in other CAM and C4 plants for carboxylation of PEP. Furthermore, we find that circadian control of key CAM pathway genes has diverged considerably in Isoetes relative to flowering plants. This suggests the existence of more evolutionary paths to CAM than previously recognized.

Metabolite analysis in Nymphaea 'Blue Bird' petals reveal the roles of flavonoids in color formation, stress amelioration, and bee orientation

In this study, metabolome of open petals (OP) and closed petals (CP) from Nymphaea 'Blue Bird' was firstly investigated. A total of 455 metabolites was identified in Nymphaea 'Blue Bird' petals, which was mainly composed of 100 flavonoids, 83 phenolic acids, 64 amino acids and derivatives, 60 lipids, 32 alkaloids, 32 organic acids, 24 nucleotides and derivatives, and 12 lignans and coumarins. By differential analysis, 192 metabolites were identified with variable importance in project ≥ 1, among which 83 and 109 metabolites were up- and down-regulated in OP group, respectively. Further analysis (Log₂ fold change ≥ 1) identified 26 and 7 metabolites exhibited significantly lower and higher contents in CP group, relative to OP group. Importantly, KEGG analysis indicated that flavonoid biosynthesis exhibited the most significant enrichment. qRT-PCR analysis indicated that the PAL, CHS, and HCDBR genes showed a significantly higher expression in OP group than in CP group. These data explain the increase of naringenin chalcone and phloretin in OP. However, there was no significant difference of total flavonoids between OP and CP groups. Considering the increase of H₂O₂ content and ultraviolet (UV) absorption peak in OP, our results implied that diurnal stressful conditions induced the degradation of flavonoids, which contributed to environmental stress amelioration. Moreover, a higher absorption peak of 360-380 nm UV was observed in the extract liquor of OP. The sensitivity maximum of the UV-photoreceptor of bees is situated around 340-380 nm UV. This suggested, as noted for the maximum absorption of dihydrokaempferol in 340-370 nm, rhythmic accumulation and loss of these differential flavonoids in Nymphaea 'Blue Bird' petals might enhance UV pattern to some degree, influencing pollinator attraction.

Biosynthesis and emission of methyl hexanoate, the major constituent of floral scent of a night-blooming water lily Victoriacruziana

Among the factors that have made flowering plants the most species-rich lineage of land plants is the interaction between flower and insect pollinators, for which floral scent plays a pivotal role. Water lilies belong to the ANA (Amborellales, Nymphaeales, and Austrobaileyales) grade of basal flowering plants. In this study, Victoria cruziana was investigated as a model night-blooming water lily for floral scent biosynthesis. Four volatile compounds, including three benzenoids and one fatty acid methyl ester methyl hexanoate, were detected from the flowers of V. cruziana during their first bloom, with methyl hexanoate accounting for 45 % of total floral volatile emission. Emission rates were largely constant before significant drop starting at the end of second bloom. To understand the molecular basis of floral scent biosynthesis in V. cruziana, particularly methyl hexanoate, a transcriptome from the whole flowers at the full-bloom stage was created and analyzed. Methyl hexanoate was hypothesized to be biosynthesized by SABATH methyltransferases. From the transcriptome, three full-length SABATH genes designated VcSABATH1-3 were identified. A full-length cDNA for each of the three VcSABATH genes was expressed in Escherichia coli to produce recombinant proteins. When tested in in vitro methyltransferase enzyme assays with different fatty acids, both VcSABATH1 and VcSABATH3 exhibited highest levels of activity with hexanoic acid to produce methyl hexanoate, with the specific activity of VcSABATH1 being about 15 % of that for VcSABATH3. VcSABATH1 and VcSABATH3 showed the highest levels of expression in stamen and pistil, respectively. In phylogenetic analysis, three VcSABATH genes clustered with other water lily SABATH methyltransferase genes including the one known for making other fatty acid methyl esters, implying both a common evolutionary origin and functional divergence. Fatty acid methyl esters are not frequent constituents of floral scents of mesangiosperms, pointing to the importance for the evolution of novel fatty acid methyltransferase for making fatty acid methyl esters in the pollination biology of water lilies.

Plastid phylogenomic insights into relationships of all flowering plant families

BACKGROUND

Flowering plants (angiosperms) are dominant components of global terrestrial ecosystems, but phylogenetic relationships at the familial level and above remain only partially resolved, greatly impeding our full understanding of their evolution and early diversification. The plastome, typically mapped as a circular genome, has been the most important molecular data source for plant phylogeny reconstruction for decades.

RESULTS

Here, we assembled by far the largest plastid dataset of angiosperms, composed of 80 genes from 4792 plastomes of 4660 species in 2024 genera representing all currently recognized families. Our phylogenetic tree (PPA II) is essentially congruent with those of previous plastid phylogenomic analyses but generally provides greater clade support. In the PPA II tree, 75% of nodes at or above the ordinal level and 78% at or above the familial level were resolved with high bootstrap support (BP ≥ 90). We obtained strong support for many interordinal and interfamilial relationships that were poorly resolved previously within the core eudicots, such as Dilleniales, Saxifragales, and Vitales being resolved as successive sisters to the remaining rosids, and Santalales, Berberidopsidales, and Caryophyllales as successive sisters to the asterids. However, the placement of magnoliids, although resolved as sister to all other Mesangiospermae, is not well supported and disagrees with topologies inferred from nuclear data. Relationships among the five major clades of Mesangiospermae remain intractable despite increased sampling, probably due to an ancient rapid radiation.

CONCLUSIONS

We provide the most comprehensive dataset of plastomes to date and a well-resolved phylogenetic tree, which together provide a strong foundation for future evolutionary studies of flowering plants.

The Melastoma dodecandrum genome and the evolution of Myrtales

Melastomataceae have abundant morphological diversity with high economic and ornamental merit in Myrtales. The phylogenetic position of Myrtales is still contested. Here, we report the first chromosome-level genome assembly of Melastoma dodecandrum in Melastomataceae. The assembled genome size was 299.81 Mb with a contig N50 value of 3.00 Mb. Genome evolution analysis indicated that M. dodecandrum, Eucalyptus grandis and Punica granatum were clustered into a clade of Myrtales and formed a sister group with the ancestor of fabids and malvids. We found that M. dodecandrum experienced four whole-genome polyploidization events: the ancient event was shared with most eudicots, one event was shared with Myrtales, and the other two events were unique to M. dodecandrum. Moreover, we identified MADS-box genes and found that the AP1-like genes expanded, and AP3-like genes might have undergone subfunctionalization. We found that the SUAR63-like genes and AG-like genes showed different expression patterns in stamens, which may be associated with heteranthery. In addition, we found that LAZY1-like genes were involved in the negative regulation of stem branching development, which may be related to its creeping features. Our study sheds new light on the evolution of Melastomataceae and Myrtales, which provides a comprehensive genetic resource for future research.

Phylogeny of gracillariid leaf-mining moths: evolution of larval behaviour inferred from phylogenomic and Sanger data

Gracillariidae is the most taxonomically diverse cosmopolitan leaf-mining moth family, consisting of nearly 2000 named species in 105 described genera, classified into eight extant subfamilies. The majority of gracillariid species are internal plant feeders as larvae, creating mines and galls in plant tissue. Despite their diversity and ecological adaptations, their phylogenetic relationships, especially among subfamilies, remain uncertain. Genomic data (83 taxa, 589 loci) were integrated with Sanger data (130 taxa, 22 loci), to reconstruct a phylogeny of Gracillariidae. Based on analyses of both datasets combined and analyzed separately, monophyly of Gracillariidae and all its subfamilies, monophyly of the clade "LAMPO" (subfamilies: Lithocolletinae, Acrocercopinae, Marmarinae, Phyllocnistinae, and Oecophyllembiinae) and relationships of its subclade "AMO" (subfamilies: Acrocercopinae, Marmarinae, and Oecophyllembiinae) were strongly supported. A sister-group relationship of Ornixolinae to the remainder of the family, and a monophyletic leaf roller lineage (Callicercops Vári + Parornichinae) + Gracillariinae, as sister to the "LAMPO" clade were supported by the most likely tree. Dating analyses indicate a mid-Cretaceous (105.3 Ma) origin of the family, followed by a rapid diversification into the nine subfamilies predating the Cretaceous-Palaeogene extinction. We hypothesize that advanced larval behaviours, such as making keeled or tentiform blotch mines, rolling leaves and galling, allowed gracillariids to better avoid larval parasitoids allowing them to further diversify. Finally, we stabilize the classification by formally re-establishing the subfamily ranks of Marmarinae stat.rev., Oecophyllembiinae stat.rev. and Parornichinae stat.rev., and erect a new subfamily, Callicercopinae Li, Ohshima and Kawahara to accommodate the enigmatic genus Callicercops.

Physalis floridana CRABS CLAW mediates neofunctionalization of GLOBOSA genes in carpel development

Floral B-function MADS-box genes, such as GLOBOSA (GLO), function in corolla and stamen organ identity specification. The functions of these genes outside these floral whorls are rarely reported. DOLL1 is a GLO gene controlling corolla and androecium organ identity. In this study we found that, in Physalis floridana double-layered-lantern 1 (doll1) mutant pollinated with wild-type pollen, fruit set was extremely low, indicating that doll1 females are dysfunctional. Stigma and style structure, stigma receptivity, pollen tube guidance, and embryo sac development were also impaired in doll1. P. floridana CRABS CLAW (PFCRC), predominantly expressed in carpels, was repressed in doll1 native carpels. Loss-of-function of PFCRC altered carpel meristem determinacy, carpel closure, and ovule number, and the resultant 'pistil' consisted of multiple spirally-arranged dorsiventral carpels occasionally with 1-2 naked ovules on the margin and trichomes at each mutated carpel tip, implying an alteration of carpel organ identity. Regulatory and genetic interactions between B-class MADS-box genes and PFCRC were revealed in a context-dependent manner in floral development. Our work reveals a new role for the B-function genes in carpel and ovule development via regulating PFCRC, providing a new understanding of genetic regulatory networks between MADS-domain and CRC transcription factors in mediating carpel organ specification, functionality, and origin.

Taxus yunnanensis genome offers insights into gymnosperm phylogeny and taxol production

Taxol, a natural product derived from Taxus, is one of the most effective natural anticancer drugs and the biosynthetic pathway of Taxol is the basis of heterologous bio-production. Here, we report a high-quality genome assembly and annotation of Taxus yunnanensis based on 10.7 Gb sequences assembled into 12 chromosomes with contig N50 and scaffold N50 of 2.89 Mb and 966.80 Mb, respectively. Phylogenomic analyses show that T. yunnanensis is most closely related to Sequoiadendron giganteum among the sampled taxa, with an estimated divergence time of 133.4-213.0 MYA. As with most gymnosperms, and unlike most angiosperms, there is no evidence of a recent whole-genome duplication in T. yunnanensis. Repetitive sequences, especially long terminal repeat retrotransposons, are prevalent in the T. yunnanensis genome, contributing to its large genome size. We further integrated genomic and transcriptomic data to unveil clusters of genes involved in Taxol synthesis, located on the chromosome 12, while gene families encoding hydroxylase in the Taxol pathway exhibited significant expansion. Our study contributes to the further elucidation of gymnosperm relationships and the Taxol biosynthetic pathway.

Transcriptome Analysis Reveals the Senescence Process Controlling the Flower Opening and Closure Rhythm in the Waterlilies (Nymphaea L.)

Most waterlily flowers open at dawn and close after noon usually for three to four days, and thereafter wilt. The short lifespan of flowers restricts the development of the flower postharvest industry. The termination of flower movements is a key event during flower aging process. However, it is still unclear when the senescence process initiates and how it terminates the movement rhythm. In this study, we observed that the opening diameter of flowers was the smallest on the fourth (last) flowering day. Subsequent transcriptome profiles generated from petals at different flowering stages showed that the multiple signaling pathways were activated at the last closure stage (Time 3, T3) of the flowers, including Ca²⁺, reactive oxygen species and far red light signaling pathways, as well as auxin, ethylene and jasmonic acid signaling pathways. Moreover, In terms of cell metabolism regulation, the genes related to hydrolase (protease, phospholipase, nuclease) were upregulated at T3 stage, indicating that petals entered the senescence stage at that time; and the genes related to water transport and cell wall modification were also differentially regulated at T3 stage, which would affect the ability of cell expand and contract, and eventually lead to petal not open after the fourth day. Collectively, our data provided a new insight into the termination of flower opening in the waterlilies, and a global understanding of the senescence process of those opening-closure rhythm flowers.

Expression and Functional Analyses of Nymphaea caerulea MADS-Box Genes Contribute to Clarify the Complex Flower Patterning of Water Lilies

Nymphaeaceae are early diverging angiosperms with large flowers characterized by showy petals and stamens not clearly whorled but presenting a gradual morphological transition from the outer elements to the inner stamens. Such flower structure makes these plant species relevant for studying flower evolution. MADS-domain transcription factors are crucial components of the molecular network that controls flower development. We therefore isolated and characterized MADS-box genes from the water lily Nymphaea caerulea. RNA-seq experiments on floral buds have been performed to obtain the transcript sequences of floral organ identity MADS-box genes. Maximum Likelihood phylogenetic analyses confirmed their belonging to specific MADS-box gene subfamilies. Their expression was quantified by RT-qPCR in all floral organs at two stages of development. Protein interactions among these transcription factors were investigated by yeast-two-hybrid assays. We found especially interesting the involvement of two different AGAMOUS-like genes (NycAG1 and NycAG2) in the water lily floral components. They were therefore functionally characterized by complementing Arabidopsis ag and shp1 shp2 mutants. The expression analysis of MADS-box genes across flower development in N. caerulea described a complex scenario made of numerous genes in numerous floral components. Their expression profiles in some cases were in line with what was expected from the ABC model of flower development and its extensions, while in other cases presented new and interesting gene expression patterns, as for instance the involvement of NycAGL6 and NycFL. Although sharing a high level of sequence similarity, the two AGAMOUS-like genes NycAG1 and NycAG2 could have undergone subfunctionalization or neofunctionalization, as only one of them could partially restore the euAG function in Arabidopsis ag-3 mutants. The hereby illustrated N. caerulea MADS-box gene expression pattern might mirror the morphological transition from the outer to the inner floral organs, and the presence of transition organs such as the petaloid stamens. This study is intended to broaden knowledge on the role and evolution of floral organ identity genes and the genetic mechanisms causing biodiversity in angiosperm flowers.

Organelle Genomes and Transcriptomes of Nymphaea Reveal the Interplay between Intron Splicing and RNA Editing

Posttranscriptional modifications, including intron splicing and RNA editing, are common processes during regulation of gene expression in plant organelle genomes. However, the intermediate products of intron-splicing, and the interplay between intron-splicing and RNA-editing were not well studied. Most organelle transcriptome analyses were based on the Illumina short reads which were unable to capture the full spectrum of transcript intermediates within an organelle. To fully investigate the intermediates during intron splicing and the underlying relationships with RNA editing, we used PacBio DNA-seq and Iso-seq, together with Illumina short reads genome and transcriptome sequencing data to assemble the chloroplast and mitochondrial genomes of Nymphaea 'Joey Tomocik' and analyze their posttranscriptional features. With the direct evidence from Iso-seq, multiple intermediates partially or fully intron-spliced were observed, and we also found that both cis- and trans-splicing introns were spliced randomly. Moreover, by using rRNA-depleted and non-Oligo(dT)-enrichment strand-specific RNA-seq data and combining direct SNP-calling and transcript-mapping methods, we identified 98 and 865 RNA-editing sites in the plastome and mitogenome of N. 'Joey Tomocik', respectively. The target codon preference, the tendency of increasing protein hydrophobicity, and the bias distribution of editing sites are similar in both organelles, suggesting their common evolutionary origin and shared editing machinery. The distribution of RNA editing sites also implies that the RNA editing sites in the intron and exon regions may splice synchronously, except those exonic sites adjacent to intron which could only be edited after being intron-spliced. Our study provides solid evidence for the multiple intermediates co-existing during intron-splicing and their interplay with RNA editing in organelle genomes of a basal angiosperm.

Insights into angiosperm evolution, floral development and chemical biosynthesis from the Aristolochia fimbriata genome

Aristolochia, a genus in the magnoliid order Piperales, has been famous for centuries for its highly specialized flowers and wide medicinal applications. Here, we present a new, high-quality genome sequence of Aristolochia fimbriata, a species that, similar to Amborella trichopoda, lacks further whole-genome duplications since the origin of extant angiosperms. As such, the A. fimbriata genome is an excellent reference for inferences of angiosperm genome evolution, enabling detection of two novel whole-genome duplications in Piperales and dating of previously reported whole-genome duplications in other magnoliids. Genomic comparisons between A. fimbriata and other angiosperms facilitated the identification of ancient genomic rearrangements suggesting the placement of magnoliids as sister to monocots, whereas phylogenetic inferences based on sequence data we compiled yielded ambiguous relationships. By identifying associated homologues and investigating their evolutionary histories and expression patterns, we revealed highly conserved floral developmental genes and their distinct downstream regulatory network that may contribute to the complex flower morphology in A. fimbriata. Finally, we elucidated the genetic basis underlying the biosynthesis of terpenoids and aristolochic acids in A. fimbriata.

Molecular Cloning and Expression Profile of Class E Genes Related to Sepal Development in Nelumbo nucifera

The lotus (Nelumbo Adans.) is an important aquatic plant with ornamental, medicinal and edible values and cultural connotations. It has single-, semi-double-, double- and thousand-petalled types of flower shape and is an ideal material for developmental research of flower doubling. The lotus is a basal eudicot species without a morphological difference between the sepals and petals and occupies a critical phylogenetic position in flowering plants. In order to investigate the genetic relationship between the sepals and petals in the lotus, the class E genes which affect sepal formation were focused on and analyzed. Here, SEPALLATA 1(NnSEP1) and its homologous genes AGAMOUS-LIKE MADS-BOXAGL9 (NnAGL9) and MADS-BOX TRANSCRIPTION FACTOR 6-like (NnMADS6-like) of the class E gene family were isolated from the flower buds of the Asian lotus (Nelumbo nucifera Gaertn.). The protein structure, subcellular localization and expression patterns of these three genes were investigated. All three genes were verified to locate in the nucleus and had typical MADS-box characteristics. NnSEP1 and NnMADS6-like were specifically expressed in the sepals, while NnAGL9 was highly expressed in the petals, suggesting that different developmental mechanisms exist in the formation of the sepals and petals in the lotus. The significant functional differences between NnSEP1, NnMADS6-like and NnAGL9 were also confirmed by a yeast two-hybrid assay. These results expand our knowledge on the class E gene family in sepal formation and will benefit fundamental research on the development of floral organs in Nelumbo.

Genome of the world's smallest flowering plant, Wolffia australiana, helps explain its specialized physiology and unique morphology

Watermeal, Wolffia australiana, is the smallest known flowering monocot and is rich in protein. Despite its great potential as a biotech crop, basic research on Wolffia is in its infancy. Here, we generated the reference genome of a species of watermeal, W. australiana, and identified the genome-wide features that may contribute to its atypical anatomy and physiology, including the absence of roots, adaxial stomata development, and anaerobic life as a turion. In addition, we found evidence of extensive genome rearrangements that may underpin the specialized aquatic lifestyle of watermeal. Analysis of the gene inventory of this intriguing species helps explain the distinct characteristics of W. australiana and its unique evolutionary trajectory.

Halim Park and Jin Hwa Park et al. report the nuclear genome sequence of the duckweed Wolffia australiana, the smallest known flowering plant. The genome assembly represents an improvement over a recently published genome and highlights genome rearrangements that may be linked to its specialized aquatic adaptations.

An update on the function and regulation of methylerythritol phosphate and mevalonate pathways and their evolutionary dynamics

Isoprenoids are among the largest and most chemically diverse classes of organic compounds in nature and are involved in the processes of photosynthesis, respiration, growth, development, and plant responses to stress. The basic building block units for isoprenoid synthesis-isopentenyl diphosphate and its isomer dimethylallyl diphosphate-are generated by the mevalonate (MVA) and methylerythritol phosphate (MEP) pathways. Here, we summarize recent advances on the roles of the MEP and MVA pathways in plant growth, development and stress responses, and attempt to define the underlying gene networks that orchestrate the MEP and MVA pathways in response to developmental or environmental cues. Through phylogenomic analysis, we also provide a new perspective on the evolution of the plant isoprenoid pathway. We conclude that the presence of the MVA pathway in plants may be associated with the transition from aquatic to subaerial and terrestrial environments, as lineages for its core components are absent in green algae. The emergence of the MVA pathway has acted as a key evolutionary event in plants that facilitated land colonization and subsequent embryo development, as well as adaptation to new and varied environments.

A chromosome-level genome assembly of rugged rose (Rosa rugosa) provides insights into its evolution, ecology, and floral characteristics

Rosa rugosa, commonly known as rugged rose, is a perennial ornamental shrub. It produces beautiful flowers with a mild fragrance and colorful seed pods. Unlike many other cultivated roses, R. rugosa adapts to a wide range of habitat types and harsh environmental conditions such as salinity, alkaline, shade, drought, high humidity, and frigid temperatures. Here, we produced and analyzed a high-quality genome sequence for R. rugosa to understand its ecology, floral characteristics and evolution. PacBio HiFi reads were initially used to construct the draft genome of R. rugosa, and then Hi-C sequencing was applied to assemble the contigs into 7 chromosomes. We obtained a 382.6 Mb genome encoding 39,704 protein-coding genes. The genome of R. rugosa appears to be conserved with no additional whole-genome duplication after the gamma whole-genome triplication (WGT), which occurred ~100 million years ago in the ancestor of core eudicots. Based on a comparative analysis of the high-quality genome assembly of R. rugosa and other high-quality Rosaceae genomes, we found a unique large inverted segment in the Chinese rose R. chinensis and a retroposition in strawberry caused by post-WGT events. We also found that floral development- and stress response signaling-related gene modules were retained after the WGT. Two MADS-box genes involved in floral development and the stress-related transcription factors DREB2A-INTERACTING PROTEIN 2 (DRIP2) and PEPTIDE TRANSPORTER 3 (PTR3) were found to be positively selected in evolution, which may have contributed to the unique ability of this plant to adapt to harsh environments. In summary, the high-quality genome sequence of R. rugosa provides a map for genetic studies and molecular breeding of this plant and enables comparative genomic studies of Rosa in the near future.

Analysis of the Coptis chinensis genome reveals the diversification of protoberberine-type alkaloids

Chinese goldthread (Coptis chinensis Franch.), a member of the Ranunculales, represents an important early-diverging eudicot lineage with diverse medicinal applications. Here, we present a high-quality chromosome-scale genome assembly and annotation of C. chinensis. Phylogenetic and comparative genomic analyses reveal the phylogenetic placement of this species and identify a single round of ancient whole-genome duplication (WGD) shared by the Ranunculaceae. We characterize genes involved in the biosynthesis of protoberberine-type alkaloids in C. chinensis. In particular, local genomic tandem duplications contribute to member amplification of a Ranunculales clade-specific gene family of the cytochrome P450 (CYP) 719. The functional versatility of a key CYP719 gene that encodes the (S)-canadine synthase enzyme involved in the berberine biosynthesis pathway may play critical roles in the diversification of other berberine-related alkaloids in C. chinensis. Our study provides insights into the genomic landscape of early-diverging eudicots and provides a valuable model genome for genetic and applied studies of Ranunculales.

Turgor regulation defect 1 proteins play a conserved role in pollen tube reproductive innovation of the angiosperms

Sexual reproduction in angiosperms is siphonogamous, and the interaction between pollen tube and pistil is critical for successful fertilization. Our previous study demonstrated that mutation of the Arabidopsis turgor regulation defect 1 (TOD1) gene leads to reduced male fertility, a result of retarded pollen tube growth in the pistil. TOD1 encodes a Golgi-localized alkaline ceramidase, a key enzyme for the production of sphingosine-1-phosphate (S1P), which is involved in the regulation of turgor pressure in plant cells. However, whether TOD1s play a conserved role in the innovation of siphonogamy is largely unknown. In this study, we provide evidence that OsTOD1, which is similar to AtTOD1, is also preferentially expressed in rice pollen grains and pollen tubes. OsTOD1 knockout results in reduced pollen tube growth potential in rice pistil. Both the OsTOD1 genomic sequence with its own promoter and the coding sequence under the AtTOD1 promoter can partially rescue the attod1 mutant phenotype. Furthermore, TOD1s from other angiosperm species can partially rescue the attod1 mutant phenotype, while TOD1s from gymnosperm species are not able to complement the attod1 mutant phenotype. Our data suggest that TOD1 acts conservatively in angiosperms, and this opens up an opportunity to dissect the role of sphingolipids in pollen tube growth in angiosperms.

Nuclear phylotranscriptomics and phylogenomics support numerous polyploidization events and hypotheses for the evolution of rhizobial nitrogen-fixing symbiosis in Fabaceae

Fabaceae are the third largest angiosperm family, with 765 genera and ∼19 500 species. They are important both economically and ecologically, and global Fabaceae crops are intensively studied in part for their nitrogen-fixing ability. However, resolution of the intrasubfamilial Fabaceae phylogeny and divergence times has remained elusive, precluding a reconstruction of the evolutionary history of symbiotic nitrogen fixation in Fabaceae. Here, we report a highly resolved phylogeny using >1500 nuclear genes from newly sequenced transcriptomes and genomes of 391 species, along with other datasets, for a total of 463 legumes spanning all 6 subfamilies and 333 of 765 genera. The subfamilies are maximally supported as monophyletic. The clade comprising subfamilies Cercidoideae and Detarioideae is sister to the remaining legumes, and Duparquetioideae and Dialioideae are successive sisters to the clade of Papilionoideae and Caesalpinioideae. Molecular clock estimation revealed an early radiation of subfamilies near the K/Pg boundary, marked by mass extinction, and subsequent divergence of most tribe-level clades within ∼15 million years. Phylogenomic analyses of thousands of gene families support 28 proposed putative whole-genome duplication/whole-genome triplication events across Fabaceae, including those at the ancestors of Fabaceae and five of the subfamilies, and further analyses supported the Fabaceae ancestral polyploidy. The evolution of rhizobial nitrogen-fixing nodulation in Fabaceae was probed by ancestral character reconstruction and phylogenetic analyses of related gene families and the results support the hypotheses of one or two switch(es) to rhizobial nodulation followed by multiple losses. Collectively, these results provide a foundation for further morphological and functional evolutionary analyses across Fabaceae.

Ancient and recent gene duplications as evolutionary drivers of the seed maturation regulators DELAY OF GERMINATION1 family genes

The DELAY OF GERMINATION1 (DOG1) family genes (DFGs) in Arabidopsis thaliana are involved in seed dormancy, reserve accumulation, and desiccation tolerance. Decoding the molecular evolution of DFGs is key to understanding how these seed programs evolved. This article demonstrates that DFGs have diverged in the four lineages DOG1, DOG1-LIKE4 (DOGL4), DOGL5 and DOGL6, whereas DOGL1, DOGL2 and DOGL3 arose separately within the DOG1 lineage. The systematic DFG nomenclature proposed in this article addresses the current issues of inconsistent DFG annotation and highlights DFG genomic synteny in angiosperms. DFG pseudogenes, or collapsed coding sequences, hidden in the genomes of early-diverging angiosperms are documented here. They suggest ancient birth and loss of DFGs over the course of angiosperm evolution. The proposed models suggest that the origin of DFG diversification dates back to the most recent common ancestor of living angiosperms. The presence of a single form of DFG in nonflowering plants is discussed. Phylogenetic analysis of gymnosperm, lycophyte, and liverwort DFGs and similar genes found in mosses and algae suggests that DFGs diverged from the TGACG motif-binding transcription factor genes before the divergence of the bryophyte lineage.

Developmental stochasticity and variation in floral phyllotaxis

Floral phyllotaxis is a relatively robust phenotype; trimerous and pentamerous arrangements are widely observed in monocots and core eudicots. Conversely, it also shows variability in some angiosperm clades such as 'ANA' grade (Amborellales, Nymphaeales, and Austrobaileyales), magnoliids, and Ranunculales. Regardless of the phylogenetic relationship, however, phyllotactic pattern formation appears to be a common process. What are the causes of the variability in floral phyllotaxis and how has the variation of floral phyllotaxis contributed to floral diversity? In this review, I summarize recent progress in studies on two related fields to develop answers to these questions. First, it is known that molecular and cellular stochasticity are inevitably found in biological systems, including plant development. Organisms deal with molecular stochasticity in several ways, such as dampening noise through gene networks or maintaining function through cellular redundancy. Recent studies on molecular and cellular stochasticity suggest that stochasticity is not always detrimental to plants and that it is also essential in development. Second, studies on vegetative and inflorescence phyllotaxis have shown that plants often exhibit variability and flexibility in phenotypes. Three types of phyllotaxis variations are observed, namely, fluctuation around the mean, transition between regular patterns, and a transient irregular organ arrangement called permutation. Computer models have demonstrated that stochasticity in the phyllotactic pattern formation plays a role in pattern transitions and irregularities. Variations are also found in the number and positioning of floral organs, although it is not known whether such variations provide any functional advantages. Two ways of diversification may be involved in angiosperm floral evolution: precise regulation of organ position and identity that leads to further specialization of organs and organ redundancy that leads to flexibility in floral phyllotaxis.

Genome-wide identification and comparative analysis of the WRKY gene family in aquatic plants and their response to abiotic stresses in giant duckweed (Spirodela polyrhiza)

WRKY is one of the largest transcription factor families across higher plant species and is involved in important biological processes and plant responses to various biotic/abiotic stresses. However, only a few investigations on WRKYs have been conducted in aquatic plants. This study first systematically analyzed the gene structure, protein properties, and phylogenetic relationship of 693 WRKYs in nine aquatic and two wetland plants at the genome-wide level. The pattern of WRKY groups in two aquatic ferns provided new evidence for the origin and evolution of WRKY genes. ARE cis-regulatory elements show an unusual high frequency in the promoter region of WRKY genes, indicating the adaptation to the aquatic habitat in aquatic plants. The WRKY gene family experienced a series of gene loss events in aquatic plants, especially group III. Further studies were conducted on the interaction network of SpWRKYs, their target genes, and non-coding RNAs. The expression profile of SpWRKYs under phosphate starvation, cold, and submergence conditions revealed that most SpWRKYs are involved in the response to abiotic stresses. Our investigations lay the foundation for further study on the mechanism of WRKYs responding to abiotic stresses in aquatic plants.

Evolution of the Subgroup 6 R2R3-MYB Genes and Their Contribution to Floral Color in the Perianth-Bearing Piperales

Flavonoids, carotenoids, betalains, and chlorophylls are the plant pigments responsible for floral color. Anthocyanins, a class of flavonoids, are largely responsible for the red, purple, pink, and blue colors. R2R3-MYB genes belonging to subgroup 6 (SG6) are the upstream regulatory factors of the anthocyanin biosynthetic pathway. The canonical members of these genes in Arabidopsis include AtMYB75, AtMYB90, AtMYB113, and AtMYB114. The Aristolochiaceae is an angiosperm lineage with diverse floral groundplans and perianth colors. Saruma henryi exhibits a biseriate perianth with green sepals and yellow petals. All other genera have sepals only, with colors ranging from green (in Lactoris) to a plethora of yellow to red and purple mixtures. Here, we isolated and reconstructed the SG6 R2R3-MYB gene lineage evolution in angiosperms with sampling emphasis in Aristolochiaceae. We found numerous species-specific duplications of this gene lineage in core eudicots and local duplications in Aristolochiaceae for Saruma and Asarum. Expression of SG6 R2R3-MYB genes examined in different developmental stages and plant organs of four Aristolochiaceae species, largely overlaps with red and purple pigments, suggesting a role in anthocyanin and flavonoid synthesis and accumulation. A directed RNA-seq analysis corroborated our RT-PCR analyses, by showing that these structural enzymes activate during perianth development in Aristolochia fimbriata and that the regulatory genes are expressed in correlation with color phenotype. Finally, the reconstruction of the flavonoid and anthocyanin metabolic pathways using predicted peptides from transcriptomic data show that all pivotal enzymes are present in the analyzed species. We conclude that the regulatory genes as well as the biosynthetic pathway are largely conserved across angiosperms. In addition, the Aristolochiaceae emerges as a remarkable group to study the genetic regulatory network for floral color, as their members exhibit an outstanding floral diversity with elaborate color patterns and the genetic complement for SG6 R2R3-MYB genes is simpler than in core eudicot model species.

Evolution-Driven Versatility of N Terminal Acetylation in Photoautotrophs

N terminal protein α-acetylation (NTA) is a pervasive protein modification that has recently attracted renewed interest. Early studies on NTA were mostly conducted in yeast and metazoans, providing a detailed portrait of the modification, which was indirectly applied to all eukaryotes. However, new findings originating from photosynthetic organisms have expanded our knowledge of this modification, revealing strong similarities as well as idiosyncratic features. Here, we review the most recent advances on NTA and its dedicated machinery in photosynthetic organisms. We discuss the cytosolic and unique plastid NTA machineries and their critical biological roles in development, stress responses, protein translocation, and stability. These new findings suggest that the multitasking plastid and cytosolic machineries evolved to support the specific needs of photoautotrophs.

Research advances in and prospects of ornamental plant genomics

The term 'ornamental plant' refers to all plants with ornamental value, which generally have beautiful flowers or special plant architectures. China is rich in ornamental plant resources and known as the "mother of gardens". Genomics is the science of studying genomes and is useful for carrying out research on genome evolution, genomic variations, gene regulation, and important biological mechanisms based on detailed genome sequence information. Due to the diversity of ornamental plants and high sequencing costs, the progress of genome research on ornamental plants has been slow for a long time. With the emergence of new sequencing technologies and a reduction in costs since the whole-genome sequencing of the first ornamental plant (Prunus mume) was completed in 2012, whole-genome sequencing of more than 69 ornamental plants has been completed in <10 years. In this review, whole-genome sequencing and resequencing of ornamental plants will be discussed. We provide analysis with regard to basic data from whole-genome studies of important ornamental plants, the regulation of important ornamental traits, and application prospects.

Fossil data support a pre-Cretaceous origin of flowering plants

Flowering plants (angiosperms) are the most diverse of all land plants, becoming abundant in the Cretaceous and achieving dominance in the Cenozoic. However, the exact timing of their origin remains a controversial topic, with molecular clocks generally placing their origin much further back in time than the oldest unequivocal fossils. To resolve this discrepancy, we developed a Bayesian method to estimate the ages of angiosperm families on the basis of the fossil record (a newly compiled dataset of ~15,000 occurrences in 198 families) and their living diversity. Our results indicate that several families originated in the Jurassic, strongly rejecting a Cretaceous origin for the group. We report a marked increase in lineage accumulation from 125 to 72 million years ago, supporting Darwin's hypothesis of a rapid Cretaceous angiosperm diversification. Our results demonstrate that a pre-Cretaceous origin of angiosperms is supported not only by molecular clock approaches but also by analyses of the fossil record that explicitly correct for incomplete sampling.

Angiosperm pollinivory in a Cretaceous beetle

Despite the crucial importance of flower-visiting insects in modern ecosystems, there is little fossil evidence on the origins of angiosperm pollination. Most reports of pollination in the Mesozoic fossil record have been based on the co-occurrence of the purported pollinators with pollen grains and assumed morphological adaptations for vectoring pollen. Here, we describe an exceptionally preserved short-winged flower beetle (Cucujoidea: Kateretidae) from mid-Cretaceous amber, Pelretes vivificus gen. et sp. nov., associated with pollen aggregations and coprolites consisting mainly of pollen, providing direct evidence of pollen-feeding in a Cretaceous beetle and confirming that diverse beetle lineages visited early angiosperms in the Cretaceous. The exquisite preservation of our fossil permits the identification of the pollen grains as Tricolpopollenites (Asteridae or Rosidae), representing a record of flower beetle pollination of a group of derived angiosperms in the Mesozoic and suggesting that potentially diverse beetle lineages visited early angiosperms by the mid-Cretaceous.

A novel predicted ADP-ribosyltransferase-like family conserved in eukaryotic evolution

The presence of many completely uncharacterized proteins, even in well-studied organisms such as humans, seriously hampers full understanding of the functioning of the living cells. ADP-ribosylation is a common post-translational modification of proteins; also nucleic acids and small molecules can be modified by the covalent attachment of ADP-ribose. This modification, important in cellular signalling and infection processes, is usually executed by enzymes from the large superfamily of ADP-ribosyltransferases (ARTs). Here, using bioinformatics approaches, we identify a novel putative ADP-ribosyltransferase family, conserved in eukaryotic evolution, with a divergent active site. The hallmark of these proteins is the ART domain nestled between flanking leucine-rich repeat (LRR) domains. LRRs are typically involved in innate immune surveillance. The novel family appears as putative novel ADP-ribosylation-related actors, most likely pseudoenzymes. Sequence divergence and lack of clearly detectable “classical” ART active site suggests the novel domains are pseudoARTs, yet atypical ART activity, or alternative enzymatic activity cannot be excluded. We propose that this family, including its human member LRRC9, may be involved in an ancient defense mechanism, with analogies to the innate immune system, and coupling pathogen detection to ADP-ribosyltransfer or other signalling mechanisms.

A high-quality Brassica napus genome reveals expansion of transposable elements, subgenome evolution and disease resistance

Rapeseed (Brassica napus L.) is a recent allotetraploid crop, which is well known for its high oil production. Here, we report a high-quality genome assembly of a typical semi-winter rapeseed cultivar, 'Zhongshuang11' (hereafter 'ZS11'), using a combination of single-molecule sequencing and chromosome conformation capture (Hi-C) techniques. Most of the high-confidence sequences (93.1%) were anchored to the individual chromosomes with a total of 19 centromeres identified, matching the exact chromosome count of B. napus. The repeat sequences in the A and C subgenomes in B. napus expanded significantly from 500 000 years ago, especially over the last 100 000 years. These young and recently amplified LTR-RTs showed dispersed chromosomal distribution but significantly preferentially clustered into centromeric regions. We exhaustively annotated the nucleotide-binding leucine-rich repeat (NLR) gene repertoire, yielding a total of 597 NLR genes in B. napus genome and 17.4% of which are paired (head-to-head arrangement). Based on the resequencing data of 991 B. napus accessions, we have identified 18 759 245 single nucleotide polymorphisms (SNPs) and detected a large number of genomic regions under selective sweep among the three major ecotype groups (winter, semi-winter and spring) in B. napus. We found 49 NLR genes and five NLR gene pairs colocated in selective sweep regions with different ecotypes, suggesting a rapid diversification of NLR genes during the domestication of B. napus. The high quality of our B. napus 'ZS11' genome assembly could serve as an important resource for the study of rapeseed genomics and reveal the genetic variations associated with important agronomic traits.

Impact of short-read sequencing on the misassembly of a plant genome

Background

Availability of plant genome sequences has led to significant advances. However, with few exceptions, the great majority of existing genome assemblies are derived from short read sequencing technologies with highly uneven read coverages indicative of sequencing and assembly issues that could significantly impact any downstream analysis of plant genomes. In tomato for example, 0.6% (5.1 Mb) and 9.7% (79.6 Mb) of short-read based assembly had significantly higher and lower coverage compared to background, respectively.

Results

To understand what the causes may be for such uneven coverage, we first established machine learning models capable of predicting genomic regions with variable coverages and found that high coverage regions tend to have higher simple sequence repeat and tandem gene densities compared to background regions. To determine if the high coverage regions were misassembled, we examined a recently available tomato long-read based assembly and found that 27.8% (1.41 Mb) of high coverage regions were potentially misassembled of duplicate sequences, compared to 1.4% in background regions. In addition, using a predictive model that can distinguish correctly and incorrectly assembled high coverage regions, we found that misassembled, high coverage regions tend to be flanked by simple sequence repeats, pseudogenes, and transposon elements.

Conclusions

Our study provides insights on the causes of variable coverage regions and a quantitative assessment of factors contributing to plant genome misassembly when using short reads and the generality of these causes and factors should be tested further in other species.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07397-5.

Up-regulated 2-alkenal reductase expression improves low-nitrogen tolerance in maize by alleviating oxidative stress

In plants, cellular lipid peroxidation is enhanced under low nitrogen (LN) stress; this increases the lipid-derived reactive carbonyl species (RCS) levels. The cellular toxicity of RCS can be reduced by various RCS-scavenging enzymes. However, the roles of these enzymes in alleviating oxidative stress and improving nutrient use efficiency (NUE) under nutrient stress remain unknown. Here, we overexpressed maize endogenous NADPH-dependent 2-alkenal reductase (ZmAER) in maize; it significantly increased the tolerance of transgenic plants (OX-AER) to LN stress. Under LN condition, the biomass, nitrogen accumulation, NUE, and leaf photosynthesis of the OX-AER plants were significantly higher than those of the wild-type (WT) plants. The leaf and root malondialdehyde and H₂ O₂ levels in the transgenic plants were significantly lower than those in WT. The expression of antioxidant enzyme-related genes ZmCAT3, ZmPOD5 and ZmPOD13 was significantly higher in the transgenic lines than in WT. Under LN stress, the nitrate reductase activity in the OX-AER leaves was significantly increased compared with that in the WT leaves. Furthermore, under LN stress, ZmNRT1.1 and ZmNRT2.5 expression was upregulated in the OX-AER plants compared with that in WT. Overall, up-regulated ZmAER expression could enhance maize's tolerance to LN stress by alleviating oxidative stress and improve NUE.

The pineapple MADS-box gene family and the evolution of early monocot flower

Unlike the flower of the model monocot rice, which has diverged greatly from the ancestral monocot flower, the pineapple (Ananas comosus) flower is more typical of monocot flowers. Here, we identified 43 pineapple genes containing MADS-box domains, including 11 type I and 32 type II genes. RNA-seq expression data generated from five pineapple floral organs (sepals, petals, stamens, pistils, and ovules) and quantitative real-time PCR revealed tissue-specific expression patterns for some genes. We found that AcAGL6 and AcFUL1 were mainly expressed in sepals and petals, suggesting their involvement in the regulation of these floral organs. A pineapple ‘ABCDE’ model was proposed based on the phylogenetic analysis and expression patterns of MADS-box genes. Unlike rice and orchid with frequent species-specific gene duplication and subsequent expression divergence, the composition and expression of the ABCDE genes were conserved in pineapple. We also found that AcSEP1/3, AcAG, AcAGL11a/b/c, and AcFUL1 were highly expressed at different stages of fruit development and have similar expression profiles, implicating these genes’ role in fruit development and ripening processes. We propose that the pineapple flower can be used as a model for studying the ancestral form of monocot flowers to investigate their development and evolutionary history.

Piper nigrum CYP719A37 Catalyzes the Decisive Methylenedioxy Bridge Formation in Piperine Biosynthesis

Black pepper (Piper nigrum) is among the world's most popular spices. Its pungent principle, piperine, has already been identified 200 years ago, yet the biosynthesis of piperine in black pepper remains largely enigmatic. In this report we analyzed the characteristic methylenedioxy bridge formation of the aromatic part of piperine by a combination of RNA-sequencing, functional expression in yeast, and LC-MS based analysis of substrate and product profiles. We identified a single cytochrome P450 transcript, specifically expressed in black pepper immature fruits. The corresponding gene was functionally expressed in yeast (Saccharomyces cerevisiae) and characterized for substrate specificity with a series of putative aromatic precursors with an aromatic vanilloid structure. Methylenedioxy bridge formation was only detected when feruperic acid (5-(4-hydroxy-3-methoxyphenyl)-2,4-pentadienoic acid) was used as a substrate, and the corresponding product was identified as piperic acid. Two alternative precursors, ferulic acid and feruperine, were not accepted. Our data provide experimental evidence that formation of the piperine methylenedioxy bridge takes place in young black pepper fruits after a currently hypothetical chain elongation of ferulic acid and before the formation of the amide bond. The partially characterized enzyme was classified as CYP719A37 and is discussed in terms of specificity, storage, and phylogenetic origin of CYP719 catalyzed reactions in magnoliids and eudicots.

Genome-wide identification and characterization of small auxin-up RNA (SAUR) gene family in plants: evolution and expression profiles during normal growth and stress response

BACKGROUND

Auxin is critical to plant growth and development, as well as stress responses. Small auxin-up RNA (SAUR) is the largest family of early auxin responsive genes in higher plants. However, the function of few SAUR genes is known owing to functional redundancy among the many family members.

RESULTS

In this study, we conducted a phylogenetic analysis using protein sequences of 795 SAURs from Anthoceros angustus, Marchantia polymorpha, Physcomitrella patens, Selaginella moellendorffii, Ginkgo biloba, Gnetum montanum, Amborella trichopoda, Arabidopsis thaliana, Oryza sativa, Zea mays, Glycine max, Medicago truncatula and Setaria italica. The phylogenetic trees showed that the SAUR proteins could be divided into 10 clades and three subfamilies, and that SAUR proteins of three bryophyte species were only located in subfamily III, which suggested that they may be ancestral. From bryophyta to anthophyta, SAUR family have appeared very large expansion. The number of SAUR gene in Fabaceae species was considerably higher than that in other plants, which may be associated with independent whole genome duplication event in the Fabaceae lineages. The phylogenetic trees also showed that SAUR genes had expanded independently monocotyledons and dicotyledons in angiosperms. Conserved motif and protein structure prediction revealed that SAUR proteins were highly conserved among higher plants, and two leucine residues in motif I were observed in almost all SAUR proteins, which suggests the residues plays a critical role in the stability and function of SAUR proteins. Expression analysis of SAUR genes using publicly available RNA-seq data from rice and soybean indicated functional similarity of members in the same clade, which was also further confirmed by qRT-PCR. Summarization of SAUR functions also showed that SAUR functions were usually consistent within a subclade.

CONCLUSIONS

This study provides insights into the evolution and function of the SAUR gene family from bryophyta to anthophyta, particularly in Fabaceae plants. Future investigation to understand the functions of SAUR family members should employ a clade as the study unit.

The RGF/GLV/CLEL Family of Short Peptides Evolved Through Lineage-Specific Losses and Diversification and Yet Conserves Its Signaling Role Between Vascular Plants and Bryophytes

Short secreted plant peptides act as key signaling molecules and control a plethora of developmental and physiological processes. The ROOT GROWTH FACTOR (RGF)/GOLVEN (GLV)/CLE-Like (CLEL) family of peptides was discovered to be involved in root development in Arabidopsis thaliana. In contrast to active research efforts, which have been revealing receptors and downstream signaling components, little attention has been paid to evolutionary processes that shaped the RGF signaling system as we know it in angiosperms today. As a first step toward understanding how RGF signaling emerged and evolved, this study aimed to elucidate the phylogenetic distribution and functional conservation of RGF-like sequences. Using publicly available, genome and transcriptome data, RGF-like sequences were searched in 27 liverworts, 22 mosses, 8 hornworts, 23 lycophytes, 23 ferns, 38 gymnosperms, and 8 angiosperms. This led to the identification of more than four hundreds of RGF-like sequences in all major extant land plant lineages except for hornworts. Sequence comparisons within and between taxonomic groups identified lineage-specific characters. Notably, one of the two major RGF subgroups, represented by A. thaliana RGF6/GLV1/CLEL6, was found only in vascular plants. This subgroup, therefore, likely emerged in a common ancestor of vascular plants after its divergence from bryophytes. In bryophytes, our results infer independent losses of RGF-like sequences in mosses and hornworts. On the other hand, a single, highly similar RGF-like sequence is conserved in liverworts, including Marchantia polymorpha, a genetically tractable model species. When constitutively expressed, the M. polymorpha RGF-like sequence (MpRGF) affected plant development and growth both in A. thaliana and M. polymorpha. This suggests that MpRGF can exert known RGF-like effects and that MpRGF is under transcriptional control so that its potent activities are precisely controlled. These data suggest that RGFs are conserved as signaling molecules in both vascular plants and bryophytes and that lineage-specific diversification has increased sequence variations of RGFs. All together, our findings form a basis for further studies into RGF peptides and their receptors, which will contribute to our understandings of how peptide signaling pathways evolve.

Functioning of PPR Proteins in Organelle RNA Metabolism and Chloroplast Biogenesis

The pentatricopeptide repeat (PPR) proteins constitute one of the largest nuclear-encoded protein families in higher plants, with over 400 members in most sequenced plant species. The molecular functions of these proteins and their physiological roles during plant growth and development have been widely studied. Generally, there is mounting evidence that PPR proteins are involved in the post-transcriptional regulation of chloroplast and/or mitochondrial genes, including RNA maturation, editing, intron splicing, transcripts' stabilization, and translation initiation. The cooperative action of RNA metabolism has profound effects on the biogenesis and functioning of both chloroplasts and mitochondria and, consequently, on the photosynthesis, respiration, and development of plants and their environmental responses. In this review, we summarize the latest research on PPR proteins, specifically how they might function in the chloroplast, by documenting their mechanism of molecular function, their corresponding RNA targets, and their specific effects upon chloroplast biogenesis and host organisms.

A reappraisal of the phylogenetic placement of the Aquilegia whole-genome duplication

The accurate placement of an ancient whole-genome duplication (WGD) in relation to the lineage divergence is important. Here, we re-investigated the Aquilegia coerulea WGD and found it is more likely lineage-specific rather than shared by all eudicots.

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

Background

Tetracentron sinense is an endemic and endangered deciduous tree. It belongs to the Trochodendrales, one of four early diverging lineages of eudicots known for having vesselless secondary wood. Sequencing and resequencing of the T. sinense genome will help us understand eudicot evolution, the genetic basis of tracheary element development, and the genetic diversity of this relict species.

Results

Here, we report a chromosome-scale assembly of the T. sinense genome. We assemble the 1.07 Gb genome sequence into 24 chromosomes and annotate 32,690 protein-coding genes. Phylogenomic analyses verify that the Trochodendrales and core eudicots are sister lineages and showed that two whole-genome duplications occurred in the Trochodendrales approximately 82 and 59 million years ago. Synteny analyses suggest that the γ event, resulting in paleohexaploidy, may have only happened in core eudicots. Interestingly, we find that vessel elements are present in T. sinense, which has two orthologs of AtVND7, the master regulator of vessel formation. T. sinense also has several key genes regulated by or regulating TsVND7.2 and their regulatory relationship resembles that in Arabidopsis thaliana. Resequencing and population genomics reveals high levels of genetic diversity of T. sinense and identifies four refugia in China.

Conclusions

The T. sinense genome provides a unique reference for inferring the early evolution of eudicots and the mechanisms underlying vessel element formation. Population genomics analysis of T. sinense reveals its genetic diversity and geographic structure with implications for conservation.

Supplementary information

The online version contains supplementary material available at 10.1186/s13059-020-02198-7.

The ancient wave of polyploidization events in flowering plants and their facilitated adaptation to environmental stress

Flowering plants, or angiosperms, consist of more than 300,000 species, far more than any other land plant lineages. The accumulated evidence indicates that multiple ancient polyploidy events occurred around 100 to 120 million years ago during the Cretaceous and drove the early diversification of four major clades of angiosperms: gamma whole-genome triplication in the common ancestor of core eudicots, tau whole-genome duplication during the early diversification of monocots, lambda whole-genome duplication during the early diversification of magnoliids, and pi whole-genome duplication in the Nymphaeales lineage. These four polyploidy events have played essential roles in the adaptive evolution and diversification of major clades of flowering plants. Here, we specifically review the current understanding of this wave of ancient whole-genome duplications and their evolutionary significance. Notably, although these ancient whole-genome duplications occurred independently, they have contributed to the expansion of many stress-related genes (e.g., heat shock transcription factors and Arabidopsis response regulators)，and these genes could have been selected for by global environmental changes in the Cretaceous. Therefore, this ancient wave of paleopolyploidy events could have significantly contributed to the adaptation of angiosperms to environmental changes, and potentially promoted the wide diversification of flowering plants.

Draft genome of the aquatic moss Fontinalis antipyretica (Fontinalaceae, Bryophyta)

Mosses comprise one of three lineages forming a sister group to extant vascular plants. Having emerged from an early split in the diversification of embryophytes, mosses may offer complementary insights into the evolution of traits following the transition to, and colonization of, land. Here, we report the draft nuclear genome of Fontinalis antipyretica (Fontinalaceae, Hypnales), a charismatic aquatic moss that is widespread in temperate regions of the Northern Hemisphere. We sequenced and de novo-assembled its genome using the 10X Genomics method. The genome comprises 385.2 Mbp, with a scaffold N50 of 45.8 Kbp. The assembly captured 87.2% of the 430 genes in the BUSCO Viridiplantae odb10 dataset. The newly generated F. antipyretica genome is the third moss genome, and the second seedless aquatic plant genome, to be sequenced and assembled to date.

A study of male fertility control in Medicago truncatula uncovers an evolutionarily conserved recruitment of two tapetal bHLH subfamilies in plant sexual reproduction

Male sterility is an important tool for plant breeding and hybrid seed production. Male-sterile mutants are largely due to an abnormal development of either the sporophytic or gametophytic anther tissues. Tapetum, a key sporophytic tissue, provides nutrients for pollen development, and its delayed degeneration induces pollen abortion. Numerous bHLH proteins have been documented to participate in the degeneration of the tapetum in angiosperms, but relatively little attention has been given to the evolution of the involved developmental pathways across the phylogeny of land plants. A combination of cellular, molecular, biochemical and evolutionary analyses was used to investigate the male fertility control in Medicago truncatula. We characterized the male-sterile mutant empty anther1 (ean1) and identified EAN1 as a tapetum-specific bHLH transcription factor necessary for tapetum degeneration. Our study uncovered an evolutionarily conserved recruitment of bHLH subfamily II and III(a + c)1 in the regulation of tapetum degeneration. EAN1 belongs to the subfamily II and specifically forms heterodimers with the subfamily III(a + c)1 members, which suggests a heterodimerization mechanism conserved in angiosperms. Our work suggested that the pathway of two tapetal-bHLH subfamilies is conserved in all land plants, and likely was established before the divergence of the spore-producing land plants.

A theoretical morphological model for quantitative description of the three-dimensional floral morphology in water lily (Nymphaea)

Water lilies (Nymphaea spp.) have diverse floral morphologies. Water lilies are not only commonly used as ornamental plants, but they are also important for understanding the diversification of basal angiosperms. Although the diversity in floral morphology of water lily provides useful information for evolutionary biology, horticulture, and horticultural science, it is difficult to describe and analyze the three-dimensional morphology of flowers. In this study, we propose a method to describe the floral morphology of water lily using a three-dimensional theoretical morphological model. The theoretical model was constructed based on three components, i.e., (1) the gradual change in size of floral organs, (2) spiral phyllotaxis, and (3) the interpolation of elevation angles, which were integrated into the model. We generated three-dimensional representation of water lily flowers and visualized theoretical morphospaces by varying each morphological parameter. The theoretical morphospace is a mathematical space of morphological spectrum generated by a theoretical morphological model. These morphospaces seems to display the large part of morphological variations of water lily. We measured morphological parameters of real flowers based on our theoretical model and display the occupation pattern of morphological parameters. We also surveyed the relation between morphological parameters and flower shape descriptions found in a catalog. In some parameters, we found breeders’ description can link to our morphological model. In addition, the relationship between the global features of floral morphology and the parameters of the theoretical model was calculated with flower silhouettes simulated with a range of parameter values and the global features of the silhouette. We used two simple indices to assess the global morphological features, which were calculated with the convex hull. The results indicated that our method can effectively provide an objective and quantitative overview of the diversity in the floral morphology of water lily.

Genome-wide identification and expression analysis of SABATH methyltransferases in tea plant (Camellia sinensis): insights into their roles in plant defense responses

SABATH methyltransferases convent plant small-molecule metabolites into volatile methyl esters, which play important roles in many biological processes and defense reactions in plants. In this study, a total of 32 SABATH genes were identified in the Camellia sinensis var. sinensis (CSS) genome, which were renamed CsSABATH1 to CsSABATH32. Genome location annotation suggested that tandem duplication was responsible for the expansion of SABATH genes in tea plant. Multiple sequence alignment and phylogenetic analysis showed that the CsSABATHs could be classified into three groups (I, II and III), which were also supported by gene structures and conserved motifs analysis. Group II contained only two CsSABATH proteins, which were closely related to PtIAMT, AtIAMT and OsIAMT. The group III SABATH genes of tea plant exhibited expansion on the CSS genome compared with Camellia sinensis var. assamica (CSA) genome. Based on RNA-seq data, the CsSABATHs exhibited tissue-specific expression patterns, and the members with high expression in buds and young leaves were also obviously upregulated after MeJA treatment. The expression of many transcription factors was significantly correlated with that of different members of the CsSABATH gene family, suggesting a potential regulatory relationship between them. Quantitative real-time PCR (qPCR) expression analysis showed that CsSABATHs could respond to exogenous JA, SA and MeSA treatments in tea plants. RNA-seq data analysis and qPCR validation suggested that CsSABATH8, 11, 16, 25, 29 and 32 might play a special role in plant defense against insect herbivory. These results provide references for evolutionary studies of the plant SABATH family and the exploration of the potential roles of CsSABATHs in tea plant defense responses.

What can lycophytes teach us about plant evolution and development? Modern perspectives on an ancient lineage

All Evo-Devo studies rely on representative sampling across the tree of interest to elucidate evolutionary trajectories through time. In land plants, genetic resources are well established in model species representing lineages including bryophytes (mosses, liverworts, and hornworts), monilophytes (ferns and allies), and seed plants (gymnosperms and flowering plants), but few resources are available for lycophytes (club mosses, spike mosses, and quillworts). Living lycophytes are a sister group to the euphyllophytes (the fern and seed plant clade), and have retained several ancestral morphological traits despite divergence from a common ancestor of vascular plants around 420 million years ago. This sister relationship offers a unique opportunity to study the conservation of traits such as sporophyte branching, vasculature, and indeterminacy, as well as the convergent evolution of traits such as leaves and roots which have evolved independently in each vascular plant lineage. To elucidate the evolution of vascular development and leaf formation, molecular studies using RNA Seq, quantitative reverse transcription polymerase chain reaction, in situ hybridisation and phylogenetics have revealed the diversification and expression patterns of KNOX, ARP, HD-ZIP, KANADI, and WOX gene families in lycophytes. However, the molecular basis of further trait evolution is not known. Here we describe morphological traits of living lycophytes and their extinct relatives, consider the molecular underpinnings of trait evolution and discuss future research required in lycophytes to understand the key evolutionary innovations enabling the growth and development of all vascular plants.

Distinct Expression and Methylation Patterns for Genes with Different Fates following a Single Whole-Genome Duplication in Flowering Plants

For most sequenced flowering plants, multiple whole-genome duplications (WGDs) are found. Duplicated genes following WGD often have different fates that can quickly disappear again, be retained for long(er) periods, or subsequently undergo small-scale duplications. However, how different expression, epigenetic regulation, and functional constraints are associated with these different gene fates following a WGD still requires further investigation due to successive WGDs in angiosperms complicating the gene trajectories. In this study, we investigate lotus (Nelumbo nucifera), an angiosperm with a single WGD during the K-pg boundary. Based on improved intraspecific-synteny identification by a chromosome-level assembly, transcriptome, and bisulfite sequencing, we explore not only the fundamental distinctions in genomic features, expression, and methylation patterns of genes with different fates after a WGD but also the factors that shape post-WGD expression divergence and expression bias between duplicates. We found that after a WGD genes that returned to single copies show the highest levels and breadth of expression, gene body methylation, and intron numbers, whereas the long-retained duplicates exhibit the highest degrees of protein-protein interactions and protein lengths and the lowest methylation in gene flanking regions. For those long-retained duplicate pairs, the degree of expression divergence correlates with their sequence divergence, degree in protein-protein interactions, and expression level, whereas their biases in expression level reflecting subgenome dominance are associated with the bias of subgenome fractionation. Overall, our study on the paleopolyploid nature of lotus highlights the impact of different functional constraints on gene fate and duplicate divergence following a single WGD in plant.

The Chimonanthus salicifolius genome provides insight into magnoliid evolution and flavonoid biosynthesis

Chimonanthus salicifolius, a member of the Calycanthaceae of magnoliids, is one of the most famous medicinal plants in Eastern China. Here, we report a chromosome-level genome assembly of C. salicifolius, comprising 820.1 Mb of genomic sequence with a contig N50 of 2.3 Mb and containing 36 651 annotated protein-coding genes. Phylogenetic analyses revealed that magnoliids were sister to the eudicots. Two rounds of ancient whole-genome duplication were inferred in the C. salicifolious genome. One is shared by Calycanthaceae after its divergence with Lauraceae, and the other is in the ancestry of Magnoliales and Laurales. Notably, long genes with > 20 kb in length were much more prevalent in the magnoliid genomes compared with other angiosperms, which could be caused by the length expansion of introns inserted by transposon elements. Homologous genes within the flavonoid pathway for C. salicifolius were identified, and correlation of the gene expression and the contents of flavonoid metabolites revealed potential critical genes involved in flavonoids biosynthesis. This study not only provides an additional whole-genome sequence from the magnoliids, but also opens the door to functional genomic research and molecular breeding of C. salicifolius.