Comparative genomics and transcriptomics of Pichia pastoris

Construction of an efficient enzyme-cell@material biocatalyst through the biofilm immobilization of Komagataella phaffii for continuous biocatalysis

The ever-growing demand for cost-effective and green biocatalytic transformations has prompted the rational design and development of robust biocatalytic tools. However, transformations are hindered by limited continuous process and enzymatic instability. Here, 10 Flo family related genes in Komagataella phaffii were systematically evaluated to assess their adhesive properties. For the first time, we identified the KpFlo11C domain of BSC1p as facilitating cell aggregation on carriers, thereby enhancing the biofilm immobilization process. Furthermore, an engineered K. phaffii strain, fixing β-galactosidase on the cell surface, was constructed by optimizing the signal peptide and gene dosage, for enhancing the efficiency of enzyme targeting and anchoring, as well as the proportion of cells displaying the enzyme. Finally, the KpFlo11C domain was overexpressed in this K. phaffii cell display system to construct the enzyme-cell@material biocatalyst, which exhibited robust continuous production of galacto-oligosaccharides (GOS) at a rate of 8.16 g/L/h in a 6-L fermenter. The development of this enzyme-cell@material biocatalyst, which offers a highly efficient, stable, low-cost, and simplified biocatalytic process, provides a direction for the application of other yeasts in large-scale industrial continuous production.

Identification and functional analysis of growth rate associated long non-coding RNAs in Komagataella phaffii

Long non-coding RNAs (lncRNAs) are a regulatory feature that have been reported to operate on both transcriptional and translational levels. With DNA-based prediction still limited, lncRNAs are most reliably identified through genome-guided mapping of RNA-Seq data. Reports of lncRNAs in yeast have been increasing in recent years and changes in their expression levels have often been associated with stressful conditions. As the transition to near zero-growth conditions likely imposes stress, we used RNA-Seq data from the non-conventional, biotechnologically established yeast Komagataella phaffii, cultivated in glucose-limited retentostats, to identify the expression of lncRNAs. Using an adapted bioinformatics pipeline, we identified 168 mostly novel lncRNAs from the K. phaffii retentostat RNA-Seq data, 36 of which demonstrate likely growth-associated expression changes. lncRNA expression levels were associated to that of possible interaction partners, in both cis and trans, suggesting potential roles in regulatory adaptations. Our analysis indicates that lncRNAs likely contribute to how K. phaffii responds to changing environmental conditions, as exemplified here by the adaptation to extremely slow growth.

Ancient origin and high diversity of zymocin-like killer toxins in the budding yeast subphylum

Zymocin is a well-characterized killer toxin secreted by some strains of the yeast Kluyveromyces lactis. It acts by cleaving a specific tRNA in sensitive recipient cells. Zymocin is encoded by a killer plasmid or virus-like element (VLE), which is a linear DNA molecule located in the cytosol. We hypothesized that a tRNA-cleaving toxin similar to zymocin may have caused the three parallel changes to the nuclear genetic code that occurred during yeast evolution, in which the codon CUG became translated as Ser or Ala instead of Leu. However, zymocin-like toxins are rare - both among species, and among strains within a species -and only four toxins of this type have previously been discovered. Here, we identified 45 zymocin-like toxin genes in Saccharomycotina, the budding yeast subphylum, using a bioinformatics strategy, and verified that many of them are toxic to Saccharomyces cerevisiae when expressed. Some of the toxin genes are located on cytosolic VLEs, whereas others are on VLE-derived DNA integrated into the nuclear genome. The toxins are extraordinarily diverse in sequence and show evidence of positive selection. Toxin genes were found in five taxonomic orders of budding yeasts, including two of the three orders that reassigned CUG codons, indicating that VLEs have been parasites of yeast species for at least 300 My and that their existence predates the genetic code changes.

Adaptation of Aglycosylated Monoclonal Antibodies for Improved Production in Komagataella phaffii

Monoclonal antibodies (mAbs) are a major class of biopharmaceuticals manufactured by well-established processes using Chinese Hamster Ovary (CHO) cells. Next-generation biomanufacturing using alternative hosts like Komagataella phaffii could improve the accessibility of these medicines, address broad societal goals for sustainability, and offer financial advantages for accelerated development of new products. Antibodies produced by K. phaffii, however, may manifest unique molecular quality attributes, like host-dependent, product-related variants, that could raise potential concerns for clinical use. We demonstrate here conservative modifications to the amino acid sequence of aglycosylated antibodies based on the human IgG1 isotype that minimize product-related variations when secreted by K. phaffii. A combination of 2-3 changes of amino acids reduced variations across six different aglycosylated versions of commercial mAbs. Expression of a modified sequence of NIST mAb in both K. phaffii and CHO cells showed comparable biophysical properties and molecular variations. These results suggest a path toward the production of high-quality mAbs that could be expressed interchangeably by either yeast or mammalian cells. Improving molecular designs of proteins to enable a range of manufacturing strategies for well-characterized biopharmaceuticals could accelerate global accessibility and innovations.

Optimizing Pichia pastoris protein secretion: Role of N-linked glycosylation on the α-mating factor secretion signal leader

The methylotrophic yeast, Pichia pastoris (P. pastoris; syn. Komagataella spp.), known for its ability to grow to high cell densities, its strong and tightly regulated promoters, and mammalian liked secretion pathway, has been widely used as a robust system to secrete heterologous proteins. The α-mating factor (MF) secretion signal leader from Saccharomyces cerevisiae (S. cerevisiae) is currently the most successfully used secretion signal sequence in the P. pastoris system. In this study, the secretion efficiency mediated by the α-MF secretion signal leaders from Komagataella pastoris (K. pastoris) and Komagataella phaffii (K. phaffii) was assessed using Enhanced Green Fluorescent Protein (EGFP) as a reporter. The results indicated that the secretion efficiency associated with the α-MF secretion signal leaders from K. pastoris and K. phaffii was notably lower in comparison to the α-MF secretion signal leader from S. cerevisiae. Further research indicated that N-linked glycosylation of the α-MF secretion signal leader enhanced the secretion of EGFP. Disruption of calnexin impaired the secretion of EGFP mediated by the N-linked glycosylated α-MF secretion signal leader, without affecting EGFP secretion mediated by the non-N-linked glycosylation α-MF secretion signal leader. The N-linked glycosylated of the α-MF secretion signal leader reduced the unfolded protein response (UPR) in the endoplasmic reticulum (ER). The enhancement of EGFP secretion by the N-linked glycosylated α-MF secretion signal leader might be achieved through the acceleration of proper folding of glycoproteins by the molecular chaperone calnexin. This study enhances the understanding of protein secretion in P. pastoris, specifically highlighting the influence of N-linked glycosylation on secretion efficiency, and could have implications for the production of recombinant proteins in bioengineering and biotechnological applications in P. pastoris.

CRISPR-Cas9 knockout screen informs efficient reduction of the Komagataella phaffii secretome

BACKGROUND

The yeast Komagataella phaffii is widely used for manufacturing recombinant proteins, but secreted titers of recombinant proteins could be improved by genetic engineering. In this study, we hypothesized that cellular resources could be redirected from production of endogenous proteins to production of recombinant proteins by deleting unneeded endogenous proteins. In non-model microorganisms such as K. phaffii, however, genetic engineering is limited by lack gene annotation and knowledge of gene essentiality.

