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

Characterization and phylogenetic analysis of the complete mitochondrial genome of Saccharomycopsis fibuligera (lindner) Klocker 1907 (saccharomycetales: saccharomycopsidaceae)

Saccharomycopsis fibuligera (Lindner) Klocker 1907 is frequently employed in the fermentation of metabolites such as citric acid, ethanol, mannitol, and pyruvate. Its heat tolerance and alcohol-producing capabilities during fermentation make it a desirable option for bread and wine production. To date, the mitochondrial genome of S. fibuligera has not been sequenced. In the present study, we obtained the full mitochondrial genome of S. fibuligera, which is 57,302 bp long and has a GC content of 24.40%. This genome contained 14 core protein-coding genes, 3 independent ORFs, 21 intronic ORFs, 25 tRNAs, and 2 rRNA genes. By utilizing the Bayesian inference phylogenetic method, we constructed phylogenetic trees for 24 Saccharomycotina fungi, which indicated that S. fibuligera is closely related to S. capsularis.

Contribution of Complex I NADH Dehydrogenase to Respiratory Energy Coupling in Glucose-Grown Cultures of Ogataea parapolymorpha

The thermotolerant yeast Ogataea parapolymorpha (formerly Hansenula polymorpha) is an industrially relevant production host that exhibits a fully respiratory sugar metabolism in aerobic batch cultures. NADH-derived electrons can enter its mitochondrial respiratory chain either via a proton-translocating complex I NADH-dehydrogenase or via three putative alternative NADH dehydrogenases. This respiratory entry point affects the amount of ATP produced per NADH/O2 consumed and therefore impacts the maximum yield of biomass and/or cellular products from a given amount of substrate. To investigate the physiological importance of complex I, a wild-type O. parapolymorpha strain and a congenic complex I-deficient mutant were grown on glucose in aerobic batch, chemostat, and retentostat cultures in bioreactors. In batch cultures, the two strains exhibited a fully respiratory metabolism and showed the same growth rates and biomass yields, indicating that, under these conditions, the contribution of NADH oxidation via complex I was negligible. Both strains also exhibited a respiratory metabolism in glucose-limited chemostat cultures, but the complex I-deficient mutant showed considerably reduced biomass yields on substrate and oxygen, consistent with a lower efficiency of respiratory energy coupling. In glucose-limited retentostat cultures at specific growth rates down to ∼0.001 h-1, both O. parapolymorpha strains showed high viability. Maintenance energy requirements at these extremely low growth rates were approximately 3-fold lower than estimated from faster-growing chemostat cultures, indicating a stringent-response-like behavior. Quantitative transcriptome and proteome analyses indicated condition-dependent expression patterns of complex I subunits and of alternative NADH dehydrogenases that were consistent with physiological observations.IMPORTANCE Since popular microbial cell factories have typically not been selected for efficient respiratory energy coupling, their ATP yields from sugar catabolism are often suboptimal. In aerobic industrial processes, suboptimal energy coupling results in reduced product yields on sugar, increased process costs for oxygen transfer, and volumetric productivity limitations due to limitations in gas transfer and cooling. This study provides insights into the contribution of mechanisms of respiratory energy coupling in the yeast cell factory Ogataea parapolymorpha under different growth conditions and provides a basis for rational improvement of energy coupling in yeast cell factories. Analysis of energy metabolism of O. parapolymorpha at extremely low specific growth rates indicated that this yeast reduces its energy requirements for cellular maintenance under extreme energy limitation. Exploration of the mechanisms for this increased energetic efficiency may contribute to an optimization of the performance of industrial processes with slow-growing eukaryotic cell factories.

The diversity of mtDNA rns introns among strains of Ophiostoma piliferum, Ophiostoma pluriannulatum and related species

Background

