Exploring contradictory phylogenetic relationships in yeasts

Fungi under Modified Atmosphere-The Effects of CO2 Stress on Cell Membranes and Description of New Yeast Stenotrophomyces fumitolerans gen. nov., sp. nov

High levels of carbon dioxide are known to inhibit the growth of microorganisms. A total of twenty strains of filamentous fungi and yeasts were isolated from habitats with enriched carbon dioxide concentration. Most strains were derived from modified atmosphere packed (MAP) food products or mofettes and were cultivated under an atmosphere of 20% CO2 and 80% O2. The influence of CO2 on fungal cell membrane fatty acid profiles was examined in this study. Major changes were the increase in linolenic acid (C18:3 cis 9, 12, 15) and, additionally in most strains, linoleic acid (C18:2 cis 9, 12) with a maximum of 24.8%, at the expense of oleic (C18:1 cis 9), palmitic (C16:0), palmitoleic (C16:1 cis 9) and stearic acid (C18:0). The degree of fatty acid unsaturation increased for all of the strains in the study, which consequently led to lower melting temperatures of the cell membranes after incubation with elevated levels of CO2, indicating fluidization of the membrane and a potential membrane malfunction. Growth was reduced in 18 out of 20 strains in laboratory experiments and a change in pigmentation was observed in several strains. Two of the isolated strains, strain WT5 and strain WR1, were found to represent a hitherto undescribed yeast for which the new genus and species Stenotrophomyces fumitolerans (MB# 849906) is proposed.

Gene tree discord, simplex plots, and statistical tests under the coalescent

A simple graphical device, the simplex plot of quartet concordance factors, is introduced to aid in the exploration of a collection of gene trees on a common set of taxa. A single plot summarizes all gene tree discord, and allows for visual comparison to the expected discord from the multispecies coalescent model (MSC) of incomplete lineage sorting on a species tree. A formal statistical procedure is described that can quantify the deviation from expectation for each subset of four taxa, suggesting when the data is not in accord with the MSC, and thus that either gene tree inference error is substantial or a more complex model such as that on a network may be required. If the collection of gene trees is in accord with the MSC, the plots reveal when substantial incomplete lineage sorting is present. Applications to both simulated and empirical multilocus data sets illustrate the insights provided.

NANUQ: a method for inferring species networks from gene trees under the coalescent model

Species networks generalize the notion of species trees to allow for hybridization or other lateral gene transfer. Under the network multispecies coalescent model, individual gene trees arising from a network can have any topology, but arise with frequencies dependent on the network structure and numerical parameters. We propose a new algorithm for statistical inference of a level-1 species network under this model, from data consisting of gene tree topologies, and provide the theoretical justification for it. The algorithm is based on an analysis of quartets displayed on gene trees, combining several statistical hypothesis tests with combinatorial ideas such as a quartet-based intertaxon distance appropriate to networks, the NeighborNet algorithm for circular split systems, and the Circular Network algorithm for constructing a splits graph.

Insights into the phylogeny of Northern Hemisphere Armillaria: Neighbor-net and Bayesian analyses of translation elongation factor 1-α gene sequences

Armillaria possesses several intriguing characteristics that have inspired wide interest in understanding phylogenetic relationships within and among species of this genus. Nuclear ribosomal DNA sequence-based analyses of Armillaria provide only limited information for phylogenetic studies among widely divergent taxa. More recent studies have shown that translation elongation factor 1-α (tef1) sequences are highly informative for phylogenetic analysis of Armillaria species within diverse global regions. This study used Neighbor-net and coalescence-based Bayesian analyses to examine phylogenetic relationships of newly determined and existing tef1 sequences derived from diverse Armillaria species from across the Northern Hemisphere, with Southern Hemisphere Armillaria species included for reference. Based on the Bayesian analysis of tef1 sequences, Armillaria species from the Northern Hemisphere are generally contained within the following four superclades, which are named according to the specific epithet of the most frequently cited species within the superclade: (i) Socialis/Tabescens (exannulate) superclade including Eurasian A. ectypa, North American A. socialis (A. tabescens), and Eurasian A. socialis (A. tabescens) clades; (ii) Mellea superclade including undescribed annulate North American Armillaria sp. (Mexico) and four separate clades of A. mellea (Europe and Iran, eastern Asia, and two groups from North America); (iii) Gallica superclade including Armillaria Nag E (Japan), multiple clades of A. gallica (Asia and Europe), A. calvescens (eastern North America), A. cepistipes (North America), A. altimontana (western USA), A. nabsnona (North America and Japan), and at least two A. gallica clades (North America); and (iv) Solidipes/Ostoyae superclade including two A. solidipes/ostoyae clades (North America), A. gemina (eastern USA), A. solidipes/ostoyae (Eurasia), A. cepistipes (Europe and Japan), A. sinapina (North America and Japan), and A. borealis (Eurasia) clade 2. Of note is that A. borealis (Eurasia) clade 1 appears basal to the Solidipes/Ostoyae and Gallica superclades. The Neighbor-net analysis showed similar phylogenetic relationships. This study further demonstrates the utility of tef1 for global phylogenetic studies of Armillaria species and provides critical insights into multiple taxonomic issues that warrant further study.

Zygosaccharomyces rouxii strains CECT 11923 and Z. rouxii CECT 10425: Two new putative hybrids?

Based on IGS-PCR RFLP polymorphism, we previously detected two Z. rouxii strains (CECT 11923 and CECT 10425) that clustered with hybrid strains (NCYC 1682, NCYC 3060 and NCYC 3061). Given the recently recognized important industrial role of hybrids, their detection is very useful. Based on the IGS1 rDNA region alignment of hybrid strains and the Z. rouxii CECT 11923 and CECT 10425, in this work, we developed a pair of Zygosaccharomyces hybrid-specific primers, HibZF/HibZR. Positive amplicons were only obtained in the Zygosaccharomyces spp. hybrids included in this study and the CECT 11923 and CECT 10425 strains analyzed here. In the present study, we applied molecular tools to highlight the nature of these strains; they are quite different from each other as well as from Z. rouxii type strain. Based on the presence of two heterologous copies of nuclear-encoded genes (SOD2 and HIS3), the sequences of divergent 5.8S-ITS rDNA, D1/D2 26S rDNA copies and, the amplification with species-specific primer for Z. rouxii and Z. pseudorouxii, we hypothesize that the CECT 11923 strain might be a hybrid strain. Whereas, CECT 10425, the sequence analysis of 5.8S-ITS rDNA and D1/D2 26S rDNA copies presented 99-100% sequence identity with Zygosaccharomyces sp. NBRC 10669 (LN849119.1) and Z. sapae ABT 301^T. Nevertheless, we discard that it could be a Z. sapae strain based on the results obtained in this study. Namely, the amplification with hybrid-specific primer designed in this study, the number of divergent copies of HIS3 (2), the fact that it only possesses one SOD2 gene and the amplification with species-specific primer for Z. pseudorouxii, therefore it could be a new species or a hybrid strain.

Multigene assessment of the species boundaries and sexual status of the basidiomycetous yeasts Cryptococcus flavescens and C. terrestris (Tremellales)

Cryptococcus flavescens and C. terrestris are phenotypically indistinguishable sister species that belong to the order Tremellales (Tremellomycetes, Basidiomycota) and which may be mistaken for C. laurentii based on phenotype. Phylogenetic separation between C. flavescens and C. terrestris was based on rDNA sequence analyses, but very little is known on their intraspecific genetic variability or propensity for sexual reproduction. We studied 59 strains from different substrates and geographic locations, and used a multilocus sequencing (MLS) approach complemented with the sequencing of mating type (MAT) genes to assess genetic variation and reexamine the boundaries of the two species, as well as their sexual status. The following five loci were chosen for MLS: the rDNA ITS-LSU region, the rDNA IGS1 spacer, and fragments of the genes encoding the largest subunit of RNA polymerase II (RPB1), the translation elongation factor 1 alpha (TEF1) and the p21-activated protein kinase (STE20). Phylogenetic network analyses confirmed the genetic separation of the two species and revealed two additional cryptic species, for which the names Cryptococcus baii and C. ruineniae are proposed. Further analyses of the data revealed a high degree of genetic heterogeneity within C. flavescens as well as evidence for recombination between lineages detected for this species. Strains of C. terrestris displayed higher levels of similarity in all analysed genes and appear to make up a single recombining group. The two MAT genes (STE3 and SXI1/SXI2) sequenced for C. flavescens strains confirmed the potential for sexual reproduction and suggest the presence of a tetrapolar mating system with a biallelic pheromone/receptor locus and a multiallelic HD locus. In C. terrestris we could only sequence STE3, which revealed a biallelic P/R locus. In spite of the strong evidence for sexual recombination in the two species, attempts at mating compatible strains of both species on culture media were unsuccessful.

Komagataella kurtzmanii sp. nov., a new sibling species of Komagataella (Pichia) pastoris based on multigene sequence analysis

A novel methanol assimilating yeast species Komagataella kurtzmanii is described using the type strain VKPM Y-727 (=KBP Y-2878 = UCD-FST 76-20 = Starmer #75-208.2 = CBS 12817 = NRRL Y-63667) isolated by W.T. Starmer from a fir flux in the Catalina Mountains, Southern AZ, USA. The new species is registered in MycoBank under MB 803919. The species was differentiated by divergence in gene sequences for D1/D2 LSU rRNA, ITS1-5.8S-ITS2, RNA polymerase subunit I, translation elongation factor-1α and mitochondrial small subunit rRNA. K. kurtzmanii differs from its phenotypically similar sibling species Komagataella pastoris, Komagataella pseudopastoris, Komagataella phaffii, Komagataella populi and Komagataella ulmi by absence of growth at 35 °C and inability to assimilate trehalose.

Networks in a large-scale phylogenetic analysis: reconstructing evolutionary history of Asparagales (Lilianae) based on four plastid genes

Phylogenetic analysis aims to produce a bifurcating tree, which disregards conflicting signals and displays only those that are present in a large proportion of the data. However, any character (or tree) conflict in a dataset allows the exploration of support for various evolutionary hypotheses. Although data-display network approaches exist, biologists cannot easily and routinely use them to compute rooted phylogenetic networks on real datasets containing hundreds of taxa. Here, we constructed an original neighbour-net for a large dataset of Asparagales to highlight the aspects of the resulting network that will be important for interpreting phylogeny. The analyses were largely conducted with new data collected for the same loci as in previous studies, but from different species accessions and greater sampling in many cases than in published analyses. The network tree summarised the majority data pattern in the characters of plastid sequences before tree building, which largely confirmed the currently recognised phylogenetic relationships. Most conflicting signals are at the base of each group along the Asparagales backbone, which helps us to establish the expectancy and advance our understanding of some difficult taxa relationships and their phylogeny. The network method should play a greater role in phylogenetic analyses than it has in the past. To advance the understanding of evolutionary history of the largest order of monocots Asparagales, absolute diversification times were estimated for family-level clades using relaxed molecular clock analyses.

Networks in a large-scale phylogenetic analysis: reconstructing evolutionary history of Asparagales (Lilianae) based on four plastid genes

Phylogenetic analysis aims to produce a bifurcating tree, which disregards conflicting signals and displays only those that are present in a large proportion of the data. However, any character (or tree) conflict in a dataset allows the exploration of support for various evolutionary hypotheses. Although data-display network approaches exist, biologists cannot easily and routinely use them to compute rooted phylogenetic networks on real datasets containing hundreds of taxa. Here, we constructed an original neighbour-net for a large dataset of Asparagales to highlight the aspects of the resulting network that will be important for interpreting phylogeny. The analyses were largely conducted with new data collected for the same loci as in previous studies, but from different species accessions and greater sampling in many cases than in published analyses. The network tree summarised the majority data pattern in the characters of plastid sequences before tree building, which largely confirmed the currently recognised phylogenetic relationships. Most conflicting signals are at the base of each group along the Asparagales backbone, which helps us to establish the expectancy and advance our understanding of some difficult taxa relationships and their phylogeny. The network method should play a greater role in phylogenetic analyses than it has in the past. To advance the understanding of evolutionary history of the largest order of monocots Asparagales, absolute diversification times were estimated for family-level clades using relaxed molecular clock analyses.

Phylogenetic origin and transcriptional regulation at the post-diauxic phase of SPI1, in Saccharomyces cerevisiae

The gene SPI1, of Saccharomyces cerevisiae, encodes a cell wall protein that is induced in several stress conditions, particularly in the postdiauxic and stationary phases of growth. It has a paralogue, SED1, which shows some common features in expression regulation and in the null mutant phenotype. In this work we have identified homologues in other species of yeasts and filamentous fungi, and we have also elucidated some aspects of the origin of SPI1, by duplication and diversification of SED1. In terms of regulation, we have found that the expression in the post-diauxic phase is regulated by genes related to the PKA pathway and stress response (MSN2/4, YAK1, POP2, SOK2, PHD1, and PHO84) and by genes involved in the PKC pathway (WSC2, PKC1, and MPK1).

Mitochondrial genomes of yeasts of the Yarrowia clade

Candida alimentaria, Candida deformans, Candida galli, and Candida phangngensis have been recently reported to be the close relatives of Yarrowia lipolytica. To explore this clade of yeasts, we sequenced the mitochondrial genome (mtDNA) of these four species and compared it with the mtDNA of Y. lipolytica. The five mtDNAs exhibit a similar architecture and a high level of similarity of protein coding sequences. Genome sizes are variable, ranging from 28 017 bp in C. phangngensis to 48 508 bp in C. galli, mainly because of the variations in intron size and number. All introns are of group I, except for a group II intron inserted in the cob gene of a single species, C. galli. Putative endonuclease coding sequences were present in most group I introns, but also twice as free-standing ORFs in C. galli. Phylogenetic relationships of the five species were explored using protein alignments. No close relative of the Yarrowia clade could be identified, but protein and rRNA gene orders were partially conserved in the mtDNA of Candida salmanticensis.

Addressing inter-gene heterogeneity in maximum likelihood phylogenomic analysis: yeasts revisited

Phylogenomic approaches to the resolution of inter-species relationships have become well established in recent years. Often these involve concatenation of many orthologous genes found in the respective genomes followed by analysis using standard phylogenetic models. Genome-scale data promise increased resolution by minimising sampling error, yet are associated with well-known but often inappropriately addressed caveats arising through data heterogeneity and model violation. These can lead to the reconstruction of highly-supported but incorrect topologies. With the aim of obtaining a species tree for 18 species within the ascomycetous yeasts, we have investigated the use of appropriate evolutionary models to address inter-gene heterogeneities and the scalability and validity of supermatrix analysis as the phylogenetic problem becomes more difficult and the number of genes analysed approaches truly phylogenomic dimensions. We have extended a widely-known early phylogenomic study of yeasts by adding additional species to increase diversity and augmenting the number of genes under analysis. We have investigated sophisticated maximum likelihood analyses, considering not only a concatenated version of the data but also partitioned models where each gene constitutes a partition and parameters are free to vary between the different partitions (thereby accounting for variation in the evolutionary processes at different loci). We find considerable increases in likelihood using these complex models, arguing for the need for appropriate models when analyzing phylogenomic data. Using these methods, we were able to reconstruct a well-supported tree for 18 ascomycetous yeasts spanning about 250 million years of evolution.

Complete sequence and analysis of the mitochondrial genome of the methylotrophic yeast Hansenula polymorpha DL-1

We determined the complete nucleotide sequence of the 41 719 bp mitochondrial genome of the methylotrophic yeast Hansenula polymorpha strain DL-1. It contains genes for three subunits of cytochrome oxidase (cox1, cox2 and cox3), three subunits of ATP synthase (atp6, atp8 and atp9), seven subunits of NADH dehydrogenase (nad1-6 and nad4L), apocytochrome b (cob), four endonuclease/maturase homologs, a ribosomal protein (rps3), large and small rRNAs and a complete set of tRNAs. The structural genes are organized in two major transcriptional units. Phylogenetic, gene content and gene order analyses revealed the close phylogenetic relationship between H. polymorpha and Brettanomyces custersianus, and support the assignment of strain DL-1 to a separate genus rather than including it in the polyphyletic genus Pichia.

Species-specific PCR primers for the rapid identification of yeasts of the genus Zygosaccharomyces

Species-specific primer pairs that produce a single band of known product size have been developed for members of the Zygosaccharomyces clade including Zygosaccharomyces bailii, Zygosaccharomyces bisporus, Zygosaccharomyces kombuchaensis, Zygosaccharomyces lentus, Zygosaccharomyces machadoi, Zygosaccharomyces mellis and Zygosaccharomyces rouxii. An existing primer pair for the provisional new species Zygosaccharomyces pseudorouxii has been confirmed as specific. The HIS3 gene, encoding imidazole-glycerolphosphate dehydratase, was used as the target gene. This housekeeping gene evolves slowly and is thus well conserved among different isolates, but shows a significant number of base pair changes between even closely related species, sufficient for species-specific primer design. The primers were tested on type and wild strains of the genus Zygosaccharomyces and on members of the Saccharomycetaceae. Sequencing of the D1/D2 region of rDNA was used to confirm the identification of all nonculture collection isolates. This approach used extracted genomic DNA, but in practice, it can be used efficiently with a rapid colony PCR protocol. The method also successfully detected known and new hybrid strains of Z. rouxii and Z. pseudorouxii. The method is rapid, robust and inexpensive. It requires little expertise by the user and is thus useful for preliminary, large-scale screens.

Mitochondrial genome from the facultative anaerobe and petite-positive yeast Dekkera bruxellensis contains the NADH dehydrogenase subunit genes

The progenitor of the Dekkera/Brettanomyces clade separated from the Saccharomyces/Kluyveromyces clade over 200 million years ago. However, within both clades, several lineages developed similar physiological traits. Both Saccharomyces cerevisiae and Dekkera bruxellensis are facultative anaerobes; in the presence of excess oxygen and sugars, they accumulate ethanol (Crabtree effect) and they both spontaneously generate respiratory-deficient mutants (petites). In order to understand the role of respiratory metabolism, the mitochondrial DNA (mtDNA) molecules of two Dekkera/Brettanomyces species were analysed. Dekkera bruxellensis mtDNA shares several properties with S. cerevisiae, such as the large genome size (76 453 bp), and the organization of the intergenic sequences consisting of spacious AT-rich regions containing a number of hairpin GC-rich cluster-like elements. In addition to a basic set of the mitochondrial genes coding for the components of cytochrome oxidase, cytochrome b, subunits of ATPase, two rRNA subunits and 25 tRNAs, D. bruxellensis also carries genes for the NADH dehydrogenase complex. Apparently, in yeast, the loss of this complex is not a precondition to develop a petite-positive, Crabtree-positive and anaerobic nature. On the other hand, mtDNA from a petite-negative Brettanomyces custersianus is much smaller (30 058 bp); it contains a similar gene set and has only short intergenic sequences.

Unifying vertical and nonvertical evolution: a stochastic ARG-based framework

Evolutionary biologists have introduced numerous statistical approaches to explore nonvertical evolution, such as horizontal gene transfer, recombination, and genomic reassortment, through collections of Markov-dependent gene trees. These tree collections allow for inference of nonvertical evolution, but only indirectly, making findings difficult to interpret and models difficult to generalize. An alternative approach to explore nonvertical evolution relies on phylogenetic networks. These networks provide a framework to model nonvertical evolution but leave unanswered questions such as the statistical significance of specific nonvertical events. In this paper, we begin to correct the shortcomings of both approaches by introducing the "stochastic model for reassortment and transfer events" (SMARTIE) drawing upon ancestral recombination graphs (ARGs). ARGs are directed graphs that allow for formal probabilistic inference on vertical speciation events and nonvertical evolutionary events. We apply SMARTIE to phylogenetic data. Because of this, we can typically infer a single most probable ARG, avoiding coarse population dynamic summary statistics. In addition, a focus on phylogenetic data suggests novel probability distributions on ARGs. To make inference with our model, we develop a reversible jump Markov chain Monte Carlo sampler to approximate the posterior distribution of SMARTIE. Using the BEAST phylogenetic software as a foundation, the sampler employs a parallel computing approach that allows for inference on large-scale data sets. To demonstrate SMARTIE, we explore 2 separate phylogenetic applications, one involving pathogenic Leptospirochete and the other Saccharomyces.

Unusual composition of a yeast chromosome arm is associated with its delayed replication

The 11.3-Mb genome of the yeast Lachancea (Saccharomyces) kluyveri displays an intriguing compositional heterogeneity: a region of approximately 1 Mb, covering almost the whole left arm of chromosome C (C-left), has an average GC content of 52.9%, which is significantly higher than the 40.4% global GC content of the rest of the genome. This region contains the MAT locus, which remains normal in composition. The excess of GC base pairs affects both coding and noncoding sequences, and thus is not due to selective pressure acting on protein sequences. It leads to a strong codon usage bias and alters the amino acid composition of the 457 proteins encoded on C-left that do not show obvious bias for functional categories, or the presence of paralogs or orthologs of essential genes of Saccharomyces cerevisiae. They share significant synteny conservation with other species of the Saccharomycetaceae, and phylogenetic analysis indicates that C-left originates from a Lachancea species. In contrast, there is a complete absence of transposable elements in C-left, whereas 18 elements per megabase are distributed across the rest of the genome. Comparative hybridization of synchronized cells using high-density genome arrays reveals that C-left is replicated later during S phase than the rest of the genome. Two possible primary causes of this major compositional heterogeneity are discussed: an ancient hybridization of two related species with very distinct GC composition, or an intrinsic mechanism, possibly associated with the loss of the silent cassettes from C-left that progressively increased the GC content and generated the delayed replication of this chromosomal arm.