Abstract
Non-homologous end joining (NHEJ)-mediated integration is effective in generating random mutagenesis to identify beneficial gene targets in the whole genome, which can significantly promote the performance of the strains. Here, a novel target leading to higher protein synthesis was identified by NHEJ-mediated integration that seriously improved fatty alcohols biosynthesis in Yarrowia lipolytica. One batch of strains transformed with fatty acyl-CoA reductase gene (FAR) showed significant differences (up to 70.53-fold) in fatty alcohol production. Whole-genome sequencing of the high-yield strain demonstrated that a new target YALI0_A00913g (“A1 gene”) was disrupted by NHEJ-mediated integration of partial carrier DNA, and reverse engineering of the A1 gene disruption (YlΔA1-FAR) recovered the fatty alcohol overproduction phenotype. Transcriptome analysis of YlΔA1-FAR strain revealed A1 disruption led to strengthened protein synthesis process that was confirmed by sfGFP gene expression, which may account for enhanced cell viability and improved biosynthesis of fatty alcohols. This study identified a novel target that facilitated synthesis capacity and provided new insights into unlocking biosynthetic potential for future genetic engineering in Y. lipolytica.
Introduction
The metabolic reactions and regulatory mechanisms that exist in organisms are not performed in isolation, thereby microbial biosynthesis of a desired product is constantly in the control of a complex intracellular network (Herrgard et al., 2006; Chen et al., 2018). In addition to the synthetic pathway, an abundance of beneficial gene targets in the genome can promote the performance of strains, such as production yield (Wang et al., 2018), toxicity tolerance (Hu et al., 2021; Li et al., 2021) or cell viability (Shiomi et al., 2013; Jiang et al., 2015). However, identification of potential targets of unknown function by rational engineering remains difficult because of our limited understanding of the underlying linkages among the metabolic processes (Ajjawi et al., 2017; Fang et al., 2021). In contrast, random mutagenesis succeeds in inducing arbitrary mutations with random distribution to produce genotypic diversity, which is an efficient approach to identify new targets and allow further unravel the inner machinery of metabolism by mapping the genotypes and phenotypes (van Opijnen et al., 2009).
Non-homologous end joining (NHEJ) is a typical pathway to repair genomic DNA double-strand breaks (DSBs; Lieber, 2010; Ranjha et al., 2018). NHEJ is active throughout the cell cycle and is considered the predominant DSB repair pathway (Bouwman and Crosetto, 2018). In human cells, NHEJ appears to repair nearly all DSBs outside of S and G2 cell cycle phases and even about 80% of DSBs within S and G2 (Beucher et al., 2009). NHEJ-mediated integration can efficiently generate random mutagenesis across the genome in Yarrowia lipolytica. Firstly, NHEJ repairs DSBs is recognized as error-prone, which tends to result in gene mutations and is rarely restored to their original DNA sequence (Burma et al., 2006; Shrivastav et al., 2008; Pannunzio et al., 2018). Secondly, NHEJ-mediated integration confers the randomness and genome-scale unbiased coverage of gene mutations in Y. lipolytica. Liu et al. revealed that 9.15 × 105 NHEJ-mediated independent integration events were relatively uniformly distributed across the genome among approximately one million colonies (Liu et al., 2022). Finally, NHEJ-mediated integration is time-saving and convenient to generate gene mutations. Compared with other mutagenesis approaches, NHEJ-mediated integration requires only linear DNA amplification without the assistance of heterologous transposons (Kumar, 2016) or a CRISPR system (Schwartz et al., 2019) and therefore the mutations are generated easily and rapidly. Up to 65 copy numbers were achieved within 96 h by Bai et al. through optimized NHEJ-mediated genome integration technology (Bai et al., 2021). Based on these advantages, NHEJ-mediated integration was utilized to screen causal mutations related to interested phenotypes. For example, 7 new targets for improving β-carotene biosynthesis or acetic acid tolerance were identified by NHEJ-mediated insertional mutagenesis in Y. lipolytica (Liu et al., 2022), which contributed to superior cell factory construction.
Yarrowia lipolytica displays a strong preference for repairing DSBs by NHEJ other than homologous recombination (HR; Kretzschmar et al., 2013; Verbeke et al., 2013). Besides the complete foreign sequences, even the fractional fragments or vector can also be spontaneously inserted into the genome, resulting in various genotypes in one batch of gene transformation (Liu et al., 2019; Longmuir et al., 2019). However, the potentially beneficial gene targets might be missed, if the researchers only select strains with phenotypes of interest without further analysis of random insertion sites (Ye et al., 2012). For example, if the whole-genome sequencing had not been performed on a high-yield eicosapentaenoic acid producing strain generated by NHEJ-mediated random insertion of pathway genes, mutation of PEX10 gene might have been overlooked (Xue et al., 2013), which was proved to be highly effective and was widely engineered to enhance lipids and terpenoids production in later works (Blazeck et al., 2014; Jin et al., 2019; Wei et al., 2021). Therefore, the incidental phenotypic differences should be valued and further investigated, which might contribute to the identification of potential targets and greatly expand the understanding of biological knowledge.
In this study, we spotted unexpected strains with fatty alcohols overproduction in the same batch after transforming the fatty acyl-CoA reductase (FAR) gene. Then, the whole-genome sequencing of the high-yield strains showed that the YALI0_A00913g (“A1 gene”) was disrupted due to the partial carrier DNA sequence was inserted into the open reading frame (ORF) by NHEJ-mediated integration. Reverse engineering of the A1 gene disruption (YlΔA1-FAR) recovered the high-yielding fatty alcohols phenotype. Transcriptome analysis of YlΔA1-FAR strain revealed that disruption of A1 gene resulted in a significant enhancement of protein synthesis, which was confirmed by characterization with green fluorescence intensity (sfGFP gene). Therefore, the disruption of new target A1 boosted the cell viability and accounted for the increased fatty alcohol production. Target identification by NHEJ-mediated integration facilitated the unlocking of biosynthetic potential and effectively increased the productivity of Y. lipolytica.
Materials and Methods
Strains, Medium and Culture Conditions
Yarrowia lipolytica strain ATCC 201249 (MATA, ura3-302, leu2-270, lys8-11, PEX17-HA) was selected as the background strain for all constructs (Gao et al., 2016), and the initial strain ATCC 201249 was provided by Professor Ying-Jin Yuan (Tianjin University, China). All derivatives constructed in this study are listed in Supplementary Table 1. The above strains of Y. lipolytica were cultured in yeast culture medium at 28°C and 250 rpm shaking speed.
TransT1 E. coli was used for plasmid construction and propagation, which was grown in Luria-Bertani broth (LB) at 37°C with constant shaking at 220 rpm. When necessary, the selective antibiotic was added (100 mg/l ampicillin or 50 mg/l kanamycin). Agar (15 g/l) was added for LB solid plate preparation. Yeast extract-peptone-dextrose (YPD) medium, consisting of 20 g/l glucose, 20 g/l peptone, and 10 g/l yeast extract, was used for transformation, activation and preculture. Agar (20 g/l) was added for YPD solid plate preparation. The rich YPD medium containing 50 g/l glucose, 20 g/l peptone, and 10 g/l yeast extract was used for fermentation of Y. lipolytica strains. The SC medium, used for screening Y. lipolytica transformants, contained 20 g/l glucose, 6.7 g/l yeast nitrogen base without amino acids, and 2 g/l complete supplement mixture (CSM) lacking uracil (SC-Ura) or leucine (SC-Leu), supplemented with uracil or leucine depending on the selection marker requirements.
For fatty alcohols flask fermentation, freshly streaked single colonies of strains were first cultivated in a 25-mL polypropylene tube with 5 ml YPD medium, and cultivate overnight at 28°C with a shaking speed of 250 rpm. After preculturing, the seed cultures were inoculated into 50 ml of fresh YPD medium with an initial OD600 of 0.2 and fermented in a 250-mL shake flask for 72 h or extend appropriately.
Construction of Plasmids and Strains
Plasmids used in this work are listed in Supplementary Table 2. Fatty acyl-CoA reductase gene from Marinobacter aquaeolei VT8 (i.e., MaFAR) and sfGFP were codon optimized and synthesized by GenScript (Nanjing, China) and listed in Supplementary Data 1. The IntX HUM were amplified from genomic DNA of Y. lipolytica and constructed the integrated plasmids pIntA, pIntB, pIntC, pIntD and pIntE through seamless cloning. And the expressed gene were amplified with the proper sticky ends, cleaved with DNA endonuclease and integrated into the corresponding integration plasmids by T4 DNA ligase. All the primers used to construct the expression cassettes and the plasmids were synthesized by Genewiz Inc. (China) and are listed in Supplementary Table 3. The genotypes of all strains were confirmed by colony polymerase chain reaction of KOD FX DNA polymerase (Toyobo Co., Ltd.; Shanghai, China).
For disrupting target genes in Y. lipolytica, the CRISPR system was used as previously described (Schwartz et al., 2016; Zhang et al., 2018). The backbone plasmid used for constructing CRISPR-Cas9 plasmids for gene knockout in Y. lipolytica was pMCS-Cen1. The synthesized gRNA was incorporated into pMCS-URA via one-step Golden Gate assembly. The corresponding plasmids were digested with restriction enzyme BamHI and HindIII, and ligated with the segment of Cas9 to form the final plasmids. Strains with Ura-marker plasmid were inoculated into YPD liquid medium and cultured at 28°C for 3–4 days, then they were screened on YPD solid medium containing 1.2 mg/ml 5-fluoroorotic acid after dilution for the loss of Ura-marker plasmid.
Transformation
Yarrowia lipolytica transformation with integrative fragments or recombinant plasmids was performed with Zymogen Frozen EZ Yeast Transformation Kit II (Zymo Research Corporation). For CRISPR plasmid transformations, the cells were transformed with plasmid and cultivated in SC-Ura liquid medium for 4 d. Then the cells were plated onto SC-Ura plates for 2 d and confirmed by sequencing. For integrative fragment transformations, approximately 2 μg of linearized DNA was used in the transformation reaction and then the cells were harvested by centrifugation at 5,000 rpm for 2 min and plated on SC agar plates without the auxotrophic compound supplemented by the fragments. Selection plates were incubated at 28°C for 2–3 days.
Extraction and Analysis of Fatty Alcohols
Extraction of fatty alcohols was similar with the previously described method (Xu et al., 2016). Briefly, cells from 0.5 ml of culture were blended with 1-Nonadecanol acid (the final concentration is about 50~500 mg/l, adjusted according to the actual output) as internal standard (IS) and extracted with 0.5 ml of ethyl acetate. The ethyl acetate extract was vortexed at a speed of 1,200 r/min for 20 min, and then centrifuged at 15000 rpm for 25 min. The supernatant organic phase was filtered through a 0.22 μm filter membrane for fatty alcohols assay using gas chromatography (GC). The extracted fatty alcohol was then determined using a Thermo Scientific TRACE 1300 GC equipped with a TG-5MS column (30 m × 0.32 mm × 0.25 μm; Thermo Scientific) and a Flame Ionization Detector (FID) operating under constant flow rate of the carrier gas (nitrogen) at 1 ml min−1. The injector temperature was 260°C; it was operated on constant pressure mode at 91 kPa. The column temperature was held at 70°C for 2 min and then increased to 290°C with the rate of 8°C/min, holding for 6 min. Individual fatty alcohol species were qualified by authentic homologous standards and quantified by comparing the peak areas with that of the internal standard using the Chromeleon 7.1 software.
Fluorescence Assay
The Y. lipolytica strains transformed with sfGFP gene were activated in Sc-Leu medium for 24 h, and then the seed cultures were inoculated into 25-mL polypropylene tube containing 5 ml fresh SC-Leu medium with an initial OD600 of 0.2. The strains were grown at 250 rpm for 48 h at 28°C. 500 ul suspensions of each polypropylene tubes were centrifuged at 5000 rpm for 2 min to remove the supernatant and the cells were washed and resuspended with the phosphate-buffered saline (PBS) buffer. Fluorescence intensity (excitation: 488 nm and emission: 530 nm) was measured using a 96-well polystyrene plates (black plate, clear bottom with lid; Corning Incorporated 3,603, United States) after dilution into the linear range of the detector by a multimode microplate reader (Tecan Infinite 200 PRO, Austria). The fluorescence intensity was normalized to cell density (OD600), which was measured using the same microplate reader. The images of sfGFP fluorescence were observed by fluorescence microscope (Olympus CX41, Japan). Experiments were performed in three biological replicates.
Whole-Genome Sequencing
Whole-genome sequencing of Y. lipolytica strains was performed with an Illumina Hiseq2000 by the Genomic Sequencing and Analysis Facility at the Beijing Genome Institute (BGI, China) using paired-End sequencing. Genomic DNA of each individual is extracted by CTAB method (Murray and Thompson, 1980), and then fragmented randomly. After electrophoresis, DNA fragments of desired length are gel purified. Adapter ligation and DNA cluster preparation are performed and subjected to sequencing. Step one, the sequencing reads are aligned onto the reference CLIB122 genome sequence using SOAPaligner,1 obtain the single reads on the alignment, and perform pre-processing on this part of the reads. The sequencing depth and coverage compared to the reference genome are calculated based on the alignment. Step two, align the sequencing reads with the insert sequence as a reference sequence, and obtain the single reads aligning with the insert sequence. Extract the intersection of the two steps, count the distribution of reads in the genome, and determine the insertion site. SOAPsnp,2 SOAPindel and SOAPsv were used to analyze the data.
RNA-Seq Transcriptomic Analysis
Yarrowia lipolytica cells cultured at 24 h on the initial logarithmic phase were collected for RNA isolation. Cells were collected by centrifugation at 5000 rpm for 2 min, washed with PBS three times and then quick-frozen by liquid nitrogen and stored at − 80°C tentatively. The RNA sample preparation, library construction and sequencing were performed at the Beijing Genome Institute (BGI, China). RNA was isolated from the cell pellet using TRIzol™ Reagent (Invitrogen, United States), and the quality was determined using Bioanalyzer 2,100 (Agilent, United States). The RNA was sheared and reverse transcribed using random primers to obtain cDNA which was used for library construction. Then, the library quality was determined using a Bioanalyzer 2100 (Agilent) analyzer, and the libraries were sequenced on the BGISEQ-500 (BGI, China) sequencing platform. All raw sequencing reads were filtered using BGI’s own filtering software SOAPnuke to remove reads containing low quality, splice contamination and unknown base N content greater than 5%. The filtered “Clean Reads” were saved in FASTQ format for data analysis to ensure reliable results. Gene expression levels were calculated for each sample using the RSEM (Li and Dewey, 2011; Langmead and Salzberg, 2012) software package. NOI Seq (Anders and Huber, 2010; Tarazona et al., 2011) method for detecting differential gene expression between groups. Finally, in-depth data analysis was performed on the BGI data analysis platform.3
Determination of ATP Level
Cellular ATP levels were determined using an Enhanced ATP Assay Kit (Beyotime S0027, Nanjing, China). Equal numbers of Y. lipolytica cells cultured at 24 h on the initial logarithmic phase were collected and lysed by an ultrasonic cell-crushing device was used for cell disruption (COLE-PARMER INSTRUMENTS CP505, United States), which were centrifuged for 5 min at 4°C and 12,000 g, and the supernatant was collected. Subsequently, 100 ul of ATP detection solution was added to a 96-well plate (Corning Incorporated 3,603, United States) and then 20 ul of the supernatant was transferred into each well. Chemiluminescence was monitored using a multimode microplate reader (Tecan Infinite 200 PRO, Austria).
Results
Incidental Obtaining of Strains With High Fatty Alcohols Producing
Fatty alcohols have gained huge market and applied value due to their broad applications in the cosmetics, detergent, surfactant, and lubricant industries (Pfleger et al., 2015; Fillet and Adrio, 2016). It has been previously proven that fatty acyl-CoA reductase from M. aquaeolei VT8 (encoded by the MaFAR gene) was the most effective enzyme to directly catalyze the reduction of fatty acyl-CoA to fatty alcohols (Figure 1A; Wahlen et al., 2009; Willis et al., 2011; Wang et al., 2016b). To establish an exogenous synthesis pathway of fatty alcohol in Y. lipolytica, the MaFAR gene was expressed by pINA1269 vector in strains ATCC 201249, and Y10301 (ΔPEX10; Zhang et al., 2019), generating strains YlA and YlP (Figures 1B,C). pINA1269 is a pBR322-based mono-copy integrative vector carrying the LEU2 gene, which is supposed to integrate into the PBR docking platform of the genome through homologous recombination (Xuan et al., 1990; Madzak et al., 2000; Nicaud et al., 2002). We picked a series of strains randomly on each plate for fermentation; however, the results were amazing. Strains with unexpected yields were found incidentally and the yield of other strains was at the general level. As shown in Figure 1D, strains from the same plate have significantly different yields and we selected one of the low-yielding strains as the control. Strain YlA-1 produced 255.3 mg/l of fatty alcohols, which was 70.53-fold higher than strain YlA-2 (3.62 mg/l), both two strains based on ATCC 201249. Similarly, YlP-1 produced 49.94-fold fatty alcohols than YlP-2. In theory, the strains obtained from the same transformation batch should not have significant differences in fatty alcohols production. But NHEJ is the main pathway for repairing DNA double strand breaks (DSB) of genome in Y. lipolytica (Richard et al., 2005; Verbeke et al., 2013), which may produce unexpected NHEJ-mediated insertions accompanied by integration of the pINA1269-FAR plasmid to generate new genotypes (Gao et al., 2014; Cui et al., 2019; Sudfeld et al., 2021).
Figure 1
Whole-Genome Sequencing Suggested Gene Disruption by NHEJ-Mediated Integration
To further investigate the potential cause of improved fatty alcohols production, we performed whole-genome sequencing of the YlA-1,2 and YlP-1,2 to confirm the genotypes (Liu et al., 2015a,b; Li and Alper, 2016). As shown in Figure 2, the FAR gene was integrated into the pBR322 locus on chromosome F as expected by comparison with the draft genome sequence of the wild-type control strain, ATCC 201249. However, additional changes of gene disruption were indeed discovered in genome of high-yield strains, which were caused by the integration of a small proportion of vector into chromosomes A and E, named A1 and E1 sites, respectively. The specific locations of disruptions are shown in Figure 2; the A1 site was located in the 118,335–118,543 region of chromosomes A, which was in the open reading frame in YALI0_A00913g gene (“A1 gene”). Similarly, the E1 site was in the 4,134,688–4,134,892 region of chromosomes E in YALI0_E34687g gene (“E1 gene”). However, these disruptions were not available in low-yield strains of YlA-2 and YlP-2 (Supplementary Figure 1; Supplementary Table 4), and we hypothesized that these additional changes led to high yield of fatty alcohols.
Figure 2
Reverse-Engineered Disruption Validation of New Target YALI0_A00913g (A1 Gene)
We further conducted reverse engineering on these targets to validate whether the disruption of unknown genes A1 and E1 can enhance the fatty alcohols production. As shown in Figure 3A, the A1 and E1 genes were deleted 2 bp (GG) and 11 bp (GAAGAGGAGTA), respectively, on the basis of the blank strain ATCC 201249 through the CRISPR/Cas system (Gao et al., 2016; Schwartz et al., 2016), obtaining the mutant strains YlΔA1 and YlΔE1. The position we selected to knock out was consistent with the site of the broken vector fragments insertion in the previous whole-genome sequencing. To reduce the occurrence of the incidental insertion by NHEJ, we selected a relatively stable pIntF vector contained 500 bp homologous to IntF (YALI0F3161413 to YALI0F3162449) flanking a gene expression cassette (Matthaus et al., 2014; Holkenbrink et al., 2018; Figure 3B). The FAR gene was expressed under the control of TEFin promoter by pIntF vector and integrated into the IntF site in strains ATCC 201249, YlΔA1 and YlΔE1, respectively. Three strains were generated, named WT-FAR, YlΔA1-FAR and YlΔE1-FAR. Fatty alcohols production for 72 h of fermentation was shown in Figure 3C, the yield of strain WT-FAR without knockout target was only 2.47 mg/l and the strain YlΔE1-FAR that knocked out the E1 target produced 4.50 mg/l of fatty alcohols. Previous studies have shown that the yield of Y. lipolytica strains expressing the MaFAR gene is only 1.5–80 mg/l (Xu et al., 2016; Wang et al., 2016a; Zhang et al., 2019), so the low-yield strains are conventional. While the yield of strain YlΔA1-FAR with knockout A1 target is as high as 810.39 mg/l (Fermentation for 90 h up to 851.8 mg/l, Supplementary Figure 2), which was even higher than the yield of previous strains optimized for metabolism (Wang et al., 2016a; Xu et al., 2016; Zhou et al., 2016), which is unconventional. It was confirmed that the disruption of A1 gene can significantly increase the yield of fatty alcohols, but the disruption of E1 gene did not reach the expected yield. Since the high fatty alcohol producing strain YlP-1 was based on Yl0301 (ΔPEX10), but we deleted the target gene E1 in ATCC 201249 instead of Y10301. This may have resulted in the YlΔE1-FAR not reaching the expected yield.
Figure 3
Transcriptome Analysis Revealed the Increased Protein Synthesis Owing to A1 Disruption
The hypothetical gene A1 (YALI0_A00913g) was weakly similar (31.4%) to the YND1 gene of Saccharomyces cerevisiae (D’Alessio et al., 2005; Maoz et al., 2005; Onuma et al., 2006; Mittelman et al., 2010) after alignment by Protein BLAST (Supplementary Figure 3). The YND1 gene encodes Golgi apyrase that is not directly related to the fatty alcohol biosynthesis pathway. In order to reveal the possible mechanism of high yield caused by the disruption of A1 gene, RNA-Seq was performed to analyze the transcriptome differences between YlΔA1-FAR and WT-FAR. As shown in Figure 4A, cell samples were taken at 24 h which was the initial logarithmic phase. Transcriptome analysis results showed that the disruption of A1 gene and excess accumulation of fatty alcohols resulted in large transcriptional changes. A total of 437 DEGs were identified in YlΔA1-FAR compared with the control, with 262 and 175 genes being upregulated and downregulated, respectively (Figure 4B). The DEGs were subjected to KEGG enrichment analysis, which were enriched in 16 pathways. Among them, the folding, sorting and degradation and translation noted by asterisks (*) are the top two gene sets with the highest proportions of upregulated genes, reaching 88.9 and 80.0%, respectively (Figure 4C). The two gene sets were obviously related to protein synthesis and processing, and the genes with significant differences in transcription level are shown in Table 1.
Figure 4
Table 1
| Functional group and gene | log2-Fold | Description |
|---|---|---|
| Translation | ||
| RNA transport | ||
| YALI0_C17567g | 2.56 | NUP98, ADAR2; nuclear pore complex protein Nup98-Nup96 (Ho et al., 2000; Krull et al., 2004) |
| YALI0_D20658g | 1.30 | UPF2, RENT2; regulator of nonsense transcripts 2 (Chamieh et al., 2008) |
| Ribosome | ||
| YALI0_B13200g | 2.66 | EIF3B; translation initiation factor 3 subunit B (Phan et al., 2001) |
| YALI0_F25399g | 1.62 | RP-S15Ae, RPS15A; small subunit ribosomal protein S15Ae (Vladimirov et al., 1996; Ikeda et al., 2017) |
| YALI0_A02497g | 1.44 | EIF4G; translation initiation factor 4G (Levy-Strumpf et al., 1997; Gingras et al., 1999) |
| YALI0_C18975g | 1.04 | EIF4G; translation initiation factor 4G (Levy-Strumpf et al., 1997; Gingras et al., 1999) |
| Aminoacyl-tRNA biosynthesis | ||
| YALI0_B13134g | 3.18 | RARS, argS; arginyl-tRNA synthetase (Mehler and Mitra, 1967; Cavarelli et al., 1998; Delagoutte et al., 2000) |
| YALI0_C15246g | 1.85 | LARS, leuS; leucyl-tRNA synthetase (Lincecum et al., 2003) |
| Folding, sorting and degradation | ||
| Protein translocation | ||
| YALI0_D08635g | 1.54 | SEC61G, SSS1, secE; protein transport protein SEC61 subunit gamma and related proteins (Falcone et al., 2011) |
| Protein processing | ||
| YALI0_E35046g | 1.74 | HSPA1s; heat shock 70 kDa protein 1/2/6/8 (Hunt and Morimoto, 1985; Mayer and Bukau, 2005; Perez-Vargas et al., 2006) |
| YALI0_E18546g | 1.64 | HSP20; HSP20 family protein (Sun et al., 2001; Liu et al., 2012) |
| YALI0_C03465g | 1.17 | HSP20; HSP20 family protein (Sun et al., 2001; Liu et al., 2012) |
| Protein degradation | ||
| YALI0_C05709g | −1.26 | TRIP12; E3 ubiquitin-protein ligase TRIP12 (Park et al., 2009) |
| YALI0_A19426g | −1.16 | PSMA2; 20S proteasome subunit alpha 2 (Coux et al., 1996; Voges et al., 1999; Smith et al., 2006) |
The differentially expressed genes in the folding, sorting and degradation and translation gene sets.
In the translation gene set, 8 genes were significantly upregulated, belonging to RNA transport, ribosome, aminoacyl-tRNA synthesis pathway. First, in the RNA transport pathway, 2 genes YALI0_C17567g and YALI0_D20658g were significantly upregulated by2.56-log2 and 1.30-log2 fold in YlΔA1-FAR strain, respectively. Among them, YALI0_C17567g gene encodes the nuclear pore complex (NPC) protein, which mediates the export of RNAs produced in the nucleus to fulfill their function in protein synthesis (Rout et al., 2000; Rout and Aitchison, 2001; Rodriguez et al., 2004). Therefore, the upregulation of YALI0_C17567g gene can facilitate the transport of RNAs out of the nucleus (Figure 5A). Then, in the ribosome pathway, 4 genes were upregulated in YlΔA1-FAR strain, there (YALI0_B13200g, YALI0_A02497g and YALI0_C18975g) of which are encoding the eukaryotic translation initiation factors (eIFs) family. The initiation of protein synthesis in eukaryotic cells is dependent on multiple eIFs and initiation is the rate-limiting step for translation under most circumstances (Pain and Standart, 2001). Upregulation of eIFs-related genes stimulated the binding of mRNA and methionyl-initiator tRNA (Met-tRNAiMet) to the ribosome (Gingras et al., 1999; Fromont-Racine et al., 2003; Valasek, 2012), thereby could enhance the translation rates (Figure 5B). Last, two genes YALI0_B13134g and YALI0_C15246g related to the aminoacyl-tRNA biosynthesis pathway were upregulated by 3.18-log2 and 1.85-log2 fold, respectively. They encode aminoacyl-tRNA synthetases (ARS), which catalyze the ligation of amino acids to their cognate tRNAs (Negrutskii and Deutscher, 1991; Cusack, 1997; Park et al., 2008; Ling et al., 2009). The upregulation of related genes facilitated the connection between amino acids and tRNA, thus promoting protein synthesis (Figure 5C).
Figure 5
In the folding, sorting and degradation gene set, genes with significant differences in transcription were mainly distributed in protein translocation, protein processing, and protein degradation pathway. First in the protein translocation pathway, the gene YALI0_D08635g encoding protein transport protein Sec61 subunit gamma was upregulation by 1.54-log2 fold (Falcone et al., 2011). The Sec61 complex mediated protein precursors crossing the endoplasmic reticulum (ER) membrane of eukaryotes, which is a vital, first committed step in the biogenesis of many proteins (Stephenson, 2005; Robson and Collinson, 2006; Rapoport, 2007). Therefore, the upregulation of YALI0_D08635g gene facilitated the precursor proteins to across the ER membrane for subsequent protein processing (Figure 5D). Then, in the protein processing pathway, 3 genes were significantly upregulated, they all belong to the heat shock protein (Hsp) family, which is a kind of molecular chaperone. The Hsp family chaperones prevent protein aggregation and catalyze polypeptide folding in the protein processing (Mayer and Bukau, 2005; Young, 2010; Ghazaei, 2017) and drive post-translational protein translocation to yeast ER (McClellan et al., 1998; Rapoport, 2007; Dudek et al., 2009). The upregulation of related genes promoted the binding of molecular chaperone to hydrophobic patches on proteins (Figure 5E). Finally, in the protein degradation pathway, the YALI0_C05709g gene encoding E3 ubiquitin-protein ligase and the YALI0_A19426g gene encoding 20S proteasome subunit alpha 2 were downregulation by 1.26-log2 and 1.16-log2 fold, respectively. E3s facilitates covalent attachment of ubiquitin to the substrate, then the ubiquitinated protein is recognized and degraded by the 26S proteasome (Johnson et al., 1995; Ciechanover et al., 2000; Park et al., 2009). Thus, downregulation of these genes might reduce protein degradation (Figure 5F).
Above all, transcriptome analysis results suggested that the protein synthesis was significantly enhanced, which potentially contributed to superior cell viability. Our results indicated that the growth of strain YlΔA1-FAR was indeed elevated (Figure 4A) and fatty alcohols production was improved (Figure 3C) compared to the control WT-FAR. As described above, the A1 gene is weakly similar to the YND1 gene of S. cerevisiae, which encodes Golgi apyrase that generally hydrolyzes both ATP and ADP (Supplementary Figure 4A). The transcript levels of the A1 gene before and after disruption are shown in Supplementary Figure 5, which indicated that the disruption did not completely disable the A1 gene but reduced expression level by 50%. The decreased expression of the A1 gene might reduce the consumption of ATP, which may facilitate the maintenance of energy (Choo et al., 2010) and provide energy for protein synthesis. We measured the intracellular ATP levels and the results are shown in Supplementary Figure 4B, the ATP content of YlΔA1-FAR was indeed higher (2.24-fold) than WT-FAR.
Verification of Improved Protein Expression by sfGFP Gene Expression
To further confirm the increased protein synthesis suggested by transcriptome analysis, we characterized protein expression before and after A1 gene disruption using superfolder green fluorescent protein (sfGFP; Pedelacq et al., 2006) as a reporter, which was work well in Y. lipolytica (Zhang et al., 2018; Bai et al., 2021). The sfGFP gene expressed by constitutive TEFin promoter in pIntF vector was integrated into the IntF site of the ATCC 201249 and YlΔA1 strains through HR (Figure 6A), generating strains of YGF01 to YGF02 (Figure 6B). Compared with the fluorescence intensity of the strain YGF01, the fluorescence intensity of the strain YGF02 that disrupts the A1 gene was increased by about 6.46 times. The fluorescence intensity of the cells was observed by fluorescence microscope, as shown in Figure 6C, the cell growth state of strain YGF02 was no different from that of strain YGF01, but the fluorescence intensity was significantly enhanced, which was consistent with the Figure 6B. These results indicated that the disruption of A1 gene truly increased protein expression, which further corroborates the transcriptome results. And the protein content of WT-FAR and YlΔA1-FAR strains was measured at 48 h, the results did support the above conclusion (Supplementary Table 5).
Figure 6
To investigate the regularity of enhanced exogenous protein expression, we next verified whether other promoters and integration sites (Figure 7A) were applicable to the increased protein expression in YlΔA1 strain. Firstly, we selected the IntX location (IntA to IntF; Matthaus et al., 2014) as the integration site, which was commonly used in Y. lipolytica. The sfGFP gene was integrated into the IntX location of strains ATCC 201249 and YlΔA1, the integration location of the resulting strains is shown in Figure 7B. After 48 h of fermentation, the results were shown in Figure 7D, when the sfGFP gene is integrated into chromosomes D and F in the YlΔA1 strain, the green fluorescence intensity was 4.30 and 6.46 times higher than ATCC201249, respectively. But the improvement was not observed when the sfGFP gene was integrated into the A and E chromosomes and the integration on the B and C chromosomes failed. Then, we selected two sets of promoters with or without intron sequences (Figure 7C) including TEFin and TEF (Tai and Stephanopoulos, 2013), FBAin and FBA (Hong et al., 2012) to express sfGFP gene in strains ATCC 201249 and YlΔA1. The fluorescence intensity of the resulting strains was shown in Figure 7E, the expression of sfGFP gene driven by the promoter with intron sequences in YlΔA1 strain was higher than ATCC 201249, but the improvement was not observed when the promoters without intron sequences. This demonstrated that disruption of the A1 gene prefers to enhance gene expression by promoters with intron sequences. And in the above results the Hp4d promoter also contained intron sequence which further proved this conclusion (Figure 1B). Overall, disruption of the A1 gene increased protein expression, which boosted cell viability and may be responsible for the improved cellular synthesis. However, the enhancement of exogenous protein expression from A1 disruption was applicable to D and F chromosomes and promoters with introns.
Figure 7
Discussion
The regulatory mechanisms of the metabolic networks are extremely complex (Nielsen and Keasling, 2016; Chen et al., 2018; Zhu et al., 2020). Therefore, a substantial set of beneficial gene targets across the genome targets that can facilitate the performance of the strains, in addition to synthetic pathways of a desired product. Multiplex automated genome engineering (MAGE; Wang et al., 2009), trackable multiplex recombineering (TRMR; Warner et al., 2010) and other similar approaches can generate genome-scale mutagenesis for the identification of new targets. However, these techniques have not yet been extensively applied to other non-model organisms. Many non-conventional yeasts and fungi display a high preference for DSBs repair by NHEJ, which is several orders of magnitude than homologous recombination (HR; Wagner and Alper, 2016). And DSBs are widely distributed throughout the genome (Dellino et al., 2019). Therefore, NHEJ-mediated integration confers the randomness and genome-scale unbiased coverage of gene mutations (Liu et al., 2022), which can efficiently generate random mutagenesis across the genome in different microorganisms. NHEJ-mediated integration sometimes can still be produced even when HR is performed to integrate genes in NHEJ-preferred strains, generating unexpected superiors phenotypes with additional insertions (Gao et al., 2014). That requires researchers to value and further investigate incidental phenotypic differences to identify genotypic alterations instead of ignoring them.
In conclusion, we identified a novel target YALI0_A00913g (“A1 gene”) that was disrupted by NHEJ-mediated insertional mutagenesis, through whole-genome sequencing an unexpectedly high fatty alcohol-producing strain obtained incidentally. Transcriptome analysis revealed that disruption of the A1 gene increased protein synthesis, which enhanced cell viability and improved product synthesis. Our research confirmed the importance of potential target excavation in the whole genome and provided a new engineering approach, which is of significant value to further stimulate the biosynthetic potential of for microorganisms.
Funding
This work was supported by the National Key Research and Development Program of China (2021YFC2104400), the National Natural Science Foundation of China (NSFC 21621004 and NSFC 22078240), the Natural Science Foundation of Tianjin City (19JCQNJC09200), and the Young Elite Scientists Sponsorship Program by Tianjin (TJSQNTJ-2018-16).
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Statements
Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found at: https://www.ncbi.nlm.nih.gov/bioproject/, PRJNA818142; https://www.ncbi.nlm.nih.gov/geo/, GSE199895.
Author contributions
YC and HS supervised the study. ML and JZ conceived and designed the experiments. ML and QB performed the experiments. ML, JZ, LF, and QB analyzed the data. YC and ML revised the manuscript. All authors contributed to the article and approved the submitted version.
Acknowledgments
We are grateful for the strain resources provided by Ying-Jin Yuan (Tianjin University, China).
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2022.898884/full#supplementary-material
References
1
AjjawiI.VerrutoJ.AquiM.SoriagaL. B.CoppersmithJ.KwokK.et al. (2017). Lipid production in Nannochloropsis gaditana is doubled by decreasing expression of a single transcriptional regulator. Nat. Biotechnol.35, 647–652. doi: 10.1038/nbt.3865
2
AndersS.HuberW. (2010). Differential expression analysis for sequence count data. Genome Biol.11:R106. doi: 10.1186/gb-2010-11-10-r106
3
BaiQ.ChengS.ZhangJ.LiM.CaoY.YuanY. (2021). Establishment of genomic library technology mediated by non-homologous end joining mechanism in Yarrowia lipolytica. Sci. China Life Sci.64, 2114–2128. doi: 10.1007/s11427-020-1885-x
4
BeucherA.BirrauxJ.TchouandongL.BartonO.ShibataA.ConradS.et al. (2009). ATM and Artemis promote homologous recombination of radiation-induced DNA double-strand breaks in G2. EMBO J.28, 3413–3427. doi: 10.1038/emboj.2009.276
5
BlazeckJ.HillA.LiuL.KnightR.MillerJ.PanA.et al. (2014). Harnessing Yarrowia lipolytica lipogenesis to create a platform for lipid and biofuel production. Nat. Commun.5:3131. doi: 10.1038/ncomms4131
6
BouwmanB. A. M.CrosettoN. (2018). Endogenous DNA double-Strand breaks during DNA transactions: emerging insights and methods for genome-wide profiling. Genes9:632. doi: 10.3390/genes9120632
7
BurmaS.ChenB. P.ChenD. J. (2006). Role of non-homologous end joining (NHEJ) in maintaining genomic integrity. DNA Repair (Amst)5, 1042–1048. doi: 10.1016/j.dnarep.2006.05.026
8
CavarelliJ.DelagoutteB.ErianiG.GangloffJ.MorasD. (1998). L-arginine recognition by yeast arginyl-tRNA synthetase. EMBO J.17, 5438–5448. doi: 10.1093/emboj/17.18.5438
9
ChamiehH.BallutL.BonneauF.Le HirH. (2008). NMD factors UPF2 and UPF3 bridge UPF1 to the exon junction complex and stimulate its RNA helicase activity. Nat. Struct. Mol. Biol.15, 85–93. doi: 10.1038/nsmb1330
10
ChenX.GaoC.GuoL.HuG.LuoQ.LiuJ.et al. (2018). DCEO biotechnology: tools to design, construct, evaluate, and optimize the metabolic pathway for biosynthesis of chemicals. Chem. Rev.118, 4–72. doi: 10.1021/acs.chemrev.6b00804
11
ChooA. Y.KimS. G.Vander HeidenM. G.MahoneyS. J.VuH.YoonS. O.et al. (2010). Glucose addiction of TSC null cells is caused by failed mTORC1-dependent balancing of metabolic demand with supply. Mol. Cell38, 487–499. doi: 10.1016/j.molcel.2010.05.007
12
CiechanoverA.OrianA.SchwartzA. L. (2000). Ubiquitin-mediated proteolysis: biological regulation via destruction. BioEssays22, 442–451. doi: 10.1002/(SICI)1521-1878(200005)22:5<442::AID-BIES6>3.0.CO;2-Q
13
CouxO.TanakaK.GoldbergA. L. (1996). Structure and functions of the 20S and 26S proteasomes. Annu. Rev. Biochem.65, 801–847. doi: 10.1146/annurev.bi.65.070196.004101
14
CuiZ.JiangX.ZhengH.QiQ.HouJ. (2019). Homology-independent genome integration enables rapid library construction for enzyme expression and pathway optimization in Yarrowia lipolytica. Biotechnol. Bioeng.116, 354–363. doi: 10.1002/bit.26863
15
CusackS. (1997). Aminoacyl-tRNA synthetases. Curr. Opin. Struct. Biol.7, 881–889. doi: 10.1016/S0959-440X(97)80161-3
16
D’AlessioC.CarameloJ. J.ParodiA. J. (2005). Absence of nucleoside diphosphatase activities in the yeast secretory pathway does not abolish nucleotide sugar-dependent protein glycosylation. J. Biol. Chem.280, 40417–40427. doi: 10.1074/jbc.M503149200
17
DelagoutteB.MorasD.CavarelliJ. (2000). tRNA aminoacylation by arginyl-tRNA synthetase: induced conformations during substrates binding. EMBO J.19, 5599–5610. doi: 10.1093/emboj/19.21.5599
18
DellinoG. I.PalluzziF.ChiarielloA. M.PiccioniR.BiancoS.FuriaL.et al. (2019). Release of paused RNA polymerase II at specific loci favors DNA double-strand-break formation and promotes cancer translocations. Nat. Genet.51, 1011–1023. doi: 10.1038/s41588-019-0421-z
19
DudekJ.BenedixJ.CappelS.GreinerM.JalalC.MullerL.et al. (2009). Functions and pathologies of BiP and its interaction partners. Cell. Mol. Life Sci.66, 1556–1569. doi: 10.1007/s00018-009-8745-y
20
FalconeD.HendersonM. P.NieuwlandH.CoughlanC. M.BrodskyJ. L.AndrewsD. W. (2011). Stability and function of the Sec61 translocation complex depends on the Sss1p tail-anchor sequence. Biochem. J.436, 291–303. doi: 10.1042/BJ20101865
21
FangL.FanJ.LuoS.ChenY.WangC.CaoY.et al. (2021). Genome-scale target identification in Escherichia coli for high-titer production of free fatty acids. Nat. Commun.12:4976. doi: 10.1038/s41467-021-25243-w
22
FilletS.AdrioJ. L. (2016). Microbial production of fatty alcohols. World J. Microbiol. Biotechnol.32:152. doi: 10.1007/s11274-016-2099-z
23
Fromont-RacineM.SengerB.SaveanuC.FasioloF. (2003). Ribosome assembly in eukaryotes. Gene313, 17–42. doi: 10.1016/S0378-1119(03)00629-2
24
GaoS.HanL.ZhuL.GeM.YangS.JiangY.et al. (2014). One-step integration of multiple genes into the oleaginous yeast Yarrowia lipolytica. Biotechnol. Lett.36, 2523–2528. doi: 10.1007/s10529-014-1634-y
25
GaoS.TongY.WenZ.ZhuL.GeM.ChenD.et al. (2016). Multiplex gene editing of the Yarrowia lipolytica genome using the CRISPR-Cas9 system. J. Ind. Microbiol. Biotechnol.43, 1085–1093. doi: 10.1007/s10295-016-1789-8
26
GhazaeiC. (2017). Role and mechanism of the Hsp70 molecular chaperone machines in bacterial pathogens. J. Med. Microbiol.66, 259–265. doi: 10.1099/jmm.0.000429
27
GingrasA. C.RaughtB.SonenbergN. (1999). eIF4 initiation factors: effectors of mRNA recruitment to ribosomes and regulators of translation. Annu. Rev. Biochem.68, 913–963. doi: 10.1146/annurev.biochem.68.1.913
28
HerrgardM. J.LeeB. S.PortnoyV.PalssonB. O. (2006). Integrated analysis of regulatory and metabolic networks reveals novel regulatory mechanisms in Saccharomyces cerevisiae. Genome Res.16, 627–635. doi: 10.1101/gr.4083206
29
HoA. K.ShenT. X.RyanK. J.KiselevaE.LevyM. A.AllenT. D.et al. (2000). Assembly and preferential localization of Nup116p on the cytoplasmic face of the nuclear pore complex by interaction with Nup82p. Mol. Cell. Biol.20, 5736–5748. doi: 10.1128/MCB.20.15.5736-5748.2000
30
HolkenbrinkC.DamM. I.KildegaardK. R.BederJ.DahlinJ.Domenech BeldaD.et al. (2018). EasyCloneYALI: CRISPR/Cas9-based synthetic toolbox for engineering of the yeast Yarrowia lipolytica. Biotechnol. J.13:e1700543. doi: 10.1002/biot.201700543
31
HongS. P.SeipJ.Walters-PollakD.RupertR.JacksonR.XueZ.et al. (2012). Engineering Yarrowia lipolytica to express secretory invertase with strong FBA1IN promoter. Yeast29, 59–72. doi: 10.1002/yea.1917
32
HuD.WangZ.HeM.MaY. (2021). Functional gene identification and corresponding tolerant mechanism of high furfural-tolerant Zymomonas mobilis strain F211. Front. Microbiol.12:736583. doi: 10.3389/fmicb.2021.736583
33
HuntC.MorimotoR. I. (1985). Conserved features of eukaryotic hsp70 genes revealed by comparison with the nucleotide sequence of human hsp70. Proc. Natl. Acad. Sci. U. S. A.82, 6455–6459. doi: 10.1073/pnas.82.19.6455
34
IkedaF.YoshidaK.TokiT.UechiT.IshidaS.NakajimaY.et al. (2017). Exome sequencing identified RPS15A as a novel causative gene for diamond-Blackfan anemia. Haematologica102, e93–e96. doi: 10.3324/haematol.2016.153932
35
JiangX. R.WangH.ShenR.ChenG. Q. (2015). Engineering the bacterial shapes for enhanced inclusion bodies accumulation. Metab. Eng.29, 227–237. doi: 10.1016/j.ymben.2015.03.017
36
JinC. C.ZhangJ. L.SongH.CaoY. X. (2019). Boosting the biosynthesis of betulinic acid and related triterpenoids in Yarrowia lipolytica via multimodular metabolic engineering. Microb. Cell Factories18:77. doi: 10.1186/s12934-019-1127-8
37
JohnsonE. S.MaP. C.OtaI. M.VarshavskyA. (1995). A proteolytic pathway that recognizes ubiquitin as a degradation signal. J. Biol. Chem.270, 17442–17456. doi: 10.1074/jbc.270.29.17442
38
KretzschmarA.OttoC.HolzM.WernerS.HubnerL.BarthG. (2013). Increased homologous integration frequency in Yarrowia lipolytica strains defective in non-homologous end-joining. Curr. Genet.59, 63–72. doi: 10.1007/s00294-013-0389-7
39
KrullS.ThybergJ.BjorkrothB.RackwitzH. R.CordesV. C. (2004). Nucleoporins as components of the nuclear pore complex core structure and Tpr as the architectural element of the nuclear basket. Mol. Biol. Cell15, 4261–4277. doi: 10.1091/mbc.e04-03-0165
40
KumarA. (2016). Multipurpose Transposon-Insertion Libraries in Yeast. Cold Spring Harb. Protoc.6, 499–502. doi: 10.1101/pdb.top080259
41
LangmeadB.SalzbergS. L. (2012). Fast gapped-read alignment with bowtie 2. Nat. Methods9, 357–359. doi: 10.1038/nmeth.1923
42
Levy-StrumpfN.DeissL. P.BerissiH.KimchiA. (1997). DAP-5, a novel homolog of eukaryotic translation initiation factor 4G isolated as a putative modulator of gamma interferon-induced programmed cell death. Mol. Cell. Biol.17, 1615–1625. doi: 10.1128/MCB.17.3.1615
43
LiH.AlperH. S. (2016). Enabling xylose utilization in Yarrowia lipolytica for lipid production. Biotechnol. J.11, 1230–1240. doi: 10.1002/biot.201600210
44
LiB.DeweyC. N. (2011). RSEM: accurate transcript quantification from RNA-Seq data with or without a reference genome. BMC Bioinform.12:323. doi: 10.1186/1471-2105-12-323
45
LiH.TangX.YangX.ZhangH. (2021). Comprehensive transcriptome and metabolome profiling reveal metabolic mechanisms of Nitraria sibirica pall. To salt stress. Sci. Rep.11:12878. doi: 10.1038/s41598-021-92317-6
46
LieberM. R. (2010). The mechanism of double-strand DNA break repair by the nonhomologous DNA end-joining pathway. Annu. Rev. Biochem.79, 181–211. doi: 10.1146/annurev.biochem.052308.093131
47
LincecumT. L.TukaloM.YaremchukA.MursinnaR. S.WilliamsA. M.SproatB. S.et al. (2003). Structural and mechanistic basis of pre- and Posttransfer editing by Leucyl-tRNA Synthetase. Mol. Cell11, 951–963. doi: 10.1016/S1097-2765(03)00098-4
48
LingJ.ReynoldsN.IbbaM. (2009). Aminoacyl-tRNA synthesis and translational quality control. Annu. Rev. Microbiol.63, 61–78. doi: 10.1146/annurev.micro.091208.073210
49
LiuI. C.ChiuS. W.LeeH. Y.LeuJ. Y. (2012). The histone deacetylase Hos2 forms an Hsp42-dependent cytoplasmic granule in quiescent yeast cells. Mol. Biol. Cell23, 1231–1242. doi: 10.1091/mbc.e11-09-0752
50
LiuY.JiangX.CuiZ.WangZ.QiQ.HouJ. (2019). Engineering the oleaginous yeast Yarrowia lipolytica for production of alpha-farnesene. Biotechnol. Biofuels12:296. doi: 10.1186/s13068-019-1636-z
51
LiuX.LiuM.ZhangJ.ChangY.CuiZ.JiB.et al. (2022). Mapping of nonhomologous end joining-mediated integration facilitates genome-scale Trackable mutagenesis in Yarrowia lipolytica. ACS Synth. Biol.11, 216–227. doi: 10.1021/acssynbio.1c00390
52
LiuL.MarkhamK.BlazeckJ.ZhouN.LeonD.OtoupalP.et al. (2015a). Surveying the lipogenesis landscape in Yarrowia lipolytica through understanding the function of a Mga2p regulatory protein mutant. Metab. Eng.31, 102–111. doi: 10.1016/j.ymben.2015.07.004
53
LiuL.PanA.SpoffordC.ZhouN.AlperH. S. (2015b). An evolutionary metabolic engineering approach for enhancing lipogenesis in Yarrowia lipolytica. Metab. Eng.29, 36–45. doi: 10.1016/j.ymben.2015.02.003
54
LongmuirS.AkhtarN.MacNeillS. A. (2019). Unexpected insertion of carrier DNA sequences into the fission yeast genome during CRISPR-Cas9 mediated gene deletion. BMC. Res. Notes12:191. doi: 10.1186/s13104-019-4228-x
55
MadzakC.TretonB.Blanchin-RolandS. (2000). Strong hybrid promoters and integrative expression/secretion vectors for quasi-constitutive expression of heterologous proteins in the yeast Yarrowia lipolytica. J. Mol. Microbiol. Biotechnol.2, 207–216.
56
MaozT.KorenR.Ben-AriI.KleinbergerT. (2005). YND1 interacts with CDC55 and is a novel mediator of E4orf4-induced toxicity. J. Biol. Chem.280, 41270–41277. doi: 10.1074/jbc.M507281200
57
MatthausF.KetelhotM.GatterM.BarthG. (2014). Production of lycopene in the non-carotenoid-producing yeast Yarrowia lipolytica. Appl. Environ. Microbiol.80, 1660–1669. doi: 10.1128/AEM.03167-13
58
MayerM. P.BukauB. (2005). Hsp70 chaperones: cellular functions and molecular mechanism. Cell. Mol. Life Sci.62, 670–684. doi: 10.1007/s00018-004-4464-6
59
McClellanA. J.EndresJ. B.VogelJ. P.PalazziD.RoseM. D.BrodskyJ. L. (1998). Specific molecular chaperone interactions and an ATP-dependent conformational change are required during posttranslational protein translocation into the yeast ER. Mol. Biol. Cell9, 3533–3545. doi: 10.1091/mbc.9.12.3533
60
MehlerA. H.MitraS. K. (1967). The activation of Arginyl transfer ribonucleic acid Synthetase by transfer ribonucleic acid. J. Biol. Chem.242, 5495–5499. doi: 10.1016/S0021-9258(18)99386-5
61
MittelmanK.ZivK.MaozT.KleinbergerT. (2010). The cytosolic tail of the Golgi apyrase Ynd1 mediates E4orf4-induced toxicity in Saccharomyces cerevisiae. PLoS One5:e15539. doi: 10.1371/journal.pone.0015539
62
MurrayM. G.ThompsonW. F. (1980). Rapid isolation of high molecular weight plant DNA. Nucleic Acids Res.8, 4321–4326. doi: 10.1093/nar/8.19.4321
63
NegrutskiiB. S.DeutscherM. P. (1991). Channeling of aminoacyl-tRNA for protein synthesis in vivo. Proc. Natl. Acad. Sci. U. S. A.88, 4991–4995. doi: 10.1073/pnas.88.11.4991
64
NicaudJ. M.MadzakC.van den BroekP.GyslerC.DubocP.NiederbergerP.et al. (2002). Protein expression and secretion in the yeast Yarrowia lipolytica. FEMS Yeast Res.2, 371–379. doi: 10.1111/j.1567-1364.2002.tb00106.x
65
NielsenJ.KeaslingJ. D. (2016). Engineering cellular metabolism. Cell164, 1185–1197. doi: 10.1016/j.cell.2016.02.004
66
OnumaK.WatanabeA.KanzakiN.KubotaT. (2006). Association kinetics of wild- and mutant-type Ynd1p in relation to quality of grown crystals. J. Phys. Chem. B110, 24876–24883. doi: 10.1021/jp0643146
67
PainV. M.StandartN. (2001). Translational control of gene expression, edited by Nahum Sonenberg, John W.B. Hershey, and Michael B. Mathews. 2000. Cold Spring Harbor, New York: cold Spring Harbor laboratory press. Cloth, 1,020 pp. $113. RNA7, 331–333. doi: 10.1017/S1355838201002540
68
PannunzioN. R.WatanabeG.LieberM. R. (2018). Nonhomologous DNA end-joining for repair of DNA double-strand breaks. J. Biol. Chem.293, 10512–10523. doi: 10.1074/jbc.TM117.000374
69
ParkS. G.SchimmelP.KimS. (2008). Aminoacyl tRNA synthetases and their connections to disease. Proc. Natl. Acad. Sci. U. S. A.105, 11043–11049. doi: 10.1073/pnas.0802862105
70
ParkY.YoonS. K.YoonJ. B. (2009). The HECT domain of TRIP12 ubiquitinates substrates of the ubiquitin fusion degradation pathway. J. Biol. Chem.284, 1540–1549. doi: 10.1074/jbc.M807554200
71
PedelacqJ. D.CabantousS.TranT.TerwilligerT. C.WaldoG. S. (2006). Engineering and characterization of a superfolder green fluorescent protein. Nat. Biotechnol.24, 79–88. doi: 10.1038/nbt1172
72
Perez-VargasJ.RomeroP.LopezS.AriasC. F. (2006). The peptide-binding and ATPase domains of recombinant hsc70 are required to interact with rotavirus and reduce its infectivity. J. Virol.80, 3322–3331. doi: 10.1128/JVI.80.7.3322-3331.2006
73
PflegerB. F.GossingM.NielsenJ. (2015). Metabolic engineering strategies for microbial synthesis of oleochemicals. Metab. Eng.29, 1–11. doi: 10.1016/j.ymben.2015.01.009
74
PhanL.SchoenfeldL. W.ValasekL.NielsenK. H.HinnebuschA. G. (2001). A subcomplex of three eIF3 subunits binds eIF1 and eIF5 and stimulates ribosome binding of mRNA and tRNA(i)met. EMBO J.20, 2954–2965. doi: 10.1093/emboj/20.11.2954
75
RanjhaL.HowardS. M.CejkaP. (2018). Main steps in DNA double-strand break repair: an introduction to homologous recombination and related processes. Chromosoma127, 187–214. doi: 10.1007/s00412-017-0658-1
76
RapoportT. A. (2007). Protein translocation across the eukaryotic endoplasmic reticulum and bacterial plasma membranes. Nature450, 663–669. doi: 10.1038/nature06384
77
RichardG. F.KerrestA.LafontaineI.DujonB. (2005). Comparative genomics of hemiascomycete yeasts: genes involved in DNA replication, repair, and recombination. Mol. Biol. Evol.22, 1011–1023. doi: 10.1093/molbev/msi083
78
RobsonA.CollinsonI. (2006). The structure of the sec complex and the problem of protein translocation. EMBO Rep.7, 1099–1103. doi: 10.1038/sj.embor.7400832
79
RodriguezM. S.DargemontC.StutzF. (2004). Nuclear export of RNA. Biol. Cell.96, 639–655. doi: 10.1016/j.biolcel.2004.04.014
80
RoutM. P.AitchisonJ. D. (2001). The nuclear pore complex as a transport machine. J. Biol. Chem.276, 16593–16596. doi: 10.1074/jbc.R100015200
81
RoutM. P.AitchisonJ. D.SupraptoA.HjertaasK.ZhaoY.ChaitB. T. (2000). The yeast nuclear pore complex: composition, architecture, and transport mechanism. J. Cell Biol.148, 635–652. doi: 10.1083/jcb.148.4.635
82
SchwartzC.ChengJ. F.EvansR.SchwartzC. A.WagnerJ. M.AnglinS.et al. (2019). Validating genome-wide CRISPR-Cas9 function improves screening in the oleaginous yeast Yarrowia lipolytica. Metab. Eng.55, 102–110. doi: 10.1016/j.ymben.2019.06.007
83
SchwartzC. M.HussainM. S.BlennerM.WheeldonI. (2016). Synthetic RNA polymerase III promoters facilitate high-efficiency CRISPR-Cas9-mediated genome editing in Yarrowia lipolytica. ACS Synth. Biol.5, 356–359. doi: 10.1021/acssynbio.5b00162
84
ShiomiD.ToyodaA.AizuT.EjimaF.FujiyamaA.ShiniT.et al. (2013). Mutations in cell elongation genes mreB, mrdA and mrdB suppress the shape defect of RodZ-deficient cells. Mol. Microbiol.87, 1029–1044. doi: 10.1111/mmi.12148
85
ShrivastavM.De HaroL. P.NickoloffJ. A. (2008). Regulation of DNA double-strand break repair pathway choice. Cell Res.18, 134–147. doi: 10.1038/cr.2007.111
86
SmithD. M.BenaroudjN.GoldbergA. (2006). Proteasomes and their associated ATPases: a destructive combination. J. Struct. Biol.156, 72–83. doi: 10.1016/j.jsb.2006.04.012
87
StephensonK. (2005). Sec-dependent protein translocation across biological membranes: evolutionary conservation of an essential protein transport pathway (review). Mol. Membr. Biol.22, 17–28. doi: 10.1080/09687860500063308
88
SudfeldC.HubacekM.FigueiredoD.NaduthodiM. I. S.van der OostJ.WijffelsR. H.et al. (2021). High-throughput insertional mutagenesis reveals novel targets for enhancing lipid accumulation in Nannochloropsis oceanica. Metab. Eng.66, 239–258. doi: 10.1016/j.ymben.2021.04.012
89
SunW.BernardC.van de CotteB.Van MontaguM.VerbruggenN. (2001). At-HSP17.6A, encoding a small heat-shock protein in Arabidopsis, can enhance osmotolerance upon overexpression. Plant J.27, 407–415. doi: 10.1046/j.1365-313X.2001.01107.x
90
TaiM.StephanopoulosG. (2013). Engineering the push and pull of lipid biosynthesis in oleaginous yeast Yarrowia lipolytica for biofuel production. Metab. Eng.15, 1–9. doi: 10.1016/j.ymben.2012.08.007
91
TarazonaS.Garcia-AlcaldeF.DopazoJ.FerrerA.ConesaA. (2011). Differential expression in RNA-seq: a matter of depth. Genome Res.21, 2213–2223. doi: 10.1101/gr.124321.111
92
ValasekL. S. (2012). ‘Ribozoomin’ – translation initiation from the perspective of the ribosome-bound eukaryotic initiation factors (eIFs). Curr. Protein Pept. Sci.13, 305–330. doi: 10.2174/138920312801619385
93
van OpijnenT.BodiK. L.CamilliA. (2009). Tn-seq: high-throughput parallel sequencing for fitness and genetic interaction studies in microorganisms. Nat. Methods6, 767–772. doi: 10.1038/nmeth.1377
94
VerbekeJ.BeopoulosA.NicaudJ. M. (2013). Efficient homologous recombination with short length flanking fragments in Ku70 deficient Yarrowia lipolytica strains. Biotechnol. Lett.35, 571–576. doi: 10.1007/s10529-012-1107-0
95
VladimirovS. N.IvanovA. V.KarpovaG. G.MusolyamovA. K.EgorovT. A.ThiedeB.et al. (1996). Characterization of the human small-ribosomal-subunit proteins by N-terminal and internal sequencing, and mass spectrometry. Eur. J. Biochem.239, 144–149. doi: 10.1111/j.1432-1033.1996.0144u.x
96
VogesD.ZwicklP.BaumeisterW. (1999). The 26S proteasome: a molecular machine designed for controlled proteolysis. Annu. Rev. Biochem.68, 1015–1068. doi: 10.1146/annurev.biochem.68.1.1015
97
WagnerJ. M.AlperH. S. (2016). Synthetic biology and molecular genetics in non-conventional yeasts: current tools and future advances. Fungal Genet. Biol.89, 126–136. doi: 10.1016/j.fgb.2015.12.001
98
WahlenB. D.OswaldW. S.SeefeldtL. C.BarneyB. M. (2009). Purification, characterization, and potential bacterial wax production role of an NADPH-dependent fatty aldehyde reductase from Marinobacter aquaeolei VT8. Appl. Environ. Microbiol.75, 2758–2764. doi: 10.1128/AEM.02578-08
99
WangH. H.IsaacsF. J.CarrP. A.SunZ. Z.XuG.ForestC. R.et al. (2009). Programming cells by multiplex genome engineering and accelerated evolution. Nature460, 894–898. doi: 10.1038/nature08187
100
WangG.ShiT.ChenT.WangX.WangY.LiuD.et al. (2018). Integrated whole-genome and transcriptome sequence analysis reveals the genetic characteristics of a riboflavin-overproducing Bacillus subtilis. Metab. Eng.48, 138–149. doi: 10.1016/j.ymben.2018.05.022
101
WangW.WeiH.KnoshaugE.Van WychenS.XuQ.HimmelM. E.et al. (2016b). Fatty alcohol production in Lipomyces starkeyi and Yarrowia lipolytica. Biotechnol. Biofuels9:227. doi: 10.1186/s13068-016-0647-2
102
WangG.XiongX.GhogareR.WangP.MengY.ChenS. (2016a). Exploring fatty alcohol-producing capability of Yarrowia lipolytica. Biotechnol. Biofuels9:107. doi: 10.1186/s13068-016-0512-3
103
WarnerJ. R.ReederP. J.Karimpour-FardA.WoodruffL. B.GillR. T. (2010). Rapid profiling of a microbial genome using mixtures of barcoded oligonucleotides. Nat. Biotechnol.28, 856–862. doi: 10.1038/nbt.1653
104
WeiL. J.CaoX.LiuJ. J.KwakS.JinY. S.WangW.et al. (2021). Increased accumulation of Squalene in engineered Yarrowia lipolytica through deletion of PEX10 and URE2. Appl. Environ. Microbiol.87:e0048121. doi: 10.1128/AEM.00481-21
105
WillisR. M.WahlenB. D.SeefeldtL. C.BarneyB. M. (2011). Characterization of a fatty acyl-CoA reductase from Marinobacter aquaeolei VT8: a bacterial enzyme catalyzing the reduction of fatty acyl-CoA to fatty alcohol. Biochemistry50, 10550–10558. doi: 10.1021/bi2008646
106
XuP.QiaoK.AhnW. S.StephanopoulosG. (2016). Engineering Yarrowia lipolytica as a platform for synthesis of drop-in transportation fuels and oleochemicals. Proc. Natl. Acad. Sci. U. S. A.113, 10848–10853. doi: 10.1073/pnas.1607295113
107
XuanJ. W.FournierP.DeclerckN.ChaslesM.GaillardinC. (1990). Overlapping reading frames at the LYS5 locus in the yeast Yarrowia lipolytica. Mol. Cell. Biol.10, 4795–4806. doi: 10.1128/MCB.10.9.4795
108
XueZ.SharpeP. L.HongS. P.YadavN. S.XieD.ShortD. R.et al. (2013). Production of omega-3 eicosapentaenoic acid by metabolic engineering of Yarrowia lipolytica. Nat. Biotechnol.31, 734–740. doi: 10.1038/nbt.2622
109
YeR. W.SharpeP. L.ZhuQ. (2012). Bioengineering of oleaginous yeast Yarrowia lipolytica for lycopene production. Methods Mol. Biol.898, 153–159. doi: 10.1007/978-1-61779-918-1_9
110
YoungJ. C. (2010). Mechanisms of the Hsp70 chaperone system. Biochem. Cell Biol.88, 291–300. doi: 10.1139/O09-175
111
ZhangJ.-L.CaoY.-X.PengY.-Z.JinC.-C.BaiQ.-Y.ZhangR.-S.et al. (2019). High production of fatty alcohols in Yarrowia lipolytica by coordination with glycolysis. SCIENCE CHINA Chem.62, 1007–1016. doi: 10.1007/s11426-019-9456-y
112
ZhangJ. L.PengY. Z.LiuD.LiuH.CaoY. X.LiB. Z.et al. (2018). Gene repression via multiplex gRNA strategy in Y. lipolytica. Microb. Cell Factories17:62. doi: 10.1186/s12934-018-0909-8
113
ZhouY. J.BuijsN. A.ZhuZ.QinJ.SiewersV.NielsenJ. (2016). Production of fatty acid-derived oleochemicals and biofuels by synthetic yeast cell factories. Nat. Commun.7:11709. doi: 10.1038/ncomms11709
114
ZhuM.LiM.ZhouW.GeG.ZhangL.JiG. (2020). Metabolomic analysis identifies Glycometabolism pathways as potential targets of Qianggan extract in hyperglycemia rats. Front. Pharmacol.11:671. doi: 10.3389/fphar.2020.00671
Summary
Keywords
new target identification, non-homologous end joining-mediated integration, Yarrowia lipolytica, fatty alcohols, RNA-Seq
Citation
Li M, Zhang J, Bai Q, Fang L, Song H and Cao Y (2022) Non-homologous End Joining-Mediated Insertional Mutagenesis Reveals a Novel Target for Enhancing Fatty Alcohols Production in Yarrowia lipolytica. Front. Microbiol. 13:898884. doi: 10.3389/fmicb.2022.898884
Received
18 March 2022
Accepted
06 April 2022
Published
25 April 2022
Volume
13 - 2022
Edited by
Haichun Gao, Zhejiang University, China
Reviewed by
Huihui Fu, Institute of Oceanology (CAS), China; Jingwen Zhou, Jiangnan University, China
Updates
Copyright
© 2022 Li, Zhang, Bai, Fang, Song and Cao.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Yingxiu Cao, caoyingxiu@tju.edu.cn
This article was submitted to Microbial Physiology and Metabolism, a section of the journal Frontiers in Microbiology
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.