Loss of RNA-dependent RNA polymerase 2 (RDR2) function causes widespread and unexpected changes in the expression of transposons, genes, and 24-nt small RNAs

RNA-directed DNA methylation mutants reduce histone methylation at the paramutated maize booster1 enhancer

Paramutation is the transfer of mitotically and meiotically heritable silencing information between two alleles. With paramutation at the maize (Zea mays) booster1 (b1) locus, the low-expressed B' epiallele heritably changes the high-expressed B-I epiallele into B' with 100% frequency. This requires specific tandem repeats and multiple components of the RNA-directed DNA methylation pathway, including the RNA-dependent RNA polymerase (encoded by mediator of paramutation1, mop1), the second-largest subunit of RNA polymerase IV and V (NRP(D/E)2a, encoded by mop2), and the largest subunit of RNA Polymerase IV (NRPD1, encoded by mop3). Mutations in mop genes prevent paramutation and release silencing at the B' epiallele. In this study, we investigated the effect of mutations in mop1, mop2, and mop3 on chromatin structure and DNA methylation at the B' epiallele, and especially the regulatory hepta-repeat 100 kb upstream of the b1 gene. Mutations in mop1 and mop3 resulted in decreased repressive histone modifications H3K9me2 and H3K27me2 at the hepta-repeat. Associated with this decrease were partial activation of the hepta-repeat enhancer function, formation of a multi-loop structure, and elevated b1 expression. In mop2 mutants, which do not show elevated b1 expression, H3K9me2, H3K27me2 and a single-loop structure like in wild-type B' were retained. Surprisingly, high CG and CHG methylation levels at the B' hepta-repeat remained in all three mutants, and CHH methylation was low in both wild type and mutants. Our results raise the possibility of MOP factors mediating RNA-directed histone methylation rather than RNA-directed DNA methylation at the b1 locus.

Potent pollen gene regulation by DNA glycosylases in maize

Although DNA methylation primarily represses transposable elements (TEs) in plants, it also represses select endosperm and pollen genes. These genes, or their cis-regulatory elements, are methylated in plant body tissues but are demethylated by DNA glycosylases (DNGs) in endosperm and pollen, enabling their transcription. Activity of either one of two DNGs, MDR1 or DNG102, is essential for pollen viability in maize. Using single-pollen mRNA sequencing on pollen segregating mutations in both genes, we identified 58 candidate DNG target genes, whose expression is strongly decreased in double mutant pollen (124-fold decrease on average). These genes account for 11.1% of the wild-type pollen polyadenylated transcriptome, but they are silent or barely detectable in the plant body. They are unusual in their tendency to lack introns but even more so in their having TE-like methylation in their coding DNA sequence. Moreover, they are strongly enriched for predicted functions in cell wall modification. While some may support development of the pollen grain cell wall, expansins and pectinases in this set of genes suggest a function in cell wall loosening to support the rapid tip growth characteristic of pollen tubes as they carry the sperm cells through maternal apoplast and extracellular matrix of the pistil. These results suggest a critical role for DNA methylation and demethylation in regulating maize genes with potential for extremely high expression in pollen but constitutive silencing elsewhere.

Loops, crosstalk, and compartmentalization: it takes many layers to regulate DNA methylation

DNA methylation is a conserved epigenetic modification associated with transposon silencing and gene regulation. The stability of this modification relies on intimate connections between DNA and histone modifications that generate self-reinforcing loops wherein the presence of one mark promotes the other. However, it is becoming increasingly clear that the efficiency of these loops is affected by cross-talk between pathways and by chromatin accessibility, which is heavily influenced by histone variants. Focusing primarily on plants, this review provides an update on the aforementioned self-reinforcing loops, highlights recent advances in understanding how DNA methylation pathways are restricted to prevent encroachment on genes, and discusses the roles of histone variants in compartmentalizing epigenetic pathways within the genome. This multilayered approach facilitates two essential, yet opposing functions, the ability to maintain heritable DNA methylation patterns while retaining the flexibility to modify these patterns during development.

Miniature Inverted-Repeat Transposable Elements: Small DNA Transposons That Have Contributed to Plant MICRORNA Gene Evolution

Angiosperms form the largest phylum within the Plantae kingdom and show remarkable genetic variation due to the considerable difference in the nuclear genome size of each species. Transposable elements (TEs), mobile DNA sequences that can amplify and change their chromosome position, account for much of the difference in nuclear genome size between individual angiosperm species. Considering the dramatic consequences of TE movement, including the complete loss of gene function, it is unsurprising that the angiosperms have developed elegant molecular strategies to control TE amplification and movement. Specifically, the RNA-directed DNA methylation (RdDM) pathway, directed by the repeat-associated small-interfering RNA (rasiRNA) class of small regulatory RNA, forms the primary line of defense to control TE activity in the angiosperms. However, the miniature inverted-repeat transposable element (MITE) species of TE has at times avoided the repressive effects imposed by the rasiRNA-directed RdDM pathway. MITE proliferation in angiosperm nuclear genomes is due to their preference to transpose within gene-rich regions, a pattern of transposition that has enabled MITEs to gain further transcriptional activity. The sequence-based properties of a MITE results in the synthesis of a noncoding RNA (ncRNA), which, after transcription, folds to form a structure that closely resembles those of the precursor transcripts of the microRNA (miRNA) class of small regulatory RNA. This shared folding structure results in a MITE-derived miRNA being processed from the MITE-transcribed ncRNA, and post-maturation, the MITE-derived miRNA can be used by the core protein machinery of the miRNA pathway to regulate the expression of protein-coding genes that harbor homologous MITE insertions. Here, we outline the considerable contribution that the MITE species of TE have made to expanding the miRNA repertoire of the angiosperms.

Genome-wide analysis of epigenetic and transcriptional changes in the pathogenesis of RGSV in rice

Rice grassy stunt virus (RGSV), a typical negative single-stranded RNA virus, invades rice and generates several disease signs, including dwarfing, tillering, and sterility. Previous research has revealed that RGSV-encoded proteins can force the host's ubiquitin-proteasome system to utilize them for viral pathogenesis. However, most of the studies were limited to a single omics level and lacked multidimensional data collection and correlation analysis on the mechanisms of RGSV-rice interactions. Here, we performed a comprehensive association analysis of genome-wide methylation sequencing, transcriptome sequencing, and histone H3K9me3 modification in RGSV-infested as well as non-infested rice leaves, and the levels of all three cytosine contexts (CG, CHG and CHH) were found to be slightly lower in RGSV-infected rice leaves than in normal rice. Large proportions of DMRs were distributed in the promoter and intergenic regions, and most DMRs were enriched in the CHH context, where the number of CHH hypo-DMRs was almost twice as high as that of hyper-DMRs. Among the genes with down-regulated expression and hypermethylation, we analyzed and identified 11 transcripts involved in fertility, plant height and tillering, and among the transcribed up-regulated and hypermethylated genes, we excavated 7 transcripts related to fertility, plant height and tillering. By analyzing the changes of histone H3K9me3 modification before and after virus infestation, we found that the distribution of H3K9me3 modification in the whole rice genome was prevalent, mainly concentrated in the gene promoter and gene body regions, which was distinctly different from the characteristics of animals. Combined with transcriptomic data, H3K9me3 mark was found to favor targeting highly expressed genes. After RGSV infection, H3K9me3 modifications in several regions of CTK and BR hormone signaling-related genes were altered, providing important targets for subsequent studies.

Soaud S, Huang Q, M. Saleh A, A. S. Abourehab M, Wan L, Cheng GT, Liu J, Ihtisham M, Noor Z, Rouf Mir R, Zhao X, Yan K, Abbas M, Li J. Natural resistance of tomato plants to Tomato yellow leaf curl virus

Tomato yellow leaf curl virus (TYLCV) is one of the most harmful afflictions in the world that affects tomato growth and production. Six regular antagonistic genes (Ty-1, Ty-2, Ty-3, Ty-4, ty-5, and Ty-6) have been transferred from wild germplasms to commercial cultivars as TYLCV protections. With Ty-1 serving as an appropriate source of TYLCV resistance, only Ty-1, Ty-2, and Ty-3 displayed substantial levels of opposition in a few strains. It has been possible to clone three TYLCV opposition genes (Ty-1/Ty-3, Ty-2, and ty-5) that target three antiviral safety mechanisms. However, it significantly impacts obtaining permanent resistance to TYLCV, trying to maintain opposition whenever possible, and spreading opposition globally. Utilizing novel methods, such as using resistance genes and identifying new resistance resources, protects against TYLCV in tomato production. To facilitate the breeders make an informed decision and testing methods for TYLCV blockage, this study highlights the portrayal of typical obstruction genes, common opposition sources, and subatomic indicators. The main goal is to provide a fictitious starting point for the identification and application of resistance genes as well as the maturation of tomato varieties that are TYLCV-resistant.

