The new phylogeny of eukaryotes

Genome content reorganization in the non-model ciliate Chilodonella uncinata: insights into nuclear architecture, DNA content, and chromosome fragmentation during macronuclear development

Ciliates are a model lineage for studies of genome architecture given their unusual genome structures. All ciliates have both somatic macronuclei (MAC) and germline micronuclei (MIC), both of which develop from a zygotic nucleus following sex (i.e., conjugation). Nuclear developmental stages are not well documented among non-model ciliates, including Chilodonella uncinata (class Phyllopharyngea), the focus of our work. Here, we characterize nuclear architecture and genome dynamics in C. uncinata by combining 4',6-diamidino-2-phenylindole (DAPI) staining and fluorescence in situ hybridization (FISH) techniques with confocal microscopy. We developed a telomere probe for staining, which alongside DAPI allows for the identification of fragmented somatic chromosomes among the total DNA in the nuclei. We quantify both total DNA and telomere-bound signals from more than 250 nuclei sampled from 116 individual cells, and analyze changes in DNA content and nuclear architecture across Chilodonella's nuclear life cycle. Specifically, we find that MAC developmental stages in the ciliate C. uncinata are different from those reported from other ciliate species. These data provide insights into nuclear dynamics during development and enrich our understanding of genome evolution in non-model ciliates.

IMPORTANCE

Ciliates are a clade of diverse single-celled eukaryotic microorganisms that contain at least one somatic macronucleus (MAC) and germline micronucleus (MIC) within each cell/organism. Ciliates rely on complex genome rearrangements to generate somatic genomes from a zygotic nucleus. However, the development of somatic nuclei has only been documented for a few model ciliate genera, including Paramecium, Tetrahymena, and Oxytricha. Here, we study the MAC developmental process in the non-model ciliate, C. uncinata. We analyze both total DNA and the generation of gene-sized somatic chromosomes using a laser scanning confocal microscope to describe C. uncinata's nuclear life cycle. We show that DNA content changes dramatically during their life cycle and in a manner that differs from previous studies on model ciliates. Our study expands knowledge of genome dynamics in ciliates and among eukaryotes more broadly.

Genome content in the non-model ciliate Chilodonella uncinata: insights into nuclear architecture, gene-sized chromosomes among the total DNA in their somatic macronuclei during their development

