Comparison of rice and Arabidopsis annotation

Cis-regulatory signatures of orthologous stress-associated bZIP transcription factors from rice, sorghum and Arabidopsis based on phylogenetic footprints

Background

The potential contribution of upstream sequence variation to the unique features of orthologous genes is just beginning to be unraveled. A core subset of stress-associated bZIP transcription factors from rice (Oryza sativa) formed ten clusters of orthologous groups (COG) with genes from the monocot sorghum (Sorghum bicolor) and dicot Arabidopsis (Arabidopsis thaliana). The total cis-regulatory information content of each stress-associated COG was examined by phylogenetic footprinting to reveal ortholog-specific, lineage-specific and species-specific conservation patterns.

Results

The most apparent pattern observed was the occurrence of spatially conserved ‘core modules’ among the COGs but not among paralogs. These core modules are comprised of various combinations of two to four putative transcription factor binding site (TFBS) classes associated with either developmental or stress-related functions. Outside the core modules are specific stress (ABA, oxidative, abiotic, biotic) or organ-associated signals, which may be functioning as ‘regulatory fine-tuners’ and further define lineage-specific and species-specific cis-regulatory signatures. Orthologous monocot and dicot promoters have distinct TFBS classes involved in disease and oxidative-regulated expression, while the orthologous rice and sorghum promoters have distinct combinations of root-specific signals, a pattern that is not particularly conserved in Arabidopsis.

Conclusions

Patterns of cis-regulatory conservation imply that each ortholog has distinct signatures, further suggesting that they are potentially unique in a regulatory context despite the presumed conservation of broad biological function during speciation. Based on the observed patterns of conservation, we postulate that core modules are likely primary determinants of basal developmental programming, which may be integrated with and further elaborated by additional intrinsic or extrinsic signals in conjunction with lineage-specific or species-specific regulatory fine-tuners. This synergy may be critical for finer-scale spatio-temporal regulation, hence unique expression profiles of homologous transcription factors from different species with distinct zones of ecological adaptation such as rice, sorghum and Arabidopsis. The patterns revealed from these comparisons set the stage for further empirical validation by functional genomics.

Genome-wide computational identification of WG/GW Argonaute-binding proteins in Arabidopsis

Domains in Arabidopsis proteins NRPE1 and SPT5-like, composed almost exclusively of repeated motifs in which only WG or GW sequences and an overall amino-acid preference are conserved, have been experimentally shown to bind multiple molecules of Argonaute (AGO) protein(s). Domain swapping between the WG/GW domains of NRPE1 and the human protein GW182 showed a conserved function. As classical sequence alignment methods are poorly-adapted to detect such weakly-conserved motifs, we have developed a tool to carry out a systematic analysis to identify genes potentially encoding AGO-binding GW/WG proteins. Here, we describe exhaustive analysis of the Arabidopsis genome for all regions potentially encoding proteins bearing WG/GW motifs and consider the possible role of some of them in AGO-dependent mechanisms. We identified 20 different candidate WG/GW genes, encoding proteins in which the predicted domains range from 92aa to 654aa. These mostly correspond to a limited number of families: RNA-binding proteins, transcription factors, glycine-rich proteins, translation initiation factors and known silencing-associated proteins such as SDE3. Recent studies have argued that the interaction between WG/GW-rich domains and AGO proteins is evolutionarily conserved. Here, we demonstrate by an in silico domain-swapping simulation between plant and mammalian WG/GW proteins that the biased amino-acid composition of the AGO-binding sites is conserved.

Transformation of a dwarf Arabidopsis mutant illustrates gibberellin hormone physiology and the function of a Green Revolution gene

The introduction of dwarfing traits into crops was a major factor in increased grain yields during the "Green Revolution." In most cases those traits were the consequence of altered synthesis or response to the gibberellin (GA) plant hormones. Our current understanding of GA synthesis and physiology has been facilitated by the characterization of mutants. To introduce concepts about GA hormone physiology and plant transformation in an undergraduate laboratory course we have used ga5, a semi-dwarf Arabidopsis mutant with reduced activity of GA 20-oxidase. In this laboratory exercise, Arabidopsis ga5 mutant plants are transformed by the floral-dip method using Agrobacterium tumefaciens carrying plasmid constructs conferring kanamycin resistance and containing the GA5 gene. Within 4 weeks, seeds of transformed plants can be easily screened by antibiotic resistance on plates. After transfer to soil the dwarf mutant plants transformed with a wild-type version of the gene show normal size. In addition to offering a visual understanding of the effect of GA on stem elongation, students learn additional techniques in this experiment, including PCR and agarose gel electrophoresis. This experiment is cost effective and can be completed within a 4-month term.

