Search

Genome-scale data sets result in an enhanced resolution of the phylogenetic inference by reducing stochastic errors. However, there is also an increase of systematic errors due to model violations, which can lead to erroneous phylogenies. Here, we explore the impact of systematic errors on the resolution of the eukaryotic phylogeny using a data set of 143 nuclear-encoded proteins from 37 species. The initial observation was that, despite the impressive amount of data, some branches had no significant statistical support. To demonstrate that this lack of resolution is due to a mutual annihilation of phylogenetic and nonphylogenetic signals, we created a series of data sets with slightly different taxon sampling. As expected, these data sets yielded strongly supported but mutually exclusive trees, thus confirming the presence of conflicting phylogenetic and nonphylogenetic signals in the original data set. To decide on the correct tree, we applied several methods expected to reduce the impact of some kinds of systematic error. Briefly, we show that (i) removing fast-evolving positions, (ii) recoding amino acids into functional categories, and (iii) using a site-heterogeneous mixture model (CAT) are three effective means of increasing the ratio of phylogenetic to nonphylogenetic signal. Finally, our results allow us to formulate guidelines for detecting and overcoming phylogenetic artefacts in genome-scale phylogenetic analyses. [ABSTRACT FROM AUTHOR]

Recent advances in resolving the tree of eukaryotes are converging on a model composed of a few large hypothetical ‘supergroups’, each comprising a diversity of primarily microbial eukaryotes (protists, or protozoa and algae). The process of resolving the tree involves the synthesis of many kinds of data, including single-gene trees, multigene analyses, and other kinds of molecular and structural characters. Here, we review the recent progress in assembling the tree of eukaryotes, describing the major evidence for each supergroup, and where gaps in our knowledge remain. We also consider other factors emerging from phylogenetic analyses and comparative genomics, in particular lateral gene transfer, and whether such factors confound our understanding of the eukaryotic tree. [Copyright &y& Elsevier]

Animal mtDNAs are typically small (≈16 kbp), circular-mapping molecules that encode 37 or fewer tightly packed genes. Here we investigate whether similarly compact mitochondrial genomes are also present in the closest unicellular relatives of animals, i.e., choanoflagellate and ichthyosporean protists. We find that the gene content and architecture of the mitochondrial genomes of the choanoflagellate Monosiga brevicollis, the ichthyosporean Amoebidium parasiticum, and Metazoa are radically different from one another. The circular-mapping choanoflagellate mtDNA with its long intergenic regions is four times as large and contains two times as many protein genes as do animal mtDNAs, whereas the ichthyosporean mitochondrial genome totals >200 kbp and consists of several hundred linear chromosomes that share elaborate terminal-specific sequence patterns. The highly peculiar organization of the ichthyosporean mtDNA raises questions about the mechanism of mitochondriai genome replication and chromosome segregation during cell division in this organism. Considering that the closest unicellular relatives of animals possess large, spacious, gene-rich mtDNAs, we posit that the distinct compaction characteristic of metazoan mitochondriai genomes occurred simultaneously with the emergence of a multicellular body plan in the animal lineage. [ABSTRACT FROM AUTHOR]

Fungi colonizing plants are gaining attention because of their ability to promote plant growth and suppress pathogens. While most studies focus on endosymbionts from grasses and legumes, the large and diverse group of ericaceous plants has been much neglected. We recently described one of the very few fungal endophytes promoting the growth of the Ericaceae Vaccinium macrocarpon (American cranberry), notably the Codinaeella isolate EC4. Here, we show that EC4 also suppresses fungal pathogens, which makes it a promising endophyte for sustainable cranberry cultivation. By dual-culture assays on agar plates, we tested the potential growth suppression (or biocontrol) of EC4 on other microbes, notably 12 pathogenic fungi and one oomycete reported to infect not only cranberry but also blueberry, strawberry, tomato plants, rose bushes and olive trees. Under greenhouse conditions, EC4 protects cranberry plantlets infected with one of the most notorious cranberry-plant pathogens, Diaporthe vaccinii , known to cause upright dieback and berry rot. The nuclear genome sequence of EC4 revealed a large arsenal of genes potentially involved in biocontrol. About ∼60 distinct clusters of genes are homologs of secondary metabolite gene clusters, some of which were shown in other fungi to synthesize nonribosomal peptides and polyketides, but in most cases, the exact compounds these clusters may produce are unknown. The EC4 genome also encodes numerous homologs of hydrolytic enzymes known to degrade fungal cell walls. About half of the nearly 250 distinct glucanases and chitinases are likely involved in biocontrol because they are predicted to be secreted outside the cell. Transcriptome analysis shows that the expression of about a quarter of the predicted secondary-metabolite gene clusters and glucan and chitin-degrading genes of EC4 is stimulated when it is co-cultured with D. vaccinii . Some of the differentially expressed EC4 genes are alternatively spliced exclusively in the presence of the pathogen, altering the proteins' domain content and subcellular localization signal, thus adding a second level of proteome adaptation in response to habitat competition. To our knowledge, this is the first report of Diaporthe -induced alternative splicing of biocontrol genes., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest. The authors declared that they were an editorial board member of Frontiers, at the time of submission. This had no impact on the peer review process and the final decision., (Copyright © 2024 Thimmappa, Salhi, Forget, Sarrasin, Bustamante Villalobos, Henrissat, Lang and Burger.)

