PlanTAPDB, a phylogeny-based resource of plant transcription-associated proteins

IT3F: a web-based tool for functional analysis of transcription factors in plants

A web-based tool, the Interspecies Transcription Factor Function Finder (IT3F), has been developed to display both evolutionary gene relationships and expression data for plant transcription factors, focussing primarily on the R2R3MYB gene subfamily for proof of concept. The graphical display of information allows users to make direct comparisons between structurally related genes and to identify those genes that are potentially orthologous, thereby assisting with their understanding of gene function. A key feature of the website is the provision of an interrogative phylogenetic tree that allows submission of new sequences corresponding to a transcription factor family or subfamily and maps their relative positions to the products of other genes on an 'existing' tree containing proteins encoded by Arabidopsis and rice genes, along with key proteins encoded by genes from other species that have been characterised functionally. In addition, a feature to select clusters of related sequences has been developed so that more detailed phylogenetic analysis can be performed to highlight potential orthologous and paralogous genes within related clusters. Arabidopsis genes that reside on duplicated regions of the genome are indicated on the tree, providing further information for interpreting gene function. An additional feature of the website allows a selected number of key Arabidopsis and rice microarray experiments to be visualised alongside the tree as a tabulated heat map of expression intensity values. Through this display, it is possible to observe relative expression levels across a whole gene family and the extent to which the expression of closely related genes within subgroups has altered since their ancestral divergence. The website is available at http://jicbio.nbi.ac.uk/IT3F/.

The role of bZIP transcription factors in green plant evolution: adaptive features emerging from four founder genes

Background

Transcription factors of the basic leucine zipper (bZIP) family control important processes in all eukaryotes. In plants, bZIPs are regulators of many central developmental and physiological processes including photomorphogenesis, leaf and seed formation, energy homeostasis, and abiotic and biotic stress responses. Here we performed a comprehensive phylogenetic analysis of bZIP genes from algae, mosses, ferns, gymnosperms and angiosperms.

Methodology/Principal Findings

We identified 13 groups of bZIP homologues in angiosperms, three more than known before, that represent 34 Possible Groups of Orthologues (PoGOs). The 34 PoGOs may correspond to the complete set of ancestral angiosperm bZIP genes that participated in the diversification of flowering plants. Homologous genes dedicated to seed-related processes and ABA-mediated stress responses originated in the common ancestor of seed plants, and three groups of homologues emerged in the angiosperm lineage, of which one group plays a role in optimizing the use of energy.

Conclusions/Significance

Our data suggest that the ancestor of green plants possessed four bZIP genes functionally involved in oxidative stress and unfolded protein responses that are bZIP-mediated processes in all eukaryotes, but also in light-dependent regulations. The four founder genes amplified and diverged significantly, generating traits that benefited the colonization of new environments.

Tobacco transcription factors: novel insights into transcriptional regulation in the Solanaceae

Tobacco (Nicotiana tabacum) is a member of the Solanaceae, one of the agronomically most important groups of flowering plants. We have performed an in silico analysis of 1.15 million gene-space sequence reads from the tobacco nuclear genome and report the detailed analysis of more than 2,500 tobacco transcription factors (TFs). The tobacco genome contains at least one member of each of the 64 well-characterized TF families identified in sequenced vascular plant genomes, indicating that evolution of the Solanaceae was not associated with the gain or loss of TF families. However, we found notable differences between tobacco and non-Solanaceae species in TF family size and evidence for both tobacco- and Solanaceae-specific subfamily expansions. Compared with TF families from sequenced plant genomes, tobacco has a higher proportion of ERF/AP2, C2H2 zinc finger, homeodomain, GRF, TCP, zinc finger homeodomain, BES, and STERILE APETALA (SAP) genes and novel subfamilies of BES, C2H2 zinc finger, SAP, and NAC genes. The novel NAC subfamily, termed TNACS, appears restricted to the Solanaceae, as they are absent from currently sequenced plant genomes but present in tomato (Solanum lycopersicum), pepper (Capsicum annuum), and potato (Solanum tuberosum). They constitute approximately 25% of NAC genes in tobacco. Based on our phylogenetic studies, we predict that many of the more than 50 tobacco group IX ERF genes are involved in jasmonate responses. Consistent with this, over two-thirds of group IX ERF genes tested showed increased mRNA levels following jasmonate treatment. Our data are a major resource for the Solanaceae and fill a void in studies of TF families across the plant kingdom.

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

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

