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

The microsporidian polar tube: origin, structure, composition, function, and application

Microsporidia are a class of obligate intracellular parasitic unicellular eukaryotes that infect a variety of hosts, even including humans. Although different species of microsporidia differ in host range and specificity, they all share a similar infection organelle, the polar tube, which is also defined as the polar filament in mature spores. In response to the appropriate environmental stimulation, the spore germinates with the polar filament everted, forming a hollow polar tube, and then the infectious cargo is transported into host cells via the polar tube. Hence, the polar tube plays a key role in microsporidian infection. Here, we review the origin, structure, composition, function, and application of the microsporidian polar tube, focusing on the origin of the polar filament, the structural differences between the polar filament and polar tube, and the characteristics of polar tube proteins. Comparing the three-dimensional structure of PTP6 homologous proteins provides new insight for the screening of additional novel polar tube proteins with low sequence similarity in microsporidia. In addition, the interaction of the polar tube with the spore wall and the host are summarized to better understand the infection mechanism of microsporidia. Due to the specificity of polar tube proteins, they are also used as the target in the diagnosis and prevention of microsporidiosis. With the present findings, we propose a future study on the polar tube of microsporidia.

Lower statistical support with larger datasets: insights from the Ochrophyta radiation

It is commonly assumed that increasing the number of characters has the potential to resolve evolutionary radiations. Here, we studied photosynthetic stramenopiles (Ochrophyta) using alignments of heterogeneous origin mitochondrion, plastid, and nucleus. Surprisingly while statistical support for the relationships between the six major Ochrophyta lineages increases when comparing the mitochondrion (6,762 sites) and plastid (21,692 sites) trees, it decreases in the nuclear (209,105 sites) tree. Statistical support is not simply related to the data set size but also to the quantity of phylogenetic signal available at each position and our ability to extract it. Here, we show that this ability for current phylogenetic methods is limited, because conflicting results were obtained when varying taxon sampling. Even though the use of a better fitting model improved signal extraction and reduced the observed conflicts, the plastid data set provided higher statistical support for the ochrophyte radiation than the larger nucleus data set. We propose that the higher support observed in the plastid tree is due to an acceleration of the evolutionary rate in one short deep internal branch, implying that more phylogenetic signal per position is available to resolve the Ochrophyta radiation in the plastid than in the nuclear data set. Our work therefore suggests that, in order to resolve radiations, beyond the obvious use of data sets with more positions, we need to continue developing models of sequence evolution that better extract the phylogenetic signal and design methods to search for genes/characters that contain more signal specifically for short internal branches.

Life Cycle of Cryptococcus neoformans

Cryptococcus neoformans is a ubiquitous environmental fungus and an opportunistic pathogen that causes fatal cryptococcal meningitis. Advances in genomics, genetics, and cellular and molecular biology of C. neoformans have dramatically improved our understanding of this important pathogen, rendering it a model organism to study eukaryotic biology and microbial pathogenesis. In light of recent progress, we describe in this review the life cycle of C. neoformans with a special emphasis on the regulation of the yeast-to-hypha transition and different modes of sexual reproduction, in addition to the impacts of the life cycle on cryptococcal populations and pathogenesis.

Analysis of the genomic basis of functional diversity in dinoflagellates using a transcriptome-based sequence similarity network

Dinoflagellates are one of the most abundant and functionally diverse groups of eukaryotes. Despite an overall scarcity of genomic information for dinoflagellates, constantly emerging high-throughput sequencing resources can be used to characterize and compare these organisms. We assembled de novo and processed 46 dinoflagellate transcriptomes and used a sequence similarity network (SSN) to compare the underlying genomic basis of functional features within the group. This approach constitutes the most comprehensive picture to date of the genomic potential of dinoflagellates. A core-predicted proteome composed of 252 connected components (CCs) of putative conserved protein domains (pCDs) was identified. Of these, 206 were novel and 16 lacked any functional annotation in public databases. Integration of functional information in our network analyses allowed investigation of pCDs specifically associated with functional traits. With respect to toxicity, sequences homologous to those of proteins found in species with toxicity potential (e.g., sxtA4 and sxtG) were not specific to known toxin-producing species. Although not fully specific to symbiosis, the most represented functions associated with proteins involved in the symbiotic trait were related to membrane processes and ion transport. Overall, our SSN approach led to identification of 45,207 and 90,794 specific and constitutive pCDs of, respectively, the toxic and symbiotic species represented in our analyses. Of these, 56% and 57%, respectively (i.e., 25,393 and 52,193 pCDs), completely lacked annotation in public databases. This stresses the extent of our lack of knowledge, while emphasizing the potential of SSNs to identify candidate pCDs for further functional genomic characterization.

Microsporidia: Obligate Intracellular Pathogens Within the Fungal Kingdom

Microsporidia are obligate intracellular pathogens related to Fungi. These organisms have a unique invasion organelle, the polar tube, which upon appropriate environmental stimulation rapidly discharges out of the spore, pierces a host cell's membrane, and serves as a conduit for sporoplasm passage into the host cell. Phylogenetic analysis suggests that microsporidia are related to the Fungi, being either a basal branch or sister group. Despite the description of microsporidia over 150 years ago, we still lack an understanding of the mechanism of invasion, including the role of various polar tube proteins, spore wall proteins, and host cell proteins in the formation and function of the invasion synapse. Recent advances in ultrastructural techniques are helping to better define the formation and functioning of the invasion synapse. Over the past 2 decades, proteomic approaches have helped define polar tube proteins and spore wall proteins as well as the importance of posttranslational modifications such as glycosylation in the functioning of these proteins, but the absence of genetic techniques for the manipulation of microsporidia has hampered research on the function of these various proteins. The study of the mechanism of invasion should provide fundamental insights into the biology of these ubiquitous intracellular pathogens that can be integrated into studies aimed at treating or controlling microsporidiosis.

The complex evolutionary dynamics of Hsp70s: a genomic and functional perspective