Compared to nuclear genomes, mitochondrial genomes (mitogenomes) are small and usually code for only a few dozen genes. Still, identifying genes and their structure can be challenging and time-consuming. Even automated tools for mitochondrial genome annotation often require manual analysis and curation by skilled experts. The most difficult steps are (i) the structural modelling of intron-containing genes; (ii) the identification and delineation of Group I and II introns; and (iii) the identification of moderately conserved, non-coding RNA (ncRNA) genes specifying 5S rRNAs, tmRNAs and RNase P RNAs. Additional challenges arise through genetic code evolution which can redefine the translational identity of both start and stop codons, thus obscuring protein-coding genes. Further, RNA editing can render gene identification difficult, if not impossible, without additional RNA sequence data. Current automated mito- and plastid-genome annotators are limited as they are typically tailored to specific eukaryotic groups. The MFannot annotator we developed is unique in its applicability to a broad taxonomic scope, its accuracy in gene model inference, and its capabilities in intron identification and classification. The pipeline leverages curated profile Hidden Markov Models (HMMs), covariance (CMs) and ERPIN models to better capture evolutionarily conserved signatures in the primary sequence (HMMs and CMs) as well as secondary structure (CMs and ERPIN). Here we formally describe MFannot, which has been available as a web-accessible service (https://megasun.bch.umontreal.ca/apps/mfannot/) to the research community for nearly 16 years. Further, we report its performance on particularly intron-rich mitogenomes and describe ongoing and future developments., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2023 Lang, Beck, Prince, Sarrasin, Rioux and Burger.)

The unicellular amoeba Acanthamoeba castellanii is ubiquitous in aquatic environments, where it preys on bacteria. The organism also hosts bacterial endosymbionts, some of which are parasitic, including human pathogens such as Chlamydia and Legionella spp. Here we report complete, high-quality genome sequences for two extensively studied A. castellanii strains, Neff and C3. Combining long- and short-read data with Hi-C, we generated near chromosome-level assemblies for both strains with 90% of the genome contained in 29 scaffolds for the Neff strain and 31 for the C3 strain. Comparative genomics revealed strain-specific functional enrichment, most notably genes related to signal transduction in the C3 strain and to viral replication in Neff. Furthermore, we characterized the spatial organization of the A. castellanii genome and showed that it is reorganized during infection by Legionella pneumophila Infection-dependent chromatin loops were found to be enriched in genes for signal transduction and phosphorylation processes. In genomic regions where chromatin organization changed during Legionella infection, we found functional enrichment for genes associated with metabolism, organelle assembly, and cytoskeleton organization. Given Legionella infection is known to alter its host's cell cycle, to exploit the host's organelles, and to modulate the host's metabolism in its favor, these changes in chromatin organization may partly be related to mechanisms of host control during Legionella infection., (© 2022 Matthey-Doret et al.; Published by Cold Spring Harbor Laboratory Press.)

Most of the described species in kingdom Fungi are contained in two phyla, the Ascomycota and the Basidiomycota (subkingdom Dikarya). As a result, our understanding of the biology of the kingdom is heavily influenced by traits observed in Dikarya, such as aerial spore dispersal and life cycles dominated by mitosis of haploid nuclei. We now appreciate that Fungi comprises numerous phylum-level lineages in addition to those of Dikarya, but the phylogeny and genetic characteristics of most of these lineages are poorly understood due to limited genome sampling. Here, we addressed major evolutionary trends in the non-Dikarya fungi by phylogenomic analysis of 69 newly generated draft genome sequences of the zoosporic (flagellated) lineages of true fungi. Our phylogeny indicated five lineages of zoosporic fungi and placed Blastocladiomycota, which has an alternation of haploid and diploid generations, as branching closer to the Dikarya than to the Chytridiomyceta. Our estimates of heterozygosity based on genome sequence data indicate that the zoosporic lineages plus the Zoopagomycota are frequently characterized by diploid-dominant life cycles. We mapped additional traits, such as ancestral cell-cycle regulators, cell-membrane- and cell-wall-associated genes, and the use of the amino acid selenocysteine on the phylogeny and found that these ancestral traits that are shared with Metazoa have been subject to extensive parallel loss across zoosporic lineages. Together, our results indicate a gradual transition in the genetics and cell biology of fungi from their ancestor and caution against assuming that traits measured in Dikarya are typical of other fungal lineages.

Fungi survive in diverse ecological niches by secreting proteins and other molecules into the environment to acquire food and interact with various biotic and abiotic stressors. Fungal secretome content is, therefore, believed to be tightly linked to fungal ecologies. We sampled 132 genomes from the early-diverging terrestrial fungal lineage zygomycetes (Mucoromycota and Zoopagomycota) and characterized their secretome composition. Our analyses revealed that phylogeny played an important role in shaping the secretome composition of zygomycete fungi with trophic mode contributing a smaller amount. Reconstruction of the evolution of secreted digestive enzymes revealed lineage-specific expansions, indicating that Mucoromycota and Zoopagomycota followed different trajectories early in their evolutionary history. We identified the presence of multiple pathogenicity-related proteins in the lineages known as saprotrophs, suggesting that either the ecologies of these fungi are incompletely known, and/or that these pathogenicity-related proteins have important functions associated with saprotrophic ecologies, both of which invite further investigation., Competing Interests: The authors declare no conflict of interest., (© 2022 The Author(s).)

Mitochondrial genomes-in particular those of fungi-often encode genes with a large number of Group I and Group II introns that are conserved at both the sequence and the RNA structure level. They provide a rich resource for the investigation of intron and gene structure, self- and protein-guided splicing mechanisms, and intron evolution. Yet, the degree of sequence conservation of introns is limited, and the primary sequence differs considerably among the distinct intron sub-groups. It makes intron identification, classification, structural modeling, and the inference of gene models a most challenging and error-prone task-frequently passed on to an "expert" for manual intervention. To reduce the need for manual curation of intron structures and mitochondrial gene models, computational methods using ERPIN sequence profiles were initially developed in 2007. Here we present a refinement of search models and alignments using the now abundant publicly available fungal mtDNA sequences. In addition, we have tested in how far members of the originally proposed sub-groups are clearly distinguished and validated by our computational approach. We confirm clearly distinct mitochondrial Group I sub-groups IA1, IA3, IB3, IC1, IC2, and ID. Yet, IB1, IB2, and IB4 ERPIN models are overlapping substantially in predictions, and are therefore combined and reported as IB. We have further explored the conversion of our ERPIN profiles into covariance models (CM). Current limitations and prospects of the CM approach will be discussed., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2022 Prince, Munoz, Filion-Bienvenue, Rioux, Sarrasin and Lang.)

