Genome-wide identification of endogenous RNA-directed DNA methylation loci associated with abundant 21-nucleotide siRNAs in Arabidopsis

Endogenous viral elements are targeted by RNA silencing pathways in banana

Endogenous banana streak virus (eBSV) integrants derived from three distinct species, present in Musa balbisiana (B) but not Musa acuminata (A) banana genomes are able to reconstitute functional episomal viruses causing banana streak disease in interspecific triploid AAB banana hybrids but not in the diploid (BB) parent line, which harbours identical eBSV loci. Here, we investigated the regulation of these eBSV. In-depth characterization of siRNAs, transcripts and methylation derived from eBSV using Illumina and bisulfite sequencing were carried out on eBSV-free Musa acuminata AAA plants and BB or AAB banana plants with eBSV. eBSV loci produce low-abundance transcripts covering most of the viral sequence and generate predominantly 24-nt siRNAs. siRNA accumulation is restricted to duplicated and inverted viral sequences present in eBSV. Both siRNA-accumulating and nonaccumulating sequences of eBSV in BB plants are heavily methylated in all three CG, CHG and CHH contexts. Our data suggest that eBSVs are controlled at the epigenetic level in BB diploids. This regulation not only prevents their awakening and systemic infection of the plant but is also probably involved in the inherent resistance of the BB plants to mealybug-transmitted viral infection. These findings are thus of relevance to other plant resources hosting integrated viruses.

DDM1-Mediated TE Silencing in Plants

Epigenetic modifications are indispensable for regulating gene bodies and TE silencing. DECREASE IN DNA METHYLATION 1 (DDM1) is a chromatin remodeller involved in histone modifications and DNA methylation. Apart from maintaining the epigenome, DDM1 also maintains key plant traits such as flowering time and heterosis. The role of DDM1 in epigenetic regulation is best characterised in plants, especially arabidopsis, rice, maize and tomato. The epigenetic changes induced by DDM1 establish the stable inheritance of many plant traits for at least eight generations, yet DDM1 does not methylate protein-coding genes. The DDM1 TE silencing mechanism is distinct and has evolved independently of other silencing pathways. Unlike the RNA-directed DNA Methylation (RdDM) pathway, DDM1 does not depend on siRNAs to enforce the heterochromatic state of TEs. Here, we review DDM1 TE silencing activity in the RdDM and non-RdDM contexts. The DDM1 TE silencing machinery is strongly associated with the histone linker H1 and histone H2A.W. While the linker histone H1 excludes the RdDM factors from methylating the heterochromatin, the histone H2A.W variant prevents TE mobility. The DDM1-H2A.W strategy alone silences nearly all the mobile TEs in the arabidopsis genome. Thus, the DDM1-directed TE silencing essentially preserves heterochromatic features and abolishes mobile threats to genome stability.

Homeobox 27, a Homeodomain Transcription Factor, Confers Tolerances to CMV by Associating with Cucumber Mosaic Virus 2b Protein

Transcription factors (TFs) play an important role in plant development; however, their role during viral infection largely remains unknown. The present study was designed to uncover the role transcription factors play in Cucumber mosaic virus (CMV) infection. During the screening of an Arabidopsis thaliana (Col-0) transcription factor library, using the CMV 2b protein as bait in the yeast two-hybrid system, the 2b protein interacted with Homeobox protein 27 (HB27). HB27 belongs to the zinc finger homeodomain family and is known to have a regulatory role in flower development, and responses to biotic and abiotic stress. The interaction between CMV 2b and HB27 proteins was further validated using in planta (bimolecular fluorescence complementation assay) and in vitro far-Western blotting (FWB) methods. In the bimolecular fluorescence complementation assay, these proteins reconstituted YFP fluorescence in the nucleus and the cytoplasmic region as small fluorescent dots. In FWB, positive interaction was detected using bait anti-MYC antibody on the target HB27-HA protein. During CMV infection, upregulation (~3-fold) of the HB27 transcript was observed at 14 days post-infection (dpi) in A. thaliana plants, and expression declined to the same as healthy plants at 21 dpi. To understand the role of the HB27 protein during CMV infection, virus accumulation was determined in HB27-overexpressing (HB27 OE) and knockout mutants. In HB27-overexpressing lines, infected plants developed mild symptoms, accumulating a lower virus titer at 21 dpi compared to wild-type plants. Additionally, knockout HB27 mutants had more severe symptoms and a higher viral accumulation than wild-type plants. These results indicate that HB27 plays an important role in the regulation of plant defense against plant virus infection.

