Transcriptional analysis of highly syntenic regions between Medicago truncatula and Glycine max using tiling microarrays.

BACKGROUND: Legumes are the third largest family of flowering plants and are unique among crop species in their ability to fix atmospheric nitrogen. As a result of recent genome sequencing efforts, legumes are now one of a few plant families with extensive genomic and transcriptomic data available in multiple species. The unprecedented complexity and impending completeness of these data create opportunities for new approaches to discovery. RESULTS: We report here a transcriptional analysis in six different organ types of syntenic regions totaling approximately 1 Mb between the legume plants barrel medic (Medicago truncatula) and soybean (Glycine max) using oligonucleotide tiling microarrays. This analysis detected transcription of over 80% of the predicted genes in both species. We also identified 499 and 660 transcriptionally active regions from barrel medic and soybean, respectively, over half of which locate outside of the predicted exons. We used the tiling array data to detect differential gene expression in the six examined organ types and found several genes that are preferentially expressed in the nodule. Further investigation revealed that some collinear genes exhibit different expression patterns between the two species. CONCLUSION: These results demonstrate the utility of genome tiling microarrays in generating transcriptomic data to complement computational annotation of the newly available legume genome sequences. The tiling microarray data was further used to quantify gene expression levels in multiple organ types of two related legume species. Further development of this method should provide a new approach to comparative genomics aimed at elucidating genome organization and transcriptional regulation.

High-resolution mapping of epigenetic modifications of the rice genome uncovers interplay between DNA methylation, histone methylation, and gene expression.

We present high-resolution maps of DNA methylation and H3K4 di- and trimethylation of two entire chromosomes and two fully sequenced centromeres in rice (Oryza sativa) shoots and cultured cells. This analysis reveals combinatorial interactions between these epigenetic modifications and chromatin structure and gene expression. Cytologically densely stained heterochromatin had less H3K4me2 and H3K4me3 and more methylated DNA than the less densely stained euchromatin, whereas centromeres had a unique epigenetic composition. Most transposable elements had highly methylated DNA but no H3K4 methylation, whereas more than half of protein-coding genes had both methylated DNA and di- and/or trimethylated H3K4. Methylation of DNA but not H3K4 was correlated with suppressed transcription. By contrast, when both DNA and H3K4 were methylated, transcription was only slightly reduced. Transcriptional activity was positively correlated with the ratio of H3K4me3/H3K4me2: genes with predominantly H3K4me3 were actively transcribed, whereas genes with predominantly H3K4me2 were transcribed at moderate levels. More protein-coding genes contained all three modifications, and more transposons contained DNA methylation in shoots than cultured cells. Differential epigenetic modifications correlated to tissue-specific expression between shoots and cultured cells. Collectively, this study provides insights into the rice epigenomes and their effect on gene expression and plant development.

In-depth query of large genomes using tiling arrays.

Identification of the transcribed regions in the newly sequenced genomes is one of the major challenges of postgenomic biology. Among different alternatives for empirical transcriptome mapping, whole-genome tiling array experiment emerged as the most comprehensive and unbiased approach. This relatively new method uses high-density oligonucleotide arrays with probes chosen uniformly from both strands of the entire genomes including all genic and intergenic regions. By hybridizing the arrays with tissue specific or pooled RNA samples, a genome-wide picture of transcription can be derived. This chapter discusses computational tools and techniques necessary to successfully conduct genome tiling array experiments.

Global identification and characterization of transcriptionally active regions in the rice genome.

