Correlation of methane production with physiological traits in Trichodesmium IMS 101 grown with methylphosphonate at different temperatures

The diazotrophic cyanobacterium Trichodesmium has been recognized as a potentially significant contributor to aerobic methane generation via several mechanisms including the utilization of methylphophonate (MPn) as a source of phosphorus. Currently, there is no information about how environmental factors regulate methane production by Trichodesmium. Here, we grew Trichodesmium IMS101 at five temperatures ranging from 16 to 31°C, and found that its methane production rates increased with rising temperatures to peak (1.028 ± 0.040 nmol CH4 μmol POC−1 day−1) at 27°C, and then declined. Its specific growth rate changed from 0.03 ± 0.01 d−1 to 0.34 ± 0.02 d−1, with the optimal growth temperature identified between 27 and 31°C. Within the tested temperature range the Q10 for the methane production rate was 4.6 ± 0.7, indicating a high sensitivity to thermal changes. In parallel, the methane production rates showed robust positive correlations with the assimilation rates of carbon, nitrogen, and phosphorus, resulting in the methane production quotients (molar ratio of carbon, nitrogen, or phosphorus assimilated to methane produced) of 227–494 for carbon, 40–128 for nitrogen, and 1.8–3.4 for phosphorus within the tested temperature range. Based on the experimental data, we estimated that the methane released from Trichodesmium can offset about 1% of its CO2 mitigation effects.