Ciliates are a model lineage for studies of genome architecture given their unusual genome structures. All ciliates have both somatic macronuclei (MAC) and germline micronuclei (MIC), both of which develop from a zygotic nucleus following sex (i.e., conjugation). Nuclear developmental stages are not as well explored among non-model ciliate genera, including Chilodonella uncinata (Class- Phyllopharyngea), the focus of our work. Here, we characterize nuclear architecture and genome dynamics in C. uncinata by combining DAPI (4',6-diamidino-2-phenylindole) staining and fluorescence in situ hybridization (FISH) techniques with confocal microscopy. We developed a telomere probe for staining alongside DAPI, which allows for the identification of fragmented somatic chromosomes among the total DNA in the nuclei. We quantify both total DNA and telomere-bound signals to explore changes in DNA content and chromosome maturation across Chilodonella's nuclear life cycle. Specifically, we find that MAC developmental stages in the ciliate C. uncinata are different than the data reported from other ciliate species. These data provide insights into nuclear dynamics during nuclear development and enrich our understanding of genome evolution in non-model ciliates.

Molecular evidence for the evolution of the eukaryotic mitochondrial arginyl-tRNA synthetase from the prokaryotic suborder Cystobacterineae

The evolutionary origin of the family of eukaryotic aminoacyl-tRNA synthetases that are essential to all living organisms is a matter of debate. In order to shed molecular light on the ancient source of arginyl-tRNA synthetase, a total of 1347 eukaryotic arginyl-tRNA synthetase sequences were mined from databases and analyzed. Their multiple sequence alignment reveals a signature sequence that is characteristic of the nuclear-encoded enzyme, which is imported into mitochondria. Using this molecular beacon, the origins of this gene can be traced to modern prokaryotes. In this way, a previous phylogenetic analysis linking Myxococcus to the emergence of the eukaryotic mitochondrial arginyl-tRNA synthetase is supported by the unique existence of the molecular signature within the suborder Cystobacterineae that includes Myxococcus.

Phylogenomic Insights into Animal Evolution

Animals make up only a small fraction of the eukaryotic tree of life, yet, from our vantage point as members of the animal kingdom, the evolution of the bewildering diversity of animal forms is endlessly fascinating. In the century following the publication of Darwin's Origin of Species, hypotheses regarding the evolution of the major branches of the animal kingdom - their relationships to each other and the evolution of their body plans - was based on a consideration of the morphological and developmental characteristics of the different animal groups. This morphology-based approach had many successes but important aspects of the evolutionary tree remained disputed. In the past three decades, molecular data, most obviously primary sequences of DNA and proteins, have provided an estimate of animal phylogeny largely independent of the morphological evolution we would ultimately like to understand. The molecular tree that has evolved over the past three decades has drastically altered our view of animal phylogeny and many aspects of the tree are no longer contentious. The focus of molecular studies on relationships between animal groups means, however, that the discipline has become somewhat divorced from the underlying biology and from the morphological characteristics whose evolution we aim to understand. Here, we consider what we currently know of animal phylogeny; what aspects we are still uncertain about and what our improved understanding of animal phylogeny can tell us about the evolution of the great diversity of animal life.

Classification and expression analyses of homeobox genes from Dictyostelium discoideum

Homeobox genes are compared between genomes in an attempt to understand the evolution of animal development. The ability of the protist, Dictyostelium discoideum, to shift between uni- and multicellularity makes this group ideal for studying the genetic changes that may have occurred during this transition. We present here the first genome-wide classification and comparative genomic analysis of the 14 homeobox genes present in D. discoideum. Based on the structural alignment of the homeodomains, they can be broadly divided into TALE and non-TALE classes. When individual homeobox genes were compared with members of known class or family, we could further classify them into 3 groups, namely, TALE, OTHER and NOVEL classes, but no HOX family was found. The 5 members of TALE class could be further divided into PBX, PKNOX, IRX and CUP families; 4 homeobox genes classified as NOVEL did not show any similarity to any known homeobox genes; while the remaining 5 were classified as OTHERS as they did show certain degree of similarity to few known homeobox genes. No unique RNA expression pattern during development of D. discoideum emerged for members of an individual group. Putative promoter analysis revealed binding sites for few homeobox transcription factors among many probable factors.

On the age of eukaryotes: evaluating evidence from fossils and molecular clocks

Our understanding of the phylogenetic relationships among eukaryotic lineages has improved dramatically over the few past decades thanks to the development of sophisticated phylogenetic methods and models of evolution, in combination with the increasing availability of sequence data for a variety of eukaryotic lineages. Concurrently, efforts have been made to infer the age of major evolutionary events along the tree of eukaryotes using fossil-calibrated molecular clock-based methods. Here, we review the progress and pitfalls in estimating the age of the last eukaryotic common ancestor (LECA) and major lineages. After reviewing previous attempts to date deep eukaryote divergences, we present the results of a Bayesian relaxed-molecular clock analysis of a large dataset (159 proteins, 85 taxa) using 19 fossil calibrations. We show that for major eukaryote groups estimated dates of divergence, as well as their credible intervals, are heavily influenced by the relaxed molecular clock models and methods used, and by the nature and treatment of fossil calibrations. Whereas the estimated age of LECA varied widely, ranging from 1007 (943-1102) Ma to 1898 (1655-2094) Ma, all analyses suggested that the eukaryotic supergroups subsequently diverged rapidly (i.e., within 300 Ma of LECA). The extreme variability of these and previously published analyses preclude definitive conclusions regarding the age of major eukaryote clades at this time. As more reliable fossil data on eukaryotes from the Proterozoic become available and improvements are made in relaxed molecular clock modeling, we may be able to date the age of extant eukaryotes more precisely.

The eukaryotic tree of life from a global phylogenomic perspective

Molecular phylogenetics has revolutionized our knowledge of the eukaryotic tree of life. With the advent of genomics, a new discipline of phylogenetics has emerged: phylogenomics. This method uses large alignments of tens to hundreds of genes to reconstruct evolutionary histories. This approach has led to the resolution of ancient and contentious relationships, notably between the building blocks of the tree (the supergroups), and allowed to place in the tree enigmatic yet important protist lineages for understanding eukaryote evolution. Here, I discuss the pros and cons of phylogenomics and review the eukaryotic supergroups in light of earlier work that laid the foundation for the current view of the tree, including the position of the root. I conclude by presenting a picture of eukaryote evolution, summarizing the most recent progress in assembling the global tree.

Epigenetics of ciliates

Research using ciliates revealed early examples of epigenetic phenomena and continues to provide novel findings. These protozoans maintain separate germline and somatic nuclei that carry transcriptionally silent and active genomes, respectively. Examining the differences in chromatin within distinct nuclei of Tetrahymena identified histone variants and established that transcriptional regulators act by modifying histones. Formation of somatic nuclei requires both transcriptional activation of silent chromatin and large-scale DNA elimination. This somatic genome remodeling is directed by homologous RNAs, acting with an RNA interference (RNAi)-related machinery. Furthermore, the content of the parental somatic genome provides a homologous template to guide this genome restructuring. The mechanisms regulating ciliate DNA rearrangements reveal the surprising power of homologous RNAs to remodel the genome and transmit information transgenerationally.

Biochemistry and evolution of anaerobic energy metabolism in eukaryotes

Major insights into the phylogenetic distribution, biochemistry, and evolutionary significance of organelles involved in ATP synthesis (energy metabolism) in eukaryotes that thrive in anaerobic environments for all or part of their life cycles have accrued in recent years. All known eukaryotic groups possess an organelle of mitochondrial origin, mapping the origin of mitochondria to the eukaryotic common ancestor, and genome sequence data are rapidly accumulating for eukaryotes that possess anaerobic mitochondria, hydrogenosomes, or mitosomes. Here we review the available biochemical data on the enzymes and pathways that eukaryotes use in anaerobic energy metabolism and summarize the metabolic end products that they generate in their anaerobic habitats, focusing on the biochemical roles that their mitochondria play in anaerobic ATP synthesis. We present metabolic maps of compartmentalized energy metabolism for 16 well-studied species. There are currently no enzymes of core anaerobic energy metabolism that are specific to any of the six eukaryotic supergroup lineages; genes present in one supergroup are also found in at least one other supergroup. The gene distribution across lineages thus reflects the presence of anaerobic energy metabolism in the eukaryote common ancestor and differential loss during the specialization of some lineages to oxic niches, just as oxphos capabilities have been differentially lost in specialization to anoxic niches and the parasitic life-style. Some facultative anaerobes have retained both aerobic and anaerobic pathways. Diversified eukaryotic lineages have retained the same enzymes of anaerobic ATP synthesis, in line with geochemical data indicating low environmental oxygen levels while eukaryotes arose and diversified.

Zygotic expression of the double-stranded RNA binding motif protein Drb2p is required for DNA elimination in the ciliate Tetrahymena thermophila

Double-stranded RNA binding motif (DSRM)-containing proteins play many roles in the regulation of gene transcription and translation, including some with tandem DSRMs that act in small RNA biogenesis. We report the characterization of the genes for double-stranded RNA binding proteins 1 and 2 (DRB1 and DRB2), two genes encoding nuclear proteins with tandem DSRMs in the ciliate Tetrahymena thermophila. Both proteins are expressed throughout growth and development but exhibit distinct peaks of expression, suggesting different biological roles. In support of this, we show that expression of DRB2 is essential for vegetative growth while DRB1 expression is not. During conjugation, Drb1p and Drb2p localize to distinct nuclear foci. Cells lacking all DRB1 copies are able to produce viable progeny, although at a reduced rate relative to wild-type cells. In contrast, cells lacking germ line DRB2 copies, which thus cannot express Drb2p zygotically, fail to produce progeny, arresting late into conjugation. This arrest phenotype is accompanied by a failure to organize the essential DNA rearrangement protein Pdd1p into DNA elimination bodies and execute DNA elimination and chromosome breakage. These results implicate zygotically expressed Drb2p in the maturation of these nuclear structures, which are necessary for reorganization of the somatic genome.

Rooting the eukaryotic tree with mitochondrial and bacterial proteins

By exploiting the large body of genome data and the considerable progress in phylogenetic methodology, recent phylogenomic studies have provided new insights into the relationships among major eukaryotic groups. However, confident placement of the eukaryotic root remains a major challenge. This is due to the large evolutionary distance separating eukaryotes from their closest relatives, the Archaea, implying a weak phylogenetic signal and strong long-branch attraction artifacts. Here, we apply a new approach to the rooting of the eukaryotic tree by using a subset of genomic information with more recent evolutionary origin-mitochondrial sequences, whose closest relatives are α-Proteobacteria. For this, we identified and assembled a data set of 42 mitochondrial proteins (mainly encoded by the nuclear genome) and performed Bayesian and maximum likelihood analyses. Taxon sampling includes the recently sequenced Thecamonas trahens, a member of the phylogenetically elusive Apusozoa. This data set confirms the relationships of several eukaryotic supergroups seen before and places the eukaryotic root between the monophyletic "unikonts" and "bikonts." We further show that T. trahens branches sister to Opisthokonta with significant statistical support and question the bikont/excavate affiliation of Malawimonas species. The mitochondrial data set developed here (to be expanded in the future) constitutes a unique alternative means in resolving deep eukaryotic relationships.

Phylogenomics-based reconstruction of protozoan species tree

We have developed a semi-automatic methodology to reconstruct the phylogenetic species tree in Protozoa, integrating different phylogenetic algorithms and programs, and demonstrating the utility of a supermatrix approach to construct phylogenomics-based trees using 31 universal orthologs (UO). The species tree obtained was formed by three major clades that were related to three groups of data: i) Species containing at least 80% of UO (25/31) in the concatenated multiple alignment or supermatrix, this clade was called C1, ii) Species containing between 50%–79% (15–24/31) of UO called C2, and iii) Species containing less than 50% (1–14/31) of UO called C3. C1 was composed by only protozoan species, C2 was composed by species related to Protozoa, and C3 was composed by some species of C1 (Protozoa) and C2 (related to Protozoa). Our phylogenomics-based methodology using a supermatrix approach proved to be reliable with protozoan genome data and using at least 25 UO, suggesting that (a) the more UO used the better, (b) using the entire UO sequence or just a conserved block of it for the supermatrix produced similar phylogenomic trees.

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.

