Phylogenomic analyses support the monophyly of Taphrinomycotina, including Schizosaccharomyces fission yeasts

Retracing the evolution of Pneumocystis species, with a focus on the human pathogen Pneumocystis jirovecii

SUMMARYEvery human being is presumed to be infected by the fungus Pneumocystis jirovecii at least once in his or her lifetime. This fungus belongs to a large group of species that appear to exclusively infect mammals, with P. jirovecii being the only one known to cause disease in humans. The mystery of P. jirovecii origin and speciation is just beginning to unravel. Here, we provide a review of the major steps of P. jirovecii evolution. The Pneumocystis genus likely originated from soil or plant-associated organisms during the period of Cretaceous ~165 million years ago and successfully shifted to mammals. The transition coincided with a substantial loss of genes, many of which are related to the synthesis of nutrients that can be scavenged from hosts or cell wall components that could be targeted by the mammalian immune system. Following the transition, the Pneumocystis genus cospeciated with mammals. Each species specialized at infecting its own host. Host specialization is presumably built at least partially upon surface glycoproteins, whose protogene was acquired prior to the genus formation. P. jirovecii appeared at ~65 million years ago, overlapping with the emergence of the first primates. P. jirovecii and its sister species P. macacae, which infects macaques nowadays, may have had overlapping host ranges in the distant past. Clues from molecular clocks suggest that P. jirovecii did not cospeciate with humans. Molecular evidence suggests that Pneumocystis speciation involved chromosomal rearrangements and the mounting of genetic barriers that inhibit gene flow among species.

Schizosaccharomyces versatilis represents a distinct evolutionary lineage of fission yeast

The fission yeast species Schizosaccharomyces japonicus is currently divided into two varieties-S. japonicus var. japonicus and S. japonicus var. versatilis. Here we examine the var. versatilis isolate CBS5679. The CBS5679 genome shows 88% identity to the reference genome of S. japonicus var. japonicus at the coding sequence level, with phylogenetic analyses suggesting that it has split from the S. japonicus lineage 25 million years ago. The CBS5679 genome contains a reciprocal translocation between chromosomes 1 and 2, together with several large inversions. The products of genes linked to the major translocation are associated with 'metabolism' and 'cellular assembly' ontology terms. We further show that CBS5679 does not generate viable progeny with the reference strain of S. japonicus. Although CBS5679 shares closer similarity to the 'type' strain of var. versatilis as compared to S. japonicus, it is not identical to the type strain, suggesting population structure within var. versatilis. We recommend that the taxonomic status of S. japonicus var. versatilis is raised, with it being treated as a separate species, Schizosaccharomyces versatilis.

Reinstatement of the fission yeast species Schizosaccharomyces versatilis Wickerham et Duprat, a sibling species of Schizosaccharomyces japonicus

Schizosaccharomyces japonicus Yukawa et Maki (1931) and Schizosaccharomyces versatilis Wickerham et Duprat (1945) have been treated as varieties of S. japonicus or as conspecific, based on various approaches including mating trials and nDNA/nDNA optical reassociation studies. However, the type strains of S. japonicus and S. versatilis differ by five substitutions (99.15% identity) and one 1-bp indel in the sequences of the D1/D2 domain of the 26S rRNA gene, and 23 substitutions (96.3% identity) and 31-bp indels in the sequences of internal transcribed spacer (ITS) of rRNA, suggesting that they may not be conspecific. To reassess their taxonomic status, we conducted mating trials and whole-genome analyses. Mating trials using the type strains showed a strong but incomplete prezygotic sterility barrier, yielding interspecies mating products at two orders of magnitude lower efficiency than intraspecies matings. These mating products, which were exclusively allodiploid hybrids, were unable to undergo the haplontic life cycle of the parents. We generated chromosome-level gap-less genome assemblies for both type strains. Whole genome sequences yielded an average nucleotide identity (ANI) of 86.4%, indicating clear separation of S. japonicus and S. versatilis. Based on these findings, we propose the reinstatement of S. versatilis as a distinct species (holotype strain: CBS 103T and ex-types: NRRL Y-1026, NBRC 1607, ATCC 9987, PYCC 7100; Mycobank no.: 847838).

A high-quality reference genome for the fission yeast Schizosaccharomyces osmophilus

Fission yeasts are an ancient group of fungal species that diverged from each other from tens to hundreds of million years ago. Among them is the preeminent model organism Schizosaccharomyces pombe, which has significantly contributed to our understandings of molecular mechanisms underlying fundamental cellular processes. The availability of the genomes of S. pombe and 3 other fission yeast species S. japonicus, S. octosporus, and S. cryophilus has enabled cross-species comparisons that provide insights into the evolution of genes, pathways, and genomes. Here, we performed genome sequencing on the type strain of the recently identified fission yeast species S. osmophilus and obtained a complete mitochondrial genome and a nuclear genome assembly with gaps only at rRNA gene arrays. A total of 5,098 protein-coding nuclear genes were annotated and orthologs for more than 95% of them were identified. Genome-based phylogenetic analysis showed that S. osmophilus is most closely related to S. octosporus and these 2 species diverged around 16 million years ago. To demonstrate the utility of this S. osmophilus reference genome, we conducted cross-species comparative analyses of centromeres, telomeres, transposons, the mating-type region, Cbp1 family proteins, and mitochondrial genomes. These analyses revealed conservation of repeat arrangements and sequence motifs in centromere cores, identified telomeric sequences composed of 2 types of repeats, delineated relationships among Tf1/sushi group retrotransposons, characterized the evolutionary origins and trajectories of Cbp1 family domesticated transposases, and discovered signs of interspecific transfer of 2 types of mitochondrial selfish elements.

Genome Comparisons of the Fission Yeasts Reveal Ancient Collinear Loci Maintained by Natural Selection

Fission yeasts have a unique life history and exhibit distinct evolutionary patterns from other yeasts. Besides, the species demonstrate stable genome structures despite the relatively fast evolution of their genomic sequences. To reveal what could be the reason for that, comparative genomic analyses were carried out. Our results provided evidence that the structural and sequence evolution of the fission yeasts were correlated. Moreover, we revealed ancestral locally collinear blocks (aLCBs), which could have been inherited from their last common ancestor. These aLCBs proved to be the most conserved regions of the genomes as the aLCBs contain almost eight genes/blocks on average in the same orientation and order across the species. Gene order of the aLCBs is mainly fission-yeast-specific but supports the idea of filamentous ancestors. Nevertheless, the sequences and gene structures within the aLCBs are as mutable as any sequences in other parts of the genomes. Although genes of certain Gene Ontology (GO) categories tend to cluster at the aLCBs, those GO enrichments are not related to biological functions or high co-expression rates, they are, rather, determined by the density of essential genes and Rec12 cleavage sites. These data and our simulations indicated that aLCBs might not only be remnants of ancestral gene order but are also maintained by natural selection.

