A large number of novel coding small open reading frames in the intergenic regions of the Arabidopsis thaliana genome are transcribed and/or under purifying selection

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).

sORF finder: a program package to identify small open reading frames with high coding potential

SUMMARY

sORF finder is a program package for identifying small open reading frames (sORFs) with high-coding potential. This application allows the identification of coding sORFs according to the nucleotide composition bias among coding sequences and the potential functional constraint at the amino acid level through evaluation of synonymous and non-synonymous substitution rates.

AVAILABILITY

Online tools and source codes are freely available at http://evolver.psc.riken.jp/.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

PLAZA: a comparative genomics resource to study gene and genome evolution in plants

The number of sequenced genomes of representatives within the green lineage is rapidly increasing. Consequently, comparative sequence analysis has significantly altered our view on the complexity of genome organization, gene function, and regulatory pathways. To explore all this genome information, a centralized infrastructure is required where all data generated by different sequencing initiatives is integrated and combined with advanced methods for data mining. Here, we describe PLAZA, an online platform for plant comparative genomics (http://bioinformatics.psb.ugent.be/plaza/). This resource integrates structural and functional annotation of published plant genomes together with a large set of interactive tools to study gene function and gene and genome evolution. Precomputed data sets cover homologous gene families, multiple sequence alignments, phylogenetic trees, intraspecies whole-genome dot plots, and genomic colinearity between species. Through the integration of high confidence Gene Ontology annotations and tree-based orthology between related species, thousands of genes lacking any functional description are functionally annotated. Advanced query systems, as well as multiple interactive visualization tools, are available through a user-friendly and intuitive Web interface. In addition, detailed documentation and tutorials introduce the different tools, while the workbench provides an efficient means to analyze user-defined gene sets through PLAZA's interface. In conclusion, PLAZA provides a comprehensible and up-to-date research environment to aid researchers in the exploration of genome information within the green plant lineage.

An archived activation tagged population of Arabidopsis thaliana to facilitate forward genetics approaches

BACKGROUND

Functional genomics tools provide researchers with the ability to apply high-throughput techniques to determine the function and interaction of a diverse range of genes. Mutagenized plant populations are one such resource that facilitate gene characterisation. They allow complex physiological responses to be correlated with the expression of single genes in planta, through either reverse genetics where target genes are mutagenized to assay the affect, or through forward genetics where populations of mutant lines are screened to identify those whose phenotype diverges from wild type for a particular trait. One limitation of these types of populations is the prevalence of gene redundancy within plant genomes, which can mask the affect of individual genes. Activation or enhancer populations, which not only provide knock-out but also dominant activation mutations, can facilitate the study of such genes.

RESULTS

We have developed a population of almost 50,000 activation tagged A. thaliana lines that have been archived as individual lines to the T3 generation. The population is an excellent tool for both reverse and forward genetic screens and has been used successfully to identify a number of novel mutants. Insertion site sequences have been generated and mapped for 15,507 lines to enable further application of the population, while providing a clear distribution of T-DNA insertions across the genome. The population is being screened for a number of biochemical and developmental phenotypes, provisional data identifying novel alleles and genes controlling steps in proanthocyanidin biosynthesis and trichome development is presented.

CONCLUSION

This publicly available population provides an additional tool for plant researcher's to assist with determining gene function for the many as yet uncharacterised genes annotated within the Arabidopsis genome sequence http://aafc-aac.usask.ca/FST. The presence of enhancer elements on the inserted T-DNA molecule allows both knock-out and dominant activation phenotypes to be identified for traits of interest.

Plant peptides in signalling: looking for new partners

A novel candidate ligand-receptor system, INFLORESCENCE DEFICIENT IN ABSCISSION (IDA) and the related receptor-like kinases (RLKs) HAESA (HAE) and HAESA-LIKE (HSL)2, has been shown to control floral abscission in Arabidopsis thaliana. Furthermore, several IDA-LIKE (IDL) proteins, which contain a conserved C-terminal domain resembling that of the CLAVATA (CLV)3-ENDOSPERM SURROUNDING REGION (ESR)-RELATED (CLE) protein family, have been shown to be partially redundant with IDA. Here, we use the genetic similarities between the IDA and CLV3 signalling systems to hypothesize that closely related peptide ligands are likely to interact with families of closely related RLKs. Guided by this hypothesis and with the aid of genetics and novel methods, ligand-receptor systems can be identified to improve our understanding of developmental processes in plants.

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.

Discovery and revision of Arabidopsis genes by proteogenomics

Gene annotation underpins genome science. Most often protein coding sequence is inferred from the genome based on transcript evidence and computational predictions. While generally correct, gene models suffer from errors in reading frame, exon border definition, and exon identification. To ascertain the error rate of Arabidopsis thaliana gene models, we isolated proteins from a sample of Arabidopsis tissues and determined the amino acid sequences of 144,079 distinct peptides by tandem mass spectrometry. The peptides corresponded to 1 or more of 3 different translations of the genome: a 6-frame translation, an exon splice-graph, and the currently annotated proteome. The majority of the peptides (126,055) resided in existing gene models (12,769 confirmed proteins), comprising 40% of annotated genes. Surprisingly, 18,024 novel peptides were found that do not correspond to annotated genes. Using the gene finding program AUGUSTUS and 5,426 novel peptides that occurred in clusters, we discovered 778 new protein-coding genes and refined the annotation of an additional 695 gene models. The remaining 13,449 novel peptides provide high quality annotation (>99% correct) for thousands of additional genes. Our observation that 18,024 of 144,079 peptides did not match current gene models suggests that 13% of the Arabidopsis proteome was incomplete due to approximately equal numbers of missing and incorrect gene models.

Unique genes in plants: specificities and conserved features throughout evolution

Background

Plant genomes contain a high proportion of duplicated genes as a result of numerous whole, segmental and local duplications. These duplications lead up to the formation of gene families, which are the usual material for many evolutionary studies. However, all characterized genomes include single-copy (unique) genes that have not received much attention. Unlike gene duplication, gene loss is not an unspecific mechanism but is rather influenced by a functional selection. In this context, we have established and used stringent criteria in order to identify suitable sets of unique genes present in plant proteomes. Comparisons of unique genes in the green phylum were used to characterize the gene and protein features exhibited by both conserved and species-specific unique genes.

Results

We identified the unique genes within both A. thaliana and O. sativa genomes and classified them according to the number of homologs in the alternative species: none (U{1:0}), one (U{1:1}) or several (U{1:m}). Regardless of the species, all the genes in these groups present some conserved characteristics, such as small average protein size and abnormal intron number. In order to understand the origin and function of unique genes, we further characterized the U{1:1} gene pairs. The possible involvement of sequence convergence in the creation of U{1:1} pairs was discarded due to the frequent conservation of intron positions. Furthermore, an orthology relationship between the two members of each U{1:1} pair was strongly supported by a high conservation in the protein sizes and transcription levels. Within the promoter of the unique conserved genes, we found a number of TATA and TELO boxes that specifically differed from their mean number in the whole genome. Many unique genes have been conserved as unique through evolution from the green alga Ostreococcus lucimarinus to higher plants. Plant unique genes may also have homologs in bacteria and we showed a link between the targeting towards plastids of proteins encoded by plant nuclear unique genes and their homology with a bacterial protein.

Conclusion

Many of the A. thaliana and O. sativa unique genes are conserved in plants for which the ancestor diverged at least 725 million years ago (MYA). Half of these genes are also present in other eukaryotic and/or prokaryotic species. Thus, our results indicate that (i) a strong negative selection pressure has conserved a number of genes as unique in genomes throughout evolution, (ii) most unique genes are subjected to a low divergence rate, (iii) they have some features observed in housekeeping genes but for most of them there is no functional annotation and (iv) they may have an ancient origin involving a possible gene transfer from ancestral chloroplasts or bacteria to the plant nucleus.

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).