Identification and characterization of lineage-specific genes within the Poaceae

Using the rice (Oryza sativa) sp. japonica genome annotation, along with genomic sequence and clustered transcript assemblies from 184 species in the plant kingdom, we have identified a set of 861 rice genes that are evolutionarily conserved among six diverse species within the Poaceae yet lack significant sequence similarity with plant species outside the Poaceae. This set of evolutionarily conserved and lineage-specific rice genes is termed conserved Poaceae-specific genes (CPSGs) to reflect the presence of significant sequence similarity across three separate Poaceae subfamilies. The vast majority of rice CPSGs (86.6%) encode proteins with no putative function or functionally characterized protein domain. For the remaining CPSGs, 8.8% encode an F-box domain-containing protein and 4.5% encode a protein with a putative function. On average, the CPSGs have fewer exons, shorter total gene length, and elevated GC content when compared with genes annotated as either transposable elements (TEs) or those genes having significant sequence similarity in a species outside the Poaceae. Multiple sequence alignments of the CPSGs with sequences from other Poaceae species show conservation across a putative domain, a novel domain, or the entire coding length of the protein. At the genome level, syntenic alignments between sorghum (Sorghum bicolor) and 103 of the 861 rice CPSGs (12.0%) could be made, demonstrating an additional level of conservation for this set of genes within the Poaceae. The extensive sequence similarity in evolutionarily distinct species within the Poaceae family and an additional screen for TE-related structural characteristics and sequence discounts these CPSGs as being misannotated TEs. Collectively, these data confirm that we have identified a specific set of genes that are highly conserved within, as well as specific to, the Poaceae.

Curated genome annotation of Oryza sativa ssp. japonica and comparative genome analysis with Arabidopsis thaliana

We present here the annotation of the complete genome of rice Oryza sativa L. ssp. japonica cultivar Nipponbare. All functional annotations for proteins and non-protein-coding RNA (npRNA) candidates were manually curated. Functions were identified or inferred in 19,969 (70%) of the proteins, and 131 possible npRNAs (including 58 antisense transcripts) were found. Almost 5000 annotated protein-coding genes were found to be disrupted in insertional mutant lines, which will accelerate future experimental validation of the annotations. The rice loci were determined by using cDNA sequences obtained from rice and other representative cereals. Our conservative estimate based on these loci and an extrapolation suggested that the gene number of rice is approximately 32,000, which is smaller than previous estimates. We conducted comparative analyses between rice and Arabidopsis thaliana and found that both genomes possessed several lineage-specific genes, which might account for the observed differences between these species, while they had similar sets of predicted functional domains among the protein sequences. A system to control translational efficiency seems to be conserved across large evolutionary distances. Moreover, the evolutionary process of protein-coding genes was examined. Our results suggest that natural selection may have played a role for duplicated genes in both species, so that duplication was suppressed or favored in a manner that depended on the function of a gene.

Functional characterization of gibberellin-regulated genes in rice using microarray system

Gibberellin (GA) is collectively referred to a group of diterpenoid acids, some of which act as plant hormones and are essential for normal plant growth and development. DNA microarray technology has become the standard tool for the parallel quantification of large numbers of messenger RNA transcripts. The power of this approach has been demonstrated in dissecting plant physiology and development, and in unraveling the underlying cellular signaling pathways. To understand the molecular mechanism by which GA regulates the growth and development of plants, with reference to the monocot model plant-rice, it is essential to identify and analyze more genes and their products at the transcription and translation levels that are regulated by GA. With the availability of draft sequences of two major rice types, indica and japonica rice, it has become possible to analyze global expression profiles of genes on a genome scale. In this review, the progress made in finding new genes in rice leaf sheath using microarray system and their characterization is discussed. It is believed that the findings made in this regard have important implications for understanding the mechanism by which GA regulates the growth and development of rice.

Decoding the rice genome

Rice cultivation is one of the most important agricultural activities on earth, with nearly 90% of it being produced in Asia. It belongs to the family of crops that includes wheat, maize and barley, and it supplies more than 50% of calories consumed by the world population. Its immense economic value and a relatively small genome size makes it a focal point for scientific investigations, so much so that four whole genome sequence drafts with varying qualities have been generated by both public and privately funded ventures. The availability of a complete and high-quality map-based sequence has provided the opportunity to study genome organization and evolution. Most importantly, the order and identity of 37,544 genes of rice have been unraveled. The sequence provides the required ingredients for functional genomics and molecular breeding programs aimed at unraveling intricate cellular processes and improving rice productivity.

