Mutations in retrotransposon AtCOPIA4 compromises resistance to Hyaloperonospora parasitica in Arabidopsis thaliana

Long-read direct RNA sequencing reveals epigenetic regulation of chimeric gene-transposon transcripts in Arabidopsis thaliana

Transposable elements (TEs) are accumulated in both intergenic and intragenic regions in plant genomes. Intragenic TEs often act as regulatory elements of associated genes and are also co-transcribed with genes, generating chimeric TE-gene transcripts. Despite the potential impact on mRNA regulation and gene function, the prevalence and transcriptional regulation of TE-gene transcripts are poorly understood. By long-read direct RNA sequencing and a dedicated bioinformatics pipeline, ParasiTE, we investigated the transcription and RNA processing of TE-gene transcripts in Arabidopsis thaliana. We identified a global production of TE-gene transcripts in thousands of A. thaliana gene loci, with TE sequences often being associated with alternative transcription start sites or transcription termination sites. The epigenetic state of intragenic TEs affects RNAPII elongation and usage of alternative poly(A) signals within TE sequences, regulating alternative TE-gene isoform production. Co-transcription and inclusion of TE-derived sequences into gene transcripts impact regulation of RNA stability and environmental responses of some loci. Our study provides insights into TE-gene interactions that contributes to mRNA regulation, transcriptome diversity, and environmental responses in plants.

Impacts of DNA methylases and demethylases on the methylation and expression of Arabidopsis ethylene signal pathway genes

Arabidopsis ethylene (ET) signal pathway plays important roles in various aspects. Cytosine DNA methylation is significant in controlling gene expression in plants. Here, we analyzed the bisulfite sequencing and mRNA sequencing data from Arabidopsis (de)methylase mutants met1, cmt3, drm1/2, ddm1, ros1-4, and rdd to investigate how DNA (de)methylases influence the DNA methylation and expression of Arabidopsis ET pathway genes. At least 32 genes are found to involved in Arabidopsis ET pathway by text mining. Among them, 14 genes are unmethylated or methylated with very low levels. ACS6 and ACS9 are conspicuously methylated within their upstream regions. The other 16 genes are predominantly methylated at the CG sites within gene body regions in wild-type plants, and mutation of MET1 resulted in almost entire elimination of the CG methylations. In addition, CG methylations within some genes are jointly maintained by MET1 and other (de)methylases. Analyses of mRNA-seq data indicated that some ET pathway genes were differentially expressed between wild-type and diverse mutants. PDF1.2, the marker gene of ET signal pathway, was found being regulated indirectly by the methylases. Eighty-two transposable elements (TEs) were identified to be associated to 15 ET pathway genes. ACS11 is found located in a heterochromatin region that contains 57 TEs, indicating its specific expression and regulation. Together, our results suggest that DNA (de)methylases are implicated in the regulation of CG methylation within gene body regions and transcriptional activity of some ET pathway genes and that maintenance of normal CG methylation is essential for ET pathway in Arabidopsis.

Widespread premature transcription termination of Arabidopsis thaliana NLR genes by the spen protein FPA

Genes involved in disease resistance are some of the fastest evolving and most diverse components of genomes. Large numbers of nucleotide-binding, leucine-rich repeat (NLR) genes are found in plant genomes and are required for disease resistance. However, NLRs can trigger autoimmunity, disrupt beneficial microbiota or reduce fitness. It is therefore crucial to understand how NLRs are controlled. Here, we show that the RNA-binding protein FPA mediates widespread premature cleavage and polyadenylation of NLR transcripts, thereby controlling their functional expression and impacting immunity. Using long-read Nanopore direct RNA sequencing, we resolved the complexity of NLR transcript processing and gene annotation. Our results uncover a co-transcriptional layer of NLR control with implications for understanding the regulatory and evolutionary dynamics of NLRs in the immune responses of plants.

The Arabidopsis PHD-finger protein EDM2 has multiple roles in balancing NLR immune receptor gene expression

