Streamlined spatial and environmental expression signatures characterize the minimalist duckweed Wolffia australiana.

Single-cell genomics permits a new resolution in the examination of molecular and cellular dynamics, allowing global, parallel assessments of cell types and cellular behaviors through development and in response to environmental circumstances, such as interaction with water and the light-dark cycle of the Earth. Here, we leverage the smallest, and possibly most structurally reduced, plant, the semiaquatic Wolffia australiana, to understand dynamics of cell expression in these contexts at the whole-plant level. We examined single-cell-resolution RNA-sequencing data and found Wolffia cells divide into four principal clusters representing the above- and below-water-situated parenchyma and epidermis. Although these tissues share transcriptomic similarity with model plants, they display distinct adaptations that Wolffia has made for the aquatic environment. Within this broad classification, discrete subspecializations are evident, with select cells showing unique transcriptomic signatures associated with developmental maturation and specialized physiologies. Assessing this simplified biological system temporally at two key time-of-day (TOD) transitions, we identify additional TOD-responsive genes previously overlooked in whole-plant transcriptomic approaches and demonstrate that the core circadian clock machinery and its downstream responses can vary in cell-specific manners, even in this simplified system. Distinctions between cell types and their responses to submergence and/or TOD are driven by expression changes of unexpectedly few genes, characterizing Wolffia as a highly streamlined organism with the majority of genes dedicated to fundamental cellular processes. Wolffia provides a unique opportunity to apply reductionist biology to elucidate signaling functions at the organismal level, for which this work provides a powerful resource.

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.

Gene expression and metabolite levels converge in the thermogenic spadix of skunk cabbage.

The inflorescence (spadix) of skunk cabbage (Symplocarpus renifolius) is strongly thermogenic and can regulate its temperature at around 23 °C even when the ambient temperature drops below freezing. To elucidate the mechanisms underlying developmentally controlled thermogenesis and thermoregulation in skunk cabbage, we conducted a comprehensive transcriptome and metabolome analysis across 3 developmental stages of spadix development. Our RNA-seq analysis revealed distinct groups of expressed genes, with selenium-binding protein 1/methanethiol oxidase (SBP1/MTO) exhibiting the highest levels in thermogenic florets. Notably, the expression of alternative oxidase (AOX) was consistently high from the prethermogenic stage through the thermogenic stage in the florets. Metabolome analysis showed that alterations in nucleotide levels correspond with the developmentally controlled and tissue-specific thermogenesis of skunk cabbage, evident by a substantial increase in AMP levels in thermogenic florets. Our study also reveals that hydrogen sulfide, a product of SBP1/MTO, inhibits cytochrome c oxidase (COX)-mediated mitochondrial respiration, while AOX-mediated respiration remains relatively unaffected. Specifically, at lower temperatures, the inhibitory effect of hydrogen sulfide on COX-mediated respiration increases, promoting a shift toward the dominance of AOX-mediated respiration. Finally, despite the differential regulation of genes and metabolites throughout spadix development, we observed a convergence of gene expression and metabolite accumulation patterns during thermogenesis. This synchrony may play a key role in developmentally regulated thermogenesis. Moreover, such convergence during the thermogenic stage in the spadix may provide a solid molecular basis for thermoregulation in skunk cabbage.

Auxin and carbohydrate control flower bud development in Anthurium andraeanum during early stage of sexual reproduction.

BACKGROUND: Flower buds of Anthurium andraeanum frequently cease to grow and abort during the early flowering stage, resulting in prolonged planting times and increased commercialization costs. Nevertheless, limited knowledge exists of the mechanism of flower development after initiation in A. andraeanum. RESULTS: In this study, the measurement of carbohydrate flow and intensity between leaves and flowers during different growth stages showed that tender leaves are strong sinks and their concomitant flowers are weak ones. This suggested that the tender leaves compete with their concomitant flower buds for carbohydrates during the early growth stages, potentially causing the abortion of the flower buds. The analysis of transcriptomic differentially expressed genes suggested that genes related to sucrose metabolism and auxin response play an important role during flower bud development. Particularly, co-expression network analysis found that AaSPL12 is a hub gene engaged in flower development by collaborating carbohydrate and auxin signals. Yeast Two Hybrid assays revealed that AaSPL12 can interact with AaARP, a protein that serves as an indicator of dormancy. Additionally, the application of exogenous IAA and sucrose can suppress the expression of AaARP, augment the transcriptional abundance of AaSPL12, and consequently expedite flower development in Anthurium andraeanum. CONCLUSIONS: Collectively, our findings indicated that the combination of auxin and sugar signals could potentially suppress the repression of AaARP protein to AaSPL12, thus advancing the development of flower buds in Anthurium andraeanum.

