EVOLUTION OF PHACUS (EUGLENOPHYCEAE) AS INFERRED FROM PELLICLE MORPHOLOGY AND SSU rDNA

This research integrates a large morphological data set into a molecular context. Nineteen pellicle characters and 62 states from 13 euglenid taxa were analyzed cladistically. The pellicle morphology of Euglena tripteris (Klebs), Lepocinclis ovata (Conrad), Phacus brachykentron (Pochmann), P. oscillans (Klebs), P. pyrum (Stein), and P. triqueter (Dujardin) is described comprehensively. These data are compared with new information on the pellicle morphology of Euglena acus (Ehrenberg), E. stellata (Mainx), and Peranema trichophorum (Stein) in addition to published data on Entosiphon sulcatum (Dujardin), Euglena gracilis (Klebs), Distigma proteus (Pringsheim), and Petalomonas cantuscygni (Cann and Pennick). Nuclear small subunit (SSU) rDNA sequences provided an independent test for establishing a robust organismal pedigree of the same taxa. A synthetic tree derived from the combined phylogenetic analyses of pellicle morphology and SSU rDNA enabled us to parsimoniously map morphological character states. This approach demonstrated the utility of pellicle morphology for inferring phylogenetic relationships of euglenids and establishing apomorphy-based clade definitions. Three robust clades with unambiguous pellicle-based apomorphies can be recognized within taxa traditionally classified as Phacus : (1) L. ovata and P. pyrum , (2) E. tripteris and P. triqueter , and (3) P. brachykentron and P. oscillans. Taxonomic concerns that emerged from these results are discussed.     

PHYLOGENY OF PHAGOTROPHIC EUGLENIDS (EUGLENOZOA): A MOLECULAR APPROACH BASED ON CULTURE MATERIAL AND ENVIRONMENTAL SAMPLES1

Molecular studies based on small subunit (SSU) rDNA sequences addressing euglenid phylogeny hitherto suffered from the lack of available data about phagotrophic species. To extend the taxon sampling, SSU rRNA genes from species of seven genera of phagotrophic euglenids were investigated. Sequence analyses revealed an increasing genetic diversity among euglenid SSU rDNA sequences compared with other well‐known eukaryotic groups, reflecting an equally broad diversity of morphological characters among euglenid phagotrophs. Phylogenetic inference using standard parsimony and likelihood approaches as well as Bayesian inference and spectral analyses revealed no clear support for euglenid monophyly. Among phagotrophs, monophyly of Petalomonas cantuscygni and Notosolenus ostium, both comprising simple ingestion apparatuses, is strongly supported. A moderately supported clade comprises phototrophic euglenids and primary osmotrophic euglenids together with phagotrophs, exhibiting a primarily flexible pellicle composed of numerous helically arranged strips and a complex ingestion apparatus with two supporting rods and four curved vanes. Comparison of molecular and morphological data is used to demonstrate the difficulties to formulate a hypothesis about how the ingestion apparatus evolved in this group.     

MOLECULAR EVOLUTION OF EUGLENOZOAN PARAXONEMAL ROD GENES par1 AND par2 COINCIDES WITH PHYLOGENETIC RECONSTRUCTION BASED ON SMALL SUBUNIT rDNA DATA1

Emergent flagella of Euglenozoa consist of two prominent structural elements: the axoneme built by microtubules with motor proteins to enable the movement of the flagellum and a highly organized protein structure of unknown function, called the paraxonemal rod (PAR), which consists of two major proteins paralleling the axoneme of euglenid and kinetoplastid emergent flagella. These flagellar structures are considered apomorphic characters of Euglenozoa. We examined the evolution of the genes par1 and par2 encoding the two major proteins, where we could show that these proteins are encoded by two very similar genes found in kinetoplastids and euglenids. The branching pattern indicated a gene duplication before the diversification into euglenids and kinetoplastids. In the clades of the genes, subtrees of euglenid and kinetoplastid monophyla arose. Both genes showed strong genetic diversity with biased GC content at taxon rather than at gene level. We also examined phylogenies inferred from PAR genes that are well in agreement with established small subunit rDNA analyses. Both showed further separation of the euglenid subtree into primary osmotrophs and a phototrophic clade, including secondarily derived osmotrophs.     