Hsp70 molecular chaperones are ubiquitous. By preventing aggregation, promoting folding, and regulating degradation, Hsp70s are major factors in the ability of cells to maintain proteostasis. Despite a wealth of functional information, little is understood about the evolutionary dynamics of Hsp70s. We undertook an analysis of Hsp70s in the fungal clade Ascomycota. Using the well-characterized 14 Hsp70s of Saccharomyces cerevisiae, we identified 491 orthologs from 53 genomes. Saccharomyces cerevisiae Hsp70s fall into seven subfamilies: four canonical-type Hsp70 chaperones (SSA, SSB, KAR, and SSC) and three atypical Hsp70s (SSE, SSZ, and LHS) that play regulatory roles, modulating the activity of canonical Hsp70 partners. Each of the 53 surveyed genomes harbored at least one member of each subfamily, and thus establishing these seven Hsp70s as units of function and evolution. Genomes of some species contained only one member of each subfamily that is only seven Hsp70s. Overall, members of each subfamily formed a monophyletic group, suggesting that each diversified from their corresponding ancestral gene present in the common ancestor of all surveyed species. However, the pattern of evolution varied across subfamilies. At one extreme, members of the SSB subfamily evolved under concerted evolution. At the other extreme, SSA and SSC subfamilies exhibited a high degree of copy number dynamics, consistent with a birth–death mode of evolution. KAR, SSE, SSZ, and LHS subfamilies evolved in a simple divergent mode with little copy number dynamics. Together, our data revealed that the evolutionary history of this highly conserved and ubiquitous protein family was surprising complex and dynamic.

Mutant HSP70 reverses autoimmune depigmentation in vitiligo

Vitiligo is an autoimmune disease characterized by destruction of melanocytes, leaving 0.5% of the population with progressive depigmentation. Current treatments offer limited efficacy. We report that modified inducible heat shock protein 70 (HSP70i) prevents T cell-mediated depigmentation. HSP70i is the molecular link between stress and the resultant immune response. We previously showed that HSP70i induces an inflammatory dendritic cell (DC) phenotype and is necessary for depigmentation in vitiligo mouse models. Here, we observed a similar DC inflammatory phenotype in vitiligo patients. In a mouse model of depigmentation, DNA vaccination with a melanocyte antigen and the carboxyl terminus of HSP70i was sufficient to drive autoimmunity. Mutational analysis of the HSP70i substrate-binding domain established the peptide QPGVLIQVYEG as invaluable for DC activation, and mutant HSP70i could not induce depigmentation. Moreover, mutant HSP70iQ435A bound human DCs and reduced their activation, as well as induced a shift from inflammatory to tolerogenic DCs in mice. HSP70iQ435A-encoding DNA applied months before spontaneous depigmentation prevented vitiligo in mice expressing a transgenic, melanocyte-reactive T cell receptor. Furthermore, use of HSP70iQ435A therapeutically in a different, rapidly depigmenting model after loss of differentiated melanocytes resulted in 76% recovery of pigmentation. Treatment also prevented relevant T cells from populating mouse skin. In addition, ex vivo treatment of human skin averted the disease-related shift from quiescent to effector T cell phenotype. Thus, HSP70iQ435A DNA delivery may offer potent treatment opportunities for vitiligo.

Diamidines versus Monoamidines as Anti-Pneumocystis Agents: An in Vivo Study

Some compounds articulated around a piperazine or an ethylenediamine linker have been evaluated in vitro to determine their activity in the presence of a 3T6 fibroblast cell line and an axenic culture of Pneumocystis carinii, respectively. The most efficient antifungal derivatives, namely N,N′-bis(benzamidine-4-yl)ethane-1,2-diamine (compound 6, a diamidine) and N-(benzamidine-4-yl)-N′-phenylethane-1,2-diamine (compound 7, a monoamidine), exhibited no cytotoxicity and were evaluated in vivo in a rat model. Only the diamidine 6 emerged as a promising hit for further studies.

The ATP-binding cassette proteins of the deep-branching protozoan parasite Trichomonas vaginalis

The ATP binding cassette (ABC) proteins are a family of membrane transporters and regulatory proteins responsible for diverse and critical cellular process in all organisms. To date, there has been no attempt to investigate this class of proteins in the infectious parasite Trichomonas vaginalis. We have utilized a combination of bioinformatics, gene sequence analysis, gene expression and confocal microscopy to investigate the ABC proteins of T. vaginalis. We demonstrate that, uniquely among eukaryotes, T. vaginalis possesses no intact full-length ABC transporters and has undergone a dramatic expansion of some ABC protein sub-families. Furthermore, we provide preliminary evidence that T. vaginalis is able to read through in-frame stop codons to express ABC transporter components from gene pairs in a head-to-tail orientation. Finally, with confocal microscopy we demonstrate the expression and endoplasmic reticulum localization of a number of T. vaginalis ABC transporters.

The parasite Trichomonas vaginalis infects in excess of 100 million people per year, and is a contributory factor to enhanced transmission rates of HIV, the causative virus in AIDS. As such, T. vaginalis infection is an important public health concern. Understanding the biology of the organism is important to determine aspects of the response to drug treatment, host:parasite interactions and so on. We have investigated an important family of proteins – the ATP binding cassette transporters – which are present in the membranes of all cells, and which contribute to a diverse spectrum of important cellular processes. The ABC transporters of T. vaginalis were identified by analysis of primary amino acid sequence data, and examined by subsequent protein and gene expression studies. Our most important conclusion is that – uniquely amongst eukaryotes - T. vaginalis has no ABC transporters capable of acting as monomers. In other words, its ABC transporters must all act by forming functional complexes with other ABC proteins. This has implications for our understanding not just of the parasite's biology, but also its evolution. In summary our analysis opens up the path for future research of individual members of the ABC protein family in T. vaginalis.

Pneumocystis carinii causes a distinctive interstitial pneumonia in immunocompetent laboratory rats that had been attributed to "rat respiratory virus"