Genome tiling microarray studies have consistently documented rich transcriptional activity beyond the annotated genes. However, systematic characterization and transcriptional profiling of the putative novel transcripts on the genome scale are still lacking. We report here the identification of 25,352 and 27,744 transcriptionally active regions (TARs) not encoded by annotated exons in the rice (Oryza. sativa) subspecies japonica and indica, respectively. The non-exonic TARs account for approximately two thirds of the total TARs detected by tiling arrays and represent transcripts likely conserved between japonica and indica. Transcription of 21,018 (83%) japonica non-exonic TARs was verified through expression profiling in 10 tissue types using a re-array in which annotated genes and TARs were each represented by five independent probes. Subsequent analyses indicate that about 80% of the japonica TARs that were not assigned to annotated exons can be assigned to various putatively functional or structural elements of the rice genome, including splice variants, uncharacterized portions of incompletely annotated genes, antisense transcripts, duplicated gene fragments, and potential non-coding RNAs. These results provide a systematic characterization of non-exonic transcripts in rice and thus expand the current view of the complexity and dynamics of the rice transcriptome.

Analysis of transcription factor HY5 genomic binding sites revealed its hierarchical role in light regulation of development.

The transcription factor LONG HYPOCOTYL5 (HY5) acts downstream of multiple families of the photoreceptors and promotes photomorphogenesis. Although it is well accepted that HY5 acts to regulate target gene expression, in vivo binding of HY5 to any of its target gene promoters has yet to be demonstrated. Here, we used a chromatin immunoprecipitation procedure to verify suspected in vivo HY5 binding sites. We demonstrated that in vivo association of HY5 with promoter targets is not altered under distinct light qualities or during light-to-dark transition. Coupled with DNA chip hybridization using a high-density 60-nucleotide oligomer microarray that contains one probe for every 500 nucleotides over the entire Arabidopsis thaliana genome, we mapped genome-wide in vivo HY5 binding sites. This analysis showed that HY5 binds preferentially to promoter regions in vivo and revealed >3000 chromosomal sites as putative HY5 binding targets. HY5 binding targets tend to be enriched in the early light-responsive genes and transcription factor genes. Our data thus support a model in which HY5 is a high hierarchical regulator of the transcriptional cascades for photomorphogenesis.

A core transcriptional network for early mesoderm development in Drosophila melanogaster.

Embryogenesis is controlled by large gene-regulatory networks, which generate spatially and temporally refined patterns of gene expression. Here, we report the characteristics of the regulatory network orchestrating early mesodermal development in the fruitfly Drosophila, where the transcription factor Twist is both necessary and sufficient to drive development. Through the integration of chromatin immunoprecipitation followed by microarray analysis (ChIP-on-chip) experiments during discrete time periods with computational approaches, we identified >2000 Twist-bound cis-regulatory modules (CRMs) and almost 500 direct target genes. Unexpectedly, Twist regulates an almost complete cassette of genes required for cell proliferation in addition to genes essential for morophogenesis and cell migration. Twist targets almost 25% of all annotated Drosophila transcription factors, which may represent the entire set of regulators necessary for the early development of this system. By combining in vivo binding data from Twist, Mef2, Tinman, and Dorsal we have constructed an initial transcriptional network of early mesoderm development. The network topology reveals extensive combinatorial binding, feed-forward regulation, and complex logical outputs as prevalent features. In addition to binary activation and repression, we suggest that Twist binds to almost all mesodermal CRMs to provide the competence to integrate inputs from more specialized transcription factors.

The transcriptome of the sea urchin embryo.

The sea urchin Strongylocentrotus purpuratus is a model organism for study of the genomic control circuitry underlying embryonic development. We examined the complete repertoire of genes expressed in the S. purpuratus embryo, up to late gastrula stage, by means of high-resolution custom tiling arrays covering the whole genome. We detected complete spliced structures even for genes known to be expressed at low levels in only a few cells. At least 11,000 to 12,000 genes are used in embryogenesis. These include most of the genes encoding transcription factors and signaling proteins, as well as some classes of general cytoskeletal and metabolic proteins, but only a minor fraction of genes encoding immune functions and sensory receptors. Thousands of small asymmetric transcripts of unknown function were also detected in intergenic regions throughout the genome. The tiling array data were used to correct and authenticate several thousand gene models during the genome annotation process.

Global identification of noncoding RNAs in Saccharomyces cerevisiae by modulating an essential RNA processing pathway.

