Subtelomere organization in the genome of the microsporidian Encephalitozoon cuniculi: patterns of repeated sequences and physicochemical signatures.

BACKGROUND: The microsporidian Encephalitozoon cuniculi is an obligate intracellular eukaryotic pathogen with a small nuclear genome (2.9 Mbp) consisting of 11 chromosomes. Although each chromosome end is known to contain a single rDNA unit, the incomplete assembly of subtelomeric regions following sequencing of the genome identified only 3 of the 22 expected rDNA units. While chromosome end assembly remains a difficult process in most eukaryotic genomes, it is of significant importance for pathogens because these regions encode factors important for virulence and host evasion. RESULTS: Here we report the first complete assembly of E. cuniculi chromosome ends, and describe a novel mosaic structure of segmental duplications (EXT repeats) in these regions. EXT repeats range in size between 3.5 and 23.8 kbp and contain four multigene families encoding membrane associated proteins. Twenty-one recombination sites were identified in the sub-terminal region of E. cuniculi chromosomes. Our analysis suggests that these sites contribute to the diversity of chromosome ends organization through Double Strand Break repair mechanisms. The region containing EXT repeats at chromosome extremities can be differentiated based on gene composition, GC content, recombination sites density and chromosome landscape. CONCLUSION: Together this study provides the complete structure of the chromosome ends of E. cuniculi GB-M1, and identifies important factors, which could play a major role in parasite diversity and host-parasite interactions. Comparison with other eukaryotic genomes suggests that terminal regions could be distinguished precisely based on gene content, genetic instability and base composition biais. The diversity of processes assciated with chromosome extremities and their biological consequences, as they are presented in the present study, emphasize the fact that great effort will be necessary in the future to characterize more carefully these regions during whole genome sequencing efforts.

Localization and functionality of microsporidian iron-sulphur cluster assembly proteins.

Microsporidia are highly specialized obligate intracellular parasites of other eukaryotes (including humans) that show extreme reduction at the molecular, cellular and biochemical level. Although microsporidia have long been considered as early branching eukaryotes that lack mitochondria, they have recently been shown to contain a tiny mitochondrial remnant called a mitosome. The function of the mitosome is unknown, because microsporidians lack the genes for canonical mitochondrial functions, such as aerobic respiration and haem biosynthesis. However, microsporidial genomes encode several components of the mitochondrial iron-sulphur (Fe-S) cluster assembly machinery. Here we provide experimental insights into the metabolic function and localization of these proteins. We cloned, functionally characterized and localized homologues of several central mitochondrial Fe-S cluster assembly components for the microsporidians Encephalitozoon cuniculi and Trachipleistophora hominis. Several microsporidial proteins can functionally replace their yeast counterparts in Fe-S protein biogenesis. In E. cuniculi, the iron (frataxin) and sulphur (cysteine desulphurase, Nfs1) donors and the scaffold protein (Isu1) co-localize with mitochondrial Hsp70 to the mitosome, consistent with it being the functional site for Fe-S cluster biosynthesis. In T. hominis, mitochondrial Hsp70 and the essential sulphur donor (Nfs1) are still in the mitosome, but surprisingly the main pools of Isu1 and frataxin are cytosolic, creating a conundrum of how these key components of Fe-S cluster biosynthesis coordinate their function. Together, our studies identify the essential biosynthetic process of Fe-S protein assembly as a key function of microsporidian mitosomes.

Sequencing of the smallest Apicomplexan genome from the human pathogen Babesia microti.

We have sequenced the genome of the emerging human pathogen Babesia microti and compared it with that of other protozoa. B. microti has the smallest nuclear genome among all Apicomplexan parasites sequenced to date with three chromosomes encoding ∼3500 polypeptides, several of which are species specific. Genome-wide phylogenetic analyses indicate that B. microti is significantly distant from all species of Babesidae and Theileridae and defines a new clade in the phylum Apicomplexa. Furthermore, unlike all other Apicomplexa, its mitochondrial genome is circular. Genome-scale reconstruction of functional networks revealed that B. microti has the minimal metabolic requirement for intraerythrocytic protozoan parasitism. B. microti multigene families differ from those of other protozoa in both the copy number and organization. Two lateral transfer events with significant metabolic implications occurred during the evolution of this parasite. The genomic sequencing of B. microti identified several targets suitable for the development of diagnostic assays and novel therapies for human babesiosis.

