The Spirodela polyrhiza genome reveals insights into its neotenous reduction fast growth and aquatic lifestyle

Genomics: To be (or not to be) a duckweed

Duckweeds are among the smallest and fastest-growing flowering plants. In a new study that combines experimental data with phylogenomic comparisons across the clade, the authors explore how changes in gene content, epigenetic pathways, and their interplay shaped the body plan, aquatic lifestyle, and clonal growth habit of this plant family.

Chromosome-level genome assembly and annotation of Anthurium amnicola

Anthurium amnicola is in the monocot family Araceae, subfamily Pothoideae and is a contributing species in Hawaii floriculture industry hybrids. To support future genetic improvements to this commodity, we sequenced and assembled the A. amnicola genome to chromosome-scale using PacBio HiFi and short-read Hi-C sequencing. A total of 98.51% of the 4.79 Gb genome is anchored into 15 chromosomes, with 99.2% gene completeness and a high LTR assembly index (LAI) score of 21.73, indicative of a complete, high-quality assembly. Annotation reveals the presence of 20,380 protein-coding genes, with 78.52% of the genome composed of repetitive sequences, predominantly long terminal repeat retrotransposons (LTR-RT). Phylogenetic analysis identified evolutionary relationships between A. amnicola and representative species in the Araceae and other plant families. This study provides the first reference genome sequence for the neotropical genus Anthurium and insights into Araceae evolution, benefiting the floriculture industry and evolutionary studies.

Insights into membrane-bound fatty acid desaturase genes in tigernut (Cyperus esculentus L.), an oil-rich tuber plant in Cyperaceae

BACKGROUND

Tigernut (Cyperus esculentus L.), an oil-rich tuber plant of the Cyperaceae family, is typical for the naturally high content of oleic acid. However, to date, genes contributing to oil composition have not been well characterized.

RESULTS

In this study, the first genome-wide analysis of tigernut genes encoding membrane-bound fatty acid desaturases (FADs), the key contributors to oil composition, is presented. According to phylogenetic analysis, ten members identified from the tigernut genome were assigned into seven out of eight evolutionary groups as defined in Arabidopsis thaliana, i.e., FAD2 (3), FAD6 (1), FAD3 (1), FAD7 (1), FAD4 (1), DES (1), and SLD (2). In contrast to the absence of an FAD5 homolog, FAD2 and SLD in tigernut were shown to have expanded via tandem and dispersed duplications, respectively. Comparison of 285 members from 29 representative plant species resulted in 11 orthogroups, where FAD2a, FAD6, FAD7, FAD3, FAD4, FAD5, DES, and SLD1 were shown to have already appeared in the ancestor of seed plants. Significantly, orthologous and syntenic analyses revealed that loss of FAD5 and expansion of SLD in tigernut are lineage-specific, occurred sometime before the radiation of core monocots, in contrast to species-specific expansion of FAD2. Moreover, though no syntenic relationship was observed between CeFAD genes, our comparative genomics analyses indicated that FAD3 and -7 are more likely to arise from segmental duplication. Structural variation and expression divergence of CeFAD genes were also observed. Gain of introns in CeFAD4, CeSLD1, and CeSLD2 was shown to be lineage-specific, occurred sometime before Cyperaceae-Juncaceae split. Tissue-specific expression analysis revealed that CeFAD2-1, CeFAD6, and CeFAD7 were constitutively expressed, whereas others were tissue-specific. Among five paralogs identified, CeFAD2-1 and CeSLD1 have evolved to be two dominant members. Putative roles of CeFAD2-1 in oil accumulation are supported by 1) exhibited an expression pattern positively associated with oil accumulation during tuber development; 2) were expressed more in tubers than their orthologs in C. rotundus. Additionally, in contrast to high expression of CrFAD3, transcript levels of CeFAD3 in tubers were fairly low, which may explain the distinct α-linolenic acid content between these two close species.

CONCLUSIONS

Our findings provide a global view of CeFAD genes, which not only highlights lineage-specific evolution of the family, but also provides valuable information for further functional analysis and genetic improvement in tigernut.

Atypical epigenetic and small RNA control of degenerated transposons and their fragments in clonally reproducing Spirodela polyrhiza

A handful of model plants have provided insight into silencing of transposable elements (TEs) through RNA-directed DNA methylation (RdDM). Guided by 24 nt long small-interfering RNAs (siRNAs), this epigenetic regulation installs DNA methylation and histone modifications like H3K9me2, which can be subsequently maintained independently of siRNAs. However, the genome of the clonally propagating duckweed Spirodela polyrhiza (Lemnaceae) has low levels of DNA methylation, very low expression of RdDM components, and near absence of 24 nt siRNAs. Moreover, some genes encoding RdDM factors, DNA methylation maintenance, and RNA silencing mechanisms are missing from the genome. Here, we investigated the distribution of TEs and their epigenetic marks in the Spirodela genome. Although abundant degenerated TEs have largely lost DNA methylation and H3K9me2 is low, they remain marked by the heterochromatin-associated H3K9me1 and H3K27me1 modifications. In contrast, we find high levels of DNA methylation and H3K9me2 in the relatively few intact TEs, which are source of 24 nt siRNAs, like RdDM-controlled TEs in other angiosperms. The data suggest that, potentially as adaptation to vegetative propagation, RdDM extent, silencing components, and targets are different from other angiosperms, preferentially focused on potentially intact TEs. It also provides evidence for heterochromatin maintenance independently of DNA methylation in flowering plants. These discoveries highlight the diversity of silencing mechanisms that exist in plants and the importance of using disparate model species to discover these mechanisms.

Integrative Analysis Provides Insights into Genes Encoding LEA_5 Domain-Containing Proteins in Tigernut (Cyperus esculentus L.)

LEA_5 domain-containing proteins constitute a small family of late embryogenesis-abundant proteins that are essential for seed desiccation tolerance and dormancy. However, their roles in non-seed storage organs such as underground tubers are largely unknown. This study presents the first genome-scale analysis of the LEA_5 family in tigernut (Cyperus esculentus L.), a Cyperaceae plant producing desiccation-tolerant tubers. Four LEA_5 genes identified from the tigernut genome are twice of two present in model plants Arabidopsis thaliana and Oryza sativa. A comparison of 86 members from 34 representative plant species revealed the monogenic origin and lineage-specific family evolution in Poales, which includes the Cyperaceae family. CeLEA5 genes belong to four out of five orthogroups identified in this study, i.e., LEA5a, LEA5b, LEA5c, and LEA5d. Whereas LEA5e is specific to eudicots, LEA5b and LEA5d appear to be Poales-specific and LEA5c is confined to families Cyperaceae and Juncaceae. Though no syntenic relationship was observed between CeLEA5 genes, comparative genomics analyses indicated that LEA5b and LEA5c are more likely to arise from LEA5a via whole-genome duplication. Additionally, local duplication, especially tandem duplication, also played a role in the family expansion in Juncus effuses, Joinvillea ascendens, and most Poaceae plants examined in this study. Structural variation (e.g., fragment insertion) and expression divergence of LEA_5 genes were also observed. Whereas LEA_5 genes in A. thaliana, O. sativa, and Zea mays were shown to be preferentially expressed in seeds/embryos, CeLEA5 genes have evolved to be predominantly expressed in tubers, exhibiting seed desiccation-like accumulation during tuber maturation. Moreover, CeLEA5 orthologs in C. rotundus showed weak expression in various stages of tuber development, which may explain the difference in tuber desiccation tolerance between these two close species. These findings highlight the lineage-specific evolution of the LEA_5 family, which facilitates further functional analysis and genetic improvement in tigernut and other species.

Insights into the Stearoyl-Acyl Carrier Protein Desaturase (SAD) Family in Tigernut (Cyperus esculentus L.), an Oil-Bearing Tuber Plant

