Expressed sequence tags from the Yukon ecotype of Thellungiella reveal that gene expression in response to cold, drought and salinity shows little overlap.

Thellungiella salsuginea (also known as T. halophila) is a close relative of Arabidopsis that is very tolerant of drought, freezing, and salinity and may be an appropriate model to identify the molecular mechanisms underlying abiotic stress tolerance in plants. We produced 6578 ESTs, which represented 3628 unique genes (unigenes), from cDNA libraries of cold-, drought-, and salinity-stressed plants from the Yukon ecotype of Thellungiella. Among the unigenes, 94.1% encoded products that were most similar in amino acid sequence to Arabidopsis and 1.5% had no match with a member of the family Brassicaceae. Unigenes from the cold library were more similar to Arabidopsis sequences than either drought- or salinity-induced sequences, indicating that latter responses may be more divergent between Thellungiella and Arabidopsis. Analysis of gene ontology using the best matched Arabidopsis locus showed that the Thellungiella unigenes represented all biological processes and all cellular components, with the highest number of sequences attributed to the chloroplast and mitochondria. Only 140 of the unigenes were found in all three abiotic stress cDNA libraries. Of these common unigenes, 70% have no known function, which demonstrates that Thellungiella can be a rich resource of genetic information about environmental responses. Some of the ESTs in this collection have low sequence similarity with those in Genbank suggesting that they may encode functions that may contribute to Thellungiella's high degree of stress tolerance when compared with Arabidopsis. Moreover, Thellungiella is a closer relative of agriculturally important Brassica spp. than Arabidopsis, which may prove valuable in transferring information to crop improvement programs.

Patterns of bacterial gene movement.

Lateral gene transfer has emerged as an important force in bacterial evolution. A substantial number of genes can be inserted into or deleted from genomes through the process of lateral transfer. In this study, we looked for atypical occurrence of genes among related organisms to detect laterally transferred genes. We have analyzed 50 bacterial complete genomes from nine groups. For each group we use a 16s rRNA phylogeny and a comparison of protein similarity to map gene insertions/deletions onto their species phylogeny. The results reveal that there is poor correlation of genes inserted, deleted, and duplicated with evolutionary branch length. In addition, the numbers of genes inserted, deleted, or duplicated within the same branch are not always correlated with each other. Nor is there any similarity within groups. For example, in the Rhizobiales group, the ratio of insertions to deletions in the evolutionary branch leading to Agrobacterium tumefaciens str. C58 (Cereon) is 0.52, but it is 39.52 for Mesorhizobium loti. Most strikingly, the number of insertions of foreign genes is much larger in the external branches of the trees. These insertions also greatly outnumber the occurrence of deletions, and yet the genome sizes of these bacteria remain roughly constant. This indicates that many of the insertions are specific to each organism and are lost before related species can evolve. Simulations of the process of insertion and deletion, tailored to each phylogeny, support this conclusion.

The closest BLAST hit is often not the nearest neighbor.

It is well known that basing phylogenetic reconstructions on uncorrected genetic distances can lead to errors in their reconstruction. Nevertheless, it is often common practice to report simply the most similar BLAST (Altschul et al. 1997) hit in genomic reports that discuss many genes (Ruepp et al. 2000; Freiberg et al. 1997). This is because BLAST hits can provide a rapid, efficient, and concise analysis of many genes at once. These hits are often interpreted to imply that the gene is most closely related to the gene or protein in the databases that returned the closest BLAST hit. Though these two may coincide, for many genes, particularly genes with few homologs, they may not be the same. There are a number of circumstances that can account for such limitations in accuracy (Eisen 2000). We stress here that genes appearing to be the most similar based on BLAST hits are often not each others closest relative phylogenetically. The extent to which this occurs depends on the availability of close relatives present in the databases. As an example we have chosen the analysis of the genomes of a crenarcheaota species Aeropyrum pernix, an organism with few close relatives fully sequenced, and Escherichia coli, an organism whose closest relative, Salmonella typhimurium, is completely sequenced.

Codon bias and base composition are poor indicators of horizontally transferred genes.