Based on previous studies, it was suspected that the mitochondrial rns gene within the Ophiostomatales is rich in introns. This study focused on a collection of strains representing Ophiostoma piliferum, Ophiostoma pluriannulatum and related species that cause blue-stain; these fungi colonize the sapwood of trees and impart a dark stain. This reduces the value of the lumber. The goal was to examine the mtDNA rns intron landscape for these important blue stain fungi in order to facilitate future annotation of mitochondrial genomes (mtDNA) and to potentially identify mtDNA introns that can encode homing endonucleases which may have applications in biotechnology.

Results

Comparative sequence analysis identified five intron insertion sites among the ophiostomatoid fungi examined. Positions mS379 and mS952 harbor group II introns, the mS379 intron encodes a reverse transcriptase, and the mS952 intron encodes a potential homing endonuclease. Positions mS569, mS1224, and mS1247 have group I introns inserted and these encode intact or eroded homing endonuclease open reading frames (ORF). Phylogenetic analysis of the intron ORFs showed that they can be found in the same insertion site in closely and distantly related species.

Conclusions

Based on the molecular markers examined (rDNA internal transcribed spacers and rns introns), strains representing O. pilifera, O. pluriannulatum and Ophiostoma novae-zelandiae could not be resolved. Phylogenetic studies suggest that introns are gained and lost and that horizontal transfer could explain the presence of related intron in distantly related fungi. With regard to the mS379 group II intron, this study shows that mitochondrial group II introns and their reverse transcriptases may also follow the life cycle previously proposed for group I introns and their homing endonucleases. This consists of intron invasion, decay of intron ORF, loss of intron, and possible reinvasion.

Electronic supplementary material

The online version of this article (doi:10.1186/s40064-016-3076-6) contains supplementary material, which is available to authorized users.

High variability of mitochondrial gene order among fungi

From their origin as an early alpha proteobacterial endosymbiont to their current state as cellular organelles, large-scale genomic reorganization has taken place in the mitochondria of all main eukaryotic lineages. So far, most studies have focused on plant and animal mitochondrial (mt) genomes (mtDNA), but fungi provide new opportunities to study highly differentiated mtDNAs. Here, we analyzed 38 complete fungal mt genomes to investigate the evolution of mtDNA gene order among fungi. In particular, we looked for evidence of nonhomologous intrachromosomal recombination and investigated the dynamics of gene rearrangements. We investigated the effect that introns, intronic open reading frames (ORFs), and repeats may have on gene order. Additionally, we asked whether the distribution of transfer RNAs (tRNAs) evolves independently to that of mt protein-coding genes. We found that fungal mt genomes display remarkable variation between and within the major fungal phyla in terms of gene order, genome size, composition of intergenic regions, and presence of repeats, introns, and associated ORFs. Our results support previous evidence for the presence of mt recombination in all fungal phyla, a process conspicuously lacking in most Metazoa. Overall, the patterns of rearrangements may be explained by the combined influences of recombination (i.e., most likely nonhomologous and intrachromosomal), accumulated repeats, especially at intergenic regions, and to a lesser extent, mobile element dynamics.

Genome sequence and analysis of methylotrophic yeast Hansenula polymorpha DL1

BACKGROUND

Hansenula polymorpha DL1 is a methylotrophic yeast, widely used in fundamental studies of methanol metabolism, peroxisome biogenesis and function, and also as a microbial cell factory for production of recombinant proteins and metabolic engineering towards the goal of high temperature ethanol production.

RESULTS

We have sequenced the 9 Mbp H. polymorpha DL1 genome and performed whole-genome analysis for the H. polymorpha transcriptome obtained from both methanol- and glucose-grown cells. RNA-seq analysis revealed the complex and dynamic character of the H. polymorpha transcriptome under the two studied conditions, identified abundant and highly unregulated expression of 40% of the genome in methanol grown cells, and revealed alternative splicing events. We have identified subtelomerically biased protein families in H. polymorpha, clusters of LTR elements at G + C-poor chromosomal loci in the middle of each of the seven H. polymorpha chromosomes, and established the evolutionary position of H. polymorpha DL1 within a separate yeast clade together with the methylotrophic yeast Pichia pastoris and the non-methylotrophic yeast Dekkera bruxellensis. Intergenome comparisons uncovered extensive gene order reshuffling between the three yeast genomes. Phylogenetic analyses enabled us to reveal patterns of evolution of methylotrophy in yeasts and filamentous fungi.