The Impact of the Opening of Producer Services on the International Competitiveness of Manufacturing Industry

Countries participating in the international division of labor each try to occupy the industrial highlands, obtain competitive advantages, and promote sustainability of economic development. Based on Porter’s Diamond model, it is widely believed that producer services are vital to support the manufacturing industry. Consequently, this paper selects samples of 55 countries and uses data from 2010 to 2017 to empirically test the impact of producer service’s opening on the export competitiveness of the manufacturing industry. The results show that the opening of producer services in a country promotes the improvement of international competitiveness of manufacturing industry, and the improvement effect is more significant in developed countries than in developing countries. Additionally, the negative impact of foreign capital access restrictions on the export competitiveness of the manufacturing industry is greater than the positive impact of service trade opening. It is also found that the restrictions on foreign capital’s opening in the financial sector have the biggest negative impact on manufacturing exports than that of other divisions of the producer service industry. To improve national competitiveness, it is suggested that the integration of service and manufacturing industries should be enhanced and to decrease FDI access to financial and transportation services restrictions.

Characterization of mitogenomes from four Mucorales species and insights into pathogenicity

BACKGROUND

Mucorales, as one major order of Zygomycetes fungi, can infect human beings and cause serious consequence. We have noticed the pathogenicity of Mucorales is closely related to energy metabolism, while mitochondria play the role of energy factories in almost all biological activities.

METHODS

Virulence of M irregularis, M hiemalis, L corymbifera and R arrhizus were verified in Galleria mellonella larvae, as well as mitochondrial gene copies analysed with RT-qPCR. Mitogenomes of the four Mucorales species were sequenced based on illumina NovaSeq technology to study their characteristic features and functional regions.

RESULTS

Variant virulence of M irregularis, M hiemalis, L corymbifera and R arrhizu were verified by clinical retrospective data and our G mellonella infection models, also copies of mitochondrial genes indicated the significant associations with pathogenicity. A total of 274.18 clean reads were generated to be assembled; the complete mitogenomes of the four Mucorales species were obtained with totally different length. After the genomes annotated and compared, M irregularis was found more similar with M hiemalis than those of L corymbifera and R arrhizus, especially the small (rrns) and large (rrnl) subunits of mitochondrial ribosomal RNA (rRNA) genes. The GC content, ncRNAs and the distribution of the SNPs and InDels were also compared, and the GC content rate of fungi seems to be related to the fungal thermal adaptability. In addition, linear mitogenomes of these four Mucorales showed diverse arrangements of orf genes and directionality of some conserved gene elements.

CONCLUSION

This study uncovered the pathogenicity variances among the four Mucorales species and the relationship between their mitogenomic features and clinical pathogenicity. Further studies like spatial structure of mitochondrial genomes and the comprehensive analysis of transcription regulation are needed.

Effects of benzo[a]pyrene toxicity on morphology and ultrastructure of Hordeum sativum

Many studies have been devoted to investigation of toxic benzo(a)pyrene (BaP) compound, but studies involving changes at the cellular level are insufficient to understand the mechanisms of polycyclic aromatic hydrocarbons (PAHs) effect on plants. To study the toxicity of BaP, a model vegetation experiment was conducted on cultivation of spring barley (Hordeum sativum distichum) on artificially polluted BaP soil at different concentrations. The article discusses the intake of BaP from the soil into the plant and its effect on the organismic and cellular levels of plant organization. The BaP content in the organs of spring barley was determined by the method of saponification. With an increase in the concentration of BaP in the soil, its content in plants also rises, which leads to inhibition of growth processes. The BaP content in the green part of Hordeum sativum increased from 0.3 µg kg^-1 in control soil up to 2.6 µg kg^-1 and 16.8 µg kg^-1 under 20 and 400 ng/g BaP applying in soil, as well as in roots: 0.9 µg kg^-1, 7.7 µg kg^-1, 42.8 µg kg^-1, respectively. Using light and electron microscopy, changes in the tissues and cells of plants were found and it was established that accumulation of BaP in plant tissues caused varying degrees of ultrastructural damage depending on the concentration of pollutant. BaP had the greatest effect on the root, significant changes were found in it both at histological and cytological levels, while changes in the leaves were observed only at the cytological level. The results provide significant information about the mechanism of action of BaP on agricultural plants.

A novel Arabidopsis phyllosphere resident Protomyces species and a re-examination of genus Protomyces based on genome sequence data

Protomyces is an understudied genus of yeast-like fungi currently defined as phytopathogens of only Umbelliferae and Compositae. Species relationships and boundaries remain controversial and molecular data are lacking. Of the 82 named Protomyces, we found few recent studies and six available cultures. We previously isolated Protomyces strains from wild Arabidopsis thaliana, a member of Brassicaceae, a family distant from accepted Protomyces hosts. We previously sequenced the genomes of all available Protomyces species, and P. arabidopsidicola sp. nov. strain C29, from Arabidopsis. Phylogenomics suggests this new species occupied a unique position in the genus. Genomic, morphological, and physiological characteristics distinguished P. arabidopsidicola sp. nov. from other Protomyces. Nuclear gene phylogenetic marker analysis suggests actin1 gene DNA sequences could be used with nuclear ribosomal DNA internal transcribed spacer sequences for rapid identification of Protomyces species. Previous studies demonstrated P. arabidopsidicola sp. nov. could persist on the Arabidopsis phyllosphere and Protomyces sequences were discovered on Arabidopsis at multiple sites in different countries. We conclude that the strain C29 represents a novel Protomyces species and propose the name of P. arabidopsidicola sp. nov. Consequently, we propose that Protomyces is not strictly associated only with the previously recognized host plants.

Novakomyces olei sp. nov., the First Member of a Novel Taphrinomycotina Lineage

Taphrinomycotina is the smallest subphylum of the phylum Ascomycota. It is an assemblage of distantly related early diverging lineages of the phylum, comprising organisms with divergent morphology and ecology; however, phylogenomic analyses support its monophyly. In this study, we report the isolation of a yeast strain, which could not be assigned to any of the currently recognised five classes of Taphrinomycotina. The strain of the novel budding species was recovered from extra virgin olive oil and characterised phenotypically by standard methods. The ultrastructure of the cell wall was investigated by transmission electron microscopy. Comparisons of barcoding DNA sequences indicated that the investigated strain is not closely related to any known organism. Tentative phylogenetic placement was achieved by maximum-likelihood analysis of the D1/D2 domain of the nuclear LSU rRNA gene. The genome of the investigated strain was sequenced, assembled, and annotated. Phylogenomic analyses placed it next to the fission Schizosaccharomyces species. To accommodate the novel species, Novakomyces olei, a novel genus Novakomyces, a novel family Novakomycetaceae, a novel order Novakomycetales, and a novel class Novakomycetes is proposed as well. Functional analysis of genes missing in N. olei in comparison to Schizosaccharomyces pombe revealed that they are biased towards biosynthesis of complex organic molecules, regulation of mRNA, and the electron transport chain. Correlating the genome content and physiology among species of Taphrinomycotina revealed some discordance between pheno- and genotype. N. olei produced ascospores in axenic culture preceded by conjugation between two cells. We confirmed that N. olei is a primary homothallic species lacking genes for different mating types.

