Sequence Analysis of the Mitochondrial Genome of Sarcophyton glaucum: Conserved Gene Order Among Octocorals

The nucleotide sequence for an 11,715-bp segment of the mitochondrial genome of the octocoral Sarcophyton glaucum is presented, completing the analysis of the entire genome for this anthozoan member of the phylum Cnidaria. The genome contained the same 13 protein-coding and 2 ribosomal RNA genes as in other animals. However, it also included an unusual mismatch repair gene homologue reported previously and codes for only a single tRNA gene. Intermediate in length compared to two other cnidarians (17,443 and 18,911 bp), this organellar genome contained the smallest amount of noncoding DNA (428, compared to 1283 and 781 nt, respectively), making it the most compact one found for the phylum to date. The mitochondrial genes of S. glaucum exhibited an identical arrangement to that found in another octocoral, Renilla kolikeri, with five protein-coding genes in the same order as has been found in insect and vertebrate mitochondrial genomes. Although gene order appears to be highly conserved among octocorals, compared to the hexacoral, Metridium senile, few similarities were found. Like other metazoan mitochondrial genomes, the A + T composition was elevated and a general bias against codons ending in G or C was observed. However, an exception to this was the infrequent use of TGA compared to TGG to code for tryptophan. This divergent codon bias is unusual but appears to be a conserved feature among two rather distantly related anthozoans. Received: 27 January 1998 / Accepted: 25 May 1998     

Jam packed genomes--a preliminary,comparative analysis of nucleomorphs

There are two ways eukaryotic cells can permanently acquire chloroplasts. They can take up a cyanobacterium and turn it into a chloroplast or they can engulf an alga that already has a chloroplast. The second method is far more common and there are at least seven major groups of protists that have obtained their chloroplasts, this way. In most cases little remains of the engulfed alga apart from its chloroplast, but in two groups, the cryptomonads and chlorarachniophytes, a small remnant nucleus of the engulfed alga is still present. These tiny nuclei, called nucleomorphs, are the smallest and most compact eukaryotic genomes known and recently the nucleomorph of the cryptomonad alga Guillardia theta, was completely sequenced (551 kilobases). The nucleomorph of the chlorarachniophyte Bigellowiella natans (380 kilobases), is also being sequenced and is about half complete. We discuss some of the similarities and differences that are emerging between these two nucleomorph genomes. Both genomes contain just three chromosomes that encode mainly housekeeping genes and a few proteins for chloroplast functions. The bulk of nucleomorph gene coding capacity, therefore, appears to be devoted to self perpetuation and creating gene and protein expression machineries to make a small number of essential chloroplast proteins. We discuss reasons why both nucleomorphs are extraordinarily compact and why their gene sequences are evolving rapidly.     

Nuclear and Nucleomorph SSU rDNA Phylogeny in the Cryptophyta and the Evolution of Cryptophyte Diversity

The plastid-bearing members of the Cryptophyta contain two functional eukaryotic genomes of different phylogenetic origin, residing in the nucleus and in the nucleomorph, respectively. These widespread and diverse protists thus offer a unique opportunity to study the coevolution of two different eukaryotic genomes within one group of organisms. In this study, the SSU rRNA genes of both genomes were PCR-amplified with specific primers and phylogenetic analyses were performed on different data sets using different evolutionary models. The results show that the composition of the principal clades obtained from the phylogenetic analyses of both genes was largely congruent, but striking differences in evolutionary rates were observed. These affected the topologies of the nuclear and nucleomorph phylogenies differently, resulting in long-branch attraction artifacts when simple evolutionary models were applied. Deletion of long-branch taxa stabilized the internal branching order in both phylogenies and resulted in a completely resolved topology in the nucleomorph phylogeny. A comparison of the tree topologies derived from SSU rDNA sequences with characters previously used in cryptophyte systematics revealed that the biliprotein type was congruent, but the type of inner periplast component incongruent, with the molecular trees. The latter is indicative of a hidden cellular dimorphism (cells with two periplast types present in a single clonal strain) of presumably widespread occurrence throughout cryptophyte diversity, which, in consequence, has far-reaching implications for cryptophyte systematics as it is practiced today.     

Identification of Protein-Tyrosine Phosphatases in Archaea

