Complete Eukaryotic Genomes Research Articles

Analysis of evolution of exon-intron structure of eukaryotic genes

The availability of multiple, complete eukaryotic genome sequences allows one to address many fundamental evolutionary questions on genome scale. One such important, long-standing problem is evolution of exon-intron structure of eukaryotic genes. Analysis of orthologous genes from completely sequenced genomes revealed numerous shared intron positions in orthologous genes from animals and plants and even between animals, plants and protists. The data on shared and lineage-specific intron positions were used as the starting point for evolutionary reconstruction with parsimony and maximum-likelihood approaches. Parsimony methods produce reconstructions with intron-rich ancestors but also infer lineage-specific, in many cases, high levels of intron loss and gain. Different probabilistic models gave opposite results, apparently depending on model parameters and assumptions, from domination of intron loss, with extremely intron-rich ancestors, to dramatic excess of gains, to the point of denying any true conservation of intron positions among deep eukaryotic lineages. Development of models with adequate, realistic parameters and assumptions seems to be crucial for obtaining more definitive estimates of intron gain and loss in different eukaryotic lineages. Many shared intron positions were detected in ancestral eukaryotic paralogues which evolved by duplication prior to the divergence of extant eukaryotic lineages. These findings indicate that numerous introns were present in eukaryotic genes already at the earliest stages of evolution of eukaryotes and are compatible with the hypothesis that the original, catastrophic intron invasion accompanied the emergence of the eukaryotic cells. Comparison of various features of old and younger introns starts shedding light on probable mechanisms of intron insertion, indicating that propagation of old introns is unlikely to be a major mechanism for origin of new ones. The existence and structure of ancestral protosplice sites were addressed by examining the context of introns inserted within codons that encode amino acids conserved in all eukaryotes and, accordingly, are not subject to selection for splicing efficiency. It was shown that introns indeed predominantly insert into or are fixed in specific protosplice sites which have the consensus sequence (A/C)AG|Gt.

Predicting Protein Function With Hierarchical Phylogenetic Profiles: the Gene3D Phylo-Tuner Method Applied to Eukaryotic Genomes

“Phylogenetic profiling” is based on the hypothesis that during evolution functionally or physically interacting genes are likely to be inherited or eliminated in a codependent manner. Creating presence–absence profiles of orthologous genes is now a common and powerful way of identifying functionally associated genes. In this approach, correctly determining orthology, as a means of identifying functional equivalence between two genes, is a critical and nontrivial step and largely explains why previous work in this area has mainly focused on using presence–absence profiles in prokaryotic species. Here, we demonstrate that eukaryotic genomes have a high proportion of multigene families whose phylogenetic profile distributions are poor in presence–absence information content. This feature makes them prone to orthology mis-assignment and unsuited to standard profile-based prediction methods. Using CATH structural domain assignments from the Gene3D database for 13 complete eukaryotic genomes, we have developed a novel modification of the phylogenetic profiling method that uses genome copy number of each domain superfamily to predict functional relationships. In our approach, superfamilies are subclustered at ten levels of sequence identity—from 30% to 100%—and phylogenetic profiles built at each level. All the profiles are compared using normalised Euclidean distances to identify those with correlated changes in their domain copy number. We demonstrate that two protein families will “auto-tune” with strong co-evolutionary signals when their profiles are compared at the similarity levels that capture their functional relationship. Our method finds functional relationships that are not detectable by the conventional presence–absence profile comparisons, and it does not require a priori any fixed criteria to define orthologous genes.

The evolutionary dynamics of eukaryotic gene order.

In eukaryotes, unlike in bacteria, gene order has typically been assumed to be random. However, the first statistically rigorous analyses of complete genomes, together with the availability of abundant gene-expression data, have forced a paradigm shift: in every complete eukaryotic genome that has been analysed so far, gene order is not random. It seems that genes that have similar and/or coordinated expression are often clustered. Here, we review this evidence and ask how such clusters evolve and how this relates to mechanisms that control gene expression.

A comprehensive evolutionary classification of proteins encoded in complete eukaryotic genomes.