Fungal mitochondrial genomes encode genes involved in crucial cellular processes, such as oxidative phosphorylation and mitochondrial translation, and the molecule has been used as a molecular marker for population genetics studies. Coccidioides immitis and C. posadasii are endemic fungal pathogens that cause coccidioidomycosis in arid regions across both American continents. To date, approximately 150 Coccidioides isolates have been sequenced to infer patterns of variation in nuclear genomes. However, less attention has been given to the mitochondrial genomes of Coccidioides. In this report, we describe the assembly and annotation of mitochondrial reference genomes for two representative strains of C. posadasii and C. immitis, as well as assess population variation among 77 selected genomes. The sizes of the circular-mapping molecules are 68.2 Kb in C. immitis and 75.1 Kb in C. posadasii. We identify 14 mitochondrial protein-coding genes common to most fungal mitochondria, which are largely syntenic across different populations and species of Coccidioides. Both Coccidioides species are characterized by a large number of group I and II introns, harboring twice the number of elements as compared to closely related Onygenales. The introns contain complete or truncated ORFs with high similarity to homing endonucleases of the LAGLIDADG and GIY-YIG families. Phylogenetic comparisons of mitochondrial and nuclear genomes show extensive phylogenetic discordance suggesting that the evolution of the two types of genetic material is not identical. This work represents the first assessment of mitochondrial genomes among isolates of both species of Coccidioides, and provides a foundation for future functional work., (© The Author(s) 2021. Published by Oxford University Press on behalf of Genetics Society of America.)

Here, we experimentally expand understanding of the reactions and enzymes involved in Acidithiobacillus thiooxidans ATCC 19377 S 0 and S 2 ⁢ O 3 2 - metabolism by developing models that integrate gene expression analyzed by RNA-Seq, solution sulfur speciation, electron microscopy and spectroscopy. The A. thiooxidans S 2 ⁢ O 3 2 - metabolism model involves the conversion of S 2 ⁢ O 3 2 - to SO 4 2 - , S 0 and S 4 ⁢ O 6 2 - , mediated by the sulfur oxidase complex (Sox), tetrathionate hydrolase (TetH), sulfide quinone reductase (Sqr), and heterodisulfate reductase (Hdr) proteins. These same proteins, with the addition of rhodanese (Rhd), were identified to convert S 0 to SO 3 2 - , S 2 ⁢ O 3 2 - and polythionates in the A. thiooxidans S 0 metabolism model. Our combined results shed light onto the important role specifically of TetH in S 2 ⁢ O 3 2 - metabolism. Also, we show that activity of Hdr proteins rather than Sdo are likely associated with S 0 oxidation. Finally, our data suggest that formation of intracellular S 2 ⁢ O 3 2 - is a critical step in S 0 metabolism, and that recycling of internally generated SO 3 2 - occurs, through comproportionating reactions that result in S 2 ⁢ O 3 2 - . Electron microscopy and spectroscopy confirmed intracellular production and storage of S 0 during growth on both S 0 and S 2 ⁢ O 3 2 - substrates., (Copyright © 2020 Camacho, Frazao, Fouillen, Nanci, Lang, Apte, Baron and Warren.)

