Comprehensive definition of genome features in Spirodela polyrhiza by high-depth physical mapping and short-read DNA sequencing strategies

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.

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.

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.

Growth, physiological parameters and DNA methylation in Spirodela polyrhiza (L.) Schleid exposed to PET micro-nanoplastic contaminated waters

The effects of polyethylene terephthalate micro-nanoplastics (PET-MNPs) were tested on the model freshwater species Spirodela polyrhiza (L.) Schleid., with focus on possible particle-induced epigenetic effects (i.e. alteration of DNA methylation status). MNPs (size ∼ 200-300 nm) were produced as water dispersions from PET bottles through repeated cycles of homogenization and used to prepare N-medium at two environmentally relevant concentrations (∼0.05 g L-1 and ∼0.1 g L-1 of MNPs). After 10 days of exposure, a reduction in fresh and dry weight was observed in treated plants, even if the average specific growth rate for both frond number and area was not altered. Impaired growth was coupled with a MNP-induced decrease of chlorophyll fluorescence parameters (i.e. ΨETo and Piabs, indicators of photochemical efficiency) and starch concentration, as well as with alterations in plant ionomic profile and oxidative status. The methylation-sensitive amplification polymorphism (MSAP) technique was used to assess possible changes in DNA methylation levels induced by plastic particles. The analysis showed unusual hypermethylation in 5'-CCGG sites that could be implicated in DNA protection from dangerous agents (i.e. reactive oxygen species) or in the formation of new epialleles. This work represents the first evidence of MNP-induced epigenetic modifications in the plant world.

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.

Diploid and tetraploid genomes of Acorus and the evolution of monocots

Monocots are a major taxon within flowering plants, have unique morphological traits, and show an extraordinary diversity in lifestyle. To improve our understanding of monocot origin and evolution, we generate chromosome-level reference genomes of the diploid Acorus gramineus and the tetraploid Ac. calamus, the only two accepted species from the family Acoraceae, which form a sister lineage to all other monocots. Comparing the genomes of Ac. gramineus and Ac. calamus, we suggest that Ac. gramineus is not a potential diploid progenitor of Ac. calamus, and Ac. calamus is an allotetraploid with two subgenomes A, and B, presenting asymmetric evolution and B subgenome dominance. Both the diploid genome of Ac. gramineus and the subgenomes A and B of Ac. calamus show clear evidence of whole-genome duplication (WGD), but Acoraceae does not seem to share an older WGD that is shared by most other monocots. We reconstruct an ancestral monocot karyotype and gene toolkit, and discuss scenarios that explain the complex history of the Acorus genome. Our analyses show that the ancestors of monocots exhibit mosaic genomic features, likely important for that appeared in early monocot evolution, providing fundamental insights into the origin, evolution, and diversification of monocots.

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.

Loss of ancestral function in duckweed roots is accompanied by progressive anatomical reduction and a re-distribution of nutrient transporters

Organ loss occurs frequently during plant and animal evolution. Sometimes, non-functional organs are retained through evolution. Vestigial organs are defined as genetically determined structures that have lost their ancestral (or salient) function.^1,2,3 Duckweeds, an aquatic monocot family, exhibit both these characteristics. They possess a uniquely simple body plan, variably across five genera, two of which are rootless. Due to the existence of closely related species with a wide diversity in rooting strategies, duckweed roots represent a powerful system for investigating vestigiality. To explore this, we employed a panel of physiological, ionomic, and transcriptomic analyses, with the main goal of elucidating the extent of vestigiality in duckweed roots. We uncovered a progressive reduction in root anatomy as genera diverge and revealed that the root has lost its salient ancestral function as an organ required for supplying nutrients to the plant. Accompanying this, nutrient transporter expression patterns have lost the stereotypical root biased localization observed in other plant species. While other examples of organ loss such as limbs in reptiles⁴ or eyes in cavefish⁵ frequently display a binary of presence/absence, duckweeds provide a unique snapshot of an organ with varying degrees of vestigialization in closely related neighbors and thus provide a unique resource for exploration of how organs behave at different stages along the process of loss.

Optimization of Molecular Methods for Detecting Duckweed-Associated Bacteria

The bacterial colonization dynamics of plants can differ between phylogenetically similar bacterial strains and in the context of complex bacterial communities. Quantitative methods that can resolve closely related bacteria within complex communities can lead to a better understanding of plant-microbe interactions. However, current methods often lack the specificity to differentiate phylogenetically similar bacterial strains. In this study, we describe molecular strategies to study duckweed-associated bacteria. We first systematically optimized a bead-beating protocol to co-isolate nucleic acids simultaneously from duckweed and bacteria. We then developed a generic fingerprinting assay to detect bacteria present in duckweed samples. To detect specific duckweed-bacterium associations, we developed a genomics-based computational pipeline to generate bacterial strain-specific primers. These strain-specific primers differentiated bacterial strains from the same genus and enabled the detection of specific duckweed-bacterium associations present in a community context. Moreover, we used these strain-specific primers to quantify the bacterial colonization of duckweed by normalization to a plant reference gene and revealed differences in colonization levels between strains from the same genus. Lastly, confocal microscopy of inoculated duckweed further supported our PCR results and showed bacterial colonization of the duckweed root-frond interface and root interior. The molecular methods introduced in this work should enable the tracking and quantification of specific plant-microbe associations within plant-microbial communities.

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.

Broad-scale factors shaping the ecological niche and geographic distribution of Spirodela polyrhiza

