Utilizing tiling microarrays for whole-genome analysis in plants

An Overview of Methodologies in Studying lncRNAs in the High-Throughput Era: When Acronyms ATTACK!

The discovery of pervasive transcription in eukaryotic genomes provided one of many surprising (and perhaps most surprising) findings of the genomic era and led to the uncovering of a large number of previously unstudied transcriptional events. This pervasive transcription leads to the production of large numbers of noncoding RNAs (ncRNAs) and thus opened the window to study these diverse, abundant transcripts of unclear relevance and unknown function. Since that discovery, recent advances in high-throughput sequencing technologies have identified a large collection of ncRNAs, from microRNAs to long noncoding RNAs (lncRNAs). Subsequent discoveries have shown that many lncRNAs play important roles in various eukaryotic processes; these discoveries have profoundly altered our understanding of the regulation of eukaryotic gene expression. Although the identification of ncRNAs has become a standard experimental approach, the functional characterization of these diverse ncRNAs remains a major challenge. In this chapter, we highlight recent progress in the methods to identify lncRNAs and the techniques to study the molecular function of these lncRNAs and the application of these techniques to the study of plant lncRNAs.

Global nucleosome positioning regulates salicylic acid mediated transcription in Arabidopsis thaliana

BACKGROUND

The nucleosome positioning regulates the gene expression and many other DNA-related processes in eukaryotes. Genome-wide mapping of nucleosome positions and correlation of genome-wide nucleosomal remodeling with the changes in the gene expression can help us understanding gene regulation on genome level.

RESULTS

In the present study, we correlate the gene expression and the genomic nucleosomal remodeling in response to salicylic acid (SA) treatment in A. thaliana. We have mapped genome-wide nucleosomes by performing tiling microarray using 146 bp mononucleosomal template DNA. The average nucleosomal coverage is approximately 346 bp per nucleosome both under the control and the SA-treated conditions. The nucleosomal coverage is more in the coding region than in the 5' regulatory regions. We observe approximately 50% nucleosomal remodeling on SA treatment where significant nucleosomal depletion and nucleosomal enrichment around the transcription start site (TSS) occur in SA induced genes and SA repressed genes respectively in response to SA treatment. Especially in the case of the SA-induced group, the nucleosomal remodeling over the minimal promoter in response to SA is especially significant in the Non-expresser of PR1 (NPR1)-dependent genes. A detailed investigation of npr1-1 mutant confirms a distinct role of NPR1 in the nucleosome remodeling over the core promoter. We have also identified several motifs for various hormonal responses; including ABRE elements in the remodeled nucleosomal regions around the promoter region in the SA regulated genes. We have further identified that the W-box and TGACG/C motif, reported to play an important role in SA-mediated induction, are enriched in nucleosome free regions (NFRs) of the promoter region of the SA induced genes.

CONCLUSIONS

This is the first study reporting genome-wide effects of SA treatment on the chromatin architecture of A. thaliana. It also reports significant role of NPR1 in genome-wide nucleosomal remodeling in response to SA.

Epigenetic landscape of germline specific genes in the sporophyte cells of Arabidopsis thaliana

In plants, the germline lineages arise in later stages of life cycle as opposed to animals where both male and female germlines are set aside early in development. This developmental divergence is associated with germline specific or preferential expression of a subset of genes that are normally repressed for the rest of plant life cycle. The gene regulatory mechanisms involved in such long-term suppression and short-term activation in plant germline remain vague. Thus, we explored the nature of epigenetic marks that are likely associated with long-term gene repression in the non-germline cells. We accessed available Arabidopsis genome-wide DNA methylation and histone modification data and queried it for epigenetic marks associated with germline genes: genes preferentially expressed in sperm cells, egg cells, synergid cells, central cells, antipodal cells or embryo sac or genes that are with enriched expression in two or more of female germline tissues. The vast majority of germline genes are associated with repression-related epigenetic histone modifications in one or more non-germline tissues, among which H3K9me2 and H3K27me3 are the most widespread repression-related marks. Interestingly, we show here that the repressive epigenetic mechanisms differ between male and female germline genes. We also highlight the diverse states of epigenetic marks in different non-germline tissues. Some germline genes also have activation-related marks in non-germline tissues, and the proportion of such genes is higher for female germline genes. Germline genes include 30 transposable element (TE) loci, to which a large number of 24-nt long small interfering RNAs were mapped, suggesting that these small RNAs take a role in suppressing them in non-germline tissues. The data presented here suggest that the majority of Arabidopsis gamete-preferentially/-enriched genes bear repressive epigenetic modifications or regulated by small RNAs.

Microarray-based ultra-high resolution discovery of genomic deletion mutations

BACKGROUND

