Identification of differentially methylated regions during vernalization revealed a role for RNA methyltransferases in bolting

Role of methylation in vernalization and photoperiod pathway: a potential flowering regulator?

Recognized as a pivotal developmental transition, flowering marks the continuation of a plant's life cycle. Vernalization and photoperiod are two major flowering pathways orchestrating numerous florigenic signals. Methylation, including histone, DNA and RNA methylation, is one of the recent foci in plant development. Considerable studies reveal that methylation seems to show an increasing potential regulatory role in plant flowering via altering relevant gene expression without altering the genetic basis. However, little has been reviewed about whether and how methylation acts on vernalization- and photoperiod-induced flowering before and after FLOWERING LOCUS C (FLC) reactivation, what role RNA methylation plays in vernalization- and photoperiod-induced flowering, how methylation participates simultaneously in both vernalization- and photoperiod-induced flowering, the heritability of methylation memory under the vernalization/photoperiod pathway, and whether and how methylation replaces vernalization/photoinduction to regulate flowering. Our review provides insight about the crosstalk among the genetic control of the flowering gene network, methylation (methyltransferases/demethylases) and external signals (cold, light, sRNA and phytohormones) in vernalization and photoperiod pathways. The existing evidence that RNA methylation may play a potential regulatory role in vernalization- and photoperiod-induced flowering has been gathered and represented for the first time. This review speculates about and discusses the possibility of substituting methylation for vernalization and photoinduction to promote flowering. Current evidence is utilized to discuss the possibility of future methylation reagents becoming flowering regulators at the molecular level.

RNA methyltransferases in plants: Breakthroughs in function and evolution

Each day it is becoming increasingly difficult not to notice the completely new, fast growing, extremely intricate and challenging world of epitranscriptomics as the understanding of RNA methylation is expanding at a hasty rate. Writers (methyltransferases), erasers (demethylases) and readers (RNA-binding proteins) are responsible for adding, removing and recognising methyl groups on RNA, respectively. Several methyltransferases identified in plants are now being investigated and recent studies have shown a connection between RNA-methyltransferases (RNA-MTases) and stress and development processes. However, compared to their animal and bacteria counterparts, the understanding of RNA methyltransferases is still incipient, particularly those located in organelles. Comparative and systematic analyses allowed the tracing of the evolution of these enzymes suggesting the existence of several methyltransferases yet to be characterised. This review outlines the functions of plant nuclear and organellar RNA-MTases in plant development and stress responses and the comparative and evolutionary discoveries made on RNA-MTases across kingdoms.

Genome-wide identification of the KNOTTED HOMEOBOX gene family and their involvement in stalk development in flowering Chinese cabbage

Gibberellin and cytokinin synergistically regulate the stalk development in flowering Chinese cabbage. KNOX proteins were reported to function as important regulators of the shoot apex to promote meristem activity by synchronously inducing CTK and suppressing GA biosynthesis, while their regulatory mechanism in the bolting and flowering is unknown. In this study, 9 BcKNOX genes were identified and mapped unevenly on 6 out of 10 flowering Chinese cabbage chromosomes. The BcKNOXs were divided into three subfamilies on the basis of sequences and gene structure. The proteins contain four conserved domains except for BcKNATM. Three BcKNOX TFs (BcKNOX1, BcKNOX3, and BcKNOX5) displayed high transcription levels on tested tissues at various stages. The major part of BcKNOX genes showed preferential expression patterns in response to low-temperature, zeatin (ZT), and GA₃ treatment, indicating that they were involved in bud differentiation and bolting. BcKNOX1 and BcKNOX5 showed high correlation level with gibberellins synthetase, and CTK metabolic genes. BcKONX1 also showed high correlation coefficients within BcRGA1 and BcRGL1 which are negative regulators of GA signaling. In addition, BcKNOX1 interacted with BcRGA1 and BcRGL1, as confirmed by yeast two-hybrid (Y2H) and biomolecular fluorescence complementation assay (BiFC). This analysis has provided useful foundation for the future functional roles' analysis of flowering Chinese cabbage KNOX genes.

