Epigenomic mapping in Arabidopsis using tiling microarrays

Quantitative trait loci from identification to exploitation for crop improvement

Advancement in the field of genetics and genomics after the discovery of Mendel's laws of inheritance has led to map the genes controlling qualitative and quantitative traits in crop plant species. Mapping of genomic regions controlling the variation of quantitatively inherited traits has become routine after the advent of different types of molecular markers. Recently, the next generation sequencing methods have accelerated the research on QTL analysis. These efforts have led to the identification of more closely linked molecular markers with gene/QTLs and also identified markers even within gene/QTL controlling the trait of interest. Efforts have also been made towards cloning gene/QTLs or identification of potential candidate genes responsible for a trait. Further new concepts like crop QTLome and QTL prioritization have accelerated precise application of QTLs for genetic improvement of complex traits. In the past years, efforts have also been made in exploitation of a number of QTL for improving grain yield or other agronomic traits in various crops through markers assisted selection leading to cultivation of these improved varieties at farmers' field. In present article, we reviewed QTLs from their identification to exploitation in plant breeding programs and also reviewed that how improved cultivars developed through introgression of QTLs have improved the yield productivity in many crops.

Analysis of Small RNA Populations Using Hybridization to DNA Tiling Arrays

Epigenetic response to stress in plants involves changes in DNA methylation, histone modifications, and expression of small noncoding RNAs (sRNA). Here we present the method of analysis of differential expression of sRNA populations using DNA tiling arrays. sRNA extracted from Arabidopsis thaliana plants exposed to pathogen elicitor or control plants were reverse-transcribed into cDNAs, and subsequently hybridized after labeling to a custom-made DNA tiling array covering Arabidopsis chromosome 4. We first designed a control experiment with eight cDNA clones corresponding to sequences located on chromosome 4 and obtained robust and specific hybridization signals. Furthermore, hybridization signals along chromosome 4 were in good agreement with sRNA abundance as previously determined by massive parallel sequence signature (MPSS) in the case of untreated plants, but differed substantially after stress treatment. These results demonstrate the utility of hybridization to DNA tiling arrays to detect major changes in sRNA abundance.

Accessing epigenetic variation in the plant methylome

Cytosine DNA methylation is the addition of a methyl group to the 5' position of a cytosine, which plays a crucial role in plant development and gene silencing. Genome-wide profiling of DNA methylation is now possible using various techniques and strategies. Using these technologies, we are beginning to elucidate the extent and impact of variation in DNA methylation between individuals and/or tissues. Here, we review the different techniques used to analyze the methylomes at the whole-genome level and their applications to better understand epigenetic variations in plants.

Insights into chromatin structure and dynamics in plants

The packaging of chromatin into the nucleus of a eukaryotic cell requires an extraordinary degree of compaction and physical organization. In recent years, it has been shown that this organization is dynamically orchestrated to regulate responses to exogenous stimuli as well as to guide complex cell-type-specific developmental programs. Gene expression is regulated by the compartmentalization of functional domains within the nucleus, by distinct nucleosome compositions accomplished via differential modifications on the histone tails and through the replacement of core histones by histone variants. In this review, we focus on these aspects of chromatin organization and discuss novel approaches such as live cell imaging and photobleaching as important tools likely to give significant insights into our understanding of the very dynamic nature of chromatin and chromatin regulatory processes. We highlight the contribution plant studies have made in this area showing the potential advantages of plants as models in understanding this fundamental aspect of biology.

Genome-wide survey of cold stress regulated alternative splicing in Arabidopsis thaliana with tiling microarray

