Annotation of the Arabidopsis genome

Expression and functional analysis of various structural domains of tobacco topoisomerase II: To understand the mechanistic insights of plant type II topoisomerases

In contrast to bacterial, yeast and animal systems, topoisomerases (topo) from plants have not been well studied. In this report, we generated four truncated topoisomerase II (Topo II) cDNA fragments encoding different functional domains of Nicotiana tabacum topo II (NtTopoII). Each of these recombinant polypeptides was expressed alone or in combination in temperature-sensitive topoisomerase II yeast mutants. Recombinant NtTopoII with truncated polypeptides fails to target the yeast nuclei and does not rescue the temperature-sensitive phenotype. In contrast complementation was achieved with the full-length NtTopoII, which localized to the yeast nucleus. These observations suggested the presence of a potent nuclear localization signal (NLS) in the extreme C-terminal 314 amino acid residues of NtTopoII that functioned effectively in the heterologous yeast system. Biochemical characterization of purified recombinant full-length and the partial NtTopoII polypeptides revealed that the ATP-binding and hydrolysis region of NtTopoIIwas located at 413 amino acid N-terminal region and this ATPase domain is functional both when it is expressed as a separate polypeptide or as part of the holoenzyme. The present findings also revealed that all NtTopoII truncated polypeptides were detrimental for in vitro supercoiled DNA relaxation and/or DNA nicking and ligation activity. Further, we discuss the possible disruption of coordinated macromolecular interface movements and the dimer interactions in truncated NtTopoII that are required for functional topoisomerase activity.

Using the Arabidopsis Information Resource (TAIR) to Find Information About Arabidopsis Genes