Cucumber RDR1s and cucumber mosaic virus suppressor protein 2b association directs host defence in cucumber plants

RNA-dependent RNA polymerases (RDRs) regulate important aspects of plant development and resistance to pathogens. The role of RDRs in virus resistance has been demonstrated using siRNA signal amplification and through the methylation of viral genomes. Cucumber (Cucumis sativus) has four RDR1 genes that are differentially induced during virus infection: CsRDR1a, CsRDR1b, and duplicated CsRDR1c1/c2. The mode of action of CsRDR1s during viral infection is unknown. Transient expression of the cucumber mosaic virus (CMV)-2b protein (the viral suppressor of RNA silencing) in cucumber protoplasts induced the expression of CsRDR1c, but not of CsRDR1a/1b. Results from the yeast two-hybrid system showed that CsRDR1 proteins interacted with CMV-2b and this was confirmed by bimolecular fluorescence complementation assays. In protoplasts, CsRDR1s localized in the cytoplasm as punctate spots. Colocalization experiments revealed that CsRDR1s and CMV-2b were uniformly dispersed throughout the cytoplasm, suggesting that CsRDR1s are redistributed as a result of interactions. Transient overexpression of individual CsRDR1a/1b genes in protoplasts reduced CMV accumulation, indicating their antiviral role. However, overexpression of CsRDR1c in protoplasts resulted in relatively higher accumulation of CMV and CMVΔ2b. In single cells, CsRDR1c enhances viral replication, leading to CMV accumulation and blocking secondary siRNA amplification of CsRDR1c by CMV-2b protein. This suggests that CMV-2b acts as both a transcription factor that induces CsRDR1c (controlling virus accumulation) and a suppressor of CsRDR1c activity.

Environmental RNAi pathways in the two-spotted spider mite

Background

RNA interference (RNAi) regulates gene expression in most multicellular organisms through binding of small RNA effectors to target transcripts. Exploiting this process is a popular strategy for genetic manipulation and has applications that includes arthropod pest control. RNAi technologies are dependent on delivery method with the most convenient likely being feeding, which is effective in some animals while others are insensitive. The two-spotted spider mite, Tetranychus urticae, is prime candidate for developing RNAi approaches due to frequent occurrence of conventional pesticide resistance. Using a sequencing-based approach, the fate of ingested RNAs was explored to identify features and conditions that affect small RNA biogenesis from external sources to better inform RNAi design.

Results

Biochemical and sequencing approaches in conjunction with extensive computational assessment were used to evaluate metabolism of ingested RNAs in T. urticae. This chelicerae arthropod shows only modest response to oral RNAi and has biogenesis pathways distinct from model organisms. Processing of synthetic and plant host RNAs ingested during feeding were evaluated to identify active substrates for spider mite RNAi pathways. Through cataloging characteristics of biochemically purified RNA from these sources, trans-acting small RNAs could be distinguished from degradation fragments and their origins documented.

Conclusions

Using a strategy that delineates small RNA processing, we found many transcripts have the potential to enter spider mite RNAi pathways, however, trans-acting RNAs appear very unstable and rare. This suggests potential RNAi pathway substrates from ingested materials are mostly degraded and infrequently converted into regulators of gene expression. Spider mites infest a variety of plants, and it would be maladaptive to generate diverse gene regulators from dietary RNAs. This study provides a framework for assessing RNAi technology in organisms where genetic and biochemical tools are absent and benefit rationale design of RNAi triggers for T.urticae.

