Evolution of a Biomass-Fermenting Bacterium To Resist Lignin Phenolics

We measured the growth of C. phytofermentans in various concentrations of 12 lignocellulosic inhibitors (see Table S1 in the supplemental material) to gain a general understanding of the relative effects of aliphatic acids, furans, and phenolics (Fig. 1; see also Fig. S1 in the supplemental material). Both aliphatic acids reduce growth; acetate (Fig. 1A) was less toxic than formate (Fig. 1B) on a mass per volume basis, but both acids had similar effects in terms of molarity (Fig. S1). At low furan concentrations, we observed normal growth rates after an extended lag phase (Fig. 1C and D), similar to other bacteria that reduce and detoxify furans (15, 16). Growth lags are proposed to be due to alcohol dehydrogenase (ADH) reducing furan, causing NADH depletion and acetaldehyde accumulation (17). Additionally, the ADH protein Cphy1179 shares 28% amino acid identity with a furfural-reducing, Zn-dependent ADH (18). However, if C. phytofermentans detoxifies furans, this mechanism is abruptly overwhelmed at concentrations above 2 g liter⁻¹ 5-HMF and 1 g liter⁻¹ furfural.

We examined the toxicities of two types of phenolic acids: hydrocinnamic acids (p-coumarate and ferulate) and hydroxybenzoic acids (vanillate and 4-hydroxybenzoic acid). We found that hydrocinnamic acids (Fig. 1E and F) are more toxic than hydroxybenzoic acids (Fig. 1G and H), which supports that the propionic group on the benzene ring in hydrocinnamic acids enhances toxicity, likely by affecting how the molecules partition into the membrane. Moreover, we found that phenolic acids are typically less toxic than the corresponding aldehydes (Fig. 1I to K) and catechol (Fig. 1L). For example, vanillate (Fig. 1G) is much less toxic than vanillin (Fig. 1I), and 4-hydroxybenzoic acid (Fig. 1H) is similarly less toxic than benzaldehyde (Fig. 1J). The enhanced toxicity of aldehydes is likely due to their reactivity, resulting in formation of adducts with nucleophilic sites on DNA, proteins, and other macromolecules (19).

We quantified genome-wide mRNA expression changes at two time (t) points (t = 0.5 h and t = 4 h) following supplementation of mid-log cultures with 2 g liter⁻¹ ferulate, which reduced growth (Fig. 2A) similar to that of the initial growth screen (Fig. 1F). Three to five million read pairs were aligned to the genome for each culture (see Table S3A in the supplemental material) to calculate gene expression levels (Table S3B). The number of differentially expressed genes (Table S3C to E) increased from 0 genes before ferulate addition (Fig. 2B) to 78 genes after 30 min (Fig. 2C) and then declined to 47 genes after 4 h (Fig. 2D). The most abundant functional categories of differentially expressed genes at a t of 0.5 h relate to the repression of energy production, coenzyme metabolism, and lipids (Fig. 2E). The coenzyme-associated genes enable siroheme biosynthesis, which is repressed in clostridia in response to redox stress (20). Lipid genes include the fab gene cluster (cphy0516-cphy0523) for fatty acid biosynthesis, which was strongly repressed at a t of 0.5 h (Fig. 2E) and recovered by a t of 4 h. While cultures continued active growth after sampling, many of the differences between cultures with or without ferulate at a t of 4 h indicate that the treatment without ferulate had depleted nutrients in the medium, triggering expression of genes to assimilate alternative carbohydrates (Fig. 2F).

Gene expression at a t of 0.5 h shows that abrupt ferulate stress induces expression of genes encoding the efflux pump cphy1055-cphy1056, which is similar to E. coli mdlAB conferring resistance to organic solvents (21). Many of the genes upregulated at a t of 0.5 h are colocated in two genomic regions. The first region encodes tad (tight adherence) cphy0029-cphy0040 genes for Flp-type type IV pilus assembly. Type IV pili are widespread in clostridia (22) for adhesion to solid substrates to form protective biofilms (23), reflecting how ferulate represses motility genes in Clostridium beijerinckii (24). The other cluster cphy1838-cphy1845 includes genes for the flavin mononucleotide (FMN)-binding protein WrbA (25) and two NADPH:FMN reductases. NADPH:FMN reductase inactivation confers ferulate resistance in C. beijerinckii by an unknown mechanism (26). While this appears to be in opposition to our data showing that NADPH:FMN reductases are upregulated by ferulate, both results support the importance of FMN-mediated oxidoreduction in ferulate resistance. This island also includes genes encoding an acetyltransferase and Cphy1845 that shares 41% amino acid identity and metal coordination with E. coli YhhW, which cleaves the plant phenolic quercetin (27). C. phytofermentans may thus upregulate genes to transform or detoxify plant phenolics, similar to some ruminal clostridia (28).