Retrotransposons: How the continuous evolutionary front shapes plant genomes for response to heat stress

Long terminal repeat retrotransposons (LTR retrotransposons) are the most abundant group of mobile genetic elements in eukaryotic genomes and are essential in organizing genomic architecture and phenotypic variations. The diverse families of retrotransposons are related to retroviruses. As retrotransposable elements are dispersed and ubiquitous, their "copy-out and paste-in" life cycle of replicative transposition leads to new genome insertions without the excision of the original element. The overall structure of retrotransposons and the domains responsible for the various phases of their replication is highly conserved in all eukaryotes. The two major superfamilies of LTR retrotransposons, Ty1/Copia and Ty3/Gypsy, are distinguished and dispersed across the chromosomes of higher plants. Members of these superfamilies can increase in copy number and are often activated by various biotic and abiotic stresses due to retrotransposition bursts. LTR retrotransposons are important drivers of species diversity and exhibit great variety in structure, size, and mechanisms of transposition, making them important putative actors in genome evolution. Additionally, LTR retrotransposons influence the gene expression patterns of adjacent genes by modulating potential small interfering RNA (siRNA) and RNA-directed DNA methylation (RdDM) pathways. Furthermore, comparative and evolutionary analysis of the most important crop genome sequences and advanced technologies have elucidated the epigenetics and structural and functional modifications driven by LTR retrotransposon during speciation. However, mechanistic insights into LTR retrotransposons remain obscure in plant development due to a lack of advancement in high throughput technologies. In this review, we focus on the key role of LTR retrotransposons response in plants during heat stress, the role of centromeric LTR retrotransposons, and the role of LTR retrotransposon markers in genome expression and evolution.

Active DNA demethylation in plants: 20 years of discovery and beyond

Maintaining proper DNA methylation levels in the genome requires active demethylation of DNA. However, removing the methyl group from a modified cytosine is chemically difficult and therefore, the underlying mechanism of demethylation had remained unclear for many years. The discovery of the first eukaryotic DNA demethylase, Arabidopsis thaliana REPRESSOR OF SILENCING 1 (ROS1), led to elucidation of the 5-methylcytosine base excision repair mechanism of active DNA demethylation. In the 20 years since ROS1 was discovered, our understanding of this active DNA demethylation pathway, as well as its regulation and biological functions in plants, has greatly expanded. These exciting developments have laid the groundwork for further dissecting the regulatory mechanisms of active DNA demethylation, with potential applications in epigenome editing to facilitate crop breeding and gene therapy.

Transposable Elements: Distribution, Polymorphism, and Climate Adaptation in Populus

Transposable elements (TEs) are a class of mobile genetic elements that make effects on shaping rapid phenotypic traits of adaptive significance. TE insertions are usually related to transcription changes of nearby genes, and thus may be subjected to purifying selection. Based on the available genome resources of Populus, we found that the composition of Helitron DNA family were highly variable and could directly influence the transcription of nearby gene expression, which are involving in stress-responsive, programmed cell death, and apoptosis pathway. Next, we indicated TEs are highly enriched in Populus trichocarpa compared with three other congeneric poplar species, especially located at untranslated regions (3'UTRs and 5'UTRs) and Helitron transposons, particularly 24-nt siRNA-targeted, are significantly associated with reduced gene expression. Additionally, we scanned a representative resequenced Populus tomentosa population, and identified 9,680 polymorphic TEs loci. More importantly, we identified a Helitron transposon located at the 3'UTR, which could reduce WRKY18 expression level. Our results highlight the importance of TE insertion events, which could regulate gene expression and drive adaptive phenotypic variation in Populus.

N. S, Pachamuthu K, Chenna S, Nair A, Shivaprasad PV. Essential role of γ-clade RNA-dependent RNA polymerases in rice development and yield-related traits is linked to their atypical polymerase activities regulating specific genomic regions

RNA-dependent RNA polymerases (RDR) generate double-stranded (ds)RNA triggers for RNA silencing across eukaryotes. Among the three clades, α-clade and β-clade members are key components of RNA silencing and mediators of stress responses across eukaryotes. However, γ-clade members are unusual in that they are represented in phylogenetically distant plants and fungi, and their functions are unknown. Using genetic, bioinformatic and biochemical methods, we show that γ-clade RDRs from Oryza sativa L. are involved in plant development as well as regulation of expression of coding and noncoding RNAs. Overexpression of γ-clade RDRs in transgenic rice and tobacco plants resulted in robust growth phenotype, whereas their silencing in rice displayed strong inhibition of growth. Small (s)RNA and RNA-seq analysis of OsRDR3 mis-expression lines suggested that it is specifically involved in the regulation of repeat-rich regions in the genome. Biochemical analysis confirmed that OsRDR3 has robust polymerase activities on both single stranded (ss)RNA and ssDNA templates similar to the activities reported for α-clade RDRs such as AtRDR6. Our results provide the first evidence of the importance of γ-clade RDRs in plant development, their atypical biochemical activities and their contribution to the regulation of gene expression.

Patterns of stability and change in the maize genome: a case study of small RNA transcriptomes in two recombinant inbred lines and their progenitors

Small RNAs (sRNAs) are epigenetic regulators of eukaryotic genes and transposable elements (TEs). Diverse sRNA expression patterns exist within a species, but how this diversity arises is not well understood. To provide a window into the dynamics of maize sRNA patterning, sRNA and mRNA transcriptomes were examined in two related Zea mays recombinant inbred lines (RILs) and their inbred parents. Analysis of these RILs revealed that most clusters of sRNA expression retained the parental sRNA expression level. However, expression states that differ from the parental allele were also observed, predominantly reflecting decreases in sRNA expression. When RIL sRNA expression differed from the parental allele, the new state was frequently similar between the two RILs, and similar to the expression state found at the allele in the other parent. Novel sRNA expression patterns, distinct from those of either parent, were rare. Additionally, examination of sRNA expression over TEs revealed one TE family, Gyma, which showed consistent enrichment for RIL sRNA expression differences compared to those found in parental alleles. These findings provide insights into how sRNA silencing might evolve over generations and suggest that further investigation into the molecular nature of sRNA trans regulators is warranted.

RNA-directed DNA methylation prevents rapid and heritable reversal of transposon silencing under heat stress in Zea mays

In large complex plant genomes, RNA-directed DNA methylation (RdDM) ensures that epigenetic silencing is maintained at the boundary between genes and flanking transposable elements. In maize, RdDM is dependent on Mediator of Paramutation1 (Mop1), a gene encoding a putative RNA dependent RNA polymerase. Here we show that although RdDM is essential for the maintenance of DNA methylation of a silenced MuDR transposon in maize, a loss of that methylation does not result in a restoration of activity. Instead, heritable maintenance of silencing is maintained by histone modifications. At one terminal inverted repeat (TIR) of this element, heritable silencing is mediated via histone H3 lysine 9 dimethylation (H3K9me2), and histone H3 lysine 27 dimethylation (H3K27me2), even in the absence of DNA methylation. At the second TIR, heritable silencing is mediated by histone H3 lysine 27 trimethylation (H3K27me3), a mark normally associated with somatically inherited gene silencing. We find that a brief exposure of high temperature in a mop1 mutant rapidly reverses both of these modifications in conjunction with a loss of transcriptional silencing. These reversals are heritable, even in mop1 wild-type progeny in which methylation is restored at both TIRs. These observations suggest that DNA methylation is neither necessary to maintain silencing, nor is it sufficient to initiate silencing once has been reversed. However, given that heritable reactivation only occurs in a mop1 mutant background, these observations suggest that DNA methylation is required to buffer the effects of environmental stress on transposable elements.

Most plant genomes are mostly transposable elements (TEs), most of which are held in check by modifications of both DNA and histones. The bulk of silenced TEs are associated with methylated DNA and histone H3 lysine 9 dimethylation (H3K9me2). In contrast, epigenetically silenced genes are often associated with histone lysine 27 trimethylation (H3K27me3). Although stress can affect each of these modifications, plants are generally competent to rapidly reset them following that stress. Here we demonstrate that although DNA methylation is not required to maintain silencing of the MuDR element, it is essential for preventing heat-induced, stable and heritable changes in both H3K9me2 and H3K27me3 at this element, and for concomitant changes in transcriptional activity. These finding suggest that RdDM acts to buffer the effects of heat on silenced transposable elements, and that a loss of DNA methylation under conditions of stress can have profound and long-lasting effects on epigenetic silencing in maize.

Locus-specific paramutation in Zea mays is maintained by a PICKLE-like chromodomain helicase DNA-binding 3 protein controlling development and male gametophyte function

