Euglena mutabilis exists in a FAB consortium with microbes that enhance cadmium tolerance.

BACKGROUND: Synthetic algal-fungal and algal-bacterial cultures have been investigated as a means to enhance the technological applications of the algae. This inclusion of other microbes has enhanced growth and improved stress tolerance of the algal culture. The goal of the current study was to investigate natural microbial consortia to gain an understanding of the occurrence and benefits of these associations in nature. The photosynthetic protist Euglena mutabilis is often found in association with other microbes in acidic environments with high heavy metal (HM) concentrations. This may suggest that microbial interactions are essential for the protist's ability to tolerate these extreme environments. Our study assessed the Cd tolerance of a natural fungal-algal-bacterial (FAB) association whereby the algae is E. mutabilis. RESULTS: This study provides the first assessment of antibiotic and antimycotic agents on an E. mutabilis culture. The results indicate that antibiotic and antimycotic applications significantly decreased the viability of E. mutabilis cells when they were also exposed to Cd. Similar antibiotic treatments of E. gracilis cultures had variable or non-significant impacts on Cd tolerance. E. gracilis also recovered better after pre-treatment with antibiotics and Cd than did E. mutabilis. The recoveries were assessed by heterotrophic growth without antibiotics or Cd. In contrast, both Euglena species displayed increased chlorophyll production upon Cd exposure. PacBio full-length amplicon sequencing and targeted Sanger sequencing identified the microbial species present in the E. mutabilis culture to be the fungus Talaromyces sp. and the bacterium Acidiphilium acidophilum. CONCLUSION: This study uncovers a possible fungal, algal, and bacterial relationship, what we refer to as a FAB consortium. The members of this consortium interact to enhance the response to Cd exposure. This results in a E. mutabilis culture that has a higher tolerance to Cd than the axenic E. gracilis. The description of this interaction provides a basis for explore the benefits of natural interactions. This will provide knowledge and direction for use when creating or maintaining FAB interactions for biotechnological purposes, including bioremediation.

The enigmatic RNase MRP of kinetoplastids.

The ribonucleoprotein RNase MRP is responsible for the processing of ribosomal RNA precursors. It is found in virtually all eukaryotes that have been examined. In the Euglenozoa, including the genera Euglena, Diplonema and kinetoplastids, MRP RNA and protein subunits have so far escaped detection using bioinformatic methods. However, we now demonstrate that the RNA component is widespread among the Euglenozoa and that these RNAs have secondary structures that conform to the structure of all other phylogenetic groups. In Euglena, we identified the same set of P/MRP protein subunits as in many other protists. However, we failed to identify any of these proteins in the kinetoplastids. This finding poses interesting questions regarding the structure and function of RNase MRP in these species.

Description of Flexiglena gen. nov. and new members of Discoplastis and Euglenaformis (Euglenida).

Environmental sampling in Poland and the United States and phylogenetic analyses based on 567 sequences of four genes (155 sequences of nuclear SSU rDNA, 139 of nuclear LSU rDNA, 135 of plastid-encoded SSU rDNA, and 138 of plastid-encoded LSU rDNA) resulted in description of the new genus Flexiglena, which has been erected by accommodating Euglena variabilis, and enriching the Discoplastis and Euglenaformis genera with five new species. Four of them have joined the Discoplastis genus, currently consisting of six representatives: D. adunca, D. angusta (=Euglena angusta), D. constricta (=Lepocinclis constricta), D. excavata (=E. excavata), D. gasterosteus (=E. gasterosteus), and D. spathirhyncha. One of them has enriched the Euglenaformis genus, currently represented by two species: Euf. chlorophoenicea (= E. chlorophoenicea) and Euf. proxima. For most studied species, the diagnostic descriptions have been emended and epitypes were designated. Furthermore, the emending of Discoplastis and Euglenaformis diagnoses was performed.

Oxygenic Phototrophs Need ζ-Carotene Isomerase (Z-ISO) for Carotene Synthesis: Functional Analysis in Arthrospira and Euglena.