RESULTS

We identified a set of endogenous secreted proteins in K. phaffii by mass spectrometry and signal peptide prediction. Our efforts to disrupt these genes were hindered by limited annotation of essential genes. To predict essential genes, therefore, we designed, transformed, and sequenced a pooled library of guide RNAs for CRISPR-Cas9-mediated knockout of all endogenous secreted proteins. We then used predicted gene essentiality to guide iterative disruptions of up to 11 non-essential genes. Engineered strains exhibited a ~20× increase in the production of human serum albumin and a twofold increase in the production of a monoclonal antibody.

CONCLUSIONS

We demonstrated that disruption of as few as six genes can increase production of recombinant proteins. Further reduction of the endogenous proteome of K. phaffii may further improve strain performance. The pooled library of secretome-targeted guides for CRISPR-Cas9 and knowledge of gene essentiality reported here will facilitate future efforts to engineer K. phaffii for production of other recombinant proteins and enzymes.

A Systematic Review of the Potential of Pichia pastoris (Komagataella phaffii) as an Alternative Host for Biologics Production

The methylotrophic yeast Pichia pastoris is garnering interest as a chassis cell factory for the manufacture of recombinant proteins because it effectively satisfies the requirements of both laboratory and industrial set up. The optimisation of P. pastoris cultivation is still necessary due to strain- and product-specific problems such as promoter strength, methanol utilisation type, and culturing conditions to realize the high yields of heterologous protein(s) of interest. Techniques integrating genetic and process engineering have been used to overcome these problems. Insight into the Pichia as an expression system utilizing MUT pathway and the development of methanol free systems are highlighted in this systematic review. Recent developments in the improved production of proteins in P. pastoris by (i) diverse genetic engineering such as codon optimization and gene dosage; (ii) cultivating tactics including co-expression of chaperones; (iii) advances in the use of the 2A peptide system, and (iv) CRISPR/Cas technologies are widely discussed. We believe that by combining these strategies, P. pastoris will become a formidable platform for the production of high value therapeutic proteins.

Genomic factors shape carbon and nitrogen metabolic niche breadth across Saccharomycotina yeasts

Organisms exhibit extensive variation in ecological niche breadth, from very narrow (specialists) to very broad (generalists). Two general paradigms have been proposed to explain this variation: (i) trade-offs between performance efficiency and breadth and (ii) the joint influence of extrinsic (environmental) and intrinsic (genomic) factors. We assembled genomic, metabolic, and ecological data from nearly all known species of the ancient fungal subphylum Saccharomycotina (1154 yeast strains from 1051 species), grown in 24 different environmental conditions, to examine niche breadth evolution. We found that large differences in the breadth of carbon utilization traits between yeasts stem from intrinsic differences in genes encoding specific metabolic pathways, but we found limited evidence for trade-offs. These comprehensive data argue that intrinsic factors shape niche breadth variation in microbes.

OPENPichia: licence-free Komagataella phaffii chassis strains and toolkit for protein expression

The industrial yeast Komagataella phaffii (formerly named Pichia pastoris) is commonly used to synthesize recombinant proteins, many of which are used as human therapeutics or in food. However, the basic strain, named NRRL Y-11430, from which all commercial hosts are derived, is not available without restrictions on its use. Comparative genome sequencing leaves little doubt that NRRL Y-11430 is derived from a K. phaffii type strain deposited in the UC Davis Phaff Yeast Strain Collection in 1954. We analysed four equivalent type strains in several culture collections and identified the NCYC 2543 strain, from which we started to develop an open-access Pichia chassis strain that anyone can use to produce recombinant proteins to industry standards. NRRL Y-11430 is readily transformable, which we found to be due to a HOC1 open-reading-frame truncation that alters cell-wall mannan. We introduced the HOC1 open-reading-frame truncation into NCYC 2543, which increased the transformability and improved secretion of some but not all of our tested proteins. We provide our genome-sequenced type strain, the hoc1tr derivative that we named OPENPichia as well as a synthetic, modular expression vector toolkit under liberal end-user distribution licences as an unencumbered OPENPichia resource for the microbial biotechnology community.

A supernumerary synthetic chromosome in Komagataella phaffii as a repository for extraneous genetic material

BACKGROUND

Komagataella phaffii (Pichia pastoris) is a methylotrophic commercially important non-conventional species of yeast that grows in a fermentor to exceptionally high densities on simple media and secretes recombinant proteins efficiently. Genetic engineering strategies are being explored in this organism to facilitate cost-effective biomanufacturing. Small, stable artificial chromosomes in K. phaffii could offer unique advantages by accommodating multiple integrations of extraneous genes and their promoters without accumulating perturbations of native chromosomes or exhausting the availability of selection markers.

RESULTS

Here, we describe a linear "nano"chromosome (of 15-25 kb) that, according to whole-genome sequencing, persists in K. phaffii over many generations with a copy number per cell of one, provided non-homologous end joining is compromised (by KU70-knockout). The nanochromosome includes a copy of the centromere from K. phaffii chromosome 3, a K. phaffii-derived autonomously replicating sequence on either side of the centromere, and a pair of K. phaffii-like telomeres. It contains, within its q arm, a landing zone in which genes of interest alternate with long (approx. 1-kb) non-coding DNA chosen to facilitate homologous recombination and serve as spacers. The landing zone can be extended along the nanochromosome, in an inch-worming mode of sequential gene integrations, accompanied by recycling of just two antibiotic-resistance markers. The nanochromosome was used to express PDI, a gene encoding protein disulfide isomerase. Co-expression with PDI allowed the production, from a genomically integrated gene, of secreted murine complement factor H, a plasma protein containing 40 disulfide bonds. As further proof-of-principle, we co-expressed, from a nanochromosome, both PDI and a gene for GFP-tagged human complement factor H under the control of PAOX1 and demonstrated that the secreted protein was active as a regulator of the complement system.

CONCLUSIONS

We have added K. phaffii to the list of organisms that can produce human proteins from genes carried on a stable, linear, artificial chromosome. We envisage using nanochromosomes as repositories for numerous extraneous genes, allowing intensive engineering of K. phaffii without compromising its genome or weakening the resulting strain.

Genetic tools for metabolic engineering of Pichia pastoris

The methylotrophic yeast Pichia pastoris (also known as Komagataella phaffii) is widely used as a yeast cell factory for producing heterologous proteins. Recently, it has gained attention for its potential in producing chemicals from inexpensive feedstocks, which requires efficient genetic engineering platforms. This review provides an overview of the current advances in developing genetic tools for metabolic engineering of P. pastoris. The topics cover promoters, terminators, plasmids, genome integration sites, and genetic editing systems, with a special focus on the development of CRISPR/Cas systems and their comparison to other genome editing tools. Additionally, this review highlights the prospects of multiplex genome integration, fine-tuning gene expression, and single-base editing systems. Overall, the aim of this review is to provide valuable insights into current genetic engineering and discuss potential directions for future efforts in developing efficient genetic tools in P. pastoris.

Effect of Biotin Starvation on Gene Expression in Komagataella phaffii Cells

