Molecular phylogeny of eukaryotes

Occurrence, diel patterns, and the influence of melatonin on the photosynthetic performance of cultured Symbiodinium

Dinoflagellata is the earliest phylum in which true circadian regulation of melatonin rhythms has been convincingly demonstrated. Here, diel profiling of melatonin in a cultured member of this phylum belonging to the genus Symbiodinium indicated that melatonin levels oscillate with significant nocturnal peaks. However, unlike in other previously studied dinoflagellate species, the diel rhythmicity of melatonin in Symbiodinium did not persist under constant dark conditions. Thus, the oscillating pattern of melatonin in Symbiodinium is presumed not to be driven by endogenous circadian control of melatonin production, but rather by changes in the daily photocycle, most likely through a mechanism involving the enhanced photo-consumption of melatonin by free radicals. Although direct interactions of melatonin with detrimental radicals have been previously studied in several basal species, including dinoflagellates, none of these investigations addressed the effects that this molecule may have on photosynthesis, a major source of radical species in unicellular algae. In the present work, real-time monitoring of oxygen evolution in Symbiodinium cultures indicated a significant decrease in photosynthesis rates upon treatment with various doses of melatonin. Analyses of chlorophyll a fluorescence and xanthophyll cycle activity confirmed this effect and further revealed that this slowdown may occur through an enhanced engagement of photoprotective mechanisms in melatonin-treated cells. These findings are of great importance as they demonstrate that in certain photoautotroph species, the interactions of elevated melatonin levels with photosynthesis may extend beyond the general purpose of antioxidant protection.

Mathematical modeling and comparison of protein size distribution in different plant, animal, fungal and microbial species reveals a negative correlation between protein size and protein number, thus providing insight into the evolution of proteomes

BACKGROUND

The sizes of proteins are relevant to their biochemical structure and for their biological function. The statistical distribution of protein lengths across a diverse set of taxa can provide hints about the evolution of proteomes.

RESULTS

Using the full genomic sequences of over 1,302 prokaryotic and 140 eukaryotic species two datasets containing 1.2 and 6.1 million proteins were generated and analyzed statistically. The lengthwise distribution of proteins can be roughly described with a gamma type or log-normal model, depending on the species. However the shape parameter of the gamma model has not a fixed value of 2, as previously suggested, but varies between 1.5 and 3 in different species. A gamma model with unrestricted shape parameter described best the distributions in ~48% of the species, whereas the log-normal distribution described better the observed protein sizes in 42% of the species. The gamma restricted function and the sum of exponentials distribution had a better fitting in only ~5% of the species. Eukaryotic proteins have an average size of 472 aa, whereas bacterial (320 aa) and archaeal (283 aa) proteins are significantly smaller (33-40% on average). Average protein sizes in different phylogenetic groups were: Alveolata (628 aa), Amoebozoa (533 aa), Fornicata (543 aa), Placozoa (453 aa), Eumetazoa (486 aa), Fungi (487 aa), Stramenopila (486 aa), Viridiplantae (392 aa). Amino acid composition is biased according to protein size. Protein length correlated negatively with %C, %M, %K, %F, %R, %W, %Y and positively with %D, %E, %Q, %S and %T. Prokaryotic proteins had a different protein size bias for %E, %G, %K and %M as compared to eukaryotes.

CONCLUSIONS

Mathematical modeling of protein length empirical distributions can be used to asses the quality of small ORFs annotation in genomic releases (detection of too many false positive small ORFs). There is a negative correlation between average protein size and total number of proteins among eukaryotes but not in prokaryotes. The %GC content is positively correlated to total protein number and protein size in prokaryotes but not in eukaryotes. Small proteins have a different amino acid bias than larger proteins. Compared to prokaryotic species, the evolution of eukaryotic proteomes was characterized by increased protein number (massive gene duplication) and substantial changes of protein size (domain addition/subtraction).

Defining DNA-based operational taxonomic units for microbial-eukaryote ecology

DNA sequence information has increasingly been used in ecological research on microbial eukaryotes. Sequence-based approaches have included studies of the total diversity of selected ecosystems, studies of the autecology of ecologically relevant species, and identification and enumeration of species of interest for human health. It is still uncommon, however, to delineate protistan species based on their genetic signatures. The reluctance to assign species-level designations based on DNA sequences is in part a consequence of the limited amount of sequence information presently available for many free-living microbial eukaryotes and in part a consequence of the problematic nature of and debate surrounding the microbial species concept. Despite the difficulties inherent in assigning species names to DNA sequences, there is a growing need to attach meaning to the burgeoning amount of sequence information entering the literature, and there is a growing desire to apply this information in ecological studies. We describe a computer-based tool that assigns DNA sequences from environmental databases to operational taxonomic units at approximately species-level distinctions. This approach provides a practical method for ecological studies of microbial eukaryotes (primarily protists) by enabling semiautomated analysis of large numbers of samples spanning great taxonomic breadth. Derivation of the algorithm was based on an analysis of complete small-subunit (18S) rRNA gene sequences and partial gene sequences obtained from the GenBank database for morphologically described protistan species. The program was tested using environmental 18S rRNA data sets for two oceanic ecosystems. A total of 388 operational taxonomic units were observed for 2,207 sequences obtained from samples collected in the western North Atlantic and eastern North Pacific oceans.

