Interpreting protein abundance in Saccharomyces cerevisiae through relational learning.

MOTIVATION: Proteomic profiles reflect the functional readout of the physiological state of an organism. An increased understanding of what controls and defines protein abundances is of high scientific interest. Saccharomyces cerevisiae is a well-studied model organism, and there is a large amount of structured knowledge on yeast systems biology in databases such as the Saccharomyces Genome Database, and highly curated genome-scale metabolic models like Yeast8. These datasets, the result of decades of experiments, are abundant in information, and adhere to semantically meaningful ontologies. RESULTS: By representing this knowledge in an expressive Datalog database we generated data descriptors using relational learning that, when combined with supervised machine learning, enables us to predict protein abundances in an explainable manner. We learnt predictive relationships between protein abundances, function and phenotype; such as α-amino acid accumulations and deviations in chronological lifespan. We further demonstrate the power of this methodology on the proteins His4 and Ilv2, connecting qualitative biological concepts to quantified abundances. AVAILABILITY AND IMPLEMENTATION: All data and processing scripts are available at the following Github repository: https://github.com/DanielBrunnsaker/ProtPredict.

Genome-scale model of Rhodotorula toruloides metabolism.

The basidiomycete red yeast Rhodotorula toruloides is a promising platform organism for production of biooils. We present rhto-GEM, the first genome-scale model (GEM) of R. toruloides metabolism, that was largely reconstructed using RAVEN toolbox. The model includes 852 genes, 2,731 reactions, and 2,277 metabolites, while lipid metabolism is described using the SLIMEr formalism allowing direct integration of lipid class and acyl chain experimental distribution data. The simulation results confirmed that the R. toruloides model provides valid growth predictions on glucose, xylose, and glycerol, while prediction of genetic engineering targets to increase production of linolenic acid, triacylglycerols, and carotenoids identified genes-some of which have previously been engineered to successfully increase production. This renders rtho-GEM valuable for future studies to improve the production of other oleochemicals of industrial relevance including value-added fatty acids and carotenoids, in addition to facilitate system-wide omics-data analysis in R. toruloides. Expanding the portfolio of GEMs for lipid-accumulating fungi contributes to both understanding of metabolic mechanisms of the oleaginous phenotype but also uncover particularities of the lipid production machinery in R. toruloides.

Interaction of Lactobacillus vini with the ethanol-producing yeasts Dekkera bruxellensis and Saccharomyces cerevisiae.

Lactobacillus vini was recently described as a contaminant in industrial ethanol fermentations and its co-occurrence with Dekkera bruxellensis was noted. We investigated the growth characteristics of L. vini in cocultivation together with either Saccharomyces cerevisiae or D. bruxellensis. Lower cell numbers of both the yeasts and L. vini as well as a decrease in ethanol and lactate formation in mixed batch cultures compared with pure cultures were noted. L. vini formed cell aggregates (flocs) in all cultivation media with different shapes in Man-Rogosa-Sharpe and yeast extract-peptone-dextrose media. Flocs' size and proportion of cells bound to flocs increased with increasing ethanol concentration. In coculture, formation of lactic acid bacteria-yeast cell aggregates consisting of a bacterial core with an outer layer of yeast cells was observed. L. vini-D. bruxellensis flocs had a bigger surface, due to cells protruding from the pseudomycelium. The involvement of mannose residues in the flocculation between L. vini and yeasts was tested. The presence of mannose induced deflocculation in a concentration-dependent manner. Less mannose was required for the deflocculation of D. bruxellensis as compared with S. cerevisiae.

Adaptation of Dekkera bruxellensis to lignocellulose-based substrate.