Plant oils rich in oleic acid (OA) are attracting considerable attention for their high nutritional value and significant industrial potential. Stearoyl-acyl carrier protein desaturases (SADs) are a class of soluble desaturases that play a key role in OA accumulation in plants. In this study, the first genome-wide characterization of the SAD gene family was conducted in tigernut (Cyperus esculentus L. var. sativus Baeck., Cyperaceae), an oil-rich tuber plant typical for its high OA content. Six SAD genes identified from the tigernut genome are comparative to seven reported in two model plants Arabidopsis thaliana and Oryza sativa, but relatively more than four were found in most Cyperaceae species examined in this study. A comparison of 161 SAD genes from 29 representative plant species reveals the monogenic origin and lineage-specific family evolution in Poales. C. esculentus SAD genes (CeSADs) were shown to constitute two evolutionary groups (i.e., FAB2 and AAD) and four out of 12 orthogroups identified in this study, i.e., FAB2a, FAB2b, FAB2c, and AAD1. Whereas FAB2a and AAD1 are widely distributed, FAB2b and FAB2c are specific to Cyperaceae, which may arise from FAB2a via tandem and dispersed duplications, respectively. Though FAB2d and AAD2 are also broadly present in monocots, they are more likely to be lost in the Cyperaceae ancestor sometime after the split with its close family, Juncaceae. In tigernut, FAB2a appears to have undergone species-specific expansion via tandem duplication. Frequent structural variation and apparent expression divergence were also observed. Though FAB2a and AAD1 usually feature two and one intron, respectively, gain of certain introns was observed in CeSAD genes, all of which have three introns. Despite recent expansion of the FAB2 group, CeFAB2-1 has evolved into the dominant member that was highly and constitutively expressed in all tested organs. Moreover, CeFAB2-1, CeAAD1, as well as CeFAB2-5 have evolved to be predominantly expressed in tubers and thus contribute to high OA accumulation. These findings highlight lineage-specific evolution of the SAD family and putative roles of CeSAD genes in tuber oil accumulation, which facilitate further functional analysis and genetic improvement in tigernut and other species.

Duckweed: a starch-hyperaccumulating plant under cultivation with a combination of nutrient limitation and elevated CO2

Introduction

The increasing global demand for starch has created an urgent need to identify more efficient and sustainable production methods. However, traditional starch sources, such as crop-based options, experience significant bottlenecks due to limitations in land use, water consumption, and the impacts of climate change. Therefore, there is a pressing need to explore and develop new sources of starch.

Methods

We develop a novel duckweed cultivation technology that combines nutrients limitation and CO2 supplementation to achieve very high starch content. In this study, we integrated whole-genome sequencing, epigenomics, transcriptomics, enzyme activity, and composition variation to elucidate the mechanisms of efficient starch accumulation in duckweed in terms of starch accumulation and carbon partitioning, regulation of the expression of genes in the starch metabolic pathway, and sucrose biosynthesis and transportation.

Results and discussion

Although Landoltia punctata exhibits dramatic gene family contraction, its starch content and productivity reached 72.2% (dry basis) and 10.4 g m-2 d-1, respectively, in 10 days, equivalent to a yield of 38.0 t ha-1 y-1, under nutrient limitation treatment with elevated CO2 levels. We also examined the mechanism of high starch accumulation in duckweed. This phenomenon is associated with the regulation of DNA methylation and transcription factors as well as the significantly upregulated transcription levels and the increased activities of key enzymes involved in starch biosynthesis. Moreover, while nitrogen redistribution was increased, sucrose biosynthesis and transportation and lignocellulose biosynthesis were reduced. These alterations led to a reduction in lignocellulose and protein contents and ultimately an increase in the accumulation of starch in the chloroplasts.

Conclusion

This work demonstrates the potential of duckweed as a highly efficient starch producer.

Insights into Genes Encoding LEA_1 Domain-Containing Proteins in Cyperus esculentus, a Desiccation-Tolerant Tuber Plant

LEA_1 domain-containing proteins constitute a class of late-embryogenesis-abundant proteins that are highly hydrophilic and predominantly accumulate in mature seeds. Though LEA_1 proteins have been proven to be essential for seed desiccation tolerance and longevity, little information is available on their roles in non-seed storage organs. In this study, a first genome-wide characterization of the LEA_1 gene family was conducted in tigernut (Cyperus esculentus L., Cyperaceae), whose underground tubers are desiccation tolerant with a moisture content of less than 6%. Five family members identified in tigernut are comparative to four to six found in seven other Cyperaceae plants, but relatively more than three reported in Arabidopsis. Further comparison of 125 members from 29 plant species supports early divergence of the LEA_1 family into two phylogenetic groups before angiosperm radiation, and gene expansion in tigernut was contributed by whole-genome duplications occurring after the split with the eudicot clade. These two phylogenetic groups could be further divided into six orthogroups in the momocot clade, five of which are present in tigernut and the remaining one is Poaceae specific. Frequent structural variation and expression divergence of paralogs were also observed. Significantly, in contrast to seed-preferential expression of LEA_1 genes in Arabidopsis, rice, and maize, transcriptional profiling and qRT-PCR analysis revealed that CeLEA1 genes have evolved to predominantly express in tubers, exhibiting a seed desiccation-like accumulation during tuber development. Moreover, CeLEA1 transcripts in tubers were shown to be considerably more than that of their orthologs in purple nutsedge, another Cyperaceae plant producing desiccation-sensitive tubers. These results imply species-specific activation and key roles of CeLEA1 genes in the acquisition of desiccation tolerance of tigernut tubers as observed in orthodox seeds. Our findings not only improve the understanding of lineage-specific evolution of the LEA_1 family, but also provide valuable information for further functional analysis and genetic improvement in tigernut.

Enzymatic Characterization of SpPAL Genes in S. polyrhiza and Overexpression of the SpPAL3

Phenylalanine ammonia-lyase (PAL, EC 4.3.1.5) catalyzes the deamination of phenylalanine, which is the initial step in the biosynthesis of phenylpropanoids. It serves as a crucial enzyme that facilitates the transfer of carbon from primary to secondary metabolism in plants. Duckweed is regarded as a promising chassis plant in synthetic biology research and application, due to its being rich in secondary metabolites and other advantages. The genes encoding PAL in Spirodela polyrhiza (L.) Schleid, the giant duckweed, were investigated in this study. Three SpPAL genes (SpPAL1-SpPAL3) were identified and cloned. All of them were successfully expressed in E. coli, and their recombinant proteins all showed PAL activity. In addition, SpPAL1 and SpPAL2 proteins could also utilize tyrosine as substrate, although the activity was low. A qRT-PCR analysis demonstrated that the expression of SpPAL3 was most pronounced in young fronds. It was found that the expression of SpPAL1 and SpPAL3 was significantly induced by MeJA treatment. Overexpression of SpPAL3 in Lemna turionifera inhibited the growth of fronds and adventitious roots in the transgenic plants, indicating the importance of SpPAL3 in duckweed besides its involvement in the secondary metabolism.

Phytohormonal regulation determines the organization pattern of shoot aerenchyma in greater duckweed (Spirodela polyrhiza)

Airspace or aerenchyma is crucial for plant development and acclimation to stresses such as hypoxia, drought, and nutritional deficiency. Although ethylene-mediated signaling cascades are known to regulate aerenchyma formation in stems and roots under hypoxic conditions, the precise mechanisms remain unclear. Moreover, the cellular dynamics underlying airspace formation in shoots are poorly understood. We investigated the stage-dependent structural dynamics of shoot aerenchyma in greater duckweed (Spirodela polyrhiza), a fast-growing aquatic herb with well-developed aerenchyma in its floating fronds. Using X-ray micro-computed tomography and histological analysis, we showed that the spatial framework of aerenchyma is established before frond volume increases, driven by cell division and expansion. The substomatal cavity connecting aerenchyma to stomata formed via programmed cell death (PCD) and was closely associated with guard cell development. Additionally, transcriptome analysis and pharmacological studies revealed that the organization of aerenchyma in greater duckweed is determined by the interplay between PCD and proliferation. This balance is governed by spatiotemporal regulation of phytohormone signaling involving ethylene, abscisic acid, and salicylic acid. Overall, our study reveals the structural dynamics and phytohormonal regulation underlying aerenchyma development in duckweed, improving our understanding of how plants establish distinct architectural arrangements. These insights hold the potential for wide-ranging application, not only in comprehending aerenchyma formation across various plant species but also in understanding how airspaces are formed within the leaves of terrestrial plants.

Covering the bases: Population genomic structure of Lemna minor and the cryptic species L. japonica in Switzerland