A prevalent and distinctive infectious interstitial pneumonia (IIP) of immunocompetent laboratory rats was suspected to be caused by a putative virus, termed rat respiratory virus, but this was never substantiated. To study this disease, 2 isolators were independently populated with rats from colonies with endemic disease, which was perpetuated by the regular addition of naive rats. After Pneumocystis was demonstrated by histopathology and polymerase chain reaction (PCR) in the lungs of rats from both isolators and an earlier bedding transmission study, the relationship between Pneumocystis and IIP was explored further by analyzing specimens from 3 contact transmission experiments, diagnostic submissions, and barrier room breeding colonies, including 1 with and 49 without IIP. Quantitative (q) PCR and immunofluorescence assay only detected Pneumocystis infection and serum antibodies in rats from experiments or colonies in which IIP was diagnosed by histopathology. In immunocompetent hosts, the Pneumocystis concentration in lungs corresponded to the severity and prevalence of IIP; seroconversion occurred when IIP developed and was followed by the concurrent clearance of Pneumocystis from lungs and resolution of disease. Experimentally infected immunodeficient RNU rats, by contrast, did not seroconvert to Pneumocystis or recover from infection. qPCR found Pneumocystis at significantly higher concentrations and much more often in lungs than in bronchial and nasal washes and failed to detect Pneumocystis in oral swabs. The sequences of a mitochondrial ribosomal large-subunit gene region for Pneumocystis from 11 distinct IIP sources were all identical to that of P. carinii. These data provide substantial evidence that P. carinii causes IIP in immunocompetent rats.

Characterization of novel Leishmania infantum recombinant proteins encoded by genes from five families with distinct capacities for serodiagnosis of canine and human visceral leishmaniasis

To expand the available panel of recombinant proteins that can be useful for identifying Leishmania-infected dogs and for diagnosing human visceral leishmaniasis (VL), we selected recombinant antigens from L. infantum, cDNA, and genomic libraries by using pools of serum samples from infected dogs and humans. The selected DNA fragments encoded homologs of a cytoplasmic heat-shock protein 70, a kinesin, a polyubiquitin, and two novel hypothetical proteins. Histidine-tagged recombinant proteins were produced after subcloning these DNA fragments and evaluated by using an enzyme-linked immunosorbent assays with panels of canine and human serum samples. The enzyme-linked immunosorbent assays with different recombinant proteins had different sensitivities (67.4-93.0% and 36.4-97.2%) and specificities (76.1-100% and 90.4-97.3%) when tested with serum samples from Leishmania-infected dogs and human patients with VL. Overall, no single recombinant antigen was sufficient to serodiagnosis all canine or human VL cases.

A tandem duplication of manganese superoxide dismutase in Nosema bombycis and its evolutionary origins

Microsporidia are a group of obligate intracellular eukaryotic parasites with small genomes. They infect animals from a wide variety of phyla, including humans. Two manganese superoxide dismutase (MnSOD) genes, designated NbMnSOD1 and NbMnSOD2, were found to be organized in a tandem array within the Nosema bombycis genome. The genes, both 678 bp in length, were found to be more similar to each other than they are to homologous genes of other Microsporidia, suggesting that the tandem duplication occurred subsequent to the development of this lineage. Reverse transcript PCR shows that mRNA for both genes is present in the spores. Analysis of the primary structure, hydrophobic cluster analysis, target signal analysis, and phylogenetic analysis all indicate that NbMnSOD1 is dimeric and targeted to the cytosol. NbMnSOD2 seems to have changed more rapidly and is under less evolutionary constraint than NbMnSOD1 suggesting that NbMnSOD2 may function under different conditions or in different tissues of its host rather than simply resulting in an increase in expression. A phylogenetic analysis of MnSOD sequences from eukaryotes, Archaea, and bacteria shows the microsporidial MnSODs to be grouped with the bacteria suggesting a possible horizontal gene transfer.

The origin and early evolution of eukaryotes in the light of phylogenomics

Phylogenomics of eukaryote supergroups suggest a highly complex last common ancestor of eukaryotes and a key role of mitochondrial endosymbiosis in the origin of eukaryotes.

Phylogenetic tree reconstruction accuracy and model fit when proportions of variable sites change across the tree

Commonly used phylogenetic models assume a homogeneous process through time in all parts of the tree. However, it is known that these models can be too simplistic as they do not account for nonhomogeneous lineage-specific properties. In particular, it is now widely recognized that as constraints on sequences evolve, the proportion and positions of variable sites can vary between lineages causing heterotachy. The extent to which this model misspecification affects tree reconstruction is still unknown. Here, we evaluate the effect of changes in the proportions and positions of variable sites on model fit and tree estimation. We consider 5 current models of nucleotide sequence evolution in a Bayesian Markov chain Monte Carlo framework as well as maximum parsimony (MP). We show that for a tree with 4 lineages where 2 nonsister taxa undergo a change in the proportion of variable sites tree reconstruction under the best-fitting model, which is chosen using a relative test, often results in the wrong tree. In this case, we found that an absolute test of model fit is a better predictor of tree estimation accuracy. We also found further evidence that MP is not immune to heterotachy. In addition, we show that increased sampling of taxa that have undergone a change in proportion and positions of variable sites is critical for accurate tree reconstruction.

EEF2 analysis challenges the monophyly of Archaeplastida and Chromalveolata

BACKGROUND

Classification of eukaryotes provides a fundamental phylogenetic framework for ecological, medical, and industrial research. In recent years eukaryotes have been classified into six major supergroups: Amoebozoa, Archaeplastida, Chromalveolata, Excavata, Opisthokonta, and Rhizaria. According to this supergroup classification, Archaeplastida and Chromalveolata each arose from a single plastid-generating endosymbiotic event involving a cyanobacterium (Archaeplastida) or red alga (Chromalveolata). Although the plastids within members of the Archaeplastida and Chromalveolata share some features, no nucleocytoplasmic synapomorphies supporting these supergroups are currently known.

METHODOLOGY/PRINCIPAL FINDINGS

This study was designed to test the validity of the Archaeplastida and Chromalveolata through the analysis of nucleus-encoded eukaryotic translation elongation factor 2 (EEF2) and cytosolic heat-shock protein of 70 kDa (HSP70) sequences generated from the glaucophyte Cyanophora paradoxa, the cryptophytes Goniomonas truncata and Guillardia theta, the katablepharid Leucocryptos marina, the rhizarian Thaumatomonas sp. and the green alga Mesostigma viride. The HSP70 phylogeny was largely unresolved except for certain well-established groups. In contrast, EEF2 phylogeny recovered many well-established eukaryotic groups and, most interestingly, revealed a well-supported clade composed of cryptophytes, katablepharids, haptophytes, rhodophytes, and Viridiplantae (green algae and land plants). This clade is further supported by the presence of a two amino acid signature within EEF2, which appears to have arisen from amino acid replacement before the common origin of these eukaryotic groups.

CONCLUSIONS/SIGNIFICANCE