Cytological, molecular and life cycle characterization of Anostracospora rigaudi n. g., n. sp. and Enterocytospora artemiae n. g., n. sp., two new microsporidian parasites infecting gut tissues of the brine shrimp Artemia.

Two new microsporidia, Anostracospora rigaudi n. g., n. sp., and Enterocytospora artemiae n. g., n. sp. infecting the intestinal epithelium of Artemia parthenogenetica Bowen and Sterling, 1978 and Artemia franciscana Kellogg, 1906 in southern France are described. Molecular analyses revealed the two species belong to a clade of microsporidian parasites that preferentially infect the intestinal epithelium of insect and crustacean hosts. These parasites are morphologically distinguishable from other gut microsporidia infecting Artemia. All life cycle stages have isolated nuclei. Fixed spores measure 1·3×0·7 µm with 5-6 polar tube coils for A. rigaudi and 1·2×0·9 µm with 4 polar tube coils for E. artemiae. Transmission of both species is horizontal, most likely through the ingestion of spores released with the faeces of infected hosts. The minute size of these species, together with their intestinal localization, makes their detection and identification difficult. We developed two species-specific molecular markers allowing each type of infection to be detected within 3-6 days post-inoculation. Using these markers, we show that the prevalence of these microsporidia ranges from 20% to 75% in natural populations. Hence, this study illustrates the usefulness of molecular approaches to study prevalent, but cryptic, infections involving microsporidian parasites of gut tissues.

First molecular characterisation of hydrogenosomes in the protozoan parasite Histomonas meleagridis.

Histomonas meleagridis is a trichomonad species that undergoes a flagellate-to-amoeba transformation during tissue invasion and causes a serious disease in gallinaceous birds (blackhead disease or histomoniasis). Living in the avian cecum, the flagellated form can be grown in vitro in the presence of an ill-defined bacterial flora. Its cytoplasm harbours numerous spherical bodies which structurally resemble hydrogenosomes. To test whether these organelles may be involved in anaerobic metabolism, we undertook the identification of H. meleagridis genes encoding some potentially conserved hydrogenosomal enzymes. The strategy was based on several PCR amplification steps using primers designed from available sequences of the phylogenetically-related human parasite Trichomonas vaginalis. We first obtained a C-terminal sequence of an iron-hydrogenase homologue (Hm_HYD) with typical active site signatures (H-cluster domain). Immunoelectron microscopy with anti-Hm_HYD polyclonal antibodies showed specific gold labelling of electron-dense organelles, thus confirming their hydrogenosomal nature. The whole genes encoding a malic enzyme (Hm_ME) and the alpha-subunit of a succinyl coenzyme A synthetase (Hm_alpha-SCS) were then identified. Short N-terminal presequences for hydrogenosomal targeting were predicted in both proteins. Anti-Hm_ME and anti-Hm_alpha-SCS antisera provided immunofluorescence staining patterns of H. meleagridis cytoplasmic granules similar to those observed with anti-Hm_HYD antiserum or mAb F5.2 known to react with T. vaginalis hydrogenosomes. Hm_ME, Hm_alpha-SCS and Hm_HYD were also detected as reactive bands on immunoblots of proteins from purified hydrogenosomes. Interestingly, anti-Hm_alpha-SCS staining of the cell surface in non-permeabilised parasites suggests a supplementary role for SCS in cytoadherence, as previously demonstrated in T. vaginalis.

The complement of protein kinases of the microsporidium Encephalitozoon cuniculi in relation to those of Saccharomyces cerevisiae and Schizosaccharomyces pombe.