Horizontal gene transfer is now recognized as an important mechanism of evolution. Several methods to detect horizontally transferred genes have been suggested. These methods are based on either nucleotide composition or the failure to find a similar gene in closely related species. Genes that evolve vertically between closely related species can be divided into those that retain homologous chromosomal positions (positional orthologs) and those that do not. By comparing open reading frames in the Escherichia coli and Salmonella typhi genomes, we identified 2,728 positional orthologs since these species split 100 MYA. A group of 1,144 novel E. coli genes were unusually diverged from their S. typhi counterparts. These novel genes included those that had been horizontally transferred into E. coli, as well as members of gene pairs that had been rearranged or deleted. Positional orthologs were used to investigate compositional methods of identifying horizontally transferred genes. A large number of E. coli genes with normal nucleotide composition have no apparent ortholog in S. typhi, and many genes of atypical composition do, in fact, have positional orthologs. A phylogenetic approach was employed to confirm selected examples of horizontal transmission among the novel groups of genes. Our analysis of 80 E. coli genes determined that a number of genes previously classified as horizontally transferred based on base composition and codon bias were native, and genes previously classified as native appeared to be horizontally transferred. Hence, atypical nucleotide composition alone is not a reliable indicator of horizontal transmission.

Evolution of simple sequence in proteins.

The proteins of Saccharomyces cerevisiae contain a high proportion of low-complexity, simple sequences. These are protein segments composed almost exclusively or largely of a single repetitive amino acid polymer and are the most commonly shared feature between proteins. We have examined a survey of other species to determine how widespread this phenomenon might be. This was done by comparing how frequently segments from one protein are present in other proteins. Any recently evolutionarily related proteins were excluded. It was found that the most commonly shared features of eukaryotic proteins were repetitive but that prokaryotes did not contain such shared, extensively redundant repeats. The proportion of eukaryotic proteins that contain a significantly repetitive fraction changes dramatically from species to species. In addition the individual amino acids present in these repeats change between species. This suggests that the primary sequence of the repeats may not be important for their function. Further tests of the yeast repeats confirmed that these repeats evolve more quickly than the remainder of the protein sequence within which they are embedded. These results show that these rapid evolving, simple sequence repeats are in fact the most commonly shared pattern between all of the genomic proteins of eukaryotes.

Enzyme evolution explained (sort of).

Sites in proteins evolve at markedly different rates; some are highly conserved, others change rapidly. We have developed a maximum likelihood method to identify regions of a protein that evolve rapidly or slowly relative to the remaining structure. We also show that solvent accessibility and distance from the catalytic site are major determinants of evolutionary rate in eubacterial isocitrate dehydrogenases. These two variables account for most of the rate heterogeneity not ascribable to stochastic effects.

Simple sequence is abundant in eukaryotic proteins.

All proteins of Saccharomyces cerevisiae have been compared to determine how frequently segments from one protein are present in other proteins. Proteins that are recently evolutionarily related were excluded. The most frequently present protein segments are long, tandem repetitions of a single amino acid. For some of these segments, up to 14% of all proteins in the genome were found to have similar peptides within them. These peptide segments may not be functional protein domains. Although they are the most common shared feature of yeast proteins, their ubiquity and simplicity argue that their probable function may be to simply serve as spacers between other protein motifs.

The mosaic nature of the eukaryotic nucleus.

The phylogenies for each of the protein-coding genes from the Methanococcus jannaschii genome were surveyed to determine the history of the major groups of life. For each gene, homologous sequences from other archaea, eucarya, and Gram-positive and Gram-negative bacteria were collected and aligned, and a phylogeny was reconstructed with a maximum-likelihood algorithm. The majority of significant phylogenies favor the eucarya and the archaca as sister groups. A smaller, but still substantial, portion of these significant phylogenies favor an eucarya/Gram-negative clade. These results indicate that support for the early history of life is not unequivocal. A chimeric origin of eukaryotes or an ancient, massive horizontal transfer of genes from Gram-negative bacteria to eucarya can explain many of the observed phylogenies.

The structural basis of molecular adaptation.

The study of molecular adaptation has long been fraught with difficulties, not the least of which is identifying out of hundreds of amino acid replacements those few directly responsible for major adaptations. Six studies are used to illustrate how phylogenies, site-directed mutagenesis, and a knowledge of protein structure combine to provide much deeper insights into the adaptive process than has hitherto been possible. Ancient genes can be reconstructed, and the phenotypes can be compared to modern proteins. Out of hundreds of amino acid replacements accumulated over billions of years those few responsible for discriminating between alternative substrates are identified. An amino acid replacement of modest effect at the molecular level causes a dramatic expansion in an ecological niche. These and other topics are creating the emerging field of "paleomolecular biochemistry."