Our EEF2 analysis strongly refutes the monophyly of the Archaeplastida and the Chromalveolata, adding to a growing body of evidence that limits the utility of these supergroups. In view of EEF2 phylogeny and other morphological evidence, we discuss the possibility of an alternative eukaryotic supergroup.

A mixed branch length model of heterotachy improves phylogenetic accuracy

Evolutionary relationships are typically inferred from molecular sequence data using a statistical model of the evolutionary process. When the model accurately reflects the underlying process, probabilistic phylogenetic methods recover the correct relationships with high accuracy. There is ample evidence, however, that models commonly used today do not adequately reflect real-world evolutionary dynamics. Virtually all contemporary models assume that relatively fast-evolving sites are fast across the entire tree, whereas slower sites always evolve at relatively slower rates. Many molecular sequences, however, exhibit site-specific changes in evolutionary rates, called "heterotachy." Here we examine the accuracy of 2 phylogenetic methods for incorporating heterotachy, the mixed branch length model--which incorporates site-specific rate changes by summing likelihoods over multiple sets of branch lengths on the same tree--and the covarion model, which uses a hidden Markov process to allow sites to switch between variable and invariable as they evolve. Under a variety of simple heterogeneous simulation conditions, the mixed model was dramatically more accurate than homotachous models, which were subject to topological biases as well as biases in branch length estimates. When data were simulated with strong versions of the types of heterotachy observed in real molecular sequences, the mixed branch length model was more accurate than homotachous techniques. Analyses of empirical data sets confirmed that the mixed branch length model can improve phylogenetic accuracy under conditions that cause homotachous models to fail. In contrast, the covarion model did not improve phylogenetic accuracy compared with homotachous models and was sometimes substantially less accurate. We conclude that a mixed branch length approach, although not the solution to all phylogenetic errors, is a valuable strategy for improving the accuracy of inferred trees.

Application of two-dimensional electrophoresis and matrix-assisted laser desorption/ionization time-of-flight mass spectrometry for proteomic analysis of the sexually transmitted parasite Trichomonas vaginalis

Trichomonas vaginalis is a sexually transmitted protozoan parasite that infects the human urogenital tract causing trichomoniasis, a worldwide disease. In this work, a fresh clinical isolate of T. vaginalis was used for study of the protein expression in this species. Two-dimensional gel electrophoresis (2-DE) and MALDI-TOF/TOF mass spectrometry (MS) were employed to create a reference map of soluble proteins in the pH range 4-7. A set of 116 proteins belonging to functional classes expressed in high and low abundance was identified by peptide mass fingerprinting and tandem MS. These identifications corresponded to 67 different proteins, suggesting that post-translational modifications are common phenomena in T. vaginalis. Identified proteins were classified into 16 groups according to biological processes. Among detected proteins we identified the major enzymes involved in both cytosolic and hydrogenosomal metabolic pathways, as well as putative protein targets for new drug design. In addition, this analysis allows validation of previous gene predictions confirming the expression of 15 hypothetical proteins. Finally, the findings here reported represent the first reference proteome map of T. vaginalis and the first steps towards the description of a comprehensive proteome map of this parasite.

A postgenomic view of the heat shock proteins in kinetoplastids

The kinetoplastids Leishmania major, Trypanosoma brucei and Trypanosoma cruzi are causative agents of a diverse spectrum of human diseases: leishmaniasis, sleeping sickness and Chagas' disease, respectively. These protozoa possess digenetic life cycles that involve development in mammalian and insect hosts. It is generally accepted that temperature is a triggering factor of the developmental programme allowing the adaptation of the parasite to the mammalian conditions. The heat shock response is a general homeostatic mechanism that protects cells from the deleterious effects of environmental stresses, such as heat. This response is universal and includes the synthesis of the heat-shock proteins (HSPs). In this review, we summarize the salient features of the different HSP families and describe their main cellular functions. In parallel, we analyse the composition of these families in kinetoplastids according to literature data and our understanding of genome sequence data. The genome sequences of these parasites have been recently completed. The HSP families described here are: HSP110, HSP104, group I chaperonins, HSP90, HSP70, HSP40 and small HSPs. All these families are widely represented in these parasites. In particular, kinetoplastids possess an unprecedented number of members of the HSP70, HSP60 and HSP40 families, suggesting key roles for these HSPs in their biology.

Phylogenetic analyses of parasites in the new millennium

Phylogenetic analysis has changed greatly in the last decade, and the most important themes in that change are reviewed here. Sequence data have become the most common source of phylogenetic information. This means that explicit models for evolutionary processes have been developed in a likelihood context, which allow more realistic data analyses. These models are becoming increasingly complex, both for nucleotides and for amino acid sequences, and so all such models need to be quantitatively assessed for each data set, to find the most appropriate one for use in any particular tree-building analysis. Bayesian analysis has been developed for tree-building and is greatly increasing in popularity. This is because a good heuristic strategy exists, which allows large data sets to be analyzed with complex evolutionary models in a practical time. Perhaps the most disappointing aspect of tree interpretation is the ongoing confusion between rooted and unrooted trees, while the effect of taxon and character sampling is often overlooked when constructing a phylogeny (especially in parasitology). The review finishes with a detailed consideration of the analysis of a multi-gene data set for several dozen taxa of Cryptosporidium (Apicomplexa), illustrating many of the theoretical and practical points highlighted in the review.

Pneumocystis oryctolagisp. nov., an uncultured fungus causing pneumonia in rabbits at weaning: review of current knowledge, and description of a new taxon on genotypic, phylogenetic and phenotypic bases

The genus Pneumocystis comprises noncultivable, highly diversified fungal pathogens dwelling in the lungs of mammals. The genus includes numerous host-species-specific species that are able to induce severe pneumonitis, especially in severely immunocompromised hosts. Pneumocystis organisms attach specifically to type-1 epithelial alveolar cells, showing a high level of subtle and efficient adaptation to the alveolar microenvironment. Pneumocystis species show little difference at the light microscopy level but DNA sequences of Pneumocystis from humans, other primates, rodents, rabbits, insectivores and other mammals present a host-species-related marked divergence. Consistently, selective infectivity could be proven by cross-infection experiments. Furthermore, phylogeny among primate Pneumocystis species was correlated with the phylogeny of their hosts. This observation suggested that cophylogeny could explain both the current distribution of pathogens in their hosts and the speciation. Thus, molecular, ultrastructural and biological differences among organisms from different mammals strengthen the view of multiple species existing within the genus Pneumocystis. The following species were subsequently described: Pneumocystis jirovecii in humans, Pneumocystis carinii and Pneumocystis wakefieldiae in rats, and Pneumocystis murina in mice. The present work focuses on Pneumocystis oryctolagi sp. nov. from Old-World rabbits. This new species has been described on the basis of both biological and phylogenetic species concepts.