An Overview of Genomics, Phylogenomics and Proteomics Approaches in Ascomycota

Fungi are among the most successful eukaryotes on Earth: they have evolved strategies to survive in the most diverse environments and stressful conditions and have been selected and exploited for multiple aims by humans. The characteristic features intrinsic of Fungi have required evolutionary changes and adaptations at deep molecular levels. Omics approaches, nowadays including genomics, metagenomics, phylogenomics, transcriptomics, metabolomics, and proteomics have enormously advanced the way to understand fungal diversity at diverse taxonomic levels, under changeable conditions and in still under-investigated environments. These approaches can be applied both on environmental communities and on individual organisms, either in nature or in axenic culture and have led the traditional morphology-based fungal systematic to increasingly implement molecular-based approaches. The advent of next-generation sequencing technologies was key to boost advances in fungal genomics and proteomics research. Much effort has also been directed towards the development of methodologies for optimal genomic DNA and protein extraction and separation. To date, the amount of proteomics investigations in Ascomycetes exceeds those carried out in any other fungal group. This is primarily due to the preponderance of their involvement in plant and animal diseases and multiple industrial applications, and therefore the need to understand the biological basis of the infectious process to develop mechanisms for biologic control, as well as to detect key proteins with roles in stress survival. Here we chose to present an overview as much comprehensive as possible of the major advances, mainly of the past decade, in the fields of genomics (including phylogenomics) and proteomics of Ascomycota, focusing particularly on those reporting on opportunistic pathogenic, extremophilic, polyextremotolerant and lichenized fungi. We also present a review of the mostly used genome sequencing technologies and methods for DNA sequence and protein analyses applied so far for fungi.

Intraspecific Diversity of Fission Yeast Mitochondrial Genomes

The fission yeast Schizosaccharomyces pombe is an important model organism, but its natural diversity and evolutionary history remain under-studied. In particular, the population genomics of the S. pombe mitochondrial genome (mitogenome) has not been thoroughly investigated. Here, we assembled the complete circular-mapping mitogenomes of 192 S. pombe isolates de novo, and found that these mitogenomes belong to 69 nonidentical sequence types ranging from 17,618 to 26,910 bp in length. Using the assembled mitogenomes, we identified 20 errors in the reference mitogenome and discovered two previously unknown mitochondrial introns. Analyzing sequence diversity of these 69 types of mitogenomes revealed two highly distinct clades, with only three mitogenomes exhibiting signs of inter-clade recombination. This diversity pattern suggests that currently available S. pombe isolates descend from two long-separated ancestral lineages. This conclusion is corroborated by the diversity pattern of the recombination-repressed K-region located between donor mating-type loci mat2 and mat3 in the nuclear genome. We estimated that the two ancestral S. pombe lineages diverged about 31 million generations ago. These findings shed new light on the evolution of S. pombe and the data sets generated in this study will facilitate future research on genome evolution.

Fungal evolution: diversity, taxonomy and phylogeny of the Fungi

The fungal kingdom comprises a hyperdiverse clade of heterotrophic eukaryotes characterized by the presence of a chitinous cell wall, the loss of phagotrophic capabilities and cell organizations that range from completely unicellular monopolar organisms to highly complex syncitial filaments that may form macroscopic structures. Fungi emerged as a 'Third Kingdom', embracing organisms that were outside the classical dichotomy of animals versus vegetals. The taxonomy of this group has a turbulent history that is only now starting to be settled with the advent of genomics and phylogenomics. We here review the current status of the phylogeny and taxonomy of fungi, providing an overview of the main defined groups. Based on current knowledge, nine phylum-level clades can be defined: Opisthosporidia, Chytridiomycota, Neocallimastigomycota, Blastocladiomycota, Zoopagomycota, Mucoromycota, Glomeromycota, Basidiomycota and Ascomycota. For each group, we discuss their main traits and their diversity, focusing on the evolutionary relationships among the main fungal clades. We also explore the diversity and phylogeny of several groups of uncertain affinities and the main phylogenetic and taxonomical controversies and hypotheses in the field.

Lipolytic Enzymes with Hydrolytic and Esterification Activities Produced by Filamentous Fungi Isolated from Decomposition Leaves in an Aquatic Environment

Microbial lipases are prominent biocatalysts able to catalyze a wide variety of reactions in aqueous and nonaqueous media. In this work, filamentous fungi isolated from leaves decomposed in an aquatic environment were screened for lipase production with hydrolytic activity and esterification. Agar plates with Tween 20 and Rhodamine B were used for selection, while submerged cultures with olive oil were subsequently used to select 38 filamentous fungi. Trichoderma harzianum, Fusarium solani, Trichoderma harzianum F5, and Penicillium sp. F36 were grown in six different culture media. F. solani presented the highest lipase production (2.37 U/mL) with esterification activity of 0.07 U/mL using medium composed of (g.L-1) KH2PO4 1.00, MgSO4 H2O 1.123, and CuSO4 0.06. Supplementation of this culture medium with organic nitrogen sources increased lipase production by 461.3% using tryptone and by 419.4% using yeast extract. Among the vegetable oils from the Amazon region, degummed cotton oil induced lipase production up to 8.14 U/mL. The lipase produced by F. solani F61 has great potential to application in conventional processes and biodiesel production by transesterification of vegetable oils, as well as food industries in the production of fatty acid esters by hydrolysis and esterification.

Production of exopolysaccharide by strains of Lactobacillus plantarum YO175 and OF101 isolated from traditional fermented cereal beverage

Lactobacillus plantarum YO175 and OF101 isolates from Nigerian traditional fermented cereal gruel 'ogi', were investigated on the basis of their capability to produce exopolysaccharide (EPS) on sucrose modified deMan Rogosa Sharpe medium (mMRS). Functional groups analysis of the EPSs produced (EPS-YO175 and EPS-OF101) by Fourier-transform infrared (FT-IR) spectroscopy revealed the presence of -OH, C=O and C-H groups. The chemical composition of EPS-YO175 and EPS-OF101 showed the presence of 87.1% and 80.62% carbohydrates and 1.21% and 1.47% protein. For maximum EPS yield, three significant factors were optimized using central composite design and response surface methodology, the predicted maximum EPS produced was 1.38 g/L and 2.19 g/L, while the experimental values were 1.36 g/L and 2.18 g/L for EPS-YO175 and EPS-OF101. The EPS samples showed strong antioxidant activities in-vitro. The scale-up of the production process of the EPS will find its potential application in food industries.