Protein-tyrosine dephosphorylation is a major mechanism in cellular regulation. A large number of protein-tyrosine phosphatases is known from Eukarya, and more recently bacterial homologues have also been identified. By employing conserved sequence patterns from both eukaryotic and bacterial protein-tyrosine phosphatases, we have identified three homologous sequences in two of the four complete archaeal genomes. Two hypothetical open reading frames in the genome of Methanococcus jannaschii (MJ0215 and MJECL20) and one in the genome of Pyrococcus horikoshii (PH1732) clearly bear all the conserved residues of this family. No homologues were found in the genomes of Archaeoglobus fulgidus and Methanobacterium thermoautotrophicum. This is the first report of protein-tyrosine phosphatase sequences in Archaea. Received: 29 April 1998 / Accepted: 27 November 1998     

Nucleomorph genomes: much ado about practically nothing

The DNA sequence of one of the smallest eukaryotic genomes has recently been finished - that of the reduced nucleus, or nucleomorph, of an algal endosymbiont that resides within a cryptomonad host cell. Its sequence promises insights into chloroplast acquisition, the constraints on genome size and the basic workings of eukaryotic cells.     

The Structure and Gene Repertoire of an Ancient Red Algal Plastid Genome

Photosynthetic eukaryotes can, according to features of their chloroplasts, be divided into two major groups: the red and the green lineage of plastid evolution. To extend the knowledge about the evolution of the red lineage we have sequenced and analyzed the chloroplast genome (cp-genome) of Cyanidium caldarium RK1, a unicellular red alga (AF022186). The analysis revealed that this genome shows several unusual structural features, such as a hypothetical hairpin structure in a gene-free region and absence of large repeat units. We provide evidence that this structural organization of the cp-genome of C. caldarium may be that of the most ancient cp-genome so far described. We also compared the cp-genome of C. caldarium to the other known cp-genomes of the red lineage. The cp-genome of C. caldarium cannot be readily aligned with that of Porphyra purpurea, a multicellular red alga, or Guillardia theta due to a displacement of a region of the cp-genome. The phylogenetic tree reveals that the secondary endosymbiosis, through which G. theta evolved, took place after the separation of the ancestors of C. caldarium and P. purpurea. We found several genes unique to the cp-genome of C. caldarium. Five of them seem to be involved in the building of bacterial cell envelopes and may be responsible for the thermotolerance of the chloroplast of this alga. Two additional genes may play a role in stabilizing the photosynthetic machinery against salt stress and detoxification of the chloroplast. Thus, these genes may be unique to the cp-genome of C. caldarium and may be required for the endurance of the extreme living conditions of this alga. Received: 3 June 2000 / Accepted: 18 July 2000     

The chlorarachniophyte: a cell with two different nuclei and two different telomeres

Chlorarachniophyte algae contain a complex chloroplast derived from the endosymbiosis of a eukaryotic alga. The reduced nucleus of the endosymbiont, the nucleomorph, is located between the inner and outer pair of membranes surrounding the chloroplast. The nucleomorph of chlorarachniophytes has previously been demonstrated to contain at least three small linear chromosomes. Here we describe cloning the end of the smallest nucleomorph chromosome which is shown to carry a telomere consisting of a tandemly repeated 7 bp sequence, TCTAGGG. Using the telomere repeat as a probe, we show that nucleomorph telomeres display typical hetero-disperse size distribution. The nucleomorph is shown to contain only three chromosomes with a haploid genome size of just 380kb. All six nucleomorph chromosome termini are identical with an rRNA cistron closely linked to the telomere. The nucleomorph chromosomes thus have relatively large inverted repeats at their ends. Chromosomes from the host nucleus are shown to have a different telomere repeat motif to that of the nucleomorph chromosomes.     

Phylogeny of Ultra-Rapidly Evolving Dinoflagellate Chloroplast Genes: A Possible Common Origin for Sporozoan and Dinoflagellate Plastids