For carotenogenesis, two biosynthetic pathways from phytoene to lycopene are known. Most bacteria and fungi require only phytoene desaturase (PDS, CrtI), whereas land plants require four enzymes: PDS (CrtP), ζ-carotene desaturase (ZDS, CrtQ), ζ-carotene isomerase (Z-ISO) and cis-carotene isomerase (CrtISO, CrtH). The gene encoding Z-ISO has been functionally identified in only two species, Arabidopsis thaliana and Zea mays, and has been little studied in other organisms. In this study, we found that the deduced amino acid sequences of Arthrospira Z-ISO and Euglena Z-ISO have 58% and 62% identity, respectively, with functional Z-ISO from Arabidopsis. We studied the function of Z-ISO genes from the cyanobacterium Arthrospira platensis and eukaryotic microalga Euglena gracilis. The Z-ISO genes of Arthrospira and Euglena were transformed into Escherichia coli strains that produced mainly 9,15,9'-tri-cis-ζ-carotene in darkness. In the resulting E. coli transformants cultured under darkness, 9,9'-di-cis-ζ-carotene was accumulated predominantly as Z-ISO in Arabidopsis. This indicates that the Z-ISO genes were involved in the isomerization of 9,15,9'-tri-cis-ζ-carotene to 9,9'-di-cis-ζ-carotene in darkness. This is the first functional analysis of Z-ISO as a ζ-carotene isomerase in cyanobacteria and eukaryotic microalgae. Green sulfur bacteria and Chloracidobacterium also use CrtP, CrtQ and CrtH for lycopene synthesis as cyanobacteria, but their genomes did not comprise Z-ISO genes. Consequently, Z-ISO is needed in oxygenic phototrophs, whereas it is not found in anoxygenic species.

The Chloroplast Genome of Euglena mutabilis-Cluster Arrangement, Intron Analysis, and Intrageneric Trends.

A comparative analysis of the chloroplast genome of Euglena mutabilis underlined a high diversity in the evolution of plastids in euglenids. Gene clusters in more derived Euglenales increased in complexity with only a few, but remarkable changes in the genus Euglena. Euglena mutabilis differed from other Euglena species in a mirror-inverted arrangement of 12 from 15 identified clusters, making it very likely that the emergence at the base of the genus Euglena, which has been considered a long branch artifact, is truly a probable position. This was corroborated by many similarities in gene arrangement and orientation with Strombomonas and Monomorphina, rendering the genome organization of E. mutabilis in certain clusters as plesiomorphic feature. By RNA analysis exact exon-intron boundaries and the type of the 77 introns identified were mostly determined unambiguously. A detailed intron study of psbC pointed at two important issues: First, the number of introns varied even between species, and no trend from few to many introns could be observed. Second, mat1 was localized in Eutreptiales exclusively in intron 1, and mat2 was not identified. With the emergence of Euglenaceae in most species, a new intron containing mat2 inserted in front of the previous intron 1 and thereby became intron 2 with mat1.

A rare case of plastid protein-coding gene duplication in the chloroplast genome of Euglena archaeoplastidiata (Euglenophyta).

Gene duplication is an important evolutionary process that allows duplicate functions to diverge, or, in some cases, allows for new functional gains. However, in contrast to the nuclear genome, gene duplications within the chloroplast are extremely rare. Here, we present the chloroplast genome of the photosynthetic protist Euglena archaeoplastidiata. Upon annotation, it was found that the chloroplast genome contained a novel tandem direct duplication that encoded a portion of RuBisCO large subunit (rbcL) followed by a complete copy of ribosomal protein L32 (rpl32), as well as the associated intergenic sequences. Analyses of the duplicated rpl32 were inconclusive regarding selective pressures, although it was found that substitutions in the duplicated region, all non-synonymous, likely had a neutral functional effect. The duplicated region did not exhibit patterns consistent with previously described mechanisms for tandem direct duplications, and demonstrated an unknown mechanism of duplication. In addition, a comparison of this chloroplast genome to other previously characterized chloroplast genomes from the same family revealed characteristics that indicated E. archaeoplastidiata was probably more closely related to taxa in the genera Monomorphina, Cryptoglena, and Euglenaria than it was to other Euglena taxa. Taken together, the chloroplast genome of E. archaeoplastidiata demonstrated multiple characteristics unique to the euglenoid world, and has justified the longstanding curiosity regarding this enigmatic taxon.

Wax Ester Fermentation and Its Application for Biofuel Production.