Plant NLR-type receptors serve as sensitive triggers of host immunity. Their expression has to be well-balanced, due to their interference with various cellular processes and dose-dependency of their defense-inducing activity. A genetic “arms race” with fast-evolving pathogenic microbes requires plants to constantly innovate their NLR repertoires. We previously showed that insertion of the COPIA-R7 retrotransposon into RPP7 co-opted the epigenetic transposon silencing signal H3K9me2 to a new function promoting expression of this Arabidopsis thaliana NLR gene. Recruitment of the histone binding protein EDM2 to COPIA-R7-associated H3K9me2 is required for optimal expression of RPP7. By profiling of genome-wide effects of EDM2, we now uncovered additional examples illustrating effects of transposons on NLR gene expression, strongly suggesting that these mobile elements can play critical roles in the rapid evolution of plant NLR genes by providing the “raw material” for gene expression mechanisms. We further found EDM2 to have a global role in NLR expression control. Besides serving as a positive regulator of RPP7 and a small number of other NLR genes, EDM2 acts as a suppressor of a multitude of additional NLR genes. We speculate that the dual functionality of EDM2 in NLR expression control arose from the need to compensate for fitness penalties caused by high expression of some NLR genes by suppression of others. Moreover, we are providing new insights into functional relationships of EDM2 with its interaction partner, the RNA binding protein EDM3/AIPP1, and its target gene IBM1, encoding an H3K9-demethylase.

We previously found the Arabidopsis thaliana PHD-finger protein EDM2 to serve as a chromatin-associated factor controlling expression of the NLR-type immune receptor gene RPP7. EDM2 binds to the transposon-silencing signal H3K9me2 and affects levels of this epigenetic mark at various loci. By genome-wide profiling of transcript- and H3K9me2-levels we now found EDM2 to have a broader role in controlling NLR gene expression. In order to mitigate fitness costs caused by its promoting effects on RPP7 expression and that of several other NLR genes, EDM2 seems to suppress expression of many additional members of this gene family. This observation is in line with multiple reports demonstrating the need for balanced expression of NLRs, which can substantially reduce overall plant fitness, but need to be present at certain minimal levels to confer sufficient immune protection. Our previous results demonstrated that the influence of EDM2 on RPP7 expression was co-opted to this immune receptor gene by the insertion of an EDM2-controlled transposon. Here, we are providing additional examples for transposon-associated effects on NLR gene expression, suggesting that these mobile elements play an important role for NLR genes by equipping members of this rapidly evolving gene family with regulatory mechanisms needed for balanced expression.

Epigenetic regulation of spurious transcription initiation in Arabidopsis

In plants, epigenetic regulation is critical for silencing transposons and maintaining proper gene expression. However, its impact on the genome-wide transcription initiation landscape remains elusive. By conducting a genome-wide analysis of transcription start sites (TSSs) using cap analysis of gene expression (CAGE) sequencing, we show that thousands of TSSs are exclusively activated in various epigenetic mutants of Arabidopsis thaliana and referred to as cryptic TSSs. Many have not been identified in previous studies, of which up to 65% are contributed by transposons. They possess similar genetic features to regular TSSs and their activation is strongly associated with the ectopic recruitment of RNAPII machinery. The activation of cryptic TSSs significantly alters transcription of nearby TSSs, including those of genes important for development and stress responses. Our study, therefore, sheds light on the role of epigenetic regulation in maintaining proper gene functions in plants by suppressing transcription from cryptic TSSs.

Epigenetic regulation can silence transposons and maintain gene expression. Here the authors survey Arabidopsis mutants defective in epigenetic regulation and show ectopic activation of thousands of cryptic TSSs and altered expression of nearby genes demonstrating the importance of suppressing spurious transcription.

Antagonistic Actions of FPA and IBM2 Regulate Transcript Processing from Genes Containing Heterochromatin