Paramutations represent directed and meiotically-heritable changes in gene regulation leading to apparent violations of Mendelian inheritance. Although the mechanism and evolutionary importance of paramutation behaviors remain largely unknown, genetic screens in maize (Zea mays) identify five components affecting 24 nucleotide RNA biogenesis as required to maintain repression of a paramutant purple plant1 (pl1) allele. Currently, the RNA polymerase IV largest subunit represents the only component also specifying proper development. Here we identify a chromodomain helicase DNA-binding 3 (CHD3) protein orthologous to Arabidopsis (Arabidopsis thaliana) PICKLE as another component maintaining both pl1 paramutation and normal somatic development but without affecting overall small RNA biogenesis. In addition, genetic tests show this protein contributes to proper male gametophyte function. The similar mutant phenotypes documented in Arabidopsis and maize implicate some evolutionarily-conserved gene regulation while developmental defects associated with the two paramutation mutants are largely distinct. Our results show that a CHD3 protein responsible for normal plant ontogeny and sperm transmission also helps maintain meiotically-heritable epigenetic regulatory variation for specific alleles. This finding implicates an intersection of RNA polymerase IV function and nucleosome positioning in the paramutation process.

Genes are switched “on” and “off” during normal development by regulating DNA accessibility within the chromosomes. How certain gene variants permanently maintain “off” states from one generation to the next remains unclear, but studies in multiple eukaryotes implicate roles for specific types of small RNAs, some of which define cytosine methylation patterns. In corn, these RNAs come from at least two RNA polymerase II-derived complexes sharing a common catalytic subunit (RPD1). Although RPD1 both controls the normal developmental switching of many genes and permanently maintains some of these “off” states across generations, how RPD1 function defines heritable DNA accessibility is unknown. We discovered that a protein (CHD3a) belonging to a group known to alter nucleosome positioning is also required to help maintain a heritable “off” state for one particular corn gene variant controlling both plant and flower color. We also found CHD3a necessary for normal plant development and sperm transmission consistent with the idea that proper nucleosome positioning defines evolutionarily-important gene expression patterns. Because both CHD3a and RPD1 maintain the heritable “off” state of a specific gene variant, their functions appear to be mechanistically linked.

The dynamic transcriptome of waxy maize (Zea mays L. sinensis Kulesh) during seed development

BACKGROUND

Waxy maize (Zea mays L. sinensis Kulesh) is a mutant of maize (Zea mays L.) with a mutation at Waxy1 (Wx1) gene locus. The seed of waxy maize has higher viscosity compared to regular maize. By now, we know little about the expression patterns of genes that involved in the seed development of waxy maize.

OBJECTIVE

By analyzing the transcriptome data during waxy maize seed development, we attempt to dig out the genes that may influence the seed development of waxy maize.

METHODS

The seeds of waxy maize inbred line SWL01 from six phases after pollination were used to do RNA-seq. Bioinformatics methods were used to analyze the expression patterns of the expressed genes, to identify the genes involved in waxy maize seed development.

RESULTS

A total of 24,546 genes including 1611 transcription factors (TFs) were detected during waxy maize seed development. Coexpression analysis of expressed genes revealed the dynamic processes of waxy maize seed development. Particularly, 2457 genes including 177 TFs were specially expressed in waxy maize seed, some of which mainly involved in the process of seed dormancy and maturation. In addition, 2681, 5686, 4491, 4386, 3669 and 4624 genes were identified to be differential expressed genes (DEGs) at six phases compared to regular maize B73, and 113 DEGs among them may be key genes that lead the difference of seed development between waxy and regular maizes in milk stage.

CONCLUSION

In summary, we elucidated the expression patterns of expressed genes during waxy maize seed development globally. A series of genes that associated with seed development were identified in our research, which may provide an important resource for functional study of waxy maize seed development to help molecular assisted breeding.

Epigenetic Regulation of ABA-Induced Transcriptional Responses in Maize

Plants are subjected to extreme environmental conditions and must adapt rapidly. The phytohormone abscisic acid (ABA) accumulates during abiotic stress, signaling transcriptional changes that trigger physiological responses. Epigenetic modifications often facilitate transcription, particularly at genes exhibiting temporal, tissue-specific and environmentally-induced expression. In maize (Zea mays), MEDIATOR OF PARAMUTATION 1 (MOP1) is required for progression of an RNA-dependent epigenetic pathway that regulates transcriptional silencing of loci genomewide. MOP1 function has been previously correlated with genomic regions adjoining particular types of transposable elements and genic regions, suggesting that this regulatory pathway functions to maintain distinct transcriptional activities within genomic spaces, and that loss of MOP1 may modify the responsiveness of some loci to other regulatory pathways. As critical regulators of gene expression, MOP1 and ABA pathways each regulate specific genes. To determine whether loss of MOP1 impacts ABA-responsive gene expression in maize, mop1-1 and Mop1 homozygous seedlings were subjected to exogenous ABA and RNA-sequencing. A total of 3,242 differentially expressed genes (DEGs) were identified in four pairwise comparisons. Overall, ABA-induced changes in gene expression were enhanced in mop1-1 homozygous plants. The highest number of DEGs were identified in ABA-induced mop1-1 mutants, including many transcription factors; this suggests combinatorial regulatory scenarios including direct and indirect transcriptional responses to genetic disruption (mop1-1) and/or stimulus-induction of a hierarchical, cascading network of responsive genes. Additionally, a modest increase in CHH methylation at putative MOP1-RdDM loci in response to ABA was observed in some genotypes, suggesting that epigenetic variation might influence environmentally-induced transcriptional responses in maize.

Extensive profiling in Arabidopsis reveals abundant polysome-associated 24-nt small RNAs including AGO5-dependent pseudogene-derived siRNAs

In a reductionist perspective, plant silencing small (s)RNAs are often classified as mediating nuclear transcriptional gene silencing (TGS) or cytosolic posttranscriptional gene silencing (PTGS). Among the PTGS diagnostics is the association of AGOs and their sRNA cargos with the translation apparatus. In Arabidopsis, this is observed for AGO1 loaded with micro(mi)RNAs and, accordingly, translational-repression (TR) is one layer of plant miRNA action. Using AGO1:miRNA-mediated TR as a paradigm, we explored, with two unrelated polysome-isolation methods, which, among the ten Arabidopsis AGOs and numerous sRNA classes, interact with translation. We found that representatives of all three AGO-clades associate with polysomes, including the TGS-effector AGO4 and stereotypical 24-nt sRNAs that normally mediate TGS of transposons/repeats. Strikingly, approximately half of these annotated 24-nt siRNAs displayed unique matches in coding regions/introns of genes, and in pseudogenes, but not in transposons/repeats commonly found in their vicinity. Protein-coding gene-derived 24-nt sRNAs correlate with gene-body methylation. Those derived from pseudogenes belong to two main clusters defined by their parental-gene expression patterns, and are vastly enriched in AGO5, itself found on polysomes. Based on their tight expression pattern in developing and mature siliques, their biogenesis, and genomic/epigenomic features of their loci-of-origin, we discuss potential roles for these hitherto unknown polysome-enriched, pseudogene-derived siRNAs.

Genetic diversity for nitrogen use efficiency in Arabidopsis thaliana

The plasticity of plant growth response to differing nitrate availability renders the identification of biomarkers difficult, but allows access to genetic factors as tools to modulate root systems to a wide range of soil conditions. Nitrogen availability is a major determinant of crop yield. While the application of fertiliser substantially increases the yield on poor soils, it also causes nitrate pollution of water resources and high costs for farmers. Increasing nitrogen use efficiency in crop plants is a necessary step to implement low-input agricultural systems. We exploited the genetic diversity present in the worldwide Arabidopsis thaliana population to study adaptive growth patterns and changes in gene expression associated with chronic low nitrate stress, to identify biomarkers associated with good plant performance under low nitrate availability. Arabidopsis accessions were grown on agar plates with limited and sufficient supply of nitrate to measure root system architecture as well as shoot and root fresh weight. Differential gene expression was determined using Affymetrix ATH1 arrays. We show that the response to differing nitrate availability is highly variable in Arabidopsis accessions. Analyses of vegetative shoot growth and root system architecture identified accession-specific reaction modes to cope with limited nitrate availability. Transcription and epigenetic factors were identified as important players in the adaption to limited nitrogen in a global gene expression analysis. Five nitrate-responsive genes emerged as possible biomarkers for NUE in Arabidopsis. The plasticity of plant growth in response to differing nitrate availability in the substrate renders the identification of morphological and molecular features as biomarkers difficult, but at the same time allows access to a multitude of genetic factors which can be used as tools to modulate and adjust root systems to a wide range of soil conditions.

High Temporal-Resolution Transcriptome Landscape of Early Maize Seed Development