Critical analysis of eukaryotic phylogeny: a case study based on the HSP70 family

Trichomonads, together with diplomonads and microsporidia, emerge at the base of the eukaryotic tree, on the basis of the small subunit rRNA phylogeny. However, phylogenies based on protein sequences such as tubulin are markedly different with these protists emerging much later. We have investigated 70 kDa heat-shock protein (HSP70), which could be a reliable phylogenetic marker. In eukaryotes, HSP70s are found in cytosol, endoplasmic reticulum, and organelles (mitochondria and chloroplasts). In Trichomonas vaginalis we identified nine different HSP70-encoding genes and sequenced three nearly complete cDNAs corresponding to cytosolic, endoplasmic reticulum, and mitochondrial-type HSP70. Phylogenies of eukaryotes were reconstructed using the classical methods while varying the number of species and characters considered. Almost all the undoubtedly monophyletic groups, defined by ultrastructural characters, were recovered. However, due to the long branch attraction phenomenon, the evolutionary rates were the main factor determining the position of species, even with the use of a close outgroup, which is an important advantage of HSP70 with respect to many other markers. Numerous variable sites are peculiar to Trichomonas and probably generated the artefactual placement of this species at the base of the eukaryotes or as the sister group of fast-evolving species. The inter-phyla relationships were not well supported and were sensitive to the reconstruction method, the number of species; and the quantity of information used. This lack of resolution could be explained by the very rapid diversification of eukaryotes, likely after the mitochondrial endosymbiosis.

Paramecium tetraurelia encodes unconventional actin containing short introns

The polymerase chain reaction was used to amplify and clone an actin gene fragment from Paramecium tetraurelia. This DNA fragment was 1,138 bp long, more than 96% of the actin coding sequence, and contained four in-frame UAA codons and two small introns located at positions unique in the actin intron catalogue. This is the first report for the phylum Ciliophora of an actin gene containing introns. The deduced amino acid sequence of this actin fragment shared 58-77% identity with other actins. When compared with rabbit alpha-muscle actin, similarities were observed mainly in subdomains 1 and 3, whereas subdomains 2 and 4 appeared to be more divergent.

Phylogenetic identification of the symbiotic hypermastigote Trichonympha agilis in the hindgut of the termite Reticulitermes speratus based on small-subunit rRNA sequence

The phylogeny of a symbiotic hypermastigote Trichonympha agilis (class Parabasalia; order Hypermastigida) in the hindgut of the lower termite Reticulitermes speratus was examined by a strategy that does not rely on cultivation. From mixed-population DNA obtained from the termite gut, small subunit (16S-like) ribosomal RNA sequences were directly amplified by the polymerase chain reaction method using primers specific for eukaryotes. Comparative sequence analysis of the clones revealed two kinds of sequences, one from the termite itself and the other from a symbiotic protist. A fluorescent-labeled oligonucleotide probe for the latter sequence was designed and used in whole-cell hybridization experiments to provide direct visual evidence that the sequence originated from a larger hypermastigote in the termite hindgut, Trichonympha agilis. According to the phylogenetic trees constructed, the hypermastigote represented one of the deepest branches of eukaryotes. The hypermastigote along with members of the order Trichomonadida formed a monophyletic lineage, indicating that this hypermastigote and trichomonads shared a recent common ancestry.

Phylogenetic position of Psalteriomonas lanterna deduced from the SSU rDNA sequence

The small subunit ribosomal DNA sequence (SSU rDNA) of the microaerophilic free-living amoeboflagellate Psalteriomonas lanterna has been sequenced and analyzed. The gene is 1,945 bp long and has a G + C content of 33.4%. Based upon ultrastructural studies, P. lanterna has been placed in the class Lyromonadea within the phylum Percolozoa Cavalier-Smith, 1991. However, based upon cytological characteristics, this microaerophilic free-living amoeboflagellate appears to be very primitive. It shares certain characteristics in common with some archezoans, i.e. it lacks mitochondria and dictyosomes but contains hydrogenosomes. Despite sharing these characteristics with the amitochondriate taxa, P. lanterna is not related to any of these taxa but instead to the Vahlkampfiidae. Therefore, we used primary sequence data and the secondary structure of the SSU rDNA gene to determine the placement P. lanterna in the phylogenetic tree. Our analyses showed that P. lanterna groups as a sister taxon to the Vahlkampfiidae but probably diverged from them quite early.

