Water lilies as emerging models for Darwin's abominable mystery

Comparative transcriptomics of seed nourishing tissues: uncovering conserved and divergent pathways in seed plants

The evolutionary and ecological success of spermatophytes is intrinsically linked to the seed habit, which provides a protective environment for the initial development of the new generation. This environment includes an ephemeral nourishing tissue that supports embryo growth. In gymnosperms this tissue originates from the asexual proliferation of the maternal megagametophyte, while in angiosperms it is a product of fertilization, and is called the endosperm. The emergence of these nourishing tissues is of profound evolutionary value, and they are also food staples for most of the world's population. Here, using Orthofinder to infer orthologue genes among newly generated and previously published datasets, we provide a comparative transcriptomic analysis of seed nourishing tissues from species of several angiosperm clades, including those of early diverging lineages, as well as of one gymnosperm. Our results show that, although the structure and composition of seed nourishing tissues has seen significant divergence along evolution, there are signatures that are conserved throughout the phylogeny. Conversely, we identified processes that are specific to species within the clades studied, and thus illustrate their functional divergence. With this, we aimed to provide a foundation for future studies on the evolutionary history of seed nourishing structures, as well as a resource for gene discovery in future functional studies.

Multi-omic dissection of the cold resistance traits of white water lily

The white water lily (Nymphaea candida), exemplifying nature's resilience, thrives in the high-altitude terrains of Xinjiang, China, serving as an ideal model for investigating cold adaptation mechanisms in aquatic plants. This study meticulously elucidates the complex cold adaptation mechanisms of the white water lily through a comprehensive and integrated methodological approach. We discovered that the water lily undergoes ecodormancy in winter, retaining high cellular viability and growth potential. During overwintering, the white water lily demonstrates effective resource reallocation, a process facilitated by morphological adjustments, thereby strengthening its resistance to cold temperatures. This enhancement is achieved particularly through the compartmentalization of large vacuoles, the accumulation of osmoregulatory substances, and an increased antioxidant capacity. We established the first exhaustive full-length transcriptome for the white water lily. A subsequent comprehensive analysis of the transcriptome, phytohormones, and metabolome uncovered a multifaceted regulatory network orchestrating cold adaptation. Our research spotlights phytohormone signaling, amino acid metabolism, and circadian rhythms as key elements in the water lily's defense against cold. The results emphasize the critical role of nitrogen metabolism, especially amino acid-related pathways, during cold stress. Metabolite profiling revealed the importance of compounds like myo-inositol and L-proline in enhancing cold tolerance. Remarkably, our study demonstrates that the white water lily notably diminishes the utilization of unsaturated fatty acids in its temperature regulation strategies. In conclusion, this research substantially enriches our understanding of the white water lily's intricate cold adaptation mechanisms, offering new perspectives on the adaptive strategies of aquatic plants and potential applications in agricultural advancement.

Diverse O-methyltransferases catalyze the biosynthesis of floral benzenoids that repel aphids from the flowers of waterlily Nymphaea prolifera

Nymphaea is a key genus of the ANA grade (Amborellales, Nymphaeales, and Austrobaileyales) of basal flowering plants, which serve as a key model to study the early evolution of floral traits. In this study, we comprehensively investigated the emission, biosynthesis, and biological function of the floral scent in a night-blossoming waterlily Nymphaea prolifera. The headspace volatile collection combined with GC-MS analysis showed that the floral scent of N. prolifera is predominately comprised by methylated benzenoids including anisole, veratrole, guaiacol, and methoxyanisole. Moreover, the emission of these floral benzenoids in N. prolifera exhibited temporal and spatial pattern with circadian rhythm and tissue specificity. By creating and mining transcriptomes of N. prolifera flowers, 12 oxygen methyltransferases (NpOMTs) were functionally identified. By in vitro enzymatic assay, NpOMT3, 6, and 7 could produce anisole and NpOMT5, 7, 9, produce guaiacol, whereas NpOMT3, 6, 9, 11 catalyzed the formation of veratrole. Methoxyanisole was identified as the universal product of all NpOMTs. Expression patterns of NpOMTs provided implication for their roles in the production of the respective benzenoids. Phylogenetic analysis of OMTs suggested a Nymphaea-specific expansion of the OMT family, indicating the evolution of lineage-specific functions. In bioassays, anisole, veratrole, and guaiacol in the floral benzenoids were revealed to play the critical role in repelling waterlily aphids. Overall, this study indicates that the basal flowering plant N. prolifera has evolved a diversity and complexity of OMT genes for the biosynthesis of methylated benzenoids that can repel insects from feeding the flowers. These findings provide new insights into the evolutional mechanism and ecological significance of the floral scent from early-diverged flowering plants.