Duckweeds, including the common duckweed Lemna minor, are increasingly used to test eco-evolutionary theories. Yet, despite its popularity and near-global distribution, the understanding of its population structure (and genetic variation therein) is still limited. It is essential that this is resolved, because of the impact genetic diversity has on experimental responses and scientific understanding. Through whole-genome sequencing, we assessed the genetic diversity and population genomic structure of 23 natural Lemna spp. populations from their natural range in Switzerland. We used two distinct analytical approaches, a reference-free kmer approach and the classical reference-based one. Two genetic clusters were identified across the described species distribution of L. minor, surprisingly corresponding to species-level divisions. The first cluster contained the targeted L. minor individuals and the second contained individuals from a cryptic species: Lemna japonica. Within the L. minor cluster, we identified a well-defined population structure with little intra-population genetic diversity (i.e., within ponds) but high inter-population diversity (i.e., between ponds). In L. japonica, the population structure was significantly weaker and genetic variation between a subset of populations was as low as within populations. This study revealed that L. japonica is more widespread than previously thought. Our findings signify that thorough genotype-to-phenotype analyses are needed in duckweed experimental ecology and evolution.

Duckweed: Beyond an Efficient Plant Model System

Duckweed (Lemnaceae) rises as a crucial model system due to its unique characteristics and wide-ranging utility. The significance of physiological research and phytoremediation highlights the intricate potential of duckweed in the current era of plant biology. Special attention to duckweed has been brought due to its distinctive features of nutrient uptake, ion transport dynamics, detoxification, intricate signaling, and stress tolerance. In addition, duckweed can alleviate environmental pollutants and enhance sustainability by participating in bioremediation processes and wastewater treatment. Furthermore, insights into the genomic complexity of Lemnaceae species and the flourishing field of transgenic development highlight the opportunities for genetic manipulation and biotechnological innovations. Novel methods for the germplasm conservation of duckweed can be adopted to preserve genetic diversity for future research endeavors and breeding programs. This review centers around prospects in duckweed research promoting interdisciplinary collaborations and technological advancements to drive its full potential as a model organism.

Population genomics and epigenomics of Spirodela polyrhiza provide insights into the evolution of facultative asexuality

Many plants are facultatively asexual, balancing short-term benefits with long-term costs of asexuality. During range expansion, natural selection likely influences the genetic controls of asexuality in these organisms. However, evidence of natural selection driving asexuality is limited, and the evolutionary consequences of asexuality on the genomic and epigenomic diversity remain controversial. We analyzed population genomes and epigenomes of Spirodela polyrhiza, (L.) Schleid., a facultatively asexual plant that flowers rarely, revealing remarkably low genomic diversity and DNA methylation levels. Within species, demographic history and the frequency of asexual reproduction jointly determined intra-specific variations of genomic diversity and DNA methylation levels. Genome-wide scans revealed that genes associated with stress adaptations, flowering and embryogenesis were under positive selection. These data are consistent with the hypothesize that natural selection can shape the evolution of asexuality during habitat expansions, which alters genomic and epigenomic diversity levels.

Mechanisms of Spirodela polyrhiza tolerance to FGD wastewater-induced heavy-metal stress: Lipidomics, transcriptomics, and functional validation

Unlike terrestrial angiosperm plants, the freshwater aquatic angiosperm duckweed (Spirodela polyrhiza) grows directly in water and has distinct responses to heavy-metal stress. Plantlets accumulate metabolites, including lipids and carbohydrates, under heavy-metal stress, but how they balance metabolite levels is unclear, and the gene networks that mediate heavy-metal stress responses remain unknown. Here, we show that heavy-metal stress induced by flue gas desulfurization (FGD) wastewater reduces chlorophyll contents, inhibits growth, reduces membrane lipid biosynthesis, and stimulates membrane lipid degradation in S. polyrhiza, leading to triacylglycerol and carbohydrate accumulation. In FGD wastewater-treated plantlets, the degraded products of monogalactosyldiacylglycerol, primarily polyunsaturated fatty acids (18:3), were incorporated into triacylglycerols. Genes involved in early fatty acid biosynthesis, β-oxidation, and lipid degradation were upregulated while genes involved in cuticular wax biosynthesis were downregulated by treatment. The transcription factor gene WRINKLED3 (SpWRI3) was upregulated in FGD wastewater-treated plantlets, and its ectopic expression increased tolerance to FGD wastewater in transgenic Arabidopsis (Arabidopsis thaliana). Transgenic Arabidopsis plants showed enhanced glutathione and lower malondialdehyde contents under stress, suggesting that SpWRI3 functions in S. polyrhiza tolerance of FGD wastewater-induced heavy-metal stress. These results provide a basis for improving heavy metal-stress tolerance in plants for industrial applications.

Contrasting patterns of 5S rDNA repeats in European and Asian ecotypes of greater duckweed, Spirodela polyrhiza (Lemnaceae)

Ribosomal DNA (rDNA) contains highly conserved, specifically organized sequences encoding ribosomal RNAs (rRNAs) separated by variable non-transcribed intergenic spacers (NTSs) and is abundant in eukaryotic genomes. These characteristics make the rDNA an informative molecular target to study genome organization, molecular evolution, and phylogenetics. In this study, we characterized the 5S rDNA repeats in the greater duckweed Spiroldela polyrhiza, a species known for its small size, rapid growth, highly conserved genome organization, and low mutation rate. Sequence analysis of at least 12 individually cloned PCR fragments containing the 5S rDNA units for each of six ecotypes that originated from Europe (Ukraine) and Asia (China) revealed two distinct types of 5S rDNA repeats containing NTSs of different lengths and nucleotide compositions. The shorter 5S rDNA repeat units had a highly homogeneous 400-bp NTS, with few ecotype- or region-specific single-nucleotide polymorphisms (SNPs). The longer 5S rDNA units had NTSs of 1056-1084 bp with characteristic intra- and inter-genomic variants due to specific SNPs and insertions/deletions of 4-15-bp DNA elements. We also detected significant variability in the ratio of short/long 5S rDNA variants between ecotypes of S. polyrhiza. The contrasting dynamics of the two types of 5S rDNA units, combined with the unusually low repeat copy number (for plants) in S. polyrhiza (46-220 copies per genome), shows that this species could serve as an excellent model for examining the mechanisms of concerted evolution and functional significance of rDNA variability.

TIP aquaporins in Cyperus esculentus: genome-wide identification, expression profiles, subcellular localizations, and interaction patterns

BACKGROUND

Tonoplast intrinsic proteins (TIPs), which typically mediate water transport across vacuolar membranes, play an essential role in plant growth, development, and stress responses. However, their characterization in tigernut (Cyperus esculentus L.), an oil-bearing tuber plant of the Cyperaceae family, is still in the infancy.

RESULTS

In this study, a first genome-wide characterization of the TIP subfamily was conducted in tigernut, resulting in ten members representing five previously defined phylogenetic groups, i.e., TIP1-5. Although the gene amounts are equal to that present in two model plants Arabidopsis and rice, the group composition and/or evolution pattern were shown to be different. Except for CeTIP1;3 that has no counterpart in both Arabidopsis and rice, complex orthologous relationships of 1:1, 1:2, 1:3, 2:1, and 2:2 were observed. Expansion of the CeTIP subfamily was contributed by whole-genome duplication (WGD), transposed, and dispersed duplications. In contrast to the recent WGD-derivation of CeTIP3;1/-3;2, synteny analyses indicated that TIP4 and - 5 are old WGD repeats of TIP2, appearing sometime before monocot-eudicot divergence. Expression analysis revealed that CeTIP genes exhibit diverse expression profiles and are subjected to developmental and diurnal fluctuation regulation. Moreover, when transiently overexpressed in tobacco leaves, CeTIP1;1 was shown to locate in the vacuolar membrane and function in homo/heteromultimer, whereas CeTIP2;1 is located in the cell membrane and only function in heteromultimer. Interestingly, CeTIP1;1 could mediate the tonoplast-localization of CeTIP2;1 via protein interaction, implying complex regulatory patterns.

CONCLUSIONS

Our findings provide a global view of CeTIP genes, which provide valuable information for further functional analysis and genetic improvement through manipulating key members in tigernut.

The unusual predominance of maintenance DNA methylation in Spirodela polyrhiza

Duckweeds are among the fastest reproducing plants, able to clonally divide at exponential rates. However, the genetic and epigenetic impact of clonality on plant genomes is poorly understood. 5-methylcytosine (5mC) is a modified base often described as necessary for the proper regulation of certain genes and transposons and for the maintenance of genome integrity in plants. However, the extent of this dogma is limited by the current phylogenetic sampling of land plant species diversity. Here we analyzed DNA methylomes, small RNAs, mRNA-seq, and H3K9me2 histone modification for Spirodela polyrhiza. S. polyrhiza has lost highly conserved genes involved in de novo methylation of DNA at sites often associated with repetitive DNA, and within genes, however, symmetrical DNA methylation and heterochromatin are maintained during cell division at certain transposons and repeats. Consequently, small RNAs that normally guide methylation to silence repetitive DNA like retrotransposons are diminished. Despite the loss of a highly conserved methylation pathway, and the reduction of small RNAs that normally target repetitive DNA, transposons have not proliferated in the genome, perhaps due in part to the rapid, clonal growth lifestyle of duckweeds.