The Arabidopsis Information Resource (TAIR; http://arabidopsis.org) is a comprehensive web resource of Arabidopsis biology for plant scientists. TAIR curates and integrates information about genes, proteins, gene function, orthologs, gene expression, mutant phenotypes, biological materials such as clones and seed stocks, genetic markers, genetic and physical maps, genome organization, images of mutant plants, protein sub-cellular localizations, publications, and the research community. The various data types are extensively interconnected and can be accessed through a variety of web-based search and display tools. This article primarily focuses on some basic methods for searching, browsing, visualizing, and analyzing information about Arabidopsis genes and genomes. Additionally, we describe how members of the community can share data via JBrowse and the Generic Online Annotation Submission Tool (GOAT) in order to make their published research more accessible and visible. © 2022 Wiley Periodicals LLC. Basic Protocol 1: TAIR homepage, sitemap, and navigation Basic Protocol 2: Finding comprehensive information about Arabidopsis genes Basic Protocol 3: Using the Arabidopsis genome browser: JBrowse Basic Protocol 4: Using the Gene Ontology annotations for gene discovery and gene function analysis Basic Protocol 5: Using gene lists to download bulk datasets Basic Protocol 6: Using TAIR's analysis tools to find short sequences and motifs Basic Protocol 7: Using the TAIR generic online annotation tool (GOAT) to submit functional annotations for Arabidopsis (or any other species) genes Basic Protocol 8: Using PhyloGenes to visualize gene families and predict functions Basic Protocol 9: Using TAIR to browse Arabidopsis literature Basic Protocol 10: Using the synteny viewer to find and display syntenic regions from Arabidopsis and other plant species.

Genome-Wide SNP discovery and genomic characterization in avocado (Persea americana Mill.)

Modern crop breeding is based on the use of genetically and phenotypically diverse plant material and, consequently, a proper understanding of population structure and genetic diversity is essential for the effective development of breeding programs. An example is avocado, a woody perennial fruit crop native to Mesoamerica with an increasing popularity worldwide. Despite its commercial success, there are important gaps in the molecular tools available to support on-going avocado breeding programs. In order to fill this gap, in this study, an avocado 'Hass' draft assembly was developed and used as reference to study 71 avocado accessions which represent the three traditionally recognized avocado horticultural races or subspecies (Mexican, Guatemalan and West Indian). An average of 5.72 M reads per individual and a total of 7,108 single nucleotide polymorphism (SNP) markers were produced for the 71 accessions analyzed. These molecular markers were used in a study of genetic diversity and population structure. The results broadly separate the accessions studied according to their botanical race in four main groups: Mexican, Guatemalan, West Indian and an additional group of Guatemalan × Mexican hybrids. The high number of SNP markers developed in this study will be a useful genomic resource for the avocado community.

Phosphatidylinositol 5 Phosphate (PI5P): From Behind the Scenes to the Front (Nuclear) Stage

Phosphatidylinositol (PI)-related signaling plays a pivotal role in many cellular aspects, including survival, cell proliferation, differentiation, DNA damage, and trafficking. PI is the core of a network of proteins represented by kinases, phosphatases, and lipases which are able to add, remove or hydrolyze PI, leading to different phosphoinositide products. Among the seven known phosphoinositides, phosphatidylinositol 5 phosphate (PI5P) was the last to be discovered. PI5P presence in cells is very low compared to other PIs. However, much evidence collected throughout the years has described the role of this mono-phosphoinositide in cell cycles, stress response, T-cell activation, and chromatin remodeling. Interestingly, PI5P has been found in different cellular compartments, including the nucleus. Here, we will review the nuclear role of PI5P, describing how it is synthesized and regulated, and how changes in the levels of this rare phosphoinositide can lead to different nuclear outputs.

Discovery of Novel Functional Centers With Rationally Designed Amino Acid Motifs

Plants are constantly exposed to environmental stresses and in part due to their sessile nature, they have evolved signal perception and adaptive strategies that are distinct from those of other eukaryotes. This is reflected at the cellular level where receptors and signalling molecules cannot be identified using standard homology-based searches querying with proteins from prokaryotes and other eukaryotes. One of the reasons for this is the complex domain architecture of receptor molecules. In order to discover hidden plant signalling molecules, we have developed a motif-based approach designed specifically for the identification of functional centers in plant molecules. This has made possible the discovery of novel components involved in signalling and stimulus-response pathways; the molecules include cyclic nucleotide cyclases, a nitric oxide sensor and a novel target for the hormone abscisic acid. Here, we describe the major steps of the method and illustrate it with recent and experimentally confirmed molecules as examples. We foresee that carefully curated search motifs supported by structural and bioinformatic assessments will uncover many more structural and functional aspects, particularly of signalling molecules.

Enhancing the GABI-Kat Arabidopsis thaliana T-DNA Insertion Mutant Database by Incorporating Araport11 Annotation

SimpleSearch provides access to a database containing information about T-DNA insertion lines of the GABI-Kat collection of Arabidopsis thaliana mutants. These mutants are an important tool for reverse genetics, and GABI-Kat is the second largest collection of such T-DNA insertion mutants. Insertion sites were deduced from flanking sequence tags (FSTs), and the database contains information about mutant plant lines as well as insertion alleles. Here, we describe improvements within the interface (available at http://www.gabi-kat.de/db/genehits.php) and with regard to the database content that have been realized in the last five years. These improvements include the integration of the Araport11 genome sequence annotation data containing the recently updated A. thaliana structural gene descriptions, an updated visualization component that displays groups of insertions with very similar insertion positions, mapped confirmation sequences, and primers. The visualization component provides a quick way to identify insertions of interest, and access to improved data about the exact structure of confirmed insertion alleles. In addition, the database content has been extended by incorporating additional insertion alleles that were detected during the confirmation process, as well as by adding new FSTs that have been produced during continued efforts to complement gaps in FST availability. Finally, the current database content regarding predicted and confirmed insertion alleles as well as primer sequences has been made available as downloadable flat files.

Triterpenoid saponin biosynthetic pathway profiling and candidate gene mining of the Ilex asprella root using RNA-Seq

Ilex asprella, which contains abundant α-amyrin type triterpenoid saponins, is an anti-influenza herbal drug widely used in south China. In this work, we first analysed the transcriptome of the I. asprella root using RNA-Seq, which provided a dataset for functional gene mining. mRNA was isolated from the total RNA of the I. asprella root and reverse-transcribed into cDNA. Then, the cDNA library was sequenced using an Illumina HiSeq™ 2000, which generated 55,028,452 clean reads. De novo assembly of these reads generated 51,865 unigenes, in which 39,269 unigenes were annotated (75.71% yield). According to the structures of the triterpenoid saponins of I. asprella, a putative biosynthetic pathway downstream of 2,3-oxidosqualene was proposed and candidate unigenes in the transcriptome data that were potentially involved in the pathway were screened using homology-based BLAST and phylogenetic analysis. Further amplification and functional analysis of these putative unigenes will provide insight into the biosynthesis of Ilex triterpenoid saponins.

A rapidly evolving genomic toolkit for Drosophila heterochromatin

Heterochromatin is the enigmatic eukaryotic genome compartment found mostly at telomeres and centromeres. Conventional approaches to sequence assembly and genetic manipulation fail in this highly repetitive, gene-sparse, and recombinationally silent DNA. In contrast, genetic and molecular analyses of euchromatin-encoded proteins that bind, remodel, and propagate heterochromatin have revealed its vital role in numerous cellular and evolutionary processes. Utilizing the 12 sequenced Drosophila genomes, Levine et al¹ took a phylogenomic approach to discover new such protein “surrogates” of heterochromatin function and evolution. This paper reported over 20 new members of what was traditionally believed to be a small and static Heterochromatin Protein 1 (HP1) gene family. The newly identified HP1 proteins are structurally diverse, lineage-restricted, and expressed primarily in the male germline. The birth and death of HP1 genes follows a “revolving door” pattern, where new HP1s appear to replace old HP1s. Here, we address alternative evolutionary models that drive this constant innovation.

Functional proteomics of barley and barley chloroplasts - strategies, methods and perspectives

Barley (Hordeum vulgare) is an important cereal grain that is used in a range of products for animal and human consumption. Crop yield and seed quality has been optimized during decades by plant breeding programs supported by biotechnology and molecular biology techniques. The recently completed whole-genome sequencing of barley revealed approximately 26,100 open reading frames, which provides a foundation for detailed molecular studies of barley by functional genomics and proteomics approaches. Such studies will provide further insights into the mechanisms of, for example, drought and stress tolerance, micronutrient utilization, and photosynthesis in barley. In the present review we present the current state of proteomics research for investigations of barley chloroplasts, i.e., the organelle that contain the photosynthetic apparatus in the plant. We describe several different proteomics strategies and discuss their applications in characterization of the barley chloroplast as well as future perspectives for functional proteomics in barley research.

Comparison of plasma membrane proteomic changes of Arabidopsis suspension-cultured cells (T87 Line) after cold and ABA treatment in association with freezing tolerance development

The plasma membrane (PM) is the primary site of freezing injury in plants. To determine global changes in PM protein profiles in association with freezing tolerance development, proteome analysis of the purified PM of Arabidopsis suspension-cultured cells (T87 line) was conducted with label-free protein quantification technology. Freezing tolerance of Arabidopsis cells at the lag growth phase (8 d old) increased after cold acclimation (CA) or ABA treatment. Proteome analysis assigned 658 proteins in the PM in total, of which 45.3% (298 proteins) were predicted to have transmembrane domains. They were classified into several functional categories, with the primary categories being proteins in transporters, signal transduction, protein destination and storage, and cell structure. After CA, 271 proteins increased and 111 proteins decreased. ABA treatment resulted in 185 increased and 56 decreased proteins. Of these, 139 increased and 49 decreased proteins were identified in common after both CA and ABA treatment. In addition, there were proteins specifically expressed in cold- (132 increased and 62 decreased) or ABA- (46 increased and 7 decreased) treated cells. Collectively, our results clearly show that (i) responses of the PM proteome to CA and ABA treatment overlap substantially but, at the same time, some proteins exhibited different response patterns in each treatment; and (ii) the majority of ABA-responsive proteins are CA-responsive proteins but not vice versa, suggesting complex interactions of CA and ABA signaling pathways in the PM proteome responses.

GABI-Kat SimpleSearch: new features of the Arabidopsis thaliana T-DNA mutant database

T-DNA insertion mutants are very valuable for reverse genetics in Arabidopsis thaliana. Several projects have generated large sequence-indexed collections of T-DNA insertion lines, of which GABI-Kat is the second largest resource worldwide. User access to the collection and its Flanking Sequence Tags (FSTs) is provided by the front end SimpleSearch (http://www.GABI-Kat.de). Several significant improvements have been implemented recently. The database now relies on the TAIRv10 genome sequence and annotation dataset. All FSTs have been newly mapped using an optimized procedure that leads to improved accuracy of insertion site predictions. A fraction of the collection with weak FST yield was re-analysed by generating new FSTs. Along with newly found predictions for older sequences about 20,000 new FSTs were included in the database. Information about groups of FSTs pointing to the same insertion site that is found in several lines but is real only in a single line are included, and many problematic FST-to-line links have been corrected using new wet-lab data. SimpleSearch currently contains data from ~71,000 lines with predicted insertions covering 62.5% of the 27,206 nuclear protein coding genes, and offers insertion allele-specific data from 9545 confirmed lines that are available from the Nottingham Arabidopsis Stock Centre.

Approaches to Fungal Genome Annotation

Fungal genome annotation is the starting point for analysis of genome content. This generally involves the application of diverse methods to identify features on a genome assembly such as protein-coding and non-coding genes, repeats and transposable elements, and pseudogenes. Here we describe tools and methods leveraged for eukaryotic genome annotation with a focus on the annotation of fungal nuclear and mitochondrial genomes. We highlight the application of the latest technologies and tools to improve the quality of predicted gene sets. The Broad Institute eukaryotic genome annotation pipeline is described as one example of how such methods and tools are integrated into a sequencing center's production genome annotation environment.

A physical map of Brassica oleracea shows complexity of chromosomal changes following recursive paleopolyploidizations

Background

Evolution of the Brassica species has been recursively affected by polyploidy events, and comparison to their relative, Arabidopsis thaliana, provides means to explore their genomic complexity.

Results

A genome-wide physical map of a rapid-cycling strain of B. oleracea was constructed by integrating high-information-content fingerprinting (HICF) of Bacterial Artificial Chromosome (BAC) clones with hybridization to sequence-tagged probes. Using 2907 contigs of two or more BACs, we performed several lines of comparative genomic analysis. Interspecific DNA synteny is much better preserved in euchromatin than heterochromatin, showing the qualitative difference in evolution of these respective genomic domains. About 67% of contigs can be aligned to the Arabidopsis genome, with 96.5% corresponding to euchromatic regions, and 3.5% (shown to contain repetitive sequences) to pericentromeric regions. Overgo probe hybridization data showed that contigs aligned to Arabidopsis euchromatin contain ~80% of low-copy-number genes, while genes with high copy number are much more frequently associated with pericentromeric regions. We identified 39 interchromosomal breakpoints during the diversification of B. oleracea and Arabidopsis thaliana, a relatively high level of genomic change since their divergence. Comparison of the B. oleracea physical map with Arabidopsis and other available eudicot genomes showed appreciable 'shadowing' produced by more ancient polyploidies, resulting in a web of relatedness among contigs which increased genomic complexity.

Conclusions

A high-resolution genetically-anchored physical map sheds light on Brassica genome organization and advances positional cloning of specific genes, and may help to validate genome sequence assembly and alignment to chromosomes.

All the physical mapping data is freely shared at a WebFPC site (http://lulu.pgml.uga.edu/fpc/WebAGCoL/brassica/WebFPC/; Temporarily password-protected: account: pgml; password: 123qwe123.

Genome-wide functional annotation and structural verification of metabolic ORFeome of Chlamydomonas reinhardtii

Background

Recent advances in the field of metabolic engineering have been expedited by the availability of genome sequences and metabolic modelling approaches. The complete sequencing of the C. reinhardtii genome has made this unicellular alga a good candidate for metabolic engineering studies; however, the annotation of the relevant genes has not been validated and the much-needed metabolic ORFeome is currently unavailable. We describe our efforts on the functional annotation of the ORF models released by the Joint Genome Institute (JGI), prediction of their subcellular localizations, and experimental verification of their structural annotation at the genome scale.

Results

We assigned enzymatic functions to the translated JGI ORF models of C. reinhardtii by reciprocal BLAST searches of the putative proteome against the UniProt and AraCyc enzyme databases. The best match for each translated ORF was identified and the EC numbers were transferred onto the ORF models. Enzymatic functional assignment was extended to the paralogs of the ORFs by clustering ORFs using BLASTCLUST.

In total, we assigned 911 enzymatic functions, including 886 EC numbers, to 1,427 transcripts. We further annotated the enzymatic ORFs by prediction of their subcellular localization. The majority of the ORFs are predicted to be compartmentalized in the cytosol and chloroplast. We verified the structure of the metabolism-related ORF models by reverse transcription-PCR of the functionally annotated ORFs. Following amplification and cloning, we carried out 454FLX and Sanger sequencing of the ORFs. Based on alignment of the 454FLX reads to the ORF predicted sequences, we obtained more than 90% coverage for more than 80% of the ORFs. In total, 1,087 ORF models were verified by 454 and Sanger sequencing methods. We obtained expression evidence for 98% of the metabolic ORFs in the algal cells grown under constant light in the presence of acetate.

Conclusions

We functionally annotated approximately 1,400 JGI predicted metabolic ORFs that can facilitate the reconstruction and refinement of a genome-scale metabolic network. The unveiling of the metabolic potential of this organism, along with structural verification of the relevant ORFs, facilitates the selection of metabolic engineering targets with applications in bioenergy and biopharmaceuticals. The ORF clones are a resource for downstream studies.

Molecular activities, biosynthesis and evolution of triterpenoid saponins

Saponins are bioactive compounds generally considered to be produced by plants to counteract pathogens and herbivores. Besides their role in plant defense, saponins are of growing interest for drug research as they are active constituents of several folk medicines and provide valuable pharmacological properties. Accordingly, much effort has been put into unraveling the modes of action of saponins, as well as in exploration of their potential for industrial processes and pharmacology. However, the exploitation of saponins for bioengineering crop plants with improved resistances against pests as well as circumvention of laborious and uneconomical extraction procedures for industrial production from plants is hampered by the lack of knowledge and availability of genes in saponin biosynthesis. Although the ability to produce saponins is rather widespread among plants, a complete synthetic pathway has not been elucidated in any single species. Current conceptions consider saponins to be derived from intermediates of the phytosterol pathway, and predominantly enzymes belonging to the multigene families of oxidosqualene cyclases (OSCs), cytochromes P450 (P450s) and family 1 UDP-glycosyltransferases (UGTs) are thought to be involved in their biosynthesis. Formation of unique structural features involves additional biosynthetical enzymes of diverse phylogenetic background. As an example of this, a serine carboxypeptidase-like acyltransferase (SCPL) was recently found to be involved in synthesis of triterpenoid saponins in oats. However, the total number of identified genes in saponin biosynthesis remains low as the complexity and diversity of these multigene families impede gene discovery based on sequence analysis and phylogeny. This review summarizes current knowledge of triterpenoid saponin biosynthesis in plants, molecular activities, evolutionary aspects and perspectives for further gene discovery.

Draft genome sequence of the oilseed species Ricinus communis

Castor bean (Ricinus communis) is an oilseed crop that belongs to the spurge (Euphorbiaceae) family, which comprises approximately 6,300 species that include cassava (Manihot esculenta), rubber tree (Hevea brasiliensis) and physic nut (Jatropha curcas). It is primarily of economic interest as a source of castor oil, used for the production of high-quality lubricants because of its high proportion of the unusual fatty acid ricinoleic acid. However, castor bean genomics is also relevant to biosecurity as the seeds contain high levels of ricin, a highly toxic, ribosome-inactivating protein. Here we report the draft genome sequence of castor bean (4.6-fold coverage), the first for a member of the Euphorbiaceae. Whereas most of the key genes involved in oil synthesis and turnover are single copy, the number of members of the ricin gene family is larger than previously thought. Comparative genomics analysis suggests the presence of an ancient hexaploidization event that is conserved across the dicotyledonous lineage.

A comparison of global gene expression measurement technologies in Arabidopsis thaliana

Microarrays and tag-based transcriptional profiling technologies represent diverse but complementary data types. We are currently conducting a comparison of high-density in situ synthesized microarrays and massively-parallel signature sequencing (MPSS) data in the model plant, Arabidopsis thaliana. The MPSS data (available at http://mpss.udel.edu/at) and the microarray data have been compiled using the same RNA source material. In this review, we outline the experimental strategy that we are using, and present preliminary data and interpretations from the transcriptional profiles of Arabidopsis leaves and roots. The preliminary data indicate that the log ratio differences of transcripts between leaves and roots measured by microarray data are in better agreement with the MPSS data than the absolute intensities measured for individual microarrays hybridized to only one of the cRNA populations. The correlation was substantially improved by focusing on a subset of genes excluding those with very low expression levels; this selection may have removed noisy data. Future reports will incorporate more than 10 tissues that have been sampled by MPSS.

FIDEL-a retrovirus-like retrotransposon and its distinct evolutionary histories in the A- and B-genome components of cultivated peanut

In this paper, we describe a Ty3-gypsy retrotransposon from allotetraploid peanut (Arachis hypogaea) and its putative diploid ancestors Arachis duranensis (A-genome) and Arachis ipaënsis (B-genome). The consensus sequence is 11,223 bp. The element, named FIDEL (Fairly long Inter-Dispersed Euchromatic LTR retrotransposon), is more frequent in the A- than in the B-genome, with copy numbers of about 3,000 (±950, A. duranensis), 820 (±480, A. ipaënsis), and 3,900 (±1,500, A. hypogaea) per haploid genome. Phylogenetic analysis of reverse transcriptase sequences showed distinct evolution of FIDEL in the ancestor species. Fluorescent in situ hybridization revealed disperse distribution in euchromatin and absence from centromeres, telomeric regions, and the nucleolar organizer region. Using paired sequences from bacterial artificial chromosomes, we showed that elements appear less likely to insert near conserved ancestral genes than near the fast evolving disease resistance gene homologs. Within the Ty3-gypsy elements, FIDEL is most closely related with the Athila/Calypso group of retrovirus-like retrotransposons. Putative transmembrane domains were identified, supporting the presence of a vestigial envelope gene. The results emphasize the importance of FIDEL in the evolution and divergence of different Arachis genomes and also may serve as an example of the role of retrotransposons in the evolution of legume genomes in general.

High-resolution melting analysis of cDNA-derived PCR amplicons for rapid and cost-effective identification of novel alleles in barley

An original method has been established for the identification of novel alleles of eukaryotic translation initiation factor 4E (eIF4E) gene, which is required for resistance to agronomically important bymoviruses, in barley germplasm. This method involves scanning for sequence variations in cDNA-derived PCR amplicons using High-resolution melting (HRM) followed by direct Sanger sequencing of only those amplicons which were predicted to carry nucleotide changes. HRM is a simple, cost-effective, rapid and high-throughput assay, which so far has only been widely used in clinical pathology for molecular diagnostic of diseases and patient genotyping. Application of HRM allowed significant reduction in the amount of expensive Sanger sequencing required for allele mining in plants. The method described here involved an investigation of total cDNA rather than genomic DNA, thus permitting the analyses of shorter (up to 300-bp) and fewer overlapping amplicons to cover the coding sequence. This strategy further reduced the allele mining costs. The sensitivity and accuracy of HRM for predicting genotypes carrying a wide range of nucleotide polymorphisms in eIF4E approached 100%. Results of the current study are promising and suggest that this method could also potentially be applied to the discovery of superior alleles controlling other important traits in barley as well in other model and crop plant species.

Establishing RNA interference as a reverse-genetic approach for gene functional analysis in protoplasts

Double-stranded (ds)RNA interference (RNAi) is widely used for functional analysis of plant genes and is achieved via generating stable transformants expressing dsRNA in planta. This study demonstrated that RNAi can also be utilized to examine gene functions in protoplasts. Because protoplasts are nongrowing cells, effective RNAi-triggered gene silencing depends not only on a depletion of gene transcripts but also on turnover rates of corresponding polypeptides. Herein, we tested if transient RNAi in protoplasts would result in the depletion of a targeted polypeptide and, because protoplasts have a limited life span, if functional assays of RNAi knockout genes would be feasible in protoplasts. We showed that protoplasts transfection with an in vitro-synthesized dsRNA against Arabidopsis (Arabidopsis thaliana) beta-glutamylcysteine synthase (ECS1), a key enzyme in the synthesis of glutathione, resulted in a 95% depletion of ECS1 transcript, a 72% decrease of ECS1 polypeptide, and a 60% drop in glutathione content. These results were comparable with those obtained upon analysis of Arabidopsis seedlings bearing the cad2-1 mutant allele of ECS1. We also improved the procedure for RNAi inactivation of several genes simultaneously. Finally, because we isolated protoplasts from tissues of 14-d-old seedlings instead of 1-month-old mature plants, the described procedure is rapid (as it only takes 20 d from seed planting to functional studies), suitable for analyzing multiple genes in parallel, and independent of cloning dsRNAs into plant expression vectors. Therefore, RNAi in protoplasts complements existing genetic tools, as it allows rapid, cost- and space-efficient initial screening and selection of genes for subsequent in planta studies.

Gene Structure Annotation at PlantGDB

The accurate identification of exons and introns that comprise a complete plant gene structure can be a time-consuming and challenging task. Novel Web-based tools facilitate the process by providing a convenient interface to current transcript evidence, and portals to relevant bioinformatics software. With a few keystrokes, the user can explore alternative transcript assemblies and, for example, select for annotation those that are clearly supported by transcript evidence and similarity to known genes. The implementation of the tool at the PlantGDB resource also allows immediate communication of the novel annotations to the community through Web display.

Genome sequencing approaches and successes

Sequence data is crucial to our understanding of crop growth and development, as differences in DNA sequence are responsible for almost all of the heritable differences between crop varieties and ecotypes. The sequence of a genome is often referred to as the genetic blueprint, and is the foundation for all additional information from the genome to the phenome. The value of DNA sequence is leading to rapid improvements in sequencing technology, increasing throughput, and reducing costs, and technological advances are accelerating with the introduction of novel approaches that are replacing the traditional Sanger-based methods. As genome sequencing becomes cheaper, it will be applied to a greater number of species with increasingly large and complex genomes. This will increase our understanding of how differences in the sequence relate to phenotypic observations, heritable traits, speciation, and evolution. Our understanding of plants will be greatly enhanced by this flow of sequence information, with direct benefit for crop improvement.

Sequencing and comparative analysis of a conserved syntenic segment in the Solanaceae

Comparative genomics is a powerful tool for gaining insight into genomic function and evolution. However, in plants, sequence data that would enable detailed comparisons of both coding and noncoding regions have been limited in availability. Here we report the generation and analysis of sequences for an unduplicated conserved syntenic segment (CSS) in the genomes of five members of the agriculturally important plant family Solanaceae. This CSS includes a 105-kb region of tomato chromosome 2 and orthologous regions of the potato, eggplant, pepper, and petunia genomes. With a total neutral divergence of 0.73-0.78 substitutions/site, these sequences are similar enough that most noncoding regions can be aligned, yet divergent enough to be informative about evolutionary dynamics and selective pressures. The CSS contains 17 distinct genes with generally conserved order and orientation, but with numerous small-scale differences between species. Our analysis indicates that the last common ancestor of these species lived approximately 27-36 million years ago, that more than one-third of short genomic segments (5-15 bp) are under selection, and that more than two-thirds of selected bases fall in noncoding regions. In addition, we identify genes under positive selection and analyze hundreds of conserved noncoding elements. This analysis provides a window into 30 million years of plant evolution in the absence of polyploidization.

Using the Arabidopsis Information Resource (TAIR) to find information about Arabidopsis genes

The Arabidopsis Information Resource (TAIR; http://www.arabidopsis.org) is a comprehensive Web resource of Arabidopsis biology for plant scientists. TAIR curates and integrates information about genes, proteins, gene expression, mutant phenotypes, biological materials such as DNA and seed stocks, genetic markers, genetic and physical maps, biochemical pathways, genome organization, images of mutant plants and protein sub-cellular localizations, publications, and the research community Data in TAIR are extensively interconnected and can be accessed through a variety of Web-based search and display tools. This unit primarily focuses on some basic methods for searching, browsing, visualizing, and analyzing information about Arabidopsis genes. Gene expression data from microarrays is a recent addition to the database and methods for accessing these data are also described. Two pattern identification programs are described for mining TAIR's unique Arabidopsis sequence data sets. We also describe how to use AraCyc for mining plant metabolic pathways.

Annotating genes of known and unknown function by large-scale coexpression analysis

About 40% of the proteins encoded in eukaryotic genomes are proteins of unknown function (PUFs). Their functional characterization remains one of the main challenges in modern biology. In this study we identified the PUF encoding genes from Arabidopsis (Arabidopsis thaliana) using a combination of sequence similarity, domain-based, and empirical approaches. Large-scale gene expression analyses of 1,310 publicly available Affymetrix chips were performed to associate the identified PUF genes with regulatory networks and biological processes of known function. To generate quality results, the study was restricted to expression sets with replicated samples. First, genome-wide clustering and gene function enrichment analysis of clusters allowed us to associate 1,541 PUF genes with tightly coexpressed genes for proteins of known function (PKFs). Over 70% of them could be assigned to more specific biological process annotations than the ones available in the current Gene Ontology release. The most highly overrepresented functional categories in the obtained clusters were ribosome assembly, photosynthesis, and cell wall pathways. Interestingly, the majority of the PUF genes appeared to be controlled by the same regulatory networks as most PKF genes, because clusters enriched in PUF genes were extremely rare. Second, large-scale analysis of differentially expressed genes was applied to identify a comprehensive set of abiotic stress-response genes. This analysis resulted in the identification of 269 PKF and 104 PUF genes that responded to a wide variety of abiotic stresses, whereas 608 PKF and 206 PUF genes responded predominantly to specific stress treatments. The provided coexpression and differentially expressed gene data represent an important resource for guiding future functional characterization experiments of PUF and PKF genes. Finally, the public Plant Gene Expression Database (http://bioweb.ucr.edu/PED) was developed as part of this project to provide efficient access and mining tools for the vast gene expression data of this study.

Metabolomics- and proteomics-assisted genome annotation and analysis of the draft metabolic network of Chlamydomonas reinhardtii

We present an integrated analysis of the molecular repertoire of Chlamydomonas reinhardtii under reference conditions. Bioinformatics annotation methods combined with GCxGC/MS-based metabolomics and LC/MS-based shotgun proteomics profiling technologies have been applied to characterize abundant proteins and metabolites, resulting in the detection of 1069 proteins and 159 metabolites. Of the measured proteins, 204 currently do not have EST sequence support; thus a significant portion of the proteomics-detected proteins provide evidence for the validity of in silico gene models. Furthermore, the generated peptide data lend support to the validity of a number of proteins currently in the proposed model stage. By integrating genomic annotation information with experimentally identified metabolites and proteins, we constructed a draft metabolic network for Chlamydomonas. Computational metabolic modeling allowed an identification of missing enzymatic links. Some experimentally detected metabolites are not producible by the currently known and annotated enzyme set, thus suggesting entry points for further targeted gene discovery or biochemical pathway research. All data sets are made available as supplementary material as well as web-accessible databases and within the functional context via the Chlamydomonas-adapted MapMan annotation platform. Information of identified peptides is also available directly via the JGI-Chlamydomonas genomic resource database (http://genome.jgi-psf.org/Chlre3/Chlre3.home.html).

Analysis of alternative splicing in plants with bioinformatics tools

Alternative splicing is a molecular mechanism utilized by a broad range of eukaryotes to extend the repertoire of functions encoded by single genes and to posttranscriptionally regulate gene expression. Recent analyses of expressed transcript sequences aligned to the complete genomes of Arabidopsis and rice indicate that alternative splicing in plants is prevalent and exhibits several features similar to other higher eukaryotes including mouse and human. This chapter reviews the computational strategies employed to study alternative splicing with bioinformatics tools and the recent findings from analyses performed on plants by applying such methods.

Fluorescent protein applications in plants

Study of plant cell biology has benefited tremendously from the use of fluorescent proteins (FPs). Development of well-established techniques in genetics, by transient expression or by Agrobacterium-mediated plant cell transformation, makes it possible to readily create material for imaging molecules tagged with FPs. Confocal microscopy of FPs is routine and, in highly scattering tissues, multiphoton microscopy improves deep imaging. The abundance of autofluorescent compounds in plants in some cases potentially interferes with FP signals, but spectral imaging is an effective tool in unmixing overlapping signals. This approach allows separate detection of DsRed and chlorophyll, DsRed and GFP, and green fluorescent protein (GFP) and yellow fluorescent protein (YFP). FPs have been targeted to most plant organelles. Free (untargeted) FPs in plant cells are not only cytoplasmic, but also go into the nucleus due to their small size. FP fluorescence is potentially unstable in acidic vacuoles. FPs have been targeted to novel compartments, including protein storage vacuoles in seeds. Endoplasmic reticulum (ER)-targeted GFP has identified novel inclusion bodies that are surprisingly dynamic. FP-tagged Rab GTPases have allowed documentation of the dynamics of membrane trafficking. Investigation of virus infections has progressed significantly with the aid of FP-tagged virus proteins. Advanced techniques are giving plant scientists the ability to quantitatively analyze the behavior of FP-tagged proteins. Fluorescence lifetime microscopy is becoming the method of choice for fluorescence resonance energy transfer (FRET) analysis of FP-tagged proteins. Fluorescence correlation spectroscopy (FCS) of FPs provides information on molecular diffusion and intermolecular interactions. Use of FPs in elucidating the behavior of plant cells has a bright future.

Disruptions in AUX1-dependent auxin influx alter hypocotyl phototropism in Arabidopsis

Phototropism represents a differential growth response by which plant organs can respond adaptively to changes in the direction of incident light to optimize leaf/stem positioning for photosynthetic light capture and root growth orientation for water/nutrient acquisition. Studies over the past few years have identified a number of components in the signaling pathway(s) leading to development of phototropic curvatures in hypocotyls. These include the phototropin photoreceptors (phot1 and phot2) that perceive directional blue-light (BL) cues and then stimulate signaling, leading to relocalization of the plant hormone auxin, as well as the auxin response factor NPH4/ARF7 that responds to changes in local auxin concentrations to directly mediate expression of genes likely encoding proteins necessary for development of phototropic curvatures. While null mutations in NPH4/ARF7 condition an aphototropic response to unidirectional BL, seedlings carrying the same mutations recover BL-dependent phototropic responsiveness if co-irradiated with red light (RL) or pre-treated with either ethylene. In the present study, we identify second-site enhancer mutations in the nph4 background that abrogate these recovery responses. One of these mutations--map1 (modifier of arf7 phenotypes 1)--was found to represent a missense allele of AUX1--a gene encoding a high-affinity auxin influx carrier previously associated with a number of root responses. Pharmacological studies and analyses of additional aux1 mutants confirmed that AUX1 functions as a modulator of hypocotyl phototropism. Moreover, we have found that the strength of dependence of hypocotyl phototropism on AUX1-mediated auxin influx is directly related to the auxin responsiveness of the seedling in question.

The cytochrome P450 enzyme CYP96A15 is the midchain alkane hydroxylase responsible for formation of secondary alcohols and ketones in stem cuticular wax of Arabidopsis

Most aerial surfaces of plants are covered by cuticular wax that is synthesized in epidermal cells. The wax mixture on the inflorescence stems of Arabidopsis (Arabidopsis thaliana) is dominated by alkanes, secondary alcohols, and ketones, all thought to be formed sequentially in the decarbonylation pathway of wax biosynthesis. Here, we used a reverse-genetic approach to identify a cytochrome P450 enzyme (CYP96A15) involved in wax biosynthesis and characterized it as a midchain alkane hydroxylase (MAH1). Stem wax of T-DNA insertional mutant alleles was found to be devoid of secondary alcohols and ketones (mah1-1) or to contain much lower levels of these components (mah1-2 and mah1-3) than wild type. All mutant lines also had increased alkane amounts, partially or fully compensating for the loss of other compound classes. In spite of the chemical variation between mutant and wild-type waxes, there were no discernible differences in the epicuticular wax crystals on the stem surfaces. Mutant stem wax phenotypes could be partially rescued by expression of wild-type MAH1 under the control of the native promoter as well as the cauliflower mosaic virus 35S promoter. Cauliflower mosaic virus 35S-driven overexpression of MAH1 led to ectopic accumulation of secondary alcohols and ketones in Arabidopsis leaf wax, where only traces of these compounds are found in the wild type. The newly formed leaf alcohols and ketones had midchain functional groups on or next to the central carbon, thus matching those compounds in wild-type stem wax. Taken together, mutant analyses and ectopic expression of MAH1 in leaves suggest that this enzyme can catalyze the hydroxylation reaction leading from alkanes to secondary alcohols and possibly also a second hydroxylation leading to the corresponding ketones. MAH1 expression was largely restricted to the expanding regions of the inflorescence stems, specifically to the epidermal pavement cells, but not in trichomes and guard cells. MAH1-green fluorescent protein fusion proteins localized to the endoplasmic reticulum, providing evidence that both intermediate and final products of the decarbonylation pathway are generated in this subcellular compartment and must subsequently be delivered to the plasma membrane for export toward the cuticle.

EuCAP, a Eukaryotic Community Annotation Package, and its application to the rice genome

BACKGROUND

Despite the improvements of tools for automated annotation of genome sequences, manual curation at the structural and functional level can provide an increased level of refinement to genome annotation. The Institute for Genomic Research Rice Genome Annotation (hereafter named the Osa1 Genome Annotation) is the product of an automated pipeline and, for this reason, will benefit from the input of biologists with expertise in rice and/or particular gene families. Leveraging knowledge from a dispersed community of scientists is a demonstrated way of improving a genome annotation. This requires tools that facilitate 1) the submission of gene annotation to an annotation project, 2) the review of the submitted models by project annotators, and 3) the incorporation of the submitted models in the ongoing annotation effort.

RESULTS

We have developed the Eukaryotic Community Annotation Package (EuCAP), an annotation tool, and have applied it to the rice genome. The primary level of curation by community annotators (CA) has been the annotation of gene families. Annotation can be submitted by email or through the EuCAP Web Tool. The CA models are aligned to the rice pseudomolecules and the coordinates of these alignments, along with functional annotation, are stored in the MySQL EuCAP Gene Model database. Web pages displaying the alignments of the CA models to the Osa1 Genome models are automatically generated from the EuCAP Gene Model database. The alignments are reviewed by the project annotators (PAs) in the context of experimental evidence. Upon approval by the PAs, the CA models, along with the corresponding functional annotations, are integrated into the Osa1 Genome Annotation. The CA annotations, grouped by family, are displayed on the Community Annotation pages of the project website http://rice.tigr.org, as well as in the Community Annotation track of the Genome Browser.

CONCLUSION

We have applied EuCAP to rice. As of July 2007, the structural and/or functional annotation of 1,094 genes representing 57 families have been deposited and integrated into the current gene set. All of the EuCAP components are open-source, thereby allowing the implementation of EuCAP for the annotation of other genomes. EuCAP is available at http://sourceforge.net/projects/eucap/.

Large-scale polymorphism of heterochromatic repeats in the DNA of Arabidopsis thaliana

BACKGROUND

The composition of the individual eukaryote's genome and its variation within a species remain poorly defined. Even for a sequenced genome such as that of the model plant Arabidopsis thaliana accession Col-0, the large arrays of heterochromatic repeats are incompletely sequenced, with gaps of uncertain size persisting in them.

RESULTS

Using geographically separate populations of A. thaliana, we assayed variation in the heterochromatic repeat arrays using two independent methods and identified significant polymorphism among them, with variation by as much as a factor of two in the centromeric 180 bp repeat, in the 45S rDNA arrays and in the Athila retroelements. In the accession with highest genome size as measured by flow cytometry, Loh-0, we found more than a two-fold increase in 5S RNA gene copies relative to Col-0; results from fluorescence in situ hybridization with 5S probes were consistent with the existence of size polymorphism between Loh-0 and Col-0 at the 5S loci. Comparative genomic hybridization results of Loh-0 and Col-0 did not support contiguous variation in copy number of protein-coding genes on the scale needed to explain their observed genome size difference. We developed a computational data model to test whether the variation we measured in the repeat fractions could account for the different genome sizes determined with flow cytometry, and found that this proposed relationship could account for about 50% of the variance in genome size among the accessions.

CONCLUSION

Our analyses are consistent with substantial repeat number polymorphism for 5S and 45S ribosomal genes among accession of A. thaliana. Differences are also suggested for centromeric and pericentromeric repeats. Our analysis also points to the difficulties in measuring the repeated fraction of the genome and suggests that independent validation of genome size should be sought in addition to flow cytometric measurements.

Arabidopsis Gene Family Profiler (aGFP)--user-oriented transcriptomic database with easy-to-use graphic interface

BACKGROUND

Microarray technologies now belong to the standard functional genomics toolbox and have undergone massive development leading to increased genome coverage, accuracy and reliability. The number of experiments exploiting microarray technology has markedly increased in recent years. In parallel with the rapid accumulation of transcriptomic data, on-line analysis tools are being introduced to simplify their use. Global statistical data analysis methods contribute to the development of overall concepts about gene expression patterns and to query and compose working hypotheses. More recently, these applications are being supplemented with more specialized products offering visualization and specific data mining tools. We present a curated gene family-oriented gene expression database, Arabidopsis Gene Family Profiler (aGFP; http://agfp.ueb.cas.cz), which gives the user access to a large collection of normalised Affymetrix ATH1 microarray datasets. The database currently contains NASC Array and AtGenExpress transcriptomic datasets for various tissues at different developmental stages of wild type plants gathered from nearly 350 gene chips.

RESULTS

The Arabidopsis GFP database has been designed as an easy-to-use tool for users needing an easily accessible resource for expression data of single genes, pre-defined gene families or custom gene sets, with the further possibility of keyword search. Arabidopsis Gene Family Profiler presents a user-friendly web interface using both graphic and text output. Data are stored at the MySQL server and individual queries are created in PHP script. The most distinguishable features of Arabidopsis Gene Family Profiler database are: 1) the presentation of normalized datasets (Affymetrix MAS algorithm and calculation of model-based gene-expression values based on the Perfect Match-only model); 2) the choice between two different normalization algorithms (Affymetrix MAS4 or MAS5 algorithms); 3) an intuitive interface; 4) an interactive "virtual plant" visualizing the spatial and developmental expression profiles of both gene families and individual genes.

CONCLUSION

Arabidopsis GFP gives users the possibility to analyze current Arabidopsis developmental transcriptomic data starting with simple global queries that can be expanded and further refined to visualize comparative and highly selective gene expression profiles.

Small cysteine-rich peptides resembling antimicrobial peptides have been under-predicted in plants

Multicellular organisms produce small cysteine-rich antimicrobial peptides as an innate defense against pathogens. While defensins, a well-known class of such peptides, are common among eukaryotes, there are other classes restricted to the plant kingdom. These include thionins, lipid transfer proteins and snakins. In earlier work, we identified several divergent classes of small putatively secreted cysteine-rich peptides (CRPs) in legumes [Graham et al. (2004)Plant Physiol. 135, 1179-97]. Here, we built sequence motif models for each of these classes of peptides, and iteratively searched for related sequences within the comprehensive UniProt protein dataset, the Institute for Genomic Research's 33 plant gene indices, and the entire genomes of the model dicot, Arabidopsis thaliana, and the model monocot and crop species, Oryza sativa (rice). Using this search strategy, we identified approximately 13,000 plant genes encoding peptides with common features: (i) an N-terminal signal peptide, (ii) a small divergent charged or polar mature peptide with conserved cysteines, (iii) a similar intron/exon structure, (iv) spatial clustering in the genomes studied, and (v) overrepresentation in expressed sequences from reproductive structures of specific taxa. The identified genes include classes of defensins, thionins, lipid transfer proteins, and snakins, plus other protease inhibitors, pollen allergens, and uncharacterized gene families. We estimate that these classes of genes account for approximately 2-3% of the gene repertoire of each model species. Although 24% of the genes identified were not annotated in the latest Arabidopsis genome releases (TIGR5, TAIR6), we confirmed expression via RT-PCR for 59% of the sequences attempted. These findings highlight limitations in current annotation procedures for small divergent peptide classes.

Comparative study of SBP-box gene family in Arabidopsis and rice

SBP-box proteins are plant-specific putative transcription factors, which contain highly conserved SBP domain and could bind specifically to promoters of the floral meristem identity gene SQUAMOSA and its orthologous genes to regulate their expressions. In this study, 17 non-redundant SBP-box genes in Arabidopsis genome and 19 in rice genome were identified by using the known SBP domain sequences as queries. The phylogenetic analysis suggested that the main characteristics of this family might have been in existence before the split of Arabidopsis and rice, and most SBP-box genes expanded in a species-specific manner after the split of monocotyledon and dicotyledon. All the SBP-box proteins were classified into 9 subgroups based on the phylogenetic tree, where each group shared similar motifs and the orders of the motifs in the same group were found almost identical. Analysis of nonsynonymous and synonymous substitution rates revealed that the SBP domain had gone through purifying selection, whereas some regions outside SBP domain had gone through positive or relaxed purifying selection. The expression patterns of the SBP-box genes were further investigated by searching against the EST database. Results showed that the Arabidopsis SBP-box genes are expressed chiefly in flowers, leaves, roots and seeds, while those in rice mainly in flowers and callus.

Transcriptional profiling of the Arabidopsis embryo

We have used laser-capture microdissection to isolate RNA from discrete tissues of globular, heart, and torpedo stage embryos of Arabidopsis (Arabidopsis thaliana). This was amplified and analyzed by DNA microarray using the Affymetrix ATH1 GeneChip, representing approximately 22,800 Arabidopsis genes. Cluster analysis showed that spatial differences in gene expression were less significant than temporal differences. Time course analysis reveals the dynamics and complexity of gene expression in both apical and basal domains of the developing embryo, with several classes of synexpressed genes identifiable. The transition from globular to heart stage is associated in particular with an up-regulation of genes involved in cell cycle control, transcriptional regulation, and energetics and metabolism. The transition from heart to torpedo stage is associated with a repression of cell cycle genes and an up-regulation of genes encoding storage proteins, and pathways of cell growth, energy, and metabolism. The torpedo stage embryo shows strong functional differentiation in the root and cotyledon, as inferred from the classes of genes expressed in these tissues. The time course of expression of the essential EMBRYO-DEFECTIVE genes shows that most are expressed at unchanging levels across all stages of embryogenesis. We show how identified genes can be used to generate cell type-specific markers and promoter activities for future application in cell biology.

Plant database resources at The Institute for Genomic Research

With the completion of the genome sequences of the model plants Arabidopsis and rice, and the continuing sequencing efforts of other economically important crop plants, an unprecedented amount of genome sequence data is now available for large-scale genomics studies and analyses, such as the identification and discovery of novel genes, comparative genomics, and functional genomics. Efficient utilization of these large data sets is critically dependent on the ease of access and organization of the data. The plant databases at The Institute for Genomic Research (TIGR) have been set up to maintain various data types including genomic sequence, annotation and analyses, expressed transcript assemblies and analyses, and gene expression profiles from microarray studies. We present here an overview of the TIGR database resources for plant genomics and describe methods to access the data.

AtEnsEMBL

A comprehensive Arabidopsis genomic resource has been developed at the Nottingham Arabidopsis Stock Centre (NASC) to support the international plant community. This browser, termed AtEnsembl, provides a detailed and user-friendly interface for accessing a wide range of Arabidopsis-based genomic information and post-genomic resources using the Ensembl browser. The resource aims to provide the broadest possible range of Ensembl features, including pointers to germplasm as well as representations of gene and protein information, links to Affymetrix gene expression data, and extensive data download capabilities.

GABI-Kat SimpleSearch: an Arabidopsis thaliana T-DNA mutant database with detailed information for confirmed insertions

Insertional mutagenesis approaches, especially by T-DNA, play important roles in gene function studies of the model plant Arabidopsis thaliana. GABI-Kat SimpleSearch (http://www.GABI-Kat.de) is a Flanking Sequence Tag (FST)-based database for T-DNA insertion mutants generated by the GABI-Kat project. Currently, the database contains >108,000 mapped FSTs from approximately 64,000 lines which cover 64% of all annotated A.thaliana protein-coding genes. The web interface allows searching for relevant insertions by gene code, keyword, line identifier, GenBank accession number of the FST, and also by BLAST. A graphic display of the genome region around the gene or the FST assists users to select insertion lines of their interests. About 3500 insertions were confirmed in the offspring of the plant from which the original FST was generated, and the seeds of these lines are available from the Nottingham Arabidopsis Stock Centre. The database now also contains additional information such as segregation data, gene-specific primers and confirmation sequences. This information not only helps users to evaluate the usefulness of the mutant lines, but also covers a big part of the molecular characterization of the insertion alleles.

MicroRNAs and other small RNAs enriched in the Arabidopsis RNA-dependent RNA polymerase-2 mutant

The Arabidopsis genome contains a highly complex and abundant population of small RNAs, and many of the endogenous siRNAs are dependent on RNA-Dependent RNA Polymerase 2 (RDR2) for their biogenesis. By analyzing an rdr2 loss-of-function mutant using two different parallel sequencing technologies, MPSS and 454, we characterized the complement of miRNAs expressed in Arabidopsis inflorescence to considerable depth. Nearly all known miRNAs were enriched in this mutant and we identified 13 new miRNAs, all of which were relatively low abundance and constitute new families. Trans-acting siRNAs (ta-siRNAs) were even more highly enriched. Computational and gel blot analyses suggested that the minimal number of miRNAs in Arabidopsis is approximately 155. The size profile of small RNAs in rdr2 reflected enrichment of 21-nt miRNAs and other classes of siRNAs like ta-siRNAs, and a significant reduction in 24-nt heterochromatic siRNAs. Other classes of small RNAs were found to be RDR2-independent, particularly those derived from long inverted repeats and a subset of tandem repeats. The small RNA populations in other Arabidopsis small RNA biogenesis mutants were also examined; a dcl2/3/4 triple mutant showed a similar pattern to rdr2, whereas dcl1-7 and rdr6 showed reductions in miRNAs and ta-siRNAs consistent with their activities in the biogenesis of these types of small RNAs. Deep sequencing of mutants provides a genetic approach for the dissection and characterization of diverse small RNA populations and the identification of low abundance miRNAs.

Genome comparison using Gene Ontology (GO) with statistical testing

BACKGROUND

Automated comparison of complete sets of genes encoded in two genomes can provide insight on the genetic basis of differences in biological traits between species. Gene ontology (GO) is used as a common vocabulary to annotate genes for comparison. Current approaches calculate the fold of unweighted or weighted differences between two species at the high-level GO functional categories. However, to ensure the reliability of the differences detected, it is important to evaluate their statistical significance. It is also useful to search for differences at all levels of GO.

RESULTS

We propose a statistical approach to find reliable differences between the complete sets of genes encoded in two genomes at all levels of GO. The genes are first assigned GO terms from BLAST searches against genes with known GO assignments, and for each GO term the abundance of genes in the two genomes is compared using a chi-squared test followed by false discovery rate (FDR) correction. We applied this method to find statistically significant differences between two cyanobacteria, Synechocystis sp. PCC6803 and Anabaena sp. PCC7120. We then studied how the set of identified differences vary when different BLAST cutoffs are used. We also studied how the results vary when only subsets of the genes were used in the comparison of human vs. mouse and that of Saccharomyces cerevisiae vs. Schizosaccharomyces pombe.

CONCLUSION

There is a surprising lack of statistical approaches for comparing complete genomes at all levels of GO. With the rapid increase of the number of sequenced genomes, we hope that the approach we proposed and tested can make valuable contribution to comparative genomics.

Systematic analysis of Arabidopsis organelles and a protein localization database for facilitating fluorescent tagging of full-length Arabidopsis proteins

Cells are organized into a complex network of subcellular compartments that are specialized for various biological functions. Subcellular location is an important attribute of protein function. To facilitate systematic elucidation of protein subcellular location, we analyzed experimentally verified protein localization data of 1,300 Arabidopsis (Arabidopsis thaliana) proteins. The 1,300 experimentally verified proteins are distributed among 40 different compartments, with most of the proteins localized to four compartments: mitochondria (36%), nucleus (28%), plastid (17%), and cytosol (13.3%). About 19% of the proteins are found in multiple compartments, in which a high proportion (36.4%) is localized to both cytosol and nucleus. Characterization of the overrepresented Gene Ontology molecular functions and biological processes suggests that the Golgi apparatus and peroxisome may play more diverse functions but are involved in more specialized processes than other compartments. To support systematic empirical determination of protein subcellular localization using a technology called fluorescent tagging of full-length proteins, we developed a database and Web application to provide preselected green fluorescent protein insertion position and primer sequences for all Arabidopsis proteins to study their subcellular localization and to store experimentally verified protein localization images, videos, and their annotations of proteins generated using the fluorescent tagging of full-length proteins technology. The database can be searched, browsed, and downloaded using a Web browser at http://aztec.stanford.edu/gfp/. The software can also be downloaded from the same Web site for local installation.

The maize genome as a model for efficient sequence analysis of large plant genomes

The genomes of flowering plants vary in size from about 0.1 to over 100 gigabase pairs (Gbp), mostly because of polyploidy and variation in the abundance of repetitive elements in intergenic regions. High-quality sequences of the relatively small genomes of Arabidopsis (0.14 Gbp) and rice (0.4 Gbp) have now been largely completed. The sequencing of plant genomes that have a more representative size (the mean for flowering plant genomes is 5.6 Gbp) has been seen as a daunting task, partly because of their size and partly because of the numerous highly conserved repeats. Nevertheless, creative strategies and powerful new tools have been generated recently in the plant genetics community, so that sequencing large plant genomes is now a realistic possibility. Maize (2.4-2.7 Gbp) will be the first gigabase-size plant genome to be sequenced using these novel approaches. Pilot studies on maize indicate that the new gene-enrichment, gene-finishing and gene-orientation technologies are efficient, robust and comprehensive. These strategies will succeed in sequencing the gene-space of large genome plants, and in locating all of these genes and adjacent sequences on the genetic and physical maps.

Generating single-copy nuclear gene data for a recent adaptive radiation

Recent adaptive radiations provide an exceptional opportunity to understand the processes of speciation and adaptation. However, reconstructing the phylogenetic history of recent and rapidly evolving clades often requires the use of multiple, independent gene genealogies. Nuclear introns are an obvious source of the necessary data but their use is often limited because degenerate primers can amplify paralogous loci. To identify PCR primers for a large number of loci in an especially rapid adaptive radiation, that of the flowering plant genus Aquilegia (Ranunculaceae), we developed an efficient method for amplifying multiple single-copy nuclear loci by sequencing a modest number of clones from a cDNA library and designing PCR primers; with one primer anchored in the 3' untranslated region (3'-UTR) and one primer in the coding region of each gene. Variation between paralogous loci evolves more quickly in 3'-UTR regions compared to adjacent exons, and therefore we achieved high specificity for isolating orthologous loci. Furthermore, we were able to identify genes containing large introns by amplifying genes from genomic DNA and comparing the PCR product size to that predicted from their cDNA sequence. In Aquilegia eight out of eleven loci were isolated with this method and six of these loci had introns. Among four genes sequenced for samples spanning the phylogenetic breadth of the genus, we found sequence variation at levels similar to that observed in ITS, further supporting the recent and rapid radiation in Aquilegia. We assessed the orthology of amplification products by phylogenetic congruence among loci, the presence of two well established phylogenetic relationships, and similarity among loci for levels of sequence variation. Higher levels of variation among samples for one locus suggest possible paralogy. Overall, this method provides an efficient means of isolating predominantly single-copy loci from both low and high-copy gene families, providing ample nuclear variation for reconstructing species-level phylogenies in non-model taxa.

AGRIS and AtRegNet. a platform to link cis-regulatory elements and transcription factors into regulatory networks

Gene regulatory pathways converge at the level of transcription, where interactions among regulatory genes and between regulators and target genes result in the establishment of spatiotemporal patterns of gene expression. The growing identification of direct target genes for key transcription factors (TFs) through traditional and high-throughput experimental approaches has facilitated the elucidation of regulatory networks at the genome level. To integrate this information into a Web-based knowledgebase, we have developed the Arabidopsis Gene Regulatory Information Server (AGRIS). AGRIS, which contains all Arabidopsis (Arabidopsis thaliana) promoter sequences, TFs, and their target genes and functions, provides the scientific community with a platform to establish regulatory networks. AGRIS currently houses three linked databases: AtcisDB (Arabidopsis thaliana cis-regulatory database), AtTFDB (Arabidopsis thaliana transcription factor database), and AtRegNet (Arabidopsis thaliana regulatory network). AtTFDB contains 1,690 Arabidopsis TFs and their sequences (protein and DNA) grouped into 50 (October 2005) families with information on available mutants in the corresponding genes. AtcisDB consists of 25,806 (September 2005) promoter sequences of annotated Arabidopsis genes with a description of putative cis-regulatory elements. AtRegNet links, in direct interactions, several hundred genes with the TFs that control their expression. The current release of AtRegNet contains a total of 187 (September 2005) direct targets for 66 TFs. AGRIS can be accessed at http://Arabidopsis.med.ohio-state.edu.