Comprehensive structural annotation of Pichia pastoris transcriptome and the response to various carbon sources using deep paired-end RNA sequencing

Engineering Yeast Peroxisomes for α-Bisabolene Production from Sole Methanol with the Aid of Proteomic Analysis

Microbial metabolic engineering provides a feasible approach to sustainably produce advanced biofuels and biochemicals from renewable feedstocks. Methanol is an ideal feedstock since it can be massively produced from CO2 through green energy, such as solar energy. However, engineering microbes to transform methanol and overproduce chemicals is challenging. Notably, the microbial production of isoprenoids from methanol is still rarely reported. Here, we extensively engineered Pichia pastoris (syn. Komagataella phaffii) for the overproduction of sesquiterpene α-bisabolene from sole methanol by optimizing the mevalonate pathway and peroxisomal compartmentalization. Furthermore, through label-free quantification (LFQ) proteomic analysis of the engineered strains, we identified the key bottlenecks in the peroxisomal targeting pathway, and overexpressing the limiting enzyme EfmvaE significantly improved α-bisabolene production to 212 mg/L with the peroxisomal pathway. The engineered strain LH122 with the optimized peroxisomal pathway produced 1.1 g/L α-bisabolene under fed-batch fermentation in shake flasks, achieving a 69% increase over that of the cytosolic pathway. This study provides a viable approach for overproducing isoprenoid from sole methanol in engineered yeast cell factories and shows that proteomic analysis can help optimize the organelle compartmentalized pathways to enhance chemical production.

Scalable protein production by Komagataella phaffii enabled by ARS plasmids and carbon source-based selection

BACKGROUND

Most recombinant Komagataella phaffii (Pichia pastoris) strains for protein production are generated by genomic integration of expression cassettes. The clonal variability in gene copy numbers, integration loci and consequently product titers limit the aptitude for high throughput applications in drug discovery, enzyme engineering or most comparative analyses of genetic elements such as promoters or secretion signals. Circular episomal plasmids with an autonomously replicating sequence (ARS), an alternative which would alleviate some of these limitations, are inherently unstable in K. phaffii. Permanent selection pressure, mostly enabled by antibiotic resistance or auxotrophy markers, is crucial for plasmid maintenance and hardly scalable for production. The establishment and use of extrachromosomal ARS plasmids with key genes of the glycerol metabolism (glycerol kinase 1, GUT1, and triosephosphate isomerase 1, TPI1) as selection markers was investigated to obtain a system with high transformation rates that can be directly used for scalable production processes in lab scale bioreactors.

RESULTS

In micro-scale deep-well plate experiments, ARS plasmids employing the Ashbya gossypii TEF1 (transcription elongation factor 1) promoter to regulate transcription of the marker gene were found to deliver high transformation efficiencies and the best performances with the reporter protein (CalB, lipase B of Candida antarctica) for both, the GUT1- and TPI1-based, marker systems. The GUT1 marker-bearing strain surpassed the reference strain with integrated expression cassette by 46% upon re-evaluation in shake flask cultures regarding CalB production, while the TPI1 system was slightly less productive compared to the control. In 5 L bioreactor methanol-free fed-batch cultivations, the episomal production system employing the GUT1 marker led to 100% increased CalB activity in the culture supernatant compared to integration construct.

CONCLUSIONS

For the first time, a scalable and methanol-independent expression system for recombinant protein production for K. phaffii using episomal expression vectors was demonstrated. Expression of the GUT1 selection marker gene of the new ARS plasmids was refined by employing the TEF1 promoter of A. gossypii. Additionally, the antibiotic-free marker toolbox for K. phaffii was expanded by the TPI1 marker system, which proved to be similarly suited for the use in episomal plasmids as well as integrative expression constructs for the purpose of recombinant protein production.

Transcriptomic Profiling at the Maternal-to-Zygotic Transition in Leech, Helobdella austinensis

The glossiphoniid leech, Helobdella austinensis, is an experimentally tractable member of the superphylum, Lophotrochozoa. Its large embryonic cells, stereotyped asymmetric cell divisions and ex vivo development capabilities makes it a favorable model for studying the molecular and cellular events of a representative spiralian. In this study, we focused on a narrow developmental time window of ~6-8 h, comprising stages just prior to and immediately following zygote deposition. Employing RNA-Seq methodology, we identified differentially expressed transcripts at this fundamental ontogenic boundary, known as the maternal-to-zygotic transition (MZT). Gene expression changes were characterized by the massive degradation of maternal RNAs (~45%) coupled with the rapid transcription of ~5000 zygotic genes (~20% of the genome) in the first mitotic cell cycle. The latter transcripts encoded a mixture of cell maintenance and regulatory proteins that predictably influence downstream developmental events.

Enabling growth-decoupled Komagataella phaffii recombinant protein production based on the methanol-free PDH promoter

The current transition towards the circular bioeconomy requires a rational development of biorefineries to sustainably fulfill the present demands. The use of Komagataella phaffii (Pichia pastoris) can meet this challenge, since it has the capability to use crude glycerol as a carbon-source, a by-product from the biodiesel industry, while producing high- and low-added value products. Recombinant protein production (RPP) using K. phaffii has often been driven either by the methanol induced AOX1 promoter (P_AOX1) and/or the constitutive GAP promoter (P_GAP). In the last years, strong efforts have been focused on developing novel expression systems that expand the toolbox variety of K. phaffii to efficiently produce diverse proteins that requires different strategies. In this work, a study was conducted towards the development of methanol-free expression system based on a heat-shock gene promoter (P_DH) using glycerol as sole carbon source. Using this promoter, the recombinant expression is strongly induced in carbon-starving conditions. The classical P_GAP was used as a benchmark, taking for both strains the lipase B from Candida antarctica (CalB) as model protein. Titer of CalB expressed under P_DH outperformed P_GAP controlled expression in shake-flask cultivations when using a slow-release continuous feeding technology, confirming that P_DH is induced under pseudo-starving conditions. This increase was also confirmed in fed-batch cultivations. Several optimization rounds were carried out for P_DH under different feeding and osmolarity conditions. In all of them the P_DH controlled process outperformed the P_GAP one in regard to CalB titer. The best P_DH approach reached 3.6-fold more specific productivity than P_GAP fed-batch at low μ. Compared to the optimum approach for P_GAP-based process, the best P_DH fed-batch strategy resulted in 2.3-fold higher titer, while the specific productivity was very similar. To summarize, P_DH is an inducible promoter that exhibited a non-coupled growth regulation showing high performance, which provides a methanol-free additional solution to the usual growth-coupled systems for RPP. Thus, this novel system emerges as a potential alternative for K. phaffii RPP bioprocess and for revaluing crude glycerol, promoting the transition towards a circular economy.

Comparative transcriptomic analysis-based identification of the regulation of foreign proteins with different stabilities expressed in Pichia pastoris

Introduction