The choice of appropriate independent variables to create models characterizing ecological niches of species is of critical importance in distributional ecology. This set of dimensions in which a niche is defined can inform about what factors limit the distributional potential of a species. We used a multistep approach to select relevant variables for modeling the ecological niche of the aquatic Spirodela polyrhiza, taking into account variability arising from using distinct algorithms, calibration areas, and spatial resolutions of variables. We found that, even after an initial selection of meaningful variables, the final set of variables selected based on statistical inference varied considerably depending on the combination of algorithm, calibration area, and spatial resolution used. However, variables representing extreme temperatures and dry periods were more consistently selected than others, despite the treatment used, highlighting their importance in shaping the distribution of this species. Other variables related to seasonality of solar radiation, summer solar radiation, and some soil proxies of nutrients in water, were selected commonly but not as frequently as the ones mentioned above. We suggest that these later variables are also important to understanding the distributional potential of the species, but that their effects may be less pronounced at the scale at which they are represented for the needs of this type of modeling. Our results suggest that an informed definition of an initial set of variables, a series of statistical steps for filtering and exploring these predictors, and model selection exercises that consider multiple sets of predictors, can improve determination of variables that shape the niche and distribution of the species, despite differences derived from factors related to data or modeling algorithms.

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.

Strategies for Intraspecific Genotyping of Duckweed: Comparison of Five Orthogonal Methods Applied to the Giant Duckweed Spirodela polyrhiza

The predominantly vegetative propagating duckweeds are of growing commercial interest. Since clonal accessions within a respective species can vary considerably with respect to their physiological as well as biochemical traits, it is critical to be able to track the clones of species of interest after their characterization. Here, we compared the efficacy of five different genotyping methods for Spirodela polyrhiza, a species with very low intraspecific sequence variations, including polymorphic NB-ARC-related loci, tubulin-gene-based polymorphism (TBP), simple sequence repeat variations (SSR), multiplexed ISSR genotyping by sequencing (MIG-seq), and low-coverage, reduced-representation genome sequencing (GBS). Four of the five approaches could distinguish 20 to 22 genotypes out of the 23 investigated clones, while TBP resolved just seven genotypes. The choice for a particular method for intraspecific genotyping can depend on the research question and the project budget, while the combination of orthogonal methods may increase the confidence and resolution for the results obtained.

Chromosome Numbers and Genome Sizes of All 36 Duckweed Species (Lemnaceae)

Usually, chromosome sets (karyotypes) and genome sizes are rather stable for distinct species and therefore of diagnostic value for taxonomy. In combination with (cyto)genomics, both features provide essential cues for genome evolution and phylogenetic relationship studies within and between taxa above the species level. We present for the first time a survey on chromosome counts and genome size measurement for one or more accessions from all 36 duckweed species and discuss the evolutionary impact and peculiarities of both parameters in duckweeds.

Spruce versus Arabidopsis: different strategies of photosynthetic acclimation to light intensity change

The acclimation of higher plants to different light intensities is associated with a reorganization of the photosynthetic apparatus. These modifications, namely, changes in the amount of peripheral antenna (LHCII) of photosystem (PS) II and changes in PSII/PSI stoichiometry, typically lead to an altered chlorophyll (Chl) a/b ratio. However, our previous studies show that in spruce, this ratio is not affected by changes in growth light intensity. The evolutionary loss of PSII antenna proteins LHCB3 and LHCB6 in the Pinaceae family is another indication that the light acclimation strategy in spruce could be different. Here we show that, unlike Arabidopsis, spruce does not modify its PSII/PSI ratio and PSII antenna size to maximize its photosynthetic performance during light acclimation. Its large PSII antenna consists of many weakly bound LHCIIs, which form effective quenching centers, even at relatively low light. This, together with sensitive photosynthetic control on the level of cytochrome b₆f complex (protecting PSI), is the crucial photoprotective mechanism in spruce. High-light acclimation of spruce involves the disruption of PSII macro-organization, reduction of the amount of both PSII and PSI core complexes, synthesis of stress proteins that bind released Chls, and formation of "locked-in" quenching centers from uncoupled LHCIIs. Such response has been previously observed in the evergreen angiosperm Monstera deliciosa exposed to high light. We suggest that, in contrast to annuals, shade-tolerant evergreen land plants have their own strategy to cope with light intensity changes and the hallmark of this strategy is a stable Chl a/b ratio.

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.

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.

Chemoenzymatic labeling of DNA methylation patterns for single-molecule epigenetic mapping

DNA methylation, specifically, methylation of cytosine (C) nucleotides at the 5-carbon position (5-mC), is the most studied and significant epigenetic modification. Here we developed a chemoenzymatic procedure to fluorescently label non-methylated cytosines in CpG context, allowing epigenetic profiling of single DNA molecules spanning hundreds of thousands of base pairs. We used a CpG methyltransferase with a synthetic S-adenosyl-l-methionine cofactor analog to transfer an azide to cytosines instead of the natural methyl group. A fluorophore was then clicked onto the DNA, reporting on the amount and position of non-methylated CpGs. We found that labeling efficiency was increased up to 2-fold by the addition of a nucleosidase, presumably by degrading the inactive by-product of the cofactor after labeling, preventing its inhibitory effect. We used the method to determine the decline in global DNA methylation in a chronic lymphocytic leukemia patient and then performed whole-genome methylation mapping of the model plant Arabidopsis thaliana. Our genome maps show high concordance with published bisulfite sequencing methylation maps. Although mapping resolution is limited by optical detection to 500–1000 bp, the labeled DNA molecules produced by this approach are hundreds of thousands of base pairs long, allowing access to long repetitive and structurally variable genomic regions.