Did trypanosomatid parasites have photosynthetic ancestors?

Some molecular phylogenies of plastid-like genes suggest that chloroplasts (the structures responsible for photosynthesis in plants and algae) might have been secondarily lost in trypanosomatid parasites. Chloroplasts are present in some euglenids, which are closely related to trypanosomatids, and it has been argued that chloroplasts arose early in the diversification of the lineage Euglenozoa, to which trypanosomatids and euglenids belong (plastids-early hypothesis). This article reviews how euglenid ultrastructural systems are functionally integrated and phylogenetically correlated. I argue that chloroplast acquisition profoundly altered the structure of certain euglenids, and that the complete absence of these modifications in other euglenozoans is most consistent with their never having had a chloroplast. Ultrastructural evidence suggests that chloroplasts arose relatively recently within a specific subgroup of euglenids and that trypanosomatids are not secondarily non-photosynthetic (plastids-recent hypothesis).     

Protein phylogenies robustly resolve the deep-level relationships within Euglenozoa

The deepest-level relationships amongst Euglenozoa remain poorly resolved, despite a rich history of morphological examination and numerous molecular phylogenetic studies of small subunit ribosomal RNA (SSU rRNA) data. We address this question using two nuclear-encoded proteins, the cytosolic isoforms of heat shock protein 90 (hsp90) and heat shock protein 70 (hsp70). For both proteins we examined sequences from the three primary groups within Euglenozoa (euglenids, diplonemids, and kinetoplastids), and from their close relatives, Heterolobosea. Maximum likelihood (ML) and ML distance analyses of these proteins support a close relationship between diplonemids and kinetoplastids to the exclusion of the euglenid Euglena gracilis. In hsp90 and combined protein analyses bootstrap support is very strong and alternative topologies are generally rejected by 'approximately unbiased' (AU) tests. This result is consistent with recent molecular biological and morphological data, but contradicts early structural accounts and many SSU rRNA analyses that favour a closer relationship between diplonemids and euglenids. However, a re-examination of an important SSU rRNA data set highlights the instability of the inferences from this marker. The protein analyses also suggest that bodonids are paraphyletic, with trypanosomatids grouping with 'clade 2' and 'clade 3' bodonids to the exclusion of 'clade 1' bodonids.     

Phylogenetic Relationships and Morphological Character Evolution of Photosynthetic Euglenids (Excavata) Inferred from Taxon‐rich Analyses of Five Genes

Photosynthetic euglenids acquired chloroplasts by secondary endosymbiosis, which resulted in changes to their mode of nutrition and affected the evolution of their morphological characters. Mapping morphological characters onto a reliable molecular tree could elucidate major trends of those changes. We analyzed nucleotide sequence data from regions of three nuclear‐encoded genes (nSSU, nLSU, hsp90), one chloroplast‐encoded gene (cpSSU) and one nuclear‐encoded chloroplast gene (psbO) to estimate phylogenetic relationships among 59 photosynthetic euglenid species. Our results were consistent with previous works; most genera were monophyletic, except for the polyphyletic genus Euglena, and the paraphyletic genus Phacus. We also analyzed character evolution in photosynthetic euglenids using our phylogenetic tree and eight morphological traits commonly used for generic and species diagnoses, including: characters corresponding to well‐defined clades, apomorphies like presence of lorica and mucilaginous stalks, and homoplastic characters like rigid cells and presence of large paramylon grains. This research indicated that pyrenoids were lost twice during the evolution of phototrophic euglenids, and that mucocysts, which only occur in the genus Euglena, evolved independently at least twice. In contrast, the evolution of cell shape and chloroplast morphology was difficult to elucidate, and could not be unambiguously reconstructed in our analyses.     