Introduction
Methane (CH 4 ) is a crucial part of the carbon cycle and a potent greenhouse gas, with a global warming potential of more than 80 times that of CO 2 over a 20-year period (IPCC, 2013).Intriguingly, methane is typically supersaturated in the ocean's top mixed layer (Lamontagne et al., 1973;Scranton and Brewer, 1977;Weber et al., 2019), a phenomenon known as "the marine methane paradox." The paradoxical nature of this phenomenon stems from the historically prevailing recognition that biogenic methane is produced exclusively by methanogenic archaea in a strictly anoxic environment, a condition incompatible with the oxic state of ocean top mixed layer (Karl and Tilbrook, 1994).
The aerobic metabolism of MPn serves as a crucial source of phosphorus and subsequently a significant source of methane (Karl et al., 2008;Repeta et al., 2016;Von Arx et al., 2023).Despite dissolved inorganic phosphate (DIP) typically being the most bioavailable form of phosphorus, its availability is commonly limited in pelagic surface waters (Sañudo-Wilhelmy et al., 2001;Dyhrman et al., 2002).It has been documented that concentrations of dissolved organic phosphorus (DOP) in the open oceans are often much higher than those of DIP (Björkman and Karl, 2003;Karl and Björkman, 2015).The primary sources of DOP in the ocean originate from biological processes that include exudation, viral lysis, autolysis and cell death, and grazing (Karl and Björkman, 2015).C-P bond phosphonates are broadly distributed in the ocean, with the nuclear magnetic resonance (NMR) spectra of ultrafiltered dissolved organic matter (DOM) revealing that phosphonates (21%) are the second predominant components, after phosphate esters (73%), of high molecular weight dissolved organic phosphorus (HMWDOP) in the Pacific Ocean, the Atlantic Ocean, and the North Sea (Clark et al., 1998;Kolowith et al., 2001;Repeta et al., 2016).At station ALOHA, MPn and its precursor 2-hydroxyethylphosphonate (2-HEP) account for approximately 20% of the HMWDOP (Repeta et al., 2016).MPn can be synthesized by microbes, such as the archaeon Nitrosopumilus maritimus (Metcalf et al., 2012) and by other abundant marine bacteria, such as Candidatus Pelagibacter sp.(Born et al., 2017).The catabolism of MPn involves active transport into the cytoplasm through the phosphonatespecific ABC transporter system integrated by the phnCDE complex (Stasi et al., 2019) and subsequent degradation by the protein complex phnGHIJK (Kamat et al., 2011;Amstrup et al., 2023).
Trichodesmium, a prominent N 2 -fixing organisms in the pelagic oceans, has long been recognized as a primary contributor to oceanic N 2 fixation (Capone et al., 1997;Bergman et al., 2013), supplying the "new" nitrogen within the euphotic zones of the tropical and subtropical regimes (Zehr and Capone, 2020).N 2 fixation and growth of Trichodesmium in the ocean is often constrained by the availability of inorganic phosphorus (Pi) (Sañudo-Wilhelmy et al., 2001;Frischkorn et al., 2018;Wang, 2022).Under Pi-limiting conditions, Trichodesmium exhibits weak competitiveness for Pi due to its lower specific affinities compared to other phytoplankton (Mccarthy and Carpenter, 1979;Sohm and Capone, 2006;Dyhrman, 2016).However, it can effectively utilize phosphate esters as a dependable source of phosphorus, supported by high levels of alkaline phosphatase (APase) activity (Mccarthy and Carpenter, 1979;Stihl et al., 2001;Sohm and Capone, 2006).Furthermore, Trichodesmium acquires phosphorus from phosphonate compounds, including MPn, through the C-P lyase pathway (Dyhrman et al., 2006), which leads to methane production as a byproduct (Karl et al., 2008).Hence, Trichodesmium is significant not only as a contributor to new nitrogen input for ocean but also as a potential source of ocean methane.However, little has been documented on the impacts of environmental drivers (e.g., temperature, P-availability) on methane production by Trichodesmium.
Trichodesmium's temperature tolerance ranges from18°C to 34°C when grown on the ingoranic phosphate, with optimal temperatures between 26-28°C (Breitbarth et al., 2007;Chappell and Webb, 2010;Fu et al., 2014).Rising temperatures impact diverse physiological and biochemical processes, including nitrogen fixation, respiration, carbon fixation, and growth (Breitbarth et al., 2007;Yvon-Durocher et al., 2012;Fu et al., 2014).As ocean warming intensifies, thermal stratification in the upper ocean stabilizes, reducing upward supply of dissolved inorganic phosphate (DIP) (Behrenfeld et al., 2006;Polovina et al., 2008;Somavilla et al., 2017;Ulloa et al., 2019).This may promote Trichodesmium to utilize more MPn, which may boost methane emission.Nevertheless, the thermal response of methane production by Trichodesmium has not been examined thus far.We hypothesize that changes in temperatures would alter levels of methane production by Trichodesmium grown on Mpn, since its utilization should correlate with the temperaturedependent assimilations of carbon, nitrogen, and phosphorus.To test this hypothesis, we employed MPn as the phosphorus source, grew Trichodesmium under various temperatures ranging from 16 to 31°C, and examined the correlations between methane production and several pivotal physiological processes, including assimilation of carbon, nitrogen, and phosphorus, as well as the specific growth rates.We found that changes in temperature affect the methane production of Trichodesmium grown on Mpn in parallel with the rates of C/N/P assimilations and growth.

Culture conditions
Trichodesmium IMS101 stock culture was grown in nitrogenfree YBCII medium (Chen et al., 1996) containing 5 μM methylphosphonate (MPn, Aladdin, CAS 993-13-5, ≥98%).The stock culture was kept at 27°C, and the light intensity was set at 110 μmol photons m −2 s −1 with a light-dark cycle of 12:12 h (light source: white LED tubes; light period: 08:00 to 20:00 local time).The cultures were run at five different temperatures (16, 20, 23, 27, 31°C), with all other environmental conditions maintained identical to those of the stock cultures.No other P source or dissolved inorganic phosphorus (DIP) was intentionally added to the cultures.The experimental cultures were kept in the exponential growth phase through regular dilutions performed every 3-15 days depending on the temperatures and growth rates.The semi-continuous cultures ensured that chlorophyll-a (Chl-a) concentrations consistently fell within the range of 0.005-0.05μg mL −1 .The experimental cultures were acclimated to their respective temperatures for more than 6 months prior to measuring physiological and biochemical parameters.

