Full-length messenger RNA sequences greatly improve genome annotation

Dual functionality of a plant U-rich intronic sequence element

In potato invertase genes, the constitutively included, 9-nucleotide (nt)-long mini-exon requires a strong branchpoint and U-rich polypyrimidine tract for inclusion. The strength of these splicing signals was demonstrated by greatly enhanced splicing of a poorly spliced intron and by their ability to support splicing of an artificial mini-exon, following their introduction. Plant introns also require a second splicing signal, UA-rich intronic elements, for efficient intron splicing. Mutation of the branchpoint caused loss of mini-exon inclusion without loss of splicing enhancement, showing that the same U-rich sequence can function as either a polypyrimidine tract or a UA-rich intronic element. The distinction between the splicing signals depended on intron context (the presence or absence of an upstream, adjacent and functional branchpoint), and on the sequence context of the U-rich elements. Polypyrimidine tracts tolerated C residues while UA-rich intronic elements tolerated As. Thus, in plant introns, U-rich splicing elements can have dual roles as either a general plant U-rich splicing signal or a polypyrimidine tract. Finally, overexpression of two different U-rich binding proteins enhanced intron recognition significantly. These results highlight the importance of co-operation between splicing signals, the importance of other nucleotides within U-rich elements for optimal binding of competing splicing factors and effects on splicing efficiency of U-rich binding proteins.

Effect of manganese toxicity on the proteome of the leaf apoplast in cowpea

Excess manganese (Mn) supply causes formation of visible brown depositions in the cell walls of leaves of cowpea (Vigna unguiculata), which consist of oxidized Mn and oxidized phenols. Because oxidation of Mn and phenolic compounds in the leaf apoplast was proposed to be catalyzed by apoplastic peroxidases (PODs), induction of these enzymes by Mn excess was investigated. POD activity increased upon prolonged Mn treatment in the leaf tissue. Simultaneously, a significant increase in the concentration of soluble apoplastic proteins in "apoplastic washing fluid" was observed. The identity of the released proteins was systematically characterized by analysis of the apoplast proteome using two-dimensional gel electrophoresis and liquid chromatography-tandem mass spectrometry. Some of the identified proteins exhibit sequence identity to acidic PODs from other plants. Several other proteins show homologies to pathogenesis-related proteins, e.g. glucanase, chitinase, and thaumatin-like proteins. Because pathogenesis-related-like proteins are known to be induced by various other abiotic and biotic stresses, a specific physiological role of these proteins in response to excess Mn supply remains to be established. The specific role of apoplastic PODs in the response of plants to Mn stress is discussed.

Functional genomics of wood quality and properties

Genomics promises to enrich the investigations of biology and biochemistry. Current advancements in genomics have major implications for genetic improvement in animals, plants, and microorganisms, and for our understanding of cell growth, development, differentiation, and communication. Significant progress has been made in the understanding of plant genomics in recent years, and the area continues to progress rapidly. Functional genomics offers enormous potential to tree improvement and the understanding of gene expression in this area of science worldwide. In this review we focus on functional genomics of wood quality and properties in trees, mainly based on progresses made in genomics study of Pinus and Populus. The aims of this review are to summarize the current status of functional genomics including: (1) Gene discovery; (2) EST and genomic sequencing; (3) From EST to functional genomics; (4) Approaches to functional analysis; (5) Engineering lignin biosynthesis; (6) Modification of cell wall biogenesis; and (7) Molecular modelling. Functional genomics has been greatly invested worldwide and will be important in identifying candidate genes whose function is critical to all aspects of plant growth, development, differentiation, and defense. Forest biotechnology industry will significantly benefit from the advent of functional genomics of wood quality and properties.

Improving the Arabidopsis genome annotation using maximal transcript alignment assemblies

The spliced alignment of expressed sequence data to genomic sequence has proven a key tool in the comprehensive annotation of genes in eukaryotic genomes. A novel algorithm was developed to assemble clusters of overlapping transcript alignments (ESTs and full-length cDNAs) into maximal alignment assemblies, thereby comprehensively incorporating all available transcript data and capturing subtle splicing variations. Complete and partial gene structures identified by this method were used to improve The Institute for Genomic Research Arabidopsis genome annotation (TIGR release v.4.0). The alignment assemblies permitted the automated modeling of several novel genes and >1000 alternative splicing variations as well as updates (including UTR annotations) to nearly half of the approximately 27 000 annotated protein coding genes. The algorithm of the Program to Assemble Spliced Alignments (PASA) tool is described, as well as the results of automated updates to Arabidopsis gene annotations.