Conifer WOX-related homeodomain transcription factors, developmental consideration and expression dynamic of WOX2 during Picea abies somatic embryogenesis

In angiosperms, the WOX family of transcription factors has important functions in meristem regulation and in control of the partitioning of developing embryos into functional domains. In this study, a putative WOX2 homologous gene was isolated from Picea abies, and its expression pattern during somatic embryo development was followed using real-time quantitative reverse transcription polymerase chain reaction (qRT-PCR). We used strategies of both absolute and relative quantification of gene expression, and benefits and disadvantages of the two methods are presented and discussed. During embryogenesis, PaWOX2 expression was highest at the earliest stages of development, but low levels were also detected in seedling tissues. No PaWOX2 expression was detected in a non-embryogenic cell culture, indicating that PaWOX2 plays a fundamental role during early somatic embryo development, and can be used as a possible marker for embryogenic potential. Additional results show that conifers, like angiosperms, contain a large number of WOX-related genes, many of them expressed during embryo development. In phylogenetic analysis based on the deduced homeodomain of retrieved pine and spruce EST sequences, no conifer WUS homolog was found. Neither did we find any homeodomain to cluster with WOX5. Interestingly, a clade including only conifer sequences derived from various tissues was resolved as sister to a Physcomitrella WOX-like gene, suggestive of the early origin of this gene family. Our results thus provide basic information for further studies of the evolution of this gene family and of their function in relation to meristem dynamics and specification of stem cells in gymnosperms.

Sequencing and analysis of the gene-rich space of cowpea

Background

Cowpea, Vigna unguiculata (L.) Walp., is one of the most important food and forage legumes in the semi-arid tropics because of its drought tolerance and ability to grow on poor quality soils. Approximately 80% of cowpea production takes place in the dry savannahs of tropical West and Central Africa, mostly by poor subsistence farmers. Despite its economic and social importance in the developing world, cowpea remains to a large extent an underexploited crop. Among the major goals of cowpea breeding and improvement programs is the stacking of desirable agronomic traits, such as disease and pest resistance and response to abiotic stresses. Implementation of marker-assisted selection and breeding programs is severely limited by a paucity of trait-linked markers and a general lack of information on gene structure and organization. With a nuclear genome size estimated at ~620 Mb, the cowpea genome is an ideal target for reduced representation sequencing.

Results

We report here the sequencing and analysis of the gene-rich, hypomethylated portion of the cowpea genome selectively cloned by methylation filtration (MF) technology. Over 250,000 gene-space sequence reads (GSRs) with an average length of 610 bp were generated, yielding ~160 Mb of sequence information. The GSRs were assembled, annotated by BLAST homology searches of four public protein annotation databases and four plant proteomes (A. thaliana, M. truncatula, O. sativa, and P. trichocarpa), and analyzed using various domain and gene modeling tools. A total of 41,260 GSR assemblies and singletons were annotated, of which 19,786 have unique GenBank accession numbers. Within the GSR dataset, 29% of the sequences were annotated using the Arabidopsis Gene Ontology (GO) with the largest categories of assigned function being catalytic activity and metabolic processes, groups that include the majority of cellular enzymes and components of amino acid, carbohydrate and lipid metabolism. A total of 5,888 GSRs had homology to genes encoding transcription factors (TFs) and transcription associated factors (TAFs) representing about 5% of the total annotated sequences in the dataset. Sixty-two (62) of the 64 well-characterized plant transcription factor (TF) gene families are represented in the cowpea GSRs, and these families are of similar size and phylogenetic organization to those characterized in other plants. The cowpea GSRs also provides a rich source of genes involved in photoperiodic control, symbiosis, and defense-related responses. Comparisons to available databases revealed that about 74% of cowpea ESTs and 70% of all legume ESTs were represented in the GSR dataset. As approximately 12% of all GSRs contain an identifiable simple-sequence repeat, the dataset is a powerful resource for the design of microsatellite markers.

Conclusion

The availability of extensive publicly available genomic data for cowpea, a non-model legume with significant importance in the developing world, represents a significant step forward in legume research. Not only does the gene space sequence enable the detailed analysis of gene structure, gene family organization and phylogenetic relationships within cowpea, but it also facilitates the characterization of syntenic relationships with other cultivated and model legumes, and will contribute to determining patterns of chromosomal evolution in the Leguminosae. The micro and macrosyntenic relationships detected between cowpea and other cultivated and model legumes should simplify the identification of informative markers for marker-assisted trait selection and map-based gene isolation necessary for cowpea improvement.