Sequencing the genomes of multiple, taxonomically diverse eukaryotes enables in-depth comparative-genomic analysis which is expected to help in reconstructing ancestral eukaryotic genomes and major events in eukaryotic evolution and in making functional predictions for currently uncharacterized conserved genes. We examined functional and evolutionary patterns in the recently constructed set of 5,873 clusters of predicted orthologs (eukaryotic orthologous groups or KOGs) from seven eukaryotic genomes: Caenorhabditis elegans, Drosophila melanogaster, Homo sapiens, Arabidopsis thaliana, Saccharomyces cerevisiae, Schizosaccharomyces pombe and Encephalitozoon cuniculi. Conservation of KOGs through the phyletic range of eukaryotes strongly correlates with their functions and with the effect of gene knockout on the organism's viability. The approximately 40% of KOGs that are represented in six or seven species are enriched in proteins responsible for housekeeping functions, particularly translation and RNA processing. These conserved KOGs are often essential for survival and might approximate the minimal set of essential eukaryotic genes. The 131 single-member, pan-eukaryotic KOGs we identified were examined in detail. For around 20 that remained uncharacterized, functions were predicted by in-depth sequence analysis and examination of genomic context. Nearly all these proteins are subunits of known or predicted multiprotein complexes, in agreement with the balance hypothesis of evolution of gene copy number. Other KOGs show a variety of phyletic patterns, which points to major contributions of lineage-specific gene loss and the 'invention' of genes new to eukaryotic evolution. Examination of the sets of KOGs lost in individual lineages reveals co-elimination of functionally connected genes. Parsimonious scenarios of eukaryotic genome evolution and gene sets for ancestral eukaryotic forms were reconstructed. The gene set of the last common ancestor of the crown group consists of 3,413 KOGs and largely includes proteins involved in genome replication and expression, and central metabolism. Only 44% of the KOGs, mostly from the reconstructed gene set of the last common ancestor of the crown group, have detectable homologs in prokaryotes; the remainder apparently evolved via duplication with divergence and invention of new genes. The KOG analysis reveals a conserved core of largely essential eukaryotic genes as well as major diversification and innovation associated with evolution of eukaryotic genomes. The results provide quantitative support for major trends of eukaryotic evolution noticed previously at the qualitative level and a basis for detailed reconstruction of evolution of eukaryotic genomes and biology of ancestral forms.

The genome reduction hypothesis and the phylogeny of eukaryotes

The first steps in eukaryotic evolution appear difficult to retrace despite the availability of an increasing amount of data. Current molecular phylogenies suggest that the eukaryotic tree would be better represented as a bush of major lineages whose order of emerge is poorly resolved. Such lack of resolution is often explained by a radiation event that would have left very little ancient signal in eukaryotic molecular markers. We suggest a complementary genomic approach that might help tackling this major issue. It rests on a hypothesis, the genome reduction hypothesis (GRH), suggesting that the divergence of major eukaryotic lineages might have been coupled with independent genomic reduction events, starting from a large and partially redundant chimerical genome. Thus, significant and coherent patterns of shared ancestral gene losses between major eukaryotic lineages might help polarizing the most basal nodes in the eukaryotic phylogeny. We propose a test for the GRH that exploits the increasing availability of complete eukaryotic genomes in public databases.

Gene loss, protein sequence divergence, gene dispensability, expression level, and interactivity are correlated in eukaryotic evolution.

Lineage-specific gene loss, to a large extent, accounts for the differences in gene repertoires between genomes, particularly among eukaryotes. We derived a parsimonious scenario of gene losses for eukaryotic orthologous groups (KOGs) from seven complete eukaryotic genomes. The scenario involves substantial gene loss in fungi, nematodes, and insects. Based on this evolutionary scenario and estimates of the divergence times between major eukaryotic phyla, we introduce a numerical measure, the propensity for gene loss (PGL). We explore the connection among the propensity of a gene to be lost in evolution (PGL value), protein sequence divergence, the effect of gene knockout on fitness, the number of protein-protein interactions, and expression level for the genes in KOGs. Significant correlations between PGL and each of these variables were detected. Genes that have a lower propensity to be lost in eukaryotic evolution accumulate fewer substitutions in their protein sequences and tend to be essential for the organism viability, tend to be highly expressed, and have many interaction partners. The dependence between PGL and gene dispensability and interactivity is much stronger than that for sequence evolution rate. Thus, propensity of a gene to be lost during evolution seems to be a direct reflection of its biological importance.