A RING-H2 zinc-finger protein gene RIE1 is essential for seed development in Arabidopsis

RING zinc-finger proteins play important roles in the regulation of development in a variety of organisms. In the plant kingdom, few genes encoding RING zinc-finger proteins have been documented with visible effects on plant growth and development. A novel gene, RIE1, encoding a RING-H2 zinc-finger protein was identified in Arabidopsis thaliana and is characterized in this paper. RIE1 encodes a predicted protein product of 359 amino acids residues with a molecular mass of 40 kDa, with a RING-H2 zinc-finger motif located at the extreme end of the C-terminus. Characterization of a Dissociation (Ds) insertion line (SGT4559) and a T-DNA insertion line (SRIE1) demonstrated that disruption of RIE1 is embryo-lethal. SGT4559 heterozygous plants produced seeds with embryo development arrested from globular to torpedo stages. Some mutant seeds were rescued by embryo culture, and the mutant (rie1) plants seemed to grow normally compared to wild-type plants, except that the mutants produced only abnormal seeds. However, RIE1 was expressed in different tissues throughout the whole plant as revealed by northern blot analysis and gene fusion assay of RIE1 promoter with the beta-glucuronidase (GUS) gene. Our results indicated that RIE1 plays an essential role in seed development.

GlimmerM, Exonomy and Unveil: three ab initio eukaryotic genefinders

We present three programs for ab initio gene prediction in eukaryotes: Exonomy, Unveil and GlimmerM. Exonomy is a 23-state Generalized Hidden Markov Model (GHMM), Unveil is a 283-state standard Hidden Markov Model (HMM) and GlimmerM is a previously-described genefinder which utilizes decision trees and Interpolated Markov Models (IMMs). All three are readily re-trainable for new organisms and have been found to perform well compared to other genefinders. Results are presented for Arabidopsis thaliana. Cases have been found where each of the genefinders outperforms each of the others, demonstrating the collective value of this ensemble of genefinders. These programs are all accessible through webservers at http://www.tigr.org/software.

GeneSeqer@PlantGDB: Gene structure prediction in plant genomes