Biodiversity of Duckweed (Lemnaceae) in Water Reservoirs of Ukraine and China Assessed by Chloroplast DNA Barcoding

Monitoring and characterizing species biodiversity is essential for germplasm preservation, academic studies, and various practical applications. Duckweeds represent a group of tiny aquatic plants that include 36 species divided into 5 genera within the Lemnaceae family. They are an important part of aquatic ecosystems worldwide, often covering large portions of the water reservoirs they inhabit, and have many potential applications, including in bioremediation, biofuels, and biomanufacturing. Here, we evaluated the biodiversity of duckweeds in Ukraine and Eastern China by characterizing specimens using the two-barcode protocol with the chloroplast atpH–atpF and psbK–psbI spacer sequences. In total, 69 Chinese and Ukrainian duckweed specimens were sequenced. The sequences were compared against sequences in the NCBI database using BLAST. We identified six species from China (Spirodela polyrhiza, Landoltia punctata, Lemna aequinoctialis, Lemna minor, Lemna turionifera, and Wolffia globosa) and six from Ukraine (S. polyrhiza, Lemna gibba, Lemna minor, Lemna trisulca, Lemna turionifera, and Wolffia arrhiza). The most common duckweed species in the samples from Ukraine were Le. minor and S. polyrhiza, accounting for 17 and 15 out of 40 specimens, respectively. The most common duckweed species in the samples from China was S. polyrhiza, accounting for 15 out of 29 specimens. La. punctata and Le. aequinoctialis were also common in China, accounting for five and four specimens, respectively. According to both atpH–atpF and psbK–psbI barcode analyses, the species identified as Le. aequinoctialis does not form a uniform taxon similar to other duckweed species, and therefore the phylogenetic status of this species requires further clarification. By monitoring duckweeds using chloroplast DNA sequencing, we not only precisely identified local species and ecotypes, but also provided background for further exploration of native varieties with diverse genetic backgrounds. These data could be useful for future conservation, breeding, and biotechnological applications.

Wolffia, a minimalist plant and synthetic biology chassis

A highly simplified species for genome engineering would facilitate rational design of a synthetic plant. A candidate species is the aquatic, non-grass monocot wolffia (Wolffia australiana) in the Lemnaceae family. Commonly known as watermeal, wolffia is a rootless ball of several thousand cells the size of a pinhead and the fastest growing plant known on Earth. Its extreme morphological reduction is coupled to transposon-mediated streamlining of its transcriptome, which represents a core set of nonredundant protein coding genes. Despite its body plan and transcriptome being highly specialized for continuous growth, wolffia retains cell types relevant to higher plants. Systems level studies with this species could enable the creation of a defined biological chassis for synthetic plant construction.

Chromosome evolution and the genetic basis of agronomically important traits in greater yam

The nutrient-rich tubers of the greater yam, Dioscorea alata L., provide food and income security for millions of people around the world. Despite its global importance, however, greater yam remains an orphan crop. Here, we address this resource gap by presenting a highly contiguous chromosome-scale genome assembly of D. alata combined with a dense genetic map derived from African breeding populations. The genome sequence reveals an ancient allotetraploidization in the Dioscorea lineage, followed by extensive genome-wide reorganization. Using the genomic tools, we find quantitative trait loci for resistance to anthracnose, a damaging fungal pathogen of yam, and several tuber quality traits. Genomic analysis of breeding lines reveals both extensive inbreeding as well as regions of extensive heterozygosity that may represent interspecific introgression during domestication. These tools and insights will enable yam breeders to unlock the potential of this staple crop and take full advantage of its adaptability to varied environments.

Intraspecific Diversity in Aquatic Ecosystems: Comparison between Spirodela polyrhiza and Lemna minor in Natural Populations of Duckweed

Samples of two duckweed species, Spirodela polyrhiza and Lemna minor, were collected around small ponds and investigated concerning the question of whether natural populations of duckweeds constitute a single clone, or whether clonal diversity exists. Amplified fragment length polymorphism was used as a molecular method to distinguish clones of the same species. Possible intraspecific diversity was evaluated by average-linkage clustering. The main criterion to distinguish one clone from another was the 95% significance level of the Jaccard dissimilarity index for replicated samples. Within natural populations of L. minor, significant intraspecific genetic differences were detected. In each of the three small ponds harbouring populations of L. minor, based on twelve samples, between four and nine distinct clones were detected. Natural populations of L. minor consist of a mixture of several clones representing intraspecific biodiversity in an aquatic ecosystem. Moreover, identical distinct clones were discovered in more than one pond, located at a distance of 1 km and 2.4 km from each other. Evidently, fronds of L. minor were transported between these different ponds. The genetic differences for S. polyrhiza, however, were below the error-threshold of the method within a pond to detect distinct clones, but were pronounced between samples of two different ponds.

The Ribosomal DNA Loci of the Ancient Monocot Pistia stratiotes L. (Araceae) Contain Different Variants of the 35S and 5S Ribosomal RNA Gene Units