Complete chloroplast 23S rRNA and psbA genes from five peridinin-containing dinoflagellates (Heterocapsa pygmaea, Heterocapsa niei, Heterocapsa rotun-data, Amphidinium carterae, and Protoceratium reticulatum) were amplified by PCR and sequenced; partial sequences were obtained from Thoracosphaera heimii and Scrippsiella trochoidea. Comparison with chloroplast 23S rRNA and psbA genes of other organisms shows that dinoflagellate chloroplast genes are the most divergent and rapidly evolving of all. Quartet puzzling, maximum likelihood, maximum parsimony, neighbor joining, and LogDet trees were constructed. Intersite rate variation and invariant sites were allowed for with quartet puzzling and neighbor joining. All psbA and 23S rRNA trees showed peridinin-containing dinoflagellate chloroplasts as monophyletic. In psbA trees they are related to those of chromists and red algae. In 23S rRNA trees, dinoflagellates are always the sisters of Sporozoa (apicomplexans); maximum likelihood analysis of Heterocapsa triquetra 16S rRNA also groups the dinoflagellate and sporozoan sequences, but the other methods were inconsistent. Thus, dinoflagellate chloroplasts may actually be related to sporozoan plastids, but the possibility of reproducible long-branch artifacts cannot be strongly ruled out. The results for all three genes fit the idea that dinoflagellate chloroplasts originated from red algae by a secondary endosymbiosis, possibly the same one as for chromists and Sporozoa. The marked disagreement between 16S rRNA trees using different phylogenetic algorithms indicates that this is a rather poor molecule for elucidating overall chloroplast phylogeny. We discuss possible reasons why both plastid and mitochondrial genomes of alveolates (Dinozoa, Sporozoa and Ciliophora) have ultra-rapid substitution rates and a proneness to unique genomic rearrangements. Received: 27 December 1999 / Accepted: 24 March 2000     

Primary Structure and Phylogenetic Relationships of a Malate Dehydrogenase Gene from Giardia lamblia

The lactate and malate dehydrogenases comprise a complex protein superfamily with multiple enzyme homologues found in eubacteria, archaebacteria, and eukaryotes. In this study we describe the sequence and phylogenetic relationships of a malate dehydrogenase (MDH) gene from the amitochondriate diplomonad protist, Giardia lamblia. Parsimony, distance, and maximum-likelihood analyses of the MDH protein family solidly position G. lamblia MDH within a eukaryote cytosolic MDH clade, to the exclusion of chloroplast, mitochondrial, and peroxisomal homologues. Furthermore, G. lamblia MDH is specifically related to a homologue from Trichomonas vaginalis. This MDH topology, together with published phylogenetic analyses of β-tubulin, chaperonin 60, valyl-tRNA synthetase, and EF-1α, suggests a sister-group relationship between diplomonads and parabasalids. Since these amitochondriate lineages contain genes encoding proteins which are characteristic of mitochondria and α-proteobacteria, their shared ancestry suggests that mitochondrial properties were lost in the common ancestor of both groups. Received: 14 September 1998 / Accepted: 29 December 1998     

The Dentin Matrix Protein 1 Gene of Prototherian and Metatherian Mammals

Mineralization of tooth dentin (the deposition of hydroxyapatite crystals in and around collagen type I fibers of the extracellular matrix) requires the involvement of several genes, among them the gene coding for the dentin matrix protein 1, DMP1. We determined the exon–intron organization of the cattle DMP1 gene and used this information to amplify by the polymerase chain reaction homologous gene fragments from the genomic DNA of two species of metatherian (marsupial) mammals and one prototherian (monotreme) species. The translated proto- and metatherian protein sequences are highly divergent from the eutherian sequences but retain the general characteristics of the DMP1 (high acidity, serine-richness, multiple glycosylation sites, and the presence of the RGD cell attachment tripeptide). They therefore appear to be functional even though, evolutionarily, teeth are in a regression phase in prototherians. It is possible, therefore, that DMP1 is also involved in other functions besides dentinogenesis. The DMP1 gene appears to evolve rapidly and apparently tolerates non-frame-shifting insertions/deletions throughout the coding sequence. Received: 22 February 1998 / Accepted: 11 July 1998     

Using Homolog Groups to Create a Whole-Genomic Tree of Free-Living Organisms: An Update