Alternative splicing plays a major role in expanding the potential informational content of eukaryotic genomes. It is an important post-transcriptional regulatory mechanism that can increase protein diversity and affect mRNA stability. Alternative splicing is often regulated in a tissue-specific and stress-responsive manner. Cold stress, which adversely affects plant growth and development, regulates the transcription and splicing of plant splicing factors. This can affect the pre-mRNA processing of many genes. To identify cold regulated alternative splicing we applied Affymetrix Arabidopsis tiling arrays to survey the transcriptome under cold treatment conditions. A novel algorithm was used for detection of statistically relevant changes in intron expression within a transcript between control and cold growth conditions. A reverse transcription polymerase chain reaction (RT-PCR) analysis of a number of randomly selected genes confirmed the changes in splicing patterns under cold stress predicted by tiling array. Our analysis revealed new types of cold responsive genes. While their expression level remains relatively unchanged under cold stress their splicing pattern shows detectable changes in the relative abundance of isoforms. The majority of cold regulated alternative splicing introduced a premature termination codon (PTC) into the transcripts creating potential targets for degradation by the nonsense mediated mRNA decay (NMD) process. A number of these genes were analyzed in NMD-defective mutants by RT-PCR and shown to evade NMD. This may result in new and truncated proteins with altered functions or dominant negative effects. The results indicate that cold affects both quantitative and qualitative aspects of gene expression.

Transcript profiling in Arabidopsis with genome tiling microarrays

Microarray technology is at present a standardized workflow for genome-wide expression analysis. Whole-genome tiling microarrays have emerged as an important platform for flexible and comprehensive expression profiling. In this chapter we describe a detailed standardized workflow for experiments assessing the transcriptome of Arabidopsis using tiling arrays and provide useful hints for critical steps from experimental design to data analysis. Although the protocol is optimized for AGRONOMICS1 arrays, it can readily be adapted to other tiling arrays. AGRONOMICS1 is the first platform that enables strand-specific expression analysis of the Arabidopsis genome with a single array. Moreover, it includes all perfect match probes from the original ATH1 array, allowing readily integration with the large existing ATH1 knowledge base. This workflow is designed for the analysis of raw data for any number of samples and it does not pose any particular hardware requirements.

Parsimonious higher-order hidden Markov models for improved array-CGH analysis with applications to Arabidopsis thaliana

Array-based comparative genomic hybridization (Array-CGH) is an important technology in molecular biology for the detection of DNA copy number polymorphisms between closely related genomes. Hidden Markov Models (HMMs) are popular tools for the analysis of Array-CGH data, but current methods are only based on first-order HMMs having constrained abilities to model spatial dependencies between measurements of closely adjacent chromosomal regions. Here, we develop parsimonious higher-order HMMs enabling the interpolation between a mixture model ignoring spatial dependencies and a higher-order HMM exhaustively modeling spatial dependencies. We apply parsimonious higher-order HMMs to the analysis of Array-CGH data of the accessions C24 and Col-0 of the model plant Arabidopsis thaliana. We compare these models against first-order HMMs and other existing methods using a reference of known deletions and sequence deviations. We find that parsimonious higher-order HMMs clearly improve the identification of these polymorphisms. Moreover, we perform a functional analysis of identified polymorphisms revealing novel details of genomic differences between C24 and Col-0. Additional model evaluations are done on widely considered Array-CGH data of human cell lines indicating that parsimonious HMMs are also well-suited for the analysis of non-plant specific data. All these results indicate that parsimonious higher-order HMMs are useful for Array-CGH analyses. An implementation of parsimonious higher-order HMMs is available as part of the open source Java library Jstacs (www.jstacs.de/index.php/PHHMM).

Array-based comparative genomics is a standard approach for the identification of DNA copy number polymorphisms between closely related genomes. The huge amounts of data produced by these experiments require efficient and accurate bioinformatics tools for the identification of copy number polymorphisms. Hidden Markov Models (HMMs) are frequently used for analyzing such data sets, but current models are based on first-order HMMs only having limited capabilities to model spatial dependencies between measurements of closely adjacent chromosomal regions. We develop parsimonious higher-order HMMs enabling the interpolation between a mixture model ignoring spatial dependencies and a higher-order HMM exhaustively modeling these dependencies to overcome this limitation. In an in-depth case study with Arabidopsis thaliana, we find that parsimonious higher-order HMMs clearly improve the identification of copy number polymorphisms in comparison to standard first-order HMMs and other frequently used methods. Functional analysis of identified polymorphisms revealed details of genomic differences between the accessions C24 and Col-0 of Arabidopsis thaliana. An additional study on human cell lines further indicates that parsimonious HMMs are well-suited for the analysis of Array-CGH data.