The Proteome Analysis database: a tool for the in silico analysis of whole proteomes.

The Proteome Analysis database (http://www.ebi.ac.uk/proteome/) has been developed by the Sequence Database Group at EBI utilizing existing resources and providing comparative analysis of the predicted protein coding sequences of the complete genomes of bacteria, archeae and eukaryotes. Three main projects are used, InterPro, CluSTr and GO Slim, to give an overview on families, domains, sites, and functions of the proteins from each of the complete genomes. Complete proteome analysis is available for a total of 89 proteome sets. A specifically designed application enables InterPro proteome comparisons for any one proteome against any other one or more of the proteomes in the database.

A Novel Family of Calmodulin-binding Transcription Activators in Multicellular Organisms

Screening of cDNA expression libraries derived from plants exposed to stress, with 35S-labeled recombinant calmodulin as a probe, revealed a new family of proteins containing a transcription activation domain and two types of DNA-binding domains designated the CG-1 domain and the transcription factor immunoglobulin domain, ankyrin repeats, and a varying number of IQ calmodulin-binding motifs. Based on domain organization and amino acid sequence comparisons, similar proteins, with the same domain organization, were identified in the genomes of other multicellular organisms including human, Drosophila, and Caenorhabditis, whereas none were found in the complete genomes of single cell eukaryotes and prokaryotes. This family of proteins was designated calmodulin-binding transcription activators (CAMTAs). Arabidopsis thaliana contains six CAMTA genes (AtCAMTA1-AtCAMTA6). The transcription activation domain of AtCAMTA1 was mapped by testing a series of protein fusions with the DNA-binding domain of the bacterial LexA transcription factor and two reporter genes fused to LexA recognition sequences in yeast cells. Two human proteins designated HsCAMTA1 and HsCAMTA2 were also shown to activate transcription in yeast using the same reporter system. Subcellular fractionation of Arabidopsis tissues revealed the presence of CAMTAs predominantly in the nucleus. Calmodulin binding assays identified a region of 25 amino acids capable of binding calmodulin with high affinity (K(d) = 1.2 nm) in the presence of calcium. We suggest that CAMTAs comprise a conserved family of transcription factors in a wide range of multicellular eukaryotes, which possibly respond to calcium signaling by direct binding of calmodulin.

Segmentation of long genomic sequences into domains with homogeneous composition with BASIO software.

We present a software system BASIO that allows one to segment a sequence into regions with homogeneous nucleotide composition at a desired length scale. The system can work with arbitrary alphabet and therefore can be applied to various (e.g. protein) sequences. Several sequences of complete genomes of eukaryotes are used to demonstrate the efficiency of the software. The BASIO suite is available for non-commercial users free of charge as a set of executables and accompanying segmentation scenarios from http://www.imb.ac.ru/compbio/basio. To obtain the source code, contact the authors.

Exon structure conservation despite low sequence similarity: a relic of dramatic events in evolution?

The evolutionary significance of introns remains a mystery. The current availability of several complete eukaryotic genomes permits new studies to probe the possible function of these peculiar genomic features. Here we investigate the degree to which gene structure (intron position, phase and length) is conserved between homologous protein domains. We find that for certain extracellular-signalling and nuclear domains, gene structures are similar even when protein sequence similarity is low or not significant and sequences can only be aligned with a knowledge of protein tertiary structure. In contrast, other domains, including most intracellular signalling modules, show little gene structure conservation. Intriguingly, many domains with conserved gene structures, such as cytokines, are involved in similar biological processes, such as the immune response. This suggests that gene structure conservation may be a record of key events in evolution, such as the origin of the vertebrate immune system or the duplication of nuclear receptors in nematodes. The results suggest ways to detect new and potentially very remote homologues, and to construct phylogenies for proteins with limited sequence similarity.

Transposable elements in the Animal Kingdom

Transposable elements (TEs) are commonly thought to be of universal occurrence in eukaryotes. Analysis of complete higher eukaryotic genomes confirms TE status as substantial genome components and provides insights into their role in shaping the genome structure of extant eukaryotes. This review addresses several recently investigated problems in transposon biology, including potential roles of promoter organization in transposon function and evolution, the ubiquity of TEs in numerous phyla of the animal kingdom, and possible connections between transposon content and the mode of reproduction.

Molecular fungal biology

Preface 1. Fungi - important organisms in history and today M. Wedde, M. Jacobs and U. Stahl 2. Fungal phylogeny M. L. Berbee and J. W. Taylor 3. The yeast Saccharomyces cerevisiae: insights from the first complete eukaryotic genome sequence H. Feldman 4. Interactions between pathway-specific and global genetic regulation and the control of pathway flux A. R. Hawkins, K. A. Wheeler, L. J. Levett, G. H. Newton and H. K. Lamb 5. Hyphal cell biology G. Robson 6. Asexual conidiation T. H. Adams and J. K. Wieser 7. Fungal cell division C. Pitt and J. Doonan 8. Sexual development of higher fungi S.-W. Chiu and D. Moore 9. Lignocellulose breakdown and utilization by fungi M. Penttila and M. Saloheimo 10. Plant disease caused by fungi - phytopathogenicity P. Bowyer 11. Animal pathogenicity S. Hosking 12. Biotechnology of filamentous fungi: applications of molecular biology D. B. Archer Index.

The COG database: a tool for genome-scale analysis of protein functions and evolution.

Rational classification of proteins encoded in sequenced genomes is critical for making the genome sequences maximally useful for functional and evolutionary studies. The database of Clusters of Orthologous Groups of proteins (COGs) is an attempt on a phylogenetic classification of the proteins encoded in 21 complete genomes of bacteria, archaea and eukaryotes (http://www. ncbi.nlm. nih.gov/COG). The COGs were constructed by applying the criterion of consistency of genome-specific best hits to the results of an exhaustive comparison of all protein sequences from these genomes. The database comprises 2091 COGs that include 56-83% of the gene products from each of the complete bacterial and archaeal genomes and approximately 35% of those from the yeast Saccharomyces cerevisiae genome. The COG database is accompanied by the COGNITOR program that is used to fit new proteins into the COGs and can be applied to functional and phylogenetic annotation of newly sequenced genomes.

Why genome analysis?

The field of comparative genomics is coming of age at an almost frightening speed. The first complete genome sequence of a cellular life form – the parasitic bacterium Haemophilus influenzae – was reported in 1995. Now, only 3.5 years later, we have about 20 complete genomes of bacteria, archaea and eukaryotes, and many more are in the pipeline. These genomes are changing the landscape of today's biology. Indeed, these days, almost all research in molecular biology, genetics and biochemistry, such as any functional study on a particular enzyme or cellular system, can, and probably should, involve a comparative-genomic aspect. That is, one cannot help asking: is my enzyme or system encoded in all genomes? Or only in a particular subset, say, in eukaryotes and archaea, but not in bacteria? Or is this system, with its distinct and elaborate design, present only in animals but its components can be found in yeast, archaea and bacteria where they perform other, perhaps still unknown functions? If I show that my enzyme is essential in the given organism, but is missing in others, what alternative enzyme(s) or strategies are they using? Can we predict an alternative solution from genome comparisons? Needless to say, one can obtain satisfactory answers to questions like this only by carefully comparing complete genomes. Indeed, it is hardly possible to ascertain that a particular gene is not present in the given genome unless the entire genome is available. In an ideal world, biologists should be able to address these and similar questions by querying an ‘intellectual’ database of genome comparisons. The reality is, of course, such that, although considerable systematic effort in genome comparison has been made by a number of research groups, not only does the dream query system not exist, but no one seems to have a good idea what it should look like. The absence of such a coherent system is a major problem but… it is also a research opportunity. A researcher equipped with an Internet-connected computer, a working knowledge of the sequence analysis techniques, and most importantly, a sustained interest in comparative genomics can make observations that, in many cases, are of interest and alue to a broad scientific community.

Genome-wide bioinformatic and molecular analysis of introns in Saccharomyces cerevisiae.

Introns have typically been discovered in an ad hoc fashion: introns are found as a gene is characterized for other reasons. As complete eukaryotic genome sequences become available, better methods for predicting RNA processing signals in raw sequence will be necessary in order to discover genes and predict their expression. Here we present a catalog of 228 yeast introns, arrived at through a combination of bioinformatic and molecular analysis. Introns annotated in the Saccharomyces Genome Database (SGD) were evaluated, questionable introns were removed after failing a test for splicing in vivo, and known introns absent from the SGD annotation were added. A novel branchpoint sequence, AAUUAAC, was identified within an annotated intron that lacks a six-of-seven match to the highly conserved branchpoint consensus UACUAAC. Analysis of the database corroborates many conclusions about pre-mRNA substrate requirements for splicing derived from experimental studies, but indicates that splicing in yeast may not be as rigidly determined by splice-site conservation as had previously been thought. Using this database and a molecular technique that directly displays the lariat intron products of spliced transcripts (intron display), we suggest that the current set of 228 introns is still not complete, and that additional intron-containing genes remain to be discovered in yeast. The database can be accessed at http://www.cse.ucsc.edu/research/compbi o/yeast_introns.html.

The translational signal database, TransTerm, is now a relational database.

TransTerm-97 contains more than 97 500 non-redundant coding-sequence initiation and termination contexts compiled from GenBank, release 101 (15-June-1997). In addition, several coding sequence parameters are available: coding sequence length, Nc, GC3, and, when it is computable, codon adaptation index (CAI). Codon usage tables and summaries of start and stop codon contexts are also included. The information covers more than 325 species and organelles, including seven complete bacterial genomes and one complete eukaryotic genome. To promote research in translational control of protein synthesis, TransTerm has been converted into a relational database to ease the process of making queries. The relational database manager, Postgresql, gives access to the database using SQL (Structured Query Language). A World Wide Web interface using forms is being completed to allow the casual user access to the database. Extensions are planned to include the full 5'-UTR, full coding sequence and 3'-UTR. TransTerm-97 is available on the World Wide Web at:http://biochem. otago.ac.nz:800/Transterm/homepage.html

Effort To Map And Sequence The Human Genome Makes Significant Progress

The Human Genome Project, an international research effort to map and sequence the genomes of humans and selected model organisms, is making significant progress toward its goals four years into its projected 15-year life. A detailed human genetic linkage map has been developed ahead of time. A physical map, consisting of overlapping pieces of DNA, is only slightly behind schedule. Base-by-base sequencing of the human genome is lagging, but sequencing of model organisms is moving along very well, with the first complete eukaryotic genome likely to be completed within two years. Human Genome Project sponsorship of a map that would show the location of expressed human genes is still in the planning stage. However, such maps have been and are being produced privately on a large scale—a state of affairs that has stirred up considerable controversy about whether the market for such data is being cornered by proprietary interests. The Human Genome Project is living ...

Two-dimensional DNA gel electrophoresis as a method for analysis of eukaryotic genome structure: Evaluation using Tetrahymena thermophila DNA

There is growing interest in mapping and analyzing complete eukaryotic genomes. Yee and Inouye (in Experimental Manipulation of Gene Expression, pp. 279–290, Academic Press, New York) demonstrated that bacterial chromosomes can be resolved into interpretable patterns of DNA fragments by means of restriction enzyme digestion and electrophoresis in two dimensions. We have begun to explore applications of this procedure to analysis of eukaryotic genomes, which are far more complex. Tetrahymena thermophila was selected as a model organism because its genome is small, roughly equivalent to that of a single human chromosome. In addition, each Tetrahymena cell contains two nuclei which differ in sequence composition and methylation. Our results demonstrate that the Tetrahymena genome can be resolved into complex patterns of fragments in two dimensions. Hybridization to Southern blots of these gels with a multiply repeated sequence probe yielded analyzable patterns of a subset of the genome. The blots reveal alterations in genome structure due to methylation and rearrangement. Future extensions of the method are discussed.