Genomic trees have been constructed based on the presence and absence of families of protein-encoding genes observed in 27 complete genomes, including genomes of 15 free-living organisms. This method does not rely on the identification of suspected orthologs in each genome, nor the specific alignment used to compare gene sequences because the protein-encoding gene families are formed by grouping any protein with a pairwise similarity score greater than a preset value. Because of this all inclusive grouping, this method is resilient to some effects of lateral gene transfer because transfers of genes are masked when the recipient genome already has a homolog (not necessarily an ortholog) of the incoming gene. Of 71 genes suspected to have been laterally transferred to the genome of Aeropyrum pernix, only approximately 7 to 15 represent genes where a lateral gene transfer appears to have generated homoplasy in our character dataset. The genomic tree of the 15 free-living taxa includes six different bacterial orders, six different archaeal orders, and two different eukaryotic kingdoms. The results are remarkably similar to results obtained by analysis of rRNA. Inclusion of the other 12 genomes resulted in a tree only broadly similar to that suggested by rRNA with at least some of the differences due to artifacts caused by the small genome size of many of these species. Very small genomes, such as those of the two Mycoplasma genomes included, fall to the base of the Bacterial domain, a result expected due to the substantial gene loss inherent to these lineages. Finally, artificial ``partial genomes' were generated by randomly selecting ORFs from the complete genomes in order to test our ability to recover the tree generated by the whole genome sequences when only partial data are available. The results indicated that partial genomic data, when sampled randomly, could robustly recover the tree generated by the whole genome sequences. Received: 30 May 2001 / Accepted: 10 October 2001     

Phylogeny and nucleomorph karyotype diversity of chlorarachniophyte algae

Chlorarachniophytes are flagellated and/or reticulopod-forming marine algae with chlorophyll a- and b-containing plastids of secondary endosymbiotic origin. They are one of only two algal groups known to possess a "nucleomorph" (i.e. the remnant nucleus of the eukaryotic endosymbiont that donated the plastid). Apart from the recently sequenced nucleomorph genome of Bigelowiella natans, little is known about the size, structure, and composition of chlorarachniophyte nucleomorph genomes. Toward the goal of better understanding nucleomorph genome diversity, as well as establishing a phylogenetic framework with which to interpret variation in chlorarachniophyte morphology, ultrastructure, and life cycle, we are studying a wide range of chlorarachniophyte strains from public culture collections and natural habitats. We have obtained 22 new chlorarachniophyte nuclear and nucleomorph 18S rRNA gene (18S rDNA) sequences and nucleomorph genome size estimates for 14 different strains. Consistent with previous studies, all of the chlorarachniophytes examined appear to possess three nucleomorph chromosomes. However, our results suggest considerable variation in nucleomorph genome size and structure, with individual chromosome sizes ranging from approximately 90 to approximately 210 kbp, and total genome sizes between approximately 330 kbp in Lotharella amoebiformis and approximately 610 kbp in unidentified chlorarachniophyte strain CCMP622. The significance of these phylogenetic and nucleomorph karyotype data is discussed.     

Selection on the codon bias of chloroplast and cyanelle genes in different plant and algal lineages

In the plant chloroplast genome the codon usage of the highly expressed psbA gene is unique and is adapted to the tRNA population, probably due to selection for translation efficiency. In this study the role of selection on codon usage in each of the fully sequenced chloroplast genomes, in addition to Chlamydomonas reinhardtii, is investigated by measuring adaptation to this pattern of codon usage. A method is developed which tests selection on each gene individually by constructing sequences with the same amino acid composition as the gene and randomly assigning codons based on the nucleotide composition of noncoding regions of that genome. The codon bias of the actual gene is then compared to a distribution of random sequences. The data indicate that within the algae selection is strong in Cyanophora paradoxa, affecting a majority of genes, of intermediate intensity in Odontella sinensis, and weaker in Porphyra purpurea and Euglena gracilis. In the plants, selection is found to be quite weak in Pinus thunbergii and the angiosperms but there is evidence that an intermediate level of selection exists in the liverwort Marchantia polymorpha. The role of selection is then further investigated in two comparative studies. It is shown that average relative codon bias is correlated with expression level and that, despite saturation levels of substitution, there is a strong correlation among the algae genomes in the degree of codon bias of homologous genes. All of these data indicate that selection for translation efficiency plays a significant role in determining the codon bias of chloroplast genes but that it acts with different intensities in different lineages. In general it is stronger in the algae than the higher plants, but within the algae Euglena is found to have several unusual features which are noted. The factors that might be responsible for this variation in intensity among the various genomes are discussed. Received: 6 June 1997 / Accepted: 24 July 1997     

A Translocated Mitochondrial Cytochrome b Pseudogene in Voles (Rodentia: Microtus)

