Efficient CRISPR/Cas12a-based genome editing toolbox for metabolic engineering in Methanococcus maripaludis

Methanococcus maripaludis is a fast-growing and genetically tractable methanogen. To become a useful host organism for the biotechnological conversion of CO2 and renewable hydrogen to fuels and value-added products, its product scope needs to be extended. Metabolic engineering requires reliable and efficient genetic tools, in particular for genome editing related to the primary metabolism that may affect cell growth. We have constructed a genome editing toolbox by utilizing Cas12a from Lachnospiraceae bacterium ND2006 (LbCas12a) in combination with the homology-directed repair machinery natively present in M. maripaludis. The toolbox enables gene knock-out with a positive rate typically above 89%, despite M. maripaludis being hyper-polyploid. We have replaced the flagellum operon (around 8.9kb) by a beta-glucuronidase gene to demonstrate a larger deletion, and to enable quantification of promotor strengths. The CRISPR/LbCas12a toolbox presented here is currently perhaps the most reliable and fastest method for genome editing in a methanogen.

Growth Coordination Between Butyrate-Oxidizing Syntrophs and Hydrogenotrophic Methanogens

Syntrophy is a thermodynamically required mutualistic cooperation between fatty acid-oxidizing bacteria and methanogens that plays the important role in organic decomposition and methanogenesis in anoxic environments. In this study, three experiments were conducted to evaluate the cell-to-cell interaction in a thermophilic coculture consisting of Syntrophothermus lipocalidus and Methanocella conradii and a mesophilic coculture consisting of Syntrophomonas wolfei and Methanococcus maripaludis. First, syntrophs and methanogens were inoculated at different initial cell ratios to evaluate the growth synchronization. The quantitative PCR analysis revealed that the organism with a lower relative abundance at the beginning always grew faster, and the cell ratio converged over time to relative constant values in both the thermophilic and mesophilic cocultures. Next, intermittent ultrasound and constant shaking treatments were used to evaluate the influence of physical disturbance on microbial aggregation in the mesophilic coculture. The fluorescence in situ hybridization and scanning electron microscopy revealed that the tendency of syntrophic aggregation was not affected by the physical disturbances, although the activity was slightly depressed. Syntrophomonas dominated in the initial microbial aggregates, which, however, did not grow until Methanococcus was attached and increased to a significant extent, indicating the local growth synchronization during the formation and maturation of syntrophic aggregates. Last, microfluidic experiments revealed that whether or not Syntrophomonas or Methanococcus was loaded first, the second organism preferred moving to the place where the first organism was located, suggesting the cell-to-cell attraction between Syntrophomonas and Methanococcus. Collectively, our study demonstrated the growth synchronization and cell-to-cell attraction between the butyrate-oxidizing bacteria and methanogens for optimizing the syntrophic cooperation.

CryoEM reveals the stochastic nature of individual ATP binding events in a group II chaperonin

Chaperonins are homo- or hetero-oligomeric complexes that use ATP binding and hydrolysis to facilitate protein folding. ATP hydrolysis exhibits both positive and negative cooperativity. The mechanism by which chaperonins coordinate ATP utilization in their multiple subunits remains unclear. Here we use cryoEM to study ATP binding in the homo-oligomeric archaeal chaperonin from Methanococcus maripaludis (MmCpn), consisting of two stacked rings composed of eight identical subunits each. Using a series of image classification steps, we obtained different structural snapshots of individual chaperonins undergoing the nucleotide binding process. We identified nucleotide-bound and free states of individual subunits in each chaperonin, allowing us to determine the ATP occupancy state of each MmCpn particle. We observe distinctive tertiary and quaternary structures reflecting variations in nucleotide occupancy and subunit conformations in each chaperonin complex. Detailed analysis of the nucleotide distribution in each MmCpn complex indicates that individual ATP binding events occur in a statistically random manner for MmCpn, both within and across the rings. Our findings illustrate the power of cryoEM to characterize a biochemical property of multi-subunit ligand binding cooperativity at the individual particle level.

Using genome comparisons of wild-type and resistant mutants of Methanococcus maripaludis to help understand mechanisms of resistance to methane inhibitors

Methane emissions from enteric fermentation in the ruminant digestive system generated by methanogenic archaea are a significant contributor to anthropogenic greenhouse gas emissions. Additionally, methane produced as an end-product of enteric fermentation is an energy loss from digested feed. To control the methane emissions from ruminants, extensive research in the last decades has been focused on developing viable enteric methane mitigation practices, particularly, using methanogen-specific inhibitors. We report here the utilization of two known inhibitors of methanogenic archaea, neomycin and chloroform, together with a recently identified inhibitor, echinomycin, to produce resistant mutants of Methanococcus maripaludis S2 and S0001. Whole-genome sequencing at high coverage (> 100-fold) was performed subsequently to investigate the potential targets of these inhibitors at the genomic level. Upon analysis of the whole-genome sequencing data, we identified mutations in a number of genetic loci pointing to potential mechanisms of inhibitor action and their underlying mechanisms of resistance.

The oligosaccharyltransferase AglB supports surface-associated growth and iron oxidation in Methanococcus maripaludis