BACKGROUND: Microsporidia, parasitic fungi-related eukaryotes infecting many cell types in a wide range of animals (including humans), represent a serious health threat in immunocompromised patients. The 2.9 Mb genome of the microsporidium Encephalitozoon cuniculi is the smallest known of any eukaryote. Eukaryotic protein kinases are a large superfamily of enzymes with crucial roles in most cellular processes, and therefore represent potential drug targets. We report here an exhaustive analysis of the E. cuniculi genomic database aimed at identifying and classifying all protein kinases of this organism with reference to the kinomes of two highly-divergent yeast species, Saccharomyces cerevisiae and Schizosaccharomyces pombe. RESULTS: A database search with a multi-level protein kinase family hidden Markov model library led to the identification of 29 conventional protein kinase sequences in the E. cuniculi genome, as well as 3 genes encoding atypical protein kinases. The microsporidian kinome presents striking differences from those of other eukaryotes, and this minimal kinome underscores the importance of conserved protein kinases involved in essential cellular processes. Approximately 30% of its kinases are predicted to regulate cell cycle progression while another approximately 28% have no identifiable homologues in model eukaryotes and are likely to reflect parasitic adaptations. E. cuniculi lacks MAP kinase cascades and almost all protein kinases that are involved in stress responses, ion homeostasis and nutrient signalling in the model fungi S. cerevisiae and S. pombe, including AMPactivated protein kinase (Snf1), previously thought to be ubiquitous in eukaryotes. A detailed database search and phylogenetic analysis of the kinomes of the two model fungi showed that the degree of homology between their kinomes of approximately 85% is much higher than that previously reported. CONCLUSION: The E. cuniculi kinome is by far the smallest eukaryotic kinome characterised to date. The difficulty in assigning clear homology relationships for nine out of the twentynine microsporidian conventional protein kinases despite its compact genome reflects the phylogenetic distance between microsporidia and other eukaryotes. Indeed, the E. cuniculi genome presents a high proportion of genes in which evolution has been accelerated by up to four-fold. There are no orthologues of the protein kinases that constitute MAP kinase pathways and many other protein kinases with roles in nutrient signalling are absent from the E. cuniculi kinome. However, orthologous kinases can nonetheless be identified that correspond to members of the yeast kinomes with roles in some of the most fundamental cellular processes. For example, E. cuniculi has clear orthologues of virtually all the major conserved protein kinases that regulate the core cell cycle machinery (Aurora, Polo, DDK, CDK and Chk1). A comprehensive comparison of the homology relationships between the budding and fission yeast kinomes indicates that, despite an estimated 800 million years of independent evolution, the two model fungi share approximately 85% of their protein kinases. This will facilitate the annotation of many of the as yet uncharacterised fission yeast kinases, and also those of novel fungal genomes.

Functional and evolutionary analysis of a eukaryotic parasitic genome.

The DNA sequences of the 11 linear chromosomes of the approximately 2.9 Mbp genome of Encephalitozoon cuniculi, an obligate intracellular parasite of mammals, include approximately 2000 putative protein-coding genes. The compactness of this genome is associated with the length reduction of various genes. Essential functions are dependent on a minimal set of genes. Phylogenetic analysis supports the hypotheses that microsporidia are related to fungi and have retained a mitochondrion-derived organelle, the mitosome.

Occurrence of the microsporidian parasite Nucleospora salmonis in four species of salmonids from the Massif Central of France.

Nucleospora salmonis (Hedrick, Groff et Baxa, 1991), an intranuclear microsporidian parasite of marine and freshwater fish, causes diseases mainly in salmonid species. Losses have been reported in stocks of salmonid fish reared in the region of Auvergne (France). The cause of chronic mortalities in the local host species raised in aquaculture and destined for supplementation of the river system Loire-Allier was examined. The presence of N. salmonis was confirmed by PCR and histology in Salmo salar L. previously and in newly investigated salmonid species, Salmo salar, Salmo trutta fario L., Thymallus thymallus (L.) and Salvelinus alpinus (L.), present in European streams. The infection by N. salmonis was consistent in all cases with characteristic symptoms of the disease in deceased or moribund fish. The small subunit ribosomal DNA from N. salmonis was partially sequenced and compared to previously characterised N. salmonis isolates. As a result, a genotype, or clonal entity, was attributed to N. salmonis among Atlantic salmon found along the Northern Atlantic coastal lines and other salmonid species co-inhabiting or co-cultivated in the Auvergne region.

Proteolytic activity in Encephalitozoon cuniculi sporogonial stages: predominance of metallopeptidases including an aminopeptidase-P-like enzyme.

A fraction enriched in spore precursor cells (sporoblasts) of the microsporidian Encephalitozoon cuniculi, an intracellular parasite of mammals, was obtained by Percoll gradient centrifugation. Soluble extracts of these cells exhibited proteolytic activity towards azocasein, with an alkaline optimum pH range (9-10). Prevalence of some metallopeptidases was supported by the stimulating effect of Ca2+, Mg2+, Mn2+ and Zn2+ ions, and inhibition by two chelating agents (EDTA and 1,10-phenanthroline), a thiol reductant (dithiothreitol) and two aminopeptidase inhibitors (bestatin and apstatin). Zymographic analysis revealed four caseinolytic bands at about 76, 70, 55 and 50 kDa. Mass spectrometry of tryptic peptides from one-dimensional gel slices identified a cytosol (leucine) aminopeptidase homologue (M17 family) in 50-kDa band and an enzyme similar to aminopeptidase P (AP-P) of cytosolic type (M24B subfamily) in 70-kDa band. Multiple sequence alignments showed conservation of critical residues for catalysis and metal binding. A long insertion in a common position was found in AP-P sequences from E. cuniculi and Nosema locustae, an insect-infecting microsporidian. The expression of cytosolic AP-P in sporogonial stages of microsporidia may suggest a key role in the attack of proline-containing peptides as a prerequisite to long-duration biosynthesis of structural proteins destined to the sporal polar tube.