CONCLUSIONS

Our results open new opportunities for in-depth understanding of many aspects of H. polymorpha life cycle, physiology and metabolism as well as genome evolution in methylotrophic yeasts and may lead to novel improvements toward the application of H. polymorpha DL-1 as a microbial cell factory.

Specific features of telomerase RNA from Hansenula polymorpha

Telomerase, a ribonucleoprotein, is responsible for the maintenance of eukaryotic genome integrity by replicating the ends of chromosomes. The core enzyme comprises the conserved protein TERT and an RNA subunit (TER) that, in contrast, displays large variations in size and structure. Here, we report the identification of the telomerase RNA from thermotolerant yeast Hansenula polymorpha (HpTER) and describe its structural features. We show further that the H. polymorpha telomerase reverse transcribes the template beyond the predicted boundary and adds a nontelomeric dT in vitro. Sequencing of the chromosomal ends revealed that this nucleotide is specifically present as a terminal nucleotide at the 3' end of telomeres. Mutational analysis of HpTER confirmed that the incorporation of dT functions to limit telomere length in this species.

Mitochondrial genome evolution in a single protoploid yeast species

Mitochondria are organelles, which play a key role in some essential functions, including respiration, metabolite biosynthesis, ion homeostasis, and apoptosis. The vast numbers of mitochondrial DNA (mtDNA) sequences of various yeast species, which have recently been published, have also helped to elucidate the structural diversity of these genomes. Although a large corpus of data are now available on the diversity of yeast species, little is known so far about the mtDNA diversity in single yeast species. To study the genetic variations occurring in the mtDNA of wild yeast isolates, we performed a genome-wide polymorphism survey on the mtDNA of 18 Lachancea kluyveri (formerly Saccharomyces kluyveri) strains. We determined the complete mt genome sequences of strains isolated from various geographical locations (in North America, Asia, and Europe) and ecological niches (Drosophila, tree exudates, soil). The mt genome of the NCYC 543 reference strain is 51,525 bp long. It contains the same core of genes as Lachancea thermotolerans, the nearest relative to L. kluyveri. To explore the mt genome variations in a single yeast species, we compared the mtDNAs of the 18 isolates. The phylogeny and population structure of L. kluyveri provide clear-cut evidence for the existence of well-defined geographically isolated lineages. Although these genomes are completely syntenic, their size and the intron content were found to vary among the isolates studied. These genomes are highly polymorphic, showing an average diversity of 28.5 SNPs/kb and 6.6 indels/kb. Analysis of the SNP and indel patterns showed the existence of a particularly high overall level of polymorphism in the intergenic regions. The dN/dS ratios obtained are consistent with purifying selection in all these genes, with the noteworthy exception of the VAR1 gene, which gave a very high ratio. These data suggest that the intergenic regions have evolved very fast in yeast mitochondrial genomes.

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.

Complete mitochondrial genome of compactin-producing fungus Penicillium solitum and comparative analysis of Trichocomaceae mitochondrial genomes

We determined the complete mitochondrial genome sequence of the compactin-producing fungus Penicillium solitum strain 20-01. The 28 601-base pair circular-mapping DNA molecule encodes a characteristic set of mitochondrial proteins and RNA genes and is intron-free. All 46 protein- and RNA-encoding genes are located on one strand and apparently transcribed in one direction. Comparative analysis of this mtDNA and previously sequenced but unannotated mitochondrial genomes of several medically and industrially important species of the Aspergillus/Penicillium group revealed their extensive similarity in terms of size, gene content and sequence, which is also reflected in the almost perfect conservation of mitochondrial gene order in Penicillium and Aspergillus. Phylogenetic analysis based on concatenated mitochondrial protein sequences confirmed the monophyletic origin of Eurotiomycetes.