The industrial yeast Pichia pastoris is widely used as a cell factory to produce proteins, chemicals and advanced biofuels. We have previously constructed P. pastoris strains that overexpress protein disulfide isomerase (PDI), which is a kind of molecular chaperone that can improve the expression of an exogenous protein when they are co-expressed. Chicken cystatin (cC) is a highly thermostable cysteine protease inhibitor and a homologous protein of human cystatin C (HCC). Wild-type cC and the two mutants, I66Q and ΔW (a truncated cC lacking the á-helix 2) represent proteins with different degrees of stability.

Methods

Wild-type cC, I66Q and ΔW were each overexpressed in P. pastoris without and with the coexpression of PDI and their extracellular levels were determined and compared. Transcriptomic profiling was performed to compare the changes in the main signaling pathways and cell components (other than endoplasmic reticulum quality control system represented by molecular chaperones) in P. pastoris in response to intracellular folding stress caused by the expression of exogenous proteins with different stabilities. Finally, hub genes hunting was also performed.

Results and discussion

The coexpression of PDI was able to increase the extracellular levels of both wild-type cC and the two mutants, indicating that overexpression of PDI could prevent the misfolding of unstable proteins or promote the degradation of the misfolded proteins to some extent. For P. pastoris cells that expressed the I66Q or ΔW mutant, GO (Gene Ontology) and KEGG (Kyoto Encyclopedia of Genes and Genomes) analyses of the common DEGs in these cells revealed a significant upregulation of the genes involved in protein processing, but a significant downregulation of the genes enriched in the Ribosome, TCA and Glycolysis/Gluconeogenesis pathways. Hub genes hunting indicated that the most downregulated ribosome protein, C4QXU7 in this case, might be an important target protein that could be manipulated to increase the expression of foreign proteins, especially proteins with a certain degree of instability.

Conclusion

These findings should shed new light on our understanding of the regulatory mechanism in yeast cells that responds to intracellular folding stress, providing valuable information for the development of a convenient platform that could improve the efficiency of heterologous protein expression in P. 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.

Comparative transcriptome and metabolome analyses reveal the methanol dissimilation pathway of Pichia pastoris

BACKGROUND

Pichia pastoris (Komagataella phaffii) is a model organism widely used for the recombinant expression of eukaryotic proteins, and it can metabolize methanol as its sole carbon and energy source. Methanol is oxidized to formaldehyde by alcohol oxidase (AOX). In the dissimilation pathway, formaldehyde is oxidized to CO₂ by formaldehyde dehydrogenase (FLD), S-hydroxymethyl glutathione hydrolase (FGH) and formate dehydrogenase (FDH).

RESULTS

The transcriptome and metabolome of P. pastoris were determined under methanol cultivation when its dissimilation pathway cut off. Firstly, Δfld and Δfgh were significantly different compared to the wild type (GS115), with a 60.98% and 23.66% reduction in biomass, respectively. The differential metabolites between GS115 and Δfld were mainly enriched in ABC transporters, amino acid biosynthesis, and protein digestion and absorption. Secondly, comparative transcriptome between knockout and wild type strains showed that oxidative phosphorylation, glycolysis and the TCA cycle were downregulated, while alcohol metabolism, proteasomes, autophagy and peroxisomes were upregulated. Interestingly, the down-regulation of the oxidative phosphorylation pathway was positively correlated with the gene order of dissimilation pathway knockdown. In addition, there were significant differences in amino acid metabolism and glutathione redox cycling that raised our concerns about formaldehyde sorption in cells.

CONCLUSIONS

This is the first time that integrity of dissimilation pathway analysis based on transcriptomics and metabolomics was carried out in Pichia pastoris. The blockage of dissimilation pathway significantly down-regulates the level of oxidative phosphorylation and weakens the methanol assimilation pathway to the point where deficiencies in energy supply and carbon fixation result in inefficient biomass accumulation and genetic replication. In addition, transcriptional upregulation of the proteasome and autophagy may be a stress response to resolve formaldehyde-induced DNA-protein crosslinking.

Engineering of Promoters for Gene Expression in Pichia pastoris

The availability of exceptionally strong and tightly regulated promoters is a key feature of Komagataella phaffii (syn. Pichia pastoris), a widely applied yeast expression system for heterologous protein production. Most commonly, the methanol-inducible promoter of the alcohol oxidase 1 gene (P_AOX1) and the constitutive promoter of the glyceraldehyde 3 phosphate dehydrogenase gene (P_GAP) have been used. Recently, also promising novel constitutive (P_GCW14), regulated (P_GTH1, P_CAT1), and bidirectional promoters (histone promoters and synthetic hybrid variants) have been reported.As natural promoters showed so far limited tunability of expression levels and regulatory profiles, various promoter engineering efforts have been undertaken for P. pastoris . P_AOX1, P_DAS2, P_GAP, and P_GCW14 have been engineered by systematic deletion studies or random mutagenesis of upstream regulatory sequences. New engineering strategies have focused on P_AOX1 core promoter modifications by random or rational approaches and transcriptional regulatory circuits to render P_AOX1 independent of methanol induction. These promoter engineering efforts in P. pastoris have resulted in improved, sequence-diversified synthetic promoter variants allowing coordinated fine-tuning of gene expression for a multitude of biotechnological applications.

Global transcriptome profiling reveals genes responding to overproduction of a small secretory, a high cysteine- and a high glycosylation-bearing protein in Yarrowia lipolytica

•

Secretion of a protein results in > 10-fold higher titer compared to its retention.

•

Overproduction of rs-Prots induces oxidative stress and detoxification response.

•

Excessive vacuolar protein degradation limits rs-Prot production in Y. lipolytica.

•

Non-classical export protein NCE102 is upregulated upon rs-Prot overloading.

•

Downregulation of cyclin CLN1 marks growth arrest in G1 under rs-Prot synthesis.

Investigation of the yeast cell’s response to recombinant secretory protein (rs-Prot) overproduction is relevant for both basic and applied research. Imbalance, overloading or stress within this process impacts the whole cell. In the present study, by using steady-state cultures and transcriptomics, we investigated the cellular response of Yarrowia lipolytica challenged with high-level expression of genes encoding proteins with significantly different biochemical characteristics: a small protein retained within the cell i) or secreted ii), a medium size secretory protein with a high number of disulfide bonds iii), or glycosylation sites iv). Extensive analysis of omics data, supported by careful manual curation, led to some anticipated observations on oxidative and unfolded protein stress (CTT1, PXMP2/4, HAC1), glycosylation (ALGs, KTRs, MNTs, MNNs), folding and translocation (SSAs, SSEs) but also generated new exciting knowledge on non-conventional protein secretion (NCE102), transcriptional regulators (FLO11, MHY1, D01353 g, RSFA, E23925g or MAF1), vacuolar proteolysis targets in Y. lipolytica (ATGs, VPSs, HSE1, PRB1, PRC1, PEP4) or growth arrest (CLN1) upon rs-Prots overproduction.

Transcriptomic analysis reveals MAPK signaling pathways affect the autolysis in baker's yeast

