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.

Retrotransposon addiction promotes centromere function via epigenetically activated small RNAs.

Retrotransposons have invaded eukaryotic centromeres in cycles of repeat expansion and purging, but the function of centromeric retrotransposons, if any, has remained unclear. In Arabidopsis, centromeric ATHILA retrotransposons give rise to epigenetically activated short interfering RNAs (easiRNAs) in mutants in DECREASE IN DNA METHYLATION1 (DDM1), which promote histone H3 lysine-9 di-methylation (H3K9me2). Here, we show that mutants which lose both DDM1 and RNA dependent RNA polymerase (RdRP) have pleiotropic developmental defects and mis-segregation of chromosome 5 during mitosis. Fertility defects are epigenetically inherited with the centromeric region of chromosome 5, and can be rescued by directing artificial small RNAs to a single family of ATHILA5 retrotransposons specifically embedded within this centromeric region. easiRNAs and H3K9me2 promote pericentromeric condensation, chromosome cohesion and proper chromosome segregation in mitosis. Insertion of ATHILA silences transcription, while simultaneously making centromere function dependent on retrotransposon small RNAs, promoting the selfish survival and spread of centromeric retrotransposons. Parallels are made with the fission yeast S. pombe, where chromosome segregation depends on RNAi, and with humans, where chromosome segregation depends on both RNAi and HELLSDDM1.

Teosinte Pollen Drive guides maize diversification and domestication by RNAi.

Meiotic drivers subvert Mendelian expectations by manipulating reproductive development to bias their own transmission. Chromosomal drive typically functions in asymmetric female meiosis, while gene drive is normally postmeiotic and typically found in males. Using single molecule and single-pollen genome sequencing, we describe Teosinte Pollen Drive, an instance of gene drive in hybrids between maize (Zea mays ssp. mays) and teosinte mexicana (Zea mays ssp. mexicana), that depends on RNA interference (RNAi). 22nt small RNAs from a non-coding RNA hairpin in mexicana depend on Dicer-Like 2 (Dcl2) and target Teosinte Drive Responder 1 (Tdr1), which encodes a lipase required for pollen viability. Dcl2, Tdr1, and the hairpin are in tight pseudolinkage on chromosome 5, but only when transmitted through the male. Introgression of mexicana into early cultivated maize is thought to have been critical to its geographical dispersal throughout the Americas, and a tightly linked inversion in mexicana spans a major domestication sweep in modern maize. A survey of maize landraces and sympatric populations of teosinte mexicana reveals correlated patterns of admixture among unlinked genes required for RNAi on at least 4 chromosomes that are also subject to gene drive in pollen from synthetic hybrids. Teosinte Pollen Drive likely played a major role in maize domestication and diversification, and offers an explanation for the widespread abundance of "self" small RNAs in the germlines of plants and animals.

Duckweeds.

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

Engineering triacylglycerol accumulation in duckweed (Lemna japonica).

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

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.

Arabidopsis retrotransposon virus-like particles and their regulation by epigenetically activated small RNA.

In Arabidopsis, LTR retrotransposons are activated by mutations in the chromatin gene DECREASE in DNA METHYLATION 1 (DDM1), giving rise to 21- to 22-nt epigenetically activated siRNA (easiRNA) that depend on RNA DEPENDENT RNA POLYMERASE 6 (RDR6). We purified virus-like particles (VLPs) from ddm1 and ddm1rdr6 mutants in which genomic RNA is reverse transcribed into complementary DNA. High-throughput short-read and long-read sequencing of VLP DNA (VLP DNA-seq) revealed a comprehensive catalog of active LTR retrotransposons without the need for mapping transposition, as well as independent of genomic copy number. Linear replication intermediates of the functionally intact COPIA element EVADE revealed multiple central polypurine tracts (cPPTs), a feature shared with HIV in which cPPTs promote nuclear localization. For one member of the ATCOPIA52 subfamily (SISYPHUS), cPPT intermediates were not observed, but abundant circular DNA indicated transposon "suicide" by auto-integration within the VLP. easiRNA targeted EVADE genomic RNA, polysome association of GYPSY (ATHILA) subgenomic RNA, and transcription via histone H3 lysine-9 dimethylation. VLP DNA-seq provides a comprehensive landscape of LTR retrotransposons and their control at transcriptional, post-transcriptional, and reverse transcriptional levels.

Epigenetic activation of meiotic recombination near Arabidopsis thaliana centromeres via loss of H3K9me2 and non-CG DNA methylation.