The early maize (Zea mays) seed undergoes several developmental stages after double fertilization to become fully differentiated within a short period of time, but the genetic control of this highly dynamic and complex developmental process remains largely unknown. Here, we report a high temporal-resolution investigation of transcriptomes using 31 samples collected at an interval of 4 or 6 h within the first six days of seed development. These time-course transcriptomes were clearly separated into four distinct groups corresponding to the stages of double fertilization, coenocyte formation, cellularization, and differentiation. A total of 22,790 expressed genes including 1415 transcription factors (TFs) were detected in early stages of maize seed development. In particular, 1093 genes including 110 TFs were specifically expressed in the seed and displayed high temporal specificity by expressing only in particular period of early seed development. There were 160, 22, 112, and 569 seed-specific genes predominantly expressed in the first 16 h after pollination, coenocyte formation, cellularization, and differentiation stage, respectively. In addition, network analysis predicted 31,256 interactions among 1317 TFs and 14,540 genes. The high temporal-resolution transcriptome atlas reported here provides an important resource for future functional study to unravel the genetic control of seed development.

RNA Interference: A Natural Immune System of Plants to Counteract Biotic Stressors

During plant-pathogen interactions, plants have to defend the living transposable elements from pathogens. In response to such elements, plants activate a variety of defense mechanisms to counteract the aggressiveness of biotic stressors. RNA interference (RNAi) is a key biological process in plants to inhibit gene expression both transcriptionally and post-transcriptionally, using three different groups of proteins to resist the virulence of pathogens. However, pathogens trigger an anti-silencing mechanism through the expression of suppressors to block host RNAi. The disruption of the silencing mechanism is a virulence strategy of pathogens to promote infection in the invaded hosts. In this review, we summarize the RNA silencing pathway, anti-silencing suppressors, and counter-defenses of plants to viral, fungal, and bacterial pathogens.

Disrupted Genome Methylation in Response to High Temperature Has Distinct Affects on Microspore Abortion and Anther Indehiscence

High-temperature (HT) stress induces male sterility, leading to yield reductions in crops. DNA methylation regulates a range of processes involved in plant development and stress responses, but its role in male sterility under HT remains unknown. Here, we investigated DNA methylation levels in cotton (Gossypium hirsutum) anthers under HT and normal temperature (NT) conditions by performing whole-genome bisulfite sequencing to investigate the regulatory roles of DNA methylation in male fertility under HT. Global disruption of DNA methylation, especially CHH methylation (where H = A, C, or T), was detected in an HT-sensitive line. Changes in the levels of 24-nucleotide small-interfering RNAs were significantly associated with DNA methylation levels. Experimental suppression of DNA methylation led to pollen sterility in the HT-sensitive line under NT conditions but did not affect the normal dehiscence of anther walls. Further transcriptome analysis showed that the expression of genes in sugar and reactive oxygen species (ROS) metabolic pathways were significantly modulated in anthers under HT, but auxin biosynthesis and signaling pathways were only slightly altered, indicating that HT disturbs sugar and ROS metabolism via disrupting DNA methylation, leading to microspore sterility. This study opens up a pathway for creating HT-tolerant cultivars using epigenetic techniques.

Redistribution of CHH Methylation and Small Interfering RNAs across the Genome of Tomato ddm1 Mutants

In plants, cytosine methylation, an epigenetic mark critical for transposon silencing, is maintained over generations by key enzymes that directly methylate DNA and is facilitated by chromatin remodelers, like DECREASE IN DNA METHYLATION1 (DDM1). Short-interfering RNAs (siRNAs) also mediate transposon DNA methylation through a process called RNA-directed DNA methylation (RdDM). In tomato (Solanum lycopersicum), siRNAs are primarily mapped to gene-rich chromosome arms, and not to pericentromeric regions as in Arabidopsis thaliana Tomato encodes two DDM1 genes. To better understand their functions and interaction with the RdDM pathway, we targeted the corresponding genes via the CRISPR/Cas9 technology, resulting in the isolation of Slddm1a and Slddm1b knockout mutants. Unlike the single mutants, Slddm1a Slddm1b double mutant plants display pleiotropic vegetative and reproductive phenotypes, associated with severe hypomethylation of the heterochromatic transposons in both the CG and CHG methylation contexts. The methylation in the CHH context increased for some heterochromatic transposons and conversely decreased for others localized in euchromatin. We found that the number of heterochromatin-associated siRNAs, including RdDM-specific small RNAs, increased significantly, likely limiting the transcriptional reactivation of transposons in Slddm1a Slddm1b Taken together, we propose that the global production of siRNAs and the CHH methylation mediated by the RdDM pathway are restricted to chromosome arms in tomato. Our data suggest that both pathways are greatly enhanced in heterochromatin when DDM1 functions are lost, at the expense of silencing mechanisms normally occurring in euchromatin.

Subtle Perturbations of the Maize Methylome Reveal Genes and Transposons Silenced by Chromomethylase or RNA-Directed DNA Methylation Pathways

DNA methylation is a chromatin modification that can provide epigenetic regulation of gene and transposon expression. Plants utilize several pathways to establish and maintain DNA methylation in specific sequence contexts. The chromomethylase (CMT) genes maintain CHG (where H = A, C or T) methylation. The RNA-directed DNA methylation (RdDM) pathway is important for CHH methylation. Transcriptome analysis was performed in a collection of Zea mays lines carrying mutant alleles for CMT or RdDM-associated genes. While the majority of the transcriptome was not affected, we identified sets of genes and transposon families sensitive to context-specific decreases in DNA methylation in mutant lines. Many of the genes that are up-regulated in CMT mutant lines have high levels of CHG methylation, while genes that are differentially expressed in RdDM mutants are enriched for having nearby mCHH islands, implicating context-specific DNA methylation in the regulation of expression for a small number of genes. Many genes regulated by CMTs exhibit natural variation for DNA methylation and transcript abundance in a panel of diverse inbred lines. Transposon families with differential expression in the mutant genotypes show few defining features, though several families up-regulated in RdDM mutants show enriched expression in endosperm tissue, highlighting the potential importance for this pathway during reproduction. Taken together, our findings suggest that while the number of genes and transposon families whose expression is reproducibly affected by mild perturbations in context-specific methylation is small, there are distinct patterns for loci impacted by RdDM and CMT mutants.

Systematic analysis of DEMETER-like DNA glycosylase genes shows lineage-specific Smi-miR7972 involved in SmDML1 regulation in Salvia miltiorrhiza

DEMETER-like DNA glycosylases (DMLs) initiate the base excision repair-dependent DNA demethylation to regulate a wide range of biological processes in plants. Six putative SmDML genes, termed SmDML1-SmDML6, were identified from the genome of S. miltiorrhiza, an emerging model plant for Traditional Chinese Medicine (TCM) studies. Integrated analysis of gene structures, sequence features, conserved domains and motifs, phylogenetic analysis and differential expression showed the conservation and divergence of SmDMLs. SmDML1, SmDML2 and SmDML4 were significantly down-regulated by the treatment of 5Aza-dC, a general DNA methylation inhibitor, suggesting involvement of SmDMLs in genome DNA methylation change. SmDML1 was predicted and experimentally validated to be target of Smi-miR7972. Computational analysis of forty whole genome sequences and almost all of RNA-seq data from Lamiids revealed that MIR7972s were only distributed in some plants of the three orders, including Lamiales, Solanales and Boraginales, and the number of MIR7972 genes varied among species. It suggests that MIR7972 genes underwent expansion and loss during the evolution of some Lamiids species. Phylogenetic analysis of MIR7972s showed closer evolutionary relationships between MIR7972s in Boraginales and Solanales in comparison with Lamiales. These results provide a valuable resource for elucidating DNA demethylation mechanism in S. miltiorrhiza.

Small RNA differential expression and regulation in Tuxpeño maize embryogenic callus induction and establishment

Somatic embryogenesis represents an alternative developmental process used to achieve genetic transformation and to approach key questions in maize development. It is known that embryogenic callus induction and plant regeneration are accompanied by microRNA expression changes. However, small RNA (sRNA) populations have not been explored during the proliferative callus subculture establishment and their impact on maintaining the dedifferentiated status and embryogenic potential is far from being completely understood. Here we globally tested the sRNA populations in explants (immature embryos), induced and established maize embryogenic callus from the Mexican cultivar VS-535, Tuxpeño landrace. We detected readjustments in 24 nt and 21-22 nt sRNAs during the embryogenic callus (EC) establishment and maintenance. A follow up on specific microRNAs (miRNAs) indicated that miRNAs related to stress response substantially increase upon the callus proliferation establishment, correlating with a reduction in some of their target levels. On the other hand, while 24 nt-long heterochromatic small interfering RNAs (hc-siRNAs) derived from transposable retroelements transiently decreased in abundance during the EC establishment, a population of 22 nt-hc-siRNAs increased. This was accompanied by reduction in transposon expression in the established callus subcultures. We conclude that stress- and development-related miRNAs are highly expressed upon maize EC callus induction and during maintenance of the subcultures, while miRNAs involved in hormone response only transiently increase during induction. In addition, the establishment of a proliferative status in embryogenic callus is accompanied by important readjustments in hc-siRNAs mapping to long tandem repeat (LTR) retrotransposons, and their expression regulation.