Plants and fungi use light and other signals to regulate development, growth, and metabolism. The fruiting bodies of the fungus Phycomyces blakesleeanus are single cells that react to environmental cues, including light, but the mechanisms are largely unknown [1]. The related fungus Mucor circinelloides is an opportunistic human pathogen that changes its mode of growth upon receipt of signals from the environment to facilitate pathogenesis [2]. Understanding how these organisms respond to environmental cues should provide insights into the mechanisms of sensory perception and signal transduction by a single eukaryotic cell, and their role in pathogenesis. We sequenced the genomes of P. blakesleeanus and M. circinelloides and show that they have been shaped by an extensive genome duplication or, most likely, a whole-genome duplication (WGD), which is rarely observed in fungi [3-6]. We show that the genome duplication has expanded gene families, including those involved in signal transduction, and that duplicated genes have specialized, as evidenced by differences in their regulation by light. The transcriptional response to light varies with the developmental stage and is still observed in a photoreceptor mutant of P. blakesleeanus. A phototropic mutant of P. blakesleeanus with a heterozygous mutation in the photoreceptor gene madA demonstrates that photosensor dosage is important for the magnitude of signal transduction. We conclude that the genome duplication provided the means to improve signal transduction for enhanced perception of environmental signals. Our results will help to understand the role of genome dynamics in the evolution of sensory perception in eukaryotes., (Copyright © 2016 Elsevier Ltd. All rights reserved.)

The large phylogenetic distance separating eukaryotic genes and their archaeal orthologs has prevented identification of the position of the eukaryotic root in phylogenomic studies. Recently, an innovative approach has been proposed to circumvent this issue: the use as phylogenetic markers of proteins that have been transferred from bacterial donor sources to eukaryotes, after their emergence from Archaea. Using this approach, two recent independent studies have built phylogenomic datasets based on bacterial sequences, leading to different predictions of the eukaryotic root. Taking advantage of additional genome sequences from the jakobid Andalucia godoyi and the two known malawimonad species (Malawimonas jakobiformis and Malawimonas californiana), we reanalyzed these two phylogenomic datasets. We show that both datasets pinpoint the same phylogenetic position of the eukaryotic root that is between "Unikonta" and "Bikonta," with malawimonad and collodictyonid lineages on the Unikonta side of the root. Our results firmly indicate that (i) the supergroup Excavata is not monophyletic and (ii) the last common ancestor of eukaryotes was a biflagellate organism. Based on our results, we propose to rename the two major eukaryotic groups Unikonta and Bikonta as Opimoda and Diphoda, respectively.

Programmed translational bypassing is a process whereby ribosomes "ignore" a substantial interval of mRNA sequence. Although discovered 25 y ago, the only experimentally confirmed example of this puzzling phenomenon is expression of the bacteriophage T4 gene 60. Bypassing requires translational blockage at a "takeoff codon" immediately upstream of a stop codon followed by a hairpin, which causes peptidyl-tRNA dissociation and reassociation with a matching "landing triplet" 50 nt downstream, where translation resumes. Here, we report 81 translational bypassing elements (byps) in mitochondria of the yeast Magnusiomyces capitatus and demonstrate in three cases, by transcript analysis and proteomics, that byps are retained in mitochondrial mRNAs but not translated. Although mitochondrial byps resemble the bypass sequence in the T4 gene 60, they utilize unused codons instead of stops for translational blockage and have relaxed matching rules for takeoff/landing sites. We detected byp-like sequences also in mtDNAs of several Saccharomycetales, indicating that byps are mobile genetic elements. These byp-like sequences lack bypassing activity and are tolerated when inserted in-frame in variable protein regions. We hypothesize that byp-like elements have the potential to contribute to evolutionary diversification of proteins by adding new domains that allow exploration of new structures and functions.

The discovery of diverse codon reassignment events has demonstrated that the canonical genetic code is not universal. Studying coding reassignment at the molecular level is critical for understanding genetic code evolution, and provides clues to genetic code manipulation in synthetic biology. Here we report a novel reassignment event in the mitochondria of Ashbya (Eremothecium) gossypii, a filamentous-growing plant pathogen related to yeast (Saccharomycetaceae). Bioinformatics studies of conserved positions in mitochondrial DNA-encoded proteins suggest that CUU and CUA codons correspond to alanine in A. gossypii, instead of leucine in the standard code or threonine in yeast mitochondria. Reassignment of CUA to Ala was confirmed at the protein level by mass spectrometry. We further demonstrate that a predicted tRNA(Ala)UAG is transcribed and accurately processed in vivo, and is responsible for Ala reassignment. Enzymatic studies reveal that tRNA(Ala)UAG is efficiently recognized by A. gossypii mitochondrial alanyl-tRNA synthetase (AgAlaRS). AlaRS typically recognizes the G3:U70 base pair of tRNA(Ala); a G3A change in Ashbya tRNA(Ala)UAG abolishes its recognition by AgAlaRS. Conversely, an A3G mutation in Saccharomyces cerevisiae tRNA(Thr)UAG confers tRNA recognition by AgAlaRS. Our work highlights the dynamic feature of natural genetic codes in mitochondria, and the relative simplicity by which tRNA identity may be switched.

Filopodia are fine actin-based cellular projections used for both environmental sensing and cell motility, and they are essential organelles for metazoan cells. In this study, we reconstruct the origin of metazoan filopodia and microvilli. We first report on the evolutionary assembly of the filopodial molecular toolkit and show that homologs of many metazoan filopodial components, including fascin and myosin X, were already present in the unicellular or colonial progenitors of metazoans. Furthermore, we find that the actin crosslinking protein fascin localizes to filopodia-like structures and microvilli in the choanoflagellate Salpingoeca rosetta. In addition, homologs of filopodial genes in the holozoan Capsaspora owczarzaki are upregulated in filopodia-bearing cells relative to those that lack them. Therefore, our findings suggest that proteins essential for metazoan filopodia and microvilli are functionally conserved in unicellular and colonial holozoans and that the last common ancestor of metazoans bore a complex and specific filopodial machinery.