Adaptation of Dekkera bruxellensis to lignocellulose hydrolysate was investigated. Cells of D. bruxellensis were grown for 72 and 192 H in batch and continuous culture, respectively (adapted cells). Cultivations in semisynthetic medium were run as controls (nonadapted cells). To test the adaptation, cells from these cultures were reinoculated in the lignocellulose medium, and growth and ethanol production characteristics were monitored. Cells adapted to lignocellulose hydrolysate had a shorter lag phase, grew faster, and produced a higher ethanol concentration as compared with nonadapted cells. A stability test showed that after cultivation in rich medium, cells partially lost the adapted phenotype but still showed faster growth and higher ethanol production as compared with nonadapted cells. Because alcohol dehydrogenase genes have been described to be involved in the adaptation to furfural in Saccharomyces cerevisiae, an analogous mechanism of adaptation to lignocelluloses hydrolysate of D. bruxellensis was hypothesized. However, gene expression analysis showed that genes homologous to S. cerevisiae ADH1 were not involved in the adaptation to lignocelluloses hydrolysate in D. bruxellensis.

AutonoMS: Automated Ion Mobility Metabolomic Fingerprinting.

Automation is dramatically changing the nature of laboratory life science. Robotic lab hardware that can perform manual operations with greater speed, endurance, and reproducibility opens an avenue for faster scientific discovery with less time spent on laborious repetitive tasks. A major bottleneck remains in integrating cutting-edge laboratory equipment into automated workflows, notably specialized analytical equipment, which is designed for human usage. Here we present AutonoMS, a platform for automatically running, processing, and analyzing high-throughput mass spectrometry experiments. AutonoMS is currently written around an ion mobility mass spectrometry (IM-MS) platform and can be adapted to additional analytical instruments and data processing flows. AutonoMS enables automated software agent-controlled end-to-end measurement and analysis runs from experimental specification files that can be produced by human users or upstream software processes. We demonstrate the use and abilities of AutonoMS in a high-throughput flow-injection ion mobility configuration with 5 s sample analysis time, processing robotically prepared chemical standards and cultured yeast samples in targeted and untargeted metabolomics applications. The platform exhibited consistency, reliability, and ease of use while eliminating the need for human intervention in the process of sample injection, data processing, and analysis. The platform paves the way toward a more fully automated mass spectrometry analysis and ultimately closed-loop laboratory workflows involving automated experimentation and analysis coupled to AI-driven experimentation utilizing cutting-edge analytical instrumentation. AutonoMS documentation is available at https://autonoms.readthedocs.io.

The influence of nitrate on the physiology of the yeast Dekkera bruxellensis grown under oxygen limitation.

A previous study showed that the use of nitrate by Dekkera bruxellensis might be an advantageous trait when ammonium is limited in sugarcane substrate for ethanol fermentation. The aim of the present work was to evaluate the influence of nitrate on the yeast physiology during cell growth in different carbon sources under oxygen limitation. If nitrate was the sole source of nitrogen, D. bruxellensis cells presented slower growth, diminished sugar consumption and growth-associated ethanol production, when compared to ammonium. These results were corroborated by the increased expression of genes involved in the pentose phosphate (PP) pathway, the tricarboxylic acid (TCA) cycle and ATP synthesis. The presence of ammonium in the mixed medium restored most parameters to the standard conditions. This work may open up a line of investigation to establish the connection between nitrate assimilation and energetic metabolism in D. bruxellensis and their influence on its fermentative capacity in oxygen-limited or oxygen-depleted conditions.

High-throughput metabolomics for the design and validation of a diauxic shift model.

Saccharomyces cerevisiae is a very well studied organism, yet ∼20% of its proteins remain poorly characterized. Moreover, recent studies seem to indicate that the pace of functional discovery is slow. Previous work has implied that the most probable path forward is via not only automation but fully autonomous systems in which active learning is applied to guide high-throughput experimentation. Development of tools and methods for these types of systems is of paramount importance. In this study we use constrained dynamical flux balance analysis (dFBA) to select ten regulatory deletant strains that are likely to have previously unexplored connections to the diauxic shift. We then analyzed these deletant strains using untargeted metabolomics, generating profiles which were then subsequently investigated to better understand the consequences of the gene deletions in the metabolic reconfiguration of the diauxic shift. We show that metabolic profiles can be utilised to not only gaining insight into cellular transformations such as the diauxic shift, but also on regulatory roles and biological consequences of regulatory gene deletion. We also conclude that untargeted metabolomics is a useful tool for guidance in high-throughput model improvement, and is a fast, sensitive and informative approach appropriate for future large-scale functional analyses of genes. Moreover, it is well-suited for automated approaches due to relative simplicity of processing and the potential to make massively high-throughput.