Comparison of the evolutionary distances among syngens and sibling species of Paramecium

The morphospecies of the genus Paramecium have several mating type groups, so-called syngens, composed of cells of complementary mating types. The Paramecium aurelia complex is composed of 15 sibling species assigned to the species from the syngen. To increase our understanding of the evolutionary relationships among syngen and sibling species of the genus Paramecium, we investigated the gene sequences of cytosol-type hsp70 from 7 syngens of Paramecium caudatum and 15 sibling species of P. aurelia. Molecular phylogenetic trees indicated that the P. aurelia complex could be divided into four lineages and separated into each sibling species. However, we did not find any obvious genetic distance among syngens of P. caudatum, and they could only be separated into two closely related groups. These results indicated that the concept of syngens in P. caudatum differs quite markedly from that of the P. aurelia complex. In addition, we also discuss the relationships among these species and other species, Paramecium jenningsi and Paramecium multimicronucleatum, which were once classified as varieties of P. aurelia.

The largest subunit of RNA polymerase II from the Glaucocystophyta: functional constraint and short-branch exclusion in deep eukaryotic phylogeny

BACKGROUND

Evolutionary analyses of the largest subunit of RNA polymerase II (RPB1) have yielded important and at times provocative results. One particularly troublesome outcome is the consistent inference of independent origins of red algae and green plants, at odds with the more widely accepted view of a monophyletic Plantae comprising all eukaryotes with primary plastids. If the hypothesis of a broader kingdom Plantae is correct, then RPB1 trees likely reflect a persistent phylogenetic artifact. To gain a better understanding of RNAP II evolution, and the presumed artifact relating to green plants and red algae, we isolated and analyzed RPB1 from representatives of Glaucocystophyta, the third eukaryotic group with primary plastids.

RESULTS

Phylogenetic analyses incorporating glaucocystophytes do not recover a monophyletic Plantae; rather they result in additional conflicts with the most widely held views on eukaryotic relationships. In particular, glaucocystophytes are recovered as sister to several amoebozoans with strong support. A detailed investigation shows that this clade can be explained by what we call "short-branch exclusion," a phylogenetic artifact integrally associated with "long-branch attraction." Other systematic discrepancies observed in RPB1 trees can be explained as phylogenetic artifacts; however, these apparent artifacts also appear in regions of the tree that support widely held views of eukaryotic evolution. In fact, most of the RPB1 tree is consistent with artifacts of rate variation among sequences and co-variation due to functional constraints related to C-terminal domain based RNAP II transcription.

CONCLUSION

Our results reveal how subtle and easily overlooked biases can dominate the overall results of molecular phylogenetic analyses of ancient eukaryotic relationships. Sources of potential phylogenetic artifact should be investigated routinely, not just when obvious "long-branch attraction" is encountered.

Shaping the mitochondrial proteome

Mitochondria are eukaryotic organelles that originated from a single bacterial endosymbiosis some 2 billion years ago. The transition from the ancestral endosymbiont to the modern mitochondrion has been accompanied by major changes in its protein content, the so-called proteome. These changes included complete loss of some bacterial pathways, amelioration of others and gain of completely new complexes of eukaryotic origin such as the ATP/ADP translocase and most of the mitochondrial protein import machinery. This renewal of proteins has been so extensive that only 14-16% of modern mitochondrial proteome has an origin that can be traced back to the bacterial endosymbiont. The rest consists of proteins of diverse origin that were eventually recruited to function in the organelle. This shaping of the proteome content reflects the transformation of mitochondria into a highly specialized organelle that, besides ATP production, comprises a variety of functions within the eukaryotic metabolism. Here we review recent advances in the fields of comparative genomics and proteomics that are throwing light on the origin and evolution of the mitochondrial proteome.

A "chimera" theory on the origin of dicyemid mesozoans: evolution driven by frequent lateral gene transfer from host to parasite

The phylogenetic status of the enigmatic dicyemid mesozoans is still uncertain. Are they primitive multicellular organisms or degenerate triploblastic animals? Presently, the latter view is accepted. A phylogenetic analysis of 18S rDNA sequences placed dicyemids within the animal clade, and this was supported by the discovery of a Hox-type gene with a lophotrochozoan signature sequence. This molecular information suggests that dicyemid mesozoans evolved from an ancestral animal degenerately. Considering their extreme simplicity, which is probably due to parasitism, they might have come from an early embryo via a radical transformation, i.e. neoteny. Irrespective of this molecular information, dicyemid mesozoans retain many protistan-like or extremely primitive features, such as tubular mitochondrial cristae, endocytic ability from the outer surface, and the absence of collagenous tissue, while they do not share noticeable synapomorphy with animals. In addition, the 5S rRNA phylogeny suggests a somewhat closer kinship with protozoan ciliates than with animals. If we accept this clear contradiction, dicyemids should be regarded as a chimera of animals and protistans. Here, we discuss the traditional theory of extreme degeneration via parasitism, and then propose a new "chimera" theory in which dicyemid mesozoans are exposed to a continual flow of genetic information via eating host tissues from the outer surface by endocytosis. Consequently, many of their intrinsic genes have been replaced by host-derived genes through lateral gene transfer (LGT), implying that LGT is a key driving force in the evolution of dicyemid mesozoans.

Phylogenomics of Eukaryotes: Impact of Missing Data on Large Alignments