Walnut N-Acetylserotonin Methyltransferase Gene Family Genome-Wide Identification and Diverse Functions Characterization During Flower Bud Development

Melatonin widely mediates multiple developmental dynamics in plants as a vital growth stimulator, stress protector, and developmental regulator. N-acetylserotonin methyltransferase (ASMT) is the key enzyme that catalyzes the final step of melatonin biosynthesis in plants and plays an essential role in the plant melatonin regulatory network. Studies of ASMT have contributed to understanding the mechanism of melatonin biosynthesis in plants. However, AMST gene is currently uncharacterized in most plants. In this study, we characterized the JrASMT gene family using bioinformatics in a melatonin-rich plant, walnut. Phylogenetic, gene structure, conserved motifs, promoter elements, interacting proteins and miRNA analyses were also performed. The expansion and differentiation of the ASMT family occurred before the onset of the plant terrestrialization. ASMT genes were more differentiated in dicotyledonous plants. Forty-six ASMT genes were distributed in clusters on 10 chromosomes of walnut. Four JrASMT genes had homologous relationships both within walnut and between species. Cis-regulatory elements showed that JrASMT was mainly induced by light and hormones, and targeted cleavage of miRNA172 and miR399 may be an important pathway to suppress JrASMT expression. Transcriptome data showed that 13 JrASMT were differentially expressed at different periods of walnut bud development. WGCNA showed that JrASMT1/10/13/23 were coexpressed with genes regulating cell fate and epigenetic modifications during early physiological differentiation of walnut female flower buds. JrASMT12/28/37/40 were highly expressed during morphological differentiation of flower buds, associated with altered stress capacity of walnut flower buds, and predicted to be involved in the regulatory network of abscisic acid, salicylic acid, and cytokinin in walnut. The qRT-PCR validated the results of differential expression analysis and further provided three JrASMT genes with different expression profiles in walnut flower bud development. Our study explored the evolutionary relationships of the plant ASMT gene family and the functional characteristics of walnut JrASMT. It provides a valuable perspective for further understanding the complex melatonin mechanisms in plant developmental regulation.

Molecular markers based on sequence variation in BoFLC1.C9 for characterizing early- and late-flowering cabbage genotypes

BACKGROUND

Cabbage (Brassica oleracea var. capitata) is popular worldwide for consumption as a leafy vegetable. Premature flowering is triggered by low temperature, and deteriorates quality of cabbage as vegetable. In general, growers prefer late-flowering varieties to assure good quality compact head. Here, we report BoFLC1.C9 as a gene with clear sequence variation between cabbage lines with different flowering times, and proposed as molecular marker to characterize early- and late-flowering cabbage lines.

RESULTS

We identified sequence variation of 67 bp insertions in intron 2, which were contributed in flowering time variation between two inbred lines through rapid down-regulation of the BoFLC1.C9 gene in early-flowering line compared to late-flowering one upon vernalization. One set of primer 'F7R7' proposed as marker, of which was explained with 83 and 80% of flowering time variation in 141 F₂ individuals and 20 commercial lines, respectively.

CONCLUSIONS

This F7R7 marker could be used as genetic tools to characterize flowering time variation and to select as well to develop early- and late-flowering cabbage cultivars.

Haplotype Variation of Flowering Time Genes of Sugar Beet and Its Wild Relatives and the Impact on Life Cycle Regimes

