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.

The genome of the Wollemi pine, a critically endangered "living fossil" unchanged since the Cretaceous, reveals extensive ancient transposon activity.

We present the genome of the living fossil, Wollemia nobilis, a southern hemisphere conifer morphologically unchanged since the Cretaceous. Presumed extinct until rediscovery in 1994, the Wollemi pine is critically endangered with less than 60 wild adults threatened by intensifying bushfires in the Blue Mountains of Australia. The 12 Gb genome is among the most contiguous large plant genomes assembled, with extremely low heterozygosity and unusual abundance of DNA transposons. Reduced representation and genome re-sequencing of individuals confirms a relictual population since the last major glacial/drying period in Australia, 120 ky BP. Small RNA and methylome sequencing reveal conservation of ancient silencing mechanisms despite the presence of thousands of active and abundant transposons, including some transferred horizontally to conifers from arthropods in the Jurassic. A retrotransposon burst 8-6 my BP coincided with population decline, possibly as an adaptation enhancing epigenetic diversity. Wollemia, like other conifers, is susceptible to Phytophthora, and a suite of defense genes, similar to those in loblolly pine, are targeted for silencing by sRNAs in leaves. The genome provides insight into the earliest seed plants, while enabling conservation efforts.

Chromatin remodeling of histone H3 variants by DDM1 underlies epigenetic inheritance of DNA methylation.

Nucleosomes block access to DNA methyltransferase, unless they are remodeled by DECREASE in DNA METHYLATION 1 (DDM1LSH/HELLS), a Snf2-like master regulator of epigenetic inheritance. We show that DDM1 promotes replacement of histone variant H3.3 by H3.1. In ddm1 mutants, DNA methylation is partly restored by loss of the H3.3 chaperone HIRA, while the H3.1 chaperone CAF-1 becomes essential. The single-particle cryo-EM structure at 3.2 Å of DDM1 with a variant nucleosome reveals engagement with histone H3.3 near residues required for assembly and with the unmodified H4 tail. An N-terminal autoinhibitory domain inhibits activity, while a disulfide bond in the helicase domain supports activity. DDM1 co-localizes with H3.1 and H3.3 during the cell cycle, and with the DNA methyltransferase MET1Dnmt1, but is blocked by H4K16 acetylation. The male germline H3.3 variant MGH3/HTR10 is resistant to remodeling by DDM1 and acts as a placeholder nucleosome in sperm cells for epigenetic inheritance.

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.

Pseudouridine guides germline small RNA transport and epigenetic inheritance.

Epigenetic modifications that arise during plant and animal development, such as DNA and histone modification, are mostly reset during gamete formation, but some are inherited from the germline including those marking imprinted genes1. Small RNAs guide these epigenetic modifications, and some are also inherited by the next generation2,3. In C. elegans, these inherited small RNAs have poly (UG) tails4, but how inherited small RNAs are distinguished in other animals and plants is unknown. Pseudouridine (Ψ) is the most abundant RNA modification but has not been explored in small RNAs. Here, we develop novel assays to detect Ψ in short RNA sequences, demonstrating its presence in mouse and Arabidopsis microRNAs and their precursors. We also detect substantial enrichment in germline small RNAs, namely epigenetically activated siRNAs (easiRNAs) in Arabidopsis pollen, and piwi-interacting piRNAs in mouse testis. In pollen, pseudouridylated easiRNAs are localized to sperm cells, and we found that PAUSED/HEN5 (PSD), the plant homolog of Exportin-t, interacts genetically with Ψ and is required for transport of easiRNAs into sperm cells from the vegetative nucleus. We further show that Exportin-t is required for the triploid block: chromosome dosage-dependent seed lethality that is epigenetically inherited from pollen. Thus, Ψ has a conserved role in marking inherited small RNAs in the germline.

Chromatin remodeling of histone H3 variants underlies epigenetic inheritance of DNA methylation.