Insights into Adaptive Regulation of the Leaf-Petiole System: Strategies for Survival of Water Lily Plants under Salt Stress

The water lily (Nymphaea tetragona) is an ancient angiosperm that belongs to the Nymphaeaceae family. As a rooted floating-leaf plant, water lilies are generally cultivated in fresh water, therefore, little is known about their survival strategies under salt stress. Long-term salt stress causes morphological changes, such as the rapid regeneration of floating leaves and a significant decrease in leaf number and surface area. We demonstrate that salt stress induces toxicity soon after treatment, but plants can adapt by regenerating floating leaves that are photosynthetically active. Transcriptome profiling revealed that ion binding was one of the most-enriched GO terms in leaf-petiole systems under salt stress. Sodium-transporter-related genes were downregulated, whereas K+ transporter genes were both up- and downregulated. These results suggest that restricting intracellular Na+ importing while maintaining balanced K+ homeostasis is an adaptive strategy for tolerating long-term salt stress. ICP-MS analysis identified the petioles and leaves as Na-hyperaccumulators, with a maximum content of over 80 g kg−1 DW under salt stress. Mapping of the Na-hyperaccumulation trait onto the phylogenetic relationships revealed that water lily plants might have a long evolutionary history from ancient marine plants, or may have undergone historical ecological events from salt to fresh water. Ammonium transporter genes involved in nitrogen metabolism were downregulated, whereas NO3−-related transporters were upregulated in both the leaves and petioles, suggesting a selective bias toward NO3− uptake under salt stress. The morphological changes we observed may be due to the reduced expression of genes related to auxin signal transduction. In conclusion, the floating leaves and submerged petioles of the water lily use a series of adaptive strategies to survive salt stress. These include the absorption and transport of ions and nutrients from the surrounding environments, and the ability to hyperaccumulate Na+. These adaptations may serve as the physiological basis for salt tolerance in water lily plants.

The complete chloroplast genome of Nymphaea atrans (Surrey Wilfrid Laurance Jacobs, 1992: Nymphaeaceae)

Nymphaea atrans belongs to the subgenus Anecphya and displays varied flower colors over successive days. Because of its excellent ornamental characteristics, this species is widely cultivated in waterscape gardens worldwide. Here, we have sequenced the complete chloroplast genome of N. atrans. The whole genome size is 160,990 bp in length with four subregions: large and small single-copy regions of 90,879 and 19,699 bp, respectively, separated by a pair of inverted repeat regions of 25,206 bp each. A total of 126 genes were annotated, including 82 coding genes, eight ribosomal RNAs, and 36 transfer RNAs. The total GC content of the complete genome was 39%. Phylogenetic analysis showed that N. atrans was closely related to N. immutabilis. In this study, we have provided the chloroplast genome of N. atrans as a valuable resource for further phylogenetic analysis of Nymphaea species.

Phylogenetic studies of magnoliids: Advances and perspectives

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

Genome-wide identification of YABBY genes in three Cymbidium species and expression patterns in C. ensifolium (Orchidaceae)