In Euglena cells under anaerobic conditions, paramylon, the storage polysaccharide, is promptly degraded and converted to wax esters. The wax esters synthesized are composed of saturated fatty acids and alcohols with chain lengths of 10-18, and the major constituents are myristic acid and myristyl alcohol. Since the anaerobic cells gain ATP through the conversion of paramylon to wax esters, the phenomenon is named "wax ester fermentation". The wax ester fermentation is quite unique in that the end products, i.e. wax esters, have relatively high molecular weights, are insoluble in water, and accumulate in the cells, in contrast to the common fermentation end products such as lactic acid and ethanol.A unique metabolic pathway involved in the wax ester fermentation is the mitochondrial fatty acid synthetic system. In this system, fatty acid are synthesized by the reversal of ß-oxidation with an exception that trans-2-enoyl-CoA reductase functions instead of acyl-CoA dehydrogenase. Therefore, acetyl-CoA is directly used as a C2 donor in this fatty acid synthesis, and the conversion of acetyl-CoA to malonyl-CoA, which requires ATP, is not necessary. Consequently, the mitochondrial fatty acid synthetic system makes possible the net gain of ATP through the synthesis of wax esters from paramylon. In addition, acetyl-CoA is provided in the anaerobic cells from pyruvate by the action of a unique enzyme, oxygen sensitive pyruvate:NADP+ oxidoreductase, instead of the common pyruvate dehydrogenase multienzyme complex.Wax esters produced by anaerobic Euglena are promising biofuels because myristic acid (C14:0) in contrast to other algal produced fatty acids, such as palmitic acid (C16:0) and stearic acid (C18:0), has a low freezing point making it suitable as a drop-in jet fuel. To improve wax ester production, the molecular mechanisms by which wax ester fermentation is regulated in response to aerobic and anaerobic conditions have been gradually elucidated by identifying individual genes related to the wax ester fermentation metabolic pathway and by comprehensive gene/protein expression analysis. In addition, expression of the cyanobacterial Calvin cycle fructose-1,6-bisphosphatase/sedohepturose-1,7-bisphosphatase, in Euglena provided photosynthesis resulting in increased paramylon accumulation enhancing wax ester production. This chapter will discuss the biochemistry of the wax ester fermentation, recent advances in our understanding of the regulation of the wax ester fermentation and genetic engineering approaches to increase production of wax esters for biofuels.

Arsenic hypertolerance in the protist Euglena mutabilis is mediated by specific transporters and functional integrity maintenance mechanisms.

Arsenic is a toxic metalloid known to cause multiple and severe cellular damages, including lipid peroxidation, protein misfolding, mutagenesis and double and single-stranded DNA breaks. Thus, exposure to this compound is lethal for most organisms but some species such as the photosynthetic protist Euglena mutabilis are able to cope with very high concentrations of this metalloid. Our comparative transcriptomic approaches performed on both an arsenic hypertolerant protist, i.e. E. mutabilis, and a more sensitive one, i.e. E. gracilis, revealed multiple mechanisms involved in arsenic tolerance. Indeed, E. mutabilis prevents efficiently the accumulation of arsenic in the cell through the expression of several transporters. More surprisingly, this protist induced the expression of active DNA reparation and protein turnover mechanisms, which allow E. mutabilis to maintain functional integrity of the cell under challenging conditions. Our observations suggest that this protist has acquired specific functions regarding arsenic and has developed an original metabolism to cope with acid mine drainages-related stresses.

The impact of elevated sulfur and nitrogen levels on cadmium tolerance in Euglena species.

Heavy metal (HM) pollution threatens human and ecosystem health. Current methods for remediating water contaminated with HMs are expensive and have limited effect. Therefore, bioremediation is being investigated as an environmentally and economically viable alternative. Freshwater protists Euglena gracilis and Euglena mutabilis were investigated for their tolerance to cadmium (Cd). A greater increase in cell numbers under Cd stress was noted for E. mutabilis but only E. gracilis showed an increase in Cd tolerance following pre-treatment with elevated concentrations of S or N. To gain insight regarding the nature of the increased tolerance RNA-sequencing was carried out on E. gracilis. This revealed transcript level changes among pretreated cells, and additional differences among cells exposed to CdCl2. Gene ontology (GO) enrichment analysis reflected changes in S and N metabolism, transmembrane transport, stress response, and physiological processes related to metal binding. Identifying these changes enhances our understanding of how these organisms adapt to HM polluted environments and allows us to target development of future pre-treatments to enhance the use of E. gracilis in bioremediation relating to heavy metals.