A microarray analysis of the rice transcriptome and its comparison to Arabidopsis

Arabidopsis and rice are the only two model plants whose finished phase genome sequence has been completed. Here we report the construction of an oligomer microarray based on the presently known and predicted gene models in the rice genome. This microarray was used to analyze the transcriptional activity of the gene models in representative rice organ types. Expression of 86% of the 41,754 known and predicted gene models was detected. A significant fraction of these expressed gene models are organized into chromosomal regions, about 100 kb in length, that exhibit a coexpression pattern. Compared with similar genome-wide surveys of the Arabidopsis transcriptome, our results indicate that similar proportions of the two genomes are expressed in their corresponding organ types. A large percentage of the rice gene models that lack significant Arabidopsis homologs are expressed. Furthermore, the expression patterns of rice and Arabidopsis best-matched homologous genes in distinct functional groups indicate dramatic differences in their degree of conservation between the two species. Thus, this initial comparative analysis reveals some basic similarities and differences between the Arabidopsis and rice transcriptomes.

Enhancing gene targeting efficiency in higher plants: rice is on the move

Meeting the challenge of routine gene targeting (GT) in higher plants is of crucial interest to researchers and plant breeders who are currently in need of a powerful tool to specifically modify a given locus in a genome. Higher plants have long been considered the last lineage resistant to targeting technology. However, a recent report described an efficient method of T-DNA-mediated targeted disruption of a non-selectable locus in rice [Terada et al., Nat Biotechnol 20: 1030-1034 (2002)]. Though this study was an obvious breakthrough, further improvement of GT frequencies may derive from a better understanding of the natural mechanisms that control homologous recombination (HR) processes. In this review, we will focus on what is known about HR and the factors which may hamper the development of routine GT by HR in higher plants. We will also present the current strategies envisaged to overcome these limitations, such as expression of recombination proteins and refinements in the design of the transformation vector.

Microarray and proteomic analysis of brassinosteroid- and gibberellin-regulated gene and protein expression in rice

Brassinosteroid (BR) and gibberellin (GA) are two groups of plant growth regulators essential for normal plant growth and development. To gain insight into the molecular mechanism by which BR and GA regulate the growth and development of plants, especially the monocot plant rice, it is necessary to identify and analyze more genes and proteins that are regulated by them. With the availability of draft sequences of two major types, japonica and indica rice, it has become possible to analyze expression changes of genes and proteins at genome scale. In this review, we summarize rice functional genomic research by using microarray and proteomic approaches and our recent research results focusing on the comparison of cDNA microarray and proteomic analyses of BR- and GA-regulated gene and protein expression in rice. We believe our findings have important implications for understanding the mechanism by which BR and GA regulate the growth and development of rice.

Identification by suppression subtractive hybridization of genes that are differentially expressed between near-isogenic maize lines in association with sugarcane mosaic virus resistance

The molecular mechanisms underlying the development and progression of sugarcane mosaic virus (SCMV) infection in maize are poorly understood. A study of differential expression was conducted to identify genes involved in resistance to SCMV. In this study, we combined suppression subtractive hybridization and macroarray hybridization to identify genes that are differently expressed in the near isogenic lines F7+ (SCMV resistant) and F7 (susceptible). Altogether, 302 differentially expressed genes were identified in four comparisons based on constitutively expressed and inducible genes, and on compatible and incompatible plant-virus interactions. Apart from genes related to metabolism, most of the functionally classified genes identified belonged to three pathogenesis-related categories: cell rescue, defense, cell death and ageing, signal transduction, and transcription. The latter three groups accounted for 56-66% of the genes classified. Some 19% (60 of 302) of the identified genes had previously been assigned to 29 bins distributed over all ten maize chromosomes. Among the mapped genes, 31% (18 of 58) were located within the Scmv2 and Scmv1 regions on chromosomes 3 and 6, respectively, which have been associated with resistance to SCMV. Promising candidate genes for Scmv1 have been identified, such as AA661457 (receptor-like kinase Xa21-binding protein 3). The implications of the genomic distribution of differentially expressed genes identified by this isogenic comparison are discussed in the context of breeding for resistance.

Rapid recent growth and divergence of rice nuclear genomes