We selected C. phytofermentans strains with increased ferulate resistance by cultivation in a GM3 automat, a dual-chamber, continuous-culture device that automates delivery of fresh medium and transfers the evolving cell suspension between twin growth chambers to prevent biofilm formation. During acclimation to increased ferulate in medium-swap mode, cell densities oscillated for 2 to 5 days because high densities triggered pulses of stressing medium (high ferulate) that reduced culture density, which in turn resulted in delivery of relaxing medium (low ferulate) that enabled recovery (Fig. 3A). As such, the ferulate-based selection in medium-swap mode is modulated by the ratio of relaxing and stressing medium. Once cell densities stabilized in the stressing medium, the growth rate at the higher ferulate concentration was improved in turbidostat mode (Fig. 3B). We initiated the growth selection with a stressing medium containing 1 g liter⁻¹ ferulate, the highest concentration at which we could establish a stable wild-type (WT) culture in the GM3. After 93 days (∼500 generations) of continuous, log-phase growth selection with incrementally higher ferulate, the culture grew with the same 3.75-h generation time in the 3 g liter⁻¹ ferulate medium as that of the WT in the absence of ferulate (Fig. 3C). Clones isolated along the growth selection are progressively more ferulate resistant in batch culture (Fig. 3D to G); while no growth was observed above 2 g liter⁻¹ ferulate in the WT strain (Fig. 3D), CFY3 clones grow robustly at the ferulate solubility limit (6 g liter⁻¹) (Fig. 3G). We assessed the ferulate resistance of 2 clones from each of the CFY1 (CFY1A and CFY1B) and CFY2 (CFY2C and CFY2D) time points and 4 clones from the CFY3 time point (CFY3E to CFY3H). The duplicate CFY1 and CFY2 clones showed similar ferulate resistance, but CFY3H is much less ferulate resistant than the 3 other clones (see Fig. S2 in the supplemental material), showing that cells in the GM3 culture are heterogeneous with respect to ferulate resistance.

We examined whether selection for ferulate resistance in glucose medium resulted in physiological changes impacting cellulose fermentation and resistance to other inhibitors. CFY3 strains degrade cellulose similar to that in the WT (Fig. 4A) and show accelerated cellulose degradation in medium supplemented with ferulate (Fig. 4B), supporting that the evolved strains are potentially improved candidates for fermentation of lignocellulose. Moreover, the evolved resistance mechanisms extend to other biomass inhibitors, as CFY3 strains are also more resistant to vanillate and acetate (see Fig. S3 in the supplemental material), albeit with considerable variability between strains. We also used mass spectrometry to investigate if ferulate was consumed or transformed in WT and CFY3 cultures, revealing that the ferulate concentration was unaltered with no products corresponding to reduced, demethoxylated, or decarboxylated ferulate (see Fig. S4 in the supplemental material). Thus, even though C. phytofermentans upregulates potential phenol-degrading enzymes in response to ferulate, the cell adapted to ferulate by reinforcing the cell or excluding this molecule rather than detoxifying it.

As the toxicity of aromatic molecules is often associated with disruption of the cell membrane, we profiled fatty acids (FAs) to determine if ferulate resistance is associated with altered membrane phospholipids (see Table S4 in the supplemental material). We found that when WT was exposed to ferulate, the plasmalogen (vinyl ether phospholipid) content in the membranes increased 18-fold. Moreover, CFY strains retained elevated plasmalogens even in the absence of ferulate (Fig. 4C). In particular, the CFY1B plasmalogen content in the medium without ferulate was 185-fold higher than that of the WT. Related clostridia similarly increase plasmalogens in response to aliphatic alcohol stress (29, 30), likely to fine tune membrane fluidity and protect from redox-mediated damage (31). The distribution of FA chain lengths in WT cells (Fig. 4D) is similar to that of other clostridia but with fewer unsaturated FAs and more cyclopropanes (32), both of which reduce membrane fluidity to protect from solvent stress (33). While the addition of ferulate had little immediate effect on the FA chains of WT cells (see Fig. S5A and B in the supplemental material), the CFY strains showed altered FAs relative to those of the WT in the absence of ferulate (Fig. 4D). The CFY1B FA profile was the most perturbed with increased hydroxylated C₁₆ and unsaturated fatty acids, largely C_18:1, which is associated with increased ethanol tolerance in E. coli (34). CFY3F shifted to branched FA (especially C₁₅) and longer chain lengths (C₁₈, C₂₀), which increase membrane rigidity (10) to potentially combat the membrane-fluidizing effects of ferulate. C. phytofermentans fatty acids are decorated with a diversity of phospho, glyco, and amino head groups (Fig. S5C). While we did not detect changes in these head groups in the WT response to ferulate or in the CFY strains, we consider it likely that they participate in the response to solvents, similar to some other bacteria (10).