Organelle genomes: going, going, gone!

The organelles of eukaryotic cells—chloroplasts and mitochondria—first arose as engulfed symbionts with their own genomes. They subsequently lost most of their genes to the nucleus, retaining a few that could not be transferred. In his Perspective, Palmer discusses recent evidence that suggests that another organelle, the hydrogenosome, is a highly modified mitochondrion that has lost all of its genetic material.

Presence of a mitochondrial-type 70-kDa heat shock protein in Trichomonas vaginalis suggests a very early mitochondrial endosymbiosis in eukaryotes

Molecular phylogenetic analyses, based mainly on ribosomal RNA, show that three amitochondriate protist lineages, diplomonads, microsporidia, and trichomonads, emerge consistently at the base of the eukaryotic tree before groups having mitochondria. This suggests that these groups could have diverged before the mitochondrial endosymbiosis. Nevertheless, since all these organisms live in anaerobic environments, the absence of mitochondria might be due to secondary loss, as demonstrated for the later emerging eukaryote Entamoeba histolytica. We have now isolated from Trichomonas vaginalis a gene encoding a chaperone protein (HSP70) that in other lineages is addressed to the mitochondrial compartment. The phylogenetic reconstruction unambiguously located this HSP70 within a large set of mitochondrial sequences, itself a sister-group of alpha-purple bacteria. In addition, the T. vaginalis protein exhibits the GDAWV sequence signature, so far exclusively found in mitochondrial HSP70 and in proteobacterial dnaK. Thus mitochondrial endosymbiosis could have occurred earlier than previously assumed. The trichomonad double membrane-bounded organelles, the hydrogenosomes, could have evolved from mitochondria.

The cloning and sequencing of a ribosomal L18 protein from an evolutionary divergent eukaryote, Trypanosoma brucei

Eukaryotic ribosomal proteins are highly conserved across widely divergent species, suggesting that strong functional constraints prevent divergence of important amino acid motifs. Using this as a basis, an evolutionary approach could be used to identify putative functional motifs. We obtained the DNA sequence of the ribosomal protein L18 from the evolutionary divergent protozoan parasite, Trypanosoma brucei. Analysis of this sequence showed that it had 46% and 43% identity with the human and yeast sequences, respectively, and 30% of amino acid residues were identical across all the species analysed. Using these data, amino acids essential to the structure and function of ribosomal protein L18 can be inferred and could provide valuable information for molecular modelling and mutational studies.

Evidence for the Heterolobosea from phylogenetic analysis of genes encoding glyceraldehyde-3-phosphate dehydrogenase

The phylogenetic relationships between major slime mould groups and the identification of their unicellular relatives has been a subject of controversy for many years. Traditionally, it has been assumed that two slime mould groups, the acrasids and the dictyostelids were related by virtue of their cellular slime mould habit; a view still endorsed by at least one current classification scheme. However, a decade ago, on the basis of detailed ultrastructural resemblances it was proposed that acrasids of the family Acrasidae were not relatives of other slime moulds but instead related to a group of mostly free-living unicellular amoebae, the Schizopyrenida. The class Heterolobosea was created to contain these organisms and has since figured in many discussions of protist evolution. We sought to test the validity of Heterolobosea by characterizing homologs of the highly conserved glycolytic enzyme glyceraldehyde-3-phosphate dehydrogenase (GAPDH) from an acrasid, Acrasis rosea; a dictyostelid, Dictyostelium discoideum; and the schizopyrenid Naegleria andersoni. Phylogenetic analysis of these and other GAPDH sequences, using maximum parsimony, neighbour-joining distance and maximum likelihood methods strongly supports the Heterolobosea hypothesis and discredits the concept of a cellular slime mould grouping. Moreover, all of our analyses place Dictyostelium discoideum as a relatively recently originating lineage, most closely related to the Metazoa, similar to other recently published phylogenies of protein-coding genes. However, GAPDH phylogenies do not show robust branching orders for most of the relationships between major groups. We propose that several of the incongruencies observed between GAPDH and other molecular phylogenies are artifacts resulting from substitutional saturation of this enzyme.

Primary structure and phylogenetic relationships of glyceraldehyde-3-phosphate dehydrogenase genes of free-living and parasitic diplomonad flagellates