The species Beta vulgaris encompasses wild and cultivated members with a broad range of phenological development. The annual life cycle is commonly found in sea beets (ssp. maritima) from Mediterranean environments which germinate, bolt, and flower within one season under long day conditions. Biennials such as the cultivated sugar beet (B. vulgaris ssp. vulgaris) as well as sea beets from northern latitudes require prolonged exposure to cold temperature over winter to acquire floral competence. Sugar beet is mainly cultivated for sugar production in Europe and is likely to have originated from sea beet. Flowering time strongly affects seed yield and yield potential and is thus a trait of high agronomic relevance. Besides environmental cues, there are complex genetic networks known to impact life cycle switch in flowering plants. In sugar beet, BTC1, BvBBX19, BvFT1, and BvFT2 are major flowering time regulators. In this study, we phenotyped plants from a diversity Beta panel encompassing cultivated and wild species from different geographical origin. Plants were grown under different day length regimes with and without vernalization. Haplotype analysis of BTC1, BvBBX19, BvFT1, and BvFT2 was performed to identify natural diversity of these genes and their impact on flowering. We found that accessions from northern latitudes flowered significantly later than those from southern latitudes. Some plants did not flower at all, indicating a strong impact of latitude of origin on life cycle. Haplotype analysis revealed a high conservation of the CCT-, REC-, BBX-, and PEBP-domains with regard to SNP occurrence. We identified sequence variation which may impact life cycle adaptation in beet. Our data endorse the importance of BTC1 in the domestication process of cultivated beets and contribute to the understanding of distribution and adaption of Beta species to different life cycle regimes in response to different environments. Moreover, our data provide a resource for haplotypes identified for the major floral regulators in beet.

Subcellular localization of Arabidopsis arogenate dehydratases suggests novel and non-enzymatic roles

Arogenate dehydratases (ADTs) catalyze the final step in phenylalanine biosynthesis in plants. The Arabidopsis thaliana genome encodes a family of six ADTs capable of decarboxylating/dehydrating arogenate into phenylalanine. Using cyan fluorescent protein (CFP)-tagged proteins, the subcellular localization patterns of all six A. thaliana ADTs were investigated in intact Nicotiana benthamiana and A. thaliana leaf cells. We show that A. thaliana ADTs localize to stroma and stromules (stroma-filled tubules) of chloroplasts. This localization pattern is consistent with the enzymatic function of ADTs as many enzymes required for amino acid biosynthesis are primarily localized to chloroplasts, and stromules are thought to increase metabolite transport from chloroplasts to other cellular compartments. Furthermore, we provide evidence that ADTs have additional, non-enzymatic roles. ADT2 localizes in a ring around the equatorial plane of chloroplasts or to a chloroplast pole, which suggests that ADT2 is a component of the chloroplast division machinery. In addition to chloroplasts, ADT5 was also found in nuclei, again suggesting a non-enzymatic role for ADT5. We also show evidence that ADT5 is transported to the nucleus via stromules. We propose that ADT2 and ADT5 are moonlighting proteins that play an enzymatic role in phenylalanine biosynthesis and a second role in chloroplast division or transcriptional regulation, respectively.

A Detailed Analysis of the BR1 Locus Suggests a New Mechanism for Bolting after Winter in Sugar Beet (Beta vulgaris L.)

Sugar beet (Beta vulgaris ssp. vulgaris) is a biennial, sucrose-storing plant, which is mainly cultivated as a spring crop and harvested in the vegetative stage before winter. For increasing beet yield, over-winter cultivation would be advantageous. However, bolting is induced after winter and drastically reduces yield. Thus, post-winter bolting control is essential for winter beet cultivation. To identify genetic factors controlling bolting after winter, a F₂ population was previously developed by crossing the sugar beet accessions BETA 1773 with reduced bolting tendency and 93161P with complete bolting after winter. For a mapping-by-sequencing analysis, pools of 26 bolting-resistant and 297 bolting F₂ plants were used. Thereby, a single continuous homozygous region of 103 kb was co-localized to the previously published BR₁ QTL for post-winter bolting resistance (Pfeiffer et al., 2014). The BR₁ locus was narrowed down to 11 candidate genes from which a homolog of the Arabidopsis CLEAVAGE AND POLYADENYLATION SPECIFICITY FACTOR 73-I (CPSF73-I) was identified as the most promising candidate. A 2 bp deletion within the BETA 1773 allele of BvCPSF73-Ia results in a truncated protein. However, the null allele of BvCPSF73-Ia might partially be compensated by a second BvCPSF73-Ib gene. This gene is located 954 bp upstream of BvCPSF73-Ia and could be responsible for the incomplete penetrance of the post-winter bolting resistance allele of BETA 1773. This result is an important milestone for breeding winter beets with complete bolting resistance after winter.

