The adaptive evolution of Euryale ferox to the aquatic environment through paleo-hexaploidization

Water wisteria genome reveals environmental adaptation and heterophylly regulation in amphibious plants

Heterophylly is a phenomenon whereby an individual plant dramatically changes leaf shape in response to the surroundings. Hygrophila difformis (Acanthaceae; water wisteria), has recently emerged as a model plant to study heterophylly because of its striking leaf shape variation in response to various environmental factors. When submerged, H. difformis often develops complex leaves, but on land it develops simple leaves. Leaf complexity is also influenced by other factors, such as light density, humidity, and temperature. Here, we sequenced and assembled the H. difformis chromosome-level genome (scaffold N50: 60.43 Mb, genome size: 871.92 Mb), which revealed 36 099 predicted protein-coding genes distributed over 15 pseudochromosomes. H. difformis diverged from its relatives during the Oligocene climate-change period and expanded gene families related to its amphibious habit. Genes related to environmental stimuli, leaf development, and other pathways were differentially expressed in submerged and terrestrial conditions, possibly modulating morphological and physiological acclimation to changing environments. We also found that auxin plays a role in H. difformis heterophylly. Finally, we discovered candidate genes that respond to different environmental conditions and elucidated the role of LATE MERISTEM IDENTITY 1 (LMI1) in heterophylly. We established H. difformis as a model for studying interconnections between environmental adaptation and morphogenesis.

Leaf unfolding and lamina biomechanics inSyngonium podophyllumandPilea peperomioides

In nature, leaves and their laminae vary in shape, appearance and unfolding behaviour. We investigated peltate leaves of two model species with peltate leaves and highly different morphology (Syngonium podophyllumandPilea peperomioides) and two distinct unfolding patterns via time-lapse recordings: we observed successive unfolding of leaf halves inS. podophyllumand simultaneous unfolding inP. peperomioides.Furthermore, we gathered relevant morphological and biomechanical data in juvenile (unfolding) and adult (fully unfolded) leaves of both species by measuring the thickness and the tensile modulus of both lamina and veins as a measure of their stiffness. InS. podophyllum, lamina and veins stiffen after unfolding, which may facilitate unfolding in the less stiff juvenile lamina. Secondary venation highly contributes to stiffness in the adult lamina ofS. podophyllum, while the lamina itself withstands tensile loads best in direction parallel to secondary veins. In contrast, the leaf ofP. peperomioideshas a higher lamina thickness and small, non-prominent venation and is equally stiff in every region and direction, although, as is the case inS. podophyllum, thickness and stiffness increase during ontogeny of leaves from juvenile to adult. It could be shown that (changes in) lamina thickness and stiffness can be well correlated with the unfolding processes of both model plants, so that we conclude that functional lamina morphology in juvenile and adult leaf stages and the ontogenetic transition while unfolding is highly dependent on biomechanical characteristics, though other factors are also taken into consideration and discussed.

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.

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

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

Characterization and Coexpression Analysis of the TIFY Family Genes in Euryale ferox Related to Leaf Development

TIFYs are plant-specific transcription factors that contain the TIFY structural domain and play an important role in plant leaf growth and development. However, the role played by TIFY in E. ferox (Euryale ferox Salisb.) leaf development has not been investigated. In this study, 23 TIFY genes were identified in E. ferox. Phylogenetic analyses of the TIFY genes showed clustering into three groups (JAZ, ZIM, and PPD). The TIFY domain was shown to be conserved. JAZ was mainly expanded via wholegenome triplication (WGT) in E. ferox. Based on analyses of the TIFY genes in nine species, we found that JAZ has a closer relationship with PPD, in addition to appearing the most recently and expanding most rapidly, leading to the rapid expansion of TIFYs in Nymphaeaceae. In addition, their different evolution types were discovered. Different gene expressions showed the distinct and corresponsive expression patterns of the EfTIFYs in different stages of tissue and leaf development. Finally, The qPCR analysis revealed that the expression of EfTIFY7.2 and EfTIFY10.1 showed an upward trend and high expression throughout leaf development. Further co-expression analysis indicated that EfTIFY7.2 might be more important for the development of E. ferox leaves. This information will be valuable when exploring the molecular mechanisms of EfTIFYs in plants.