Yeast autolysis refers to the process in which cells degrade and release intracellular contents under specific conditions by endogenous enzymes such as proteases, nucleases and lipid enzymes. Protein-rich baker's yeast is widely used to produce yeast extract in food industry, however, the molecular mechanism related to baker's yeast autolysis is still unclear. In this study, RNA-seq technology and biochemical analysis were performed to analyze the autolysis processes in baker's yeast. The differentially expressed genes (DEGs), 27 autolysis-related euKaryotic Ortholog Groups (KOG) and three types of autolysis-induced Gene Ontology (GO) were identified and analyzed in detail. A total of 143 Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways under autolysis were also assigned. Interestingly, the DEGs were significantly enriched in the mitogen-activated protein kinase (MAPK) signaling pathways and metabolic pathways, and key genes MID2, MTL1, SLT2, PTP2, HKR1 and GPD1 may play important roles in autolysis. Further quantitative PCR was performed to verify the expression pattern in baker's yeast autolysis. Together, all these results indicated that MAPK pathways might play an essential role during autolysis process through inhibiting the metabolism and disrupting cell wall in baker's yeast. This result may provide important clues for the in-depth interpretation of the yeast autolysis mechanism.

Transcriptome analysis on fructose as the sole carbon source enhancing perylenequinones production of endophytic fungus Shiraia sp. Slf14

Perylenequinones (PQ), a class of naturally occurring polypeptides, are widely used as a clinical drug for treating skin diseases and as a photodynamic therapy against cancers and viruses. In this study, the effects of different carbon sources on PQ biosynthesis by Shiraia sp. Slf14 were compared, and the underlying molecular mechanism of fructose as the sole carbon to enhance PQ production was investigated by transcriptome analysis. The results indicated that fructose enhanced PQ yield to 1753.64 mg/L, which was 1.73-fold higher than that obtained with glucose. Comparative transcriptome analysis demonstrated that most of the upregulated genes were related to transport systems, energy and central carbon metabolism in Shiraia sp. Slf14 cultured in fructose. The genes involved in glycolysis and pentose phosphate pathways, and encoding citrate synthase, ATP-citrate lyase, and acetyl-CoA carboxylase were substantially upregulated, resulting in increased overall acetyl-CoA and malonyl-CoA production. However, genes involved in gluconeogenesis, glyoxylate cycle pathway, and fatty acid synthesis were significantly downregulated, resulting in higher acetyl-CoA influx for PQ formation. In particular, the putative PQ biosynthetic cluster was upregulated in Shiraia sp. Slf14 cultured in fructose, leading to a significant increase in PQ production. The results of real-time qRT-PCR and related enzyme activities were also consistent with those of transcriptome analysis. These findings provide a remarkable insight into the underlying mechanism of PQ biosynthesis and pave the way for improvements in PQ production by Shiraia sp. Slf14.

Transcriptomic Analysis of Pichia pastoris (Komagataella phaffii) GS115 During Heterologous Protein Production Using a High-Cell-Density Fed-Batch Cultivation Strategy

Pichia pastoris (Komagataella phaffii) is a methylotrophic yeast that is widely used in industry as a host system for heterologous protein expression. Heterologous gene expression is typically facilitated by strongly inducible promoters derived from methanol utilization genes or constitutive glycolytic promoters. However, protein production is usually accomplished by a fed-batch induction process, which is known to negatively affect cell physiology, resulting in limited protein yields and quality. To assess how yields of exogenous proteins can be increased and to further understand the physiological response of P. pastoris to the carbon conversion of glycerol and methanol, as well as the continuous induction of methanol, we analyzed recombinant protein production in a 10,000-L fed-batch culture. Furthermore, we investigated gene expression during the yeast cell culture phase, glycerol feed phase, glycerol-methanol mixture feed (GM) phase, and at different time points following methanol induction using RNA-Seq. We report that the addition of the GM phase may help to alleviate the adverse effects of methanol addition (alone) on P. pastoris cells. Secondly, enhanced upregulation of the mitogen-activated protein kinase (MAPK) signaling pathway was observed in P. pastoris following methanol induction. The MAPK signaling pathway may be related to P. pastoris cell growth and may regulate the alcohol oxidase1 (AOX1) promoter via regulatory factors activated by methanol-mediated stimulation. Thirdly, the unfolded protein response (UPR) and ER-associated degradation (ERAD) pathways were not significantly upregulated during the methanol induction period. These results imply that the presence of unfolded or misfolded phytase protein did not represent a serious problem in our study. Finally, the upregulation of the autophagy pathway during the methanol induction phase may be related to the degradation of damaged peroxisomes but not to the production of phytase. This work describes the metabolic characteristics of P. pastoris during heterologous protein production under high-cell-density fed-batch cultivation. We believe that the results of this study will aid further in-depth studies of P. pastoris heterologous protein expression, regulation, and secretory mechanisms.

Screening for functional IRESes using α-complementation system of β-galactosidase in Pichia pastoris

BACKGROUND

Pichia pastoris is becoming a promising chassis cell for metabolic engineering and synthetic biology after its whole genome and transcriptome sequenced. However, the current systems for multigene co-expression in P. pastoris are not efficient. The internal ribosome entry site (IRES) has an ability to recruit the ribosome to initiate protein synthesis by cap-independent translation manner. This study seeks to screen IRES sequences that are functional in P. pastoris, which will allow P. pastoris to express multiple proteins in a single mRNA and increase its efficacy as a platform for metabolic engineering and synthetic biology.

RESULTS

In order to efficiently screen the IRES sequences, we first set out to create a screening system using LacZ gene. Due to the cryptic transcription of the LacZ gene, we established the α-complementation system of β-galactosidase in P. pastoris with the optimum length of the α-complementing peptide at ~ 92 amino acids. The optimal α-complementing peptide was then used as the second reporter to screen IRESes in the engineered GS115 expressing the corresponding ω-peptide. A total of 34 reported IRESes were screened. After ruling out all false positive or negative IRESes, only seven IRESes were functional in P. pastoris, which were from TEV, PVY, RhPV, TRV, KSHV, crTMV viruses and the 5'-UTR of the YAP1 gene of S. cerevisiae.

CONCLUSIONS

We showed here that α-complementation also works in P. pastoris and it can be used in a variety of in vivo studies. The functional IRESes screened in this study can be used to introduce multiple genes into P. pastoris via a prokaryotic-like polycistronic manner, which provided new efficient tools for metabolic engineering and synthetic biology researches in P. pastoris.

Deletion of Gcw13 represses autophagy in Pichia pastoris cells grown in methanol medium with sufficient amino acids

OBJECTIVE

The purpose of this article is to study the underlying cause of the induction of autophagy in Pichia pastoris cells grown in amino acid-rich methanol medium during methanol adaptation.

RESULTS

Autophagy was induced in P. pastoris GS115 when cells were grown in amino acid-rich methanol medium. Transcriptome analysis revealed that genes involved in amino acid biosynthesis were upregulated. The deletion of Gcw13, a GPI-anchored protein that plays a role in the endocytosis of the general amino acid permease Gap1, resulted in the inhibition of autophagy, the activation of TORC1 and an increase in the uptake of glutamine and asparagine in methanol-grown cells.