Most eukaryotic lineages belong to one of a few major groups. However, several protistan lineages have not yet been robustly placed in any of these groups. Both the breviates and apusomonads are two such lineages that appear to be related to the Amoebozoa and Opisthokonta (i.e. the 'unikonts' or Amorphea); however, their precise phylogenetic positions remain unclear. Here, we describe a novel microaerophilic breviate, Pygsuia biforma gen. nov. sp. nov., isolated from a hypoxic estuarine sediment. Ultrastructurally, this species resembles the breviate genera Breviata and Subulatomonas but has two cell morphologies, adherent and swimming. Phylogenetic analyses of the small sub-unit rRNA gene show that Pygsuia is the sister to the other breviates. We constructed a 159-protein supermatrix, including orthologues identified in RNA-seq data from Pygsuia. Phylogenomic analyses of this dataset show that breviates, apusomonads and Opisthokonta form a strongly supported major eukaryotic grouping we name the Obazoa. Although some phylogenetic methods disagree, the balance of evidence suggests that the breviate lineage forms the deepest branch within Obazoa. We also found transcripts encoding a nearly complete integrin adhesome from Pygsuia, indicating that this protein complex involved in metazoan multicellularity may have evolved earlier in eukaryote evolution than previously thought.

Sponges (phylum Porifera) are a large and ancient group of morphologically simple but ecologically important aquatic animals. Although their body plan and lifestyle are relatively uniform, sponges show extensive molecular and genetic diversity. In particular, mitochondrial genomes from three of the four previously studied classes of Porifera (Demospongiae, Hexactinellida, and Homoscleromorpha) have distinct gene contents, genome organizations, and evolutionary rates. Here, we report the mitochondrial genome of Clathrina clathrus (Calcinea, Clathrinidae), a representative of the fourth poriferan class, the Calcarea, which proves to be the most unusual. Clathrina clathrus mitochondrial DNA (mtDNA) consists of six linear chromosomes 7.6-9.4 kb in size and encodes at least 37 genes: 13 protein codings, 2 ribosomal RNAs (rRNAs), and 24 transfer RNAs (tRNAs). Protein genes include atp9, which has now been found in all major sponge lineages, but no atp8. Our analyses further reveal the presence of a novel genetic code that involves unique reassignments of the UAG codons from termination to tyrosine and of the CGN codons from arginine to glycine. Clathrina clathrus mitochondrial rRNAs are encoded in three (srRNA) and ≥6 (lrRNA) fragments distributed out of order and on several chromosomes. The encoded tRNAs contain multiple mismatches in the aminoacyl acceptor stems that are repaired posttranscriptionally by 3'-end RNA editing. Although our analysis does not resolve the phylogenetic position of calcareous sponges, likely due to their high rates of mitochondrial sequence evolution, it confirms mtDNA as a promising marker for population studies in this group. The combination of unusual mitochondrial features in C. clathrus redefines the extremes of mtDNA evolution in animals and further argues against the idea of a "typical animal mtDNA."

Gigaspora rosea is a member of the arbuscular mycorrhizal fungi (AMF; Glomeromycota) and a distant relative of Glomus species that are beneficial to plant growth. To allow for a better understanding of Glomeromycota, we have sequenced the mitochondrial DNA of G. rosea. A comparison with Glomus mitochondrial genomes reveals that Glomeromycota undergo insertion and loss of mitochondrial plasmid-related sequences and exhibit considerable variation in introns. The gene order between the two species is almost completely reshuffled. Furthermore, Gigaspora has fragmented cox1 and rns genes, and an unorthodox initiator tRNA that is tailored to decoding frequent UUG initiation codons. For the fragmented cox1 gene, we provide evidence that its RNA is joined via group I-mediated trans-splicing, whereas rns RNA remains in pieces. According to our model, the two cox1 precursor RNA pieces are brought together by flanking cox1 exon sequences that form a group I intron structure, potentially in conjunction with the nad5 intron 3 sequence. Finally, we present analyses that address the controversial phylogenetic association of Glomeromycota within fungi. According to our results, Glomeromycota are not a separate group of paraphyletic zygomycetes but branch together with Mortierellales, potentially also Harpellales.

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.

Initial steps in the synthesis of functional tRNAs require 5'- and 3'-processing of precursor tRNAs (pre-tRNAs), which in yeast mitochondria are achieved by two endonucleases, RNase P and RNase Z. In this study, using a combination of detergent-free Blue Native Gel Electrophoresis, proteomics and in vitro testing of pre-tRNA maturation, we reveal the physical association of these plus other mitochondrial activities in a large, stable complex of 136 proteins. It contains a total of seven proteins involved in RNA processing including RNase P and RNase Z, five out of six subunits of the mitochondrial RNA degradosome, components of the fatty acid synthesis pathway, translation, metabolism and protein folding. At the RNA level, there are the small and large rRNA subunits and RNase P RNA. Surprisingly, this complex is absent in an oar1Δ deletion mutant of the type II fatty acid synthesis pathway, supporting a recently published functional link between pre-tRNA processing and the FAS II pathway--apparently by integration into a large complex as we demonstrate here. Finally, the question of mt-RNase P localization within mitochondria was investigated, by GFP-tracing of a known protein subunit (Rpm2p). We find that about equal fractions of RNase P are soluble versus membrane-attached.