Genomic analysis of Nypa fruticans elucidates its intertidal adaptations and early palm evolution

Nypa fruticans (Wurmb), a mangrove palm species with origins dating back to the Late Cretaceous period, is a unique species for investigating long-term adaptation strategies to intertidal environments and the early evolution of palms. Here, we present a chromosome-level genome sequence and assembly for N. fruticans. We integrated the genomes of N. fruticans and other palm family members for a comparative genomic analysis, which confirmed that the common ancestor of all palms experienced a whole-genome duplication event around 89 million years ago, shaping the distinctive characteristics observed in this clade. We also inferred a low mutation rate for the N. fruticans genome, which underwent strong purifying selection and evolved slowly, thus contributing to its stability over a long evolutionary period. Moreover, ancient duplicates were preferentially retained, with critical genes having experienced positive selection, enhancing waterlogging tolerance in N. fruticans. Furthermore, we discovered that the pseudogenization of Early Methionine-labelled 1 (EM1) and EM6 in N. fruticans underly its crypto-vivipary characteristics, reflecting its intertidal adaptation. Our study provides valuable genomic insights into the evolutionary history, genome stability, and adaptive evolution of the mangrove palm. Our results also shed light on the long-term adaptation of this species and contribute to our understanding of the evolutionary dynamics in the palm family.

Genetic Mechanism of Non-Targeted-Site Resistance to Diquat in Spirodela polyrhiza

Understanding non-target-site resistance (NTSR) to herbicides represents a pressing challenge as NTSR is widespread in many weeds. Using giant duckweed (Spirodela polyrhiza) as a model, we systematically investigated genetic and molecular mechanisms of diquat resistance, which can only be achieved via NTSR. Quantifying the diquat resistance of 138 genotypes, we revealed an 8.5-fold difference in resistance levels between the most resistant and most susceptible genotypes. Further experiments suggested that diquat uptake and antioxidant-related processes jointly contributed to diquat resistance in S. polyrhiza. Using a genome-wide association approach, we identified several candidate genes, including a homolog of dienelactone hydrolase, that are associated with diquat resistance in S. polyrhiza. Together, these results provide new insights into the mechanisms and evolution of NTSR in plants.

The "Duckweed Dip": Aquatic Spirodela polyrhiza Plants Can Efficiently Uptake Dissolved, DNA-Wrapped Carbon Nanotubes from Their Environment for Transient Gene Expression

Duckweeds (Lemnaceae) are aquatic nongrass monocots that are the smallest and fastest-growing flowering plants in the world. While having simplified morphologies, relatively small genomes, and many other ideal traits for emerging applications in plant biotechnology, duckweeds have been largely overlooked in this era of synthetic biology. Here, we report that Greater Duckweed (Spirodela polyrhiza), when simply incubated in a solution containing plasmid-wrapped carbon nanotubes (DNA-CNTs), can directly uptake the DNA-CNTs from their growth media with high efficiency and that transgenes encoded within the plasmids are expressed by the plants─without the usual need for large doses of nanomaterials or agrobacterium to be directly infiltrated into plant tissue. This process, called the "duckweed dip", represents a streamlined, "hands-off" tool for transgene delivery to a higher plant that we expect will enhance the throughput of duckweed engineering and help to realize duckweed's potential as a powerhouse for plant synthetic biology.

Transcriptomic and physiological analysis of the response of Spirodela polyrrhiza to sodium nitroprusside

BACKGROUND

Spirodela polyrrhiza is a simple floating aquatic plant with great potential in synthetic biology. Sodium nitroprusside (SNP) stimulates plant development and increases the biomass and flavonoid content in some plants. However, the molecular mechanism of SNP action is still unclear.

RESULTS

To determine the effect of SNP on growth and metabolic flux in S. polyrrhiza, the plants were treated with different concentrations of SNP. Our results showed an inhibition of growth, an increase in starch, soluble protein, and flavonoid contents, and enhanced antioxidant enzyme activity in plants after 0.025 mM SNP treatment. Differentially expressed transcripts were analysed in S. polyrrhiza after 0.025 mM SNP treatment. A total of 2776 differentially expressed genes (1425 upregulated and 1351 downregulated) were identified. The expression of some genes related to flavonoid biosynthesis and NO biosynthesis was upregulated, while the expression of some photosynthesis-related genes was downregulated. Moreover, SNP stress also significantly influenced the expression of transcription factors (TFs), such as ERF, BHLH, NAC, and WRKY TFs.

CONCLUSIONS

Taken together, these findings provide novel insights into the mechanisms of underlying the SNP stress response in S. polyrrhiza and show that the metabolic flux of fixed CO2 is redirected into the starch synthesis and flavonoid biosynthesis pathways after SNP treatment.

Induced responses contribute to rapid adaptation of Spirodela polyrhiza to herbivory by Lymnaea stagnalis

Herbivory-induced responses in plants are typical examples of phenotypic plasticity, and their evolution is thought to be driven by herbivory. However, direct evidence of the role of induced responses in plant adaptive evolution to herbivores is scarce. Here, we experimentally evolve populations of an aquatic plant (Spirodela polyrhiza, giant duckweed) and its native herbivore (Lymnaea stagnalis, freshwater snail), testing whether herbivory drives rapid adaptive evolution in plant populations using a combination of bioassays, pool-sequencing, metabolite analyses, and amplicon metagenomics. We show that snail herbivory drove rapid phenotypic changes, increased herbivory resistance, and altered genotype frequencies in the plant populations. Additional bioassays suggest that evolutionary changes of induced responses contributed to the rapid increase of plant resistance to herbivory. This study provides direct evidence that herbivory-induced responses in plants can be subjected to selection and have an adaptive role by increasing resistance to herbivores.

Evidence for the role of sound on the growth and signal response in duckweed

Sound vibration, an external mechanical force, has been proven to modulate plant growth and development like rain, wind, and vibration. However, the role of sound on plants, especially on signal response, has been usually neglected in research. Herein, we investigated the growth state, gene expression, and signal response in duckweed treated with soft music. The protein content in duckweed after music treatment for 7 days was about 1.6 times that in duckweed without music treatment. Additionally, the potential maximum photochemical efficiency of photosystem II (Fv/Fm) ratio in duckweed treated with music was 0.78, which was significantly higher in comparison with the control group (P < .01). Interestingly, music promoted the Glu and Ca signaling response. To further explore the global molecular mechanism, we performed transcriptome analysis and the library preparations were sequenced on an Illumina Hiseq platform. A total of 1296 differentially expressed genes (DEGs) were found for all these investigated genes in duckweed treated with music compared to the control group. Among these, up-regulation of the expression of metabolism-related genes related to glycolysis, cell wall biosynthesis, oxidative phosphorylation, and pentose phosphate pathways were found. Overall, these results provided a molecular basis to music-triggered signal response, transcriptomic, and growth changes in duckweed, which also highlighted the potential of music as an environmentally friendly stimulus to promote improved protein production in duckweed.

Characteristics of Amorphophallus konjac as indicated by its genome

Amorphophallus konjac, belonging to the genus Amorphophallus of the Araceae family, is an economically important crop widely used in health products and biomaterials. In the present work, we performed the whole-genome assembly of A. konjac based on the NovaSeq platform sequence data. The final genome assembly was 4.58 Gb with a scaffold N50 of 3212 bp. The genome includes 39,421 protein-coding genes, and 71.75% of the assemblies were repetitive sequences. Comparative genomic analysis showed 1647 gene families have expanded and 2685 contracted in the A. konjac genome. Likewise, genome evolution analysis indicated that A. konjac underwent whole-genome duplication, possibly contributing to the expansion of certain gene families. Furthermore, we identified many candidate genes involved in the tuber formation and development, cellulose and lignification synthesis. The genome of A. konjac obtained in this work provides a valuable resource for the further study of the genetics, genomics, and breeding of this economically important crop, as well as for evolutionary studies of Araceae family.

Mechanisms of metabolic adaptation in the duckweed Lemna gibba: an integrated metabolic, transcriptomic and flux analysis

BACKGROUND