Resolving the relationships between Metazoa and other eukaryotic groups as well as between metazoan phyla is central to the understanding of the origin and evolution of animals. The current view is based on limited data sets, either a single gene with many species (e.g., ribosomal RNA) or many genes but with only a few species. Because a reliable phylogenetic inference simultaneously requires numerous genes and numerous species, we assembled a very large data set containing 129 orthologous proteins ( approximately 30,000 aligned amino acid positions) for 36 eukaryotic species. Included in the alignments are data from the choanoflagellate Monosiga ovata, obtained through the sequencing of about 1,000 cDNAs. We provide conclusive support for choanoflagellates as the closest relative of animals and for fungi as the second closest. The monophyly of Plantae and chromalveolates was recovered but without strong statistical support. Within animals, in contrast to the monophyly of Coelomata observed in several recent large-scale analyses, we recovered a paraphyletic Coelamata, with nematodes and platyhelminths nested within. To include a diverse sample of organisms, data from EST projects were used for several species, resulting in a large amount of missing data in our alignment (about 25%). By using different approaches, we verify that the inferred phylogeny is not sensitive to these missing data. Therefore, this large data set provides a reliable phylogenetic framework for studying eukaryotic and animal evolution and will be easily extendable when large amounts of sequence information become available from a broader taxonomic range.

LUCA, or just a conserved Archaeon?: Comments on Xue et al. (2003)

In their recent paper, Xue et al. used an unusual technique of rooting the universal phylogenetic tree, which resulted in positioning of the last universal common ancestor within Archaea. The present paper brings some criticisms on the methods and results achieved.

Oral microbial heat-shock proteins and their potential contributions to infections

The oral cavity is a complex ecosystem in which several hundred microbial species normally cohabit harmoniously. However, under certain special conditions, the growth of some micro-organisms with a pathogenic potential is promoted, leading to infections such as dental caries, periodontal disease, and stomatitis. The physiology and pathogenic properties of micro-organisms are influenced by modifications in environmental conditions that lead to the synthesis of specific proteins known as the heat-shock proteins (HSPs). HSPs are families of highly conserved proteins whose main role is to allow micro-organisms to survive under stress conditions. HSPs act as molecular chaperones in the assembly and folding of proteins, and as proteases when damaged or toxic proteins have to be degraded. Several pathological functions have been associated with these proteins. Many HSPs of oral micro-organisms, particularly periodontopathogens, have been identified, and some of their properties-including location, cytotoxicity, and amino acid sequence homology with other HSPs-have been reported. Since these proteins are immunodominant antigens in many human pathogens, studies have recently focused on the potential contributions of HSPs to oral diseases. The cytotoxicity of some bacterial HSPs may contribute to tissue destruction, whereas the presence of common epitopes in host proteins and microbial HSPs may lead to autoimmune responses. Here, we review the current knowledge regarding HSPs produced by oral micro-organisms and discuss their possible contributions to the pathogenesis of oral infections.

Searching factors causing implausible non-monophyly: ssu rDNA phylogeny of Isopoda Asellota (Crustacea: Peracarida) and faster evolution in marine than in freshwater habitats

This contribution addresses two questions: which alignment patterns are causing non-monophyly of the Asellota and what is the phylogenetic history of this group? The Asellota are small benthic crustaceans occurring in most aquatic habitats. In view of the complex morphological apomorphies known for this group, monophyly of the Asellota has never been questioned. Using ssu rDNA sequences of outgroups and of 16 asellote species from fresh water, littoral marine habitats and from deep-sea localities, the early divergence between the lineages in fresh water and in the ocean, and the monophyly of the deep-sea taxon Munnopsidae are confirmed. Relative substitution rates of freshwater species are much lower than in other isopod species, rates being highest in some littoral marine genera (Carpias and Jaera). Furthermore, more sequence sites are variable in marine than in freshwater species, the latter conserve outgroup character states. Monophyly is recovered with parsimony methods, but not with distance and maximum likelihood analyses, which tear apart the marine from the freshwater species. The information content of alignments was studied with spectra of supporting positions. The scarcity of signal (=apomorphic nucleotides) supporting monophyly of the Asellota is attributed to a short stem-line of this group or to erosion of signal in fast evolving marine species. Parametric boostrapping in combination with spectra indicates that a tree model cannot explain the data and that monophyly of the Asellota should not be rejected even though many topologies do not recover this taxon.

Conflicting phylogenetic signals at the base of the metazoan tree

A phylogenetic framework is essential for under-standing the origin and evolution of metazoan development. Despite a number of recent molecular studies and a rich fossil record of sponges and cnidarians, the evolutionary relationships of the early branching metazoan groups to each other and to a putative outgroup, the choanoflagellates, remain uncertain. This situation may be the result of the limited amount of phylogenetic information found in single genes and the small number of relevant taxa surveyed. To alleviate the effect of these analytical factors in the phylogenetic recons-truction of early branching metazoan lineages, we cloned multiple protein-coding genes from two choanoflagellates and diverse sponges, cnidarians, and a ctenophore. Comparisons of sequences for alpha-tubulin, beta-tubulin, elongation factor 2, HSP90, and HSP70 robustly support the hypothesis that choanoflagellates are closely affiliated with animals. However, analyses of single and concatenated amino acid sequences fail to resolve the relationships either between early branching metazoan groups or between Metazoa and choano-flagellates. We demonstrate that variable rates of evolution among lineages, sensitivity of the analyses to taxon selection, and conflicts in the phylogenetic signal contained in different amino acid sequences obscure the phylogenetic associations among the early branching Metazoa. These factors raise concerns about the ability to resolve the phylogenetic history of animals with molecular sequences. A consensus view of animal evolution may require investigations of genome-scale characters.

Mitochondrial connection to the origin of the eukaryotic cell

Phylogenetic evidence is presented that primitively amitochondriate eukaryotes containing the nucleus, cytoskeleton, and endomembrane system may have never existed. Instead, the primary host for the mitochondrial progenitor may have been a chimeric prokaryote, created by fusion between an archaebacterium and a eubacterium, in which eubacterial energy metabolism (glycolysis and fermentation) was retained. A Rickettsia-like intracellular symbiont, suggested to be the last common ancestor of the family Rickettsiaceae and mitochondria, may have penetrated such a host (pro-eukaryote), surrounded by a single membrane, due to tightly membrane-associated phospholipase activity, as do present-day rickettsiae. The relatively rapid evolutionary conversion of the invader into an organelle may have occurred in a safe milieu via numerous, often dramatic, changes involving both partners, which resulted in successful coupling of the host glycolysis and the symbiont respiration. Establishment of a potent energy-generating organelle made it possible, through rapid dramatic changes, to develop genuine eukaryotic elements. Such sequential, or converging, global events could fill the gap between prokaryotes and eukaryotes known as major evolutionary discontinuity.

