Survival strategies of aerobic methanotrophs under hypoxia in methanogenic lake sediments

Background Microbial methane oxidation, methanotrophy, plays a crucial role in mitigating the release of the potent greenhouse gas methane from aquatic systems. While aerobic methanotrophy is a well-established process in oxygen-rich environments, emerging evidence suggests their activity in hypoxic conditions. However, the adaptability of these methanotrophs to such environments has remained poorly understood. Here, we explored the genetic adaptability of aerobic methanotrophs to hypoxia in the methanogenic sediments of Lake Kinneret (LK). These LK methanogenic sediments, situated below the oxidic and sulfidic zones, were previously characterized by methane oxidation coupled with iron reduction via the involvement of aerobic methanotrophs. Results In order to explore the adaptation of the methanotrophs to hypoxia, we conducted two experiments using LK sediments as inoculum: (i) an aerobic "classical" methanotrophic enrichment with ambient air employing DNA stable isotope probing (DNA-SIP) and (ii) hypoxic methanotrophic enrichment with repeated spiking of 1% oxygen. Analysis of 16S rRNA gene amplicons revealed the enrichment of Methylococcales methanotrophs, being up to a third of the enriched community. Methylobacter, Methylogaea, and Methylomonas were prominent in the aerobic experiment, while hypoxic conditions enriched primarily Methylomonas. Using metagenomics sequencing of DNA extracted from these experiments, we curated five Methylococcales metagenome-assembled genomes (MAGs) and evaluated the genetic basis for their survival in hypoxic environments. A comparative analysis with an additional 62 Methylococcales genomes from various environments highlighted several core genetic adaptations to hypoxia found in most examined Methylococcales genomes, including high-affinity cytochrome oxidases, oxygen-binding proteins, fermentation-based methane oxidation, motility, and glycogen use. We also found that some Methylococcales, including LK Methylococcales, may denitrify, while metals and humic substances may also serve as electron acceptors alternative to oxygen. Outer membrane multi-heme cytochromes and riboflavin were identified as potential mediators for the utilization of metals and humic material. These diverse mechanisms suggest the ability of methanotrophs to thrive in ecological niches previously thought inhospitable for their growth. Conclusions Our study sheds light on the ability of enriched Methylococcales methanotrophs from methanogenic LK sediments to survive under hypoxia. Genomic analysis revealed a spectrum of genetic capabilities, potentially enabling these methanotrophs to function. The identified mechanisms, such as those enabling the use of alternative electron acceptors, expand our understanding of methanotroph resilience in diverse ecological settings. These findings contribute to the broader knowledge of microbial methane oxidation and have implications for understanding and potential contribution methanotrophs may have in mitigating methane emissions in various environmental conditions. Supplementary Information The online version contains supplementary material available at 10.1186/s40793-024-00586-1.


