- 1ERA Chair in Gas Fermentation Technologies, Institute of Technology, University of Tartu, Tartu, Estonia
- 2LanzaTech Inc., Skokie, IL, United States
Gas fermentation has emerged as a sustainable route to produce fuels and chemicals by recycling inexpensive one-carbon (C1) feedstocks from gaseous and solid waste using gas-fermenting microbes. Currently, acetogens that utilise the Wood-Ljungdahl pathway to convert carbon oxides (CO and CO2) into valuable products are the most advanced biocatalysts for gas fermentation. However, our understanding of the functionalities of the genes involved in the C1-fixing gene cluster and its closely-linked genes is incomplete. Here, we investigate the role of two genes with unclear functions—hypothetical protein (hp; LABRINI_07945) and CooT nickel binding protein (nbp; LABRINI_07950)—directly adjacent and expressed at similar levels to the C1-fixing gene cluster in the gas-fermenting model-acetogen Clostridium autoethanogenum. Targeted deletion of either the hp or nbp gene using CRISPR/nCas9, and phenotypic characterisation in heterotrophic and autotrophic batch and autotrophic bioreactor continuous cultures revealed significant growth defects and altered by-product profiles for both ∆hp and ∆nbp strains. Variable effects of gene deletion on autotrophic batch growth on rich or minimal media suggest that both genes affect the utilisation of complex nutrients. Autotrophic chemostat cultures showed lower acetate and ethanol production rates and higher carbon flux to CO2 and biomass for both deletion strains. Additionally, proteome analysis revealed that disruption of either gene affects the expression of proteins of the C1-fixing gene cluster and ethanol synthesis pathways. Our work contributes to a better understanding of genotype-phenotype relationships in acetogens and offers engineering targets to improve carbon fixation efficiency in gas fermentation.
1 Introduction
The high concentration of greenhouse gases (GHG) in the atmosphere is leading to potentially irreversible climate change threatening Earth’s sustainability. Fossil fuel combustion for energy, heat, and transportation is widely recognised as the primary source of GHG (mainly CO2) from human activities (Sixth Assessment Report IPCC, 2021). Likewise, increasing population leads to higher volume of generated waste [e.g., plastic waste and municipal solid waste (MSW)] that is negatively impacting ecosystems. Although decarbonising energy generation is important to lowering GHG emissions, we also need to adopt sustainable technologies for the production of fuels and chemicals as these will mostly stay carbon-based for the foreseeable future. The advancement of such technologies will be vital in achieving the United Nations Sustainable Development Goal (SDG) 7—Affordable and Clean Energy and SDG 13—Climate Action by the target year 2050.
Gas fermentation has emerged as an attractive route for sustainable production of low-carbon fuels and chemicals from gaseous one-carbon (C1) waste feedstocks (e.g., industrial waste gases (CO2, CO, and CH4) and syngas [from gasified biomass or MSW (CO, H2, and CO2)] using gas-fermenting microbes (Liew et al., 2016; Redl et al., 2017; Fackler et al., 2021; Pavan et al., 2022). Anaerobic acetogens are the most advanced biocatalysts for gas fermentation as they can naturally use gas as their sole energy and carbon source (Drake et al., 2006). Acetogens use the most energy-efficient CO2 fixation pathway, the Wood-Ljungdahl pathway (WLP), to fix carbon oxides (CO2 and/or CO) into many metabolic products via acetyl-CoA. Notably, the model-acetogen Clostridium autoethanogenum is being used as a commercial-scale cell factory for ethanol production (Köpke and Simpson, 2020) and has been demonstrated at pilot scale to produce two valuable industrial solvents: acetone and isopropanol (Liew et al., 2022).
The WLP is considered to be the first biochemical pathway on Earth (Russell & Martin, 2004) and it plays a vital role in global carbon cycles by converting ∼20% of CO2 on Earth each year to billions of tons of acetate (Drake et al., 2006; Ljungdahl, 2009). The genes encoding the enzymes involved in the WLP belong to a large C1-fixing gene cluster (Figure 1) that has been annotated as CAETHG_1606–1621 (Brown et al., 2014). The authors proposed that the cluster contains 16 genes, including the WLP central enzyme—the bifunctional CO dehydrogenase/acetyl-CoA synthase (CODH/ACS; CAETHG_1620-21/1608) complex—that catalyses CO oxidation/CO2 reduction and acetyl-CoA synthesis (Ragsdale and Pierce, 2008) and is essential for autotrophic growth (Liew et al., 2016). The cluster also encodes the pathway’s most abundant but slowest enzyme in C. autoethanogenum (Valgepea et al., 2022): the formyl-THF ligase (Fhs; CAETHG_1618). Though the biochemistry of the WLP is well-described (Ragsdale, 2008; Ragsdale and Pierce, 2008; Can et al., 2014), our understanding of the functionalities of genes within the proposed C1-fixing gene cluster is incomplete. For instance, functions of genes CAETHG_1607 and 1613 are unclear in C. autoethanogenum. In fact, the C), our understanding of the functionalities of genes within the proposed C1-fixing gene cluster is incomplete. For instance, functions of genes CAETHG_1607 and 1613 are unclear in C. autoethanogenum. In fact, the C), our understanding of the functionalities of genes within the proposed C1-fixing gene cluster is incomplete. For instance, functions of genes CAETHG_1607 and 1613 are unclear in C. autoethanogenum. In fact, the C1-fixing gene cluster might be larger as several adjacent genes with unknown functions show similarly high gene and protein expression levels (Valgepea et al., 2017a-fixing gene cluster might be larger as several adjacent genes with unknown functions show similarly high gene and protein expression levels (Valgepea et al., 2017a-fixing gene cluster might be larger as several adjacent genes with unknown functions show similarly high gene and protein expression levels (Valgepea et al., 2017a; 2022). Two of such genes, found right downstream of the currently annotated cluster, are hypothetical protein (hp; CAETHG_1604) and CooT family nickel binding protein (nbp; CAETHG_1605). Both nbp and hp genes are expressed highly in C. autoethanogenum and their homologs accompany the WLP genes of other commonly studied acetogens like Clostridium ljungdahlii, Clostridium drakei, Clostridium carboxidivorans, Eubacterium limosum, as well as thermophilic acetogens Acetobacterium woodii and Thermoanaerobacter kivui. Despite the clear genomic conservation of the hp and nbp genes among acetogens, their function has not yet been described. While nickel binding is predicted for Nbp based on gene sequence annotation, and histidine-rich amino acid sequence of the Hp protein suggests a metal-binding activity, the validity of these functional predictions remains to be demonstrated.
FIGURE 1. C1-fixing gene cluster in Clostridium autoethanogenum. The gene cluster contains 16 genes, including five genes with unconfirmed biochemical functions (CAETHG_1606; 1607; 1612; 1613; 1619). Gene names and annotations based on Valgepea et al., 2022. Figure created using BioRender.com.
Determination of gene functionalities in acetogens through genotype-phenotype studies has been limited mostly due to slow, labour-intensive, and challenging workflows for constructing genetically modified strains (Pavan et al., 2022). However, gene-editing tools for acetogens have rapidly improved recently (Jin et al., 2020; Bourgade et al., 2021; Pavan et al., 2022). We thus aimed to take advantage of these developments to shed light on the functionalities of the hp (CAETHG_1604) and the CooT family nickel binding gene nbp (CAETHG_1605), referred to as LABRINI_07945 for hp and LABRINI_07950 for nbp based on the C. autoethanogenum base strain “LAbrini” used in this work (Ingelman et al., 2023).
Here, we used CRISPR/nCas9 in C. autoethanogenum for the targeted deletion of two genes adjacent to the C1-fixing gene cluster with unclear functions: named hp (LABRINI_07945) and nbp (LABRINI_07950). Phenotypic characterisation in heterotrophic and autotrophic batch and autotrophic bioreactor continuous cultures revealed significant growth defects and altered by-product profiles for both deletion strains compared to the base strain LAbrini. Additionally, proteomics showed that disruption of either gene affects the expression of proteins of the C1-fixing gene cluster and ethanol synthesis pathways. Our findings contribute to a better understanding of genotype-phenotype relationships which is needed to advance the metabolic engineering of acetogen factories for gas fermentation.
2 Materials and methods
2.1 Bacterial strains, growth media and cultivation conditions
A derivate of C. autoethanogenum DSM 10061 strain deposited in the German Collection of Microorganisms and Cell Cultures (DSMZ) named “LAbrini” (Ingelman et al., 2023) was stored as a glycerol stock at −80°C and used as the base strain for genetic engineering. Escherichia coli strains NEB Turbo and NEB Express (New England BioLabs) were used for cloning, plasmid assembly, and propagation. For transformant selection, the cultivation of E. coli was done in Lysogeny Broth (LB) medium supplemented with ampicillin (100 μg/mL).
For C. autoethanogenum cultivation on solid medium, YTF agar containing 10 g/L yeast extract (YE), 16 g/L tryptone, and 5 g/L fructose (Liew et al., 2017) was used and supplemented with 4 μg/mL clarithromycin if needed (see below). For liquid batch cultures, cells were grown in chemically defined PETC-MES medium (Valgepea et al., 2017b) with or without YE (1.5 g/L) and supplemented with 0.4 g/L of cysteine-HCl∙H2O as the reducing agent while the carbon source was either 5 g/L of fructose or syngas (50% CO, 20% H2, 20% CO2, and 10% Ar; AS Eesti AGA). The PETC-MES medium contained per litre: 1 g NH4Cl, 0.1 g KCl, 0.2 g MgSO4·7H2O, 0.2 g KH2PO4, 0.02 g CaCl2, 20 g 2-(N-morpholino)ethanesulfonic acid (MES), 0.25 g sodium acetate, 0.05 g Fe(SO4)·7H2O, 0.05 g nitriolotriacetic acid, 0.5 mL of 2 g/L resazurin, 10 mL trace metal solution (TMS) and 10 mL of Wolfe’s vitamin solution. The TMS solution contained per litre: 2 g nitrilotriacetic acid, 1 g MnSO4·H2O, 0.8 g Fe(SO4)·7H2O, 0.2 g CoCl2·6H2O, 0.2 mg ZnSO4·7H2O, 0.02 g CuCl2·2H2O, 0.02 g NaMoO4·2H2O, 0.03 g Na2SeO3·5H2O, 0.02 g NiCl2·6H2O and 0.02 g Na2WO4·2H2O. The Wolfe’s vitamin solution contained per litre: 2 mg biotin, 2 mg folic acid, 10 mg pyridoxine hydrochloride, 5 mg thiamine-HCl, 5 mg riboflavin, 5 mg nicotinic acid, 5 mg calcium pantothenate, 0.1 mg vitamin B12, 5 mg 4-aminobenzoic acid and 5 mg thioctic acid. The medium was prepared anaerobically, and the pH was adjusted to 5.7. Batch fermentations were performed in 250 mL Schott bottles with 50 mL liquid volume in biological triplicates incubated horizontally at 37°C with orbital shaking at 120 RPM under strictly anaerobic conditions unless stated otherwise. For autotrophic batch cultures, bottle headspace was pressurised to 140 kPa with syngas. Cultures were sampled during the exponential growth phase for growth characterisation and extracellular metabolome analysis to determine the maximum specific growth rate (μmax) and production yields of acetate, ethanol, and 2,3-butanediol. At least five time points were used for μmax calculation resulting in correlation coefficients R2 ≥ 0.98 between cultivation time and natural logarithm (ln) of the culture’s optical density (OD600). Product yields (mmol/gram of dry cell weight [DCW]) were estimated using at least three time points by linear regression between the respective product concentration (mmol/L) and biomass concentration (gDCW/L) with R2 ≥ 0.92. Product yields were estimated across time-points that showed no consumption of any of the three quantified products (i.e., products did not serve as substrates).
Continuous chemostat fermentations were carried out as described in detail before (Valgepea et al., 2017a). Shortly, all three C. autoethanogenum strains were grown in bioreactors on syngas under strictly anaerobic conditions at 37°C and pH 5 in chemically defined medium (without YE) at a dilution rate of 1 day−1. Chemostat continuous cultures were performed in 1.4 L Multifors bioreactors (Infors AG) at a working volume of 750 mL connected to a Hiden HPR-20-QIC mass spectrometer (Hiden Analytical) for online high-resolution off-gas analysis. Steady-state results were collected after OD, gas uptake and production rates had been stable for at least 3–5 working volumes.
2.2 Genetic engineering
2.2.1 Genomic and plasmid DNA manipulations
Routine C. autoethanogenum genomic DNA extraction for PCR diagnostics was done with PureLink Genomic DNA extraction kit (Invitrogen™, Thermo Fischer Scientific), while plasmid DNA was extracted from E. coli using FavorGen FavorPrep™ Plasmid DNA Extraction Mini Kit (Favorgen Biotech Corp). All purified genomic and plasmid DNA fragments were quantified using the NanoDrop™ 1000 Spectrophotometer (Thermo Fischer Scientific). PCR for amplification of DNA fragments for sequencing and cloning, routine screening, and analytical procedures was performed using the Phusion™ High-Fidelity DNA Polymerase (Thermo Fischer Scientific). All primers used in this study (see Supplementary Table S1) were designed using the NetPrimer software (PREMIER Biosoft International) and synthesised by Integrated DNA Technologies. PCR products and DNA fragments were purified from agarose gels with FavorPrep Gel/PCR Purification kit (Favorgen Biotech Corp).
2.2.2 Assembly of CRISPR/nCas9 plasmids
A CRISPR/nCas9 gene editing plasmid for C. autoethanogenum was constructed based on the pNickClos2.0 plasmid (Addgene plasmid #73228) (Li et al., 2016) with the following modifications to the backbone sequence: origin of transfer and the required traJ coding sequence from pFX01 (Xia et al., 2020) was added to allow potential transfer of the resulting plasmids into recipient cells via conjugation. Additionally, one of the two NotI restriction sites was removed to enable convenient cloning of the homology arms (HAs) using NotI and XhoI sites.
To specify the target site for the nCas9 enzyme, single guide RNAs (sgRNAs) were designed using CRISPR RGEN Tools (Park et al., 2015). We used the Basic Local Alignment Search Tool (Altschul et al., 1990) to cross-check that the selected spacer sequences would lead to unique cuts in the genome, thus minimising off-target cleavages. Inverse PCR (iPCR) was employed to clone the unique sgRNA for each target gene into the backbone plasmid to create intermediate plasmids: pGFT021 for LABRINI_07945 (hp) and pGFT022 for LABRINI_07950 (nbp). The two HAs were obtained by PCR amplification of the 1 kb regions upstream and downstream of each target gene from the C. autoethanogenum genome using various primer pairs (Supplementary Table S1). The derived 5′ and 3′ HA fragments were fused using overlap extension PCR as described by Bryksin and Matsumura (2013) and moved into the appropriate sgRNA-containing vectors using the same technique or via restriction cloning. The resulting plasmids containing all necessary parts for CRISPR/nCas9-mediated gene deletion were named pGFT036 for hp and pGFT048 for nbp KO. See Supplementary File S1 (pGFT036) and Supplementary File S2 (pGFT048) for plasmid maps. Plasmids constructed in this study are listed in Supplementary Table S2.
2.2.3 Plasmid transformation into C. autoethanogenum using electroporation
Before electroporation, plasmids were propagated and isolated from the E. coli NEB Express strain to lose the Dcm methylation pattern, which has been shown to be targeted by the native type IV restriction system of C. autoethanogenum (Woods et al., 2019). Preparation of electrocompetent C. autoethanogenum cells and transformation using electroporation was performed as previously described (Leang et al., 2013). Electroporated cells were recovered in 10 mL YTF medium for 18–24 h at 37°C with 120 RPM of agitation. Recovery cultures were mixed with 1.5% molten agar YTF supplemented with clarithromycin (4 μg/mL). Transformation plates were incubated inside Oxoid™ AnaeroJars™ (Thermo Fischer Scientific) at 37°C until colonies became visible (in less than 10 days).
2.2.4 Transformant colony screening and plasmid curing
Colonies were first inoculated into liquid YTF medium containing clarithromycin and then plated on selective YTF agar to isolate individual transformant colonies. Next, transformant colonies were randomly chosen for colony PCR screening to verify C. autoethanogenum (primers acsB-F and -R), plasmid presence (primers oriT-F and -R), and successful gene deletion. Primers hp-F, hp-R and, nbp-F, nbp-R were used to confirm deletions of hp and nbp genes, respectively. PCR products obtained with the gene deletion control primers were purified and verified by Sanger sequencing.
Plasmid-carrying C. autoethanogenum colonies with target gene deletions were then subcultured three or more times in 10 mL YTF medium without the antibiotic before plating on non-selective YTF agar to isolate individual colonies. Next, colonies were randomly chosen and patch-plated on both YTF agar without and with clarithromycin (4 μg/mL). Plasmid loss in colonies that failed to grow in the presence of clarithromycin was confirmed by colony PCR using primers oriT-F and oriT-R as above. In total, 10 out of 10 colonies for hp gene deletion and 1 out of 30 colonies for nbp gene deletion had lost their plasmid. The plasmid-cured colonies were then inoculated into liquid YTF medium for outgrowth and preparation of glycerol stocks, stored as Δhp and Δnbp.
2.3 Analytical methods
2.3.1 Biomass concentration analysis
Biomass concentration (gDCW/L) was estimated by measuring culture OD at 600 nm and using the correlation coefficient of 0.23 between OD and DCW using the methodology established before (Peebo et al., 2014).
2.3.2 Extracellular metabolome analysis
Extracellular metabolome analysis from filtered broth samples was performed as described before (Valgepea et al., 2017b). We note that cells produced 2R,3R-butanediol.
2.3.3 Bioreactor off-gas analysis
Bioreactor off-gas analysis for the determination of specific gas uptake (CO and H2) and production rates (CO2, ethanol) (mmol/gDCW/h) have been described before (Valgepea et al., 2017a). The Faraday Cup detector monitored the intensities of H2, CO, ethanol, H2S, Ar, and CO2 at 2, 14, 31, 34, 40, and 44 amu, respectively.
2.3.4 Carbon balance analysis
Carbon recoveries and balances were determined as described before (Valgepea et al., 2017a). Briefly, carbon balancing of substrate carbon (CO, cysteine) between growth products (acetate, ethanol, 2,3-BDO, CO2, biomass) was calculated as C-mol fraction of each product from total C-mol carbon uptake. Carbon recoveries were calculated as the fraction of summed C-mol products from total C-mol substrates. Ethanol stripping and the total soluble CO2 fraction in culture broth were also considered to achieve more accurate carbon balancing.
2.4 Proteome analysis
2.4.1 Sample preparation
Proteome analysis was carried out for four biological replicate cultures for each of the three strains: LAbrini, Δhp, and Δnbp. C. autoethanogenum chemostat cultures were sampled for proteomics by immediately pelleting 2 mL of the culture by centrifugation (20,000 × g for 1 min at 4°C) and stored at −80°C until analysis. Cell pellets were suspended with 10 volumes (relative to pellet volume) of chaotrope-based lysis buffer (6 M guanidine-HCl, 100 mM Tris-HCl pH 8.0, 50 mM dithiothreitol), heated at 95°C for 10 min and sonicated with the Bioruptor sonicator (Diagenode) for 15 cycles (30 s ON, 60 s OFF; “High” setting) at 4°C. Next, 0.5 volumes (relative to lysate volume) of Silibeads TypZY-s 0.4–0.6 mm (Sigmund Lindner) beads were added to the lysate and bead beating was carried out with the FastPrep24 device (MP Biomedicals) by 2 × 40 s at 6 m/s. Next, samples were centrifuged at 17,000 × g for 5 min, and the supernatant was separated from the beads and transferred to a new tube. A small aliquot of the lysate was then precipitated with trichloroacetic acid-deoxycholate (TCA-DOC) precipitation and the protein concentration was determined with the Micro BCA assay (Thermo Fisher Scientific). Next, 20 µg of lysate protein was alkylated with 100 mM chloroacetamide by incubating for 1 h in the dark at room temperature. Samples were then processed, and proteins were digested to peptides by the SP3 protocol as described elsewhere (Hughes et al., 2018).
2.4.2 LC-MS/MS analysis
Sample order for LC-MS/MS analysis was randomised to avoid batch effects and blanks were flanked to minimise carry-over. 1 μg of peptides for each sample was injected to an Ultimate 3500 RSLCnano system (Dionex) using a 0.3 × 5 mm trap-column (5 µm C18 particles, Dionex) and an in-house packed (3 µm C18 particles, Dr Maisch) analytical 50 cm × 75 µm emitter-column (New Objective). Columns were operated at 45°C. Peptides were eluted at 300 nL/min with an 8%–45% B 90 min gradient (buffer B: 80% acetonitrile + 0.1% formic acid, buffer A: 0.1% formic acid) to a Q Exactive HF (Thermo Fisher Scientific) mass spectrometer (MS) using a nano-electrospray source (spray voltage of 2.5 kV in positive mode). The MS was operated using a data-independent acquisition (DIA) approach (Gillet et al., 2012) with variable isolation windows over a range of 400–1,100 m/z. Briefly, one full range 400–1,100 m/z MS1 scan was collected at a resolution setting of 60,000 (max ion injection time of 60 ms, max of 3e6 ions), followed by 25 overlappings (overlap of 1 m/z) DIA isolation windows as follows (start‒end m/z, isolation width): 399–472, 25; 471–510, 20; 509–608, 15; 607–703, 20; 702–790, 30; 789–868, 40; 867–966, 50; and 965–1101, 136. Each DIA scan was collected at a resolution setting of 30,000 (max ion injection time of 41 ms, max of 3e6 ions), the default charge state was set to +3 and normalised collision energy to 27. The overall method cycle time was approximately 2.3 s.
2.4.3 DIA MS data analysis and differential protein expression analysis
DIA MS data was analysed using version 1.8 of the software suite DIA-NN (Demichev et al., 2019) with default settings. The spectral library was generated in silico using smart profiling from the protein sequence database of C. autoethanogenum LAbrini (NCBI Genbank CP110420) (Ingelman et al., 2023) with a precursor of 1% false discovery rate (FDR) filter. The following settings were used for analysis: full trypsin specificity with 1 missed cleavage allowed for peptides with a length of 7–30 AAs; fixed modifications of cysteine carbamidomethylation and methionine N-terminal excision; precursor charge range 1–4 and m/z range 400–1,100; fragment m/z range 200–1,800; double-pass mode for neural network classifier; Robust LC (high precision) for quantification strategy; and RT-dependent cross-run normalization. Label-free protein quantification and normalisation were performed within the DIA-NN workflow using the MaxLFQ method (Cox et al., 2014). We confidently quantitated 28,752 peptides and 2,127 proteins across all samples, and 27,358 peptides and 2,107 proteins on average within each sample after the removal of shared peptides from the analysis.
Protein expression fold-changes (FC) with q-values between Δhp and LAbrini, and Δnbp and LAbrini were determined using the software Perseus (Tyanova et al., 2016) with Student’s T-test. Only proteins with at least two peptides in all 12 samples (1,835) analysed were used to ensure higher quantification accuracy. Proteins were considered to be differentially expressed by an FC > 1.5 and a q-value < 0.05 after FDR correction (Benjamini and Hochberg, 1995). Differentially expressed proteins are presented in Supplementary Table S3. Proteomics data have been deposited to the ProteomeXchange Consortium (http://proteomecentral.proteomexchange.org) via the PRIDE partner repository (Perez-Riverol et al., 2022) with the data set identifier PXD040151.
3 Results and discussion
3.1 Single-gene deletion of hp and nbp in C. autoethanogenum using CRISPR/nCas9 knock-out plasmids
To study the role of two genes adjacent to the C1-fixing gene cluster (Figure 1) with unclear functions, LABRINI_07945 (hp) and LABRINI_07950 (nbp), each gene was targeted for deletion using CRISPR/nCas9. Plasmids pGFT036 and pGFT048 carrying required components (Figure 2A) for hp and nbp deletion, respectively, were transferred into C. autoethanogenum by electroporation that yields transformation efficiencies up to 200 CFU/μg DNA in our laboratory. We obtained 40 and 2 colonies for hp and nbp deletion, respectively, and colony PCR screening of 6 colonies for hp and 2 for nbp deletion confirmed C. autoethanogenum genotype and the presence of our relatively large plasmid (∼12.5 kb) in all colonies. The gene deletion control PCR showed that all the hp colonies screened had a mixed base strain/knock-out (KO) genotype evidenced by the presence of both the 2444-bp fragment of the base strain allele and the predicted 2103-bp fragment of the KO allele (Figure 2B). Six rounds of liquid selective sub-culturing of the mixed colonies was needed to obtain 1 colony (from 20 screened) with a clean hp KO genotype (Figure 2C). In contrast, the 2 nbp colonies (nbp-C1 and nbp-C2) obtained showed a 1444-bp fragment of the desired KO allele, i.e., a clean KO genotype (Figure 2B). Sanger sequencing confirmed both hp and nbp gene deletions in selected colonies and the presence of a PstI site, proving that the provided plasmid-based DNA repair template, in the form of the fused 1 kb HA interspaced with the PstI site, was used during homologous recombination to repair the nick from nCas9 (Figure 2C). Successful deletion of both genes was also confirmed by proteome analysis (see below). Finally, the obtained Δhp and Δnbp cells were cured of the plasmid before proceeding with phenotypic characterisation.
FIGURE 2. Successful deletion of target genes using CRISPR/nCas9 KO plasmids in C. autoethanogenum. (A) Representative plasmid map and an overview of CRISPR/nCas9 mediated gene deletion mechanism. The yellow star indicates the genomic recognition site of the Cas9-sgRNA complex. (B) PCR Screening of transformant colonies using control primers to amplify gene regions. (C) Confirmation of hp (LABRINI_07945) and nbp (LABRINI_07950) deletion via Sanger Sequencing. PC, positive control (non-transformed C. autoethanogenum); NC, negative control (no template); LAbrini, C. autoethanogenum base strain allele; KO, knock-out/deletion allele; SSB, single-strand break. Figure partly created using BioRender.com.
3.2 Deletion of hp and nbp genes confer growth defects in rich-medium batch cultures
Next, we grew Δhp and Δnbp strains in batch cultures on PETC-MES medium with YE (i.e., rich medium) in heterotrophic (fructose) and autotrophic (syngas: CO+H2+CO2) conditions to investigate potential phenotypic effects of hp and nbp deletions. Interestingly, the maximum specific growth rate (µmax) of both mutant strains was significantly impaired compared to the C. autoethanogenum base strain “LAbrini” on both fructose and syngas (Figure 3A). The µmax for the Δhp strain was 0.08 ± 0.001 (average ± standard deviation) and 0.05 ± 0.001 on fructose and syngas, respectively, while for the Δnbp strain, the µmax was 0.07 ± 0.002 on fructose and 0.06 ± 0.001 on syngas (Figure 3A). The LAbrini base strain showed µmax values of 0.10 ± 0.001 on fructose and 0.08 ± 0.001 on syngas. Thus, µmax for deletion strains were substantially lower compared to base strain: for Δhp, 21% and 37% lower on fructose and syngas and for Δnbp, 34% and 26% lower on fructose and syngas, respectively (Figure 3A). Interestingly, Δhp had a more pronounced growth defect on syngas, while Δnbp was affected more on fructose. The µmax values for both growth conditions were also statistically different between deletion strains. In addition to the lower µmax, Δhp strain exhibited a longer lag phase and slightly lower peak OD on syngas compared to the base strain LAbrini (Figure 3B). The growth profile and peak OD of the Δnbp strain, however, were similar to LAbrini despite the lower µmax values (Figure 3C). Altogether, our results demonstrate that the deletion of either the hp or the nbp gene confers significant growth defects for hetero- and autotrophic rich medium batch cultures, suggesting a notable role for both genes in the physiology of C. autoethanogenum.
FIGURE 3. Batch growth of Δhp, Δnbp, and LAbrini strains on PETC-MES with YE. (A) Maximum specific growth rate (µmax) in heterotrophic (fructose) and autotrophic (syngas) growth conditions. Data are average ± standard deviation between three biological replicates. Asterisk denotes values statistically different according to t-test (p-value ≤ 0.01). (B, C) Autotrophic growth profiles on syngas. LAbrini, C. autoethanogenum base strain; YE, yeast extract.
3.3 Deletion of hp and nbp show no effects on growth in minimal medium autotrophic batch cultures
We also investigated the growth of all three strains in PETC-MES medium without YE to assess the effects of hp and nbp deletions on autotrophic growth on syngas in minimal medium batch cultures. Intriguingly, there was no difference in the µmax values between the three strains: 0.08 ± 0.004, 0.08 ± 0.005, and 0.08 ± 0.003 for the Δhp, Δnbp, and LAbrini strains, respectively (Figure 4). Although one would expect slower growth in minimal medium, Δhp and Δnbp strains grew faster in minimal medium compared to rich medium, whereas LAbrini showed similarly high µmax values for autotrophic growth of acetogens on both media. Like the µmax data, the three strains did not show significant variations in peak OD and growth curves (i.e., lag phase and entry into stationary phase) (Figure 4B). Thus, both genes seem to play a role in C. autoethanogenum for utilisation of complex nutrients due to significant growth defects in rich medium with no phenotypic effects in minimal medium for Δhp and Δnbp strains.
FIGURE 4. Autotrophic batch growth of Δhp, Δnbp, and LAbrini strains in PETC-MES without YE. (A) Maximum specific growth rates (µmax) on syngas. Data are average ± standard deviation between two biological replicates. (B) Growth profiles on syngas. LAbrini, C. autoethanogenum base strain; YE, yeast extract.
3.4 Deletion of hp and nbp alter growth by-product yields in autotrophic batch cultures
We determined the production yields (mmol/gDCW) of C. autoethanogenum growth by-products — ethanol, acetate, and 2,3-butanediol (2,3-BDO) — during the above-described autotrophic batch cultures to check if these are affected by the deletion of either the hp or the nbp gene. Remarkably, the Δhp strain showed a >5-fold higher acetate yield (550 ± 45 mmol/gDCW) relative to LAbrini (92 ± 7 mmol/gDCW) during growth in rich medium (Figure 5A). The production of 2,3-BDO was also slightly higher for Δhp compared to LAbrini (10 ± 1 vs. 7 ± 1 mmol/gDCW). While the ethanol yield was not different for Δhp when compared to LAbrini (96 ± 5 vs. 102 ± 13 mmol/gDCW), ethanol production did not coincide with acetate production in the Δhp strain (data not shown). Namely, Δhp was producing acetate in the early exponential growth phase and started producing ethanol only in the late exponential phase.
FIGURE 5. Growth by-product yields and acetate-to-ethanol ratios of Δhp, Δnbp, and LAbrini strains in autotrophic batch cultures with and without YE. (A, B) Growth by-product yields. (C) Acetate/ethanol (Cmol). Data are average ± standard deviation between three biological replicates on PETC-MES + YE or two biological replicates on PETC-MES -YE. Asterisk denotes values statistically different according to t-test: *, p-value < 0.05; **, p-value < 0.01. LAbrini, C. autoethanogenum base strain; YE, yeast extract; gDCW, gram of dry cell weight; 2,3-BDO, 2,3-butanediol.
Also, the Δnbp strain showed significantly different product yields compared to LAbrini during growth in rich medium: 41% higher acetate (130 ± 0.2 vs. 92 ± 7 mmol/gDCW), 44% higher 2,3-BDO (10 ± 1 vs. 7 ± 1 mmol/gDCW), and 33% lower ethanol compared to LAbrini (65 ± 8 vs. 96 ± 6 to mmol/gDCW) (Figure 5A). In addition, both acetate and ethanol production yields were different between the two deletion strains. The substantially higher acetate production yield of the Δhp strain might indicate that deletion of the hp gene improves the uptake of YE components and their catabolism to acetate. However, this did not lead to faster growth, thus, acetate production was potentially not linked to ATP production.
We observed no difference in the product yields between the Δhp strain and LAbrini on minimal medium (Figure 5B), which is not surprising as growth was also comparable (Figure 4). However, despite similar growth for Δnbp and LAbrini strains on minimal medium (Figure 4), the Δnbp strain showed an 86% lower ethanol yield than LAbrini (9 ± 6 vs. 65 ± 1 to mmol/gDCW) while acetate and 2,3-BDO yields were similar (Figure 5B). Only the ethanol production yield was different between the deletion strains on minimal medium. Interestingly, while LAbrini maintained the acetate-to-ethanol ratio at ∼1 on both media, the ratio was higher for the Δhp strain (∼5 vs. ∼1) and lower for the Δnbp strain (∼2 vs. ∼11) on rich medium compared to minimal medium during autotrophic growth (Figure 5C). These ratios are in the range of those previously seen in autotrophic batch and continuous bioreactor cultures of C. autoethanogenum (Valgepea et al., 2017a; Valgepea et al., 2018; Ingelman et al., 2023). Our results imply that both hp and nbp genes play a notable role in carbon distribution between the two major by-products — acetate and ethanol — during autotrophic batch growth of C. autoethanogenum.
3.5 Deletion of nbp does not affect autotrophic growth in nickel-supplemented minimal medium batch cultures
Since Nbp is annotated as a nickel binding protein and deletion of nbp strongly affected growth in rich medium cultures (Figures 3, 5), we hypothesised that the effects could arise from a potentially lower capacity of the Δnbp strain to catabolise the higher nickel levels in rich medium, i.e., YE contains nickel. We thus tested the effect of nickel supplementation on autotrophic batch growth of Δnbp and LAbrini strains in PETC-MES without YE. We added 0.6 mg/L NiCl2·6H2O to the 0.2 mg/L in PETC-MES (0.8 mg/L final concentration) to facilitate potential effects on growth as a previous study showed variable influence of nickel in the range of 0.1–0.4 mg/L (NiCl2·6H2O) on heterotrophic growth of C. ljungdahlii (Mann et al., 2022).
We observed no difference between µmax values of Δnbp and LAbrini strains growing on syngas and in PETC-MES medium without YE with supplemented nickel (Figure 6A). These values were not statistically significantly (p-value > 0.05) different from respective values during growth without additional nickel (Figure 4A). Likewise, by-product yields between the two strains with nickel supplementation did not differ (Figure 6B). Both acetate and ethanol production yields with nickel addition did, however, differ from yields observed during growth without extra nickel (Figure 5B). Effect of nickel levels on acetate and ethanol production is also seen during heterotrophic growth of C. ljungdahlii (Mann et al., 2022). Our results showed that deletion of nbp did not affect autotrophic growth in nickel-supplemented minimal medium batch cultures. As several enzymes in the WLP are dependent on nickel, the potential role of nbp in C. autoethanogenum nickel metabolism remains to be further elucidated.
FIGURE 6. Effect of nickel supplementation on autotrophic (syngas) batch growth and growth by-product yields of Δnbp and LAbrini strains in PETC-MES without YE. (A) Maximum specific growth rates (µmax). (B) Growth by-product yields. Data are average ± standard deviation between three biological replicates. LAbrini, C. autoethanogenum base strain; YE, yeast extract; gDCW, gram of dry cell weight; 2,3-BDO, 2,3-butanediol.
3.6 Deletion of hp and nbp affect growth characteristics of steady-state chemostat cultures
We also compared the performance of the mutant strains to the LAbrini base strain in steady-state bioreactor continuous cultures as these yield high-quality physiological data (Adamberg et al., 2015; Valgepea et al., 2018) and are also the industrially relevant fermentation mode for acetogens (Köpke and Simpson, 2020). Cells were grown in minimal medium in syngas chemostats at a dilution rate of 1 day-1 and steady-state data at biomass concentrations ∼1.5 gDCW/L were collected for by-product production (ethanol, acetate, 2,3-BDO) and gas uptake and production (CO, H2, CO2) analysis.
While acetate production in minimal medium batch cultures was similar across all three strains, both Δhp and Δnbp strains showed lower specific acetate production rates (qace; mmol/gDCW/day) in chemostats (Figure 7A). Specific ethanol production rates (qEtOH) were also lower for both deletion strains compared to LAbrini (Figure 7A), whereas ethanol production differed only for the Δnbp strain in batch cultures (Figure 5). Production of 2,3-BDO (q2,3-BDO) was similar across the three strains in chemostats (Figure 7A). These steady-state data are consistent with above-described batch data by showing that both hp and nbp genes affect carbon distribution during autotrophic growth of the model-acetogen C. autoethanogenum. Still, differences in bioreactor chemostat and bottle batch data demonstrate the importance of testing strains in the industrially relevant fermentation mode for acetogens (Köpke and Simpson, 2020).
FIGURE 7. Steady-state growth characteristics of Δhp, Δnbp, and LAbrini strains in syngas-grown chemostats. (A) Specific by-product production rates (mmol/gDCW/day). (B) Specific gas uptake (CO, H2) and production (CO2) rates. Asterisk denotes values statistically different according to t-test (p-value < 0.05). (C) Carbon balances. Carbon recoveries were normalised to 100% to have a fair comparison of carbon distributions between different strains. Data are average ± standard deviation between four biological replicates. LAbrini, C. autoethanogenum base strain; EtOH, ethanol; 2,3-BDO, 2,3-butanediol; Ace, acetate; gDCW, gram of dry cell weight; q, specific rate.
Bioreactor off-gas analysis revealed that the specific CO uptake rates (qCO; mmol/gDCW/day) were slightly lower for both deletion strains compared to LAbrini (873 ± 19): by 8% for the Δhp strain (805 ± 19) and by 5% for the Δnbp strain (830 ± 15) (Figure 7B). At the same time, we detected no differences in the specific H2 uptake rates (qH2; mmol/gDCW/day) and CO2 production rates (qCO2; mmol/gDCW/day) among the three strains (Figure 7B). The slightly higher qCO2/qCO and qH2/qCO for both deletion strains indicates that cells needed to increase the generation of reducing power in the form of reduced ferredoxin as a response to gene deletion since both CO and H2 oxidation supply reduced ferredoxin (Bertsch and Müller, 2015; Pavan et al., 2022). Quantification of carbon flows from substrates to products showed that both Δhp and Δnbp strains diverted more carbon to biomass and CO2 with a concomitant lower carbon flux to ethanol and acetate (Figure 7C). The slightly lower acetate-to-ethanol ratios for Δhp (1.2) and Δnbp (1.1) compared to LAbrini (1.3) are consistent with higher qCO2/qCO for deletion strains since ethanol production generally results in a greater loss of carbon from CO as CO2 (Bertsch and Müller, 2015).
3.7 Proteome analysis reveals intriguing effects of deletion of hp and nbp genes
We used proteomics to determine the effects of hp and nbp gene deletions on protein expression in the above-described chemostat cultures of C. autoethanogenum. Our analysis quantified 93 differentially expressed proteins (fold-change > 1.5 and a q-value < 0.05) between Δhp and LAbrini, and 77 between Δnbp and LAbrini (Figure 8; Supplementary Table S3). Proteomics also confirmed the deletion of each gene in respective strains as no peptides for either protein in corresponding strains were detected. Additionally, the data showed no effect of deletion of either gene on expression of the protein encoded by the other gene, i.e., hp deletion on Nbp expression or nbp deletion on Hp expression. Intriguingly, deletion of either gene increased the expression of the methyltransferase (AcsE; 07970) and a component of the correnoid iron sulphur protein (AcsD; 07980) (Table 1) that are linked to a critical activity of the WLP: synthesis of acetyl-CoA. AcsE and AcsD supply the methyl group to the carbon monoxide dehydrogenase (AcsA; 08025) and acetyl-CoA synthase (AcsB; 07965) CODH/ACS enzyme complex (Lemaire and Wagner, 2021), that is essential for autotrophic growth of C. autoethanogenum (Liew et al., 2016). At the same time, the expression of AcsB was repressed in both deletion strains (Table 1). These are all among the top 10 most abundant proteins for autotrophic growth of C. autoethanogenum, except AcsA (Valgepea et al., 2022). These results show that the two genes adjacent to the C1-fixing cluster studied here—hp and nbp—affect the expression of key proteins for C1-fixation. Notably, the effect of hp and nbp on C1-fixation extends beyond the WLP as the expression of the most abundant CODH in C. autoethanogenum (Valgepea et al., 2022)—CooS1; 15015—was repressed in both mutant strains (Table 1). Further work is needed to understand the interplay between the hp, nbp, and the affected enzymes with a central role in the autotrophic growth of C. autoethanogenum.
FIGURE 8. Volcano plots showing differentially expressed proteins in Δhp and Δnbp strains compared to LAbrini. Differentially expressed proteins are with fold-change > 1.5 and a q-value < 0.05. Blue markers, downregulated proteins; red markers, upregulated proteins; grey, proteins not expressed differentially. See Supplementary Table S3 for data.
Deletion of our target genes also affected the expression of proteins involved in the synthesis of the most attractive native product of C. autoethanogenum: ethanol. In contrast to the lowered ethanol production in the deletion strains, proteome analysis revealed that expression of alcohol dehydrogenases (19700 and 18005) was elevated in both mutant strains, including the most abundant alcohol dehydrogenase in C. autoethanogenum (Valgepea et al., 2022): Adh4, 09125 (Table 1). Similarly, increased expression of both aldehyde ferredoxin oxidoreductases (AOR1, 00445 and AOR2, 00495), that couple ethanol synthesis to ATP production and the bifunctional acetaldehyde-CoA/alcohol dehydrogenase AdhE1 (18700) were measured in deletion strains. Such seeming discrepancies between either the transcripts or proteins linked to ethanol synthesis and the measured ethanol production have been detected in C. autoethanogenum before (Valgepea et al., 2017a; Valgepea et al., 2018). Thus, more systemic and in-depth investigations are required to decipher the genotype-phenotype relationships for ethanol synthesis in acetogens.
Among other differentially expressed proteins, it was unexpected to see repression of 13 ribosomal proteins in the Δhp strain (1.7- to 2.3-fold, q < 0.05; Supplementary Table S3). In addition, several oxidoreductases (05235, 14800, and 18000; ∼2-fold, q < 0.03) and sigma-54-related transcriptional regulators (14025 and 00450; ∼3.2- and 1.6-fold, q < 0.05) were upregulated in the deletion strains. Also noteworthy is the ∼1.8-fold (q < 0.01) upregulation of a PocR ligand-binding domain-containing protein (16380) in both Δhp and Δnbp strains, as we recently identified a missense mutation in the PocR gene of a CO-evolved C. autoethanogenum strain obtained via adaptive laboratory evolution (Ingelman et al., 2023). Moreover, PocR has an alternative annotation in C. autoethanogenum as a histidine kinase (Humphreys et al., 2015) that regulates sporulation and metabolism in Clostridia (Xin et al., 2020).
4 Conclusion
In summary, our study shows that CRISPR/nCas9 can effectively be used for targeted gene deletion in the industrially relevant acetogen C. autoethanogenum and provides valuable insights to accelerate the rational genetic engineering of acetogen cell factories. Additionally, the characterisation of the constructed deletion strains offers valuable clues to the in vivo functionalities of the hp and nbp genes that seem to influence both growth and product patterns of C. autoethanogenum during autotrophic growth. Furthermore, our proteomics results shed light on enzymes that potentially play a key role in ethanol production. Future work will need to confirm the biochemical function of closely linked genes in the WLP and identify possible targets for improving carbon fixation efficiencies.
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 in the article/Supplementary Material.
Author contributions
Conceptualization: UJN, KR, and KV; Methodology: UJN, KR, KS, AH, MK, and KV; Formal analysis: UJN, KR, LD, PP, and KV; Investigation: UJN, KR, LD, PP, and KV; Resources: KV; Writing—Original Draft: UJN, KR, and KV; Writing—Review and Editing: UJN, KR, LD, PP, KS, AH, MK, and KV; Supervision: KV; Project Administration: KV; Funding Acquisition: KV. All authors listed have made a substantial, direct, and intellectual contribution to the work and approved it for publication.
Funding
This work was funded by the European Union’s Horizon 2020 research and innovation programme under grant agreement N810755 and the Estonian Research Council’s grant agreement PSG289.
Acknowledgments
We thank Bastian Molitor and Peng-Fei Xia for providing us the pFX01 plasmid.
Conflict of interest
LanzaTech has interest in commercial gas fermentation with C. autoethanogenum. AH and MK are employees of LanzaTech.
The remaining 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.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fbioe.2023.1167892/full#supplementary-material
References
Adamberg, K., Valgepea, K., and Vilu, R. (2015). Advanced continuous cultivation methods for systems microbiology. Microbiol 161, 1707–1719. doi:10.1099/mic.0.000146
Altschul, S. F., Gish, W., Miller, W., Myers, E. W., and Lipman, D. J. (1990). Basic local alignment search tool. J. Mol. Biol. 215, 403–410. doi:10.1016/S0022-2836(05)80360-2
Benjamini, Y., and Hochberg, Y. (1995). Controlling the false discovery rate: A practical and powerful approach to multiple testing. J. R. Stat. Soc. B 57, 289–300. doi:10.1111/j.2517-6161.1995.tb02031.x
Bertsch, J., and Müller, V. (2015). Bioenergetic constraints for conversion of syngas to biofuels in acetogenic bacteria. Biotechnol. Biofuels 8, 210. doi:10.1186/s13068-015-0393-x
Bourgade, B., Minton, N. P., and Islam, M. A. (2021). Genetic and metabolic engineering challenges of C1-gas fermenting acetogenic chassis organisms. FEMS Microbiol. Rev. 45, fuab008. fuab008. doi:10.1093/femsre/fuab008
Brown, S. D., Nagaraju, S., Utturkar, S., De Tissera, S., Segovia, S., Mitchell, W., et al. (2014). Comparison of single-molecule sequencing and hybrid approaches for finishing the genome of Clostridium autoethanogenum and analysis of CRISPR systems in industrial relevant Clostridia. Biotechnol. Biofuels 7, 40. doi:10.1186/1754-6834-7-40
Bryksin, A., and Matsumura, I. (2013). “Overlap extension PCR cloning,” in Methods in molecular Biology. Editors K. Polizzi, and C. Kontoravdi (Totowa, NJ: Humana Press), 31–42. doi:10.1007/978-1-62703-625-2_4/COVER
Can, M., Armstrong, F. A., and Ragsdale, S. W. (2014). Structure, function, and mechanism of the nickel metalloenzymes. Chem. Rev. 8, 4149. doi:10.1021/cr400461p
Cox, J., Hein, M. Y., Luber, C. A., Paron, I., Nagaraj, N., and Mann, M. (2014). Accurate proteome-wide label-free quantification by delayed normalization and maximal peptide ratio extraction, termed MaxLFQ. Mol. Cell. Proteomics 13, 2513–2526. doi:10.1074/MCP.M113.031591
Demichev, V., Messner, C. B., Vernardis, S. I., Lilley, K. S., and Ralser, M. (2019). DIA-NN: Neural networks and interference correction enable deep proteome coverage in high throughput. Nat. Methods 171, 41–44. doi:10.1038/s41592-019-0638-x
Drake, H. L., Küsel, K., and Matthies, C. (2006). “Acetogenic Prokaryotes,” in Prokaryotes (Berlin, Germany: Springer), 354–420.
Fackler, N., Heijstra, B. D., Rasor, B. J., Brown, H., Martin, J., Ni, Z., et al. (2021). Stepping on the gas to a circular economy: Accelerating development of carbon-negative chemical production from gas fermentation. Annu. Rev. Chem. Biomol. Eng. 12, 439–470. doi:10.1146/annurev-chembioeng-120120-021122
Gillet, L. C., Navarro, P., Tate, S., Röst, H., Selevsek, N., Reiter, L., et al. (2012). Targeted data extraction of the MS/MS spectra generated by data-independent acquisition: A new concept for consistent and accurate proteome analysis. Mol. Cell. Proteomics 11, O111.016717. doi:10.1074/mcp.O111.016717
Hughes, C. S., Moggridge, S., Müller, T., Sorensen, P. H., Morin, G. B., and Krijgsveld, J. (2018). Single-pot, solid-phase-enhanced sample preparation for proteomics experiments. Nat. Protoc. 141, 68–85. doi:10.1038/s41596-018-0082-x
Humphreys, C. M., McLean, S., Schatschneider, S., Millat, T., Henstra, A. M., Annan, F. J., et al. (2015). Whole genome sequence and manual annotation of Clostridium autoethanogenum, an industrially relevant bacterium. BMC Genomics 16, 1085. doi:10.1186/s12864-015-2287-5
Ingelman, H., Heffernan, J. K., Harris, A., Brown, S. D., Yar Saqib, A., Pinheiro, M. J., et al. (2023). Autotrophic adaptive laboratory evolution of the acetogen Clostridium autoethanogenum delivers the gas-fermenting strain LAbrini with superior growth, products, and robustness. bioRxiv. doi:10.1101/2023.01.28.526018
IPCC (2021). Climate change 2021: The physical science basis. Contribution of working group I to the sixth assessment Report of the intergovernmental panel on climate change. Available at: https://report.ipcc.ch/ar6/wg1/IPCC_AR6_WGI_FullReport.pdf.
Jin, S., Bae, J., Song, Y., Pearcy, N., Shin, J., Kang, S., et al. (2020). Synthetic Biology on acetogenic bacteria for highly efficient conversion of C1 gases to biochemicals. Int. J. Mol. Sci. 21, 7639. doi:10.3390/ijms21207639
Köpke, M., and Simpson, S. D. (2020). Pollution to products: Recycling of ‘above ground’ carbon by gas fermentation. Curr. Opin. Biotechnol. 65, 180–189. doi:10.1016/j.copbio.2020.02.017
Leang, C., Ueki, T., Nevin, K. P., and Lovley, D. R. (2013). A genetic system for Clostridium ljungdahlii: A chassis for autotrophic production of biocommodities and a model homoacetogen. Appl. Environ. Microbiol. 79, 1102–1109. doi:10.1128/AEM.02891-12
Lemaire, O. N., and Wagner, T. (2021). Gas channel rerouting in a primordial enzyme: Structural insights of the carbon-monoxide dehydrogenase/acetyl-CoA synthase complex from the acetogen Clostridium autoethanogenum. Biochim. Biophys. Acta - Bioenerg. 1862, 148330. doi:10.1016/j.bbabio.2020.148330
Li, Q., Chen, J., Minton, N. P., Zhang, Y., Wen, Z., Liu, J., et al. (2016). CRISPR-based genome editing and expression control systems in Clostridium acetobutylicum and Clostridium beijerinckii. Biotechnol. J. 11, 961–972. doi:10.1002/BIOT.201600053
Liew, F. E., Nogle, R., Abdalla, T., Rasor, B. J., Canter, C., Jensen, R. O., et al. (2022). Carbon-negative production of acetone and isopropanol by gas fermentation at industrial pilot scale. Nat. Biotechnol. 40, 335–344. doi:10.1038/s41587-021-01195-w
Liew, F., Henstra, A. M., Köpke, M., Winzer, K., Simpson, S. D., and Minton, N. P. (2017). Metabolic engineering of Clostridium autoethanogenum for selective alcohol production. Metab. Eng. 40, 104–114. doi:10.1016/j.ymben.2017.01.007
Liew, F., Henstra, A. M., Winzer, K., Köpke, M., Simpson, S. D., and Minton, N. P. (2016). Insights into CO2 fixation pathway of Clostridium autoethanogenum by targeted mutagenesis. MBio 7, 004277–e516. doi:10.1128/mBio.00427-16
Ljungdahl, L. G. (2009). A life with acetogens, thermophiles, and cellulolytic anaerobes. Annu. Rev. Microbiol. 63, 1–25. doi:10.1146/annurev.micro.091208.073617
Mann, M., Effert, D., Kottenhahn, P., Hüser, A., Philipps, G., Jennewein, S., et al. (2022). Impact of different trace elements on metabolic routes during heterotrophic growth of C. ljungdahlii investigated through online measurement of the carbon dioxide transfer rate. Biotechnol. Prog. 38, e3263. doi:10.1002/BTPR.3263
Park, J., Bae, S., and Kim, J. S. (2015). Cas-designer: A web-based tool for choice of CRISPR-cas9 target sites. Bioinformatics 31, 4014–4016. doi:10.1093/BIOINFORMATICS/BTV537
Pavan, M., Reinmets, K., Garg, S., Mueller, A. P., Marcellin, E., Köpke, M., et al. (2022). Advances in systems metabolic engineering of autotrophic carbon oxide-fixing biocatalysts towards a circular economy. Metab. Eng. 71, 117–141. doi:10.1016/j.ymben.2022.01.015
Peebo, K., Valgepea, K., Nahku, R., Riis, G., Õun, M., Adamberg, K., et al. (2014). Coordinated activation of PTA-ACS and TCA cycles strongly reduces overflow metabolism of acetate in Escherichia coli. Appl. Microbiol. Biotechnol. 98, 5131–5143. doi:10.1007/s00253-014-5613-y
Perez-Riverol, Y., Bai, J., Bandla, C., García-Seisdedos, D., Hewapathirana, S., Kamatchinathan, S., et al. (2022). The PRIDE database resources in 2022: A hub for mass spectrometry-based proteomics evidences. Nucleic Acids Res. 50, D543–D552. doi:10.1093/NAR/GKAB1038
Ragsdale, S. W. (2008). Enzymology of the Wood-Ljungdahl pathway of acetogenesis. Ann. N. Y. Acad. Sci. 1125, 129–136. doi:10.1196/annals.1419.015
Ragsdale, S. W., and Pierce, E. (2008). Acetogenesis and the wood–ljungdahl pathway of CO2 fixation. Biochim. Biophys. Acta 1784, 1873–1898. doi:10.1016/j.bbapap.2008.08.012
Redl, S., Sukumara, S., Ploeger, T., Wu, L., Jensen, T. Ø., Nielsen, A. T., et al. (2017). Thermodynamics and economic feasibility of acetone production from syngas using the thermophilic production host Moorella thermoacetica. Biotechnol. Biofuels 10, 150. doi:10.1186/s13068-017-0827-8
Russell, M. J., and Martin, W. (2004). The rocky roots of the acetyl-CoA pathway. Trends biochem. Sci. 29, 358–363. doi:10.1016/j.tibs.2004.05.007
Tyanova, S., Temu, T., Sinitcyn, P., Carlson, A., Hein, M. Y., Geiger, T., et al. (2016). The Perseus computational platform for comprehensive analysis of (prote)omics data. Nat. Methods 13, 731–740. doi:10.1038/nmeth.3901
Valgepea, K., de Souza Pinto Lemgruber, R., Abdalla, T., Binos, S., Takemori, N., Takemori, A., et al. (2018). H2 drives metabolic rearrangements in gas-fermenting Clostridium autoethanogenum. Biotechnol. Biofuels 11, 55. doi:10.1186/s13068-018-1052-9
Valgepea, K., de Souza Pinto Lemgruber, R., Meaghan, K., Palfreyman, R. W., Abdalla, T., Heijstra, B. D., et al. (2017a). Maintenance of ATP homeostasis triggers metabolic shifts in gas-fermenting acetogens. Cell Syst. 4, 505–515.e5. doi:10.1016/j.cels.2017.04.008
Valgepea, K., Loi, K. Q., Behrendorff, J. B., Lemgruber, R. de S. P., Plan, M., Hodson, M. P., et al. (2017b). Arginine deiminase pathway provides ATP and boosts growth of the gas-fermenting acetogen Clostridium autoethanogenum. Metab. Eng. 41, 202–211. doi:10.1016/j.ymben.2017.04.007
Valgepea, K., Talbo, G., Takemori, N., Takemori, A., Ludwig, C., Mahamkali, V., et al. (2022). Absolute proteome quantification in the gas-fermenting acetogen Clostridium autoethanogenum. mSystems 7, e0002622–22. doi:10.1128/msystems.00026-22
Woods, C., Humphreys, C. M., Rodrigues, R. M., Ingle, P., Rowe, P., Henstra, A. M., et al. (2019). A novel conjugal donor strain for improved DNA transfer into Clostridium spp. Anaerobe 59, 184–191. doi:10.1016/j.anaerobe.2019.06.020
Xia, P.-F., Casini, I., Schulz, S., Klask, C.-M., Angenent, L. T., and Molitor, B. (2020). Reprogramming acetogenic bacteria with CRISPR-targeted base editing via deamination. ACS Synth. Biol. 9, 2162–2171. doi:10.1021/acssynbio.0c00226
Keywords: gas fermentation, acetogen, CRISPR/Cas9, genetic engineering, metabolomics, proteomics, chemostat, syngas
Citation: Nwaokorie UJ, Reinmets K, de Lima LA, Pawar PR, Shaikh KM, Harris A, Köpke M and Valgepea K (2023) Deletion of genes linked to the C1-fixing gene cluster affects growth, by-products, and proteome of Clostridium autoethanogenum. Front. Bioeng. Biotechnol. 11:1167892. doi: 10.3389/fbioe.2023.1167892
Received: 16 February 2023; Accepted: 03 May 2023;
Published: 15 May 2023.
Edited by:
Isabelle Meynial-Salles, Institut National des Sciences Appliquées de Toulouse (INSA), FranceReviewed by:
Byung-Kwan Cho, Korea Advanced Institute of Science and Technology (KAIST), Republic of KoreaArmin Ehrenreich, Technical University of Munich, Germany
Copyright © 2023 Nwaokorie, Reinmets, de Lima, Pawar, Shaikh, Harris, Köpke and Valgepea. 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: Kaspar Valgepea, kaspar.valgepea@ut.ee
†These authors have contributed equally to this work