The freshwater plant water lettuce (Pistia stratiotes L.) grows in warm climatic zones and is used for phytoremediation and biomass production. P. stratiotes belongs to the Araceae, an ecologically and structurally diverse early monocot family, but the phylogenetic relationships among Araceae members are poorly understood. Ribosomal DNAs (rDNAs), including the 35S and 5S rDNA, encode the RNA components of ribosomes and are widely used in phylogenetic and evolutionary studies of various plant taxa. Here, we comprehensively characterized the chromosomal locations and molecular organization of 35S and 5S rDNA genes in water lettuce using karyological and molecular methods. Fluorescence in situ hybridization revealed a single location for the 35S and 5S rDNA loci, each on a different pair of the species' 28 chromosomes. Molecular cloning and nucleotide sequencing of 35S rDNA of P. stratiotes, the first representative Araceae sensu stricto in which such a study was performed, displayed typical structural characteristics. The full-length repeat showed high sequence conservation of the regions producing the 18S, 5.8S, and 25S rRNAs and divergence of the internal transcribed spacers ITS1 and ITS2 as well as the large intergenic spacer (IGS). Alignments of the deduced sequence of 18S rDNA with the sequences available for other Araceae and representatives of other clades were used for phylogenetic analysis. Examination of 11 IGS sequences revealed significant intra-genomic length variability due to variation in subrepeat number, with four types of units detected within the 35S rDNA locus of the P. stratiotes genome (estimated size 407 Mb/1C). Similarly, the 5S rDNA locus harbors gene units comprising a conserved 119-bp sequence encoding 5S rRNA and two types of non-transcribed spacer (NTS) sequences. Type I was classified into four subtypes, which apparently originated via progressive loss of subrepeats within the duplicated NTS region containing the 3' part of the 5S rRNA gene. The minor Type II NTS is shorter than Type I and differs in nucleotide composition. Some DNA clones containing two or three consecutive 5S rDNA repeats harbored 5S rDNA genes with different types of NTSs, confirming the mosaic composition of the 5S rDNA locus.

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.

Personal Perspectives on Plant Ribosomal RNA Genes Research: From Precursor-rRNA to Molecular Evolution

The history of rDNA research started almost 90 years ago when the geneticist, Barbara McClintock observed that in interphase nuclei of maize the nucleolus was formed in association with a specific region normally located near the end of a chromosome, which she called the nucleolar organizer region (NOR). Cytologists in the twentieth century recognized the nucleolus as a common structure in all eukaryotic cells, using both light and electron microscopy and biochemical and genetic studies identified ribosomes as the subcellular sites of protein synthesis. In the mid- to late 1960s, the synthesis of nuclear-encoded rRNA was the only system in multicellular organisms where transcripts of known function could be isolated, and their synthesis and processing could be studied. Cytogenetic observations of NOR regions with altered structure in plant interspecific hybrids and detailed knowledge of structure and function of rDNA were prerequisites for studies of nucleolar dominance, epistatic interactions of rDNA loci, and epigenetic silencing. In this article, we focus on the early rDNA research in plants, performed mainly at the dawn of molecular biology in the 60 to 80-ties of the last century which presented a prequel to the modern genomic era. We discuss - from a personal view - the topics such as synthesis of rRNA precursor (35S pre-rRNA in plants), processing, and the organization of 35S and 5S rDNA. Cloning and sequencing led to the observation that the transcribed and processed regions of the rRNA genes vary enormously, even between populations and species, in comparison with the more conserved regions coding for the mature rRNAs. Epigenetic phenomena and the impact of hybridization and allopolyploidy on rDNA expression and homogenization are discussed. This historical view of scientific progress and achievements sets the scene for the other articles highlighting the immense progress in rDNA research published in this special issue of Frontiers in Plant Science on "Molecular organization, evolution, and function of ribosomal DNA."

The Dynamics of NO3- and NH4+ Uptake in Duckweed Are Coordinated with the Expression of Major Nitrogen Assimilation Genes

Duckweed plants play important roles in aquatic ecosystems worldwide. They rapidly accumulate biomass and have potential uses in bioremediation of water polluted by fertilizer runoff or other chemicals. Here we studied the assimilation of two major sources of inorganic nitrogen, nitrate (NO3- ) and ammonium (NH4+), in six duckweed species: Spirodela polyrhiza, Landoltia punctata, Lemna aequinoctialis, Lemna turionifera, Lemna minor, and Wolffia globosa. All six duckweed species preferred NH4+ over NO3- and started using NO3- only when NH4+ was depleted. Using the available genome sequence, we analyzed the molecular structure and expression of eight key nitrogen assimilation genes in S. polyrhiza. The expression of genes encoding nitrate reductase and nitrite reductase increased about 10-fold when NO3- was supplied and decreased when NH4+ was supplied. NO3- and NH4+ induced the glutamine synthetase (GS) genes GS1;2 and the GS2 by 2- to 5-fold, respectively, but repressed GS1;1 and GS1;3. NH4+ and NO3- upregulated the genes encoding ferredoxin- and NADH-dependent glutamate synthases (Fd-GOGAT and NADH-GOGAT). A survey of nitrogen assimilation gene promoters suggested complex regulation, with major roles for NRE-like and GAATC/GATTC cis-elements, TATA-based enhancers, GA/CTn repeats, and G-quadruplex structures. These results will inform efforts to improve bioremediation and nitrogen use efficiency.

High-quality reference genome sequences of two coconut cultivars provide insights into evolution of monocot chromosomes and differentiation of fiber content and plant height

BACKGROUND

Coconut is an important tropical oil and fruit crop whose evolutionary position renders it a fantastic species for the investigation of the evolution of monocot chromosomes and the subsequent differentiation of ancient plants.

RESULTS

