Sequence analysis of a 40-kb Arabidopsis thaliana genomic region located at the top of chromosome 1

Analysis of antisense expression by whole genome tiling microarrays and siRNAs suggests mis-annotation of Arabidopsis orphan protein-coding genes

Background

MicroRNAs (miRNAs) and trans-acting small-interfering RNAs (tasi-RNAs) are small (20–22 nt long) RNAs (smRNAs) generated from hairpin secondary structures or antisense transcripts, respectively, that regulate gene expression by Watson-Crick pairing to a target mRNA and altering expression by mechanisms related to RNA interference. The high sequence homology of plant miRNAs to their targets has been the mainstay of miRNA prediction algorithms, which are limited in their predictive power for other kingdoms because miRNA complementarity is less conserved yet transitive processes (production of antisense smRNAs) are active in eukaryotes. We hypothesize that antisense transcription and associated smRNAs are biomarkers which can be computationally modeled for gene discovery.

Principal Findings

We explored rice (Oryza sativa) sense and antisense gene expression in publicly available whole genome tiling array transcriptome data and sequenced smRNA libraries (as well as C. elegans) and found evidence of transitivity of MIRNA genes similar to that found in Arabidopsis. Statistical analysis of antisense transcript abundances, presence of antisense ESTs, and association with smRNAs suggests several hundred Arabidopsis ‘orphan’ hypothetical genes are non-coding RNAs. Consistent with this hypothesis, we found novel Arabidopsis homologues of some MIRNA genes on the antisense strand of previously annotated protein-coding genes. A Support Vector Machine (SVM) was applied using thermodynamic energy of binding plus novel expression features of sense/antisense transcription topology and siRNA abundances to build a prediction model of miRNA targets. The SVM when trained on targets could predict the “ancient” (deeply conserved) class of validated Arabidopsis MIRNA genes with an accuracy of 84%, and 76% for “new” rapidly-evolving MIRNA genes.

Conclusions

Antisense and smRNA expression features and computational methods may identify novel MIRNA genes and other non-coding RNAs in plants and potentially other kingdoms, which can provide insight into antisense transcription, miRNA evolution, and post-transcriptional gene regulation.

A natural antisense transcript of the Petunia hybrida Sho gene suggests a role for an antisense mechanism in cytokinin regulation

The Sho gene from Petunia hybrida encodes an enzyme responsible for the synthesis of plant cytokinins. The 3' region of the Sho gene contains a promoter in the opposite orientation that produces a partially overlapping antisense transcript. Although Sho expression varies significantly in individual cell types, the sense and antisense transcript levels maintain a stable ratio in most tissue types. In reporter lines for the antisense promoter, we observed a change in antisense promoter activity in newly formed tissue that had been induced by prolonged culture on cytokinins or following decapitation. We interpret these data as a reflection of tissue-specific threshold levels for activation of the antisense transcript. In all tissue types tested, we detect a pool of antisense RNA of approximately 35 nt, which derives from the region where Sho sense and antisense transcripts overlap. We detect a second pool of putative dsRNA breakdown products of approximately 24 nt in all tissues tested, except roots, which are the main source of cytokinin synthesis. Our data suggest that antisense transcription can be activated in a tissue-specific manner to adjust local cytokinin synthesis via degradation of Sho dsRNA. We therefore propose that, in addition to cytokinin transport and inactivation, regulation of local cytokinin synthesis via antisense transcription represents yet another device for the complex control of local cytokinin levels in plants.

Arabidopsis transcriptome profiling indicates that multiple regulatory pathways are activated during cold acclimation in addition to the CBF cold response pathway