TOBFAC: the database of tobacco transcription factors

BACKGROUND

Regulation of gene expression at the level of transcription is a major control point in many biological processes. Transcription factors (TFs) can activate and/or repress the transcriptional rate of target genes and vascular plant genomes devote approximately 7% of their coding capacity to TFs. Global analysis of TFs has only been performed for three complete higher plant genomes - Arabidopsis (Arabidopsis thaliana), poplar (Populus trichocarpa) and rice (Oryza sativa). Presently, no large-scale analysis of TFs has been made from a member of the Solanaceae, one of the most important families of vascular plants. To fill this void, we have analysed tobacco (Nicotiana tabacum) TFs using a dataset of 1,159,022 gene-space sequence reads (GSRs) obtained by methylation filtering of the tobacco genome. An analytical pipeline was developed to isolate TF sequences from the GSR data set. This involved multiple (typically 10-15) independent searches with different versions of the TF family-defining domain(s) (normally the DNA-binding domain) followed by assembly into contigs and verification. Our analysis revealed that tobacco contains a minimum of 2,513 TFs representing all of the 64 well-characterised plant TF families. The number of TFs in tobacco is higher than previously reported for Arabidopsis and rice.

RESULTS

TOBFAC: the database of tobacco transcription factors, is an integrative database that provides a portal to sequence and phylogeny data for the identified TFs, together with a large quantity of other data concerning TFs in tobacco. The database contains an individual page dedicated to each of the 64 TF families. These contain background information, domain architecture via Pfam links, a list of all sequences and an assessment of the minimum number of TFs in this family in tobacco. Downloadable phylogenetic trees of the major families are provided along with detailed information on the bioinformatic pipeline that was used to find all family members. TOBFAC also contains EST data, a list of published tobacco TFs and a list of papers concerning tobacco TFs. The sequences and annotation data are stored in relational tables using a PostgrelSQL relational database management system. The data processing and analysis pipelines used the Perl programming language. The web interface was implemented in JavaScript and Perl CGI running on an Apache web server. The computationally intensive data processing and analysis pipelines were run on an Apple XServe cluster with more than 20 nodes.

CONCLUSION

TOBFAC is an expandable knowledgebase of tobacco TFs with data currently available for over 2,513 TFs from 64 gene families. TOBFAC integrates available sequence information, phylogenetic analysis, and EST data with published reports on tobacco TF function. The database provides a major resource for the study of gene expression in tobacco and the Solanaceae and helps to fill a current gap in studies of TF families across the plant kingdom. TOBFAC is publicly accessible at http://compsysbio.achs.virginia.edu/tobfac/.

The evolution of plant regulatory networks: what Arabidopsis cannot say for itself

Genetic and molecular analyses in the dicot model plant Arabidopsis thaliana have begun to shed some light on regulatory networks in plants. However, comparisons with other species are necessary to validate networks identified in model species on the evolutionary scale. Many key regulatory proteins are encoded by members of transcription factor gene families. Orthologous genes can be identified by phylogenetic reconstructions based on conserved protein domains and functionally substantiated by gene expression patterns and mutant analyses. Recent comparative analyses of different pathways involved in shoot meristem development reveal not only conservation from basal land plants to angiosperms but also evolutionary freedom for significant adaptations in the course of plant speciation.

PlantTFDB: a comprehensive plant transcription factor database