Here, we report the assembly and annotation of reference-grade genomes of Cn. tall and Cn. dwarf, whose genome sizes are 2.40 Gb and 2.39 Gb, respectively. The comparative analysis reveals that the two coconut subspecies diverge about 2-8 Mya while the conserved Arecaceae-specific whole-genome duplication (ω WGD) occurs approximately 47-53 Mya. It additionally allows us to reconstruct the ancestral karyotypes of the ten ancient monocot chromosomes and the evolutionary trajectories of the 16 modern coconut chromosomes. Fiber synthesis genes in Cn. tall, related to lignin and cellulose synthesis, are found at a higher copy number and expression level than dwarf coconuts. Integrated multi-omics analysis reveals that the difference in coconut plant height is the result of altered gibberellin metabolism, with both the GA20ox copy number and a single-nucleotide change in the promoter together leading to the difference in plant height between Cn. tall and Cn. dwarf.

CONCLUSION

We provide high-quality coconut genomes and reveal the genetic basis of trait differences between two coconuts through multi-omics analysis. We also reveal that the selection of plant height has been targeted for the same gene for millions of years, not only in natural selection of ancient plant as illustrated in coconut, but also for artificial selection in cultivated crops such as rice and maize.

Genomic analysis of the polyamine biosynthesis pathway in duckweed Spirodela polyrhiza L.: presence of the arginine decarboxylase pathway, absence of the ornithine decarboxylase pathway, and response to abiotic stresses

Identification of the polyamine biosynthetic pathway genes in duckweed S. polyrhiza reveals presence of prokaryotic as well as land plant-type ADC pathway but absence of ODC encoding genes. Their differential gene expression and transcript abundance is shown modulated by exogenous methyl jasmonate, salinity, and acidic pH. Genetic components encoding for polyamine (PA) biosynthetic pathway are known in several land plant species; however, little is known about them in aquatic plants. We utilized recently sequenced three duckweed (Spirodela polyrhiza) genome assemblies to map PA biosynthetic pathway genes in S. polyrhiza. PA biosynthesis in most higher plants except for Arabidopsis involves two pathways, via arginine decarboxylase (ADC) and ornithine decarboxylase (ODC). ADC-mediated PA biosynthetic pathway genes, namely, one arginase (SpARG1), two arginine decarboxylases (SpADC1, SpADC2), one agmatine iminohydrolase/deiminase (SpAIH), one N-carbamoyl putrescine amidase (SpCPA), three S-adenosylmethionine decarboxylases (SpSAMDc1, 2, 3), one spermidine synthase (SpSPDS1) and one spermine synthase (SpSPMS1) in S. polyrhiza genome were identified here. However, no locus was found for ODC pathway genes in this duckweed. Hidden Markov Model protein domain analysis established that SpADC1 is a prokaryotic/biodegradative type ADC and its molecular phylogenic classification fell in a separate prokaryotic origin ADC clade with SpADC2 as a biosynthetic type of arginine decarboxylase. However, thermospermine synthase (t-SPMS)/Aculis5 genes were not found present. Instead, one of the annotated SPDS may also function as SPMS, since it was found associated with the SPMS phylogenetic clade along with known SPMS genes. Moreover, we demonstrate that S. polyrhiza PA biosynthetic gene transcripts are differentially expressed in response to unfavorable conditions, such as exogenously added salt, methyl jasmonate, or acidic pH environment as well as in extreme temperature regimes. Thus, S. polyrhiza genome encodes for complete polyamine biosynthesis pathway and the genes are transcriptionally active in response to changing environmental conditions suggesting an important role of polyamines in this aquatic plant.

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.

Characterization of Frond and Flower Development and Identification of FT and FD Genes From Duckweed Lemna aequinoctialis Nd

Duckweeds (Araceae: Lemnoideae) are aquatic monocotyledonous plants that are characterized by their small size, rapid growth, and wide distribution. Developmental processes regulating the formation of their small leaf-like structures, called fronds, and tiny flowers are not well characterized. In many plant species, flowering is promoted by the florigen activation complex, whose major components are florigen FLOWERING LOCUS T (FT) protein and transcription factor FD protein. How this complex is regulated at the molecular level during duckweed flowering is also not well understood. In this study, we characterized the course of developmental changes during frond development and flower formation in Lemna aequinoctialis Nd, a short-day plant. Detailed observations of frond and flower development revealed that cell proliferation in the early stages of frond development is active as can be seen in the separate regions corresponding to two budding pouches in the proximal region of the mother frond. L. aequinoctialis produces two stamens of different lengths with the longer stamen growing more rapidly. Using high-throughput RNA sequencing (RNA-seq) and de novo assembly of transcripts from plants induced to flower, we identified the L. aequinoctialis FT and FD genes, whose products in other angiosperms form a transcriptional complex to promote flowering. We characterized the protein-protein interaction of duckweed FT and FD in yeast and examined the functions of the two gene products by overexpression in Arabidopsis. We found that L. aequinoctialis FTL1 promotes flowering, whereas FTL2 suppresses flowering.

Mosaic Arrangement of the 5S rDNA in the Aquatic Plant Landoltia punctata (Lemnaceae)