CONCLUSIONS

Our results demonstrated that the autophagy induced in P. pastoris cells grown in amino acid-rich methanol medium was nitrogen source independent and may be due to a Gcw13-dependent decrease in amino acid uptake during methanol adaptation.

Statistically Designed Medium Reveals Interactions between Metabolism and Genetic Information Processing for Production of Stable Human Serum Albumin in Pichia pastoris

Human serum albumin (HSA), sourced from human serum, has been an important therapeutic protein for several decades. Pichia pastoris is strongly considered as an expression platform, but proteolytic degradation of recombinant HSA in the culture filtrate remains a major bottleneck for use of this system. In this study, we have reported the development of a medium that minimized proteolytic degradation across different copy number constructs. A synthetic codon-optimized copy of HSA was cloned downstream of α-factor secretory signal sequence and expressed in P. pastoris under the control of Alcohol oxidase 1 promoter. A two-copy expression cassette was also prepared. Culture conditions and medium components were identified and optimized using statistical tools to develop a medium that supported stable production of HSA. Comparative analysis of transcriptome data obtained by cultivation on optimized and unoptimized medium indicated upregulation of genes involved in methanol metabolism, alternate nitrogen assimilation, and DNA transcription, whereas enzymes of translation and secretion were downregulated. Several new genes were identified that could serve as possible targets for strain engineering of this yeast.

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.

Enhancing the efficiency of the Pichia pastoris AOX1 promoter via the synthetic positive feedback circuit of transcription factor Mxr1

Background

The methanol-regulated AOX1 promoter (P_AOX1) is the most widely used promoter in the production of recombinant proteins in the methylotrophic yeast Pichia pastoris. However, as the tight regulation and methanol dependence of P_AOX1 restricts its application, it is necessary to develop a flexible induction system to avoid the problems of methanol without losing the advantages of P_AOX1. The availability of synthetic biology tools enables researchers to reprogram the cellular behaviour of P. pastoris to achieve this goal.

Results

The characteristics of P_AOX1 are highly related to the expression profile of methanol expression regulator 1 (Mxr1). In this study, we applied a biologically inspired strategy to reprogram regulatory networks in P. pastoris. A reprogrammed P. pastoris was constructed by inserting a synthetic positive feedback circuit of Mxr1 driven by a weak AOX2 promoter (P_AOX2). This novel approach enhanced P_AOX1 efficiency by providing extra Mxr1 and generated switchable Mxr1 expression to allow P_AOX1 to be induced under glycerol starvation or carbon-free conditions. Additionally, the inhibitory effect of glycerol on P_AOX1 was retained because the synthetic circuit was not activated in response to glycerol. Using green fluorescent protein as a demonstration, this reprogrammed P. pastoris strain displayed stronger fluorescence intensity than non-reprogrammed cells under both methanol induction and glycerol starvation. Moreover, with single-chain variable fragment (scFv) as the model protein, increases in extracellular scFv productivity of 98 and 269% were observed in Mxr1-reprogrammed cells under methanol induction and glycerol starvation, respectively, compared to productivity in non-reprogrammed cells under methanol induction.

Conclusions

We successfully demonstrate that the synthetic positive feedback circuit of Mxr1 enhances recombinant protein production efficiency in P. pastoris and create a methanol-free induction system to eliminate the potential risks of methanol.

Electronic supplementary material

The online version of this article (10.1186/s12896-018-0492-4) contains supplementary material, which is available to authorized users.

Presence of protein production enhancers results in significantly higher methanol-induced protein production in Pichia pastoris

Background

The yeast Komagataella phaffii, better known as Pichia pastoris, is a commonly used host for recombinant protein production. Here expression vectors are reported that address the different steps of the transcription–translation–secretion pathway of heterologous protein production.

Results

Transcription and translation enhancing elements were introduced in an expression cassette for the production of recombinant Aspergillus niger feruloyl esterase A. The yield was increased by threefold as compared to the yield without these elements. Multiple copy strains were selected using a zeocin resistance marker in the expression cassette and showed another sixfold higher yield. Modification of the C-terminal amino acid sequence of the secretion signal did not significantly improve the production yield. Similar data were obtained for the production of another protein, recombinant human interleukin 8. Upscaling to fed-batch fermentation conditions resulted in a twofold increase for reference strains, while for strains with enhancing elements a tenfold improvement was observed.

Conclusions

Pichia pastoris is used for recombinant protein production in industrial fermentations. By addressing the transcription and translation of mRNA coding for recombinant protein, significant yield improvement was obtained. The yield improvement obtained under microscale conditions was maintained under fed-batch fermentation conditions. These data demonstrate the potential of these expression vectors for large scale application as improved production of proteins has major implications on the economics and sustainability of biocatalyst dependent production processes e.g. for the production of pharmaceuticals and for the bioconversions of complex molecules.

Investigating the effect of carbon source on rabies virus glycoprotein production in Pichia pastoris by a transcriptomic approach

Several factors affect protein expression in Pichia pastoris, one among them is the carbon source. In this work, we studied the effect of this factor on the expression level of rabies virus glycoprotein (RABV-G) in two recombinant clones harboring seven copies of the gene of interest. The expression was driven either by the constitutive glyceraldehyde-3-phosphate dehydrogenase (GAP) promoter or the inducible alcohol oxidase1 (AOX1) promoter. Clones were compared in terms of cell physiology and carbon source metabolism. The transcription levels of 16 key genes involved in the central metabolic pathway, the methanol catabolism, and the oxidative stress were investigated in both clones. Cell size, as a parameter reflecting cell physiological changes, was also monitored. Our results showed that when glucose was used as the sole carbon source, large cells were obtained. Transcript levels of the genes of the central metabolic pathway were also upregulated, whereas antioxidative gene transcript levels were low. By contrast, the use of methanol as a carbon source generated small cells and a shift in carbon metabolism toward the dissimilatory pathway by the upregulation of formaldehyde dehydrogenase gene and the downregulation of those of the central metabolic. These observations are in favor of the use of glucose to enhance the expression of RABV-G in P. pastoris.

Increased dosage of AOX1 promoter-regulated expression cassettes leads to transcription attenuation of the methanol metabolism in Pichia pastoris

The methanol-regulated alcohol oxidase promoter (P_AOX1) of Pichia pastoris is one of the strongest promoters for heterologous gene expression in this methylotrophic yeast. Although increasing gene dosage is one of the most common strategies to increase recombinant protein productivities, the increase of gene dosage of Rhizopus oryzae lipase (ROL) in P. pastoris has been previously shown to reduce cell growth, lipase production and substrate consumption in high-copy strains. To better assess that physiological response, transcriptomics analysis was performed of a subset of strains with 1 to 15 ROL copies. The macroscopic physiological parameters confirm that growth yield and carbon uptake rate are gene dosage dependent, and were supported by the transcriptomic data, showing the impact of increased dosage of AOX1 promoter-regulated expression cassettes on P. pastoris physiology under steady methanolic growth conditions. Remarkably, increased number of cassettes led to transcription attenuation of the methanol metabolism and peroxisome biogenesis in P. pastoris, concomitant with reduced secretion levels of the heterologous product. Moreover, our data also point to a block in ROL mRNA translation in the higher ROL-copies constructs, while the low productivities of multi-copy strains under steady growth conditions do not appear to be directly related to UPR and ERAD induction.