Most microbial organisms grow as surface-attached communities known as biofilms. However, the mechanisms whereby methanogenic archaea grow attached to surfaces have remained understudied. Here, we show that the oligosaccharyltransferase AglB is essential for growth of Methanococcus maripaludis strain JJ on glass or metal surfaces. AglB glycosylates several cellular structures such as pili, archaella, and the cell surface layer (S-layer). We show that the S-layer of strain JJ, but not strain S2, is a glycoprotein, that only strain JJ was capable of growth on surfaces, and that deletion of aglB blocked S-layer glycosylation and abolished surface-associated growth. A strain JJ mutant lacking structural components of the type IV-like pilus did not have a growth defect under any conditions tested, while a mutant lacking the pre-flagellin peptidase (Δ flaK ) was only defective for surface growth when formate was provided as the sole electron donor. Finally, for strains that are capable of Fe 0 oxidation, we show that deletion of aglB decreases the rate of anaerobic Fe 0 oxidation, presumably due to decreased association of biomass with the Fe 0 surface. Together, these data provide an initial characterization of surface-associated growth in a member of the methanogenic archaea. Importance Methanogenic archaea are responsible for producing the majority of methane on Earth and catalyze the terminal reactions in the degradation of organic matter in anoxic environments. Methanogens often grow as biofilms associated with surfaces or partner organisms; however, the molecular details of surface-associated growth remain uncharacterized. We have found evidence that glycosylation of the cell surface layer is essential for growth of M. maripaludis on surfaces and can enhance rates of anaerobic iron corrosion. These results provide insight into the physiology of surface-associated methanogenic organisms and highlight the importance of surface association to anaerobic iron corrosion.

An alternative resource allocation strategy in the chemolithoautotrophic archaeon Methanococcus maripaludis

Most microorganisms in nature spend the majority of time in a state of slow or zero growth and slow metabolism under limited energy or nutrient flux rather than growing at maximum rates. Yet, most of our knowledge has been derived from studies on fast-growing bacteria. Here, we systematically characterized the physiology of the methanogenic archaeon Methanococcus maripaludis during slow growth. M. maripaludis was grown in continuous culture under energy (formate)-limiting conditions at different dilution rates ranging from 0.09 to 0.002 h−1, the latter corresponding to 1% of its maximum growth rate under laboratory conditions (0.23 h−1). While the specific rate of methanogenesis correlated with growth rate as expected, the fraction of cellular energy used for maintenance increased and the maintenance energy per biomass decreased at slower growth. Notably, proteome allocation between catabolic and anabolic pathways was invariant with growth rate. Unexpectedly, cells maintained their maximum methanogenesis capacity over a wide range of growth rates, except for the lowest rates tested. Cell size, cellular DNA, RNA, and protein content as well as ribosome numbers also were largely invariant with growth rate. A reduced protein synthesis rate during slow growth was achieved by a reduction in ribosome activity rather than via the number of cellular ribosomes. Our data revealed a resource allocation strategy of a methanogenic archaeon during energy limitation that is fundamentally different from commonly studied versatile chemoheterotrophic bacteria such as E. coli.

Perfect Match Genomic Landscape strategy: Refinement and customization of reference genomes

When addressing a genomic question, having a reliable and adequate reference genome is of utmost importance. This drives the necessity to refine and customize reference genomes (RGs). Our laboratory has recently developed a strategy, the Perfect Match Genomic Landscape (PMGL), to detect variation between genomes [K. Palacios-Flores et al.. Genetics 208, 1631–1641 (2018)]. The PMGL is precise and sensitive and, in contrast to most currently used algorithms, is nonstatistical in nature. Here we demonstrate the power of PMGL to refine and customize RGs. As a proof-of-concept, we refined different versions of the Saccharomyces cerevisiae RG. We applied the automatic PMGL pipeline to refine the genomes of microorganisms belonging to the three domains of life: the archaea Methanococcus maripaludis and Pyrococcus furiosus; the bacteria Escherichia coli, Staphylococcus aureus, and Bacillus subtilis; and the eukarya Schizosaccharomyces pombe, Aspergillus oryzae, and several strains of Saccharomyces paradoxus. We analyzed the reference genome of the virus SARS-CoV-2 and previously published viral genomes from patients’ samples with COVID-19. We performed a mutation-accumulation experiment in E. coli and show that the PMGL strategy can detect specific mutations generated at any desired step of the whole procedure. We propose that PMGL can be used as a final step for the refinement and customization of any haploid genome, independently of the strategies and algorithms used in its assembly.

Synergistic epistasis enhances the co-operativity of mutualistic interspecies interactions

Early evolution of mutualism is characterized by big and predictable adaptive changes, including the specialization of interacting partners, such as through deleterious mutations in genes not required for metabolic cross-feeding. We sought to investigate whether these early mutations improve cooperativity by manifesting in synergistic epistasis between genomes of the mutually interacting species. Specifically, we have characterized evolutionary trajectories of syntrophic interactions of Desulfovibrio vulgaris (Dv) with Methanococcus maripaludis (Mm) by longitudinally monitoring mutations accumulated over 1000 generations of nine independently evolved communities with analysis of the genotypic structure of one community down to the single-cell level. We discovered extensive parallelism across communities despite considerable variance in their evolutionary trajectories and the perseverance within many evolution lines of a rare lineage of Dv that retained sulfate-respiration (SR+) capability, which is not required for metabolic cross-feeding. An in-depth investigation revealed that synergistic epistasis across pairings of Dv and Mm genotypes had enhanced cooperativity within SR− and SR+ assemblages, enabling their coexistence within the same community. Thus, our findings demonstrate that cooperativity of a mutualism can improve through synergistic epistasis between genomes of the interacting species, enabling the coexistence of mutualistic assemblages of generalists and their specialized variants.