Duckweeds are a group of monocotyledonous aquatic plants in the Araceae superfamily, represented by 37 species divided into five genera. Duckweeds are the fastest growing flowering plants and are distributed around the globe; moreover, these plants have multiple applications, including biomass production, wastewater remediation, and making pharmaceutical proteins. Dotted duckweed (Landoltia punctata), the sole species in genus Landoltia, is one of the most resilient duckweed species. The ribosomal DNA (rDNA) encodes the RNA components of ribosomes and represents a significant part of plant genomes but has not been comprehensively studied in duckweeds. Here, we characterized the 5S rDNA genes in L. punctata by cloning and sequencing 25 PCR fragments containing the 5S rDNA repeats. No length variation was detected in the 5S rDNA gene sequence, whereas the nontranscribed spacer (NTS) varied from 151 to 524 bp. The NTS variants were grouped into two major classes, which differed both in nucleotide sequence and the type and arrangement of the spacer subrepeats. The dominant class I NTS, with a characteristic 12-bp TC-rich sequence present in 3-18 copies, was classified into four subclasses, whereas the minor class II NTS, with shorter, 9-bp nucleotide repeats, was represented by two identical sequences. In addition to these diverse subrepeats, class I and class II NTSs differed in their representation of cis-elements and the patterns of predicted G-quadruplex structures, which may influence the transcription of the 5S rDNA. Similar to related duckweed species in the genus Spirodela, L. punctata has a relatively low rDNA copy number, but in contrast to Spirodela and the majority of other plants, the arrangement of the 5S rDNA units demonstrated an unusual, heterogeneous pattern in L. punctata, as revealed by analyzing clones containing double 5S rDNA neighboring units. Our findings may further stimulate the research on the evolution of the plant rDNA and discussion of the molecular forces driving homogenization of rDNA repeats in concerted evolution.

In silico analysis of glycosyltransferase 2 family genes in duckweed (Spirodela polyrhiza) and its role in salt stress tolerance

Plant glycosyltransferase 2 (GT2) family genes are involved in plant abiotic stress tolerance. However, the roles of GT2 genes in the abiotic resistance in freshwater plants are largely unknown. We identified seven GT2 genes in duckweed, remarkably more than those in the genomes of Arabidopsis thaliana, Oryza sativa, Amborella trichopoda, Nymphaea tetragona, Persea americana, Zostera marina, and Ginkgo biloba, suggesting a significant expansion of this family in the duckweed genome. Phylogeny resolved the GT2 family into two major clades. Six duckweed genes formed an independent subclade in Clade I, and the other was clustered in Clade II. Gene structure and protein domain analysis showed that the lengths of the seven duckweed GT2 genes were varied, and the majority of GT2 genes harbored two conserved domains, PF04722.12 and PF00535.25. The expression of all Clade I duckweed GT2 genes was elevated at 0 h after salt treatment, suggesting a common role of these genes in rapid response to salt stress. The gene Sp01g00794 was highly expressed at 12 and 24 h after salt treatment, indicating its association with salt stress resilience. Overall, these results are essential for studies on the molecular mechanisms in stress response and resistance in aquatic plants.

Genome-Wide Identification of the Nramp Gene Family in Spirodela polyrhiza and Expression Analysis under Cadmium Stress

Natural resistance-associated macrophage proteins (Nramps) are specific metal transporters in plants with different functions among various species. The evolutionary and functional information of the Nramp gene family in Spirodela polyrhiza has not been previously reported in detail. To identify the Nramp genes in S. polyrhiza, we performed genome-wide identification, characterization, classification, and cis-elements analysis among 22 species with 138 amino acid sequences. We also conducted chromosomal localization and analyzed the synteny relationship, promoter, subcellular localization, and expression patterns in S. polyrhiza. β-Glucuronidase staining indicated that SpNramp1 and SpNramp3 mainly accumulated in the root and joint between mother and daughter frond. Moreover, SpNramp1 was also widely displayed in the frond. SpNramp2 was intensively distributed in the root and frond. Quantitative real-time PCR results proved that the SpNramp gene expression level was influenced by Cd stress, especially in response to Fe or Mn deficiency. The study provides detailed information on the SpNramp gene family and their distribution and expression, laying a beneficial foundation for functional research.

Improved Spirodela polyrhiza genome and proteomic analyses reveal a conserved chromosomal structure with high abundance of chloroplastic proteins favoring energy production

Duckweeds are a monophyletic group of rapidly reproducing aquatic monocots in the Lemnaceae family. Given their clonal, exponentially fast reproduction, a key question is whether genome structure is conserved across the species in the absence of meiotic recombination. Here, we studied the genome and proteome of Spirodela polyrhiza, or greater duckweed, which has the largest body plan yet the smallest genome size in the family (1C=150 Mb). Using Oxford Nanopore sequencing combined with Hi-C scaffolding, we generated a highly contiguous, chromosome-scale assembly of S. polyrhiza line Sp7498 (Sp7498_HiC). Both the Sp7498_HiC and Sp9509 genome assemblies reveal large chromosomal misorientations relative to a recent PacBio assembly of Sp7498, highlighting the need for orthogonal long-range scaffolding techniques such as Hi-C and BioNano optical mapping. Shotgun proteomics of Sp7498 verified the expression of ~2250 proteins and revealed a high abundance of proteins involved in photosynthesis and carbohydrate metabolism among other functions. In addition, a strong increase in chloroplast proteins was observed that correlated to chloroplast density. This Sp7498_HiC genome was generated cheaply and quickly with a single Oxford Nanopore MinION flow cell and one Hi-C library in a classroom setting. Combining these data with a mass spectrometry-generated proteome illustrates the utility of duckweed as a model for genomics- and proteomics-based education.

Flowering and Seed Production across the Lemnaceae