Genome structure and gene content in protist mitochondrial DNAs.

Although the collection of completely sequenced mitochondrial genomes is expanding rapidly, only recently has a phylogenetically broad representation of mtDNA sequences from protists (mostly unicellular eukaryotes) become available. This review surveys the 23 complete protist mtDNA sequences that have been determined to date, commenting on such aspects as mitochondrial genome structure, gene content, ribosomal RNA, introns, transfer RNAs and the genetic code and phylogenetic implications. We also illustrate the utility of a comparative genomics approach to gene identification by providing evidence that orfB in plant and protist mtDNAs is the homolog of atp8 , the gene in animal and fungal mtDNA that encodes subunit 8 of the F0portion of mitochondrial ATP synthase. Although several protist mtDNAs, like those of animals and most fungi, are seen to be highly derived, others appear to be have retained a number of features of the ancestral, proto-mitochondrial genome. Some of these ancestral features are also shared with plant mtDNA, although the latter have evidently expanded considerably in size, if not in gene content, in the course of evolution. Comparative analysis of protist mtDNAs is providing a new perspective on mtDNA evolution: how the original mitochondrial genome was organized, what genes it contained, and in what ways it must have changed in different eukaryotic phyla.

Searching for substitution rate heterogeneity.

A permutation method for detecting regional substitution rate heterogeneity in DNA sequences is described. An estimated phylogeny for the sequences is required. The method is likelihood based and searches through all regions of a DNA sequence to find any regions that have an unusually large (or small) number of substitutions in comparison to an optimal substitution rate for the entire sequence. Likelihood ratio tests do not achieve proper statistical significance levels due to the large number of tests performed. Empirical corrections based on permutations are suggested and shown to give correct statistical accuracy. The ability of the method to detect rate heterogeneity is good even when the region is small and with less than twice the rate of evolution in the remainder of the sequence.

An ancestral mitochondrial DNA resembling a eubacterial genome in miniature.

Mitochondria, organelles specialized in energy conservation reactions in eukaryotic cells, have evolved from eubacteria-like endosymbionts whose closest known relatives are the rickettsial group of alpha-proteobacteria. Because characterized mitochondrial genomes vary markedly in structure, it has been impossible to infer from them the initial form of the proto-mitochondrial genome. This would require the identification of minimally derived mitochondrial DNAs that better reflect the ancestral state. Here we describe such a primitive mitochondrial genome, in the freshwater protozoon Reclinomonas americana. This protist displays ultrastructural characteristics that ally it with the retortamonads, a protozoan group that lacks mitochondria. R. americana mtDNA (69,034 base pairs) contains the largest collection of genes (97) so far identified in any mtDNA, including genes for 5S ribosomal RNA, the RNA component of RNase P, and at least 18 proteins not previously known to be encoded in mitochondria. Most surprising are four genes specifying a multisubunit, eubacterial-type RNA polymerase. Features of gene content together with eubacterial characteristics of genome organization and expression not found before in mitochondrial genomes indicate that R. americana mtDNA more closely resembles the ancestral proto-mitochondrial genome than any other mtDNA investigated to date.

Protein engineering reveals ancient adaptive replacements in isocitrate dehydrogenase.

Evolutionary analysis indicates that eubacterial NADP-dependent isocitrate dehydrogenases (EC 1.1.1.42) first evolved from an NAD-dependent precursor about 3.5 billion years ago. Selection in favor of utilizing NADP was probably a result of niche expansion during growth on acetate, where isocitrate dehydrogenase provides 90% of the NADPH necessary for biosynthesis. Amino acids responsible for differing coenzyme specificities were identified from x-ray crystallographic structures of Escherichia coli isocitrate dehydrogenase and the distantly related Thermus thermophilus NAD-dependent isopropylmalate dehydrogenase. Site-directed mutagenesis at sites lining the coenzyme binding pockets has been used to invert the coenzyme specificities of both enzymes. Reconstructed ancestral sequences indicate that these replacements are ancestral. Hence the adaptive history of molecular evolution is amenable to experimental investigation.