Genomics and evolution of Pneumocystis species

The genus Pneumocystis comprises highly diversified fungal species that cause severe pneumonia in individuals with a deficient immune system. These fungi infect exclusively mammals and present a strict host species specificity. These species have co-diverged with their hosts for long periods of time (> 100 MYA). Details of their biology and evolution are fragmentary mainly because of a lack of an established long-term culture system. Recent genomic advances have unlocked new areas of research and allow new hypotheses to be tested. We review here new findings of the genomic studies in relation with the evolutionary trajectory of these fungi and discuss the impact of genomic data analysis in the context of the population genetics. The combination of slow genome decay and limited expansion of specific gene families and introns reflect intimate interactions of these species with their hosts. The evolutionary adaptation of these organisms is profoundly influenced by their population structure, which in turn is determined by intrinsic features such as their self-fertilizing mating system, high host specificity, long generation times, and transmission mode. Essential key questions concerning their adaptation and speciation remain to be answered. The next cornerstone will consist in the establishment of a long-term culture system and genetic manipulation that should allow unravelling the driving forces of Pneumocystis species evolution.

Phylogenetic taxon definitions for Fungi, Dikarya, Ascomycota and Basidiomycota

Phylogenetic taxon definitions (PTDs) are explicit, phylogeny-based statements that specify clades. PTDs are central to the system of rank-free classification that is governed by the PhyloCode, but they can also be used to clarify the meanings of ranked names. We present PTDs for four major groups: Fungi, Dikarya, Ascomycota, and Basidiomycota.

Advances in Fungal Phylogenomics and Their Impact on Fungal Systematics

In the past decade, advances in next-generation sequencing technologies and bioinformatic pipelines for phylogenomic analysis have led to remarkable progress in fungal systematics and taxonomy. A number of long-standing questions have been addressed using comparative analysis of genome sequence data, resulting in robust multigene phylogenies. These have added to, and often surpassed traditional morphology or single-gene phylogenetic methods. In this chapter, we provide a brief history of fungal systematics and highlight some examples to demonstrate the impact of phylogenomics on this field. We conclude by discussing some of the challenges and promises in fungal biology posed by the ongoing genomics revolution.

The mitochondrial genome of the plant-pathogenic fungus Stemphylium lycopersici uncovers a dynamic structure due to repetitive and mobile elements

Stemphylium lycopersici (Pleosporales) is a plant-pathogenic fungus that has been associated with a broad range of plant-hosts worldwide. It is one of the causative agents of gray leaf spot disease in tomato and pepper. The aim of this work was to characterize the mitochondrial genome of S. lycopersici CIDEFI-216, to use it to trace taxonomic relationships with other fungal taxa and to get insights into the evolutionary history of this phytopathogen. The complete mitochondrial genome was assembled into a circular double-stranded DNA molecule of 75,911 bp that harbors a set of 37 protein-coding genes, 2 rRNA genes (rns and rnl) and 28 tRNA genes, which are transcribed from both sense and antisense strands. Remarkably, its gene repertoire lacks both atp8 and atp9, contains a free-standing gene for the ribosomal protein S3 (rps3) and includes 13 genes with homing endonuclease domains that are mostly located within its 15 group I introns. Strikingly, subunits 1 and 2 of cytochrome oxidase are encoded by a single continuous open reading frame (ORF). A comparative mitogenomic analysis revealed the large extent of structural rearrangements among representatives of Pleosporales, showing the plasticity of their mitochondrial genomes. Finally, an exhaustive phylogenetic analysis of the subphylum Pezizomycotina based on mitochondrial data reconstructed their relationships in concordance with several studies based on nuclear data. This is the first report of a mitochondrial genome belonging to a representative of the family Pleosporaceae.

Changes in nucleosome formation at gene promoters in the archiascomycetous yeast Saitoella complicata

We measured the degree of nucleosome formation at the gene promoters in trichostatin A-treated (1, 2, and 3 µg/mL) cells of the archiascomycete Saitoella complicata and those in enlarged S. complicata cells after zymolyase treatment. TSA-treated and enlarged cells showed similar changes in nucleosome occupancy in five out of six positions in the gene promoters. These results suggest that changes in nucleosome formation at the gene promoters could serve as stress response mechanisms elicited in response to spheroplast (zymolyase treatment) and TSA treatment. In addition, we demonstrated that changes in nucleosome position occurred mainly in cells treated with 1 µg/mL TSA, whereas cells treated with 2 and 3 µg/mL TSA did not exhibit significant changes in the degree of nucleosome formation.

Reconstructing the Backbone of the Saccharomycotina Yeast Phylogeny Using Genome-Scale Data

Understanding the phylogenetic relationships among the yeasts of the subphylum Saccharomycotina is a prerequisite for understanding the evolution of their metabolisms and ecological lifestyles. In the last two decades, the use of rDNA and multilocus data sets has greatly advanced our understanding of the yeast phylogeny, but many deep relationships remain unsupported. In contrast, phylogenomic analyses have involved relatively few taxa and lineages that were often selected with limited considerations for covering the breadth of yeast biodiversity. Here we used genome sequence data from 86 publicly available yeast genomes representing nine of the 11 known major lineages and 10 nonyeast fungal outgroups to generate a 1233-gene, 96-taxon data matrix. Species phylogenies reconstructed using two different methods (concatenation and coalescence) and two data matrices (amino acids or the first two codon positions) yielded identical and highly supported relationships between the nine major lineages. Aside from the lineage comprised by the family Pichiaceae, all other lineages were monophyletic. Most interrelationships among yeast species were robust across the two methods and data matrices. However, eight of the 93 internodes conflicted between analyses or data sets, including the placements of: the clade defined by species that have reassigned the CUG codon to encode serine, instead of leucine; the clade defined by a whole genome duplication; and the species Ascoidea rubescens These phylogenomic analyses provide a robust roadmap for future comparative work across the yeast subphylum in the disciplines of taxonomy, molecular genetics, evolutionary biology, ecology, and biotechnology. To further this end, we have also provided a BLAST server to query the 86 Saccharomycotina genomes, which can be found at http://y1000plus.org/blast.

Phylogenetic Resolution of Deep Eukaryotic and Fungal Relationships Using Highly Conserved Low-Copy Nuclear Genes