A non-cell-autonomous circadian rhythm of bioluminescence reporter activities in individual duckweed cells.

The circadian clock is responsible for the temporal regulation of various physiological processes in plants. Individual cells contain a circadian oscillator consisting of a clock gene circuit that coordinates physiological rhythms within the plant body in an orderly manner. The coordination of time information has been studied from the perspective of cell-cell local coupling and long-distance communication between tissues based on the view that the behavior of circadian oscillators represents physiological rhythms. Here, we report the cellular circadian rhythm of bioluminescence reporters that are not governed by the clock gene circuit in expressing cells. We detected cellular bioluminescence rhythms with different free-running periods in the same cells using a dual-color bioluminescence monitoring system in duckweed (Lemna minor) transfected with Arabidopsis CIRCADIAN CLOCK ASSOCIATED 1::luciferace+ (AtCCA1::LUC+) and Cauliflower mosaic virus 35S::modified click-beetle red-color luciferase (CaMV35S::PtRLUC) reporters. Co-transfection experiments with the two reporters and a clock gene-overexpressing effector revealed that the AtCCA1::LUC+ rhythm, but not the CaMV35S::PtRLUC rhythm, was altered in cells with a dysfunctional clock gene circuit. This indicated that the AtCCA1::LUC+ rhythm is a direct output of the cellular circadian oscillator, whereas the CaMV35S::PtRLUC rhythm is not. After plasmolysis, the CaMV35S::PtRLUC rhythm disappeared, whereas the AtCCA1::LUC+ rhythm persisted. This suggests that the CaMV35S::PtRLUC bioluminescence has a symplast/apoplast-mediated circadian rhythm generated at the organismal level. The CaMV35S::PtRLUC-type bioluminescence rhythm was also observed when other bioluminescence reporters were expressed. These results reveal that the plant circadian system consists of both cell-autonomous and noncell-autonomous rhythms that are unaffected by cellular oscillators.

Characterization of the cryptic interspecific hybrid Lemna×mediterranea by an integrated approach provides new insights into duckweed diversity.

Lemnaceae taxonomy is challenged by the particular morphology of these tiny free-floating angiosperms. Although molecular taxonomy has helped clarify the phylogenetic history of this family, some inconsistency with morphological data leads to frequent misclassifications in the genus Lemna. Recently, the finding that Lemna japonica is an interspecific hybrid between Lemna minor and Lemna turionifera provided a clear explanation for one such taxonomic question. Here we demonstrated that L. minor is also capable of hybridizing with Lemna gibba, generating a cryptic but widespread taxon in the Mediterranean area. The nothotaxon Lemna ×mediterranea is described and compared with clones of the putative parental species L. minor and L. gibba. Genetic analysis by nuclear and plastid markers, as well as genome size measurement, revealed that two different cytotypes, diploid and triploid, originated by at least two independent hybridization events. Despite high overall similarity, morphometrical, physiological, and biochemical analyses showed an intermediate position of L. ×mediterranea between its parental species in most qualitative and quantitative characters, and also separation of the two hybrid cytotypes by some criteria. These data provide evidence that hybridization and polyploidization, driving forces of terrestrial plant evolution, contribute to duckweed genetic diversity and may have shaped the phylogenetic history of these mainly asexual, aquatic plants.

Return of the Lemnaceae: duckweed as a model plant system in the genomics and postgenomics era.

The aquatic Lemnaceae family, commonly called duckweed, comprises some of the smallest and fastest growing angiosperms known on Earth. Their tiny size, rapid growth by clonal propagation, and facile uptake of labeled compounds from the media were attractive features that made them a well-known model for plant biology from 1950 to 1990. Interest in duckweed has steadily regained momentum over the past decade, driven in part by the growing need to identify alternative plants from traditional agricultural crops that can help tackle urgent societal challenges, such as climate change and rapid population expansion. Propelled by rapid advances in genomic technologies, recent studies with duckweed again highlight the potential of these small plants to enable discoveries in diverse fields from ecology to chronobiology. Building on established community resources, duckweed is reemerging as a platform to study plant processes at the systems level and to translate knowledge gained for field deployment to address some of society's pressing needs. This review details the anatomy, development, physiology, and molecular characteristics of the Lemnaceae to introduce them to the broader plant research community. We highlight recent research enabled by Lemnaceae to demonstrate how these plants can be used for quantitative studies of complex processes and for revealing potentially novel strategies in plant defense and genome maintenance.