Substitution rate variation in closely related rodent species.

The existence of evolutionary rate variation has previously been demonstrated between different orders, different species and even between different regions of the same gene. To examine rate variation between closely related species of rodents we have sequenced the adenine phosphorybosyltransferase (APRT) gene from Mus spicilegus, Mus pahari, Mastomys hildebrandtii, Stochomys longicaudatus and Gerbillus campestris and compared these sequences with the previously published Mus musculus, Rattus norvegicus and Mesocricetus auratus APRT sequences. The alignment of these eight rodent APRT sequences reveals two large insertions within the introns: an insertion with sequence similar to a B1 repetitive element is found within Mastomys and an insertion with sequence similar to a B2 repetitive element is found within M. pahari. A phylogeny for the rodent APRTs agrees with the previously published rodent phylogeny based on other molecular and morphological data. The relative rate test which is often used to test for variation in rates of evolution in different lineages is shown here to be sensitive to the choice of outgroup and therefore should be used with great caution. This sensitivity is detectable only with closely related species and results from the prevalence of homoplastic substitutions. Rate variation is demonstrated within the APRT exons and introns and between the rodent species (with the most significant difference being a rate difference in M. spicilegus). In addition, some third codon positions are shown to be more prone to substitution than others. This clearly demonstrates that even between very closely related species there is ample evidence of major differences in rates of evolution among species, among regions of the gene and among different positions within the gene. We also demonstrate that standard methods of analysis might not detect this variation.

Evolutionary rate variation within Mus APRT.

Rodents are thought to have relatively high rates of evolution, twice as fast as the rates for mammals in other orders. However, the uniformly high rates of evolution inferred for the order Rodentia from Mus musculus and Rattus norvegicus are not consistently found for other rodent species. Using a maximum likelihood phylogenetic algorithm (DNAML), we show here that Mus spicilegus has a fivefold different rate of evolution in 1100 bp around the adenine phosphoribosyltransferase gene (APRT) since its divergence from a common ancestor with Mus musculus. A greater than threefold difference in rates is also found in a comparison of the number of evolutionary events directly detected from the APRT sequences of these two closely related Mus species. The evolutionary events can be directly detected, since M. spicilegus, M. musculus, and the four rodent outgroup species used to determine the ancestral sequence are so closely related. One of the major differences between M. spicilegus and M. musculus that might affect evolutionary rate is the degree of commensalism with man. The Mus species therefore provide a useful model for testing various hypotheses for the causes of rate variations between genes, and possibly, between lineages.

The origin of the eukaryotic cell.

Molecular sequence data are beginning to provide important insights into the evolutionary origin of eukaryotic cells. Global phylogenies of numerous protein sequences indicate that the eukaryotic cell nucleus is a chimera, which has received major contributions from both a Gram-negative eubacterium and an archaebacterium. Recent studies also indicate that the formation of the nuclear envelope and the endoplasmic reticulum was accompanied by duplication of genes for the molecular chaperone proteins (e.g. hsp70, hsp90), which facilitate protein transport across membranes. Based on these observations, it is suggested that the ancestral eukaryotic cell arose by a unique endosymbiotic event involving engulfment of an eocyte archaebacterium by a Gram-negative eubacterial host.

Protein-based phylogenies support a chimeric origin for the eukaryotic genome.

The phylogenetic position of the archaebacteria and the place of eukaryotes in the history of life remain a question of debate. Recent studies based on some protein-sequence data have obtained unusual phylogenies for these organisms. We therefore collected the protein sequences that were available with representatives from each of the major forms of life: the gram-negative bacteria, gram-positive bacteria, archaebacteria, and eukaryotes. Monophyletic, unrooted phylogenies based on these sequence data show that seven of 24 proteins yield a significant gram-positive-archaebacteria clade/gram-negative-eukaryotic clade. The phylogenies for these seven proteins cannot be explained by the traditional three-way split of the eukaryotes, archaebacteria, and eubacteria. Nine of the 24 proteins yield the traditional gram-positive-gram-negative clade/archaebacteria-eukaryotic clade. The remaining eight proteins give phylogenies that cannot be statistically distinguished. These results support the hypothesis of a chimeric origin for the eukaryotic cell nucleus formed from the fusion of an archaebacteria and a gram-negative bacteria.