Alpha-tubulin and small subunit rRNA phylogenies of peritrichs are congruent and do not support the clustering of mobilids and sessilids (Ciliophora, Oligohymenophorea)

Peritrich ciliates have been traditionally subdivided into two orders, Sessilida and Mobilida within the subclass Peritrichia. However, all the existing small subunit (SSU) rRNA phylogenetic trees showed that the sessilids and mobilids did not branch together. To shed some light on this disagreement, we tested whether or not the classic Peritrichia is a monophyletic group by assessing the reliability of the SSU rRNA phylogeny in terms of congruency with alpha-tubulin phylogeny. For this purpose, we obtained 10 partial alpha-tubulin sequences from peritrichs and built phylogenetic trees based on alpha-tubulin nucleotide and amino acid data. A phylogenetic tree from the alpha-tubulin and SSU rRNA genes in combination was also constructed and compared with that from the SSU rRNA gene using a similar species sampling. Our results show that the mobilids and sessilids are consistently separated in all trees, which reinforces the idea that the peritrichs do not constitute a monophyletic group. However, in all alpha-tubulin gene trees, the urceolariids and trichodiniids do not group together, suggested mobilids may not be a monophyletic group.

Molecular Data are Transforming Hypotheses on the Origin and Diversification of Eukaryotes

The explosion of molecular data has transformed hypotheses on both the origin of eukaryotes and the structure of the eukaryotic tree of life. Early ideas about the evolution of eukaryotes arose through analyses of morphology by light microscopy and later electron microscopy. Though such studies have proven powerful at resolving more recent events, theories on origins and diversification of eukaryotic life have been substantially revised in light of analyses of molecular data including gene and, increasingly, whole genome sequences. By combining these approaches, progress has been made in elucidating both the origin and diversification of eukaryotes. Yet many aspects of the evolution of eukaryotic life remain to be illuminated.

Can comprehensive background knowledge be incorporated into substitution models to improve phylogenetic analyses? A case study on major arthropod relationships

BACKGROUND

Whenever different data sets arrive at conflicting phylogenetic hypotheses, only testable causal explanations of sources of errors in at least one of the data sets allow us to critically choose among the conflicting hypotheses of relationships. The large (28S) and small (18S) subunit rRNAs are among the most popular markers for studies of deep phylogenies. However, some nodes supported by this data are suspected of being artifacts caused by peculiarities of the evolution of these molecules. Arthropod phylogeny is an especially controversial subject dotted with conflicting hypotheses which are dependent on data set and method of reconstruction. We assume that phylogenetic analyses based on these genes can be improved further i) by enlarging the taxon sample and ii) employing more realistic models of sequence evolution incorporating non-stationary substitution processes and iii) considering covariation and pairing of sites in rRNA-genes.

RESULTS

We analyzed a large set of arthropod sequences, applied new tools for quality control of data prior to tree reconstruction, and increased the biological realism of substitution models. Although the split-decomposition network indicated a high noise content in the data set, our measures were able to both improve the analyses and give causal explanations for some incongruities mentioned from analyses of rRNA sequences. However, misleading effects did not completely disappear.

CONCLUSION

Analyses of data sets that result in ambiguous phylogenetic hypotheses demand for methods, which do not only filter stochastic noise, but likewise allow to differentiate phylogenetic signal from systematic biases. Such methods can only rely on our findings regarding the evolution of the analyzed data. Analyses on independent data sets then are crucial to test the plausibility of the results. Our approach can easily be extended to genomic data, as well, whereby layers of quality assessment are set up applicable to phylogenetic reconstructions in general.