Background
Methane, a potent greenhouse gas, has more than threefold increased its atmospheric concentrations since the pre-industrial era.Around half of global methane emissions arise from natural inland waters production, including lake sediments, wetlands, rivers, and reservoirs [1].Methane is consumed naturally by aerobic and anaerobic methanotrophs.Aerobic methanotrophy is observed in oxic environments and performed by aerobic bacterial methanotrophs, employing oxygen for methane activation and catalyze the oxidation of methane to methanol through the action of methane monooxygenase (MMO) [2,3].Anaerobic oxidation of methane (AOM) can be linked to sulfate reduction through anaerobic archaeal methanotrophs (ANMEs) and sulfate-reducing bacteria [4].This process effectively prevents the release of up to 90% of the produced methane in marine sediments, and can be efficient also in rich sulfate freshwater systems [5].In sulfate-depleted environments such as most freshwater sediments, AOM can consume over 50% of the produced methane [6] and be coupled to other electron acceptors such as nitrate, nitrite, metal oxides, and humic substances [7][8][9][10][11].
Microorganisms participating in AOM have been identified not only as ANMEs but surprisingly also as those considered aerobic bacterial methanotrophs.These bacteria were observed to be active and involved in methane oxidation beneath the chemocline in the anoxic hypolimnion [12,13] and in certain freshwater lake sediments [14][15][16][17].Certain methanotrophs employ unique mechanisms to obtain oxygen, such as the disproportionation of various molecules.For instance, Methylomirabilis (NC10) produces oxygen through the disproportionation of nitrite, utilizing it subsequently for the oxidation of methane [18].Additionally, specific alphaproteobacterial methanotrophs utilize methanobactins to generate both oxygen and hydrogen through the disproportionation of water.[19].Some bacterial methanotrophs exhibit versatility by utilizing alternative electron acceptors.Methylococcales bacteria Methylomonas denitrificans and Methylocaldum sp. have been experimentally shown to couple denitrification with methane oxidation [20][21][22], while others, such as Methylomonas, Methylosinus and Methylococcus capsulatus demonstrate methane oxidation coupled to iron reduction [23][24][25].The adaptation mechanisms, that enable such a switch, remain poorly understood.
Here we explore the remarkable genetic adaptability of Methylococcales to hypoxia in methanogenic sediments of Lake Kinneret (LK, Sea of Galilee), where our previous studies confirmed methane oxidation coupled to iron reduction (Fe-AOM) beneath the sulfate reduction zone in the iron rich methanogenic zone [26].The mediation of this Fe-AOM process was proposed to involve archaea methanogens and bacterial methanotrophs [27][28][29].Indepth analyses, including isotopes of specific fatty acid lipids, quantification of the functional gene pmoA, and metagenomic analysis, confirmed the involvement of Methylococcales-like methanotrophs in methane oxidation [15,27,28].Aerobic methanotrophy was also shown to boost iron reduction in these sediments [25], but the potential microbial mechanisms that may allow the methanotrophs to survive under hypoxia and stimulate iron reduction have remained unclear.
We hypothesized that LK methanogenic sediments harbor Methylococcales with unique strategies for survival in oxygen-limited conditions, and this research aimed to identify the key mechnisems enabling it.In order to explore this adaptation to hypoxia we conducted two experiments using LK methanogenic sediments as inoculum: (i) an aerobic "classic" methanotrophic enrichment with ambient air employing DNA stable isotope probing (DNA-SIP) and (ii) hypoxic methanotrophic enrichment with repeated spiking of 1% oxygen.Using metagenomics sequencing of DNA extracted from these experiments, we evaluated the genetic basis for their survival by comparing Methylococcales metagenome-assembled genomes (MAGs) from LK to Methylococcales genomes from diverse lineages and environments.including the.We investigated the genomic potential associated with various survival mechanisms including oxygen acquisition strategies, metabolic versatility and use of alternative electron acceptors.

Sampling site
LK is a monomictic lake in northern Israel, stratified usually between March and December, leading to anoxic hypolimnion [30].The seasonal changes lead to variations in the geochemical porewater profiles in the sediments, primarily of methane, iron, and sulfur [26].In the mixed period, the oxic-anoxic boundary reaches the sediment-water interface, affecting chemical profiles within the sediment.Sulfate in the sediment is depleted and potential contribution methanotrophs may have in mitigating methane emissions in various environmental conditions.
at around a depth of 10-20 cm depending on the stratification period and is followed by iron reduction.Methane concentration profiles generally increase with depth up to saturation levels, however, methane concentration profiles also hint at the presence of a "deep sink" correlated with an increase in concentrations of reduced iron [15,26].We collected sediment as inoculum for both experiments from the deepest point of LK (station A) using a gravity corer, as previously described by Bar-Or et al. [27].The specific sediment used for the experiments originated from a depth of 25-40 cm below the sediment surface.

Study design overview
To achieve this study's aims, we designed and conducted two experiments intended mainly to enrich methanotrophs.Experiment 1 was a "classical" enrichment with ambient air, combined with DNA-Stable Isotope Probing (DNA-SIP).This experiment focused on identifying the active methanotrophs thriving under oxic conditions (Fig. 1A).In the second experiment, we stimulated methanotroph growth under hypoxic conditions (Fig. 1B).Following enrichments, we used amplicon sequencing of the 16S rRNA gene to study microbial diversity.Using metagenomic sequencing of representative samples, we curated Methylococcales MAGs.We investigated functions involved in adaptation to anoxia in these MAGs, alongside other Methylococcales genomes, using various bioinformatic platforms (Fig. 1C).