Methylotrophic yeast Komagataella phaffii is widely used in biotechnology for recombinant protein production. Due to the practical significance of these yeasts, it is extremely important to properly select cultivation conditions and optimize the media composition. In this study the effect of biotin starvation on gene expression in K. phaffii at transcriptome level was investigated. It was demonstrated, that the response of K. phaffii cell to biotin deficiency strongly depends on the carbon source in the medium. In the media containing glycerol, biotin deficiency led to activation of the genes involved in biotin metabolism, glyoxylate cycle, and synthesis of acetyl-CoA in cytoplasm, as well as repression of the genes, involved in lipo- and gluconeogenesis. In the methanol-containing media, biotin deficiency primarily led to repression of the genes, involved in protein synthesis, and activation of cell response to oxidative stress.

Engineering Komagataella phaffii to biosynthesize cordycepin from methanol which drives global metabolic alterations at the transcription level

Cordycepin has the potential to be an alternative to the disputed herbicide glyphosate. However, current laborious and time-consuming production strategies at low yields based on Cordyceps militaris lead to extremely high cost and restrict its application in the field of agriculture. In this study, Komagataella phaffii (syn. Pichia pastoris) was engineered to biosynthesize cordycepin from methanol, which could be converted from CO2. Combined with fermentation optimization, cordycepin content in broth reached as high as 2.68 ± 0.04 g/L within 168 h, around 15.95 mg/(L·h) in productivity. Additionally, a deaminated product of cordycepin was identified at neutral or weakly alkaline starting pH during fermentation. Transcriptome analysis found the yeast producing cordycepin was experiencing severe inhibition in methanol assimilation and peroxisome biogenesis, responsible for delayed growth and decreased carbon flux to pentose phosphate pathway (PPP) which led to lack of precursor supply. Amino acid interconversion and disruption in RNA metabolism were also due to accumulation of cordycepin. The study provided a unique platform for the manufacture of cordycepin based on the emerging non-conventional yeast and gave practical strategies for further optimization of the microbial cell factory.

Advances in Komagataella phaffii Engineering for the Production of Renewable Chemicals and Proteins

The need for a more sustainable society has prompted the development of bio-based processes to produce fuels, chemicals, and materials in substitution for fossil-based ones. In this context, microorganisms have been employed to convert renewable carbon sources into various products. The methylotrophic yeast Komagataella phaffii has been extensively used in the production of heterologous proteins. More recently, it has been explored as a host organism to produce various chemicals through new metabolic engineering and synthetic biology tools. This review first summarizes Komagataella taxonomy and diversity and then highlights the recent approaches in cell engineering to produce renewable chemicals and proteins. Finally, strategies to optimize and develop new fermentative processes using K. phaffii as a cell factory are presented and discussed. The yeast K. phaffii shows an outstanding performance for renewable chemicals and protein production due to its ability to metabolize different carbon sources and the availability of engineering tools. Indeed, it has been employed in producing alcohols, carboxylic acids, proteins, and other compounds using different carbon sources, including glycerol, glucose, xylose, methanol, and even CO2.

Isolation and evaluation of strong endogenous promoters for the heterologous expression of proteins in Pichia pastoris

BACKGROUND

The heterologous expression of biosynthetic pathway genes for pharmaceutical or fine chemical production usually requires to express more than one gene in the host cells. In eukaryotes, the pathway flux is typically balanced by controlling the transcript levels of the genes involved. It is difficult to balance the stoichiometric fine-tuning of the reaction steps of the pathway by acting on one or two promoters. Furthermore, the promoter used should not be identical to avoid loss of inserted genes by recombination or dilute its transcription factors.

RESULTS

Based on RNA-seq data, 18 candidate genes with the highest transcription levels at three carbon sources (glucose, glycerol and methanol) were selected and their promoter regions were isolated from GS115 genome. The performance of these promoters on the level of protein production was evaluated using LacZ and EGFP genes as the reporters, respectively. These isolated promoters all exhibited activity to express LacZ gene. Using LacZ as a reporter, of the 18 promoter candidates, 9 promoters showed higher expression levels for the reporter compare to pGAP, a strong promoter widely used for constitutive expression of heterologous proteins in Pichia pastoris. These promoters with high expression levels were further employed to evaluate secreted expression using EGFP as a reporter. 6 promoters exhibited stronger protein expression compare to pGAP. Interestingly, the protein expression driven by pFDH1 was slightly higher than that of commonly used pAOX1 at methanol, and methanol-induced expression of pFDH1 was not repressed by glycerol.

CONCLUSION

The various promoters identified in this study could be used for heterologous expression of biosynthetic pathway genes for pharmaceutical or fine chemical production. the methanol-induced pFDH1 that is not repressed by glycerol is an attractive alternative to pAOX1 and may provide a novel way to produce heterologous proteins in Pichia pastoris.

Comparative proteomics analysis of Pichia pastoris cultivating in glucose and methanol

The methylotrophic yeast Pichia pastoris (syn. Komagataella phaffii) has been extensively engineered for protein production, and is attracting attention as a chassis cell for methanol biotransformation toward production of small molecules. However, the relatively unclear methanol metabolism hampers the metabolic rewiring to improve the biosynthetic efficiency. We here performed a label-free quantitative proteomic analysis of Pichia pastoris when cultivated in minimal media containing methanol and glucose, respectively. There were 243, 158 up-regulated proteins and 244, 304 down-regulated proteins in log and stationary phase, respectively, when cultivated in methanol medium compared with that of glucose medium. Peroxisome enrichment further improved the characterization of more differentially expressed proteins (481 proteins in log phase and 524 proteins in stationary phase). We demonstrated the transaldolase isoenzyme (Tal2, Protein ID: C4R244) was highly up-regulated in methanol medium cultivation, which plays an important role in methanol utilization. Our work provides important information for understanding methanol metabolism in methyltrophic yeast and will help to engineer methanol biotransformation in P. pastoris.

Genotypic and phenotypic diversity among Komagataella species reveals a hidden pathway for xylose utilization

BACKGROUND

The yeast genus Komagataella currently consists of seven methylotrophic species isolated from tree environments. Well-characterized strains of K. phaffii and K. pastoris are important hosts for biotechnological applications, but the potential of other species from the genus remains largely unexplored. In this study, we characterized 25 natural isolates from all seven described Komagataella species to identify interesting traits and provide a comprehensive overview of the genotypic and phenotypic diversity available within this genus.

RESULTS

Growth tests on different carbon sources and in the presence of stressors at two different temperatures allowed us to identify strains with differences in tolerance to high pH, high temperature, and growth on xylose. As Komagataella species are generally not considered xylose-utilizing yeasts, xylose assimilation was characterized in detail. Growth assays, enzyme activity measurements and 13C labeling confirmed the ability of K. phaffii to utilize D-xylose via the oxidoreductase pathway. In addition, we performed long-read whole-genome sequencing to generate genome assemblies of all Komagataella species type strains and additional K. phaffii and K. pastoris isolates for comparative analysis. All sequenced genomes have a similar size and share 83-99% average sequence identity. Genome structure analysis showed that K. pastoris and K. ulmi share the same rearrangements in difference to K. phaffii, while the genome structure of K. kurtzmanii is similar to K. phaffii. The genomes of the other, more distant species showed a larger number of structural differences. Moreover, we used the newly assembled genomes to identify putative orthologs of important xylose-related genes in the different Komagataella species.

CONCLUSIONS

