Mengxu Li
Jinlai Zhang1,2
Lixia Fang
Hao Song
Yingxiu Cao- 1Frontier Science Center for Synthetic Biology and Key Laboratory of Systems Bioengineering (Ministry of Education), School of Chemical Engineering and Technology, Tianjin University, Tianjin, China
- 2Key Laboratory of Systems Bioengineering (Ministry of Education), Tianjin University, Tianjin, China
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. Strains with different fatty alcohols yields were found in the same batch of transformed of FAR gene. (A) Schematic diagram of the fatty alcohol biosynthesis pathway engineering strategies in Yarrowia lipolytica. FAR, fatty acyl-CoA reductase gene from Marinobacter aquaeolei VT8; PEX10, the peroxisome synthesis gene. Fatty acyl-CoA reductase converts fatty acyl-CoA to fatty alcohols. The gene PEX10 was deleted to prevent the degradation of fatty alcohols. (B) Schematic diagram of the plasmid pINA1269-FAR used to transform Y. lipolytica. The FAR gene was driven by the Hp4d promoter and the plasmid was linearized by NotI and transformed into different basic strains and two transformants were randomly selected from each plate for fermentation (C). (C) Schematic diagram of establishing exogenous synthesis pathway of fatty alcohols. (D) Fatty alcohols production of strains from the same batch of transformants. All strains were cultured in rich YPD medium. Samples were taken at 72 h. Data are presented as mean ± SD of three biological replicates.
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. Whole-genome sequencing determined gene disruption in fatty alcohol high-yielding strains. YlA and YlP strains were applied to the whole-genome sequencing. The FAR gene marked in orange was integrated into the chromosome F as expected. The gene disruption caused by the integration of broken vector fragments into chromosomes A and E was marked with red vertical line, named A1 and E1 sites, respectively. The A1 site was located in the 118,335–118,543 region of chromosomes A, and the E1 site was in the 4,134,688–4,134,892 region of chromosomes E.
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. Reverse engineering of the candidate targets derived from WGS to identify beneficial gene for enhanced fatty alcohol production. (A) CRISPR/Cas-mediated knockout of candidate gene targets in wild-type Y. lipolytica strain ATCC 201249. The deletion of bases at A1 and E1 sites in ATCC 201249 genome were marked with red and strikethrough. The knockout location is contained within the location of the broken vector fragment insertion. (B) Schematic diagram of the plasmid pIntF-FAR used to transform Y. lipolytica. The FAR gene is driven by the TEFin promoter and flanked by 2 homology arms of IntF location. The plasmid was linearized by NotI and transformed into different basic strains and two transformants were randomly selected from each plate for fermentation. (C) The effect of A1 and E1 sites disruption on fatty alcohols production. All strains were cultured in rich YPD medium. Values are the mean of three biological replicates ± standard deviation (n = 3) after 72 h.
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. Transcriptome analysis of YlΔA1-FAR and WT-FAR. (A) Growth rate of YlΔA1-FAR strain with high fatty alcohols production compared with WT-FAR under rich medium conditions. All cultivations were performed in rich YPD medium for 108 h. Transcriptome analysis (RNA-Seq) was performed on the samples at 24 h (pointed by arrow). Results are presented as mean ± s. d. of three biological replicates. (B) Fold changes (log2) of genes differentially expressed (adjusted p < 0.05, |Fold change| > 1) in YlΔA1-FAR cells relative to WT-FAR cells. Genes at the positive side of the x axis are upregulated (highlighted in red), and genes at the negative side are downregulated (highlighted in blue). (C) KEGG pathways enrichment analysis. Gene sets were defined by KEGG pathways, which showed the percentage of genes that were either up- (red) or down-regulated (blue).
Table 1. 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. Changes in transcript levels of the protein synthesis process. The (A) RNA transport, (B) Ribosome and (C) Aminoacyl-tRNA synthesis pathway belong to the translation gene set, the (D) protein translocation, and (E) protein processing and (F) protein degradation pathway belong to the folding, sorting and degradation gene set. Among the six pathways, the up- and down-regulated pathways are marked with red and blue arrows, respectively. The boxes indicate fold-change in gene expression of YlΔA1-FAR compared to WT-FAR control. Red boxes indicate upregulated gene expression, and blue boxes indicate downregulation gene expression (p < 0.05).
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. Verification of improved protein expression by sfGFP gene expression. (A) Schematic diagram of the sfGFP insertion cassette used to transform Y. lipolytica. The sfGFP gene was driven by the TEFin promoter and stopped by the LIP1 terminator, flanking by two 500 bp homologous arms to IntF site. (B) The A1 gene disruption increased protein expression. Characterization of protein expression by green fluorescence intensity. All strains were cultured in SC-Leu medium. Values are the mean of three biological replicates ± standard deviation (n = 3) after 48 h. (C) Fluorescence images of the strains YGF01 and YGF02. The images of sfGFP fluorescence were observed by a fluorescence microscope (Olympus CX41, Japan).
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. Additional regularities for enhanced expression of exogenous protein. (A) Schematic diagram of transformed fragments with different promoters and homologous arms. (B) Schematic diagram of the distribution of IntX integration sites on the chromosome in Y. lipolytica. (C) Schematic diagram of the promoters with or without intron sequences. (D) The impact of different integration sites on protein expression from A1 disruption. Characterization of protein expression by green fluorescence intensity. The octothorpe (#) indicated the integration failure of sfGFP gene on chromosome B and C. (E) The impact of different promoters on protein expression from A1 disruption. Characterization of protein expression by green fluorescence intensity. All strains were cultured in SC-Leu medium. Values are the mean of three biological replicates standard deviation (n = 3) after 48 h.
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.
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.
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).
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.
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.
Acknowledgments
We are grateful for the strain resources provided by Ying-Jin Yuan (Tianjin University, China).
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
Footnotes
References
Ajjawi, I., Verruto, J., Aqui, M., Soriaga, L. B., Coppersmith, J., Kwok, K., 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
Anders, S., and Huber, W. (2010). Differential expression analysis for sequence count data. Genome Biol. 11:R106. doi: 10.1186/gb-2010-11-10-r106
Bai, Q., Cheng, S., Zhang, J., Li, M., Cao, Y., and Yuan, Y. (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
Beucher, A., Birraux, J., Tchouandong, L., Barton, O., Shibata, A., Conrad, S., 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
Blazeck, J., Hill, A., Liu, L., Knight, R., Miller, J., Pan, A., et al. (2014). Harnessing Yarrowia lipolytica lipogenesis to create a platform for lipid and biofuel production. Nat. Commun. 5:3131. doi: 10.1038/ncomms4131
Bouwman, B. A. M., and Crosetto, N. (2018). Endogenous DNA double-Strand breaks during DNA transactions: emerging insights and methods for genome-wide profiling. Genes 9:632. doi: 10.3390/genes9120632
Burma, S., Chen, B. P., and Chen, D. 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
Cavarelli, J., Delagoutte, B., Eriani, G., Gangloff, J., and Moras, D. (1998). L-arginine recognition by yeast arginyl-tRNA synthetase. EMBO J. 17, 5438–5448. doi: 10.1093/emboj/17.18.5438
Chamieh, H., Ballut, L., Bonneau, F., and Le Hir, H. (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
Chen, X., Gao, C., Guo, L., Hu, G., Luo, Q., Liu, J., 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
Choo, A. Y., Kim, S. G., Vander Heiden, M. G., Mahoney, S. J., Vu, H., Yoon, S. O., et al. (2010). Glucose addiction of TSC null cells is caused by failed mTORC1-dependent balancing of metabolic demand with supply. Mol. Cell 38, 487–499. doi: 10.1016/j.molcel.2010.05.007
Ciechanover, A., Orian, A., and Schwartz, A. L. (2000). Ubiquitin-mediated proteolysis: biological regulation via destruction. BioEssays 22, 442–451. doi: 10.1002/(SICI)1521-1878(200005)22:5<442::AID-BIES6>3.0.CO;2-Q
Coux, O., Tanaka, K., and Goldberg, A. 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
Cui, Z., Jiang, X., Zheng, H., Qi, Q., and Hou, J. (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
Cusack, S. (1997). Aminoacyl-tRNA synthetases. Curr. Opin. Struct. Biol. 7, 881–889. doi: 10.1016/S0959-440X(97)80161-3
D’Alessio, C., Caramelo, J. J., and Parodi, A. 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
Delagoutte, B., Moras, D., and Cavarelli, J. (2000). tRNA aminoacylation by arginyl-tRNA synthetase: induced conformations during substrates binding. EMBO J. 19, 5599–5610. doi: 10.1093/emboj/19.21.5599
Dellino, G. I., Palluzzi, F., Chiariello, A. M., Piccioni, R., Bianco, S., Furia, L., 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
Dudek, J., Benedix, J., Cappel, S., Greiner, M., Jalal, C., Muller, L., 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
Falcone, D., Henderson, M. P., Nieuwland, H., Coughlan, C. M., Brodsky, J. L., and Andrews, D. 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
Fang, L., Fan, J., Luo, S., Chen, Y., Wang, C., Cao, Y., 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
Fillet, S., and Adrio, J. L. (2016). Microbial production of fatty alcohols. World J. Microbiol. Biotechnol. 32:152. doi: 10.1007/s11274-016-2099-z
Fromont-Racine, M., Senger, B., Saveanu, C., and Fasiolo, F. (2003). Ribosome assembly in eukaryotes. Gene 313, 17–42. doi: 10.1016/S0378-1119(03)00629-2
Gao, S., Han, L., Zhu, L., Ge, M., Yang, S., Jiang, Y., 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
Gao, S., Tong, Y., Wen, Z., Zhu, L., Ge, M., Chen, D., 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
Ghazaei, C. (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
Gingras, A. C., Raught, B., and Sonenberg, N. (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
Herrgard, M. J., Lee, B. S., Portnoy, V., and Palsson, B. 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
Ho, A. K., Shen, T. X., Ryan, K. J., Kiseleva, E., Levy, M. A., Allen, T. 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
Holkenbrink, C., Dam, M. I., Kildegaard, K. R., Beder, J., Dahlin, J., Domenech Belda, D., 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
Hong, S. P., Seip, J., Walters-Pollak, D., Rupert, R., Jackson, R., Xue, Z., et al. (2012). Engineering Yarrowia lipolytica to express secretory invertase with strong FBA1IN promoter. Yeast 29, 59–72. doi: 10.1002/yea.1917
Hu, D., Wang, Z., He, M., and Ma, Y. (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
Hunt, C., and Morimoto, R. 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
Ikeda, F., Yoshida, K., Toki, T., Uechi, T., Ishida, S., Nakajima, Y., et al. (2017). Exome sequencing identified RPS15A as a novel causative gene for diamond-Blackfan anemia. Haematologica 102, e93–e96. doi: 10.3324/haematol.2016.153932
Jiang, X. R., Wang, H., Shen, R., and Chen, G. Q. (2015). Engineering the bacterial shapes for enhanced inclusion bodies accumulation. Metab. Eng. 29, 227–237. doi: 10.1016/j.ymben.2015.03.017
Jin, C. C., Zhang, J. L., Song, H., and Cao, Y. X. (2019). Boosting the biosynthesis of betulinic acid and related triterpenoids in Yarrowia lipolytica via multimodular metabolic engineering. Microb. Cell Factories 18:77. doi: 10.1186/s12934-019-1127-8
Johnson, E. S., Ma, P. C., Ota, I. M., and Varshavsky, A. (1995). A proteolytic pathway that recognizes ubiquitin as a degradation signal. J. Biol. Chem. 270, 17442–17456. doi: 10.1074/jbc.270.29.17442
Kretzschmar, A., Otto, C., Holz, M., Werner, S., Hubner, L., and Barth, G. (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
Krull, S., Thyberg, J., Bjorkroth, B., Rackwitz, H. R., and Cordes, V. C. (2004). Nucleoporins as components of the nuclear pore complex core structure and Tpr as the architectural element of the nuclear basket. Mol. Biol. Cell 15, 4261–4277. doi: 10.1091/mbc.e04-03-0165
Kumar, A. (2016). Multipurpose Transposon-Insertion Libraries in Yeast. Cold Spring Harb. Protoc. 6, 499–502. doi: 10.1101/pdb.top080259
Langmead, B., and Salzberg, S. L. (2012). Fast gapped-read alignment with bowtie 2. Nat. Methods 9, 357–359. doi: 10.1038/nmeth.1923
Levy-Strumpf, N., Deiss, L. P., Berissi, H., and Kimchi, A. (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
Li, H., and Alper, H. S. (2016). Enabling xylose utilization in Yarrowia lipolytica for lipid production. Biotechnol. J. 11, 1230–1240. doi: 10.1002/biot.201600210
Li, B., and Dewey, C. 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
Li, H., Tang, X., Yang, X., and Zhang, H. (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
Lieber, M. 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
Lincecum, T. L., Tukalo, M., Yaremchuk, A., Mursinna, R. S., Williams, A. M., Sproat, B. S., et al. (2003). Structural and mechanistic basis of pre- and Posttransfer editing by Leucyl-tRNA Synthetase. Mol. Cell 11, 951–963. doi: 10.1016/S1097-2765(03)00098-4
Ling, J., Reynolds, N., and Ibba, M. (2009). Aminoacyl-tRNA synthesis and translational quality control. Annu. Rev. Microbiol. 63, 61–78. doi: 10.1146/annurev.micro.091208.073210
Liu, I. C., Chiu, S. W., Lee, H. Y., and Leu, J. Y. (2012). The histone deacetylase Hos2 forms an Hsp42-dependent cytoplasmic granule in quiescent yeast cells. Mol. Biol. Cell 23, 1231–1242. doi: 10.1091/mbc.e11-09-0752
Liu, Y., Jiang, X., Cui, Z., Wang, Z., Qi, Q., and Hou, J. (2019). Engineering the oleaginous yeast Yarrowia lipolytica for production of alpha-farnesene. Biotechnol. Biofuels 12:296. doi: 10.1186/s13068-019-1636-z
Liu, X., Liu, M., Zhang, J., Chang, Y., Cui, Z., Ji, B., 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
Liu, L., Markham, K., Blazeck, J., Zhou, N., Leon, D., Otoupal, P., 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
Liu, L., Pan, A., Spofford, C., Zhou, N., and Alper, H. 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
Longmuir, S., Akhtar, N., and MacNeill, S. A. (2019). Unexpected insertion of carrier DNA sequences into the fission yeast genome during CRISPR-Cas9 mediated gene deletion. BMC. Res. Notes 12:191. doi: 10.1186/s13104-019-4228-x
Madzak, C., Treton, B., and Blanchin-Roland, S. (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.
Maoz, T., Koren, R., Ben-Ari, I., and Kleinberger, T. (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
Matthaus, F., Ketelhot, M., Gatter, M., and Barth, G. (2014). Production of lycopene in the non-carotenoid-producing yeast Yarrowia lipolytica. Appl. Environ. Microbiol. 80, 1660–1669. doi: 10.1128/AEM.03167-13
Mayer, M. P., and Bukau, B. (2005). Hsp70 chaperones: cellular functions and molecular mechanism. Cell. Mol. Life Sci. 62, 670–684. doi: 10.1007/s00018-004-4464-6
McClellan, A. J., Endres, J. B., Vogel, J. P., Palazzi, D., Rose, M. D., and Brodsky, J. L. (1998). Specific molecular chaperone interactions and an ATP-dependent conformational change are required during posttranslational protein translocation into the yeast ER. Mol. Biol. Cell 9, 3533–3545. doi: 10.1091/mbc.9.12.3533
Mehler, A. H., and Mitra, S. 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
Mittelman, K., Ziv, K., Maoz, T., and Kleinberger, T. (2010). The cytosolic tail of the Golgi apyrase Ynd1 mediates E4orf4-induced toxicity in Saccharomyces cerevisiae. PLoS One 5:e15539. doi: 10.1371/journal.pone.0015539
Murray, M. G., and Thompson, W. F. (1980). Rapid isolation of high molecular weight plant DNA. Nucleic Acids Res. 8, 4321–4326. doi: 10.1093/nar/8.19.4321
Negrutskii, B. S., and Deutscher, M. 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
Nicaud, J. M., Madzak, C., van den Broek, P., Gysler, C., Duboc, P., Niederberger, P., 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
Nielsen, J., and Keasling, J. D. (2016). Engineering cellular metabolism. Cell 164, 1185–1197. doi: 10.1016/j.cell.2016.02.004
Onuma, K., Watanabe, A., Kanzaki, N., and Kubota, T. (2006). Association kinetics of wild- and mutant-type Ynd1p in relation to quality of grown crystals. J. Phys. Chem. B 110, 24876–24883. doi: 10.1021/jp0643146
Pain, V. M., and Standart, N. (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. RNA 7, 331–333. doi: 10.1017/S1355838201002540
Pannunzio, N. R., Watanabe, G., and Lieber, M. 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
Park, S. G., Schimmel, P., and Kim, S. (2008). Aminoacyl tRNA synthetases and their connections to disease. Proc. Natl. Acad. Sci. U. S. A. 105, 11043–11049. doi: 10.1073/pnas.0802862105
Park, Y., Yoon, S. K., and Yoon, J. 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
Pedelacq, J. D., Cabantous, S., Tran, T., Terwilliger, T. C., and Waldo, G. S. (2006). Engineering and characterization of a superfolder green fluorescent protein. Nat. Biotechnol. 24, 79–88. doi: 10.1038/nbt1172
Perez-Vargas, J., Romero, P., Lopez, S., and Arias, C. 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
Pfleger, B. F., Gossing, M., and Nielsen, J. (2015). Metabolic engineering strategies for microbial synthesis of oleochemicals. Metab. Eng. 29, 1–11. doi: 10.1016/j.ymben.2015.01.009
Phan, L., Schoenfeld, L. W., Valasek, L., Nielsen, K. H., and Hinnebusch, A. 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
Ranjha, L., Howard, S. M., and Cejka, P. (2018). Main steps in DNA double-strand break repair: an introduction to homologous recombination and related processes. Chromosoma 127, 187–214. doi: 10.1007/s00412-017-0658-1
Rapoport, T. A. (2007). Protein translocation across the eukaryotic endoplasmic reticulum and bacterial plasma membranes. Nature 450, 663–669. doi: 10.1038/nature06384
Richard, G. F., Kerrest, A., Lafontaine, I., and Dujon, B. (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
Robson, A., and Collinson, I. (2006). The structure of the sec complex and the problem of protein translocation. EMBO Rep. 7, 1099–1103. doi: 10.1038/sj.embor.7400832
Rodriguez, M. S., Dargemont, C., and Stutz, F. (2004). Nuclear export of RNA. Biol. Cell. 96, 639–655. doi: 10.1016/j.biolcel.2004.04.014
Rout, M. P., and Aitchison, J. D. (2001). The nuclear pore complex as a transport machine. J. Biol. Chem. 276, 16593–16596. doi: 10.1074/jbc.R100015200
Rout, M. P., Aitchison, J. D., Suprapto, A., Hjertaas, K., Zhao, Y., and Chait, B. 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
Schwartz, C., Cheng, J. F., Evans, R., Schwartz, C. A., Wagner, J. M., Anglin, S., 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
Schwartz, C. M., Hussain, M. S., Blenner, M., and Wheeldon, I. (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
Shiomi, D., Toyoda, A., Aizu, T., Ejima, F., Fujiyama, A., Shini, T., 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
Shrivastav, M., De Haro, L. P., and Nickoloff, J. A. (2008). Regulation of DNA double-strand break repair pathway choice. Cell Res. 18, 134–147. doi: 10.1038/cr.2007.111
Smith, D. M., Benaroudj, N., and Goldberg, A. (2006). Proteasomes and their associated ATPases: a destructive combination. J. Struct. Biol. 156, 72–83. doi: 10.1016/j.jsb.2006.04.012
Stephenson, K. (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
Sudfeld, C., Hubacek, M., Figueiredo, D., Naduthodi, M. I. S., van der Oost, J., Wijffels, R. 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
Sun, W., Bernard, C., van de Cotte, B., Van Montagu, M., and Verbruggen, N. (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
Tai, M., and Stephanopoulos, G. (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
Tarazona, S., Garcia-Alcalde, F., Dopazo, J., Ferrer, A., and Conesa, A. (2011). Differential expression in RNA-seq: a matter of depth. Genome Res. 21, 2213–2223. doi: 10.1101/gr.124321.111
Valasek, L. 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
van Opijnen, T., Bodi, K. L., and Camilli, A. (2009). Tn-seq: high-throughput parallel sequencing for fitness and genetic interaction studies in microorganisms. Nat. Methods 6, 767–772. doi: 10.1038/nmeth.1377
Verbeke, J., Beopoulos, A., and Nicaud, J. 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
Vladimirov, S. N., Ivanov, A. V., Karpova, G. G., Musolyamov, A. K., Egorov, T. A., Thiede, B., 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
Voges, D., Zwickl, P., and Baumeister, W. (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
Wagner, J. M., and Alper, H. 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
Wahlen, B. D., Oswald, W. S., Seefeldt, L. C., and Barney, B. 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
Wang, H. H., Isaacs, F. J., Carr, P. A., Sun, Z. Z., Xu, G., Forest, C. R., et al. (2009). Programming cells by multiplex genome engineering and accelerated evolution. Nature 460, 894–898. doi: 10.1038/nature08187
Wang, G., Shi, T., Chen, T., Wang, X., Wang, Y., Liu, D., 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
Wang, W., Wei, H., Knoshaug, E., Van Wychen, S., Xu, Q., Himmel, M. E., et al. (2016b). Fatty alcohol production in Lipomyces starkeyi and Yarrowia lipolytica. Biotechnol. Biofuels 9:227. doi: 10.1186/s13068-016-0647-2
Wang, G., Xiong, X., Ghogare, R., Wang, P., Meng, Y., and Chen, S. (2016a). Exploring fatty alcohol-producing capability of Yarrowia lipolytica. Biotechnol. Biofuels 9:107. doi: 10.1186/s13068-016-0512-3
Warner, J. R., Reeder, P. J., Karimpour-Fard, A., Woodruff, L. B., and Gill, R. T. (2010). Rapid profiling of a microbial genome using mixtures of barcoded oligonucleotides. Nat. Biotechnol. 28, 856–862. doi: 10.1038/nbt.1653
Wei, L. J., Cao, X., Liu, J. J., Kwak, S., Jin, Y. S., Wang, W., 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
Willis, R. M., Wahlen, B. D., Seefeldt, L. C., and Barney, B. 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. Biochemistry 50, 10550–10558. doi: 10.1021/bi2008646
Xu, P., Qiao, K., Ahn, W. S., and Stephanopoulos, G. (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
Xuan, J. W., Fournier, P., Declerck, N., Chasles, M., and Gaillardin, C. (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
Xue, Z., Sharpe, P. L., Hong, S. P., Yadav, N. S., Xie, D., Short, D. 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
Ye, R. W., Sharpe, P. L., and Zhu, Q. (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
Young, J. C. (2010). Mechanisms of the Hsp70 chaperone system. Biochem. Cell Biol. 88, 291–300. doi: 10.1139/O09-175
Zhang, J.-L., Cao, Y.-X., Peng, Y.-Z., Jin, C.-C., Bai, Q.-Y., Zhang, R.-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
Zhang, J. L., Peng, Y. Z., Liu, D., Liu, H., Cao, Y. X., Li, B. Z., et al. (2018). Gene repression via multiplex gRNA strategy in Y. lipolytica. Microb. Cell Factories 17:62. doi: 10.1186/s12934-018-0909-8
Zhou, Y. J., Buijs, N. A., Zhu, Z., Qin, J., Siewers, V., and Nielsen, J. (2016). Production of fatty acid-derived oleochemicals and biofuels by synthetic yeast cell factories. Nat. Commun. 7:11709. doi: 10.1038/ncomms11709
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
Edited by:
Haichun Gao, Zhejiang University, ChinaReviewed by:
Huihui Fu, Institute of Oceanology (CAS), ChinaJingwen Zhou, Jiangnan University, China
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