Repressive epigenetic marks, such as DNA and histone methylation, are sometimes located within introns. In Arabidopsis (Arabidopsis thaliana), INCREASE IN BONSAI METHYLATION2 (IBM2), an RNA-binding protein containing a bromo-adjacent homology domain, is required to process functional transcript isoforms of genes carrying intronic heterochromatin. In a genetic screen for suppressors of the ibm2 mutation, we identified FPA, an RNA-binding protein that promotes use of proximal polyadenylation sites in genes targeted by IBM2, including IBM1 encoding an essential H3K9 histone demethylase and the disease resistance gene RECOGNITION OF PERONOSPORA PARASITICA7 Both IBM2 and FPA are involved in the processing of their common mRNA targets: Transcription of IBM2 target genes is restored when FPA is mutated in ibm2 and impaired in transgenic plants overexpressing FPA By contrast, transposons targeted by IBM2 and localized outside introns are not under this antagonistic control. The DNA methylation patterns of some genes and transposons are modified in fpa plants, including the large intron of IBM1, but these changes are rather limited and reversed when the mutant is complemented, indicating that FPA has a restricted role in mediating silencing. These data reveal a complex regulation by IBM2 and FPA pathways in processing mRNAs of genes bearing heterochromatic marks.

siRNAs regulate DNA methylation and interfere with gene and lncRNA expression in the heterozygous polyploid switchgrass

BACKGROUND

Understanding the DNA methylome and its relationship with non-coding RNAs, including microRNAs (miRNAs) and long non-coding RNAs (lncRNAs), is essential for elucidating the molecular mechanisms underlying key biological processes in plants. Few studies have examined the functional roles of the DNA methylome in grass species with highly heterozygous polyploid genomes.

RESULTS

We performed genome-wide DNA methylation profiling in the tetraploid switchgrass (Panicum virgatum L.) cultivar 'Alamo' using bisulfite sequencing. Single-base-resolution methylation patterns were observed in switchgrass leaf and root tissues, which allowed for characterization of the relationship between DNA methylation and mRNA, miRNA, and lncRNA populations. The results of this study revealed that siRNAs positively regulate DNA methylation of the mCHH sites surrounding genes, and that DNA methylation interferes with gene and lncRNA expression in switchgrass. Ninety-six genes covered by differentially methylated regions (DMRs) were annotated by GO analysis as being involved in stimulus-related processes. Functionally, 82% (79/96) of these genes were found to be hypomethylated in switchgrass root tissue. Sequencing analysis of lncRNAs identified two lncRNAs that are potential precursors of miRNAs, which are predicted to target genes that function in cellulose biosynthesis, stress regulation, and stem and root development.

CONCLUSIONS

This study characterized the DNA methylome in switchgrass and elucidated its relevance to gene and non-coding RNAs. These results provide valuable genomic resources and references that will aid further epigenetic research in this important biofuel crop.

Dynamic Roles for Small RNAs and DNA Methylation during Ovule and Fiber Development in Allotetraploid Cotton

DNA methylation is essential for plant and animal development. In plants, methylation occurs at CG, CHG, and CHH (H = A, C or T) sites via distinct pathways. Cotton is an allotetraploid consisting of two progenitor genomes. Each cotton fiber is a rapidly-elongating cell derived from the ovule epidermis, but the molecular basis for this developmental transition is unknown. Here we analyzed methylome, transcriptome, and small RNAome and revealed distinct changes in CHH methylation during ovule and fiber development. In ovules, CHH hypermethylation in promoters correlated positively with siRNAs, inducing RNA-dependent DNA methylation (RdDM), and up-regulation of ovule-preferred genes. In fibers, the ovule-derived cells generated additional heterochromatic CHH hypermethylation independent of RdDM, which repressed transposable elements (TEs) and nearby genes including fiber-related genes. Furthermore, CHG and CHH methylation in genic regions contributed to homoeolog expression bias in ovules and fibers. Inhibiting DNA methylation using 5-aza-2'-deoxycytidine in cultured ovules has reduced fiber cell number and length, suggesting a potential role for DNA methylation in fiber development. Thus, RdDM-dependent methylation in promoters and RdDM-independent methylation in TEs and nearby genes could act as a double-lock feedback mechanism to mediate gene and TE expression, potentiating the transition from epidermal to fiber cells during ovule and seed development.