Secondary siRNAs in Plants: Biosynthesis, Various Functions, and Applications in Virology

The major components of RNA silencing include both transitive and systemic small RNAs, which are technically called secondary sRNAs. Double-stranded RNAs trigger systemic silencing pathways to negatively regulate gene expression. The secondary siRNAs generated as a result of transitive silencing also play a substantial role in gene silencing especially in antiviral defense. In this review, we first describe the discovery and pathways of transitivity with emphasis on RNA-dependent RNA polymerases followed by description on the short range and systemic spread of silencing. We also provide an in-depth view on the various size classes of secondary siRNAs and their different roles in RNA silencing including their categorization based on their biogenesis. The other regulatory roles of secondary siRNAs in transgene silencing, virus-induced gene silencing, transitivity, and trans-species transfer have also been detailed. The possible implications and applications of systemic silencing and the different gene silencing tools developed are also described. The details on mobility and roles of secondary siRNAs derived from viral genome in plant defense against the respective viruses are presented. This entails the description of other compatible plant-virus interactions and the corresponding small RNAs that determine recovery from disease symptoms, exclusion of viruses from shoot meristems, and natural resistance. The last section presents an overview on the usefulness of RNA silencing for management of viral infections in crop plants.

Recent advances in understanding plant antiviral RNAi and viral suppressors of RNAi

Molecular plant-virus interactions provide an excellent model to understanding host antiviral immunity and viral counter-defense mechanisms. The primary antiviral defense is triggered inside the infected plant cell by virus-derived small-interfering RNAs, which guide homology-dependent RNA interference (RNAi) and/or RNA-directed DNA methylation (RdDM) to target RNA and DNA viruses. In counter-defense, plant viruses have independently evolved viral suppressors of RNAi (VSRs) to specifically antagonize antiviral RNAi. Recent studies have shown that plant antiviral responses are regulated by endogenous small silencing RNAs, RNA decay and autophagy and that some known VSRs of plant RNA and DNA viruses also target these newly recognized defense responses to promote infection. This review focuses on these recent advances that have revealed multilayered regulation of plant-virus interactions.

Epigenetic Variation at a Genomic Locus Affecting Biomass Accumulation under Low Nitrogen in Arabidopsis thaliana

Nitrogen (N) is a macronutrient determining crop yield. The application of N fertilisers can substantially increase the yield, but excess use also causes the nitrate pollution of water resources and increases production costs. Increasing N use efficiency (NUE) in crop plants is an important step to implement low-input agricultural systems. We used Arabidopsis thaliana as model system to investigate the natural genetic diversity in traits related to NUE. Natural variation was used to study adaptive growth patterns and changes in gene expression associated with limited nitrate availability. A genome-wide association study revealed an association of eight SNP markers on Chromosome 1 with shoot growth under limited N. The identified linkage disequilibrium (LD) interval includes the DNA sequences of three cysteine/histidine-rich C1 domain proteins in tandem orientation. These genes differ in promoter structure, methylation pattern and expression level among accessions, correlating with growth performance under N deficiency. Our results suggest the involvement of epigenetic regulation in the expression of NUE-related traits.

The salt-induced transcription factor GmMYB84 confers salinity tolerance in soybean