The evolution of the duckweed ionome mirrors losses in structural complexity.

BACKGROUND AND AIMS: The duckweeds (Lemnaceae) consist of 36 species exhibiting impressive phenotypic variation, including the progressive evolutionary loss of a fundamental plant organ, the root. Loss of roots and reduction of vascular tissues in recently derived taxa occur in concert with genome expansions of ≤14-fold. Given the paired loss of roots and reduction in structural complexity in derived taxa, we focus on the evolution of the ionome (whole-plant elemental contents) in the context of these fundamental changes in body plan. We expect that progressive vestigiality and eventual loss of roots might have both adaptive and maladaptive consequences that are hitherto unknown. METHODS: We quantified the ionomes of 34 accessions in 21 species across all duckweed genera, spanning 70 Myr in this rapidly cycling plant (doubling times are as rapid as ~24 h). We related both micro- and macroevolutionary ionome contrasts to body plan remodelling and showed nimble microevolutionary shifts in elemental accumulation and exclusion in novel accessions. KEY RESULTS: We observed a robust directional trend in calcium and magnesium levels, decreasing from the ancestral representative Spirodela genus towards the derived rootless Wolffia, with the latter also accumulating cadmium. We also identified abundant within-species variation and hyperaccumulators of specific elements, with this extensive variation at the fine (as opposed to broad) scale. CONCLUSIONS: These data underscore the impact of root loss and reveal the very fine scale of microevolutionary variation in hyperaccumulation and exclusion of a wide range of elements. Broadly, they might point to trade-offs not well recognized in ionomes.

The immediate metabolomic effects of whole-genome duplication in the greater duckweed, Spirodela polyrhiza.

PREMISE: In plants, whole-genome duplication (WGD) is a common mutation with profound evolutionary potential. Given the costs associated with a superfluous genome copy, polyploid establishment is enigmatic. However, in the right environment, immediate phenotypic changes following WGD can facilitate establishment. Metabolite abundances are the direct output of the cell's regulatory network and determine much of the impact of environmental and genetic change on the phenotype. While it is well known that an increase in the bulk amount of genetic material can increase cell size, the impact of gene dosage multiplication on the metabolome remains largely unknown. METHODS: We used untargeted metabolomics on four genetically distinct diploid-neoautotetraploid pairs of the greater duckweed, Spirodela polyrhiza, to investigate how WGD affects metabolite abundances per cell and per biomass. RESULTS: Autopolyploidy increased metabolite levels per cell, but the response of individual metabolites varied considerably. However, the impact on metabolite level per biomass was restricted because the increased cell size reduced the metabolite concentration per cell. Nevertheless, we detected both quantitative and qualitative effects of WGD on the metabolome. Many effects were strain-specific, but some were shared by all four strains. CONCLUSIONS: The nature and impact of metabolic changes after WGD depended strongly on the genotype. Dosage effects have the potential to alter the plant metabolome qualitatively and quantitatively, but were largely balanced out by the reduction in metabolite concentration due to an increase in cell size in this species.

Actinomycetospora lemnae sp. nov., A Novel Actinobacterium Isolated from Lemna aequinoctialis Able to Enhance Duckweed Growth.

Duckweed-associated actinobacteria are co-existing microbes that affect duckweed growth and adaptation. In this study, we aimed to report a novel actinobacterium species and explore its ability to enhance duckweed growth. Strain DW7H6T was isolated from duckweed, Lemna aequinoctialis. Phylogenetic analysis based on its 16S rRNA gene sequence revealed that the strain was most closely related to Actinomycetospora straminea IY07-55T (99.0%), Actinomycetospora chibensis TT04-21T (98.9%), Actinomycetospora lutea TT00-04T (98.8%) and Actinomycetospora callitridis CAP 335T (98.4%). Chemotaxonomic and morphological characteristics of strain DW7H6T were consistent with members of the genus Actinomycetospora, while average nucleotide identity (ANI) and digital DNA-DNA hybridization (dDDH) between the draft genomes of this strain and its closely related type strains were below the proposed threshold values used for species discrimination. Based on chemotaxonomic, phylogenetic, phenotypic, and genomic evidence obtained, we describe a novel Actinomycetospora species, for which the name Actinomycetospora lemnae sp. nov. is proposed. The type strain is DW7H6T (TBRC 15165T, NBRC 115294T). Additionally, the duckweed-associated actinobacterium strain DW7H6T was able to enhance duckweed growth when compared to the control, in which the number of fronds and biomass dry weight were increased by up to 1.4 and 1.3 fold, respectively. Moreover, several plant-associated gene features in the genome of strain DW7H6T potentially involved in plant-microbe interactions were identified.