Epigenetic inheritance refers to the faithful replication of DNA methylation and histone modification independent of DNA sequence. Nucleosomes block access to DNA methyltransferases, unless they are remodeled by DECREASE IN DNA METHYLATION1 (DDM1 Lsh/HELLS ), a Snf2-like master regulator of epigenetic inheritance. We show that DDM1 activity results in replacement of the transcriptional histone variant H3.3 for the replicative variant H3.1 during the cell cycle. In ddm1 mutants, DNA methylation can be restored by loss of the H3.3 chaperone HIRA, while the H3.1 chaperone CAF-1 becomes essential. The single-particle cryo-EM structure at 3.2 Å of DDM1 with a variant nucleosome reveals direct engagement at SHL2 with histone H3.3 at or near variant residues required for assembly, as well as with the deacetylated H4 tail. An N-terminal autoinhibitory domain binds H2A variants to allow remodeling, while a disulfide bond in the helicase domain is essential for activity in vivo and in vitro . We show that differential remodeling of H3 and H2A variants in vitro reflects preferential deposition in vivo . DDM1 co-localizes with H3.1 and H3.3 during the cell cycle, and with the DNA methyltransferase MET1 Dnmt1 . DDM1 localization to the chromosome is blocked by H4K16 acetylation, which accumulates at DDM1 targets in ddm1 mutants, as does the sperm cell specific H3.3 variant MGH3 in pollen, which acts as a placeholder nucleosome in the germline and contributes to epigenetic inheritance.

Oncogene-like addiction to aneuploidy in human cancers.

Most cancers exhibit aneuploidy, but its functional significance in tumor development is controversial. Here, we describe ReDACT (Restoring Disomy in Aneuploid cells using CRISPR Targeting), a set of chromosome engineering tools that allow us to eliminate specific aneuploidies from cancer genomes. Using ReDACT, we created a panel of isogenic cells that have or lack common aneuploidies, and we demonstrate that trisomy of chromosome 1q is required for malignant growth in cancers harboring this alteration. Mechanistically, gaining chromosome 1q increases the expression of MDM4 and suppresses p53 signaling, and we show that TP53 mutations are mutually exclusive with 1q aneuploidy in human cancers. Thus, tumor cells can be dependent on specific aneuploidies, raising the possibility that these "aneuploidy addictions" could be targeted as a therapeutic strategy.

Oncogene-like addiction to aneuploidy in human cancers.

Most cancers exhibit aneuploidy, but its functional significance in tumor development is controversial. Here, we describe ReDACT (Restoring Disomy in Aneuploid cells using CRISPR Targeting), a set of chromosome engineering tools that allow us to eliminate specific aneuploidies from cancer genomes. Using ReDACT, we created a panel of isogenic cells that have or lack common aneuploidies, and we demonstrate that trisomy of chromosome 1q is required for malignant growth in cancers harboring this alteration. Mechanistically, gaining chromosome 1q increases the expression of MDM4 and suppresses TP53 signaling, and we show that TP53 mutations are mutually-exclusive with 1q aneuploidy in human cancers. Thus, specific aneuploidies play essential roles in tumorigenesis, raising the possibility that targeting these "aneuploidy addictions" could represent a novel approach for cancer treatment.

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.

AGO104 is a RdDM effector of paramutation at the maize b1 locus.

Although paramutation has been well-studied at a few hallmark loci involved in anthocyanin biosynthesis in maize, the cellular and molecular mechanisms underlying the phenomenon remain largely unknown. Previously described actors of paramutation encode components of the RNA-directed DNA-methylation (RdDM) pathway that participate in the biogenesis of 24-nucleotide small interfering RNAs (24-nt siRNAs) and long non-coding RNAs. In this study, we uncover an ARGONAUTE (AGO) protein as an effector of the RdDM pathway that is in charge of guiding 24-nt siRNAs to their DNA target to create de novo DNA methylation. We combined immunoprecipitation, small RNA sequencing and reverse genetics to, first, validate AGO104 as a member of the RdDM effector complex and, then, investigate its role in paramutation. We found that AGO104 binds 24-nt siRNAs involved in RdDM, including those required for paramutation at the b1 locus. We also show that the ago104-5 mutation causes a partial reversion of the paramutation phenotype at the b1 locus, revealed by intermediate pigmentation levels in stem tissues. Therefore, our results place AGO104 as a new member of the RdDM effector complex that plays a role in paramutation at the b1 locus in maize.

The genetic and epigenetic landscape of the Arabidopsis centromeres.