Overlapping RdDM and non-RdDM mechanisms work together to maintain somatic repression of a paramutagenic epiallele of maize pericarp color1

Allelic variation at the Zea mays (maize) pericarp color1 (p1) gene has been attributed to epigenetic gene regulation. A p1 distal enhancer, 5.2 kb upstream of the transcriptional start site, has demonstrated variation in DNA methylation in different p1 alleles/epialleles. In addition, DNA methylation of sequences within the 3’ end of intron 2 also plays a role in tissue-specific expression of p1 alleles. We show here a direct evidence for small RNAs’ involvement in regulating p1 that has not been demonstrated previously. The role of mediator of paramutation1 (mop1) was tested in the maintenance of somatic silencing at distinct p1 alleles: the non-paramutagenic P1-wr allele and paramutagenic P1-rr’ epiallele. The mop1-1 mutation gradually relieves the silenced phenotype after multiple generations of exposure; P1-wr;mop1-1 plants display a loss of 24-nt small RNAs and DNA methylation in the 3’ end of the intron 2, a region close to a Stowaway transposon. In addition, a MULE sequence within the proximal promoter of P1-wr shows depletion of 24nt siRNAs in mop1-1 plants. Release of silencing was not correlated with small RNAs at the distal enhancer region of the P1-wr allele. We found that the somatic silencing of the paramutagenic P1-rr’ is correlated with significantly reduced H3K9me2 in the distal enhancer of P1-rr’; mop1-1 plants, while symmetric DNA methylation is not significantly different. This study highlights that the epigenetic regulation of p1 alleles is controlled both via RdDM as well as non-RdDM mechanisms.

Maize RNA PolIV affects the expression of genes with nearby TE insertions and has a genome-wide repressive impact on transcription

BACKGROUND

RNA-directed DNA methylation (RdDM) is a plant-specific epigenetic process that relies on the RNA polymerase IV (Pol IV) for the production of 24 nucleotide small interfering RNAs (siRNA) that guide the cytosine methylation and silencing of genes and transposons. Zea mays RPD1/RMR6 gene encodes the largest subunit of Pol IV and is required for normal plant development, paramutation, transcriptional repression of certain transposable elements (TEs) and transcriptional regulation of specific alleles.

RESULTS

In this study we applied a total RNA-Seq approach to compare the B73 and rpd1/rmr6 leaf transcriptomes. Although previous studies indicated that loss of siRNAs production in RdDM mutants provokes a strong loss of CHH DNA methylation but not massive gene or TEs transcriptional activation in both Arabidopsis and maize, our total RNA-Seq analysis of rpd1/rmr6 transcriptome reveals that loss of Pol IV activity causes a global increase in the transcribed fraction of the maize genome. Our results point to the genes with nearby TE insertions as being the most strongly affected by Pol IV-mediated gene silencing. TEs modulation of nearby gene expression is linked to alternative methylation profiles on gene flanking regions, and these profiles are strictly dependent on specific characteristics of the TE member inserted. Although Pol IV is essential for the biogenesis of siRNAs, the genes with associated siRNA loci are less affected by the pol IV mutation.

CONCLUSIONS

This deep and integrated analysis of gene expression, TEs distribution, smallRNA targeting and DNA methylation levels, reveals that loss of Pol IV activity globally affects genome regulation, pointing at TEs as modulator of nearby gene expression and indicating the existence of multiple level epigenetic silencing mechanisms. Our results also suggest a predominant role of the Pol IV-mediated RdDM pathway in genome dominance regulation, and subgenome stability and evolution in maize.

Altered nucleosome positions in maize haplotypes and mutants of a subset of SWI/SNF-like proteins

Chromatin remodelers alter DNA-histone interactions in eukaryotic organisms and have been well characterized in yeast and Arabidopsis. While there are maize proteins with similar domains as known remodelers, the ability of the maize proteins to alter nucleosome position has not been reported. Mutant alleles of several maize proteins (RMR1, CHR101, CHR106, CHR127, and CHR156) with similar functional domains to known chromatin remodelers were identified. Altered gene expression of Chr101, Chr106, Chr127, and Chr156 was demonstrated in plants homozygous for the mutant alleles. These mutant genotypes were subjected to nucleosome position analysis to determine whether misregulation of putative maize chromatin proteins would lead to altered DNA-histone interactions. Nucleosome position changes were observed in plants homozygous for chr101, chr106, chr127, and chr156 mutant alleles, suggesting that CHR101, CHR106, CHR127, and CHR156 may affect chromatin structure. The role of RNA polymerases in altering DNA-histone interactions was also tested. Changes in nucleosome position were demonstrated in homozygous mop2-1 individuals. These changes were demonstrated at the b1 tandem repeats and at newly identified loci. Additionally, differential DNA-histone interactions and altered gene expression of putative chromatin remodelers were demonstrated between different maize haplotypes.

Genome-wide Nucleosome Occupancy and Organization Modulates the Plasticity of Gene Transcriptional Status in Maize

Nucleosomes are fundamental units of chromatin that play critical roles in gene regulation by modulating DNA accessibility. However, their roles in regulating tissue-specific gene transcription are poorly understood. Here, we present genome-wide nucleosome maps of maize shoot and endosperm generated by sequencing the micrococcal nuclease digested nucleosomal DNA. The changes of gene transcriptional status between shoot and endosperm were accompanied by preferential nucleosome loss from the promoters and shifts in the first nucleosome downstream of the transcriptional start sites (+1 nucleosome) and upstream of transcriptional termination sites (-1 nucleosome). Intrinsically DNA-encoded nucleosome organization was largely associated with the capacity of a gene to alter its transcriptional status among different tissues. Compared with tissue-specific genes, constitutively expressed genes showed more pronounced 5' and 3' nucleosome-depleted regions as well as further +1 nucleosome to transcriptional start sites and -1 nucleosome to transcriptional termination sites. Moreover, nucleosome organization was more highly correlated with the plasticity of gene transcriptional status than with its expression level when examined using in vivo and predicted nucleosome data. In addition, the translational efficiencies of tissue-specific genes appeared to be greater than those of constitutively expressed genes. Taken together, our results indicate that intrinsically DNA-encoded nucleosome organization is important, beyond its role in regulating gene expression levels, in determining the plasticity of gene transcriptional status.

Arabidopsis proteins with a transposon-related domain act in gene silencing

Transposable elements (TEs) are prevalent in most eukaryotes, and host genomes have devised silencing strategies to rein in TE activity. One of these, transcriptional silencing, is generally associated with DNA methylation and short interfering RNAs. Here we show that the Arabidopsis genes MAIL1 and MAIN define an alternative silencing pathway independent of DNA methylation and short interfering RNAs. Mutants for MAIL1 or MAIN exhibit release of silencing and appear to show impaired condensation of pericentromeric heterochromatin. Phylogenetic analysis suggests not only that MAIL1 and MAIN encode a retrotransposon-related plant mobile domain, but also that host plant mobile domains were captured by DNA transposons during plant evolution. Our results reveal a role for Arabidopsis proteins with a transposon-related domain in gene silencing.

Analysis of small RNAs revealed differential expressions during pollen and embryo sac development in autotetraploid rice

Background

Partial pollen and embryo sac sterilities are the two main reasons for low fertility in autotetraploid rice. Our previous study revealed that small RNAs changes may associate with pollen fertility in autotetraploid rice. However, knowledge on comparative analysis between the development of pollen and embryo sac by small RNAs in autotetraploid rice is still unknown. In the present study, WE-CLSM (whole-mount eosin B-staining confocal laser scanning microscopy) and high-throughput sequencing technology was employed to examine the cytological variations and to analyze small RNAs changes during pollen and embryo sac development in autotetraploid rice compared with its diploid counterpart.

Results

A total of 321 and 368 differentially expressed miRNAs (DEM) were detected during pollen and embryo sac development in autotetraploid rice, respectively. Gene Ontology enrichment analysis on the targets of DEM associated with embryo sac and pollen development revealed 30 prominent functional gene classes, such as cell differentiation and signal transduction during embryo sac development, while only 7 prominent functional gene classes, such as flower development and transcription factor activity, were detected during pollen development in autotetraploid rice. The expression levels of 39 DEM, which revealed interaction with meiosis-related genes, showed opposite expression patterns during pollen and embryo sac development. Of these DEM, osa-miR1436_L + 3_1ss5CT and osa-miR167h-3p were associated with the female meiosis, while osa-miR159a.1 and osa-MIR159a-p5 were related with the male meiosis. 21 nt-phasiRNAs were detected during both pollen and embryo sac development, while 24 nt-phasiRNAs were found only in pollen development, which displayed down-regulation in autotetraploid compared to diploid rice and their spatial-temporal expression patterns were similar to osa-miR2275d. 24 nt TEs-siRNAs were found to be up-regulated in embryo sac but down-regulated in pollen development.