Experiment 1-ambient air enrichment and DNA stable isotope fingerprinting
Sediment samples collected from LK in September 2019 were used to enrich methanotrophs and explore the active methane-utilizing organisms.Approximately 5 g of homogenized sediment was added to 300 ml sealed serum vials with black rubber septa.The vials were prefilled with modified DNMS medium [31] (1:10 w/v).Prior to inoculation, the headspace of each vial was purged with nitrogen gas for a brief period (e.g., 1 h) to flush out any methane gas remaining in the sediment.Following this, the headspace was replaced with ambient air.
Fig. 1 Experimental layout.Experimental design-sediment cores were collected from the methanotrophic zone in Lake Kineret and used to initiate two enrichment experiments.A Experiment 1: Incubated with ambient air, representing a "classical" enrichment approach coupled with DNA-SIP analysis.B Experiment 2: Hypoxic enrichment by repeatedly introducing small amounts (spikes) of oxygen (1%) to maintain low-oxygen conditions.Extracted DNA was used for 16S rRNA amplicon sequencing, and representing samples were also used to assemble and bin Methylococcales MAGs.Genome-based metabolic comparative analyses were performed using different bioinformatic platforms C Subsequently, either 12 C-methane (control) or 13 C-methane (used for DNA-SIP) was introduced at a concentration of 20% (v/v) as the sole carbon source.Incubations proceeded at 25 °C in the dark on a 120 rpm rotary shaker for a total of 21 days.Triplicate incubations were performed for each treatment, with detailed information on sampling times provided in the Supplementary Information (SI, Table S1).A separate 5 g sediment sample (time-zero control) was preserved in a 15 ml Falcon tube at − 20 °C for subsequent DNA extraction.
After incubation, sediment samples were carefully collected and processed for DNA extraction.Briefly, samples were centrifuged, supernatants discarded, and the resulting pellets were stored at − 20 °C.DNA was then extracted from 0.5 g aliquots of each sample (time-zero, 12 C-methane, and 13 C-methane) using a PowerSoil kit (Qiagen) following the manufacturer's instructions.To achieve sufficient DNA concentration for DNA-SIP analysis, seven aliquots from each sample were extracted and pooled.DNA quality and quantity were assessed using a NanoDrop Spectrophotometer (Thermo Fisher Scientific).
For DNA-SIP, samples were further processed through CsCl density gradient using a Beckman Coulter ultracentrifuge and a VTi 65.2 rotor, operating at 44,100 rpm at 20 °C for 65 h, as previously outlined [32].Twelve fractions were collected from each sample, and their DNA content and their density was determined using an AR200 digital hand-held refractometer (Reichert, Buffalo, NY, USA).These DNA fractions were recovered, purified by PEG-6000 precipitation, resuspended in 30 μl TE buffer, and the DNA quantity was validated using Fluorometer (Qubit, Invitrogen).