Members of the YABBY gene family play significant roles in lamina development in cotyledons, floral organs, and other lateral organs. The Orchidaceae family is one of the largest angiosperm groups. Some YABBYs have been reported in Orchidaceae. However, the function of YABBY genes in Cymbidium is currently unknown. In this study, 24 YABBY genes were identified in Cymbidium ensifolium, C. goeringii, and C. sinense. We analyzed the conserved domains and motifs, the phylogenetic relationships, chromosome distribution, collinear correlation, and cis-elements of these three species. We also analyzed expression patterns of C. ensifolium and C. goeringii. Phylogenetic relationships analysis indicated that 24 YABBY genes were clustered in four groups, INO, CRC/DL, YAB2, and YAB3/FIL. For most YABBY genes, the zinc finger domain was located near the N-terminus and the helix-loop-helix domain (YABBY domain) near the C-terminus. Chromosomal location analysis results suggested that only C. goeringii YABBY has tandem repeat genes. Almost all the YABBY genes displayed corresponding one-to-one relationships in the syntenic relationships analysis. Cis-elements analysis indicated that most elements were clustered in light-responsive elements, followed by MeJA-responsive elements. Expression patterns showed that YAB2 genes have high expression in floral organs. RT-qPCR analysis showed high expression of CeYAB3 in lip, petal, and in the gynostemium. CeCRC and CeYAB2.2 were highly expressed in gynostemium. These findings provide valuable information of YABBY genes in Cymbidium species and the function in Orchidaceae.

Phylogenomics as an effective approach to untangle cross-species hybridization event: A case study in the family Nymphaeaceae

Hybridization is common and considered as an important evolutionary force to increase intraspecific genetic diversity. Detecting hybridization events is crucial for understanding the evolutionary history of species and further improving molecular breeding. The studies on identifying hybridization events through the phylogenomic approach are still limited. We proposed the conception and method of identifying allopolyploidy events by phylogenomics. The reconciliation and summary of nuclear multi-labeled gene family trees were adopted to untangle hybridization events from next-generation data in our novel phylogenomic approach. Given horticulturalists’ relatively clear cultivated crossbreeding history, the water lily family is a suitable case for examining recent allopolyploidy events. Here, we reconstructed and confirmed the well-resolved nuclear phylogeny for the Nymphaeales family in the context of geological time as a framework for identifying hybridization signals. We successfully identified two possible allopolyploidy events with the parental lineages for the hybrids in the family Nymphaeaceae based on summarization from multi-labeled gene family trees of Nymphaeales. The lineages where species Nymphaea colorata and Nymphaea caerulea are located may be the progenitors of horticultural cultivated species Nymphaea ‘midnight’ and Nymphaea ‘Woods blue goddess’. The proposed hybridization hypothesis is also supported by horticultural breeding records. Our methodology can be widely applied to identify hybridization events and theoretically facilitate the genome breeding design of hybrid plants.

Differential efficacy of water lily cultivars in phytoremediation of eutrophic water contaminated with phosphorus and nitrogen

Excessive inputs of phosphorus (P) and nitrogen (N) trigger eutrophication of the water bodies, which promotes the undesirable growth of algal bloom and deterioration of the water quality, and aquatic biodiversity. Macrophytes provide an environmentally benign and economically viable paradigm for the ecological restoration of eutrophic aquatic ecosystems. Water lily (Nymphaea) is largely used as ornamental plant for landscaping, and it has been documented that water lily possesses the potentiality in mitigating polluted aquatic environments. In the present study, water lily cultivars Nymphaea Texas Dawn (NTD), Nymphaea Colorado (NC), Nymphaea Madame Wilfron Gonnère (NMWG), and Nymphaea Sunshine Princess (NSP) were investigated for their potency in alleviating the eutrophication. The concentrations of total P and total N were significantly higher in the leaves of NC and NSP compared with NTD and NMWG. Therefore, NC and NSP were selected for subsequent studies to decipher their recuperation efficacy on eutrophic waters at different growth stages. NC and NSP significantly reduced the concentrations of eutrophication indicators i.e., total P, NH₄⁺-N, and chemical oxygen demand in different gradients of the simulated eutrophic water in a growth-dependent manner. On the contrary, NC and NSP triggered a significant increase in the concentration of dissolved oxygen particularly at the seedling stage. Notably, the concentrations of total P (shoot and root) and total N (root) were relatively higher in NSP than NC. The study thus revealed a growth-dependent differential efficacy of NSP and NC in mitigating the different eutrophic waters.