A Genetic Trap in Yeast for Inhibitors of SARS-CoV-2 Main Protease.

The ongoing COVID-19 pandemic urges searches for antiviral agents that can block infection or ameliorate its symptoms. Using dissimilar search strategies for new antivirals will improve our overall chances of finding effective treatments. Here, we have established an experimental platform for screening of small molecule inhibitors of the SARS-CoV-2 main protease in Saccharomyces cerevisiae cells, genetically engineered to enhance cellular uptake of small molecules in the environment. The system consists of a fusion of the Escherichia coli toxin MazF and its antitoxin MazE, with insertion of a protease cleavage site in the linker peptide connecting the MazE and MazF moieties. Expression of the viral protease confers cleavage of the MazEF fusion, releasing the MazF toxin from its antitoxin, resulting in growth inhibition. In the presence of a small molecule inhibiting the protease, cleavage is blocked and the MazF toxin remains inhibited, promoting growth. The system thus allows positive selection for inhibitors. The engineered yeast strain is tagged with a fluorescent marker protein, allowing precise monitoring of its growth in the presence or absence of inhibitor. We detect an established main protease inhibitor by a robust growth increase, discernible down to 1 µM. The system is suitable for robotized large-scale screens. It allows in vivo evaluation of drug candidates and is rapidly adaptable for new variants of the protease with deviant site specificities. IMPORTANCE The COVID-19 pandemic may continue for several years before vaccination campaigns can put an end to it globally. Thus, the need for discovery of new antiviral drug candidates will remain. We have engineered a system in yeast cells for the detection of small molecule inhibitors of one attractive drug target of SARS-CoV-2, its main protease, which is required for viral replication. The ability to detect inhibitors in live cells brings the advantage that only compounds capable of entering the cell and remain stable there will score in the system. Moreover, because of its design in yeast cells, the system is rapidly adaptable for tuning the detection level and eventual modification of the protease cleavage site in the case of future mutant variants of the SARS-CoV-2 main protease or even for other proteases.

Assembly and Analysis of the Genome Sequence of the Yeast Brettanomyces naardenensis CBS 7540.

Brettanomyces naardenensis is a spoilage yeast with potential for biotechnological applications for production of innovative beverages with low alcohol content and high attenuation degree. Here, we present the first annotated genome of B. naardenensis CBS 7540. The genome of B. naardenensis CBS 7540 was assembled into 76 contigs, totaling 11,283,072 nucleotides. In total, 5168 protein-coding sequences were annotated. The study provides functional genome annotation, phylogenetic analysis, and discusses genetic determinants behind notable stress tolerance and biotechnological potential of B. naardenensis.

Identification and characterisation of two high-affinity glucose transporters from the spoilage yeast Brettanomyces bruxellensis.

The yeast Brettanomyces bruxellensis (syn. Dekkera bruxellensis) is an emerging and undesirable contaminant in industrial low-sugar ethanol fermentations that employ the yeast Saccharomyces cerevisiae. High-affinity glucose import in B. bruxellensis has been proposed to be the mechanism by which this yeast can outcompete S. cerevisiae. The present study describes the characterization of two B. bruxellensis genes (BHT1 and BHT3) believed to encode putative high-affinity glucose transporters. In vitro-generated transcripts of both genes as well as the S. cerevisiae HXT7 high-affinity glucose transporter were injected into Xenopus laevis oocytes and subsequent glucose uptake rates were assayed using 14C-labelled glucose. At 0.1 mM glucose, Bht1p was shown to transport glucose five times faster than Hxt7p. pH affected the rate of glucose transport by Bht1p and Bht3p, indicating an active glucose transport mechanism that involves proton symport. These results suggest a possible role for BHT1 and BHT3 in the competitive ability of B. bruxellensis.