Transcription factor activation and DNA methylation are important plant responses to abiotic stress. Here, we established that the salinity stress-induced expression of the soybean (Glycine max) transcription factor-encoding gene GmMYB84 relies on DNA methylation. The level of DNA methylation at sequences 690 nt to 950 nt upstream of the GmMYB84 transcription initiation codon was markedly reduced in plants exposed to salinity stress, resulting in a higher abundance of transcripts. When challenged with salinity stress, plants constitutively expressing GmMYB84 outperformed untransformed plants with respect to their germination rate, primary root elongation, proline accumulation, antioxidant enzyme activity, membrane integrity, and K+ levels. Arabidopsis thaliana plants heterologously expressing GmMYB84 were more tolerant to salt stress and exhibited higher germination rates than the wild type. Electrophoretic mobility shift assays revealed that GmMYB84 binds to the cis-regulatory sequences of GmAKT1, the homolog of ARABIDOPSIS K+ TRANSPORTER 1 (AKT1). Thus, DNA methylation modulates the salinity stress-induced expression of the soybean transcription factor-encoding gene GmMYB84 and thereby confers salinity stress tolerance.

Crosstalk between epigenetic silencing and infection by tobacco rattle virus in Arabidopsis

DNA methylation is an important epigenetic mechanism for controlling innate immunity against microbial pathogens in plants. Little is known, however, about the manner in which viral infections interact with DNA methylation pathways. Here we investigate the crosstalk between epigenetic silencing and viral infections in Arabidopsis inflorescences. We found that tobacco rattle virus (TRV) causes changes in the expression of key transcriptional gene silencing factors with RNA-directed DNA methylation activities that coincide with changes in methylation at the whole genome level. Viral susceptibility/resistance was altered in DNA (de)methylation-deficient mutants, suggesting that DNA methylation is an important regulatory system controlling TRV proliferation. We further show that several transposable elements (TEs) underwent transcriptional activation during TRV infection, and that TE regulation likely involved both DNA methylation-dependent and -independent mechanisms. We identified a cluster of disease resistance genes regulated by DNA methylation in infected plants that were enriched for TEs in their promoters. Interestingly, TEs and nearby resistance genes were co-regulated in TRV-infected DNA (de)methylation mutants. Our study shows that DNA methylation contributes to modulate the outcome of viral infections in Arabidopsis, and opens up new possibilities for exploring the role of TE regulation in antiviral defence.

Roles of Small RNAs in Virus-Plant Interactions

Small RNAs (sRNAs), including microRNAs (miRNAs) and short interfering RNAs (siRNAs), are non-coding but powerful RNA molecules of 20–30 nucleotides in length. sRNAs play crucial regulatory roles in diverse plant biological processes. Recently, many studies on sRNAs have been reported. We summarize new findings of sRNAs in virus-plant interactions to accelerate the function analysis of sRNAs. The main content of this review article includes three parts: virus-responsive sRNAs, function analysis of sRNAs in virus pathogenicity or host resistance, and some sRNAs-mediated underlying mechanisms in virus-plant interactions. New findings of sRNAs deepen our understanding about sRNAs’ roles, which might contribute to the design of novel control measures against plant viruses.

CMV2b-Dependent Regulation of Host Defense Pathways in the Context of Viral Infection

RNA silencing (or RNA interference, RNAi) plays direct roles in plant host defenses against viruses. Viruses encode suppressors of RNAi (VSRs) to counteract host antiviral defenses. The generation of transgenic plants expressing VSRs facilitates the understanding of the mechanisms of VSR-mediated interference with the endogenous silencing pathway. However, studying VSRs independent of other viral components simplifies the complex roles of VSRs during natural viral infection. While suppression of transgene silencing by the VSR 2b protein encoded by cucumber mosaic virus (CMV) requires 2b-small RNA (sRNA) binding activity, suppression of host antiviral defenses requires the binding activity of both sRNAs and AGOs proteins. This study, aimed to understand the functions of 2b in the context of CMV infection; thus, we performed genome-wide analyses of differential DNA methylation regions among wild-type CMV-infected, CMVΔ2b-infected, and 2b-transgenic Arabidopsis plants. These analyses, together with transcriptome sequencing and RT-qPCR analyses, show that while the majority of induced genes in 2b-transgenic plants were involved in extensive metabolic pathways, CMV-infection 2b-dependent induced genes were enriched in plant immunity pathways, including salicylic acid (SA) signaling. Together with infection with CMV mutants that expressed the 2b functional domains of sRNA or AGO binding, our data demonstrate that CMV-accelerated SA signaling depends on 2b-sRNA binding activity which is also responsible for virulence.