Growth and Nutritional Quality of Lemnaceae Viewed Comparatively in an Ecological and Evolutionary Context

This review focuses on recently characterized traits of the aquatic floating plant Lemna with an emphasis on its capacity to combine rapid growth with the accumulation of high levels of the essential human micronutrient zeaxanthin due to an unusual pigment composition not seen in other fast-growing plants. In addition, Lemna's response to elevated CO₂ was evaluated in the context of the source-sink balance between plant sugar production and consumption. These and other traits of Lemnaceae are compared with those of other floating aquatic plants as well as terrestrial plants adapted to different environments. It was concluded that the unique features of aquatic plants reflect adaptations to the freshwater environment, including rapid growth, high productivity, and exceptionally strong accumulation of high-quality vegetative storage protein and human antioxidant micronutrients. It was further concluded that the insensitivity of growth rate to environmental conditions and plant source-sink imbalance may allow duckweeds to take advantage of elevated atmospheric CO₂ levels via particularly strong stimulation of biomass production and only minor declines in the growth of new tissue. It is proposed that declines in nutritional quality under elevated CO₂ (due to regulatory adjustments in photosynthetic metabolism) may be mitigated by plant-microbe interaction, for which duckweeds have a high propensity.

The complete chloroplast genome of Nymphaea thermarum (Nymphaeaceae) from Rwanda, Africa

Nymphaea thermarum is classified in the Nymphaeaceae, and is the smallest water lily in the world. It has been extinct its native environment and needs urgent protection. Here, we report and characterize the complete chloroplast genome of N. thermarum. The total length of the chloroplast genome is 159,849 bp and the GC content is 39.2% (A: 30.1%, C: 20.0%, G: 19.2%, T: 30.8%). The chloroplast genome consists of 8 rRNA, 37 tRNA, and 85 protein-coding genes. Phylogenetic analysis of N. thermarum fully resolved this taxon in a clade with Nymphaea capensis. The chloroplast genome of N. thermarum provides scientific guidance for its conservation genetics and also contributes genome resources for the phylogenetic relationship of Nymphaea.

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.

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.

Genome-Wide Identification of the YABBY Gene Family in Seven Species of Magnoliids and Expression Analysis in Litsea

The YABBY gene family, specific to seed plants, encodes a class of transcription factors in the lamina maintenance and development of lateral organs. Magnoliids are sisters to the clade-containing eudicots and monocots, which have rapidly diversified among the common ancestors of these three lineages. However, prior to this study, information on the function of the YABBY genes in magnoliids was extremely limited to the third major clades and the early diverging lineage of Mesangiospermae. In this study, the sum of 55 YABBY genes including five genes in INO, six in CRC, eight in YAB2, 22 in YAB5, and 14 in FIL clade were identified from seven magnoliid plants. Sequence analysis showed that all encoded YABBY protein sequences possess the highly conserved YABBY domain and C2C2 zinc-finger domain. Gene and protein structure analysis indicates that a certain number of exons were highly conserved and similar in the same class, and YABBY genes encode proteins of 71-392 amino acids and an open reading frame of 216-1179 bp in magnoliids. Additionally, the predicted molecular weight and isoelectric point of YABBY proteins in three species ranged from 7689.93 to 43578.13 and from 5.33 to 9.87, respectively. Meanwhile, the YABBY gene homolog expression of Litsea was detected at a temporal and spatial level during various developmental stages of leaf and reproductive tissues. This research could provide a brief overview of YABBY gene family evolution and its differential expression in magnoliids. Therefore, this comprehensive diversification analysis would provide a new insight into further understanding of the function of genes in seven magnoliids.

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.

A Critical Cross-Species Comparison of Pollen from Nelumbo nucifera Gaertn. vs. Nymphaea lotus L. for Authentication of Thai Medicinal Herbal Tea