Support for the Coelomata clade of animals from a rigorous analysis of the pattern of intron conservation

Many intron positions are conserved in varying subsets of eukaryotic genomes and, consequently, comprise a potentially informative class of phylogenetic characters. Roy and Gilbert developed a method of phylogenetic reconstruction using the patterns of intron presence-absence in eukaryotic genes and, applying this method to the analysis of animal phylogeny, obtained support for an Ecdysozoa clade (Roy SW, Gilbert W. 2005. Resolution of a deep animal divergence by the pattern of intron conservation. Proc Natl Acad Sci USA. 102:4403-4408). The critical assumption in the method was the independence of intron loss in different branches of the phylogenetic tree. Here, this assumption is refuted by showing that the branch-specific intron loss rates are strongly correlated. We show that different tree topologies are obtained, in each case with a significant statistical support, when different subsets of intron positions are analyzed. The analysis of the conserved intron positions supports the Coelomata topology, that is, a clade comprised of arthropods and chordates, whereas the analysis of more variable intron positions favors the Ecdysozoa topology, that is, a clade of arthropods and nematodes. We show, however, that the support for Ecdysozoa is fully explained by parallel loss of introns in nematodes and arthropods, a factor that does not contribute to the analysis of the conserved introns. The developed procedure for the identification and analysis of conserved introns and other characters with minimal or no homoplasy is expected to be useful for resolving many hard phylogenetic problems.

The Biological Big Bang model for the major transitions in evolution

BACKGROUND

Major transitions in biological evolution show the same pattern of sudden emergence of diverse forms at a new level of complexity. The relationships between major groups within an emergent new class of biological entities are hard to decipher and do not seem to fit the tree pattern that, following Darwin's original proposal, remains the dominant description of biological evolution. The cases in point include the origin of complex RNA molecules and protein folds; major groups of viruses; archaea and bacteria, and the principal lineages within each of these prokaryotic domains; eukaryotic supergroups; and animal phyla. In each of these pivotal nexuses in life's history, the principal "types" seem to appear rapidly and fully equipped with the signature features of the respective new level of biological organization. No intermediate "grades" or intermediate forms between different types are detectable. Usually, this pattern is attributed to cladogenesis compressed in time, combined with the inevitable erosion of the phylogenetic signal.

HYPOTHESIS

I propose that most or all major evolutionary transitions that show the "explosive" pattern of emergence of new types of biological entities correspond to a boundary between two qualitatively distinct evolutionary phases. The first, inflationary phase is characterized by extremely rapid evolution driven by various processes of genetic information exchange, such as horizontal gene transfer, recombination, fusion, fission, and spread of mobile elements. These processes give rise to a vast diversity of forms from which the main classes of entities at the new level of complexity emerge independently, through a sampling process. In the second phase, evolution dramatically slows down, the respective process of genetic information exchange tapers off, and multiple lineages of the new type of entities emerge, each of them evolving in a tree-like fashion from that point on. This biphasic model of evolution incorporates the previously developed concepts of the emergence of protein folds by recombination of small structural units and origin of viruses and cells from a pre-cellular compartmentalized pool of recombining genetic elements. The model is extended to encompass other major transitions. It is proposed that bacterial and archaeal phyla emerged independently from two distinct populations of primordial cells that, originally, possessed leaky membranes, which made the cells prone to rampant gene exchange; and that the eukaryotic supergroups emerged through distinct, secondary endosymbiotic events (as opposed to the primary, mitochondrial endosymbiosis). This biphasic model of evolution is substantially analogous to the scenario of the origin of universes in the eternal inflation version of modern cosmology. Under this model, universes like ours emerge in the infinite multiverse when the eternal process of exponential expansion, known as inflation, ceases in a particular region as a result of false vacuum decay, a first order phase transition process. The result is the nucleation of a new universe, which is traditionally denoted Big Bang, although this scenario is radically different from the Big Bang of the traditional model of an expanding universe. Hence I denote the phase transitions at the end of each inflationary epoch in the history of life Biological Big Bangs (BBB).

CONCLUSION

A Biological Big Bang (BBB) model is proposed for the major transitions in life's evolution. According to this model, each transition is a BBB such that new classes of biological entities emerge at the end of a rapid phase of evolution (inflation) that is characterized by extensive exchange of genetic information which takes distinct forms for different BBBs. The major types of new forms emerge independently, via a sampling process, from the pool of recombining entities of the preceding generation. This process is envisaged as being qualitatively different from tree-pattern cladogenesis.

Intron dynamics in ribosomal protein genes

The role of spliceosomal introns in eukaryotic genomes remains obscure. A large scale analysis of intron presence/absence patterns in many gene families and species is a necessary step to clarify the role of these introns. In this analysis, we used a maximum likelihood method to reconstruct the evolution of 2,961 introns in a dataset of 76 ribosomal protein genes from 22 eukaryotes and validated the results by a maximum parsimony method. Our results show that the trends of intron gain and loss differed across species in a given kingdom but appeared to be consistent within subphyla. Most subphyla in the dataset diverged around 1 billion years ago, when the “Big Bang” radiation occurred. We speculate that spliceosomal introns may play a role in the explosion of many eukaryotes at the Big Bang radiation.

A type-1 metacaspase from Acanthamoeba castellanii

The complete sequence of a type-1 metacaspase from Acanthamoeba castellanii is reported comprising 478 amino acids. The metacaspase was recovered from an expression library using sera specific for membrane components implicated in stimulating encystation. A central domain of 155 amino acid residues contains the Cys/His catalytic dyad and is the most conserved region containing at least 30 amino acid identities in all metacaspases. The Acanthamoeba castellanii metacaspase has the most proline-rich N-terminus so far reported in type-1 metacaspases with over 40 prolines in the first 150 residues. Ala-Pro-Pro is present 11 times. Phylogenies constructed using only the conserved proteolytic domains or the complete sequences show identical branching patterns, differing only in the rates of change.

Phylum-wide analysis of SSU rDNA reveals deep phylogenetic relationships among nematodes and accelerated evolution toward crown Clades