Evidence for Transitional Stages in the Evolution of Euglenid Group II Introns and Twintrons in the Monomorphina aenigmatica Plastid Genome

Background

Photosynthetic euglenids acquired their plastid by secondary endosymbiosis of a prasinophyte-like green alga. But unlike its prasinophyte counterparts, the plastid genome of the euglenid Euglena gracilis is riddled with introns that interrupt almost every protein-encoding gene. The atypical group II introns and twintrons (introns-within-introns) found in the E. gracilis plastid have been hypothesized to have been acquired late in the evolution of euglenids, implying that massive numbers of introns may be lacking in other taxa. This late emergence was recently corroborated by the plastid genome sequences of the two basal euglenids, Eutreptiella gymnastica and Eutreptia viridis, which were found to contain fewer introns.

Methodology/Principal Findings

To gain further insights into the proliferation of introns in euglenid plastids, we have characterized the complete plastid genome sequence of Monomorphina aenigmatica, a freshwater species occupying an intermediate phylogenetic position between early and late branching euglenids. The M. aenigmatica UTEX 1284 plastid genome (74,746 bp, 70.6% A+T, 87 genes) contains 53 intron insertion sites, of which 41 were found to be shared with other euglenids including 12 of the 15 twintron insertion sites reported in E. gracilis.

Conclusions

The pattern of insertion sites suggests an ongoing but uneven process of intron gain in the lineage, with perhaps a minimum of two bursts of rapid intron proliferation. We also identified several sites that represent intermediates in the process of twintron evolution, where the external intron is in place, but not the internal one, offering a glimpse into how these convoluted molecular contraptions originate.     

Ultrastructure and Molecular Phylogenetic Position of Heteronema scaphurum: A Eukaryovorous Euglenid with a Cytoproct

Euglenids comprise a distinct clade of flagellates with diverse modes of nutrition, including phagotrophy, osmotrophy and phototrophy. Much of the previous research on euglenids has focused on phototrophic species because of their ecological abundance and significance as indicators for the health of aquatic ecosystems. Although largely understudied, phagotrophic species probably represent the majority of euglenid diversity. Phagotrophic euglenids tend to be either bacterivorous or eukaryovorous and use an elaborate feeding apparatus for capturing prey cells. We characterized the ultrastructure and molecular phylogenetic position of Heteronema scaphurum, a eukaryovorous euglenid collected in freshwater. This species was equipped with a distinct cytoproct through which waste products were eliminated in the form of faecal pellets; a cytoproct has not been reported in any other member of the Euglenida. Heteronema scaphurum also had a novel predatory mode of feeding. The euglenid ensnared and corralled several green algal prey cells (i.e. Chlamydomonas) with hook‐like flagella covered in mucous before engulfing the bundle of prey cells whole. Molecular phylogenetic analyses inferred from small subunit rDNA sequences placed this species with other eukaryovorous euglenids, which was consistent with ultrastructural features associated with the feeding apparatus, flagellar apparatus, extrusomes, and pellicle.     

Unusual mitochondrial genome structures throughout the Euglenozoa

Mitochondrial DNA of Kinetoplastea is composed of different chromosomes, the maxicircle (bearing 'regular' genes) and numerous minicircles (specifying guide RNAs involved in RNA editing). In trypanosomes [Kinetoplastea], DNA circles are compacted into a single dense body, the kinetoplast. This report addresses the question whether multi-chromosome mitochondrial genomes and compacted chromosome organization are restricted to Kinetoplastea or rather occur throughout Euglenozoa, i.e., Kinetoplastea, Euglenida and Diplonemea. To this end, we investigated the diplonemid Rhynchopus euleeides and the euglenids Petalomonas cantuscygni, Peranema trichophorum and Entosiphon sulcatum, using light and electron microscopy and molecular techniques. Our findings together with previously published data show that multi-chromosome mitochondrial genomes prevail across Euglenozoa, while kinetoplast-like mtDNA packaging is confined to trypanosomes.     