By employing the nuclear DNA of the African rice Oryza glaberrima as a reference genome, the timing, natures, mechanisms, and specificities of recent sequence evolution in the indica and japonica subspecies of Oryza sativa were identified. The data indicate that the genome sizes of both indica and japonica have increased substantially, >2% and >6%, respectively, since their divergence from a common ancestor, mainly because of the amplification of LTR-retrotransposons. However, losses of all classes of DNA sequence through unequal homologous recombination and illegitimate recombination have attenuated the growth of the rice genome. Small deletions have been particularly frequent throughout the genome. In >1 Mb of orthologous regions that we analyzed, no cases of complete gene acquisition or loss from either indica or japonica were found, nor was any example of precise transposon excision detected. The sequences between genes were observed to have a very high rate of divergence, indicating a molecular clock for transposable elements that is at least 2-fold more rapid than synonymous base substitutions within genes. We found that regions prone to frequent insertions and deletions also exhibit higher levels of point mutation. These results indicate a highly dynamic rice genome with competing processes for the generation and removal of genetic variation.

Structural and functional analysis of rice genome

Rice is an excellent system for plant genomics as it represents a modest size genome of 430 Mb. It feeds more than half the population of the world. Draft sequences of the rice genome, derived by whole-genome shotgun approach at relatively low coverage (4-6 X), were published and the International Rice Genome Sequencing Project (IRGSP) declared high quality (>10 X), genetically anchored, phase 2 level sequence in 2002. In addition, phase 3 level finished sequence of chromosomes 1, 4 and 10 (out of 12 chromosomes of rice) has already been reported by scientists from IRGSP consortium. Various estimates of genes in rice place the number at >50,000. Already, over 28,000 full-length cDNAs have been sequenced, most of which map to genetically anchored genome sequence. Such information is very useful in revealing novel features of macro- and micro-level synteny of rice genome with other cereals. Microarray analysis is unraveling the identity of rice genes expressing in temporal and spatial manner and should help target candidate genes useful for improving traits of agronomic importance. Simultaneously, functional analysis of rice genome has been initiated by marker-based characterization of useful genes and employing functional knock-outs created by mutation or gene tagging. Integration of this enormous information is expected to catalyze tremendous activity on basic and applied aspects of rice genomics.

Whole genome sequence comparisons and "full-length" cDNA sequences: a combined approach to evaluate and improve Arabidopsis genome annotation

To evaluate the existing annotation of the Arabidopsis genome further, we generated a collection of evolutionary conserved regions (ecores) between Arabidopsis and rice. The ecore analysis provides evidence that the gene catalog of Arabidopsis is not yet complete, and that a number of these annotations require re-examination. To improve the Arabidopsis genome annotation further, we used a novel "full-length" enriched cDNA collection prepared from several tissues. An additional 1931 genes were covered by new "full-length" cDNA sequences, raising the number of annotated genes with a corresponding "full-length" cDNA sequence to about 14,000. Detailed comparisons between these "full-length" cDNA sequences and annotated genes show that this resource is very helpful in determining the correct structure of genes, in particular, those not yet supported by "full-length" cDNAs. In addition, a total of 326 genomic regions not included previously in the Arabidopsis genome annotation were detected by this cDNA resource, providing clues for new gene discovery. Because, as expected, the two data sets only partially overlap, their combination produces very useful information for improving the Arabidopsis genome annotation.

MIPS Arabidopsis thaliana Database (MAtDB): an integrated biological knowledge resource for plant genomics

Arabidopsis thaliana is the most widely studied model plant. Functional genomics is intensively underway in many laboratories worldwide. Beyond the basic annotation of the primary sequence data, the annotated genetic elements of Arabidopsis must be linked to diverse biological data and higher order information such as metabolic or regulatory pathways. The MIPS Arabidopsis thaliana database MAtDB aims to provide a comprehensive resource for Arabidopsis as a genome model that serves as a primary reference for research in plants and is suitable for transfer of knowledge to other plants, especially crops. The genome sequence as a common backbone serves as a scaffold for the integration of data, while, in a complementary effort, these data are enhanced through the application of state-of-the-art bioinformatics tools. This information is visualized on a genome-wide and a gene-by-gene basis with access both for web users and applications. This report updates the information given in a previous report and provides an outlook on further developments. The MAtDB web interface can be accessed at http://mips.gsf.de/proj/thal/db.

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.

Expressed sequence tags: alternative or complement to whole genome sequences?

Over three million sequences from approximately 200 plant species have been deposited in the publicly available plant expressed sequence tag (EST) sequence databases. Many of the ESTs have been sequenced as an alternative to complete genome sequencing or as a substrate for cDNA array-based expression analyses. This creates a formidable resource from both biodiversity and gene-discovery standpoints. Bioinformatics-based sequence analysis tools have extended the scope of EST analysis into the fields of proteomics, marker development and genome annotation. Although EST collections are certainly no substitute for a whole genome scaffold, this "poor man's genome" resource forms the core foundations for various genome-scale experiments within the as yet unsequenceable plant genomes.