Oligonucleotide microarray-based comparative genomic hybridization (CGH) offers an attractive possible route for the rapid and cost-effective genome-wide discovery of deletion mutations. CGH typically involves comparison of the hybridization intensities of genomic DNA samples with microarray chip representations of entire genomes, and has widespread potential application in experimental research and medical diagnostics. However, the power to detect small deletions is low.

RESULTS

Here we use a graduated series of Arabidopsis thaliana genomic deletion mutations (of sizes ranging from 4 bp to ~5 kb) to optimize CGH-based genomic deletion detection. We show that the power to detect smaller deletions (4, 28 and 104 bp) depends upon oligonucleotide density (essentially the number of genome-representative oligonucleotides on the microarray chip), and determine the oligonucleotide spacings necessary to guarantee detection of deletions of specified size.

CONCLUSIONS

Our findings will enhance a wide range of research and clinical applications, and in particular will aid in the discovery of genomic deletions in the absence of a priori knowledge of their existence.

Mapping genomic features of tiling microarray data by TileMapper

The recent revolution of genomics techniques has allowed the detection of various sequence features and biological variations on whole-genome scale. However, these high-resolution data present significant challenges for experimental biologists to understand and analyze. The conventional way is to use genome browsers to locate and visualize regions of interest. But it lacks user-friendly data mining functionality. Here we present a protocol that allows rapid annotation of genomic coordinate data by using TileMapper. Interesting biological annotations from large-scale genomic data, such as transcriptome analysis, chromatin immunoprecipitation on chip, or methyl-DNA immunoprecipitation (MeDIP) studies generated from the tiling microarrays and other platforms, could be analyzed without requiring computational skills. The outputs are saved in tabulated format, which permit flexible and simple processing in spreadsheet software, or to be exported to other pipelines for subsequent analysis.

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.

A brief introduction to tiling microarrays: principles, concepts, and applications

Technological achievements have always contributed to the advancement of biomedical research. It has never been more so than in recent times, when the development and application of innovative cutting-edge technologies have transformed biology into a data-rich quantitative science. This stunning revolution in biology primarily ensued from the emergence of microarrays over two decades ago. The completion of whole-genome sequencing projects and the advance in microarray manufacturing technologies enabled the development of tiling microarrays, which gave unprecedented genomic coverage. Since their first description, several types of application of tiling arrays have emerged, each aiming to tackle a different biological problem. Although numerous algorithms have already been developed to analyze microarray data, new method development is still needed not only for better performance but also for integration of available microarray data sets, which without doubt constitute one of the largest collections of biological data ever generated. In this chapter we first introduce the principles behind the emergence and the development of tiling microarrays, and then discuss with some examples how they are used to investigate different biological problems.

Comparative co-expression analysis in plant biology

The analysis of gene expression data generated by high-throughput microarray transcript profiling experiments has shown that transcriptionally coordinated genes are often functionally related. Based on large-scale expression compendia grouping multiple experiments, this guilt-by-association principle has been applied to study modular gene programmes, identify cis-regulatory elements or predict functions for unknown genes in different model plants. Recently, several studies have demonstrated how, through the integration of gene homology and expression information, correlated gene expression patterns can be compared between species. The incorporation of detailed functional annotations as well as experimental data describing protein-protein interactions, phenotypes or tissue specific expression, provides an invaluable source of information to identify conserved gene modules and translate biological knowledge from model organisms to crops. In this review, we describe the different steps required to systematically compare expression data across species. Apart from the technical challenges to compute and display expression networks from multiple species, some future applications of plant comparative transcriptomics are highlighted.

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.

DNA methylation screening and analysis

DNA methylation is an epigenetic form of gene regulation that is universally important throughout the life course, especially during in utero and postnatal development. DNA methylation aids in cell cycle regulation and cellular differentiation processes. Previous studies have demonstrated that DNA methylation profiles may be altered by diet and the environment, and that these profiles are especially vulnerable during development. Thus, it is important to understand the role of DNA methylation in developmental governance and subsequent disease progression. A variety of molecular methods exist to assay for global, gene-specific, and epigenome-wide methylation. Here we describe these methods and discuss their relative strengths and limitations.

A 3,000-loci transcription map of chromosome 3B unravels the structural and functional features of gene islands in hexaploid wheat