Assessment of otocephalan and protacanthopterygian concepts in the light of multiple molecular phylogenies

The rise of cladistics in ichthyology has dramatically improved our knowledge of teleostean basal interrelationships. However, some questions have remained open, among them the reliability of the Otocephala, a clade grouping clupeomorphs and ostariophysans, and the relationships of the Esocoidei. These two questions have been investigated in the light of new DNA sequences (from 28S and rhodopsin genes) and sequences from data banks (cytochrome b, 12-16S, 18S, MLL and RAG1). The ability of each of these markers to resolve basal teleostean interrelationships is assessed, and the cytochrome b was not found appropriate. Practical (i.e. different taxonomic samplings) and epistemological grounds led us to perform multiple separated phylogenetic analyses, in order to estimate the reliability of the above clades from their repeatability among trees from independent sequence data. The Otocephala are found monophyletic from most of the datasets; otherwise, they are not significantly contradicted from the others, which exhibit unresolved relationships. We conclude that the evidence provided here favours the sister-group relationship of clupeomorphs and ostariophysans. Morphological evidence including fossils is discussed, concluding that morphological works have not yet provided sufficient data to support this group. Salmonids and esocoids are found sister-groups from every molecular dataset in which these groups were sampled. Based on these convincing results, the Protacanthopterygii of Johnson and Patterson [1] are redefined, including the Esocoidei.

Cytokinesis arrest and nuclear fission in low density populations of trichomonad protozoan

Cell growth of anaerobic protozoan Tritrichomonas foetus was analyzed. This protozoan usually proliferates in extremely high density, but protozoan parasites were dispersed uniformly in F-bouillon medium and cell division stopped temporarily. However, nuclear fission continued and giant polynucleated cells formed. Later, cell division resumed and cells returned to normal form. In conditioned medium, cytokinesis of the dispersed parasites did not stop. Results indicated that T. foetus cells secreted an extracellular factor that influenced cytokinesis.

Ancient phylogenetic relationships

Traditional views on deep evolutionary events have been seriously challenged over the last few years, following the identification of major pitfalls affecting molecular phylogeny reconstruction. Here we describe the principally encountered artifacts, notably long branch attraction, and their causes (i.e., difference in evolutionary rates, mutational saturation, compositional biases). Additional difficulties due to phenomena of biological nature (i.e., lateral gene transfer, recombination, hidden paralogy) are also discussed. Moreover, contrary to common beliefs, we show that the use of rare genomic events can also be misleading and should be treated with the same caution as standard molecular phylogeny. The universal tree of life, as described in most textbooks, is partly affected by tree reconstruction artifacts, e.g. (i) the bacterial rooting of the universal tree of life; (ii) the early emergence of amitochondriate lineages in eukaryotic phylogenies; and (iii) the position of hyperthermophilic taxa in bacterial phylogenies. We present an alternative view of this tree, based on recent evidence obtained from reanalyses of ancient data sets and from novel analyses of large combination of genes.

Pneumocystis carinii f. sp. hominis is not infectious for SCID mice

The infectious power of Pneumocystis carinii f. sp. hominis was explored by inoculating SCID mice intranasally with either P. carinii f. sp. hominis or P. carinii f. sp. muris isolates. Only mice inoculated with mouse parasites developed Pneumocystis pneumonia, as assessed by microscopy and PCR. These results suggest that humans do not contract pneumocystosis from animals.

The chaperonin genes of jakobid and jakobid-like flagellates: implications for eukaryotic evolution

The jakobids are free-living mitochondriate protists that share ultrastructural features with certain amitochondriate groups and possess the most bacterial-like mitochondrial genomes described thus far. Jakobids belong to a diverse group of mitochondriate and amitochondriate eukaryotes, the excavate taxa. The relationships among the various excavate taxa and their relationships to other putative deep-branching protist groups are largely unknown. With the hope of clarifying these issues, we have isolated the cytosolic chaperonin CCTalpha gene from the jakobid Reclinomonas americana (strains 50394 and 50283), the jakobid-like malawimonad Malawimonas jakobiformis, two heteroloboseans (Acrasis rosea and Naegleria gruberi), a euglenozoan (Trypanosoma brucei), and a parabasalid (Monocercomonas sp.). We also amplified the CCTdelta gene from M. jakobiformis. The Reclinomonas and Malawimonas sequences presented here are among the first nuclear protein-coding genes to be described from these organisms. Unlike other putative early diverging protist lineages, a high density of spliceosomal introns was found in the jakobid and malawimonad CCTs-similar to that observed in vertebrate protein-coding genes. An analysis of intron positions in CCT genes from protists, plants, animals, and fungi suggests that many of the intron-sparse or intron-lacking protist lineages may not be primitively so but have lost spliceosomal introns during their evolutionary history. In phylogenetic trees constructed from CCTalpha protein sequences, R. americana (but not M. jakobiformis) shows a weak but consistent affinity for the Heterolobosea and Euglenozoa.

The analysis of 100 genes supports the grouping of three highly divergent amoebae: Dictyostelium, Entamoeba, and Mastigamoeba

The phylogenetic relationships of amoebae are poorly resolved. To address this difficult question, we have sequenced 1,280 expressed sequence tags from Mastigamoeba balamuthi and assembled a large data set containing 123 genes for representatives of three phenotypically highly divergent major amoeboid lineages: Pelobionta, Entamoebidae, and Mycetozoa. Phylogenetic reconstruction was performed on approximately 25,000 aa positions for 30 species by using maximum-likelihood approaches. All well-established eukaryotic groups were recovered with high statistical support, validating our approach. Interestingly, the three amoeboid lineages strongly clustered together in agreement with the Conosa hypothesis [as defined by T. Cavalier-Smith (1998) Biol. Rev. Cambridge Philos. Soc. 73, 203-266]. Two amitochondriate amoebae, the free-living Mastigamoeba and the human parasite Entamoeba, formed a significant sister group to the exclusion of the mycetozoan Dictyostelium. This result suggested that a part of the reductive process in the evolution of Entamoeba (e.g., loss of typical mitochondria) occurred in its free-living ancestors. Applying this inexpensive expressed sequence tag approach to many other lineages will surely improve our understanding of eukaryotic evolution.

Heterotachy, an important process of protein evolution