Additive inheritance of histone modifications in Arabidopsis thaliana intra-specific hybrids

Plant genomes are earmarked with defined patterns of chromatin marks. Little is known about the stability of these epigenomes when related, but distinct genomes are brought together by intra-species hybridization. Arabidopsis thaliana accessions and their reciprocal hybrids were used as a model system to investigate the dynamics of histone modification patterns. The genome-wide distribution of histone modifications H3K4me2 and H3K27me3 in the inbred parental accessions Col-0, C24 and Cvi and their hybrid offspring was compared by chromatin immunoprecipitation in combination with genome tiling array hybridization. The analysis revealed that, in addition to DNA sequence polymorphisms, chromatin modification variations exist among accessions of A. thaliana. The range of these variations was higher for H3K27me3 (typically a repressive mark) than for H3K4me2 (typically an active mark). H3K4me2 and H3K27me3 were rather stable in response to intra-species hybridization, with mainly additive inheritance in hybrid offspring. In conclusion, intra-species hybridization does not result in gross changes to chromatin modifications.

The epigenome and plant development

The epigenomic regulation of chromatin structure and genome stability is essential for the interpretation of genetic information and ultimately the determination of phenotype. High-resolution maps of plant epigenomes have been obtained through a combination of chromatin technologies and genomic tiling microarrays and through high-throughput sequencing-based approaches. The transcriptomic activity of a plant at a certain stage of development is controlled by genome-wide combinatorial interactions of epigenetic modifications. Tissue- or environment-specific epigenomes are established during plant development. Epigenomic reprogramming triggered by the activation and movement of small RNAs is important for plant gametogenesis. Genome-wide loss of DNA methylation in the endosperm and the accompanying endosperm-specific gene expression during seed development provide a genomic insight into epigenetic regulation of gene imprinting in plants. Global changes of histone modifications during plant responses to different light environments play an important regulatory role in a sophisticated light-regulated transcriptional network. Epigenomic natural variation that developed during evolution is important for phenotypic diversity and can potentially contribute to the molecular mechanisms of complex biological phenomena such as heterosis in plants.

Chlorine ions but not sodium ions alter genome stability of Arabidopsis thaliana

Various environmental stresses influence plant genome stability. Most of these stresses, such as ionizing radiation, heavy metals and organic chemicals, represent potent DNA-damaging agents. Here, we show that exposure to NaCl, the stress that is not thought to cause direct DNA damage, results in an increase in the level of strand breaks and homologous recombination rates (RRs) in Arabidopsis thaliana plants. The effect of salt stress on the RR was found to be primarily associated with Cl(-) ions, since exposure of plants to Na(2)SO(4) did not increase the RR, whereas exposure to MgCl(2) resulted in an increase. Changes in the number of strand breaks and in the RR were also paralleled by transcriptional activation of AtRad51 and down-regulation of AtKu70. The progeny of exposed plants exhibited higher RRs, higher expression of AtRad51, lower expression of AtKu70, higher tolerance to salt and methyl methane sulfate (MMS) stresses, as well as a higher increase in RR upon further exposure to stress. Our experiments showed that NaCl is a genotoxic stress that leads to somatic and transgenerational changes in recombination rates, and these changes are primarily triggered by exposure to Cl(-) ions.

Arabidopsis tiling array analysis to identify the stress-responsive genes

Plants respond and adapt to drought, cold, and high-salinity stresses. Stress-inducible gene products function in the stress response and tolerance in plants. Using cDNA microarrays and oligonucleotide microarrays, stress-inducible genes have been identified in various plant species so far. Recently, tiling array technology has become a powerful tool for the whole-genome transcriptome analysis. We applied the Arabidopsis Affymetrix tiling arrays to study the whole-genome transcriptome under drought, cold, and high-salinity stresses and identified a large number of drought, cold, and high-salinity stress-inducible genes and transcriptional units (TUs).