Noncoding RNAs (ncRNAs) perform essential cellular tasks and play key regulatory roles in all organisms. Although several new ncRNAs in yeast were recently discovered by individual studies, to our knowledge no comprehensive empirical search has been conducted. We demonstrate a powerful and versatile method for global identification of previously undescribed ncRNAs by modulating an essential RNA processing pathway through the depletion of a key ribonucleoprotein enzyme component, and monitoring differential transcriptional activities with genome tiling arrays during the time course of the ribonucleoprotein depletion. The entire Saccharomyces cerevisiae genome was scanned during cell growth decay regulated by promoter-mediated depletion of Rpp1, an essential and functionally conserved protein component of the RNase P enzyme. In addition to most verified genes and ncRNAs, expression was detected in 98 antisense and intergenic regions, 74 that were further confirmed to contain previously undescribed RNAs. A class of ncRNAs, located antisense to coding regions of verified protein-coding genes, is discussed in this article. One member, HRA1, is likely involved in 18S rRNA maturation.

Genome-wide transcription analyses in rice using tiling microarrays.

Sequencing and computational annotation revealed several features, including high gene numbers, unusual composition of the predicted genes and a large number of genes lacking homology to known genes, that distinguish the rice (Oryza sativa) genome from that of other fully sequenced model species. We report here a full-genome transcription analysis of the indica rice subspecies using high-density oligonucleotide tiling microarrays. Our results provided expression data support for the existence of 35,970 (81.9%) annotated gene models and identified 5,464 unique transcribed intergenic regions that share similar compositional properties with the annotated exons and have significant homology to other plant proteins. Elucidating and mapping of all transcribed regions revealed an association between global transcription and cytological chromosome features, and an overall similarity of transcriptional activity between duplicated segments of the genome. Collectively, our results provide the first whole-genome transcription map useful for further understanding the rice genome.

A pilot study of transcription unit analysis in rice using oligonucleotide tiling-path microarray.

As the international efforts to sequence the rice genome are completed, an immediate challenge and opportunity is to comprehensively and accurately define all transcription units in the rice genome. Here we describe a strategy of using high-density oligonucleotide tiling-path microarrays to map transcription of the japonica rice genome. In a pilot experiment to test this approach, one array representing the reverse strand of the last 11.2 Mb sequence of chromosome 10 was analyzed in detail based on a mathematical model developed in this study. Analysis of the array data detected 77% of the reference gene models in a mixture of four RNA populations. Moreover, significant transcriptional activities were found in many of the previously annotated intergenic regions. These preliminary results demonstrate the utility of genome tiling microarrays in evaluating annotated rice gene models and in identifying novel transcription units that will facilitate rice genome annotation.

Identification of transcribed sequences in Arabidopsis thaliana by using high-resolution genome tiling arrays.

Using a maskless photolithography method, we produced DNA oligonucleotide microarrays with probe sequences tiled throughout the genome of the plant Arabidopsis thaliana. RNA expression was determined for the complete nuclear, mitochondrial, and chloroplast genomes by tiling 5 million 36-mer probes. These probes were hybridized to labeled mRNA isolated from liquid grown T87 cells, an undifferentiated Arabidopsis cell culture line. Transcripts were detected from at least 60% of the nearly 26,330 annotated genes, which included 151 predicted genes that were not identified previously by a similar genome-wide hybridization study on four different cell lines. In comparison with previously published results with 25-mer tiling arrays produced by chromium masking-based photolithography technique, 36-mer oligonucleotide probes were found to be more useful in identifying intron-exon boundaries. Using two-dimensional HPLC tandem mass spectrometry, a small-scale proteomic analysis was performed with the same cells. A large amount of strongly hybridizing RNA was found in regions "antisense" to known genes. Similarity of antisense activities between the 25-mer and 36-mer data sets suggests that it is a reproducible and inherent property of the experiments. Transcription activities were also detected for many of the intergenic regions and the small RNAs, including tRNA, small nuclear RNA, small nucleolar RNA, and microRNA. Expression of tRNAs correlates with genome-wide amino acid usage.