We sequenced the genomes of eight CFY1 to CFY3 clones, giving between 106- and 705-fold coverage (see Table S5A in the supplemental material) to identify DNA variants relative to the wild type (Table S5B). Seven single-nucleotide variations (SNVs) and short insertions/deletions (indels) are present in all of the CFY genomes (Fig. 5A), which likely fixed in the population during an early selective sweep. These variants caused nonsynonymous changes in 5 proteins, including a homolog of Cap5F (Cphy3503), a protein for biosynthesis of capsular polysaccharides (35) that is associated with biofilm formation (36) and stress resistance (37). Strains subsequently accrued strain-specific mutations consistent with the population exploring alternative mutational pathways to improve ferulate resistance, particularly by modifying sensor kinases that can transduce signals associated with ferulate stress, fatty acid biosynthesis, and the surface layer (S-layer) (Fig. 5A). For example, the CFY1 and CFY2 strains incurred coding variants in 3 genes putatively encoding fatty acid biosynthesis proteins: Cphy3113 for anaerobic synthesis of unbranched fatty acids (38), the fatty acid dehydratase FabZ (Cphy0520), and the reductase FabV (Cphy1286) for the final step in fatty acid elongation. The genomes of CFY3E to CFY3G (high resistance) and CFY3H (low resistance) differ by variants in Cphy3510, the most highly expressed protein in the proteome that is proposed to form the S-layer (39). The S-layer is a protein lattice that provides mechanical stabilization, sites for extracellular protein attachment, and a selective barrier for molecules (40).

Intergenic changes that arose in the CFY genomes affect the expression levels of adjacent genes. For example, a 15-bp sequence between the first two genes of the ABC glucose transporter operon (cphy2241-cphy2243) was duplicated in the CFY3 strains (see Fig. S6A and B in the supplemental material). The repeated sequence forms an inverted repeat (IR) similar to repeated extragenic palindrome (REP) sequences, a widespread mechanism in bacteria to tune gene expression by modulating the stability of different mRNA segments within an operon (41). Duplication of this putative REP increases the mRNA secondary structure of the cphy2243-cphy2242 intergenic region (Fig. S6B), supporting functions similar to those of REP that increase expression by forming stable stem-loop structures that protect mRNA from ribonucleases (42, 43). Similarly, we found that mRNA expression of the two genes downstream of the insertion was elevated (Fig. 5B), which may have increased fitness because the GM3 growth selections were done in glucose medium. The mRNA expression of genes in two colocated operons with upstream point mutations was upregulated in the CFY strains (Fig. 5C). The A-to-G transition upstream of cphy1464 created a TG dinucleotide 2 bp upstream of the Pribnow hexamer (Fig. S6A) that enhances transcription in other bacteria (44, 45) and is present in the consensus −10 promoter sequence in C. phytofermentans (46). We propose that the upregulated operons cphy1459-cphy1461 and cphy1464-cphy1465 either enable increased production of malonyl coenzyme A (malonyl-CoA) for fatty acid biosynthesis or neutralize intracellular pH in response to ferulic acid stress through production of ammonium and lactate (7) and bicarbonate buffering (Fig. S6C).