High-density linkage map construction and QTL analysis for earliness-related traits in Gossypium hirsutum L

Background

Gossypium hirsutum L., or upland cotton, is an important renewable resource for textile fiber. To enhance understanding of the genetic basis of cotton earliness, we constructed an intra-specific recombinant inbred line population (RIL) containing 137 lines, and performed linkage map construction and quantitative trait locus (QTL) mapping.

Results

Using restriction-site associated DNA sequencing, a genetic map composed of 6,434 loci, including 6,295 single nucleotide polymorphisms and 139 simple sequence repeat loci, was developed from RIL population. This map spanned 4,071.98 cM, with an average distance of 0.63 cM between adjacent markers. A total of 247 QTLs for six earliness-related traits were detected in 6 consecutive years. In addition, 55 QTL coincidence regions representing more than 60 % of total QTLs were found on 22 chromosomes, which indicated that several earliness-related traits might be simultaneously improved. Fine-mapping of a 2-Mb region on chromosome D3 associated with five stable QTLs between Marker25958 and Marker25963 revealed that lines containing alleles derived from CCRI36 in this region exhibited smaller phenotypes and earlier maturity. One candidate gene (EMF2) was predicted and validated by quantitative real-time PCR in early-, medium- and late-maturing cultivars from 3- to 6-leaf stages, with highest expression level in early-maturing cultivar, CCRI74, lowest expression level in late-maturing cultivar, Bomian1.

Conclusions

We developed an SNP-based genetic map, and this map is the first high-density genetic map for short-season cotton and has the potential to provide deeper insights into earliness. Cotton earliness-related QTLs and QTL coincidence regions will provide useful materials for QTL fine mapping, gene positional cloning and MAS. And the gene, EMF2, is promising for further study.

Electronic supplementary material

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

Refined Pichia pastoris reference genome sequence

Strains of the species Komagataella phaffii are the most frequently used "Pichia pastoris" strains employed for recombinant protein production as well as studies on peroxisome biogenesis, autophagy and secretory pathway analyses. Genome sequencing of several different P. pastoris strains has provided the foundation for understanding these cellular functions in recent genomics, transcriptomics and proteomics experiments. This experimentation has identified mistakes, gaps and incorrectly annotated open reading frames in the previously published draft genome sequences. Here, a refined reference genome is presented, generated with genome and transcriptome sequencing data from multiple P. pastoris strains. Twelve major sequence gaps from 20 to 6000 base pairs were closed and 5111 out of 5256 putative open reading frames were manually curated and confirmed by RNA-seq and published LC-MS/MS data, including the addition of new open reading frames (ORFs) and a reduction in the number of spliced genes from 797 to 571. One chromosomal fragment of 76kbp between two previous gaps on chromosome 1 and another 134kbp fragment at the end of chromosome 4, as well as several shorter fragments needed re-orientation. In total more than 500 positions in the genome have been corrected. This reference genome is presented with new chromosomal numbering, positioning ribosomal repeats at the distal ends of the four chromosomes, and includes predicted chromosomal centromeres as well as the sequence of two linear cytoplasmic plasmids of 13.1 and 9.5kbp found in some strains of P. pastoris.

Centromeres of the Yeast Komagataella phaffii (Pichia pastoris) Have a Simple Inverted-Repeat Structure

Centromere organization has evolved dramatically in one clade of fungi, the Saccharomycotina. These yeasts have lost the ability to make normal eukaryotic heterochromatin with histone H3K9 methylation, which is a major component of pericentromeric regions in other eukaryotes. Following this loss, several different types of centromere emerged, including two types of sequence-defined (“point”) centromeres, and the epigenetically defined “small regional” centromeres of Candida albicans. Here we report that centromeres of the methylotrophic yeast Komagataella phaffii (formerly called Pichia pastoris) are structurally defined. Each of its four centromeres consists of a 2-kb inverted repeat (IR) flanking a 1-kb central core (mid) region. The four centromeres are unrelated in sequence. CenH3 (Cse4) binds strongly to the cores, with a decreasing gradient along the IRs. This mode of organization resembles Schizosaccharomyces pombe centromeres but is much more compact and lacks the extensive flanking heterochromatic otr repeats. Different isolates of K. phaffii show polymorphism for the orientation of the mid regions, due to recombination in the IRs. CEN4 is located within a 138-kb region that changes orientation during mating-type switching, but switching does not induce recombination of centromeric IRs. Our results demonstrate that evolutionary transitions in centromere organization have occurred in multiple yeast clades.

Comparative genomics and transcriptomics of Pichia pastoris

Background

Pichia pastoris has emerged as an important alternative host for producing recombinant biopharmaceuticals, owing to its high cultivation density, low host cell protein burden, and the development of strains with humanized glycosylation. Despite its demonstrated utility, relatively little strain engineering has been performed to improve Pichia, due in part to the limited number and inconsistent frameworks of reported genomes and transcriptomes. Furthermore, the co-mingling of genomic, transcriptomic and fermentation data collected about Komagataella pastoris and Komagataella phaffii, the two strains co-branded as Pichia, has generated confusion about host performance for these genetically distinct species. Generation of comparative high-quality genomes and transcriptomes will enable meaningful comparisons between the organisms, and potentially inform distinct biotechnological utilies for each species.

Results

Here, we present a comprehensive and standardized comparative analysis of the genomic features of the three most commonly used strains comprising the tradename Pichia: K. pastoris wild-type, K. phaffii wild-type, and K. phaffii GS115. We used a combination of long-read (PacBio) and short-read (Illumina) sequencing technologies to achieve over 1000X coverage of each genome. Construction of individual genomes was then performed using as few as seven individual contigs to create gap-free assemblies. We found substantial syntenic rearrangements between the species and characterized a linear plasmid present in K. phaffii. Comparative analyses between K. phaffii genomes enabled the characterization of the mutational landscape of the GS115 strain. We identified and examined 35 non-synonomous coding mutations present in GS115, many of which are likely to impact strain performance. Additionally, we investigated transcriptomic profiles of gene expression for both species during cultivation on various carbon sources. We observed that the most highly transcribed genes in both organisms were consistently highly expressed in all three carbon sources examined. We also observed selective expression of certain genes in each carbon source, including many sequences not previously reported as promoters for expression of heterologous proteins in yeasts.

Conclusions

Our studies establish a foundation for understanding critical relationships between genome structure, cultivation conditions and gene expression. The resources we report here will inform and facilitate rational, organism-wide strain engineering for improved utility as a host for protein production.