Plants in the family Lemnaceae are aquatic monocots and the smallest, simplest, and fastest growing angiosperms. Their small size, the smallest family member is 0.5 mm and the largest is 2.0 cm, as well as their diverse morphologies make these plants ideal for laboratory studies. Their rapid growth rate is partially due to the family's neotenous lifestyle, where instead of maturing and producing flowers, the plants remain in a juvenile state and continuously bud asexually. Maturation and flowering in the wild are rare in most family members. To promote further research on these unique plants, we have optimized laboratory flowering protocols for 3 of the 5 genera: Spirodela; Lemna; and Wolffia in the Lemnaceae. Duckweeds were widely used in the past for research on flowering, hormone and amino acid biosynthesis, the photosynthetic apparatus, and phytoremediation due to their aqueous lifestyle and ease of aseptic culture. There is a recent renaissance in interest in growing these plants as non-lignified biomass sources for fuel production, and as a resource-efficient complete protein source. The genome sequences of several Lemnaceae family members have become available, providing a foundation for genetic improvement of these plants as crops. The protocols for maximizing flowering described herein are based on screens testing daylength, a variety of media, supplementation with salicylic acid or ethylenediamine-N,N'-bis(2-hydroxyphenylacetic acid) (EDDHA), as well as various culture vessels for effects on flowering of verified Lemnaceae strains available from the Rutgers Duckweed Stock Cooperative.

Limitation of current probe design for oligo-cross-FISH, exemplified by chromosome evolution studies in duckweeds

Duckweeds represent a small, free-floating aquatic family (Lemnaceae) of the monocot order Alismatales with the fastest growth rate among flowering plants. They comprise five genera (Spirodela, Landoltia, Lemna, Wolffiella, and Wolffia) varying in genome size and chromosome number. Spirodela polyrhiza had the first sequenced duckweed genome. Cytogenetic maps are available for both species of the genus Spirodela (S. polyrhiza and S. intermedia). However, elucidation of chromosome homeology and evolutionary chromosome rearrangements by cross-FISH using Spirodela BAC probes to species of other duckweed genera has not been successful so far. We investigated the potential of chromosome-specific oligo-FISH probes to address these topics. We designed oligo-FISH probes specific for one S. intermedia and one S. polyrhiza chromosome (Fig. 1a). Our results show that these oligo-probes cross-hybridize with the homeologous regions of the other congeneric species, but are not suitable to uncover chromosomal homeology across duckweeds genera. This is most likely due to too low sequence similarity between the investigated genera and/or too low probe density on the target genomes. Finally, we suggest genus-specific design of oligo-probes to elucidate chromosome evolution across duckweed genera.

Genome and time-of-day transcriptome of Wolffia australiana link morphological minimization with gene loss and less growth control

Rootless plants in the genus Wolffia are some of the fastest growing known plants on Earth. Wolffia have a reduced body plan, primarily multiplying through a budding type of asexual reproduction. Here, we generated draft reference genomes for Wolffia australiana (Benth.) Hartog & Plas, which has the smallest genome size in the genus at 357 Mb and has a reduced set of predicted protein-coding genes at about 15,000. Comparison between multiple high-quality draft genome sequences from W. australiana clones confirmed loss of several hundred genes that are highly conserved among flowering plants, including genes involved in root developmental and light signaling pathways. Wolffia has also lost most of the conserved nucleotide-binding leucine-rich repeat (NLR) genes that are known to be involved in innate immunity, as well as those involved in terpene biosynthesis, while having a significant overrepresentation of genes in the sphingolipid pathways that may signify an alternative defense system. Diurnal expression analysis revealed that only 13% of Wolffia genes are expressed in a time-of-day (TOD) fashion, which is less than the typical ∼40% found in several model plants under the same condition. In contrast to the model plants Arabidopsis and rice, many of the pathways associated with multicellular and developmental processes are not under TOD control in W. australiana, where genes that cycle the conditions tested predominantly have carbon processing and chloroplast-related functions. The Wolffia genome and TOD expression data set thus provide insight into the interplay between a streamlined plant body plan and optimized growth.

Research Progress of a Potential Bioreactor: Duckweed

Recently, plant bioreactors have flourished into an exciting area of synthetic biology because of their product safety, inexpensive production cost, and easy scale-up. Duckweed is the smallest and fastest-growing aquatic plant, and has advantages including simple processing and the ability to grow high biomass in smaller areas. Therefore, duckweed could be used as a new potential bioreactor for biological products such as vaccines, antibodies, pharmaceutical proteins, and industrial enzymes. Duckweed has made a breakthrough in biosynthesis as a chassis plant and is being utilized for the production of plenty of biological products or bio-derivatives with multiple uses and high values. This review summarizes the latest progress on genetic background, genetic transformation system, and bioreactor development of duckweed, and provides insights for further exploration and application of duckweed.

Duckweeds: their utilization, metabolites and cultivation

Duckweeds are floating plants of the family Lemnaceae, comprising 5 genera and 36 species. They typically live in ponds or lakes and are found worldwide, except the polar regions. There are two duckweed subfamilies-namely Lemnoidea and Wolffioideae, with 15 and 21 species, respectively. Additionally, they have characteristic reproduction methods. Several metabolites have also been reported in various duckweeds. Duckweeds have a wide range of adaptive capabilities and are particularly suitable for experiments requiring high productivity because of their speedy growth and reproduction rates. Duckweeds have been studied for their use as food/feed resources and pharmaceuticals, as well as for phytoremediation and industrial applications. Because there are numerous duckweed species, culture conditions should be optimized for industrial applications. Here, we review and summarize studies on duckweed species and their utilization, metabolites, and cultivation methods to support the extended application of duckweeds in future.

Identification, Phylogeny, and Comparative Expression of the Lipoxygenase Gene Family of the Aquatic Duckweed, Spirodela polyrhiza, during Growth and in Response to Methyl Jasmonate and Salt