The putative chitin deacetylase of Encephalitozoon cuniculi: a surface protein implicated in microsporidian spore-wall formation.

Microsporidia are fungal-like unicellular eukaryotes which develop as obligate intracellular parasites. They differentiate into resistant spores that are protected by a thick cell wall composed of glycoproteins and chitin. Despite an extensive description of the fibrillar structure of this wall, very little is known about its protein components and deposit mechanisms. In this study on the human pathogen Encephalitozoon cuniculi, we identify by mass spectrometry the target of polyclonal antibodies previously raised against a 33-kDa protein located at the outer face of the parasite plasma membrane. This 254-amino acid protein is encoded by the ECU11_0510 open reading frame and presents two isoforms of 33 and 55 kDa. Sequence analysis supports an assignment to the polysaccharide deacetylase family with a suspected chitin deacetylase activity (EcCDA). As demonstrated by TEM studies, EcCDA is present at the plasma membrane of the early stages of E. cuniculi life-cycle. At the sporoblast stage, the enzyme accumulates especially in paramural bodies which are convolutions of the plasma membrane opened to the wall. The identification of an EcCDA homologue in the insect parasite Antonospora locustae (ex Nosema locustae) suggests a widespread distribution of this enzyme among Microsporidia. This characterization of a new microsporidian surface protein creates new perspectives to understand spore wall formation and spore resistance.

Description of a xenoma-inducing microsporidian, Microgemma tincae n. sp., parasite of the teleost fish Symphodus tinca from Tunisian coasts.

A xenoma-inducing microsporidian species was found to infect the liver of the teleost fish, peacock wrasse Symphodus (Crenilabrus) tinca. Minimal estimates of the prevalence of the parasite in fishes caught along Tunisian coasts were as high as 43 % for Bizerte samples (over 2 yr) and 72% for Monastir samples (over 3 yr). Developmental stages were dispersed within a xenoma structure that was bounded only by the plasma membrane of the hypertrophic host cell. Ultrastructural features support allocation to the genus Microgemma Ralphs and Matthews, 1986. Meronts were multinucleate plasmodia and were surrounded by rough endoplasmic reticulum (RER) of the host cell. Merogonic plasmodia developed into sporogonic plasmodia, with loss of the RER interface. Sporogony was polysporoblastic. Ovocylindrical spores (3.6 x 1.2 microm) harbored a lamellar polaroplast and a polar tube that was coiled 9 times. Spore features and host specificity led us to propose a new species, Microgemma tincae. The conversion of M. tincae xenomas into well-visible cyst structures or granulomas reflected an efficient host response involving the infiltration of phagocytic cells, degradation of various parasite stages and formation of a thick fibrous wall. The small subunit rDNA gene of M. tincae was partially sequenced. Phylogenetic analysis confirms the placement within the family Tetramicriidae represented by the genera Tetramicra and Microgemma.

Post-genomics of microsporidia, with emphasis on a model of minimal eukaryotic proteome: a review.

The genome sequence of the microsporidian parasite Encephalitozoon cuniculi Levaditi, Nicolau et Schoen, 1923 contains about 2,000 genes that are representative of a non-redundant potential proteome composed of 1,909 protein chains. The purpose of this review is to relate some advances in the characterisation of this proteome through bioinformatics and experimental approaches. The reduced diversity of the set of E. cuniculi proteins is perceptible in all the compilations of predicted domains, orthologs, families and superfamilies, available in several public databases. The phyletic patterns of orthologs for seven eukaryotic organisms support an extensive gene loss in the fungal clade, with additional deletions in E. cuniculi. Most microsporidial orthologs are the smallest ones among eukaryotes, justifying an interest in the use of these compacted proteins to better discriminate between essential and non-essential regions. The three components of the E. cuniculi mRNA capping apparatus have been especially well characterized and the three-dimensional structure of the cap methyltransferase has been elucidated following the crystallisation of the microsporidial enzyme Ecm1. So far, our mass spectrometry-based analyses of the E. cuniculi spore proteome has led to the identification of about 170 proteins, one-quarter of these having no clearly predicted function. Immunocytochemical studies are in progress to determine the subcellular localisation of microsporidia-specific proteins. Post-translational modifications such as phosphorylation and glycosylation are expected to be soon explored.