Chl-a and specific growth rate
Chl-a concentration was determined using the spectrophotometric method.The cells were filtered onto GF/F filters and subsequently extracted in pure methanol overnight at 4°C in the dark.After extraction, the samples were centrifuged at 12,000 g for 4 min.The resulting supernatant was then scanned for absorbance across the wavelength range of 250-800 nm using a spectrophotometer (Cary 60, Agilent, CA, United States).Chl-a concentration was calculated using the following formula (Ritchie, 2006): where OD 665 and OD 750 were absorbance at wavelengths 665 nm and 750 nm, respectively.The following equation was used to calculate the specific growth rate ( ) µ based on the Chl-a concentration: where m 2 and m 1 are the Chl-a values at time t 2 and t 1 , respectively.
The following equation modified by Norberg (2004) was utilized to fit the thermal growth curve of Trichodesmium: where the specific growth rate (μ) is a function of temperature (T).In this equation, the coefficient w represents the thermal niche width, whereas the explicit biological significance of coefficients a, b, and z remain unspecific.Collectively, these four coefficients can be used to derive both the maximum growth rate and the optimum growth temperature (T opt ):

Chl-a fluorescence
The effective photosynthetic quantum yield of photosystem II (Φ II ) and relative electron transport rate (rETR) were measured by Multi-color PAM (Multi-color PAM, Walz).Samples were acclimated to white actinic light with photon flux intensities similar to the growth conditions for 2 min to measure F s , and then F m ′ was measured using a saturating pulse (8,000 μmol photons m −2 s −1 , 0.8 s) to obtain the effective quantum yield (Φ II ) as follows (Genty et al., 1989): Subsequently, rapid light curves (RLC) were measured at 11 actinic light levels [E], from 5 to 2,904 μmol photons m −2 s −1 with each light exposure lasting for 30 s. Relative electron transport rates (rETR) were calculated as follows (Ralph and Gademann, 2005): The RLC was fitted by the following equation (Supplementary Figure S1) (Eilers and Peeters, 1988), where a, b, and c are the fitting coefficients.These three coefficients were used to derive the photosynthetic light-harvesting efficiency (α), maximum relative electron transport rate (rETR max ):

Methane production
A Cavity Ring-Down Spectroscopy gas analyzer (Picarro Model G2308, CA, United States) was used to measure methane concentrations.A 350 mL sample (V l ) was placed within a polycarbonate (PC) bottle, leaving a 270 mL headspace (V g ).The PC bottle was sealed using a silicone stopper to ensure an airtight condition.The silicone stopper was equipped with two three-way valves for sampling.During the incubation of the samples at each of the temperatures (16, 20, 23, 27, 31°C) for methane measurement, 200 mL of gas in the bottle was replaced with sterile air every 12 h.The methane content in the headspace gas was shaken to equilibrate the dissolved gas with the headspace before subsequent measurement using the Picarro analyzer according to Lenhart et al. (2016).The total methane production rates (b CH 4 ) were calculated based on the modified formula (Johnson et al., 1990): where C g indicates the methane concentration in the headspace (nmol L −1 ), C M represents the methane concentration in the headspace after replacing with the sterile air, C air denotes the methane concentration of the sterile air, and C g1 and C g2 signify the methane concentrations at two time points, t 1 and t 2 .Additionally, β is the Benson coefficient.V g represents the headspace volume, while V l corresponds to the volume of the culture medium.Methane production rates specific to Trichodesmium were determined based on the rates obtained from samples extracted from fractions containing heterotrophic bacteria.Several control experiments were Zou et al. 10.3389/fmicb.2024.1396369Frontiers in Microbiology 04 frontiersin.orgconducted to assess the potential influence of heterotrophic bacteria and to account for systematic errors.To test the bacterial contribution to methane production we filtered YBCII media through a 1.2 μm polycarbonate membrane, which lacked Trichodesmium but retained the heterotrophic bacteria in the cultures.It should be noted that the bacteria attached to the filaments of Trichodesmium might not be filtered off into the medium.