Inference of evolutionary relationships between nematodes is severely hampered by their conserved morphology, the high frequency of homoplasy, and the scarcity of phylum-wide molecular data. To study the origin of nematode radiation and to unravel the phylogenetic relationships between distantly related species, 339 nearly full-length small-subunit rDNA sequences were analyzed from a diverse range of nematodes. Bayesian inference revealed a backbone comprising 12 consecutive dichotomies that subdivided the phylum Nematoda into 12 clades. The most basal clade is dominated by the subclass Enoplia, and members of the order Triplonchida occupy positions most close to the common ancestor of the nematodes. Crown Clades 8-12, a group formerly indicated as "Secernentea" that includes Caenorhabditis elegans and virtually all major plant and animal parasites, show significantly higher nucleotide substitution rates than the more basal Clades 1-7. Accelerated substitution rates are associated with parasitic lifestyles (Clades 8 and 12) or short generation times (Clades 9-11). The relatively high substitution rates in the distal clades resulted in numerous autapomorphies that allow in most cases DNA barcode-based species identification. Teratocephalus, a genus comprising terrestrial bacterivores, was shown to be most close to the starting point of Secernentean radiation. Notably, fungal feeding nematodes were exclusively found basal to or as sister taxon next to the 3 groups of plant parasitic nematodes, namely, Trichodoridae, Longidoridae, and Tylenchomorpha. The exclusive common presence of fungivorous and plant parasitic nematodes supports a long-standing hypothesis that states that plant parasitic nematodes arose from fungivorous ancestors.

An Empirical Assessment of Long-Branch Attraction Artefacts in Deep Eukaryotic Phylogenomics

In the context of exponential growing molecular databases, it becomes increasingly easy to assemble large multigene data sets for phylogenomic studies. The expected increase of resolution due to the reduction of the sampling (stochastic) error is becoming a reality. However, the impact of systematic biases will also become more apparent or even dominant. We have chosen to study the case of the long-branch attraction artefact (LBA) using real instead of simulated sequences. Two fast-evolving eukaryotic lineages, whose evolutionary positions are well established, microsporidia and the nucleomorph of cryptophytes, were chosen as model species. A large data set was assembled (44 species, 133 genes, and 24,294 amino acid positions) and the resulting rooted eukaryotic phylogeny (using a distant archaeal outgroup) is positively misled by an LBA artefact despite the use of a maximum likelihood-based tree reconstruction method with a complex model of sequence evolution. When the fastest evolving proteins from the fast lineages are progressively removed (up to 90%), the bootstrap support for the apparently artefactual basal placement decreases to virtually 0%, and conversely only the expected placement, among all the possible locations of the fast-evolving species, receives increasing support that eventually converges to 100%. The percentage of removal of the fastest evolving proteins constitutes a reliable estimate of the sensitivity of phylogenetic inference to LBA. This protocol confirms that both a rich species sampling (especially the presence of a species that is closely related to the fast-evolving lineage) and a probabilistic method with a complex model are important to overcome the LBA artefact. Finally, we observed that phylogenetic inference methods perform strikingly better with simulated as opposed to real data, and suggest that testing the reliability of phylogenetic inference methods with simulated data leads to overconfidence in their performance. Although phylogenomic studies can be affected by systematic biases, the possibility of discarding a large amount of data containing most of the nonphylogenetic signal allows recovering a phylogeny that is less affected by systematic biases, while maintaining a high statistical support.

Comparison of the Yeast Proteome to Other Fungal Genomes to Find Core Fungal Genes

The purpose of this research was to search for evolutionarily conserved fungal sequences to test the hypothesis that fungi have a set of core genes that are not found in other organisms, as these genes may indicate what makes fungi different from other organisms. By comparing 6355 predicted or known yeast (Saccharomyces cerevisiae) genes to the genomes of 13 other fungi using Standalone TBLASTN at an e-value <1E-5, a list of 3340 yeast genes was obtained with homologs present in at least 12 of 14 fungal genomes. By comparing these common fungal genes to complete genomes of animals (Fugu rubripes, Caenorhabditis elegans), plants (Arabidopsis thaliana, Oryza sativa), and bacteria (Agrobacterium tumefaciens, Xylella fastidiosa), a list of common fungal genes with homologs in these plants, animals, and bacteria was produced (938 genes), as well as a list of exclusively fungal genes without homologs in these other genomes (60 genes). To ensure that the 60 genes were exclusively fungal, these were compared using TBLASTN to the major sequence databases at GenBank: NR (nonredundant), EST (expressed sequence tags), GSS (genome survey sequences), and HTGS (unfinished high-throughput genome sequences). This resulted in 17 yeast genes with homologs in other fungal genomes, but without known homologs in other organisms. These 17 core, fungal genes were not found to differ from other yeast genes in GC content or codon usage patterns. More intensive study is required of these 17 genes and other common fungal genes to discover unique features of fungi compared to other organisms.

The diversity of dolichol-linked precursors to Asn-linked glycans likely results from secondary loss of sets of glycosyltransferases

The vast majority of eukaryotes (fungi, plants, animals, slime mold, and euglena) synthesize Asn-linked glycans (Alg) by means of a lipid-linked precursor dolichol-PP-GlcNAc2Man9Glc3. Knowledge of this pathway is important because defects in the glycosyltransferases (Alg1-Alg12 and others not yet identified), which make dolichol-PP-glycans, lead to numerous congenital disorders of glycosylation. Here we used bioinformatic and experimental methods to characterize Alg glycosyltransferases and dolichol-PP-glycans of diverse protists, including many human pathogens, with the following major conclusions. First, it is demonstrated that common ancestry is a useful method of predicting the Alg glycosyltransferase inventory of each eukaryote. Second, in the vast majority of cases, this inventory accurately predicts the dolichol-PP-glycans observed. Third, Alg glycosyltransferases are missing in sets from each organism (e.g., all of the glycosyltransferases that add glucose and mannose are absent from Giardia and Plasmodium). Fourth, dolichol-PP-GlcNAc2Man5 (present in Entamoeba and Trichomonas) and dolichol-PP- and N-linked GlcNAc2 (present in Giardia) have not been identified previously in wild-type organisms. Finally, the present diversity of protist and fungal dolichol-PP-linked glycans appears to result from secondary loss of glycosyltransferases from a common ancestor that contained the complete set of Alg glycosyltransferases.