Conclusion

The above results not only provide the small RNAs changes during four landmark stages of pollen and embryo sac development in autotetraploid rice but also have identified specifically expressed miRNAs, especially meiosis-related miRNAs, pollen-specific-24 nt-phasiRNAs and TEs-siRNAs in autotetraploid rice. Together, these findings provide a foundation for understanding the effect of polyploidy on small RNAs expression patterns during pollen and embryo sac development that may lead to different abnormalities in autotetraploid rice.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3526-8) contains supplementary material, which is available to authorized users.

Epigenetic regulation and epigenomic landscape in rice

Epigenetic regulation has been implicated in the control of complex agronomic traits in rice (Oryza sativa), a staple food crop and model monocot plant. Recent advances in high-throughput sequencing and the moderately complex genome of rice have made it possible to study epigenetic regulation in rice on a genome-wide scale. This review discusses recent advances in our understanding of epigenetic regulation in rice, with an emphasis on the roles of key epigenetic regulators, the epigenomic landscape, epigenetic variation, transposon repression, and plant development.

The Tomato Spotted Wilt Virus Genome Is Processed Differentially in its Plant Host Arachis hypogaea and its Thrips Vector Frankliniella fusca

Thrips-transmitted tospoviruses are economically important viruses affecting a wide range of field and horticultural crops worldwide. Tomato spotted wilt virus (TSWV) is the type member of the Tospovirus genus with a broad host range of more than 900 plant species. Interactions between these viruses and their plant hosts and insect vectors via RNAi pathways are likely a key determinant of pathogenicity. The current investigation, for the first time, compares biogenesis of small RNAs between the plant host and insect vector in the presence or absence of TSWV. Unique viral small interfering RNA (vsiRNA) profiles are evident for Arachis hypogaea (peanut) and Frankliniella fusca (thrips vector) following infection with TSWV. Differences between vsiRNA profiles for these plant and insect species, such as the relative abundance of 21 and 22 nt vsiRNAs and locations of alignment hotspots, reflect the diverse siRNA biosynthesis pathways of their respective kingdoms. The presence of unique vsiRNAs in F. fusca samples indicates that vsiRNA generation takes place within the thrips, and not solely through uptake via feeding on vsiRNAs produced in infected A. hypogaea. The study also shows key vsiRNA profile differences for TSWV among plant families, which are evident in the case of A. hypogaea, a legume, and members of Solanaceae (S. lycopersicum and Nicotiana benthamiana). Distinctively, overall small RNA (sRNA) biogenesis in A. hypogaea is markedly affected with an absence of the 24 nt sRNAs in TSWV-infected plants, possibly leading to wide-spread molecular and phenotypic perturbations specific to this species. These findings add significant information on the host-virus-vector interaction in terms of RNAi pathways and may lead to better crop and vector specific control strategies.

Mechanistic and evolutionary questions about epigenetic conflicts between transposable elements and their plant hosts

Transposable elements (TEs) constitute the majority of plant genomes, but most are epigenetically inactivated by their host. Research over the last decade has elucidated many of the molecular components that are required for TE silencing. In contrast, the evolutionary dynamics between TEs and silencing pathways are less clear. Here, we discuss current information about these dynamics from both mechanistic and evolutionary perspectives. We highlight new evidence that palindromic sequences within TEs may act as signals for host recognition and that cis-regulatory regions of TEs may be sites of ongoing arms races with host defenses. We also discuss patterns of TE aging after they are silenced; while there is not yet a consensus, it appears that TEs are removed more rapidly near genes, such that older TE insertions tend to be farther from genes. We conclude by discussing the energetic costs for maintaining silencing pathways, which appear to be substantive. The maintenance of silencing pathways across many species suggests that epigenetic emergencies are frequent.

Genome-Wide Characterization of Maize Small RNA Loci and Their Regulation in the required to maintain repression6-1 (rmr6-1) Mutant and Long-Term Abiotic Stresses

Endogenous small RNAs (sRNAs) contribute to gene regulation and genome homeostasis, but their activities and functions are incompletely known. The maize genome has a high number of transposable elements (TEs; almost 85%), some of which spawn abundant sRNAs. We performed sRNA and total RNA sequencing from control and abiotically stressed B73 wild-type plants and rmr6-1 mutants. RMR6 encodes the largest subunit of the RNA polymerase IV complex and is responsible for accumulation of most 24-nucleotide (nt) small interfering RNAs (siRNAs). We identified novel MIRNA loci and verified miR399 target conservation in maize. RMR6-dependent 23-24 nt siRNA loci were specifically enriched in the upstream region of the most highly expressed genes. Most genes misregulated in rmr6-1 did not show a significant correlation with loss of flanking siRNAs, but we identified one gene supporting existing models of direct gene regulation by TE-derived siRNAs. Long-term drought correlated with changes of miRNA and sRNA accumulation, in particular inducing down-regulation of a set of sRNA loci in the wild-typeleaf.

Modeling the Morphometric Evolution of the Maize Shoot Apical Meristem

The maize (Zea mays subsp. mays L.) shoot apical meristem (SAM) is a self-replenishing pool of stem cells that produces all above-ground plant tissues. Improvements in image acquisition and processing techniques have allowed high-throughput, quantitative genetic analyses of SAM morphology. As with other large-scale phenotyping efforts, meaningful descriptions of genetic architecture depend on the collection of relevant measures. In this study, we tested two quantitative image processing methods to describe SAM morphology within the genus Zea, represented by 33 wild relatives of maize and 841 lines from a domesticated maize by wild teosinte progenitor (MxT) backcross population, along with previously reported data from several hundred diverse maize inbred lines. Approximating the MxT SAM as a paraboloid derived eight parabolic estimators of SAM morphology that identified highly overlapping quantitative trait loci (QTL) on eight chromosomes, which implicated previously identified SAM morphology candidate genes along with new QTL for SAM morphological variation. Using a Fourier-transform related method of comprehensive shape analysis, we detected cryptic SAM shape variation that identified QTL on six chromosomes. We found that Fourier transform shape descriptors and parabolic estimation measures are highly correlated and identified similar QTL. Analysis of shoot apex contours from 73 anciently diverged plant taxa further suggested that parabolic shape may be a universal feature of plant SAMs, regardless of evolutionary clade. Future high-throughput examinations of SAM morphology may benefit from the ease of acquisition and phenotypic fidelity of modeling the SAM as a paraboloid.

RNA-directed DNA methylation enforces boundaries between heterochromatin and euchromatin in the maize genome

The maize genome is relatively large (∼ 2.3 Gb) and has a complex organization of interspersed genes and transposable elements, which necessitates frequent boundaries between different types of chromatin. The examination of maize genes and conserved noncoding sequences revealed that many of these are flanked by regions of elevated asymmetric CHH (where H is A, C, or T) methylation (termed mCHH islands). These mCHH islands are quite short (∼ 100 bp), are enriched near active genes, and often occur at the edge of the transposon that is located nearest to genes. The analysis of DNA methylation in other sequence contexts and several chromatin modifications revealed that mCHH islands mark the transition from heterochromatin-associated modifications to euchromatin-associated modifications. The presence of an mCHH island is fairly consistent in several distinct tissues that were surveyed but shows some variation among different haplotypes. The presence of insertion/deletions in promoters often influences the presence and position of an mCHH island. The mCHH islands are dependent upon RNA-directed DNA methylation activities and are lost in mop1 and mop3 mutants, but the nearby genes rarely exhibit altered expression levels. Instead, loss of an mCHH island is often accompanied by additional loss of DNA methylation in CG and CHG contexts associated with heterochromatin in nearby transposons. This suggests that mCHH islands and RNA-directed DNA methylation near maize genes may act to preserve the silencing of transposons from activity of nearby genes.

Small RNA and degradome sequencing reveals microRNAs and their targets involved in tomato pedicel abscission

We constructed small RNA and degradome sequencing libraries to identify miRNAs and targets involved in tomato pedicel abscission, and confirmed their roles via quantitative real-time PCR. MicroRNAs (miRNAs) are endogenous small RNAs which play crucial negatively regulatory roles at both the transcriptional and post-transcriptional levels in plants; however, limited knowledge is available on the expression profiles of miRNAs and their target genes during tomato pedicel abscission. Taking advantage of small RNA (sRNA) and degradome sequencing technology, a total of 56 known and 11 novel candidate miRNAs targeting 223 mRNA genes were confirmed during pedicel abscission. Gene ontology annotation and KEGG pathway analysis showed that these target genes were significantly enriched in intracellular, membrane-bounded organelle-related biological processes as well as in metabolic, plant-pathogen interaction and hormone signaling pathways. We screened 17 miRNA/target pairs for further analysis and performed quantitative real-time PCR to identify the roles. Cluster analysis of selected miRNAs revealed that the expression profiles of miRNAs varied in different stages of abscission and could be impacted by ethylene treatment. In the present study, the correlations between miRNAs and targets suggested a complex regulatory network of miRNA-mediated target interaction during pedicel abscission. Additionally, the expression profiles of miRNAs and their targets changed by ethylene might be a considerable reason why ethylene promotes pedicel abscission. Our study provides new insights into the expression and regulatory profiles of miRNAs during tomato pedicel abscission.