Many of the eukaryotic phylogenomic analyses published to date were based on alignments of hundreds to thousands of genes. Frequently, in such analyses, the most realistic evolutionary models currently available are often used to minimize the impact of systematic error. However, controversy remains over whether or not idiosyncratic gene family dynamics (i.e., gene duplications and losses) and incorrect orthology assignments are always appropriately taken into account. In this paper, we present an innovative strategy for overcoming orthology assignment problems. Rather than identifying and eliminating genes with paralogy problems, we have constructed a data set comprised exclusively of conserved single-copy protein domains that, unlike most of the commonly used phylogenomic data sets, should be less confounded by orthology miss-assignments. To evaluate the power of this approach, we performed maximum likelihood and Bayesian analyses to infer the evolutionary relationships within the opisthokonts (which includes Metazoa, Fungi, and related unicellular lineages). We used this approach to test 1) whether Filasterea and Ichthyosporea form a clade, 2) the interrelationships of early-branching metazoans, and 3) the relationships among early-branching fungi. We also assessed the impact of some methods that are known to minimize systematic error, including reducing the distance between the outgroup and ingroup taxa or using the CAT evolutionary model. Overall, our analyses support the Filozoa hypothesis in which Ichthyosporea are the first holozoan lineage to emerge followed by Filasterea, Choanoflagellata, and Metazoa. Blastocladiomycota appears as a lineage separate from Chytridiomycota, although this result is not strongly supported. These results represent independent tests of previous phylogenetic hypotheses, highlighting the importance of sophisticated approaches for orthology assignment in phylogenomic analyses.

The standard genetic code is used by most living organisms, yet deviations have been observed in many genomes, suggesting that the genetic code has been evolving. In certain yeast mitochondria, CUN codons are reassigned from leucine to threonine, which requires an unusual tRNA(Thr) with an enlarged 8-nt anticodon loop ( ). To trace its evolutionary origin we performed a comprehensive phylogenetic analysis which revealed that evolved from yeast mitochondrial tRNA(His). To understand this tRNA identity change, we performed mutational and biochemical experiments. We show that Saccharomyces cerevisiae mitochondrial threonyl-tRNA synthetase (MST1) could attach threonine to both and the regular , but not to the wild-type tRNA(His). A loss of the first nucleotide (G(-1)) in tRNA(His) converts it to a substrate for MST1 with a K(m) value (0.7 μM) comparable to that of (0.3 μM), and addition of G(-1) to allows efficient histidylation by histidyl-tRNA synthetase. We also show that MST1 from Candida albicans, a yeast in which CUN codons remain assigned to leucine, could not threonylate , suggesting that MST1 has coevolved with . Our work provides the first clear example of a recent recoding event caused by alloacceptor tRNA gene recruitment.

Mitochondrial genome diversity in closely related species provides an excellent platform for investigation of chromosome architecture and its evolution by means of comparative genomics. In this study, we determined the complete mitochondrial DNA sequences of eight Candida species and analyzed their molecular architectures. Our survey revealed a puzzling variability of genome architecture, including circular- and linear-mapping and multipartite linear forms. We propose that the arrangement of large inverted repeats identified in these genomes plays a crucial role in alterations of their molecular architectures. In specific arrangements, the inverted repeats appear to function as resolution elements, allowing genome conversion among different topologies, eventually leading to genome fragmentation into multiple linear DNA molecules. We suggest that molecular transactions generating linear mitochondrial DNA molecules with defined telomeric structures may parallel the evolutionary emergence of linear chromosomes and multipartite genomes in general and may provide clues for the origin of telomeres and pathways implicated in their maintenance.

According to the chromalveolate hypothesis (Cavalier-Smith T. 1999. Principles of protein and lipid targeting in secondary symbiogenesis: euglenoid, dinoflagellate, and sporozoan plastid origins and the eukaryote family tree. J Eukaryot Microbiol 46:347-366), the four eukaryotic groups with chlorophyll c-containing plastids originate from a single photosynthetic ancestor, which acquired its plastids by secondary endosymbiosis with a red alga. So far, molecular phylogenies have failed to either support or disprove this view. Here, we devise a phylogenomic falsification of the chromalveolate hypothesis that estimates signal strength across the three genomic compartments: If the four chlorophyll c-containing lineages indeed derive from a single photosynthetic ancestor, then similar amounts of plastid, mitochondrial, and nuclear sequences should allow to recover their monophyly. Our results refute this prediction, with statistical support levels too different to be explained by evolutionary rate variation, phylogenetic artifacts, or endosymbiotic gene transfer. Therefore, we reject the chromalveolate hypothesis as falsified in favor of more complex evolutionary scenarios involving multiple higher order eukaryote-eukaryote endosymbioses.