Mitochondria of Protists

Over the past several decades, our knowledge of the origin and evolution of mitochondria has been greatly advanced by determination of complete mitochondrial genome sequences. Among the most informative mitochondrial genomes have been those of protists (primarily unicellular eukaryotes), some of which harbor the most gene-rich and most eubacteria-like mitochondrial DNAs (mtDNAs) known. Comparison of mtDNA sequence data has provided insights into the radically diverse trends in mitochondrial genome evolution exhibited by different phylogenetically coherent groupings of eukaryotes, and has allowed us to pinpoint specific protist relatives of the multicellular eukaryotic lineages (animals, plants, and fungi). This comparative genomics approach has also revealed unique and fascinating aspects of mitochondrial gene expression, highlighting the mitochondrion as an evolutionary playground par excellence.

Identification of a Giardia krr1 Homolog Gene and the Secondarily Anucleolate Condition of Giaridia lamblia

Giaridia lamblia was long considered to be one of the most primitive eukaryotes and to lie close to the transition between prokaryotes and eukaryotes, but several supporting features, such as lack of mitochondrion and Golgi, have been challenged recently. It was also reported previously that G. lamblia lacked nucleolus, which is the site of pre-rRNA processing and ribosomal assembling in the other eukaryotic cells. Here, we report the identification of the yeast homolog gene, krr1, in the anucleolate eukaryote, G. lamblia. The krr1 gene, encoding one of the pre-rRNA processing proteins in yeast, is actively transcribed in G. lamblia. The deduced protein sequence of G. lamblia krr1 is highly similar to yeast KRR1p that contains a single-KH domain. Our database searches indicated that krr1 genes actually present in diverse eukaryotes and also seem to present in Archaea. However, only the eukaryotic homologs, including that of G. lamblia, have the single-KH domain, which contains the conserved motif KR(K)R. Fibrillarin, another important pre-rRNA processing protein has also been identified previously in G. lamblia. Moreover, our database search shows that nearly half of the other nucleolus-localized protein genes of eukaryotic cells also have their homologs in Giardia. Therefore, we suggest that a common mechanism of pre-RNA processing may operate in the anucleolate eukaryote G. lamblia and in the other eukaryotes and that like the case of "lack of mitochondrion," "lack of nucleolus" may not be a primitive feature, but a secondarily evolutionary condition of the parasite.

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.

A divergent ADP/ATP carrier in the hydrogenosomes of Trichomonas gallinae argues for an independent origin of these organelles

The evolution of mitochondrial ADP and ATP exchanging proteins (AACs) highlights a key event in the evolution of the eukaryotic cell, as ATP exporting carriers were indispensable in establishing the role of mitochondria as ATP-generating cellular organelles. Hydrogenosomes, i.e. ATP- and hydrogen-generating organelles of certain anaerobic unicellular eukaryotes, are believed to have evolved from the same ancestral endosymbiont that gave rise to present day mitochondria. Notably, the hydrogenosomes of the parasitic anaerobic flagellate Trichomonas seemed to be deficient in mitochondrial-type AACs. Instead, HMP 31, a different member of the mitochondrial carrier family (MCF) with a hitherto unknown function, is abundant in the hydrogenosomal membranes of Trichomonas vaginalis. Here we show that the homologous HMP 31 of closely related Trichomonas gallinae specifically transports ADP and ATP with high efficiency, as do genuine mitochondrial AACs. However, phylogenetic analysis and its resistance against bongkrekic acid (BKA, an efficient inhibitor of mitochondrial-type AACs) identify HMP 31 as a member of the mitochondrial carrier family that is distinct from all mitochondrial and hydrogenosomal AACs studied so far. Thus, our data support the hypothesis that the various hydrogenosomes evolved repeatedly and independently.

Endosomes, glycosomes, and glycosylphosphatidylinositol catabolism in Leishmania major

Glycosylphosphatidylinositols (GPIs) serve as membrane anchors of polysaccharides and proteins in the protozoan parasite Leishmania major. Free GPIs that are not attached to macromolecules are present in L. major as intermediates of protein-GPI and polysaccharide-GPI synthesis or as terminal glycolipids. The importance of the intracellular location of GPIs in vivo for functions of the glycolipids is not appreciated. To examine the roles of intracellular free GPI pools for attachment to polypeptide, a GPI-specific phospholipase C (GPI-PLCp) from Trypanosoma brucei was used to probe trafficking of GPI pools inside L. major. The locations of GPIs were determined, and their catabolism by GPI-PLCp was analyzed with respect to the intracellular location of the enzyme. GPIs accumulated on the endo-lysosomal system, where GPI-PLCp was also detected. A peptide motif [CS][CS]-x(0,2)-G-x(1)-C-x(2,3)-S-x(3)-L formed part of an endosome targeting signal for GPI-PLCp. Mutations of the endosome targeting motif caused GPI-PLCp to associate with glycosomes (peroxisomes). Endosomal GPI-PLCp caused a deficiency of protein-GPI in L. major, whereas glycosomal GPI-PLCp failed to produce the GPI deficiency. We surmise that (i) endo-lysosomal GPIs are important for biogenesis of GPI-anchored proteins in L. major; (ii) sequestration of GPI-PLCp to glycosomes protects free protein-GPIs from cleavage by the phospholipase. In T. brucei, protein-GPIs are concentrated at the endoplasmic reticulum, separated from GPI-PLCp. These observations support a model in which glycosome sequestration of a catabolic GPI-PLCp preserves free protein-GPIs in vivo.

Common evolutionary origin of mitochondrial and rickettsial respiratory chains