Paramutation in evolution, population genetics and breeding

Paramutation is a fascinating phenomenon in which directed allelic interactions result in heritable changes in the state of an allele. Paramutation has been carefully characterized at a handful of loci but the prevalence of paramutable/paramutagenic alleles is not well characterized within genomes or populations. In order to consider the role of paramutation in evolutionary processes and plant breeding, we focused on several questions. First, what causes certain alleles to become subject to paramutation? While paramutation clearly involves epigenetic regulation it is also true that only certain alleles defined by genetic sequences are able to participate in paramutation. Second, what is the prevalence of paramutation? There are only a handful of well-documented examples of paramutation. However, there is growing evidence that many loci may undergo changes in chromatin state or expression that are similar to changes observed as a result of paramutation. Third, how will paramutation events be inherited in natural or artificial populations? Many factors, including stability of epigenetic state, mating style and ploidy, may influence the prevalence of paramutation states within populations. Developing a clear understanding of the mechanisms and frequency of paramutation in crop plant genomes will facilitate new opportunities in genetic manipulation, and will also enhance plant breeding programs and our understanding of genome evolution.

Overexpression of Arabidopsis thaliana ERI, the homolog of C. elegans Enhancer of RNAinterference, leads to enhanced growth

Organisms adopt a wide range of strategies to adapt to change. Gene silencing describes the ability of organisms to modulate the expression of susceptible genes at certain times at the transcriptional or the translational level. In all known eukaryotic organisms 21-nt long short interfering RNAs (siRNAs) are the effector molecules of post-transcriptional gene silencing (PTGS), while 24-nt long siRNAs are involved in PTGS in plants. Mutant studies in Caenorhabditis elegans lead to the identification of the enzyme ERI (Enhancer of RNAinterference) with enhanced PTGS. Although the genes involved in growth vigor and growth rate are still unknown, it becomes clearer that the population of small RNAs plays a role in the very early phase of plant development. To pinpoint the link between growth and siRNAs, the expression of Arabidopsis uni-gene Enhancer of RNAi (ERI) homolog from C. elegans was modulated. Increased degradation of small RNAs was achieved by ectopic AtERI overexpression in planta. Based on global small RNA analysis, AtERI overexpression affects mainly the population of 21 mers, excluding miRNAs. To identify target genes, AtERI gain-of-function mutants were analyzed, and differentially abundant small RNAs were identified. Plants with an elevated level of AtERI were bigger in all three light intensities analyzed, indicating an inhibitory function of particular small RNAs in plant growth, with differences in relative growth rates depending on developmental stage and light intensity. Understanding the role of these siRNAs could open new avenues for enhancing plant growth.

Methylation-sensitive expression of a DNA demethylase gene serves as an epigenetic rheostat

Genomes must balance active suppression of transposable elements (TEs) with the need to maintain gene expression. In Arabidopsis, euchromatic TEs are targeted by RNA-directed DNA methylation (RdDM). Conversely, active DNA demethylation prevents accumulation of methylation at genes proximal to these TEs. It is unknown how a cellular balance between methylation and demethylation activities is achieved. Here we show that both RdDM and DNA demethylation are highly active at a TE proximal to the major DNA demethylase gene ROS1. Unexpectedly, and in contrast to most other genomic targets, expression of ROS1 is promoted by DNA methylation and antagonized by DNA demethylation. We demonstrate that inducing methylation in the ROS1 proximal region is sufficient to restore ROS1 expression in an RdDM mutant. Additionally, methylation-sensitive expression of ROS1 is conserved in other species, suggesting it is adaptive. We propose that the ROS1 locus functions as an epigenetic rheostat, tuning the level of demethylase activity in response to methylation alterations, thus ensuring epigenomic stability.

Organisms must adapt to dynamic and variable internal and external environments. Maintaining homeostasis in core biological processes is crucial to minimizing the deleterious consequences of environmental fluctuations. Genomes are also dynamic and variable, and must be robust against stresses, including the invasion of genomic parasites, such as transposable elements (TEs). In this work we present the discovery of an epigenetic rheostat in plants that maintains homeostasis in levels of DNA methylation. DNA methylation typically silences transcription of TEs. Because there is positive feedback between existing and de novo DNA methylation, it is critical that methylation is not allowed to spread and potentially silence transcription of genes. To maintain homeostasis, methylation promotes the production of a demethylase enzyme that removes methylation from gene-proximal regions. The demethylation of genes is therefore always maintained in concert with the levels of methylation suppressing TEs. In addition, this DNA demethylating enzyme also represses its own production in a negative feedback loop. Together, these feedback mechanisms shed new light on how the conflict between gene expression and genome defense is maintained in homeostasis. The presence of this rheostat in multiple species suggests it is an evolutionary conserved adaptation.

Regulation of a maize HD-ZIP IV transcription factor by a non-conventional RDR2-dependent small RNA

Small non-coding RNAs are versatile riboregulators that control gene expression at the transcriptional or post-transcriptional level, governing many facets of plant development. Here we present evidence for the existence of a 24 nt small RNA (named small1) that is complementary to the 3' UTR of OCL1 (Outer Cell Layer1), the founding member of the maize HD-ZIP IV gene family encoding plant-specific transcription factors that are mainly involved in epidermis differentiation and specialization. The biogenesis of small1 depends on DICER-like 3 (DCL3), RNA-dependent RNA polymerase 2 (RDR2) and RNA polymerase IV, components that are usually required for RNA-dependent DNA-methylation. Unexpectedly, GFP sensor experiments in transient and stable transformation systems revealed that small1 may regulate its target at the post-transcriptional level, mainly through translational repression. This translational repression is attenuated in an rdr2 mutant background in which small1 does not accumulate. Our experiments further showed the possible involvement of a secondary stem-loop structure present in the 3' UTR of OCL1 for efficient target repression, suggesting the existence of several regulatory mechanisms affecting OCL1 mRNA stability and translation.

Spatiotemporally dynamic, cell-type-dependent premeiotic and meiotic phasiRNAs in maize anthers

Maize anthers, the male reproductive floral organs, express two classes of phased small-interfering RNAs (phasiRNAs). PhasiRNA precursors are transcribed by RNA polymerase II and map to low-copy, intergenic regions similar to PIWI-interacting RNAs (piRNAs) in mammalian testis. From 10 sequential cohorts of staged maize anthers plus mature pollen we find that 21-nt phased siRNAs from 463 loci appear abruptly after germinal and initial somatic cell fate specification and then diminish, whereas 24-nt phasiRNAs from 176 loci coordinately accumulate during meiosis and persist as anther somatic cells mature and haploid gametophytes differentiate into pollen. Male-sterile ocl4 anthers defective in epidermal signaling lack 21-nt phasiRNAs. Male-sterile mutants with subepidermal defects--mac1 (excess meiocytes), ms23 (defective pretapetal cells), and msca1 (no normal soma or meiocytes)--lack 24-nt phasiRNAs. ameiotic1 mutants (normal soma, no meiosis) accumulate both 21-nt and 24-nt phasiRNAs, ruling out meiotic cells as a source or regulator of phasiRNA biogenesis. By in situ hybridization, miR2118 triggers of 21-nt phasiRNA biogenesis localize to epidermis; however, 21-PHAS precursors and 21-nt phasiRNAs are abundant subepidermally. The miR2275 trigger, 24-PHAS precursors, and 24-nt phasiRNAs all accumulate preferentially in tapetum and meiocytes. Therefore, each phasiRNA type exhibits independent spatiotemporal regulation with 21-nt premeiotic phasiRNAs dependent on epidermal and 24-nt meiotic phasiRNAs dependent on tapetal cell differentiation. Maize phasiRNAs and mammalian piRNAs illustrate putative convergent evolution of small RNAs in male reproduction.

Nascent transcription affected by RNA polymerase IV in Zea mays

All eukaryotes use three DNA-dependent RNA polymerases (RNAPs) to create cellular RNAs from DNA templates. Plants have additional RNAPs related to Pol II, but their evolutionary role(s) remain largely unknown. Zea mays (maize) RNA polymerase D1 (RPD1), the largest subunit of RNA polymerase IV (Pol IV), is required for normal plant development, paramutation, transcriptional repression of certain transposable elements (TEs), and transcriptional regulation of specific alleles. Here, we define the nascent transcriptomes of rpd1 mutant and wild-type (WT) seedlings using global run-on sequencing (GRO-seq) to identify the broader targets of RPD1-based regulation. Comparisons of WT and rpd1 mutant GRO-seq profiles indicate that Pol IV globally affects transcription at both transcriptional start sites and immediately downstream of polyadenylation addition sites. We found no evidence of divergent transcription from gene promoters as seen in mammalian GRO-seq profiles. Statistical comparisons identify genes and TEs whose transcription is affected by RPD1. Most examples of significant increases in genic antisense transcription appear to be initiated by 3'-proximal long terminal repeat retrotransposons. These results indicate that maize Pol IV specifies Pol II-based transcriptional regulation for specific regions of the maize genome including genes having developmental significance.