Electronic supplementary material

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

Transcriptomic Analysis of the Endangered Neritid Species Clithon retropictus: De Novo Assembly, Functional Annotation, and Marker Discovery

An aquatic gastropod belonging to the family Neritidae, Clithon retropictus is listed as an endangered class II species in South Korea. The lack of information on its genomic background limits the ability to obtain functional data resources and inhibits informed conservation planning for this species. In the present study, the transcriptomic sequencing and de novo assembly of C. retropictus generated a total of 241,696,750 high-quality reads. These assembled to 282,838 unigenes with mean and N50 lengths of 736.9 and 1201 base pairs, respectively. Of these, 125,616 unigenes were subjected to annotation analysis with known proteins in Protostome DB, COG, GO, and KEGG protein databases (BLASTX; E ≤ 0.00001) and with known nucleotides in the Unigene database (BLASTN; E ≤ 0.00001). The GO analysis indicated that cellular process, cell, and catalytic activity are the predominant GO terms in the biological process, cellular component, and molecular function categories, respectively. In addition, 2093 unigenes were distributed in 107 different KEGG pathways. Furthermore, 49,280 simple sequence repeats were identified in the unigenes (>1 kilobase sequences). This is the first report on the identification of transcriptomic and microsatellite resources for C. retropictus, which opens up the possibility of exploring traits related to the adaptation and acclimatization of this species.

Epistasis together with partial dominance, over-dominance and QTL by environment interactions contribute to yield heterosis in upland cotton

QTL mapping based on backcross and RIL populations suggests that epistasis together with partial dominance, over-dominance and their environmental interactions of QTLs play an important role in yield heterosis in upland cotton. A backcross population (BC) was constructed to explore the genetic basis of heterosis in upland cotton (Gossypium hirsutum L.). For yield and yield components, recombinant inbred line (RIL) and BC populations were evaluated simultaneously at three different locations. A total of 35 and 30 quantitative trait loci (QTLs) were detected based on the RILs and BC data, respectively. Six (16.7 %) additive QTLs, 19 (52.8 %) partial dominant QTLs and 11 (30.6 %) over-dominant QTLs were detected by single-locus analysis using composite interval mapping in BC population. QTLs detected for mid-parent heterosis (MPH) were mostly related to those detected in the BC population. No significant correlation was found between marker heterozygosity and performance. It indicated that heterozygosity was not always favorable for performance. Two-locus analysis revealed 46, 25 and 12 QTLs with main effects (M-QTLs), and 55, 63 and 33 QTLs involved in digenic interactions (E-QTLs) were detected for yield and yield components in RIL, BC and MPH, respectively. A large number of M-QTLs and E-QTLs showed QTL by environment interactions (QEs) in three environments. These results suggest that epistasis together with partial dominance, over-dominance and QEs all contribute to yield heterosis in upland cotton.

Internal Initiation of Translation of mRNA in the Methylotrophic Yeast Hansenula polymorpha

Besides regular cap-dependent translation of mRNA, eukaryotes exploit internal initiation of translation driven by internal ribosome entry sites (IRESs). It is supposed that internal initiation provides translation of cellular mRNAs under stress conditions where the cap-dependent initiation is reduced. A number of IRESs have been characterized in mammalian mRNAs, but only a few examples are known in lower eukaryotes, particularly in yeasts. Here we identified two IRESs in the thermotolerant methylotrophic yeast Hansenula polymorpha DL-1. These sites are located in 5'-untranslated regions of genes HPODL_02249 and HPODL_04025 encoding a hypothetical membrane protein and actin-binding protein, respectively. In Saccharomyces cerevisiae cells, both IRESs drive expression of a second gene of a bicistronic mRNA, as well as translation of hairpin-containing monocistronic mRNA. The possibility of spurious splicing or presence of a cryptic promoter in the IRES sequences was ruled out, indicating that expression of a second gene of a bicistronic mRNA was IRES-dependent. We evaluated IRES activity of both elements and found that under normal physiological conditions its contribution to the overall translation of the respective mRNAs in yeast cells is about 0.3-0.4%. Therefore, these results suggest that the IRES-dependent translation initiation mechanism exists in Hansenula polymorpha.

A Toolbox of Diverse Promoters Related to Methanol Utilization: Functionally Verified Parts for Heterologous Pathway Expression in Pichia pastoris

The heterologous expression of biosynthetic pathways for pharmaceutical or fine chemical production requires suitable expression hosts and vectors. In eukaryotes, the pathway flux is typically balanced by stoichiometric fine-tuning of reaction steps by varying the transcript levels of the genes involved. Regulated (inducible) promoters are desirable to allow a separation of pathway expression from cell growth. Ideally, the promoter sequences used should not be identical to avoid loss by recombination. The methylotrophic yeast Pichia pastoris is a commonly used protein production host, and single genes have been expressed at high levels using the methanol-inducible, strong, and tightly regulated promoter of the alcohol oxidase 1 gene (PAOX1). Here, we have studied the regulation of the P. pastoris methanol utilization (MUT) pathway to identify a useful set of promoters that (i) allow high coexpression and (ii) differ in DNA sequence to increase genetic stability. We noticed a pronounced involvement of the pentose phosphate pathway (PPP) and genes involved in the defense of reactive oxygen species (ROS), providing strong promoters that, in part, even outperform PAOX1 and offer novel regulatory profiles. We have applied these tightly regulated promoters together with novel terminators as useful tools for the expression of a heterologous biosynthetic pathway. With the synthetic biology toolbox presented here, P. pastoris is now equipped with one of the largest sets of strong and co-regulated promoters of any microbe, moving it from a protein production host to a general industrial biotechnology host.

Multi-omics data driven analysis establishes reference codon biases for synthetic gene design in microbial and mammalian cells

In this study, we analyzed multi-omics data and subsets thereof to establish reference codon usage biases for codon optimization in synthetic gene design. Specifically, publicly available genomic, transcriptomic, proteomic and translatomic data for microbial and mammalian expression hosts, Escherichia coli, Saccharomyces cerevisiae, Pichia pastoris and Chinese hamster ovary (CHO) cells, were compiled to derive their individual codon and codon pair frequencies. Then, host dependent and -omics specific codon biases were generated and compared by principal component analysis and hierarchical clustering. Interestingly, our results indicated the similar codon bias patterns of the highly expressed transcripts, highly abundant proteins, and efficiently translated mRNA in microbial cells, despite the general lack of correlation between mRNA and protein expression levels. However, for CHO cells, the codon bias patterns among various -omics subsets are not distinguishable, forming one cluster. Thus, we further investigated the effect of different input codon biases on codon optimized sequences using the codon context (CC) and individual codon usage (ICU) design parameters, via in silico case study on the expression of human IFNγ sequence in CHO cells. The results supported that CC is more robust design parameter than ICU for improved heterologous gene design.

Systems-level organization of yeast methylotrophic lifestyle

Background