Acyl-CoA dehydrogenases (ACADs), which are key enzymes in fatty acid and amino acid catabolism, form a large, pan-taxonomic protein family with at least 13 distinct subfamilies. Yet most reported ACAD members have no subfamily assigned, and little is known about the taxonomic distribution and evolution of the subfamilies. In completely sequenced genomes from approximately 210 species (eukaryotes, bacteria and archaea), we detect ACAD subfamilies by rigorous ortholog identification combining sequence similarity search with phylogeny. We then construct taxonomic subfamily-distribution profiles and build phylogenetic trees with orthologous proteins. Subfamily profiles provide unparalleled insight into the organisms' energy sources based on genome sequence alone and further predict enzyme substrate specificity, thus generating explicit working hypotheses for targeted biochemical experimentation. Eukaryotic ACAD subfamilies are traditionally considered as mitochondrial proteins, but we found evidence that in fungi one subfamily is located in peroxisomes and participates in a distinct beta-oxidation pathway. Finally, we discern horizontal transfer, duplication, loss and secondary acquisition of ACAD genes during evolution of this family. Through these unorthodox expansion strategies, the ACAD family is proficient in utilizing a large range of fatty acids and amino acids-strategies that could have shaped the evolutionary history of many other ancient protein families.

Sequencing of cDNA libraries is an efficient and inexpensive approach to analyze the protein-coding portion of a genome. It is frequently used for surveying the genomes of poorly studied eukaryotes, and is particularly useful for species that are not easily amenable to genome sequencing, because they are nonaxenic and/or difficult to cultivate. In this chapter, we describe protocols that have been applied successfully to construct and normalize a variety of cDNA libraries from many different species of free-living protists and fungi, and that require only small quantities of cell material.

The organelle genome database GOBASE, now in its 21st release (June 2008), contains all published mitochondrion-encoded sequences (approximately 913,000) and chloroplast-encoded sequences (approximately 250,000) from a wide range of eukaryotic taxa. For all sequences, information on related genes, exons, introns, gene products and taxonomy is available, as well as selected genome maps and RNA secondary structures. Recent major enhancements to database functionality include: (i) addition of an interface for RNA editing data, with substitutions, insertions and deletions displayed using multiple alignments; (ii) addition of medically relevant information, such as haplotypes, SNPs and associated disease states, to human mitochondrial sequence data; (iii) addition of fully reannotated genome sequences for Escherichia coli and Nostoc sp., for reference and comparison; and (iv) a number of interface enhancements, such as the availability of both genomic and gene-coding sequence downloads, and a more sophisticated literature reference search functionality with links to PubMed where available. Future projects include the transfer of GOBASE features to NCBI/GenBank, allowing long-term preservation of accumulated expert information. The GOBASE database can be found at http://gobase.bcm.umontreal.ca/. Queries about custom and large-scale data retrievals should be addressed to gobase@bch.umontreal.ca.

Several morphologically dissimilar ascomycete fungi including Schizosaccharomyces, Taphrina, Saitoella, Pneumocystis, and Neolecta have been grouped into the taxon Taphrinomycotina (Archiascomycota or Archiascomycotina), originally based on rRNA phylogeny. These analyses lack statistically significant support for the monophyly of this grouping, and although confirmed by more recent multigene analyses, this topology is contradicted by mitochondrial phylogenies. To resolve this inconsistency, we have assembled phylogenomic mitochondrial and nuclear data sets from four distantly related taphrinomycotina taxa: Schizosaccharomyces pombe, Pneumocystis carinii, Saitoella complicata, and Taphrina deformans. Our phylogenomic analyses based on nuclear data (113 proteins) conclusively support the monophyly of Taphrinomycotina, diverging as a sister group to Saccharomycotina + Pezizomycotina. However, despite the improved taxon sampling, Taphrinomycotina continue to be paraphyletic with the mitochondrial data set (13 proteins): Schizosaccharomyces species associate with budding yeasts (Saccharomycotina) and the other Taphrinomycotina group as a sister group to Saccharomycotina + Pezizomycotina. Yet, as Schizosaccharomyces and Saccharomycotina species are fast evolving, the mitochondrial phylogeny may be influenced by a long-branch attraction (LBA) artifact. After removal of fast-evolving sequence positions from the mitochondrial data set, we recover the monophyly of Taphrinomycotina. Our combined results suggest that Taphrinomycotina is a legitimate taxon, that this group of species diverges as a sister group to Saccharomycotina + Pezizomycotina, and that phylogenetic positioning of yeasts and fission yeasts with mitochondrial data is plagued by a strong LBA artifact.

This chapter outlines key considerations for constructing and implementing an EST database. Instead of showing the technological details step by step, emphasis is put on the design of an EST database suited to the specific needs of EST projects and how to choose the most suitable tools. Using TBestDB as an example, we illustrate the essential factors to be considered for database construction and the steps for data population and annotation. This process employs technologies such as PostgreSQL, Perl, and PHP to build the database and interface, and tools such as AutoFACT for data processing and annotation. We discuss these in comparison to other available technologies and tools, and explain the reasons for our choices.

The evolution of multicellular animals (Metazoa) from their unicellular ancestors was a key transition that was accompanied by the emergence and diversification of gene families associated with multicellularity. To clarify the timing and order of specific events in this transition, we conducted expressed sequence tag surveys on 4 putative protistan relatives of Metazoa including the choanoflagellate Monosiga ovata, the ichthyosporeans Sphaeroforma arctica and Amoebidium parasiticum, and the amoeba Capsaspora owczarzaki, and 2 members of Amoebozoa, Acanthamoeba castellanii and Mastigamoeba balamuthi. We find that homologs of genes involved in metazoan multicellularity exist in several of these unicellular organisms, including 1 encoding a membrane-associated guanylate kinase with an inverted arrangement of protein-protein interaction domains (MAGI) in Capsaspora. In Metazoa, MAGI regulates tight junctions involved in cell-cell communication. By phylogenomic analyses of genes encoded in nuclear and mitochondrial genomes, we show that the choanoflagellates are the closest relatives of the Metazoa, followed by the Capsaspora and Ichthyosporea lineages, although the branching order between the latter 2 groups remains unclear. Understanding the function of "metazoan-specific" proteins we have identified in these protists will clarify the evolutionary steps that led to the emergence of the Metazoa.