Chromosomal composition of the genome in the monomorphic diplokaryotic microsporidium Paranosema grylli: analysis by two-dimensional pulsed-field gel electrophoresis.

The molecular karyotype of Paranosema grylli Sokolova, Seleznev, Dolgikh et Issi, 1994, a monomorphic diplokaryotic microsporidium, comprises numerous bright and faint bands of nonstoichiometric staining intensity. Restriction analysis of chromosomal DNAs by "karyotype and restriction display" 2-D PFGE has demonstrated that the complexity of molecular karyotype of P. grylli is related to the pronounced length polymorphism of-homologous chromosomes. The background of this phenomenon is discussed in the context of ploidy state, reproductive strategy and population structure in this microsporidium. We propose that the remarkable size variation between homologous chromosomes in P. grylli may be a consequence of ectopic recombination at the chromosome extremities.

Fine analysis of the Pneumocystis carinii f. sp. carinii genome by two-dimensional pulsed-field gel electrophoresis.

Pneumocystis carinii is a general designation for a group of unusual unicellular fungal parasites responsible of pneumopathy in animal hosts. Divided into several subgroups termed the 'special forms', P. carinii is prone to an extensive karyotype variation. In previous studies, the nuclear genome of these organisms has been considered to be haploid and a set of 16 chromosomes has been assigned to P. carinii f. sp. carinii, a special form known to infect rats. We report the analysis of the genome of an isolate representative of the karyotype 1 of this special form, using two-dimensional pulsed-field gel electrophoresis procedures. The 'karyotype and restriction display' (KARD) fingerprints indicated the presence of 17 different chromosomes. The haploid genome size was estimated to be 8.4 Mbp. Some homologous chromosomes were distinguished on the basis of a single restriction fragment length polymorphism, which raises the possibility of a diploid nucleus. A restriction map of the chromosome 15, characterized by two homologues with a size difference of 7 kb, was constructed. Hybridization data indicated that insertion/deletion events may have occurred within subtelomeric regions which carry genes encoding the major surface glycoprotein (MSG) of Pneumocystis.

The microsporidian polar tube: evidence for a third polar tube protein (PTP3) in Encephalitozoon cuniculi.

The invasion strategy used by microsporidia is primarily related to spore germination. Small differentiated spores of these fungi-related parasites inject their contents into target cells through the lumen of a rapidly extruded polar tube, as a prerequisite to obligate intracellular development. Previous studies in Encephalitozoon species that infect mammals have identified two major antigenic polar tube proteins (PTP1 and PTP2) which are predicted to contribute to the high tensile strength of the polar tube via an assembly process dependent on disulfide linkages. By immunoscreening of a cDNA library, we found that a novel PTP is encoded by a single transcription unit (3990 bp) located on the chromosome XI of E. cuniculi. PTP3 is predicted to be synthesized as a 1256-amino acid precursor with a cleavable signal peptide. The mature protein lacks cysteine residue and its large acidic core is flanked by highly basic N- and C-terminal regions. Immunolocalization data indicated that PTP3 is involved in the sporoblast-to-spore polar tube biogenesis. A transcriptional up-regulation during sporogony is supported by a strong increase in the relative amount of Ecptp mRNAs within host cells sampled at late post-infection times. To begin to explore polar tube-associated protein interactions, spore proteins were extracted in the presence of SDS and dithiothreitol then incubated with a chemical cross-linker (DSP or sulfo-EGS). A large multimeric complex was formed and shown to contain PTP1, PTP2 and PTP3 with a few other proteins. PTP3 is hypothesized to play a role in the control of the polar tube extrusion as part of a specific response to ionic stimuli.

Comparative analysis of sequences encoding ABC systems in the genome of the microsporidian Encephalitozoon cuniculi.