Adaptive function can be imparted by structural changes to the genome resulting from recombination and transposition of insertion sequences (IS elements). C. phytofermentans encodes 31 IS elements (see Table S6 in the supplemental material), including 12 ISL3 comprised of 2 isoforms—8 ISL3-1 elements and 4 ISL3-2 elements. IS elements inactivate genes through their transposition and act as the substrates for homologous recombination. In addition, the IRs of all 12 ISL3 contain a 5′-TTGACA-3′ sequence matching an outward-facing, consensus −35 box from this organism (46) (Table S6), suggesting that ISL3 could activate expression of adjacent genes (47). An ISL3-1 was precisely deleted in all CFY genomes as evidenced by reduced read coverage (Fig. 6A) as well as BioNano optical mapping and Sanger sequencing (see Fig. S7A and B in the supplemental material), showing that ISL3 are active in C. phytofermentans. Further, the CFY2C and CFY2D strains share a 333-kb duplication from cphy2178 to cphy2461 (276 genes) (Fig. 6B), which we showed by PCR exists as a tandem duplication joined by a novel ISL3-2 insertion (Fig. 6C; Fig. S7C). We did not observe any extrachromosomal DNA by pulsed-field gel in the CFY2C strain, supporting that the cphy2178-ISL3-cphy2461 fragment did not excise as a circular molecule. Further, BioNano sequencing of DNA molecules greater than 400 kb spanning the junctions of the duplicated region localizes the rearrangement as a genomic, tandem duplication (Fig. 6C). We quantified the relative abundance of cells bearing this duplication by quantitative PCR of the cphy2178-ISL3-cphy2461 fragment (Fig. 6D). The duplication arose between day 13 and 20 and overtook the population to comprise 68% of cells by day 40; supporting it was the subject of positive selection. Subsequently, this variant declined in the population, representing 1% of cells at day 63, as it was gradually replaced by mutants with higher fitness; it was not present in any of the CFY3 genomes.

C. phytofermentans ISDg (ATCC 700394) was cultured anaerobically in GS2 medium (50). The growth of batch cultures containing inhibitors (see Table S1 in the supplemental material) was measured in 100-well microtiter plates (Bioscreen 9502550) containing 400 μl of GS2 medium with 3 g liter⁻¹ glucose supplemented with a given inhibitor neutralized to pH 7. Wells were inoculated with 1:10 volume cells grown to log phase in a medium without inhibitor. The plates were sealed in the anaerobic chamber (2% H₂, 98% N₂) by press-fitting adhesive sheets (Qiagen 1018104) (51) and incubated at 37°C in a Thermo Scientific Bioscreen C. The cell densities (optical density at 600 nm [OD₆₀₀]) were measured every 15 min with 30 s of shaking before each reading. Cellulose cultures were inoculated into GS2 containing 10 g liter⁻¹ cellulose (0.5- by 5-cm strips of Whatman filter paper 1001-090, >98% cellulose content). Cellulose degradation was measured as the dry mass of cellulose remaining in culture by collecting the remaining cellulose on 11-μm-pore-size filters by vacuum filtration and drying it overnight at 65°C (52).

Ferulate-resistant C. phytofermentans clones were selected using a GM3 automat (14), a dual-chamber continuous-culture device that maintained anaerobic conditions by flushing cultures with 100% N₂ gas. A 50-ml culture was maintained at 30°C with optical density readings every 30 s and was transferred between growth chambers every 12 h to clean the empty chamber with 5 N sodium hydroxide. Cells were acclimated to increased ferulate using medium-swap mode, a modified chemostat (6-h generation time) with dilutions every 30 min of stressing medium (high ferulate) if the cell density exceeded the density threshold (measured as an OD₈₈₀ of 30, which is equivalent to an OD₆₀₀ of 0.4) and relaxing medium (low ferulate) otherwise. Once cell densities stabilized at a constant cell density in the stressing medium for 24 h, the GM3 was run as a turbidostat using the stressing medium until the culture reattained a 3.75-h generation time, similar to that of the WT strain in medium without ferulate. In turbidostat mode, 20% of the culture volume was replaced with fresh medium each time the cell density reached an OD₈₈₀ of 30. Initially, the stressing medium contained 1 g liter⁻¹ ferulate, the highest concentration at which a WT culture can be established in the GM3, and the relaxing medium lacked ferulate. The medium-swap/turbidostat approach was iterated by incrementing the stress medium by 0.5 g liter⁻¹ ferulate and replacing the relaxing medium with the previous stressing medium. Samples from the GM3 culture were plated to isolate colonies, called CFY1, CFY2, and CFY3 clones, at the end of the turbidostat selections in 1, 2, and 3 g liter⁻¹ ferulate, respectively.