Duckweeds are small, rapidly growing aquatic flowering plants. Due to their ability for biomass production at high rates they represent promising candidates for biofuel feedstocks. Duckweeds are also excellent model organisms because they can be maintained in well-defined liquid media, usually reproduce asexually, and because genomic resources are becoming increasingly available. To demonstrate the utility of duckweed for integrated metabolic studies, we examined the metabolic adaptation of growing Lemna gibba cultures to different nutritional conditions.

RESULTS

To establish a framework for quantitative metabolic research in duckweeds we derived a central carbon metabolism network model of Lemna gibba based on its draft genome. Lemna gibba fronds were grown with nitrate or glutamine as nitrogen source. The two conditions were compared by quantification of growth kinetics, metabolite levels, transcript abundance, as well as by 13C-metabolic flux analysis. While growing with glutamine, the fronds grew 1.4 times faster and accumulated more protein and less cell wall components compared to plants grown on nitrate. Characterization of photomixotrophic growth by 13C-metabolic flux analysis showed that, under both metabolic growth conditions, the Calvin-Benson-Bassham cycle and the oxidative pentose-phosphate pathway are highly active, creating a futile cycle with net ATP consumption. Depending on the nitrogen source, substantial reorganization of fluxes around the tricarboxylic acid cycle took place, leading to differential formation of the biosynthetic precursors of the Asp and Gln families of proteinogenic amino acids. Despite the substantial reorganization of fluxes around the tricarboxylic acid cycle, flux changes could largely not be associated with changes in transcripts.

CONCLUSIONS

Through integrated analysis of growth rate, biomass composition, metabolite levels, and metabolic flux, we show that Lemna gibba is an excellent system for quantitative metabolic studies in plants. Our study showed that Lemna gibba adjusts to different nitrogen sources by reorganizing central metabolism. The observed disconnect between gene expression regulation and metabolism underscores the importance of metabolic flux analysis as a tool in such studies.

Lemnaceae as Novel Crop Candidates for CO2 Sequestration and Additional Applications

Atmospheric carbon dioxide (CO2) is projected to be twice as high as the pre-industrial level by 2050. This review briefly highlights key responses of terrestrial plants to elevated CO2 and compares these with the responses of aquatic floating plants of the family Lemnaceae (duckweeds). Duckweeds are efficient at removing CO2 from the atmosphere, which we discuss in the context of their exceptionally high growth rates and capacity for starch storage in green tissue. In contrast to cultivation of terrestrial crops, duckweeds do not contribute to CO2 release from soils. We briefly review how this potential for contributions to stabilizing atmospheric CO2 levels is paired with multiple additional applications and services of duckweeds. These additional roles include wastewater phytoremediation, feedstock for biofuel production, and superior nutritional quality (for humans and livestock), while requiring minimal space and input of light and fertilizer. We, furthermore, elaborate on other environmental factors, such as nutrient availability, light supply, and the presence of a microbiome, that impact the response of duckweed to elevated CO2. Under a combination of elevated CO2 with low nutrient availability and moderate light supply, duckweeds' microbiome helps maintain CO2 sequestration and relative growth rate. When incident light intensity increases (in the presence of elevated CO2), the microbiome minimizes negative feedback on photosynthesis from increased sugar accumulation. In addition, duckweed shows a clear propensity for absorption of ammonium over nitrate, accepting ammonium from their endogenous N2-fixing Rhizobium symbionts, and production of large amounts of vegetative storage protein. Finally, cultivation of duckweed could be further optimized using hydroponic vertical farms where nutrients and water are recirculated, saving both resources, space, and energy to produce high-value products.

The "Duckweed Dip": Aquatic Spirodela polyrhiza Plants Can Efficiently Uptake Dissolved, DNA-Wrapped Carbon Nanotubes from Their Environment for Transient Gene Expression

Duckweeds (Lemnaceae) are aquatic non-grass monocots that are the smallest and fastest-growing flowering plants in the world. While having simplified morphologies, relatively small genomes, and many other ideal traits for emerging applications in plant biotechnology, duckweeds have been largely overlooked in this era of synthetic biology. Here, we report that Greater Duckweed (Spirodela polyrhiza), when simply incubated in a solution containing plasmid-wrapped carbon nanotubes (DNA-CNTs), can directly up-take the DNA-CNTs from their growth media with high efficiency and that transgenes encoded within the plasmids are expressed by the plants-without the usual need for large doses of nanomaterials or agrobacterium to be directly infiltrated into plant tissue. This process, called the "duckweed dip", represents a streamlined, 'hands-off' tool for transgene delivery to a higher plant that we expect will enhance the throughput of duckweed engineering and help to realize duckweed's potential as a powerhouse for plant synthetic biology. (148 words).

Photosynthesis in the Biomass Model Species Lemna minor Displays Plant-Conserved and Species-Specific Features

Lemnaceae are small freshwater plants with extraordinary high growth rates. We aimed to test whether this correlates with a more efficient photosynthesis, the primary energy source for growth. To this end, we compared photosynthesis properties of the duckweed Lemna minor and the terrestrial model plant Arabidopsis thaliana. Chlorophyll fluorescence analyses revealed high similarity in principle photosynthesis characteristics; however, Lemna exhibited a more effective light energy transfer into photochemistry and more stable photosynthesis parameters especially under high light intensities. Western immunoblot analyses of representative photosynthesis proteins suggested potential post-translational modifications in Lemna proteins that are possibly connected to this. Phospho-threonine phosphorylation patterns of thylakoid membrane proteins displayed a few differences between the two species. However, phosphorylation-dependent processes in Lemna such as photosystem II antenna association and the recovery from high-light-induced photoinhibition were not different from responses known from terrestrial plants. We thus hypothesize that molecular differences in Lemna photosynthesis proteins are associated with yet unidentified mechanisms that improve photosynthesis and growth efficiencies. We also developed a high-magnification video imaging approach for Lemna multiplication which is useful to assess the impact of external factors on Lemna photosynthesis and growth.

Sixth International Conference on Duckweed Research and Applications Presents Lemnaceae as a Model Plant System in the Genomics and Postgenomics Era

The 6th International Conference on Duckweed Research and Applications (6th ICDRA) was organized at the Institute of Plant Genetics and Crop Plant Research (IPK) located in Gatersleben, Germany, from 29 May to 1 June 2022. The growing community of duckweed research and application specialists was noted with participants from 21 different countries including an increased share of newly integrated young researchers. The four-day conference focused on diverse aspects of basic and applied research together with practical applications of these tiny aquatic plants that could have an enormous potential for biomass production.

Genomics of turions from the Greater Duckweed reveal its pathways for dormancy and re-emergence strategy

Over 15 families of aquatic plants are known to use a strategy of developmental switching upon environmental stress to produce dormant propagules called turions. However, few molecular details for turion biology have been elucidated due to the difficulties in isolating high-quality nucleic acids from this tissue. We successfully developed a new protocol to isolate high-quality transcripts and carried out RNA-seq analysis of mature turions from the Greater Duckweed Spirodela polyrhiza. Comparison of turion transcriptomes to that of fronds, the actively growing leaf-like tissue, were carried out. Bioinformatic analysis of high confidence, differentially expressed transcripts between frond and mature turion tissues revealed major pathways related to stress tolerance, starch and lipid metabolism, and dormancy that are mobilized to reprogram frond meristems for turion differentiation. We identified the key genes that are likely to drive starch and lipid accumulation during turion formation, as well as those in pathways for starch and lipid utilization upon turion germination. Comparison of genome-wide cytosine methylation levels also revealed evidence for epigenetic changes in the formation of turion tissues. Similarities between turions and seeds provide evidence that key regulators for seed maturation and germination were retooled for their function in turion biology.

Metabolic flexibility during a trophic transition reveals the phenotypic plasticity of greater duckweed (Spirodela polyrhiza 7498)

The greater duckweed (Spirodela polyrhiza 7498) exhibits trophic diversity (photoautotrophic, heterotrophic, photoheterotrophic, and mixotrophic growth) depending on the availability of exogenous organic carbon sources and light. Here, we show that the ability to transition between various trophic growth conditions is an advantageous trait, providing great phenotypic plasticity and metabolic flexibility in S. polyrhiza 7498. By comparing S. polyrhiza 7498 growth characteristics, metabolic acclimation, and cellular ultrastructure across these trophic modes, we show that mixotrophy decreases photosynthetic performance and relieves the CO2 limitation of photosynthesis by enhancing the CO2 supply through the active respiration pathway. Proteomic and metabolomic analyses corroborated that S. polyrhiza 7498 increases its intracellular CO2 and decreases reactive oxygen species under mixotrophic and heterotrophic conditions, which substantially suppressed the wasteful photorespiration and oxidative-damage pathways. As a consequence, mixotrophy resulted in a higher biomass yield than the sum of photoautotrophy and heterotrophy. Our work provides a basis for using trophic transitions in S. polyrhiza 7498 for the enhanced accumulation of value-added products.