Analysis of small RNA populations using hybridization to DNA tiling arrays

Small RNA (sRNA) populations extracted from Arabidopsis plants submitted or not to biotic stress, were reverse-transcribed into cDNAs, and these were subsequently hybridized after labelling to a custom-made DNA tiling array covering Arabidopsis chromosome 4. We first designed a control experiment with eight cDNA clones corresponding to sequences located on chromosome 4 and obtained robust and specific hybridization signals. Furthermore, hybridization signals along chromosome 4 were in good agreement with sRNA abundance as previously determined by Massive Parallel Sequence Signature (MPSS) in the case of untreated plants, but differed substantially after stress treatment. These results demonstrate the utility of hybridization to DNA tiling arrays to detect major changes in small RNA populations.

A transposon-induced epigenetic change leads to sex determination in melon

Sex determination in plants leads to the development of unisexual flowers from an originally bisexual floral meristem. This mechanism results in the enhancement of outcrossing and promotes genetic variability, the consequences of which are advantageous to the evolution of a species. In melon, sexual forms are controlled by identity of the alleles at the andromonoecious (a) and gynoecious (g) loci. We previously showed that the a gene encodes an ethylene biosynthesis enzyme, CmACS-7, that represses stamen development in female flowers. Here we show that the transition from male to female flowers in gynoecious lines results from epigenetic changes in the promoter of a transcription factor, CmWIP1. This natural and heritable epigenetic change resulted from the insertion of a transposon, which is required for initiation and maintenance of the spreading of DNA methylation to the CmWIP1 promoter. Expression of CmWIP1 leads to carpel abortion, resulting in the development of unisexual male flowers. Moreover, we show that CmWIP1 indirectly represses the expression of the andromonoecious gene, CmACS-7, to allow stamen development. Together our data indicate a model in which CmACS-7 and CmWIP1 interact to control the development of male, female and hermaphrodite flowers in melon.

Next-generation sequencing reveals complex relationships between the epigenome and transcriptome in maize

Epigenetic modifications and small RNAs play an important role in gene regulation. Here, we discuss results of our Solexa/Illumina 1G sequencing-based survey of DNA methylation, activating and repressive histone modifications, small RNAs and mRNA in the maize genome. We analyze tissue-specific epigenetic patterns, discuss antagonistic relationships between repressive epigenetic marks and highlight synergistic relationships between activating histone modifications. We discuss our observation that small RNAs show a tissue-specific distribution in maize. Whereas 24-nucleotide long small interfering RNAs (siRNAs) accumulated preferentially in shoots, 21-nucleotide long micro RNAs (miRNAs) were the most abundant group in roots, which follows the transcript level of mop1. Furthermore, we discuss the possibility that a novel class of 22-nucleotide siRNAs might originate from long double-stranded RNAs in an RNA-dependent RNA polymerase (RdRP)-independent manner. This supports the intriguing possibility that maize possesses at least two distinct pathways to generate siRNAs, one of which relies on RdRP and a second one that might be RdRP-independent.

Genome-wide and organ-specific landscapes of epigenetic modifications and their relationships to mRNA and small RNA transcriptomes in maize

Maize (Zea mays) has an exceptionally complex genome with a rich history in both epigenetics and evolution. We report genomic landscapes of representative epigenetic modifications and their relationships to mRNA and small RNA (smRNA) transcriptomes in maize shoots and roots. The epigenetic patterns differed dramatically between genes and transposable elements, and two repressive marks (H3K27me3 and DNA methylation) were usually mutually exclusive. We found an organ-specific distribution of canonical microRNAs (miRNAs) and endogenous small interfering RNAs (siRNAs), indicative of their tissue-specific biogenesis. Furthermore, we observed that a decreasing level of mop1 led to a concomitant decrease of 24-nucleotide siRNAs relative to 21-nucleotide miRNAs in a tissue-specific manner. A group of 22-nucleotide siRNAs may originate from long-hairpin double-stranded RNAs and preferentially target gene-coding regions. Additionally, a class of miRNA-like smRNAs, whose putative precursors can form short hairpins, potentially targets genes in trans. In summary, our data provide a critical analysis of the maize epigenome and its relationships to mRNA and smRNA transcriptomes.