Broad taxonomic sampling of mitochondrial cytochrome c oxidase subunit I does not solve the relationships between Rotifera and Acanthocephala

The phylogenetic relationships within Syndermata (Acanthocephala + Rotifera) are still unresolved. Cladistic morphological analyses support monophyly of Rotifera and Eurotatoria (Bdelloidea + Monogononta), while molecular phylogenies of 18S, 28S, COI, hsp82 and EST propose different topologies, with at least six contrasting scenarios. All these phylogenies are characterized by poor taxon sampling; thus, our aim is to solve the relationships within Syndermata sampling as many sequences as possible from one single locus. We reconstructed phylogenetic relationship using more than 1000 sequences of COI. We performed Maximum Likelihood and Bayesian phylogenetic reconstructions on amino acid alignments, using either Gnathostomulida or Platyhelminthes as an outgroup, and then we performed SH tests to provide confidence on the best phylogenetic hypotheses. All four major clades (Acanthocephala, Bdelloidea, Monogononta and Seisonidea) are always highly supported. The basal relationship among the four clades is not consistently resolved by any of the phylogenetic reconstructions; nevertheless, there is a strong support for a clade of Acanthocephala + Bdelloidea from the SH tests, in agreement with other phylogenies from ribosomal genes and EST analyses.     

Utility of NADP-dependent isocitrate dehydrogenase for species-level evolutionary inference in angiosperm phylogeny: a case study in Saltugilia

NADP-dependent isocitrate dehydrogenase is a low-copy nuclear gene family. We have sequenced two regions from an idh gene (idhB) near the 3' terminal end. The first fragment encodes 4 exons and 3 introns and is between approximately 600 and 950 bp in length. The second fragment includes three additional exons and introns and is between approximately 1200 and 1500 bp in length. The phylogenetic utility of the two sequence regions was evaluated in Polemoniaceae with a focus on Saltugilia, an incipient species complex that lacks phylogenetic resolution among these same taxa based on nuclear ribosomal ITS and chloroplast trnL. Multiple sequences from several individuals, multiple individuals from several populations, and multiple populations from all Saltugilia species were sampled to evaluate the taxonomic level at which idhB was useful as a phylogenetic marker in this clade. Phylogenies based on idhB sequences were compared with topological resolution and clade composition in ITS and trnL phylogenies. Phylogenies based on idhB and idhB in combination with ITS and trnL are better resolved than any other phylogenies for Saltugilia published to date, and character evolution within Saltugilia is explored.     

Genomic organization of hsp90 gene family in Arabidopsis

We have isolated six members of the hsp90 gene family from Arabidopsis thaliana. Three genes designated hsp81.2, 81.3 and 81.4 are clustered within a 15 kb genomic region while two of these are 1.5 kb apart in a head-to-head orientation. The deduced amino acid sequence shows that the members can be divided into two types. The hsp81.1, 81.2, 81.3 and 81.4 genes comprise the cytosolic hsp90 type having few introns. However, the hsp88.1 and 89.1 genes comprising the organelle type are composed of 18 or 19 introns. Sequence comparison showed there is high homology among the cytosolic members while there is less homology among the organelle members. The expression of the hsp90 genes and mRNA accumulation in plants and calli is very low at control temperatures and is strongly induced by heat-shock. Arsenite stress strongly stimulates the expression of this gene family.     

Ultrastructural and immunocytochemical investigation of paramylon combined with new 18S rDNA-based secondary structure analysis clarifies phylogenetic affiliation of <Emphasis Type="Italic">Entosiphon sulcatum</Emphasis> (Euglenida: Euglenozoa)