Complete nucleotide sequences have been established for two genes (gap1 and gap2) coding for glyceraldehyde-3-phosphate dehydrogenase (GAPDH, EC 1.2.1.12) homologs in the diplomonad Giardia lamblia. In addition, almost complete sequences of the GAPDH open reading frames were obtained from PCR products for two free-living diplomonad species, Trepomonas agillis and Hexamita inflata, and a parasite of Atlantic salmon, an as yet unnamed species with morphological affinities to Spironucleus. Giardia lamblia gap1 and the genes from the three other diplomonad species show high similarity to each other and to other glycolytic GAPDH genes. All amino-acyl residues known to be highly conserved in this enzyme are also conserved in these sequences. Giardia lamblia gap2 gene is more divergent and its putative translation reveals the presence of a cysteine and serine-rich insertion resembling a metal binding finger. This motif has not yet been noted in other GAPDH molecules. All sequences contain an S-loop signature with characteristics close to those of eukaryotes. In phylogenetic reconstructions based on the derived amino acid sequences with neighbor-joining, parsimony and maximum-likelihood methods the four typical GAPDH sequences of diplomonads cluster into a single clade. Within this clade, G. lambia gap1 shares a common ancestor with the rest of the genes. The latter are more closely related to each other, indicating an early separation of the lineage leading to the genus Giardia from the lineage encompassing the morphologically less differentiated genera, Trepomonas, Hexamita and that of the unnamed species. This result is discordant with the orthogonal evolution of diplomonads suggested on the basis of comparative morphology. In neighbor-joining reconstructions G. lamblia gap2 occupies a variable position, due to its great divergence. In parsimony and maximum likelihood analysis however, it shares a most recent common ancestor with the typical G. lamblia gap1 gene, suggesting that it diverged after the separation of the Giardia lineage. The position of the diplomonad clade in broader phylogenetic reconstructions is firmly within the typical cytosolic glycolytic representatives of GAPDH of eukaryotes.

Culture and phylogenetic characterization of Trichomitus trypanoides Duboscq & Grassè 1924, n. comb.: a trichomonad flagellate isolated from the hindgut of the termite Reticulitermes santonensis Feytaud

A trichomonad flagellate strain R1 was isolated from the hindgut contents of the termite Reticulitermes santonensis Feytaud. The flagellate was cultivated at 28 degrees C in anaerobic medium containing yeast extract, minerals and vitamins. The isolate fed on living bacteria. It showed the typical morphological and ultrastructural features of the trichomonads, closely resembling Trichomitus trypanoides. In order to determine its phylogenetic position the small subunit ribosomal DNA (SSU rDNA) of the flagellate was amplified in vitro using the polymerase chain reaction (PCR), cloned in a plasmid vector and sequenced. Comparison of the obtained sequence with so far available SSU rRNA/rDNA sequences showed strongest similarity (89%) to the sequence of Tritrichomonas foetus. The phylogenetic analysis with parsimony and distance matrix methods placed Trichomitus trypanoides strain R1 near by the root of the phylogenetically so far analyzed eukaryotic organisms. This confirms that termites harbour hindgut symbionts, which originate from very early evolved eukaryotes.

The evolutionary expansion of the trypanosomatid flagellates

The trypanosomatids combine a relatively uniform morphology with ability to parasitise a very diverse range of hosts including animals, plants and other protists. Along with their sister family, the biflagellate bodonids, they are set apart from other eukaryotes by distinctive organisational features, such as the kinetoplast-mitochondrion and RNA editing, isolation of glycolysis enzymes in the glycosome, use of the flagellar pocket for molecular traffic into and out of the cell, a unique method of generating cortical microtubules, and bizarre nuclear organisation. These features testify to the antiquity and isolation of the kinetoplast-bearing flagellates (Kinetoplastida). Molecular sequencing techniques (especially small subunit ribosomal RNA gene sequencing) are now radically reshaping previous ideas on the phylogeny of these organisms. The idea that the monogenetic (MG) trypanosomatids gave rise to the digenetic (DG) genera is losing ground to a view that, after the bodonids, the African trypanosomes (DG) represent the most ancient lineage, followed by Trypanosoma cruzi (DG), then Blastocrithidia (MG), Herpetomonas (MG) and Phytomonas (DG), with Leptomonas (MG), Crithidia (MG), Leishmania (DG) and Endotrypanum (DG) forming the crown of the evolutionary tree. Vast genetic distances (12% divergence) separate T. brucei and T. cruzi, while the Leishmania species are separated by very short distances (less than 1% divergence). These phylogenetic conclusions are supported by studies on RNA editing and on the nature of the parasite surface. The trypanosomatids seem to be able to adapt with ease their energy metabolism to the availability of substrates and oxygen, and this may give them the ability to institute new life cycles if host behaviour patterns allow. Sexual processes, though present in at least some trypanosomatids, may have played only a minor part in generating diversity during trypanosomatid evolution. On the other hand, the development of altruistic behaviour on the part of some life cycle stages may be a hitherto unconsidered way of maximising fitness in this group. It is concluded that, owing to organisational constraints, the trypanosomatids can undergo substantial molecular variation while registering very little in the way of morphological change.