A comprehensive and reliable eukaryotic tree of life is important for many aspects of biological studies from comparative developmental and physiological analyses to translational medicine and agriculture. Both gene-rich and taxon-rich approaches are effective strategies to improve phylogenetic accuracy and are greatly facilitated by marker genes that are universally distributed, well conserved, and orthologous among divergent eukaryotes. In this article, we report the identification of 943 low-copy eukaryotic genes and we show that many of these genes are promising tools in resolving eukaryotic phylogenies, despite the challenges of determining deep eukaryotic relationships. As a case study, we demonstrate that smaller subsets of ∼20 and 52 genes could resolve controversial relationships among widely divergent taxa and provide strong support for deep relationships such as the monophyly and branching order of several eukaryotic supergroups. In addition, the use of these genes resulted in fungal phylogenies that are congruent with previous phylogenomic studies that used much larger datasets, and successfully resolved several difficult relationships (e.g., forming a highly supported clade with Microsporidia, Mitosporidium and Rozella sister to other fungi). We propose that these genes are excellent for both gene-rich and taxon-rich analyses and can be applied at multiple taxonomic levels and facilitate a more complete understanding of the eukaryotic tree of life.

Scaling up discovery of hidden diversity in fungi: impacts of barcoding approaches

The fungal kingdom is a hyperdiverse group of multicellular eukaryotes with profound impacts on human society and ecosystem function. The challenge of documenting and describing fungal diversity is exacerbated by their typically cryptic nature, their ability to produce seemingly unrelated morphologies from a single individual and their similarity in appearance to distantly related taxa. This multiplicity of hurdles resulted in the early adoption of DNA-based comparisons to study fungal diversity, including linking curated DNA sequence data to expertly identified voucher specimens. DNA-barcoding approaches in fungi were first applied in specimen-based studies for identification and discovery of taxonomic diversity, but are now widely deployed for community characterization based on sequencing of environmental samples. Collectively, fungal barcoding approaches have yielded important advances across biological scales and research applications, from taxonomic, ecological, industrial and health perspectives. A major outstanding issue is the growing problem of 'sequences without names' that are somewhat uncoupled from the traditional framework of fungal classification based on morphology and preserved specimens. This review summarizes some of the most significant impacts of fungal barcoding, its limitations, and progress towards the challenge of effective utilization of the exponentially growing volume of data gathered from high-throughput sequencing technologies.This article is part of the themed issue 'From DNA barcodes to biomes'.

A Survey of the Gene Repertoire of Gigaspora rosea Unravels Conserved Features among Glomeromycota for Obligate Biotrophy

Arbuscular mycorrhizal (AM) fungi are a diverse group of soil fungi (Glomeromycota) that form the most ancient mutualistic association termed AM symbiosis with a majority of land plants, improving their nutrition uptake and resistance to stresses. In contrast to their great ecological implications, the knowledge of the molecular biological mechanisms involved is still scant, partly due to the limited genomic resources available. Here, we describe the gene repertoire of a new AM fungus Gigaspora rosea (Diversisporales). Among the 86332 non-redundant virtual transcripts assembled, 15346 presented similarities with proteins in the Refseq database and 10175 were assigned with GO terms. KOG and Interpro domain annotations clearly showed an enrichment of genes involved in signal transduction in G. rosea. KEGG pathway analysis indicates that most primary metabolic processes are active in G. rosea. However, as for Rhizophagus irregularis, several metabolic genes were not found, including the fatty acid synthase (FAS) gene. This finding supports the hypothesis that AM fungi depend on the lipids produced by their hosts. Furthermore, the presence of a large number of transporters and 100s of secreted proteins, together with the reduced number of plant cell wall degrading enzymes could be interpreted as an evolutionary adaptation to its mutualistic obligate biotrophy. The detection of meiosis-related genes suggests that G. rosea might use a cryptic sexual process. Lastly, a phylogeny of basal fungi clearly shows Glomeromycota as a sister clade to Mucoromycotina, not only to the Mucorales or Mortierellales. The characterization of the gene repertoire from an AM fungal species belonging to the order of Diversisporales and its comparison with the gene sets of R. irregularis (Glomerales) and Gigaspora margarita (Diversisporales), reveal that AM fungi share several features linked to mutualistic obligate biotrophy. This work contributes to lay the foundation for forthcoming studies into the genomics of Diversisporales, and also illuminates the utility of comparing gene repertoires of species from Diversisporales and other clades of Glomeromycota to gain more insights into the genetics and evolution of this fungal group.

Draft Genome Sequence of the Archiascomycetous Yeast Saitoella complicata

The draft genome sequence of the archiasomycetous yeast Saitoella complicata was determined. The assembly of newly and previously sequenced data sets resulted in 104 contigs (total of 14.1 Mbp; N₅₀, 239 kbp). On the newly assembled genome, a total of 6,933 protein-coding sequences (7,119 transcripts, including alternative splicing forms) were identified.

The early diverging ascomycetous budding yeast Saitoella complicata has three histone deacetylases belonging to the Clr6, Hos2, and Rpd3 lineages

We sequenced the genomic DNA and the transcribed RNA of the ascomycetous budding yeast Saitoella complicata, which belongs to the earliest lineage (Taphrinomycotina) of ascomycetes. We found 3 protein-coding regions similar to Clr6 of Schizosaccharomyces (a member of Taphrinomycotina). Clr6 has a structure similar to that of Rpd3 and Hos2 of Saccharomyces. These proteins belong to the class 1 histone deacetylase (HDAC) family. The phylogenetic tree showed that the Clr6, Hos2, and Rpd3 lineages are separated in fungal HDACs. Basidiomycetes have 3 proteins belonging to the Clr6, Hos2, and Rpd3 lineages. On the other hand, whereas ascomycetes except for Schizosaccharomyces have the Hos2 and Rpd3 homologs, and lack the Clr6 homolog, Schizosaccharomyces has the Clr6 and Hos2 homologs, and lacks the Rpd3 homolog. Interestingly, Pneumocystis and Saitoella have 3 proteins belonging to the Clr6, Hos2, and Rpd3 lineages. Thus, these fungi are the first ascomycete found to possess all 3 types. Our findings indicated that Taphrinomycotina has conserved the Clr6 protein, suggesting that the ancestor of Dikarya (ascomycetes and basidiomycetes) had 3 proteins belonging to the Clr6, Hos2, and Rpd3 lineages. During ascomycete evolution, Pezizomycotina and Saccharomycotina appear to have lost their Clr6 homologs and Schizosaccharomyces to have lost its Rpd3 homolog.

Mitochondrial genome evolution in yeasts: an all-encompassing view