"Bau Luang" or Nelumbo nucifera Gaertn. is an aquatic medicinal herb that has been used as a component of traditional medicines, medicinal products, and herbal tea for good health, particularly in Asia. The stamen of N. nucifera is an important part of this medicinal plant that is used in the form of dried and/or powdered stamens for herbal tea as well as the main ingredient of some traditional remedies. However, there is another aquatic herb called "Bau Sai" or Nymphaea lotus L. that is distributed in similar locations. Living plants of these two aquatic species may be classified according to their morphology, but the dried and powdered stamens of these two medicinal species are difficult to distinguish. The major reason of adulteration is the higher price of Bau Luang stamen. As a result, various methods of authentication, such as pollen micromorphology evaluation using scanning electron microscopy (SEM) analysis, bioinformatics analysis of two nuclear and plastic DNA markers, phytochemical stamen profiling, and Fourier transform infrared (FTIR) analysis of stamen plant material authentication from Bau Luang and Bau Sai, have been used in this present research in order to avoid some adulteration and/or misuse between the dried stamens of Bau Luang and Bau Sai. These results showed that the micro-morphology of pollen (size of pollen grain, number of apertures, and surface ornamentation) from the SEM analysis, some phytochemical compounds and the FTIR sporopollenin-to-protein ratio signal analysis are potential tools for authentication and identification of these two medicinal plants from their dried-stamen materials. This model of investigation may also be used to distinguish dried plant material from other problematic plant groups.

Water lily (Nymphaea thermarum) genome reveals variable genomic signatures of ancient vascular cambium losses

For more than 225 million y, all seed plants were woody trees, shrubs, or vines. Shortly after the origin of angiosperms ∼140 million y ago (MYA), the Nymphaeales (water lilies) became one of the first lineages to deviate from their ancestral, woody habit by losing the vascular cambium, the meristematic population of cells that produces secondary xylem (wood) and phloem. Many of the genes and gene families that regulate differentiation of secondary tissues also regulate the differentiation of primary xylem and phloem, which are produced by apical meristems and retained in nearly all seed plants. Here, we sequenced and assembled a draft genome of the water lily Nymphaea thermarum, an emerging system for the study of early flowering plant evolution, and compared it to genomes from other cambium-bearing and cambium-less lineages (e.g., monocots and Nelumbo). This revealed lineage-specific patterns of gene loss and divergence. Nymphaea is characterized by a significant contraction of the HD-ZIP III transcription factors, specifically loss of REVOLUTA, which influences cambial activity in other angiosperms. We also found the Nymphaea and monocot copies of cambium-associated CLE signaling peptides display unique substitutions at otherwise highly conserved amino acids. Nelumbo displays no obvious divergence in cambium-associated genes. The divergent genomic signatures of convergent loss of vascular cambium reveals that even pleiotropic genes can exhibit unique divergence patterns in association with independent events of trait loss. Our results shed light on the evolution of herbaceousness-one of the key biological innovations associated with the earliest phases of angiosperm evolution.

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.

Evolution of development of pollen performance

With the origin of pollination in ancient seed plants, the male gametophyte ("pollen") began to evolve a new and unique life history stage, the progamic phase, a post-pollination period in which pollen sexual maturation occurs in interaction with sporophyte-derived tissues. Pollen performance traits mediate the timing of the fertilization process, often in competition with other pollen, via the speed of pollen germination, sperm development, and pollen tube growth. Studies of pollen development rarely address the issue of performance or its evolution, which involves linking variation in developmental rates to relative fitness within populations or to adaptations on a macroevolutionary scale. Modifications to the pollen tube pathway and changes in the intensity of pollen competition affect the direction and strength of selection on pollen performance. Hence, pollen developmental evolution is always contextual-it involves both the population biology of pollen reaching stigmas and the co-evolution of sporophytic traits, such as the pollen tube pathway and mating system. For most species, performance evolution generally reflects a wandering history of periods of directional selection and relaxed selection, channeled by developmental limitations, a pattern that favors the accumulation of diversity and redundancy in developmental mechanisms and the genetic machinery. Developmental biologists are focused on finding universal mechanisms that underlie pollen function, and these are largely mechanisms that have evolved through their effects on performance. Here, we suggest ways in which studies of pollen performance or function could progress by cross-fertilization between the "evo" and "devo" fields.