Reproductive phasiRNAs in grasses are compositionally distinct from other classes of small RNAs

Little is known about the characteristics and function of reproductive phased, secondary, small interfering RNAs (phasiRNAs) in the Poaceae, despite the availability of significant genomic resources, experimental data, and a growing number of computational tools. We utilized machine-learning methods to identify sequence-based and positional features that distinguish phasiRNAs in rice and maize from other small RNAs (sRNAs). We developed Random Forest classifiers that can distinguish reproductive phasiRNAs from other sRNAs in complex sets of sequencing data, utilizing sequence-based (k-mers) and features describing position-specific sequence biases. The classification performance attained is > 80% in accuracy, sensitivity, specificity, and positive predicted value. Feature selection identified important features in both ends of phasiRNAs. We demonstrated that phasiRNAs have strand specificity and position-specific nucleotide biases potentially influencing AGO sorting; we also predicted targets to infer functions of phasiRNAs, and computationally assessed their sequence characteristics relative to other sRNAs. Our results demonstrate that machine-learning methods effectively identify phasiRNAs despite the lack of characteristic features typically present in precursor loci of other small RNAs, such as sequence conservation or structural motifs. The 5'-end features we identified provide insights into AGO-phasiRNA interactions. We describe a hypothetical model of competition for AGO loading between phasiRNAs of different nucleotide compositions.

Transposon-associated epigenetic silencing during Pleurotus ostreatus life cycle

Transposable elements constitute an important fraction of eukaryotic genomes. Given their mutagenic potential, host-genomes have evolved epigenetic defense mechanisms to limit their expansion. In fungi, epigenetic modifications have been widely studied in ascomycetes, although we lack a global picture of the epigenetic landscape in basidiomycetes. In this study, we analysed the genome-wide epigenetic and transcriptional patterns of the white-rot basidiomycete Pleurotus ostreatus throughout its life cycle. Our results performed by using high-throughput sequencing analyses revealed that strain-specific DNA methylation profiles are primarily involved in the repression of transposon activity and suggest that 21 nt small RNAs play a key role in transposon silencing. Furthermore, we provide evidence that transposon-associated DNA methylation, but not sRNA production, is directly involved in the silencing of genes surrounded by transposons. Remarkably, we found that nucleus-specific methylation levels varied in dikaryotic strains sharing identical genetic complement but different subculture conditions. Finally, we identified key genes activated in the fruiting process through the comparative analysis of transcriptomes. This study provides an integrated picture of epigenetic defense mechanisms leading to the transcriptional silencing of transposons and surrounding genes in basidiomycetes. Moreover, our findings suggest that transcriptional but not methylation reprogramming triggers fruitbody development in P. ostreatus.

Epigenetic Changes in the Regulation of Nicotiana tabacum Response to Cucumber Mosaic Virus Infection and Symptom Recovery through Single-Base Resolution Methylomes

Plants have evolved multiple mechanisms to respond to viral infection. These responses have been studied in detail at the level of host immune response and antiviral RNA silencing (RNAi). However, the possibility of epigenetic reprogramming has not been thoroughly investigated. Here, we identified the role of DNA methylation during viral infection and performed reduced representation bisulfite sequencing (RRBS) on tissues of Cucumber mosaic virus (CMV)-infected Nicotiana tabacum at various developmental stages. Differential methylated regions are enriched with CHH sequence contexts, 80% of which are located on the gene body to regulate gene expression in a temporal style. The methylated genes depressed by methyltransferase inhibition largely overlapped with methylated genes in response to viral invasion. Activation in the argonaute protein and depression in methyl donor synthase revealed the important role of dynamic methylation changes in modulating viral clearance and resistance signaling. Methylation-expression relationships were found to be required for the immune response and cellular components are necessary for the proper defense response to infection and symptom recovery.