Mitochondria are important organelles that harbor their own genomes encoding a key set of proteins that ensure respiration and provide the eukaryotic cell with energy. Recent advances in high-throughput sequencing technologies present a unique opportunity to explore mitochondrial (mt) genome evolution. The Saccharomycotina yeasts have proven to be the leading organisms for mt comparative and population genomics. In fact, the explosion of complete yeast mt genome sequences has allowed for a broader view of the mt diversity across this incredibly diverse subphylum, both within and between closely related species. Here, we summarize the present state of yeast mitogenomics, including the currently available data and what it reveals concerning the diversity of content, organization, structure and evolution of mt genomes.

Comparative genomics of Taphrina fungi causing varying degrees of tumorous deformity in plants

Taphrina fungi are biotrophic plant pathogens that cause plant deformity diseases. We sequenced the genomes of four Taphrina species—Taphrina wiesneri, T. deformans, T. flavorubra, and T. populina—which parasitize Prunus, Cerasus, and Populus hosts with varying severity of disease symptoms. High levels of gene synteny within Taphrina species were observed, and our comparative analysis further revealed that these fungi may utilize multiple strategies in coping with the host environment that are also found in some specialized dimorphic species. These include species-specific aneuploidy and clusters of highly diverged secreted proteins located at subtelomeres. We also identified species differences in plant hormone biosynthesis pathways, which may contribute to varying degree of disease symptoms. The genomes provide a rich resource for investigation into Taphrina biology and evolutionary studies across the basal ascomycetes clade.

Complementation of a manganese-dependent superoxide dismutase-deficient yeast strain with Pneumocystis carinii sod2 gene

Manganese-dependent superoxide dismutase (MnSOD) is one of the key enzymes involved in the cellular defense against oxidative stress. Previously, the Pneumocystis carinii sod2 gene (Pcsod2) was isolated and characterized. Based on protein sequence comparison, Pcsod2 was suggested to encode a putative MnSOD protein likely to be targeted into the mitochondrion. In this work, the Pcsod2 was cloned and expressed as a recombinant protein in EG110 Saccharomyces cerevisiae strain lacking the MnSOD-coding gene (Scsod2) in order to investigate the function and subcellular localization of P. carinii MnSOD (PcMnSOD). The Pcsod2 gene was amplified by PCR and cloned into the pYES2.1/V5-His-TOPO(®) expression vector. The recombinant construct was then transformed into EG110 strain. Once its expression had been induced, PcMnSOD was able to complement the growth defect of EG110 yeast cells that had been exposed to the redox-cycling compound menadione. N-term sequencing of the PcMnSOD protein allowed identifying the cleavage site of a mitochondrial targeting peptide. Immune-colocalization of PcMnSOD and yeast CoxIV further confirmed the mitochondrial localization of the PcMnSOD. Heterologous expression of PcMnSOD in yeast indicates that Pcsod2 encodes an active MnSOD, targeted to the yeast mitochondrion that allows the yeast cells to grow in the presence of reactive oxygen species (ROS).

Comparative genomics of pneumocystis species suggests the absence of genes for myo-inositol synthesis and reliance on inositol transport and metabolism

In the context of deciphering the metabolic strategies of the obligate pathogenic fungi in the genus Pneumocystis, the genomes of three species (P. carinii, P. murina, and P. jirovecii) were compared among themselves and with the free-living, phylogenetically related fission yeast (Schizosaccharomyces pombe). The underrepresentation of amino acid metabolism pathways compared to those in S. pombe, as well as the incomplete steroid biosynthesis pathway, were confirmed for P. carinii and P. jirovecii and extended to P. murina. All three Pneumocystis species showed overrepresentation of the inositol phosphate metabolism pathway compared to that in the fission yeast. In addition to those known in S. pombe, four genes, encoding inositol-polyphosphate multikinase (EC 2.7.1.151), inositol-pentakisphosphate 2-kinase (EC 2.7.1.158), phosphoinositide 5-phosphatase (EC 3.1.3.36), and inositol-1,4-bisphosphate 1-phosphatase (EC 3.1.3.57), were identified in the two rodent Pneumocystis genomes, P. carinii and P. murina. The P. jirovecii genome appeared to contain three of these genes but lacked phosphoinositide 5-phosphatase. Notably, two genes encoding enzymes essential for myo-inositol synthesis, inositol-1-phosphate synthase (INO1) and inositol monophosphatase (INM1), were absent from all three genomes, suggesting that Pneumocystis species are inositol auxotrophs. In keeping with the need to acquire exogenous inositol, two genes with products homologous to fungal inositol transporters, ITR1 and ITR2, were identified in P. carinii and P. murina, while P. jirovecii contained only the ITR1 homolog. The ITR and inositol metabolism genes in P. murina and P. carinii were expressed during fulminant infection as determined by reverse transcriptase real-time PCR of cDNA from infected lung tissue. Supplementation of in vitro culture with inositol yielded significant improvement of the viability of P. carinii for days 7 through 14.

IMPORTANCE

Microbes in the genus Pneumocystis are obligate pathogenic fungi that reside in mammalian lungs and cause Pneumocystis pneumonia in hosts with weakened immune systems. These fungal infections are not responsive to standard antifungal therapy. A long-term in vitro culture system is not available for these fungi, impeding the study of their biology and genetics and new drug development. Given that all genomes of the Pneumocystis species analyzed lack the genes for inositol synthesis and contain inositol transporters, Pneumocystis fungi, like S. pombe, appear to be inositol auxotrophs. Inositol is important for the pathogenesis, virulence, and mating processes in Candida albicans and Cryptococcus neoformans, suggesting similar importance within the Pneumocystis species as well. This is the first report to (i) characterize genes in the inositol phosphate metabolism and transport pathways in Pneumocystis species and (ii) identify inositol as a supplement that improved the viability of P. carinii in in vitro culture.

A dually located multi-HMG-box protein of Aspergillus nidulans has a crucial role in conidial and ascospore germination

Seven HMG-box proteins of Aspergillus nidulans have been identified in the genomic databases. Three of these have the characteristics of non-specific DNA-binding proteins. One of these, AN1267 (HmbB), comprises one canonical HMG-box in its C-terminus and upstream of the canonical box two structurally related boxes, to be called Shadow-HMG-boxes. This protein defines, together with the Podospora anserina mtHMG1, a clade of proteins present in the Pezizomycotina, with orthologues in some of the Taphrinomycotina. HmbB localizes primarily to the mitochondria but occasionally in nuclei. The deletion of the cognate gene results in a number of pleiotropic effects, including those on hyphal morphology, sensitivity to oxidative stress, absence of sterigmatocystin production and changes in the profile of conidial metabolites. The most striking phenotype of deletion strains is a dramatic decrease in conidial and ascospore viability. We show that this is most likely due to the protein being essential to maintain mitochondrial DNA in spores.

Use of gene sequence analyses and genome comparisons for yeast systematics