Defining gene and QTL networks

Current technologies for high-throughput molecular profiling of large numbers of genetically different individuals offer great potential for elucidating the genotype-to-phenotype relationship. Variation in molecular and phenotypic traits can be correlated to DNA sequence variation using the methods of quantitative trait locus (QTL) mapping. In addition, the correlation structure in the molecular and phenotypic traits can be informative for inferring the underlying molecular networks. For this, new methods are emerging to distinguish among causality, reactivity, or independence of traits based upon logic involving underlying QTL. These methods are becoming increasingly popular in plant genetic studies as well as in studies on many other organisms.

Finding the fifth base: genome-wide sequencing of cytosine methylation

Complete sequences of myriad eukaryotic genomes, including several human genomes, are now available, and recent dramatic developments in DNA sequencing technology are opening the floodgates to vast volumes of sequence data. Yet, despite knowing for several decades that a significant proportion of cytosines in the genomes of plants and animals are present in the form of methylcytosine, until very recently the precise locations of these modified bases have never been accurately mapped throughout a eukaryotic genome. Advanced "next-generation" DNA sequencing technologies are now enabling the global mapping of this epigenetic modification at single-base resolution, providing new insights into the regulation and dynamics of DNA methylation in genomes.

Epigenomic consequences of immortalized plant cell suspension culture

Plant cells grown in culture exhibit genetic and epigenetic instability. Using a combination of chromatin immunoprecipitation and DNA methylation profiling on tiling microarrays, we have mapped the location and abundance of histone and DNA modifications in a continuously proliferating, dedifferentiated cell suspension culture of Arabidopsis. We have found that euchromatin becomes hypermethylated in culture and that a small percentage of the hypermethylated genes become associated with heterochromatic marks. In contrast, the heterochromatin undergoes dramatic and very precise DNA hypomethylation with transcriptional activation of specific transposable elements (TEs) in culture. High throughput sequencing of small interfering RNA (siRNA) revealed that TEs activated in culture have increased levels of 21-nucleotide (nt) siRNA, sometimes at the expense of the 24-nt siRNA class. In contrast, TEs that remain silent, which match the predominant 24-nt siRNA class, do not change significantly in their siRNA profiles. These results implicate RNA interference and chromatin modification in epigenetic restructuring of the genome following the activation of TEs in immortalized cell culture.

Whole-genome microarrays: applications and technical issues

DNA microarrays have become a mainstream tool in experimental plant biology. The constant improvements in the technological platforms have enabled the development of the tiling DNA microarrays that cover the whole genome, which in turn catalyzed the wide variety of creative applications of such microarrays in the areas as diverse as global studies of genetic variation, DNA-binding proteins, DNA methylation, and chromatin and transcriptome dynamics. This chapter attempts to summarize such applications as well as discusses some technical and strategic issues that are particular to the use of tiling microarrays.

Variation for homoeologous gene silencing in hexaploid wheat

The absence of expression of individual members of a homoeologous set of genes in a polyploid is a well-established phenomenon. However, the extent to which such 'homoeologous silencing' can vary between individual genotypes within a species is unexplored. We have used the single-strand conformation polymorphism assay to identify homoeologue non-expression at 15 single-copy genes across a panel of 16 wheat varieties, representative of the genetic diversity present in modern northern European winter wheat (Triticum aestivum). There was no evidence for any homoeologous silencing at seven of the fifteen genes, but in the remaining eight, at least one of the three homoeologues varied qualitatively for expression in either the root or the seedling leaf. The identity of the non-expressed homoeologue was generally consistent, but when the expression profiles of eight informative genes were compared, only two varieties shared the same pattern of silencing. A small-scale study suggested that silencing patterns were largely inherited across self-pollinated generations, and some evidence is presented for the epigenetic segregation of these patterns in a population bred from parents having contrasting silencing profiles. Epigenetic variation exerts a significant effect on phenotype, so given the ubiquity and variability in homoeologous silencing observed in wheat, we suggest that it is likely to play a considerable role in generating phenotypic variation. Thus epigenetic profiling may need to be incorporated as part of the analytical tool kit for predictive wheat breeding.