Chromosomal genome assembly of the ethanol production strain CBS 11270 indicates a highly dynamic genome structure in the yeast species Brettanomyces bruxellensis.

Here, we present the genome of the industrial ethanol production strain Brettanomyces bruxellensis CBS 11270. The nuclear genome was found to be diploid, containing four chromosomes with sizes of ranging from 2.2 to 4.0 Mbp. A 75 Kbp mitochondrial genome was also identified. Comparing the homologous chromosomes, we detected that 0.32% of nucleotides were polymorphic, i.e. formed single nucleotide polymorphisms (SNPs), 40.6% of them were found in coding regions (i.e. 0.13% of all nucleotides formed SNPs and were in coding regions). In addition, 8,538 indels were found. The total number of protein coding genes was 4897, of them, 4,284 were annotated on chromosomes; and the mitochondrial genome contained 18 protein coding genes. Additionally, 595 genes, which were annotated, were on contigs not associated with chromosomes. A number of genes was duplicated, most of them as tandem repeats, including a six-gene cluster located on chromosome 3. There were also examples of interchromosomal gene duplications, including a duplication of a six-gene cluster, which was found on both chromosomes 1 and 4. Gene copy number analysis suggested loss of heterozygosity for 372 genes. This may reflect adaptation to relatively harsh but constant conditions of continuous fermentation. Analysis of gene topology showed that most of these losses occurred in clusters of more than one gene, the largest cluster comprising 33 genes. Comparative analysis against the wine isolate CBS 2499 revealed 88,534 SNPs and 8,133 indels. Moreover, when the scaffolds of the CBS 2499 genome assembly were aligned against the chromosomes of CBS 11270, many of them aligned completely, some have chunks aligned to different chromosomes, and some were in fact rearranged. Our findings indicate a highly dynamic genome within the species B. bruxellensis and a tendency towards reduction of gene number in long-term continuous cultivation.

Proteome analysis of xylose metabolism in Rhodotorula toruloides during lipid production.

BACKGROUND: Rhodotorula toruloides is a promising platform organism for production of lipids from lignocellulosic substrates. Little is known about the metabolic aspects of lipid production from the lignocellolosic sugar xylose by oleaginous yeasts in general and R. toruloides in particular. This study presents the first proteome analysis of the metabolism of R. toruloides during conversion of xylose to lipids. RESULTS: Rhodotorula toruloides cultivated on either glucose or xylose was subjected to comparative analysis of its growth dynamics, lipid composition, fatty acid profiles and proteome. The maximum growth and sugar uptake rate of glucose-grown R. toruloides cells were almost twice that of xylose-grown cells. Cultivation on xylose medium resulted in a lower final biomass yield although final cellular lipid content was similar between glucose- and xylose-grown cells. Analysis of lipid classes revealed the presence of monoacylglycerol in the early exponential growth phase as well as a high proportion of free fatty acids. Carbon source-specific changes in lipid profiles were only observed at early exponential growth phase, where C18 fatty acids were more saturated in xylose-grown cells. Proteins involved in sugar transport, initial steps of xylose assimilation and NADPH regeneration were among the proteins whose levels increased the most in xylose-grown cells across all time points. The levels of enzymes involved in the mevalonate pathway, phospholipid biosynthesis and amino acids biosynthesis differed in response to carbon source. In addition, xylose-grown cells contained higher levels of enzymes involved in peroxisomal beta-oxidation and oxidative stress response compared to cells cultivated on glucose. CONCLUSIONS: The results obtained in the present study suggest that sugar import is the limiting step during xylose conversion by R. toruloides into lipids. NADPH appeared to be regenerated primarily through pentose phosphate pathway although it may also involve malic enzyme as well as alcohol and aldehyde dehydrogenases. Increases in enzyme levels of both fatty acid biosynthesis and beta-oxidation in xylose-grown cells was predicted to result in a futile cycle. The results presented here are valuable for the development of lipid production processes employing R. toruloides on xylose-containing substrates.