The complete mitochondrial genome of the early flowering plant Nymphaea colorata is highly repetitive with low recombination

BACKGROUND

Mitochondrial genomes of flowering plants (angiosperms) are highly dynamic in genome structure. The mitogenome of the earliest angiosperm Amborella is remarkable in carrying rampant foreign DNAs, in contrast to Liriodendron, the other only known early angiosperm mitogenome that is described as 'fossilized'. The distinctive features observed in the two early flowering plant mitogenomes add to the current confusions of what early flowering plants look like. Expanded sampling would provide more details in understanding the mitogenomic evolution of early angiosperms. Here we report the complete mitochondrial genome of water lily Nymphaea colorata from Nymphaeales, one of the three orders of the earliest angiosperms.

RESULTS

Assembly of data from Pac-Bio long-read sequencing yielded a circular mitochondria chromosome of 617,195 bp with an average depth of 601×. The genome encoded 41 protein coding genes, 20 tRNA and three rRNA genes with 25 group II introns disrupting 10 protein coding genes. Nearly half of the genome is composed of repeated sequences, which contributed substantially to the intron size expansion, making the gross intron length of the Nymphaea mitochondrial genome one of the longest among angiosperms, including an 11.4-Kb intron in cox2, which is the longest organellar intron reported to date in plants. Nevertheless, repeat mediated homologous recombination is unexpectedly low in Nymphaea evidenced by 74 recombined reads detected from ten recombinationally active repeat pairs among 886,982 repeat pairs examined. Extensive gene order changes were detected in the three early angiosperm mitogenomes, i.e. 38 or 44 events of inversions and translocations are needed to reconcile the mitogenome of Nymphaea with Amborella or Liriodendron, respectively. In contrast to Amborella with six genome equivalents of foreign mitochondrial DNA, not a single horizontal gene transfer event was observed in the Nymphaea mitogenome.

CONCLUSIONS

The Nymphaea mitogenome resembles the other available early angiosperm mitogenomes by a similarly rich 64-coding gene set, and many conserved gene clusters, whereas stands out by its highly repetitive nature and resultant remarkable intron expansions. The low recombination level in Nymphaea provides evidence for the predominant master conformation in vivo with a highly substoichiometric set of rearranged molecules.

Molecular breeding of water lily: engineering cold stress tolerance into tropical water lily

Water lilies (order Nymphaeales) are rich in both economic and cultural values. They grow into aquatic herbs, and are divided into two ecological types: tropical and hardy. Although tropical water lilies have more ornamental and medicinal values compared to the hardy water lily, the study and utilization of tropical water lilies in both landscaping and pharmaceutical use is greatly hindered due to their limited planting area. Tropical water lilies cannot survive the winter in areas beyond 24.3°N latitude. Here, the transgenic pipeline through the pollen-tube pathway was generated for water lily for the first time. To improve cold stress tolerance of tropical water lilies, a gene encoding choline oxidase (CodA) driven by a cold stress-inducible promoter was transformed into a tropical water lily through the pollen-tube transformation. Six independent transgenic lines were tested for survival rate during two winter seasons from 2015 to 2017 in Hangzhou (30.3°N latitude). PCR and southern blot detection revealed that the CodA gene had been integrated into the genome. Reverse transcription PCR showed that CodA gene was induced after cold stress treatment, and further quantitative real-time PCR revealed different expressions among six 4 lines and line 3 had the highest expression. Multiple physiological experiments showed that after cold stress treatment, both the conductivity and malondialdehyde (MDA) levels from transgenic plants were significantly lower than those of non-transgenic plants, whereas the content of betaine and the activity of superoxide dismutase, catalase, and peroxidase were higher than those from non-transgenic plants. These results suggest that expression of exogenous CodA gene significantly improved the cold stress tolerance of tropical water lilies through a wide range of physiological alterations. Our results currently expanded a six-latitude cultivating area of the tropical water lilies. These results not only illuminate the bright future for water lily breeding but will also facilitate the functional genomic studies.