Arabidopsis transcriptome analysis under drought, cold, high-salinity and ABA treatment conditions using a tiling array

Plants respond and adapt to drought, cold and high-salinity stresses in order to survive. In this study, we applied Arabidopsis Affymetrix tiling arrays to study the whole genome transcriptome under drought, cold, high-salinity and ABA treatment conditions. The bioinformatic analysis using the tiling array data showed that 7,719 non-AGI transcriptional units (TUs) exist in the unannotated "intergenic" regions of Arabidopsis genome. These include 1,275 and 181 TUs that are induced and downregulated, respectively, by the stress or ABA treatments. Most of the non-AGI TUs are hypothetical non-protein-coding RNAs. About 80% of the non-AGI TUs belong to pairs of the fully overlapping sense-antisense transcripts (fSATs). Significant linear correlation between the expression ratios (treated/untreated) of the sense TUs and the ratios of the antisense TUs was observed in the SATs of AGI code/non-AGI TU. We studied the biogenesis mechanisms of the stress- or ABA-inducible antisense RNAs and found that the expression of sense TUs is necessary for the stress- or ABA-inducible expression of the antisense TUs in the fSATs (AGI code/non-AGI TU).

Comprehensive high-throughput arrays for relative methylation (CHARM)

This study was originally conceived to test in a rigorous way the specificity of three major approaches to high-throughput array-based DNA methylation analysis: (1) MeDIP, or methylated DNA immunoprecipitation, an example of antibody-mediated methyl-specific fractionation; (2) HELP, or HpaII tiny fragment enrichment by ligation-mediated PCR, an example of differential amplification of methylated DNA; and (3) fractionation by McrBC, an enzyme that cuts most methylated DNA. These results were validated using 1466 Illumina methylation probes on the GoldenGate methylation assay and further resolved discrepancies among the methods through quantitative methylation pyrosequencing analysis. While all three methods provide useful information, there were significant limitations to each, specifically bias toward CpG islands in MeDIP, relatively incomplete coverage in HELP, and location imprecision in McrBC. However, we found that with an original array design strategy using tiling arrays and statistical procedures that average information from neighboring genomic locations, much improved specificity and sensitivity could be achieved, e.g., approximately 100% sensitivity at 90% specificity with McrBC. We term this approach "comprehensive high-throughput arrays for relative methylation" (CHARM). While this approach was applied to McrBC analysis, the array design and computational algorithms are fractionation method-independent and make this a simple, general, relatively inexpensive tool suitable for genome-wide analysis, and in which individual samples can be assayed reliably at very high density, allowing locus-level genome-wide epigenetic discrimination of individuals, not just groups of samples. Furthermore, unlike the other approaches, CHARM is highly quantitative, a substantial advantage in application to the study of human disease.

Utilizing tiling microarrays for whole-genome analysis in plants

The recent explosion in available genome sequence data has ushered in an era in which analysis of a whole genome can be performed in a single experiment. While DNA microarrays have long been the established technology for measuring gene expression levels, standard expression arrays use relatively few probes for each gene and are typically biased toward known and predicted gene structures. Recently, with the availability of complete genome sequences for many organisms, very-high-density oligonucleotide-based microarrays that span the entire genome have emerged as the preferred platform for genomic analysis. Whole-genome tiling microarrays can be employed for a myriad of purposes, including empirical annotation of the transcriptome, chromatin immunoprecipitation-chip studies, analysis of alternative RNA splicing, characterization of the methylation state of cytosine bases throughout a genome (methylome), and DNA polymorphism discovery. Here, we review several applications of whole-genome technology to obtain a variety of genomic-scale information in plants.

Epigenetic natural variation in Arabidopsis thaliana