Centromeres attach chromosomes to spindle microtubules during cell division and, despite this conserved role, show paradoxically rapid evolution and are typified by complex repeats. We used long-read sequencing to generate the Col-CEN Arabidopsis thaliana genome assembly that resolves all five centromeres. The centromeres consist of megabase-scale tandemly repeated satellite arrays, which support CENTROMERE SPECIFIC HISTONE H3 (CENH3) occupancy and are densely DNA methylated, with satellite variants private to each chromosome. CENH3 preferentially occupies satellites that show the least amount of divergence and occur in higher-order repeats. The centromeres are invaded by ATHILA retrotransposons, which disrupt genetic and epigenetic organization. Centromeric crossover recombination is suppressed, yet low levels of meiotic DNA double-strand breaks occur that are regulated by DNA methylation. We propose that Arabidopsis centromeres are evolving through cycles of satellite homogenization and retrotransposon-driven diversification.

Regulation of retrotransposition in Arabidopsis.

Plant genomes are largely comprised of retrotransposons which can replicate through 'copy and paste' mechanisms. Long terminal repeat (LTR) retrotransposons are the major class of retrotransposons in plant species, and importantly they broadly affect the expression of nearby genes. Although most LTR retrotransposons are non-functional, active retrotranspositions have been reported in plant species or mutants under normal growth condition and environmental stresses. With the well-defined reference genome and numerous mutant alleles, Arabidopsis studies have significantly expanded our understanding of retrotransposon regulation. Active LTR retrotransposon loci produce virus-like particles to perform reverse transcription, and their complementary DNA can be inserted into new genomic loci. Due to the detrimental consequences of retrotransposition, plants like animals, have developed transcriptional and post-transcriptional silencing mechanisms. Recently several different genome-wide techniques have been developed to understand LTR retrotransposition in Arabidopsis and different plant species. Transposome, methylome, transcriptome, translatome and small RNA sequencing data have revealed how host silencing mechanisms can affect multiple steps of retrotransposition. These recent advances shed light on future mechanistic studies of retrotransposition as well as retrotransposon diversity.

Small RNAs guide histone methylation in Arabidopsis embryos.

Epigenetic reprogramming occurs during gametogenesis as well as during embryogenesis to reset the genome for early development. In flowering plants, many heterochromatic marks are maintained in sperm, but asymmetric DNA methylation is mostly lost. Asymmetric DNA methylation is dependent on small RNA but the re-establishment of silencing in embryo is not well understood. Here we demonstrate that small RNAs direct the histone H3 lysine 9 dimethylation during Arabidopsis thaliana embryonic development, together with asymmetric DNA methylation. This de novo silencing mechanism depends on the catalytic domain of SUVH9, a Su(Var)3-9 homolog thought to be catalytically inactive.

Circadian immunity, sunrise time and the seasonality of respiratory infections.

The innate and adaptive immune response are regulated by biological clocks, and circulating lymphocytes are lowest at sunrise. Accordingly, severity of disease in mouse models is highly dependent on the time of day of viral infection. Here, we explore whether circadian immunity contributes significantly to seasonality of respiratory viruses, including influenza and SARS-CoV-2. Susceptibility-Infection-Recovery-Susceptibility (SIRS) models of influenza and SIRS-derived models of COVID-19 suggest that local sunrise time is a better predictor of the basic reproductive number (R0) than climate, even when day length is taken into account. Moreover, these models predict a window of susceptibility when local sunrise time corresponds to the morning commute and contact rate is expected to be high. Counterfactual modeling suggests that retaining daylight savings time in the fall would reduce the length of this window, and substantially reduce seasonal waves of respiratory infections.

Male fertility in Arabidopsis requires active DNA demethylation of genes that control pollen tube function.

Active DNA demethylation is required for sexual reproduction in plants but the molecular determinants underlying this epigenetic control are not known. Here, we show in Arabidopsis thaliana that the DNA glycosylases DEMETER (DME) and REPRESSOR OF SILENCING 1 (ROS1) act semi-redundantly in the vegetative cell of pollen to demethylate DNA and ensure proper pollen tube progression. Moreover, we identify six pollen-specific genes with increased DNA methylation as well as reduced expression in dme and dme;ros1. We further show that for four of these genes, reinstalling their expression individually in mutant pollen is sufficient to improve male fertility. Our findings demonstrate an essential role of active DNA demethylation in regulating genes involved in pollen function.

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.