DNA methylation signatures of duplicate gene evolution in angiosperms

Gene duplication is a source of evolutionary novelty. DNA methylation may play a role in the evolution of duplicate genes (paralogs) through its association with gene expression. While this relationship has been examined to varying extents in a few individual species, the generalizability of these results at either a broad phylogenetic scale with species of differing duplication histories or across a population remains unknown. We applied a comparative epigenomics approach to 43 angiosperm species across the phylogeny and a population of 928 Arabidopsis (Arabidopsis thaliana) accessions, examining the association of DNA methylation with paralog evolution. Genic DNA methylation was differentially associated with duplication type, the age of duplication, sequence evolution, and gene expression. Whole genome duplicates were typically enriched for CG-only gene-body methylated or unmethylated genes, while single-gene duplications were typically enriched for non-CG methylated or unmethylated genes. Non-CG methylation, in particular, was characteristic of more recent single-gene duplicates. Core angiosperm gene families differentiated into those which preferentially retain paralogs and 'duplication-resistant' families, which convergently reverted to singletons following duplication. Duplication-resistant families that still have paralogous copies were, uncharacteristically for core angiosperm genes, enriched for non-CG methylation. Non-CG methylated paralogs had higher rates of sequence evolution, higher frequency of presence-absence variation, and more limited expression. This suggests that silencing by non-CG methylation may be important to maintaining dosage following duplication and be a precursor to fractionation. Our results indicate that genic methylation marks differing evolutionary trajectories and fates between paralogous genes and have a role in maintaining dosage following duplication.

Phylotranscriptomic analyses reveal multiple whole-genome duplication events, the history of diversification and adaptations in the Araceae

BACKGROUND AND AIMS

The Araceae are one of the most diverse monocot families with numerous morphological and ecological novelties. Plastid and mitochondrial genes have been used to investigate the phylogeny and to interpret shifts in the pollination biology and biogeography of the Araceae. In contrast, the role of whole-genome duplication (WGD) in the evolution of eight subfamilies remains unclear.

METHODS

New transcriptomes or low-depth whole-genome sequences of 65 species were generated through Illumina sequencing. We reconstructed the phylogenetic relationships of Araceae using concatenated and species tree methods, and then estimated the age of major clades using TreePL. We inferred the WGD events by Ks and gene tree methods. We investigated the diversification patterns applying time-dependent and trait-dependent models. The expansions of gene families and functional enrichments were analysed using CAFE and InterProScan.

KEY RESULTS

Gymnostachydoideae was the earliest diverging lineage followed successively by Orontioideae, Lemnoideae and Lasioideae. In turn, they were followed by the clade of 'bisexual climbers' comprised of Pothoideae and Monsteroideae, which was resolved as the sister to the unisexual flowers clade of Zamioculcadoideae and Aroideae. A special WGD event ψ (psi) shared by the True-Araceae clade occurred in the Early Cretaceous. Net diversification rates first declined and then increased through time in the Araceae. The best diversification rate shift along the stem lineage of the True-Araceae clade was detected, and net diversification rates were enhanced following the ψ-WGD. Functional enrichment analyses revealed that some genes, such as those encoding heat shock proteins, glycosyl hydrolase and cytochrome P450, expanded within the True-Araceae clade.

CONCLUSIONS

Our results improve our understanding of aroid phylogeny using the large number of single-/low-copy nuclear genes. In contrast to the Proto-Araceae group and the lemnoid clade adaption to aquatic environments, our analyses of WGD, diversification and functional enrichment indicated that WGD may play a more important role in the evolution of adaptations to tropical, terrestrial environments in the True-Araceae clade. These insights provide us with new resources to interpret the evolution of the Araceae.

Duckweeds

Mateo-Elizalde et al. introduce duckweeds, a family of freshwater plants.

Engineering triacylglycerol accumulation in duckweed (Lemna japonica)

Duckweeds are amongst the fastest growing of higher plants, making them attractive high-biomass targets for biofuel feedstock production. Their fronds have high rates of fatty acid synthesis to meet the demand for new membranes, but triacylglycerols (TAG) only accumulate to very low levels. Here we report on the engineering of Lemna japonica for the synthesis and accumulation of TAG in its fronds. This was achieved by expression of an estradiol-inducible cyan fluorescent protein-Arabidopsis WRINKLED1 fusion protein (CFP-AtWRI1), strong constitutive expression of a mouse diacylglycerol:acyl-CoA acyltransferase2 (MmDGAT), and a sesame oleosin variant (SiOLE(*)). Individual expression of each gene increased TAG accumulation by 1- to 7-fold relative to controls, while expression of pairs of these genes increased TAG by 7- to 45-fold. In uninduced transgenics containing all three genes, TAG accumulation increased by 45-fold to 3.6% of dry weight (DW) without severely impacting growth, and by 108-fold to 8.7% of DW after incubation on medium containing 100 μm estradiol for 4 days. TAG accumulation was accompanied by an increase in total fatty acids of up to three-fold to approximately 15% of DW. Lipid droplets from fronds of all transgenic lines were visible by confocal microscopy of BODIPY-stained fronds. At a conservative 12 tonnes (dry matter) per acre and 10% (DW) TAG, duckweed could produce 350 gallons of oil/acre/year, approximately seven-fold the yield of soybean, and similar to that of oil palm. These findings provide the foundation for optimizing TAG accumulation in duckweed and present a new opportunity for producing biofuels and lipidic bioproducts.

Duckweeds for Phytoremediation of Polluted Water

Tiny aquatic plants from the Lemnaceae family, commonly known as duckweeds, are often regarded as detrimental to the environment because of their ability to quickly populate and cover the surfaces of bodies of water. Due to their rapid vegetative propagation, duckweeds have one of the fastest growth rates among flowering plants and can accumulate large amounts of biomass in relatively short time periods. Due to the high yield of valuable biomass and ease of harvest, duckweeds can be used as feedstock for biofuels, animal feed, and other applications. Thanks to their efficient absorption of nitrogen- and phosphate-containing pollutants, duckweeds play an important role in the restorative ecology of water reservoirs. Moreover, compared to other species, duckweed species and ecotypes demonstrate exceptionally high adaptivity to a variety of environmental factors; indeed, duckweeds remove and convert many contaminants, such as nitrogen, into plant biomass. The global distribution of duckweeds and their tolerance of ammonia, heavy metals, other pollutants, and stresses are the major factors highlighting their potential for use in purifying agricultural, municipal, and some industrial wastewater. In summary, duckweeds are a powerful tool for bioremediation that can reduce anthropogenic pollution in aquatic ecosystems and prevent water eutrophication in a simple, inexpensive ecologically friendly way. Here we review the potential for using duckweeds in phytoremediation of several major water pollutants: mineral nitrogen and phosphorus, various organic chemicals, and heavy metals.

Ammonium Uptake, Mediated by Ammonium Transporters, Mitigates Manganese Toxicity in Duckweed, Spirodela polyrhiza

Nitrogen is an essential nutrient that affects all aspects of the growth, development and metabolic responses of plants. Here we investigated the influence of the two major sources of inorganic nitrogen, nitrate and ammonium, on the toxicity caused by excess of Mn in great duckweed, Spirodela polyrhiza. The revealed alleviating effect of ammonium on Mn-mediated toxicity, was complemented by detailed molecular, biochemical and evolutionary characterization of the species ammonium transporters (AMTs). Four genes encoding AMTs in S. polyrhiza, were classified as SpAMT1;1, SpAMT1;2, SpAMT1;3 and SpAMT2. Functional testing of the expressed proteins in yeast and Xenopus oocytes clearly demonstrated activity of SpAMT1;1 and SpAMT1;3 in transporting ammonium. Transcripts of all SpAMT genes were detected in duckweed fronds grown in cultivation medium, containing a physiological or 50-fold elevated concentration of Mn at the background of nitrogen or a mixture of nitrate and ammonium. Each gene demonstrated an individual expression pattern, revealed by RT-qPCR. Revealing the mitigating effect of ammonium uptake on manganese toxicity in aquatic duckweed S. polyrhiza, the study presents a comprehensive analysis of the transporters involved in the uptake of ammonium, shedding a new light on the interactions between the mechanisms of heavy metal toxicity and the regulation of the plant nitrogen metabolism.