Cytosine methylation of repetitive sequences is widespread in plant genomes, occurring in both symmetric (CpG and CpNpG) as well as asymmetric sequence contexts. We used the methylation-dependent restriction enzyme McrBC to profile methylated DNA using tiling microarrays of Arabidopsis Chromosome 4 in two distinct ecotypes, Columbia and Landsberg erecta. We also used comparative genome hybridization to profile copy number polymorphisms. Repeated sequences and transposable elements (TEs), especially long terminal repeat retrotransposons, are densely methylated, but one third of genes also have low but detectable methylation in their transcribed regions. While TEs are almost always methylated, genic methylation is highly polymorphic, with half of all methylated genes being methylated in only one of the two ecotypes. A survey of loci in 96 Arabidopsis accessions revealed a similar degree of methylation polymorphism. Within-gene methylation is heritable, but is lost at a high frequency in segregating F(2) families. Promoter methylation is rare, and gene expression is not generally affected by differences in DNA methylation. Small interfering RNA are preferentially associated with methylated TEs, but not with methylated genes, indicating that most genic methylation is not guided by small interfering RNA. This may account for the instability of gene methylation, if occasional failure of maintenance methylation cannot be restored by other means.

Arabidopsis TFL2/LHP1 specifically associates with genes marked by trimethylation of histone H3 lysine 27

TERMINAL FLOWER 2/LIKE HETEROCHROMATIN PROTEIN 1 (TFL2/LHP1) is the only Arabidopsis protein with overall sequence similarity to the HETEROCHROMATIN PROTEIN 1 (HP1) family of metazoans and S. pombe. TFL2/LHP1 represses transcription of numerous genes, including the flowering-time genes FLOWERING LOCUS T (FT) and FLOWERING LOCUS C (FLC), as well as the floral organ identity genes AGAMOUS (AG) and APETALA 3 (AP3). These genes are also regulated by proteins of the Polycomb repressive complex 2 (PRC2), and it has been proposed that TFL2/LHP1 represents a potential stabilizing factor of PRC2 activity. Here we show by chromatin immunoprecipitation and hybridization to an Arabidopsis Chromosome 4 tiling array (ChIP-chip) that TFL2/LHP1 associates with hundreds of small domains, almost all of which correspond to genes located within euchromatin. We investigated the chromatin marks to which TFL2/LHP1 binds and show that, in vitro, TFL2/LHP1 binds to histone H3 di- or tri-methylated at lysine 9 (H3K9me2 or H3K9me3), the marks recognized by HP1, and to histone H3 trimethylated at lysine 27 (H3K27me3), the mark deposited by PRC2. However, in vivo TFL2/LHP1 association with chromatin occurs almost exclusively and co-extensively with domains marked by H3K27me3, but not H3K9me2 or -3. Moreover, the distribution of H3K27me3 is unaffected in lhp1 mutant plants, indicating that unlike PRC2 components, TFL2/LHP1 is not involved in the deposition of this mark. Rather, our data suggest that TFL2/LHP1 recognizes specifically H3K27me3 in vivo as part of a mechanism that represses the expression of many genes targeted by PRC2.

Stable repression of gene expression is an important aspect of the developmental programs of higher organisms. In plants and animals, DNA is organized within chromatin, which contains at its core a set of evolutionarily conserved proteins called histones. These proteins can be modified for example by methylation or acetylation of lysines or phosphorylation of serines. Specific combinations of these histone modifications are interpreted by other chromatin proteins and thereby play essential roles in gene regulation. One such potential effector of the histone code in the flowering plant Arabidopsis is TERMINAL FLOWER 2/LIKE HETEROCHROMATIN PROTEIN 1 (TFL2/LHP1). Here we present highly detailed “epigenomic” maps that establish that TFL2/LHP1 associates with a subset of Arabidopsis genes that are marked by tri-methylation of Lysine 27 of histone H3. In plants and animals, an evolutionarily conserved complex called PRC2 deposits this mark. In Drosophila and mammals this modified histone is then read by another complex, called PRC1, to maintain the stable repression of genes. In Arabidopsis however, no PRC1 complex exists, and our results provide evidence that TFL2/LHP1 may fulfill a related function.