Genome-wide transcriptional analysis of flagellar regeneration in Chlamydomonas reinhardtii identifies orthologs of ciliary disease genes.

The important role that cilia and flagella play in human disease creates an urgent need to identify genes involved in ciliary assembly and function. The strong and specific induction of flagellar-coding genes during flagellar regeneration in Chlamydomonas reinhardtii suggests that transcriptional profiling of such cells would reveal new flagella-related genes. We have conducted a genome-wide analysis of RNA transcript levels during flagellar regeneration in Chlamydomonas by using maskless photolithography method-produced DNA oligonucleotide microarrays with unique probe sequences for all exons of the 19,803 predicted genes. This analysis represents previously uncharacterized whole-genome transcriptional activity profiling study in this important model organism. Analysis of strongly induced genes reveals a large set of known flagellar components and also identifies a number of important disease-related proteins as being involved with cilia and flagella, including the zebrafish polycystic kidney genes Qilin, Reptin, and Pontin, as well as the testis-expressed tubby-like protein TULP2.

Global identification of human transcribed sequences with genome tiling arrays.

Elucidating the transcribed regions of the genome constitutes a fundamental aspect of human biology, yet this remains an outstanding problem. To comprehensively identify coding sequences, we constructed a series of high-density oligonucleotide tiling arrays representing sense and antisense strands of the entire nonrepetitive sequence of the human genome. Transcribed sequences were located across the genome via hybridization to complementary DNA samples, reverse-transcribed from polyadenylated RNA obtained from human liver tissue. In addition to identifying many known and predicted genes, we found 10,595 transcribed sequences not detected by other methods. A large fraction of these are located in intergenic regions distal from previously annotated genes and exhibit significant homology to other mammalian proteins.

A gene expression map for the euchromatic genome of Drosophila melanogaster.

We used a maskless photolithography method to produce DNA oligonucleotide microarrays with unique probe sequences tiled throughout the genome of Drosophila melanogaster and across predicted splice junctions. RNA expression of protein coding and nonprotein coding sequences was determined for each major stage of the life cycle, including adult males and females. We detected transcriptional activity for 93% of annotated genes and RNA expression for 41% of the probes in intronic and intergenic sequences. Comparison to genome-wide RNA interference data and to gene annotations revealed distinguishable levels of expression for different classes of genes and higher levels of expression for genes with essential cellular functions. Differential splicing was observed in about 40% of predicted genes, and 5440 previously unknown splice forms were detected. Genes within conserved regions of synteny with D. pseudoobscura had highly correlated expression; these regions ranged in length from 10 to 900 kilobase pairs. The expressed intergenic and intronic sequences are more likely to be evolutionarily conserved than nonexpressed ones, and about 15% of them appear to be developmentally regulated. Our results provide a draft expression map for the entire nonrepetitive genome, which reveals a much more extensive and diverse set of expressed sequences than was previously predicted.

Characterization of synthetic DNA bar codes in Saccharomyces cerevisiae gene-deletion strains.

Incorporation of strain-specific synthetic DNA tags into yeast Saccharomyces cerevisiae gene-deletion strains has enabled identification of gene functions by massively parallel growth rate analysis. However, it is important to confirm the sequences of these tags, because mutations introduced during construction could lead to significant errors in hybridization performance. To validate this experimental system, we sequenced 11,812 synthetic 20-mer molecular bar codes and adjacent sequences (>1.8 megabases synthetic DNA) by pyrosequencing and Sanger methods. At least 31% of the genome-integrated 20-mer tags contain differences from those originally synthesized. However, these mutations result in anomalous hybridization in only a small subset of strains, and the sequence information enables redesign of hybridization probes for arrays. The robust performance of the yeast gene-deletion dual oligonucleotide bar-code design in array hybridization validates the use of molecular bar codes in living cells for tracking their growth phenotype.