Carbon and nitrogen assimilation
In this study carbon and nitrogen assimilation rates were determined by assessing changes over time for particulate organic carbon (POC) and particulate organic nitrogen (PON) respectively.Briefly, samples for POC and PON measurements were taken at 0, 12, and 24 h after the start of the light period.The changes in particulate organic nitrogen (PON) and particulate organic carbon (POC) were then analyzed to determine the assimilation rate of carbon and nitrogen.Upon sampling, cells were filtered onto pre-combusted (450°C, 4 h) GF/F filters and rinsed with 100 mL fresh nitrogen-free YBCII.Subsequently, the filters were acidified for 24 h in HCl fumes and then dried for 24 h to remove unassimilated inorganic carbon.An elemental analyzer (Vario EL cube, Elementary, Germany) was used to quantify POC and PON.Changes in POC or PON between samples taken at 12 h and 0 h represented the assimilation during the light period, and changes between samples taken at 24 h and 0 h represented the daily assimilation.The results were comparable to the POC and PON production rates, which were calculated as POC or PON content (nmol Chl a −1 ) × specific growth rate μ (d −1 , Supplementary Figure S2) (Tong et al., 2019).
We acknowledge that the POC changes effectively represent integrated assimilation of inorganic carbon and recycling of organic carbon leaked into the media by Trichodesmium.Similarly, the nitrogen assimilation rates, as determined in this study, are indicative of the sum of both N 2 -fixation and the recycling of biogenic nitrogen leaked into the media by Trichodesmium.

Dissolved and particulate phosphorus measurement
To determine dissolved inorganic phosphorus (DIP), the samples were filtered through a 0.22 μm cellulose acetate membrane and were then analyzed using an auto-analyzer (AA3, Seal, Germany) at room temperature.For the measurement of particulate phosphorus (PP), we followed the Solórzano method (Solórzano and Sharp, 1980).In brief, the samples were filtered on pre-combusted (450°C, 4 h) GF/F filters (25 mm, Whatman, United States), and rinsed with 100 mL phosphorus-free YBCII artificial seawater.Subsequently, the filters were soaked with 0.017 M MgSO 4 , dried at 95°C, and then baked for 2 h at 450°C.Before measurement, the samples were hydrolyzed with acid (0.2 M HCl) at 80°C for 30 min.These procedures convert PP to DIP, which was subsequently quantified with an autoanalyzer (see above).The assimilation rates of phosphorus during a daily cycle were calculated using the PP contents at 0 h, 12 h, and 24 h.