The GeneSeqer@PlantGDB Web server (http://www.plantgdb.org/cgi-bin/GeneSeqer.cgi) provides a gene structure prediction tool tailored for applications to plant genomic sequences. Predictions are based on spliced alignment with source-native ESTs and full-length cDNAs or non-native probes derived from putative homologous genes. The tool is illustrated with applications to refinement of current gene structure annotation and de novo annotation of draft genomic sequences. The service should facilitate expert annotation as a community effort by providing convenient access to all public plant sequences via the PlantGDB database, a simple four-step protocol for spliced alignment and visually appealing displays of the predicted gene structures in addition to detailed sequence alignments.

Computational discovery of internal micro-exons

Very short exons, also known as micro-exons, occur in large numbers in some eukaryotic genomes. Existing annotation tools have a limited ability to recognize these short sequences, which range in length up to 25 bp. Here, we describe a computational method for the identification of micro-exons using near-perfect alignments between cDNA and genomic DNA sequences. Using this method, we detected 319 micro-exons in 4 complete genomes, of which 224 were previously unknown, human (170), the nematode Caenorhabditis elegans (4), the fruit fly Drosophila melanogaster (14), and the mustard plant Arabidopsis thaliana (36). Comparison of our computational method with popular cDNA alignment programs shows that the new algorithm is both efficient and accurate. The algorithm also aids in the discovery of micro-exon-skipping events and cross-species micro-exon conservation.

Refined annotation of the Arabidopsis genome by complete expressed sequence tag mapping

Expressed sequence tags (ESTs) currently encompass more entries in the public databases than any other form of sequence data. Thus, EST data sets provide a vast resource for gene identification and expression profiling. We have mapped the complete set of 176,915 publicly available Arabidopsis EST sequences onto the Arabidopsis genome using GeneSeqer, a spliced alignment program incorporating sequence similarity and splice site scoring. About 96% of the available ESTs could be properly aligned with a genomic locus, with the remaining ESTs deriving from organelle genomes and non-Arabidopsis sources or displaying insufficient sequence quality for alignment. The mapping provides verified sets of EST clusters for evaluation of EST clustering programs. Analysis of the spliced alignments suggests corrections to current gene structure annotation and provides examples of alternative and non-canonical pre-mRNA splicing. All results of this study were parsed into a database and are accessible via a flexible Web interface at http://www.plantgdb.org/AtGDB/.

Comparison of rice and Arabidopsis annotation

Several versions of the rice genome were published in 2002, providing a first overview of the genome content of this model monocot. At the same time, the genome of the model dicot, Arabidopsis thaliana, reached a new level of annotation as thousands of full-length cDNA sequences were integrated with the genome sequence.

Genome-wide analysis of NBS-LRR-encoding genes in Arabidopsis

The Arabidopsis genome contains approximately 200 genes that encode proteins with similarity to the nucleotide binding site and other domains characteristic of plant resistance proteins. Through a reiterative process of sequence analysis and reannotation, we identified 149 NBS-LRR-encoding genes in the Arabidopsis (ecotype Columbia) genomic sequence. Fifty-six of these genes were corrected from earlier annotations. At least 12 are predicted to be pseudogenes. As described previously, two distinct groups of sequences were identified: those that encoded an N-terminal domain with Toll/Interleukin-1 Receptor homology (TIR-NBS-LRR, or TNL), and those that encoded an N-terminal coiled-coil motif (CC-NBS-LRR, or CNL). The encoded proteins are distinct from the 58 predicted adapter proteins in the previously described TIR-X, TIR-NBS, and CC-NBS groups. Classification based on protein domains, intron positions, sequence conservation, and genome distribution defined four subgroups of CNL proteins, eight subgroups of TNL proteins, and a pair of divergent NL proteins that lack a defined N-terminal motif. CNL proteins generally were encoded in single exons, although two subclasses were identified that contained introns in unique positions. TNL proteins were encoded in modular exons, with conserved intron positions separating distinct protein domains. Conserved motifs were identified in the LRRs of both CNL and TNL proteins. In contrast to CNL proteins, TNL proteins contained large and variable C-terminal domains. The extant distribution and diversity of the NBS-LRR sequences has been generated by extensive duplication and ectopic rearrangements that involved segmental duplications as well as microscale events. The observed diversity of these NBS-LRR proteins indicates the variety of recognition molecules available in an individual genotype to detect diverse biotic challenges.

Auxin signaling in Arabidopsis leaf vascular development

A number of observations have implicated auxin in the formation of vascular tissues in plant organs. These include vascular strand formation in response to local auxin application, the effects of impaired auxin transport on vascular patterns and suggestive phenotypes of Arabidopsis auxin response mutants. In this study, we have used molecular markers to visualize auxin response patterns in developing Arabidopsis leaves as well as Arabidopsis mutants and transgenic plants to trace pathways of auxin signal transduction controlling the expression of early procambial genes. We show that in young Arabidopsis leaf primordia, molecular auxin response patterns presage sites of procambial differentiation. This is the case not only in normal development but also upon experimental manipulation of auxin transport suggesting that local auxin signals are instrumental in patterning Arabidopsis leaf vasculature. We further found that the activity of the Arabidopsis gene MONOPTEROS, which is required for proper vascular differentiation, is also essential in a spectrum of auxin responses, which include the regulation of rapidly auxin-inducible AUX/IAA genes, and discovered the tissue-specific vascular expression profile of the class I homeodomain-leucine zipper gene, AtHB20. Interestingly, MONOPTEROS activity is a limiting factor in the expression of AtHB8 and AtHB20, two genes encoding transcriptional regulators expressed early in procambial development. Our observations connect general auxin signaling with early controls of vascular differentiation and suggest molecular mechanisms for auxin signaling in patterned cell differentiation.

Powdery mildew-induced Mla mRNAs are alternatively spliced and contain multiple upstream open reading frames

In barley (Hordeum vulgare), the Mla13 powdery mildew resistance gene confers Rar1-dependent, AvrMla13-specific resistance to Blumeria graminis f. sp. hordei (Bgh). We have identified cDNA and genomic copies of Mla13 and used this coiled-coil nucleotide-binding site leucine-rich repeat protein-encoding gene as a model for the regulation of host resistance to obligate biotrophic fungi in cereals. We demonstrate quantitatively that a rapid increase in the accumulation of Mla transcripts and transcripts of the Mla-signaling genes, Rar1 and Sgt1, is triggered between 16 and 20 h post inoculation, the same time frame that haustoria of avirulent Bgh make contact with the host cell plasma membrane. An abundance of Mla13 cDNAs revealed five classes of transcript leader regions containing two alternatively spliced introns and up to three upstream open reading frames (uORFs). Alternative splicing of introns in the transcript leader region results in a different number of uORFs and variability in the size of uORF2. These results indicate that regulation of Mla transcript accumulation is not constitutive and that induction is coordinately controlled by recognition-specific factors. The sudden increase in specific transcript levels could account for the rapid defense response phenotype conferred by Mla6 and Mla13.

Update on chloroplast proteomics

Currently, relatively few proteomics studies of chloroplast have been published, but the field has just started emerging and is likely to develop more rapidly in the future. While the complex membrane structure of the chloroplast makes it difficult to study its entire proteome by global approaches, proteomics has considerably increased our knowledge of the proteins of single compartments such as, for instance, the envelope and the thylakoid lumen. Proteomics has also succeeded in the subunit characterisation of select protein complexes such as the ribosomes and the cytochrome b (6)f complex. In addition, proteomics was successfully applied to find new potential target pathways for thioredoxin-mediated signal transduction. In this review, we present an overview of the latest developments in the field of chloroplast proteomics and discuss their impact on photosynthesis research. In addition, we summarise the current state of research in proteomics of the photosynthetic cyanobactrium Synechocystis sp. PCC 6803.

CATMA: a complete Arabidopsis GST database

The Complete Arabidopsis Transcriptome Micro Array (CATMA) database contains gene sequence tag (GST) and gene model sequences for over 70% of the predicted genes in the Arabidopsis thaliana genome as well as primer sequences for GST amplification and a wide range of supplementary information. All CATMA GST sequences are specific to the gene for which they were designed, and all gene models were predicted from a complete reannotation of the genome using uniform parameters. The database is searchable by sequence name, sequence homology or direct SQL query, and is available through the CATMA website at http://www.catma.org/.

Cloning and sequencing of cDNAs for hypothetical genes from chromosome 2 of Arabidopsis

About 25% of the genes in the fully sequenced and annotated Arabidopsis genome have structures that are predicted solely by computer algorithms with no support from either nucleic acid or protein homologs from other species or expressed sequence matches from Arabidopsis. These are referred to as "hypothetical genes." On chromosome 2, sequenced by The Institute for Genomic Research, there are approximately 800 hypothetical genes among a total of approximately 4,100 genes. To test their expression under various growth conditions and in specific tissues, we used six cDNA populations prepared from cold-treated, heat-treated, and pathogen (Xanthomonas campestris pv campestris)-infected plants, callus, roots, and young seedlings. To date, 169 hypothetical genes were tested, and 138 of them are found to be expressed in one or more of the six cDNA populations. By sequencing multiple clones from each 5'- and 3'-rapid amplification of cDNA ends (RACE) product and assembling the sequences, we generated full-length sequences for 16 of these genes. For 14 genes, there was one full-length assembly that precisely supported the intron-exon boundaries of their gene predictions, adding only 5'- and 3'-untranslated region sequences. However, for three of these genes, the other assemblies represent additional exons and alternatively spliced or unspliced introns. For the remaining two genes, the cDNA sequences reveal major differences with predicted gene structures. In addition, a total of six genes displayed more than one polyadenylation site. These data will be used to update gene models in The Institute for Genomic Research annotation database ATH1.

Assessing the impact of comparative genomic sequence data on the functional annotation of the Drosophila genome

BACKGROUND

It is widely accepted that comparative sequence data can aid the functional annotation of genome sequences; however, the most informative species and features of genome evolution for comparison remain to be determined.

RESULTS

We analyzed conservation in eight genomic regions (apterous, even-skipped, fushi tarazu, twist, and Rhodopsins 1, 2, 3 and 4) from four Drosophila species (D. erecta, D. pseudoobscura, D. willistoni, and D. littoralis) covering more than 500 kb of the D. melanogaster genome. All D. melanogaster genes (and 78-82% of coding exons) identified in divergent species such as D. pseudoobscura show evidence of functional constraint. Addition of a third species can reveal functional constraint in otherwise non-significant pairwise exon comparisons. Microsynteny is largely conserved, with rearrangement breakpoints, novel transposable element insertions, and gene transpositions occurring in similar numbers. Rates of amino-acid substitution are higher in uncharacterized genes relative to genes that have previously been studied. Conserved non-coding sequences (CNCSs) tend to be spatially clustered with conserved spacing between CNCSs, and clusters of CNCSs can be used to predict enhancer sequences.

CONCLUSIONS

Our results provide the basis for choosing species whose genome sequences would be most useful in aiding the functional annotation of coding and cis-regulatory sequences in Drosophila. Furthermore, this work shows how decoding the spatial organization of conserved sequences, such as the clustering of CNCSs, can complement efforts to annotate eukaryotic genomes on the basis of sequence conservation alone.

An integrated computational pipeline and database to support whole-genome sequence annotation

We describe here our experience in annotating the Drosophila melanogaster genome sequence, in the course of which we developed several new open-source software tools and a database schema to support large-scale genome annotation. We have developed these into an integrated and reusable software system for whole-genome annotation. The key contributions to overall annotation quality are the marshalling of high-quality sequences for alignments and the design of a system with an adaptable and expandable flexible architecture.

Apollo: a sequence annotation editor

The well-established inaccuracy of purely computational methods for annotating genome sequences necessitates an interactive tool to allow biological experts to refine these approximations by viewing and independently evaluating the data supporting each annotation. Apollo was developed to meet this need, enabling curators to inspect genome annotations closely and edit them. FlyBase biologists successfully used Apollo to annotate the Drosophila melanogaster genome and it is increasingly being used as a starting point for the development of customized annotation editing tools for other genome projects.

Computational analysis of core promoters in the Drosophila genome

BACKGROUND

The core promoter, a region of about 100 base-pairs flanking the transcription start site (TSS), serves as the recognition site for the basal transcription apparatus. Drosophila TSSs have generally been mapped by individual experiments; the low number of accurately mapped TSSs has limited analysis of promoter sequence motifs and the training of computational prediction tools.

RESULTS

We identified TSS candidates for about 2,000 Drosophila genes by aligning 5' expressed sequence tags (ESTs) from cap-trapped cDNA libraries to the genome, while applying stringent criteria concerning coverage and 5'-end distribution. Examination of the sequences flanking these TSSs revealed the presence of well-known core promoter motifs such as the TATA box, the initiator and the downstream promoter element (DPE). We also define, and assess the distribution of, several new motifs prevalent in core promoters, including what appears to be a variant DPE motif. Among the prevalent motifs is the DNA-replication-related element DRE, recently shown to be part of the recognition site for the TBP-related factor TRF2. Our TSS set was then used to retrain the computational promoter predictor McPromoter, allowing us to improve the recognition performance to over 50% sensitivity and 40% specificity. We compare these computational results to promoter prediction in vertebrates.

CONCLUSIONS

There are relatively few recognizable binding sites for previously known general transcription factors in Drosophila core promoters. However, we identified several new motifs enriched in promoter regions. We were also able to significantly improve the performance of computational TSS prediction in Drosophila.

Annotation of the Drosophila melanogaster euchromatic genome: a systematic review

BACKGROUND

The recent completion of the Drosophila melanogaster genomic sequence to high quality and the availability of a greatly expanded set of Drosophila cDNA sequences, aligning to 78% of the predicted euchromatic genes, afforded FlyBase the opportunity to significantly improve genomic annotations. We made the annotation process more rigorous by inspecting each gene visually, utilizing a comprehensive set of curation rules, requiring traceable evidence for each gene model, and comparing each predicted peptide to SWISS-PROT and TrEMBL sequences.

RESULTS

Although the number of predicted protein-coding genes in Drosophila remains essentially unchanged, the revised annotation significantly improves gene models, resulting in structural changes to 85% of the transcripts and 45% of the predicted proteins. We annotated transposable elements and non-protein-coding RNAs as new features, and extended the annotation of untranslated (UTR) sequences and alternative transcripts to include more than 70% and 20% of genes, respectively. Finally, cDNA sequence provided evidence for dicistronic transcripts, neighboring genes with overlapping UTRs on the same DNA sequence strand, alternatively spliced genes that encode distinct, non-overlapping peptides, and numerous nested genes.

CONCLUSIONS

Identification of so many unusual gene models not only suggests that some mechanisms for gene regulation are more prevalent than previously believed, but also underscores the complex challenges of eukaryotic gene prediction. At present, experimental data and human curation remain essential to generate high-quality genome annotations.