Microsporidia are amitochondriate eukaryotic microbes with fungal affinities and a common status of obligate intracellular parasites. A set of 13 potential genes encoding ATP-binding cassette (ABC) systems was identified in the fully sequenced genome of Encephalitozoon cuniculi. Our analyses of multiple alignments, phylogenetic trees and conserved motifs support a distribution of E. cuniculi ABC systems within only four subfamilies. Six half transporters are homologous to the yeast ATM1 mitochondrial protein, a finding which is in agreement with the hypothesis of a cryptic mitochondrion-derived compartment playing a role in the synthesis and transport of Fe-S clusters. Five half transporters are similar to the human ABCG1 and ABCG2 proteins, involved in regulation of lipid trafficking and anthracyclin resistance respectively. Two proteins with duplicated ABC domains are clearly candidate to non-transport ABC systems: the first is homologous to mammalian RNase L inhibitor and the second to the yeast translation initiation regulator GCN20. An unusual feature of ABC systems in E. cuniculi is the lack of homologs of P-glycoprotein and other ABC transporters which are involved in multiple drug resistance in a large number of eukaryotic microorganisms.

The genes encoding cAMP-dependent protein kinase catalytic subunit homologues of the microsporidia Encephalitozoon intestinalis and E. cuniculi: molecular characterisation and phylogenetic analysis.

A gene encoding a protein kinase was identified by homology-based PCR amplification in Encephalitozoon intestinalis, a microsporidian parasite pathogenic to humans, and its orthologue has been identified by database mining in the genome of the related species E. cuniculi, whose sequence has been recently published. Phylogenetic analysis revealed that the proteins encoded by these genes are homologues of the cAMP-dependent protein kinase catalytic subunits (PKAc). Southern blot analysis indicated that the EiPKAc gene is present in two copies in the E. intestinalis genome, whereas the E. cuniculi orthologue (EcPKAc) is a single copy gene. RT-PCR data showed that the EiPKAc gene is expressed in at least one of the intracellular stages during infection of the mammalian host cell by E. intestinalis.

Spraguea lophii (Microsporidia) parasite of the teleost fish, Lophius piscatorius from Tunisian coasts: evidence for an extensive chromosome length polymorphism.

A microsporidian of the genus Spraguea was found parasitizing the nervous tissues of Lophius piscatorius collected from various localities in the Mediterranean coastal areas of Tunisia. The tissue localization, the infection focus aspect and sporal dimorphism are characteristics of Spraguea lophii species. Molecular data based on partial sequence of SSUrRNA encoding gene shows few nucleotide polymorphisms, compared to all described Spraguea isolates. Molecular karyotype obtained on pulsed field gel electrophoresis (1D-PFGE) shows a profile with 14 stained bands in the range of 230-880 kbp and a genome size estimated to 6.700 kbp. The rare cutter endonuclease MluI KARD 2-D-PFGE fingerprint shows an extensive chromosome length polymorphism, but the number of chromosome is unchanged and consists of 15 different molecules. The extensive chromosome length polymorphism is associated to a reduced number of genetic events.

Genome sequence of the stramenopile Blastocystis, a human anaerobic parasite.

BACKGROUND: Blastocystis is a highly prevalent anaerobic eukaryotic parasite of humans and animals that is associated with various gastrointestinal and extraintestinal disorders. Epidemiological studies have identified different subtypes but no one subtype has been definitively correlated with disease. RESULTS: Here we report the 18.8 Mb genome sequence of a Blastocystis subtype 7 isolate, which is the smallest stramenopile genome sequenced to date. The genome is highly compact and contains intriguing rearrangements. Comparisons with other available stramenopile genomes (plant pathogenic oomycete and diatom genomes) revealed effector proteins potentially involved in the adaptation to the intestinal environment, which were likely acquired via horizontal gene transfer. Moreover, Blastocystis living in anaerobic conditions harbors mitochondria-like organelles. An incomplete oxidative phosphorylation chain, a partial Krebs cycle, amino acid and fatty acid metabolisms and an iron-sulfur cluster assembly are all predicted to occur in these organelles. Predicted secretory proteins possess putative activities that may alter host physiology, such as proteases, protease-inhibitors, immunophilins and glycosyltransferases. This parasite also possesses the enzymatic machinery to tolerate oxidative bursts resulting from its own metabolism or induced by the host immune system. CONCLUSIONS: This study provides insights into the genome architecture of this unusual stramenopile. It also proposes candidate genes with which to study the physiopathology of this parasite and thus may lead to further investigations into Blastocystis-host interactions.