Log-phase cultures (OD₆₀₀ of 0.8) of WT C. phytofermentans ISDg were diluted with 1 volume medium either lacking ferulate or containing 4 g liter⁻¹ ferulate (2 g liter⁻¹ final concentration). Samples for RNA and cell densities were taken from duplicate cultures for each treatment immediately before ferulate addition and 0.5 and 4 h afterward. Total RNA was extracted using TRI reagent (Sigma 93289), and 20 μg RNA was treated with 4 U Turbo DNase (Ambion AM2238) for 30 min at 37°C. RNA was purified by Zymoclean (Zymo Research R1015) to capture RNAs of >200 bp. Five micrograms of total RNA was depleted of rRNA by Ribo-Zero (Illumina MRZMB126), yielding 200 to 400 ng RNA, and purified by Zymo Concentrator-5 (total capture) into 10 μl of water. cDNA libraries were prepared from 100 ng RNA using the TruSeq stranded mRNA kit (Illumina 15031047) and sequenced on an Illumina HiSeq 2000 sequencer with paired-end 150-bp reads. Reads were aligned to the C. phytofermentans ISDg genome (NCBI accession number NC_010001.1 ) using Bowtie 2 (53). Gene expression was calculated as reads per kilobase of gene per million (RPKM) using the easyRNASeq bioconductor package (54). Differential expression was defined as a greater than 4-fold change in expression and a DESeq (55) P value of <0.01 after Bonferroni correction for multiple testing of the 3,902 genes in the genome.