A chromosome-scale genome assembly of Artemisia argyi reveals unbiased subgenome evolution and key contributions of gene duplication to volatile terpenoid diversity

Artemisia argyi Lévl. et Vant., a perennial Artemisia herb with an intense fragrance, is widely used in traditional medicine in China and many other Asian countries. Here, we present a chromosome-scale genome assembly of A. argyi comprising 3.89 Gb assembled into 17 pseudochromosomes. Phylogenetic and comparative genomic analyses revealed that A. argyi underwent a recent lineage-specific whole-genome duplication (WGD) event after divergence from Artemisia annua, resulting in two subgenomes. We deciphered the diploid ancestral genome of A. argyi, and unbiased subgenome evolution was observed. The recent WGD led to a large number of duplicated genes in the A. argyi genome. Expansion of the terpene synthase (TPS) gene family through various types of gene duplication may have greatly contributed to the diversity of volatile terpenoids in A. argyi. In particular, we identified a typical germacrene D synthase gene cluster within the expanded TPS gene family. The entire biosynthetic pathways of germacrenes, (+)-borneol, and (+)-camphor were elucidated in A. argyi. In addition, partial deletion of the amorpha-4,11-diene synthase (ADS) gene and loss of function of ADS homologs may have resulted in the lack of artemisinin production in A. argyi. Our study provides new insights into the genome evolution of Artemisia and lays a foundation for further improvement of the quality of this important medicinal plant.

Molecular identification of phenylalanine ammonia lyase-encoding genes EfPALs and EfPAL2-interacting transcription factors in Euryale ferox

Flavonoids are one of the most important secondary metabolites in plants, and phenylalanine ammonia-lyase (PAL) is the first rate-limiting enzyme for their biosynthesis. However, detailed information on the regulation of PAL in plants is still little. In this study, PAL in E. ferox was identified and functionally analyzed, and its upstream regulatory network was investigated. Through genome-wide identification, we obtained 12 putative PAL genes from E. ferox. Phylogenetic tree and synteny analysis revealed that PAL in E. ferox was expanded and mostly preserved. Subsequently, enzyme activity assays demonstrated that EfPAL1 and EfPAL2 both catalyzed the production of cinnamic acid from phenylalanine only, with EfPAL2 exhibiting a superior enzyme activity. Overexpression of EfPAL1 and EfPAL2 in Arabidopsis thaliana, respectively, both enhanced the biosynthesis of flavonoids. Furthermore, two transcription factors, EfZAT11 and EfHY5, were identified by yeast one-hybrid library assays as binding to the promoter of EfPAL2, and further luciferase (LUC) activity analysis indicated that EfZAT11 promoted the expression of EfPAL2, while EfHY5 repressed the expression of EfPAL2. These results suggested that EfZAT11 and EfHY5 positively and negatively regulate flavonoid biosynthesis, respectively. Subcellular localization revealed that EfZAT11 and EfHY5 were localized in the nucleus. Our findings clarified the key EfPAL1 and EfPAL2 of flavonoid biosynthesis in E. ferox and established the upstream regulatory network of EfPAL2, which would provide novel information for the study of flavonoid biosynthesis mechanism.

Combined analysis of the transcriptome and metabolome provides insights into the fleshy stem expansion mechanism in stem lettuce