Some yeasts have evolved a methylotrophic lifestyle enabling them to utilize the single carbon compound methanol as a carbon and energy source. Among them, Pichia pastoris (syn. Komagataella sp.) is frequently used for the production of heterologous proteins and also serves as a model organism for organelle research. Our current knowledge of methylotrophic lifestyle mainly derives from sophisticated biochemical studies which identified many key methanol utilization enzymes such as alcohol oxidase and dihydroxyacetone synthase and their localization to the peroxisomes. C1 assimilation is supposed to involve the pentose phosphate pathway, but details of these reactions are not known to date.

Results

In this work we analyzed the regulation patterns of 5,354 genes, 575 proteins, 141 metabolites, and fluxes through 39 reactions of P. pastoris comparing growth on glucose and on a methanol/glycerol mixed medium, respectively. Contrary to previous assumptions, we found that the entire methanol assimilation pathway is localized to peroxisomes rather than employing part of the cytosolic pentose phosphate pathway for xylulose-5-phosphate regeneration. For this purpose, P. pastoris (and presumably also other methylotrophic yeasts) have evolved a duplicated methanol inducible enzyme set targeted to peroxisomes. This compartmentalized cyclic C1 assimilation process termed xylose-monophosphate cycle resembles the principle of the Calvin cycle and uses sedoheptulose-1,7-bisphosphate as intermediate. The strong induction of alcohol oxidase, dihydroxyacetone synthase, formaldehyde and formate dehydrogenase, and catalase leads to high demand of their cofactors riboflavin, thiamine, nicotinamide, and heme, respectively, which is reflected in strong up-regulation of the respective synthesis pathways on methanol. Methanol-grown cells have a higher protein but lower free amino acid content, which can be attributed to the high drain towards methanol metabolic enzymes and their cofactors. In context with up-regulation of many amino acid biosynthesis genes or proteins, this visualizes an increased flux towards amino acid and protein synthesis which is reflected also in increased levels of transcripts and/or proteins related to ribosome biogenesis and translation.

Conclusions

Taken together, our work illustrates how concerted interpretation of multiple levels of systems biology data can contribute to elucidation of yet unknown cellular pathways and revolutionize our understanding of cellular biology.

Electronic supplementary material

The online version of this article (doi:10.1186/s12915-015-0186-5) contains supplementary material, which is available to authorized users.

The selenium content of SEPP1 versus selenium requirements in vertebrates

Selenoprotein P (SEPP1) distributes selenium (Se) throughout the body via the circulatory system. For vertebrates, the Se content of SEPP1 varies from 7 to 18 Se atoms depending on the species, but the reason for this variation remains unclear. Herein we provide evidence that vertebrate SEPP1 Sec content correlates positively with Se requirements. As the Se content of full length SEPP1 is genetically determined, this presents a unique case where a nutrient requirement can be predicted based on genomic sequence information.

Translational arrest due to cytoplasmic redox stress delays adaptation to growth on methanol and heterologous protein expression in a typical fed-batch culture of Pichia pastoris

Results

We have followed a typical fed-batch induction regime for heterologous protein production under the control of the AOX1 promoter using both microarray and metabolomic analysis. The genetic constructs involved 1 and 3 copies of the TRY1 gene, encoding human trypsinogen. In small-scale laboratory cultures, expression of the 3 copy-number construct induced the unfolded protein response (UPR) sufficiently that titres of extracellular trypsinogen were lower in the 3-copy construct than with the 1-copy construct. In the fed-batch-culture, a similar pattern was observed, with higher expression from the 1-copy construct, but in this case there was no significant induction of UPR with the 3-copy strain. Analysis of the microarray and metabolomic information indicates that the 3-copy strain was undergoing cytoplasmic redox stress at the point of induction with methanol. In this Crabtree-negative yeast, this redox stress appeared to delay the adaptation to growth on methanol and supressed heterologous protein production, probably due to a block in translation.

Conclusion

Although redox imbalance as a result of artificially imposed hypoxia has previously been described, this is the first time that it has been characterised as a result of a transient metabolic imbalance and shown to involve a stress response which can lead to translational arrest. Without detailed analysis of the underlying processes it could easily have been mis-interpreted as secretion stress, transmitted through the UPR.

Pichia pastoris regulates its gene-specific response to different carbon sources at the transcriptional, rather than the translational, level

Background

The methylotrophic, Crabtree-negative yeast Pichia pastoris is widely used as a heterologous protein production host. Strong inducible promoters derived from methanol utilization genes or constitutive glycolytic promoters are typically used to drive gene expression. Notably, genes involved in methanol utilization are not only repressed by the presence of glucose, but also by glycerol. This unusual regulatory behavior prompted us to study the regulation of carbon substrate utilization in different bioprocess conditions on a genome wide scale.

Results

We performed microarray analysis on the total mRNA population as well as mRNA that had been fractionated according to ribosome occupancy. Translationally quiescent mRNAs were defined as being associated with single ribosomes (monosomes) and highly-translated mRNAs with multiple ribosomes (polysomes). We found that despite their lower growth rates, global translation was most active in methanol-grown P. pastoris cells, followed by excess glycerol- or glucose-grown cells. Transcript-specific translational responses were found to be minimal, while extensive transcriptional regulation was observed for cells grown on different carbon sources. Due to their respiratory metabolism, cells grown in excess glucose or glycerol had very similar expression profiles. Genes subject to glucose repression were mainly involved in the metabolism of alternative carbon sources including the control of glycerol uptake and metabolism. Peroxisomal and methanol utilization genes were confirmed to be subject to carbon substrate repression in excess glucose or glycerol, but were found to be strongly de-repressed in limiting glucose-conditions (as are often applied in fed batch cultivations) in addition to induction by methanol.

Conclusions

P. pastoris cells grown in excess glycerol or glucose have similar transcript profiles in contrast to S. cerevisiae cells, in which the transcriptional response to these carbon sources is very different. The main response to different growth conditions in P. pastoris is transcriptional; translational regulation was not transcript-specific. The high proportion of mRNAs associated with polysomes in methanol-grown cells is a major finding of this study; it reveals that high productivity during methanol induction is directly linked to the growth condition and not only to promoter strength.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1393-8) contains supplementary material, which is available to authorized users.

Development of chromosome-specific markers with high polymorphism for allotetraploid cotton based on genome-wide characterization of simple sequence repeats in diploid cottons (Gossypium arboreum L. and Gossypium raimondii Ulbrich)

Background

Tetraploid cotton contains two sets of homologous chromosomes, the At- and Dt-subgenomes. Consequently, many markers in cotton were mapped to multiple positions during linkage genetic map construction, posing a challenge to anchoring linkage groups and mapping economically-important genes to particular chromosomes. Chromosome-specific markers could solve this problem. Recently, the genomes of two diploid species were sequenced whose progenitors were putative contributors of the At- and Dt-subgenomes to tetraploid cotton. These sequences provide a powerful tool for developing chromosome-specific markers given the high level of synteny among tetraploid and diploid cotton genomes. In this study, simple sequence repeats (SSRs) on each chromosome in the two diploid genomes were characterized. Chromosome-specific SSRs were developed by comparative analysis and proved to distinguish chromosomes.