Genomes were sequenced for clones isolated from the GM3 samples: CFY1 (2 clones, CFY1A and CFY1B), CFY2 (2 clones, CFY2C and CFY2D), and CFY3 (4 clones, CFY3E to CFY3H). Genomic DNA (15 to 20 μg) was extracted from 3-ml cultures using the Sigma GenElute bacterial genomic DNA kit (NA2110). DNA (100 to 250 ng) was fragmented by Covaris E220 (Covaris, Inc., Woburn, MA, USA) to a 600-bp mean fragment size. The DNA was end-repaired, 3′ A-tailed, and ligated to Illumina compatible adapters using the NEBNext DNA sample prep master mix set 1 (New England BioLabs E6040). Ligation products were purified with 1 volume solid-phase reversible immobilization (SPRI) beads (Beckman Coulter A63880) and amplified by 12 cycles of PCR using the Kapa Hifi HotStart NGS library amplification kit (Kapa Biosystems KK2611) with P5/P7 primers. PCR products were purified (0.8 volume SPRI beads) and run on a 2% agarose gel, and DNA (700 to 800 bp) was excised and purified using the NucleoSpin extract II DNA purification kit (Macherey-Nagel 740609). cDNA libraries were sequenced using 300-bp paired-end reads on an Illumina MiSeq instrument. Reads were quality filtered by Picard (https://github.com/broadinstitute/picard ) and aligned to the C. phytofermentans ISDg reference genome (NCBI accession number NC_010001.1 ) using Bowtie 2 (53). Sequence variants (single nucleotide polymorphisms [SNPs], indels) in the CFY strains relative to the reference genome were identified using GATK (56) as described previously (57). Structural variations were detected using the breseq split-read analysis tool (58). Insertion sequences (IS) were identified using ISfinder (59).

High-molecular-weight DNA of strain CFY2C was extracted in agar plugs, which were solubilized with 0.4 U of GELase (Epicentre G09200) and dialyzed for 45 min. DNA was treated using IrysPrep reagent kit (BioNano Genomics) to prepare nicked, labeled, repaired, and stained (NLRS) DNA. Briefly, 300 ng of DNA was nicked with 10 U Nt.BspQI (NEB R0644S) for 2 h at 37°C. Nicked DNA was incubated for 1 h at 72°C with fluorescently labeled Alexa 546-dUTP and Taq polymerase (NEB M0273). Nicks were ligated using Taq ligase (NEB M0208) with deoxynucleoside triphosphates (dNTPs). DNA was counterstained with YOYO-1 (Life Technologies). NLRS DNA was loaded into IrysChips (BioNano Genomics), and data were collected on the Irys instrument (BioNano Genomics) until reaching ≥1,000-fold coverage of molecules ≥100 kb.

CFY2C DNA molecules were filtered using BioNano IrysView software (version 2.5.1) retaining molecules ≥100 kb with at least 6 label sites, yielding 32,359 molecules with an N₅₀ of 172 kb. The NCBI assembly (NCBI accession number NC_010001 ) was in silico digested with BspQI (5′-GCTCTTC-3′) and used to align and assemble CFY2C molecules using the BioNano assembly pipeline (Pipeline version 4618; RefAligner and Assembler version 4704) with the parameters used for small genomes. Molecules ≥400 kb were aligned against the NCBI assembly with a tandem duplication of bp 2689393 to 3023191 joined by an ISL3-2 element in order to identify molecules spanning the duplicated zone.

We measured mRNA expression by quantitative reverse transcription-PCR (qRT-PCR) as described previously (60). Briefly, RNA was extracted from log-phase WT and CFY3E cultures as for RNA-seq. RNA was reverse transcribed (Applied Biosystems 4368814), and mRNA expression was quantified by quantitative PCR (qPCR) (KAPA KK4621) with primers in Table S2 in the supplemental material. Expression values are means of triplicate measurements of duplicate cultures calculated by the threshold cycle method as 2^−ΔCT (61), normalized to 16S rRNA levels and multiplied by a scaling factor of 10⁶. To calculate the relative abundance of DNA variants in the GM3 cultures, genomic DNA was extracted as for genome sequencing from samples taken directly from the GM3 at 8 time points. The abundance of a DNA variant was measured by qPCR (KAPA KK4621) relative to 16S (primers in Table S2). The abundance of the DNA variant in the mixed population was calculated as 2^−ΔΔCT relative to a CFY clone that bears the variant in 100% of cells.

Ferulate concentrations were compared in WT and CFY3E cultures after 5 days of growth in GS2 medium containing 6 g liter⁻¹ ferulate by liquid chromatography-electrospray ionization-mass spectrometry (LC-ESI-MS) and tandem mass spectrometry (MS/MS) using a Dionex TCC-3000RS chromatographic system (Thermo Fisher Scientific) coupled to an Orbitrap Elite mass spectrometer (Thermo Electron Corporation) equipped with a heated electrospray ionization (HESI) source. High-pressure liquid chromatography (HPLC) separation was performed on a 5-μm, 4.6- by 150-mm SeQuant ZIC-pHILIC column (Merck) at 40°C with a flow rate of 0.5 ml min⁻¹ and a mobile phase of 10 mM (NH₄)₂CO₃ pH 9.9 (phase A) and acetonitrile (phase B). Elution was conducted using the following gradient conditions: 2 min at 80% phase B, 20-min linear gradient from 80% to 40% of phase B, 8 min at 40% phase B, 5-min increase to 80% phase B, and 15 min of 80% phase B. The mass spectrometer was operated in ESI negative ion mode using a −4.5 kV ion spray, a 275°C capillary temperature, and a mass resolution of 60,000. Sheath gas, auxiliary gas, and sweep gas flow rates were set to 60, 10, and 2 arbitrary units, respectively. Mass spectra were analyzed using Xcalibur version 2.2 (Thermo Fisher Scientific).

Cellular fatty acids were analyzed in WT and CFY3E cultures with and without ferulate and CFY1B, CFY2C, CFY2D, and CFY3F cultures without ferulate. Late log-phase cells were collected by centrifugation from cultures grown in medium either with 2 g liter⁻¹ ferulate or lacking ferulate. Fatty acid methyl esters (FAME) were obtained from 100 mg of cells by saponification, methylation, and extraction (62) and were identified using the DSMZ identification service (Braunschweig, Germany). Briefly, FAME mixtures were separated using the Sherlock microbial identification system (MIS) (Microbial ID, USA): an Agilent model 6890N gas chromatograph with a 5% phenyl-methyl silicone capillary column (0.2 mm by 25 m), a flame ionization detector, and an automatic sampler (Agilent model 7683A). Peaks were integrated and fatty acid names and percentages calculated using Sherlock MIS standard software (Microbial ID, USA). Plasmalogens were quantified as dimethyl acetyl fatty acids. Polar lipids were extracted from 100 mg cells using a chloroform–methanol–0.3% aqueous NaCl mixture (1:2:0.8, vol/vol/vol) by stirring overnight. Cells were centrifuged, and the polar lipids were recovered in the chloroform phase by adjusting the chloroform–methanol–0.3% aqueous NaCl mixture to 1:1:0.9 (vol/vol/vol). Polar lipids were resolved by two-dimensional (2D) silica gel thin-layer chromatography; dimension 1 was chloroform–methanol–water (65:25:4, vol/vol/vol), and dimension 2 was chloroform–methanol–acetic acid–water (80:12:15:4, vol/vol/vol/vol). Total lipids were detected using molybdatophosphoric acid, and specific functional groups were identified using spray reagents specific for defined functional groups (63).

All data underlying the findings are fully available without restriction. Sequencing files in FASTQ format are available in the European Nucleotide Archive under study accession numbers ERP018602 (RNA-seq) and ERP018603 (whole-genome sequencing).