By characterizing the phenotypes of 25 natural Komagataella isolates, we could identify strains with improved growth on different relevant carbon sources and stress conditions. Our data on the phenotypic and genotypic diversity will provide the basis for the use of so-far neglected Komagataella strains with interesting characteristics and the elucidation of the genetic determinants of improved growth and stress tolerance for targeted strain improvement.

Burden Imposed by Heterologous Protein Production in Two Major Industrial Yeast Cell Factories: Identifying Sources and Mitigation Strategies

Production of heterologous proteins, especially biopharmaceuticals and industrial enzymes, in living cell factories consumes cellular resources. Such resources are reallocated from normal cellular processes toward production of the heterologous protein that is often of no benefit to the host cell. This competition for resources is a burden to host cells, has a negative impact on cell fitness, and may consequently trigger stress responses. Importantly, this often causes a reduction in final protein titers. Engineering strategies to generate more burden resilient production strains offer sustainable opportunities to increase production and profitability for this growing billion-dollar global industry. We review recently reported impacts of burden derived from resource competition in two commonly used protein-producing yeast cell factories: Saccharomyces cerevisiae and Komagataella phaffii (syn. Pichia pastoris). We dissect possible sources of burden in these organisms, from aspects related to genetic engineering to protein translation and export of soluble protein. We also summarize advances as well as challenges for cell factory design to mitigate burden and increase overall heterologous protein production from metabolic engineering, systems biology, and synthetic biology perspectives. Lastly, future profiling and engineering strategies are highlighted that may lead to constructing robust burden-resistant cell factories. This includes incorporation of systems-level data into mathematical models for rational design and engineering dynamical regulation circuits in production strains.

Strains and Molecular Tools for Recombinant Protein Production in Pichia pastoris

Within the last two decades, the methylotrophic yeast Pichia pastoris (Komagataella phaffii) has become an important alternative to E. coli or mammalian cell lines for the production of recombinant proteins. Easy handling, strong promoters, and high cell density cultivations as well as the capability of posttranslational modifications are some of the major benefits of this yeast. The high secretion capacity and low level of endogenously secreted proteins further promoted the rapid development of a versatile Pichia pastoris toolbox. This chapter reviews common and new "Pichia tools" and their specific features. Special focus is given to expression strains, such as different methanol utilization, protease-deficient or glycoengineered strains, combined with application highlights. Different promoters and signal sequences are also discussed.

SuMO-Fil: Supervised multi-omic filtering prior to performing network analysis

Multi-omic analyses that integrate many high-dimensional datasets often present significant deficiencies in statistical power and require time consuming computations to execute the analytical methods. We present SuMO-Fil to remedy against these issues which is a pre-processing method for Supervised Multi-Omic Filtering that removes variables or features considered to be irrelevant noise. SuMO-Fil is intended to be performed prior to downstream analyses that detect supervised gene networks in sparse settings. We accomplish this by implementing variable filters based on low similarity across the datasets in conjunction with low similarity with the outcome. This approach can improve accuracy, as well as reduce run times for a variety of computationally expensive downstream analyses. This method has applications in a setting where the downstream analysis may include sparse canonical correlation analysis. Filtering methods specifically for cluster and network analysis are introduced and compared by simulating modular networks with known statistical properties. The SuMO-Fil method performs favorably by eliminating non-network features while maintaining important biological signal under a variety of different signal settings as compared to popular filtering techniques based on low means or low variances. We show that the speed and accuracy of methods such as supervised sparse canonical correlation are increased after using SuMO-Fil, thus greatly improving the scalability of these approaches.

Engineered yeast genomes accurately assembled from pure and mixed samples

Yeast whole genome sequencing (WGS) lacks end-to-end workflows that identify genetic engineering. Here we present Prymetime, a tool that assembles yeast plasmids and chromosomes and annotates genetic engineering sequences. It is a hybrid workflow-it uses short and long reads as inputs to perform separate linear and circular assembly steps. This structure is necessary to accurately resolve genetic engineering sequences in plasmids and the genome. We show this by assembling diverse engineered yeasts, in some cases revealing unintended deletions and integrations. Furthermore, the resulting whole genomes are high quality, although the underlying assembly software does not consistently resolve highly repetitive genome features. Finally, we assemble plasmids and genome integrations from metagenomic sequencing, even with 1 engineered cell in 1000. This work is a blueprint for building WGS workflows and establishes WGS-based identification of yeast genetic engineering.

Translational landscape and protein biogenesis demands of the early secretory pathway in Komagataella phaffii

Background

Eukaryotes use distinct networks of biogenesis factors to synthesize, fold, monitor, traffic, and secrete proteins. During heterologous expression, saturation of any of these networks may bottleneck titer and yield. To understand the flux through various routes into the early secretory pathway, we quantified the global and membrane-associated translatomes of Komagataella phaffii.

Results

By coupling Ribo-seq with long-read mRNA sequencing, we generated a new annotation of protein-encoding genes. By using Ribo-seq with subcellular fractionation, we quantified demands on co- and posttranslational translocation pathways. During exponential growth in rich media, protein components of the cell-wall represent the greatest number of nascent chains entering the ER. Transcripts encoding the transmembrane protein PMA1 sequester more ribosomes at the ER membrane than any others. Comparison to Saccharomyces cerevisiae reveals conservation in the resources allocated by gene ontology, but variation in the diversity of gene products entering the secretory pathway.

Conclusion

A subset of host proteins, particularly cell-wall components, impose the greatest biosynthetic demands in the early secretory pathway. These proteins are potential targets in strain engineering aimed at alleviating bottlenecks during heterologous protein production.

Komagataella phaffii as Emerging Model Organism in Fundamental Research

Komagataella phaffii (Pichia pastoris) is one of the most extensively applied yeast species in pharmaceutical and biotechnological industries, and, therefore, also called the biotech yeast. However, thanks to more advanced strain engineering techniques, it recently started to gain attention as model organism in fundamental research. So far, the most studied model yeast is its distant cousin, Saccharomyces cerevisiae. While these data are of great importance, they limit our knowledge to one organism only. Since the divergence of the two species 250 million years ago, K. phaffii appears to have evolved less rapidly than S. cerevisiae, which is why it remains more characteristic of the common ancient yeast ancestors and shares more features with metazoan cells. This makes K. phaffii a valuable model organism for research on eukaryotic molecular cell biology, a potential we are only beginning to fully exploit. As methylotrophic yeast, K. phaffii has the intriguing property of being able to efficiently assimilate methanol as a sole source of carbon and energy. Therefore, major efforts have been made using K. phaffii as model organism to study methanol assimilation, peroxisome biogenesis and pexophagy. Other research topics covered in this review range from yeast genetics including mating and sporulation behavior to other cellular processes such as protein secretion, lipid biosynthesis and cell wall biogenesis. In this review article, we compare data obtained from K. phaffii with S. cerevisiae and other yeasts whenever relevant, elucidate major differences, and, most importantly, highlight the big potential of using K. phaffii in fundamental research.

Towards Genomic Criteria for Delineating Fungal Species