Dirigent protein subfamily function and structure in terrestrial plant phenol metabolism

Aquatic plant transition to land, and subsequent terrestrial plant species diversification, was accompanied by the emergence and massive elaboration of plant phenol chemo-diversity. Concomitantly, dirigent protein (DP) and dirigent-like protein subfamilies, derived from large multigene families, emerged and became extensively diversified. DP biochemical functions as gateway entry points into new and diverse plant phenol skeletal types then markedly expanded. DPs have at least eight non-uniformly distributed subfamilies, with different DP subfamily members of known biochemical/physiological function now implicated as gateway entries to lignan, lignin, aromatic diterpenoid, pterocarpan and isoflavene pathways. While some other DP subfamily members have jacalin domains, both these and indeed the majority of DPs throughout the plant kingdom await discovery of their biochemical roles. Methods and approaches were developed to discover DP biochemical function as gateway entry points to distinct plant phenol skeletal types in land plants. Various DP 3D X-ray structural determinations enabled structure-based comparative sequence analysis and modeling to understand similarities and differences among the different DP subfamilies. We consider that the core DP β-barrel fold and associated characteristics are likely common to all DPs, with several residues conserved and nearly invariant. There is also considerable variation in residue composition and topography of the putative substrate binding pockets, as well as substantial differences in several loops, such as the β1-β2 loop. All DPs likely bind and stabilize quinone methide intermediates, while guiding distinctive regio- and/or stereo-chemical entry into Nature's chemo-diverse land plant phenol metabolic classes.

Examination of the Metallothionein Gene Family in Greater Duckweed Spirodela polyrhiza

Duckweeds are aquatic plants that proliferate rapidly in a wide range of freshwaters, and they are regarded as a potential source of sustainable biomass for various applications and the cost-effective bioremediation of heavy metal pollutants. To understand the cellular and molecular basis that underlies the high metal tolerance and accumulation capacity of duckweeds, we examined the forms and transcript profiles of the metallothionein (MT) gene family in the model duckweed Spirodela polyrhiza, whose genome has been completely sequenced. Four S. polyrhiza MT-like genes were identified and annotated as SpMT2a, SpMT2b, SpMT3, and SpMT4. All except SpMT2b showed high sequence homology including the conserved cysteine residues with the previously described MTs from flowering plants. The S. polyrhiza genome appears to lack the root-specific Type 1 MT. The transcripts of SpMT2a, SpMT2b, and SpMT3 could be detected in the vegetative whole-plant tissues. The transcript abundance of SpMT2a was upregulated several-fold in response to cadmium stress, and the heterologous expression of SpMT2a conferred copper and cadmium tolerance to the metal-sensitive ∆cup1 strain of Saccharomyces cerevisiae. Based on these results, we proposed that SpMT2a may play an important role in the metal detoxification mechanism of duckweed.

Characterization of defense responses against bacterial pathogens in duckweeds lacking EDS1

ENHANCED DISEASE SUSCEPTIBILITY 1 (EDS1) mediates the induction of defense responses against pathogens in most angiosperms. However, it has recently been shown that a few species have lost EDS1. It is unknown how defense against disease unfolds and evolves in the absence of EDS1. We utilize duckweeds; a collection of aquatic species that lack EDS1, to investigate this question. We established duckweed-Pseudomonas pathosystems and used growth curves and microscopy to characterize pathogen-induced responses. Through comparative genomics and transcriptomics, we show that the copy number of infection-associated genes and the infection-induced transcriptional responses of duckweeds differ from other model species. Pathogen defense in duckweeds has evolved along different trajectories than in other plants, including genomic and transcriptional reprogramming. Specifically, the miAMP1 domain-containing proteins, which are absent in Arabidopsis, showed pathogen responsive upregulation in duckweeds. Despite such divergence between Arabidopsis and duckweed species, we found conservation of upregulation of certain genes and the role of hormones in response to disease. Our work highlights the importance of expanding the pool of model species to study defense responses that have evolved in the plant kingdom independent of EDS1.

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

BACKGROUND AND AIMS

Cycads are regarded as an ancient lineage of living seed plants, and hold important clues to understand the early evolutionary trends of seed plants. The molecular phylogeny and spatio-temporal diversification of one of the species-rich genera of cycads, Macrozamia, have not been well reconstructed.

METHODS

We analysed a transcriptome dataset of 4740 single-copy nuclear genes (SCGs) of 39 Macrozamia species and two outgroup taxa. Based on concatenated (maximum parsimony, maximum likelihood) and multispecies coalescent analyses, we first establish a well-resolved phylogenetic tree of Macrozamia. To identify cyto-nuclear incongruence, the plastid protein coding genes (PCGs) from transcriptome data are extracted using the software HybPiper. Furthermore, we explore the biogeographical history of the genus and shed light on the pattern of floristic exchange between three distinct areas of Australia. Six key diagnostic characters are traced on the phylogenetic framework using two comparative methods, and infra-generic classification is investigated.

KEY RESULTS

The tree topologies of concatenated and multi-species coalescent analyses of SCGs are mostly congruent with a few conflicting nodes, while those from plastid PCGs show poorly supported relationships. The genus contains three major clades that correspond to their distinct distributional areas in Australia. The crown group of Macrozamia is estimated to around 11.80 Ma, with a major expansion in the last 5-6 Myr. Six morphological characters show homoplasy, and the traditional phenetic sectional division of the genus is inconsistent with this current phylogeny.

CONCLUSIONS

This first detailed phylogenetic investigation of Macrozamia demonstrates promising prospects of SCGs in resolving phylogenetic relationships within cycads. Our study suggests that Macrozamia, once widely distributed in Australia, underwent major extinctions because of fluctuating climatic conditions such as cooling and mesic biome disappearance in the past. The current close placement of morphologically distinct species in the phylogenetic tree may be related to neotenic events that occurred in the genus.

Expansion and contraction of small RNA and methylation machinery throughout plant evolution

The revolution in sequencing has created a wealth of plant genomes that can be mined to understand the evolution of biological complexity. Complexity is often driven by gene duplication, which allows paralogs to specialize in an activity of the ancestral gene or acquire novel functions. Angiosperms encode a variety of gene silencing pathways that share related machinery for small RNA biosynthesis and function. Recent phylogenetic analysis of these gene families plots the expansion, specialization, and occasional contraction of this core machinery. This analysis reveals the ancient origin of RNA-directed DNA Methylation in early land plants, or possibly their algal ancestors, as well as ongoing duplications that evolve novel small RNA pathways.

Evolution of the DEHYDRATION-RESPONSIVE ELEMENT-BINDING PROTEIN subfamily in green plants

Adapting to unfavorable environments is a necessary step in plant terrestrialization and radiation. The dehydration-responsive element-binding (DREB) protein subfamily plays a pivotal role in plant abiotic stress regulation. However, relationships between the origin and expansion of the DREB subfamily and adaptive evolution of land plants are still being elucidated. Here, we constructed the evolutionary history of the DREB subfamily by compiling APETALA2/ethylene-responsive element-binding protein superfamily genes from 169 representative species of green plants. Through extensive phylogenetic analyses and comparative genomic analysis, our results revealed that the DREB subfamily diverged from the ethylene-responsive factor (ERF) subfamily in the common ancestor of Zygnemophyceae and Embryophyta during the colonization of land by plants, followed by expansions to form three different ancient archetypal genes in Zygnemophyceae species, designated as groups archetype-I, archetype-II/III, and archetype-IV. Four large-scale expansions paralleling the evolution of land plants led to the nine-subgroup divergence of group archetype-II/III in angiosperms, and five whole-genome duplications during Brassicaceae and Poaceae radiation shaped the diversity of subgroup IIb-1. We identified a Poaceae-specific gene in subgroup IIb-1, ERF014, remaining in a Poaceae-specific microsynteny block and co-evolving with a small heat shock protein cluster. Expression analyses demonstrated that heat acclimation may have driven the neofunctionalization of ERF014s in Pooideae by engaging in the conserved heat-responsive module in Poaceae. This study provides insights into lineage-specific expansion and neofunctionalization in the DREB subfamily, together with evolutionary information valuable for future functional studies of plant stress biology.