Comprehensive phylogenetic analysis of the subunits of respiratory chain was carried out using a variety of mitochondrial and bacterial sequences including those from all unfinished alpha-proteobacterial genomes known to date. Maximum likelihood, neighbor-joining, and maximum parsimony consensus trees, based on four proton-translocating complexes, placed mitochondria as a sister group to the order Rickettsiales of obligate endosymbiotic bacteria to the exclusion of free-living alpha-proteobacteria. Thus, phylogenetic relationship of most eukaryotic respiratory enzymes conforms to canonical pattern of mitochondrial ancestry, prior established in analyses of ribosomal RNAs, which are encoded by residual mitochondrial genomes. These data suggest that mitochondria may have derived from a reduced intracellular bacterium and that respiration may be the only evolutionary novelty brought into eukaryotes by mitochondrial endosymbiont.

Phylogenetic affinity of a Giardia lamblia cysteine desulfurase conforms to canonical pattern of mitochondrial ancestry

Among a few potential archezoan groups, only the Metamonada (diplomonads, retortamonads, and oxymonads) still retain the status of amitochondriate protists that diverged before the acquisition or retention of mitochondria. Indeed, finding that diplomonad genomes harbor a gene encoding a mitochondrial type chaperonin 60, the most compelling evidence for their secondarily amitochondriate nature, may be interpreted as an acquisition of this important general chaperone during some transient alpha-proteobacterial endosymbiosis. Recently published data on the cysteine desulfurase IscS demonstrated an alpha-proteobacterial origin of mitochondrial enzymes including a diplomonad Giardia lamblia homolog. An extended phylogenetic analysis of IscS is reported here that revealed a full canonical pattern of mitochondrial ancestry for the giardial enzyme. The above canonical pattern, a sister group relationship of mitochondria and rickettsiae exclusive of free-living alpha-proteobacteria, was robustly confirmed by a comprehensive analysis of Cob and Cox1 subunits of the respiratory chain encoded by resident mitochondrial genes. Given that Fe-S cluster assembly involving IscS represents an essential mitochondrial function, these data strongly suggest that diplomonads once harbored bona fide mitochondria.

Giardia lamblia RNA polymerase II: amanitin-resistant transcription

Giardia lamblia is an early branching eukaryote, and although distinctly eukaryotic in its cell and molecular biology, transcription and translation in G. lamblia demonstrate important differences from these processes in higher eukaryotes. The cyclic octapeptide amanitin is a relatively selective inhibitor of eukaryotic RNA polymerase II (RNAP II) and is commonly used to study RNAP II transcription. Therefore, we measured the sensitivity of G. lamblia RNAP II transcription to alpha-amanitin and found that unlike most other eukaryotes, RNAP II transcription in Giardia is resistant to 1 mg/ml amanitin. In contrast, 50 microg/ml amanitin inhibits 85% of RNAP III transcription activity using leucyl-tRNA as a template. To better understand transcription in G. lamblia, we identified 10 of the 12 known eukaryotic rpb subunits, including all 10 subunits that are required for viability in Saccharomyces cerevisiae. The amanitin motif (amanitin binding site) of Rpb1 from G. lamblia has amino acid substitutions at six highly conserved sites that have been associated with amanitin resistance in other organisms. These observations of amanitin resistance of Giardia RNA polymerase II support previous proposals of the mechanism of amanitin resistance in other organisms and provide a molecular framework for the development of novel drugs with selective activity against G. lamblia.

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.

On the origins of cells: a hypothesis for the evolutionary transitions from abiotic geochemistry to chemoautotrophic prokaryotes, and from prokaryotes to nucleated cells

All life is organized as cells. Physical compartmentation from the environment and self-organization of self-contained redox reactions are the most conserved attributes of living things, hence inorganic matter with such attributes would be life's most likely forebear. We propose that life evolved in structured iron monosulphide precipitates in a seepage site hydrothermal mound at a redox, pH and temperature gradient between sulphide-rich hydrothermal fluid and iron(II)-containing waters of the Hadean ocean floor. The naturally arising, three-dimensional compartmentation observed within fossilized seepage-site metal sulphide precipitates indicates that these inorganic compartments were the precursors of cell walls and membranes found in free-living prokaryotes. The known capability of FeS and NiS to catalyse the synthesis of the acetyl-methylsulphide from carbon monoxide and methylsulphide, constituents of hydrothermal fluid, indicates that pre-biotic syntheses occurred at the inner surfaces of these metal-sulphide-walled compartments, which furthermore restrained reacted products from diffusion into the ocean, providing sufficient concentrations of reactants to forge the transition from geochemistry to biochemistry. The chemistry of what is known as the RNA-world could have taken place within these naturally forming, catalyticwalled compartments to give rise to replicating systems. Sufficient concentrations of precursors to support replication would have been synthesized in situ geochemically and biogeochemically, with FeS (and NiS) centres playing the central catalytic role. The universal ancestor we infer was not a free-living cell, but rather was confined to the naturally chemiosmotic, FeS compartments within which the synthesis of its constituents occurred. The first free-living cells are suggested to have been eubacterial and archaebacterial chemoautotrophs that emerged more than 3.8 Gyr ago from their inorganic confines. We propose that the emergence of these prokaryotic lineages from inorganic confines occurred independently, facilitated by the independent origins of membrane-lipid biosynthesis: isoprenoid ether membranes in the archaebacterial and fatty acid ester membranes in the eubacterial lineage. The eukaryotes, all of which are ancestrally heterotrophs and possess eubacterial lipids, are suggested to have arisen ca. 2 Gyr ago through symbiosis involving an autotrophic archaebacterial host and a heterotrophic eubacterial symbiont, the common ancestor of mitochondria and hydrogenosomes. The attributes shared by all prokaryotes are viewed as inheritances from their confined universal ancestor. The attributes that distinguish eubacteria and archaebacteria, yet are uniform within the groups, are viewed as relics of their phase of differentiation after divergence from the non-free-living universal ancestor and before the origin of the free-living chemoautotrophic lifestyle. The attributes shared by eukaryotes with eubacteria and archaebacteria, respectively, are viewed as inheritances via symbiosis. The attributes unique to eukaryotes are viewed as inventions specific to their lineage. The origin of the eukaryotic endomembrane system and nuclear membrane are suggested to be the fortuitous result of the expression of genes for eubacterial membrane lipid synthesis by an archaebacterial genetic apparatus in a compartment that was not fully prepared to accommodate such compounds, resulting in vesicles of eubacterial lipids that accumulated in the cytosol around their site of synthesis. Under these premises, the most ancient divide in the living world is that between eubacteria and archaebacteria, yet the steepest evolutionary grade is that between prokaryotes and eukaryotes.