Cotton is the world’s largest source of renewable textile fiber and is an allotetraploid crop consisting of two progenitor genomes. In plants, de novo CHH (H = A, T, or C) methylation depends on RNA-directed DNA methylation (RdDM) and CHROMOMETHYLASE2 (CMT2)-mediated pathways. The biological significance of the two pathways is largely unknown. Here we show dynamic roles of these two pathways in ovule and fiber development. RdDM-dependent CHH methylation is linked to gene activation in ovules, and additional CMT2-dependent methylation leads to silencing of transposons and nearby genes in fibers. Moreover, DNA methylation affects expression bias of homoeologous genes and fiber development. These findings provide novel insights into epigenetic regulation of organ development and polyploid evolution.

Genome size diversity in angiosperms and its influence on gene space

Genome size varies c. 2400-fold in angiosperms (flowering plants), although the range of genome size is skewed towards small genomes, with a mean genome size of 1C=5.7Gb. One of the most crucial factors governing genome size in angiosperms is the relative amount and activity of repetitive elements. Recently, there have been new insights into how these repeats, previously discarded as 'junk' DNA, can have a significant impact on gene space (i.e. the part of the genome comprising all the genes and gene-related DNA). Here we review these new findings and explore in what ways genome size itself plays a role in influencing how repeats impact genome dynamics and gene space, including gene expression.

Epigenetic regulation and functional exaptation of transposable elements in higher plants

Transposable elements (TEs) are mobile genetic elements that can proliferate in their host genomes. Because of their robust amplification, TEs have long been considered 'selfish DNA', harmful insertions that can threaten host genome integrity. The idea of TEs as junk DNA comes from analysis of epigenetic silencing of their mobility in plants and animals. This idea contrasts with McClintock's characterization of TEs as 'controlling elements'. Emerging studies on the regulatory functions of TEs in plant genomes have updated McClintock's characterization, indicating exaptation of TEs for genetic regulation. In this review, we summarize recent progress in TE silencing, particularly in Arabidopsis and rice, and show that TEs provide an abundant, natural source of regulation for the host genome.

How important are transposons for plant evolution?

For decades, transposable elements have been known to produce a wide variety of changes in plant gene expression and function. This has led to the idea that transposable element activity has played a key part in adaptive plant evolution. This Review describes the kinds of changes that transposable elements can cause, discusses evidence that those changes have contributed to plant evolution and suggests future strategies for determining the extent to which these changes have in fact contributed to plant adaptation and evolution. Recent advances in genomics and phenomics for a range of plant species, particularly crops, have begun to allow the systematic assessment of these questions.

A systems approach for identifying resistance factors to Rice stripe virus

Rice stripe virus (RSV) causes disease that can severely affect the productivity of rice (Oryza sativa). Several RSV-resistant cultivars have been developed. However, host factors conferring RSV resistance in these cultivars are still elusive. Here, we present a systems approach for identifying potential rice resistance factors. We developed two near-isogenic lines (NIL), RSV-resistant NIL22 and RSV-susceptible NIL37, and performed gene expression profiling of the two lines in RSV-infected and RSV-uninfected conditions. We identified 237 differentially expressed genes (DEG) between NIL22 and NIL37. By integrating with known quantitative trait loci (QTL), we selected 11 DEG located within the RSV resistance QTL as RSV resistance factor candidates. Furthermore, we identified 417 DEG between RSV-infected and RSV-uninfected conditions. Using an interaction network-based method, we selected 20 DEG highly interacting with the two sets of DEG as RSV resistance factor candidates. Among the 31 candidates, we selected the final set of 21 potential RSV resistance factors whose differential expression was confirmed in the independent samples using real-time reverse-transcription polymerase chain reaction. Finally, we reconstructed a network model delineating potential association of the 21 selected factors with resistance-related processes. In summary, our approach, based on gene expression profiling, revealed potential host resistance factors and a network model describing their relationships with resistance-related processes, which can be further validated in detailed experiments.