Many plants, including Arabidopsis, increase in freezing tolerance in response to low, nonfreezing temperatures, a phenomenon known as cold acclimation. Previous studies established that cold acclimation involves rapid expression of the CBF transcriptional activators (also known as DREB1 proteins) in response to low temperature followed by induction of the CBF regulon (CBF-targeted genes), which contributes to an increase in freezing tolerance. Here, we present the results of transcriptome-profiling experiments indicating the existence of multiple low-temperature regulatory pathways in addition to the CBF cold response pathway. The transcript levels of approximately 8000 genes were determined at multiple times after plants were transferred from warm to cold temperature and in warm-grown plants that constitutively expressed CBF1, CBF2, or CBF3. A total of 306 genes were identified as being cold responsive, with transcripts for 218 genes increasing and those for 88 genes decreasing threefold or more at one or more time points during the 7-day experiment. These results indicate that extensive downregulation of gene expression occurs during cold acclimation. Of the cold-responsive genes, 48 encode known or putative transcription factors. Two of these, RAP2.1 and RAP2.6, were activated by CBF expression and thus presumably control subregulons of the CBF regulon. Transcriptome comparisons indicated that only 12% of the cold-responsive genes are certain members of the CBF regulon. Moreover, at least 28% of the cold-responsive genes were not regulated by the CBF transcription factors, including 15 encoding known or putative transcription factors, indicating that these cold-responsive genes are members of different low-temperature regulons. Significantly, CBF expression at warm temperatures repressed the expression of eight genes that also were downregulated by low temperature, indicating that in addition to gene induction, gene repression is likely to play an integral role in cold acclimation.

The ANGUSTIFOLIA gene of Arabidopsis, a plant CtBP gene, regulates leaf-cell expansion, the arrangement of cortical microtubules in leaf cells and expression of a gene involved in cell-wall formation

We previously showed that the ANGUSTIFOLIA (AN) gene regulates the width of leaves of Arabidopsis thaliana, by controlling the polar elongation of leaf cells. In the present study, we found that the abnormal arrangement of cortical microtubules (MTs) in an leaf cells appeared to account entirely for the abnormal shape of the cells. It suggested that the AN gene might regulate the polarity of cell growth by controlling the arrangement of cortical MTs. We cloned the AN gene using a map-based strategy and identified it as the first member of the CtBP family to be found in plants. Wild-type AN cDNA reversed the narrow-leaved phenotype and the abnormal arrangement of cortical MTs of the an-1 mutation. In the animal kingdom, CtBPs self-associate and act as co-repressors of transcription. The AN protein can also self-associate in the yeast two-hybrid system. Furthermore, microarray analysis suggested that the AN gene might regulate the expression of certain genes, e.g. the gene involved in formation of cell walls, MERI5. A discussion of the molecular mechanisms involved in the leaf shape regulation is presented based on our observations.

Identification and analysis of Arabidopsis expressed sequence tags characteristic of non-coding RNAs

Sequencing of the Arabidopsis genome has led to the identification of thousands of new putative genes based on the predicted proteins they encode. Genes encoding tRNAs, ribosomal RNAs, and small nucleolar RNAs have also been annotated; however, a potentially important class of genes has largely escaped previous annotation efforts. These genes correspond to RNAs that lack significant open reading frames and encode RNA as their final product. Accumulating evidence indicates that such "non-coding RNAs" (ncRNAs) can play critical roles in a wide range of cellular processes, including chromosomal silencing, transcriptional regulation, developmental control, and responses to stress. Approximately 15 putative Arabidopsis ncRNAs have been reported in the literature or have been annotated. Although several have homologs in other plant species, all appear to be plant specific, with the exception of signal recognition particle RNA. Conversely, none of the ncRNAs reported from yeast or animal systems have homologs in Arabidopsis or other plants. To identify additional genes that are likely to encode ncRNAs, we used computational tools to filter protein-coding genes from genes corresponding to 20,000 expressed sequence tag clones. Using this strategy, we identified 19 clones with characteristics of ncRNAs, nine putative peptide-coding RNAs with open reading frames smaller than 100 amino acids, and 11 that could not be differentiated between the two categories. Again, none of these clones had homologs outside the plant kingdom, suggesting that most Arabidopsis ncRNAs are likely plant specific. These data indicate that ncRNAs represent a significant and underdeveloped aspect of Arabidopsis genomics that deserves further study.