The EPIP peptide of INFLORESCENCE DEFICIENT IN ABSCISSION is sufficient to induce abscission in arabidopsis through the receptor-like kinases HAESA and HAESA-LIKE2

In Arabidopsis thaliana, the final step of floral organ abscission is regulated by INFLORESCENCE DEFICIENT IN ABSCISSION (IDA): ida mutants fail to abscise floral organs, and plants overexpressing IDA display earlier abscission. We show that five IDA-LIKE (IDL) genes are expressed in different tissues, but plants overexpressing these genes have phenotypes similar to IDA-overexpressing plants, suggesting functional redundancy. IDA/IDL proteins have N-terminal signal peptides and a C-terminal conserved motif (extended PIP [EPIP]) at the C terminus (EPIP-C). IDA can, similar to CLAVATA3, be processed by an activity from cauliflower meristems. The EPIP-C of IDA and IDL1 replaced IDA function in vivo, when the signal peptide was present. In addition, synthetic IDA and IDL1 EPIP peptides rescued ida and induced early floral abscission in wild-type flowers. The EPIP-C of the other IDL proteins could partially substitute for IDA function. Similarly to ida, a double mutant between the receptor-like kinases (RLKs) HAESA (HAE) and HAESA-LIKE2 (HSL2) displays nonabscising flowers. Neither overexpression of IDA nor synthetic EPIP or EPIP-C peptides could rescue the hae hsl2 abscission deficiency. We propose that IDA and the IDL proteins constitute a family of putative ligands that act through RLKs to regulate different events during plant development.