A full-length cytochrome b pseudogene was found in rodents; it has apparently been translocated from a mitochondrion to the nuclear genome in the subfamily Arvicolinae. The pseudogene (ψcytb) differed from its mitochondrial counterpart at 201 of 1143 sites (17.6%) and by four indels. Cumulative evidence suggests that the pseudogene has been translocated to the nucleus. Phylogenetic reconstruction indicates that the pseudogene arose before the diversification of M. arvalis/M. rossiaemeridionalis from M. oeconomus, but after the divergence of the peromyscine/sigmodontine/arvicoline clades some ∼10 MYA. Published rates of divergence between mitochondrial genes and their nuclear pseudogenes suggest that the translocation of this mitochondrial gene to the nuclear genome occurred some 6 MYA, in agreement with the phylogenetic evidence. Received: 16 January 1998 / Accepted: 18 July 1998     

Gene Order Breakpoint Evidence in Animal Mitochondrial Phylogeny

Multiple genome rearrangement methodology facilitates the inference of animal phylogeny from gene orders on the mitochondrial genome. The breakpoint distance is preferable to other, highly correlated but computationally more difficult, genomic distances when applied to these data. A number of theories of metazoan evolution are compared to phylogenies reconstructed by ancestral genome optimization, using a minimal total breakpoints criterion. The notion of unambiguously reconstructed segments is introduced as a way of extracting the invariant aspects of multiple solutions for a given ancestral genome; this enables a detailed reconstruction of the evolution of non-tRNA mitochondrial gene order. Received: 15 July 1998 / Accepted: 5 March 1999     

Heat Shock Protein 70 Family: Multiple Sequence Comparisons, Function, and Evolution

The heat shock protein 70 kDa sequences (HSP70) are of great importance as molecular chaperones in protein folding and transport. They are abundant under conditions of cellular stress. They are highly conserved in all domains of life: Archaea, eubacteria, eukaryotes, and organelles (mitochondria, chloroplasts). A multiple alignment of a large collection of these sequences was obtained employing our symmetric-iterative ITERALIGN program (Brocchieri and Karlin 1998). Assessments of conservation are interpreted in evolutionary terms and with respect to functional implications. Many archaeal sequences (methanogens and halophiles) tend to align best with the Gram-positive sequences. These two groups also miss a signature segment [about 25 amino acids (aa) long] present in all other HSP70 species (Gupta and Golding 1993). We observed a second signature sequence of about 4 aa absent from all eukaryotic homologues, significantly aligned in all prokaryotic sequences. Consensus sequences were developed for eight groups [Archaea, Gram-positive, proteobacterial Gram-negative, singular bacteria, mitochondria, plastids, eukaryotic endoplasmic reticulum (ER) isoforms, eukaryotic cytoplasmic isoforms]. All group consensus comparisons tend to summarize better the alignments than do the individual sequence comparisons. The global individual consensus ``matches' 87% with the consensus of consensuses sequence. A functional analysis of the global consensus identifies a (new) highly significant mixed charge cluster proximal to the carboxyl terminus of the sequence highlighting the hypercharge run EEDKKRRER (one-letter aa code used). The individual Archaea and Gram-positive sequences contain a corresponding significant mixed charge cluster in the location of the charge cluster of the consensus sequence. In contrast, the four Gram-negative proteobacterial sequences of the alignment do not have a charge cluster (even at the 5% significance level). All eukaryotic HSP70 sequences have the analogous charge cluster. Strikingly, several of the eukaryotic isoforms show multiple mixed charged clusters. These clusters were interpreted with supporting data related to HSP70 activity in facilitating chaperone, transport, and secretion function. We observed that the consensus contains only a single tryptophan residue and a single conserved cysteine. This is interpreted with respect to the target rule for disaggregating misfolded proteins. The mitochondrial HSP70 connections to bacterial HSP70 are analyzed, suggesting a polyphyletic split of Trypanosoma and Leishmania protist mitochondrial (Mt) homologues separated from Mt-animal/fungal/plant homologues. Moreover, the HSP70 sequences from the amitochondrial Entamoeba histolytica and Trichomonas vaginalis species were analyzed. The E. histolytica HSP70 is most similar to the higher eukaryotic cytoplasmic sequences, with significantly weaker alignments to ER sequences and much diminished matching to all eubacterial, mitochondrial, and chloroplast sequences. This appears to be at variance with the hypothesis that E. histolytica rather recently lost its mitochondrial organelle. T. vaginalis contains two HSP70 sequences, one Mt-like and the second similar to eukaryotic cytoplasmic sequences suggesting two diverse origins. Received: 29 January 1998 / Accepted: 14 May 1998