Eukaryotic centromeres contain the kinetochore, which connects chromosomes to the spindle allowing segregation. During meiosis, centromeres are suppressed for inter-homolog crossover, as recombination in these regions can cause chromosome missegregation and aneuploidy. Plant centromeres are surrounded by transposon-dense pericentromeric heterochromatin that is epigenetically silenced by histone 3 lysine 9 dimethylation (H3K9me2), and DNA methylation in CG and non-CG sequence contexts. However, the role of these chromatin modifications in control of meiotic recombination in the pericentromeres is not fully understood. Here, we show that disruption of Arabidopsis thaliana H3K9me2 and non-CG DNA methylation pathways, for example, via mutation of the H3K9 methyltransferase genes KYP/SUVH4 SUVH5 SUVH6, or the CHG DNA methyltransferase gene CMT3, increases meiotic recombination in proximity to the centromeres. Using immunocytological detection of MLH1 foci and genotyping by sequencing of recombinant plants, we observe that H3K9me2 and non-CG DNA methylation pathway mutants show increased pericentromeric crossovers. Increased pericentromeric recombination in H3K9me2/non-CG mutants occurs in hybrid and inbred backgrounds and likely involves contributions from both the interfering and noninterfering crossover repair pathways. We also show that meiotic DNA double-strand breaks (DSBs) increase in H3K9me2/non-CG mutants within the pericentromeres, via purification and sequencing of SPO11-1-oligonucleotides. Therefore, H3K9me2 and non-CG DNA methylation exert a repressive effect on both meiotic DSB and crossover formation in plant pericentromeric heterochromatin. Our results may account for selection of enhancer trap Dissociation (Ds) transposons into the CMT3 gene by recombination with proximal transposon launch-pads.

The Conserved RNA Binding Cyclophilin, Rct1, Regulates Small RNA Biogenesis and Splicing Independent of Heterochromatin Assembly.

RNAi factors and their catalytic activities are essential for heterochromatin assembly in S. pombe. This has led to the idea that siRNAs can promote H3K9 methylation by recruiting the cryptic loci regulator complex (CLRC), also known as recombination in K complex (RIKC), to the nucleation site. The conserved RNA-binding protein Rct1 (AtCyp59/SIG-7) interacts with splicing factors and RNA polymerase II. Here we show that Rct1 promotes processing of pericentromeric transcripts into siRNAs via the RNA recognition motif. Surprisingly, loss of siRNA in rct1 mutants has no effect on H3K9 di- or tri-methylation, resembling other splicing mutants, suggesting that post-transcriptional gene silencing per se is not required to maintain heterochromatin. Splicing of the Argonaute gene is also defective in rct1 mutants and contributes to loss of silencing but not to loss of siRNA. Our results suggest that Rct1 guides transcripts to the RNAi machinery by promoting splicing of elongating non-coding transcripts.

Survey of the total fatty acid and triacylglycerol composition and content of 30 duckweed species and cloning of a Δ6-desaturase responsible for the production of γ-linolenic and stearidonic acids in Lemna gibba.

BACKGROUND: Duckweeds, i.e., members of the Lemnoideae family, are amongst the smallest aquatic flowering plants. Their high growth rate, aquatic habit and suitability for bio-remediation make them strong candidates for biomass production. Duckweeds have been studied for their potential as feedstocks for bioethanol production; however, less is known about their ability to accumulate reduced carbon as fatty acids (FA) and oil. RESULTS: Total FA profiles of thirty duckweed species were analysed to assess the natural diversity within the Lemnoideae. Total FA content varied between 4.6% and 14.2% of dry weight whereas triacylglycerol (TAG) levels varied between 0.02% and 0.15% of dry weight. Three FA, 16:0 (palmitic), 18:2Δ9,12 (Linoleic acid, or LN) and 18:3Δ9,12,15 (α-linolenic acid, or ALA) comprise more than 80% of total duckweed FA. Seven Lemna and two Wolffiela species also accumulate polyunsaturated FA containing Δ6-double bonds, i.e., GLA and SDA. Relative to total FA, TAG is enriched in saturated FA and deficient in polyunsaturated FA, and only five Lemna species accumulate Δ6-FA in their TAG. A putative Δ6-desaturase designated LgDes, with homology to a family of front-end Δ6-FA and Δ8-spingolipid desaturases, was identified in the assembled DNA sequence of Lemna gibba. Expression of a synthetic LgDes gene in Nicotiana benthamiana resulted in the accumulation of GLA and SDA, confirming it specifies a Δ6-desaturase. CONCLUSIONS: Total accumulation of FA varies three-fold across the 30 species of Lemnoideae surveyed. Nine species contain GLA and SDA which are synthesized by a Δ6 front-end desaturase, but FA composition is otherwise similar. TAG accumulates up to 0.15% of total dry weight, comparable to levels found in the leaves of terrestrial plants. Polyunsaturated FA is underrepresented in TAG, and the Δ6-FA GLA and SDA are found in the TAG of only five of the nine Lemna species that produce them. When present, GLA is enriched and SDA diminished relative to their abundance in the total FA pool.