Loss of Small-RNA-Directed DNA Methylation in the Plant Cell Cycle Promotes Germline Reprogramming and Somaclonal Variation.

5-methyl cytosine is widespread in plant genomes in both CG and non-CG contexts. During replication, hemi-methylation on parental DNA strands guides symmetric CG methylation on nascent strands, but non-CG methylation requires modified histones and small RNA guides. Here, we used immortalized Arabidopsis cell suspensions to sort replicating nuclei and determine genome-wide cytosine methylation dynamics during the plant cell cycle. We find that symmetric mCG and mCHG are selectively retained in actively dividing cells in culture, whereas mCHH is depleted. mCG becomes transiently asymmetric during S phase but is rapidly restored in G2, whereas mCHG remains asymmetric throughout the cell cycle. Hundreds of loci gain ectopic CHG methylation, as well as 24-nt small interfering RNAs (siRNAs) and histone H3 lysine dimethylation (H3K9me2), without gaining CHH methylation. This suggests that spontaneous epialleles that arise in plant cell cultures are stably maintained by siRNA and H3K9me2 independent of the canonical RNA-directed DNA methylation (RdDM) pathway. In contrast, loci that fail to produce siRNA may be targeted for demethylation when the cell cycle arrests. Comparative analysis with methylomes of various tissues and cell types suggests that loss of small-RNA-directed non-CG methylation during DNA replication promotes germline reprogramming and epigenetic variation in plants propagated as clones.

Comparing DNA replication programs reveals large timing shifts at centromeres of endocycling cells in maize roots.

Plant cells undergo two types of cell cycles-the mitotic cycle in which DNA replication is coupled to mitosis, and the endocycle in which DNA replication occurs in the absence of cell division. To investigate DNA replication programs in these two types of cell cycles, we pulse labeled intact root tips of maize (Zea mays) with 5-ethynyl-2'-deoxyuridine (EdU) and used flow sorting of nuclei to examine DNA replication timing (RT) during the transition from a mitotic cycle to an endocycle. Comparison of the sequence-based RT profiles showed that most regions of the maize genome replicate at the same time during S phase in mitotic and endocycling cells, despite the need to replicate twice as much DNA in the endocycle and the fact that endocycling is typically associated with cell differentiation. However, regions collectively corresponding to 2% of the genome displayed significant changes in timing between the two types of cell cycles. The majority of these regions are small with a median size of 135 kb, shift to a later RT in the endocycle, and are enriched for genes expressed in the root tip. We found larger regions that shifted RT in centromeres of seven of the ten maize chromosomes. These regions covered the majority of the previously defined functional centromere, which ranged between 1 and 2 Mb in size in the reference genome. They replicate mainly during mid S phase in mitotic cells but primarily in late S phase of the endocycle. In contrast, the immediately adjacent pericentromere sequences are primarily late replicating in both cell cycles. Analysis of CENH3 enrichment levels in 8C vs 2C nuclei suggested that there is only a partial replacement of CENH3 nucleosomes after endocycle replication is complete. The shift to later replication of centromeres and possible reduction in CENH3 enrichment after endocycle replication is consistent with a hypothesis that centromeres are inactivated when their function is no longer needed.

Targeted reprogramming of H3K27me3 resets epigenetic memory in plant paternal chromatin.

Epigenetic marks are reprogrammed in the gametes to reset genomic potential in the next generation. In mammals, paternal chromatin is extensively reprogrammed through the global erasure of DNA methylation and the exchange of histones with protamines1,2. Precisely how the paternal epigenome is reprogrammed in flowering plants has remained unclear since DNA is not demethylated and histones are retained in sperm3,4. Here, we describe a multi-layered mechanism by which H3K27me3 is globally lost from histone-based sperm chromatin in Arabidopsis. This mechanism involves the silencing of H3K27me3 writers, activity of H3K27me3 erasers and deposition of a sperm-specific histone, H3.10 (ref. 5), which we show is immune to lysine 27 methylation. The loss of H3K27me3 facilitates the transcription of genes essential for spermatogenesis and pre-configures sperm with a chromatin state that forecasts gene expression in the next generation. Thus, plants have evolved a specific mechanism to simultaneously differentiate male gametes and reprogram the paternal epigenome.