Deciphering the epitranscriptome: A green perspective

The advent of high-throughput sequencing technologies coupled with new detection methods of RNA modifications has enabled investigation of a new layer of gene regulation - the epitranscriptome. With over 100 known RNA modifications, understanding the repertoire of RNA modifications is a huge undertaking. This review summarizes what is known about RNA modifications with an emphasis on discoveries in plants. RNA ribose modifications, base methylations and pseudouridylation are required for normal development in Arabidopsis, as mutations in the enzymes modifying them have diverse effects on plant development and stress responses. These modifications can regulate RNA structure, turnover and translation. Transfer RNA and ribosomal RNA modifications have been mapped extensively and their functions investigated in many organisms, including plants. Recent work exploring the locations, functions and targeting of N⁶ -methyladenosine (m⁶ A), 5-methylcytosine (m⁵ C), pseudouridine (Ψ), and additional modifications in mRNAs and ncRNAs are highlighted, as well as those previously known on tRNAs and rRNAs. Many questions remain as to the exact mechanisms of targeting and functions of specific modified sites and whether these modifications have distinct functions in the different classes of RNAs.

Epigenomics and bolting tolerance in sugar beet genotypes

In sugar beet (Beta vulgaris altissima), bolting tolerance is an essential agronomic trait reflecting the bolting response of genotypes after vernalization. Genes involved in induction of sugar beet bolting have now been identified, and evidence suggests that epigenetic factors are involved in their control. Indeed, the time course and amplitude of DNA methylation variations in the shoot apical meristem have been shown to be critical in inducing sugar beet bolting, and a few functional targets of DNA methylation during vernalization have been identified. However, molecular mechanisms controlling bolting tolerance levels among genotypes are still poorly understood. Here, gene expression and DNA methylation profiles were compared in shoot apical meristems of three bolting-resistant and three bolting-sensitive genotypes after vernalization. Using Cot fractionation followed by 454 sequencing of the isolated low-copy DNA, 6231 contigs were obtained that were used along with public sugar beet DNA sequences to design custom Agilent microarrays for expression (56k) and methylation (244k) analyses. A total of 169 differentially expressed genes and 111 differentially methylated regions were identified between resistant and sensitive vernalized genotypes. Fourteen sequences were both differentially expressed and differentially methylated, with a negative correlation between their methylation and expression levels. Genes involved in cold perception, phytohormone signalling, and flowering induction were over-represented and collectively represent an integrative gene network from environmental perception to bolting induction. Altogether, the data suggest that the genotype-dependent control of DNA methylation and expression of an integrative gene network participate in bolting tolerance in sugar beet, opening up perspectives for crop improvement.

How to Isolate a Plant's Hypomethylome in One Shot

Genome assembly remains a challenge for large and/or complex plant genomes due to their abundant repetitive regions resulting in studies focusing on gene space instead of the whole genome. Thus, DNA enrichment strategies facilitate the assembly by increasing the coverage and simultaneously reducing the complexity of the whole genome. In this paper we provide an easy, fast, and cost-effective variant of MRE-seq to obtain a plant's hypomethylome by an optimized methyl filtration protocol followed by next generation sequencing. The method is demonstrated on three plant species with knowingly large and/or complex (polyploid) genomes: Oryza sativa, Picea abies, and Crocus sativus. The identified hypomethylomes show clear enrichment for genes and their flanking regions and clear reduction of transposable elements. Additionally, genomic sequences around genes are captured including regulatory elements in introns and up- and downstream flanks. High similarity of the results obtained by a de novo assembly approach with a reference based mapping in rice supports the applicability for studying and understanding the genomes of nonmodel organisms. Hence we show the high potential of MRE-seq in a wide range of scenarios for the direct analysis of methylation differences, for example, between ecotypes, individuals, within or across species harbouring large, and complex genomes.