Thermal dependence of metabolic processes
The thermal dependence of methane production and the related assimilation of carbon, nitrogen and phosphorus was analyzed using the Boltzmann-Arrhenius equation (Padfield et al., 2016): where b(T) represents the metabolic rate at temperature T (Kelvin, K), k Boltzmann's constant (8.62 × 10 −5 eV K −1 ), b(T c ) the rate of metabolism normalized to an arbitrary reference temperature, T c = 25°C, and E a is the activation energy (in electron volts, eV) for the metabolic process.Supra-optimal temperatures could deviate the metabolic rate from the Boltzmann-Arrhenius equation.Therefore, corresponding values, usually those observed at 31°C, were excluded from this analysis.
The sensitivity of methane production and assimilation of carbon, nitrogen and phosphorus to temperature changes (Q 10 ) from 16°C to 27°C was assessed by the following model (Van't Hoff and Lehfeldt, 1899): where Rate 1 and Rate 2 indicate metabolic rates at 16°C (T 1 ) and 27°C (T 2 ), respectively.

Statistical analysis
The data were provided as the means of three replicates (independent cultures) with standard deviation (SD) (n = 3).To examine the statistical differences between treatments, one-way ANOVA and Tukey's test were used.The Brown-Forsythe test and Shapiro-Wilk test were used to check data homoscedasticity and normality, respectively.

Specific growth rates, ratio of Chl-a to POC, ratio of POC to PON and Chl-a fluorescence
The specific growth rate of Trichodesmium increased with temperature (Figure 1; one-way ANOVA, p < 0.001) and was sensitive to temperature changes, with a Q 10 value (the rate increase fold for every 10-degree rise in the temperature) of growth rate for temperatures ranged from 16°C to 27°C reaching 8.6 ± 2.3.As the culture temperature increased from 16°C to 27°C, the specific growth rate exhibited a tenfold increase from 0.03 ± 0.01 d −1 to 0.34 ± 0.02 d −1 (Figure 1, Tukey's test, p = 0.0001).Increasing the growth temperature to 31°C did not significantly change the specific growth rate (0.32 ± 0.05 d −1 , Figure 1, Tukey's test, p = 0.999).According to the Norberg equation, the optimal growth temperature (T opt ) was 29.2°C, corresponding to a maximum growth rate of 0.38 d −1 .
The ratio of Chl-a to POC increased with the rise in temperatures (one-way ANOVA, p < 0.0001).The Chl a: POC was 68 ± 1 ng Chl a μmol POC −1 at the lowest growth temperature of 16°C, and increased by 100% to 136 ± 3 ng Chl a μmol POC −1 at 31°C (Table 1; Tukey's test, p < 0.0001).
In contrast, the ratios of POC to PON were negatively correlated with temperatures (Table 1; one-way ANOVA, p = 0.0004).When the growth temperature increased from 16°C to 31°C, the ratios of POC to PON decreased by 15% from 6.60 ± 0.24 to 5.64 ± 0.26 (Table 1; Tukey's test, p = 0.004).

Methane production
Methane was produced during both the day and the night, exhibiting a typical temperature response curve (Figure 2, one-way ANOVA, p < 0.0001) and demonstrated sensitivity to temperature changes, with a Q 10 value for temperature ranged from 16°C to 27°C, reaching 4.6 ± 0.7.At 16°C, the methane production rates during the light period and daily cycle were 0.120 ± 0.050 nmol CH 4 μm POC −1 day −1 and 0.198 ± 0.042 nmol CH 4 μmol POC −1 day −1 , respectively (Figure 2).When the growth temperature increased to 27°C, the rate of methane production reached the maximum, with light-period and daily-cycle values increasing by 407% to 0.608 ± 0.018 nmol CH 4 μmol POC −1 day −1 (Tukey's test, p < 0.0001) and by 419% to 1.028 ± 0.040 nmol CH 4 μmol POC −1 day −1 (Tukey's test, p < 0.0001), respectively (Figure 2).The activation energy (E a ) for methane production was 1.09 ± 0.13 eV for light period and 1.08 ± 0.08 eV for daily-cycle.Notably, the amount of methane production in each treatment increased with incubation time (Supplementary Figure S3A) and showed a positive correlation with biomass (Supplementary Figure S4).Methane production was also observed in bacterial controls, and the production rate increased with rising temperature (Supplementary Figure S3B).Trichodesmium dominated the methane production, with heterotrophic bacterial contribution to the total production being 19.7% at 16°C, 1.5% at 20°C, 3.8% at 23°C, 1.9% at 27°C, and 11.4% at 31°C, respectively (Supplementary Figure S3A).Methane production by Trichodesmium was derived by subtracting the bacterial production from the total production, with the possible contribution of the attached heterotrophic bacteria to the filaments of Trichodesmium being ignored.

Carbon, nitrogen, and phosphorus assimilation
Carbon, nitrogen, and phosphorus assimilation rates increased with higher growth temperatures (Figure 3, one-way ANOVA, p < 0.0001, p < 0.0001, and p < 0.0001), demonstrating a high sensitivity to temperature changes.The Q 10 values for temperature ranged from 16°C to 27°C for POC, PON and POP production were 4.9 ± 0.7, 2.6 ± 0.5, and 2.6 ± 0.3, respectively.When the temperature was raised from 16°C to 31°C, the light-period and daily rates of carbon assimilation increased by 577% from 81 ± 34 to 548 ± 128 nmol C nmol POC −1 day −1 (Tukey's test, p < 0.0001) and  The rETRmax and α were derived from the rapid light curves (Supplementary Figure S1).Values represent the means ± SD of independent biological replicates (n = 3).
The methane production correlated positively with rates of C, N, and P assimilation, specific growth and relative electron transfer (Figures 4, 5).The daily methane production increased with higher carbon, nitrogen, and phosphorus assimilation (Figure 4).The correlation coefficients between the daily methane production and the assimilation of carbon, nitrogen and phosphorus were 0.83, 0.74, and 0.77, respectively (Figure 4, p < 0.0001, p < 0.0001, p < 0.0001).The methane production quotients (MPQ), ratio of the carbon, nitrogen or phosphorus assimilation rates to the methane production rates, were 253-494 for carbon, 40-128 for nitrogen, and 1.3-3.4 for phosphorus, respectively (Table 2).
The relationship between metabolic activity and growth rates to methane production rates, which was established based on the positive correlations of methane production with the specific growth rate (μ) and relative electron transport rate (rETR), provided the correlation coefficients (R-square) of 0.91 and 0.87, respectively (Figure 5, p < 0.0001, p < 0.0001).Under different growth temperatures, the ratio of specific growth rate to methane production rate ranged from 0.17 to 0.41 and that of rETR to methane production rate ranged from 118 to 215 (Table 2).

Discussion
The thermal dependence of key metabolic processes, including N 2 fixation, photosynthetic CO 2 fixation and specific growth rate, have been widely explored in Trichodesmium (Mulholland and Bernhardt, 2005;Breitbarth et al., 2007;Boatman et al., 2017).These studies have contributed to our understanding of its current distribution in natural environments and provided insights to predicting its future behavior under influence of climate change (Breitbarth et al., 2007;Fu et al., 2014;Jiang et al., 2018;Yi et al., 2020).This study presents the first report on how temperature impacts methane production in Trichodesmium and the relationship The effect of growth temperature on methane production by Trichodesmium IMS101.The cultures were acclimated to the temperatures for 6 months, and the methane production rates were normalized to particulate organic carbon (POC).The hollow circles represent methane production during the daytime (12 h), whereas the solid ones represent the diel production, summing that of light and dark periods.Values represent the means ± SD of independent biological replicates (n = 3). of methane production to the assimilation of carbon, nitrogen and phosphorus, although it has been known to release methane for almost two decades (Karl et al., 2008;Beversdorf et al., 2010;Repeta et al., 2016).Within the experimental temperature range of 16 to 31°C, the daily methane production increased with growth temperature and saturated at 27°C, with a corresponding value of 1.028 ± 0.040 nmol CH 4 μmol POC −1 day −1 (Figure 2).Obvious linear positive correlations were detected between methane production and the assimilations of carbon, nitrogen and phosphorus (Figure 4).
Unlike the assimilation of carbon and nitrogen, which predominantly occurred during the light period (Figure 3), substantial methane production and phosphorus assimilation were observed during dark periods as well (Figure 2).In Trichodesmium, photosynthetic CO 2 fixation and biological N 2 fixation provide the substrates for other physiological and biochemical processes (Berman-Frank, 2001;Milligan, 2007).Moreover, the circadian clock of Trichodesmium confines the expression of genes related to CO 2 and N 2 fixation to light period (Ohki et al., 1992;Chen et al., 1996Chen et al., , 1998;;Dodd et al., 2005;Masotti et al., 2007;Haydon et al., 2013).For MPn uptake and C-P cleavage, the required energy is directly provided by ATP (Kamat et al., 2011;Stasi et al., 2019).During the night, mitochondrial respiration likely served as the primary source of ATP, which consumes the organic carbon stored during the daytime period (Figure 3A).Within the temperature range of 16-31°C, the specific growth, the methane production rate and the assimilation rates of carbon, nitrogen and phosphorus exhibited high sensitivity to thermal changes with Q10 greater than 2 (Supplementary Table S1).The methane production of Trichodesmium saturated at 27°C and decreased when the temperature exceeded 31°C, indicating a threshold for warming to promote methane production.Methane production in the associated heterotrophic bacteria (control) also increased with rising temperatures  Methane production rates of Trichodesmium IMS101 as a function of specific growth rates (A) and relative electron transport rates (rETR) (B).The rETR, were derived from the rapid light curves (Supplementary Figure S1), and represent rates measured at the acclimation light level.The correlation of methane production and the specific growth rates was based on Figures 1, 2 and peaked at 31°C Figure S3), suggesting that future ocean warming may promote the growth of these bacteria and potentially stimulate degradation of organic matters with a possibility to indirectly affect MPn-based methane production by the diazotroph.
While the specific growth rates of Trichodesmium grown with the MPn as the main phosphorus source (Figure 1) were comparable to that of Trichodesmium grown on inorganic phosphate (Breitbarth et al., 2007;Jiang et al., 2018;Yi et al., 2020), Trichodesmium grown on MPn demonstrated resilience to temperature as low as 16°C, contrasting sharply with the constraints observed in Trichodesmium cultured with inorganic phosphate.Specifically, Trichodesmium IMS101 did not grow below 20°C when supplied with inorganic phosphate (Breitbarth et al., 2007).In studies involving three other strains of Trichodesmium erythraeum (KO4-20, RLI, 21-75), the minimum growth temperature was reported to be 18°C (Fu et al., 2014).In the oceans, Trichodesmium is predominately found in tropical and subtropical waters with temperatures exceeding 20°C (Laroche and Breitbarth, 2005;Luo et al., 2012).Although its presence in higher latitudes with colder waters is occasionally reported (Díez et al., 2012;Rees et al., 2016;Sabeur et al., 2016), it had been thought that these Trichodesmium cells were merely transported by ocean currents and were unable to establish a sustained presence (Laroche and Breitbarth, 2005).Our results provide an additional interpretation for Trichodesmium's persistence in higher latitudes where temperatures drop below 20°C.When MPn was utilized as the primary phosphorus source, Trichodesmium demonstrated the capacity to acclimate and grow at temperatures as low as 16°C, possibly sustaining a seed community in colder regions (Shaika et al., 2023).However, the underlined mechanisms for the tolerance of the low temperature in Trichodesmium grown with Mpn need to be explored in future studies.
Based on the significant positive correlations between methane production and the assimilation of carbon, nitrogen, and phosphorus, as well as with the specific growth rate (μ) and relative electron transport rate (rETR) (Figures 4, 5), we established a series of methane production quotients (MPQ), representing the ratios between key metabolic rates and methane production rates (Table 2).Given that Trichodesmium's capacity for carbon and nitrogen fixation has been extensively studied over the past few decades (Sañudo-Wilhelmy et al., 2001;Mulholland et al., 2006;Bergman et al., 2013), we applied the methane production quotients for carbon and nitrogen to estimate Trichodesmium's counteractive roles in mitigating greenhouse effect.Based on the MPQ (Table 2), for every 227-494 nmol of CO 2 assimilated by Trichodesmium, 1 nmol of CH 4 is emitted.
Given that the global warming potential of CH 4 is about 80 times that of CO 2 (Collins et al., 2013), and taking into account that the contribution of MPn accounts for 5-9.8% (Clark et al., 1998;Repeta et al., 2016;Sosa et al., 2020) of the other phosphorus sources in the oceanic areas where Trichodesmium thrives, the CH 4 emission via MPn demethylation from Trichodesmium can offset 0.25-1.12% of its CO 2 mitigation effects.This calculation does not take the oxidation of methane during diffusion processes into account (Hakemian and Rosenzweig, 2007).As Trichodesmium inhabit the upper 40 meters and form extensive sea surface blooms, it is plausible to assume that methane produced by Trichodesmium will be released to the atmosphere.Furthermore, methane production in Trichodesmium is mainly coupled with the utilization of MPn.Investigating the concentrations of MPn in different waters would promote more accurate estimates on a global scale.
Under the influence of global warming, concurrent ocean warming (Collins et al., 2013) is predicted to shoal thermal stratification within the upper layers of the oceans (Somavilla et al., 2017;Frischkorn et al., 2018;Ulloa et al., 2019) and reduce the upward transport of dissolved inorganic phosphate (DIP).Consequently, Trichodesmium may increasingly resort to utilizing MPn, potentially leading to enhanced methane production.Additional drivers may impact methane production.The optimal growth temperature changes in Fe-replete versus Fe-deplete cells (Jiang et al., 2018).The Fe and MPn availability may further affect methane production and requires future study.
TABLE 2 The quotients of carbon (C as : CH 4pro ), nitrogen (N as : CH 4pro ) and phosphorus (P as : CH 4pro ) assimilations to methane production rate, and the ratio of specific growth rate (μ: CH 4pro ) and relative electron transfer rate (rETR: CH 4pro ) to methane production rate in Trichodesmium IMS101 after having acclimated for 6 months to specific temperatures ranging from 16 to 31°C.Values represent the means ± SD of independent biological replicates (n = 3).Zou et al. 10.3389/fmicb.2024.1396369Frontiers in Microbiology frontiersin.org The author(s) declare that financial support was received for the research, authorship, and/or publication of this article.This study was supported by National Natural Science Foundation of China (No. 42361144840) to KG and the joint NSFC-Israel Science Foundation (ISF) grant No. 3051/23 to IB-F.

FIGURE 1
FIGURE 1 Specific growth rates of Trichodesmium IMS101 during the exponential stage after having acclimated at each of the specific temperatures [ranging from 16 to 31°C for 6 months].The dotted line represents the 95% prediction bands of the best-fit line.Values represent the means ± SD of independent biological replicates (n = 3).

FIGURE 3
FIGURE 3Carbon, nitrogen, and phosphorus assimilation rates (A-C) of Trichodesmuim IMS 101 normalized to particular organic carbon (POC) during the light period (hollow circle) or the diel cycle (solid circle) after being acclimated for 6 months to specific temperatures.Values represent the means ± SD of independent biological replicates (n = 3).

FIGURE 4
FIGURE 4Methane production rates of Trichodesmium IMS101 correlated with assimilation rates of carbon (A), nitrogen (B), and phosphorus (C), which were derived from Figures 2, 3.The hollow symbols and dotted line represent the correlations for the light period, while the solid symbols and lines represent the daily values.The numbers next to the different symbols indicate the growth temperatures (°C).Each data point represents an independent biological replicate, and different symbols represent different biological replicates.

FIGURE 5
FIGURE 5 . The numbers next to the symbols indicate the growth-acclimation temperatures (°C).Each data point represents an independent biological replicate, and different symbols represent different biological replicates.

TABLE 1
Chlorophyll a content normalized to cellular particular organic carbon (Chl a: POC), ratio of POC to PON (POC: PON), effective quantum yield (Φ II ), photosynthetic light-gathering efficiency (α), maximum relative electron transport rate (rETR max ) of Trichodesmium IMS101 after having acclimated for 6 months to specific temperatures ranging from 16 to 31°C.