Two lineages of Lemna aequinoctialis (Araceae, Lemnoideae) based on physiology, morphology, and phylogeny.

Lemna aequinoctialis Welw. is a widely spread species that has diverse physiological and molecular properties. Flower characteristics are important factors in deducing taxonomical status; however, owing to the rarity of flowering observations in Lemna, studying them has been a prolonged challenge. In this study, physiological and morphological analyses were conducted by inducing flowering, and molecular analysis was done based on the two chloroplast DNA loci (matK, atpF-atpH intergeneric spacer) of L. aequinoctialis sensu Landolt (1986) from 70 strains found in 70 localities in Japan, Korea, Thailand, and the US. In total, 752 flowering fronds from 13 strains were observed based on axenic conditions. Two different trends in flower organ development-protogyny and adichogamy-were detected in these strains. Their physiological traits were divided into two groups, showing different morphological features based on frond thickness, root cap, and anther sizes. Molecular analysis showed two lineages corresponding to two physiological groups. These were identified as L. aequinoctialis sensu Beppu et al. (1985) and L. aoukikusa Beppu et Murata based on the description of the nomenclature of L. aoukikusa. These were concluded as independent taxa and can be treated as different species. Furthermore, the distribution of L. aoukikusa is not only limited to Japan.

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.

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.

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.

The genome and preliminary single-nuclei transcriptome of Lemna minuta reveals mechanisms of invasiveness.

The ability to trace every cell in some model organisms has led to the fundamental understanding of development and cellular function. However, in plants the complexity of cell number, organ size, and developmental time makes this a challenge even in the diminutive model plant Arabidopsis (Arabidopsis thaliana). Duckweed, basal nongrass aquatic monocots, provide an opportunity to follow every cell of an entire plant due to their small size, reduced body plan, and fast clonal growth habit. Here we present a chromosome-resolved genome for the highly invasive Lesser Duckweed (Lemna minuta) and generate a preliminary cell atlas leveraging low cell coverage single-nuclei sequencing. We resolved the 360 megabase genome into 21 chromosomes, revealing a core nonredundant gene set with only the ancient tau whole-genome duplication shared with all monocots, and paralog expansion as a result of tandem duplications related to phytoremediation. Leveraging SMARTseq2 single-nuclei sequencing, which provided higher gene coverage yet lower cell count, we profiled 269 nuclei covering 36.9% (8,457) of the L. minuta transcriptome. Since molecular validation was not possible in this nonmodel plant, we leveraged gene orthology with model organism single-cell expression datasets, gene ontology, and cell trajectory analysis to define putative cell types. We found that the tissue that we computationally defined as mesophyll expressed high levels of elemental transport genes consistent with this tissue playing a role in L. minuta wastewater detoxification. The L. minuta genome and preliminary cell map provide a paradigm to decipher developmental genes and pathways for an entire plant.

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.

Target sequence data shed new light on the infrafamilial classification of Araceae.

PREMISE: Recent phylogenetic studies of the Araceae have confirmed the position of the duckweeds nested within the aroids, and the monophyly of a clade containing all the unisexual flowered aroids plus the bisexual-flowered Calla palustris. The main objective of the present study was to better resolve the deep phylogenetic relationships among the main lineages within the family, particularly the relationships between the eight currently recognized subfamilies. We also aimed to confirm the phylogenetic position of the enigmatic genus Calla in relation to the long-debated evolutionary transition between bisexual and unisexual flowers in the family. METHODS: Nuclear DNA sequence data were generated for 128 species across 111 genera (78%) of Araceae using target sequence capture and the Angiosperms 353 universal probe set. RESULTS: The phylogenomic data confirmed the monophyly of the eight Araceae subfamilies, but the phylogenetic position of subfamily Lasioideae remains uncertain. The genus Calla is included in subfamily Aroideae, which has also been expanded to include Zamioculcadoideae. The tribe Aglaonemateae is newly defined to include the genera Aglaonema and Boycea. CONCLUSIONS: Our results strongly suggest that new research on African genera (Callopsis, Nephthytis, and Anubias) and Calla will be important for understanding the early evolution of the Aroideae. Also of particular interest are the phylogenetic positions of the isolated genera Montrichardia, Zantedeschia, and Anchomanes, which remain only moderately supported here.