Genome evolution in maize: from genomes back to genes

Maize occupies dual roles as both (a) one of the big-three grain species (along with rice and wheat) responsible for providing more than half of the calories consumed around the world, and (b) a model system for plant genetics and cytogenetics dating back to the origin of the field of genetics in the early twentieth century. The long history of genetic investigation in this species combined with modern genomic and quantitative genetic data has provided particular insight into the characteristics of genes linked to phenotypes and how these genes differ from many other sequences in plant genomes that are not easily distinguishable based on molecular data alone. These recent results suggest that the number of genes in plants that make significant contributions to phenotype may be lower than the number of genes defined by current molecular criteria, and also indicate that syntenic conservation has been underemphasized as a marker for gene function.

Detection of Pol IV/RDR2-dependent transcripts at the genomic scale in Arabidopsis reveals features and regulation of siRNA biogenesis

Twenty-four-nucleotide small interfering (si)RNAs are central players in RNA-directed DNA methylation (RdDM), a process that establishes and maintains DNA methylation at transposable elements to ensure genome stability in plants. The plant-specific RNA polymerase IV (Pol IV) is required for siRNA biogenesis and is believed to transcribe RdDM loci to produce primary transcripts that are converted to double-stranded RNAs (dsRNAs) by RDR2 to serve as siRNA precursors. Yet, no such siRNA precursor transcripts have ever been reported. Here, through genome-wide profiling of RNAs in genotypes that compromise the processing of siRNA precursors, we were able to identify Pol IV/RDR2-dependent transcripts from tens of thousands of loci. We show that Pol IV/RDR2-dependent transcripts correspond to both DNA strands, whereas the RNA polymerase II (Pol II)-dependent transcripts produced upon derepression of the loci are derived primarily from one strand. We also show that Pol IV/RDR2-dependent transcripts have a 5′ monophosphate, lack a poly(A) tail at the 3′ end, and contain no introns; these features distinguish them from Pol II-dependent transcripts. Like Pol II-transcribed genic regions, Pol IV-transcribed regions are flanked by A/T-rich sequences depleted in nucleosomes, which highlights similarities in Pol II- and Pol IV-mediated transcription. Computational analysis of siRNA abundance from various mutants reveals differences in the regulation of siRNA biogenesis at two types of loci that undergo CHH methylation via two different DNA methyltransferases. These findings begin to reveal features of Pol IV/RDR2-mediated transcription at the heart of genome stability in plants.

Differential nuclease sensitivity profiling of chromatin reveals biochemical footprints coupled to gene expression and functional DNA elements in maize

The eukaryotic genome is organized into nucleosomes, the fundamental units of chromatin. The positions of nucleosomes on DNA regulate protein-DNA interactions and in turn influence DNA-templated events. Despite the increasing number of genome-wide maps of nucleosome position, how global changes in gene expression relate to changes in nucleosome position is poorly understood. We show that in nucleosome occupancy mapping experiments in maize (Zea mays), particular genomic regions are highly susceptible to variation introduced by differences in the extent to which chromatin is digested with micrococcal nuclease (MNase). We exploited this digestion-linked variation to identify protein footprints that are hypersensitive to MNase digestion, an approach we term differential nuclease sensitivity profiling (DNS-chip). Hypersensitive footprints were enriched at the 5' and 3' ends of genes, associated with gene expression levels, and significantly overlapped with conserved noncoding sequences and the binding sites of the transcription factor KNOTTED1. We also found that the tissue-specific regulation of gene expression was linked to tissue-specific hypersensitive footprints. These results reveal biochemical features of nucleosome organization that correlate with gene expression levels and colocalize with functional DNA elements. This approach to chromatin profiling should be broadly applicable to other species and should shed light on the relationships among chromatin organization, protein-DNA interactions, and genome regulation.

Dynamic transcriptome landscape of maize embryo and endosperm development

Maize (Zea mays) is an excellent cereal model for research on seed development because of its relatively large size for both embryo and endosperm. Despite the importance of seed in agriculture, the genome-wide transcriptome pattern throughout seed development has not been well characterized. Using high-throughput RNA sequencing, we developed a spatiotemporal transcriptome atlas of B73 maize seed development based on 53 samples from fertilization to maturity for embryo, endosperm, and whole seed tissues. A total of 26,105 genes were found to be involved in programming seed development, including 1,614 transcription factors. Global comparisons of gene expression highlighted the fundamental transcriptomic reprogramming and the phases of development. Coexpression analysis provided further insight into the dynamic reprogramming of the transcriptome by revealing functional transitions during maturation. Combined with the published nonseed high-throughput RNA sequencing data, we identified 91 transcription factors and 1,167 other seed-specific genes, which should help elucidate key mechanisms and regulatory networks that underlie seed development. In addition, correlation of gene expression with the pattern of DNA methylation revealed that hypomethylation of the gene body region should be an important factor for the expressional activation of seed-specific genes, especially for extremely highly expressed genes such as zeins. This study provides a valuable resource for understanding the genetic control of seed development of monocotyledon plants.

Diverse gene-silencing mechanisms with distinct requirements for RNA polymerase subunits in Zea mays

In Zea mays, transcriptional regulation of the b1 (booster1) gene requires a distal enhancer and MEDIATOR OF PARAMUTATION1 (MOP1), MOP2, and MOP3 proteins orthologous to Arabidopsis components of the RNA-dependent DNA methylation pathway. We compared the genetic requirements for MOP1, MOP2, and MOP3 for endogenous gene silencing by two hairpin transgenes with inverted repeats of the a1 (anthocyaninless1) gene promoter (a1pIR) and the b1 gene enhancer (b1IR), respectively. The a1pIR transgene induced silencing of endogenous A1 in mop1-1 and mop3-1, but not in Mop2-1 homozygous plants. This finding suggests that transgene-derived small interfering RNAs (siRNAs) circumvented the requirement for MOP1, a predicted RNA-dependent RNA polymerase, and MOP3, the predicted largest subunit of RNA polymerase IV (Pol IV). Because the Arabidopsis protein orthologous to MOP2 is the second largest subunit of Pol IV and V, our results may indicate that hairpin-induced siRNAs cannot bypass the requirement for the predicted scaffolding activity of Pol V. In contrast to a1pIR, the b1IR transgene silenced endogenous B1 in all three homozygous mutant genotypes--mop1-1, Mop2-1, and mop3-1--suggesting that transgene mediated b1 silencing did not involve MOP2-containing Pol V complexes. Based on the combined results for a1, b1, and three previously described loci, we propose a speculative hypothesis of locus-specific deployment of Pol II, MOP2-containing Pol V, or alternative versions of Pol V with second largest subunits other than MOP2 to explain the mechanistic differences in silencing at specific loci, including one example associated with paramutation.

Arabidopsis and maize RidA proteins preempt reactive enamine/imine damage to branched-chain amino acid biosynthesis in plastids

RidA (for Reactive Intermediate Deaminase A) proteins are ubiquitous, yet their function in eukaryotes is unclear. It is known that deleting Salmonella enterica ridA causes Ser sensitivity and that S. enterica RidA and its homologs from other organisms hydrolyze the enamine/imine intermediates that Thr dehydratase forms from Ser or Thr. In S. enterica, the Ser-derived enamine/imine inactivates a branched-chain aminotransferase; RidA prevents this damage. Arabidopsis thaliana and maize (Zea mays) have a RidA homolog that is predicted to be plastidial. Expression of either homolog complemented the Ser sensitivity of the S. enterica ridA mutant. The purified proteins hydrolyzed the enamines/imines formed by Thr dehydratase from Ser or Thr and protected the Arabidopsis plastidial branched-chain aminotransferase BCAT3 from inactivation by the Ser-derived enamine/imine. In vitro chloroplast import assays and in vivo localization of green fluorescent protein fusions showed that Arabidopsis RidA and Thr dehydratase are chloroplast targeted. Disrupting Arabidopsis RidA reduced root growth and raised the root and shoot levels of the branched-chain amino acid biosynthesis intermediate 2-oxobutanoate; Ser treatment exacerbated these effects in roots. Supplying Ile reversed the root growth defect. These results indicate that plastidial RidA proteins can preempt damage to BCAT3 and Ile biosynthesis by hydrolyzing the Ser-derived enamine/imine product of Thr dehydratase.