Detection, identification and classification of yeasts have undergone a major transformation in the past decade and a half following application of gene sequence analyses and genome comparisons. Development of a database (barcode) of easily determined gene sequences from domains 1 and 2 (D1/D2) of large subunit rRNA and from the internal transcribed spacer (ITS) now permits many laboratories to identify species accurately and this has led to a doubling in the number of known species of yeasts over the past decade. Phylogenetic analysis of gene sequences has resulted in major revision of yeast systematics, resulting in redefinition of nearly all genera. Future work calls for application of genomics to refine our understanding of the species concept and to provide a better understanding of the boundaries of genera and higher levels of classification. This increased understanding of phylogeny is expected to allow prediction of the genetic potential of various clades and species for biotechnological applications and adaptation to environmental changes.

Sequencing and comparison of the mitochondrial COI gene from isolates of Arbuscular Mycorrhizal Fungi belonging to Gigasporaceae and Glomeraceae families

Arbuscular Mycorrhizal Fungi (AMF) are well known for their ecological importance and their positive influence on plants. The genetics and phylogeny of this group of fungi have long been debated. Nuclear markers are the main tools used for phylogenetic analyses, but they have sometimes proved difficult to use because of their extreme variability. Therefore, the attention of researchers has been moving towards other genomic markers, in particular those from the mitochondrial DNA. In this study, 46 sequences of different AMF isolates belonging to two main clades Gigasporaceae and Glomeraceae have been obtained from the mitochondrial gene coding for the Cytochrome c Oxidase I (COI), representing the largest dataset to date of AMF COI sequences. A very low level of divergence was recorded in the COI sequences from the Gigasporaceae, which could reflect either a slow rate of evolution or a more recent evolutionary divergence of this group. On the other hand, the COI sequence divergence between Gigasporaceae and Glomeraceae was high, with synonymous divergence reaching saturated levels. This work also showed the difficulty in developing valuable mitochondrial markers able to effectively distinguish all Glomeromycota species, especially those belonging to Gigasporaceae, yet it represents a first step towards the development of a full mtDNA-based dataset which can be used for further phylogenetic investigations of this fungal phylum.

UBL1 of Fusarium verticillioides links the N-end rule pathway to extracellular sensing and plant pathogenesis

Fusarium verticillioides produces fumonisin mycotoxins during colonization of maize. Currently, molecular mechanisms underlying responsiveness of F.verticillioides to extracellular cues during pathogenesis are poorly understood. In this study, insertional mutants were created and screened to identify genes involved in responses to extracellular starch. In one mutant, the restriction enzyme-mediated integration cassette disrupted a gene (UBL1) encoding a UBR-Box/RING domain E3 ubiquitin ligase involved in the N-end rule pathway. Disruption of UBL1 in F.verticillioides (Δubl1) influenced conidiation, hyphal morphology, pigmentation and amylolysis. Disruption of UBL1 also impaired kernel colonization, but the ratio of fumonisin B1 per unit growth was not significantly reduced. The inability of a Δubl1 mutant to recognize an N-end rule degron confirmed involvement of UBL1 in the N-end rule pathway. Additionally, Ubl1 physically interacted with two G protein α subunits of F.verticillioides, thus implicating UBL1 in G protein-mediated sensing of the external environment. Furthermore, deletion of the UBL1 orthologue in F.graminearum reduced virulence on wheat and maize, thus indicating that UBL1 has a broader role in virulence among Fusarium species. This study provides the first linkage between the N-end rule pathway and fungal pathogenesis, and illustrates a new mechanism through which fungi respond to the external environment.

Toward genome-enabled mycology

Genome-enabled mycology is a rapidly expanding field that is characterized by the pervasive use of genome-scale data and associated computational tools in all aspects of fungal biology. Genome-enabled mycology is integrative and often requires teams of researchers with diverse skills in organismal mycology, bioinformatics and molecular biology. This issue of Mycologia presents the first complete fungal genomes in the history of the journal, reflecting the ongoing transformation of mycology into a genome-enabled science. Here, we consider the prospects for genome-enabled mycology and the technical and social challenges that will need to be overcome to grow the database of complete fungal genomes and enable all fungal biologists to make use of the new data.

Phylogenomic analyses of nuclear genes reveal the evolutionary relationships within the BEP clade and the evidence of positive selection in Poaceae

BEP clade of the grass family (Poaceae) is composed of three subfamilies, i.e. Bambusoideae, Ehrhartoideae, and Pooideae. Controversies on the phylogenetic relationships among three subfamilies still persist in spite of great efforts. However, previous evidence was mainly provided from plastid genes with only a few nuclear genes utilized. Given different evolutionary histories recorded by plastid and nuclear genes, it is indispensable to uncover their relationships based on nuclear genes. Here, eleven species with whole-sequenced genome and six species with transcriptomic data were included in this study. A total of 121 one-to-one orthologous groups (OGs) were identified and phylogenetic trees were reconstructed by different tree-building methods. Genes which might have undergone positive selection and played important roles in adaptive evolution were also investigated from 314 and 173 one-to-one OGs in two bamboo species and 14 grass species, respectively. Our results support the ((B, P) E) topology with high supporting values. Besides, our findings also indicate that 24 and nine orthologs with statistically significant evidence of positive selection are mainly involved in abiotic and biotic stress response, reproduction and development, plant metabolism and enzyme etc. from two bamboo species and 14 grass species, respectively. In summary, this study demonstrates the power of phylogenomic approach to shed lights on the evolutionary relationships within the BEP clade, and offers valuable insights into adaptive evolution of the grass family.

Genome sequencing of the plant pathogen Taphrina deformans, the causal agent of peach leaf curl

Taphrina deformans is a fungus responsible for peach leaf curl, an important plant disease. It is phylogenetically assigned to the Taphrinomycotina subphylum, which includes the fission yeast and the mammalian pathogens of the genus Pneumocystis. We describe here the genome of T. deformans in the light of its dual plant-saprophytic/plant-parasitic lifestyle. The 13.3-Mb genome contains few identifiable repeated elements (ca. 1.5%) and a relatively high GC content (49.5%). A total of 5,735 protein-coding genes were identified, among which 83% share similarities with other fungi. Adaptation to the plant host seems reflected in the genome, since the genome carries genes involved in plant cell wall degradation (e.g., cellulases and cutinases), secondary metabolism, the hallmark glyoxylate cycle, detoxification, and sterol biosynthesis, as well as genes involved in the biosynthesis of plant hormones. Genes involved in lipid metabolism may play a role in its virulence. Several locus candidates for putative MAT cassettes and sex-related genes akin to those of Schizosaccharomyces pombe were identified. A mating-type-switching mechanism similar to that found in ascomycetous yeasts could be in effect. Taken together, the findings are consistent with the alternate saprophytic and parasitic-pathogenic lifestyles of T. deformans.

IMPORTANCE