Results

A total of 200,744 and 142,409 SSRs were detected on the 13 chromosomes of Gossypium arboreum L. and Gossypium raimondii Ulbrich, respectively. Chromosome-specific SSRs were obtained by comparing SSR flanking sequences from each chromosome with those from the other 25 chromosomes. The average was 7,996 per chromosome. To confirm their chromosome specificity, these SSRs were used to distinguish two homologous chromosomes in tetraploid cotton through linkage group construction. The chromosome-specific SSRs and previously-reported chromosome markers were grouped together, and no marker mapped to another homologous chromosome, proving that the chromosome-specific SSRs were unique and could distinguish homologous chromosomes in tetraploid cotton. Because longer dinucleotide AT-rich repeats were the most polymorphic in previous reports, the SSRs on each chromosome were sorted by motif type and repeat length for convenient selection. The primer sequences of all chromosome-specific SSRs were also made publicly available.

Conclusion

Chromosome-specific SSRs are efficient tools for chromosome identification by anchoring linkage groups to particular chromosomes during genetic mapping and are especially useful in mapping of qualitative-trait genes or quantitative trait loci with just a few markers. The SSRs reported here will facilitate a number of genetic and genomic studies in cotton, including construction of high-density genetic maps, positional gene cloning, fingerprinting, and genetic diversity and comparative evolutionary analyses among Gossypium species.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1265-2) contains supplementary material, which is available to authorized users.

Effect of nitrogen source and inorganic phosphate concentration on methanol utilization and PEX genes expression in Pichia pastoris

Methylotrophic yeast Pichia pastoris has proved to be especially useful for production of various heterologous proteins. In biotechnology it is very important to maintain the balance between high levels of heterologous gene expression and cell viability. Decisive understanding of gene regulation mechanisms is essential for reaching this goal. In this study, we investigated the effect of different nitrogen sources and phosphate concentration in media on methanol utilization. It was shown that expression levels of main genes, which are involved in methanol utilization (MUT genes) and in functioning of peroxisomes (PEX genes), are maximal when ammonium sulphate is used as a nitrogen source. Expression of these genes is decreased in media with poor nitrogen sources, such as proline. Addition of rapamycin to the media completely removed repression of AOX1 promoter in media with proline, which allows proposing that Tor-kinase is involved in establishing of nitrogen regulation of this gene. It was also shown that MUT genes expression levels get higher, when the phosphate concentration in media is increased.

Mating-type switching by chromosomal inversion in methylotrophic yeasts suggests an origin for the three-locus Saccharomyces cerevisiae system

Saccharomyces cerevisiae has a complex system for switching the mating type of haploid cells, requiring the genome to have three mating-type (MAT)-like loci and a mechanism for silencing two of them. How this system originated is unknown, because the three-locus system is present throughout the family Saccharomycetaceae, whereas species in the sister Candida clade have only one locus and do not switch. Here we show that yeasts in a third clade, the methylotrophs, have a simpler two-locus switching system based on reversible inversion of a section of chromosome with MATa genes at one end and MATalpha genes at the other end. In Hansenula polymorpha the 19-kb invertible region lies beside a centromere so that, depending on the orientation, either MATa or MATalpha is silenced by centromeric chromatin. In Pichia pastoris, the orientation of a 138-kb invertible region puts either MATa or MATalpha beside a telomere and represses transcription of MATa2 or MATalpha2. Both species are homothallic, and inversion of their MAT regions can be induced by crossing two strains of the same mating type. The three-locus system of S. cerevisiae, which uses a nonconservative mechanism to replace DNA at MAT, likely evolved from a conservative two-locus system that swapped genes between expression and nonexpression sites by inversion. The increasing complexity of the switching apparatus, with three loci, donor bias, and cell lineage tracking, can be explained by continuous selection to increase sporulation ability in young colonies. Our results provide an evolutionary context for the diversity of switching and silencing mechanisms.

Characterization of natural antisense transcript, sclerotia development and secondary metabolism by strand-specific RNA sequencing of Aspergillus flavus

Aspergillus flavus has received much attention owing to its severe impact on agriculture and fermented products induced by aflatoxin. Sclerotia morphogenesis is an important process related to A. flavus reproduction and aflatoxin biosynthesis. In order to obtain an extensive transcriptome profile of A. flavus and provide a comprehensive understanding of these physiological processes, the isolated mRNA of A. flavus CA43 cultures was subjected to high-throughput strand-specific RNA sequencing (ssRNA-seq). Our ssRNA-seq data profiled widespread transcription across the A. flavus genome, quantified vast transcripts (73% of total genes) and annotated precise transcript structures, including untranslated regions, upstream open reading frames (ORFs), alternative splicing variants and novel transcripts. We propose natural antisense transcripts in A. flavus might regulate gene expression mainly on the post-transcriptional level. This regulation might be relevant to tune biological processes such as aflatoxin biosynthesis and sclerotia development. Gene Ontology annotation of differentially expressed genes between the mycelia and sclerotia cultures indicated sclerotia development was related closely to A. flavus reproduction. Additionally, we have established the transcriptional profile of aflatoxin biosynthesis and its regulation model. We identified potential genes linking sclerotia development and aflatoxin biosynthesis. These genes could be used as targets for controlled regulation of aflatoxigenic strains of A. flavus.

GC-rich DNA elements enable replication origin activity in the methylotrophic yeast Pichia pastoris

The well-studied DNA replication origins of the model budding and fission yeasts are A/T-rich elements. However, unlike their yeast counterparts, both plant and metazoan origins are G/C-rich and are associated with transcription start sites. Here we show that an industrially important methylotrophic budding yeast, Pichia pastoris, simultaneously employs at least two types of replication origins--a G/C-rich type associated with transcription start sites and an A/T-rich type more reminiscent of typical budding and fission yeast origins. We used a suite of massively parallel sequencing tools to map and dissect P. pastoris origins comprehensively, to measure their replication dynamics, and to assay the global positioning of nucleosomes across the genome. Our results suggest that some functional overlap exists between promoter sequences and G/C-rich replication origins in P. pastoris and imply an evolutionary bifurcation of the modes of replication initiation.

Genome sequence and analysis of methylotrophic yeast Hansenula polymorpha DL1

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

Identification and characterization of P GCW14 : a novel, strong constitutive promoter of Pichia pastoris

The available promoters in the Pichia pastoris expression platform are still limited. We selected and identified a novel strong constitutive promoter, P GCW14 , and tested its promoter activity using enhanced green fluorescent protein (EGFP) as a reporter. Potential promoter regions of P GCW14 were cloned upstream of the EGFP gene and promoter activity was analyzed by measuring fluorescence intensity. P GCW14 exhibited significantly stronger promoter activity than the classic strong constitutive promoters P TEF1 and P GAP under various carbon sources, suggesting that P GCW14 is a strong and constitutive promoter. Hence, P GCW14 can be used as a promoter for high-level expression of heterologous proteins.