Because of functional constraints, substitution rates vary among the positions of a protein but are usually assumed to be constant at a given site during evolution. The distribution of the rates across the sequence positions generally fits a Gamma distribution. Models of sequence evolution were accordingly designed and led to improved phylogenetic reconstruction. However, it has been convincingly demonstrated that the evolutionary rate of a given position is not always constant throughout time. We called such within-site rate variations heterotachy (for "different speed" in Greek). Yet, heterotachy was found among homologous sequences of distantly related organisms, often with different functions. In such cases, the functional constraints are likely different, which would explain the different distribution of variable sites. To evaluate the importance of heterotachy, we focused on amino acid sequences of mitochondrial cytochrome b, for which the function is likely the same in all vertebrates. Using 2,038 sequences, we demonstrate that 95% of the variable positions are heterotachous, i.e., underwent dramatic variations of substitution rate among vertebrate lineages. Heterotachy even occurs at small evolutionary scale, and in these cases it is very unlikely to be related to functional changes. Since a large number of sequences are required to efficiently detect heterotachy, the extent of this phenomenon could not be estimated for all proteins yet. It could be as large as for cytochrome b, since this protein is not a peculiar case. The observations made here open several new avenues of research, such as the understanding of the evolution of functional constraints or the improvement of phylogenetic reconstruction methods.

Phylogenetic position of Blastocystis hominis and of stramenopiles inferred from multiple molecular sequence data

Blastocystis hominis, a parasite of the human intestine, has recently been positioned within stramenopiles by the small subunit rRNA phylogeny. To further confirm its phylogenetic position using multiple molecular sequence data, we determined the nucleotide sequences putatively encoding small subunit ribosomal RNA, cytosolic-type 70-kDa heat shock protein, translation elongation factor 2, and the non-catalytic 'B' subunit of vacuolar ATPase of B. hominis (HE87-1 strain). Moreover, we determined the translation elongation factor 2 sequence of an apicomplexan parasite, Plasmodium falciparum, that belongs to alveolates. The maximum likelihood analyses of small subunit rRNA and cytosolic-type 70-kDa heat shock protein clearly demonstrated that B. hominis (HE87-1 strain) is positioned within stramenopiles, being congruent with the previous small subunit rRNA analysis, including the sequences of B. hominis (Nand strain) and a Blastocystis isolate from guinea pig. Although no clear resolution among major eukaryotic groups was obtained by the individual phylogenies based on the four molecules analyzed here, a combined analysis of various molecules, including these, clearly indicated that Blastocystis/stramenopiles are the closest relatives of alveolates.

Maximum-likelihood phylogenetic analysis under a covarion-like model

Here, a model allowing covarion-like evolution of DNA sequences is introduced. In contrast to standard representation of the distribution of evolutionary rates, this model allows the site-specific rate to vary between lineages. This is achieved by adding as few as two parameters to the widely used among-site rate variation model, namely, (1) the proportion of sites undergoing rate changes and (2) the rate of rate change. This model is implemented in the likelihood framework, allowing parameter estimation, comparison of models, and tree reconstruction. An application to ribosomal RNA sequences suggests that covarions (i.e., site-specific rate changes) play an important role in the evolution of these molecules. Neglecting them results in a severe underestimate of the variance of rates across sites. It has, however, little influence on the estimation of ancestral G+C contents obtained from a nonhomogeneous model, or on the resulting inferences about the evolution of thermophyly. This theoretical effort should be useful for the study of protein adaptation, which presumably proceeds in a typical covarion-like manner.

Giardia lamblia expresses a proteobacterial-like DnaK homolog

We identified a novel gene encoding molecular chaperone HSP70 in the amitochondriate parasite Giardia lamblia. The predicted protein is similar to bacterial DnaK and mitochondrial HSP70s. The gene is transcribed and translated at a constant level during trophozoite growth and encystation. Alignment of the sequence with a data set of cytosolic, endoplasmic reticulum (ER), mitochondrial, and DnaK HSP70 homologs indicated that the sequence was extremely divergent and contained insertions unique to giardial HSP70s. Phylogenetic analyses demonstrated that this sequence was distinct from the cytosolic and ER forms and was most similar to proteobacterial and mitochondrial DnaKs. However, a specific relationship with the alpha proteobacterial and mitochondrial sequences was not strongly supported by phylogenetic analyses of this data set, in contrast to similar analyses of cpn60. These data neither confirm nor reject the possibility that this gene is a relic of secondary mitochondrial loss; they leave open the possibility that it was acquired in a separate endosymbiotic event.

Microsporidia: emerging pathogenic protists

Microsporidia are eukaryotic spore forming obligate intracellular protozoan parasites first recognized over 100 years ago. These organisms infect all of the major animal groups and are now recognized as opportunistic pathogens of humans. Microsporidian spores are common in the environment and microsporidia pathogenic to humans have been found in water supplies. The genera Nosema, Vittaforma, Brachiola, Pleistophora, Encephalitozoon, Enterocytozoon, Septata (reclassified to Encephalitozoon) and Trachipleistophora have been found in human infections. These organisms have the smallest known eukaryotic genomes. Microsporidian ribosomal RNA sequences have proven useful as diagnostic tools as well as for phylogenetic analysis. Recent phylogenetic analysis suggests that Microsporidia are related to the fungi. These organisms are defined by the presence of a unique invasion organelle consisting of a single polar tube that coils around the interior of the spore. All microsporidia exhibit the same response to stimuli, that is, the polar tube discharges from the anterior pole of the spore in an explosive reaction. If the polar tube is discharged next to a cell, it can pierce the cell and transfer its sporoplasm into the cell. A technique was developed for the purification of polar tube proteins (PTPs) using differential extraction followed by reverse phase HPLC. This method was used to purify the PTPs from Glugea americanus, Encephalitozoon cuniculi, Enc. hellem and Enc. intestinalis. These PTPs demonstrate conserved characteristics such as solubility, hydrophobicity, mass, proline content and immunologic epitopes. The major PTP gene from Enc. cuniculi and Enc. hellem has been cloned and expressed in vitro. The gene sequences support the importance of ER and in the formation of the polar tube as suggested by morphologic studies. Analysis of the cloned proteins also indicates that secondary structural characteristics are conserved. These characteristics are probably important in the function of this protein during the eversion/assembly of the polar tube and in providing elasticity and resiliency for sporoplasm passage.