To improve our understanding of the organization and regulation of the wheat (Triticum aestivum) gene space, we established a transcription map of a wheat chromosome (3B) by hybridizing a newly developed wheat expression microarray with bacterial artificial chromosome pools from a new version of the 3B physical map as well as with cDNA probes derived from 15 RNA samples. Mapping data for almost 3,000 genes showed that the gene space spans the whole chromosome 3B with a 2-fold increase of gene density toward the telomeres due to an increase in the number of genes in islands. Comparative analyses with rice (Oryza sativa) and Brachypodium distachyon revealed that these gene islands are composed mainly of genes likely originating from interchromosomal gene duplications. Gene Ontology and expression profile analyses for the 3,000 genes located along the chromosome revealed that the gene islands are enriched significantly in genes sharing the same function or expression profile, thereby suggesting that genes in islands acquired shared regulation during evolution. Only a small fraction of these clusters of cofunctional and coexpressed genes was conserved with rice and B. distachyon, indicating a recent origin. Finally, genes with the same expression profiles in remote islands (coregulation islands) were identified suggesting long-distance regulation of gene expression along the chromosomes in wheat.

Identification of alternative splicing events regulated by an Arabidopsis serine/arginine-like protein, atSR45a, in response to high-light stress using a tiling array

We have demonstrated that an Arabidopsis serine/arginine rich-like protein, atSR45a, interacts with other splicing factors and its expression is markedly induced by high-light stress, suggesting the involvement of atSR45a in the regulation of stress-responsive alternative splicing. A whole-genome tiling array identified the alternative splicing of genes regulated by atSR45a by comparing gene expression profiles in wild-type and knockout atSR45a (KO-sr45a) plants under high-light stress. The expression levels of genomic regions within 217 genes were significantly altered in the KO-sr45a plants compared with the wild-type plants. Many genes encoded factors involved in signal transduction, cell cycle and DNA processing, protein fate and transcription. A semi-quantitative reverse transcription-PCR (RT-PCR) analysis confirmed changes in the transcript levels and/or alternative splicing efficiency under high-light stress in 18 genes, suggesting that atSR45a affects directly or indirectly not only alternative splicing efficiency but also the transcription of these target genes. Changes in the expression of atSR45a in response to high-light stress temporally correlated with changes in the alternative splicing efficiency and transcript levels of three and one target genes, respectively. Sequencing of the alternatively spliced variants of three target genes showed that atSR45a suppresses the splicing efficiency of intron retention-type alternative splicing events. These findings indicated the importance of atSR45a to the diversification of the transcriptome under high-light stress.

Internet Resources for Gene Expression Analysis in Arabidopsis thaliana

The number of online databases and web-tools for gene expression analysis in Arabidopsis thaliana has increased tremendously during the last years. These resources permit the database-assisted identification of putative cis-regulatory DNA sequences, their binding proteins, and the determination of common cis-regulatory motifs in coregulated genes. DNA binding proteins may be predicted by the type of cis-regulatory motif. Further questions of combinatorial control based on the interaction of DNA binding proteins and the colocalization of cis-regulatory motifs can be addressed. The database-assisted spatial and temporal expression analysis of DNA binding proteins and their target genes may help to further refine experimental approaches. Signal transduction pathways upstream of regulated genes are not yet fully accessible in databases mainly because they need to be manually annotated. This review focuses on the use of the AthaMap and PathoPlant^® databases for gene expression regulation analysis and discusses similar and complementary online databases and web-tools. Online databases are helpful for the development of working hypothesis and for designing subsequent experiments.

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.

Digital gene expression signatures for maize development

Genome-wide expression signatures detect specific perturbations in developmental programs and contribute to functional resolution of key regulatory networks. In maize (Zea mays) inflorescences, mutations in the RAMOSA (RA) genes affect the determinacy of axillary meristems and thus alter branching patterns, an important agronomic trait. In this work, we developed and tested a framework for analysis of tag-based, digital gene expression profiles using Illumina's high-throughput sequencing technology and the newly assembled B73 maize reference genome. We also used a mutation in the RA3 gene to identify putative expression signatures specific to stem cell fate in axillary meristem determinacy. The RA3 gene encodes a trehalose-6-phosphate phosphatase and may act at the interface between developmental and metabolic processes. Deep sequencing of digital gene expression libraries, representing three biological replicate ear samples from wild-type and ra3 plants, generated 27 million 20- to 21-nucleotide reads with frequencies spanning 4 orders of magnitude. Unique sequence tags were anchored to 3'-ends of individual transcripts by DpnII and NlaIII digests, which were multiplexed during sequencing. We mapped 86% of nonredundant signature tags to the maize genome, which associated with 37,117 gene models and unannotated regions of expression. In total, 66% of genes were detected by at least nine reads in immature maize ears. We used comparative genomics to leverage existing information from Arabidopsis (Arabidopsis thaliana) and rice (Oryza sativa) in functional analyses of differentially expressed maize genes. Results from this study provide a basis for the analysis of short-read expression data in maize and resolved specific expression signatures that will help define mechanisms of action for the RA3 gene.

Seasonal and developmental timing of flowering