Experiment 2-hypoxic enrichments
This experiment attempted to enrich methanotrophs in LK sediments with low oxygen conditions following the methods of Vigderovich et al. [25].Sediment samples collected in March 2021 set as the starting material for the experiment.Each vial contained 7g of methanogenic sediments from LK, 50 mg of hematite (10 mM final concentration) and anoxic filtered porewater (1:3 w/v), previously extracted from the same depth in a 60 ml serum vials sealed with black rubber septa.Hematite, an iron oxide mineral, can potentially serve as an alternative electron acceptor for methanotrophs under anoxic conditions [25].The anoxic porewater further ensured the initial absence of oxygen within the vials.To establish these anoxic conditions, each vial underwent continuous nitrogen purging for about an hour.Following this, the vials were stored in the dark at 4 °C for 5 days.After the anoxic period, the enrichment process was initiated.Vials were vigorously mixed (vortexed) and flushed again with nitrogen for 15 min to ensure the removal of any residual oxygen.Subsequently, 1.5 ml of air was introduced into each vial, resulting in a final headspace oxygen concentration of 1% (v/v).Two treatment groups were established CH 4 + O 2 treatment and N 2 + O 2 as Control.For the CH 4 + O 2 treatment, eleven vials received an additional 1.5 ml of 99.95% methane gas, resulting in a final headspace oxygen concentration of 5% (v/v) allowing methanotrophs to potentially uptake methane as an energy source.For the N 2 + O 2 Control, triplicate vials were supplemented with an additional 1.5 ml of 99.999% nitrogen gas instead of methane, resulting in a final headspace oxygen concentration of 5% (v/v), serving as a control for the effect of methane availability.The enrichment incubations proceeded at 25 °C in the dark with the vials inverted to minimize headspace gas exchange.Oxygen consumption within the vials was monitored throughout the experiment using PSt6 sensors in designated vials.Vials were sacrificed at various time intervals, detailed in the Supplementary Information (Table S2) and the experiment lasted for 40 days.These samples were promptly stored at -80 °C for subsequent DNA/RNA extraction.DNA was successfully extracted from all samples using a PowerSoil kit (Qiagen) following the manufacturer's instructions.However, RNA concentrations using Quant-ittTM RiboGreen RNA kit (ThermoFisher) were found to be below the detection limit.

DNA library preparation and sequencing
A total of thirty-seven samples underwent 16S rRNA gene sequencing, while three samples were also allocated for metagenome analysis.Comprehensive details regarding the samples utilized in the DNA-SIP and hypoxic experiments can be found in Tables S1 and S2, respectively.Sequencing of the V4 region of the 16S rRNA gene was executed using the primer pair 515f-806r [33,34], and the procedure was conducted via Illumina sequencing at Hylabs, Rehovot, Israel.
Metagenomic libraries were generated from three distinct DNA samples: one originating from the DNA-SIP experiment and two from the hypoxic experiment (refer to Tables S1 and S2 for detailed descriptions, respectively).Subsequently, the libraries underwent sequencing at Novogene, Singapore, with each sample producing 65-80 million 2 × 150 bp paired-end reads using Illumina NovaSeq.The library construction process employed the NEBNext ® Ultra ™ II DNA Library Prep Kit.[33,34]