How big is the iceberg of which organellar genes in nuclear genomes are but the tip?

As more and more complete bacterial and archaeal genome sequences become available, the role of lateral gene transfer (LGT) in shaping them becomes more and more clear. Over the long term, it may be the dominant force, affecting most genes in most prokaryotes. We review the history of LGT, suggesting reasons why its prevalence and impact were so long dismissed. We discuss various methods purporting to measure the extent of LGT, and evidence for and against the notion that there is a core of never-exchanged genes shared by all genomes, from which we can deduce the "true" organismal tree. We also consider evidence for, and implications of, LGT between prokaryotes and phagocytic eukaryotes.

Autochthonous eukaryotic diversity in hydrothermal sediment and experimental microcolonizers at the Mid-Atlantic Ridge

The diversity and mode of life of microbial eukaryotes in hydrothermal systems is very poorly known. We carried out a molecular survey based on 18S ribosomal RNA genes of eukaryotes present in different hydrothermal niches at the Mid-Atlantic Ridge. These included metal-rich and rare-earth-element-rich hydrothermal sediments of the Rainbow site, fluid-seawater mixing regions, and colonization devices (microcolonizers) containing organic, iron-rich, and porous mineral substrates that were exposed for 15 days to a fluid source. We identified considerable phylogenetic diversity, both at kingdom level and within kinetoplastids and alveolates. None of our sequences affiliates to photosynthesizing lineages, suggesting that we are targeting only autochthonous deep-sea communities. Although sediment harbored most phylogenetic diversity, microcolonizers predominantly contained bodonids and ciliates, indicating that these protists pioneer the colonization process. Given the large variety of divergent lineages detected within the alveolates in deep-sea plankton, hydrothermal sediments, and vents, alveolates seem to dominate the deep ocean in terms of diversity. Compared with data from the Pacific Guaymas basin, some protist lineages seem ubiquitous in hydrothermal areas, whereas others, notably kinetoplastid lineages, very abundant and diverse in our samples, so far have been detected only in Atlantic systems. Unexpectedly, although alvinellid polychaetes are considered endemic of Pacific vents, we detected alvinellid-related sequences at the fluid-seawater interface and in microcolonizers. This finding can boost further studies on deep-sea vent animal biology and biogeography.

Intron gain and loss in the evolution of the conserved eukaryotic recombination machinery

Intron conservation, intron gain or loss and putative intron sliding events were determined for a set of three genes (SPO11, MRE11 and DMC1) involved in basic aspects of recombination in eukaryotes. These are ancient genes and present in nearly all of the major kingdoms. MRE11 is of bacterial origin and can be found in all kingdoms. DMC1 is a specialized homolog of the bacterial RecA protein, whereas the SPO11 gene is of archaebacterial origin. Only unique homologs of SPO11 are found in animals and fungi whereas three distantly related SPO11 copies are present in plant genomes. A comparison of the respective intron positions and phases of all genes was performed, demonstrating that a quarter of the intron positions were perfectly conserved over more than 1 000 000 000 years. Regarding the remaining three quarters of the introns we found insertions to be about three times more frequent than deletions. Aligning the introns of the three different SPO11 homologs of Arabidopsis thaliana we propose a conclusive model of their evolution. We postulate that at least one duplication event occurred shortly after the divergence of plants from animals and fungi and that a respective homolog has been retained in a protist group, the apicomplexa.

Evolution of collagens

The extracellular matrix is often defined as the substance that gives multicellular organisms (from plants to vertebrates) their structural integrity, and is intimately involved in their development. Although the general functions of extracellular matrices are comparable, their compositions are quite distinct. One of the specific components of metazoan extracellular matrices is collagen, which is present in organisms ranging from sponges to humans. By comparing data obtained in diploblastic, protostomic, and deuterostomic animals, we have attempted to trace the evolution of collagens and collagen-like proteins. Moreover, the collagen story is closely involved with the emergence and evolution of metazoa. The collagen triple helix is one of numerous modules that arose during the metazoan radiation which permit the formation of large multimodular proteins. One of the advantages of this module is its involvement in oligomerization, in which it acts as a structural organizer that is not only relatively resistant to proteases but also permits the creation of multivalent supramolecular networks.

Characterization and evolutionary analysis of a large polygalacturonase gene family in the oomycete plant pathogen Phytophthora cinnamomi

Polygalacturonases (PGs) are secreted by fungal pathogens during saprophytic and parasitic growth, and their degradation of pectin in the plant cell wall is believed to play a major role in tissue invasion and maceration. In this study, PG activity was demonstrated in culture filtrates of the oomycete plant pathogen, Phytophthora cinnamomi. A P. cinnamomi pg gene fragment amplified using degenerate primers based on conserved regions in fungal and plant PGs was used to isolate 17 complete P. cinnamomi pg genes and pseudogenes from a genomic library and partial sequence for another two genes. Gel blotting of genomic DNA indicated that there may be even more pg genes in the P. cinnamomi genome. P. cinnamomi pg gene sequences were expressed in PG-deficient yeast and found to confer PG activity, thereby confirming their functional identity. The predicted mature P. cinnamomi PGs fall into subgroups that exhibit large differences in the extent of N-glycosylation, isoelectric points, and N- and C-terminal structure. Evidence for birth-and-death and reticulate evolution in the P. cinnamomi pg gene family was obtained, and some codons for surface exposed residues in the P. cinnamomi PGs were shown to have been subject to diversifying selection. Contrary to accepted phylogenies for other proteins, phylogenetic analysis of the P. cinnamomi PGs revealed a closer relationship with PGs from true fungi than with those from plants.