The discussion of fungal species delineation has yet to reach a consensus, despite the advancements in technology, which helped modernise traditional approaches. In particular, the phylogenetic species concept was one of the tools that has been used with considerable success across the fungal kingdom. The fast rise of fungal genomics provides an unprecedented opportunity to expand measuring the relatedness of fungal strains to the level of whole genomes. However, the use of genomic information in taxonomy has only just begun, and few methodological guidelines have been suggested so far. Here, a simple approach of computationally measuring genomic distances and their use as a standard for species delineation is investigated. A fixed threshold genomic distance calculated by the quick and easy-to-use tools Mash and Dashing proved to be an unexpectedly widely applicable and robust criterion for determining whether two genomes belong to the same or to different species. The accuracy of species delineation in an uncurated dataset of GenBank fungal genomes was close to 90%—and exceeded 90% with minimal curation. As expected, the discriminative power of this approach was lower at higher taxonomic ranks, but still significantly larger than zero. Simple instructions for calculation of a genomic distance between two genomes and species similarity thresholds at different k-mer sizes are suggested. The calculation of genomic distance is identified as a powerful approach for delineating fungal species and is proposed—not as the only criterion—but as an additional tool in the versatile toolbox of fungal taxonomy.

Identifying Improved Sites for Heterologous Gene Integration Using ATAC-seq

Constructing efficient cellular factories often requires integration of heterologous pathways for synthesis of novel compounds and improved cellular productivity. Few genomic sites are routinely used, however, for efficient integration and expression of heterologous genes, especially in nonmodel hosts. Here, a data-guided framework for informing suitable integration sites for heterologous genes based on ATAC-seq was developed in the nonmodel yeast Komagataella phaffii. Single-copy GFP constructs were integrated using CRISPR/Cas9 into 38 intergenic regions (IGRs) to evaluate the effects of IGR size, intensity of ATAC-seq peaks, and orientation and expression of adjacent genes. Only the intensity of accessibility peaks was observed to have a significant effect, with higher expression observed from IGRs with low- to moderate-intensity peaks than from high-intensity peaks. This effect diminished for tandem, multicopy integrations, suggesting that the additional copies of exogenous sequence buffered the transcriptional unit of the transgene against effects from endogenous sequence context. The approach developed from these results should provide a basis for nominating suitable IGRs in other eukaryotic hosts from an annotated genome and ATAC-seq data.

Microbe Profile: Komagataella phaffii: a methanol devouring biotech yeast formerly known as Pichia pastoris

Methylotrophic yeasts of the genus Komagataella are abundantly found in tree exudates. Their ability to utilize methanol as carbon and energy source relies on an assimilation pathway localized in largely expanded peroxisomes, and a cytosolic methanol dissimilation pathway. Other substrates like glucose or glycerol are readily utilized as well. Komagataella yeasts usually grow as haploid cells and are secondary homothallic as they can switch mating type. Upon mating diploid cells sporulate readily, forming asci with four haploid spores. Their ability to secrete high amounts of heterologous proteins made them interesting for biotechnology, which expands today also to other products of primary and secondary metabolism.

The Methylotroph Gene Order Browser (MGOB) reveals conserved synteny and ancestral centromere locations in the yeast family Pichiaceae