Assembly and Comparative Analysis of the Complete Mitochondrial Genome of Two Species of Calla Lilies (Zantedeschia, Araceae).

In this study, the mitochondrial genomes of two calla species, Zantedeschia aethiopica Spreng. and Zantedeschia odorata Perry., were assembled and compared for the first time. The Z. aethiopica mt genome was assembled into a single circular chromosome, measuring 675,575 bp in length with a 45.85% GC content. In contrast, the Z. odorata mt genome consisted of bicyclic chromosomes (chromosomes 1 and 2), measuring 719,764 bp and exhibiting a 45.79% GC content. Both mitogenomes harbored similar gene compositions, with 56 and 58 genes identified in Z. aethiopica and Z. odorata, respectively. Analyses of codon usage, sequence repeats, gene migration from chloroplast to mitochondrial, and RNA editing were conducted for both Z. aethiopica and Z. odorata mt genomes. Phylogenetic examination based on the mt genomes of these two species and 30 other taxa provided insights into their evolutionary relationships. Additionally, the core genes in the gynoecium, stamens, and mature pollen grains of the Z. aethiopica mt genome were investigated, which revealed maternal mitochondrial inheritance in this species. In summary, this study offers valuable genomic resources for future research on mitogenome evolution and the molecular breeding of calla lily.

Comparative Physiological and Transcriptome Analyses of Tolerant and Susceptible Cultivars Reveal the Molecular Mechanism of Cold Tolerance in Anthurium andraeanum.

Anthurium andraeanum is a tropical ornamental flower. The cost of Anthurium production is higher under low temperature (non-freezing) conditions; therefore, it is important to increase its cold tolerance. However, the molecular mechanisms underlying the response of Anthurium to cold stress remain elusive. In this study, comparative physiological and transcriptome sequencing analyses of two cultivars with contrasting cold tolerances were conducted to evaluate the cold stress response at the flowering stage. The activities of superoxide dismutase and peroxidase and the contents of proline, soluble sugar, and malondialdehyde increased under cold stress in the leaves of the cold tolerant cultivar Elegang (E) and cold susceptible cultivar Menghuang (MH), while the soluble protein content decreased in MH and increased in E. Using RNA sequencing, 24,695 differentially expressed genes (DEGs) were identified from comparisons between cultivars under the same conditions or between the treatment and control groups of a single cultivar, 9132 of which were common cold-responsive DEGs. Heat-shock proteins and pectinesterases were upregulated in E and downregulated in MH, indicating that these proteins are essential for Anthurium cold tolerance. Furthermore, four modules related to cold treatment were obtained by weighted gene co-expression network analysis. The expression of the top 20 hub genes in these modules was induced by cold stress in E or MH, suggesting they might be crucial contributors to cold tolerance. DEGs were significantly enriched in plant hormone signal transduction pathways, trehalose metabolism, and ribosomal proteins, suggesting these processes play important roles in Anthurium's cold stress response. This study provides a basis for elucidating the mechanism of cold tolerance in A. andraeanum and potential targets for molecular breeding.

Transcriptomic and Functional Analyses of Two Cadmium Hyper-Enriched Duckweed Strains Reveal Putative Cadmium Tolerance Mechanisms.

Cadmium (Cd) is one of the most toxic metals in the environment and exerts deleterious effects on plant growth and production. Duckweed has been reported as a promising candidate for Cd phytoremediation. In this study, the growth, Cd enrichment, and antioxidant enzyme activity of duckweed were investigated. We found that both high-Cd-tolerance duckweed (HCD) and low-Cd-tolerance duckweed (LCD) strains exposed to Cd were hyper-enriched with Cd. To further explore the underlying molecular mechanisms, a genome-wide transcriptome analysis was performed. The results showed that the growth rate, chlorophyll content, and antioxidant enzyme activities of duckweed were significantly affected by Cd stress and differed between the two strains. In the genome-wide transcriptome analysis, the RNA-seq library generated 544,347,670 clean reads, and 1608 and 2045 differentially expressed genes were identified between HCD and LCD, respectively. The antioxidant system was significantly expressed during ribosomal biosynthesis in HCD but not in LCD. Fatty acid metabolism and ethanol production were significantly increased in LCD. Alpha-linolenic acid metabolism likely plays an important role in Cd detoxification in duckweed. These findings contribute to the understanding of Cd tolerance mechanisms in hyperaccumulator plants and lay the foundation for future phytoremediation studies.