DNA Microarray Analyses in Higher Plants

DNA microarrays were originally devised and described as a convenient technology for the global analysis of plant gene expression. Over the past decade, their use has expanded enormously to cover all kingdoms of living organisms. At the same time, the scope of applications of microarrays has increased beyond expression analyses, with plant genomics playing a leadership role in the on-going development of this technology. As the field has matured, the rate-limiting step has moved from that of the technical process of data generation to that of data analysis. We currently face major problems in dealing with the accumulating datasets, not simply with respect to how to archive, access, and process the huge amounts of data that have been and are being produced, but also in determining the relative quality of the different datasets. A major recognized concern is the appropriate use of statistical design in microarray experiments, without which the datasets are rendered useless. A vigorous area of current research involves the development of novel statistical tools specifically for microarray experiments. This article describes, in a necessarily selective manner, the types of platforms currently employed in microarray research and provides an overview of recent activities using these platforms in plant biology.

An RNA-dependent RNA polymerase is required for paramutation in maize

Paramutation is an allele-dependent transfer of epigenetic information, which results in the heritable silencing of one allele by another. Paramutation at the b1 locus in maize is mediated by unique tandem repeats that communicate in trans to establish and maintain meiotically heritable transcriptional silencing. The mop1 (mediator of paramutation1) gene is required for paramutation, and mop1 mutations reactivate silenced Mutator elements. Plants carrying mutations in the mop1 gene also stochastically exhibit pleiotropic developmental phenotypes. Here we report the map-based cloning of mop1, an RNA-dependent RNA polymerase gene (RDRP), most similar to the RDRP in plants that is associated with the production of short interfering RNA (siRNA) targeting chromatin. Nuclear run-on assays reveal that the tandem repeats required for b1 paramutation are transcribed from both strands, but siRNAs were not detected. We propose that the mop1 RDRP is required to maintain a threshold level of repeat RNA, which functions in trans to establish and maintain the heritable chromatin states associated with paramutation.

Whole-Genome Microarray in Arabidopsis Facilitates Global Analysis of Retained Introns

Alternative splicing (AS) is an important post-transcriptional regulatory mechanism that can increase protein diversity and affect mRNA stability. Different types of AS have been observed; these include exon skipping, alternative donor or acceptor site and intron retention. In humans, exon skipping is the most common type while intron retention is rare. In contrast, in Arabidopsis, intron retention is the most prevalent AS type (approximately 40%). Here we show that direct transcript expression analysis using high-density oligonucleotide-based whole-genome microarrays (WGAs) is particularly amenable for assessing global intron retention in Arabidopsis. By applying a novel algorithm retained introns are detected in 8% of the transcripts examined. A sampling of 14 transcripts showed that 86% can be confirmed by RT-PCR. This rate of detection predicts an overall total AS rate of 20% for Arabidopsis compared with 10-22% based on EST/cDNA-based analysis. These findings will facilitate monitoring constitutive and dynamic whole-genome splicing on the next generation WGA slides.

The Arabidopsis genome: A foundation for plant research

The sequence of the first plant genome was completed and published at the end of 2000. This spawned a series of large-scale projects aimed at discovering the functions of the 25,000+ genes identified in Arabidopsis thaliana (Arabidopsis). This review summarizes progress made in the past five years and speculates about future developments in Arabidopsis research and its implications for crop science. The provision of large populations of gene disruption lines to the research community has greatly accelerated the impact of genomics on many areas of plant science. The tools and community organization required for plant integrative and systems biology approaches are now ready to accomplish the next big step in plant biology--the integration of knowledge and modeling of biological processes. In the future, plant science will continue to be enriched by the alignment of high-quality basic research (generally conducted in Arabidopsis), with strategic objectives in crop plants. The sequence and analysis of an increasing number of crop plant genomes enhance this alignment and provide new insights into genome evolution and crop plant domestication.