The yeast family Pichiaceae, also known as the 'methylotrophs clade', is a relatively little studied group of yeasts despite its economic and clinical relevance. To explore the genome evolution and synteny relationships within this family, we developed the Methylotroph Gene Order Browser (MGOB, http://mgob.ucd.ie) similar to our previous gene order browsers for other yeast families. The dataset contains genome sequences from nine Pichiaceae species, including our recent reference sequence of Pichia kudriavzevii. As an example, we demonstrate the conservation of synteny around the MOX1 locus among species both containing and lacking the MOX1 gene for methanol assimilation. We found ancient clusters of genes that are conserved as adjacent between Pichiaceae and Saccharomycetaceae. Surprisingly, we found evidence that the locations of some centromeres have been conserved among Pichiaceae species, and between Pichiaceae and Saccharomycetaceae, even though the centromeres fall into different structural categories-point centromeres, inverted repeats and retrotransposon cluster centromeres.

Comparative genome-scale analysis of Pichia pastoris variants informs selection of an optimal base strain

Komagataella phaffii, also known as Pichia pastoris, is a common host for the production of biologics and enzymes, due to fast growth, high productivity, and advancements in host engineering. Several K. phaffii variants are commonly used as interchangeable base strains, which confounds efforts to improve this host. In this study, genomic and transcriptomic analyses of Y-11430 (CBS7435), GS115, X-33, and eight other variants enabled a comparative assessment of the relative fitness of these hosts for recombinant protein expression. Cell wall integrity explained the majority of the variation among strains, impacting transformation efficiency, growth, methanol metabolism, and secretion of heterologous proteins. Y-11430 exhibited the highest activity of genes involved in methanol utilization, up to two-fold higher transcription of heterologous genes, and robust growth. With a more permeable cell wall, X-33 displayed a six-fold higher transformation efficiency and up to 1.2-fold higher titers than Y-11430. X-33 also shared nearly all mutations, and a defective variant of HIS4, with GS115, precluding robust growth. Transferring two beneficial mutations identified in X-33 into Y-11430 resulted in an optimized base strain that provided up to four-fold higher transformation efficiency and three-fold higher protein titers, while retaining robust growth. The approach employed here to assess unique banked variants in a species and then transfer key beneficial variants into a base strain should also facilitate rational assessment of a broad set of other recombinant hosts.

Host-Informed Expression of CRISPR Guide RNA for Genomic Engineering in Komagataella phaffii

There is growing interest in the use of nonmodel microorganisms as hosts for biopharmaceutical manufacturing. These hosts require genomic engineering to meet clinically relevant product qualities and titers, but the adaptation of tools for editing genomes, such as CRISPR-Cas9, has been slow for poorly characterized hosts. Specifically, a lack of biochemical characterization of RNA polymerase III transcription has hindered reliable expression of guide RNAs in new hosts. Here, we present a sequencing-based strategy for the design of host-specific cassettes for modular, reliable, expression of guide RNAs. Using this strategy, we achieved up to 95% gene editing efficiency in the methylotrophic yeast Komagataella phaffii. We applied this approach for the rapid, multiplexed engineering of a complex phenotype, achieving humanized product glycosylation in two sequential steps of engineering. Reliable extension of simple gene editing tools to nonmodel manufacturing hosts will enable rapid engineering of manufacturing strains tuned for specific product profiles and potentially decrease the costs and timelines for process development.

Genomic diversity and meiotic recombination among isolates of the biotech yeast Komagataella phaffii (Pichia pastoris)

BACKGROUND

Komagataella phaffii is a yeast widely used in the pharmaceutical and biotechnology industries, and is one of the two species that were previously called Pichia pastoris. However, almost all laboratory work on K. phaffii has utilized strains derived from a single natural isolate, CBS7435. There is little information about the sequence diversity of K. phaffii or the genetic properties of this species.

RESULTS

We sequenced the genomes of all the known isolates of K. phaffii. We made a genetic cross between derivatives of two isolates that differ at 44,000 single nucleotide polymorphism sites, and used this cross to analyze the rate and landscape of meiotic recombination. We conducted tetrad analysis by making use of the property that K. phaffii haploids do not mate in rich media, which enabled us to isolate and sequence the four types of haploid cell that are present in the colony that forms when a tetra-type ascus germinates.

CONCLUSIONS

We found that only four distinct natural isolates of K. phaffii exist in public yeast culture collections. The meiotic recombination rate in K. phaffii is approximately 3.5 times lower than in Saccharomyces cerevisiae, with an average of 25 crossovers per meiosis. Recombination is suppressed, and genetic diversity among natural isolates is low, in a region around centromeres that is much larger than the centromeres themselves. Our work lays a foundation for future quantitative trait locus analysis in K. phaffii.

Annotation of a hypothetical protein coding gene PAS_chr2-2_0152 containing Lysine Methyl transferase SMYD domain from Komagataella phaffii GS115

The methylotrophic yeast Komagataella phaffii is an industrial workhorse yeast species that has been widely used in biotechnology industries for recombinant protein production. Genome sequencing of this yeast in 2009 have enabled scientists to assign and characterize functions to most of its proteins while few hypothetical proteins remain uncharacterized. Therefore, it is of interest to characterize the hypothetical protein coding gene PAS_chr2-2_0152 as SET containing the ZNF-MYND (SMYD) domain. They share a homology with other methylotrophic and non-methylotrophic yeast species together with known SMYD proteins of Homo sapiens, with conserved distinctive SMYD domain patterns. A homology model is developed using the crystal structure of human histone-lysine methyl transferase smyd3 as template. These data points to that the hypothetical protein is a potential histones and non-histone lysine methyl transferase regulating cell cycle, chromatin remodeling, DNA damage response, homologous recombination and transcription in Komagataella phaffii. Data also suggests the evolutionary syntenic conservation of DNA damage regulator (RFX) and lysine methyl transferase (SMYD) genes in some yeast lineages, pointing to a conserved role requiring further confirmation.

Transcriptomic Analysis of the Influence of Methanol Assimilation on the Gene Expression in the Recombinant Pichia pastoris Producing Hirudin Variant 3

Hirudin and its variants, as strong inhibitors against thrombin, are present in the saliva of leeches and are recognized as potent anticoagulants. However, their yield is far from the clinical requirement up to now. In this study, the production of hirudin variant 3 (HV3) was successfully realized by cultivating the recombinant Pichia pastoris GS115/pPIC9K-hv3 under the regulation of the promoter of AOX1 encoding alcohol oxidase (AOX). The antithrombin activity in the fermentation broth reached the maximum value of 5000 ATU/mL. To explore an effective strategy for improving HV3 production in the future, we investigated the influence of methanol assimilation on the general gene expression in this recombinant by transcriptomic study. The results showed that methanol was partially oxidized into CO₂, and the rest was converted into glycerone-P which subsequently entered into central carbon metabolism, energy metabolism, and amino acid biosynthesis. However, the later metabolic processes were almost all down-regulated. Therefore, we propose that the up-regulated central carbon metabolism, energy, and amino acid metabolism should be beneficial for methanol assimilation, which would accordingly improve the production of HV3.

Deletion of DNA ligase IV homolog confers higher gene targeting efficiency on homologous recombination in Komagataella phaffii

The non-conventional yeast Komagataella phaffii, formerly Pichia pastoris, is a popular host for recombinant protein production. The relatively lower gene targeting efficiency observed in this species occurs due to high levels of non-homologous recombination activity. In the current study, we explored the function of the K. phaffii homolog of DNA ligase IV (Dnl4p) by creating a DNL4-disrupted strain. To assess the roles of non-homologous end joining (NHEJ)-related proteins in this species, strains deleted for either or both genes encoding Dnl4p or the telomeric Ku complex subunit (Ku70p) were generated. These deletions were constructed by either of two distinct marker-recycling methods (yielding either a seamless gene deletion or a Cre-loxP-mediated gene deletion). The resulting dnl4- and/or ku70-deleted K. phaffii strains were used to evaluate gene targeting efficiency in gene knock-out and gene knock-in experiments. The Dnl4p-defective strain showed improved gene targeting efficiency for homologous recombination compared to the wild-type and Ku70p-deffective strains. The dnl4 ku70 double knock-out strain exhibited a further improvement in gene targeting efficiency. Thus, the K. phaffii dnl4 and dnl4 ku70 deletion strains are expected to serve as useful platforms for functional analysis and strain development in this species.

Coverage-Versus-Length Plots, a Simple Quality Control Step for de Novo Yeast Genome Sequence Assemblies

Illumina sequencing has revolutionized yeast genomics, with prices for commercial draft genome sequencing now below $200. The popular SPAdes assembler makes it simple to generate a de novo genome assembly for any yeast species. However, whereas making genome assemblies has become routine, understanding what they contain is still challenging. Here, we show how graphing the information that SPAdes provides about the length and coverage of each scaffold can be used to investigate the nature of an assembly, and to diagnose possible problems. Scaffolds derived from mitochondrial DNA, ribosomal DNA, and yeast plasmids can be identified by their high coverage. Contaminating data, such as cross-contamination from other samples in a multiplex sequencing run, can be identified by its low coverage. Scaffolds derived from the bacteriophage PhiX174 and Lambda DNAs that are frequently used as molecular standards in Illumina protocols can also be detected. Assemblies of yeast genomes with high heterozygosity, such as interspecies hybrids, often contain two types of scaffold: regions of the genome where the two alleles assembled into two separate scaffolds and each has a coverage level C, and regions where the two alleles co-assembled (collapsed) into a single scaffold that has a coverage level 2C. Visualizing the data with Coverage-vs.-Length (CVL) plots, which can be done using Microsoft Excel or Google Sheets, provides a simple method to understand the structure of a genome assembly and detect aberrant scaffolds or contigs. We provide a Python script that allows assemblies to be filtered to remove contaminants identified in CVL plots.

Challenging the workhorse: Comparative analysis of eukaryotic micro-organisms for expressing monoclonal antibodies

For commercial protein therapeutics, Chinese hamster ovary (CHO) cells have an established history of safety, proven capability to express a wide range of therapeutic proteins and high volumetric productivities. Expanding global markets for therapeutic proteins and increasing concerns for broadened access of these medicines has catalyzed consideration of alternative approaches to this platform. Reaching these objectives likely will require an order of magnitude increase in volumetric productivity and a corresponding reduction in the costs of manufacture. For CHO-based manufacturing, achieving this combination of targeted improvements presents challenges. Based on a holistic analysis, the choice of host cells was identified as the single most influential factor for both increasing productivity and decreasing costs. Here we evaluated eight wild-type eukaryotic micro-organisms with prior histories of recombinant protein expression. The evaluation focused on assessing the potential of each host, and their corresponding phyla, with respect to key attributes relevant for manufacturing, namely (a) growth rates in industry-relevant media, (b) adaptability to modern techniques for genome editing, and (c) initial characterization of product quality. These characterizations showed that multiple organisms may be suitable for production with appropriate engineering and development and highlighted that yeast in general present advantages for rapid genome engineering and development cycles.

Relative Abundance of Transcripts ( RATs): Identifying differential isoform abundance from RNA-seq

The biological importance of changes in RNA expression is reflected by the wide variety of tools available to characterise these changes from RNA-seq data. Several tools exist for detecting differential transcript isoform usage (DTU) from aligned or assembled RNA-seq data, but few exist for DTU detection from alignment-free RNA-seq quantifications. We present the RATs, an R package that identifies DTU transcriptome-wide directly from transcript abundance estimates. RATs is unique in applying bootstrapping to estimate the reliability of detected DTU events and shows good performance at all replication levels (median false positive fraction < 0.05). We compare RATs to two existing DTU tools, DRIM-Seq & SUPPA2, using two publicly available simulated RNA-seq datasets and a published human RNA-seq dataset, in which 248 genes have been previously identified as displaying significant DTU. RATs with default threshold values on the simulated Human data has a sensitivity of 0.55, a Matthews correlation coefficient of 0.71 and a false discovery rate (FDR) of 0.04, outperforming both other tools. Applying the same thresholds for SUPPA2 results in a higher sensitivity (0.61) but poorer FDR performance (0.33). RATs and DRIM-seq use different methods for measuring DTU effect-sizes complicating the comparison of results between these tools, however, for a likelihood-ratio threshold of 30, DRIM-Seq has similar FDR performance to RATs (0.06), but worse sensitivity (0.47). These differences persist for the simulated drosophila dataset. On the published human RNA-seq dataset the greatest agreement between the tools tested is 53%, observed between RATs and SUPPA2. The bootstrapping quality filter in RATs is responsible for removing the majority of DTU events called by SUPPA2 that are not reported by RATs. All methods, including the previously published qRT-PCR of three of the 248 detected DTU events, were found to be sensitive to annotation differences between Ensembl v60 and v87.

Established and Upcoming Yeast Expression Systems

Yeast was the first microorganism used by mankind for biotransformation of feedstock that laid the foundations of industrial biotechnology. Long historical use, vast amount of data, and experience paved the way for Saccharomyces cerevisiae as a first yeast cell factory, and still it is an important expression platform as being the production host for several large volume products. Continuing special needs of each targeted product and different requirements of bioprocess operations have led to identification of different yeast expression systems. Modern bioprocess engineering and advances in omics technology, i.e., genomics, transcriptomics, proteomics, secretomics, and interactomics, allow the design of novel genetic tools with fine-tuned characteristics to be used for research and industrial applications. This chapter focuses on established and upcoming yeast expression platforms that have exceptional characteristics, such as the ability to utilize a broad range of carbon sources or remarkable resistance to various stress conditions. Besides the conventional yeast S. cerevisiae, established yeast expression systems including the methylotrophic yeasts Pichia pastoris and Hansenula polymorpha, the dimorphic yeasts Arxula adeninivorans and Yarrowia lipolytica, the lactose-utilizing yeast Kluyveromyces lactis, the fission yeast Schizosaccharomyces pombe, and upcoming yeast platforms, namely, Kluyveromyces marxianus, Candida utilis, and Zygosaccharomyces bailii, are compiled with special emphasis on their genetic toolbox for recombinant protein production.

Neutron scattering in the biological sciences: progress and prospects

The scattering of neutrons can be used to provide information on the structure and dynamics of biological systems on multiple length and time scales. Pursuant to a National Science Foundation-funded workshop in February 2018, recent developments in this field are reviewed here, as well as future prospects that can be expected given recent advances in sources, instrumentation and computational power and methods. Crystallography, solution scattering, dynamics, membranes, labeling and imaging are examined. For the extraction of maximum information, the incorporation of judicious specific deuterium labeling, the integration of several types of experiment, and interpretation using high-performance computer simulation models are often found to be particularly powerful.

Structural and Functional Characterization of PA14/Flo5-Like Adhesins From Komagataella pastoris

Cell–cell and cell-substrate based adhesion of yeasts are major determinants of their adoption of different life styles. Genome-mining of ascomycetous GPI-anchored cell wall proteins with lectin-like PA14 domains identified a unique class of putative adhesins in the clade of methylotrophic Komagataella yeasts, many of which are known to colonize plants and insects involving yet unknown adhesion mechanisms. Here, we report the functional and structural analysis of two of its members: KpFlo1 (=Cea1), that is highly specific for terminal N-acetylglucosamine moieties, and KpFlo2, which represents an orphan lectin with intact binding site but unknown specificity. Crystal structures of the Cea1 adhesion domain complexed to N-acetylglucosamine and N,N′-diacetylchitobiose reveal a Ca²⁺-dependent binding mode that differs from other members of the PA14/Flo5 adhesin family. Heterologous expression of Cea1A in Saccharomyces cerevisiae promotes cellular adhesion to non-reducing ends of non-crystalline chitin. Overall, our data suggest that high-affinity recognition of β-GlcNAc-capped glycans by Cea1 enable Komagataella species to interact with surface cues present in fungi and insects.

On-demand manufacturing of clinical-quality biopharmaceuticals

Conventional manufacturing of protein biopharmaceuticals in centralized, large-scale, single-product facilities is not well-suited to the agile production of drugs for small patient populations or individuals. Previous solutions for small-scale manufacturing are limited in both process reproducibility and product quality, owing to their complicated means of protein expression and purification. We describe an automated, benchtop, multiproduct manufacturing system, called Integrated Scalable Cyto-Technology (InSCyT), for the end-to-end production of hundreds to thousands of doses of clinical-quality protein biologics in about 3 d. Unlike previous systems, InSCyT includes fully integrated modules for sustained production, efficient purification without the use of affinity tags, and formulation to a final dosage form of recombinant biopharmaceuticals. We demonstrate that InSCyT can accelerate process development from sequence to purified drug in 12 weeks. We used integrated design to produce human growth hormone, interferon α-2b and granulocyte colony-stimulating factor with highly similar processes on this system and show that their purity and potency are comparable to those of marketed reference products.

Identification and characterization of the Komagataella phaffii mating pheromone genes

The methylotrophic yeast Komagataella phaffii (Pichia pastoris) is a haploid yeast that is able to form diploid cells by mating once nitrogen becomes limiting. Activation of the mating response requires the secretion of a- and α-factor pheromones, which bind to G-protein coupled receptors on cells of opposite mating type. In K. phaffii, the genes coding for the α-factor (MFα), the pheromone surface receptors and the conserved a-factor biogenesis pathway have been annotated previously. Initial homology-based search failed to identify potential a-factor genes (MFA). By using transcriptome data of heterothallic strains under mating conditions, we found two K. phaffiia-factor genes. Deletion of both MFA genes prevented mating of a-type cells. MFA single mutants were still able to mate and activate the mating response pathway in α-type cells. A reporter assay was used to confirm the biological activity of synthetic a- and α-factor peptides. The identification of the a-factor genes enabled the first characterization of the role and regulation of the mating pheromone genes and the response of K. phaffii to synthetic pheromones and will help to gain a better understanding of the mating behavior of K. phaffii.

Metabolic engineering of Pichia pastoris

Besides its use for efficient production of recombinant proteins the methylotrophic yeast Pichia pastoris (syn. Komagataella spp.) has been increasingly employed as a platform to produce metabolites of varying origin. We summarize here the impressive methodological developments of the last years to model and analyze the metabolism of P. pastoris, and to engineer its genome and metabolic pathways. Efficient methods to insert, modify or delete genes via homologous recombination and CRISPR/Cas9, supported by modular cloning techniques, have been reported. An outstanding early example of metabolic engineering in P. pastoris was the humanization of protein glycosylation. More recently the cell metabolism was engineered also to enhance the productivity of heterologous proteins. The last few years have seen an increased number of metabolic pathway design and engineering in P. pastoris, mainly towards the production of complex (secondary) metabolites. In this review, we discuss the potential role of P. pastoris as a platform for metabolic engineering, its strengths, and major requirements for future developments of chassis strains based on synthetic biology principles.

Genome-wide analysis of long non-coding RNAs in Pichia pastoris during stress by RNA sequencing

Long non-coding RNAs play significant roles in many biological processes. The roles of lncRNAs in Pichia pastoris remain unclear. In this work, we focused on the identification of lncRNAs in P. pastoris and exploration of their potential roles in stress response to PLA₂ overexpression and methanol induction. By strand specific RNA sequencing, 208 novel long non-coding RNAs were identified and analyzed. Bioinformatic analysis showed potential trans-target genes and cis-regulated genes of 39 differential lncRNAs. Functional annotation and sequence motif analysis indicated that lncRNAs participate in pathways related to methanol degradation and production of the recombinant protein. The differential expression of lncRNAs was validated by qRT-PCR. Lastly, the potential functions of three lncRNAs were evaluated by knockdown of their expression and analysis of the expression levels of target genes. Our study identifies novel lncRNAs in P. pastoris induced during use as a bioreactor, facilitating future functional research.

Effect of Plasmid Design and Type of Integration Event on Recombinant Protein Expression in Pichia pastoris

Pichia pastoris (syn. Komagataella phaffii) is one of the most common eukaryotic expression systems for heterologous protein production. Expression cassettes are typically integrated in the genome to obtain stable expression strains. In contrast to Saccharomyces cerevisiae, where short overhangs are sufficient to target highly specific integration, long overhangs are more efficient in P. pastoris and ectopic integration of foreign DNA can occur. Here, we aimed to elucidate the influence of ectopic integration by high-throughput screening of >700 transformants and whole-genome sequencing of 27 transformants. Different vector designs and linearization approaches were used to mimic the most common integration events targeted in P. pastoris Fluorescence of an enhanced green fluorescent protein (eGFP) reporter protein was highly uniform among transformants when the expression cassettes were correctly integrated in the targeted locus. Surprisingly, most nonspecifically integrated transformants showed highly uniform expression that was comparable to specific integration, suggesting that nonspecific integration does not necessarily influence expression. However, a few clones (<10%) harboring ectopically integrated cassettes showed a greater variation spanning a 25-fold range, surpassing specifically integrated reference strains up to 6-fold. High-expression strains showed a correlation between increased gene copy numbers and high reporter protein fluorescence levels. Our results suggest that for comparing expression levels between strains, the integration locus can be neglected as long as a sufficient numbers of transformed strains are compared. For expression optimization of highly expressible proteins, increasing copy number appears to be the dominant positive influence rather than the integration locus, genomic rearrangements, deletions, or single-nucleotide polymorphisms (SNPs).IMPORTANCE Yeasts are commonly used as biotechnological production hosts for proteins and metabolites. In the yeast Saccharomyces cerevisiae, expression cassettes carrying foreign genes integrate highly specifically at the targeted sites in the genome. In contrast, cassettes often integrate at random genomic positions in nonconventional yeasts, such as Pichia pastoris (syn. Komagataella phaffii). Hence, cells from the same transformation event often behave differently, with significant clonal variation necessitating the screening of large numbers of strains. The importance of this study is that we systematically investigated the influence of integration events in more than 700 strains. Our findings provide novel insight into clonal variation in P. pastoris and, thus, how to avoid pitfalls and obtain reliable results. The underlying mechanisms may also play a role in other yeasts and hence could be generally relevant for recombinant yeast protein production strains.

Effect of Pyruvate Decarboxylase Knockout on Product Distribution Using Pichia pastoris (Komagataella phaffii) Engineered for Lactic Acid Production

Lactic acid is the monomer unit of the bioplastic poly-lactic acid (PLA). One candidate organism for lactic acid production is Pichia pastoris, a yeast widely used for heterologous protein production. Nevertheless, this yeast has a poor fermentative capability that can be modulated by controlling oxygen levels. In a previous study, lactate dehydrogenase (LDH) activity was introduced into P. pastoris, enabling this yeast to produce lactic acid. The present study aimed to increase the flow of pyruvate towards the production of lactic acid in P. pastoris. To this end, a strain designated GLp was constructed by inserting the bovine lactic acid dehydrogenase gene (LDHb) concomitantly with the interruption of the gene encoding pyruvate decarboxylase (PDC). Aerobic fermentation, followed by micro-aerophilic culture two-phase fermentations, showed that the GLp strain achieved a lactic acid yield of 0.65 g/g. The distribution of fermentation products demonstrated that the acetate titer was reduced by 20% in the GLp strain with a concomitant increase in arabitol production: arabitol increased from 0.025 g/g to 0.174 g/g when compared to the GS115 strain. Taken together, the results show a significant potential for P. pastoris in producing lactic acid. Moreover, for the first time, physiological data regarding co-product formation have indicated the redox balance limitations of this yeast.

Creation of Stable Heterothallic Strains of Komagataella phaffii Enables Dissection of Mating Gene Regulation

The methylotrophic yeast Komagataella phaffii (Pichia pastoris) is homothallic and has been reported to switch mating type by an ancient inversion mechanism. Two mating-type (MAT) loci include homologs of the MATa and MATα transcription factor genes, with the expression from one locus downregulated by telomere position effects. However, not much is known about mating gene regulation, since the mixture of mating types complicates detailed investigations. In this study, we developed K. phaffii strains with stable mating types by deletion of the inverted-repeat region required for mating-type switching. These heterothallic strains retain their ability to mate with cells of the opposite mating type and were used to further elucidate mating gene regulation. Functional analysis of MAT mutant strains revealed the essential role of MATa2 and MATα1 in diploid cell formation. Disruption of MATa1 or MATα2 did not affect mating; however, in diploid cells, both genes are required for sporulation and the repression of shmoo formation. The heterothallic strains generated in this study allowed the first detailed characterization of mating gene regulation in K. phaffii They will be a valuable tool for further studies investigating cell-type-specific behavior and will enable in-depth genetic analyses and strain hybridization in this industrially relevant yeast species.

Integrated Bottom-Up and Top-Down Liquid Chromatography-Mass Spectrometry for Characterization of Recombinant Human Growth Hormone Degradation Products

With the advent of biosimilars to the U.S. market, it is important to have better analytical tools to ensure product quality from batch to batch. In addition, the recent popularity of using a continuous process for production of biopharmaceuticals, the traditional bottom-up method, alone for product characterization and quality analysis is no longer sufficient. Bottom-up method requires large amounts of material for analysis and is labor-intensive and time-consuming. Additionally, in this analysis, digestion of the protein with enzymes such as trypsin could induce artifacts and modifications which would increase the complexity of the analysis. On the other hand, a top-down method requires a minimum amount of sample and allows for analysis of the intact protein mass and sequence generated from fragmentation within the instrument. However, fragmentation usually occurs at the N-terminal and C-terminal ends of the protein with less internal fragmentation. Herein, we combine the use of the complementary techniques, a top-down and bottom-up method, for the characterization of human growth hormone degradation products. Notably, our approach required small amounts of sample, which is a requirement due to the sample constraints of small scale manufacturing. Using this approach, we were able to characterize various protein variants, including post-translational modifications such as oxidation and deamidation, residual leader sequence, and proteolytic cleavage. Thus, we were able to highlight the complementarity of top-down and bottom-up approaches, which achieved the characterization of a wide range of product variants in samples of human growth hormone secreted from Pichia pastoris.