At-TAX: a whole genome tiling array resource for developmental expression analysis and transcript identification in Arabidopsis thaliana

Gene expression maps for model organisms, including Arabidopsis thaliana, have typically been created using gene-centric expression arrays. Here, we describe a comprehensive expression atlas, Arabidopsis thaliana Tiling Array Express (At-TAX), which is based on whole-genome tiling arrays. We demonstrate that tiling arrays are accurate tools for gene expression analysis and identified more than 1,000 unannotated transcribed regions. Visualizations of gene expression estimates, transcribed regions, and tiling probe measurements are accessible online at the At-TAX homepage.

Lilliputians get into the limelight: novel class of small peptide genes in morphogenesis

Generally, bioactive small peptides are derived from precursors with signal sequences at their N-terminal ends, which undergo modification and proteolysis through a secretory pathway. By contrast, small peptides encoded in short open reading frames (sORF) lack signaling sequences and therefore are released into the cytoplasm, which may result in their having functions distinct from those of secreted peptides. Several small peptides encoded by sORF are involved in the morphogenesis of multicellular organisms. POLARIS, ROTUNDIFOLIA4, and Enod40 are plant peptides that are involved, respectively, in root formation, leaf shape control, and cortical cell division during nodule formation. Brick1/HSPC300 is an evolutionarily conserved component of the actin reorganization complex. polished rice/tarsal-less and mille-pattes encode related small peptides that are required for epithelial morphogenesis in Drosophila and segmentation in Tribolium. There are only a few known examples of small peptides encoded by sORF, and their molecular functions are still largely obscure. Nevertheless, an increasing number of sORF genes is being identified, and further research should reveal their roles in novel molecular mechanisms underlying developmental events.

Silencing a candidate nematode effector gene corresponding to the tomato resistance gene Mi-1 leads to acquisition of virulence

The Mi-1 gene in tomato confers effective resistance against several species of root-knot nematode, including Meloidogyne javanica. A strain of M. javanica that can reproduce on tomato with Mi-1 was obtained from a culture of an avirulent strain after greenhouse selection. DNA blots and amplified fragment length polymorphism (AFLP) analysis indicated that the two nematode strains are closely related. Expression patterns visualized as cDNA AFLPs were nearly identical except for a cDNA fragment, Cg-1, that was present in the avirulent strain but not in the virulent strain. DNA blots showed that Cg-1 corresponds to a member of a small gene family with one or more copies missing in the virulent strain compared with the avirulent strain. Except for the presence of a histone stem loop near the 3' end of the transcript, Cg-1 shows no similarity to other sequences in GenBank. The longest open reading frame is 32 amino acids and initiates at the fourth AUG in the predicted transcript. When nematode juveniles of the Mi-1-avirulent strain were soaked in dsRNA corresponding to part of the predicted Cg-1 transcript, they produced progeny that were virulent on tomato carrying the Mi-1 gene, strongly suggesting that Cg-1 is required in the nematode for Mi-1-mediated resistance.

The draft genome of the transgenic tropical fruit tree papaya (Carica papaya Linnaeus)

Papaya, a fruit crop cultivated in tropical and subtropical regions, is known for its nutritional benefits and medicinal applications. Here we report a 3x draft genome sequence of 'SunUp' papaya, the first commercial virus-resistant transgenic fruit tree to be sequenced. The papaya genome is three times the size of the Arabidopsis genome, but contains fewer genes, including significantly fewer disease-resistance gene analogues. Comparison of the five sequenced genomes suggests a minimal angiosperm gene set of 13,311. A lack of recent genome duplication, atypical of other angiosperm genomes sequenced so far, may account for the smaller papaya gene number in most functional groups. Nonetheless, striking amplifications in gene number within particular functional groups suggest roles in the evolution of tree-like habit, deposition and remobilization of starch reserves, attraction of seed dispersal agents, and adaptation to tropical daylengths. Transgenesis at three locations is closely associated with chloroplast insertions into the nuclear genome, and with topoisomerase I recognition sites. Papaya offers numerous advantages as a system for fruit-tree functional genomics, and this draft genome sequence provides the foundation for revealing the basis of Carica's distinguishing morpho-physiological, medicinal and nutritional properties.

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.