Conservation of tRNA and rRNA 5-methylcytosine in the kingdom Plantae

BACKGROUND

Post-transcriptional methylation of RNA cytosine residues to 5-methylcytosine (m(5)C) is an important modification that regulates RNA metabolism and occurs in both eukaryotes and prokaryotes. Yet, to date, no transcriptome-wide identification of m(5)C sites has been undertaken in plants. Plants provide a unique comparative system for investigating the origin and evolution of m(5)C as they contain three different genomes, the nucleus, mitochondria and chloroplast. Here we use bisulfite conversion of RNA combined with high-throughput IIlumina sequencing (RBS-seq) to identify single-nucleotide resolution of m(5)C sites in non-coding ribosomal RNAs and transfer RNAs of all three sub-cellular transcriptomes across six diverse species that included, the single-celled algae Nannochloropsis oculata, the macro algae Caulerpa taxifolia and multi-cellular higher plants Arabidopsis thaliana, Brassica rapa, Triticum durum and Ginkgo biloba.

RESULTS

Using the plant model Arabidopsis thaliana, we identified a total of 39 highly methylated m(5)C sites in predicted structural positions of nuclear tRNAs and 7 m(5)C sites in rRNAs from nuclear, chloroplast and mitochondrial transcriptomes. Both the nucleotide position and percent methylation of tRNAs and rRNAs m(5)C sites were conserved across all species analysed, from single celled algae N. oculata to multicellular plants. Interestingly the mitochondrial and chloroplast encoded tRNAs were devoid of m(5)C in A. thaliana and this is generally conserved across Plantae. This suggests independent evolution of organelle methylation in animals and plants, as animal mitochondrial tRNAs have m(5)C sites. Here we characterize 5 members of the RNA 5-methylcytosine family in Arabidopsis and extend the functional characterization of TRDMT1 and NOP2A/OLI2. We demonstrate that nuclear tRNA methylation requires two evolutionarily conserved methyltransferases, TRDMT1 and TRM4B. trdmt1 trm4b double mutants are hypersensitive to the antibiotic hygromycin B, demonstrating the function of tRNA methylation in regulating translation. Additionally we demonstrate that nuclear large subunit 25S rRNA methylation requires the conserved RNA methyltransferase NSUN5. Our results also suggest functional redundancy of at least two of the NOP2 paralogs in Arabidopsis.

CONCLUSIONS

Our data demonstrates widespread occurrence and conservation of non-coding RNA methylation in the kingdom Plantae, suggesting important and highly conserved roles of this post-transcriptional modification.

The FLC-like gene BvFL1 is not a major regulator of vernalization response in biennial beets

Many plant species in temperate climate regions require vernalization over winter to initiate flowering. Flowering Locus C (FLC) and FLC-like genes are key regulators of vernalization requirement and growth habit in winter-annual and perennial Brassicaceae. In the biennial crop species Beta vulgaris ssp. vulgaris in the evolutionarily distant Caryophyllales clade of core eudicots growth habit and bolting time are controlled by the vernalization and photoperiod response gene BTC1 and the downstream BvFT1-BvFT2 module. B. vulgaris also contains a vernalization-responsive FLC homolog (BvFL1). Here, to further elucidate the regulation of vernalization response and growth habit in beet, we functionally characterized BvFL1 by RNAi and over-expression in transgenic plants. BvFL1 RNAi neither eliminated the requirement for vernalization of biennial beets nor had a major effect on bolting time after vernalization. Over-expression of BvFL1 resulted in a moderate late-bolting phenotype, with bolting after vernalization being delayed by approximately 1 week. By contrast, RNAi-induced down-regulation of the BvFT1-BvFT2 module led to a strong delay in bolting after vernalization by several weeks. The data demonstrate for the first time that an FLC homolog does not play a major role in the control of vernalization response in a dicot species outside the Brassicaceae.