ZFWD: a novel subfamily of plant proteins containing a C3H zinc finger and seven WD40 repeats

We describe a new subfamily of WD repeat proteins characterised by the presence of a C3H zinc finger at the N-terminal part of the protein associated with seven WD40 repeats. We have identified four members of this subfamily in Arabidopsis thaliana, one of them with associated expressed sequence tags (ESTs). We have also identified homologous ESTs in rice, cotton, maize, poplar, pine tree and the ice plant. We do not observe animal homologues, suggesting that this subfamily could be specific for plants. Our data suggest an important role for these proteins. Based on the high sequence conservation within the conserved domains, we suggest that these proteins could have a regulatory function.

The sense of naturally transcribed antisense RNAs in plants

Naturally occurring antisense transcripts are well documented in mammals and prokaryotes but little is known about their existence and effects in plants. Generally, antisense RNAs are believed to control gene expression negatively by annealing to the complementary sequences of the sense transcript. The resulting double-stranded RNAs are thought either to affect RNA stability, transcription and/or translation directly, or to generate a signal for gene silencing and defense against viruses.

Isolation and analysis of cDNA within a 300 kb Arabidopsis thaliana genomic region located around the 100 map unit of chromosome 1

In order to analyze the organization of genes located at the 100 map unit of chromosome 1, we screened cDNAs hybridized with approximately 300kb contiguous DNA using four P1 clones and one YAC clone. A total of 40 kinds of cDNA were isolated, and their entire sequences were determined. A comparison with the GenBank/EMBL database indicated that three of the cDNAs have been found in Arabidopsis, and that similar sequences to 18 of the cDNAs had been detected in Arabidopsis or other organisms. cDNAs were aligned on a physical map of the contiguous DNA, and the transcriptional direction of each cDNA was determined. This contiguous DNA contains a large direct repeat, which contains five genes. In addition, identical or very similar sequences to two cDNAs are located in a narrow region. Thus, a total of 50 genes were identified, and the gene density was revealed to be approximately one gene every 6kb. In addition, cDNA sequencing revealed the existence of unusual transcripts. A sequence of seven cDNAs seemed to have no significant open reading frames. Furthermore, the existence of antisense RNA and the possibility of alternative splicing were also revealed.

Classification of Arabidopsis thaliana gene sequences: clustering of coding sequences into two groups according to codon usage improves gene prediction

While genomic sequences are accumulating, finding the location of the genes remains a major issue that can be solved only for about a half of them by homology searches. Prediction methods are thus required, but unfortunately are not fully satisfying. Most prediction methods implicitly assume a unique model for genes. This is an oversimplification as demonstrated by the possibility to group coding sequences into several classes in Escherichia coli and other genomes. As no classification existed for Arabidopsis thaliana, we classified genes according to the statistical features of their coding sequences. A clustering algorithm using a codon usage model was developed and applied to coding sequences from A. thaliana, E. coli, and a mixture of both. By using it, Arabidopsis sequences were clustered into two classes. The CU1 and CU2 classes differed essentially by the choice of pyrimidine bases at the codon silent sites: CU2 genes often use C whereas CU1 genes prefer T. This classification discriminated the Arabidopsis genes according to their expressiveness, highly expressed genes being clustered in CU2 and genes expected to have a lower expression, such as the regulatory genes, in CU1. The algorithm separated the sequences of the Escherichia-Arabidopsis mixed data set into five classes according to the species, except for one class. This mixed class contained 89 % Arabidopsis genes from CU1 and 11 % E. coli genes, mostly horizontally transferred. Interestingly, most genes encoding organelle-targeted proteins, except the photosynthetic and photoassimilatory ones, were clustered in CU1. By tailoring the GeneMark CDS prediction algorithm to the observed coding sequence classes, its quality of prediction was greatly improved. Similar improvement can be expected with other prediction systems.