Overexpression of the wheat expansin gene TaEXPA2 improves oxidative stress tolerance in transgenic Arabidopsis plants

Expansins play an important role in plant stress tolerance. In a previous study, we cloned the wheat expansin gene TaEXPA2. Here, we analyze its involvement in oxidative stress tolerance. First, we observed that the expression of TaEXPA2 in wheat seedlings was upregulated during H₂O₂ stress. Then, we assembled a TaEXPA2 gene expression vector, transformed it to Arabidopsis, and obtained transgenic plants overexpressing TaEXPA2 (labeled OE). When exposed to H₂O₂, both OE and wild-type (Col) plants were damaged by oxidative stress, as indicated by decolored leaves and increased malondialdehyde (MDA) content. Damage in OE plants was less severe than in Col plants (WT), and this was accompanied by higher activity of cell wall peroxidase (POD) enzymes, including soluble POD, ionically bound POD, and covalently bound POD. The expansin activities of the OE plants were also higher than WT under oxidative stress. We further obtained the Arabidopsis mutant atexpa2 (AtEXPA2 is homologous to TaEXPA2), and found that the antioxidant ability of atexpa2 was lower than that in Col plants, accompanied by depressed activity of POD enzymes and expansins in cell walls. We transformed wheat TaEXPA2 to atexpa2 and obtained plants (labeled Rs) capable of recovering the antioxidant capacity. Oxidative stress tolerance in Rs plants was higher than that of Col plants, and the Rs plants also had higher levels of cell wall POD enzyme and expansin activity. Finally, we identified 13 POD genes in Arabidopsis thaliana and analyzed their expression patterns using quantitative real-time PCR. The expression of 4 of these genes (AtPOD31, AtPOD33, AtPOD34 and AtPOD71) was significantly upregulated during exposure to H₂O₂. We speculate that the 4 genes upregulated by H₂O₂ treatment are involved in the increased activity of POD in the cell wall. We suggest that TaEXPA2 may regulate antioxidant capacity in plants by regulating the activity of cell wall peroxidase.

Phytopathogen-induced changes to plant methylomes

DNA methylation is a dynamic and reversible type of epigenetic mark that contributes to cellular physiology by affecting transcription activity, transposon mobility and genome stability. When plants are infected with pathogens, plant DNA methylation patterns can change, indicating an epigenetic interplay between plant host and pathogen. In most cases methylation can change susceptibility. While DNA hypomethylation appears to be a common phenomenon during the susceptible interaction, the levels and patterns of hypomethylation in transposable elements and genic regions may mediate distinct responses against various plant pathogens. The effect of DNA methylation on the plant immune response and other cellular activities and molecular functions is established by localized differential DNA methylation via cis-regulatory mechanisms as well as through trans-acting mechanisms. Understanding the epigenetic differences that control the phenotypic variations between susceptible and resistant interactions should facilitate the identification of new sources of resistance mediated by epigenetic mechanisms, which can be exploited to endow pathogen resistance to crops.

Regulation of genome-wide DNA methylation by mobile small RNAs

Contents Summary 540 I. Introduction 540 II. There are different types of sRNA mobility 541 III. Mechanisms of sRNA movement 541 IV. Long-distance, shoot-root, mobile siRNAs influence DNA methylation in recipient tissues 541 V. Classes of interactions between shoot-root mobile siRNAs and DNA methylation 542 VI. Loci targeted directly and indirectly by shoot-root mobile siRNAs are associated with different histone modifications 543 VII. Is mobile siRNA-regulated DNA methylation important in specific tissues or under specific conditions? 543 VIII. Mobile sRNAs can be used to modify plant traits 544 IX. Conclusions 544 Acknowledgements 544 References 544 SUMMARY: RNA-directed DNA methylation (RdDM) at cytosine residues regulates gene expression, silences transposable elements and influences genome stability. The mechanisms responsible for RdDM are guided to target loci by small RNAs (sRNAs) that can move within plants cell to cell and long distance. Here we discuss recent advances in the understanding of interactions between mobile sRNAs and DNA methylation. We describe the mechanisms of sRNA movement, the differences between known classes of mobile sRNA-DNA methylation interactions and the limits of current knowledge. Finally, we discuss potential applications of mobile sRNAs in modifying plant traits.