Transcription factors (TFs) play key roles in controlling gene expression. Systematic identification and annotation of TFs, followed by construction of TF databases may serve as useful resources for studying the function and evolution of transcription factors. We developed a comprehensive plant transcription factor database PlantTFDB (http://planttfdb.cbi.pku.edu.cn), which contains 26 402 TFs predicted from 22 species, including five model organisms with available whole genome sequence and 17 plants with available EST sequences. To provide comprehensive information for those putative TFs, we made extensive annotation at both family and gene levels. A brief introduction and key references were presented for each family. Functional domain information and cross-references to various well-known public databases were available for each identified TF. In addition, we predicted putative orthologs of those TFs among the 22 species. PlantTFDB has a simple interface to allow users to search the database by IDs or free texts, to make sequence similarity search against TFs of all or individual species, and to download TF sequences for local analysis.

Moss bioreactors producing improved biopharmaceuticals

Plants may serve as superior production systems for complex recombinant pharmaceuticals. Current strategies for improving plant-based systems include the development of large-scale production facilities as well as the optimisation of protein modifications. While post-translational modifications of plant proteins generally resemble those of mammalian proteins, certain plant-specific protein-linked sugars are immunogenic in humans, a fact that restricts the use of plants in biopharmaceutical production so far. The moss Physcomitrella patens was developed as a contained tissue culture system for recombinant protein production in photo-bioreactors. By targeted gene replacements, moss strains were created with non-immunogenic humanised glycan patterns. These were proven to be superior to currently used mammalian cell lines in producing antibodies with enhanced effectiveness.

Sulfate assimilation in basal land plants - what does genomic sequencing tell us?

Sulfate assimilation is a pathway providing reduced sulfur for the synthesis of cysteine, methionine, co-enzymes such as iron-sulfur centres, thiamine, lipoic acid, or Coenzyme A, and many secondary metabolites, e.g., glucosinolates or alliins. The pathway is relatively well understood in flowering plants, but very little information exists on sulfate assimilation in basal land plants. Since the finding of a putative 3'-phosphoadenosine 5'-phosphosulfate reductase in PHYSCOMITRELLA PATENS, an enigmatic enzyme thought to exist in fungi and some bacteria only, it has been evident that sulfur metabolism in lower plants may substantially differ from seed plant models. The genomic sequencing of two basal plant species, the Bryophyte PHYSCOMITRELLA PATENS, and the Lycophyte SELAGINELLA MOELLENDORFFII, opens up the possibility to search for differences between lower and higher plants at the genomic level. Here we describe the similarities and differences in the organisation of the sulfate assimilation pathway between basal and advanced land plants derived from genome comparisons of these two species with ARABIDOPSIS THALIANA and ORYZA SATIVA, two seed plants with sequenced genomes. We found differences in the number of genes encoding sulfate transporters, adenosine 5'-phosphosulfate reductase, and sulfite reductase between the lower and higher plants. The consequences for regulation of the pathway and evolution of sulfate assimilation in plants are discussed.

Comparative analysis of the SBP-box gene families in P. patens and seed plants

To come to a better understanding of the evolution and function of the SBP-box transcription factor family in plants, we identified, isolated and characterized 13 of its members from the moss Physcomitrella patens. For the majority of the moss SBP-box genes, clear orthologous relationships with family members of flowering plants could be established by phylogenetic analysis based on the conserved DNA-binding SBP-domain, as well as additional synapomorphic molecular characters. The P. patens SBP-box genes cluster in four separable groups. One of these consists exclusively of moss genes; the three others are shared with family members of Arabidopsis and rice. Besides the family defining DNA-binding SBP-domain, other features can be found conserved between moss and other plant SBP-domain proteins. An AHA-like motif conserved from the unicellular alga Chlamydomonas reinhardtii to flowering plants, was found able to promote transcription in a heterologous yeast system. The conservation of a functional microRNA response element in the mRNA of three of the moss SBP-box genes supports the idea of an ancient origin of microRNA dependent regulation of SBP-box gene family members. As our current knowledge concerning the roles of SBP-box genes in plant development is scarce and the model system P. patens allows targeted mutation, the material we isolated and characterized will be helpful to generate the mutant phenotypes necessary to further elucidate these roles.

Dating the early evolution of plants: detection and molecular clock analyses of orthologs

Orthologs generally are under selective pressure against loss of function, while paralogs usually accumulate mutations and finally die or deviate in terms of function or regulation. Most ortholog detection methods contaminate the resulting datasets with a substantial amount of paralogs. Therefore we aimed to implement a straightforward method that allows the detection of ortholog clusters with a reduced amount of paralogs from completely sequenced genomes. The described cross-species expansion of the reciprocal best BLAST hit method is a time-effective method for ortholog detection, which results in 68% truly orthologous clusters and the procedure specifically enriches single-copy orthologs. The detection of true orthologs can provide a phylogenetic toolkit to better understand evolutionary processes. In a study across six photosynthetic eukaryotes, nuclear genes of putative mitochondrial origin were shown to be over-represented among single copy orthologs. These orthologs are involved in fundamental biological processes like amino acid metabolism or translation. Molecular clock analyses based on this dataset yielded divergence time estimates for the red/green algae (1,142 MYA), green algae/land plant (725 MYA), mosses/seed plant (496 MYA), gymno-/angiosperm (385 MYA) and monocotyledons/core eudicotyledons (301 MYA) divergence times.