As a stem variety of lettuce, the fleshy stem is the main product organ of stem lettuce. The molecular mechanism of fleshy stem expansion in stem lettuce is a complex biological process. In the study, the material accumulation, gene expression, and morphogenesis during fleshy stem expansion process were analyzed by the comparative analysis of metabolome, transcriptome and the anatomical studies. The anatomical studies showed that the occurrence and activity of vascular cambium mainly led to the development of fleshy stems; and the volume of pith cells gradually increased and arranged tightly during the expansion process. A total of 822 differential metabolites and 9,383 differentially expressed genes (DEGs) were identified by the metabolomics and transcriptomics analyses, respectively. These changes significantly enriched in sugar synthesis, glycolysis, and plant hormone anabolism. The expression profiles of genes in the sugar metabolic pathway gradually increased in fleshy stem expansion process. But the sucrose content was the highest in the early stage of fleshy stem expansion, other sugars such as fructose and glucose content increased during fleshy stem expansion process. Plant hormones, including IAA, GA, CTK, and JA, depicted important roles at different stem expansion stages. A total of 1,805 DEGs were identified as transcription factors, such as MYB, bHLH, and bZIP, indicating that these transcription factor families might regulate the fleshy stems expansion in lettuce. In addition, the expression patterns identified by qRT-PCR were consistent with the expression abundance identified by the transcriptome data. The important genes and metabolites identified in the lettuce fleshy stem expansion process will provide important information for the further molecular mechanism study of lettuce fleshy stem growth and development.

Euryale Small Auxin Up RNA62 promotes cell elongation and seed size by altering the distribution of indole-3-acetic acid under the light

Euryale (Euryale ferox Salisb.) is an aquatic crop used as both food and drug in Asia, but its utilization is seriously limited due to low yield. Previously, we hypothesized that Euryale small auxin up RNAs (EuSAURs) regulate seed size, but the underlying biological functions and molecular mechanisms remain unclear. Here, we observed that the hybrid Euryale lines (HL) generate larger seeds with higher indole-3-acetic acid (IAA) concentrations than those in the North Gordon Euryale (WT). Histological analysis suggested that a larger ovary in HL is attributed to longer cells around. Overexpression of EuSAUR62 in rice (Oryza sativa L.) resulted in larger glumes and grains and increased the length of glume cells. Immunofluorescence and protein interaction assays revealed that EuSAUR62 modulates IAA accumulation around the rice ovary by interacting with the rice PIN-FORMED 9, an auxin efflux carrier protein. Euryale basic region/leucine zipper 55 (EubZIP55), which was highly expressed in HL, directly binds to the EuSAUR62 promoter and activated the expression of EuSAUR62. Constant light increased the expression of both EubZIP55 and EuSAUR62 with auxin-mediated hook curvature in HL seedlings. Overall, we proposed that EuSAUR62 is a molecular bridge between light and IAA and plays a crucial role in regulating the size of the Euryale seed.

EfABI4 Transcription Factor Is Involved in the Regulation of Starch Biosynthesis in Euryale ferox Salisb Seeds

Starch is the final product of photosynthesis and the main storage form in plants. Studies have shown that there is a close synergistic regulatory relationship between ABA signal transduction and starch biosynthesis. In this study, we employed RNA sequencing (RNA-Seq) to investigate transcriptomic changes of the Euryale ferox seeds treated by exogenous ABA. The differentially expressed genes engaged in the "Starch and sucrose" and "TCA cycle" pathway. Furthermore, the key transcription factor EfABI4 in ABA signaling pathway and the key genes of starch biosynthesis (EfDBE1, EfSBE2, EfSS1, EfSS2, EfSS3, EfSS4 and EfGBSS1) were significantly up-regulated. Further, the Euryale ferox plant was treated with ABA, it was found that the total starch content of Euryale ferox seeds at different development stages was significantly higher than that of the control, and the key genes of starch synthesis in Euryale ferox seeds were also significantly up-regulated. Finally, yeast one-hybrid and dual luciferase assay proved that EfABI4 can promote the expression of EfSS1 by directly binding to its promoter. Subcellular localization results showed that EfABI4 protein was located at the nucleus and EfSS1 protein was located in the cytomembrane. These findings revealed that ABA promotes starch synthesis and accumulation by mediating EfABI4 to directly promote EfSS1 gene expression, which is helpful for understanding starch synthesis in seeds.