De novo assembly of Dekkera bruxellensis: a multi technology approach using short and long-read sequencing and optical mapping.

BACKGROUND: It remains a challenge to perform de novo assembly using next-generation sequencing (NGS). Despite the availability of multiple sequencing technologies and tools (e.g., assemblers) it is still difficult to assemble new genomes at chromosome resolution (i.e., one sequence per chromosome). Obtaining high quality draft assemblies is extremely important in the case of yeast genomes to better characterise major events in their evolutionary history. The aim of this work is two-fold: on the one hand we want to show how combining different and somewhat complementary technologies is key to improving assembly quality and correctness, and on the other hand we present a de novo assembly pipeline we believe to be beneficial to core facility bioinformaticians. To demonstrate both the effectiveness of combining technologies and the simplicity of the pipeline, here we present the results obtained using the Dekkera bruxellensis genome. METHODS: In this work we used short-read Illumina data and long-read PacBio data combined with the extreme long-range information from OpGen optical maps in the task of de novo genome assembly and finishing. Moreover, we developed NouGAT, a semi-automated pipeline for read-preprocessing, de novo assembly and assembly evaluation, which was instrumental for this work. RESULTS: We obtained a high quality draft assembly of a yeast genome, resolved on a chromosomal level. Furthermore, this assembly was corrected for mis-assembly errors as demonstrated by resolving a large collapsed repeat and by receiving higher scores by assembly evaluation tools. With the inclusion of PacBio data we were able to fill about 5 % of the optical mapped genome not covered by the Illumina data.

Enhanced ethanol production from wheat straw by integrated storage and pre-treatment (ISP).

Integrated storage and pre-treatment (ISP) combines biopreservation of moist material under airtight conditions and pre-treatment. Moist wheat straw was inoculated with the biocontrol yeast Wickerhamomyces anomalus, the xylan degrading yeast Scheffersomyces stipitis or a co-culture of both. The samples and non-inoculated controls were stored at 4 or 15 °C. The non-inoculated controls were heavily contaminated with moulds, in contrast to the samples inoculated with W. anomalus or S. stipitis. These two yeasts were able to grow on wheat straw as sole source of nutrients. When ethanol was produced from moist wheat straw stored for four weeks at 4 °C with S. stipitis, an up to 40% enhanced yield (final yield 0.15 g ethanol per g straw dry weight) was obtained compared to a dry sample (0.107 g/g). In all other moist samples, stored for four weeks at 4 °C or 15 °C, 6-35% higher yields were obtained. Thus, energy efficient bio-preservation can improve the pre-treatment efficiency for lignocellulose biomass, which is a critical bottleneck in its conversion to biofuels.

Transcriptome of the alternative ethanol production strain Dekkera bruxellensis CBS 11270 in sugar limited, low oxygen cultivation.

Dekkera bruxellensis can outcompete Saccharomyces cerevisiae in environments with low sugar concentrations. It is usually regarded as a spoilage yeast but has lately been identified as an alternative ethanol production organism. In this study, global gene expression in the industrial isolate D. bruxellensis CBS 11270 under oxygen and glucose limitation was investigated by whole transcriptome sequencing using the AB SOLiD technology. Among other observations, we noted expression of respiratory complex I NADH-ubiquinone reductase although D. bruxellensis is a Crabtree positive yeast. The observed higher expression of NADH-generating enzymes compared to NAD(+)-generating enzymes might be the reason for the previously observed NADH imbalance and resulting Custer effect in D. bruxellensis. Low expression of genes involved in glycerol production is probably the molecular basis for high efficiency of D. bruxellensis metabolism under nutrient limitation. No D. bruxellensis homologs to the genes involved in the final reactions of glycerol biosynthesis were detected. A high number of expressed sugar transporter genes is consistent with the hypothesis that the competitiveness of D. bruxellensis is due to a higher affinity for the limiting substrate.