The phylogeny of phagotrophic euglenids is widely based on nuclear small subunit ribosomal DNA (18S rDNA) sequence data, but most analyses suffer from weakness in statistical support regarding the “connecting backbone” between monophyletic clades. Moreover, the position of Entosiphon has remained unclear. Testing the 18S rDNA capability for a phylogeny of phagotrophic euglenids, we isolated sequences of Peranema sp. and Ploeotia edaphica and utilized secondary structure data as a prerequisite for recognition of homologous positions. We found a unique, clade-specific nucleotide substitution in the deduced 18S rDNA helix 44. Since our 18S rDNA phylogenies could only in part resolve positions of phagotrophic lineages, but did not verify that of Entosiphon, we investigated the phagotrophic key taxa Peranema trichophorum, Petalomonas cantuscygni, Ploeotia costata, and Entosiphon sulcatum ultrastructurally. Additionally, we explored the presence or absence of the euglenid reserve carbohydrate paramylon by specific staining with monoclonal anti-β-1,3-glucan antibodies. Paramylon was found to be clearly present in P. trichophorum and E. sulcatum, but was absent in Pt. cantuscygni and Pl. costata. Combined results of our molecular, ultrastructural, and immunocytochemical investigations suggest that Entosiphon sulcatum is the sister taxon of a monophyletic euglenid crown clade, characterized by a helical pellicle, which we propose to rename. This phylogenetic affiliation is confirmed by a clade-specific primary absence of the unique nucleotide substitution in helix 44 and by the common presence of paramylon.     

DNA barcoding in autotrophic euglenids: evaluation of COI and 18s rDNA

Autotrophic euglenids (Euglenophyceae) are a common and abundant group of microbial eukaryotes in freshwater habitats. They have a limited number of features, which can be observed using light microscopy, thus species identification is often problematic. Establishing a barcode for this group is therefore an important step toward the molecular identification of autotrophic euglenids. Based on the literature, we selected verified species and used a plethora of available methods to validate two molecular markers: COI and 18S rDNA (the whole sequence and three fragments separately) as potential DNA barcodes. Analyses of the COI gene were performed based on the data set of 43 sequences (42 obtained in this study) representing 24 species and the COI gene was discarded as a DNA barcode mainly due to a lack of universal primer sites. For 18S rDNA analyses we used a data set containing 263 sequences belonging to 86 taxonomically verified species. We demonstrated that the whole 18S rDNA is too long to be a useful marker, but from the three shorter analyzed variable regions we recommend variable regions V2V3 and V4 of 18S rDNA as autotrophic euglenid barcodes due to their high efficiency (above 95% and 90%, respectively).     

Phylogenetic position of Rhynchopus sp. and Diplonema ambulator as indicated by analyses of euglenozoan small subunit ribosomal DNA

The taxa Rhynchopus Skuja and Diplonema Griessmann were first described as remarkable protists with euglenid affinities. Later on, the placement of Diplonema within the Euglenozoa was confirmed by molecular data. For this study two new sequences were added to the euglenozoan data set. The uncertainly placed Rhynchopus can be identified as a close relative to Diplonema by small subunit ribosomal DNA (SSU rDNA) analysis. The new sequence of Diplonema ambulator is in close relationship to two other Diplonema species. Our molecular analyses clearly support the monophyly of the diplonemids comprising Rhynchopus and Diplonema. Yet the topology at the base of the euglenozoan tree remains unresolved, and especially the monophyly of the euglenids is arguable. SSU rDNA sequence analyses suggest that significantly different GC contents, high mutational saturation in the euglenids, and different evolutionary rates in the euglenozoan clades make it difficult to identify any sister group to the diplonemids.     

Evolution of Distorted Pellicle Patterns in Rigid Photosynthetic Euglenids (Phacus Dujardin)