Were class C iron-containing superoxide dismutases of trypanosomatid parasites initially imported into a complex plastid? A hypothesis based on analyses of their N-terminal targeting signals.

Trypanosomatid parasites possess 2 distinct iron-containing superoxide dismutases (Fe-SODs) designated SODA and SODC, both of which are targeted to their mitochondria. In contrast to SODAs that carry typical mitochondrial transit peptides, SODCs have highly unusual mitochondrial targeting signals. Our analyses clearly show that these pre-sequences are bipartite possessing a signal peptide-like domain followed by a transit peptide-like domain. Consequently, they resemble N-terminal extensions of proteins targeted to multi-membrane plastids, suggesting that trypanosomatids once contained a eukaryotic alga-derived plastid. Further support for this hypothesis comes from striking similarities in length, hydropathy profile, and amino acid composition of SODC pre-sequences to those of Euglena and dinoflagellate plastid proteins. To account for these data, we propose that the Trypanosomatidae initially possessed a gene encoding a mitochondrial Fe-SOD with a classical mitochondrial transit peptide. Before or after plastid acquisition, a gene duplication event gave rise to SODA and SODC. In a subsequent evolutionary step a signal peptide was linked to SODC, enabling its import into the plastid. When the trypanosomatid plastid subsequently was lost, natural selection favoured adaptation of the SODC N-terminal signal as a mitochondrial transit peptide and re-targeting to the mitochondrion.

Differences in planktonic microbial communities associated with three types of macrophyte stands in a shallow lake.

Little is known about how various substances from living and decomposing aquatic macrophytes affect the horizontal patterns of planktonic bacterial communities. Study sites were located within Lake Kolon, which is a freshwater marsh and can be characterised by open-water sites and small ponds with different macrovegetation (Phragmites australis, Nymphea alba and Utricularia vulgaris). Our aim was to reveal the impact of these macrophytes on the composition of the planktonic microbial communities using comparative analysis of environmental parameters, microscopy and pyrosequencing data. Bacterial 16S rRNA gene sequences were dominated by members of phyla Proteobacteria (36%-72%), Bacteroidetes (12%-33%) and Actinobacteria (5%-26%), but in the anoxic sample the ratio of Chlorobi (54%) was also remarkable. In the phytoplankton community, Cryptomonas sp., Dinobryon divergens, Euglena acus and chrysoflagellates had the highest proportion. Despite the similarities in most of the measured environmental parameters, the inner ponds had different bacterial and algal communities, suggesting that the presence and quality of macrophytes directly and indirectly controlled the composition of microbial plankton.

Identification and comparative analysis of the chloroplast alpha-subunit gene of DNA-dependent RNA polymerase from seven Euglena species.

When the sequence of the Euglena gracilis chloroplast genome was reported in 1993 the alpha-subunit gene (rpoA) of RNA polymerase appeared to be missing, based on a comparison of all putative reading frames to the then known rpoA loci. Since there has been a large increase in known rpoA sequences, the question of a Euglena chloroplast rpoA gene was re-examined. A previously described unknown reading frame of 161 codons was found to be part of an rpoA gene split by a single group III intron. This rpoA gene, which is highly variable from species to species, was then isolated and characterized in five other euglenoid species, Euglena anabaena, Euglena granulata, Euglena myxocylindracea, Euglena stellata and Euglena viridis, and in the Astasia longa plastid genome. All seven Euglena rpoA genes have either one or three group III introns. The rpoA gene products in Euglena spp. appear to be the most variable in this gene family when compared to the rpoA gene in other species of bacteria, algae and plants. Additionally, Euglena rpoA proteins lack a C-terminal domain required for interaction with some regulatory proteins, a feature shared only with some chlorophyte green algae. The E.gracilis rpoA gene is the distal cistron of a multigene cluster that includes genes for carbohydrate biosynthesis, photosynthetic electron transport, an antenna complex and ribosomal proteins. This study provides new insights into the transcription system of euglenoid plastids, the organization of the plastid genome, group III intron evolution and euglenoid phylogeny.

Analysis of intergenic spacer transcripts suggests 'read-around' transcription of the extrachromosomal circular rDNA in Euglena gracilis.