Genome-wide identification of nitrate transporter genes from Spirodela polyrhiza and characterization of SpNRT1.1 function in plant development

Nitrate transporter (NRT) genes that participate in nitrate transport and distribution are indispensable for plant growth, development, and stress tolerance. Spirodela polyrhiza has the smallest genome among monocotyledon plants, and it has strong nitrate absorbance and phytoremediation abilities. However, the evolutionary history, expression patterns, and functions of the NRT gene family in S. polyrhiza are not well understood. Here, we identified 29 NRT members in the S. polyrhiza genome. Gene structure and phylogeny analyses showed that S. polyrhiza nitrate transporter (SpNRTs) genes were divided into eight clades without gene expansion compared with that in Arabidopsis. Transcriptomic analysis showed that SpNRT genes have spatiotemporal expression patterns and respond to abiotic stress. Functional analysis revealed that in S. polyrhiza, SpNRT1.1 expression was strongly induced by treatment with nitrate and ammonium. Overexpression of SpNRT1.1 significantly repressed primary root length, and the number and total length of lateral roots. This was more pronounced in high ammonium concentration medium. Overexpressed SpNRT1.1 in Arabidopsis significantly improved biomass and delayed flowering time, indicating that the nitrate transport ability of SpNRT1.1 differs from AtNRT1.1. In conclusion, our results provide valuable information about the evolution of the NRT family in higher plants and the function of SpNRT1.1.

Evolutionary analysis and functional characterization of BZR1 gene family in celery revealed their conserved roles in brassinosteroid signaling

BACKGROUND

Brassinosteroids (BRs) are a group of essential steroid hormones involved in diverse developmental and physiological processes in plants. The Brassinazole-resistant 1 (BZR1) transcription factors are key components of BR signaling and integrate a wide range of internal and environmental signals to coordinate plant development, growth, and resistance to abiotic and biotic stresses. Although the BZR1 family has been fully studied in Arabidopsis, celery BZR1 family genes remain largely unknown.

RESULTS

Nine BZR1 genes were identified in the celery genome, and categorized into four classes based on phylogenetic and gene structure analyses. All the BZR1 proteins shared a typical bHLH (basic helix-loop-helix) domain that is highly conserved across the whole family in Arabidopsis, grape, lettuce, ginseng, and three Apiaceae species. Both duplications and losses of the BZR1 gene family were detected during the shaping of the celery genome. Whole-genome duplication (WGD) or segmental duplication contributed 55.56% of the BZR1 genes expansion, and the γ as well as celery-ω polyploidization events made a considerable contribution to the production of the BZR1 paralogs in celery. Four AgBZR1 members (AgBZR1.1, AgBZR1.3, AgBZR1.5, and AgBZR1.9), which were localized both in the nucleus and cytoplasm, exhibit transcription activation activity in yeast. AgBZR1.5 overexpression transgenic plants in Arabidopsis showed curled leaves with bent, long petioles and constitutive BR-responsive phenotypes. Furthermore, the AgBZR1 genes possessed divergent expression patterns with some overlaps in roots, petioles, and leaves, suggesting an extensive involvement of AgBZR1s in the developmental processes in celery with both functional redundancy and divergence.

CONCLUSIONS

Our results not only demonstrated that AgBZR1 played a conserved role in BR signaling but also suggested that AgBZR1 might be extensively involved in plant developmental processes in celery. The findings lay the foundation for further study on the molecular mechanism of the AgBZR1s in regulating the agronomic traits and environmental adaptation of celery, and provide insights for future BR-related genetic breeding of celery and other Apiaceae crops.

Phenotypic plasticity promotes species coexistence

Ecological explanations for species coexistence assume that species' traits, and therefore the differences between species, are fixed on short timescales. However, species' traits are not fixed, but can instead change rapidly as a consequence of phenotypic plasticity. Here we use a combined experimental-theoretical approach to demonstrate that plasticity in response to interspecific competition between two aquatic plants allows for species coexistence where competitive exclusion is otherwise predicted to occur. Our results show that rapid trait changes in response to a shift in the competitive environment can promote coexistence in a way that is not captured by common measures of niche differentiation.

Metabolome and transcriptome analyses of the flavonoid biosynthetic pathway for the efficient accumulation of anthocyanins and other flavonoids in a new duckweed variety (68-red)

Duckweed is a kind of aquatic plant with the characteristics of high nutritional value and medicinal benefits. However, most researches focused on the natural germplasms. The underlying metabolic pathway remains to be systematically elaborated in duckweed. In our laboratory, one reddish-purple mutant with high-flavonoids was screened from a mutant library of Spirodela polyrhiza 6068, named 68-red. The content of anthocyanins and proanthocyanidins in 68-red mutant increased by 563.47% and 231.19%, respectively, compared to wild type. It is interesting that cynaroside and orientin content were significantly increased, in contrast, apigetrin and vitexin were decreased in 68-red mutant. Considering this, metabolome and transcriptome were employed to explore the flavonoids biosynthetic pathway. Here, a total of 734 metabolites were identified in the wild type and 68-red mutant. Among which, cyanidin-3-O-glucoside, cyanidin-3-O-galactoside, pelargonidin-3-O-glucoside and pelargonidin-3-O-(6″-O-malonyl)glucoside were significantly accumulated, which were positively correlated with deep reddish-purple of 68-red mutant. In addition, proanthocyanidins (B1, B2, B3, B4, C1, C2), flavonoid and its glycosides (11 luteolin and its glycosides, 14 quercetin and its glycosides, 14 kaempferol and its glycosides, 2 apigenin glycosides) were significantly accumulated, 2 apigenin glycosides were down-regulated in 68-red mutant. The transcriptome data and qRT-PCR indicated that 16 enzyme genes in flavonoids biosynthetic pathway (PAL, C4H, CHSs, F3H, ANS, ANR, F3'Hs, DFRs, LAR, GT1, BZ1) were significantly up-regulated in 68-red mutant. Correlation analysis found that three copies of F3'H gene play important roles in the synthesis of anthocyanins, luteolin and apigenin glycosides. In conclusion, the 68-red mutant is a high quality germplasm resources for food and medical industry. Metabolome and transcriptome provide new insight for exploring the enzyme genes and functional metabolites in duckweed.

Identification and expression analysis of GARP superfamily genes in response to nitrogen and phosphorus stress in Spirodela polyrhiza

BACKGROUND

GARP transcription factors perform critical roles in plant development and response to environmental stimulus, especially in the phosphorus (P) and nitrogen (N) sensing and uptake. Spirodela polyrhiza (giant duckweed) is widely used for phytoremediation and biomass production due to its rapid growth and efficient N and P removal capacities. However, there has not yet been a comprehensive analysis of the GRAP gene family in S. polyrhiza.

RESULTS

We conducted a comprehensive study of GRAP superfamily genes in S. polyrhiza. First, we investigated 35 SpGARP genes which have been classified into three groups based on their gene structures, conserved motifs, and phylogenetic relationship. Then, we identified the duplication events, performed the synteny analysis, and calculated the K_a/K_s ratio in these SpGARP genes. The regulatory and co-expression networks of SpGARPs were further constructed using cis-acting element analysis and weighted correlation network analysis (WGCNA). Finally, the expression pattern of SpGARP genes were analyzed using RNA-seq data and qRT-PCR, and several NIGT1 transcription factors were found to be involved in both N and P starvation responses.

CONCLUSIONS

The study provides insight into the evolution and function of GARP superfamily in S. polyrhiza, and lays the foundation for the further functional verification of SpGARP genes.

Pervasive genome duplications across the plant tree of life and their links to major evolutionary innovations and transitions

Whole-genome duplication (WGD) has occurred repeatedly during plant evolution and diversification, providing genetic layers for evolving new functions and phenotypes. Advances in long-read sequencing technologies have enabled sequencing and assembly of over 1000 plant genomes spanning nearly 800 species, in which a large set of ancient WGDs has been uncovered. Here, we review the recently reported WGDs that occurred in major plant lineages and key evolutionary positions, and highlight their contributions to morphological innovation and adaptive evolution. Current gaps and challenges in integrating enormous volumes of sequenced plant genomes, accurately inferring WGDs, and developing web-based analysis tools are emphasized. Looking to the future, ambitious genome sequencing projects and global efforts may substantially recapitulate the plant tree of life based on broader sampling of phylogenetic diversity, reveal much of the timetable of ancient WGDs, and address the biological significance of WGDs in plant adaptation and radiation.