The coordination of the timing of flowering with seasonal and development cues is a critical life-history trait that has been shaped by evolution to maximize reproductive success. Decades of studying many plant species have revealed several of the fascinating systems that plants have evolved to control flowering time: such as the perception of day length in leaves, which leads to the production of a mobile signal, florigen, that promotes flowering at the shoot apical meristem; the vernalization process in which exposure to prolonged cold results in meristem competence to flower; and the juvenile to adult phase transition. Arabidopsis research has contributed greatly to understanding these systems at a molecular level.

The Listeria transcriptional landscape from saprophytism to virulence

The bacterium Listeria monocytogenes is ubiquitous in the environment and can lead to severe food-borne infections. It has recently emerged as a multifaceted model in pathogenesis. However, how this bacterium switches from a saprophyte to a pathogen is largely unknown. Here, using tiling arrays and RNAs from wild-type and mutant bacteria grown in vitro, ex vivo and in vivo, we have analysed the transcription of its entire genome. We provide the complete Listeria operon map and have uncovered far more diverse types of RNAs than expected: in addition to 50 small RNAs (<500 nucleotides), at least two of which are involved in virulence in mice, we have identified antisense RNAs covering several open-reading frames and long overlapping 5' and 3' untranslated regions. We discovered that riboswitches can act as terminators for upstream genes. When Listeria reaches the host intestinal lumen, an extensive transcriptional reshaping occurs with a SigB-mediated activation of virulence genes. In contrast, in the blood, PrfA controls transcription of virulence genes. Remarkably, several non-coding RNAs absent in the non-pathogenic species Listeria innocua exhibit the same expression patterns as the virulence genes. Together, our data unravel successive and coordinated global transcriptional changes during infection and point to previously unknown regulatory mechanisms in bacteria.

Formalin can alter the intracellular localization of some transcription factors in Saccharomyces cerevisiae

Indirect immunofluorescence (IF) microscopy is a frequently used method to determine intracellular protein localization. It is especially useful for low abundance proteins, for example the GATA-factors (Gln3, Gat1) which activate nitrogen catabolite repression (NCR)-sensitive transcription. Limiting nitrogen or treating cells with Tor pathway inhibitor, rapamycin, elicits nuclear GATA-factor localization and increased NCR-sensitive transcription, whereas excess nitrogen restricts these proteins to the cytoplasm and decreases transcription. The initial step of the IF procedure is formalin-fixation that quenches cellular activity and fixes protein locations via cross-linking. We find that under some conditions, formalin itself can influence GATA-factor localization. With low formalin (0.8% or 1.6%), Gat1-Myc(13) became more nuclear, and with higher concentrations (5.6%), it became more cytoplasmic. Gln3-Myc(13) localization, on the other hand, did not respond to low formalin, but became more cytoplasmic at the higher concentration. Interestingly, the high concentration of formalin had no demonstrable effect when the GATA factors were completely nuclear, i.e. after rapamycin (Gat1-Myc(13)) or Msx (Gln3-Myc(13)) treatment. These effects are most likely elicited by polyoxymethylene glycols, which significantly increase the osmolarity of the medium (0.5-2). We suggest that varying degrees of osmotic stress and transcription factor movement in response to it can occur after the beginning of fixation but before proteins become immobilized.

Plant pathogenesis-related (PR) proteins: a focus on PR peptides

The novel classes of plant pathogenesis-related (PR) proteins identified during the last decade also include novel peptide families. This review specifically focuses on these pathogenesis-related peptides, including proteinase inhibitors (PR-6 family), plant defensins (PR-12 family), thionins (PR-13 family) and lipid transfer proteins (PR-14 family). For each family of PR peptides, the general features concerning occurrence, expression and possible functions of their members are described. Next, more specifically the occurrence of each PR peptide family in the model plant Arabidopsis thaliana is discussed. Single-gene studies performed on particular gene members of a PR peptide family are reported. In addition, expression data of yet undescribed gene members of that particular PR peptide family are presented by consultation of publicly available micro-array databases. Finally an update is provided on the potential role of these PR peptides in A. thaliana, with a focus on their possible involvement in plant defense.

Use, tolerance and avoidance of amplified RNA silencing by plants

In plants and several other organisms, the effects of RNA silencing can be amplified by the action of cellular RNA-DEPENDENT RNA POLYMERASES (RDRs). These enzymes were primarily studied for their role in antiviral defense in plants, but it is becoming increasingly apparent that they also have important endogenous functions, including the control of chromatin structure and the regulation of cellular gene expression. Recent evidence suggests that endogenous RDR activities intercept several RNA quality control pathways that normally prevent or restrain widespread amplification of silencing, which is likely to be detrimental. Plants appear, however, to have evolved sophisticated measures to tolerate or exploit amplified silencing under specific biological circumstances.