We report here the sequence of the 1743 bp intergenic spacer (IGS) that separates the 3'-end of the large subunit ribosomal RNA (rRNA) gene from the 5'-end of the small subunit (SSU) rRNA gene in the circular, extrachromosomal ribosomal DNA (rDNA) of Euglena gracilis. The IGS contains a 277 nt stretch of sequence that is related to a sequence found in ITS 1, an internal transcribed spacer between the SSU and 5.8S rRNA genes. Primer extension analysis of IGS transcripts identified three abundant reverse transcriptase stops that may be analogous to the transcription initiation site (TIS) and two processing sites (A' and A0) that are found in this region in other eukaryotes. Features that could influence processing at these sites include an imperfect palindrome near site A0 and a sequence near site A' that could potentially base pair with U3 small nucleolar RNA. Our identification of the TIS (verified by mung bean nuclease analysis) is considered tentative because we also detected low-abundance transcripts upstream of this site throughout the entire IGS. This result suggests the possibility of 'read-around' transcription, i.e. transcription that proceeds multiple times around the rDNA circle without termination.

Recent horizontal intron transfer to a chloroplast genome.

Evidence is presented for the recent, horizontal transfer of a self-splicing, homing group II intron from a cyanobacteria to the chloroplast genome of Euglena myxocylindracea. The psbA gene of E.myxocylindracea was found to contain a single 2566 nt group II intron with a gene in domain 4 for a 575 amino acid maturase. The predicted secondary structure and tertiary interactions of the group II intron, as well as the derived maturase primary sequence, most closely resemble the homing intron of the cyanobacterium Calothrix and the rnl introns of Porphyra purpurea mitochondria, while being only distantly related to all other Euglena plastid introns and maturases. All main functional domains of the intron-encoded proteins of known homing introns are conserved, including reverse transcriptase domains 1-7, the zinc finger domain and domain X. The close relationship with cyanobacterial introns was confirmed by phylogenetic analysis. Both the full-length psbA intron and a Delta-maturase variant self-splice in vitro in two independent assays. The psbA intron is the first example of a self-splicing chloroplast group II intron from any organism. These results support the conclusion that the psbA intron is the result of a recent horizontal transfer into the E.myxocylindracea chloroplast genome from a cyanobacterial donor and should prompt a reconsideration of horizontal transfer mechanisms to account for the origin of other chloroplast genetic elements.

Solar Radiation Stress in Natural Acidophilic Biofilms of Euglena mutabilis Revealed by Metatranscriptomics and PAM Fluorometry.

The daily photosynthetic performance of a natural biofilm of the extreme acidophilic Euglena mutabilis from Río Tinto (SW, Spain) under full solar radiation was analyzed by means of pulse amplitude-modulated (PAM) fluorescence measurements and metatrascriptomic analysis. Natural E. mutabilis biofilms undergo large-scale transcriptomic reprogramming during midday due to a dynamic photoinhibition and solar radiation stress. Photoinhibition is due to UV radiation and not to light intensity, as revealed by PAM fluorometry analysis. In order to minimize the negative effects of solar radiation, our data supports the presence of a circadian rhythm in this euglenophyte that increases their opportunity to survive. Differential gene expression throughout the day (at 12:00, 20:00 and night) was monitored by massive Illumina parallel sequencing of metatranscriptomic libraries. The transcription pattern was altered in genes involved in Photosystem II stability and repair, UV damaged DNA repair, non-photochemical quenching and oxidative stress, supporting the photoinhibition detected by PAM fluorometry at midday.

Gene Discovery for Synthetic Biology: Exploring the Novel Natural Product Biosynthetic Capacity of Eukaryotic Microalgae.

Eukaryotic microalgae are an incredibly diverse group of organisms whose sole unifying feature is their ability to photosynthesize. They are known for producing a range of potent toxins, which can build up during harmful algal blooms causing damage to ecosystems and fisheries. Genome sequencing is lagging behind in these organisms because of their genetic complexity, but transcriptome sequencing is beginning to make up for this deficit. As more sequence data becomes available, it is apparent that eukaryotic microalgae possess a range of complex natural product biosynthesis capabilities. Some of the genes concerned are responsible for the biosynthesis of known toxins, but there are many more for which we do not know the products. Bioinformatic and analytical techniques have been developed for natural product discovery in bacteria and these approaches can be used to extract information about the products synthesized by algae. Recent analyses suggest that eukaryotic microalgae produce many complex natural products that remain to be discovered.

Localizing Proteins in Fixed Giardia lamblia and Live Cultured Mammalian Cells by Confocal Fluorescence Microscopy.

Confocal fluorescence microscopy and electron microscopy (EM) are complementary methods for studying the intracellular localization of proteins. Confocal fluorescence microscopy provides a rapid and technically simple method to identify the organelle in which a protein localizes but only EM can identify the suborganellular compartment in which that protein is present. Confocal fluorescence microscopy, however, can provide information not obtainable by EM but required to understand the dynamics and interactions of specific proteins. In addition, confocal fluorescence microscopy of cells transfected with a construct encoding a protein of interest fused to a fluorescent protein tag allows live cell studies of the subcellular localization of that protein and the monitoring in real time of its trafficking. Immunostaining methods for confocal fluorescence microscopy are also faster and less involved than those for EM allowing rapid optimization of the antibody dilution needed and a determination of whether protein antigenicity is maintained under fixation conditions used for EM immunogold labeling. This chapter details a method to determine by confocal fluorescence microscopy the intracellular localization of a protein by transfecting the organism of interest, in this case Giardia lamblia, with the cDNA encoding the protein of interest and then processing these organisms for double label immunofluorescence staining after chemical fixation. Also presented is a method to identify the organelle targeting information in the presequence of a precursor protein, in this case the presequence of the precursor to the Euglena light harvesting chlorophyll a/b binding protein of photosystem II precursor (pLHCPII), using live cell imaging of mammalian COS7 cells transiently transfected with a plasmid encoding a pLHCPII presequence fluorescent protein fusion and stained with organelle-specific fluorescent dyes.

Hybridization of bean, spinach, maize and Euglena chloroplast transfer RNAs with homologous and heterologous chloroplast DNAs. An approach to the study of homology between chloroplast tRNAs from various species.

Chloroplast tRNAs from two dicotyledons (spinach and bean), a monocotyledon (maize) and a green alga (Euglena) have been fractionated by two-dimensional gel electrophoresis. The individual tRNAs have been identified, albeled with 125I or 32P, and used in tRNA-DNA hybridization experiments. Spinach chloroplast tRNAs hybridize as well, and maize chloroplast tRNAs almost as well as bean chloroplast tRNAs to bean chloroplast DNA, thus suggesting a high degree of homology between the chloroplast tRNAs from the two dicotyledons and between the tRNAs from the two dicotyledons and those of the monocotyledon. But Euglena total chloroplast tRNA hybridizes very poorly to bean chloroplast DNA, and among the 14 individual tRNAs tested, only one, Euglena chloroplast tRNAPhe, hybridizes to both maize and bean chloroplast DNAs, which is in good agreement with the fact that Euglena and bean chloroplast tRNAsPhe have almost identical primary structures.

Euglenoid flagellates: a multifaceted biotechnology platform.

Euglenoid flagellates are mainly fresh water protists growing in highly diverse environments making them well-suited for a multiplicity of biotechnology applications. Phototrophic euglenids possesses complex chloroplasts of green algal origin bounded by three membranes. Euglena nuclear and plastid genome organization, gene structure and gene expression are distinctly different from other organisms. Our observations on the model organism Euglena gracilis indicate that transcription of both the plastid and nuclear genome is insensitive to environmental changes and that gene expression is regulated mainly at the post-transcriptional level. Euglena plastids have been proposed as a site for the production of proteins and value added metabolites of biotechnological interest. Euglena has been shown to be a suitable protist species to be used for production of several compounds that are used in the production of cosmeceuticals and nutraceuticals, such as α-tocopherol, wax esters, polyunsaturated fatty acids, biotin and tyrosine. The storage polysaccharide, paramylon, has immunostimulatory properties and has shown a promise for biomaterials production. Euglena biomass can be used as a nutritional supplement in aquaculture and in animal feed. Diverse applications of Euglena in environmental biotechnology include ecotoxicological risk assessment, heavy metal bioremediation, bioremediation of industrial wastewater and contaminated water.

Gene duplication and the evolution of photosynthetic reaction center proteins.

We investigate the evolutionary relationships between photosynthetic reaction center proteins (D1, D2, L and M) and demonstrate that the pattern of nucleotide substitution in these is more complicated than has been assumed in previous phylogenetic analyses. We show that there are serious violations of methodological assumptions in previous published studies. We conclude that there is equal support for hypotheses indicating (i) a single gene duplication of an ancestral reaction center protein followed by diversification and (ii) two independent gene duplications giving rise to proteins in oxygenic and anoxygenic systems.