Methods for discovering genomic loci exhibiting complex patterns of differential methylation

BACKGROUND

Cytosine methylation is widespread in most eukaryotic genomes and is known to play a substantial role in various regulatory pathways. Unmethylated cytosines may be converted to uracil through the addition of sodium bisulphite, allowing genome-wide quantification of cytosine methylation via high-throughput sequencing. The data thus acquired allows the discovery of methylation 'loci'; contiguous regions of methylation consistently methylated across biological replicates. The mapping of these loci allows for associations with other genomic factors to be identified, and for analyses of differential methylation to take place.

RESULTS

The segmentSeq R package is extended to identify methylation loci from high-throughput sequencing data from multiple experimental conditions. A statistical model is then developed that accounts for biological replication and variable rates of non-conversion of cytosines in each sample to compute posterior likelihoods of methylation at each locus within an empirical Bayesian framework. The same model is used as a basis for analysis of differential methylation between multiple experimental conditions with the baySeq R package. We demonstrate the capability of this method to analyse complex data sets in an analysis of data derived from multiple Dicer-like mutants in Arabidopsis. This reveals several novel behaviours at distinct sets of loci in response to loss of one or more of the Dicer-like proteins that indicate an antagonistic relationship between the Dicer-like proteins at at least some methylation loci. Finally, we show in simulation studies that this approach can be significantly more powerful in the detection of differential methylation than many existing methods in data derived from both mammalian and plant systems.

CONCLUSIONS

The methods developed here make it possible to analyse high-throughput sequencing of the methylome of any given organism under a diverse set of experimental conditions. The methods are able to identify methylation loci and evaluate the likelihood that a region is truly methylated under any given experimental condition, allowing for downstream analyses that characterise differences between methylated and non-methylated regions of the genome. Futhermore, diverse patterns of differential methylation may also be characterised from these data.

Extensive transcriptomic and epigenomic remodelling occurs during Arabidopsis thaliana germination

BACKGROUND

Seed germination involves progression from complete metabolic dormancy to a highly active, growing seedling. Many factors regulate germination and these interact extensively, forming a complex network of inputs that control the seed-to-seedling transition. Our understanding of the direct regulation of gene expression and the dynamic changes in the epigenome and small RNAs during germination is limited. The interactions between genome, transcriptome and epigenome must be revealed in order to identify the regulatory mechanisms that control seed germination.

RESULTS

We present an integrated analysis of high-resolution RNA sequencing, small RNA sequencing and MethylC sequencing over ten developmental time points in Arabidopsis thaliana seeds, finding extensive transcriptomic and epigenomic transformations associated with seed germination. We identify previously unannotated loci from which messenger RNAs are expressed transiently during germination and find widespread alternative splicing and divergent isoform abundance of genes involved in RNA processing and splicing. We generate the first dynamic transcription factor network model of germination, identifying known and novel regulatory factors. Expression of both microRNA and short interfering RNA loci changes significantly during germination, particularly between the seed and the post-germinative seedling. These are associated with changes in gene expression and large-scale demethylation observed towards the end of germination, as the epigenome transitions from an embryo-like to a vegetative seedling state.

CONCLUSIONS

This study reveals the complex dynamics and interactions of the transcriptome and epigenome during seed germination, including the extensive remodelling of the seed DNA methylome from an embryo-like to vegetative-like state during the seed-to-seedling transition. Data are available for exploration in a user-friendly browser at https://jbrowse.latrobe.edu.au/germination_epigenome .