Traditionally, phylogenetic analyses over many genes combine data into a contiguous block. Under this concatenated model, all genes are assumed to evolve at the same rate. However, it is clear that genes evolve at very different rates and that accounting for this rate heterogeneity is important if we are to accurately infer phylogenies from heterogeneous multigene data sets. There remain open questions regarding how best to incorporate gene rate parameters into phylogenetic models and which properties of real data correlate with improved fit over the concatenated model. In this study, two methods of accounting for gene rate heterogeneity are compared: the n-parameter method, which allows for each of the n gene partitions to have a gene rate parameter, and the alpha-parameter method, which fits a distribution to the gene rates. Results demonstrate that the n-parameter method is both computationally faster and in general provides a better fit over the concatenated model than the alpha-parameter method. Furthermore, improved model fit over the concatenated model is highly correlated with the presence of a gene with a slow relative rate of evolution.

In phylogenetic analyses with combined multigene or multiprotein data sets, accounting for differing evolutionary dynamics at different loci is essential for accurate tree prediction. Existing maximum likelihood (ML) and Bayesian approaches are computationally intensive. We present an alternative approach that is orders of magnitude faster. The method, Distance Rates (DistR), estimates rates based upon distances derived from gene/protein sequence data. Simulation studies indicate that this technique is accurate compared with other methods and robust to missing sequence data. The DistR method was applied to a fungal mitochondrial data set, and the rate estimates compared well to those obtained using existing ML and Bayesian approaches. Inclusion of the protein rates estimated from the DistR method into the ML calculation of trees as a branch length multiplier resulted in a significantly improved fit as measured by the Akaike Information Criterion (AIC). Furthermore, bootstrap support for the ML topology was significantly greater when protein rates were used, and some evident errors in the concatenated ML tree topology (i.e., without protein rates) were corrected. [Bayesian credible intervals; DistR method; multigene phylogeny; PHYML; rate heterogeneity.].

Phylogenetic relationships among the metazoan phyla are the subject of an ongoing controversy. Analysis of mitochondrial gene arrangements is a powerful tool to investigate these relationships; however, its previous application outside of individual animal phyla has been hampered by the lack of informative out-group data. To address this shortcoming, we determined complete mitochondrial DNA sequences for the demosponges Geodia neptuni and Tethya actinia, two representatives of the most basal animal phylum, the Porifera. With sponges as an outgroup, we investigated phylogenetic relationships of nine bilaterian phyla using both breakpoint analysis of global mitochondrial gene arrangements and maximum parsimony analysis of mitochondrial gene adjacencies. Our results provide strong support for a group that includes protostome (but not deuterostome) coelomate, pseudocoelomate, and acoelomate animals, thus clearly rejecting the Coelomata hypothesis. Two other groups of bilaterian animals, Lophotrochozoa and Ambulacraria, are also supported by our analyses. However, due to the remarkable stability of mitochondrial gene arrangements in Deuterostomia and the Ecdysozoa, conclusions on their evolutionary history cannot be drawn.

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.

We have sequenced the mitochondrial DNA (mtDNA) of Hyaloraphidium curvatum, an organism previously classified as a colorless green alga but now recognized as a lower fungus based on molecular data. The 29.97-kbp mitochondrial chromosome is maintained as a monomeric, linear molecule with identical, inverted repeats (1.43 kbp) at both ends, a rare genome architecture in mitochondria. The genome encodes only 14 known mitochondrial proteins, 7 tRNAs, the large subunit rRNA and small subunit rRNA (SSU rRNA), and 3 ORFs. The SSU rRNA is encoded in two gene pieces that are located 8 kbp apart on the mtDNA. Scrambled and fragmented mitochondrial rRNAs are well known from green algae and alveolate protists but are unprecedented in fungi. Protein genes code for apocytochrome b; cytochrome oxidase 1, 2, and 3, NADH dehydrogenase 1, 2, 3, 4, 4L, 5, and 6, and ATP synthase 6, 8, and 9 subunits, and several of these genes are organized in operon-like clusters. The set of seven mitochondrially encoded tRNAs is insufficient to recognize all codons that occur in the mitochondrial protein genes. When taking into account the pronounced codon bias, at least 16 nuclear-encoded tRNAs are assumed to be imported into the mitochondria. Three of the seven predicted mitochondria-encoded tRNA sequences carry mispairings in the first three positions of the acceptor stem. This strongly suggests that these tRNAs are edited by a mechanism similar to the one seen in the fungus Spizellomyces punctatus and the rhizopod amoeba Acanthamoeba castellanii. Our phylogenetic analysis confirms with overwhelming support that H. curvatum is a member of the chytridiomycete fungi, specifically related to the Monoblepharidales.