Bioinformatics
For 16S rRNA gene amplicons, demultiplexed paired-end reads were analyzed within QIIME2 V2020.6 pipeline [35].By applying the DADA2 pipeline [36], implemented in QIIME2, reads were truncated according to their quality plots, chimeras were removed, and reads were merged and grouped into amplicon sequence variants (ASVs).Taxonomy was assigned to ASVs by Silva 138 99% classifier [37].Beta diversity was visualized by Principal Coordinate Analysis (PCoA), in which the dissimilarity between samples was estimated with a Bray-Curtis distance matrix using the distance and ordination functions in the R package phyloseq [38].
Phylogenetic trees were generated with the BV-BRC platform [53].Marker proteins that are universally conserved across the bacterial domain were extracted from genomes [54].A hundred single-copy markers that were present in all genomes (See Table S5 for gene list) analyzed in this study were used for alignment with MUS-CLE [55].The randomized accelerated maximum likelihood (RAxML version 8.2.11) tree was calculated [56].Final representation of the tree was curated using Itol version 6.9 (https:// itol.embl.de/).

Results and discussion
Ambient oxygen selects for the enrichment of aerobic methanotrophs DNA-SIP experiments utilizing labeled 13 C-methane resulted in the isolation of 13 C-enriched DNA observed in fractions with a density ranging from 1.72 to 1.70 g ml −1 , as illustrated in the labeled fraction (Fig. 2a).Analysis of the relative abundance of 16S rRNA genes revealed bacterial communities that differed between the labeled and unlabeled fractions, as evidenced by Principal Coordinates Analysis (PCoA) based on Bray-Curtis dissimilarity (Fig. 2b).Axis 1 and 2 explained 43.9% and 17.4% of the variation, respectively, indicating distinct communities in the labeled and unlabeled fractions, while 12 C-methane-fed communities clustered together.The use of the 12 C-methane-fed sample as a control accounted for differences in the guanine/cytosine (GC) content of DNA samples.
Despite Methylococcales comprising less than 1% of the initial microbial population at time zero, enrichment during the experiment resulted in Methylococcales accounting for approximately one-third of the microbial community (refer to Figure S1).In the fractions fed with labeled 13 C-methane, Methylobacter emerged as the predominant taxon with a relative read abundance ranging from 18 to 31%, alongside the enrichment of other methanotrophs (Fig. 2c).This finding aligns with previous studies that identified Methylobacter as a dominant and active species in diverse natural environments [57], including oxic [58] and anoxic freshwater lake sediments [59], anoxic water columns [60], and wetlands [61].Other enriched Methylococcales species in this experiment included Methylomonas (5-8%) and Methylogaea (up to 3%) (Fig. 2c).In contrast, within the 12 C-methane-enriched cultures, Methylobacter exhibited a relative abundance of 6-10%, while the relative abundance of Methylomonas and Methylogaea varied significantly among biological triplicates, ranging from 4% to up to 20% (refer to Figure S1).

Hypoxic conditions select for Methylomonas
We observed a notable adaptation of microbial communities to periodic spiking with 1% oxygen and methane, leading to a significant increase in oxygen consumption (0.69 ± 0.04 mg L −1 h −1 ), in contrast to controls where the oxygen consumption was markedly lower (0.02 ± 0.01 mg L −1 h −1 ) (Fig. 3a).Microbial communities in cultures deprived of methane exhibited clustering with those from non-enriched samples (time zero), indicating that the microbial community in these samples remained  relatively consistent despite exposure to oxygen.These findings suggest that methane serves as the electron donor for oxygen respiration, especially considering the limited availability of alternative electron donors in these methanogenic sediments (Fig. 3a).
Further analysis of the microbial community revealed the specific selection of microorganisms through periodic spiking with 1% oxygen (Fig. 3b), predominantly enriching for Methylomonas (Fig. 3c).
These outcomes align with previous studies, reinforcing the notion that Methylomonas may exhibit better adaptation to low oxygen conditions, such as, better uptake of oxygen and methane, compared to Methylobacter [62].

Methylococcales lineages dominate enrichment cultures
In this investigation, methanotrophs were selectively enriched from methanogenic LK sediments through two distinct enrichment experiments.The first experiment involved the use of ambient air coupled with DNA-stable isotope probing (DNA-SIP), while the second experiment included repeated injections of 1% oxygen to simulate hypoxic conditions.From both experiments, we identified five Methylococcales metagenome-assembled genomes (MAGs), detailed in Table 1 for genome statistics.The metagenomic analysis uncovered the enrichment of Methylomonas, Methylogaea, and Methylobacter (MAGs 1-5 respectively).The closest relatives to these lineages were Methylomonas sp.ZR1 [63], Methylogaea oryzae strain E10 [64] and Methylobacter tundripaludum strain OWC-G53F [63] (see Fig. 4 and Table 2 for ANI values).For the hypoxic enrichments, two MAGs, namely, Methylomonas LK_4 and LK_5, which are related to Methylomonas sp.strain FW.007 [65].A detailed genomic comparison between all LK Methylococcales MAGs highlighting functions of interest is given below and summarized in Table 2.

Metabolic reconstruction of novel Lake Kinneret methanotrophs
Genome-base metabolic reconstruction of the five MAGs affirmed the characteristic metabolic framework of Methylococcales methanotrophs (Table 2).The prediction of methane oxidation through the particulate methane monooxygenase was consistent across all LK Methylococcales lineages.With the exception of Methylogaea LK_2, which exclusively harbored the pxmABC genes encoding an alternative form of particulate methane monooxygenase potentially functioning under oxygen-limiting conditions [20], the remaining LK Methylococcales lineages were predicted to possess one copy of the pmoCAB operon, encoding the canonical particulate methane monooxygenase.The genomes of Methylomonas LK_4 and LK_5 contained both pxmABC and pmoCAB operons, potentially broadening their range of affinities (See figure on next page.)Fig. 4 Phylogenomic analysis and metabolic profiling.A phylogenomic tree along with the metabolic presence-absence profile of 67 Methylococcales genomes.A comprehensive list of the proteins is available in Table S7 and predicted OMC proteins are available in Table S8.Additional information regarding presence-absence of kye genes particulate methane monooxygenase pmoCAB operon and pxmABC operon, soluble methane monooxygenase (mmoXYBZDC), lanthanide-dependent methanol dehydrogenases (xoxF), and methanol dehydrogenase (mxaF) is available in Table S9.Notably, the five LK Methylococcales are highlighted in bold.The phylogenetic tree is built on 100 genes (refer to Table S5), and the taxa clustering percentage is based on 100 bootstrap resamples, consistently yielding values of 98 or higher (specific values not shown) Fig. 4 (See legend on previous page.) for methane and oxygen [66].In addition, none of the genomes of the LK Methylococcales lineages encoded the soluble methane oxygenase.The lanthanide-dependent xoxF [67] was predicted to catalyze methanol oxidation in all LK Methylococcales lineages except Methylogaea, whereas the calcium-dependent mxaFI [68] was predicted only in the genome of LK Methylomonas and likely within the Methylogaea LK_2 genome (only mxaI was found).The potential secretion of methanol by methanotrophs into the environment, with subsequent consumption by syntrophic partners, indicative of cross-feeding between methanotrophs and methylotrophs [69], was considered highly likely.This inference is supported by the elevated relative abundance of Methylotenera in our enrichment experiments (refer to Figures S1 and  S2 for microbial community relative abundance).

Adaptations to hypoxic conditions in Methylococcales
To assess whether the genomic adaptations to hypoxic conditions are unique to LK Methylococcales, we compared their genomes and with an additional 62 Methylococcales genomes from diverse environments (Fig. 4).These genomes represent Methylococcales from groundwater, contaminated rivers, sewage systems, lake sediments, rice fields, volcano mud, and others (Supplementary Table S3).
This comparison revealed a core set of mechanisms shared by most Methylococcales, including those from LK potentially enabling them to function effectively in hypoxic conditions.These include (i) Enhanced oxygen usage by cytochrome bd ubiquinol oxidase (found in 88% of analyzed genomes), allowing efficient oxygen respiration potentially sustaining growth at ≤ 3 nM molecular oxygen [72].(ii) Enhancing respiration under hypoxia by oxygen-binding hemerythrin, by increasing the activity of pMMO [73] (present all of the analyzed genomes).(iii) Flagella-mediated motility across gradients to optimize oxygen and methane availability (present in 95% of the analyzed genomes) [74].(iv) Alternative electron acceptors including nitric oxide reduction, as well as riboflavins (found in 70% and 89% of the genomes, respectively).These last electron acceptors, are soluble secreted electron shuttles mediating extracellular electron transfer (EET) [75][76][77], which can be reversibly oxidized and reduced, carrying electrons between cells and insoluble electron acceptors such as manganese and iron oxides over large distances [71,75,78].Manganese and iron oxides are highly abundant in methanogenic sediments [79].Methylococcales can store carbohydrates and use alternative metabolic pathways to provide energy under oxygen limitation.Genes for glycogen synthesis and degradation were found in all of the genomes, likely allowing Methylococcales to conserve resources during periods of limited nutrients [80,81].Additional metabolic alternative is the ability to produce fermentation products like succinate and acetate was found in 66% and 51% of genomes, respectively.The occurrence of alcohol dehydrogenases (present in 91% of genomes) and bidirectional hydrogenases (present in 99% of genomes) indicates potential for alcohol and hydrogen production [60,82].
While many metabolic mechanisms were shared among Methylococcales, we identified several less abundant mechanisms that were present in LK Methylococcales.include (i) production of lactate (found in 21% of the genomes), predicted in Methylobacter LK_1 and Methylogaea LK_2 MAGs; (ii) Methane-dependent denitrification serves as a link between the carbon and nitrogen cycles [20,83].While nitrous oxide reduction was found as a common trait, respiratory nitrate reductase and nitric oxide-forming nitrite reductase (with respective genes identified in 36% and 27% of the genomes) were observed as less prevalent mechanisms.However, despite being relatively uncommon among the Methylococcales, these mechanisms were found more abundant in LK Methylococcales and exhibited in Methylobacter LK_1, Methylomonas LK_3, and LK_5 MAGs; (iii) Outer membrane cytochromes (OMCs, found in 62% of the genomes) needed to reduce iron [24,59,84], were found in Methylomonas LK_3, LK_4 and LK_5.Methylococcales appear to lack electrically conductive pili (e-pili) that can support EET [84], as > 120 amino acid long PilA proteins in Methylococcales were longer than the canonical 60-90 amino acid-long e-pili [85].
Additional manual investigation of OMCs using the InterPro database confirmed that at least some of these sequences belong to the multiheme cytochrome superfamily with at least one copy of predicted as extracellular OMC (found in 55% of the genomes), and fewer OMCs were predicted as membrane-bound (found in 10% of the genomes).Others were either unrelated to OMCs or predicted to be hydroxylamine oxidoreductase or cytochromes c-552 involved in ammonia oxidation and nitrite reduction.Among LK Methylococcales, only Methylomonas LK_3 encoded a membrane-associated OMC, hinting at the possibility of EET in this organism, but this strategy to cope with oxygen limitation is not widespread in LK.

Conclusions
Our study highlights the potential prevalence and diversity of adaptive strategies utilized by methanotrophic bacteria in low-oxygen environments, specifically within LK sediments.These LK sediments, located 20 cm below the sediment water interface and 20 m below the hypolimnion, were previously shown to be involve Methylococcales in methane oxidation and stimulation of iron reduction.We hypothesized that this environment harbors Methylococcales with unique strategies for survival.Through enrichment experiments, comparative metagenomics, and genomic analyses of diverse Methylococcales lineages, we propose several potential mechanisms enabling these organisms to thrive under oxygen-limited conditions.Our findings demonstrate that most Methylococcales, including those from LK, possess a set of traits enabling their survival in hypoxic environments: spanning effective usage of trace oxygen, motility for reaching optimal oxygen concentrations, glycogen storage, alternative energy generation using fermentation, and the use of alternative electron acceptors and possible EET mediated by riboflavins.Some less widespread functions employed by LK Methylococcales include lactate productionn, methane-dependent denitrification, and EET via outer membrane cytochromes (OMCs).Further studies are necessary to validate these potential strategies.Experimental examination of these potential adaptations may unveil additional layers of complexity in the survival strategies.Additionally, investigating whether these adaptations are expressed under hypoxic conditions is essential for understanding their role in hypoxia tolerance.This knowledge may contribute to a more comprehensive understanding of global methane cycling.

Fig. 2
Fig. 2 Aerobic enrichment experiment employing DNA-SIP.A The average DNA concentration (ng/μl) for each fraction (n = 3 biological replicates) in the 13 C-methane-fed cultures (blue) and the 12 C-methane-fed cultures (colorless).Error bars indicate the standard error (standard deviation/n) of the DNA concentration, while the error bars for density are smaller than the symbol.B Principal Component Analysis (PCoA) based on Bray-Curtis dissimilarity illustrating microbial diversity in labeled DNA fractions (circles), unlabeled DNA fractions (triangles), and time zero samples that were not fractionated (squares) for both 13 C-methane-fed cultures (in blue) and 12 C-methane-fed cultures (colorless).C The average relative abundance of dominant Methylococcales in the labeled DNA fractions of the 13 C-methane-fed cultures

Fig. 3
Fig. 3 Hypoxic enrichment experiment.A In vitro monitoring of oxygen levels (%) in bottles exposed to O 2 + CH 4 (blue) and O 2 + N 2 (green).B Principal Component Analysis illustrating microbial diversity exposed to O 2 + CH 4 (circles), O 2 + N 2 (triangles), and time zero samples without enrichment (squares).These are color-coded based on treatment, with blue, green, and gray representing O 2 + CH 4 , O 2 + N 2 , and time zero, respectively.C The average relative abundance of dominant Methylococcales spp.enriched in the O 2 + CH 4 treatment

Table 1
Genome assembly statisticsa mb-million base pairs, b ANI-Average nucleotide identity

Table 2
List of presence/ absence functionsDetailed description of all predicted proteins used to construct this table are listed in TableS6