Peach leaf curl is an important plant disease which causes significant losses of fruit production. We report here the genome sequence of the causative agent of the disease, the fungus Taphrina deformans. The genome carries characteristic genes that are important for the plant infection process. These include (i) proteases that allow degradation of the plant tissues; (ii) secondary metabolites which are products favoring interaction of the fungus with the environment, including the host; (iii) hormones that are responsible for the symptom of severely distorted leaves on the host; and (iv) drug detoxification enzymes that confer resistance to fungicides. The availability of the genome allows the design of new drug targets as well as the elaboration of specific management strategies to fight the disease.

Examining new phylogenetic markers to uncover the evolutionary history of early-diverging fungi: comparing MCM7, TSR1 and rRNA genes for single- and multi-gene analyses of the Kickxellomycotina

The recently recognised protein-coding genes MCM7 and TSR1 have shown significant promise for phylogenetic resolution within the Ascomycota and Basidiomycota, but have remained unexamined within other fungal groups (except for Mucorales). We designed and tested primers to amplify these genes across early-diverging fungal clades, with emphasis on the Kickxellomycotina, zygomycetous fungi with characteristic flared septal walls forming pores with lenticular plugs. Phylogenetic tree resolution and congruence with MCM7 and TSR1 were compared against those inferred with nuclear small (SSU) and large subunit (LSU) rRNA genes. We also combined MCM7 and TSR1 data with the rDNA data to create 3- and 4-gene trees of the Kickxellomycotina that help to resolve evolutionary relationships among and within the core clades of this subphylum. Phylogenetic inference suggests that Barbatospora, Orphella, Ramicandelaber and Spiromyces may represent unique lineages. It is suggested that these markers may be more broadly useful for phylogenetic studies among other groups of early-diverging fungi.

Alternative splicing and subfunctionalization generates functional diversity in fungal proteomes

Alternative splicing is commonly used by the Metazoa to generate more than one protein from a gene. However, such diversification of the proteome by alternative splicing is much rarer in fungi. We describe here an ancient fungal alternative splicing event in which these two proteins are generated from a single alternatively spliced ancestral SKI7/HBS1 gene retained in many species in both the Ascomycota and Basidiomycota. While the ability to express two proteins from a single SKI7/HBS1 gene is conserved in many fungi, the exact mechanism by which they achieve this varies. The alternative splicing was lost in Saccharomyces cerevisiae following the whole-genome duplication event as these two genes subfunctionalized into the present functionally distinct HBS1 and SKI7 genes. When expressed in yeast, the single gene from Lachancea kluyveri generates two functionally distinct proteins. Expression of one of these proteins complements hbs1, but not ski7 mutations, while the other protein complements ski7, but not hbs1. This is the first known case of subfunctionalization by loss of alternative splicing in yeast. By coincidence, the ancestral alternatively spliced gene was also duplicated in Schizosaccharomyces pombe with subsequent subfunctionalization and loss of splicing. Similar subfunctionalization by loss of alternative splicing in fungi also explains the presence of two PTC7 genes in the budding yeast Tetrapisispora blattae, suggesting that this is a common mechanism to preserve duplicate alternatively spliced genes.

The role of duplicated genes in originating new functions is an important question in evolution. Almost all species have duplicated genes that carry out similar but not identical functions. Similar proteins that perform different functions can also be generated when one gene generates multiple mRNAs by alternative splicing that are translated into multiple similar proteins. This alternative splicing is prevalent in animal cells, but much rarer in fungi. Here we show that most fungi use alternative splicing to make a Ski7 protein and a Hbs1 protein from the same gene. Two fungi, budding yeast and fission yeast, have been much better characterized than other fungi, and co-incidentally they both have duplicated this alternatively spliced gene, resulting in two similar genes that are no longer alternatively spliced. Finally, we describe another example where two duplicate genes replace one alternatively spliced gene, suggesting that this is a common mechanism to divide functions among duplicate genes.

Sequencing and characterization of the complete mitochondrial genomes of three Pneumocystis species provide new insights into divergence between human and rodent Pneumocystis

Pneumocystis jirovecii is an important opportunistic pathogen associated with AIDS and other immunodeficient conditions. Currently, very little is known about its nuclear and mitochondrial genomes. In this study, we sequenced the complete mitochondrial genome (mtDNA) of this organism and its closely related species Pneumocystis carinii and Pneumocystis murina by a combination of sequencing technologies. Our study shows that P. carinii and P. murina mtDNA share a nearly identical number and order of genes in a linear configuration, whereas P. jirovecii has a circular mtDNA containing nearly the same set of genes but in a different order. Detailed studies of the mtDNA terminal structures of P. murina and P. carinii suggest a unique replication mechanism for linear mtDNA. Phylogenetic analysis supports a close association of Pneumocystis species with Taphrina, Saitoella, and Schizosaccharomyces, and divergence within Pneumocystis species, with P. murina and P. carinii being more closely related to each other than either is to P. jirovecii. Comparative analysis of four complete P. jirovecii mtDNA sequences in this study and previously reported mtDNA sequences for diagnosing and genotyping suggests that the current diagnostic and typing methods can be improved using the complete mtDNA data. The availability of the complete P. jirovecii mtDNA also opens the possibility of identifying new therapeutic targets.

The common ancestral core of vertebrate and fungal telomerase RNAs

Telomerase is a ribonucleoprotein with an intrinsic telomerase RNA (TER) component. Within yeasts, TER is remarkably large and presents little similarity in secondary structure to vertebrate or ciliate TERs. To better understand the evolution of fungal telomerase, we identified 74 TERs from Pezizomycotina and Taphrinomycotina subphyla, sister clades to budding yeasts. We initially identified TER from Neurospora crassa using a novel deep-sequencing-based approach, and homologous TER sequences from available fungal genome databases by computational searches. Remarkably, TERs from these non-yeast fungi have many attributes in common with vertebrate TERs. Comparative phylogenetic analysis of highly conserved regions within Pezizomycotina TERs revealed two core domains nearly identical in secondary structure to the pseudoknot and CR4/5 within vertebrate TERs. We then analyzed N. crassa and Schizosaccharomyces pombe telomerase reconstituted in vitro, and showed that the two RNA core domains in both systems can reconstitute activity in trans as two separate RNA fragments. Furthermore, the primer-extension pulse-chase analysis affirmed that the reconstituted N. crassa telomerase synthesizes TTAGGG repeats with high processivity, a common attribute of vertebrate telomerase. Overall, this study reveals the common ancestral cores of vertebrate and fungal TERs, and provides insights into the molecular evolution of fungal TER structure and function.

Group I intron-mediated trans-splicing in mitochondria of Gigaspora rosea and a robust phylogenetic affiliation of arbuscular mycorrhizal fungi with Mortierellales

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.