The Arabidopsis Information Resource (TAIR): gene structure and function annotation

The Arabidopsis Information Resource (TAIR, http://arabidopsis.org) is the model organism database for the fully sequenced and intensively studied model plant Arabidopsis thaliana. Data in TAIR is derived in large part from manual curation of the Arabidopsis research literature and direct submissions from the research community. New developments at TAIR include the addition of the GBrowse genome viewer to the TAIR site, a redesigned home page, navigation structure and portal pages to make the site more intuitive and easier to use, the launch of several TAIR web services and a new genome annotation release (TAIR7) in April 2007. A combination of manual and computational methods were used to generate this release, which contains 27,029 protein-coding genes, 3889 pseudogenes or transposable elements and 1123 ncRNAs (32,041 genes in all, 37,019 gene models). A total of 681 new genes and 1002 new splice variants were added. Overall, 10,098 loci (one-third of all loci from the previous TAIR6 release) were updated for the TAIR7 release.

The Arabidopsis Information Resource (TAIR): gene structure and function annotation

The Arabidopsis Information Resource (TAIR, http://arabidopsis.org) is the model organism database for the fully sequenced and intensively studied model plant Arabidopsis thaliana. Data in TAIR is derived in large part from manual curation of the Arabidopsis research literature and direct submissions from the research community. New developments at TAIR include the addition of the GBrowse genome viewer to the TAIR site, a redesigned home page, navigation structure and portal pages to make the site more intuitive and easier to use, the launch of several TAIR web services and a new genome annotation release (TAIR7) in April 2007. A combination of manual and computational methods were used to generate this release, which contains 27 029 protein-coding genes, 3889 pseudogenes or transposable elements and 1123 ncRNAs (32 041 genes in all, 37 019 gene models). A total of 681 new genes and 1002 new splice variants were added. Overall, 10 098 loci (one-third of all loci from the previous TAIR6 release) were updated for the TAIR7 release.

Genome-wide expression profiling of the Arabidopsis female gametophyte identifies families of small, secreted proteins

The female gametophyte of flowering plants, the embryo sac, develops within the diploid (sporophytic) tissue of the ovule. While embryo sac–expressed genes are known to be required at multiple stages of the fertilization process, the set of embryo sac–expressed genes has remained poorly defined. In particular, the set of genes responsible for mediating intracellular communication between the embryo sac and the male gametophyte, the pollen grain, is unknown. We used high-throughput cDNA sequencing and whole-genome tiling arrays to compare gene expression in wild-type ovules to that in dif1 ovules, which entirely lack embryo sacs, and myb98 ovules, which are impaired in pollen tube attraction. We identified nearly 400 genes that are downregulated in dif1 ovules. Seventy-eight percent of these embryo sac–dependent genes were predicted to encode for secreted proteins, and 60% belonged to multigenic families. Our results define a large number of candidate extracellular signaling molecules that may act during embryo sac development or fertilization; less than half of these are represented on the widely used ATH1 expression array. In particular, we found that 37 out of 40 genes encoding Domain of Unknown Function 784 (DUF784) domains require the synergid-specific transcription factor MYB98 for expression. Several DUF784 genes were transcribed in synergid cells of the embryo sac, implicating the DUF784 gene family in mediating late stages of embryo sac development or interactions with pollen tubes. The coexpression of highly similar proteins suggests a high degree of functional redundancy among embryo sac genes.

During the sexual reproduction of flowering plants, a pollen tube delivers sperm cells to a specialized group of cells known as the embryo sac, which contains the egg cell. It is known that embryo sacs are active participants in guiding the growth of pollen tubes, in facilitating fertilization, and in initiating seed development. However, the genes responsible for the complex biology of embryo sacs are poorly understood. The authors use two recently developed technologies, whole-genome tiling microarrays and high-throughput cDNA sequencing, to identify hundreds of genes expressed in embryo sacs of Arabidopsis thaliana. Most embryo sac–dependent genes have no known function, and include entire families of related genes that are only expressed in embryo sacs. Furthermore, most embryo sac–dependent genes encode small proteins that are potentially secreted from their cells of origin, suggesting that they may act as intracellular signals or to modify the extracellular matrix during fertilization or embryo sac development. These results illustrate the extent to which our understanding of plant sexual reproduction is limited and identifies hundreds of candidate genes for future studies investigating the molecular biology of the embryo sac.