Lipoxygenases (LOXs) (EC 1.13.11.12) catalyze the oxygenation of fatty acids and produce oxylipins, including the plant hormone jasmonic acid (JA) and its methyl ester, methyl jasmonate (MeJA). Little information is available about the LOX gene family in aquatic plants. We identified a novel LOX gene family comprising nine LOX genes in the aquatic plant Spirodela polyrhiza (greater duckweed). The reduced anatomy of S. polyrhiza did not lead to a reduction in LOX family genes. The 13-LOX subfamily, with seven genes, predominates, while the 9-LOX subfamily is reduced to two genes, an opposite trend from known LOX families of other plant species. As the 13-LOX subfamily is associated with the synthesis of JA/MeJA, its predominance in the Spirodela genome raises the possibility of a higher requirement for the hormone in the aquatic plant. JA-/MeJA-based feedback regulation during culture aging as well as the induction of LOX gene family members within 6 h of salt exposure are demonstrated.

Chromosome-scale genome assembly for the duckweed Spirodela intermedia, integrating cytogenetic maps, PacBio and Oxford Nanopore libraries

Duckweeds are small, free-floating, morphologically highly reduced organisms belonging to the monocot order Alismatales. They display the most rapid growth among flowering plants, vary ~ 14-fold in genome size and comprise five genera. Spirodela is the phylogenetically oldest genus with only two mainly asexually propagating species: S. polyrhiza (2n = 40; 160 Mbp/1C) and S. intermedia (2n = 36; 160 Mbp/1C). This study combined comparative cytogenetics and de novo genome assembly based on PacBio, Illumina and Oxford Nanopore (ON) reads to obtain the first genome reference for S. intermedia and to compare its genomic features with those of the sister species S. polyrhiza. Both species' genomes revealed little more than 20,000 putative protein-coding genes, very low rDNA copy numbers and a low amount of repetitive sequences, mainly Ty3/gypsy retroelements. The detection of a few new small chromosome rearrangements between both Spirodela species refined the karyotype and the chromosomal sequence assignment for S. intermedia.

Estimation of the SNP Mutation Rate in Two Vegetatively Propagating Species of Duckweed

Mutation rate estimates for vegetatively reproducing organisms are rare, despite their frequent occurrence across the tree of life. Here we report mutation rate estimates in two vegetatively reproducing duckweed species, Lemna minor and Spirodela polyrhiza We use a modified approach to estimating mutation rates by taking into account the reduction in mutation detection power that occurs when new individuals are produced from multiple cell lineages. We estimate an extremely low per generation mutation rate in both species of duckweed and note that allelic coverage at de novo mutation sites is very skewed. We also find no substantial difference in mutation rate between mutation accumulation lines propagated under benign conditions and those grown under salt stress. Finally, we discuss the implications of interpreting mutation rate estimates in vegetatively propagating organisms.

Taro Genome Assembly and Linkage Map Reveal QTLs for Resistance to Taro Leaf Blight

Taro (Colocasia esculenta) is a food staple widely cultivated in the humid tropics of Asia, Africa, Pacific and the Caribbean. One of the greatest threats to taro production is Taro Leaf Blight caused by the oomycete pathogen Phytophthora colocasiae Here we describe a de novo taro genome assembly and use it to analyze sequence data from a Taro Leaf Blight resistant mapping population. The genome was assembled from linked-read sequences (10x Genomics; ∼60x coverage) and gap-filled and scaffolded with contigs assembled from Oxford Nanopore Technology long-reads and linkage map results. The haploid assembly was 2.45 Gb total, with a maximum contig length of 38 Mb and scaffold N50 of 317,420 bp. A comparison of family-level (Araceae) genome features reveals the repeat content of taro to be 82%, >3.5x greater than in great duckweed (Spirodela polyrhiza), 23%. Both genomes recovered a similar percent of Benchmarking Universal Single-copy Orthologs, 80% and 84%, based on a 3,236 gene database for monocot plants. A greater number of nucleotide-binding leucine-rich repeat disease resistance genes were present in genomes of taro than the duckweed, ∼391 vs. ∼70 (∼182 and ∼46 complete). The mapping population data revealed 16 major linkage groups with 520 markers, and 10 quantitative trait loci (QTL) significantly associated with Taro Leaf Blight disease resistance. The genome sequence of taro enhances our understanding of resistance to TLB, and provides markers that may accelerate breeding programs. This genome project may provide a template for developing genomic resources in other understudied plant species.

Intraspecific variations in cadmium tolerance and phytoaccumulation in giant duckweed (Spirodela polyrhiza)

Duckweeds are widely recognized for the heavy metal phytoremediation. However, the intraspecific variations in biological responses of duckweeds to heavy metal remain largely unknown. Here, the toxicity and phytoaccumulation of cadmium (Cd) were synchronously evaluated in 30 accessions of giant duckweed (Spirodela polyrhiza) collected from different provenances in Southern China. Exposure to 1 μM Cd decreased relative growth rates of dry weight, fronds number and fronds area, as well as photosynthetic pigment contents, while it increased H₂O₂ accumulation, lipid peroxidation and activities of anti-oxidant enzymes in the majority of accessions. Cd treatment led to remarkable Cd accumulation but little changes in the starch content in giant duckweed. The biological responses to Cd varied among the accessions. Further correlation analysis indicated that growth traits and Cd concentration were positively correlated with Cd accumulation, while the contents of chlorophyll, H₂O₂ and MDA were negatively associated with Cd accumulation. Our results proved the great intraspecific variation in Cd tolerance of giant duckweed, suggesting a valuable natural resource for Cd phytoremediation. Moreover, different mechanisms may be exploited by S. polyrhiza for phytoaccumulation, but growth maintenance, Cd uptake and antioxidative enzyme-independent ROS-scavenging under Cd exposure are the common mechanisms contributing to Cd accumulation ability.