ABSTRACT. Members of the euglenid genus Phacus are morphologically differentiated from other photosynthetic species by the presence of a rigid cytoskeleton (pellicle) and predominantly dorsoventrally flattened, leaf-shaped cells. In order to better understand the evolutionary history of this lineage, we used scanning electron microscopy to examine patterns of pellicle strips in Phacus acuminatus, Phacus longicauda var. tortus, Phacus triqueter, Phacus segretii, Phacus pleuronectes, Phacus similis, Phacus pusillus, Phacus orbicularis, Phacus warszewiczii , and Discoplastis spathirhyncha , a putative close relative of Phacus and Lepocinclis . Our observations showed that while the earliest diverging species in our analyses, namely P. warszewiczii , has three whorls of exponential reduction, most members of Phacus have clustered patterns of posterior strip reduction that are bilaterally symmetrical distortions of the radially symmetrical "whorled" patterns found in other photosynthetic euglenids. Comparative morphology, interpreted within the context of molecular phylogenetic analyses of combined nuclear small subunit rDNA and partial nuclear large subunit rDNA sequences, demonstrates that clustered patterns of posterior strip reduction arose after the divergence of Phacus from other photosynthetic euglenids and are the result of developmental processes that govern individual strip length. Clustered patterns of pellicle strips in Phacus do not appear to be adaptively significant themselves; they evolved in association with the origin of cell flattening and cell rigidity, which may be adaptations to a planktonic lifestyle.     

Use of the nuclear gene glyceraldehyde 3-phosphate dehydrogenase for phylogeny reconstruction of recently diverged lineages in Mitthyridium (Musci: Calymperaceae)

A portion of the nuclear gene glyceraldehyde 3-phosphate dehydrogenase (gpd) was sequenced in 26 representatives of the paleotropical moss, Mitthyridium, and a group of 20 outgroup taxa to assess its utility for phylogenetic reconstruction compared with the better understood chloroplast markers, rps4 and trnL. Primers based on plant and fungal sequences were designed to amplify gpd in plants universally with the exclusion of fungal contaminants. The piece amplified spanned 4 introns and 3 of 9 exons, based on comparisons with complete sequence from Arabidopsis. Size variation in gpd ranged from 891 to 1007 bp, in part attributable to 6 indels of variable length found within the introns. Intron 6 contributed most of the length variation and contained a variable purine-repeat motif of possible use as a microsatellite. Phylogenetic analyses of the full gpd amplicon yielded well-resolved trees that were in nearly full accord with the trees derived from the cpDNA partitions for analyses of both the ingroup and ingroup + outgroup taxon sets. Pairwise nucleotide substitution rates of gpd were as much as 2.2 times higher than those in rps4 and 2.8 times higher than in trnL. Excision of the introns left suitable numbers of parsimony informative characters and demonstrated that the full gpd amplicon could be compartmentalized to provide resolution for both shallow and deep phylogenetic branches. Exons of gpd were found to behave in a clock-like fashion for the 26 ingroup taxa and select outgroups. In general, gpd was found to hold great promise not only for improving resolution of chloroplast-derived phylogenies, but also for phylogenetic reconstruction of recent, diversifying lineages.     

Structure and Expression of a Heat-Shock Protein 83 Gene of Pharbitis nil

Four cDNA clones representing mRNAs whose levels were affected by a photoperiod that induces flowering in Pharbitis nil were isolated by a differential hybridization screening procedure. The level of mRNAs represented by three clones (12L, 15L, and 17L) increased following a photoperiod that induces flowering and that represented by the fourth clone (clone 27) increased under conditions in which flowering was inhibited. DNA sequence analysis showed that one cDNA, clone 17L, is homologous to members of the 83- to 90-kD heat-shock protein (hsp) gene family. The corresponding gene, hsp83A, was isolated and its DNA sequence was determined. hsp83A encodes a protein that exhibits 70% amino acid identity with Drosophila melanogaster HSP83. The P. nil hsp83A gene contains two introns within the coding region. hsp83A mRNA was not detectable in cotyledons of plants grown in continuous light, but its level increased transiently following a 14-h dark period and reached a maximum 2 h after the lights were turned on. A dramatic increase in the level of hsp83A mRNA was also found 2 h after an end-of-day dark treatment. Genomic Southern blot analysis demonstrated that the P. nil hsp83-90 gene family consists of at least six members, one of which appears to be constitutively expressed in the light.